** This test verifies energy type damage evolution with exponential
** 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 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 COH3D6
*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=COH3D6, ELSET=SPEC1_MODE1
11, 14, 13, 11, 18, 17, 15
12, 12, 11, 13, 16, 15, 17
*ELEMENT, TYPE=COH3D6, ELSET=SPEC1_MODE2
21, 24, 23, 21, 28, 27, 25
22, 22, 21, 23, 26, 25, 27
*ELEMENT, TYPE=COH3D6, ELSET=SPEC1_MODE3
31, 34, 33, 31, 38, 37, 35
32, 32, 31, 33, 36, 35, 37
*ELSET, ELSET=EALL_TAB
SPEC1_MODE1, SPEC1_MODE2, SPEC1_MODE3
**
*COHESIVE SECTION, RESPONSE=TRACTION SEPARATION, ELSET=EALL_TAB, MATERIAL=COH_TAB, CONTROLS=CONT
*MATERIAL, NAME=COH_TAB
*ELASTIC,TYPE=TRACTION
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=QUADS
8000., 2000., 2000.
*DAMAGE EVOLUTION, TYPE=ENERGY, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=EXPONENTIAL
24.4864, 0.0
5.8083, 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=COH3D6, ELSET=SPEC2_MODE1
41, 44, 43, 41, 48, 47, 45
42, 42, 41, 43, 46, 45, 47
*ELEMENT, TYPE=COH3D6, ELSET=SPEC2_MODE2
51, 54, 53, 51, 58, 57, 55
52, 52, 51, 53, 56, 55, 57
*ELEMENT, TYPE=COH3D6, ELSET=SPEC2_MODE3
61, 64, 63, 61, 68, 67, 65
62, 62, 61, 63, 66, 65, 67
*ELSET, ELSET=EALL_POWER
SPEC2_MODE1, SPEC2_MODE2, SPEC2_MODE3
**
*COHESIVE SECTION, RESPONSE=TRACTION SEPARATION, ELSET=EALL_POWER, MATERIAL=COH_POWER, CONTROLS=CONT
*MATERIAL, NAME=COH_POWER
*ELASTIC,TYPE=TRACTION
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=QUADS
8000., 2000., 2000.
*DAMAGE EVOLUTION, TYPE=ENERGY, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=POWER, SOFTENING=EXPONENTIAL, POWER=1.5
24.4864, 5.8083, 5.8083
*DENSITY
100.
**#####################################
**TEST CASE VII: PURE TENSILE (MODE 1)
**TEST CASE VIII: SIMPLE SHEAR (MODE 2)
**TEST CASE IX: SIMPLE SHEAR (MODE 3)  
**#####################################
*SYSTEM
0, 4, 0
*NODE, NSET=SPEC3_MODE1
71, 0., 0., 0.
72, 1., 0., 0.
73, 1., 0., 1.
74, 0., 0., 1.
75, 0., 1., 0.
76, 1., 1., 0.
77, 1., 1., 1.
78, 0., 1., 1.
*SYSTEM
2, 4, 0
*NODE, NSET=SPEC3_MODE2
81, 0., 0., 0.
82, 1., 0., 0.
83, 1., 0., 1.
84, 0., 0., 1.
85, 0., 1., 0.
86, 1., 1., 0.
87, 1., 1., 1.
88, 0., 1., 1.
*SYSTEM
4, 4, 0
*NODE, NSET=SPEC3_MODE3
91, 0., 0., 0.
92, 1., 0., 0.
93, 1., 0., 1.
94, 0., 0., 1.
95, 0., 1., 0.
96, 1., 1., 0.
97, 1., 1., 1.
98, 0., 1., 1.
*NSET, NSET=SPEC3_MODE1_TOP
75, 76, 77, 78
*NSET, NSET=SPEC3_MODE1_BOT
71, 72, 73, 74
*NSET, NSET=SPEC3_MODE2_TOP
85, 86, 87, 88
*NSET, NSET=SPEC3_MODE2_BOT
81, 82, 83, 84
*NSET, NSET=SPEC3_MODE3_TOP
95, 96, 97, 98
*NSET, NSET=SPEC3_MODE3_BOT
91, 92, 93, 94
*ELEMENT, TYPE=COH3D6, ELSET=SPEC3_MODE1
71, 74, 73, 71, 78, 77, 75
72, 72, 71, 73, 76, 75, 77
*ELEMENT, TYPE=COH3D6, ELSET=SPEC3_MODE2
81, 84, 83, 81, 88, 87, 85
82, 82, 81, 83, 86, 85, 87
*ELEMENT, TYPE=COH3D6, ELSET=SPEC3_MODE3
91, 94, 93, 91, 98, 97, 95
92, 92, 91, 93, 96, 95, 97
*ELSET, ELSET=EALL_BK
SPEC3_MODE1, SPEC3_MODE2, SPEC3_MODE3
**
*COHESIVE SECTION, RESPONSE=TRACTION SEPARATION, ELSET=EALL_BK, MATERIAL=COH_BK, CONTROLS=CONT
*MATERIAL, NAME=COH_BK
*ELASTIC,TYPE=TRACTION
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=QUADS
8000., 2000., 2000.
*DAMAGE EVOLUTION, TYPE=ENERGY, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=BK, SOFTENING=EXPONENTIAL, POWER=1.5
24.4864, 5.8083, 5.8083
*DENSITY
100.
**
*NSET, NSET=NALL_BOT
SPEC1_MODE1_BOT, SPEC1_MODE2_BOT, SPEC1_MODE3_BOT, SPEC2_MODE1_BOT, SPEC2_MODE2_BOT, SPEC2_MODE3_BOT,
SPEC3_MODE1_BOT, SPEC3_MODE2_BOT, SPEC3_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
SPEC3_MODE1_TOP, 1, 1
SPEC3_MODE1_TOP, 3, 3
SPEC3_MODE2_TOP, 2, 3
SPEC3_MODE3_TOP, 1, 2
*ELSET, ELSET=QA_TEST_EALL
EALL_TAB, EALL_POWER, EALL_BK
*NSET, NSET=QA_TEST_NALL_TOP
SPEC1_MODE1_TOP, SPEC1_MODE2_TOP, SPEC1_MODE3_TOP, SPEC2_MODE1_TOP, SPEC2_MODE2_TOP, SPEC2_MODE3_TOP,
SPEC3_MODE1_TOP, SPEC3_MODE2_TOP, SPEC3_MODE3_TOP
*SECTION CONTROLS, NAME=CONT, MAX DEGRADATION=0.9999, ELEMENT DELETION=YES
**#####################################
**DYNAMIC, EXPLICIT  STEP
**#####################################
*STEP
*DYNAMIC, EXPLICIT
, 1.0
*BOUNDARY, TYPE=VELOCITY
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
SPEC3_MODE1_TOP, 2, 2, 0.01
SPEC3_MODE2_TOP, 1, 1, 0.01
SPEC3_MODE3_TOP, 3, 3, 0.01
*OUTPUT, FIELD, NUMBER INTERVAL=20, TIME MARKS=YES
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
S, SDEG, LE, NE, DMICRT, STATUS
MMIXDMI,MMIXDME
*NODE OUTPUT, NSET=QA_TEST_NALL_TOP
U, RF
*OUTPUT, HISTORY, TIME INTERVAL=0.01
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
S, SDEG, LE, NE, DMICRT, STATUS
MMIXDMI,MMIXDME
*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, ELSET=SPEC3_MODE1
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC3_MODE2
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT, ELSET=SPEC3_MODE3
ALLSE, ALLKE, ALLDMD, ALLIE
*ENERGY OUTPUT
ETOTAL, ALLWK
*END STEP