*HEADING
DRUCKER-PRAGER/LINEAR ELASTICITY, UNIAXIAL COMPRESSION TEST, 3D
*NODE,NSET=ALLN
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.
*ELEMENT,TYPE=C3D8,ELSET=ALLE
1,1,2,3,4,5,6,7,8
*Distribution Table, name=distTableElastic
 modulus, ratio
*Distribution, name=distElastic, table=distTableElastic, location=element
 , 200.e3, 0.25
 alle, 300.0e3, 0.3
*SOLID SECTION,ELSET=ALLE,MATERIAL=ALLE
*MATERIAL,NAME=ALLE
*ELASTIC
distElastic,
*DRUCKER PRAGER
40.,1.,40.
*DRUCKER PRAGER HARDENING
6.E3,0.
9.E3,.02
11.E3,.063333
12.E3,.11
12.E3,1.
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
*STEP,INC=20
*STATIC,DIRECT
1.,20.
*BOUNDARY
7,3,,-.2
5,3,,-.2
6,3,,-.2
8,3,,-.2
*EL PRINT
S, 
SINV, 
E, 
PE, 
EE, 
*NODE PRINT
U, 
RF, 
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,E,PE
*END STEP