*HEADING
  Damage test for COHESIVE BEHAVIOR, displacement type
  damage evolution with exponential softening (exponent=2) considering
  both energy and traction 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
*SYSTEM
*NODE, NSET=SPEC1_MODE2
*SYSTEM
*NODE, NSET=SPEC1_MODE3
*NSET, NSET=SPEC1_MODE1_TOP
*NSET, NSET=SPEC1_MODE1_BOT
*NSET, NSET=SPEC1_MODE2_TOP
*NSET, NSET=SPEC1_MODE2_BOT
*NSET, NSET=SPEC1_MODE3_TOP
*NSET, NSET=SPEC1_MODE3_BOT
*ELEMENT, TYPE=C3D8R, ELSET=SPEC1_MODE1
*ELEMENT, TYPE=C3D8R, ELSET=SPEC1_MODE2
*ELEMENT, TYPE=C3D8R, ELSET=SPEC1_MODE3
*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
*DAMAGE INITIATION, CRITERION=MAXU
0.002, 0.001, 0.001
*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
*NODE, NSET=SPEC2_MODE1
*SYSTEM
*NODE, NSET=SPEC2_MODE2
*SYSTEM
*NODE, NSET=SPEC2_MODE3
*NSET, NSET=SPEC2_MODE1_TOP
*NSET, NSET=SPEC2_MODE1_BOT
*NSET, NSET=SPEC2_MODE2_TOP
*NSET, NSET=SPEC2_MODE2_BOT
*NSET, NSET=SPEC2_MODE3_TOP
*NSET, NSET=SPEC2_MODE3_BOT
*ELEMENT, TYPE=C3D8R, ELSET=SPEC2_MODE1
*ELEMENT, TYPE=C3D8R, ELSET=SPEC2_MODE2
*ELEMENT, TYPE=C3D8R, ELSET=SPEC2_MODE3
*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
*DAMAGE INITIATION, CRITERION=MAXU
0.002, 0.001, 0.001
*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
*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,CSMAXUCRT,
**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