*HEADING
: DEBONDING OF AN INTERFACE VIA COD CRITERION
** complete model.
*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, NSET=NTOP
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_ELEMENTS-1
1,6,20,20,2,1000,1000
**
***
*NCOPY, OLD SET=NTOP, NEW SET=NBOT,REFLECT=LINE,CHANGE NUMBER=3000
-1.0,0.0,0.0, 15.0,0.0,0.0
*ELCOPY,ELEMENT SHIFT=2000,OLD SET=SOLID_ELEMENTS-1, NEW SET=SOLID_ELEMENTS-2,
SHIFT NODES=3000,REFLECT
** element redefined to get correct node definitions
*ELEMENT, TYPE=CPE4, ELSET=SOL2
2101,3101,4101,4121,121
*ELSET, ELSET=SOLIDS
SOLID_ELEMENTS-1,SOLID_ELEMENTS-2
*SOLID SECTION,ELSET=SOLIDS,MAT=ELASTIC
1., 
*MATERIAL,NAME=ELASTIC
*ELASTIC
1.E6,.3
**
*** Contact definitions
**
*ELSET, ELSET=TOP
1,21,41,61,81
*surface, NAME=B-SURFACE
TOP, S1
*ELSET, ELSET=BOT
2001,2021,2041,2061,2081
*surface, NAME=A-SURFACE
BOT, S4
*SURFACE INTERACTION, NAME=BOTH
1.0, 
*CONTACT PAIR,INTERACTION=BOTH,SMALL SLIDING
B-SURFACE, A-SURFACE
*NSET,NSET=NODE_NUM_41
41,3041
*NSET,NSET=NODE_NUM_61
61,3061
*NSET,NSET=NODE_NUM_81
81,3081
*NSET,NSET=NODE_NUM_101
101,3101
*NSET,NSET=NREF
1,
***
**** Nodeset defining the bonded portion of the surface
***
*NSET, NSET=BOND
41,61,81,101
*INITIAL CONDITIONS, TYPE=CONTACT
B-SURFACE,A-SURFACE,BOND
**
***
*STEP,INC=10
*STATIC
8.0E-04,0.00287
*BOUNDARY
2001,2,2,1.
2021,2,2,1.
5001,2,2,-1.
5021,2,2,-1.
121,1,2
*DEBOND,SLAVE=B-SURFACE,MASTER=A-SURFACE,TIME INC=.0002, FREQUENCY=1,OUTPUT=BOTH
0.,1.
.0002,0.,
*FRACTURE CRITERION,TYPE=COD,TOLERANCE=0.05,DISTANCE=2.0
0.3499, 2.0
0.3654, 4.0
0.4235, 6.0
0.3142, 8.0
*ENERGY PRINT
*EL PRINT
S,
E,
*CONTACT PRINT
CSTRESS,CDISP
DBT,DBSF,DBS
*CONTACT FILE
CSTRESS,CDISP
DBT,DBSF,DBS
*EL FILE,ELSET=SOLIDS
S, E
*EL PRINT,ELSET=SOLIDS
S,E
*CONTOUR INTEGRAL,CONT=2,OUTPUT=BOTH
NODE_NUM_41,1.0,
*CONTOUR INTEGRAL,CONT=2,OUTPUT=BOTH
NODE_NUM_61,1.0,
*CONTOUR INTEGRAL,CONT=2,OUTPUT=BOTH
NODE_NUM_81,1.0,
*CONTOUR INTEGRAL,CONT=2,OUTPUT=BOTH
NODE_NUM_101,1.0,
*END STEP