anchor.diff:  
   ! Parameters
   
  +! ===> START INPUT DATA
  +
   Vf = ! [-] Fiber volume fraction
   
   t = 1             ! [mm] 2t = thickness of the element
  
positive.code: 
  /SOLU               ! Enter the solution processor
  
  ANTYPE,0			! Analysis type,static
  
negative.code: GL = ! [MPa] UD in-plane shear modulus
  GT = ! [MPa] UD transverse shear modulus
  
  ! ===> END INPUT DATA
  
  L = t/rhon        ! [mm] length of the RVE
  a = atRatio*t     ! [mm] 2a = crack length
  
  