*HEADING
 Mohr-Coulomb, C3D8 element
 Uniaxial Compression with Static Linear Perturbations
 and Multiple Load Cases
*NODE,NSET=ALLN
 1, 0., 0., 0.
 2, 1., 0., 0.
 3, 1., 1., 0.
 4, 0., 1., 0.
 5, 0., 0., 1.
 6, 1., 0., 1.
 7, 1., 1., 1.
 8, 0., 1., 1.
*ELEMENT,TYPE=C3D8,ELSET=ALLE
 1,1,2,3,4,5,6,7,8
*Distribution Table, name=distTableElastic
 modulus, ratio
*Distribution, name=distElastic, table=distTableElastic, location=element
 , 1.0, 0.0
 ALLE,  300.E3,.3
*SOLID SECTION,ELSET=ALLE,MATERIAL=ALLE
*MATERIAL,NAME=ALLE
*ELASTIC
 distElastic,
*MOHR COULOMB
 40.,40.
*MOHR COULOMB HARDENING
  6.E3, 0.
  9.E3, 0.02
 11.E3, 0.063333
 12.E3, 0.11
 12.E3, 1.
*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
 4,3
 7,3,,-.1
 5,3,,-.1
 6,3,,-.1
 8,3,,-.1
*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
 4,3
 7,3,,-.01
 5,3,,-.01
 6,3,,-.01
 8,3,,-.01
*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
 4,3
 7,3,,-.02
 5,3,,-.02
 6,3,,-.02
 8,3,,-.02
*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
 8,3
 3,3,,-.02
 1,3,,-.02
 2,3,,-.02
 4,3,,-.02
*END LOAD CASE
*EL PRINT
 S, 
 SINV, 
 E, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
 S, 
 SINV, 
 E, 
*NODE OUTPUT
 U, 
 RT, 
*END STEP