** This test verifies energy type damage evolution with linear
** softening considering energy mode mix measure.
** 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 shear (Mode 2) 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 COH2D4
*PREPRINT, MODEL=YES
**#####################################
**TEST CASE I: PURE TENSILE (MODE 1)
**TEST CASE II: SIMPLE SHEAR (MODE 2)
**#####################################
*NODE, NSET=SPEC1_MODE1
11, 0., 0.
12, 1., 0.
13, 1., 1.
14, 0., 1.
*SYSTEM
2, 0, 0
*NODE, NSET=SPEC1_MODE2
21, 0., 0.
22, 1., 0.
23, 1., 1.
24, 0., 1.
*NSET, NSET=SPEC1_MODE1_TOP
13, 14
*NSET, NSET=SPEC1_MODE1_BOT
11, 12
*NSET, NSET=SPEC1_MODE2_TOP
23, 24
*NSET, NSET=SPEC1_MODE2_BOT
21, 22
*ELEMENT, TYPE=COH2D4, ELSET=SPEC1_MODE1
11, 11, 12, 13, 14
*ELEMENT, TYPE=COH2D4, ELSET=SPEC1_MODE2
21, 21, 22, 23, 24
*ELSET, ELSET=EALL_TAB
SPEC1_MODE1, SPEC1_MODE2
**
*COHESIVE SECTION, RESPONSE=TRACTION SEPARATION, ELSET=EALL_TAB, MATERIAL=COH_TAB
*MATERIAL, NAME=COH_TAB
*ELASTIC,TYPE=TRACTION
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=QUADS
8000., 2000.
*DAMAGE EVOLUTION, TYPE=ENERGY, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=LINEAR
32.0, 0.0
8.0, 1.0
*DENSITY
100.
**
**#####################################
**TEST CASE III: PURE TENSILE (MODE 1)
**TEST CASE IV: SIMPLE SHEAR (MODE 2)
**#####################################
*SYSTEM
0, 2, 0
*NODE, NSET=SPEC2_MODE1
31, 0., 0.
32, 1., 0.
33, 1., 1.
34, 0., 1.
*SYSTEM
2, 2, 0
*NODE, NSET=SPEC2_MODE2
41, 0., 0.
42, 1., 0.
43, 1., 1.
44, 0., 1.
*NSET, NSET=SPEC2_MODE1_TOP
33, 34
*NSET, NSET=SPEC2_MODE1_BOT
31, 32
*NSET, NSET=SPEC2_MODE2_TOP
43, 44
*NSET, NSET=SPEC2_MODE2_BOT
41, 42
*ELEMENT, TYPE=COH2D4, ELSET=SPEC2_MODE1
31, 31, 32, 33, 34
*ELEMENT, TYPE=COH2D4, ELSET=SPEC2_MODE2
41, 41, 42, 43, 44
*ELSET, ELSET=EALL_POWER
SPEC2_MODE1, SPEC2_MODE2
**
*COHESIVE SECTION, RESPONSE=TRACTION SEPARATION, ELSET=EALL_POWER, MATERIAL=COH_POWER
*MATERIAL, NAME=COH_POWER
*ELASTIC,TYPE=TRACTION
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=QUADS
8000., 2000.
*DAMAGE EVOLUTION, TYPE=ENERGY, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=POWER, SOFTENING=LINEAR, POWER=1.5
32.0, 8.0, 8.0
*DENSITY
100.
**#####################################
**TEST CASE V: PURE TENSILE (MODE 1)
**TEST CASE VI: SIMPLE SHEAR (MODE 2)
**#####################################
*SYSTEM
0, 4, 0
*NODE, NSET=SPEC3_MODE1
51, 0., 0.
52, 1., 0.
53, 1., 1.
54, 0., 1.
*SYSTEM
2, 4, 0
*NODE, NSET=SPEC3_MODE2
61, 0., 0.
62, 1., 0.
63, 1., 1.
64, 0., 1.
*NSET, NSET=SPEC3_MODE1_TOP
53, 54
*NSET, NSET=SPEC3_MODE1_BOT
51, 52
*NSET, NSET=SPEC3_MODE2_TOP
63, 64
*NSET, NSET=SPEC3_MODE2_BOT
61, 62
*ELEMENT, TYPE=COH2D4, ELSET=SPEC3_MODE1
51, 51, 52, 53, 54
*ELEMENT, TYPE=COH2D4, ELSET=SPEC3_MODE2
61, 61, 62, 63, 64
*ELSET, ELSET=EALL_BK
SPEC3_MODE1, SPEC3_MODE2
**
*COHESIVE SECTION, RESPONSE=TRACTION SEPARATION, ELSET=EALL_BK, MATERIAL=COH_BK
*MATERIAL, NAME=COH_BK
*ELASTIC,TYPE=TRACTION
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=QUADS
8000., 2000.
*DAMAGE EVOLUTION, TYPE=ENERGY, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=BK, SOFTENING=LINEAR, POWER=1.5
32.0, 8.0, 8.0
*DENSITY
100.
**
*NSET, NSET=NALL_BOT
SPEC1_MODE1_BOT, SPEC1_MODE2_BOT, SPEC2_MODE1_BOT, SPEC2_MODE2_BOT, SPEC3_MODE1_BOT, SPEC3_MODE2_BOT
*BOUNDARY
NALL_BOT, 1, 2
SPEC1_MODE1_TOP, 1, 1
SPEC1_MODE2_TOP, 2, 2
SPEC2_MODE1_TOP, 1, 1
SPEC2_MODE2_TOP, 2, 2
SPEC3_MODE1_TOP, 1, 1
SPEC3_MODE2_TOP, 2, 2
*ELSET, ELSET=QA_TEST_EALL
EALL_TAB, EALL_POWER, EALL_BK
*NSET, NSET=QA_TEST_NALL_TOP
SPEC1_MODE1_TOP, SPEC1_MODE2_TOP, SPEC2_MODE1_TOP, SPEC2_MODE2_TOP, SPEC3_MODE1_TOP, SPEC3_MODE2_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
SPEC2_MODE1_TOP, 2, 2, 0.01
SPEC2_MODE2_TOP, 1, 1, 0.01
SPEC3_MODE1_TOP, 2, 2, 0.01
SPEC3_MODE2_TOP, 1, 1, 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=SPEC2_MODE1
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC2_MODE2
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC3_MODE1
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC3_MODE2
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT
ETOTAL, ALLWK
*END STEP