LNG儲(chǔ)罐流固耦合工況下的地震響應(yīng)分析

張?jiān)品澹?杜展祥， 李騰飛， 滕振超

（東北石油大學(xué)，黑龍江大慶 163318）

中國(guó)包括管道天然氣和液化天然氣（LNG）在內(nèi)的天然氣進(jìn)口在2021年同比增長(zhǎng)19.9%，達(dá)到創(chuàng)紀(jì)錄的1.216 億t，中國(guó)已成為全球天然氣需求增長(zhǎng)最旺盛的國(guó)家之一. 相較于石油和煤炭而言，天然氣更加環(huán)保，在天然氣的進(jìn)出口和使用過(guò)程中，儲(chǔ)存液化天然氣的裝置顯得尤為重要，液化天然氣主要儲(chǔ)存在LNG儲(chǔ)罐中，其安全性、穩(wěn)定性和耐久性一直是人們的重點(diǎn)關(guān)注對(duì)象. 地震是常見(jiàn)的自然災(zāi)害，會(huì)對(duì)LNG儲(chǔ)罐的安全性產(chǎn)生巨大危害. 當(dāng)?shù)卣鸢l(fā)生時(shí)，LNG 儲(chǔ)罐在內(nèi)部LNG 液體作用下會(huì)產(chǎn)生變形或運(yùn)動(dòng)，而儲(chǔ)罐的變形或運(yùn)動(dòng)又會(huì)反過(guò)來(lái)影響內(nèi)部LNG 液體的運(yùn)動(dòng). 在進(jìn)行地震響應(yīng)分析時(shí)往往需要考慮流固耦合作用，進(jìn)行耦合時(shí)兩者會(huì)有交界面的存在，并且在此交界面上兩者的參數(shù)傳遞是雙向的［1］.

關(guān)于儲(chǔ)罐的流固耦合，國(guó)外的研究起源較早，早期對(duì)儲(chǔ)罐模型進(jìn)行研究時(shí)，往往對(duì)儲(chǔ)罐罐壁的彈性進(jìn)行忽略，假定儲(chǔ)罐罐壁為剛性進(jìn)行計(jì)算. 此種計(jì)算模型由Hoskins和Jacobsen［2］提出，被稱為質(zhì)量-彈簧系統(tǒng)模型，后來(lái)被改進(jìn)：Housner［3］分析了流體容器在水平加速度作用下產(chǎn)生的水動(dòng)力壓力，將承受水平地表激振時(shí)，儲(chǔ)罐內(nèi)液體的振動(dòng)響應(yīng)分為了沖擊運(yùn)動(dòng)和對(duì)流運(yùn)動(dòng)，做沖擊運(yùn)動(dòng)的液體質(zhì)量稱為剛性沖擊質(zhì)量部分，其余液體則做獨(dú)立的自由晃動(dòng)，其質(zhì)量稱為對(duì)流質(zhì)量部分. Sharari等［4］利用有限元分析軟件ABAQUS對(duì)1994年北嶺地震和1995年神戶地震影響下的LNG儲(chǔ)罐進(jìn)行了非線性動(dòng)力分析，得出儲(chǔ)罐在不同基礎(chǔ)類型下的地震響應(yīng)情況. Subba和Gorla［5］采用有限元程序ALGOR將熱剖面與結(jié)構(gòu)設(shè)計(jì)進(jìn)行了耦合分析. Kim等［6］考慮了流體-結(jié)構(gòu)-油之間的相互作用提出了兩種簡(jiǎn)化的液化天然氣儲(chǔ)罐地震反應(yīng)分析模型. Lahlou和Rachid［7］采用ANSYS有限元軟件考慮了LNG儲(chǔ)罐泄露工況，對(duì)低溫下的LNG儲(chǔ)罐溫度場(chǎng)進(jìn)行了分析研究，得出其相關(guān)分布情況，然后對(duì)儲(chǔ)罐各節(jié)點(diǎn)的溫度和熱應(yīng)變等進(jìn)行了計(jì)算.

國(guó)內(nèi)研究中，王大鈞等［8］對(duì)殼-液耦合系統(tǒng)進(jìn)行了研究，得出低頻大幅重力波的產(chǎn)生是由于液體在柱形彈性容器中受到高頻激勵(lì)的作用，并且對(duì)中國(guó)古文物龍洗現(xiàn)象的產(chǎn)生進(jìn)行了介紹. 王暉等［9］對(duì)儲(chǔ)液容器流固耦合模態(tài)的變化進(jìn)行了研究，在研究過(guò)程中采用了強(qiáng)耦合的研究方法，探究了儲(chǔ)罐固有頻率的影響因素：儲(chǔ)罐內(nèi)液體的深度和儲(chǔ)罐結(jié)構(gòu)的剛度. 張?jiān)品宓龋?0］利用ANSYS 有限元軟件建立了內(nèi)罐泄露條件下圓柱形LNG 儲(chǔ)罐混凝土外墻的有限元模型，將動(dòng)液壓力等效為附加質(zhì)量，在半液和滿液狀態(tài)下，求得外墻自振頻率環(huán)向波變化曲線. 袁朝慶等［11］利用有限元軟件ADINA 建立了泄露工況下外罐空罐和滿罐兩種有限元模型，并對(duì)儲(chǔ)罐施加EI-Centro 地震波，考慮液固耦合作用，探究其地震響應(yīng)情況，通過(guò)分析得出儲(chǔ)罐內(nèi)部液體會(huì)對(duì)儲(chǔ)罐的地震響應(yīng)產(chǎn)生影響，且滿罐地震作用下的加速度和位移均比空罐的大.翟希梅等［12］采用直接耦合法對(duì)LNG 儲(chǔ)罐進(jìn)行了精細(xì)化建模，考慮了四種工況分別為：空罐、正常工作時(shí)滿液位、滿液位泄露、半液位泄露，探究鋼制內(nèi)罐和預(yù)應(yīng)力混凝土外罐的振動(dòng)特性，然后對(duì)液體和罐體的相互作用進(jìn)行分析，得出液體對(duì)LNG 儲(chǔ)罐罐體振動(dòng)特性的影響規(guī)律. 潘德濤［13］建立了LNG 儲(chǔ)罐外罐空罐和內(nèi)罐泄露工況下外罐滿罐兩種有限元模型，得出了兩種不同工況的儲(chǔ)罐在水平地震作用下的位移、加速度、環(huán)向應(yīng)力等. 蘇娟等［14］借助于ANSYS 有限元軟件，采用熱-固間接耦合方法，研究了儲(chǔ)罐泄露工況下其溫度場(chǎng)分布. 呂克克［15］利用ANSYS 對(duì)LNG 全容型儲(chǔ)罐進(jìn)行了三維有限元建模，對(duì)多種工況下LNG 儲(chǔ)罐外罐的溫度場(chǎng)和應(yīng)力場(chǎng)分布進(jìn)行了研究. 魏新［16］利用ANSYS 有限元軟件對(duì)二維LNG 儲(chǔ)罐結(jié)構(gòu)進(jìn)行了精細(xì)化建模，在建模過(guò)程中對(duì)所有保冷層也進(jìn)行了細(xì)致化的建模，分析了LNG儲(chǔ)罐在正常工作，內(nèi)罐泄露以及臨近火災(zāi)三種常見(jiàn)工況下的溫度場(chǎng)和應(yīng)力場(chǎng)不同工況下的溫度、應(yīng)力和位移變化情況. 張世凱［17］以16萬(wàn)m3LNG全容式儲(chǔ)罐為研究對(duì)象，考慮液固耦合的作用，運(yùn)用時(shí)程分析的方法研究了不同液位高度、不同地震波、不同作用方向下的地震響應(yīng)情況，并對(duì)樁-土結(jié)構(gòu)在靜力和動(dòng)力作用下的應(yīng)力和變形情況進(jìn)行了研究. 尚潤(rùn)丹［18］基于流固耦合原理，對(duì)在近場(chǎng)地震波作用下LNG 儲(chǔ)罐隔震和不隔震兩種不同儲(chǔ)罐結(jié)構(gòu)的地震響應(yīng)進(jìn)行研究，并對(duì)LNG 儲(chǔ)罐各工況下的地震響應(yīng)特點(diǎn)進(jìn)行了分析總結(jié)，探究了儲(chǔ)罐隔震前后加速度、應(yīng)力和位移的相關(guān)變化情況. 陳泓宇［19］考慮了LNG 儲(chǔ)罐流固耦合作用，運(yùn)用時(shí)程分析法，采用有限元軟件ANSYS 對(duì)LNG 儲(chǔ)罐混凝土外罐、鋼制內(nèi)罐和穹頂三個(gè)部分分別進(jìn)行了有限元分析，得出LNG 儲(chǔ)罐各主要部位的應(yīng)力應(yīng)變和破壞形式. 以上的分析且相關(guān)研究主要集中在其振動(dòng)特性分析，振動(dòng)周期求解以及求地震作用下的位移、加速度、環(huán)向應(yīng)力等方面，很少對(duì)地震作用下LNG 儲(chǔ)罐的流固耦合綜合來(lái)進(jìn)行考慮，本文將圍繞此方向展開(kāi)研究.

1 有限元模型的建立

本文的研究對(duì)象為某實(shí)際工程中容量為160 000 m3的大型全容式LNG 儲(chǔ)罐，其由預(yù)應(yīng)力混凝土外罐、9%鎳鋼內(nèi)罐，以及內(nèi)罐和外罐之間的保溫層構(gòu)成. 預(yù)應(yīng)力混凝土外罐內(nèi)徑為82 m，高度為38.6 m，厚度為800 mm，內(nèi)罐外徑為80 m，內(nèi)罐和外罐之間為保溫層，保溫層主要由兩部分構(gòu)成：200 mm 厚的玻璃纖維氈層、800 mm厚的膨脹珍珠巖層. 混凝土穹頂跨度為82 m，矢高約為10.98 m. 儲(chǔ)罐罐底采用隔熱結(jié)構(gòu)，從上到下依次為：450 mm厚的泡沫玻璃磚層，100 mm厚的混凝土層. 其內(nèi)部結(jié)構(gòu)具體情況如圖1所示.

圖1 LNG儲(chǔ)罐內(nèi)部結(jié)構(gòu)示意圖Fig.1 Internal structure diagram of LNG storage tank

2 流-固耦合有限元模型的驗(yàn)證

根據(jù)項(xiàng)忠權(quán)等研究油管的動(dòng)力特性提出的意見(jiàn)，計(jì)算LNG儲(chǔ)罐外罐空罐狀態(tài)下采用的模態(tài)頻率：

式中：f1為梁剪切振動(dòng)基本周期；α為彎曲變形的影響系數(shù)；ξ1為截面變形影響系數(shù)；ρs為外罐壁所用材料的密度；H為外罐壁的高度；μ為罐壁材料的泊松比；G為罐壁材料的剪切模量；R為罐體內(nèi)半徑.

利用ABAQUS 有限元軟件得其前四階模態(tài)云圖如下圖2所示.對(duì)空罐模型的頻率進(jìn)行提取，得表1.

圖2 LNG儲(chǔ)罐空罐模態(tài)云圖Fig.2 Modal cloud diagram of empty LNG tank

表1 前十階模態(tài)頻率表Tab.1 Modal frequencies of first ten orders

空罐利用項(xiàng)忠權(quán)法［20］有限元模型的驗(yàn)證結(jié)果如表2所示.

表2 空罐有限元模型驗(yàn)證Tab.2 Verification of empty tank finite element model

3 LNG儲(chǔ)罐流-固耦合地震響應(yīng)分析

3.1 LNG儲(chǔ)罐地震波和相關(guān)參數(shù)的選取

3.1.1 地震波的選取和輸入

為了更加全面真實(shí)地反映儲(chǔ)罐在地震作用下真實(shí)受力情況，綜合考慮建筑物的場(chǎng)地類別以及地震設(shè)計(jì)分組，選取兩種實(shí)際自然下的強(qiáng)震記錄和一組人工模擬的加速度時(shí)程曲線來(lái)對(duì)其地震響應(yīng)情況進(jìn)行分析.

根據(jù)文獻(xiàn)［21］，并結(jié)合LNG儲(chǔ)罐結(jié)構(gòu)的基本情況，擬采用地震波的相關(guān)參數(shù)如下：

1）EI-Centro波，間隔0.02 s，持續(xù)時(shí)間53.76 s，加速度峰值出現(xiàn)在第2.14 s，峰值為341.7 cm/s2，場(chǎng)地土屬Ⅱ～Ⅲ類.

2）Taft 波，間隔0.02 s，持續(xù)時(shí)間54.40 s，加速度峰值出現(xiàn)在第3.72 s，峰值為175.9 cm/s2，場(chǎng)地土屬Ⅱ類.

3）SHM2人工波，間隔0.02 s，持續(xù)時(shí)間78.64 s，加速度峰值出現(xiàn)在第13 s，峰值為35 cm/s2，場(chǎng)地土屬Ⅳ類.

三類地震波加速度特性和時(shí)程曲線分別如表3、圖3 所示. 根據(jù)抗震設(shè)計(jì)規(guī)范推薦的相關(guān)方法，本文中輸入的EI-Centro地震波、Taft波和SHM2人工波，持續(xù)時(shí)間均為25 s.

表3 三類地震波加速度特性Tab.3 Acceleration characteristics of three kinds of seismic waves

圖3 各地震波加速度時(shí)程曲線Fig.3 Time history curves of local seismic wave accelerations

3.1.2 地震波的施加

按抗震設(shè)防烈度為8度設(shè)防，設(shè)計(jì)基本地震加速度值為0.20 g，對(duì)地震波進(jìn)行調(diào)整，然后將其按圖4方式施加于罐底，并將儲(chǔ)罐罐壁劃分為0～13共14個(gè)節(jié)點(diǎn)，對(duì)各個(gè)節(jié)點(diǎn)的地震響應(yīng)情況進(jìn)行提取，如圖5所示.

圖4 X向地震波施加方式示意圖Fig.4 Schematic diagram of X-direction seismic wave application mode

圖5 LNG儲(chǔ)罐外罐節(jié)點(diǎn)示意圖（單位：m）Fig.5 Schematic diagram of outer tank nodes of LNG storage tank

3.2 EI-Centro地震波下的地震響應(yīng)

3.2.1 加速度時(shí)程分析

當(dāng)探究液固耦合對(duì)加速度時(shí)程分析的影響時(shí)，采用附加質(zhì)量法，對(duì)四種不同溫度工況分別考慮液固耦合的作用，得其各工況下的加速度時(shí)程曲線如圖6所示，然后對(duì)不同工況下的加速度時(shí)程曲線提取加速度峰值，如表4所示.

表4 EI-Centro波液固耦合工況下的加速度峰值提取表Tab.4 Acceleration peak extraction table under liquid-solid coupling condition of EI-Centro wave

圖6 EI-Centro波液固耦合工況下的加速度時(shí)程曲線Fig.6 Acceleration time history curves of EI Centro wave under liquid-solid coupling condition

對(duì)同一節(jié)點(diǎn)橫向?qū)Ρ葴囟雀淖兊墓r時(shí)，發(fā)現(xiàn)溫度改變對(duì)加速度峰值影響不大. 當(dāng)對(duì)流固耦合進(jìn)行考慮時(shí)，對(duì)中部節(jié)點(diǎn)的加速度峰值和上部節(jié)點(diǎn)的加速度峰值進(jìn)行比較時(shí)，可以得出上部節(jié)點(diǎn)的加速度峰值明顯大于中部節(jié)點(diǎn)的加速度峰值，分析原因主要是由于罐底與承臺(tái)間有約束，而罐頂相當(dāng)于自由端.

3.2.2 應(yīng)力時(shí)程分析

當(dāng)考慮流固耦合對(duì)應(yīng)力時(shí)程分析的影響時(shí)，需控制溫度相同，取相同溫度工況為283.15 K，僅改變是否考慮流固耦合這一條件，得出應(yīng)力沿罐壁高度的變化如圖7所示.

從圖7 中可以直觀看出不考慮流固耦合時(shí)的應(yīng)力明顯小于考慮流固耦合時(shí)的應(yīng)力，當(dāng)不考慮流固耦合時(shí)其應(yīng)力最大值為4.83×107Pa，當(dāng)考慮流固耦合時(shí)其應(yīng)力最大值為12×107Pa，應(yīng)力最大值相差7.17×107Pa，相差很大. 當(dāng)不考慮流固耦合時(shí)其應(yīng)力最小值為0.644×107Pa，當(dāng)考慮流固耦合時(shí)其應(yīng)力最小值為1.73×107Pa，應(yīng)力最小值兩者相差1.086×107Pa，其最小值兩者相差較大.

圖7 EI-Centro地震波下應(yīng)力沿罐壁高度的變化情況Fig.7 Changes of stresses along the tank wall height under EI-Centro wave

3.2.3 位移時(shí)程分析

當(dāng)研究流固耦合對(duì)位移時(shí)程分析的影響時(shí)，控制溫度條件恒定，取相同溫度工況為283.15 K，僅對(duì)是否考慮流固耦合這一條件作出改變，得出位移沿高度的變化情況如圖8所示.

圖8 中可以清晰明了地看出流固耦合對(duì)位移時(shí)程分析的影響，LNG儲(chǔ)罐位移最大處主要集中在儲(chǔ)罐罐壁的中部位置，且是否考慮流固耦合對(duì)儲(chǔ)罐位移的影響較大. 當(dāng)不考慮流固耦合時(shí)，283.15 K溫度工況下其位移最大值為22.344 6 mm，位移最小值為22.286 4 mm. 當(dāng)考慮流固耦合時(shí)，283.15 K溫度工況下其位移最大值為22.488 6 mm，位移最小值為22.297 6 mm. 兩者進(jìn)行比較，其最大位移之差為0.144 mm，差別較大.

圖8 EI-Centro地震波下位移沿罐壁高度的變化情況Fig.8 Changes of displacements along the tank wall height under EI-Centro wave

通過(guò)上述分析及上表可以了解到，LNG 儲(chǔ)罐在X向的EI-Centro 波作用下，最大應(yīng)力一般出現(xiàn)在儲(chǔ)罐底部，因此在地震頻發(fā)地區(qū)需要對(duì)LNG 儲(chǔ)罐罐底進(jìn)行加強(qiáng)，防止儲(chǔ)罐在地震作用下發(fā)生震害，危及人民生命財(cái)產(chǎn)安全. 在地震發(fā)生時(shí)，LNG儲(chǔ)罐由于內(nèi)部液體與儲(chǔ)罐的相互作用將明顯導(dǎo)致儲(chǔ)罐應(yīng)力增大，儲(chǔ)罐的最大位移一般出現(xiàn)在儲(chǔ)罐中部，LNG儲(chǔ)罐在地震作用下由于內(nèi)部液體與罐體的相互作用將會(huì)導(dǎo)致位移明顯增大.

3.3 Taft地震波下的地震響應(yīng)

3.3.1 加速度時(shí)程分析

當(dāng)探究液固耦合對(duì)加速度時(shí)程分析的影響時(shí)，與EI-Centro地震波類似，采用附加質(zhì)量法，對(duì)四種不同溫度工況分別考慮液固耦合的作用，得出273.15、283.15、293.15、303.15 K四種不同溫度工況下液固耦合的加速度時(shí)程曲線如圖9所示，然后對(duì)不同工況下的加速度時(shí)程曲線提取加速度峰值，如表5所示.

表5 Taft波液固耦合工況下的加速度峰值提取表Tab.5 Peak acceleration extraction table under Taft wave liquid-solid coupling condition

考慮中部節(jié)點(diǎn)和頂部節(jié)點(diǎn)的加速度時(shí)，可以明顯看出頂部節(jié)點(diǎn)13的加速度峰值無(wú)論在不同溫度工況下還是在是否考慮流固耦合工況下均明顯大于中部節(jié)點(diǎn)7的加速度峰值，對(duì)同一節(jié)點(diǎn)橫向?qū)Ρ葴囟雀淖兊墓r時(shí)，發(fā)現(xiàn)溫度改變對(duì)加速度峰值影響不大，與上述EI-Centro地震波下加速度的時(shí)程分析相類似.

3.3.2 應(yīng)力時(shí)程分析

當(dāng)考慮流固耦合對(duì)應(yīng)力時(shí)程分析的影響時(shí)，需控制溫度相同，取相同溫度工況為283.15 K，僅改變是否考慮流固耦合這一條件，得出應(yīng)力沿罐壁高度的變化如圖10所示.

當(dāng)對(duì)流固耦合工況進(jìn)行考慮時(shí)，應(yīng)力沿罐壁高度相關(guān)分布情況如上所示，可以明顯看到考慮流固耦合時(shí)應(yīng)力遠(yuǎn)大于沒(méi)有考慮流固耦合時(shí)的應(yīng)力.

3.3.3 位移時(shí)程分析

當(dāng)考慮流固耦合工況時(shí)，控制溫度條件不變，選取恒定溫度工況為283.15 K，僅就是否考慮流固耦合進(jìn)行研究，結(jié)果如圖11所示.

圖11 展示了是否考慮流固耦合對(duì)位移時(shí)程分析的影響，得出位移沿罐壁高度的相關(guān)變化情況. 當(dāng)考慮流固耦合時(shí)最大位移為1.700 61 mm，最小位移1.623 82 mm，不考慮流固耦合時(shí)最大位移為1.684 71 mm，最小位移為1.170 71 mm，位移最小值相差0.453 11 mm.

圖11 Taft波下位移沿罐壁高度的變化情況Fig.11 Changes of displacements along the tank wall height under Taft wave

由圖10 可知，對(duì)LNG 儲(chǔ)罐外罐施加X(jué)向Taft 地震波，儲(chǔ)罐最大應(yīng)力位于儲(chǔ)罐罐底和罐壁與穹頂?shù)倪B接處，需要對(duì)這些薄弱部位進(jìn)行加強(qiáng). 而是否考慮流固耦合時(shí)兩者位移差為0.015 90 mm，因此在對(duì)地震響應(yīng)進(jìn)行分析時(shí)應(yīng)著重考慮流固耦合的影響.

圖10 Taft波下應(yīng)力沿罐壁高度的變化Fig.10 Changes of stresses along the tank wall height under Taft wave

3.4 SHM2人工波下的地震響應(yīng)

3.4.1 加速度時(shí)程分析

與EI-Centro地震波和Taft地震波類似，當(dāng)探究液固耦合對(duì)加速度時(shí)程分析的影響時(shí)，采用附加質(zhì)量法，對(duì)273.15、283.15、293.15、303.15 K四種不同溫度工況分別考慮液固耦合的作用，得出四種不同溫度工況下液固耦合的加速度時(shí)程曲線如圖12所示，然后對(duì)不同工況下的加速度時(shí)程曲線提取加速度峰值（表6）.

表6 SHM2人工波液固耦合工況下的加速度峰值提取表Tab.6 Peak acceleration extraction table of SHM2 artificial wave under liquid-solid coupling condition

圖12 SHM2人工波液固耦合工況下的加速度時(shí)程曲線Fig.12 Acceleration time history curves of SHM2 artificial wave under liquid-solid coupling condition

對(duì)于頂部節(jié)點(diǎn)而言，其加速度峰值相較于中部節(jié)點(diǎn)7而言，明顯大于中部節(jié)點(diǎn)7的加速度峰值，這種現(xiàn)象出現(xiàn)的原因主要是因?yàn)長(zhǎng)NG儲(chǔ)罐外罐壁底部受到承臺(tái)的約束，導(dǎo)致外罐壁底部的加速度幾乎為零，而外罐壁中部相對(duì)于底部而言受承臺(tái)約束的影響較小，因此其加速度比罐壁底部大，穹頂無(wú)約束相當(dāng)于自由端，

地震作用下的加速度最大，因此頂部節(jié)點(diǎn)13的加速度峰值要大于中部節(jié)點(diǎn)7的加速度峰值. 與上述EI-Centro地震波下加速度的時(shí)程分析相類似，對(duì)同一節(jié)點(diǎn)橫向?qū)Ρ葴囟雀淖兊墓r時(shí)，同樣發(fā)現(xiàn)溫度改變對(duì)加速度峰值影響不大.

3.4.2 應(yīng)力時(shí)程分析

當(dāng)探究流固耦合對(duì)應(yīng)力時(shí)程分析的影響時(shí)，需控制溫度相同，取相同溫度工況為283.15 K，僅改變是否考慮流固耦合這一條件，得出罐壁應(yīng)力沿高度的變化（圖13）.

通過(guò)圖13可以，直觀看出是否考慮流固耦合工況下應(yīng)力沿罐壁高度的相關(guān)變化情況，LNG儲(chǔ)罐罐壁考慮流固耦合時(shí)其應(yīng)力明顯大于罐壁不考慮流固耦合時(shí)的應(yīng)力. 當(dāng)不考慮流固耦合時(shí)其最大應(yīng)力為10.7×107Pa，當(dāng)考慮流固耦合時(shí)其最大應(yīng)力為31.9×107Pa，兩者相差21.2×107Pa，考慮流固耦合時(shí)的應(yīng)力值幾乎為不考慮流固耦合應(yīng)力值的三倍，差別很大. 對(duì)于其最小應(yīng)力值而言，當(dāng)不考慮流固耦合時(shí)其最小應(yīng)力值為2.97×107Pa，而考慮流固耦合時(shí)最小應(yīng)力值為11.7×107Pa，兩者相差8.73×107Pa.

圖13 SHM2人工波下應(yīng)力沿罐壁高度的變化情況Fig.13 Changes of stresses along the tank wall height under SHM2 artificial wave

3.4.3 位移時(shí)程分析

當(dāng)考慮流固耦合對(duì)位移時(shí)程分析的影響時(shí)，首先需要控制溫度條件恒定，取相同溫度工況為283.15 K，僅對(duì)是否考慮流固耦合這一條件作出改變，得出位移沿高度的變化情況（圖14）.

從圖14 中可以直觀看出，當(dāng)僅考慮流固耦合時(shí)，其最大位移位于罐壁的頂部，而不考慮罐體內(nèi)部液體的作用時(shí)，罐壁的最大位移位于罐壁的底部. 當(dāng)不考慮流固耦合時(shí)罐壁的最大位移為25.307 mm，當(dāng)考慮流固耦合時(shí)罐壁的最大位移為26.253 mm，兩者相差0.946 mm，差別很大. 當(dāng)不考慮流固耦合時(shí)罐壁的最小位移為24.971 mm，當(dāng)考慮流固耦合時(shí)罐壁的最小位移為25.308 mm，兩者相差0.337 mm，差別較大.

圖14 SHM2人工波下位移沿罐壁高度的變化情況Fig.14 Changes of displacements along the tank wall height under SHM2 artificial wave

通過(guò)上述分析及上表可知，LNG儲(chǔ)罐在X向SHM2人工波下的最大應(yīng)力一般位于儲(chǔ)罐罐壁的底部和頂部，因此在設(shè)計(jì)和建造的過(guò)程中應(yīng)重點(diǎn)加強(qiáng)這些薄弱部位，防止儲(chǔ)罐在地震作用下產(chǎn)生損害. 對(duì)于位移而言，改變流固耦合時(shí)其位移變化為0.946 mm，位移變化較大，在地震發(fā)生時(shí)，由于罐內(nèi)液體與罐壁的相互作用，將導(dǎo)致罐壁的應(yīng)力和位移明顯增大.

4 結(jié)論

本章利用附加質(zhì)量的方法探究了流固耦合下的地震響應(yīng)情況，得出LNG儲(chǔ)罐混凝土外罐在EI-Centro地震波、Taft地震波及SHM2人工波下的加速度時(shí)程曲線，然后研究了同一溫度下是否考慮流固耦合工況下的應(yīng)力時(shí)程曲線和位移時(shí)程曲線，得出不同工況下應(yīng)力和位移沿罐壁高度的相關(guān)變化情況，并對(duì)結(jié)果進(jìn)行了分析.

1）通過(guò)LNG儲(chǔ)罐外罐在地震作用下的時(shí)程分析可得，內(nèi)罐在正常工作工況下，外罐加速度、位移、應(yīng)力在溫度改變時(shí)變化并不明顯. 但當(dāng)內(nèi)罐液體泄露，外罐在地震作用下，內(nèi)部液體與外罐罐壁發(fā)生液固耦合作用，產(chǎn)生較大的沖擊力，使加速度、位移、應(yīng)力均有一定程度的增大，液固耦合的影響對(duì)儲(chǔ)罐地震響應(yīng)的影響大.

2）對(duì)外罐加速度進(jìn)行時(shí)程分析時(shí)，可知外罐的加速度峰值從底部到頂部逐漸增大，主要是由于外罐底部與承臺(tái)相連接，在地震發(fā)生時(shí)會(huì)抑制罐體的晃動(dòng)，而頂部沒(méi)有約束，相當(dāng)于自由端，因此在地震作用下會(huì)有較大的加速度峰值產(chǎn)生.

3）當(dāng)對(duì)外罐的應(yīng)力進(jìn)行時(shí)程分析時(shí)，可以得出外罐的最大應(yīng)力幾乎全部位于外罐的底部，少部分位于穹頂和外罐的連接處，這與儲(chǔ)罐的破壞形式相符，產(chǎn)生的主要原因是內(nèi)部液體對(duì)外罐壁的沖擊壓力，因此在對(duì)LNG儲(chǔ)罐設(shè)計(jì)和建造的過(guò)程中需要對(duì)薄弱部位進(jìn)行加強(qiáng)，保障儲(chǔ)罐的安全性和穩(wěn)定性.

4）通過(guò)對(duì)三種地震波分別進(jìn)行加速度、位移、應(yīng)力時(shí)程分析可知不同波形地震波下的地震響應(yīng)是有所差別的，在SHM2 人工波地震動(dòng)作用下結(jié)構(gòu)響應(yīng)最為明顯，而在另外兩條地震波作用下結(jié)構(gòu)響應(yīng)相對(duì)弱一些，因此在對(duì)結(jié)構(gòu)進(jìn)行地震響應(yīng)分析時(shí)，需對(duì)結(jié)構(gòu)進(jìn)行多種地震動(dòng)下的參數(shù)化分析，保證儲(chǔ)罐在多種地震波下安全工作.