*HEADING
   3D Double cantilever beam test model 
   - use vcct based low cycle fatigue to model crack growth
   orthogonal mesh, short run  
**
*parameter
** User input data defining the interface properties
** Fracture toughness (N/mm):
 GIc   = 1.2
 GIIc  = 6.0
 GIIIc = 6.0
** B-K parameter:
 eta = 1.75
 modeMixLaw=3
** mixed mode parameter
 am = 1.0
 an=1.5
 ao=1.0
** Damage and tolerance parameters
damv=0
tol=0.1
width=12.7 
*NODE
  1,   0.0, 0.0
 11,   0.0, 1.0
 51,   0.0, 0.0
 61,   0.0, 1.0
** 
**
 9001,   4.5, 0.0
 9011,   4.5, 1.0
 9051,   4.5, 0.0
 9061,   4.5, 1.0
**
18001,   9.0, 0.0
18011,   9.0, 1.0
18051,   9.0, 0.0
18061,   9.0, 1.0
**
*NGEN,NSET=NS1
    1 , 9001, 100
   51 , 9051, 100
 9001 ,18001, 100
 9051 ,18051, 100
*NGEN,NSET=NS2
   11 , 9011, 100
   61 , 9061, 100
 9011 ,18011, 100
 9061 ,18061, 100
**
*NFILL
NS1,NS2,10,1
**
**
*element,type=S4
20001, 1, 101, 102, 2
20051,51, 151, 152,52
**
*ELGEN,ELSET=lower_elements
20001, 180,100,100, 10, 1,1 
*ELGEN,ELSET=upper_elements
20051, 180,100,100, 10, 1,1
**
**
*Shell Section, elset=upper_elements, material=MAT_upper,Orientation=Or1
0.2,
** 
*Shell Section, elset=lower_elements, material=MAT_lower,Orientation=Or1
0.2,
** 
*Orientation,Name=Or1
1.0, 0.0, 0.0, 0.0, 1.0, 0.0
3, 0.0 
**
*MATERIAL,NAME=MAT_upper
*ELASTIC,TYPE=ENGINEERING CONSTANTS
8000000.0,8000000.0,8000000.0, 0.0, 0.0, 0.0,4000000.0,4000000.0,
4000000.0, 
**
*MATERIAL,NAME=MAT_lower
*ELASTIC,TYPE=ENGINEERING CONSTANTS
8000000.0,8000000.0,8000000.0, 0.0, 0.0, 0.0,4000000.0,4000000.0,
4000000.0, 
*elset,elset=slave,generate
20101,37901,100
20102,37902,100
20103,37903,100
20104,37904,100
20105,37905,100
20106,37906,100
20107,37907,100
20108,37908,100
20109,37909,100
20110,37910,100
*surface,type=element,name=slave
slave,SPOS
*surface,type=element,name=master
upper_elements,SNEG
*ELSET,ELSET=ETIED,GENERATE
29101,37901,100
29102,37902,100
29103,37903,100
29104,37904,100
29105,37905,100
29106,37906,100
29107,37907,100
29108,37908,100
29109,37909,100
29110,37910,100
*NSET,NSET=NTIED,ELSET=ETIED
*initial conditions, type=contact
slave, master,ntied
*Contact Pair, interaction=SURFS, small sliding, type=NODE TO SURFACE, no thickness
slave, master
*Surface Interaction, name=SURFS
1.,
*Friction, slip tolerance=0.005
0.,
*clearance,slave=slave,master=master,value=1.0E-8
*NSET,NSET=FIX,GENERATE
**18001, 18011
18051, 18061
*BOUNDARY  
FIX,1,6
**
**
*MPC
BEAM, 1, 6
BEAM, 2, 6
BEAM, 3, 6
BEAM, 4, 6
BEAM, 5, 6
BEAM, 7, 6
BEAM, 8, 6
BEAM, 9, 6
BEAM,10, 6
BEAM,11, 6
**
BEAM, 51, 56
BEAM, 52, 56
BEAM, 53, 56
BEAM, 54, 56
BEAM, 55, 56
BEAM, 57, 56
BEAM, 58, 56
BEAM, 59, 56
BEAM, 60, 56
BEAM, 61, 56
**
**
*NSET, NSET=LP
 6, 56
**
***************************************************************
** STEP: Step-1
**
*amplitude,name=test
0,0,0.5,1,1,0
*amplitude,name=test2
0,0,0.5,1,1,0.0001
*STEP, INC=1000
*STATIC,direct
0.02,1
*debond, slave=slave, master=master
*FRACTURE CRITERION,TYPE=VCCT,MIXED MODE BEHAVIOR=REEDER,TOLERANCE=<tol>
<GIc>, <GIIc>, <GIIIc>, <eta>
*BOUNDARY,amp=test2
**
 6,3,3, -0.05
56,3,3,  0.05
**
**
** OUTPUT REQUESTS
** 
**contact controls,automatic tolerances
*output, field,freq=10
*node output
u
*Output, history,freq=10
*NODE output,NSET=LP
 U3, RF3
*End Step
*STEP, INC=1000
*direct cyclic,fatigue
0.01,1,,,49,49,,20
3,6,400000,1.1
*debond, slave=slave, master=master
*FRACTURE CRITERION,TYPE=fatigue,MIXED MODE BEHAVIOR=REEDER
0.5,-1.44,4.8768E-6,1.15,,,1.2,6.0,
6.0,1.75
*controls,type=direct cyclic
,100,1E-3,100,1E-3
**contact controls,automatic tolerances
*BOUNDARY,amp=test
**
 6,3,3, -0.05
56,3,3,  0.05
**
**
*output, field,freq=100
*node output
u
*contact output,slave=slave,master=master
bdstat, enrrt
*Output, history,freq=50
*NODE output,NSET=LP
 U3, RF3
*End Step