anchor.diff:  tRatio = 1        ! [-]  ratio of bounding ply thickness to main ply
   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
   
   EL = ! [MPa] UD longitudinal Young's modulus
   ET = ! [MPa] UD transverse Young's modulus
  
positive.code: 
  tBPly = tRatio*(2*t) ! [mm] thickness of the bounding ply
  tTotal = t + tBPly   ! [mm] thickness of the bounding ply
  
  ! Create Geometry
  
  
negative.code: MSHKEY, 0
  AMESH, 3, 4, 1
  
  FINISH              ! Finish pre-processing
  
  /SOLU               ! Enter the solution processor
  