*HEADING
 THREE-DIMENSIONAL  TEST [COH3D8] 
 -- 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,nset=nall
1,0.,0.,-1.0
2,1.,0.,-1.0
3,1.,1.,-1.0
4,0.,1.,-1.0
10,0.,0.,0.
20,1.,0.,0.
30,1.,1.,0.
40,0.,1.,0.
50,0.,0.,0.01325E-3
60,1.,0.,0.01325E-3
70,1.,1.,0.01325E-3
80,0.,1.,0.01325E-3
5 ,0.,0.,1.0
6 ,1.,0.,1.0
7 ,1.,1.,1.0
8 ,0.,1.,1.0
*ELEMENT, TYPE=C3D8I,ELSET=SOLID
10,1,2,3,4,10,20,30,40
20,50,60,70,80,5,6,7,8
*SOLID SECTION, ELSET=SOLID, MAT=MAT1
*MATERIAL, NAME=MAT1
*ELASTIC
200E9,0.3
*ELEMENT,TYPE=COH3D8, ELSET=COHESIVE
1,10,20,30,40,50,60,70,80
*COHESIVE SECTION,ELSET=COHESIVE,MATERIAL=ADHESIVE,RESPONSE=TRACTION SEPARATION,
THICKNESS=GEOMETRY,CONTROLS=CTR
*SECTION CONTROLS,ELEMENT DELETION=YES,NAME=CTR,VISCOSITY=5.E-3
*MATERIAL,NAME=ADHESIVE
*ELASTIC,TYPE=COUPLED TRACTION
<EMOD>,<EMOD>,<EMOD>,<EK12>,<EK13>,<EK23>
*DAMAGE INITIATION,CRITERION=MAXS
<ULTI>, <ULTII>,<ULTII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=POWER LAW,POWER=2.0,
DEGRADATION=MULTIPLICATIVE
<GIC>, <GIIC>,<GIIC>
*DAMAGE INITIATION,CRITERION=QUADS
<ULTI>, <ULTII>,<ULTII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=BK,POWER=<ETA>
<GIC>, <GIIC>,<GIIC>
*DAMAGE INITIATION,CRITERION=MAXE
<NSTRNI>,<NSTRNII>,<NSTRNII>
**<DISI>, <DISII>, <DISII>
*DAMAGE EVOLUTION,TYPE=DISPLACEMENT
1.995E-4
*DAMAGE INITIATION,CRITERION=QUADE
<NSTRNI>,<NSTRNII>,<NSTRNII>
**<DISI>, <DISII>, <DISII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=BK,POWER=1.0,
DEGRADATION=MULTIPLICATIVE
<GIC>, <GIIC>,<GIIC>
**
*NSET,NSET=BOT
1,2,3,4
*NSET,NSET=TOP
5,6,7,8
*BOUNDARY
BOT,1
BOT,3
nall,2
*STEP,NLGEOM=YES,INC=1000
*STATIC
 0.01, 1.0, 1.E-8, 0.05
*BOUNDARY
TOP,1,1,0.00002
**TOP,2,2,0.00003
TOP,3,3,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
SDEG
*EL PRINT
SDEG
STATUS
DMICRT
*END STEP