*HEADING
  Ductile criterion for damage initiation with Drucker-Prager plasticity
*PREPRINT, MODEL=YES
*NODE,NSET=COH3D8
11, 0., 0., 0.
12, 0., 0., 1.
13, 1., 0., 1.
14, 1., 0., 0.
15, 0., 1., 0.
16, 0., 1., 1.
17, 1., 1., 1.
18, 1., 1., 0.
*SYSTEM
0., 2., 0.
*NODE,NSET=COH3D6
21, 0., 0., 0.
22, 0., 0., 1.
23, 1., 0., 1.
24, 1., 0., 0.
25, 0., 1., 0.
26, 0., 1., 1.
27, 1., 1., 1.
28, 1., 1., 0.
*SYSTEM
0., 4., 0.
*NODE,NSET=COH2D4
31, 0., 0., 0.
32, 1., 0., 0.
33, 1., 1., 0.
34, 0., 1., 0.
*SYSTEM
0., 6., 0.
*NODE,NSET=COHAX4
41, 1., 0., 0.
42, 2., 0., 0.
43, 2., 1., 0.
44, 1., 1., 0.
*ELEMENT,TYPE=COH3D8, ELSET=COH3D8
11, 11, 12, 13, 14, 15, 16, 17, 18
*ELEMENT,TYPE=COH3D6, ELSET=COH3D6
21, 21, 22, 23, 25, 26, 27
22, 24, 21, 23, 28, 25, 27
*ELEMENT,TYPE=COH2D4, ELSET=COH2D4
31, 31, 32, 33, 34
*ELEMENT,TYPE=COHAX4, ELSET=COHAX4
41, 41, 42, 43, 44
*NSET, NSET=NALL
COH3D8, COH3D6, COH2D4, COHAX4
*ELSET, ELSET=EALL
COH3D8, COH3D6, COH2D4, COHAX4
*COHESIVE SECTION, THICKNESS=SPECIFIED, MATERIAL=COH, ELSET=EALL, RESPONSE=CONTINUUM, CONTROLS=CONT
*MATERIAL,NAME=COH
*DENSITY
10.
*ELASTIC
200.E3,0.3
*DRUCKER PRAGER
40., 1., 40.
*DRUCKER PRAGER HARDENING
200., 0.
220., 0.0009
*DAMAGE INITIATION, CRITERION=DUCTILE
0.0004
*DAMAGE EVOLUTION,TYPE=DISPLACEMENT,SOFTENING=LINEAR
0.0026
*SECTION CONTROLS, NAME=CONT, MAX DEGRADATION=0.95, ELEMENT DELETION=YES
*NSET, NSET=QA_TEST_NALL
NALL
*ELSET, ELSET=QA_TEST_EALL
EALL
*NSET, NSET=TOP_3D
15, 16, 17, 18,
25, 26, 27, 28
*NSET, NSET=BOTTOM_3D
11, 12, 13, 14, 
21, 22, 23, 24
*NSET, NSET=TOP_2D
33, 34, 
43, 44
*NSET, NSET=BOTTOM_2D
31, 32, 
41, 42
*BOUNDARY
 TOP_3D, 1, 1
 TOP_3D, 3, 3
 BOTTOM_3D, 1, 3
 TOP_2D, 1, 1
 BOTTOM_2D, 1, 2
*STEP
*DYNAMIC, EXPLICIT
,1.
*BOUNDARY, TYPE=VELOCITY
TOP_3D, 2, 2, 0.01
TOP_2D, 2, 2, 0.01
*OUTPUT,FIELD, NUMBER INTERVAL=20, TIME MARKS=YES
*ELEMENT OUTPUT,ELSET=QA_TEST_EALL
S, LE, NE, PE, PEEQ,
STATUS,SDEG,DMICRT
*NODE OUTPUT,NSET=QA_TEST_NALL
U, RF
*OUTPUT,HISTORY, TIME INTERVAL=0.01
*ELEMENT OUTPUT,ELSET=QA_TEST_EALL
S, LE, NE, PE, PEEQ,
STATUS,SDEG,DMICRT
*ENERGY OUTPUT, VARIABLE=PRESELECT
*END STEP