陳嚴(yán)飛 馬尚 董紹華 婁方宇 倪恒 趙永濤 孫秉才

1中國(guó)石油大學(xué)（北京），油氣管道輸送安全國(guó)家工程實(shí)驗(yàn)室/城市油氣輸配技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室

2大連理工大學(xué)，工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室

3中國(guó)石油集團(tuán)安全環(huán)保技術(shù)研究院

為滿足國(guó)家戰(zhàn)略儲(chǔ)備需要，我國(guó)持續(xù)建設(shè)大型國(guó)家原油儲(chǔ)備庫(kù)。目前，大型原油儲(chǔ)備基地大多建于沿海地帶，且大部分儲(chǔ)罐都是在填海造陸形成的軟土上建成，軟土土質(zhì)較差、強(qiáng)度低、壓縮性大，導(dǎo)致儲(chǔ)罐地基易發(fā)生沉降，對(duì)儲(chǔ)罐安全造成不利影響。此外，在沿海地帶，儲(chǔ)罐易受到強(qiáng)風(fēng)影響，承受較大的風(fēng)載荷，大型儲(chǔ)罐屬于薄殼結(jié)構(gòu)，對(duì)外界風(fēng)載荷作用較為敏感，容易出現(xiàn)屈曲失穩(wěn)現(xiàn)象，風(fēng)載荷略微增大也會(huì)使儲(chǔ)罐罐壁發(fā)生較大變形，嚴(yán)重時(shí)會(huì)引起罐壁撕裂，導(dǎo)致泄漏事故。因此，對(duì)于沿海地帶的大型儲(chǔ)罐的地基沉降問(wèn)題及風(fēng)致屈曲問(wèn)題需要重視。

目前，針對(duì)大型儲(chǔ)罐地基沉降和大型儲(chǔ)罐風(fēng)致屈曲的研究已取得了諸多成果。陳嚴(yán)飛等人考慮儲(chǔ)罐敏感參數(shù)提出了罐頂徑向位移預(yù)測(cè)方法[1]。MAHER、PURDY、HOLROYD 等學(xué)者通過(guò)理論研究、風(fēng)洞試驗(yàn)等方法，分別提出了風(fēng)載荷在儲(chǔ)罐罐壁上的分布規(guī)律[2-4]。郭煥良通過(guò)風(fēng)洞試驗(yàn)，研究了大型儲(chǔ)罐的風(fēng)載荷特征以及罐組的風(fēng)壓干擾效應(yīng)[5]。浙江大學(xué)趙陽(yáng)教授課題組的林寅通過(guò)設(shè)計(jì)實(shí)施風(fēng)洞試驗(yàn)，獲得了敞口儲(chǔ)罐、平頂儲(chǔ)罐和穹頂儲(chǔ)罐的風(fēng)載荷數(shù)據(jù)，并對(duì)儲(chǔ)罐罐壁平均風(fēng)載荷體型系數(shù)進(jìn)行了擬合，為工程設(shè)計(jì)提供了參考[6]。魏化中采用數(shù)值計(jì)算方法，對(duì)四種不同容積的大型儲(chǔ)罐進(jìn)行了風(fēng)致屈曲分析，并在此基礎(chǔ)上研究了抗風(fēng)結(jié)構(gòu)（抗風(fēng)圈、加強(qiáng)圈）對(duì)儲(chǔ)罐抗風(fēng)能力的影響，結(jié)果表明抗風(fēng)結(jié)構(gòu)的存在可顯著提高儲(chǔ)罐抵抗風(fēng)載荷的性能[7]。張煒等人采用數(shù)值計(jì)算方法，對(duì)104m3儲(chǔ)罐的風(fēng)致屈曲性能進(jìn)行了研究，結(jié)果表明儲(chǔ)罐對(duì)風(fēng)載荷作用十分敏感，在此基礎(chǔ)上還研究了包邊角鋼對(duì)儲(chǔ)罐抗風(fēng)能力的影響，結(jié)果表明有包邊角鋼的儲(chǔ)罐臨界屈曲載荷遠(yuǎn)大于無(wú)包邊角鋼的儲(chǔ)罐[8]。宮建國(guó)等人采用數(shù)值計(jì)算方法，針對(duì)10×104m3儲(chǔ)罐，分別研究其在靜態(tài)風(fēng)壓和動(dòng)態(tài)風(fēng)壓下的結(jié)構(gòu)響應(yīng)[9]。王喆采用數(shù)值計(jì)算方法，研究了無(wú)液狀態(tài)下不同高徑比、徑厚比的大型儲(chǔ)罐的抗風(fēng)性能，并與有液狀態(tài)下儲(chǔ)罐的抗風(fēng)性能進(jìn)行了對(duì)比，結(jié)果表明儲(chǔ)液的存在可以增強(qiáng)儲(chǔ)罐抵抗風(fēng)載荷的能力，此外還研究了地基傾斜沉降條件下儲(chǔ)罐的抗風(fēng)性能，這也是少數(shù)考慮地基沉降與風(fēng)載荷耦合的研究之一[10]。

學(xué)者們?cè)谘芯繒r(shí)很少有考慮風(fēng)載荷和沉降載荷兩種因素耦合的情況，即地基沉降發(fā)生后大型儲(chǔ)罐的風(fēng)致屈曲問(wèn)題。實(shí)際上，建于沿海地帶的大型儲(chǔ)罐自投入運(yùn)行起便會(huì)開(kāi)始發(fā)生地基沉降，且隨著時(shí)間推移，地基沉降不斷增大，儲(chǔ)罐罐壁會(huì)產(chǎn)生不同程度的變形，若此時(shí)罐壁再承受較大的風(fēng)載荷，便有可能出現(xiàn)屈曲失穩(wěn)現(xiàn)象，影響儲(chǔ)罐的安全運(yùn)行。因此，對(duì)儲(chǔ)罐運(yùn)行時(shí)承受的風(fēng)載荷進(jìn)行計(jì)算，建立了儲(chǔ)罐有限元模型，對(duì)未發(fā)生地基沉降的儲(chǔ)罐進(jìn)行風(fēng)致屈曲計(jì)算，得到風(fēng)壓作用下儲(chǔ)罐的臨界屈曲載荷；然后對(duì)實(shí)測(cè)地基沉降進(jìn)行了傅里葉分解，獲得罐周沉降；最后改變儲(chǔ)罐的風(fēng)壓駐點(diǎn)周向角，研究了各因素對(duì)儲(chǔ)罐風(fēng)致屈曲的影響，可為地基不均勻沉降下儲(chǔ)罐風(fēng)致屈曲研究提供參考。

1 風(fēng)載荷

建于沿海地區(qū)的大型儲(chǔ)罐會(huì)受到風(fēng)載荷作用。實(shí)際作用于儲(chǔ)罐罐壁表面的風(fēng)壓在儲(chǔ)罐高度及罐周方向上都按一定的規(guī)律分布。大量風(fēng)洞實(shí)驗(yàn)的結(jié)果表明，風(fēng)壓在儲(chǔ)罐高度方向上的分布變化不大[11]，而在罐周方向上的風(fēng)壓分布可用傅里葉余弦函數(shù)表示，其形式為：

式中：p(θ)為周向角θ處的風(fēng)壓，Pa；p0為基本風(fēng)壓，Pa；i為傅里葉余弦函數(shù)的階數(shù)；a0、ai為各階傅里葉余弦函數(shù)的幅值，也稱為風(fēng)壓系數(shù)，屬于經(jīng)驗(yàn)值，一般由試驗(yàn)或觀測(cè)得到；n為傅里葉余弦函數(shù)的最高階數(shù)。

RISH、GORENC、PIRCHER 等通過(guò)研究，分別對(duì)風(fēng)壓系數(shù)進(jìn)行了定義[12-14]，國(guó)外部分標(biāo)準(zhǔn)規(guī)范也對(duì)風(fēng)壓系數(shù)進(jìn)行了定義，如EN 1993-4-1 等，具體的風(fēng)壓系數(shù)如表1 所示。按照表1 中各風(fēng)壓系數(shù)計(jì)算得出風(fēng)壓在罐周方向上的分布情況（圖1）。

表1 傅里葉風(fēng)壓系數(shù)Tab.1 Fourier coefficient of wind pressure

圖1 風(fēng)壓在罐周方向上的分布Fig.1 Circumferential distribution of wind pressure on the storage tank

由圖1 可知，由各種風(fēng)壓系數(shù)確定的風(fēng)壓分布都是關(guān)于周向角0°對(duì)稱的，且在周向角0°處即風(fēng)壓駐點(diǎn)處取得最大值，駐點(diǎn)附近的儲(chǔ)罐罐壁屬于迎風(fēng)面；周向角-30°至+30°范圍內(nèi)的儲(chǔ)罐罐壁承受風(fēng)的正壓，該范圍內(nèi)的儲(chǔ)罐罐壁易發(fā)生風(fēng)致屈曲，而周向角±30°至±180°范圍內(nèi)的儲(chǔ)罐罐壁承受風(fēng)的負(fù)壓，此范圍內(nèi)的儲(chǔ)罐罐壁不易發(fā)生風(fēng)致屈曲；負(fù)風(fēng)壓約在周向角±60°至±90°范圍內(nèi)取得最大值。

對(duì)比各組風(fēng)壓系數(shù)可知，只有EN 1993-4-1 考慮了儲(chǔ)罐高度與直徑的影響，其他各組風(fēng)壓系數(shù)均未考慮儲(chǔ)罐結(jié)構(gòu)特征對(duì)罐周方向上風(fēng)壓分布的影響，因此在后續(xù)研究中，選擇EN 1993-4-1 定義的風(fēng)壓系數(shù)計(jì)算罐周方向上的風(fēng)壓分布。

2 風(fēng)載荷單獨(dú)作用下的儲(chǔ)罐屈曲計(jì)算

2.1 儲(chǔ)罐有限元模型

采用我國(guó)某10×104m3大型外浮頂儲(chǔ)罐的結(jié)構(gòu)尺寸數(shù)據(jù)建立儲(chǔ)罐有限元模型（圖2）。所選儲(chǔ)罐直徑為80 m，高度為21.7 m，底板直徑為80.4 m，厚度為20 mm，儲(chǔ)罐所用板材鋼的彈性模量為2.06×1011Pa，泊松比為0.3。在儲(chǔ)罐罐壁頂部設(shè)有包邊角鋼，尺寸為100 mm×12 mm。有限元模型的邊界條件為：儲(chǔ)罐為空罐；約束儲(chǔ)罐底板的徑向位移及切向位移，防止儲(chǔ)罐發(fā)生剛體運(yùn)動(dòng)。

圖2 儲(chǔ)罐有限元模型Fig.2 Finite element model of the storage tank

SCHMIDT 等人曾在其研究中搭建過(guò)小尺寸儲(chǔ)罐風(fēng)致屈曲實(shí)驗(yàn)?zāi)Ｐ?，并獲得了相關(guān)實(shí)驗(yàn)數(shù)據(jù)[15]。采用SCHMIDT 等人實(shí)驗(yàn)中的儲(chǔ)罐結(jié)構(gòu)參數(shù)及邊界條件建立有限元模型，將有限元計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比（表2）。由表2 可知，有限元計(jì)算結(jié)果與SCHMIDT 實(shí)驗(yàn)結(jié)果十分接近，相對(duì)誤差僅為1.49%。因此，可以認(rèn)為建立儲(chǔ)罐有限元模型、施加邊界條件、計(jì)算風(fēng)致屈曲時(shí)采用的方法合理準(zhǔn)確，能夠用于后續(xù)的分析計(jì)算及研究。

表2 有限元結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比Tab.2 Comparation of results of finite element model and experiment

2.2 儲(chǔ)罐風(fēng)致屈曲計(jì)算

對(duì)未發(fā)生地基沉降的儲(chǔ)罐，首先需確定風(fēng)壓系數(shù)。選擇EN 1993-4-1 定義的風(fēng)壓系數(shù)用于風(fēng)致屈曲計(jì)算（表3）。

表3 風(fēng)壓系數(shù)Tab.3 Wind pressure coefficient

采用弧長(zhǎng)法對(duì)未發(fā)生地基沉降儲(chǔ)罐的風(fēng)致屈曲進(jìn)行計(jì)算，記錄每個(gè)載荷子步的計(jì)算結(jié)果，得到儲(chǔ)罐在整個(gè)加載過(guò)程中的載荷-位移曲線（圖3）。

圖3 儲(chǔ)罐載荷-位移曲線Fig.3 Load-displacement curve of the storage tank

從圖3 可以看出：隨著載荷的增大，初始儲(chǔ)罐徑向位移與載荷呈線性關(guān)系，當(dāng)載荷增大至一定值時(shí)，儲(chǔ)罐徑向位移與載荷由原本的線性關(guān)系轉(zhuǎn)變?yōu)榉蔷€性關(guān)系，即儲(chǔ)罐發(fā)生了屈曲失穩(wěn)，此時(shí)的載荷大小即為儲(chǔ)罐的臨界屈曲載荷。儲(chǔ)罐發(fā)生屈曲失穩(wěn)后，儲(chǔ)罐失去原本的承載性能，因此儲(chǔ)罐在實(shí)際運(yùn)行中承受的風(fēng)載荷應(yīng)小于儲(chǔ)罐的臨界屈曲載荷，保證儲(chǔ)罐不會(huì)發(fā)生屈曲失穩(wěn)現(xiàn)象。根據(jù)儲(chǔ)罐載荷-位移曲線，得到未發(fā)生地基沉降儲(chǔ)罐風(fēng)致屈曲的臨界屈曲載荷為8183.6 Pa。

選擇儲(chǔ)罐發(fā)生屈曲失穩(wěn)前（p0=7 884 Pa）、發(fā)生屈曲失穩(wěn)時(shí)（p0=8 183.6 Pa）、發(fā)生屈曲失穩(wěn)后載荷繼續(xù)增大（p0=8 213.7 Pa）、發(fā)生屈曲失穩(wěn)后載荷出現(xiàn)減?。╬0=8 090 Pa）四種狀態(tài)的徑向位移云圖進(jìn)行對(duì)比，觀察儲(chǔ)罐屈曲失穩(wěn)過(guò)程中的變形情況（圖4）。

由圖4 可知：未發(fā)生地基沉降的儲(chǔ)罐因風(fēng)載荷導(dǎo)致屈曲時(shí)，屈曲發(fā)生在第四層罐壁處，且位于風(fēng)壓駐點(diǎn)附近；發(fā)生屈曲后屈曲位置的徑向位移不斷增大；屈曲形態(tài)表現(xiàn)為局部罐壁發(fā)生內(nèi)凹，使罐壁出現(xiàn)波浪形的“褶皺”。

圖4 儲(chǔ)罐徑向位移云圖Fig.4 Cloud map of radial displacement of the storage tank

3 地基沉降及風(fēng)載荷耦合條件下儲(chǔ)罐屈曲計(jì)算

3.1 沉降參數(shù)

選擇我國(guó)某10×104m3大型外浮頂儲(chǔ)罐的實(shí)測(cè)地基沉降數(shù)據(jù)用于計(jì)算，對(duì)該儲(chǔ)罐的實(shí)測(cè)地基沉降數(shù)據(jù)進(jìn)行傅里葉分解，將各階諧波進(jìn)行相加，即可得到罐周沉降與儲(chǔ)罐周向角的關(guān)系式，利用此關(guān)系式作出罐周沉降的曲線（圖5）。

圖5 儲(chǔ)罐罐周沉降Fig.5 Circumferential settlement of the storage tank

3.2 風(fēng)壓駐點(diǎn)周向角變化時(shí)儲(chǔ)罐的風(fēng)致屈曲

地基發(fā)生不均勻沉降后，儲(chǔ)罐罐壁各處均會(huì)產(chǎn)生不同程度的變形，當(dāng)風(fēng)從不同方向吹向儲(chǔ)罐時(shí)，承受風(fēng)壓作用的罐壁結(jié)構(gòu)并不一樣，從而對(duì)儲(chǔ)罐抵抗風(fēng)致屈曲的性能產(chǎn)生影響，使儲(chǔ)罐風(fēng)致屈曲的計(jì)算結(jié)果出現(xiàn)差異。以風(fēng)壓駐點(diǎn)周向角0°為起點(diǎn)，改變風(fēng)壓駐點(diǎn)周向角，每隔20°進(jìn)行一次計(jì)算，獲得不同風(fēng)壓駐點(diǎn)周向角下儲(chǔ)罐的臨界屈曲載荷及屈曲形態(tài)（表4）。

由表4 可知：不同風(fēng)壓駐點(diǎn)周向角下儲(chǔ)罐的臨界屈曲載荷存在明顯的差異，儲(chǔ)罐地基發(fā)生不均勻沉降后，罐壁各處發(fā)生不同程度的變形，進(jìn)而對(duì)儲(chǔ)罐抵抗風(fēng)致屈曲的性能產(chǎn)生或增強(qiáng)或削弱的作用。

表4 不同風(fēng)壓駐點(diǎn)周向角下儲(chǔ)罐臨界屈曲載荷Tab.4 Critical buckling load of the storage tank under various circumferential angles of wind stagnation point

為了研究地基不均勻沉降下儲(chǔ)罐的風(fēng)致屈曲機(jī)理，考察了不同風(fēng)壓駐點(diǎn)周向角下的地基沉降與儲(chǔ)罐臨界屈曲載荷之間的關(guān)系，以及地基沉降引起的罐壁變形與儲(chǔ)罐臨界屈曲載荷之間的關(guān)系。該儲(chǔ)罐不同風(fēng)壓駐點(diǎn)周向角下的地基沉降量如圖6a 所示，臨界屈曲載荷如圖6b 所示，不同風(fēng)壓駐點(diǎn)周向角下儲(chǔ)罐局部屈曲處因沉降導(dǎo)致的徑向位移如圖6c 所示。

由圖6 可知：不同風(fēng)壓駐點(diǎn)周向角下的儲(chǔ)罐地基沉降量與臨界屈曲載荷之間并不存在明顯的單調(diào)關(guān)系，且儲(chǔ)罐局部屈曲處因沉降導(dǎo)致的徑向位移與臨界屈曲載荷之間亦不存在明顯的單調(diào)關(guān)系。因此，不同風(fēng)壓駐點(diǎn)周向角下的地基沉降和儲(chǔ)罐局部屈曲處因沉降導(dǎo)致的徑向位移與儲(chǔ)罐臨界屈曲載荷之間皆無(wú)直接關(guān)系。

圖6 地基沉降、儲(chǔ)罐臨界屈曲載荷與徑向位移Fig.6 Foundation settlement，critical buckling load and radial displacement of the storage tank

根據(jù)以上結(jié)論可以推測(cè)，罐壁變形對(duì)儲(chǔ)罐臨界屈曲載荷的影響并不取決于某一處的變形大小，而是取決于該處附近整體的變形趨勢(shì)。將不同風(fēng)壓駐點(diǎn)周向角下的計(jì)算結(jié)果按照儲(chǔ)罐臨界屈曲載荷是否大于未沉降儲(chǔ)罐進(jìn)行劃分，并分別進(jìn)行詳細(xì)分析。

3.2.1 儲(chǔ)罐臨界屈曲載荷大于未沉降儲(chǔ)罐的情況

當(dāng)風(fēng)壓駐點(diǎn)周向角分別為-100、-80、-40、20、80、100、140、160 和180°時(shí)，地基不均勻沉降下儲(chǔ)罐的臨界屈曲載荷大于未沉降儲(chǔ)罐。根據(jù)以上風(fēng)壓駐點(diǎn)周向角下的計(jì)算結(jié)果，分別選擇以風(fēng)壓駐點(diǎn)周向角為中心±20°范圍內(nèi)、儲(chǔ)罐局部屈曲失穩(wěn)處同一高度附近的罐壁，獲得其因地基沉降導(dǎo)致的徑向位移，并由此得到該范圍內(nèi)罐壁的變形情況（圖7）。

由圖7 可知，在儲(chǔ)罐臨界屈曲載荷大于未沉降儲(chǔ)罐的風(fēng)壓駐點(diǎn)周向角下，僅圖7b、圖7f 和圖7h的風(fēng)壓駐點(diǎn)位于儲(chǔ)罐罐壁內(nèi)凹最大處附近，其余情況下風(fēng)壓駐點(diǎn)都位于儲(chǔ)罐罐壁外凸最大處附近。再將儲(chǔ)罐屈曲失穩(wěn)處的高度及周向角考慮在內(nèi)，進(jìn)行綜合對(duì)比，結(jié)果見(jiàn)表5。

表5 各計(jì)算工況屈曲及變形對(duì)比(儲(chǔ)罐臨界屈曲載荷大于未沉降儲(chǔ)罐)Tab.5 Comparison of buckling and deformation under various calculation conditions(critical buckling load of the tank is greater than that of the unsettled tank)

圖7 儲(chǔ)罐局部屈曲處附近罐壁因沉降導(dǎo)致的徑向位移(儲(chǔ)罐臨界屈曲載荷大于未沉降儲(chǔ)罐)Fig.7 Radial displacement of the tank wall due to settlement near the local buckling area of the storage tank(critical buckling load of the tank is greater than that of the unsettled tank)

由圖7 及表5 可以得出：當(dāng)風(fēng)壓駐點(diǎn)周向角為-80、100 和160°時(shí)，風(fēng)壓駐點(diǎn)位于罐壁內(nèi)凹最大處附近，不利于罐壁抵抗風(fēng)壓，但屈曲發(fā)生在兩加強(qiáng)圈之間，加強(qiáng)圈可提高罐壁抵抗風(fēng)壓的性能，兩種因素共同作用，最后使臨界屈曲載荷有所提高；當(dāng)風(fēng)壓駐點(diǎn)周向角為-100°時(shí)，風(fēng)壓駐點(diǎn)位于罐壁外凸最大處附近，同時(shí)屈曲發(fā)生在兩加強(qiáng)圈之間，兩種因素共同作用使得臨界屈曲載荷顯著增大；當(dāng)風(fēng)壓駐點(diǎn)周向角為-40、20、80、140 和180°時(shí)，雖然屈曲發(fā)生于第四層罐壁，加強(qiáng)圈的作用很小，但風(fēng)壓駐點(diǎn)位于罐壁外凸最大處附近，故臨界屈曲載荷亦有所增大。此外，無(wú)論風(fēng)壓駐點(diǎn)附近罐壁變形是外凸還是內(nèi)凹，屈曲都傾向于發(fā)生在罐壁相對(duì)凹陷的部位，當(dāng)風(fēng)壓駐點(diǎn)附近罐壁變形為外凸時(shí)，屈曲不會(huì)發(fā)生在風(fēng)壓駐點(diǎn)處，而是傾向于發(fā)生在罐壁的相對(duì)凹陷部位，屈曲失穩(wěn)處會(huì)不同程度地遠(yuǎn)離風(fēng)壓駐點(diǎn)；當(dāng)風(fēng)壓駐點(diǎn)附近罐壁變形為內(nèi)凹時(shí)，由于該處已經(jīng)是罐壁的相對(duì)凹陷部位，屈曲傾向于發(fā)生于此部位，因此屈曲失穩(wěn)處距離風(fēng)壓駐點(diǎn)往往不會(huì)太遠(yuǎn)。

3.2.2 儲(chǔ)罐臨界屈曲載荷小于未沉降儲(chǔ)罐的情況

當(dāng)風(fēng)壓駐點(diǎn)周向角分別為-160、-140、-120、-60、-20、0、40、60 和120°時(shí)，地基不均勻沉降下儲(chǔ)罐的臨界屈曲載荷小于未沉降儲(chǔ)罐。根據(jù)以上風(fēng)壓駐點(diǎn)周向角下的計(jì)算結(jié)果，分別選擇以風(fēng)壓駐點(diǎn)周向角為中心±20°范圍內(nèi)、儲(chǔ)罐局部屈曲失穩(wěn)處同一高度附近的罐壁，獲得其因地基沉降導(dǎo)致的徑向位移，并由此得到該范圍內(nèi)罐壁的變形情況（圖8）。

由圖8 可知：在儲(chǔ)罐臨界屈曲載荷小于未沉降儲(chǔ)罐的風(fēng)壓駐點(diǎn)周向角下，所有風(fēng)壓駐點(diǎn)都位于儲(chǔ)罐罐壁內(nèi)凹最大處附近。再將儲(chǔ)罐屈曲失穩(wěn)處的高度及周向角考慮在內(nèi)，進(jìn)行綜合對(duì)比，結(jié)果如表6所示。

圖8 儲(chǔ)罐局部屈曲處附近罐壁因沉降導(dǎo)致的徑向位移(儲(chǔ)罐臨界屈曲載荷小于未沉降儲(chǔ)罐)Fig.8 Radial displacement of the tank wall due to settlement near the local buckling area of the storage tank(critical buckling load of the tank is less than that of the unsettled tank)

由表6 可知：當(dāng)風(fēng)壓駐點(diǎn)周向角為60°時(shí)，儲(chǔ)罐局部屈曲失穩(wěn)處的罐壁高度為12.6m，即儲(chǔ)罐屈曲處位于儲(chǔ)罐自底向上第一個(gè)加強(qiáng)圈與第二個(gè)加強(qiáng)圈之間，此時(shí)的風(fēng)壓駐點(diǎn)位于罐壁內(nèi)凹最大處附近，而風(fēng)壓作用會(huì)使儲(chǔ)罐罐壁發(fā)生內(nèi)凹直至屈曲。因此地基沉降引起的罐壁內(nèi)凹不僅無(wú)法抵消風(fēng)壓的作用，反而會(huì)減弱罐壁抵抗風(fēng)壓的性能，進(jìn)而使儲(chǔ)罐臨界屈曲載荷減小，以上兩種因素綜合作用，最終使儲(chǔ)罐的臨界屈曲載荷相對(duì)未沉降儲(chǔ)罐略有減小。當(dāng)風(fēng)壓駐點(diǎn)周向角為-160、-140、-120、-60、-20、0、40 和120°時(shí)，儲(chǔ)罐局部屈曲失穩(wěn)處位于第四層罐壁，加強(qiáng)圈對(duì)罐壁抵抗風(fēng)壓的性能貢獻(xiàn)很小。

表6 各計(jì)算工況屈曲及變形對(duì)比(儲(chǔ)罐臨界屈曲載荷小于未沉降儲(chǔ)罐)Tab.6 Comparison of buckling and deformation under various calculation conditions(critical buckling load of the tank is less than that of the unsettled tank)

此外，由圖8 及表6 可以得出：對(duì)于儲(chǔ)罐臨界屈曲載荷小于未沉降儲(chǔ)罐的情況，風(fēng)壓駐點(diǎn)附近的罐壁變形均表現(xiàn)為內(nèi)凹，而屈曲也都傾向于發(fā)生在罐壁相對(duì)凹陷部位。

4 結(jié)論

儲(chǔ)罐地基不均勻沉降及儲(chǔ)罐承受的風(fēng)載荷皆會(huì)影響儲(chǔ)罐安全運(yùn)行。通過(guò)建立大型儲(chǔ)罐數(shù)值仿真模型，研究了儲(chǔ)罐風(fēng)致屈曲及地基不均勻沉降下儲(chǔ)罐風(fēng)致屈曲的有關(guān)規(guī)律。通過(guò)對(duì)不同風(fēng)壓駐點(diǎn)周向角下儲(chǔ)罐的風(fēng)致屈曲進(jìn)行計(jì)算分析，得到以下結(jié)論：不同周向角下的地基沉降與儲(chǔ)罐臨界屈曲載荷之間不存在直接聯(lián)系，儲(chǔ)罐局部屈曲處因地基沉降導(dǎo)致的徑向位移與儲(chǔ)罐臨界屈曲載荷之間亦不存在直接聯(lián)系；儲(chǔ)罐臨界屈曲載荷的變化由儲(chǔ)罐發(fā)生風(fēng)致屈曲的位置及風(fēng)壓駐點(diǎn)附近罐壁變形情況共同決定；地基沉降發(fā)生后，部分風(fēng)壓駐點(diǎn)周向角下儲(chǔ)罐發(fā)生風(fēng)致屈曲的位置會(huì)發(fā)生改變，由第四層罐壁變?yōu)閮?chǔ)罐由下至上第一及第二個(gè)加強(qiáng)圈之間，此因會(huì)使儲(chǔ)罐臨界屈曲載荷增大，同時(shí)改變儲(chǔ)罐的屈曲形態(tài)；當(dāng)風(fēng)壓駐點(diǎn)附近罐壁因地基沉降引起的變形表現(xiàn)為外凸時(shí)，儲(chǔ)罐臨界屈曲載荷增大，而當(dāng)風(fēng)壓駐點(diǎn)附近罐壁因地基沉降引起的變形表現(xiàn)為內(nèi)凹時(shí)，儲(chǔ)罐臨界屈曲載荷減??；無(wú)論風(fēng)壓駐點(diǎn)附近罐壁因地基沉降引起的罐壁變形表現(xiàn)為外凸還是內(nèi)凹，都不會(huì)改變儲(chǔ)罐的屈曲形態(tài)，并且儲(chǔ)罐的風(fēng)致屈曲都傾向于發(fā)生在罐壁相對(duì)凹陷部位。