*HEADING
 Limit Load Studies
 Mohr-Coulomb, Plane Strain, Dilatant flow 
*RESTART,WRITE,FREQUENCY=1
*NODE
   1,   0.
   7,  60.
  13, 180.
  15, 228.
  19, 348.
  20, 696.
 801,   0., 144.
 807,  60., 144.
 813, 180., 144.
 815, 228., 144.
 819, 348., 144.
 820, 696., 144.
*NGEN,NSET=BASE
 1,7
 7,13
 13,15
 15,19
*NSET,NSET=F1
 801,
*NSET,NSET=F2,GENERATE
 802,807
*NGEN,NSET=CENTER
 1,801,100
*NGEN,NSET=TOP
 801,807
 807,813
 813,815
 815,819
*NFILL
 BASE,TOP,8,100
*NGEN,NSET=FAR
 20,820,200
*ELEMENT,TYPE=CPE8R
 1,1,3,203,201,2,103,202,101
*ELGEN,ELSET=ALL
 1,4,200,1,9,2,10
*ELSET,ELSET=PRINTELS
 1,2,3,4
*SOLID SECTION,ELSET=ALL,MATERIAL=A1
*MATERIAL,NAME=A1
*ELASTIC
 30000.,0.3
** Dilation angles larger than 17.831 deg. my result in negative plastic work
** under high confining stress. 
*MOHR COULOMB
 20.,17.83
*MOHR COULOMB HARDENING
 10.,
*ELEMENT,TYPE=CINPE5R
 101,219,19,20,220,119
*ELGEN,ELSET=FAR
 101,4,200,1
*SOLID SECTION,ELSET=FAR,MATERIAL=A2
*MATERIAL,NAME=A2
*ELASTIC
 30000.,0.3
*EQUATION
 2, 
 F2,2,1.,801,2,-1.
*BOUNDARY
 CENTER,1
 F2,1
 BASE,1,2
*STEP,INC=50,UNSYMM=YES
 PRESCRIBE DISPLACEMENT
*STATIC
.025,1.,,.1
*BOUNDARY
 801,2,,-5.0
*MONITOR,NODE=801,DOF=2
*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