** This test verifies displacement type damage evolution with exponential
** softening (exponent=2) considering both energy and traction mode mix measures.
** The following have been incorporated in the test:
**	1. Uncoupled modes => *ELASTIC, TYPE=TRACTION
**	2. DAMAGE INITIATION CRITERION = QUADS
**	3. Deformation in each test case is considered either purely normal
**         (Mode 1) or purely in the first (Mode 2) or second (Mode 3) shear 
**         direction.
**	4. No temperature dependence
**	5. Isotropic shear behavior 
**	6. Thickness direction = Global 2 direction
**	7. Thickness = 1.0 => Nominal strain = Separation
*HEADING
  Damage test for COH3D8
*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=QUADS
8000., 2000., 2000.
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, MODE MIX RATIO=TRACTION, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=EXPONENTIAL
0.006, 2.0, 0.0
0.007, 2.0, 1.0
*DENSITY
100.
**
**#####################################
**TEST CASE IV: PURE TENSILE (MODE 1)
**TEST CASE V: SIMPLE SHEAR (MODE 2)
**TEST CASE VI: SIMPLE SHEAR (MODE 3)  
**#####################################
*SYSTEM
0, 2, 0
*NODE, NSET=SPEC2_MODE1
41, 0., 0., 0.
42, 1., 0., 0.
43, 1., 0., 1.
44, 0., 0., 1.
45, 0., 1., 0.
46, 1., 1., 0.
47, 1., 1., 1.
48, 0., 1., 1.
*SYSTEM
2, 2, 0
*NODE, NSET=SPEC2_MODE2
51, 0., 0., 0.
52, 1., 0., 0.
53, 1., 0., 1.
54, 0., 0., 1.
55, 0., 1., 0.
56, 1., 1., 0.
57, 1., 1., 1.
58, 0., 1., 1.
*SYSTEM
4, 2, 0
*NODE, NSET=SPEC2_MODE3
61, 0., 0., 0.
62, 1., 0., 0.
63, 1., 0., 1.
64, 0., 0., 1.
65, 0., 1., 0.
66, 1., 1., 0.
67, 1., 1., 1.
68, 0., 1., 1.
*NSET, NSET=SPEC2_MODE1_TOP
45, 46, 47, 48
*NSET, NSET=SPEC2_MODE1_BOT
41, 42, 43, 44
*NSET, NSET=SPEC2_MODE2_TOP
55, 56, 57, 58
*NSET, NSET=SPEC2_MODE2_BOT
51, 52, 53, 54
*NSET, NSET=SPEC2_MODE3_TOP
65, 66, 67, 68
*NSET, NSET=SPEC2_MODE3_BOT
61, 62, 63, 64
*ELEMENT, TYPE=COH3D8, ELSET=SPEC2_MODE1
41, 41, 44, 43, 42, 45, 48, 47, 46
*ELEMENT, TYPE=COH3D8, ELSET=SPEC2_MODE2
51, 51, 54, 53, 52, 55, 58, 57, 56
*ELEMENT, TYPE=COH3D8, ELSET=SPEC2_MODE3
61, 61, 64, 63, 62, 65, 68, 67, 66
*ELSET, ELSET=EALL_ENER
SPEC2_MODE1, SPEC2_MODE2, SPEC2_MODE3
**
*COHESIVE SECTION, RESPONSE=TRACTION SEPARATION, ELSET=EALL_ENER, MATERIAL=COH_ENER
*MATERIAL, NAME=COH_ENER
*ELASTIC,TYPE=TRACTION
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=QUADS
8000., 2000., 2000.
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=EXPONENTIAL
0.006, 2.0, 0.0
0.007, 2.0, 1.0
*DENSITY
100.
**
*NSET, NSET=NALL_BOT
SPEC1_MODE1_BOT, SPEC1_MODE2_BOT, SPEC1_MODE3_BOT, SPEC2_MODE1_BOT, SPEC2_MODE2_BOT, SPEC2_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
SPEC2_MODE1_TOP, 1, 1
SPEC2_MODE1_TOP, 3, 3
SPEC2_MODE2_TOP, 2, 3
SPEC2_MODE3_TOP, 1, 2
*ELSET, ELSET=QA_TEST_EALL
EALL_DISP, EALL_ENER
*NSET, NSET=QA_TEST_NALL_TOP
SPEC1_MODE1_TOP, SPEC1_MODE2_TOP, SPEC1_MODE3_TOP, SPEC2_MODE1_TOP, SPEC2_MODE2_TOP, SPEC2_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
SPEC2_MODE1_TOP, 2, 2, 0.01
SPEC2_MODE2_TOP, 1, 1, 0.01
SPEC2_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
*ENERGY OUTPUT, ELSET=SPEC1_MODE1
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC1_MODE2
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC1_MODE3
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC2_MODE1
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC2_MODE2
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC2_MODE3
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT
ETOTAL, ALLWK
*END STEP