蔣友寶,謝銘武,廖 強
(1.長沙理工大學 土木與建筑學院,長沙410004;2.東南大學 混凝土及預應力混凝土結(jié)構(gòu)教育部重點實驗室,南京210096;3.湖南湘創(chuàng)建材科技有限公司,長沙410000)
鋼拱結(jié)構(gòu)因具有良好的力學性能及優(yōu)美流暢的建筑造型,已在工程結(jié)構(gòu)中得到了廣泛的應用,如大跨度的公共建筑、橋梁等。由于鋼拱受力一般較鋼壓彎構(gòu)件復雜一些,因此很多學者對其平面內(nèi)穩(wěn)定承載力設計方法進行了深入研究。Komatsu等[1]較早研究了兩鉸拋物線鋼拱的極限承載力。隨后,Kuranishi等[2]提出了一個針對兩鉸鋼拱的承載力計算公式,該公式由兩段不同的彎矩和軸力相關方程組成。Pi等[3]分析了鋼拱在全跨荷載、半跨荷載、拱頂集中荷載和拱四分點集中荷載單獨作用下的平面內(nèi)穩(wěn)定極限承載力,并基于澳大利亞結(jié)構(gòu)設計規(guī)范的柱子曲線給出了相應的設計公式。中國學者在鋼拱平面內(nèi)穩(wěn)定承載力設計方法方面也開展了較多研究。郭彥林等[4-6]采用試驗研究與彈塑性有限元分析相結(jié)合的思路,對不同截面形式的拋物線鋼拱平面內(nèi)穩(wěn)定承載力進行了分析;得到了殘余應力對壓彎圓弧鋼拱平面內(nèi)穩(wěn)定承載力的影響較小的結(jié)論;在此基礎上提出了相應的設計方法。此外,童根樹等[7]考慮實際工程中鋼拱水平支座的非剛性性質(zhì),給出了水平彈性支承圓弧鋼拱的平面內(nèi)穩(wěn)定承載力設計方法。目前,中國已制定了相應的拱形鋼結(jié)構(gòu)技術(shù)規(guī)程[8]。但需說明的是,現(xiàn)行規(guī)程中的鋼拱平面內(nèi)承載力驗算式主要是在單一荷載工況下鋼拱承載力分析結(jié)果的基礎上得到的[9];且承載力分析時主要是按一致缺陷(按第1階屈曲模態(tài)確定的缺陷)分布來考慮缺陷的不利影響,而已有較多的研究[10-11]表明,一致缺陷分布并非結(jié)構(gòu)的最不利缺陷,甚至可能會高估結(jié)構(gòu)的承載力。例如,蔣友寶[10]對一正交拱式結(jié)構(gòu)在不同缺陷分布情形下的承載力進行了分析,結(jié)果表明當缺陷分布按前8階屈曲模態(tài)分別考慮時,第8階屈曲模態(tài)對應的缺陷最不利,其承載力最低,而按第1階屈曲模態(tài)對應的缺陷考慮的承載力則要高出約20%。設計鋼拱結(jié)構(gòu)時一般需考慮全跨荷載和半跨荷載的作用組合,其中半跨分布荷載可為雪荷載[12]等。蔣友寶等[13-14]研究了全跨和半跨荷載組合下鋼拱結(jié)構(gòu)荷載效應函數(shù)模型及承載力可靠性設計方法。結(jié)果初步表明,由于半跨荷載對鋼拱承載性能的不利影響,導致按現(xiàn)行荷載分項系數(shù)設計的鋼拱結(jié)構(gòu)的安全性可能偏低。本文以全跨和半跨荷載組合作用下的兩鉸拋物線鋼拱結(jié)構(gòu)為研究對象,對比分析了不同截面和矢跨比下現(xiàn)行規(guī)程方法、文獻[2]建議方法、考慮二階效應的強度設計方法、基于完善和一致缺陷拱模型的有限元分析方法的承載力計算結(jié)果的差異,得到了一些有益的結(jié)論。
鋼拱一般承受軸壓力和平面內(nèi)彎矩的共同作用,現(xiàn)行規(guī)程[8]給出的實腹式截面圓弧及拋物線鋼拱的平面內(nèi)穩(wěn)定承載力計算公式為
其中:N和M 為拱中截面的最大軸力和最大彎矩;φ為軸心受壓拱的平面內(nèi)穩(wěn)定系數(shù);A為截面面積;f為鋼材抗壓和抗拉屈服強度;γx為截面塑性發(fā)展系數(shù),可按現(xiàn)行鋼結(jié)構(gòu)設計規(guī)范[15]的規(guī)定取值;Wx為拱軸線平面內(nèi)彎曲的毛截面模量;α為一與支承條件、截面形式有關的系數(shù),可按規(guī)程中表6.2.3確定,對于兩鉸拋物線箱型截面拱,其取值為0.83。
該方法選用鋼拱四分點處的截面彎矩M和軸力N來驗算承載力。當N/Ny>ncr時,驗算式為
當N/Ny≤ncr時,驗算式為
式中:ncr=(1-bmcr-amcr2)/c;a、b、c、α、β 和 mcr等系數(shù)的計算式參見文獻[2];而Ny和My分別為鋼拱截面受壓屈服軸力和邊緣屈服彎矩,可分別計算為Ny=Af和My=Wxf。
設有一兩鉸拋物線鋼拱結(jié)構(gòu),其跨度為60m,彈性模量Es=206GPa,按Q345鋼材考慮,屈服強度值為315MPa,矢跨比f1=0.15,箱型截面尺寸b×h×t1×t2(寬度×高度×翼緣厚度×腹板厚度)??紤]3種截面尺寸,見表1。
表1 不同截面尺寸及對應的長細比
其中長細比λ按式(4)計算。
其中:S為拱軸弧線長度;ix為鋼拱平面內(nèi)彎曲的截面回轉(zhuǎn)半徑;fy為鋼材屈服強度標準值。
對于3種截面的鋼拱結(jié)構(gòu),規(guī)程方法和文獻[2]方法對應的彎矩軸力相關方程見圖1。可知,在鋼拱截面給定的情形下,當彎矩值大于某一很小的數(shù)值時,規(guī)程方法的N-M 相關曲線在文獻[2]建議方法的N-M 相關曲線的上方,即此時按文獻[2]建議方法設計的鋼拱驗算截面所能承受的內(nèi)力要小一些;但當彎矩很小時,規(guī)程方法的N-M相關曲線卻反而位于文獻[2]建議方法的N-M相關曲線的下方,即此時按規(guī)程方法設計的鋼拱驗算截面所能承受的內(nèi)力要小一些。
圖1 規(guī)程方法和文獻[2]方法的N-M 相關曲線比較
而兩種方法選用的驗算截面并不相同,規(guī)程方法采用內(nèi)力最大值截面,文獻[2]方法采用四分點截面。顯然規(guī)程方法采用的內(nèi)力更為不利,因此當彎矩很小時按規(guī)程方法設計所能承受的極限荷載較小,要保守一些,見4.2節(jié);而當彎矩較大時,兩種設計方法的保守程度還需進一步深入分析方能確定,具體見4.3節(jié)。
對于承受全跨荷載g和半跨荷載q的兩鉸拋物線鋼拱結(jié)構(gòu),見圖2。當采用一階分析、內(nèi)力線性組合的原則來求得截面軸力和彎矩值時,一般有
其中:c1(q/g)和c2(q/g)分別表示加載比例為q/g時彎矩和軸力對應的荷載效應系數(shù);F為加載幅值。
圖2 全跨和半跨荷載組合下兩鉸鋼拱結(jié)構(gòu)計算模型
由于規(guī)程方法和文獻[2]方法的承載力驗算式均是按一階分析內(nèi)力來考慮的,因此將式(5)和式(6)分別代入式(1)、式(2)和式(3)即可求得加載比例為q/g時的極限荷載值。
該方法是指在建立鋼拱有限元模型時,按照拱軸線線形和最低階屈曲模態(tài)分布的變形來確定鋼拱節(jié)點的位置。文中一致缺陷拱模型中考慮的最大缺陷幅值為S/500。
選用ANSYS軟件進行分析,為提高分析精度,采用能同時考慮幾何非線性和材料非線性特性的beam188單元來模擬鋼拱,其中鋼材按理想彈塑性材料考慮。由于按中國建筑鋼材考慮的殘余應力對壓彎鋼拱平面內(nèi)承載力的影響較小[6],因而分析時暫不考慮殘余應力的影響。
考慮到一致缺陷拱模型的承載力并非最低,如文獻[2]建議的承載力驗算式(2)和式(3)即是在完善拱模型承載力計算結(jié)果的基礎上得到的。
作為對比,文中還建立了完善拱的有限元模型。有限元分析仍選用ANSYS軟件,其建模、求解思路與一致缺陷拱模型相同。
德國鋼結(jié)構(gòu)規(guī)范DIN18800[16]建議了一種鋼拱平面內(nèi)穩(wěn)定承載力設計方法,即對鋼拱進行二階彈性分析,并使計入等效幾何缺陷影響的最大應力不超過0.9f。
這種方法實際上是一種強度設計方法[6]。考慮到殘余應力的影響不大,且分析時尚未考慮截面進入塑性,因而其最大應力可放寬至1.0f??紤]二階效應的強度設計方法對應的承載力驗算式為
其中:σ為截面正應力;N2nd和M2nd為考慮二階效應的截面軸力和彎矩,其余變量同式(1)。
由于一致缺陷可能并非最不利,因此下文采用該方法計算極限承載力時按完善鋼拱考慮。
此處對比各種計算方法的差異時,半跨荷載先考慮為半跨滿布荷載。對于兩鉸拋物線鋼拱,其承載力的影響因素主要有矢跨比和截面等參數(shù)。為此,以1.3節(jié)中的截面2鋼拱結(jié)構(gòu)為基準模型,通過變化截面和矢跨比參數(shù),得到多個鋼拱模型。然后采用5種方法:規(guī)程方法、文獻[2]方法、基于一致缺陷拱模型的有限元方法、基于完善拱模型的有限元方法和考慮二階效應的強度設計方法分別計算其極限荷載。
對于表1中3種不同截面的鋼拱結(jié)構(gòu),其在不同荷載組合比例下的極限荷載計算結(jié)果見圖3。
由圖3可知,在各種截面下不同方法計算得到的極限荷載結(jié)果差異較大。但總體來說,極限荷載計算結(jié)果大致可分為3類:規(guī)程方法和基于一致缺陷拱模型的有限元方法得到的極限承載力較為接近,為第1類;考慮二階效應的強度設計方法與文獻[2]建議方法的計算結(jié)果較為接近,為第2類;基于完善拱模型的有限元方法得到的計算結(jié)果與上兩類結(jié)果有不小的差異,可歸為第3類。
圖3 不同截面下按各種方法計算的極限荷載對比
當半跨滿布荷載占高比重時,即q/g大于某一臨界值r0時,第1類計算結(jié)果最高,第3類計算結(jié)果次之(截面3稍微不滿足),第2類計算結(jié)果最低。對于截面1、截面2和截面3,這一臨界值r0均約為0.20。而當全跨荷載占高比重時,即q/g小于某一值時,第3類計算結(jié)果最高,第2類計算結(jié)果次之,第1類計算結(jié)果最低。
在基準模型中拱矢跨比f1為0.15,當其變化為0.1和0.2時,其極限荷載的計算結(jié)果見圖4。
從圖4中可以看出,5種方法得到的鋼拱結(jié)構(gòu)極限荷載的變化趨勢與矢跨比f1=0.15時的情形(圖3(b))基本相同。但對于矢跨比f1=0.2時的情形,臨界值r0約為0.10。
圖4 不同矢跨比下按各種方法計算的極限荷載對比
依據(jù)上述不同截面、矢跨比下的兩鉸拋物線鋼拱結(jié)構(gòu)極限荷載計算結(jié)果,可知5種方法得到的計算結(jié)果雖有較大差異,但大致可分為3類。
由于精確地確定鋼拱最不利缺陷分布較為困難,因此當無法得知何種缺陷分布較為不利時,基于完善拱模型的有限元方法的計算結(jié)果可作為一個較好的比較參考值。而參數(shù)分析結(jié)果表明:當半跨滿布荷載占高比重時,規(guī)程方法得到的極限荷載值會大于同荷載比例下的參考值,最高可超出約20%。這表明此時規(guī)程方法計算得到的極限荷載值偏于不保守,因而按規(guī)程方法設計鋼拱結(jié)構(gòu)將偏于不安全。
另外,當全跨荷載占高比重時,如僅有全跨荷載(q/g=0),規(guī)程方法得到的極限荷載值最低,這表明此時按規(guī)程方法設計鋼拱結(jié)構(gòu)將偏于保守。
需說明的是,規(guī)程方法和文獻[2]方法給出的平面內(nèi)穩(wěn)定承載力計算式也是在大量鋼拱結(jié)構(gòu)承載力有限元分析結(jié)果的基礎上得到的。因此文中討論的5種方法都可歸為有限元方法。
由于影響鋼拱結(jié)構(gòu)極限荷載值的因素有初始缺陷、殘余應力和塑性發(fā)展程度等,而全面考慮這些因素精確求解鋼拱結(jié)構(gòu)極限荷載值的難度非常大,因此分析時往往會根據(jù)實際需要而作出一些簡化。上述5種方法對這些因素的考慮情形見表2??梢姴煌椒▽@些因素的考慮方式均不完全相同,因此得到的計算結(jié)果必然有較大差異。
表2 不同方法考慮的影響因素對比
已有的較多研究[6,10-11]表明,初始缺陷對壓彎受力結(jié)構(gòu)的承載力影響較大。因此根據(jù)這一特點,可首先將上述5種方法歸為兩類,即不考慮缺陷的文獻[2]方法、基于完善拱模型的有限元方法和考慮二階效應的強度設計方法為一類;考慮初始缺陷的規(guī)程方法和基于一致缺陷拱模型的有限元方法為一類。再考慮殘余應力和塑性發(fā)展對結(jié)構(gòu)承載力的影響,其中殘余應力為不利因素而塑性發(fā)展為有利因素;顯然若只考慮兩個因素的其中之一,則要么會過低估計結(jié)構(gòu)承載力,要么會過高估計結(jié)構(gòu)承載力,即導致承載力計算結(jié)果會有一個較大的偏差;從這個角度考慮可將不考慮殘余應力而僅考慮塑性發(fā)展的基于完善拱模型的有限元方法的計算結(jié)果單獨列選出來作為一類。因此最終可判定5種方法的計算結(jié)果大致可分為3類。
規(guī)程方法較基于完善拱模型的有限元方法保守的原因較為簡單,即當全跨荷載占高比重時不考慮初始缺陷的完善拱結(jié)構(gòu)的受力狀況接近于純壓受力狀態(tài),彎矩很小,這是一種較理想的情形,顯然高估了鋼拱結(jié)構(gòu)的極限荷載。因而以下主要分析全跨荷載占高比重時規(guī)程方法較文獻[2]方法保守的原因。
規(guī)程方法和文獻[2]方法都是基于彎矩軸力相關方程來計算極限荷載值的。區(qū)別在于,規(guī)程方法的承載力計算式(1)是按各截面最大軸力和彎矩考慮的,而文獻[2]方法卻是按l/4截面軸力和彎矩考慮的,因此即使對于同一鋼拱結(jié)構(gòu),兩種方法計算承載力時所用的荷載效應系數(shù)值也會不同。為此定義兩個系數(shù)rM和rN,其計算式分別為
可見這兩個系數(shù)反映了截面最大內(nèi)力與l/4截面內(nèi)力的比值。表3給出了f1=0.15的鋼拱結(jié)構(gòu)基準模型和將f1變化為0.1和0.2情形時rM和rN的計算結(jié)果。
由表3可知,鋼拱結(jié)構(gòu)最大軸力和l/4截面軸力的效應系數(shù)比值隨荷載比例q/g的變化較明顯,即rN可取區(qū)間[1.05,1.37]內(nèi)的多個分散的相應值。而最大彎矩和l/4截面彎矩的效應系數(shù)比值則隨q/g的變化較不明顯,基本上僅取兩個值,即q/g=0時取較大值,約為1.33,其余荷載比例值下均約為1.0。
表3 rM和rN隨荷載比例的變化
可見,當全跨荷載占高比重時,按最大彎矩和軸力考慮的規(guī)程方法對應的設計內(nèi)力(尤其是彎矩)取值較文獻[2]方法大。但需指出的是,全跨荷載占高比重時,總體上鋼拱結(jié)構(gòu)受力將以受壓為主,總彎矩將較小,這樣規(guī)程方法對應的設計抗力取值將較?。ㄒ妶D1)。綜合設計內(nèi)力和抗力兩個方面,可知采用規(guī)程方法設計鋼拱結(jié)構(gòu)將較文獻[2]方法保守較多。
采用的分析思路為:在完善拱模型上施加規(guī)程方法計算所得的極限荷載,然后考慮幾何非線性效應進行有限元分析,檢驗完善拱結(jié)構(gòu)中應力的大小。以文中兩鉸拋物線鋼拱基準模型為例,圖5給出了不同q/g時完善拱模型的應力分布情形。
可看出,在規(guī)程方法計算所得的極限荷載作用下,當q/g在[0.2,∞]內(nèi)取值時,完善拱結(jié)構(gòu)的較多區(qū)域內(nèi)的應力比均超過1.0較多,最高達1.4(q/g=0.5和q/g=1.0)。這表明當施加的荷載值尚未達到規(guī)程方法計算的極限荷載時,完善拱截面已開始進入塑性,因而半跨滿布荷載占高比重時規(guī)程方法計算結(jié)果要偏于不保守。同時,由于規(guī)程方法是按一致缺陷模態(tài)來考慮的,這亦表明半跨滿布荷載占高比重時完善拱結(jié)構(gòu)的受力較一致缺陷拱結(jié)構(gòu)不利一些。
圖5 規(guī)程方法求得的極限荷載下完善拱的應力分布
由于半跨滿布荷載占高比重的情形涵蓋了大部分荷載組合比例值,因此按現(xiàn)行規(guī)程設計兩鉸拋物線鋼拱承載力會偏于不安全的現(xiàn)象應引起重視。
為進一步明確不同計算方法的差異,繼續(xù)考察半跨荷載為按最不利影響線布置的情形。
一般來說,鋼拱結(jié)構(gòu)受較大彎矩時較為不利,因此為確定最不利荷載組合下的半跨荷載,分別選用l/8、l/4、3l/8、l/2截面作為控制截面,以截面的彎矩影響線作為確定最不利半跨荷載位置的依據(jù)。各截面的彎矩影響線如圖6所示。
圖6 控制截面彎矩影響線
可知,當半跨荷載按影響線加載時,l/4截面的彎矩最大(對于等截面鋼拱,此時較為不利),對應的半跨荷載布置區(qū)域為0~0.43l。
為明確比較重點,此處僅對比規(guī)程方法、完善結(jié)構(gòu)有限元方法和一致缺陷結(jié)構(gòu)有限元方法的差異。當半跨荷載按上述最不利荷載布置時,對于表1中三種不同截面的鋼拱結(jié)構(gòu),其在不同荷載組合比例下的極限荷載計算結(jié)果如圖7所示。
圖7 不同截面下按各種方法計算的極限荷載對比
可見,與半跨滿布荷載情形類似,考慮半跨荷載的最不利布置后,規(guī)程方法計算結(jié)果與一致缺陷有限元分析結(jié)果仍較為接近;全跨荷載占高比重時規(guī)程結(jié)果仍會保守較多;半跨滿布荷載占高比重時規(guī)程結(jié)果仍然較不保守。
但需說明的是,當鋼拱長細比較小時,如圖7(c),半跨不利荷載占高比重情形下規(guī)程結(jié)果較不保守的程度會變得不明顯一些。另外,由于上述半跨不利荷載(0~0.43l)布置與半跨滿布荷載(0~0.5l)布置較為接近,因而按半跨不利荷載情形計算得到的極限荷載值僅略低于半跨滿布荷載下的對應值。
考慮全跨荷載與半跨荷載的組合,對兩鉸拋物線鋼拱承載力計算的多種方法進行了對比分析,得出以下結(jié)論:
1)按規(guī)程方法計算得到兩鉸拋物線鋼拱結(jié)構(gòu)極限荷載值與考慮一致缺陷、幾何非線性和材料非線性影響,但不考慮殘余應力影響的有限元分析結(jié)果較為接近。
2)當全跨荷載占高比重時,按規(guī)程方法得到的極限荷載值在文中分析參數(shù)范圍內(nèi)均較保守。
3)當半跨荷載占高比重時,按規(guī)程方法設計兩鉸拋物線鋼拱結(jié)構(gòu)會偏于不安全。原因主要是此時一致缺陷分布并非最不利的缺陷形式,其對應的極限荷載值甚至會高于完善拱模型的極限荷載值。
[1]Komatsu S,Shinke T.Ultimate load carrying capacity of steel arches [J].Journal of Structural Division,1977,103(12):2323-2326.
[2]Kuranishi S,Yabuki T.Ultimate strength design criterion for two-h(huán)inged steel arch structures [J].Structural Engineering/Earthquake Engineering Proceedings of Japan Society of Civil Engineers,1984,1(2):115-123.
[3]Pi Y L,Trahair N S.In-plane buckling and design of steel arches [J].Journal of Structural Engineering,ASCE,1999,125(11):1291-1298.
[4]林冰,郭彥林.純壓拋物線拱平面內(nèi)穩(wěn)定性及設計方法研究[J].建筑結(jié)構(gòu)學報,2009,30(3):103-111.Lin B,Guo Y L.In-plane stability behavior and application of parabolic arches under pure compression [J].Journal of Building Structures,2009,30(3):103-111.
[5]郭彥林,林冰,郭宇飛.焊接工字形截面拋物線拱平面內(nèi)穩(wěn)定性試驗研究[J].建筑結(jié)構(gòu)學報,2009,30(3):95-102.Guo Y L,Lin B,Guo Y F.Experimental study on in-plane stability of welded I-section steel parabolic arches [J].Journal of Building Structures,2009,30(3):95-102.
[6]郭彥林,林冰,郭宇飛.壓彎圓弧拱平面內(nèi)穩(wěn)定承載力設計方法的理論與試驗研究[J].土木工程學報,2011,44(3):8-15.Guo Y L,Lin B,Guo Y F.Theoretical and experimental studies of in-plane stability design of circular arches subjected to axial force and moment [J].China Civil Engineering Journal,2011,44(3):8-15.
[7]楊洋,童根樹.水平彈性支承圓弧鋼拱的平面內(nèi)極限承載力研究[J].工程力學,2012,29(3):45-54.Yang Y,Tong G S.Study for in-plane ultimate strength of steel circular arches with horizontal spring supports[J].Engineering Mechanics,2012,29(3):45-54.
[8]JGJ/T 249—2011拱形鋼結(jié)構(gòu)技術(shù)規(guī)程[S].
[9]郭彥林,竇超.鋼拱結(jié)構(gòu)設計理論與我國鋼拱結(jié)構(gòu)技術(shù)規(guī)程[J].鋼結(jié)構(gòu),2009,24(5):59-70.Guo Y L,Dou C.Design theory and Chinese technical specification for steel arch structure[J].Steel Construction,2009,24(5):59-70.
[10]蔣友寶.正交拱式結(jié)構(gòu)的穩(wěn)定分析[J].建筑科學,2008,24(9):14-18.Jiang Y B.Stability analysis of the ortho-arched structure[J].Building Science,2008,24(9):14-18.
[11]盧家森,張其林.球面網(wǎng)殼最不利幾何缺陷的凸集和概率模型[J].工程力學,2013,30(7):100-104,121.Lu J S,Zhang J L.Convex and probabilistic models of uncertainties in initial geometric imperfections of latticed shells[J].Engineering Mechanics,2013,30(7):100-104,121.
[12]Hochstenbach F M,Irwin P A,Gamble S L.Parametric studies of unbalanced snow loads on arched roofs[C]//Blandford G E.Structures 2004:Building on the Past,Securing the Future,Nashville,Tennessee,United States,ASCE,2004:1-11.
[13]蔣友寶,賀藝華.全跨和半跨荷載組合下荷載效應函數(shù)模型研究[J].工程力學,2011,28(2):123-128.Jiang Y B,He Y H.Research on load effect function model under full-span and half-span loading [J].Engineering Mechanics,2011,28(2):123-128.
[14]蔣友寶,廖國宇,謝銘武.鋼筋混凝土框架柱和輕鋼拱結(jié)構(gòu)失效方程復雜特性與設計可靠度[J].建筑結(jié)構(gòu)學報,2014,35(4):192-198.Jiang Y B,Liao G Y,Xie M W.Complex failure function and design reliability for RC frame columns and lightweight steel arch structures [J].Journal of Building Structures,2014,35(4):192-198.
[15]GB 50017—2003鋼結(jié)構(gòu)設計規(guī)范[S].
[16]DIN-18800-2-1990Structural steelwork,analysis of safety against buckling of linear members and frames[S].1990.