*HEADING
CRACK PROPAGATION ALONG AN INTERFACE VIA CRITICAL STRESS CRITERION
** Centre crack specimen:complete model.
*RESTART,WRITE,FREQ=1
*NODE
1,
25,24.,
201,,10.
225,24.,10.
999, 20.0,0.0
*NGEN, NSET=BOT1
1, 25
*NGEN,NSET=TOP1
201,225
*NFILL, NSET=N1
BOT1,TOP1,4,50
*ELEMENT,TYPE=CPE4
1,1,3,103,101
*ELGEN,ELSET=SOL1
1,12,2,1,2,100,100
*NCOPY, OLD SET=N1,NEW SET=N2,REFLECT=LINE,CHANGE NUMBER=1000
-1.0, 0.0, 0.0, 30.0,0.0, 0.0
*ELCOPY,ELEMENT SHIFT=1000,OLD SET=SOL1,NEW SET=SOL2,REFLECT,SHIFT NODES=1000
*ELEMENT, TYPE=CPE4, ELSET=SOL2
1001, 1, 1101,1103,1003
1012,1023,1123,1125,25
*ELSET,ELSET=SOLID
SOL1,SOL2
*NSET, NSET=TOP2, GEN
1201,1225
*SOLID SECTION,ELSET=SOLID,MAT=ELASTIC
1., 
*MATERIAL,NAME=ELASTIC
*ELASTIC
1.E6,.3
**
*** Contact definitions.
**
*ELSET, ELSET=TOPSL, GEN
2,11
*ELSET, ELSET=BOTMS, GEN
1002, 1011
*ELSET, ELSET=TOPDL, GEN
101,112
*ELSET,ELSET=BOTDL, GEN
1101, 1112
*surface, NAME=TOPSL
TOPSL, S1
*surface, NAME=BOTMS
BOTMS, S4
*SURFACE INTERACTION, NAME=BOTH
1.0, 
*CONTACT PAIR,INTERACTION=BOTH,SMALL SLIDING
TOPSL,BOTMS
**
** Nodeset definitions for bonding
*NSET,NSET=BOND,GEN
1,11
15,25
*INITIAL CONDITIONS,TYPE=CONTACT
TOPSL,BOTMS,BOND
*BOUNDARY
1, 2
25, 2
213, 1
1213, 1
**
**** step definitioN
**
*STEP,INC=20
*STATIC
8.0E-04, 0.001313
*DEBOND,SLAVE=TOPSL,MASTER=BOTMS,TIME INC=1.E-04
0.,1.
.0003,0.
*FRACTURE CRITERION,TYPE=CRITICAL STRESS,TOLERANCE=0.01,DISTANCE=0.0
30000, 1.E12, 1.E12
*BOUNDARY
TOP1, 2, 2, 0.2363
TOP2, 2, 2, -0.2363
*ENERGY PRINT
*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
S, E
*END STEP
**
***  step 2: 
**
*STEP,INC=20
*STATIC
1.0E-04, 0.001187
*DEBOND,SLAVE=TOPSL,MASTER=BOTMS,TIME INC=1.E-04
0.,1.
.0003,0.
*FRACTURE CRITERION,TYPE=CRITICAL STRESS,TOLERANCE=0.01,DISTANCE=0.0
59063, 1.E12, 1.E12
*BOUNDARY
TOP1, 2, 2, 0.45
TOP2, 2, 2, -0.45
*END STEP