*HEADING
 : NAFEMS BENCHMARK EXAMPLE 3DNLG-10
 ELASTO-PLASTIC BEHAVIOUR OF A STIFFENED CYLINDRICAL PANEL UNDER COMPRESSIVE END LOAD
 TEST 9-DIGIT LABELS
*PREPRINT,ECHO=YES,HIST=YES,MODEL=YES
*NODE,NSET=PANEL,SYSTEM=C
900000001, 400.000000, -9.000000, 0.000000
900000002, 400.000000, -6.750000, 0.000000
900000003, 400.000000, -4.500000, 0.000000
900000004, 400.000000, -2.250000, 0.000000
900000005, 400.000000, 0.000000, 0.000000
900000006, 400.000000, 2.250000, 0.000000
900000007, 400.000000, 4.500000, 0.000000
900000008, 400.000000, 6.750000, 0.000000
900000009, 400.000000, 9.000000, 0.000000
900000011, 399.929993, -9.000000, 10.000000
900000012, 399.947418, -6.750000, 10.000000
900000013, 399.964935, -4.500000, 10.000000
900000014, 399.982452, -2.250000, 10.000000
900000015, 400.000000, 0.000000, 10.000000
900000016, 400.017548, 2.250000, 10.000000
900000017, 400.035065, 4.500000, 10.000000
900000018, 400.052582, 6.750000, 10.000000
900000019, 400.070007, 9.000000, 10.000000
900000021, 399.861023, -9.000000, 20.000000
900000022, 399.895294, -6.750000, 20.000000
900000023, 399.929993, -4.500000, 20.000000
900000024, 399.964935, -2.250000, 20.000000
900000025, 400.000000, 0.000000, 20.000000
900000026, 400.035065, 2.250000, 20.000000
900000027, 400.070007, 4.500000, 20.000000
900000028, 400.104706, 6.750000, 20.000000
900000029, 400.138977, 9.000000, 20.000000
900000031, 399.794189, -9.000000, 30.000000
900000032, 399.844086, -6.750000, 30.000000
900000033, 399.895294, -4.500000, 30.000000
900000034, 399.947418, -2.250000, 30.000000
900000035, 400.000000, 0.000000, 30.000000
900000036, 400.052582, 2.250000, 30.000000
900000037, 400.104706, 4.500000, 30.000000
900000038, 400.155914, 6.750000, 30.000000
900000039, 400.205811, 9.000000, 30.000000
900000041, 399.730469, -9.000000, 40.000000
900000042, 399.794189, -6.750000, 40.000000
900000043, 399.861023, -4.500000, 40.000000
900000044, 399.929993, -2.250000, 40.000000
900000045, 400.000000, 0.000000, 40.000000
900000046, 400.070007, 2.250000, 40.000000
900000047, 400.138977, 4.500000, 40.000000
900000048, 400.205811, 6.750000, 40.000000
900000049, 400.269531, 9.000000, 40.000000
900000051, 399.670837, -9.000000, 50.000000
900000052, 399.746063, -6.750000, 50.000000
900000053, 399.827271, -4.500000, 50.000000
900000054, 399.912598, -2.250000, 50.000000
900000055, 400.000000, 0.000000, 50.000000
900000056, 400.087402, 2.250000, 50.000000
900000057, 400.172729, 4.500000, 50.000000
900000058, 400.253937, 6.750000, 50.000000
900000059, 400.329163, 9.000000, 50.000000
900000061, 399.616241, -9.000000, 60.000000
900000062, 399.700104, -6.750000, 60.000000
900000063, 399.794189, -4.500000, 60.000000
900000064, 399.895294, -2.250000, 60.000000
900000065, 400.000000, 0.000000, 60.000000
900000066, 400.104706, 2.250000, 60.000000
900000067, 400.205811, 4.500000, 60.000000
900000068, 400.299896, 6.750000, 60.000000
900000069, 400.383759, 9.000000, 60.000000
900000071, 399.567444, -9.000000, 70.000000
900000072, 399.656708, -6.750000, 70.000000
900000073, 399.761871, -4.500000, 70.000000
900000074, 399.878113, -2.250000, 70.000000
900000075, 400.000000, 0.000000, 70.000000
900000076, 400.121887, 2.250000, 70.000000
900000077, 400.238129, 4.500000, 70.000000
900000078, 400.343292, 6.750000, 70.000000
900000079, 400.432556, 9.000000, 70.000000
900000081, 399.525238, -9.000000, 80.000000
900000082, 399.616241, -6.750000, 80.000000
900000083, 399.730469, -4.500000, 80.000000
900000084, 399.861023, -2.250000, 80.000000
900000085, 400.000000, 0.000000, 80.000000
900000086, 400.138977, 2.250000, 80.000000
900000087, 400.269531, 4.500000, 80.000000
900000088, 400.383759, 6.750000, 80.000000
900000089, 400.474762, 9.000000, 80.000000
900000091, 399.490265, -9.000000, 90.000000
900000092, 399.579071, -6.750000, 90.000000
900000093, 399.700104, -4.500000, 90.000000
900000094, 399.844086, -2.250000, 90.000000
900000095, 400.000000, 0.000000, 90.000000
900000096, 400.155914, 2.250000, 90.000000
900000097, 400.299896, 4.500000, 90.000000
900000098, 400.420929, 6.750000, 90.000000
900000099, 400.509735, 9.000000, 90.000000
900000101, 399.463013, -9.000000, 100.000000
900000102, 399.545502, -6.750000, 100.000000
900000103, 399.670837, -4.500000, 100.000000
900000104, 399.827271, -2.250000, 100.000000
900000105, 400.000000, 0.000000, 100.000000
900000106, 400.172729, 2.250000, 100.000000
900000107, 400.329163, 4.500000, 100.000000
900000108, 400.454498, 6.750000, 100.000000
900000109, 400.536987, 9.000000, 100.000000
900000111, 399.443939, -9.000000, 110.000000
900000112, 399.515808, -6.750000, 110.000000
900000113, 399.642853, -4.500000, 110.000000
900000114, 399.810638, -2.250000, 110.000000
900000115, 400.000000, 0.000000, 110.000000
900000116, 400.189362, 2.250000, 110.000000
900000117, 400.357147, 4.500000, 110.000000
900000118, 400.484192, 6.750000, 110.000000
900000119, 400.556061, 9.000000, 110.000000
900000121, 399.433319, -9.000000, 120.000000
900000122, 399.490265, -6.750000, 120.000000
900000123, 399.616241, -4.500000, 120.000000
900000124, 399.794189, -2.250000, 120.000000
900000125, 400.000000, 0.000000, 120.000000
900000126, 400.205811, 2.250000, 120.000000
900000127, 400.383759, 4.500000, 120.000000
900000128, 400.509735, 6.750000, 120.000000
900000129, 400.566681, 9.000000, 120.000000
900000131, 399.431305, -9.000000, 130.000000
900000132, 399.469086, -6.750000, 130.000000
900000133, 399.591064, -4.500000, 130.000000
900000134, 399.777924, -2.250000, 130.000000
900000135, 400.000000, 0.000000, 130.000000
900000136, 400.222076, 2.250000, 130.000000
900000137, 400.408936, 4.500000, 130.000000
900000138, 400.530914, 6.750000, 130.000000
900000139, 400.568695, 9.000000, 130.000000
900000141, 399.437958, -9.000000, 140.000000
900000142, 399.452454, -6.750000, 140.000000
900000143, 399.567444, -4.500000, 140.000000
900000144, 399.761871, -2.250000, 140.000000
900000145, 400.000000, 0.000000, 140.000000
900000146, 400.238129, 2.250000, 140.000000
900000147, 400.432556, 4.500000, 140.000000
900000148, 400.547546, 6.750000, 140.000000
900000149, 400.562042, 9.000000, 140.000000
900000151, 399.453125, -9.000000, 150.000000
900000152, 399.440491, -6.750000, 150.000000
900000153, 399.545502, -4.500000, 150.000000
900000154, 399.746063, -2.250000, 150.000000
900000155, 400.000000, 0.000000, 150.000000
900000156, 400.253937, 2.250000, 150.000000
900000157, 400.454498, 4.500000, 150.000000
900000158, 400.559509, 6.750000, 150.000000
900000159, 400.546875, 9.000000, 150.000000
900000161, 399.476624, -9.000000, 160.000000
900000162, 399.433319, -6.750000, 160.000000
900000163, 399.525238, -4.500000, 160.000000
900000164, 399.730469, -2.250000, 160.000000
900000165, 400.000000, 0.000000, 160.000000
900000166, 400.269531, 2.250000, 160.000000
900000167, 400.474762, 4.500000, 160.000000
900000168, 400.566681, 6.750000, 160.000000
900000169, 400.523376, 9.000000, 160.000000
900000171, 399.508057, -9.000000, 170.000000
900000172, 399.431000, -6.750000, 170.000000
900000173, 399.506805, -4.500000, 170.000000
900000174, 399.715149, -2.250000, 170.000000
900000175, 400.000000, 0.000000, 170.000000
900000176, 400.284851, 2.250000, 170.000000
900000177, 400.493195, 4.500000, 170.000000
900000178, 400.569000, 6.750000, 170.000000
900000179, 400.491943, 9.000000, 170.000000
900000181, 399.546967, -9.000000, 180.000000
900000182, 399.433563, -6.750000, 180.000000
900000183, 399.490265, -4.500000, 180.000000
900000184, 399.700104, -2.250000, 180.000000
900000185, 400.000000, 0.000000, 180.000000
900000186, 400.299896, 2.250000, 180.000000
900000187, 400.509735, 4.500000, 180.000000
900000188, 400.566437, 6.750000, 180.000000
900000189, 400.453033, 9.000000, 180.000000
900000191, 399.592804, -9.000000, 190.000000
900000192, 399.440948, -6.750000, 190.000000
900000193, 399.475647, -4.500000, 190.000000
900000194, 399.685333, -2.250000, 190.000000
900000195, 400.000000, 0.000000, 190.000000
900000196, 400.314667, 2.250000, 190.000000
900000197, 400.524353, 4.500000, 190.000000
900000198, 400.559052, 6.750000, 190.000000
900000199, 400.407196, 9.000000, 190.000000
900000201, 399.644806, -9.000000, 200.000000
900000202, 399.453125, -6.750000, 200.000000
900000203, 399.463013, -4.500000, 200.000000
900000204, 399.670837, -2.250000, 200.000000
900000205, 400.000000, 0.000000, 200.000000
900000206, 400.329163, 2.250000, 200.000000
900000207, 400.536987, 4.500000, 200.000000
900000208, 400.546875, 6.750000, 200.000000
900000209, 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
900001001, 0.000000, 2.500000, 0.000000
900001002, 0.000000, 5.000000, 0.000000
900001003, 0.000000, 7.500000, 0.000000
900001004, 0.000000, 10.000000, 0.000000
900001011, 10.000000, 2.500000, -0.009007
900001012, 10.000000, 5.000000, -0.016644
900001013, 10.000000, 7.500000, -0.021746
900001014, 10.000000, 10.000000, -0.023538
900001021, 20.000000, 2.500000, -0.017959
900001022, 20.000000, 5.000000, -0.033185
900001023, 20.000000, 7.500000, -0.043358
900001024, 20.000000, 10.000000, -0.046930
900001031, 30.000000, 2.500000, -0.026801
900001032, 30.000000, 5.000000, -0.049521
900001033, 30.000000, 7.500000, -0.064703
900001034, 30.000000, 10.000000, -0.070034
900001041, 40.000000, 2.500000, -0.035477
900001042, 40.000000, 5.000000, -0.065552
900001043, 40.000000, 7.500000, -0.085648
900001044, 40.000000, 10.000000, -0.092705
900001051, 50.000000, 2.500000, -0.043934
900001052, 50.000000, 5.000000, -0.081179
900001053, 50.000000, 7.500000, -0.106066
900001054, 50.000000, 10.000000, -0.114805
900001061, 60.000000, 2.500000, -0.052120
900001062, 60.000000, 5.000000, -0.096306
900001063, 60.000000, 7.500000, -0.125830
900001064, 60.000000, 10.000000, -0.136197
900001071, 70.000000, 2.500000, -0.059985
900001072, 70.000000, 5.000000, -0.110839
900001073, 70.000000, 7.500000, -0.144818
900001074, 70.000000, 10.000000, -0.156750
900001081, 80.000000, 2.500000, -0.067481
900001082, 80.000000, 5.000000, -0.124688
900001083, 80.000000, 7.500000, -0.162913
900001084, 80.000000, 10.000000, -0.176336
900001091, 90.000000, 2.500000, -0.074560
900001092, 90.000000, 5.000000, -0.137769
900001093, 90.000000, 7.500000, -0.180004
900001094, 90.000000, 10.000000, -0.194834
900001101, 100.000000, 2.500000, -0.081179
900001102, 100.000000, 5.000000, -0.150000
900001103, 100.000000, 7.500000, -0.195984
900001104, 100.000000, 10.000000, -0.212132
900001111, 110.000000, 2.500000, -0.087298
900001112, 110.000000, 5.000000, -0.161306
900001113, 110.000000, 7.500000, -0.210757
900001114, 110.000000, 10.000000, -0.228122
900001121, 120.000000, 2.500000, -0.092879
900001122, 120.000000, 5.000000, -0.171618
900001123, 120.000000, 7.500000, -0.224230
900001124, 120.000000, 10.000000, -0.242705
900001131, 130.000000, 2.500000, -0.097887
900001132, 130.000000, 5.000000, -0.180872
900001133, 130.000000, 7.500000, -0.236321
900001134, 130.000000, 10.000000, -0.255792
900001141, 140.000000, 2.500000, -0.102292
900001142, 140.000000, 5.000000, -0.189011
900001143, 140.000000, 7.500000, -0.246955
900001144, 140.000000, 10.000000, -0.267302
900001151, 150.000000, 2.500000, -0.106066
900001152, 150.000000, 5.000000, -0.195984
900001153, 150.000000, 7.500000, -0.256066
900001154, 150.000000, 10.000000, -0.277164
900001161, 160.000000, 2.500000, -0.109186
900001162, 160.000000, 5.000000, -0.201750
900001163, 160.000000, 7.500000, -0.263599
900001164, 160.000000, 10.000000, -0.285317
900001171, 170.000000, 2.500000, -0.111633
900001172, 170.000000, 5.000000, -0.206271
900001173, 170.000000, 7.500000, -0.269506
900001174, 170.000000, 10.000000, -0.291711
900001181, 180.000000, 2.500000, -0.113392
900001182, 180.000000, 5.000000, -0.209520
900001183, 180.000000, 7.500000, -0.273752
900001184, 180.000000, 10.000000, -0.296306
900001191, 190.000000, 2.500000, -0.114451
900001192, 190.000000, 5.000000, -0.211478
900001193, 190.000000, 7.500000, -0.276309
900001194, 190.000000, 10.000000, -0.299075
900001201, 200.000000, 2.500000, -0.114805
900001202, 200.000000, 5.000000, -0.212132
900001203, 200.000000, 7.500000, -0.277164
900001204, 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=S4R5,ELSET=PANEL
900000001,900000001,900000011,900000012,900000002
*ELGEN,ELSET=PANEL
900000001,8,1,1,20,10,8
*ELEMENT,TYPE=S4R5,ELSET=STIFF
900000201,900000005,900000015,900001011,900001001
900000221,900001001,900001011,900001012,900001002
*ELGEN,ELSET=STIFF
900000201,20,10,1
900000221,20,10,1,3,1,20
*NSET,NSET=EDGEA,GEN
900000001,900000201,10
900000009,900000209,10
*NSET,NSET=ENDAP
900000001,900000002,900000003,900000004,900000005,900000006,900000007,900000008,900000009
*NSET,NSET=ENDAS
900001001,900001002,900001003,900001004
*NSET,NSET=SYMMP
900000201,900000202,900000203,900000204,900000205,900000206,900000207,900000208,900000209
*NSET,NSET=SYMMS
900001201,900001202,900001203,900001204
*NSET,NSET=NODE3
900000005, 
*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
 900000001,900000008,1
 900000201,900000261,20
 201,261,20
*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,,900000005,3,0.2
*MONITOR,NODE=900000005,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