*HEADING
 Mohr-Coulomb, CPE4R element
 Triaxial Stress (inhomogeneous)
*IMPORT, UPDATE=NO,STATE=YES,step=1
ALLE
*BOUNDARY
 1,PINNED
 2,2
 4,1
*STEP,INC=1
*STATIC,DIRECT
 1.0e-5,1.0e-5
*DLOAD
 1,P2,5000.
 1,P3,5000.
*EL PRINT
 S, 
 SINV, 
 E, 
 PE, 
 EE, 
*NODE PRINT
 U, 
 RF, 
*EL FILE,FREQ=10
 S,E,PE
*output,field
*node output
 u
*element output
 press,mises,peeq,s
*output,history
*element output,elset=alle
 press,mises,peeq,s
*END STEP
*STEP,INC=20
*STATIC,DIRECT
 1.,20.
*BOUNDARY
 3,2,,-.4
 4,2,,-.4
 2,1,,.1
 3,1,,.5
*EL PRINT
 S, 
 SINV, 
 E, 
 PE, 
 EE, 
*NODE PRINT
 U, 
 RF, 
*EL FILE,FREQ=10
 S,E,PE
*output,field
*node output
 u
*element output
 press,mises,peeq,s
*output,history
*element output,elset=alle
 press,mises,peeq,s
*END STEP