*HEADING
DRUCKER-PRAGER/LINEAR ELASTICITY, CRITERION=EXP W/ TEST DATA, UC TEST, 3D
WITH STATIC LINEAR PERTURBATIONS
*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
*SOLID SECTION,ELSET=ALLE,MATERIAL=ALLE
*MATERIAL,NAME=ALLE
*ELASTIC
300.E3,.3
*DRUCKER PRAGER,SHEAR CRITERION=EXPONENT,TEST DATA 
,,,40.0
*TRIAXIAL TEST DATA, B=1.0
0,3377.199
0,-6000.00
*DRUCKER PRAGER HARDENING
6.E3,0.
9.E3,.02
11.E3,.063333
12.E3,.11
12.E3,1.
*STEP,INC=10
STEP 1 - NONLINEAR
*STATIC,DIRECT
1.,10.
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,3,,-.1
5,3,,-.1
6,3,,-.1
8,3,,-.1
*EL PRINT
S, 
SINV, 
E, 
PE, 
EE, 
*NODE PRINT
U, 
RF, 
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,E,PE
*END STEP
*STEP,PERTURBATION,INC=1
STEP 2 - LINEAR PERTURBATION
*STATIC,DIRECT
.1,.1
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,3,,-.01
5,3,,-.01
6,3,,-.01
8,3,,-.01
*EL PRINT
S, 
SINV, 
E, 
*NODE PRINT
U, 
RF, 
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,E
*END STEP
*STEP,INC=10
STEP 3 - REST OF NONLINEAR
*STATIC,DIRECT
1.,10.
*BOUNDARY
7,3,,-.2
5,3,,-.2
6,3,,-.2
8,3,,-.2
*END STEP
*STEP,PERTURBATION,INC=1
STEP 4 - LINEAR PERTURBATION
*STATIC,DIRECT
.1,.1
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,3,,-.02
5,3,,-.02
6,3,,-.02
8,3,,-.02
*EL PRINT
S, 
SINV, 
E, 
*NODE PRINT
U, 
RF, 
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,E
*END STEP