anchor.diff:  tTotal = t + tBPly   ! [mm] thickness of the bounding ply
   
   elSize = daOvera*a ! [mm] size of element in refined region close to crack tip
  +
  +appliedDisp = epsx*L ! [mm] applied displacement
   
   ! Create Geometry
   
  
positive.code: 
  ! Parameters
  
  Vf = ! [-] Fiber volume fraction
  
  t = 1             ! [mm] 2t = thickness of the element
  
negative.code: atRatio = 0.1     ! [-]  ratio of crack length to main ply thickness
  rhon = 0.01       ! [-]  normalized crack density
  daOvera = 0.05    ! [-]  ratio of crack increment (i.e. crack tip element size) to crack length
  epsx = 0.01       ! [-]  applied strain
  
  EL = ! [MPa] UD longitudinal Young's modulus
  ET = ! [MPa] UD transverse Young's modulus
  