*HEADING
  Damage test for COHESIVE BEHAVIOR for energy type damage evolution 
  with exponential softening (exponent=2) considering energy mode mix 
  measures. CRITERION=MAXU
*PREPRINT, MODEL=YES, HISTORY=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.
111, 0., 1., 0.
112, 1., 1., 0.
113, 1., 1, 1.
114, 0., 1., 1.
115, 0., 2., 0.
116, 1., 2., 0.
117, 1., 2., 1.
118, 0., 2., 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.
121, 0., 1., 0.
122, 1., 1., 0.
123, 1., 1., 1.
124, 0., 1., 1.
125, 0., 2., 0.
126, 1., 2., 0.
127, 1., 2., 1.
128, 0., 2., 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.
131, 0., 1., 0.
132, 1., 1., 0.
133, 1., 1., 1.
134, 0., 1., 1.
135, 0., 2., 0.
136, 1., 2., 0.
137, 1., 2., 1.
138, 0., 2., 1.
*NSET, NSET=SPEC1_MODE1_TOP
115, 116, 117, 118
*NSET, NSET=SPEC1_MODE1_BOT
11, 12, 13, 14
*NSET, NSET=SPEC1_MODE2_TOP
125, 126, 127, 128
*NSET, NSET=SPEC1_MODE2_BOT
21, 22, 23, 24
*NSET, NSET=SPEC1_MODE3_TOP
135, 136, 137, 138
*NSET, NSET=SPEC1_MODE3_BOT
31, 32, 33, 34
*ELEMENT, TYPE=C3D8R, ELSET=SPEC1_MODE1
11, 11, 14, 13, 12, 15, 18, 17, 16
111, 111, 114, 113, 112, 115, 118, 117, 116
*ELEMENT, TYPE=C3D8, ELSET=SPEC1_MODE2
21, 21, 24, 23, 22, 25, 28, 27, 26
121, 121, 124, 123, 122, 125, 128, 127, 126
*ELEMENT, TYPE=C3D8, ELSET=SPEC1_MODE3
31, 31, 34, 33, 32, 35, 38, 37, 36
131, 131, 134, 133, 132, 135, 138, 137, 136
*ELSET, ELSET=EALL_TAB
SPEC1_MODE1, SPEC1_MODE2, SPEC1_MODE3
*SURFACE, NAME=SPEC1_MODE1_SURF1
11,S2
*SURFACE, NAME=SPEC1_MODE1_SURF2
111,S1
*SURFACE, NAME=SPEC1_MODE2_SURF1
21,S2
*SURFACE, NAME=SPEC1_MODE2_SURF2
121,S1
*SURFACE, NAME=SPEC1_MODE3_SURF1
31,S2
*SURFACE, NAME=SPEC1_MODE3_SURF2
131,S1
**
*SURFACE INTERACTION,NAME=COH_TAB
*COHESIVE BEHAVIOR, Eligibility = Current Contacts 
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=MAXU
0.002,0.001,0.001
*DAMAGE EVOLUTION, TYPE=ENERGY, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=EXPONENTIAL
24.4864, 0.0
5.8083, 1.0
*DAMAGE STABILIZATION
5.E-5
*SOLID SECTION, ELSET=EALL_TAB, MATERIAL=MAT
*MATERIAL, NAME=MAT
*ELASTIC
4.E+6,
*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., 1., 0.
42, 1., 1., 0.
43, 1., 1., 1.
44, 0., 1., 1.
45, 0., 2., 0.
46, 1., 2., 0.
47, 1., 2., 1.
48, 0., 2., 1.
141, 0., 2., 0.
142, 1., 2., 0.
143, 1., 2., 1.
144, 0., 2., 1.
145, 0., 3., 0.
146, 1., 3., 0.
147, 1., 3., 1.
148, 0., 3., 1.
*SYSTEM
2, 2, 0
*NODE, NSET=SPEC2_MODE2
51, 0., 1., 0.
52, 1., 1., 0.
53, 1., 1., 1.
54, 0., 1., 1.
55, 0., 2., 0.
56, 1., 2., 0.
57, 1., 2., 1.
58, 0., 2., 1.
151, 0., 2., 0.
152, 1., 2., 0.
153, 1., 2., 1.
154, 0., 2., 1.
155, 0., 3., 0.
156, 1., 3., 0.
157, 1., 3., 1.
158, 0., 3., 1.
*SYSTEM
4, 2, 0
*NODE, NSET=SPEC2_MODE3
61, 0., 1., 0.
62, 1., 1., 0.
63, 1., 1., 1.
64, 0., 1., 1.
65, 0., 2., 0.
66, 1., 2., 0.
67, 1., 2., 1.
68, 0., 2., 1.
161, 0., 2., 0.
162, 1., 2., 0.
163, 1., 2., 1.
164, 0., 2., 1.
165, 0., 3., 0.
166, 1., 3., 0.
167, 1., 3., 1.
168, 0., 3., 1.
*NSET, NSET=SPEC2_MODE1_TOP
145, 146, 147, 148
*NSET, NSET=SPEC2_MODE1_BOT
41, 42, 43, 44
*NSET, NSET=SPEC2_MODE2_TOP
155, 156, 157, 158
*NSET, NSET=SPEC2_MODE2_BOT
51, 52, 53, 54
*NSET, NSET=SPEC2_MODE3_TOP
165, 166, 167, 168
*NSET, NSET=SPEC2_MODE3_BOT
61, 62, 63, 64
*ELEMENT, TYPE=C3D8R, ELSET=SPEC2_MODE1
41, 41, 44, 43, 42, 45, 48, 47, 46
141, 141, 144, 143, 142, 145, 148, 147, 146
*ELEMENT, TYPE=C3D8, ELSET=SPEC2_MODE2
51, 51, 54, 53, 52, 55, 58, 57, 56
151, 151, 154, 153, 152, 155, 158, 157, 156
*ELEMENT, TYPE=C3D8, ELSET=SPEC2_MODE3
61, 61, 64, 63, 62, 65, 68, 67, 66
161, 161, 164, 163, 162, 165, 168, 167, 166
*ELSET, ELSET=EALL_POWER
SPEC2_MODE1, SPEC2_MODE2, SPEC2_MODE3
*SURFACE, NAME=SPEC2_MODE1_SURF1
41,S2
*SURFACE, NAME=SPEC2_MODE1_SURF2
141,S1
*SURFACE, NAME=SPEC2_MODE2_SURF1
51,S2
*SURFACE, NAME=SPEC2_MODE2_SURF2
151,S1
*SURFACE, NAME=SPEC2_MODE3_SURF1
61,S2
*SURFACE, NAME=SPEC2_MODE3_SURF2
161,S1
**
*SURFACE INTERACTION, NAME=COH_POWER
*COHESIVE BEHAVIOR, Eligibility = Current Contacts 
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=MAXU
0.002,0.001,0.001
*DAMAGE EVOLUTION, TYPE=ENERGY, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=POWER, SOFTENING=EXPONENTIAL,POWER=1.5
24.4864, 5.8083, 5.8083
*DAMAGE STABILIZATION
5.E-5
*SOLID SECTION, ELSET=EALL_POWER, MATERIAL=MAT1
*MATERIAL, NAME=MAT1
*ELASTIC
4.E+6,
*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., 2., 0.
72, 1., 2., 0.
73, 1., 2., 1.
74, 0., 2., 1.
75, 0., 3., 0.
76, 1., 3., 0.
77, 1., 3., 1.
78, 0., 3., 1.
171, 0., 3., 0.
172, 1., 3., 0.
173, 1., 3., 1.
174, 0., 3., 1.
175, 0., 4., 0.
176, 1., 4., 0.
177, 1., 4., 1.
178, 0., 4., 1.
*SYSTEM
2, 4, 0
*NODE, NSET=SPEC3_MODE2
81, 0., 2., 0.
82, 1., 2., 0.
83, 1., 2., 1.
84, 0., 2., 1.
85, 0., 3., 0.
86, 1., 3., 0.
87, 1., 3., 1.
88, 0., 3., 1.
181, 0., 3., 0.
182, 1., 3., 0.
183, 1., 3., 1.
184, 0., 3., 1.
185, 0., 4., 0.
186, 1., 4., 0.
187, 1., 4., 1.
188, 0., 4., 1.
*SYSTEM
4, 4, 0
*NODE, NSET=SPEC3_MODE3
91, 0., 2., 0.
92, 1., 2., 0.
93, 1., 2., 1.
94, 0., 2., 1.
95, 0., 3., 0.
96, 1., 3., 0.
97, 1., 3., 1.
98, 0., 3., 1.
191, 0., 3., 0.
192, 1., 3., 0.
193, 1., 3., 1.
194, 0., 3., 1.
195, 0., 4., 0.
196, 1., 4., 0.
197, 1., 4., 1.
198, 0., 4., 1.
*NSET, NSET=SPEC3_MODE1_TOP
175, 176, 177, 178
*NSET, NSET=SPEC3_MODE1_BOT
71, 72, 73, 74
*NSET, NSET=SPEC3_MODE2_TOP
185, 186, 187, 188
*NSET, NSET=SPEC3_MODE2_BOT
81, 82, 83, 84
*NSET, NSET=SPEC3_MODE3_TOP
195, 196, 197, 198
*NSET, NSET=SPEC3_MODE3_BOT
91, 92, 93, 94
*ELEMENT, TYPE=C3D8R, ELSET=SPEC3_MODE1
71, 71, 74, 73, 72, 75, 78, 77, 76
171, 171, 174, 173, 172, 175, 178, 177, 176
*ELEMENT, TYPE=C3D8, ELSET=SPEC3_MODE2
81, 81, 84, 83, 82, 85, 88, 87, 86
181, 181, 184, 183, 182, 185, 188, 187, 186
*ELEMENT, TYPE=C3D8, ELSET=SPEC3_MODE3
91, 91, 94, 93, 92, 95, 98, 97, 96
191, 191, 194, 193, 192, 195, 198, 197, 196
*ELSET, ELSET=EALL_BK
SPEC3_MODE1, SPEC3_MODE2, SPEC3_MODE3
**
*SURFACE, NAME=SPEC3_MODE1_SURF1
71,S2
*SURFACE, NAME=SPEC3_MODE1_SURF2
171,S1
*SURFACE, NAME=SPEC3_MODE2_SURF1
81,S2
*SURFACE, NAME=SPEC3_MODE2_SURF2
181,S1
*SURFACE, NAME=SPEC3_MODE3_SURF1
91,S2
*SURFACE, NAME=SPEC3_MODE3_SURF2
191,S1
**
*SURFACE INTERACTION, NAME=COH_BK
*COHESIVE BEHAVIOR, Eligibility = Current Contacts 
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=MAXU
0.002,0.001,0.001
*DAMAGE EVOLUTION, TYPE=ENERGY, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=BK, SOFTENING=EXPONENTIAL, POWER=1.5
24.4864, 5.8083, 5.8083
*DAMAGE STABILIZATION
5.E-5
*SOLID SECTION, ELSET=EALL_BK, MATERIAL=MAT2
*MATERIAL, NAME=MAT2
*ELASTIC
4.E+6,
*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
SPEC1_MODE1, SPEC1_MODE2, SPEC1_MODE3, SPEC2_MODE1, SPEC2_MODE2, SPEC2_MODE3,
SPEC3_MODE1, SPEC3_MODE2, SPEC3_MODE3
*CONTACT PAIR, INTERACTION=COH_TAB
SPEC1_MODE1_SURF1,SPEC1_MODE1_SURF2
SPEC1_MODE2_SURF1,SPEC1_MODE2_SURF2
SPEC1_MODE3_SURF1,SPEC1_MODE3_SURF2
*CONTACT PAIR, INTERACTION=COH_POWER
SPEC2_MODE1_SURF1,SPEC2_MODE1_SURF2
SPEC2_MODE2_SURF1,SPEC2_MODE2_SURF2
SPEC2_MODE3_SURF1,SPEC2_MODE3_SURF2
*CONTACT PAIR, INTERACTION=COH_BK
SPEC3_MODE1_SURF1,SPEC3_MODE1_SURF2
SPEC3_MODE2_SURF1,SPEC3_MODE2_SURF2
SPEC3_MODE3_SURF1,SPEC3_MODE3_SURF2
*AMPLITUDE,NAME=AMP,DEFINITION=SMOOTH STEP
0.0,0.0,1.0,1.0,2.0,1.0
*STEP, inc=1000
*STATIC
0.005, 1.0,,0.005
*BOUNDARY, TYPE=DISPLACEMENT, AMPLITUDE=AMP
SPEC1_MODE1_TOP, 2, 2, 0.02
SPEC1_MODE2_TOP, 1, 1, 0.02
SPEC1_MODE3_TOP, 3, 3, 0.02
SPEC2_MODE1_TOP, 2, 2, 0.02
SPEC2_MODE2_TOP, 1, 1, 0.02
SPEC2_MODE3_TOP, 3, 3, 0.02
SPEC3_MODE1_TOP, 2, 2, 0.02
SPEC3_MODE2_TOP, 1, 1, 0.02
SPEC3_MODE3_TOP, 3, 3, 0.02
*OUTPUT, FIELD, FREQ=10
***OUTPUT, FIELD
**NODE OUTPUT, NSET=QA_TEST_NALL
**U, 
*CONTACT OUTPUT
CSDMG,CSMAXUCRT
**CFORCE,CSTRESS
**OUTPUT, HISTORY
**ENERGY OUTPUT
**ETOTAL, ALLWK
** These output produce constant QC diffs on frame timestamps in parallel runs
***OUTPUT, HISTORY
***CONTACT OUTPUT, SECONDARY=SPEC1_MODE1_SURF1,MAIN=SPEC1_MODE1_SURF2
**CFT,
***CONTACT OUTPUT, SECONDARY=SPEC1_MODE2_SURF1,MAIN=SPEC1_MODE2_SURF2
**CFT,
***CONTACT OUTPUT, SECONDARY=SPEC1_MODE3_SURF1,MAIN=SPEC1_MODE3_SURF2
**CFT,
***CONTACT OUTPUT, SECONDARY=SPEC2_MODE1_SURF1,MAIN=SPEC2_MODE1_SURF2
**CFT,
****CONTACT OUTPUT, SECONDARY=SPEC2_MODE2_SURF1,MAIN=SPEC2_MODE2_SURF2
****CFT,
****CONTACT OUTPUT, SECONDARY=SPEC2_MODE3_SURF1,MAIN=SPEC2_MODE3_SURF2
****CFT,
***CONTACT OUTPUT, SECONDARY=SPEC3_MODE1_SURF1,MAIN=SPEC3_MODE1_SURF2
**CFT,
****CONTACT OUTPUT, SECONDARY=SPEC3_MODE2_SURF1,MAIN=SPEC3_MODE2_SURF2
****CFT,
****CONTACT OUTPUT, SECONDARY=SPEC3_MODE3_SURF1,MAIN=SPEC3_MODE3_SURF2
****CFT,
*END STEP