** This test verifies displacement type damage evolution with tabular
** 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 = 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 1 direction
**	7. Thickness = 1.0 => Nominal strain = Separation
** Test cases I - II will give the same results as:
**                COHAX4_qdS_damDisp_softLin_xpl.inp
** Test cases III - IV will give the same results as:
**                COHAX4_qdS_damDisp_softExp_xpl.inp
*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=QUADS
8000., 2000.
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, MODE MIX RATIO=TRACTION, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=TABULAR
0.0000, 0.000, 0.0 
0.4444, 0.001, 0.0
0.6667, 0.002, 0.0
0.8000, 0.003, 0.0
0.8889, 0.004, 0.0
0.9524, 0.005, 0.0
1.0000, 0.006, 0.0
0.0000, 0.000, 1.0 
0.5714, 0.001, 1.0
0.7619, 0.002, 1.0
0.8571, 0.003, 1.0
0.9143, 0.004, 1.0
0.9524, 0.005, 1.0
0.9796, 0.006, 1.0
1.0000, 0.007, 1.0
*DENSITY
100.
**
*MATERIAL, NAME=COH_ENER
*ELASTIC,TYPE=TRACTION
4.0E6, 2.0E6, 2.0E6
*DAMAGE INITIATION, CRITERION=QUADS
8000., 2000.
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=TABULAR
0.0000, 0.000, 0.0
0.5519, 0.001, 0.0
0.7814, 0.002, 0.0
0.8924, 0.003, 0.0
0.9506, 0.004, 0.0
0.9823, 0.005, 0.0
1.0000, 0.006, 0.0
0.0000, 0.000, 1.0
0.6437, 0.001, 1.0
0.8345, 0.002, 1.0
0.9164, 0.003, 1.0
0.9575, 0.004, 1.0
0.9799, 0.005, 1.0
0.9926, 0.006, 1.0
1.0000, 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
**#####################################
**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
*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
ETOTAL, ALLWK
*END STEP