*HEADING
 : NAFEMS BENCHMARK EXAMPLE 3DNLG-10
 ELASTO-PLASTIC BEHAVIOUR OF A STIFFENED CYLINDRICAL PANEL UNDER COMPRESSIVE END LOAD
*NODE,NSET=PANEL,SYSTEM=C
    1,400.000000, -9.000000,  0.000000
    2,400.000000, -4.500000,  0.000000
    3,400.000000,  0.000000,  0.000000
    4,400.000000,  4.500000,  0.000000
    5,400.000000,  9.000000,  0.000000
   11,399.929993, -9.000000, 10.000000
   13,400.000000,  0.000000, 10.000000
   15,400.070007,  9.000000, 10.000000
   21,399.861023, -9.000000, 20.000000
   22,399.929993, -4.500000, 20.000000
   23,400.000000,  0.000000, 20.000000
   24,400.070007,  4.500000, 20.000000
   25,400.138977,  9.000000, 20.000000
   31,399.794189, -9.000000, 30.000000
   33,400.000000,  0.000000, 30.000000
   35,400.205811,  9.000000, 30.000000
   41,399.730469, -9.000000, 40.000000
   42,399.861023, -4.500000, 40.000000
   43,400.000000,  0.000000, 40.000000
   44,400.138977,  4.500000, 40.000000
   45,400.269531,  9.000000, 40.000000
   51,399.670837, -9.000000, 50.000000
   53,400.000000,  0.000000, 50.000000
   55,400.329163,  9.000000, 50.000000
   61,399.616241, -9.000000, 60.000000
   62,399.794189, -4.500000, 60.000000
   63,400.000000,  0.000000, 60.000000
   64,400.205811,  4.500000, 60.000000
   65,400.383759,  9.000000, 60.000000
   71,399.567444, -9.000000, 70.000000
   73,400.000000,  0.000000, 70.000000
   75,400.432556,  9.000000, 70.000000
   81,399.525238, -9.000000, 80.000000
   82,399.730469, -4.500000, 80.000000
   83,400.000000,  0.000000, 80.000000
   84,400.269531,  4.500000, 80.000000
   85,400.474762,  9.000000, 80.000000
   91,399.490265, -9.000000, 90.000000
   93,400.000000,  0.000000, 90.000000
   95,400.509735,  9.000000, 90.000000
  101,399.463013, -9.000000,100.000000
  102,399.670837, -4.500000,100.000000
  103,400.000000,  0.000000,100.000000
  104,400.329163,  4.500000,100.000000
  105,400.536987,  9.000000,100.000000
  111,399.443939, -9.000000,110.000000
  113,400.000000,  0.000000,110.000000
  115,400.556061,  9.000000,110.000000
  121,399.433319, -9.000000,120.000000
  122,399.616241, -4.500000,120.000000
  123,400.000000,  0.000000,120.000000
  124,400.383759,  4.500000,120.000000
  125,400.566681,  9.000000,120.000000
  131,399.431305, -9.000000,130.000000
  133,400.000000,  0.000000,130.000000
  135,400.568695,  9.000000,130.000000
  141,399.437958, -9.000000,140.000000
  142,399.567444, -4.500000,140.000000
  143,400.000000,  0.000000,140.000000
  144,400.432556,  4.500000,140.000000
  145,400.562042,  9.000000,140.000000
  151,399.453125, -9.000000,150.000000
  153,400.000000,  0.000000,150.000000
  155,400.546875,  9.000000,150.000000
  161,399.476624, -9.000000,160.000000
  162,399.525238, -4.500000,160.000000
  163,400.000000,  0.000000,160.000000
  164,400.474762,  4.500000,160.000000
  165,400.523376,  9.000000,160.000000
  171,399.508057, -9.000000,170.000000
  173,400.000000,  0.000000,170.000000
  175,400.491943,  9.000000,170.000000
  181,399.546967, -9.000000,180.000000
  182,399.490265, -4.500000,180.000000
  183,400.000000,  0.000000,180.000000
  184,400.509735,  4.500000,180.000000
  185,400.453033,  9.000000,180.000000
  191,399.592804, -9.000000,190.000000
  193,400.000000,  0.000000,190.000000
  195,400.407196,  9.000000,190.000000
  201,399.644806, -9.000000,200.000000
  202,399.463013, -4.500000,200.000000
  203,400.000000,  0.000000,200.000000
  204,400.536987,  4.500000,200.000000
  205,400.355194,  9.000000,200.000000
**
** Transform panel nodes to a cylindrical system.
**
*TRANSFORM,NSET=PANEL,TYPE=C
0,0,0, 0,0,1
**
** Define a new co-ordinate system to use in the definition
** of the stiffener.
**
*SYSTEM
400,0.0,0.0,400,0.0,200
200,0.0,0.0
**
** Define stiffener nodes in rectangular co-ordinates.
**
*NODE,NSET=STIFF,SYSTEM=R
 1001,  0.000000,  2.500000,  0.000000
 1002,  0.000000,  5.000000,  0.000000
 1003,  0.000000,  7.500000,  0.000000
 1004,  0.000000, 10.000000,  0.000000
 1012, 10.000000,  5.000000, -0.016644
 1014, 10.000000, 10.000000, -0.023538
 1021, 20.000000,  2.500000, -0.017959
 1022, 20.000000,  5.000000, -0.033185
 1023, 20.000000,  7.500000, -0.043358
 1024, 20.000000, 10.000000, -0.046930
 1032, 30.000000,  5.000000, -0.049521
 1034, 30.000000, 10.000000, -0.070034
 1041, 40.000000,  2.500000, -0.035477
 1042, 40.000000,  5.000000, -0.065552
 1043, 40.000000,  7.500000, -0.085648
 1044, 40.000000, 10.000000, -0.092705
 1052, 50.000000,  5.000000, -0.081179
 1054, 50.000000, 10.000000, -0.114805
 1061, 60.000000,  2.500000, -0.052120
 1062, 60.000000,  5.000000, -0.096306
 1063, 60.000000,  7.500000, -0.125830
 1064, 60.000000, 10.000000, -0.136197
 1072, 70.000000,  5.000000, -0.110839
 1074, 70.000000, 10.000000, -0.156749
 1081, 80.000000,  2.500000, -0.067481
 1082, 80.000000,  5.000000, -0.124688
 1083, 80.000000,  7.500000, -0.162913
 1084, 80.000000, 10.000000, -0.176335
 1092, 90.000000,  5.000000, -0.137769
 1094, 90.000000, 10.000000, -0.194834
 1101,100.000000,  2.500000, -0.081179
 1102,100.000000,  5.000000, -0.150000
 1103,100.000000,  7.500000, -0.195984
 1104,100.000000, 10.000000, -0.212132
 1112,110.000000,  5.000000, -0.161306
 1114,110.000000, 10.000000, -0.228122
 1121,120.000000,  2.500000, -0.092879
 1122,120.000000,  5.000000, -0.171618
 1123,120.000000,  7.500000, -0.224230
 1124,120.000000, 10.000000, -0.242705
 1132,130.000000,  5.000000, -0.180872
 1134,130.000000, 10.000000, -0.255792
 1141,140.000000,  2.500000, -0.102292
 1142,140.000000,  5.000000, -0.189011
 1143,140.000000,  7.500000, -0.246955
 1144,140.000000, 10.000000, -0.267302
 1152,150.000000,  5.000000, -0.195984
 1154,150.000000, 10.000000, -0.277164
 1161,160.000000,  2.500000, -0.109186
 1162,160.000000,  5.000000, -0.201749
 1163,160.000000,  7.500000, -0.263598
 1164,160.000000, 10.000000, -0.285317
 1172,170.000000,  5.000000, -0.206271
 1174,170.000000, 10.000000, -0.291711
 1181,180.000000,  2.500000, -0.113391
 1182,180.000000,  5.000000, -0.209520
 1183,180.000000,  7.500000, -0.273751
 1184,180.000000, 10.000000, -0.296306
 1192,190.000000,  5.000000, -0.211478
 1194,190.000000, 10.000000, -0.299075
 1201,200.000000,  2.500000, -0.114805
 1202,200.000000,  5.000000, -0.212132
 1203,200.000000,  7.500000, -0.277164
 1204,200.000000, 10.000000, -0.300000
**
** Transform the stiffener nodes.
**
*TRANSFORM,NSET=STIFF,TYPE=R
0.0,0.0,100,-100,0.0,0,0
*ELEMENT,TYPE=S8R5,ELSET=PANEL
1,1,21,23,3,11,22,13,2
*ELGEN,ELSET=PANEL
1,2,2,1,10,20,10,1
*ELEMENT,TYPE=S8R5,ELSET=STIFF
101,3,23,1022,1002,13,1021,1012,1001
111,1002,1022,1024,1004,1012,1023,1014,1003
*ELGEN,ELSET=STIFF
101,10,20
111,10,20
*NSET,NSET=EDGEA,GEN
1,201,10
5,205,10
*NSET,NSET=ENDAP
1,2,3,4,5
*NSET,NSET=ENDAS
1001,1002,1003,1004
*NSET,NSET=SYMMP
201,202,203,204,205
*NSET,NSET=SYMMS
1201,1202,1203,1204
*NSET,NSET=NODE3
3, 
*SHELL SECTION, ELSET=PANEL, MATERIAL=MAT1
1.0,9
*SHELL SECTION, ELSET=STIFF, MATERIAL=MAT1
1.0,9
*MATERIAL, NAME=MAT1
*ELASTIC, TYPE=ISO
210000., 0.3
*PLASTIC
350,0.0
*BOUNDARY
EDGEA,2,2,0
EDGEA,4,4,0
EDGEA,6,6,0
ENDAP,1,2,0
ENDAP,4,4,0
ENDAP,6,6,0
ENDAS,2,4,0
SYMMP,3,5,0
SYMMS,1,1,0
SYMMS,5,6,0
*ELSET, ELSET=LOAD, GENERATE
 1,2,1
 101,111,10
*SURFACE, NAME=SURF_LOAD
 LOAD, E4
**
**
** Initially this job is linear - until plasticity
** occurs. 
** In step 1 apply load using a normal static 
** step. 
** In step to apply the remainder of the load
** using the RIKS algorithm
**
**
**
*RESTART,WRITE,F=10
*STEP, AMPLITUDE=RAMP, INC=10000, NLGEOM
*STATIC,RIKS
0.05, 1.,,0.1,,3,3,0.2
*MONITOR,NODE=3,DOF=3
*CONTROLS,PARAM=FIELD,FIELD=DISPLACEMENT
.01,1.
*CONTROLS,PARAM=FIELD,FIELD=ROTATION
.01,1.
*DSLOAD, CONSTANT RESULTANT=YES, FOLLOWER=NO
 SURF_LOAD, EDLD, 191., 0.,0.,1.
*OUTPUT, HISTORY
*NODE OUTPUT, NSET=NODE3
 U3,
*NODE PRINT,F=0
*EL PRINT,F=0
*NODE FILE,F=1,NSET=NODE3
U, 
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=NODE3
U,
*NODE FILE,F=1,NSET=ENDAP 
CF, 
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=ENDAP
CF,
*NODE FILE,F=1,NSET=ENDAS
CF, 
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=ENDAS
CF,
*ENDSTEP