*HEADING
PDEBCODCO1: CRACK PROPAGATION AND COALESCENCE: CRACK LENGTH VS TIME CRITERION
** symmetry conditions taken into account.
*RESTART,WRITE,FREQ=1
*NODE
1,
25,24.,
201,,10.
225,24.,10.
999, 20.0,0.0
*NGEN, NSET=BOT
1, 25
*NGEN,NSET=TOP
201,225
*NFILL
BOT,TOP,4,50
*NODE
13,12.0,0.00001
*ELEMENT,TYPE=CPE4
1,1,3,103,101
*ELGEN,ELSET=SOLID
1,12,2,1,2,100,100
*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, 25.1, 0.0
**
*** Contact definitions.
**
*ELSET, ELSET=BOT, GEN
1, 11
*ELSET, ELSET=TOP, GEN
101, 112
*surface, NAME=BSURF
BOT, S1
*SURFACE INTERACTION, NAME=BOTH
1.0, 
*CONTACT PAIR,INTERACTION=BOTH,SMALL SLIDING
BSURF, ASURF
**
** Nodeset definitions for bonding
**
*NSET,NSET=BOND,GEN
5,11
15,23
*INITIAL CONDITIONS,TYPE=CONTACT, NORMAL
BSURF,ASURF,BOND
*NSET, NSET=NREF
1, 
*BOUNDARY
25, 1, 2
999, 1, 2
999, 6
**
*** step definitions
**
*STEP,INC=50
*STATIC
0.001,0.009
*DEBOND,SLAVE=BSURF,MASTER=ASURF,TIME INC=3.5E-04
0.,1.
.0009,0.
*FRACTURE CRITERION,TYPE=CRACK LENGTH,TOLERANCE=0.01,NSET=NREF, SYMMETRY
 2.1586E-03, 4.0
 2.8853E-03, 6.0
 3.5658E-03, 8.0
 4.2725E-03, 10.0
 5.4874E-03, 14.0
 6.3624E-03, 16.0
 7.2374E-03, 18.0
 8.1124E-03, 20.0
 8.9874E-03, 22.0
*BOUNDARY
201, 2, 2, 0.18
203, 2, 2, 0.18
*ENERGY PRINT
*CONTACT PRINT
CSTRESS,CDISP
DBT,DBSF,DBS
*EL PRINT,ELSET=SOLID
S, E
*OUTPUT,FIELD
*CONTACT OUTPUT, VARIABLE=ALL
CSTRESS,CDISP
*CONTACT OUTPUT
DBT,DBSF,DBS
*OUTPUT,FIELD
*ELEMENT OUTPUT
S, E
*END STEP