/*
             OO+O+O                                                           
            O+++OO                                                            
            O+OO+O                                                            
           O+O+O+O                                                            
           OOO+OO                                                             
           OO+O+O                                                             
          O+O+OO                                                              
          OO+O+O                                                              
         O+O+O+O                                                              
         OO+O+O                                                               
        O+O+O+O                                                               
        O+OO+O                                                                
       O+O+OO                                                                 
       OO+O+O                                                                 
     OOOOOOO+O                                                                
  O+O+OOO+OO+OO                                                               
OOOOOOOOOOOOOO+OO                                                             
OOOOOOOO+OOOO+OOO                                                             
OOOOOO+OO+OOOOOOO                                                             
OOOOO       OOOOO                                                             
 O+OO       OOOO                                                              
  OOO       OOO                                                               

*/

problem description
title="CST Sample Problem (Logan #8.22)" nodes=196 elements=307

nodes
1  x=0 y=50 z=0 constraint=free
2  x=20 y=10 z=0
3  x=40 y=0 z=0
4  x=40 y=50 z=0 constraint=fixed
5  x=55 y=70 z=0 constraint=free
6  x=95 y=70 z=0
7  x=115 y=50 z=0
8  x=115 y=0 z=0 constraint=fixed
9  x=135 y=10 z=0 constraint=free
10  x=155 y=50 z=0
11  x=155 y=90 z=0
12  x=115 y=160 z=0
13  x=170 y=404 z=0
14  x=122 y=404 z=0 force=pushme
15  x=67 y=160 z=0
16  x=0 y=90 z=0
17  x=26.2595 y=117.435 z=0
18  x=30 y=50 z=0
19  x=20 y=50 z=0
20  x=10 y=50 z=0
21  x=5 y=40 z=0
22  x=10 y=30 z=0
23  x=15 y=20 z=0
24  x=30 y=5 z=0
25  x=40 y=10 z=0
26  x=40 y=20 z=0
27  x=40 y=30 z=0
28  x=40 y=40 z=0
29  x=49.2519 y=79.4871 z=0
30  x=43.5038 y=88.9741 z=0
31  x=37.7557 y=98.4612 z=0
32  x=32.0076 y=107.948 z=0
33  x=19.6947 y=110.577 z=0
34  x=13.1298 y=103.718 z=0
35  x=6.56489 y=96.8588 z=0
36  x=0 y=80 z=0
37  x=0 y=70 z=0
38  x=0 y=60 z=0
39  x=47.5 y=60 z=0
40  x=99 y=160 z=0
41  x=83 y=160 z=0
42  x=58.8519 y=151.487 z=0
43  x=50.7038 y=142.974 z=0
44  x=42.5557 y=134.461 z=0
45  x=34.4076 y=125.948 z=0
46  x=68.3333 y=70 z=0
47  x=81.6667 y=70 z=0
48  x=107 y=74 z=0
49  x=119 y=78 z=0
50  x=131 y=82 z=0
51  x=143 y=86 z=0
52  x=149.286 y=100 z=0
53  x=143.571 y=110 z=0
54  x=137.857 y=120 z=0
55  x=132.143 y=130 z=0
56  x=126.429 y=140 z=0
57  x=120.714 y=150 z=0
58  x=125 y=50 z=0
59  x=135 y=50 z=0
60  x=145 y=50 z=0
61  x=155 y=63.3333 z=0
62  x=155 y=76.6667 z=0
63  x=105 y=60 z=0
64  x=115 y=40 z=0
65  x=115 y=30 z=0
66  x=115 y=20 z=0
67  x=115 y=10 z=0
68  x=125 y=5 z=0
69  x=140 y=20 z=0
70  x=145 y=30 z=0
71  x=150 y=40 z=0
72  x=118.929 y=177.429 z=0
73  x=122.857 y=194.857 z=0
74  x=126.786 y=212.286 z=0
75  x=130.714 y=229.714 z=0
76  x=134.643 y=247.143 z=0
77  x=138.571 y=264.571 z=0
78  x=142.5 y=282 z=0
79  x=146.429 y=299.429 z=0
80  x=150.357 y=316.857 z=0
81  x=154.286 y=334.286 z=0
82  x=158.214 y=351.714 z=0
83  x=162.143 y=369.143 z=0
84  x=166.071 y=386.571 z=0
85  x=154 y=404 z=0
86  x=138 y=404 z=0
87  x=118.071 y=386.571 z=0
88  x=114.143 y=369.143 z=0
89  x=110.214 y=351.714 z=0
90  x=106.286 y=334.286 z=0
91  x=102.357 y=316.857 z=0
92  x=98.4286 y=299.429 z=0
93  x=94.5 y=282 z=0
94  x=90.5714 y=264.571 z=0
95  x=86.6429 y=247.143 z=0
96  x=82.7143 y=229.714 z=0
97  x=78.7857 y=212.286 z=0
98  x=74.8571 y=194.857 z=0
99  x=70.9286 y=177.429 z=0
100  x=15.8281 y=41.1347 z=0
101  x=21.1837 y=32.5994 z=0
102  x=29.506 y=38.2678 z=0
103  x=23.4558 y=21.9638 z=0
104  x=31.778 y=27.6323 z=0
105  x=30.8991 y=14.8505 z=0
106  x=24.9546 y=101.66 z=0
107  x=14.9486 y=93.5215 z=0
108  x=25.3499 y=89.7387 z=0
109  x=35.7513 y=85.9558 z=0
110  x=11.0391 y=83.1662 z=0
111  x=21.4404 y=79.3834 z=0
112  x=31.8418 y=75.6006 z=0
113  x=42.2431 y=71.8177 z=0
114  x=12.8673 y=70.7242 z=0
115  x=23.2687 y=66.9414 z=0
116  x=33.67 y=63.1585 z=0
117  x=11.8878 y=59.3033 z=0
118  x=40.8546 y=110.492 z=0
119  x=43.2864 y=122.454 z=0
120  x=50.9941 y=99.0945 z=0
121  x=53.4259 y=111.057 z=0
122  x=55.8578 y=123.019 z=0
123  x=58.2896 y=134.981 z=0
124  x=59.4604 y=79.4672 z=0
125  x=61.8922 y=91.4292 z=0
126  x=64.3241 y=103.391 z=0
127  x=66.7559 y=115.353 z=0
128  x=69.1878 y=127.316 z=0
129  x=71.6196 y=139.278 z=0
130  x=74.0514 y=151.24 z=0
131  x=71.845 y=79.1137 z=0
132  x=74.2769 y=91.0758 z=0
133  x=76.7087 y=103.038 z=0
134  x=79.1405 y=115 z=0
135  x=81.5724 y=126.962 z=0
136  x=84.0042 y=138.924 z=0
137  x=86.4361 y=150.886 z=0
138  x=84.3015 y=79.1137 z=0
139  x=86.7333 y=91.0758 z=0
140  x=89.1652 y=103.038 z=0
141  x=91.597 y=115 z=0
142  x=94.0289 y=126.962 z=0
143  x=96.4607 y=138.924 z=0
144  x=98.8925 y=150.886 z=0
145  x=97.8941 y=84.7022 z=0
146  x=100.326 y=96.6643 z=0
147  x=102.758 y=108.626 z=0
148  x=105.19 y=120.588 z=0
149  x=107.621 y=132.551 z=0
150  x=110.053 y=144.513 z=0
151  x=110.385 y=84.8709 z=0
152  x=112.817 y=96.833 z=0
153  x=115.249 y=108.795 z=0
154  x=117.68 y=120.757 z=0
155  x=120.112 y=132.719 z=0
156  x=124.092 y=91.0207 z=0
157  x=126.523 y=102.983 z=0
158  x=128.955 y=114.945 z=0
159  x=137.798 y=97.1705 z=0
160  x=146.214 y=77.4792 z=0
161  x=135.111 y=72.1423 z=0
162  x=138.975 y=60.5485 z=0
163  x=124.805 y=64.4131 z=0
164  x=111.384 y=66.0303 z=0
165  x=139.172 y=41.1347 z=0
166  x=125.494 y=38.2678 z=0
167  x=133.816 y=32.5994 z=0
168  x=123.222 y=27.6323 z=0
169  x=131.544 y=21.9638 z=0
170  x=124.101 y=14.8505 z=0
171  x=149.305 y=390.219 z=0
172  x=131.304 y=377.944 z=0
173  x=147.729 y=372.527 z=0
174  x=127.482 y=360.993 z=0
175  x=143.908 y=355.575 z=0
176  x=123.661 y=344.041 z=0
177  x=140.087 y=338.623 z=0
178  x=119.84 y=327.089 z=0
179  x=136.265 y=321.671 z=0
180  x=116.019 y=310.137 z=0
181  x=132.444 y=304.719 z=0
182  x=112.198 y=293.185 z=0
183  x=128.623 y=287.767 z=0
184  x=108.377 y=276.233 z=0
185  x=124.802 y=270.815 z=0
186  x=104.556 y=259.281 z=0
187  x=120.981 y=253.863 z=0
188  x=100.735 y=242.329 z=0
189  x=117.16 y=236.911 z=0
190  x=96.9134 y=225.377 z=0
191  x=113.339 y=219.959 z=0
192  x=93.0922 y=208.425 z=0
193  x=109.518 y=203.007 z=0
194  x=89.2711 y=191.473 z=0
195  x=105.696 y=186.056 z=0
196  x=87.6948 y=173.781 z=0

CSTPlaneStress elements
1  nodes=[100,101,102]   material=steel    
2  nodes=[101,103,104]   
3  nodes=[102,101,104]   
4  nodes=[104,103,105]   
5  nodes=[20,100,19]   
6  nodes=[1,21,20]   
7  nodes=[21,100,20]   
8  nodes=[22,100,21]   
9  nodes=[22,101,100]   
10  nodes=[23,101,22]   
11  nodes=[23,103,101]   
12  nodes=[2,103,23]   
13  nodes=[2,105,103]   
14  nodes=[24,105,2]   
15  nodes=[3,25,24]   
16  nodes=[25,105,24]   
17  nodes=[26,105,25]   
18  nodes=[26,104,105]   
19  nodes=[27,104,26]   
20  nodes=[27,102,104]   
21  nodes=[28,102,27]   
22  nodes=[4,18,28]   
23  nodes=[18,102,28]   
24  nodes=[19,102,18]   
25  nodes=[19,100,102]   
26  nodes=[106,107,108]   
27  nodes=[107,110,111]   
28  nodes=[108,107,111]   
29  nodes=[108,111,112]   
30  nodes=[109,108,112]   
31  nodes=[109,112,113]   
32  nodes=[111,110,114]   
33  nodes=[111,114,115]   
34  nodes=[112,111,115]   
35  nodes=[112,115,116]   
36  nodes=[113,112,116]   
37  nodes=[115,114,117]   
38  nodes=[17,33,32]   
39  nodes=[33,106,32]   
40  nodes=[34,106,33]   
41  nodes=[34,107,106]   
42  nodes=[35,107,34]   
43  nodes=[35,110,107]   
44  nodes=[16,110,35]   
45  nodes=[36,110,16]   
46  nodes=[36,114,110]   
47  nodes=[37,114,36]   
48  nodes=[37,117,114]   
49  nodes=[38,117,37]   
50  nodes=[1,20,38]   
51  nodes=[20,117,38]   
52  nodes=[19,117,20]   
53  nodes=[19,115,117]   
54  nodes=[18,115,19]   
55  nodes=[18,116,115]   
56  nodes=[4,116,18]   
57  nodes=[39,116,4]   
58  nodes=[39,113,116]   
59  nodes=[5,113,39]   
60  nodes=[29,113,5]   
61  nodes=[29,109,113]   
62  nodes=[30,109,29]   
63  nodes=[31,109,30]   
64  nodes=[31,108,109]   
65  nodes=[31,106,108]   
66  nodes=[32,106,31]   
67  nodes=[118,120,121]   
68  nodes=[119,118,121]   
69  nodes=[119,121,122]   
70  nodes=[120,125,126]   
71  nodes=[121,120,126]   
72  nodes=[121,126,127]   
73  nodes=[122,121,127]   
74  nodes=[122,127,128]   
75  nodes=[123,122,128]   
76  nodes=[123,128,129]   
77  nodes=[125,124,131]   
78  nodes=[125,131,132]   
79  nodes=[126,125,132]   
80  nodes=[126,132,133]   
81  nodes=[127,126,133]   
82  nodes=[127,133,134]   
83  nodes=[128,127,134]   
84  nodes=[128,134,135]   
85  nodes=[129,128,135]   
86  nodes=[129,135,136]   
87  nodes=[130,129,136]   
88  nodes=[130,136,137]   
89  nodes=[132,131,138]   
90  nodes=[132,138,139]   
91  nodes=[133,132,139]   
92  nodes=[133,139,140]   
93  nodes=[134,133,140]   
94  nodes=[134,140,141]   
95  nodes=[135,134,141]   
96  nodes=[135,141,142]   
97  nodes=[136,135,142]   
98  nodes=[136,142,143]   
99  nodes=[137,136,143]   
100  nodes=[137,143,144]   
101  nodes=[139,138,145]   
102  nodes=[139,145,146]   
103  nodes=[140,139,146]   
104  nodes=[140,146,147]   
105  nodes=[141,140,147]   
106  nodes=[141,147,148]   
107  nodes=[142,141,148]   
108  nodes=[142,148,149]   
109  nodes=[143,142,149]   
110  nodes=[143,149,150]   
111  nodes=[144,143,150]   
112  nodes=[146,145,151]   
113  nodes=[146,151,152]   
114  nodes=[147,146,152]   
115  nodes=[147,152,153]   
116  nodes=[148,147,153]   
117  nodes=[148,153,154]   
118  nodes=[149,148,154]   
119  nodes=[149,154,155]   
120  nodes=[150,149,155]   
121  nodes=[152,151,156]   
122  nodes=[152,156,157]   
123  nodes=[153,152,157]   
124  nodes=[153,157,158]   
125  nodes=[154,153,158]   
126  nodes=[157,156,159]   
127  nodes=[31,118,32]   
128  nodes=[31,120,118]   
129  nodes=[30,120,31]   
130  nodes=[30,125,120]   
131  nodes=[29,125,30]   
132  nodes=[29,124,125]   
133  nodes=[5,124,29]   
134  nodes=[46,124,5]   
135  nodes=[46,131,124]   
136  nodes=[47,131,46]   
137  nodes=[47,138,131]   
138  nodes=[6,138,47]   
139  nodes=[6,145,138]   
140  nodes=[48,145,6]   
141  nodes=[48,151,145]   
142  nodes=[49,151,48]   
143  nodes=[49,156,151]   
144  nodes=[50,156,49]   
145  nodes=[50,159,156]   
146  nodes=[51,159,50]   
147  nodes=[11,52,51]   
148  nodes=[52,159,51]   
149  nodes=[53,159,52]   
150  nodes=[53,157,159]   
151  nodes=[53,158,157]   
152  nodes=[54,158,53]   
153  nodes=[55,158,54]   
154  nodes=[55,154,158]   
155  nodes=[55,155,154]   
156  nodes=[56,155,55]   
157  nodes=[56,150,155]   
158  nodes=[57,150,56]   
159  nodes=[12,150,57]   
160  nodes=[12,144,150]   
161  nodes=[40,144,12]   
162  nodes=[40,137,144]   
163  nodes=[41,137,40]   
164  nodes=[41,130,137]   
165  nodes=[15,130,41]   
166  nodes=[42,130,15]   
167  nodes=[42,129,130]   
168  nodes=[42,123,129]   
169  nodes=[43,123,42]   
170  nodes=[44,123,43]   
171  nodes=[44,122,123]   
172  nodes=[44,119,122]   
173  nodes=[45,119,44]   
174  nodes=[45,118,119]   
175  nodes=[17,118,45]   
176  nodes=[32,118,17]   
177  nodes=[162,161,163]   
178  nodes=[11,160,62]   
179  nodes=[51,160,11]   
180  nodes=[50,160,51]   
181  nodes=[50,161,160]   
182  nodes=[49,161,50]   
183  nodes=[49,163,161]   
184  nodes=[49,164,163]   
185  nodes=[48,164,49]   
186  nodes=[6,63,48]   
187  nodes=[63,164,48]   
188  nodes=[7,164,63]   
189  nodes=[7,163,164]   
190  nodes=[58,163,7]   
191  nodes=[59,163,58]   
192  nodes=[59,162,163]   
193  nodes=[60,162,59]   
194  nodes=[10,61,60]   
195  nodes=[61,162,60]   
196  nodes=[61,160,162]   
197  nodes=[162,160,161]   
198  nodes=[62,160,61]   
199  nodes=[165,166,167]   
200  nodes=[167,166,168]   
201  nodes=[167,168,169]   
202  nodes=[169,168,170]   
203  nodes=[59,166,165]   
204  nodes=[58,166,59]   
205  nodes=[7,64,58]   
206  nodes=[64,166,58]   
207  nodes=[65,166,64]   
208  nodes=[65,168,166]   
209  nodes=[66,168,65]   
210  nodes=[66,170,168]   
211  nodes=[67,170,66]   
212  nodes=[8,68,67]   
213  nodes=[68,170,67]   
214  nodes=[9,170,68]   
215  nodes=[9,169,170]   
216  nodes=[69,169,9]   
217  nodes=[69,167,169]   
218  nodes=[70,167,69]   
219  nodes=[70,165,167]   
220  nodes=[71,165,70]   
221  nodes=[10,60,71]   
222  nodes=[60,165,71]   
223  nodes=[59,165,60]   
224  nodes=[171,172,173]   
225  nodes=[173,172,174]   
226  nodes=[173,174,175]   
227  nodes=[175,174,176]   
228  nodes=[175,176,177]   
229  nodes=[177,176,178]   
230  nodes=[177,178,179]   
231  nodes=[179,178,180]   
232  nodes=[179,180,181]   
233  nodes=[181,180,182]   
234  nodes=[181,182,183]   
235  nodes=[183,182,184]   
236  nodes=[183,184,185]   
237  nodes=[185,184,186]   
238  nodes=[185,186,187]   
239  nodes=[187,186,188]   
240  nodes=[187,188,189]   
241  nodes=[189,188,190]   
242  nodes=[189,190,191]   
243  nodes=[191,190,192]   
244  nodes=[191,192,193]   
245  nodes=[193,192,194]   
246  nodes=[193,194,195]   
247  nodes=[195,194,196]   
248  nodes=[86,171,85]   
249  nodes=[86,172,171]   
250  nodes=[14,87,86]   
251  nodes=[87,172,86]   
252  nodes=[88,172,87]   
253  nodes=[88,174,172]   
254  nodes=[89,174,88]   
255  nodes=[89,176,174]   
256  nodes=[90,176,89]   
257  nodes=[90,178,176]   
258  nodes=[91,178,90]   
259  nodes=[91,180,178]   
260  nodes=[92,180,91]   
261  nodes=[92,182,180]   
262  nodes=[93,182,92]   
263  nodes=[93,184,182]   
264  nodes=[94,184,93]   
265  nodes=[94,186,184]   
266  nodes=[95,186,94]   
267  nodes=[95,188,186]   
268  nodes=[96,188,95]   
269  nodes=[96,190,188]   
270  nodes=[97,190,96]   
271  nodes=[97,192,190]   
272  nodes=[98,192,97]   
273  nodes=[98,194,192]   
274  nodes=[99,194,98]   
275  nodes=[99,196,194]   
276  nodes=[15,41,99]   
277  nodes=[41,196,99]   
278  nodes=[40,196,41]   
279  nodes=[40,195,196]   
280  nodes=[12,72,40]   
281  nodes=[72,195,40]   
282  nodes=[73,195,72]   
283  nodes=[73,193,195]   
284  nodes=[74,193,73]   
285  nodes=[74,191,193]   
286  nodes=[75,191,74]   
287  nodes=[75,189,191]   
288  nodes=[76,189,75]   
289  nodes=[76,187,189]   
290  nodes=[77,187,76]   
291  nodes=[77,185,187]   
292  nodes=[78,185,77]   
293  nodes=[78,183,185]   
294  nodes=[79,183,78]   
295  nodes=[79,181,183]   
296  nodes=[80,181,79]   
297  nodes=[80,179,181]   
298  nodes=[81,179,80]   
299  nodes=[81,177,179]   
300  nodes=[82,177,81]   
301  nodes=[82,175,177]   
302  nodes=[83,175,82]   
303  nodes=[83,173,175]   
304  nodes=[84,173,83]   
305  nodes=[84,171,173]   
306  nodes=[13,85,84]   
307  nodes=[85,171,84]   

material properties
steel  E=200000 A=0 Ix=0 Iy=0 Iz=0 J=0 G=0 nu=0.25 t=10 rho=0 kappa=0

constraints
free  tx=u ty=u tz=u rx=u ry=u rz=u
fixed  tx=c ty=c tz=u rx=u ry=u rz=u

forces
pushme  Fx=700 Fy=0 Fz=0 Mx=0 My=0 Mz=0

end