*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, -6.750000, 0.000000
3, 400.000000, -4.500000, 0.000000
4, 400.000000, -2.250000, 0.000000
5, 400.000000, 0.000000, 0.000000
6, 400.000000, 2.250000, 0.000000
7, 400.000000, 4.500000, 0.000000
8, 400.000000, 6.750000, 0.000000
9, 400.000000, 9.000000, 0.000000
11, 399.929993, -9.000000, 10.000000
12, 399.947418, -6.750000, 10.000000
13, 399.964935, -4.500000, 10.000000
14, 399.982452, -2.250000, 10.000000
15, 400.000000, 0.000000, 10.000000
16, 400.017548, 2.250000, 10.000000
17, 400.035065, 4.500000, 10.000000
18, 400.052582, 6.750000, 10.000000
19, 400.070007, 9.000000, 10.000000
21, 399.861023, -9.000000, 20.000000
22, 399.895294, -6.750000, 20.000000
23, 399.929993, -4.500000, 20.000000
24, 399.964935, -2.250000, 20.000000
25, 400.000000, 0.000000, 20.000000
26, 400.035065, 2.250000, 20.000000
27, 400.070007, 4.500000, 20.000000
28, 400.104706, 6.750000, 20.000000
29, 400.138977, 9.000000, 20.000000
31, 399.794189, -9.000000, 30.000000
32, 399.844086, -6.750000, 30.000000
33, 399.895294, -4.500000, 30.000000
34, 399.947418, -2.250000, 30.000000
35, 400.000000, 0.000000, 30.000000
36, 400.052582, 2.250000, 30.000000
37, 400.104706, 4.500000, 30.000000
38, 400.155914, 6.750000, 30.000000
39, 400.205811, 9.000000, 30.000000
41, 399.730469, -9.000000, 40.000000
42, 399.794189, -6.750000, 40.000000
43, 399.861023, -4.500000, 40.000000
44, 399.929993, -2.250000, 40.000000
45, 400.000000, 0.000000, 40.000000
46, 400.070007, 2.250000, 40.000000
47, 400.138977, 4.500000, 40.000000
48, 400.205811, 6.750000, 40.000000
49, 400.269531, 9.000000, 40.000000
51, 399.670837, -9.000000, 50.000000
52, 399.746063, -6.750000, 50.000000
53, 399.827271, -4.500000, 50.000000
54, 399.912598, -2.250000, 50.000000
55, 400.000000, 0.000000, 50.000000
56, 400.087402, 2.250000, 50.000000
57, 400.172729, 4.500000, 50.000000
58, 400.253937, 6.750000, 50.000000
59, 400.329163, 9.000000, 50.000000
61, 399.616241, -9.000000, 60.000000
62, 399.700104, -6.750000, 60.000000
63, 399.794189, -4.500000, 60.000000
64, 399.895294, -2.250000, 60.000000
65, 400.000000, 0.000000, 60.000000
66, 400.104706, 2.250000, 60.000000
67, 400.205811, 4.500000, 60.000000
68, 400.299896, 6.750000, 60.000000
69, 400.383759, 9.000000, 60.000000
71, 399.567444, -9.000000, 70.000000
72, 399.656708, -6.750000, 70.000000
73, 399.761871, -4.500000, 70.000000
74, 399.878113, -2.250000, 70.000000
75, 400.000000, 0.000000, 70.000000
76, 400.121887, 2.250000, 70.000000
77, 400.238129, 4.500000, 70.000000
78, 400.343292, 6.750000, 70.000000
79, 400.432556, 9.000000, 70.000000
81, 399.525238, -9.000000, 80.000000
82, 399.616241, -6.750000, 80.000000
83, 399.730469, -4.500000, 80.000000
84, 399.861023, -2.250000, 80.000000
85, 400.000000, 0.000000, 80.000000
86, 400.138977, 2.250000, 80.000000
87, 400.269531, 4.500000, 80.000000
88, 400.383759, 6.750000, 80.000000
89, 400.474762, 9.000000, 80.000000
91, 399.490265, -9.000000, 90.000000
92, 399.579071, -6.750000, 90.000000
93, 399.700104, -4.500000, 90.000000
94, 399.844086, -2.250000, 90.000000
95, 400.000000, 0.000000, 90.000000
96, 400.155914, 2.250000, 90.000000
97, 400.299896, 4.500000, 90.000000
98, 400.420929, 6.750000, 90.000000
99, 400.509735, 9.000000, 90.000000
101, 399.463013, -9.000000, 100.000000
102, 399.545502, -6.750000, 100.000000
103, 399.670837, -4.500000, 100.000000
104, 399.827271, -2.250000, 100.000000
105, 400.000000, 0.000000, 100.000000
106, 400.172729, 2.250000, 100.000000
107, 400.329163, 4.500000, 100.000000
108, 400.454498, 6.750000, 100.000000
109, 400.536987, 9.000000, 100.000000
111, 399.443939, -9.000000, 110.000000
112, 399.515808, -6.750000, 110.000000
113, 399.642853, -4.500000, 110.000000
114, 399.810638, -2.250000, 110.000000
115, 400.000000, 0.000000, 110.000000
116, 400.189362, 2.250000, 110.000000
117, 400.357147, 4.500000, 110.000000
118, 400.484192, 6.750000, 110.000000
119, 400.556061, 9.000000, 110.000000
121, 399.433319, -9.000000, 120.000000
122, 399.490265, -6.750000, 120.000000
123, 399.616241, -4.500000, 120.000000
124, 399.794189, -2.250000, 120.000000
125, 400.000000, 0.000000, 120.000000
126, 400.205811, 2.250000, 120.000000
127, 400.383759, 4.500000, 120.000000
128, 400.509735, 6.750000, 120.000000
129, 400.566681, 9.000000, 120.000000
131, 399.431305, -9.000000, 130.000000
132, 399.469086, -6.750000, 130.000000
133, 399.591064, -4.500000, 130.000000
134, 399.777924, -2.250000, 130.000000
135, 400.000000, 0.000000, 130.000000
136, 400.222076, 2.250000, 130.000000
137, 400.408936, 4.500000, 130.000000
138, 400.530914, 6.750000, 130.000000
139, 400.568695, 9.000000, 130.000000
141, 399.437958, -9.000000, 140.000000
142, 399.452454, -6.750000, 140.000000
143, 399.567444, -4.500000, 140.000000
144, 399.761871, -2.250000, 140.000000
145, 400.000000, 0.000000, 140.000000
146, 400.238129, 2.250000, 140.000000
147, 400.432556, 4.500000, 140.000000
148, 400.547546, 6.750000, 140.000000
149, 400.562042, 9.000000, 140.000000
151, 399.453125, -9.000000, 150.000000
152, 399.440491, -6.750000, 150.000000
153, 399.545502, -4.500000, 150.000000
154, 399.746063, -2.250000, 150.000000
155, 400.000000, 0.000000, 150.000000
156, 400.253937, 2.250000, 150.000000
157, 400.454498, 4.500000, 150.000000
158, 400.559509, 6.750000, 150.000000
159, 400.546875, 9.000000, 150.000000
161, 399.476624, -9.000000, 160.000000
162, 399.433319, -6.750000, 160.000000
163, 399.525238, -4.500000, 160.000000
164, 399.730469, -2.250000, 160.000000
165, 400.000000, 0.000000, 160.000000
166, 400.269531, 2.250000, 160.000000
167, 400.474762, 4.500000, 160.000000
168, 400.566681, 6.750000, 160.000000
169, 400.523376, 9.000000, 160.000000
171, 399.508057, -9.000000, 170.000000
172, 399.431000, -6.750000, 170.000000
173, 399.506805, -4.500000, 170.000000
174, 399.715149, -2.250000, 170.000000
175, 400.000000, 0.000000, 170.000000
176, 400.284851, 2.250000, 170.000000
177, 400.493195, 4.500000, 170.000000
178, 400.569000, 6.750000, 170.000000
179, 400.491943, 9.000000, 170.000000
181, 399.546967, -9.000000, 180.000000
182, 399.433563, -6.750000, 180.000000
183, 399.490265, -4.500000, 180.000000
184, 399.700104, -2.250000, 180.000000
185, 400.000000, 0.000000, 180.000000
186, 400.299896, 2.250000, 180.000000
187, 400.509735, 4.500000, 180.000000
188, 400.566437, 6.750000, 180.000000
189, 400.453033, 9.000000, 180.000000
191, 399.592804, -9.000000, 190.000000
192, 399.440948, -6.750000, 190.000000
193, 399.475647, -4.500000, 190.000000
194, 399.685333, -2.250000, 190.000000
195, 400.000000, 0.000000, 190.000000
196, 400.314667, 2.250000, 190.000000
197, 400.524353, 4.500000, 190.000000
198, 400.559052, 6.750000, 190.000000
199, 400.407196, 9.000000, 190.000000
201, 399.644806, -9.000000, 200.000000
202, 399.453125, -6.750000, 200.000000
203, 399.463013, -4.500000, 200.000000
204, 399.670837, -2.250000, 200.000000
205, 400.000000, 0.000000, 200.000000
206, 400.329163, 2.250000, 200.000000
207, 400.536987, 4.500000, 200.000000
208, 400.546875, 6.750000, 200.000000
209, 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
1011, 10.000000, 2.500000, -0.009007
1012, 10.000000, 5.000000, -0.016644
1013, 10.000000, 7.500000, -0.021746
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
1031, 30.000000, 2.500000, -0.026801
1032, 30.000000, 5.000000, -0.049521
1033, 30.000000, 7.500000, -0.064703
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
1051, 50.000000, 2.500000, -0.043934
1052, 50.000000, 5.000000, -0.081179
1053, 50.000000, 7.500000, -0.106066
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
1071, 70.000000, 2.500000, -0.059985
1072, 70.000000, 5.000000, -0.110839
1073, 70.000000, 7.500000, -0.144818
1074, 70.000000, 10.000000, -0.156750
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.176336
1091, 90.000000, 2.500000, -0.074560
1092, 90.000000, 5.000000, -0.137769
1093, 90.000000, 7.500000, -0.180004
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
1111, 110.000000, 2.500000, -0.087298
1112, 110.000000, 5.000000, -0.161306
1113, 110.000000, 7.500000, -0.210757
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
1131, 130.000000, 2.500000, -0.097887
1132, 130.000000, 5.000000, -0.180872
1133, 130.000000, 7.500000, -0.236321
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
1151, 150.000000, 2.500000, -0.106066
1152, 150.000000, 5.000000, -0.195984
1153, 150.000000, 7.500000, -0.256066
1154, 150.000000, 10.000000, -0.277164
1161, 160.000000, 2.500000, -0.109186
1162, 160.000000, 5.000000, -0.201750
1163, 160.000000, 7.500000, -0.263599
1164, 160.000000, 10.000000, -0.285317
1171, 170.000000, 2.500000, -0.111633
1172, 170.000000, 5.000000, -0.206271
1173, 170.000000, 7.500000, -0.269506
1174, 170.000000, 10.000000, -0.291711
1181, 180.000000, 2.500000, -0.113392
1182, 180.000000, 5.000000, -0.209520
1183, 180.000000, 7.500000, -0.273752
1184, 180.000000, 10.000000, -0.296306
1191, 190.000000, 2.500000, -0.114451
1192, 190.000000, 5.000000, -0.211478
1193, 190.000000, 7.500000, -0.276309
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=S8R,ELSET=PANEL
1,1,21,23,3,11,22,13,2
*ELGEN,ELSET=PANEL
1,4,2,1,10,20,4
*ELEMENT,TYPE=S8R,ELSET=STIFF
101,5,25,1022,1002,15,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
9,209,10
*NSET,NSET=ENDAP
1,2,3,4,5,6,7,8,9
*NSET,NSET=ENDAS
1001,1002,1003,1004
*NSET,NSET=SYMMP
201,202,203,204,205,206,207,208,209
*NSET,NSET=SYMMS
1201,1202,1203,1204
*NSET,NSET=NODE3
5, 
*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,4,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,,5,3,0.2
*MONITOR,NODE=5,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