管友海，賈娟娟，林 楠，蘇 寧

（1. 中國石油大學(華東)，山東 青島 266580; 2. 中國鐵建十六局集團有限公司，北京100018）

大型LNG儲罐液固耦合模態(tài)分析

管友海1，賈娟娟1，林 楠2，蘇 寧1

（1. 中國石油大學(華東)，山東 青島 266580; 2. 中國鐵建十六局集團有限公司，北京100018）

利用ADIAN軟件中的Lanczos特征值算法，以山東青島市董家口一座16萬m3大型預應力混凝土儲罐為研究對象，考慮液固耦合效應，對LNG儲罐內外罐進行模態(tài)分析，比較了考慮液固耦合的理論計算結果與有限元計算的結果，證明了ADINA內計算的正確性及有效性。對比分析外罐對罐液固耦合和液體晃動的影響。得出液固耦合作用下，配筋方式及預應力鋼筋間距對大型LNG混凝土外罐的梁式振動頻率的影響；外罐自振頻率跟泄露液體深度的曲線關系；對于LNG儲罐的優(yōu)化設計提供了設計依據(jù)，對地震作用下儲罐的動態(tài)分析提供數(shù)據(jù)，為研究LNG儲罐的抗震設計具有明顯的指導性意義。

LNG儲罐; 模態(tài)分析; 液固耦合; 液體泄漏

隨著國內經濟的快速發(fā)展以及人們環(huán)保意識的增強，天然氣作為一種清潔能源應用越來越廣泛，據(jù)國際權威機構預測，到2020年，天然氣將超過石油成為能源結構中的第一組成部分。LNG預應力混凝土儲罐是國際上廣泛使用的儲存液化天然氣的容器。由于地震活動比較頻繁，LNG儲罐在地震作用下一旦發(fā)生破壞，就會產生一系列的次生災害，如火災、爆炸等，做好儲罐的抗震工作尤為重要。

山東省青島市LNG儲罐液位于膠南市的董家口港，是 國家“十一五”液化天然氣規(guī)劃的重點項目，規(guī)劃建設3座16萬m3儲罐，所以根據(jù)青島市董家口具體的地質條件以及相應的抗震烈度，對LNG儲罐進行模態(tài)分析是非常必要的。本文結合實際，考慮液固耦合作用分析預應力混凝土儲罐內外罐的模態(tài)，液固耦合模態(tài)分析是結構反應譜分析的基礎，且儲罐液固耦合的模態(tài)分析能為地震作用下儲罐的動態(tài)分析提供精確的數(shù)據(jù)，因此對地震作用下的LNG儲罐進行液固耦合模態(tài)分析意義重大。

1 儲罐有限元模型的建立

以山東省青島市董家口16萬m3大型預應力混凝土儲罐為例，建立分析模型，在模型的建立過程中，單位取用國際單位制。儲罐采用直角坐標系，外罐采用三維實體單元，內罐采用等參單元，液體采用三維流體單元[1,2]。其模型參數(shù)如表1所示。

將大型 LNG 儲罐底部視為固定端，全方位自由度約束，即三個平動自由度和三個轉動自由度，不考慮地基對罐體動力反應的影響；不考慮罐頂?shù)恼駝?，將罐頂豎向振動約束[3,4]；考慮重力作用 。LNG儲罐有限元模型如圖1所示。

(1)設立資產清算組。公安機關對非法集資行為立案偵查后，如果發(fā)現(xiàn)涉案企業(yè)或個人有資不抵債的情況，可以報請法院以裁定的方式委托“防非組”組建非法集資案件資產清算組。該清算組隸屬于“防非組”，成員包括“防非組”的全部或部分人員以及一定比例的律師、會計師、集資參與人的代表等。該清算組行使類似于資產管理人的職權，參照我國《破產法》規(guī)定的清算程序，負責非法集資犯罪涉案財物的保管、清理、估價和分配。該清算組具有民事主體資格，可以依法進行必要的民事活動。

表1 模型參數(shù)Table 1 Model parameter

圖1 LNG儲罐有限元模型Fig.1 LNG storage tank finite model

2 儲罐模態(tài)分析方法的理論計算值

4.3 內罐的模態(tài)分析

對于大型LNG內儲罐的罐計算，主要應用《鋼制常壓立式圓筒形儲罐抗震鑒定標（SHT3026-2005）》[6]內。按照規(guī)范，罐的固液耦連基本自振頻率：

R為內罐半徑； ρL為LNG液體密度； E為內罐材料的彈性模量； δ3為位于罐壁1/3處的罐壁名義厚度,m）。

內罐固液耦連基本自振頻率為2.398 Hz；周期為0.417 s。

1.2.3 脫落、剔除或終止標準 因不良反應和(或)合并癥病情惡化不能完成治療者，因療效不佳而中斷或未按要求配合完成評分及療效評價者。

其中， cth為雙曲余切函數(shù)； g為重力加速。

內罐液體晃動的基本頻率為0.1027Hz；周期為9.736 3 s.

3 大型LNG儲罐內外罐動力特性計算結果

在正常情況下，大型LNG內儲罐的罐和外罐的模態(tài)分析是相互影響的。下表為儲罐考慮液固耦合內罐的頻率與液體晃動頻率以及內外罐頻率相互影響結果。

通過表2可以看出，內外罐液固耦合頻率，晃動頻率和外罐頻率與理論分析的結果相吻合，內罐和外罐單獨計算的頻率要大于整體計算的結果，說明質量的增加會降低頻率；外罐對內罐液固耦合的頻率影響較大，對內罐液體晃動幾乎無影響。

丁振軍表示，盡管近年來市場經歷著一些改變，但是歐化農業(yè)堅持不變的是產品的品牌、生產線、生產團隊、生產控制，特別是產品的品質不會改變。

表2 儲罐內外動力特性計算結果Table 2 Dynamic characteristics results of inner and outer tank

4 大型LNG儲罐的模態(tài)分析

4.1 外罐的模態(tài)分析

全介質自承式光纜(即ADSS光纜)的張力承載元件主要為紡綸纖維，具備造價低、安全性高、易維護等優(yōu)點，施工和運維與電網關聯(lián)度小，可在輸電線路帶電情況下進行操作。

通過表3的分析，不同的配筋方式和有無預應力鋼筋對大型LNG混凝土外罐的梁式振動頻率的影響很小，其相對誤差不超過0.01%，對于預應力筋的優(yōu)化和設計有非常重要的指導意義。并且理論計算值與有限元計算值得相對誤差不超過4%，說明有限元計算的準確性。

表3 外罐前三階梁式振動Table 3 The first three beam vibration of outer tank

大型LNG儲罐外罐的結構及邊界條件都是對稱的，在模態(tài)分析中會出現(xiàn)振型相似，但相位不同的情況，同時大型LNG儲罐的外罐這種殼體結構，在其環(huán)向極易產生多波環(huán)向效應。LNG儲罐內罐液固耦合和外罐的振動形式都是以 cosnθ型梁式振動為主[7]，且在抗震設計中只是需要梁式震動，即中n=1的情況。外罐的前三階振型如表3所示，分析不同的配筋方式和有無預應力鋼筋對大型LNG混凝土外罐的梁式振動頻率的影響。

在傳統(tǒng)教學中英語教師在施教過程中往往采用填鴨式的教學，教師關注學生的筆試多一點而輕視學生的聽說能力培養(yǎng)。在新課程理念的指引下，小學英語教師應試著從傳統(tǒng)的教學形式中解脫出來，改進教學方式和方法，提高學生的綜合素質。但這些愿望的實現(xiàn)需要相關的政策支持與指導。

液體泄漏的情況下，將動液壓力等效為與外墻具有相同加速度的附加質量[8,9]，不考慮表面重力波的影響[10]。圖3為內罐液體泄露條件下不考慮內罐的影響，大型LNG儲罐的混凝土外罐梁式振動的第一階頻率與泄露液體深度的變化曲線。其中h/H為泄露液體與最大泄露液體高度的比值。由圖2可以看出，外罐自振頻率隨泄露液體深度的增加而減小，當內罐液體完全泄露時，液體對大型LNG儲罐的外罐的頻率影響最大。當泄露液體深度低于35%時，頻率變化較慢；當超過35%后，頻率隨泄露液體高度的增加下降較快。所以當液體發(fā)生泄漏時，外罐的頻率增大，周期減小，如果此時發(fā)生地震作用外罐的地震反應可能會增大，加上液體泄漏，可能會造成嚴重的后果。所以設計時應嚴格控制液體的泄漏，避免災害的發(fā)生。

4.2 液體泄漏時外罐的頻率

對于大型LNG混凝土外罐的計算，主要基于鄧克利類型逼近的方法[5]， 由此理論，經計算得到大型LNG混凝土外罐的自振頻率為6.905 5 Hz；周期為0.144 8 s。

按照規(guī)范，內罐的晃動基本自振周期：

大型LNG儲罐的內罐的振動形式也是以多波環(huán)向效應為主。當內罐出現(xiàn)cosnθ（n≥2） 的環(huán)向振型時，由于儲液結構物的振動而產生的動液壓力基本上是在cosθ 的壓力分布基礎上迭加cosnθ（n≥2）的壓力值，并且 cosθ（即n=1）的壓力值占主要成分。所以梁式振動在后期的時程分析及動液壓力分析中有極其重要的意義，下表為內罐的前三階頻率。

圖 2 外罐固有頻率與泄露液體深度關系曲線Fig.2 Curve relationship between natural frequency and leaking depth of liquid

通過表4可以得出，液固耦合振動頻率與內罐液面晃動頻率的有限元計算值與理論值相差不多，說明有限元計算的正確性和模型的正確性。

主要采用BA樓控控制的方式，通過硬接點的形式接入DDC模塊實現(xiàn)，要了解風機盤管控制原理。風機盤管的電氣設計情況（見圖2）不同于通常的動力設備，風機盤管的功率較小，一般不單獨設置控制箱，電氣設計類似照明設計，單獨設置回路，且同時并接多個風機盤管設備，設置的數(shù)量與回路功率有關，只要保持功率在合理范圍內即可。

表4 內罐前三階梁式振動Table 4 The first three beam vibration of inner tank

圖3 儲液深度與內罐的第一振型頻率影響曲線Fig.3 Curve relationship between the first mode frequency of the inner tank and the liquid storage depth

在正常情況下，內罐最高設計液面高度（相對于內罐底面）為34.73 m，正常操作液面高度為34.26 m，最低報警液面高度為2.3 m。儲液深度對內罐的第一振型頻率影響曲線如圖3所示，其中h/H為儲液深度體與最高的設計液面高度的比值，從圖中可以看出，儲液深度越大，第一階頻率越小。

4 結 論

本文通過ADINA有限元軟件模擬，得到山東省青島市董家口LNG接收站16萬m3LNG儲罐的液固耦合的模態(tài)，并對其進行了分析，得出以下結論：

（1）不同的配筋方式和有無預應力鋼筋對大型LNG混凝土外罐的梁式振動頻率的影響很小，并且理論計算值與有限元計算值得相對誤差不超過4%，說明有限元計算的準確性。

2.1.3 金芪降糖片供試品溶液的制備 精密稱取金芪降糖片粉末0.5 g，甲醇定容至50 mL容量瓶，超聲提取1 h，稱質量，補足質量，進樣前0.45 μm微孔濾膜過濾。

（2）外罐自振頻率隨泄露液體深度的增加而減小，當內罐液體完全泄露時，液體對大型LNG儲罐的外罐的頻率影響最大。

（3）內罐和外罐單獨計算的頻率要大于整體計算的結果，說明質量的增加會降低頻率；外罐對內罐液固耦合的頻率影響較大，對內罐液體晃動幾乎無影響。

以上結論為青島市LNG接收站儲罐的優(yōu)化設計提供了依據(jù)，為儲罐的抗震設計提供精確的數(shù)據(jù)，對于儲罐的抗震設計有明顯的指導意義。

［1］馮衛(wèi)民，宋立，肖光宇.基于AD1N A的壓力管道流固耦合分析[J].武漢大學學報，2009,2( 2) : 264-267.

［2］陳建平，周儒榮，韋忠煊.用結構有限元程序進行液體晃動動力分析[J].中國機械工程，2003,14(19):1631-1633.

［3］王翠翠.考慮液固荊合儲液罐非線性地震反應分析[D].哈爾濱:中國地震局工程力學研究所，2008.

［4］袁朝慶，孫建剛，徐松.儲罐抗震計算模型的簡化[J].壓力容器,2008,25(4) :16-20.

［5］Zhang Ruifu, Weng Dagen and Ren Xiaosong，Seismic analysis of a LNG storage tank isolated by a multiple friction pendulum system[J].EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION，2011,10(2).

［6］SHT3026-2005鋼質常壓立式圓筒形儲罐抗震鑒定標準[S].2005.

［7］張營,孫建剛,崔利富等.16萬方LNG儲罐的動特性分析［J］．大連民族學院院報,2010,12(5):464-466.

［8］劉傳卿，圓柱形消化池靜動態(tài)數(shù)值分析[D]，天津:天津大學,2005.

［9］HOUSNER G W. The dynamic behavior of water tanks [J]. Bulletin of the Seismological Society of America, 1963, 53(2): 381-387.

［10］FISHER D. Dynamic fluid effects in liquid-lled flexible cylindrical tanks [J]. Earthquake Engineering and Structural Dynamics,1979, 7:587-601.

Liquid-Solid Coupling Modal Analysis of LNG Storage Tanks Under Earthquake

GUAN You-hai1, JIA Juan-juan1, LIN Nan2, SU Ning1
（1. China University of Petroleum (Huadong), Shandong Qingdao 266580; 2. China Railway 16th Bureau Group Co.,Ltd., Beijing 100018）

Taking a 1.6×105m3large prestressed concrete tank in Dongjiakou, Qingdao, Shandong as the research object, Lanczos eigenvalue algorithm in ADIAN was used to carry out modal analysis of inner and outer tank of LNG storage tanks considering the liquid-solid coupling effect. The calculated result was compared with the result calculated by finite element method. The correctness and effectiveness of Lanczos eigenvalue algorithm in ADIAN were proved. Effects of outer tank on inner tank liquid-solid coupling and liquid sloshing were analyzed. The results show that, under liquid-solid coupling, reinforcement scheme and spacing of the prestressed tendon have effect on the beam vibration frequency of large LNG outer tank; outer tank natural vibration frequency has curve relationship with leaking liquid depth. This paper provides the design basis for the optimization design of LNG storage tank and data for dynamic analysis of the LNG storage tanks under earthquake, and has the guiding significance to the anti-seismic design of LNG storage tanks.

LNG storage tanks; Modal analysis; Liquid-solid coupling; Liquid leakage

TQ 051

A

1671-0460（2015）01-0148-03

國家自然科學基金項目（51308545）和山東省自然科學基金項目（ZR2012EEL23）共同資助。

2014-08-05

管友海（1975-），男，山東臨沂人，副教授，博士，2006年畢業(yè)于中國海洋大學環(huán)境工程專業(yè)，研究方向：結構抗震與城市防災信息化。E-mail：guanyouhai@upc.edu.cn。

賈娟娟（1990-），女，研究方向：結構抗震與城市防災信息化。E-mail：jiajjdemail@163.com。