*HEADING
 LIMIT LOAD STUDIES, DRUCKER PRAGER, 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
*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
*STEP,INC=50
  PRESCRIBE DISPLACEMENT
*STATIC
.025,1.,,.1
*BOUNDARY
801,2,,-5.0
*MONITOR,NODE=801,DOF=2
*CONTROLS,ANALYSIS=DISCONTINUOUS
*CONTROLS, PARAMETERS=FIELD
,0.0,
*EL PRINT,ELSET=PRINTELS,FREQUENCY=10
S,E
SINV,
ENER,
E,IE
PE,
*NODE PRINT,FREQUENCY=5
U,RF
*NODE PRINT,NSET=F1
U,RF
*EL FILE,ELSET=PRINTELS,FREQUENCY=10
S,
SINV,
ENER,
IE,
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT,ELSET=PRINTELS
S,
SINV,
ENER,
IE,
*NODE FILE,NSET=F1
U,RF
*OUTPUT,FIELD
*NODE OUTPUT,NSET=F1
U,RF
*END STEP