*HEADING
DRUCKER-PRAGER/LINEAR ELASTICITY, SHEAR CRITERION=EXP, BIAXIAL TENSION, 3D
*PREPRINT,MODEL=YES
*RESTART,WRITE
*NODE,NSET=ALL
1,0.,0.,0.
2,1.,0.,0.
3,1.,1.,0.
4,0.,1.,0.
5,0.,0.,1.
6,1.,0.,1.
7,1.,1.,1.
8,0.,1.,1.
**
11,0.,0.,2.
12,1.,0.,2.
13,1.,1.,2.
14,0.,1.,2.
15,0.,0.,3.
16,1.,0.,3.
17,1.,1.,3.
18,0.,1.,3.
**
21,0.,0.,4.
22,1.,0.,4.
23,1.,1.,4.
24,0.,1.,4.
25,0.,0.,5.
26,1.,0.,5.
27,1.,1.,5.
28,0.,1.,5.
*NSET,NSET=FACE1
1,2,3,4
11,12,13,14
21,22,23,24
*NSET,NSET=FACE2
5,6,7,8
15,16,17,18
25,26,27,28
*NSET,NSET=FACE3
1,2,5,6
11,12,15,16
21,22,25,26
*NSET,NSET=FACE4
2,3,6,7
12,13,16,17
22,23,26,27
*NSET,NSET=FACE5
3,4,7,8
13,14,17,18
23,24,27,28
*NSET,NSET=FACE6
4,1,8,5
14,11,18,15
24,21,28,25
*ELEMENT,TYPE=C3D8R,ELSET=LIN
1,1,2,3,4,5,6,7,8
*ELEMENT,TYPE=C3D8R,ELSET=HYP
2,11,12,13,14,15,16,17,18
*ELEMENT,TYPE=C3D8R,ELSET=EXP
3,21,22,23,24,25,26,27,28
*ELSET, ELSET=ALLE
 LIN,HYP,EXP
*SOLID SECTION,ELSET=LIN,MATERIAL=DPLIN
*SOLID SECTION,ELSET=HYP,MATERIAL=DPHYP
*SOLID SECTION,ELSET=EXP,MATERIAL=DPEXP
**
*MATERIAL,NAME=DPLIN
*DENSITY
1.0E-2
*ELASTIC
300.E3,.3
*DRUCKER PRAGER
40.,1.,40.
*DRUCKER PRAGER HARDENING
6.E3,0.
9.E3,.02
11.E3,.063333
12.E3,.11
12.E3,1.
**
*MATERIAL,NAME=DPHYP
*DENSITY
1.0E-2
*ELASTIC
300.E3,.3
*DRUCKER PRAGER,SHEAR CRITERION=HYPERBOLIC
40.,5150.0,,40.
*DRUCKER PRAGER HARDENING 
6.E3,0.
9.E3,.02
11.E3,.063333
12.E3,.11
12.E3,1.
**
*MATERIAL,NAME=DPEXP
*DENSITY
1.0E-2
*ELASTIC
300.E3,.3
*DRUCKER PRAGER,SHEAR CRITERION=EXPONENT
1.19175,1.0,,40.0
*DRUCKER PRAGER HARDENING 
6.E3,0.
9.E3,.02
11.E3,.063333
12.E3,.11
12.E3,1.
**
*STEP,NLGEOM,INC=20
UNIAXIAL TENSION
*STATIC
0.1,1.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,.02
*EL PRINT,F=10
S, 
E, 
PEEQ,
*NODE PRINT,F=10
U,RF
*OUTPUT,FIELD,FREQ=10
*NODE OUTPUT
U,RF
*ELEMENT OUTPUT
S,PEEQ,PE
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=ALLE
S,PEEQ,PE,PENER
*END STEP