*HEADING
*PREPRINT, MODEL=YES
*NODE,NSET=ALL
 11,
 12,1.
 13,1.,1.,
 14,0.,1.,
 15,0.,0.,1.
 16,1.,0.,1.
 17,1.,1.,1.
 18,0.,1.,1.
*NSET,NSET=FACE11
 11,12,13,14
*NSET,NSET=FACE12
 15,16,17,18
*NSET,NSET=FACE13
 11,12,15,16
*NSET,NSET=FACE14
 13,16,17
*NSET,NSET=FACE15
 13,14,17,18
*NSET,NSET=FACE16
 14,11,18,15
*NSET, NSET=NDISP1
 12
*ELEMENT,TYPE=C3D8R,ELSET=C3D
 1,11,12,13,14,15,16,17,18
**
*ELSET, ELSET=EALL
 C3D
*SOLID SECTION,ELSET=EALL,MATERIAL=low_density_foam
 1.0,
*MATERIAL,NAME=low_density_foam
*DENSITY
 5.0000e-010,
*HYPERFOAM, N=3
-1048.43,0.3025,532.20,0.3958,517.027,0.2135,0.2,0.2,
0.2
*MULLINS EFFECT
 1.75,  0.3,  0.6
*EQUATION
 2,
 FACE14,1,1,NDISP1,1,-1
*NSET, NSET=NDISP
 NDISP1
*BOUNDARY
 FACE11,3
 FACE13,2
 FACE16,1
*STEP,NLGEOM
 UNIAXIAL LOADING
*STATIC
5.0e-4,1.0e-2,,5.0e-4
*BOUNDARY
 NDISP,1,1,0.5
*OUTPUT,FIELD, VARIABLE=PRESELECT
*OUTPUT, HISTORY, VARIABLE=PRESELECT
*ELEMENT OUTPUT, ELSET=EALL
 S11, LE11, NE11, ENER
*NODE OUT, NSET=NDISP
 RF1
*ENERGY OUTPUT
 ALLDMD
*EL PRINT
DGP,
SENER,
*NODE PRINT,NSET=NDISP
U1,RF1
*END STEP
*STEP,NLGEOM
 UNLOADING
*STATIC
5.0e-4,1.0e-2,,5.0e-4
*BOUNDARY,OP=NEW
 FACE11,3
 FACE13,2
 FACE16,1
 NDISP,1,1,0.
*END STEP
*STEP,NLGEOM
 UNIAXIAL LOADING
*STATIC
5.0e-4,1.0e-2,,5.0e-4
*BOUNDARY
 NDISP,1,1,1.0
*END STEP
*STEP,NLGEOM
 UNLOADING
*STATIC
5.0e-4,1.0e-2,,5.0e-4
*BOUNDARY,OP=NEW
 FACE11,3
 FACE13,2
 FACE16,1
 NDISP,1,1,0.
*END STEP