*HEADING
 Limit Load Studies
 Mohr-Coulomb, Plane Strain, Dilatant flow 
*RESTART,WRITE
*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=QA_TEST_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=CPE4R
 1,1,3,203,201
*ELGEN,ELSET=ALL
 1,4,200,1,9,2,10
*ELSET,ELSET=QA_TEST_PRINTELS
 1,2,3,4
*SOLID SECTION,ELSET=ALL,MATERIAL=A1
*MATERIAL,NAME=A1
*density 
 1.0e-6
*ELASTIC
 30000.,0.3
** Dilation angles larger than 17.831 deg. may result in negative plastic work
** under high confining stress. 
*MOHR COULOMB
 20.,17.83
*MOHR COULOMB HARDENING
 10.,
*ELEMENT,TYPE=CINPE4
 101,219,19,20,220
*ELGEN,ELSET=FAR
 101,4,200,1
*SOLID SECTION,ELSET=FAR,MATERIAL=A2
*MATERIAL,NAME=A2
*density 
 1.0e-6
*ELASTIC
 30000.,0.3
*EQUATION
 2, 
 F2,2,1.,801,2,-1.
*BOUNDARY
 CENTER,1
 F2,1
 BASE,1,2
*STEP
 PRESCRIBE DISPLACEMENT
*dynamic,explicit
 ,1.
*BOUNDARY, type=velocity
 801,2,,-5.0
*MONITOR,NODE=801,DOF=2
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT,ELSET=QA_TEST_PRINTELS
S,
SINV,
ENER,
IE,
*OUTPUT,FIELD
*NODE OUTPUT,NSET=QA_TEST_F1
U,RF
*END STEP