*HEADING
UNIAXIAL STRETCHING OF A RUBBER SHEET W/ A CIRC. HOLE - VISCOELASTIC
MARLOW STRAIN ENERGY FUNCTION USING UNIAXIAL TEST DATA INPUT 
VISCOELASTIC RELAXATION AFTER LOADING TO 1/2 STRESS LEVEL
*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=BIDERMAN
*HYPERELASTIC,MARLOW,MODULI=INSTANTANEOUS
*UNIAXIAL TEST DATA,SMOOTH
  0.0   ,  0.0
 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
*VISCOELASTIC,TIME=PRONY
0.25,,5.
0.25,,10.
*SOLID SECTION,MATERIAL=BIDERMAN,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
*EL FILE,FREQUENCY=99
S,E
*ENERGY FILE
ALLSE
*OUTPUT, HISTORY
*NODE OUTPUT,NSET=ONE
U1,RF1
*END STEP
*STEP,NLGEOM,INC=120
*VISCO,CETOL=0.1
1.,100.
*BOUNDARY
1601,1,,20.
*EL PRINT,ELSET=TOPE,FREQUENCY=10
S,E
*EL PRINT,ELSET=LEFTE,FREQUENCY=10
S,E
*NODE PRINT,NSET=ONE,FREQUENCY=10
U1,RF1
*NODE PRINT,NSET=LEFT,FREQUENCY=10
U2,RF1
*NODE PRINT,NSET=BOTTOM,FREQUENCY=10
U1,RF2
*NODE FILE,NSET=ONE,FREQUENCY=10
U,RF
*EL FILE,FREQUENCY=20
S,E,CE
*ENERGY FILE,FREQUENCY=10
ALLSE
*END STEP