*HEADING
   2D Single leg bend Delamination model 
   tf=0.8 ply ts=24 ply (4% zpins not)
   use vcct based low cycle fatigue
*parameter
** Fracture toughness:
 GIc   = 5.63
 GIIc  = 9.305
 GIIIc = 0.0
** Spatially varying properties
 GIc1=1.505
 GIc2=1.6686
 GIc3=1.8279
 GIc4=2.14222
 GIc5=2.4436	
 GIc6=2.732
 GIc7=3.011
 GIc8=3.2744
 GIc9=3.524
 GIc10=3.7607
 GIc11=3.98443
 GIc12=4.1952
 GIc13=4.39308	
 GIc14=4.57799
 GIc15=4.74997
 GIc16=4.909
 GIc17=5.0551
 GIc18=5.188
 GIc19=5.30848
 GIc20=5.4157
 GIc21=5.51
 GIc22=5.59	
 GIc23=5.63
** B-K parameter:
 modeMixLaw=1
 eta=1.75
** mixed mode parameter
 am = 0.0
 an=0.0
 ao=0.0
damv=0
**5.0E-5
tol=0.2 
** width in the plane strain direction
 width =1.0
** All Nodes
*NODE
100000, 0.00, 0.0000, 0.0 
100024, 0.00, 0.1992, 0.0
200024, 0.00, 0.1992, 0.0
200048, 0.00, 0.3984, 0.0
**
160000, 6.00, 0.0000, 0.0 
160024, 6.00, 0.1992, 0.0
260024, 6.00, 0.1992, 0.0
260048, 6.00, 0.3984, 0.0
**
*NGEN,NSET=FL
100000, 100024, 1
*NGEN,NSET=FR
160000, 160024, 1
*NFILL,NSET=NFLANGE
FL,FR,200,300
**
*NGEN,NSET=SL
200024, 200048, 1
*NGEN,NSET=SR
260024, 260048, 1
*NFILL,NSET=NSKIN
SL,SR,200,300
**
** FLANGE ELEMENTS
*ELEMENT, TYPE=CPE4I
700017,100016, 100316, 100317, 100017, 
*ELGEN
700017,8,1,1
**
*ELGEN, ELSET=EFLANGE_0
**700017,160,300,300
700018,160,300,300
**700019,160,300,300
700020,160,300,300
700021,160,300,300
**700022,160,300,300
700023,160,300,300
**700024,160,300,300
**
*ELGEN, ELSET=EFLANGE_45
700017,160,300,300
**700018,160,300,300
700019,160,300,300
**700020,160,300,300
**700021,160,300,300
700022,160,300,300
**700023,160,300,300
700024,160,300,300
**
** SKIN ELEMENTS
*ELEMENT, TYPE=CPE4I
800025,200024, 200324, 200325, 200025, 
*ELGEN
800025, 24,1,1
**
*ELGEN, ELSET=ESKIN_0
**800025,200,300,300
800026,200,300,300
**800027,200,300,300
800028,200,300,300
**800029,200,300,300
800030,200,300,300
**800031,200,300,300
800032,200,300,300
**800033,200,300,300
800034,200,300,300
**800035,200,300,300
800036,200,300,300
800037,200,300,300
**800038,200,300,300
800039,200,300,300
**800040,200,300,300
800041,200,300,300
**800042,200,300,300
800043,200,300,300
**800044,200,300,300
800045,200,300,300
**800046,200,300,300
800047,200,300,300
**800048,200,300,300
**
*ELGEN, ELSET=ESKIN_45
800025,200,300,300
**800026,200,300,300
800027,200,300,300
**800028,200,300,300
800029,200,300,300
**800030,200,300,300
800031,200,300,300
**800032,200,300,300
800033,200,300,300
**800034,200,300,300
800035,200,300,300
**800036,200,300,300
**800037,200,300,300
800038,200,300,300
**800039,200,300,300
800040,200,300,300
**800041,200,300,300
800042,200,300,300
**800043,200,300,300
800044,200,300,300
**800045,200,300,300
800046,200,300,300
**800047,200,300,300
800048,200,300,300
**
*BOUNDARY
200024, 1
200024, 2
260024, 2
**
*ORIENTATION, NAME=OID1
 1,0,0, 0,1,0
*SOLID SECTION,ELSET=ESKIN_0, MATERIAL=PW8552_0, ORIENTATION=OID1
*SOLID SECTION,ELSET=ESKIN_45, MATERIAL=PW8552_0, ORIENTATION=OID1
*SOLID SECTION,ELSET=EFLANGE_0, MATERIAL=PW8552_0, ORIENTATION=OID1
*SOLID SECTION,ELSET=EFLANGE_45, MATERIAL=PW8552_0, ORIENTATION=OID1
**
** TYPE=ENGINEERING CONSTANTS
** E1, E2, E3, Nu12, Nu13, Nu23, G12, G13
** G23
**
*MATERIAL,NAME=PW8552_0
*ELASTIC
  9100000.0,0.33
**
** Fracture elements (Rcurve) and MPC TIE included
*MPC
TIE, 100024, 200024
TIE, 100324, 200324
TIE, 100624, 200624
TIE, 100924, 200924
TIE, 101224, 201224
TIE, 101524, 201524
TIE, 101824, 201824
TIE, 102124, 202124
TIE, 102424, 202424
TIE, 102724, 202724
TIE, 103024, 203024
TIE, 103324, 203324
TIE, 103624, 203624
TIE, 103924, 203924
TIE, 104224, 204224
TIE, 104524, 204524
TIE, 104824, 204824
TIE, 105124, 205124
TIE, 105424, 205424
TIE, 105724, 205724
TIE, 106024, 206024
TIE, 106324, 206324
TIE, 106624, 206624
TIE, 106924, 206924
TIE, 107224, 207224
TIE, 107524, 207524
TIE, 107824, 207824
TIE, 108124, 208124
TIE, 108424, 208424
TIE, 108724, 208724
TIE, 109024, 209024
TIE, 109324, 209324
TIE, 109624, 209624
TIE, 109924, 209924
TIE, 110224, 210224
TIE, 110524, 210524
TIE, 110824, 210824
TIE, 111124, 211124
TIE, 111424, 211424
TIE, 111724, 211724
TIE, 112024, 212024
TIE, 112324, 212324
TIE, 112624, 212624
TIE, 112924, 212924
TIE, 113224, 213224
TIE, 113524, 213524
TIE, 113824, 213824
TIE, 114124, 214124
TIE, 114424, 214424
TIE, 114724, 214724
TIE, 115024, 215024
TIE, 115324, 215324
TIE, 115624, 215624
TIE, 115924, 215924
TIE, 116224, 216224
TIE, 116524, 216524
TIE, 116824, 216824
TIE, 117124, 217124
TIE, 117424, 217424
TIE, 117724, 217724
TIE, 118024, 218024
TIE, 118324, 218324
TIE, 118624, 218624
TIE, 118924, 218924
TIE, 119224, 219224
TIE, 119524, 219524
TIE, 119824, 219824
**
*nset, nset=bond, generate
120124, 145324, 300
*elset, elset=top, generate
820125, 859725, 300
*elset, elset=bot, generate
720124, 747724, 300
*surface, type=element, name=topsurf
top, s1
*surface, type=element, name=botsurf
bot, s3
*contact pair, interaction=crack, adjust=bond,small sliding
botsurf, topsurf
*surface interaction, name=crack
1.0
*initial conditions, type=contact
botsurf, topsurf, bond
**
*NSET, NSET=LP
 230024
*nset,nset=L145324
145324
*NODAL ENERGY RATE
bond,<GIc>, <GIIc>,<GIIIc>
145324,<GIc1>,<GIIc>,<GIIIc>
145024,<GIc2>,<GIIc>,<GIIIc>
144724,<GIc3>,<GIIc>,<GIIIc>
144424,<GIc4>,<GIIc>,<GIIIc>
144124,<GIc5>,<GIIc>,<GIIIc>
143824,<GIc6>,<GIIc>,<GIIIc>
143524,<GIc7>,<GIIc>,<GIIIc>
143224,<GIc8>,<GIIc>,<GIIIc>
142924,<GIc9>,<GIIc>,<GIIIc>
142624,<GIc10>,<GIIc>,<GIIIc>
142324,<GIc11>,<GIIc>,<GIIIc>
142024,<GIc12>,<GIIc>,<GIIIc>
141724,<GIc13>,<GIIc>,<GIIIc>
141424,<GIc14>,<GIIc>,<GIIIc>
141124,<GIc15>,<GIIc>,<GIIIc>
140824,<GIc16>,<GIIc>,<GIIIc>
140524,<GIc17>,<GIIc>,<GIIIc>
140224,<GIc18>,<GIIc>,<GIIIc>
139924,<GIc19>,<GIIc>,<GIIIc>
139624,<GIc20>,<GIIc>,<GIIIc>
139324,<GIc21>,<GIIc>,<GIIIc>
139024,<GIc22>,<GIIc>,<GIIIc>
138724,<GIc23>,<GIIc>,<GIIIc>
****** Step #1 *************
*amplitude,name=test
0,0,0.5,1,1,0
*amplitude,name=test2
0,0,0.5,1,1,0.0001
*STEP, NLGEOM, INC=5000
*STATIC,direct
0.01, 1,1.e-12
**
*controls, analysis=discontinuous
*controls, parameter=time incrementation
 , , , , , , , 30,
*CONTROLS, PARAMETER=FIELD,FIELD=DISPLACEMENT
,1.0
**contact controls,automatic tolerances
*BOUNDARY,amp=test2
230024,2,2,-0.12
*debond, slave=botsurf, master=topsurf
*FRACTURE CRITERION,TYPE=VCCT,NODAL ENERGY RATE
, , , <eta>
*output, field,freq=1
*element output
s
*node output
u,
*contact output,SLAVE=botsurf, MASTER=topsurf
cstress
cdisp
*output, history,freq=5
*node output, nset=LP
u2,
rf2
*END STEP
*STEP,INC=5000
*direct cyclic,fatigue
0.01,1,,,25,25,,20
3,6,20000,1.1
*debond, slave=botsurf, master=topsurf
*FRACTURE CRITERION,TYPE=fatigue,NODAL ENERGY RATE
0.5,-1.44,4.8768E-6,1.15,,0.95,,,
,1.75
*controls,type=direct cyclic
,100,1E-3,100,1E-3
**contact controls,automatic tolerances
*BOUNDARY,amp=test
230024,2,2,-0.12
*output,field, freq=50
*element output
s
*node output
u
*contact output,SLAVE=botsurf,MASTER=topsurf
bdstat
*output, history,freq=5
*node output,nset=LP
u2
rf2
*END STEP