*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
*nset, nset=ngrp1
 11,13,14
*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, test data input
**1.75,0.3,0.6
*uniaxial test data
1.104,0.5
1.005,0.475
0.914,0.45
0.8291,0.425
0.6775,0.375
0.610,0.35
0.5474,0.325
0.4892,0.3
0.4349,0.275
0.3843,0.25
0.3367,0.225
0.2919,0.2
0.2496,0.175
0.2094,0.15
0.1711,0.125
0.1344,0.1
0.099174,0.075
0.065122,0.05
0.032106,0.025
0.,0.
*biaxial test data
1.162,0.5
1.042,0.475
0.9319,0.45
0.8323,0.425
0.7424,0.4
0.6617,0.375
0.5891,0.35
0.5238,0.325
0.4649,0.3
0.4114,0.275
0.3626,0.25
0.3177,0.225
0.2760,0.2
0.2369,0.175
0.1999,0.15
0.1646,0.125
0.1305,0.1
0.097381,0.075
0.064801,0.05
0.032453,0.025
0.,0.
*planar test data
2.278,1.0
2.007,0.95
1.763,0.9
1.546,0.85
1.355,0.8
1.188,0.75
1.043,0.7
0.9162,0.65
0.8046,0.6
0.7055,0.55
0.6167,0.5
0.5360,0.45
0.4619,0.4
0.3931,0.35
0.3285,0.3
0.2675,0.25
0.2094,0.2
0.1539,0.15
0.1006,0.1
0.049390,0.05
0.,0.
*equation
 2,
 face14,1,1,ndisp1,1,-1
*equation
 2,
 ngrp1,3,1,ndisp1,3,-1
*nset, nset=ndisp
 ndisp1
*boundary
 face12,3
 ndisp1,3
 face13,2
 face16,1
*step 
 planar loading
*dynamic,explicit
 ,1.0e-2
*boundary, type=velocity
 ndisp,1,1,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
*end step
*step 
 planar loading
*dynamic,explicit
 ,1.0e-2
*boundary, type=velocity
 ndisp,1,1,1.0e2
*end step
*step 
 unloading
*dynamic,explicit
 ,1.0e-2
*boundary, type=velocity
 ndisp,1,1,-1.0e2
*end step