*HEADING
 NAFEMS BENCHMARK EXAMPLE 3DNLG-9
 LARGE ELASTIC DEFLECTION OF PINCHED HEMISPHERICAL SHELL [CSS8]
*PREPRINT,MODEL=YES
*NODE
    1, 10.0,    0.0,    0.0
 1001,  9.96,   0.0,    0.0
   17,  0.0,   10.0,    0.0
 1017,  0.0,    9.96,   0.0
  801,  3.0901, 0.0,    9.5105
 1801,  3.0778, 0.0,    9.4725
  817,  0.0,    3.0901, 9.5105
 1817,  0.0,    3.0778, 9.4725
*NGEN,LINE=C,NSET=A
   1,  17,1,,0.0,0.0,0.0
1001,1017,1,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=D
 801, 817,1,,0.0,0.0,9.5105
1801,1817,1,,0.0,0.0,9.4725
*NGEN,LINE=C,NSET=B
   1, 801,50,,0.0,0.0,0.0
1001,1801,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=C
  17, 817,50,,0.0,0.0,0.0
1017,1817,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=E
   2, 802,50,,0.0,0.0,0.0
1002,1802,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=F
   3, 803,50,,0.0,0.0,0.0
1003,1803,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=G
   4,804,50,,0.0,0.0,0.0
1004,1804,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=H
   5, 805,50,,0.0,0.0,0.0
1005,1805,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=I
   6, 806,50,,0.0,0.0,0.0
1006,1806,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=J
   7, 807,50,,0.0,0.0,0.0
1007,1807,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=K
   8, 808,50,,0.0,0.0,0.0
1008,1808,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=L
   9, 809,50,,0.0,0.0,0.0
1009,1809,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=M
  10,810,50,,0.0,0.0,0.0
1010,1810,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=N
  11, 811,50,,0.0,0.0,0.0
1011,1811,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=O
  12, 812,50,,0.0,0.0,0.0
1012,1812,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=P
  13, 813,50,,0.0,0.0,0.0
1013,1813,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=Q
  14, 814,50,,0.0,0.0,0.0
1014,1814,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=R
  15, 815,50,,0.0,0.0,0.0
1015,1815,50,,0.0,0.0,0.0
*NGEN,LINE=C,NSET=S
  16, 816,50,,0.0,0.0,0.0
1016,1816,50,,0.0,0.0,0.0
*NSET,NSET=NALL
A,B,C,E,F,G,H,I,J,
K,L,M,N,O,P,Q,R,S
*NSET,NSET=LOADED
1,17
*ELEMENT,TYPE=CSS8
1,1001,1002,1052,1051,1,2,52,51
*ELGEN,ELSET=ELALL
1,16,1,1,16,50,20
*SOLID SECTION,ELSET=ELALL,MATERIAL=A
*MATERIAL,NAME=A
*ELASTIC
6.825E7,0.3
*BOUNDARY
B,YSYMM
C,XSYMM
1,3
*NSET, NSET=NCORNER
 1, 1001, 17, 1017
*RESTART,WRITE,FREQ=99
*STEP,INC=50,NLGEOM
*STATIC
 0.1,1.0,,0.1
*CLOAD
    1,1,-20.0
 1001,1,-20.0
   17,2,+20.0
 1017,2,+20.0
*EL PRINT,FRE=0
*NODE PRINT,FRE=100,NSET=NCORNER
 U, 
 CF, 
*NODE FILE,FRE=100,NSET=NALL
 U, 
 CF, 
 RF, 
*OUTPUT,FIELD,FREQ=100
*NODE OUTPUT,NSET=NALL
 U,
 CF,
 RF,
*END STEP
*STEP,INC=50,NLGEOM
*STATIC
 0.1,1.0,,0.1
*CLOAD
    1,1,-30.0
 1001,1,-30.0
   17,2,+30.0
 1017,2,+30.0
*END STEP
*STEP,INC=50,NLGEOM
*STATIC
 0.1,1.0,,0.1
*CLOAD
    1,1,-50.0
 1001,1,-50.0
   17,2,50.0
 1017,2,50.0
*END STEP