*HEADING
 CRACK GROWTH IN A THREE POINT BEND SPECIMEN:  
 COD CRITERION.   
**
** Input file defining the model data.
**
*INCLUDE, INPUT=crackgrowth_model.inp
**      
** step  1: load the specimen
**      
*STEP,NLGEOM    
*STATIC 
0.001, 1.0      
**
*** load application    
**      
*BOUNDARY       
9997, 1, 1, 1.0 
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=0
*CONTACT PRINT,SLAVE=DBDSL,MASTER=DBDMS 
DBT, DBSF, DBS  
*CONTACT FILE,FREQUENCY=10,SLAVE=DBDSL,
MASTER=DBDMS
DBT, DBSF, DBS  
*PRINT, CONTACT=YES     
*END STEP       
**
*** step 2: Debond the first node.
**
*STEP,NLGEOM,INC=1
*STATIC
0.005, 0.005, 0.00001
*BOUNDARY
9997, 1, 1, 1.005
*DEBOND,SLAVE=DBDSL,MASTER=DBDMS,
TIME INCREMENT=.005,FREQUENCY=1,OUTPUT=BOTH
0.0   , 1.0
0.005, 0.0
*FRACTURE CRITERION,TYPE=COD,TOLERANCE=0.01,
DISTANCE=2.36,SYMMETRY
2.446, 0.028
*CONTROLS, PARAMETERS=FIELD, FIELD=DISPLACEMENT
 , 0.1
*CONTACT FILE,SLAVE=DBDSL,MASTER=DBDMS
DBT,DBSF,DBS
*CONTACT PRINT,SLAVE=DBDSL,MASTER=DBDMS
DBT,DBSF,DBS
*PRINT, CONTACT=YES
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=0
*END STEP
**
** step 3: Load to the point where the second
** node is about to debond.
**
*STEP,NLGEOM
*STATIC
0.001, 0.0644, 0.00001
*BOUNDARY
9997, 1, 1, 1.0694
*CONTACT FILE,FREQUENCY=10,SLAVE=DBDSL,
MASTER=DBDMS
DBT,DBSF,DBS
*CONTACT PRINT,SLAVE=DBDSL,MASTER=DBDMS
DBT,DBSF,DBS
*PRINT, CONTACT=YES
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=0
*END STEP
***
***** step 4: Debond the second node.
***
*STEP,NLGEOM
*STATIC
1.0E-06, 1.E-06
*BOUNDARY
9997, 1, 1, 1.069401
*DEBOND,SLAVE=DBDSL,MASTER=DBDMS,
TIME INCREMENT=.001,FREQUENCY=2,OUTPUT=BOTH
0.0   , 1.0
0.005, 0.0
*FRACTURE CRITERION,TYPE=COD,TOLERANCE=0.01,
DISTANCE=2.389,SYMMETRY
2.475, 0.06
*CONTACT FILE,SLAVE=DBDSL,MASTER=DBDMS
DBT,DBSF,DBS
*CONTACT PRINT,SLAVE=DBDSL,MASTER=DBDMS
DBT,DBSF,DBS
*PRINT, CONTACT=YES
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=0
*END STEP
**
** step 5: Load to the point where the third
** node is about to debond
**
*STEP,NLGEOM
*STATIC
0.001, 0.085, 0.00001
*BOUNDARY
9997, 1, 1, 1.1544
*CONTACT FILE,FREQUENCY=10,SLAVE=DBDSL,
MASTER=DBDMS
DBT,DBSF,DBS
*CONTACT PRINT,SLAVE=DBDSL,MASTER=DBDMS
DBT,DBSF,DBS
*PRINT, CONTACT=YES
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=0
*END STEP
***
***** step 6: Debond the third node
***
*STEP,NLGEOM
*STATIC
1.0E-06, 1.E-06
*BOUNDARY
9997, 1, 1, 1.154401
*DEBOND,SLAVE=DBDSL,MASTER=DBDMS,
TIME INCREMENT=.001,FREQUENCY=2,OUTPUT=BOTH
0.0   , 1.0
0.005, 0.0
*FRACTURE CRITERION,TYPE=COD,TOLERANCE=0.01,
DISTANCE=2.426,SYMMETRY
2.518, 0.0934
*CONTACT FILE,SLAVE=DBDSL,MASTER=DBDMS
DBT,DBSF,DBS
*CONTACT PRINT,SLAVE=DBDSL,MASTER=DBDMS
DBT,DBSF,DBS
*PRINT, CONTACT=YES
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=0
*END STEP
**
*** step 7: Debond the rest of the nodes.
**
**
** Use the following lines if the analysis is
** to run through to completion
**STEP,NLGEOM,INC=500
**STATIC
**0.001, 1.752, 0.00001
**BOUNDARY
**9997, 1, 1, 2.951
**
*STEP,NLGEOM
*STATIC
0.001, 0.40, 0.00001
*BOUNDARY
9997, 1, 1, 1.40
*DEBOND,SLAVE=DBDSL,MASTER=DBDMS,
TIME INCREMENT=.001,FREQUENCY=2,OUTPUT=BOTH
0.0   , 1.0
0.005, 0.0
*FRACTURE CRITERION,TYPE=COD,TOLERANCE=0.01,
DISTANCE=0.04,SYMMETRY
0.023, 0.262
0.021, 1.037
0.019, 1.66
0.012, 2.18
*CONTACT FILE,FREQUENCY=10,SLAVE=DBDSL,
MASTER=DBDMS
DBT,DBSF,DBS
*CONTACT PRINT,SLAVE=DBDSL,MASTER=DBDMS
DBT,DBSF,DBS
*PRINT, CONTACT=YES
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=0
*END STEP