黃 思，易天坤，歐晨希，林冠堂

(1.華南理工大學(xué) 機械與汽車工程學(xué)院， 廣州 510641； 2.廣東省特種設(shè)備檢測研究院珠海檢測院， 廣東 珠海 519002)

0 引言

球罐被廣泛應(yīng)用于石化行業(yè)，主要用于儲存液化石油氣[1]。當(dāng)施加到結(jié)構(gòu)上的載荷達(dá)到某一值時，若繼續(xù)增加微小的載荷，結(jié)構(gòu)的位移會發(fā)生很大的改變，稱結(jié)構(gòu)發(fā)生了屈曲或者失穩(wěn)。由于受到安裝施工和土壤等環(huán)境的影響[2]，球罐支柱將產(chǎn)生沉降和偏移。球罐支柱的沉降和偏移一旦達(dá)到一定程度時，球罐可能因局部受力不均而發(fā)生屈曲變形和破壞。

國內(nèi)外學(xué)者相繼開展了圓柱罐和球罐在沉降下的靜力學(xué)和屈曲研究。靜力學(xué)研究方面，周蓓蓓等[3]采用有限元法研究了支柱沉降形式對球罐總體應(yīng)力分布的影響；Yang等[4]采用有限元法計算了支柱不均勻沉降下的球罐應(yīng)力和支柱位移；Gao等[5-6]研究了球罐支柱沉降對球殼和支柱應(yīng)力的影響；房江祥[7]運用ABAQUS軟件模擬了1.5 m3球罐的支柱沉降對其拉桿應(yīng)力和極限載荷的影響；Chen等[8]通過ANSYS軟件建立了圓柱罐的有限元模型，研究了圓柱罐在不均勻沉降下的變形行為。屈曲研究方面，Jonaidi等[9-10]采用有限元與實驗驗證相結(jié)合的方法，研究了圓柱罐在諧波沉降下的屈曲響應(yīng)；Godoy等[11-12]通過實驗方法研究了局部沉降下圓柱罐的屈曲響應(yīng)；范海貴[13]采用現(xiàn)場實測、理論推導(dǎo)和有限元數(shù)值模擬相結(jié)合的方法，研究了圓柱罐在沉降和軸向載荷聯(lián)合作用下的屈曲行為；趙陽等[14]分別研究了整體沉降和局部沉降下圓柱罐的屈曲變形行為。

綜上所述，現(xiàn)有研究主要是針對圓柱罐和球罐在沉降下的靜力學(xué)和屈曲研究，還未對球罐支柱在偏移下的靜力學(xué)和屈曲進(jìn)行過研究。因此，根據(jù)現(xiàn)場檢測的球罐支柱偏移數(shù)據(jù)，開展球罐在支柱偏移下的屈曲研究，求解計算球罐支柱偏移的極限值。該極限值可為現(xiàn)場監(jiān)測數(shù)據(jù)提供理論參考和作為實時監(jiān)測的警戒值，為預(yù)防球罐在支柱偏移下發(fā)生屈曲失效而引發(fā)重大安全事故具有重要意義。

1 計算模型

圖1是根據(jù)表1[15]球罐的設(shè)計參數(shù)所建立的球罐三維模型，其材料特性參數(shù)見表2。

表1 丙烷球罐設(shè)計參數(shù)

圖1 球罐三維模型示意圖

表2 球罐主要結(jié)構(gòu)的材料特性參數(shù)

2 屈曲響應(yīng)分析

2.1 網(wǎng)格劃分

運用Workbench平臺中Mesh對圖1的計算模型進(jìn)行網(wǎng)格劃分，為避免網(wǎng)格數(shù)對計算結(jié)果的影響，需考察計算域在不同網(wǎng)格數(shù)下的計算結(jié)果?？紤]到本文僅對球罐進(jìn)行屈曲變形的研究，圖2給出了通過Workbench平臺計算得到的球罐的最大總位移xmax隨網(wǎng)格數(shù)的變化規(guī)律。

從圖2可知，當(dāng)球罐的網(wǎng)格數(shù)達(dá)到85萬之后，球罐最大總位移xmax值基本穩(wěn)定。考慮到計算精度和成本，選用988 554為球罐計算域的網(wǎng)格數(shù)，模型網(wǎng)格劃分結(jié)果見圖3。

圖2 網(wǎng)格無關(guān)性檢驗曲線

圖3 球罐網(wǎng)格劃分結(jié)果示意圖

2.2 載荷與約束條件

考慮球罐充裝率φ為90%的工況，球罐的載荷有重力G、內(nèi)壁壓力p。其中壓力p隨液下深度h的變化規(guī)律為：

(1)

式中：po為球罐的工作壓力，為1.42 MPa；ρl為液態(tài)丙烷的密度，為474.04 kg/m3；g為重力加速度，取9.81 m/s2；H為充液高度，其中H與φ有如下關(guān)系[16]：

(2)

式中：R為球罐的內(nèi)半徑。

圖4給出了12根支柱編號的示意圖。分別單獨研究單根支柱沿著周向偏移和徑向偏移情況。周向偏移是指支柱底部沿著圖4以O(shè)為圓心的大圓圓周方向的偏移，定義順時針方向為正，逆時針方向為負(fù)。徑向偏移是指支柱底部沿著以O(shè)為圓心的大圓半徑方向的偏移，定義向內(nèi)偏移為正，向外偏移為負(fù)。根據(jù)現(xiàn)場檢測的支柱偏移數(shù)據(jù)，支柱的偏移量一般在30 mm以內(nèi)。根據(jù)周向偏移的對稱性，對于周向偏移僅考慮正偏移的情況。因此，設(shè)置的單根支柱的周向偏移量Wt和徑向偏移量Wr見表3、4。

圖4 支柱編號示意圖

表3 單根支柱周向偏移量

2.3 屈曲分析原理

按照臨界狀態(tài)的特性，屈曲可分為跳躍屈曲、極值屈曲和分支屈曲，分別對應(yīng)圖5中的(a)、(b)和(c)所示的載荷P—位移x曲線[17]。圖5中的載荷P為廣義載荷，可以是力載荷或是位移載荷等，本文研究的載荷為球罐支柱偏移的位移載荷。

圖5 不同屈曲特性的載荷P-位移x曲線

當(dāng)球罐支柱的偏移達(dá)到一定水平時，球罐在自重G、內(nèi)壁壓力p和支柱偏移Wr和Wt等作用下將發(fā)生整體或局部彎曲大變形。因此，需對球罐開展屈曲分析，確定其支柱偏移的極限值。針對球罐的屈曲分析可采用“弧長法”?！盎￠L法”最初是由Risks等提出[18]，后經(jīng)Roma等[19]逐步改進(jìn)而完善?；￠L法是目前結(jié)構(gòu)屈曲分析中效率較高的一種方法。

弧長法中非線性平衡方程組一般可以表示為：

K·x=F

(3)

式中：矩陣K為切向剛度矩陣，向量x為節(jié)點位移，向量F為外載荷?，F(xiàn)將方程寫成增量形式，并引入載荷因子λ，則平衡方程：

K·Δx=ΔλF

(4)

此時平衡方程共有n+1個未知量：Δx1,Δx2,…,Δxn,λ，而方程數(shù)只有n個，因此需再補充一約束條件：

Δλ2FTF+ΔxTΔx=Δl2

(5)

弧長法的求解流程如圖6(a)所示，通過該流程的求解得到完整的λ-xmax曲線。如圖6(b)所示，弧長法是利用平衡方程和弧長約束條件尋找平衡點的過程。由圖6(b)、表3和表4中的球罐支柱的偏移是逐步加載的。當(dāng)支柱加載的偏移量超過C點所對應(yīng)的極限值時，球罐某些部位將會突然喪失原有的幾何形狀，發(fā)生大幅的屈曲變形，即C點后球罐發(fā)生了屈曲?；诨￠L法的數(shù)值模擬方法，首先分別單獨研究球罐在單根支柱底部周向偏移和徑向偏移下的屈曲響應(yīng)過程，然后根據(jù)現(xiàn)場檢測的支柱偏移數(shù)據(jù)，開展球罐在支柱偏移下的屈曲研究。

表4 單根支柱徑向偏移量

圖6 弧長法迭代求解過程

3 支柱偏移的屈曲分析

在應(yīng)用Workbench平臺對球罐進(jìn)行屈曲響應(yīng)分析時，分別單獨將表3和表4中球罐支柱的偏移逐步加載。載荷因子λ表示球罐屈曲響應(yīng)過程中支柱偏移的加載比例，例如λ=0.5時，球罐支柱加載的偏移分別為0.5Wr和0.5Wt。當(dāng)λ從0增大到1時，意味著球罐支柱的偏移全部加載完成。

圖7是將表3中的數(shù)據(jù)逐步加載并求解得到的單根支柱底部周向偏移的λ-xmax曲線，其中xmax為最大總位移，λ為載荷因子。圖8為該狀態(tài)下的曲線中A、B、C和D四點處的變形云圖，其中A、B、C和D四點對應(yīng)的λ值分別為0.05、0.15、0.25和0.55。

圖7 單根支柱底部周向偏移下的λ-xmax曲線

圖8 單根支柱底部周向偏移下球罐屈曲響應(yīng)的變形云圖(變形縮放系數(shù)=20)

對照圖7與圖5可知，該球罐在支柱底部周向偏移下的屈曲響應(yīng)符合分支屈曲的特點，在C點前球罐的變形隨支柱底部周向偏移的增大呈線性變化，說明球罐處于平衡狀態(tài)，還未發(fā)生屈曲。C點時球罐的變形云圖如圖8(c)所示，對應(yīng)的最大總位移為7.5 mm(在支柱1的底部)。當(dāng)支柱底部周向偏移超過C點對應(yīng)的極限值(λ>0.25，即Wt>7.5 mm)，球罐局部位移突變發(fā)生了屈曲，如圖8(d)所示，此時最大總位移為102.8 mm，屈曲的位置在支柱1與支柱2之間的拉桿處。

圖9是將表4中的數(shù)據(jù)逐步加載并求解得到的在單根支柱底部徑向偏移下的λ-xmax曲線。圖10為單根支柱底部徑向內(nèi)偏移的曲線中A、B、C和D四點處的變形云圖，其中A、B、C和D四點對應(yīng)的λ值分別為0.25、0.5、0.75和0.9。

圖9 單根支柱底部徑向偏移下的λ-xmax曲線

圖10 單根支柱底部徑向內(nèi)偏移下球罐屈曲響應(yīng)的變形云圖(變形縮放系數(shù)=20)

從圖9可知，支柱底部徑向外偏移時未出現(xiàn)位移突變，表明支柱底部徑向外偏移在現(xiàn)有數(shù)據(jù)范圍內(nèi)還未使球罐發(fā)生屈曲。將圖9與圖5進(jìn)行對照可知，該球罐在支柱底部徑向內(nèi)偏移下的屈曲響應(yīng)符合分支屈曲的特點。在C點前球罐的變形隨支柱徑向偏移的增大呈線性變化，這時球罐處于平衡狀態(tài)未發(fā)生屈曲。C點時球罐的變形云圖如圖10(c)所示，對應(yīng)的最大總位移為23.1 mm。當(dāng)支柱徑向內(nèi)偏移超過C點時對應(yīng)的極限值(λ>0.75，即Wr>22.5 mm)，球罐局部發(fā)生了屈曲，如圖10(d)所示，此時最大總位移為39.4 mm，屈曲的位置在支柱1與支柱12之間的拉桿處。

4 檢測支柱偏移的屈曲分析

全站儀具有高精度、高穩(wěn)定性和高可靠性技術(shù)等優(yōu)點，極其適合用于球罐支柱偏移位移的測量。對現(xiàn)場球罐支柱的偏移數(shù)據(jù)的測量工作采用徠卡TCR1202+R1000電子全站儀，其測量的精度為1 mm。表5是用該儀器在現(xiàn)場測量的支柱偏移數(shù)據(jù)。

表5 現(xiàn)場檢測的球罐支柱偏移值 mm

為了模擬球罐在檢測支柱偏移下的屈曲響應(yīng)過程，將表5中的球罐支柱的偏移逐步加載，在Workbench平臺運用弧長法進(jìn)行求解計算。圖11是通過求解得到的球罐在檢測支柱偏移下的λ-xmax曲線。圖12為該情形下的曲線中A、B、C和D四點處的變形云圖，其中A、B、C和D四點對應(yīng)的λ值分別為0.1、0.3、0.5和0.7。

圖11 檢測支柱偏移下的λ-xmax曲線

圖12 檢測支柱偏移下球罐屈曲響應(yīng)的變形云圖(變形縮放系數(shù)=20)

由圖11可知，球罐的變形隨支柱偏移的增大呈線性變化，直到支柱的偏移加載完成都未出現(xiàn)明顯的位移突變，說明球罐在該情形下并未發(fā)生屈曲。分析其原因，發(fā)現(xiàn)所測球罐支柱底部的徑向偏移數(shù)據(jù)均是30 mm以內(nèi)的外偏移，由數(shù)值求解結(jié)果可知球罐在30 mm內(nèi)的徑向外偏移不發(fā)生屈曲。同時球罐支柱底部的周向偏移僅有支柱4、7和8超過了數(shù)值計算結(jié)果的極限值7.5 mm，加之支柱的相互約束使支柱之間的周向偏移實際未超過極限值。

5 結(jié)論

1) 在所研究的偏移量范圍內(nèi)支柱底部徑向外偏移未出現(xiàn)屈曲。支柱底部徑向內(nèi)偏移的極限值比底部周向偏移的極限值大，支柱底部徑向內(nèi)偏移的屈曲變形量比底部周向偏移的變形量小，說明球罐支柱底部在周向偏移時更容易屈曲，在支柱底部徑向外偏移最不易屈曲。

2) 球罐的屈曲變形最大的位置位于拉桿處，屈曲變形較大的位置是偏移的支柱與球殼的連接處。該球罐在支柱底部徑向內(nèi)偏移和周向偏移下的屈曲響應(yīng)均符合分支屈曲特點。

3) 根據(jù)現(xiàn)場檢測的球罐支柱的偏移數(shù)據(jù)，求解發(fā)現(xiàn)該球罐未發(fā)生屈曲變形，仍屬于小幅彈性變形。