*HEADING
 Verification for coupled traction and multiple damage criteria
*NODE, NSET=NALL
 1, 0., 0.
 2, 1., 0.
 3, 1., 1.
 4, 0., 1.
*ELEMENT, TYPE=COH2D4, ELSET=EALL
1, 1, 2, 3, 4
*COHESIVE SECTION, MATERIAL=COH, ELSET=EALL, RESPONSE=TRACTION SEPARATION, THICKNESS=GEOMETRY, CONTROLS=CONT
*MATERIAL, NAME=COH
*DENSITY
 10.
*ELASTIC, TYPE=COUPLED TRACTION
4.E5, 2.E5, 2.E5, 1.E5, 2.E5, 3.E5
*DAMAGE INITIATION, CRITERION=MAXS
8000., 2000., 2000.
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR, DEGRADATION=MULTIPLICATIVE
0.06
*DAMAGE INITIATION, CRITERION=MAXE
0.02, 0.01, 0.01
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR, DEGRADATION=MAXIMUM
0.06
*DAMAGE INITIATION, CRITERION=QUADS
8000., 2000., 2000.
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR, DEGRADATION=MAXIMUM
0.06
*DAMAGE INITIATION, CRITERION=QUADE
0.02, 0.01, 0.01
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR, DEGRADATION=MULTIPLICATIVE
0.06
*SECTION CONTROLS, NAME=CONT, ELEMENT DELETION=YES
*NSET, NSET=TOP
3, 4, 
*NSET, NSET=BOT
1, 2
*BOUNDARY
BOT, 1, 2
*STEP
*DYNAMIC, EXPLICIT
, 1.
*BOUNDARY, TYPE=VELOCITY
TOP, 1, 1, 0.10
TOP, 2, 2, 0.10
*BULK VISCOSITY
0.0, 0.0
*OUTPUT, FIELD, NUMBER INTERVAL=20, TIME MARKS=YES
*ELEMENT OUTPUT
S, LE, SDEG, DMICRT, STATUS
*NODE OUTPUT
U,
*OUTPUT, HISTORY, TIME INTERVAL=0.01
*ELEMENT OUTPUT, ELSET=EALL
S, LE, SDEG, DMICRT, STATUS
*NODE OUTPUT, NSET=NALL
U,
*ENERGY OUTPUT, VARIABLE=PRESELECT, ELSET=EALL
*END STEP