*HEADING
  Ductile criterion for damage initiation with Mises 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, 44, 41, 42, 43
*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
*PLASTIC
200.,0.
220.,.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, 
*NSET, NSET=TOP_AX
42, 43
*NSET, NSET=BOTTOM_2D
31, 32, 41, 44
*BOUNDARY
 TOP_3D, 1, 1
 TOP_3D, 3, 3
 BOTTOM_3D, 1, 3
 TOP_2D, 1, 1
 TOP_AX, 2, 2
 BOTTOM_2D, 1, 2
*STEP
*DYNAMIC, EXPLICIT
,1.
*BOUNDARY, TYPE=VELOCITY
TOP_3D, 2, 2, 0.01
TOP_2D, 2, 2, 0.01
TOP_AX, 1, 1, 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