*HEADING
PDEBCRGRNLG: DEBONDING OF AN INTERFACE VIA CRACK LENGTH VS TIME CRITERION
*RESTART,WRITE,FREQ=1
*NODE
1,
121,12.,
2001,,10.
2121,12.,10.
2,0.,
122,12.
999, 12.0, 0.0
*NGEN,NSET=BOT
1,121,10
*NGEN,NSET=BOT2
11,121,10
*NGEN,NSET=TOP
2001,2121,10
*NFILL
BOT,TOP,4,500
*NSET, NSET=TERM2,GEN
41,101,20
*ELEMENT,TYPE=CPE4
1,1,21,1021,1001,11,521,1011,501
*ELGEN,ELSET=SOLID
1,6,20,20,2,1000,1000
*SOLID SECTION,ELSET=SOLID,MAT=ELASTIC
1., 
*MATERIAL,NAME=ELASTIC
*ELASTIC
1.E6,.3
**
*** Contact definitions.
**
*rigid body,analytical surface=asurf, REF NODE=999
*surface,TYPE=SEGMENTS,NAME=ASURF
START,-1.0, 0.0
LINE, 10.1, 0.0
*ELSET, ELSET=BOT
1,21,41,61,81
*surface, NAME=BSURF
BOT, S1
*SURFACE INTERACTION, NAME=BOTH
1.0, 
*CONTACT PAIR,INTERACTION=BOTH
BSURF, ASURF
*NSET,NSET=N41
41,
*NSET,NSET=N61
61,
*NSET,NSET=N81
81,
*NSET,NSET=N101
101,
*NSET,NSET=NREF
1,
***
**** Nodeset defining the bonded portion of the surface
***
*NSET, NSET=BOND
41,61,81,101
*INITIAL CONDITIONS, TYPE=CONTACT, NORMAL
BSURF,ASURF,BOND
**
***
*STEP,INC=20, NLGEOM
*STATIC
8.0E-04,0.00287
*BOUNDARY
2001,2,2,1.
2021,2,2,1.
121,1,2
999, 1, 2
999, 6
*DEBOND,SLAVE=BSURF,MASTER=ASURF,TIME INC=.0002, FREQUENCY=1,OUTPUT=BOTH
0.,1.
.0002,0.,
*FRACTURE CRITERION,TYPE=CRACK LENGTH,TOLERANCE=0.01,NSET=NREF
1.0106E-03, 4.0
1.2106E-03, 6.0
1.7106E-03,  8.0
1.8797E-03, 10.1
*ENERGY PRINT
*EL PRINT
S,
E,
*CONTACT PRINT
CSTRESS,CDISP
DBT,DBSF,DBS
*CONTACT FILE
CSTRESS,CDISP
DBT,DBSF,DBS
*EL FILE,ELSET=SOLID
S, E
*EL PRINT,ELSET=SOLID
S,E
*CONTOUR INTEGRAL,SYMM,CONT=2,OUTPUT=BOTH
N41,1.0,
*CONTOUR INTEGRAL,SYMM,CONT=2,OUTPUT=BOTH
N61,1.0,
*CONTOUR INTEGRAL,SYMM,CONT=2,OUTPUT=BOTH
N81,1.0,
*CONTOUR INTEGRAL,SYMM,CONT=2,OUTPUT=BOTH
N101,1.0,
*END STEP