*HEADING
UNIAXIAL STRETCHING OF A RUBBER SHEET W/ CIRCULAR HOLE - OGDEN - CPS4 
STRAIN ENERGY FUNCTION USING TEST DATA INPUT (TRELOAR'S EXPERIMENTAL TEST DATA)
*RESTART,WRITE,FREQUENCY=10
*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
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=CPS4
1,1,201,203,3
*ELGEN,ELSET=ALLE
1,8,200,1,16,2,8
*ELSET,ELSET=TOPE,GENERATE
72,128,8
*ELSET,ELSET=LEFTE
121,128,8
*BOUNDARY
BOTTOM,2
TOP,2
LEFT,1
*EQUATION
2, 
RMPC,1,1.,1601,1,-1.
*MATERIAL,NAME=OGDEN
*HYPERELASTIC,OGDEN,N=3,TEST DATA INPUT
** TRELOAR'S EXPERIMENTAL DATA
*UNIAXIAL TEST DATA
 22.0543,  0.1338
 34.6574,  0.2675
 44.1101,  0.3567
 59.8636,  0.6242
 75.6171,  0.8917
 85.0684,  1.1592
 97.6714,  1.4268
126.0280,  2.0510
151.2328,  2.5860
176.4390,  3.0318
230.0003,  3.7898
283.5617,  4.3694
333.9726,  4.8153
390.6858,  5.1720
441.0967,  5.4395
491.5077,  5.7070
545.0690,  5.9299
598.6304,  6.0637
649.0413,  6.1975
702.6027,  6.3312
756.1654,  6.4650
809.7267,  6.5541
913.6990,  6.6433
*BIAXIAL TEST DATA
 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
*PLANAR TEST DATA
  8.5339,  0.0690
 22.7570,  0.1034
 34.1354,  0.1724
 47.7896,  0.2828
 59.7370,  0.4276
 85.3386,  0.8483
110.9402,  1.3862
136.5418,  2.0000
158.1609,  2.4897
184.3314,  3.0345
211.6397,  3.4483
235.8190,  3.7793
258.8604,  4.0621
*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
*END STEP