*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
*ELEMENT, TYPE=C3D8I,ELSET=SOLID
*SOLID SECTION, ELSET=SOLID, MAT=MAT1
*MATERIAL, NAME=MAT1
*ELASTIC
200E9,0.3
*ELEMENT,TYPE=COH3D8, ELSET=COHESIVE
*COHESIVE SECTION,ELSET=COHESIVE,MATERIAL=ADHESIVE,RESPONSE=TRACTION SEPARATION,
THICKNESS=GEOMETRY,CONTROLS=CTR
*SECTION CONTROLS,ELEMENT DELETION=YES,NAME=CTR,VISCOSITY=1.E-4
*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
*NSET,NSET=TOP
*BOUNDARY
BOT,1
BOT,3
nall,2
*STEP,NLGEOM=YES,INC=1000
*STATIC
 0.0001, 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
*EL PRINT
SDEG
STATUS
DMICRT
*END STEP