*HEADING
UNIAXIAL STRETCHING OF A RUBBER SHEET      
W/ CIRCULAR HOLE - MARLOW STRAIN ENERGY FUNCTION USING BIAXIAL TEST DATA INPUT
(TRELOAR'S EXPERIMENTAL TEST DATA)
*NODE
1,0.25,0.
33,0.,0.25
1601,3.25,0.
1617,3.25,3.25
1633,0.,3.25
10000,0.,0.
*NGEN,LINE=C
1,33,1,10000
*NGEN
1601,1617,1
1617,1633,1
*NSET,NSET=DOWN,GENERATE
1,33,1
*NSET,NSET=UP,GENERATE
1601,1633,1
*NFILL, BIAS=0.75, TWO STEP
DOWN,UP,8,200
*NSET,NSET=TOP,GENERATE
1617,1633,2
*NSET,NSET=BOTTOM,GENERATE
1,1601,200
*NSET,NSET=LEFT,GENERATE
33,1633,200
*NSET,NSET=RMPC,GENERATE
1603,1617,2
*NSET,NSET=ONE
1601, 
*ELEMENT,TYPE=CPS8R
1,1,401,405,5,201,403,205,3
*ELGEN,ELSET=ALLE
1,4,400,1,8,4,4
*ELSET,ELSET=TOPE
20,24,28,32
*ELSET,ELSET=LEFTE
29,30,31,32
*BOUNDARY
BOTTOM,2
TOP,2
LEFT,1
*EQUATION
2, 
RMPC,1,1.,1601,1,-1.
*MATERIAL,NAME=OGDEN
*HYPERELASTIC,MARLOW
** TRELOAR'S EXPERIMENTAL DATA
*BIAXIAL TEST DATA,SMOOTH
  0.0   ,  0.0
 13.3470,  0.0200
 22.6147,  0.0600
 34.2592,  0.1100
 37.2930,  0.1400
 47.2776,  0.2000
 62.9770,  0.3100
 73.7183,  0.4200
 93.9066,  0.6800
110.6472,  0.9400
139.1830,  1.4900
179.7103,  2.0300
208.8008,  2.4300
247.4819,  2.7500
285.9668,  3.0700
319.3114,  3.2600
350.6421,  3.4500
*SOLID SECTION,MATERIAL=OGDEN,ELSET=ALLE
0.079, 
*STEP,NLGEOM,INC=100
*STATIC
.05,1.
*BOUNDARY
1601,1,,20.
*EL PRINT,ELSET=TOPE,FREQUENCY=1
S,E
*EL PRINT,ELSET=LEFTE,FREQUENCY=5
S,E
*NODE PRINT,NSET=ONE,FREQUENCY=1
U1,RF1
*NODE PRINT,NSET=LEFT,FREQUENCY=1
U2,RF1
*NODE PRINT,NSET=BOTTOM,FREQUENCY=1
U1,RF2
*NODE FILE,NSET=ONE,FREQUENCY=1
U,RF
*ENERGY FILE
ALLSE
*OUTPUT, HISTORY
*NODE OUTPUT,NSET=ONE
U1,RF1
*END STEP