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: GL = ! [MPa] UD in-plane shear modulus
  GT = ! [MPa] UD transverse shear modulus
  
  L = t/rhon        ! [mm] length of the RVE
  a = atRatio*t     ! [mm] 2a = crack length
  
  
negative.code: 
  /SOLU               ! Enter the solution processor
  
  ANTYPE,0            ! Analysis type,static
  
  
  ! Define Displacement Constraints on Lines   (dl command)
  ! DL, LINE, AREA, Lab, Value1, Value2
  DL, 1, ,SYMM
  DL, 2, ,SYMM
  DL, 6, ,SYMM
  DL, 7, ,SYMM
  DL, 3, ,UX,appliedDisp
  DL, 4, ,UY,appliedDisp
  
  SOLVE                ! Solve the problem
  
  FINISH               ! Finish the solution processor
  
  SAVE                 ! Save your work to the database
  