*HEADING
PRESSURIZING OF A FLAT RUBBER PLATE (C3D8IH) - OGDEN HYPERELASTICITY  
WITH UNIFORM PRESTRESS (RADIAL AND CIRCUMFERENTIAL)
*NODE
1,0.,0.,0.
11,7.5,0.0,0.0
21,0.0,.125,0.0
31,7.5,.125,0.0
41,0.0,.25,0.0
51,7.5,.25,0.0
*NGEN,LINE=P,NSET=PLANE1
1,11,1,,2.142857143
21,31,1,,2.142857143,.125
41,51,1,,2.142857143,.25,.0
*NCOPY, CHANGE NUMBER=100, OLD SET=PLANE1, NEW SET=PLANE2, SHIFT
0.,0.,0.,
0.,0.,0.,0.,1.,0.,-1.0
*NODE
1000,0,0,0
1001,0,1,0
*MPC
CYCLSYM, plane2,plane1,1000,1001
*NSET,NSET=CENTER
1,21,41
*ELEMENT,TYPE=C3D8IH,ELSET=ONE
1, 1,2,22,21,101,102,122,121
*ELGEN,ELSET=ALL
1, 2,20,1, 10,1,2
*ELSET,ELSET=BOTEL,GENERATE
1,19,2
*BOUNDARY
PLANE1,3
31,1
31,2
CENTER,1
*EQUATION
3, 
51,1,1.0,11,1,1.0,31,1,-2.
3, 
51,2,1.0,11,2,1.0,31,2,-2.
*ELSET,ELSET=ONE
1, 
*SOLID SECTION,MATERIAL=A1,ORIENTATION=CYL,ELSET=ALL
*MATERIAL,NAME=A1
*HYPERELASTIC,OGDEN,N=2
160.,2.,40.,-2.
*ORIENTATION,SYSTEM=CYLINDRICAL,NAME=CYL
0.,0.,0.,0.,1.,0.
*INITIAL CONDITIONS,TYPE=STRESS
ALL, 50.,0.,50.
*RESTART,WRITE,FREQUENCY=2
*STEP,NLGEOM, UNSYMM=YES
 MAINTAIN PRESTRESS - THIS STEP IS FOR VERIFICATION PURPOSES OF INITIAL
 STRESSES IN HYPERELASTICITY
*STATIC
*BOUNDARY
CENTER,XSYMM
31,1,2
*EL PRINT,FREQUENCY=5,ELSET=ALL
COORD , E
DG, 
LOADS, 
S,MISES
ENER, 
*EL PRINT,FREQUENCY=1,ELSET=ONE,SUMMARY=NO
COORD , E
DG, 
LOADS, 
S,MISES
ENER, 
*NODE FILE,NSET=CENTER
U ,RF
*EL FILE,FREQUENCY=2,ELSET=ONE
S, 
LOADS, 
*NODE PRINT
U,RF
*END STEP
*STEP,NLGEOM,INC=12, UNSYMM=YES
*STATIC,RIKS
.01,1.0,.00001,1.0,1.0,1,1,10000.
*DLOAD
BOTEL,P3,200.
*MONITOR,NODE=21,DOF=2
*END STEP