*heading
*preprint, model=yes
*node,nset=all
*nset,nset=face11
*nset,nset=face12
*nset,nset=face13
*nset,nset=face14
*nset,nset=face15
*nset,nset=face16
*nset, nset=ndisp1
*nset, nset=ngrp1
*element,type=c3d8r,elset=c3d
**
*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, test data input
**1.75, 0.3, 0.6
*biaxial test data
*biaxial test data
*equation
 face14,1,1,ndisp1,1,-1
*equation
 ngrp1,3,1,ndisp1,3,-1
*nset, nset=ndisp
 ndisp1
*boundary,type=velocity
 face12,3
 face13,2
 face16,1
*step
 biaxial loading
*dynamic,explicit
 ,1.0e-2
*boundary, type=velocity
 ndisp,1,1,0.5e2
 ndisp,3,3,-0.5e2
*output,field
*element output
 s, le
*node out
 rf, u
*output, history, variable=preselect
*element output, elset=eall
 s11, le11, ne11, ener
*node out, nset=ndisp
 rf1,rf3
*end step
*step
 unloading
*dynamic,explicit
 ,1.0e-2
*boundary, type=velocity
 ndisp,1,1,-0.5e2
 ndisp,3,3,0.5e2
*end step
*step
 biaxial loading
*dynamic,explicit
 ,1.0e-2
*boundary, type=velocity
 ndisp,1,1,1.0e2
 ndisp,3,3,-1.0e2
*end step
*step
 unloading
*dynamic,explicit
 ,1.0e-2
*boundary,type=velocity
 ndisp,1,1,-1.0e2
 ndisp,3,3,1.0e2
*end step