*HEADING
 NAFEMS BENCHMARK EXAMPLE 3DNLG-4, ELEMENT S4R
 LATERAL TORSIONAL BUCKLING OF AN ELASTIC CANTILEVER SUBJECTED TO TRANSVERSE END LOAD
 CONSERVATIVE LOAD CASE 
*NODE
    1,      0.000,      0.000,      0.000
    2,     10.000,      0.000,      0.000
    3,     20.000,      0.000,      0.000
    4,     30.000,      0.000,      0.000
    5,     40.000,      0.001,      0.000
    6,     50.000,      0.002,      0.000
    7,     60.000,      0.003,      0.000
    8,     70.000,      0.004,      0.000
    9,     80.000,      0.006,      0.000
   10,     90.000,      0.008,      0.000
   11,    100.000,      0.009,      0.000
  101,      0.000,      0.000,      0.833
  102,     10.000,      0.000,      0.833
  103,     20.000,      0.000,      0.833
  104,     30.000,      0.000,      0.833
  105,     40.000,      0.001,      0.833
  106,     50.000,      0.002,      0.833
  107,     60.000,      0.003,      0.833
  108,     70.000,      0.005,      0.833
  109,     80.000,      0.006,      0.833
  110,     90.000,      0.008,      0.833
  111,    100.000,      0.010,      0.833
  201,      0.000,      0.000,      1.667
  202,     10.000,      0.000,      1.667
  203,     20.000,      0.000,      1.667
  204,     30.000,      0.000,      1.667
  205,     40.000,      0.001,      1.667
  206,     50.000,      0.002,      1.667
  207,     60.000,      0.003,      1.667
  208,     70.000,      0.005,      1.667
  209,     80.000,      0.006,      1.667
  210,     90.000,      0.008,      1.667
  211,    100.000,      0.010,      1.667
  301,      0.000,      0.000,      2.500
  302,     10.000,      0.000,      2.500
  303,     20.000,      0.000,      2.500
  304,     30.000,      0.001,      2.500
  305,     40.000,      0.001,      2.500
  306,     50.000,      0.002,      2.500
  307,     60.000,      0.003,      2.500
  308,     70.000,      0.005,      2.500
  309,     80.000,      0.006,      2.500
  310,     90.000,      0.008,      2.500
  311,    100.000,      0.010,      2.500
  401,      0.000,      0.000,      3.333
  402,     10.000,      0.000,      3.333
  403,     20.000,      0.000,      3.333
  404,     30.000,      0.001,      3.333
  405,     40.000,      0.001,      3.333
  406,     50.000,      0.002,      3.333
  407,     60.000,      0.003,      3.333
  408,     70.000,      0.005,      3.333
  409,     80.000,      0.006,      3.333
  410,     90.000,      0.008,      3.333
  411,    100.000,      0.010,      3.333
  501,      0.000,      0.000,      4.167
  502,     10.000,      0.000,      4.167
  503,     20.000,      0.000,      4.167
  504,     30.000,      0.001,      4.167
  505,     40.000,      0.001,      4.167
  506,     50.000,      0.002,      4.167
  507,     60.000,      0.003,      4.167
  508,     70.000,      0.005,      4.167
  509,     80.000,      0.006,      4.167
  510,     90.000,      0.008,      4.167
  511,    100.000,      0.010,      4.167
  601,      0.000,      0.000,      5.000
  602,     10.000,      0.000,      5.000
  603,     20.000,      0.000,      5.000
  604,     30.000,      0.001,      5.000
  605,     40.000,      0.001,      5.000
  606,     50.000,      0.002,      5.000
  607,     60.000,      0.004,      5.000
  608,     70.000,      0.005,      5.000
  609,     80.000,      0.007,      5.000
  610,     90.000,      0.008,      5.000
  611,    100.000,      0.010,      5.000
*NSET,NSET=FIX,GEN
1,601,100
*NSET,NSET=N311
311, 
*ELEMENT,TYPE=S4R
1,1,2,102,101
*ELGEN,ELSET=SHELL
1,10,1,1,6,100,100
*SHELL SECTION,ELSET=SHELL,MATERIAL=MAT
0.2, 
*MATERIAL,NAME=MAT
*ELASTIC
1.0E4,0.0
*BOUNDARY
FIX,1,6
*STEP,NLGEOM,INC=50
*STATIC,RIKS
0.0005,1.0,1E-4,0.001,0.02,311,2,0.065
*CLOAD
311,3,-1.0
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=1,NSET=N311
U2,CF3
*NODE FILE,FREQ=1,NSET=N311
U, 
CF, 
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=N311
U,
CF,
*END STEP