*HEADING
 LIMIT LOAD STUDIES, DRUCKER PRAGER CREEP AND PLASTICITY, DILATANT FLOW -
*RESTART,WRITE,FREQUENCY=1
*NODE
*NGEN,NSET=BASE
*NSET,NSET=F1
*NSET,NSET=F2,GENERATE
*NGEN,NSET=CENTER
*NGEN,NSET=TOP
*NFILL
BASE,TOP,8,100
*NGEN,NSET=FAR
*ELEMENT,TYPE=CPE8R
*ELGEN,ELSET=ALL
*ELSET,ELSET=PRINTELS
*SOLID SECTION,ELSET=ALL,MATERIAL=A1
*MATERIAL,NAME= A1
*ELASTIC
30000.,0.3
*DRUCKER PRAGER,SHEAR CRITERION=LINEAR
*DRUCKER PRAGER HARDENING
*DRUCKER PRAGER CREEP, LAW=TIME
*ELEMENT,TYPE=CINPE5R
*ELGEN,ELSET=FAR
*SOLID SECTION,ELSET=FAR,MATERIAL=A2
*MATERIAL,NAME= A2
*ELASTIC
30000.,0.3
*EQUATION
F2,2,1.,801,2,-1.
*BOUNDARY
CENTER,1
F2,1
BASE,1,2
*AMPLITUDE, NAME=RAMP
0.,0.,1.,1.
*STEP,INC=50, UNSYMM=YES
  PRESCRIBE DISPLACEMENT
*VISCO, CETOL=0.01
.025,1.,,.1
*BOUNDARY, AMP=RAMP
801,2,,-5.0
*MONITOR,NODE=801,DOF=2
*CONTROLS,ANALYSIS=DISCONTINUOUS
*EL PRINT,ELSET=PRINTELS,FREQUENCY=10
S,E
SINV,
ENER,
E,IE
PE,
CE,
*NODE PRINT,FREQUENCY=5
U,RF
*NODE PRINT,NSET=F1
U,RF
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT,ELSET=PRINTELS
S,
SINV,
ENER,
IE,
PE,
CE,
*OUTPUT,FIELD
*NODE OUTPUT,NSET=F1
U,RF
*END STEP