*HEADING
SLUMP TEST LINEAR DRUCKER PRAGER                  
(ANG FRIC=20); BETA=30.136; d=1.85E-03 (based on mohr data and 
 plane strain)
 PSI=30.136 (ASSOCIATED FLOW)
*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=CAX4
1,1,2,12,11
*ELGEN,ELSET=ALLE
1,5,1,1,10,10,10
*ELSET,ELSET=OUTE,GENERATE
1,5,1
*SOLID SECTION,ELSET=ALLE,MATERIAL=CONC
*MATERIAL,NAME=CONC
*ELASTIC
2.25,.125
*DRUCKER PRAGER, SHEAR CRITERION=LINEAR
30.164, 1. ,30.164
*DRUCKER PRAGER HARDENING, TYPE=SHEAR
1.844E-03,
*DENSITY
.1,
*BOUNDARY
CENT,1
BOT,2
*RESTART,WRITE,FREQUENCY=100
*STEP,INC=100,NLGEOM
DIMENSIONLESS SLUMP TEST UNDER GRAVITY LOAD
*STATIC
.02,1.0,,.02
*DLOAD
ALLE,GRAV,0.6666,0,-1.,0.
*CONTROLS,ANALYSIS=DISCONTINUOUS
*CONTROLS,PARAMETERS=FIELD
.01,10.,,5.E-5
*NODE PRINT,FREQUENCY=0
*EL PRINT,FREQUENCY=0
*NODE FILE,FREQUENCY=1,NSET=ALLN
U, 
*NODE FILE,FREQUENCY=1,NSET=BOT
RF, 
*EL FILE,FREQUENCY=1,ELSET=OUTE
S,E, SINV, PE
LOADS, 
*MONITOR, NODE=101, DOF=2
*OUTPUT,FIELD, OP=ADD,FREQUENCY=999
*NODE OUTPUT,NSET=ALLN
U, 
*NODE OUTPUT,NSET=BOT
RF, 
*ELEMENT OUTPUT,ELSET=OUTE
S,E, SINV, PE
*OUTPUT,HISTORY,OP=ADD,FREQUENCY=1
*NODE OUTPUT,NSET=ALLN
U, 
*NODE OUTPUT,NSET=BOT
RF, 
*ELEMENT OUTPUT,ELSET=OUTE
S,E, SINV, PE
*END STEP