*HEADING
: DEBONDING OF AN INTERFACE VIA CRITICAL STRESS CRITERION
** complete model.
**preprint,history=yes,model=yes
*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=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=CPE4, ELSET=SOL2
2101,3101,4101,4121,121
*ELSET, ELSET=SOLID
SOL1,SOL2
*SOLID SECTION,ELSET=SOLID,MAT=ELASTIC
1., 
*MATERIAL,NAME=ELASTIC
*ELASTIC
1.E6,.3
**
*** Contact definitions
**
*ELSET, ELSET=TOP
1,21,41,61,81
*surface, NAME=BSURF
TOP, S1
*ELSET, ELSET=BOT
2001,2021,2041,2061,2081
*surface, NAME=ASURF
BOT, S4
*SURFACE INTERACTION, NAME=BOTH
1.0, 
*CONTACT PAIR,INTERACTION=BOTH,SMALL SLIDING
BSURF, ASURF
*NSET,NSET=N41
41,
3041,
*NSET,NSET=N61
61,
3061,
*NSET,NSET=N81
81,
3081,
*NSET,NSET=N101
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
BSURF,ASURF,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=BSURF,MASTER=ASURF,TIME INC=.0002, FREQUENCY=1,OUTPUT=BOTH
0.,1.
.0002,0.,
*FRACTURE CRITERION,TYPE=CRITICAL STRESS,TOLERANCE=0.05,DISTANCE=0.0
4.5E+4,1.E12,1.E12
*ENERGY PRINT
*EL PRINT
S,
E,
*CONTACT PRINT
CSTRESS,CDISP
DBT,DBSF,DBS
*EL PRINT,ELSET=SOLID
S,E
*CONTOUR INTEGRAL,CONT=2,OUTPUT=BOTH
N41,1.0,
*CONTOUR INTEGRAL,CONT=2,OUTPUT=BOTH
N61,1.0,
*CONTOUR INTEGRAL,CONT=2,OUTPUT=BOTH
N81,1.0,
*CONTOUR INTEGRAL,CONT=2,OUTPUT=BOTH
N101,1.0,
*OUTPUT,FIELD
*CONTACT OUTPUT, VARIABLE=ALL
CSTRESS,CDISP
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=SOLID
S, E
**NODE OUTPUT
**U,
*END STEP