*HEADING
 Test *COHESIVE BEHAVIOR for damage evolution with mode mix dependence
 TYPE=ENERGY,MODE MIX RATIO=TRACTION
*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=C3D8R, 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,6
*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.
*MATERIAL, NAME=MAT
*ELASTIC
2.E6
*DENSITY
1.,
*ELSET, ELSET=QA_TEST_EALL
SOLID,
*AMPLITUDE, NAME=DISP, DEFINITION=SMOOTH STEP
0., 0., 1., 1.
*FILTER, NAME=SMOOTH, TYPE=BUTTERWORTH
300.,4
*STEP, NLGEOM=YES
*DYNAMIC, EXPLICIT
,1.
*BOUNDARY, TYPE=DISPLACEMENT, AMPLITUDE=DISP
UPPER,1,3,1.
*CONTACT
*CONTACT INCLUSIONS
BOT,TOP
*CONTACT PROPERTY ASSIGNMENT
BOT,TOP,COHESIVE
*OUTPUT, FIELD, NUMBER=4, TIMEMARKS=YES
*NODE OUTPUT
U,
*CONTACT OUTPUT
CSDMG,CSMAXSCRT
***CFORCE, CSTRESS
*OUTPUT, HISTORY, TIME INTERVAL=0.01
*ENERGY OUTPUT
ETOTAL, ALLWK
*OUTPUT, HISTORY, FREQ=10, FILTER=SMOOTH
*CONTACT OUTPUT,SURFACE=BOT
CFT,
*END STEP