*HEADING
 NAFEMS BENCHMARK EXAMPLE 3DNLG-7, 16x16 mesh of ELEMENT SC8R
 LARGE DEFLECTION ELASTIC RESPONSE OF A HINGED SPHERICAL SHELL UNDER UNIFORM
 PRESSURE LOADING
*PREPRINT, MODEL=YES
*NODE, NSET=NALL
*ELEMENT, TYPE=SC8R, ELSET=ELALL
*NSET,NSET=CORNER1
*NSET,NSET=CORNER2
*NSET,NSET=CORNER3
*NSET,NSET=CORNER4
*NSET,NSET=CORNERS
CORNER1, CORNER2, CORNER3, CORNER4,
*NSET, NSET=SET1_TOP
*NSET, NSET=SET1_BOT
*NSET, NSET=SET1
SET1_BOT,
SET1_TOP,
CORNER1,
CORNER2,
*NSET, NSET=SET2_TOP
*NSET, NSET=SET2_BOT
*NSET, NSET=SET2
SET2_BOT,
SET2_TOP,
CORNER2,
CORNER3,
*NSET, NSET=SET3_TOP
*NSET, NSET=SET3_BOT
*NSET, NSET=SET3
SET3_BOT,
SET3_TOP,
CORNER3,
CORNER4,
*NSET, NSET=SET4_TOP
*NSET, NSET=SET4_BOT
*NSET, NSET=SET4
SET4_TOP,
SET4_BOT,
CORNER4,
CORNER1,
*NSET,NSET=CENTRE
*NSET, NSET=EDGE
SET1, SET2, SET3, SET4,
*SHELL SECTION,ELSET=ELALL,MAT=MAT, POISSON=ELASTIC
100,
*MATERIAL,NAME=MAT
*ELASTIC
69.0,0.3
*BOUNDARY
CORNERS,1,3
*EQUATION
SET1_BOT,1,1.0,SET1_TOP,1,1.0
SET1_BOT,2,1.0,SET1_TOP,2,1.0
SET1_BOT,3,1.0,SET1_TOP,3,1.0
SET2_BOT,1,1.0,SET2_TOP,1,1.0
SET2_BOT,2,1.0,SET2_TOP,2,1.0
SET2_BOT,3,1.0,SET2_TOP,3,1.0
SET3_BOT,1,1.0,SET3_TOP,1,1.0
SET3_BOT,2,1.0,SET3_TOP,2,1.0
SET3_BOT,3,1.0,SET3_TOP,3,1.0
SET4_BOT,1,1.0,SET4_TOP,1,1.0
SET4_BOT,2,1.0,SET4_TOP,2,1.0
SET4_BOT,3,1.0,SET4_TOP,3,1.0
*RESTART,WRITE,FREQ=999
*STEP,INC=100,NLGEOM
 UNIFORM PRESSURE LOADING - RIKS
*STATIC,RIKS
0.08,1.0,,0.05,1.0
*DLOAD
** put half of the load on top and bottom
ELALL,P1,-0.05
ELALL,P2, 0.05
*EL PRINT,FRE=1
LOADS,
*NODE PRINT,FRE=1,NSET=NALL
U3,
*OUTPUT,FIELD,FREQ=1
*ELEMENT OUTPUT
LOADS,
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=NALL
U,
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT,NSET=CENTRE
U,
*END STEP