*HEADING
 RATE DEPENDENT FOAM PLASTICITY, UNIAXIAL COMPRESSION TEST, C3D8
 AND MULTIPLE LOAD CASE TEST
 TEMPERATURE DEPENDENT RATE DEPENDENT PARAMETERS
*NODE,NSET=FACE1
 1, 0., 0., 0.
 2, 1., 0., 0.
 3, 1., 1., 0.
 4, 0., 1., 0.
*NODE,NSET=FACE2
 5, 0., 0., 1.
 6, 1., 0., 1.
 7, 1., 1., 1.
 8, 0., 1., 1.
*NSET,NSET=FACE3
 1, 2, 5, 6
*NSET,NSET=FACE4
 2, 3, 6, 7
*NSET,NSET=FACE5
 3, 4, 7, 8
*NSET,NSET=FACE6
 1, 4, 5, 8
*NSET,NSET=ALLN
 FACE1, FACE2
*NSET,NSET=OUT
 7,
*ELEMENT,TYPE=C3D8,ELSET=FOAM
 1, 1, 2, 3, 4, 5, 6, 7, 8
*Distribution Table, name=distTableElastic
 modulus, ratio
*Distribution, name=distElastic, table=distTableElastic, location=element
 , 10.e3, 0.25
 FOAM, 3.E6, 0.2
*SOLID SECTION,ELSET=FOAM,MATERIAL=FOAM
*MATERIAL,NAME=FOAM
*ELASTIC
 distElastic
*CRUSHABLE FOAM
 1.1, 0.1
*CRUSHABLE FOAM HARDENING
 2.200E5,  0.0
 2.465E5,  0.1
 2.729E5,  0.2
 2.990E5,  0.3
 3.245E5,  0.4
 3.493E5,  0.5
 3.733E5,  0.6
 3.962E5,  0.7
 4.180E5,  0.8
 4.387E5,  0.9
 4.583E5,  1.0
 4.938E5,  1.2
 5.248E5,  1.4
 5.515E5,  1.6
 5.743E5,  1.8
 5.936E5,  2.0
 6.294E5,  2.5
 6.520E5,  3.0
 6.833E5,  5.0
 6.883E5, 10.0
*RATE DEPENDENT
  9, 0.9, 10.
 11, 1.1, 20.
**
*INITIAL CONDITIONS,TYPE=TEMPERATURE
 ALLN, 15.
*BOUNDARY
 FACE1, 3
 FACE3, 2
 FACE6, 1
**
*STEP,INC=20,NLGEOM,UNSYMM=YES
*STATIC,DIRECT
 1.E-2, 20.E-2
*BOUNDARY
 FACE2, 3, 3, -0.6
**
*EL PRINT
 S, 
 LE, 
 PE, 
*NODE PRINT
 U, 
 RF, 
*EL FILE
 S, 
 LE, 
 PE, 
*NODE FILE
 U, 
 RF, 
**
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
 S, 
 PE,
 LE,
*OUTPUT,HISTORY,FRE=1
*ELEMENT OUTPUT,ELSET=FOAM
 MISES,
 PRESS,
 PEEQ,
 LE33,
*NODE OUTPUT,NSET=OUT
 U3,
 RF3,
**
*END STEP
**
*STEP,Perturbation
*STATIC
*LOAD CASE,name=compression
*BOUNDARY
 FACE2, 3, 3, -0.6
*END LOAD CASE
*LOAD CASE,name=tension
*BOUNDARY
 FACE2, 3, 3,0.3
*END LOAD CASE
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
 S, 
 PE,
 LE,
*NODE OUTPUT
U,
RF
*END STEP