*HEADING
 Test *COHESIVE BEHAVIOR for damage evolution with mode mix dependence
 TYPE=ENERGY,MODE MIX RATIO=TRACTION
*PREPRINT, MODEL=YES, HISTORY=YES
*NODE
1, 0., 0., 0.
2, 1., 0., 0.
3, 1., 1., 0.
4, 0., 1., 0.
5, 0., 0., 1.
6, 1., 0., 1.
7, 1., 1., 1.
8, 0., 1., 1.
101, 0., 0., 1.
102, 1., 0., 1.
103, 1., 1., 1.
104, 0., 1., 1.
105, 0., 0., 2.
106, 1., 0., 2.
107, 1., 1., 2.
108, 0., 1., 2.
*NSET, NSET=UPPER
105,106,107,108
*NSET, NSET=LOWER
1,2,3,4
*ELEMENT, TYPE=C3D8, ELSET=SOLID
1, 1,2,3,4,5,6,7,8
2, 101,102,103,104,105,106,107,108
*SOLID SECTION, ELSET=SOLID, MATERIAL=MAT
1.0,
*SURFACE, NAME=BOT
1,S2
*SURFACE, NAME=TOP
2,S1
*BOUNDARY
LOWER, 1,3
*SURFACE INTERACTION, NAME=COHESIVE
*COHESIVE BEHAVIOR, Eligibility = Current Contacts 
2.E6, 1.E6, 1.E6
*DAMAGE INITIATION, CRITERION=MAXS
1.1E5, 1.E5, 1.E5
*DAMAGE EVOLUTION,TYPE=ENERGY,MODE MIX RATIO=TRACTION, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=LINEAR
122103.659214, 0., 0.
30000., 1., 0.
20000., 1., 1.
*DAMAGE STABILIZATION
1.E-4
*MATERIAL, NAME=MAT
*ELASTIC
2.E6
*DENSITY
1.,
*ELSET, ELSET=QA_TEST_EALL
SOLID,
*AMPLITUDE, NAME=DISP, DEFINITION=SMOOTH STEP
0., 0., 1., 1.
*CONTACT PAIR, INTERACTION=COHESIVE
BOT,TOP
*STEP, NLGEOM=YES, INC=1000
*STATIC
0.01, 1.0,,0.01
*BOUNDARY, TYPE=DISPLACEMENT, AMPLITUDE=DISP
UPPER,1,3,1.
*OUTPUT, FIELD, FREQ=4
*NODE OUTPUT
U,
*CONTACT OUTPUT
CSDMG,CSMAXSCRT
***CFORCE, CSTRESS
*OUTPUT, HISTORY
*ENERGY OUTPUT
ETOTAL, ALLWK
*OUTPUT, HISTORY, FREQ=4
*CONTACT OUTPUT,SECONDARY=BOT,MAIN=TOP
CFT,CSDMG,CSMAXSCRT
*END STEP