*HEADING
 Mohr-Coulomb, C3D8 element
 Uniaxial Compression with Static Linear Perturbations
 and Multiple Load Cases
*NODE,NSET=ALLN
*ELEMENT,TYPE=C3D8,ELSET=ALLE
*SOLID SECTION,ELSET=ALLE,MATERIAL=ALLE
*MATERIAL,NAME=ALLE
*ELASTIC
 300.E3,.3
*MOHR COULOMB
*MOHR COULOMB HARDENING
*STEP,INC=10
 Step 1 - Nonlinear
*STATIC,DIRECT
 1.,10.
*BOUNDARY
 1,PINNED
 2,2
 5,2
 6,2
 4,1
 5,1
 8,1
 2,3
 3,3
*EL PRINT
 S,
 SINV,
 E,
 PE,
 EE,
*NODE PRINT
 U,
 RF,
*END STEP
*STEP,PERTURBATION,INC=1
 Step 2 - Linear Perturbation
*STATIC,DIRECT
 .1,.1
*BOUNDARY
 1,PINNED
 2,2
 5,2
 6,2
 4,1
 5,1
 8,1
 2,3
 3,3
*EL PRINT
 S,
 SINV,
 E,
*NODE PRINT
 U,
 RF,
*END STEP
*STEP,INC=10
 Step 3 - Rest of Nonlinear
*STATIC,DIRECT
 1.,10.
*BOUNDARY
 7,3,,-.2
 5,3,,-.2
 6,3,,-.2
 8,3,,-.2
*END STEP
*STEP,PERTURBATION,INC=1
 Step 4 - MULTIPLE LOAD CASES
*STATIC,DIRECT
 .1,.1
*LOAD CASE,NAME="UNIAXIAL COMPRESSION FACE2"
*BOUNDARY,OP=NEW
 1,PINNED
 2,2
 5,2
 6,2
 4,1
 5,1
 8,1
 2,3
 3,3
*END LOAD CASE
*LOAD CASE,NAME= "UNIAXIAL COMPRESSION FACE1"
*BOUNDARY,OP=NEW
 5,PINNED
 6,2
 1,2
 2,2
 8,1
 1,1
 4,1
 6,3
 7,3
*END LOAD CASE
*EL PRINT
 S,
 SINV,
 E,
*NODE PRINT
 U,
 RF,
*OUTPUT,FIELD
*ELEMENT OUTPUT
 S,
 SINV,
 E,
*NODE OUTPUT
 U,
 RF,
*END STEP