*HEADING
  Cohesive elements with Drucker-Prager plasticity material model
*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
*MATERIAL,NAME=COH
*DENSITY
1000., 
*ELASTIC
2.E7,0.3
*DRUCKER PRAGER
40., 1.0, 40.
*DRUCKER PRAGER HARDENING
4.E4,0.0
*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
*STATIC
0.01, 1., 1.0E-5, 0.01
*BOUNDARY
TOP_3D, 2, 2, 0.01
TOP_2D, 2, 2, 0.01
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=QA_TEST_EALL
S, LE, NE, EE, PE, ER, SINV
*NODE OUTPUT,NSET=QA_TEST_NALL
U, RF
*OUTPUT,HISTORY
*ENERGY OUTPUT, VARIABLE=PRESELECT
*END STEP