譚福穎 喬 玲 韓曉林

(東南大學(xué)土木工程學(xué)院,南京 210096)(東南大學(xué)江蘇省工程力學(xué)分析重點(diǎn)實(shí)驗(yàn)室,南京 210096)

殼體結(jié)構(gòu)在航空航天工程和土木工程中得到了廣泛應(yīng)用,而圓柱殼屈曲問(wèn)題一直是殼體研究中最為關(guān)注的課題之一,這類(lèi)結(jié)構(gòu)在承受荷載時(shí)大多在強(qiáng)度破壞前就已發(fā)生失穩(wěn)破壞,因此研究影響薄壁殼體穩(wěn)定性的因素一直受到研究者們的重視.

對(duì)于廣義梁理論,國(guó)內(nèi)外學(xué)者進(jìn)行了大量研究.1989年德國(guó)學(xué)者Schardt[1]首次提出廣義梁理論,此理論是在經(jīng)典的Vlasov梁理論的基礎(chǔ)上進(jìn)行拓展,通過(guò)研究局部坐標(biāo)系下平面彎曲截面下變形,從而分析全局域和局部域下薄壁棱柱形矩形結(jié)構(gòu)構(gòu)件的屈曲行為.由于Schardt的研究成果主要以德文的形式發(fā)表,在一定程度上影響了該理論的傳播和應(yīng)用.Davies等[2]最先以英文的方式將廣義梁理論進(jìn)行闡述,并應(yīng)用該理論對(duì)冷彎薄壁型鋼(主要是C形和Z形截面)的相關(guān)屈曲進(jìn)行了系統(tǒng)的理論分析和試驗(yàn)驗(yàn)證.文獻(xiàn)[3-7]將廣義梁理論擴(kuò)展到多種材料形式、考慮剪切扭轉(zhuǎn)變形、幾何非線(xiàn)性、一階靜力和二階穩(wěn)定以及任意分支開(kāi)口截面下的薄壁桿件分析中,并給出了廣義梁有限元模型以及利用數(shù)值計(jì)算方法求解出了廣義梁方程.我國(guó)學(xué)者李開(kāi)禧[8]1990年提出過(guò)類(lèi)似的理論,研究在中線(xiàn)為直線(xiàn)的假定條件下,利用圖解法求解薄壁桿件的截面變形問(wèn)題.與國(guó)外相比,國(guó)內(nèi)對(duì)于廣義梁理論這種充分考慮截面變形梁理論的研究仍然滯后,國(guó)內(nèi)相關(guān)的研究鮮有報(bào)道.

本文將廣義梁理論引入并推廣到薄壁圓形截面柱殼的穩(wěn)定性分析中,開(kāi)展了對(duì)中空薄壁圓柱殼屈曲的研究,從一階線(xiàn)性分析拓展到各階屈曲分析,通過(guò)獲得屈曲應(yīng)力和與之對(duì)應(yīng)的臨界軸長(zhǎng)表達(dá)式,研究了不同長(zhǎng)細(xì)比下臨界應(yīng)力隨軸長(zhǎng)、厚度的變化關(guān)系,為薄壁圓柱殼穩(wěn)定性分析提供了一個(gè)新方法.

1 薄壁圓柱殼穩(wěn)定性理論分析

1.1 幾何方程

如圖1所示,全局坐標(biāo)系XYZ中,圓柱殼半徑為r,厚度為t,軸長(zhǎng)為L(zhǎng).在中性面上建立一個(gè)局部坐標(biāo)系x,θ和z.其中軸向x∈[0，L],圓周向θ∈[0，2π],徑向z∈[-t/2，+t/2].分別用u,v,w表示局部坐標(biāo)系下點(diǎn)的各向位移.

圖1 全局坐標(biāo)系及局部坐標(biāo)系下圓柱殼參數(shù)示意圖

基于Love-Kirchhoff 假設(shè)和薄殼理論,應(yīng)變位移關(guān)系為[9]

(1)

由經(jīng)典梁理論可知,中性橫截面上任一點(diǎn)位移可表示為截面位移與截面翹曲的組合,因此u,v,w可表示成

(2)

式中,x,θ為獨(dú)立的自變量和位移振幅函數(shù);φk(x)表示截面翹曲程度;位移函數(shù)uk(θ)表示在同一截面上的截面軸向位移.將式(2)代入式(1),根據(jù)線(xiàn)性應(yīng)變和非線(xiàn)性應(yīng)變的定義,將應(yīng)變劃分為與截面位移u,v,w成線(xiàn)性關(guān)系的線(xiàn)性應(yīng)變,以及與截面位移u,v,w成非線(xiàn)性關(guān)系的非線(xiàn)性應(yīng)變2部分.整理得

(3)

(4)

1.2 能量變分方程

薄膜應(yīng)變能U可通過(guò)由薄膜內(nèi)力在中面薄膜變形中所做內(nèi)功計(jì)算而得.取微分單元體,由其薄膜力的元功可得能量變分方程為

(5)

(6)

其中,Q11=Q22=E/(1-ν2),Q12=νQ11,Q33=G，E為彈性模量,v為泊松比.

(7)

(8)

殼體穩(wěn)定性分析的首要問(wèn)題是臨界荷載問(wèn)題,而該問(wèn)題是由前屈曲平衡性態(tài)來(lái)決定的.當(dāng)殼體進(jìn)入屈曲狀態(tài),會(huì)出現(xiàn)非線(xiàn)性薄膜應(yīng)力,使之偏離到屈曲后附加狀態(tài),即非線(xiàn)性狀態(tài),但附加狀態(tài)偏離屈曲狀態(tài)的位移無(wú)限小,因此,非線(xiàn)性項(xiàng)應(yīng)變可以忽略不計(jì).故而,在忽略非線(xiàn)性應(yīng)變前提下,將上述各式代入到能量變分方程(5)中，可得

(9)

由以上分析可知,能量泛函δU是一個(gè)只與軸向u(x,θ)有關(guān)的函數(shù),求解u(x,θ)的表達(dá)式即可得到各未知參量.由于u(x,θ)可以分解為2個(gè)正交函數(shù)uk(θ)和φk(x)的乘積,因此問(wèn)題歸結(jié)為求解uk(θ)和φk(x)的表達(dá)式.

1.3 各階屈曲模態(tài)

根據(jù)uk(θ)和φk(x)的正交性可知,位移函數(shù)中u(x,θ)中位移振幅函數(shù)φk(x)表示截面翹曲程度,uk(θ)表示在同一截面上的截面位移,因此獲得uk(θ)即可得到殼體橫截面屈曲模態(tài).

由廣義梁理論可知,若要求解能量泛函δU,矩陣Cik,Bik必須對(duì)角化.矩陣Cik,Bik對(duì)角化，即為uk,uk,θθ,uk,θθθθ滿(mǎn)足獨(dú)立正交[11]條件,可表示為

(10)

由正交條件可知，uk(θ)可用三角周期函數(shù)表示為

uk(θ)=rsinkθ,uk(θ)=rcoskθ

(11)

對(duì)應(yīng)于一個(gè)給定的k值,uk有正弦和余弦函數(shù)2種表達(dá)形式,且兩者相互正交.由工程實(shí)際可知,相同的屈曲載荷下由于位移分叉存在2種形態(tài)類(lèi)似的屈曲模態(tài).引入同一屈曲波數(shù)m可對(duì)應(yīng)于不同的階數(shù)k值,得到2種相似的屈曲模態(tài).因此,當(dāng)表達(dá)式設(shè)為正弦函數(shù)時(shí),k=2m;設(shè)為余弦函數(shù)時(shí),k=2m+1.因此,筒殼各階截面位移表達(dá)式為

(12a)

(12b)

圖2為薄壁圓柱殼2階～13階模態(tài).由圖可知,在同一屈曲載荷下由于位移分叉存在2種形態(tài)類(lèi)似的屈曲模態(tài),如2階、3階同為屈曲模態(tài),4階、5階同為扭轉(zhuǎn)模態(tài)等,也再次證明所假設(shè)的正交對(duì)稱(chēng)位移函數(shù)是可取的.

此外,在實(shí)際工程中需要考慮軸向變形和前屈曲問(wèn)題時(shí),必須要得到其軸向變形模態(tài).對(duì)比軸向壓縮模態(tài)的力學(xué)性質(zhì)不難發(fā)現(xiàn),在發(fā)生壓縮模態(tài)時(shí),只存在軸向位移.因此有

ue=1,ve=0,we=0

(13)

式中,下標(biāo)e表示軸向伸縮模態(tài)下的位移分量.結(jié)合式(12)可以得到,當(dāng)m=0時(shí)，第1階軸壓模態(tài)為

m=0,k=1,u1=1,v1=0,w1=0

(14)

圖2 薄壁圓柱殼各階屈曲模態(tài)

柱殼屈曲時(shí)會(huì)產(chǎn)生此模態(tài).因此,各階模態(tài)還必須包含一個(gè)軸對(duì)稱(chēng)模態(tài)(用下標(biāo)a表示),其各位移分量為

ua=0,va=0,wa=1

(15)

ut=0,vt=r,wt=0

(16)

變形圖如圖3(c)所示.至此,通過(guò)分析獲得了所有可能的屈曲模態(tài)和相應(yīng)的位移函數(shù)uk(θ)的表達(dá)式.因此,下面的研究目標(biāo)是位移振幅函數(shù)φk(x)的表達(dá)式.

圖3 變形圖

1.4 特征值函數(shù)

在廣義梁理論中,通過(guò)特征值求臨界屈曲應(yīng)力[10]有Galerkin法和有限元法.本文為了求解薄壁圓柱殼的能量變分方程,采取Galerkin方法[11],由邊界條件決定位移振幅函數(shù)φk(x)的表達(dá)式,從而獲得能量變分方程的特征值函數(shù),進(jìn)而求得臨界屈曲應(yīng)力.

1) 邊界條件為簡(jiǎn)支、特征值函數(shù)為三角函數(shù)形式.其位移振幅函數(shù)可以設(shè)成

(17)

式中,dk為第k階模態(tài)的振幅(可由具體問(wèn)題給出);n為縱向半波數(shù);L為結(jié)構(gòu)的軸向長(zhǎng)度.

2) 邊界條件為非簡(jiǎn)支.位移振幅函數(shù)可以設(shè)成[12]

(18)

本文主要研究均勻壓縮下的簡(jiǎn)支薄壁圓柱殼的穩(wěn)定性問(wèn)題,因此,變分方程(9)變形為

(19)

整理后,特征值方程為

式中,λb為分叉屈曲下的特征值,并且只與軸向坐標(biāo)L有關(guān).將位移振幅函數(shù)代入式(20),分叉屈曲特征值λb為

(21)

(22)

(23)

(24)

(25)

2 薄壁圓柱殼穩(wěn)定性數(shù)值仿真

通過(guò)理論推導(dǎo)已得到臨界應(yīng)力的表達(dá)式,因此,為了驗(yàn)證此方法的適用性和正確性,下文將針對(duì)具體算例,采用理論方法與數(shù)值解法進(jìn)行對(duì)比,研究軸長(zhǎng)和厚度對(duì)屈曲應(yīng)力的影響.

2.1 構(gòu)件尺寸比

當(dāng)構(gòu)件半徑和厚度為定量,隨著軸長(zhǎng)的不斷增長(zhǎng),構(gòu)件可劃分為極短柱殼、中短柱殼和細(xì)長(zhǎng)柱殼3類(lèi).本文以均勻壓縮圓柱殼為研究對(duì)象,已知圓柱殼中心截面半徑r=100 mm,厚度t=1 mm.采用鋁合金材料,彈性模量E=68 GPa,泊松比ν=0.33,柱頂受到1 600 N的均勻軸向荷載,兩端簡(jiǎn)支.

由文獻(xiàn)[14]可知,當(dāng)L≤17.37 mm時(shí),構(gòu)件屬于極短柱殼;當(dāng)L>17.37 mm時(shí),構(gòu)件屬于中短柱殼；當(dāng)L>1 004.30 mm時(shí),構(gòu)件屬于細(xì)長(zhǎng)柱殼.

2.2 圓柱殼軸長(zhǎng)與屈曲應(yīng)力的關(guān)系

基于以上理論分析,簡(jiǎn)支邊界條件下的正弦振幅函數(shù)為φk(x),分別令縱向半屈曲波數(shù)n=1,2,得到屈曲應(yīng)力σb和軸長(zhǎng)L的關(guān)系如圖4所示.

分析圖4中n=1的情況可知:

1) 當(dāng)軸長(zhǎng)L<17.37 mm時(shí),主要的屈服模態(tài)以軸對(duì)稱(chēng)模態(tài)的形式出現(xiàn),此時(shí)只有徑向位移,屈曲模態(tài)如圖3(b)所示.由構(gòu)件細(xì)長(zhǎng)比可知,構(gòu)件為極短圓柱殼.極短圓柱殼下臨界屈曲應(yīng)力值和對(duì)應(yīng)的臨界軸長(zhǎng)為[15-16]

圖4 屈曲應(yīng)力σb和軸長(zhǎng)L關(guān)系圖

2) 當(dāng)軸長(zhǎng)17.37

3) 當(dāng)軸長(zhǎng)L>1 000 mm,構(gòu)件為細(xì)長(zhǎng)圓柱殼,屈曲應(yīng)力σb隨著軸長(zhǎng)L的增長(zhǎng)不斷減小,主要以撓曲模態(tài)(m=1)的形式發(fā)生屈曲.

結(jié)合n=1和n=2的情況可得:

1) 由圖4可知當(dāng)縱向半屈曲波數(shù)n=2和n=1,屈曲應(yīng)力σb隨著軸長(zhǎng)L的變化趨勢(shì)完全一致,因此,n=2與n=1相比,相當(dāng)于整體圖形在橫坐標(biāo)方向平移一段距離.

(2) 當(dāng)L<100 mm,屈曲應(yīng)力存在2個(gè)相等極小值點(diǎn),分別對(duì)應(yīng)于軸對(duì)稱(chēng)模態(tài);當(dāng)1001 600 mm時(shí),構(gòu)件只以撓曲模態(tài)(m=1)的形式發(fā)生屈曲.

2.3 壁厚與屈曲應(yīng)力的關(guān)系

在2.2節(jié)中已研究了由圓柱殼軸長(zhǎng)L引起的屈曲應(yīng)力σb變化趨勢(shì),而影響圓柱殼穩(wěn)定性的因素往往還有構(gòu)件的厚度t.由于本文研究對(duì)象為薄壁圓柱殼,厚度t取為1～5 mm,其他參數(shù)不變,邊界條件仍為兩端簡(jiǎn)支.此外,采用有限元軟件建模,得到相應(yīng)條件下圓柱殼的屈曲承載力值,與本文方法對(duì)比結(jié)果如表1所示.

表1 3種柱殼隨厚度變化的屈曲承載力 N/mm

由表1可知,隨著厚度的增加,屈曲承載力整體呈增大趨勢(shì).在一定范圍內(nèi),增大厚度能夠提高結(jié)構(gòu)的穩(wěn)定性.通過(guò)本文推導(dǎo)方法所得的理論解和有限元計(jì)算結(jié)果進(jìn)行對(duì)比可知,極短和中短圓柱殼屈曲承載力誤差在2%～3%之間,細(xì)長(zhǎng)柱殼的誤差稍微偏大,但也小于5%,其主要原因可能是由于細(xì)長(zhǎng)柱殼屈曲受圓周波數(shù)的影響.極短柱殼屈曲承載力增長(zhǎng)幅度為三者中最大,此時(shí)厚度對(duì)于提高結(jié)構(gòu)穩(wěn)定性的作用最為明顯.對(duì)比于細(xì)長(zhǎng)柱,厚度對(duì)細(xì)長(zhǎng)柱的屈曲承載力的影響相對(duì)較小.

3 結(jié)語(yǔ)

基于廣義梁理論研究了薄壁圓柱殼穩(wěn)定性分析方法,并采用該方法研究了軸長(zhǎng)和厚度對(duì)屈曲應(yīng)力的影響.在數(shù)值仿真中,極短圓柱殼主要以軸對(duì)稱(chēng)模態(tài)的形式發(fā)生屈服,存在極小值.中短圓柱殼的屈曲主要發(fā)生在圓周波數(shù)m在2～10區(qū)段內(nèi),屈曲應(yīng)力σb隨著軸長(zhǎng)L的增長(zhǎng)整體呈現(xiàn)下降的趨勢(shì),其中也存在幾個(gè)局部極小的屈曲應(yīng)力.細(xì)長(zhǎng)圓柱殼主要以撓曲模態(tài)(m=1)的形式發(fā)生屈曲,屈曲應(yīng)力σb隨著軸長(zhǎng)L的增長(zhǎng)不斷減小.另外,在厚度變化的情況下,厚度對(duì)極短圓柱殼的影響為三者中最大,其次為中短圓柱殼,最后為細(xì)長(zhǎng)圓柱殼.但總體來(lái)說(shuō),一定范圍內(nèi),隨著厚度的增加,屈曲承載力都相應(yīng)增加.因此,在一定范圍內(nèi),通過(guò)增加厚度提高構(gòu)件的屈曲承載能力也是顯而易見(jiàn)的.通過(guò)與已有文獻(xiàn)結(jié)論、有限元分析的結(jié)果對(duì)比,證明了將廣義梁理論推廣到圓形截面進(jìn)行屈曲分析是可行的,而且計(jì)算精度較高,與已有結(jié)果保持較好的一致性,這對(duì)薄壁圓柱殼穩(wěn)定性分析提供了一個(gè)新的方法途徑.

)

[1] Schardt R.Verallgemeinertetechnischebiegetheorie[M]. Berlin: Springer, 1989.

[2] Davies J M, Leach P, Heinz D. Second-order generalized beam theory [J].JournalofConstructionalSteelResearch, 1994,31(2/3): 221-241.

[3] Silvestre N. Generalised beam theory to analyze the buckling behavior of circular cylindrical shells and tubes [J].Thin-WalledStructures, 2007,45(2): 185-198.

[4] Ranzi G, Luogo A. A new approach for thin-walled member analysis in the framework of GBT [J].Thin-WalledStructures, 2011,49(11): 1404-1414.

[5] Dinis P B, Camotim D, Silvestre N. GBT formulation to analyse the buckling behaviour of thin-walled members with arbitrarily ‘branched’ open cross-sections [J].Thin-WalledStructures, 2006,44(1): 20-38.

[6] Basaglia C, Camotim D, Silvestre N. Global buckling analysis of plane and space thin-walled frames in the context of GBT [J].Thin-WalledStructures, 2008,46(1): 79-101.

[7] Basaglia C, Camotim D, Silvestre N. Post-buckling analysis of thin-walled steel frames using generalised beam theory (GBT) [J].Thin-WalledStructures, 2013,62: 229-242.

[8] 李開(kāi)禧. 彈性薄壁桿件的翹曲[M]. 北京:中國(guó)建筑工業(yè)出版社,1990.

[9] Ju G T, Kyriakides S. Bifurcation and localization instabilities in cylindrical shells under bending-Ⅱ： predictions [J].InternationalJournalofSolidsandStructures, 1992,29(9): 1143-1171.

[10] Camotim D, Silvestre N, Goncalves R,et al. GBT-based structural analysis of thin-walled members: overview, recent progress and future developments [J].EngineeringStructures,MechanicsandConstruction,SMCD, 2006,14-17: 187-204.

[11] Silvestre N, Camotim D. GBT buckling analysis of pultruded FRP lipped channel members [J].ComputerandStructures, 2003,81(18/19): 1889-1904.

[12] Silvestre N, Camotim D. Distortional buckling formulae for cold formed steel C and Z-section members: part Ⅰ-derivation [J].Thin-WalledStructures, 2004,42(11): 1567-1597.

[13] Singer J, Arbocz J, Weller T.Bucklingexperiments-experimentalmethodsinbucklingofthin-walledstructures[M]. New York, USA: Wiley, 2002.

[14] 陳興華,龍連春.軸壓薄壁圓柱殼彈性失穩(wěn)承載力分類(lèi)及對(duì)比[C]//北京力學(xué)會(huì)第15屆學(xué)術(shù)年會(huì)論文摘要集.北京,2009:215-216.

Chen Xinghua, Long Lianchun. Thin-walled cylindrical shells’ elastic buckling bearing capacity of classification and comparison with axial compressive [C]//MechanicsofBeijingtheFifteenthAnnualMeetingofAbstracts. Beijing, 2009:215-216.(in Chinese)

[15] Timoshenko S P, Gere J M.Theoryofelasticstability[M]. New York, USA: McGraw-Hill, 1961.

[16] Brush D O, Almroth B O.Bucklingofbars,platesandshells[M]. New York, USA: McGraw-Hill, 1975.