*HEADING
 Mohr-Coulomb, C3D8 element
 Uniaxial Compression with temperature dependence 
*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
 30.,30.,10.
 50.,50.,30.
*MOHR COULOMB HARDENING
  6.E3, 0.
  9.E3, 0.02
 11.E3, 0.063333
 12.E3, 0.11
 12.E3, 1.
*BOUNDARY
 1,PINNED
 2,2
 5,2
 6,2
 4,1
 5,1
 8,1
 2,3
 3,3
 4,3
*INITIAL CONDITIONS,TYPE=TEMPERATURE
 ALLN,20.
*STEP,INC=20
*STATIC,DIRECT
 1.,20.
*BOUNDARY
 7,3,,-.2
 5,3,,-.2
 6,3,,-.2
 8,3,,-.2
*EL PRINT
 S, 
 SINV, 
 E, 
 PE, 
 EE, 
*NODE PRINT
 U, 
 RF, 
*EL FILE,FREQ=10
 S,E,PE
*END STEP