*HEADING
 HYPERELASTIC MARLOW 
 TRELOAR'S EXPERIMENTAL DATA
 BIAXIAL TEST DATA ONLY
*RESTART,WRITE,FREQUENCY=5
*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
** Since the S11 output is Cauchy or true stress, we need to
** determine the nominal stress for post-processing.
** Nodes 3,6,7 are tied to node 2 in dof 1 so that:
** Nominal stress (dof 1) = RF1 (@ node 2) / Original area
** (w/c is 1 x 1 = 1)
2, 
FACE42,1,1, 2,1,-1
**
*NSET,NSET=FACE22
5,6,7
*EQUATION
2,
FACE22,3,1, 8,3,-1
**
*ELEMENT,TYPE=C3D8RH,ELSET=ONE
1,1,2,3,4,5,6,7,8
*SOLID SECTION,ELSET=ONE,MATERIAL=TREL
*MATERIAL,NAME=TREL
*HYPERELASTIC,MARLOW
*PLANAR TEST DATA
  0.0,     0.0
  0.6000,  0.0690
  1.6000,  0.1034
  2.4000,  0.1724
  3.3600,  0.2828
  4.2000,  0.4276
  6.0000,  0.8483
  7.8000,  1.3862
  9.6000,  2.0000
 11.1200,  2.4897
 12.9600,  3.0345
 14.8800,  3.4483
 16.5800,  3.7793
 18.2000,  4.0621
*STEP,NLGEOM,INC=50
Step 1: Uniaxial Tension
*STATIC,DIRECT
.25,7.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,7.
*ENERGY PRINT
*EL PRINT,FREQUENCY=5
S, 
E, 
*NODE PRINT,FREQUENCY=5
U,RF
*NODE FILE,FREQUENCY=2
U,RF
*OUTPUT,FIELD,FREQUENCY=0
*NODE OUTPUT,NSET=ALL
U,RF
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=ALL
U,RF
*END STEP
*STEP,NLGEOM,INC=20
Step 2: Unload
*STATIC,DIRECT
.5,7.
*BOUNDARY,OP=MOD
FACE4,1
*END STEP
*STEP,NLGEOM,INC=50
Step 3: Planar Tension (Pure Shear)
*STATIC,DIRECT
.25,4.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE5,2 
FACE6,1
FACE4,1,1,4.
*END STEP
*STEP,NLGEOM,INC=20
Step 4: Unload
*STATIC,DIRECT
.5,4.
*BOUNDARY,OP=MOD
FACE4,1
*END STEP
*STEP,NLGEOM,INC=50
Step 5: Biaxial Tension
*STATIC
.2,4.,,0.25
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,4.
FACE5,2,2,4.
*END STEP
*STEP,NLGEOM,INC=20
Step 6: Unload
*STATIC
.1,4.,,0.5
*BOUNDARY,OP=MOD
FACE4,1
FACE5,2
*END STEP