anchor.diff:  
   ! ===> START INPUT DATA
   
  -Vf = 0.000079! [-] Fiber volume fraction
  -
  +L = 10.0          ! [um]  length of the RVE
  +W = 10.0          ! [um]  width of the RVE
   Rf = 1.0          ! [um]  fiber's radius
   theta0 = 0.0      ! [deg] angular position of debond
   deltatheta = 10.0 ! [deg] 2*deltatheta = debond size
   deltaphi = 5.0    ! [deg] refined area size
   delta = 0.05      ! [deg] angular size of element at crack
  +aOverW = 0.1      ! [-] ratio of debond length to RVE's width
   
   epsx = 0.01       ! [-]  applied strain
   uniP = 0.0        ! [-]  uniform pressure applied to crack face
  
positive.code: 
  ! ===> END INPUT DATA
  
  L = 2*t/rhon        ! [mm] length of the RVE
  
  appliedDisp = epsx*L ! [mm] applied displacement
  
  
negative.code: 
  ! Points
  
  K, 1, -L, W, 0.0   ! SW corner, external face with debond
  K, 2, L, W, 0.0    ! SE corner, external face with debond
  K, 3, L, W, L      ! NE corner, external face with debond
  K, 4, -L, W, L     ! NW corner, external face with debond
  
  K, 5, -L, 0.0, 0.0   ! SW corner, external face without debond
  K, 6, L, 0.0, 0.0    ! SE corner, external face without debond
  K, 7, L, 0.0, L      ! NE corner, external face without debond
  K, 8, -L, 0.0, L     ! NW corner, external face without debond
  
  
  
  ! Lines
  
  L, 1, 2            !1
  
  
  ! Areas
  
  