** Test to verify damage evolution independent of mode mix. 
*HEADING
  Damage evolution independent of mode mix
*PREPRINT, MODEL=YES
**#####################################
**TEST CASE I: PURE TENSILE (MODE 1)
**TEST CASE II: SIMPLE SHEAR (MODE 2)
**TEST CASE III: SIMPLE SHEAR (MODE 3)  
**#####################################
*NODE, NSET=SPEC1_MODE1
11, 0., 0., 0.
12, 1., 0., 0.
13, 1., 0., 1.
14, 0., 0., 1.
15, 0., 1., 0.
16, 1., 1., 0.
17, 1., 1., 1.
18, 0., 1., 1.
*SYSTEM
2, 0, 0
*NODE, NSET=SPEC1_MODE2
21, 0., 0., 0.
22, 1., 0., 0.
23, 1., 0., 1.
24, 0., 0., 1.
25, 0., 1., 0.
26, 1., 1., 0.
27, 1., 1., 1.
28, 0., 1., 1.
*SYSTEM
4, 0, 0
*NODE, NSET=SPEC1_MODE3
31, 0., 0., 0.
32, 1., 0., 0.
33, 1., 0., 1.
34, 0., 0., 1.
35, 0., 1., 0.
36, 1., 1., 0.
37, 1., 1., 1.
38, 0., 1., 1.
*NSET, NSET=SPEC1_MODE1_TOP
15, 16, 17, 18
*NSET, NSET=SPEC1_MODE1_BOT
11, 12, 13, 14
*NSET, NSET=SPEC1_MODE2_TOP
25, 26, 27, 28
*NSET, NSET=SPEC1_MODE2_BOT
21, 22, 23, 24
*NSET, NSET=SPEC1_MODE3_TOP
35, 36, 37, 38
*NSET, NSET=SPEC1_MODE3_BOT
31, 32, 33, 34
*ELEMENT, TYPE=COH3D8, ELSET=SPEC1_MODE1
11, 11, 14, 13, 12, 15, 18, 17, 16
*ELEMENT, TYPE=COH3D8, ELSET=SPEC1_MODE2
21, 21, 24, 23, 22, 25, 28, 27, 26
*ELEMENT, TYPE=COH3D8, ELSET=SPEC1_MODE3
31, 31, 34, 33, 32, 35, 38, 37, 36
*ELSET, ELSET=EALL_DISP
SPEC1_MODE1, SPEC1_MODE2, SPEC1_MODE3
**
*COHESIVE SECTION, RESPONSE=TRACTION SEPARATION, ELSET=EALL_DISP, MATERIAL=COH_DISP
*MATERIAL, NAME=COH_DISP
*ELASTIC,TYPE=TRACTION
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=MAXS
8000., 2000., 2000.
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR
0.006
*DENSITY
100.
**
*NSET, NSET=NALL_BOT
SPEC1_MODE1_BOT, SPEC1_MODE2_BOT, SPEC1_MODE3_BOT
*BOUNDARY
NALL_BOT, 1, 3
SPEC1_MODE1_TOP, 1, 1
SPEC1_MODE1_TOP, 3, 3
SPEC1_MODE2_TOP, 2, 3
SPEC1_MODE3_TOP, 1, 2
*ELSET, ELSET=QA_TEST_EALL
EALL_DISP
*NSET, NSET=QA_TEST_NALL_TOP
SPEC1_MODE1_TOP, SPEC1_MODE2_TOP, SPEC1_MODE3_TOP
**#####################################
**STATIC STEP
**#####################################
*STEP, INC=1000
*STATIC
0.001, 1.0, 1.0E-5, 0.01
*BOUNDARY
SPEC1_MODE1_TOP, 2, 2, 0.01
SPEC1_MODE2_TOP, 1, 1, 0.01
SPEC1_MODE3_TOP, 3, 3, 0.01
*OUTPUT, FIELD
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
S, E, SDEG, DMICRT, STATUS
*NODE OUTPUT, NSET=QA_TEST_NALL_TOP
U, RF
*OUTPUT, HISTORY
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
S, E, SDEG, DMICRT, STATUS
*END STEP