** This test verifies displacement type damage evolution with linear
** softening 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 = MAXE
**	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 1 direction
**	7. Thickness = 1.0 => Nominal strain = Separation
*HEADING
  Damage test for COHAX4
*PREPRINT, MODEL=YES
**
*PART, NAME=SPEC1
*NODE
11, 1., 0.
12, 2., 0.
13, 2., 1.
14, 1., 1.
*NSET, NSET=FRONT
12, 13
*NSET, NSET=BACK
11, 14
*ELEMENT, TYPE=COHAX4, ELSET=EALL
11, 11, 12, 13, 14
*COHESIVE SECTION, RESPONSE=TRACTION SEPARATION, ELSET=EALL, MATERIAL=COH_DISP, STACK DIRECTION=1
*END PART
*PART, NAME=SPEC2
*NODE
11, 1., 0.
12, 2., 0.
13, 2., 1.
14, 1., 1.
*NSET, NSET=FRONT
12, 13
*NSET, NSET=BACK
11, 14
*ELEMENT, TYPE=COHAX4, ELSET=EALL
11, 11, 12, 13, 14
*COHESIVE SECTION, RESPONSE=TRACTION SEPARATION, ELSET=EALL, MATERIAL=COH_ENER, STACK DIRECTION=1
*END PART
**
**###########################################################
**INSTANCE I: PURE TENSILE (MODE 1); MODE MIX BASED ON TRACTION
**INSTANCE II: SIMPLE SHEAR (MODE 2); MODE MIX BASED ON TRACTION
**INSTANCE III: PURE TENSILE (MODE 1); MODE MIX BASED ON ENERGY
**INSTANCE IV: SIMPLE SHEAR (MODE 2); MODE MIX BASED ON ENERGY
**###########################################################
*ASSEMBLY, NAME=SPEC_ASSY
*INSTANCE, NAME=SPEC1_MODE1, PART=SPEC1
0., 0., 0.
*END INSTANCE
*INSTANCE, NAME=SPEC1_MODE2, PART=SPEC1
0., 2., 0.
*END INSTANCE
*INSTANCE, NAME=SPEC2_MODE1, PART=SPEC2
0., 4., 0.
*END INSTANCE
*INSTANCE, NAME=SPEC2_MODE2, PART=SPEC2
0., 6., 0.
*END INSTANCE
*NSET, NSET=NALL_BACK
SPEC1_MODE1.BACK, SPEC1_MODE2.BACK, SPEC2_MODE1.BACK, SPEC2_MODE2.BACK
*NSET, NSET=QA_TEST_NALL_FRONT
SPEC1_MODE1.FRONT, SPEC1_MODE2.FRONT, SPEC2_MODE1.FRONT, SPEC2_MODE2.FRONT
*ELSET, ELSET=QA_TEST_EALL
SPEC1_MODE1.EALL, SPEC1_MODE2.EALL, SPEC2_MODE1.EALL, SPEC2_MODE2.EALL
*END ASSEMBLY
**
*MATERIAL, NAME=COH_DISP
*ELASTIC,TYPE=TRACTION
4.0E6, 2.0E6, 2.0E6
*DAMAGE INITIATION, CRITERION=MAXE
0.002, 0.001
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, MODE MIX RATIO=TRACTION, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=LINEAR
0.006, 0.0
0.007, 1.0
*DENSITY
100.
**
*MATERIAL, NAME=COH_ENER
*ELASTIC,TYPE=TRACTION
4.0E6, 2.0E6, 2.0E6
*DAMAGE INITIATION, CRITERION=MAXE
0.002, 0.001
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=LINEAR
0.006, 0.0
0.007, 1.0
*DENSITY
100.
**
*BOUNDARY
NALL_BACK, 1, 2
SPEC1_MODE1.FRONT, 2, 2
SPEC1_MODE2.FRONT, 1, 1
SPEC2_MODE1.FRONT, 2, 2
SPEC2_MODE2.FRONT, 1, 1
**#####################################
**STATIC STEP
**#####################################
*STEP, INC=1000
*STATIC
0.001, 1.0, 1.0E-5, 0.01
*BOUNDARY
SPEC1_MODE1.FRONT, 1, 1, 0.01
SPEC1_MODE2.FRONT, 2, 2, 0.01
SPEC2_MODE1.FRONT, 1, 1, 0.01
SPEC2_MODE2.FRONT, 2, 2, 0.01
*OUTPUT, FIELD
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
S, E, SDEG, DMICRT, STATUS
*NODE OUTPUT, NSET=QA_TEST_NALL_FRONT
U, RF
*OUTPUT, HISTORY
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
S, E, SDEG, DMICRT, STATUS
*ENERGY OUTPUT, ELSET=SPEC1_MODE1.EALL
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC1_MODE2.EALL
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC2_MODE1.EALL
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC2_MODE2.EALL
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT
ETOTAL, ALLWK
*END STEP