*HEADING
 HYPERFOAM TEST - UNIAXIAL COMPRESSION - UNI + SS DATA - C3D8R - N=2 
*NODE,NSET=ALL
1,
2,1.
3,1.,1.,
4,0.,1.,
5,0.,0.,1.
6,1.,0.,1.
7,1.,1.,1.
8,0.,1.,1.
*NSET,NSET=FACE1
1,2,3,4
*NSET,NSET=FACE2
5,6,7,8
*NSET,NSET=FACE3                                                       
1,2,5,6
*NSET,NSET=FACE4
2, 
*NSET,NSET=FACE42
3,6,7
*NSET,NSET=FACE5
3,4,7,8
*NSET,NSET=FACE6
4,1,8,5
*EQUATION
** The stress output S is Cauchy or true stress and strain output E is 
** accumulated strain. For postprocessing purposes (to compare with test data), 
** nominal stress and strain values are needed and we can do the following:
** Nodes 3,6,7 (NSET=FACE42) are tied to node 2 in dof 1 so that:
** Nominal stress (dof 1) = RF1 (@ node 2) / Original area (1 x 1 = 1)
** Nominal strain can be requested with the NE output identifier.
2, 
FACE42,1,1, 2,1,-1
*ELEMENT,TYPE=C3D8R,ELSET=ONE
 1, 1,2,3,4,5,6,7,8
*SOLID SECTION,ELSET=ONE,MATERIAL=FOAM
*MATERIAL,NAME=FOAM
** To change the strain energy function order, change the N parameter.
*HYPERFOAM,N=2,TEST DATA INPUT,POISSON=0.
*UNIAXIAL TEST DATA
** nominal stress, nominal strain
 -.0217,   -.05
 -.0317,   -.10
 -.0367,   -.15
 -.0402,   -.20                              
 -.0433,   -.25
 -.0467,   -.30
 -.0504,   -.35
 -.0542,   -.40
 -.0604,   -.45
 -.0668,   -.50
 -.0759,   -.55
 -.0909,   -.60
 -.1083,   -.65
 -.1410,   -.70
 -.1933,   -.75
 -.2896,   -.80
** To fit uniaxial test data only, remove the *SIMPLE SHEAR TEST DATA card
** and the stress-strain data below.
*SIMPLE SHEAR TEST DATA      
** nom. shear stress, nom. shear strain, nom. transverse stress
  .0140, .08, .0046
  .0334, .16, .0166
  .0533, .24, .0366
  .0853, .32, .0573
  .1280, .40, .0817
  .1653, .48, .1098
  .2080, .56, .1394
  .2560, .64, .1666
  .2987, .72, .1904
*RESTART,WRITE,FREQUENCY=5
*STEP,NLGEOM,INC=20
Step 1: Uniaxial Compression to 80% strain
*STATIC
.05,.80,.05,.05
*BOUNDARY
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,-.80
*ENERGY PRINT
*EL PRINT,SUMMARY=NO,FREQUENCY=5
S, 
E, 
LE, 
NE, 
DG, 
DGP, 
*EL FILE,FREQUENCY=5
S,E,LE,NE,DG,DGP
EP,LEP,NEP
*NODE PRINT
U,RF
*NODE FILE,FREQUENCY=2
U,RF
**
*OUTPUT,FIELD,FREQUENCY=5
*NODE OUTPUT
U, 
RF, 
*ELEMENT OUTPUT
S, 
E, 
LE, 
NE, 
*OUTPUT,HISTORY,FREQUENCY=1
*ELEMENT OUTPUT,ELSET=ONE
S,E,NE
*NODE OUTPUT,NSET=ALL
U, 
RF, 
*END STEP