*HEADING
: DEBONDING OF AN INTERFACE VIA CRITICAL STRESS CRITERION
** complete model.
*RESTART,WRITE,FREQ=1
*NODE
1,2
121,14.,
2001,2.0,10.
2121,14.,10.
2,2.0,
122,14.
999, 14.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=CAX4
1,1,21,1021,1001,11,521,1011,501
*ELGEN,ELSET=SOL1
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=SOL1, NEW SET=SOL2,SHIFT NODES=3000,REFLECT
** element redefined to get correct node definitions
*ELEMENT, TYPE=CAX4, ELSET=SOL2
2001,1,4001,4021,3021
*ELSET, ELSET=SOLID
SOL1,SOL2
*SOLID SECTION,ELSET=SOLID,MAT=ELASTIC
1., 
*MATERIAL,NAME=ELASTIC
*ELASTIC
5.E7,.3
**
*** Contact definitions
**
*ELSET, ELSET=TOP
21,41,61,81,101
*surface, NAME=BSURF
TOP, S1
*ELSET, ELSET=BOT
2021,2041,2061,2081,2101
*surface, NAME=ASURF
BOT, S4
*SURFACE INTERACTION, NAME=BOTH
1.0, 
*CONTACT PAIR,INTERACTION=BOTH,SMALL SLIDING
BSURF, ASURF
*NSET,NSET=NREF
1,
***
**** Nodeset defining the bonded portion of the surface
***
*NSET, NSET=BOND
21,41,61,81
*INITIAL CONDITIONS, TYPE=CONTACT
BSURF,ASURF,BOND
*ELSET,ELSET=TOPDL, GEN
1001, 1101,20
*ELSET,ELSET=BOTDL, GEN
3001,3101, 20
*NSET, NSET=AXIS
1,1001,2001,4001,5001
**
***
*BOUNDARY
1,1, 2
**
*STEP,INC=20
*STATIC
0.0008, 0.002
*DLOAD
TOPDL, P3, -5.E05
BOTDL, P2, -5.E05
*DEBOND,SLAVE=BSURF,MASTER=ASURF,TIME INC=.0001, FREQUENCY=1,OUTPUT=BOTH
0.,1.
.0005,0.,
*FRACTURE CRITERION,TYPE=CRITICAL STRESS,TOLERANCE=0.01,DISTANCE=0.0
5.E+5, 1.E12, 1.E12
*ENERGY PRINT
*EL PRINT
S,
E,
*CONTACT PRINT
CSTRESS,CDISP
DBT,DBSF,DBS
*EL PRINT,ELSET=SOLID
S,E
*OUTPUT,FIELD
*CONTACT OUTPUT, VARIABLE=ALL
CSTRESS,CDISP
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=SOLID
S, E
*END STEP