*HEADING
 NAFEMS BENCHMARK EXAMPLE 3DNLG-2  :- ELEMENT S4R
 ELASTIC LARGE DEFLECTION RESPONSE OF PEAR-SHAPED CYLINDER UNDER END SHORTENING
*RESTART,WRITE,F=99
*NODE,NSET=ONE
    1,      0.000,      0.000,      0.000
    2,      0.088,      0.000,      0.000
    3,      0.177,      0.000,      0.000
    4,      0.265,      0.000,      0.000
    5,      0.354,      0.000,      0.000
    6,      0.442,      0.000,      0.000
    7,      0.530,      0.000,      0.000
    8,      0.619,      0.000,      0.000
    9,      0.707,      0.000,      0.000
   10,      0.825,      0.007,      0.000
   11,      0.941,      0.028,      0.000
   12,      1.053,      0.062,      0.000
   13,      1.161,      0.109,      0.000
   14,      1.263,      0.169,      0.000
   15,      1.357,      0.240,      0.000
   16,      1.441,      0.321,      0.000
   17,      1.516,      0.412,      0.000
   18,      1.580,      0.511,      0.000
   19,      1.631,      0.617,      0.000
   20,      1.670,      0.729,      0.000
   21,      1.695,      0.844,      0.000
   22,      1.706,      0.961,      0.000
   23,      1.704,      1.079,      0.000
   24,      1.688,      1.195,      0.000
   25,      1.658,      1.309,      0.000
   26,      1.615,      1.419,      0.000
   27,      1.560,      1.523,      0.000
   28,      1.492,      1.619,      0.000
   29,      1.414,      1.707,      0.000
   30,      1.344,      1.778,      0.000
   31,      1.273,      1.849,      0.000
   32,      1.202,      1.919,      0.000
   33,      1.131,      1.990,      0.000
   34,      1.061,      2.061,      0.000
   35,      0.990,      2.131,      0.000
   36,      0.919,      2.202,      0.000
   37,      0.849,      2.273,      0.000
   38,      0.778,      2.343,      0.000
   39,      0.707,      2.414,      0.000
   40,      0.649,      2.467,      0.000
   41,      0.588,      2.516,      0.000
   42,      0.522,      2.560,      0.000
   43,      0.454,      2.598,      0.000
   44,      0.383,      2.631,      0.000
   45,      0.309,      2.658,      0.000
   46,      0.233,      2.680,      0.000
   47,      0.156,      2.695,      0.000
   48,      0.078,      2.704,      0.000
   49,      0.000,      2.707,      0.000
*NCOPY,CHANGENUMBER=1000,OLDSET=ONE,SHIFT,NEWSET=TWO
0.0,0.0,0.40
0.0,0.0,0.0,0.0,0.0,1.0,0.0
*NFILL,NSET=ALL
ONE,TWO,10,100
*ELEMENT,TYPE=S4R
1,1,2,102,101
*ELGEN,ELSET=ALL
1,48,1,1,10,100,100
*SHELL SECTION,ELSET=ALL,MATERIAL=MAT1
0.01, 
*MATERIAL,NAME=MAT1
*ELASTIC
1.0E7,0.3
******************
*NSET,NSET=ENDA,GEN
1,49,1
*NSET,NSET=ENDB,GEN
1001,1049,1
*NSET,NSET=CENTRE,GEN
1,1001,100
49,1049,100
*NSET,NSET=DEPEND,GEN
1002,1049,1
*NSET,NSET=INDEP
1001, 
*NSET,NSET=NODE1
1, 
*BOUNDARY
ENDA,ZSYMM
CENTRE,XSYMM
ENDB,1,2
*EQUATION
2, 
**DEPEND,2,1.0,INDEP,2,-1.0
DEPEND,3,1.0,1001,3,-1.0
*STEP,NLGEOM,INC=200
*STATIC,RIKS
0.2,1.0,,5.0,,1,2,-0.032
*BOUNDARY
**ENDB,3,,-0.0001
INDEP,3,,-0.0001
*NODE PRINT,F=0
*EL PRINT,F=0
** for plotting
*NODE FILE,FRE=1,NSET=INDEP
RF, 
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=INDEP
RF,
*NODE FILE,FRE=1,NSET=NODE1
U, 
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=NODE1
U,
*NODE PRINT,FRE=1,NSET=INDEP
RF, 
*NODE PRINT,FRE=1,NSET=NODE1
U2, 
*END STEP