本站小編為你精心準備了肋環形加勁淺球殼屈曲及穩定承載力分析參考范文，愿這些范文能點燃您思維的火花，激發您的寫作靈感。歡迎深入閱讀并收藏。

《安徽工業大學學報》2016年第2期

摘要：

選取36種結構形式的基本模型，采用弧長法對以第一階屈曲模態作為初始形態的結構進行非線性屈曲分析，探討邊界條件、初始缺陷、矢跨比和截面尺寸對結構屈曲性能的影響。通過分析曲率半徑、等效薄膜厚度和等效彎曲厚度與臨界屈曲荷載之間的量化關系，采用最小二乘法擬合肋環型加勁淺球殼臨界失穩式，分析并選取安全系數，得到穩定容許承載力標準值計算式。

關鍵詞：

氣柜頂蓋；肋環形加勁淺球殼；臨界屈曲荷載；有效厚度；穩定承載力

大容量的鋼結構儲罐隨鋼鐵廠副產二次能源合理利用的要求應運而生。然而，在實際應用中，大型儲罐的頂蓋穩定性理論研究欠缺，指導設計的相關規范較少，致使在已建成的儲罐中，時有頂蓋失穩的狀況發生[1-2]。環肋形加勁淺球殼是各種儲罐拱頂常采用的形式之一，因其用鋼量較省、減輕了頂蓋的荷載而倍受青睞。此類結構的整體穩定性問題仍是研究熱點。采用有限元模型的簡化和計算模式使問題的研究有了一定進展。如加大主梁在鋪板面內的抗彎模量建立簡化模型，以一階極限荷載值作為結構發生失穩的臨界荷載值[3]；將40倍到60倍鋪板厚度的板寬計入次肋剛度的計算屈曲荷載的模型[4]；考慮鋪板局部屈曲后強度，利用幾何非線性弧長法進行全過程分析，研究結構失穩的原因和屈曲過程[5]。本課題以華東某10萬m3一段橡膠膜密封氣柜的頂蓋為背景，建立36種結構模型，采用有限元模擬與理論分析相結合的方法，研究矢跨比、邊界條件、初始缺陷等因素對加勁殼的屈曲性能的影響，確定有效薄膜厚度和彎曲厚度的取值，得出相應設計參數，擬合穩定容許承載力計算式。

1有限元模型的選取與計算

1.1頂蓋結構

華東某10萬m3肋環型氣柜的頂蓋由中心環梁、邊環梁、經向主肋、環向主肋、經向次肋、頂板組成的淺球形加勁球殼(圖1)。經向主肋的數量為32根，截面為槽鋼[18a；環向主肋為9道，沿經向圓弧等距離分布，外側4道的截面為槽鋼[18a，內側5道的截面為槽鋼[14a。中心環梁的截面為槽鋼[18a。經向主肋與環向主肋構成的骨架構造單元中設有經向次肋，截面為L75×6角鋼，數量由外環帶向內環帶逐步減少；蓋板采用4mm鋼板。其他模型幾何參數如下：球殼曲率半徑為80810mm，矢高為6210mm，跨度為62132mm，矢跨比為1/10，中心角為22.60°，中心環梁半徑為1200mm。

1.2有限元模型的選取

根據工程實例，本研究模型保持結構的跨度、主次肋的分布形式和數目不變，改變模型的矢跨比和構件的截面尺寸，得到以矢跨比分組的36種結構形式的基本模型，見表1。其中A，B，C，D組分別對應矢跨比為1/9，1/10，1/11，1/12，即曲率半徑為73384，80810，88230，95975mm的模型。經向次肋均為角鋼L75×6。

1.3有限元模型的臨界屈曲荷載計算

在對上述基本模型進行有限元分析的過程中，每個模型的邊界條件分別取為固接和鉸接，結構的初始缺陷分別取L/200，L/300,L/500,L/800,L/1000。參照一致缺陷模態法[6-8]，先進行特征值屈曲分析，得到第一階屈曲模態，以第一階屈曲模態作為結構的初始形態，定義初始缺陷值；然后采用弧長法對該結構進行非線性屈曲分析，得到頂蓋結構的整體失穩時臨界屈曲荷載的計算結果。選取經向主肋的豎向位移最大的節點進行研究，以1/300跨度為初始缺陷值、邊界條件為固接和鉸接所計算的屈曲臨界荷載值分別為qcr,F和qcr,P。將36種模型的屈曲臨界荷載值列入表1。其中，模型A1初始缺陷為L/300的固接和鉸接的屈曲變形圖見圖2。

2肋環型加勁淺球殼的屈曲影響因素分析

2.1邊界條件

氣柜頂蓋邊緣設有T型截面的邊環梁，經向主肋、次肋和板與邊環梁可靠連接。在肋環型加勁淺球殼模型的分析過程中，采用固接和鉸接兩種邊界條件來計算頂蓋屈曲臨界荷載，對不同矢跨比結構的兩種計算結果進行比較，見圖3。通過對比可知，固接時的屈曲荷載與鉸接時的屈曲荷載最大平均偏差率約為7.1%，因此固接和鉸接時結構的屈曲荷載值基本持平。分析其原因在于：大跨度的殼狀結構，支座之間的距離較大，臨界屈曲波長較小，它仍以特定的臨界屈曲波長發生失穩[9]，臨界屈曲荷載與塌陷區域處的截面屬性直接相關，與遠處的支座形式間接相關。因此，固接和鉸接形式的支座對屈曲臨界荷載的影響較小。

2.2初始缺陷

鑒于邊界條件對屈曲荷載影響較小，故選取固接情況下，初始缺陷值分別為結構跨度L的1/1000，1/800，1/500，1/300和1/200等5種的計算結果，繪制各組結構的屈曲臨界荷載隨初始缺陷變化的關系見圖4。由圖4可知：隨著結構初始缺陷的增大，結構的屈曲臨界荷載呈現出減小的趨勢；但是總體上，初始缺陷的變化對結構的屈曲臨界荷載值的影響較小；當初始缺陷超過L/300后有一定影響，屈曲荷載的偏差率增大。因此，L/300的初始缺陷可作為結構安裝誤差的限定值，以控制初始缺陷對結構屈曲荷載的影響。

2.3矢跨比和截面尺寸

基于上述分析結論，邊界條件和初始缺陷對模型整體失穩時的屈曲荷載值影響較小，按矢跨比分類，選取在邊界條件為固結、初始缺陷為L/300時的模型非線性屈曲計算結果，繪制各組結構屈曲節點的荷載-位移變化曲線圖5。由圖5可知：結構失穩時，隨矢跨比的降低，最大變形值增大，最大臨界屈曲荷載值降低；結構的臨界荷載值隨主肋截面尺寸的增大而增大；在曲線的直線上升段，結構的位移值很小，表明結構處于彈性工作階段，其剛度較大；達到極值點之前，結構的位移增值比荷載的增值漸漸加大，表明結構處于彈塑性工作階段，雖然剛度逐漸退化，但仍處于穩定平衡狀態；曲線下降段，維持平衡的荷載減小，結構變形加快，表明結構的剛度矩陣為非正定，結構的剛度退化明顯，呈現脆性破壞特征，處于不穩定平衡狀態，最終結構出現較大范圍塌陷，發生整體失穩。因此，頂蓋結構的失穩模態為極值點失穩，該極限點對應的荷載可視為結構整體失穩的臨界屈曲荷載值，即穩定極限承載力。

3穩定容許承載力的計算公式擬合

對于各向同性的球面殼狀結構，引入分離剛度的概念，以等效彎曲厚度tB和等效薄膜厚度tm分別用于描述構件的慣性矩和面積對結構穩定承載力的貢獻，彈性臨界法向荷載表示為[10]qcr=CEæèçöø÷tmR2æèçöø÷tBtm32（1）式中：C為系數；E為鋼材的彈性模量；R為球殼的曲率半徑。

3.1肋環形加勁淺球殼的等效彎曲厚度和等效薄膜厚度

肋環形加勁淺球殼的每個計算單元在經向和環向的截面屬性不同。由圖6可以得出一個計算單元（陰影部分）經向和環向的組合截面屬性：lφ，lθ是經向和環向的長度；Aφ，Aθ是經向和環向的截面面積；Iφ，Iθ是經向和環向的截面慣性矩。由于肋環加勁球殼的薄膜厚度和彎曲厚度在經向和環向是不同的，借鑒文獻[11-12]的等效厚度擬合方式，取頂板、經向和環向的計算厚度在組合截面的權重均為1，得到相應的等效厚度：tm=31tφA+1t+1tθA(2)tB=t3φI+t3+t3θI33(3)式中：t為頂板厚度；tφA=Aφlφ，tθA=Aθlθ為經向、環向等效薄膜厚度；tφI=12Iφlφ，tθI=12Iθlθ為經向、環向等效彎曲厚度。

3.2穩定承載力公式的擬合

針對36種基本模型的計算結果，選取初始缺陷為L/300、邊界條件為固接的臨界屈曲荷載計算結果(見表1)，分析等效效薄膜厚度和等彎曲厚度與臨界屈曲荷載的量化關系，進而擬合肋環型加勁淺球殼整體穩定性的極限承載力計算公式。參照式(1)，可擬定肋環型加勁淺球殼的穩定極限承載力表達式為qcr=αER-2tβmtγB(4)式中：α為系數；β為等效薄膜厚度的指數系數；γ為等效彎曲厚度的指數系數。因為式(4)為非線性多元函數，根據表1中的有限元計算結果，qcr,F和該位移最大點所在計算單元的截面屬性R，tm和tB計算值，引入MATLAB中的基于Gauss-Newton算法的最小二乘擬合函數nlinfit函數，對式(4)進行擬合，得出未知系數α，β和γ，從而得到穩定極限承載力計算式為qcr=0.5748ER-2t0.5385mt1.4615B(5)結構的極限承載力是結構完全破壞前所能承受的外荷載的最大能力。具有一定安全度和不引起超出使用范圍的變形的外荷載是容許承載力，它兼顧強度和變形兩方面的要求。因此，以結構的極限荷載作為容許承載力是不安全的，安全系數K的選擇至關重要[13-14]。

針對本文所分析的肋環形加勁淺球殼結構的安全系數K值確定，應該考慮以下因素：以系數k1=1.4考慮外部作用(施工荷載和半跨荷載)的不確定性可能帶來的不利影響因素；以系數k2=1.4考慮結構抗力的不確定性(如部分構件的部分截面在結構達到臨界點前進入彈塑性狀態，呈現脆性破壞特征)可能帶來的不利影響因素；以系數k3=1.3來考慮結構工作條件中的其他不利因素(如低溫環境)可能帶來的不利影響因素；以系數k4=2來考慮結構穩定性分析中可能的不精確性(如網格尺寸偏大造成模型剛度增大、板和肋的截面剛度相差很大造成變形不協調)可能帶來的不利影響因素。因此，對于按彈塑性全過程分析求得的穩定極限承載力，其安全系數K應取為1.4×1.4×1.3×2≈5。將穩定極限承載力計算式(5)除以安全系數5，則得肋環形加勁淺球殼結構的穩定容許承載力標準值簡化算式qks=0.115ER-2t0.54mt1.46B(6)

4結論

1)邊界條件采用固接與鉸接對頂蓋結構整體失穩的臨界屈曲荷載值影響較小。

2)初始缺陷的變化對結構的屈曲臨界荷載值的有一定影響，但影響較小。L/300的初始缺陷可作為結構安裝誤差的限定值，以控制初始缺陷對結構屈曲荷載的影響。

3)矢跨比和截面尺寸對臨界屈曲荷載值影響明顯。臨界屈曲荷載值隨矢跨比的降低而降低，隨截面尺寸的增大而增大。

4)初步設計時，式(6)可以用來估算肋環型加勁淺球殼的容許承載力。

參考文獻：

[1]李靜波.尺寸超規范的帶肋拱頂油罐的設計[J].山東化工，2007,36(6):23-24.

[2]黃文霞.渣油罐拱頂失穩分析及修復探討[J].石油化工設備技術，2008,29(1):25-28.

[3]周偉.新型煤氣柜的強度和穩定性分析[D].重慶：重慶大學，2004.

[4]常安果，周圣凱，龔宗宜.環肋型球面殼煤氣柜頂蓋穩定分析[J].鋼結構，2012,7(27):1-10.

[5]黃斌，毛文筠.大型煤氣柜頂蓋的幾何非線性全過程分析[J].武漢理工大學學報，2006,28(4):63-66.

[6]陳昕，沈世釗.網殼結構的幾何非線性分析[J].土木工程學報，1990,23(3):47-57.

[7]鄧長根，徐忠根.缺陷結構穩定分析的一致缺陷模態近似法[J].鋼結構，2012(增刊):209-216.

[8]陳紹蕃.鋼結構穩定設計指南[M].3版.北京：中國建筑工業出版社，2013.

[11]中國石油天然氣集團公司.立式圓筒形鋼制焊接油罐設計規范：GB50341—2014[S].北京：中國計劃出版社，2014.

[12]中國建筑科學研究院.鋼筋混凝土薄殼結構設計規程：JGJ22—2012[S].北京：中國建筑工業出版社，2012.

[13]中國建筑科學研究院.網殼技術規程：JGJ61—2003[S].北京：中國建筑工業出版社，2003.

[14]中國建筑科學研究院.空間網格結構技術規程：JGJ7—2010[S].北京：中國建筑工業出版社，2010.

作者:賈冬云 王濤 楊軍 單位:安徽工業大學建筑工程學院