*HEADING
  Damage test for COHESIVE BEHAVIOR with displacement type damage 
  evolution with exponential softening (exponent=2) considering 
  both energy and traction mode mix measures. CRITERION=MAXS
*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=C3D8R, ELSET=SPEC1_MODE2
21, 21, 24, 23, 22, 25, 28, 27, 26
121, 121, 124, 123, 122, 125, 128, 127, 126
*ELEMENT, TYPE=C3D8R, ELSET=SPEC1_MODE3
31, 31, 34, 33, 32, 35, 38, 37, 36
131, 131, 134, 133, 132, 135, 138, 137, 136
*ELSET, ELSET=EALL_DISP
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_DISP
*COHESIVE BEHAVIOR, Eligibility = Current Contacts 
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=QUADS
8000., 2000., 2000.
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, MODE MIX RATIO=TRACTION, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=EXPONENTIAL
0.006, 2.0, 0.0
0.007, 2.0, 1.0
*SOLID SECTION, ELSET=EALL_DISP, 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=C3D8R, ELSET=SPEC2_MODE2
51, 51, 54, 53, 52, 55, 58, 57, 56
151, 151, 154, 153, 152, 155, 158, 157, 156
*ELEMENT, TYPE=C3D8R, ELSET=SPEC2_MODE3
61, 61, 64, 63, 62, 65, 68, 67, 66
161, 161, 164, 163, 162, 165, 168, 167, 166
*ELSET, ELSET=EALL_ENER
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_ENER
*COHESIVE BEHAVIOR, Eligibility = Current Contacts 
4.0E6,2.0E6,2.0E6
*DAMAGE INITIATION, CRITERION=QUADS
8000., 2000., 2000.
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, MODE MIX RATIO=ENERGY, MIXED MODE BEHAVIOR=TABULAR, SOFTENING=EXPONENTIAL
0.006, 2.0, 0.0
0.007, 2.0, 1.0
*SOLID SECTION, ELSET=EALL_ENER, MATERIAL=MAT1
*MATERIAL, NAME=MAT1
*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
*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
*ELSET, ELSET=QA_TEST_EALL
EALL_DISP, EALL_ENER
*NSET, NSET=QA_TEST_NALL
SPEC1_MODE1, SPEC1_MODE2, SPEC1_MODE3, SPEC2_MODE1, SPEC2_MODE2, SPEC2_MODE3
**#####################################
**DYNAMIC, EXPLICIT  STEP
**#####################################
*AMPLITUDE,NAME=AMP,DEFINITION=SMOOTH STEP
0.0,0.0,1.0,1.0,2.0,1.0
*FILTER,NAME=SMOOTH,TYPE=BUTTERWORTH
50.,4
*STEP
*DYNAMIC, EXPLICIT
, 1.0
*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
*CONTACT
*CONTACT INCLUSIONS
SPEC1_MODE1_SURF1,SPEC1_MODE1_SURF2
SPEC1_MODE2_SURF1,SPEC1_MODE2_SURF2
SPEC1_MODE3_SURF1,SPEC1_MODE3_SURF2
SPEC2_MODE1_SURF1,SPEC2_MODE1_SURF2
SPEC2_MODE2_SURF1,SPEC2_MODE2_SURF2
SPEC2_MODE3_SURF1,SPEC2_MODE3_SURF2
*CONTACT PROPERTY ASSIGNMENT
SPEC1_MODE1_SURF1,SPEC1_MODE1_SURF2,COH_DISP
SPEC1_MODE2_SURF1,SPEC1_MODE2_SURF2,COH_DISP
SPEC1_MODE3_SURF1,SPEC1_MODE3_SURF2,COH_DISP
SPEC2_MODE1_SURF1,SPEC2_MODE1_SURF2,COH_ENER
SPEC2_MODE2_SURF1,SPEC2_MODE2_SURF2,COH_ENER
SPEC2_MODE3_SURF1,SPEC2_MODE3_SURF2,COH_ENER
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIME MARKS=YES
***OUTPUT, FIELD, NUMBER INTERVAL=200, TIME MARKS=YES
*NODE OUTPUT, NSET=QA_TEST_NALL
U, 
*CONTACT OUTPUT
CSDMG,CSQUADSCRT
**CFORCE,CSTRESS
*OUTPUT, HISTORY, TIME INTERVAL=0.01
*ENERGY OUTPUT
ETOTAL, ALLWK
*OUTPUT, HISTORY, FREQ=5, FILTER=SMOOTH
*CONTACT OUTPUT, SURFACE=SPEC1_MODE1_SURF1
CFT,
*CONTACT OUTPUT, SURFACE=SPEC1_MODE2_SURF1
CFT,
*CONTACT OUTPUT, SURFACE=SPEC1_MODE3_SURF1
CFT,
*CONTACT OUTPUT, SURFACE=SPEC2_MODE1_SURF1
CFT,
*CONTACT OUTPUT, SURFACE=SPEC2_MODE2_SURF1
CFT,
*CONTACT OUTPUT, SURFACE=SPEC2_MODE3_SURF1
CFT,
*END STEP