*HEADING
 THREE-DIMENSIONAL  TEST [COH2D4] 
 -- Using coupled traction elasticity and multi damage models
*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES
*PARAMETER
 ULTI   =  3.30E6
 ULTII  =  7.0E6
** DISI = 5.0E-8
** DISII= 2.5E-8
 NSTRNI  = 3.77358490566E-3
 NSTRNII = 1.88679245283E-3 
 EMOD = 8.5E8
 EK12 = 1.E7
 EK13 = 2.E7
 EK23 = 3.E7
 GIC   = 0.33E3
 GIIC  = 0.80E3
 ETA = 2.284
 THICK = 0.01325E-3
*NODE
 1,0.,-1.0
 2,1.,-1.0
10,0., 0.
20,1., 0.
50,0. ,0.01325E-3
60,1. ,0.01325E-3
5 ,0. ,1.0
6 ,1. ,1.0
*ELEMENT, TYPE=CPE4I,ELSET=SOLID
10,1,2,20,10
20,50,60,6,5
*SOLID SECTION, ELSET=SOLID, MAT=MAT1
*MATERIAL, NAME=MAT1
*ELASTIC
200E9,0.3
*ELEMENT,TYPE=COH2D4, ELSET=COHESIVE
1,10,20,60,50
*COHESIVE SECTION,ELSET=COHESIVE,MATERIAL=ADHESIVE,
RESPONSE=TRACTION SEPARATION,THICKNESS=GEOMETRY,CONTROLS=CTR
*SECTION CONTROLS,ELEMENT DELETION=YES,NAME=CTR,VISCOSITY=5.E-4
*MATERIAL,NAME=ADHESIVE
*ELASTIC,TYPE=COUPLED TRACTION
<EMOD>,<EMOD>,<EMOD>,<EK12>,<EK13>,<EK23>
*DAMAGE INITIATION,CRITERION=MAXS
<ULTI>, <ULTII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=POWER LAW,POWER=2.0,
DEGRADATION=MULTIPLICATIVE
<GIC>, <GIIC>
*DAMAGE INITIATION,CRITERION=QUADS
<ULTI>, <ULTII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=BK,POWER=<ETA>
<GIC>, <GIIC>
*DAMAGE INITIATION,CRITERION=MAXE
<NSTRNI>,<NSTRNII>
** <DISI>, <DISII>
*DAMAGE EVOLUTION,TYPE=DISPLACEMENT
1.995E-4
*DAMAGE INITIATION,CRITERION=QUADE
<NSTRNI>,<NSTRNII>
**<DISI>, <DISII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=BK,POWER=1.0,
DEGRADATION=MULTIPLICATIVE
<GIC>, <GIIC>
**
*NSET,NSET=BOT
1,2
*NSET,NSET=TOP
5,6
*BOUNDARY
BOT,1,2
*STEP,NLGEOM=YES,INC=1000
*STATIC
 0.01, 1.0,, 0.005
*BOUNDARY
TOP,1,1,0.00002
TOP,2,2,0.00004
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
STATUS
DMICRT
SDEG
*NODE OUTPUT
U,
RF,
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT, NSET=TOP
U,
RF
*ELEMENT OUTPUT,ELSET=COHESIVE
S,
E,
STATUS
DMICRT
*EL PRINT
SDEG
STATUS
DMICRT
*END STEP