*HEADING
   3D Single leg bend Delamination model 
   VCCT based low cycle fatigue
**
*parameter
** User input data defining the interface properties
** Fracture toughness
 GIc   = 2.0
 GIIc  = 2.0
 GIIIc = 30.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
**5.0e-5
tol=0.4
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
20001, 180,100,100, 10, 1,1 
*ELGEN,ELSET=upper
20051, 180,100,100, 10, 1,1
**
**
*Shell Section, elset=upper, material=MAT_upper,Orientation=Or1,OFFSET=-0.5
0.066,
** 
*Shell Section, elset=lower, material=MAT_lower,Orientation=Or1,OFFSET=0.5
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=secondary,generate
20001,37901,100
20002,37902,100
20003,37903,100
20004,37904,100
20005,37905,100
20006,37906,100
20007,37907,100
20008,37908,100
20009,37909,100
20010,37910,100
*surface,type=element,name=secondary
secondary,SPOS
*surface,type=element,name=main
upper,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
secondary, main,ntied
*Contact Pair, interaction=SURFS, small sliding
secondary, main
*Surface Interaction, name=SURFS
1.,
*Friction, slip tolerance=0.005
0.,
*clearance,secondary=secondary,main=main,value=3.0E-9
*MPC
BEAM, 12051, 12056
BEAM, 12052, 12056
BEAM, 12053, 12056
BEAM, 12054, 12056
BEAM, 12055, 12056
BEAM, 12057, 12056
BEAM, 12058, 12056
BEAM, 12059, 12056
BEAM, 12060, 12056
BEAM, 12061, 12056
**
**********************************************************************************
** Boundary constraints
**
*NSET, NSET=RL, GENERATE
 6051, 6061, 1
*NSET, NSET=RR, GENERATE
18051, 18061, 1
*NSET, NSET=NOSLIDE
 6051, 18051
*NSET, NSET=NOSLIP
 6051
*NSET, NSET=LP
 12056
**
*BOUNDARY
 RL, 3
 RR, 3
 NOSLIDE, 2
 NOSLIP, 1
*amplitude,name=test
0,0,0.5,1,1,0
** STEP: Step-1
**
*nset,nset=tip
9101,9102,9103,9104,9105,9106,9107,9108,
9109,9110,9111
*TIME POINTS, NAME=TimePointName,GEN
0,1,0.25
*STEP, INC=1000,nlgeom,unsymm=yes
*STATIC
0.1,1
**contact controls,automatic tolerances
*debond, secondary=secondary, main=main
*FRACTURE CRITERION,TYPE=VCCT,MIXED MODE BEHAVIOR=REEDER,TOLERANCE=<tol>
<GIc>, <GIIc>, <GIIIc>, <eta>
*BOUNDARY,amp=test
LP, 3,3, -0.035
*output, field,freq=5
*element output
s
**s,e
*node output
u,
*contact output,secondary=secondary,main=main
dbt
*Output, history,freq=5
*NODE output,NSET=LP
 U3, RF3
*End Step
*STEP, INC=1000,nlgeom,unsymm=yes
*fatigue,TIME POINTS=TimePointName
0.1,1,1E-9,0.25
3,6,400000,1.1,0.1
**contact controls,stabilize,SECONDARY=secondary,MAIN=main,stabilize=5E-6
**contact controls,automatic tolerances
*debond, secondary=secondary, main=main,debonding force=ramp
*FRACTURE CRITERION,TYPE=fatigue,MIXED MODE BEHAVIOR=Tabular
-6.46186,-1, 0,0
-4.16186, 1,0,0
-6.46186,-1, 0.5,0
-4.16186, 1,0.5,0
-6.46186, -1,1.0,0
-4.16186,1, 1.0,0
*BOUNDARY,amp=test
LP, 3,3, -0.035
*output, field,TIME POINTS=TimePointName
*node output
u
*element output
s
*contact output,secondary=secondary,main=main
bdstat, enrrt,cycle
*Output, history,freq=1
*NODE output,NSET=LP
 U3, RF3
*End Step