*HEADING
 HYPERFOAM TEST - SIMPLE SHEAR - 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=FACE4
2, 
*NSET,NSET=FACE42
3,6,7
*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 dofs 1 and 3 so that:
**  dof 3: shearing direction; dof 1: transverse direction
** Nom. shear stress (dof 3) = RF3 (@ node 2) / Original area (1 x 1 = 1)
** Nom. transv. stress (dof 1) = RF1 (@ node 2) / Original area (1 x 1 = 1)
** Nominal shear strain may be requested with the NE output identifier.
2, 
FACE42,1,1, 2,1,-1
2, 
FACE42,3,1, 2,3,-1
*ELEMENT,TYPE=C3D8R,ELSET=ONE
 1, 1,2,3,4,5,6,7,8
*SOLID SECTION,ELSET=ONE,MATERIAL=FOAM
*MATERIAL,NAME=FOAM
*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: Simple Shear to 72% strain
*STATIC
.08,.72,.08,.08
*BOUNDARY
FACE6,1,3
FACE4,1,2
FACE4,3,3,.72
*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