*HEADING
PRESSURIZING OF A FLAT RUBBER PLATE (S4R)                 
*NODE
 1, 2.652,2.652
 5, 5.3033, 5.3033
 4401, 0.0, 0.0
  401, 3.75, 0.0
 4001, 0.0 , 3.75
 4005, 0.0, 7.5
 405, 7.5,0.0
*NGEN,LINE=C,NSET=C1
 5,405,100
*NGEN,LINE=C,NSET=C2
 5,4005,1000
*NGEN,NSET=M1
 1,401,100
*NGEN,NSET=M2
 1,4001,1000
*NFILL
 M1,C1,4,1
 M2,C2,4,1
*NGEN
 4001,4401,100
 101, 4101,1000
 201, 4201,1000
 301, 4301,1000
 401, 4401,1000
*ELEMENT,TYPE=S4R
 1,1,101,102,2
*ELGEN,ELSET=EALL
 1,4,100,1,4,1,4
*ELEMENT,TYPE=S4R
 101,1,2,1002,1001
*ELGEN,ELSET=EALL
 101,4,1,1,4,1000,4   
*ELEMENT,TYPE=S4R
 1001, 1001, 1101, 101,1
*ELGEN,ELSET=EALL
 1001, 4, 100,1,4,1000,4
*SHELL SECTION,MATERIAL=M1,ELSET=EALL,CONTROL=HGC
 0.5,
*SECTION CONTROLS,NAME=HGC,HOURGLASS=STIFFNESS
*MATERIAL,NAME=M1
*HYPERELASTIC,OGDEN,N=2
160.,2.,40.,-2.
*NSET,NSET=XAXIS,GENERATE
 401,405
 401,4401,1000
*NSET,NSET=YAXIS,GENERATE
 4001,4005
 4001,4401,100
*NSET,NSET=HINGE
 C1,C2
*BOUNDARY
 XAXIS,YSYMM
 YAXIS,XSYMM
 HINGE,1,3
*RESTART,WRITE,FREQUENCY=2
*ELSET,ELSET=E1016
 1016,
*NSET,NSET=N4401
 4401,
*STEP,NLGEOM,INC=18
*STATIC,RIKS
.01,1.0,.00001,1.0,1.0,4401,1,10000.
*DLOAD
EALL,P,2.E+2
*MONITOR,NODE=4401,DOF=3
*EL PRINT,FREQUENCY=2
 S11,S22,S12
 E11,E22,E12,
 STH,
 DG,
*EL FILE,FREQUENCY=2,ELSET=E1016
LOADS, 
S , 
*OUTPUT,FIELD
*ELEMENT OUTPUT
LOADS,
ENER,
*OUTPUT,FIELD,FREQ=2
*ELEMENT OUTPUT,ELSET=E1016
LOADS,
S,
*EL FILE,ELSET=E1016
STH, 
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=E1016
STH,
*NODE FILE,NSET=N4401
U , 
*OUTPUT,FIELD
*NODE OUTPUT,NSET=N4401
U,
*END STEP