** This test verifies energy type damage evolution with linear
** softening considering energy 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_TAB, 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_POWER, STACK DIRECTION=1
*END PART
*PART, NAME=SPEC3
*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_BK, STACK DIRECTION=1
*END PART
**
**###########################################################
**INSTANCE I: PURE TENSILE (MODE 1); MODE MIX IN TABULAR FORM
**INSTANCE II: SIMPLE SHEAR (MODE 2); MODE MIX IN TABULAR FORM 
**INSTANCE III: PURE TENSILE (MODE 1); MODE MIX BASED ON POWER LAW
**INSTANCE IV: SIMPLE SHEAR (MODE 2); MODE MIX BASED ON POWER LAW
**INSTANCE V: PURE TENSILE (MODE 1); MODE MIX BASED ON BK LAW
**INSTANCE VI: SIMPLE SHEAR (MODE 2); MODE MIX BASED ON BK LAW
**###########################################################
*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
*INSTANCE, NAME=SPEC3_MODE1, PART=SPEC3
0., 8., 0.
*END INSTANCE
*INSTANCE, NAME=SPEC3_MODE2, PART=SPEC3
0., 10., 0.
*END INSTANCE
*NSET, NSET=NALL_BACK
SPEC1_MODE1.BACK, SPEC1_MODE2.BACK, SPEC2_MODE1.BACK, SPEC2_MODE2.BACK, SPEC3_MODE1.BACK, SPEC3_MODE2.BACK
*NSET, NSET=QA_TEST_NALL_FRONT
SPEC1_MODE1.FRONT, SPEC1_MODE2.FRONT, SPEC2_MODE1.FRONT, SPEC2_MODE2.FRONT, SPEC3_MODE1.FRONT, SPEC3_MODE2.FRONT
*ELSET, ELSET=QA_TEST_EALL
SPEC1_MODE1.EALL, SPEC1_MODE2.EALL, SPEC2_MODE1.EALL, SPEC2_MODE2.EALL, SPEC3_MODE1.EALL, SPEC3_MODE2.EALL
*END ASSEMBLY
**
*MATERIAL, NAME=COH_TAB
*ELASTIC,TYPE=TRACTION
4.0E6, 2.0E6, 2.0E6
*DAMAGE INITIATION, CRITERION=MAXE
0.002, 0.001
*DAMAGE EVOLUTION, TYPE=ENERGY, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=LINEAR
32.0, 0.0
8.0, 1.0
*DENSITY
100.
**
*MATERIAL, NAME=COH_POWER
*ELASTIC,TYPE=TRACTION
4.0E6, 2.0E6, 2.0E6
*DAMAGE INITIATION, CRITERION=MAXE
0.002, 0.001
*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.
**
*MATERIAL, NAME=COH_BK
*ELASTIC,TYPE=TRACTION
4.0E6, 2.0E6, 2.0E6
*DAMAGE INITIATION, CRITERION=MAXE
0.002, 0.001
*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.
**
*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
SPEC3_MODE1.FRONT, 2, 2
SPEC3_MODE2.FRONT, 1, 1
**#####################################
**DYNAMIC, EXPLICIT  STEP
**#####################################
*STEP
*DYNAMIC, EXPLICIT
, 1.0
*BOUNDARY, TYPE=VELOCITY
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
SPEC3_MODE1.FRONT, 1, 1, 0.01
SPEC3_MODE2.FRONT, 2, 2, 0.01
*OUTPUT, FIELD, NUMBER INTERVAL=20, TIME MARKS=YES
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
S, SDEG, LE, NE, DMICRT, STATUS
*NODE OUTPUT, NSET=QA_TEST_NALL_FRONT
U, RF
*OUTPUT, HISTORY, TIME INTERVAL=0.01
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
S, SDEG, LE, NE, 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, ELSET=SPEC3_MODE1.EALL
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC3_MODE2.EALL
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT
ETOTAL, ALLWK
*END STEP