*HEADING
SLUMP TEST, EXPONENTIAL DRUCKER PRAGER   
(ANG FRIC=5) PSI=8.574;  d=1.9899e-3 (based on mohr data and plane strain)
*NODE
1,0.,0.
6,.1,0.
101,0.,.3
106,.05,.3
*NGEN,NSET=BOT
1,6,1
*NGEN,NSET=TOP
101,106,1
*NFILL,NSET=ALLN
BOT,TOP,10,10
*NSET,NSET=CENT,GENERATE
1,101,10
*ELEMENT,TYPE=CAX4R
1,1,2,12,11
*ELGEN,ELSET=ALLE
1,5,1,1,10,10,10
*ELSET,ELSET=OUTE,GENERATE
1,5,1
*ELSET,ELSET=QA_TEST
OUTE
*SOLID SECTION,ELSET=ALLE,MATERIAL=CONC
*MATERIAL,NAME=CONC
*ELASTIC
2.25,.125
*DRUCKER PRAGER, SHEAR CRITERION=EXPONENT, ECCENTRICITY=0.1
6.63254, 1.0,, 8.574
*DRUCKER PRAGER HARDENING, TYPE=SHEAR
1.989E-3,
*DENSITY
.1,
*BOUNDARY
CENT,1
BOT,2
*AMPLITUDE, NAME=AMP, DEFINITION=SMOOTH STEP
0.0, 0.0, 30.0, 1.0
*STEP, NLGEOM=YES
DIMENSIONLESS SLUMP TEST UNDER GRAVITY LOAD
*DYNAMIC, EXPLICIT
 ,30.0
*DLOAD, amplitude=amp
ALLE,GRAV,0.6666,0,-1.,0.
*OUTPUT,FIELD
*NODE OUTPUT,NSET=ALLN
U, 
*NODE OUTPUT,NSET=BOT
RF, 
*ELEMENT OUTPUT
S,LE,PE
*ELEMENT OUTPUT,ELSET=QA_TEST
PE
*OUTPUT,HISTORY,VARIABLE=PRESELECT
*NODE OUTPUT,NSET=ALLN
U, 
*NODE OUTPUT,NSET=BOT
RF, 
*END STEP