劉婷婷，阮詩倫*，尹江洲，岳前進，謝 彬

(1. 大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室，遼寧 大連 116024；2. 中海油研究總院，北京 100027)

浮式液化天然氣液貨圍護系統(tǒng)的失效模式分析

劉婷婷1，阮詩倫1*，尹江洲1，岳前進1，謝 彬2

(1. 大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室，遼寧 大連 116024；2. 中海油研究總院，北京 100027)

針對浮式液化天然氣生產(chǎn)儲卸裝置(FLNG)中液貨圍護系統(tǒng)因晃蕩引起的結(jié)構(gòu)安全問題，以主流艙型GTT NO96薄膜型液艙圍護系統(tǒng)為對象，分別建立了圍護系統(tǒng)截面整體和局部有限元模型。采用有限元數(shù)值仿真方法系統(tǒng)地分析了準(zhǔn)靜力加載下液貨圍護系統(tǒng)可能的失效模式，得到了不同失效模式下的危險區(qū)域及對應(yīng)的極限荷載，為FLNG液艙圍護系統(tǒng)結(jié)構(gòu)安全設(shè)計提供參考。結(jié)果表明：當(dāng)前設(shè)計中的GTT NO96型液艙的主要失效模式是發(fā)生在次層木箱頂板與兩層木箱加強擋板交叉點處的強度失效，失效載荷為0.38 MPa。

圍護系統(tǒng)； 失效模式； 浮式液化天然氣生產(chǎn)儲卸裝置； 晃蕩荷載

0 引 言

我國海洋天然氣資源十分豐富，但以深海氣田、邊際小氣田和低品位天然氣資源為主，因此為了節(jié)約開采成本、提高天然氣開采效率，需要發(fā)展浮式液化天然氣生產(chǎn)儲卸裝置(FLNG)。圍護系統(tǒng)是FLNG中低溫液化天然氣(LNG)儲存的主要裝置。要實現(xiàn)隔絕熱傳遞及防止LNG泄漏的功能，其結(jié)構(gòu)設(shè)計關(guān)系整個FLNG裝備安全與否。FLNG對圍護系統(tǒng)的特殊要求包括更大的艙容、可提供足夠甲板空間安裝液化工藝系統(tǒng)模塊以及適應(yīng)海上長期作業(yè)。現(xiàn)有FLNG設(shè)計所采用的圍護系統(tǒng)是法國Gaztransport & Technigaz(GTT)公司的GTT系列薄膜型液艙以及日本IHI公司的SPB獨立型液艙。從國內(nèi)工程建造和使用經(jīng)驗上看，GTT NO96薄膜型液艙圍護系統(tǒng)目前仍舊是FLNG國產(chǎn)化設(shè)計的最佳選擇之一。

由于形狀和艙壁設(shè)計原因，GTT薄膜型液艙內(nèi)的LNG晃蕩荷載非常顯著，在晃蕩沖擊荷載作用下，結(jié)構(gòu)可能發(fā)生破壞，使低溫LNG液體泄漏，從而導(dǎo)致整個FLNG的災(zāi)難性后果。因此對FLNG液艙圍護系統(tǒng)的結(jié)構(gòu)失效模式開展研究是有效進行結(jié)構(gòu)設(shè)計的重要前提。相關(guān)的船級社規(guī)范對LNG運輸船液艙圍護系統(tǒng)的失效模式進行了闡述并提出失效準(zhǔn)則和結(jié)構(gòu)分析方法[1-4]，但尚無FLNG液艙結(jié)構(gòu)設(shè)計規(guī)范。Hirotsugu等[5]分別測試了靜態(tài)和動態(tài)荷載作用下GTT NO96液艙的力學(xué)特性，為數(shù)值模擬提供了有價值的參考。Arswendy等[6-8]分別采用有限元數(shù)值模擬方法研究了GTT NO96型液艙的剪切失效特性和強度極限，并探討了數(shù)值分析方法。

目前已有的液艙圍護系統(tǒng)結(jié)構(gòu)研究缺乏從結(jié)構(gòu)設(shè)計角度出發(fā)、根據(jù)FLNG特點對失效模式進行的系統(tǒng)分析。圍護系統(tǒng)截面分析結(jié)果可以為液艙整體結(jié)構(gòu)設(shè)計提供材料特性和整體剛度參數(shù)；而有限元方法能夠快速提取和分析結(jié)構(gòu)設(shè)計變量，全面評估結(jié)構(gòu)特性，并大量節(jié)約設(shè)計成本，更好地對結(jié)構(gòu)外形和材料進行優(yōu)化。本文對GTT NO96薄膜型液艙圍護系統(tǒng)的結(jié)構(gòu)力學(xué)特性進行分析，采用有限元方法模擬木箱截面結(jié)構(gòu)靜力失效并對失效模式進行系統(tǒng)分析，以期為FLNG液艙圍護系統(tǒng)結(jié)構(gòu)的安全設(shè)計提供一定參考。

1 FLNG液艙圍護系統(tǒng)的結(jié)構(gòu)特點

薄膜型液艙是當(dāng)前LNG運輸船的主流艙型，在船型性能方面優(yōu)于MOSS型和SPB型。按隔熱材料和施工方式的不同，代表性的薄膜型液艙主要有GTT NO96、MarkIII以及CS1等，其中，GTT NO96液艙采用木箱隔熱，結(jié)構(gòu)如圖1所示。其主體結(jié)構(gòu)由兩層36%殷瓦鋼金屬層和兩個獨立的木箱隔熱層組成。殷瓦鋼構(gòu)成的內(nèi)壁結(jié)構(gòu)厚度僅為0.7 mm，因此圍護系統(tǒng)的主要承載結(jié)構(gòu)為木箱隔熱層。木箱的隔熱功能通過填充密實的珍珠巖粉末實現(xiàn)；采用特殊螺栓設(shè)計連接將木箱固定，次層隔熱木箱底板與內(nèi)部船殼通過樹脂繩達到平整度設(shè)計要求；木箱由不同厚度層合板釘裝而成，包括上頂板、下底板、側(cè)面封板以及用于加強木箱截面彎曲剛度的加肋擋板。層合板幾何參數(shù)設(shè)計關(guān)系到木箱截面結(jié)構(gòu)強度，在液艙整體設(shè)計中根據(jù)荷載分布在薄弱環(huán)節(jié)進行加強設(shè)計。殷瓦鋼采用舌型設(shè)計與木箱隔熱層連接，然后高標(biāo)準(zhǔn)焊接在隔熱木箱頂板防止LNG液體泄漏。在FLNG運行過程中，隔熱層中充滿了氮氣并裝有傳感器對溫度和壓力進行監(jiān)測。

圖1 GTT NO96薄膜型圍護系統(tǒng)結(jié)構(gòu)示意圖Fig. 1 Schematic drawing of GTT NO96 containment system structure

2 液艙圍護系統(tǒng)截面有限元模型

根據(jù)結(jié)構(gòu)特點，GTT NO96型圍護系統(tǒng)截面的有限元模型需要進行一定的簡化，以提高分析效率。其中，兩層殷瓦鋼對整體結(jié)構(gòu)法向剛度的影響可以忽略；樹脂繩作用是保證平整度，對木箱力學(xué)特性影響較小；珍珠巖粉末在靜力加載中對木箱極限強度沒有影響[5]。木箱結(jié)構(gòu)由不同厚度、不同尺寸的層合板組成，幾何參數(shù)可由文獻[3]獲得。目前FLNG液艙的木箱設(shè)計基本沿用的是LNG運輸船設(shè)計，其中層合板厚度、擋板高度和數(shù)量等參數(shù)設(shè)計決定了木箱結(jié)構(gòu)強度。GTT公司提出了四種木箱類型，分別是標(biāo)準(zhǔn)型、加強型、超級加強型和終極加強型，通過對液艙內(nèi)木箱排布的設(shè)計可以針對艙內(nèi)荷載作用分布調(diào)整液艙整體強度。本文數(shù)值計算所采用的木箱為標(biāo)準(zhǔn)型。

首先采用殼體單元建立木箱截面的整體有限元模型，如圖2所示。首層隔熱木箱與次層隔熱木箱之間可以細(xì)微運動，可用LINK 8桿單元進行連接模擬。對于層合板，采用SHELL163各向同性殼單元模擬整體木箱截面。該單元可以分析薄的及中等厚度的板殼結(jié)構(gòu)，支持線性、大變形和變厚度非線性分析，密度為710 kg/m3，彈性模量為7 000 MPa，泊松比為0.17。次層隔熱木箱直接固定在船體內(nèi)殼，在木箱底板法線方向施加位移約束，木箱四角均被螺栓直接固定在內(nèi)殼上，在底板的四角施加全約束；另外，根據(jù)實際情況，木箱截面?zhèn)冗吔绲氖讓幽鞠涞装搴痛螌幽鞠漤敯宀捎门c外圍封板相同的法向位移約束[9]。

圖2 GTT NO96型液艙木箱有限元模型和幾何尺寸Fig. 2 Finite element model and geometry for box structure segment of GTT NO96 containment system

由以上殼單元模型可以得到層合板彎曲應(yīng)力和擋板屈曲模擬，但無法直接反映沿板厚分布的剪切應(yīng)力以及兩層木箱擋板交叉點的壓應(yīng)力分布，不能全面預(yù)測包括剪切失效在內(nèi)的木箱結(jié)構(gòu)失效模式。采用實體單元的有限元模型就可以直接獲得剪應(yīng)力，反映木箱局部應(yīng)力響應(yīng)。為了避免交點處應(yīng)力集中，分別截取首層和次層木箱內(nèi)擋板間中點建立實體有限元模型，如圖3所示。

圖3 GTT NO96型液艙木箱局部有限元模型尺寸Fig. 3 Geometry of local finite element model of GTT NO96 containment system

在局部有限元模型中，采用SOLID185實體單元對木箱頂板、擋板和底板進行模擬，在局部模型的截面施加對稱邊界條件，其他約束條件與整體木箱模型相同。

隔熱木箱是GTT NO96型液艙的主要承載結(jié)構(gòu)，艙內(nèi)LNG晃蕩所產(chǎn)生的沖擊荷載直接作用于隔熱層并傳遞到船殼?；问帥_擊荷載具有很強的隨機性，通常采用模型試驗方法進行測量。為了降低非線性對相似率的影響，反映出沖擊荷載對液艙結(jié)構(gòu)的作用，壓力傳感器設(shè)計為可測量到1 m2面積內(nèi)平均壓力的點陣布局[2]，因此在木箱結(jié)構(gòu)分析中就可以認(rèn)為沖擊荷載均布于木箱頂板表面。

3 FLNG液艙圍護系統(tǒng)失效模式分析

3.1 液艙圍護系統(tǒng)失效模式概述

圍護系統(tǒng)結(jié)構(gòu)失效模式研究是進行結(jié)構(gòu)設(shè)計的重點，在其基礎(chǔ)上能夠更有針對性地進行液艙結(jié)構(gòu)設(shè)計。許多設(shè)計規(guī)范及報告都對GTT NO96型液艙可能發(fā)生的破壞位置和形式進行了闡述[1,3]。主要失效模式可以總結(jié)為：(1)木箱結(jié)構(gòu)強度失效，即木箱結(jié)構(gòu)應(yīng)力超過強度極限；(2)屈曲失效，可能發(fā)生在兩層木箱封板及擋板位置，是超過屈曲壓力而產(chǎn)生的結(jié)構(gòu)失穩(wěn)；(3)剪切失效，擋板對頂板及底板會產(chǎn)生剪切作用，同時擋板與封板的交叉點也會發(fā)生剪切失效；(4)殷瓦鋼失效，厚度為0.7 mm的殷瓦鋼緊貼于平整的木箱頂板表面，一旦頂板發(fā)生斷裂破壞，殷瓦鋼發(fā)生破損，將會導(dǎo)致低溫LNG泄漏，造成災(zāi)難性后果。

3.2 木箱強度失效分析

為了研究木箱結(jié)構(gòu)達到極限強度時的晃蕩荷載以及最先發(fā)生強度破壞的位置，在整體殼體單元木箱有限元模型頂部加載持續(xù)增加的準(zhǔn)靜態(tài)均布荷載，一直達到結(jié)構(gòu)強度極限，如圖4所示。層合板材料強度極限為20 MPa，此時的均布壓力為0.38 MPa。

圖4 達到強度極限時木箱的結(jié)構(gòu)應(yīng)力狀態(tài)Fig. 4 Von Mises stress distribution of box structure at extreme state

通過比較，應(yīng)力最大位置發(fā)生在次層木箱頂板與首層和次層加強擋板的交叉位置處，如圖5(a)所示，此時次層木箱頂板由于壓力產(chǎn)生結(jié)構(gòu)壓潰失效。同時，首層木箱頂板以及木箱外圍擋板也會形成較大應(yīng)力，可能發(fā)生彎曲失效。

圖5 強度極限時木箱層合板應(yīng)力狀態(tài)Fig. 5 Von Mises stress distribution of the primary and secondary boxes at extreme state

3.3 木箱屈曲失效分析

GTT NO96隔熱木箱采用了擋板進行結(jié)構(gòu)加強，因而也帶來了擋板及封板屈曲失效問題。采用特征值屈曲分析方法能夠初步預(yù)測發(fā)生靜態(tài)屈曲時結(jié)構(gòu)達到屈曲極限時的晃蕩荷載，如圖6所示。此時，屈曲臨界應(yīng)力為1.48 MPa，最大位移為9.98 mm，出現(xiàn)在截面次層木箱外層封板位置。美國船級社于2006年進行了大型LNG運輸船圍護結(jié)構(gòu)強度試驗[2]。試驗結(jié)果表明，當(dāng)箱型結(jié)構(gòu)承受垂向壓力時，下層木箱外圍封板首先發(fā)生屈曲失效，進而進入塑性變形并崩潰。該試驗?zāi)軌蝌炞C屈曲失效發(fā)生位置，但由于強度失效位置是木箱內(nèi)部層合板，難以進行觀察測量，因此強度失效及屈曲失效的定量分析還需進行深入研究和比較。

圖6 GTT NO96木箱結(jié)構(gòu)屈曲失效狀態(tài)Fig. 6 Bulking failure of GTT NO96 box structure

3.4 木箱剪切失效分析

如上所述，由于殼單元有限元模型的限制，需要采用實體單元構(gòu)建的木箱局部模型分析結(jié)構(gòu)剪切破壞狀態(tài)及具體位置。在木箱頂部施加持續(xù)增加的準(zhǔn)靜態(tài)均布荷載直到剪切應(yīng)力達到極限強度。當(dāng)均布壓力達到0.53 MPa時，圍護系統(tǒng)截面結(jié)構(gòu)達到許用剪切應(yīng)力2.8 MPa[2]。如圖7所示，首先達到許用剪切應(yīng)力的是次層木箱頂板。

圖7 達到剪切極限時結(jié)構(gòu)的剪應(yīng)力狀態(tài)Fig. 7 Shearing stress distribution of box structure at extreme state

3.5 殷瓦鋼失效分析

殷瓦鋼薄膜是液艙圍護系統(tǒng)的主要部分，構(gòu)成了裝載低溫LNG的完整界面，并與絕熱層一起保證船體結(jié)構(gòu)不受低溫侵害而產(chǎn)生脆性破壞[10]。殷瓦鋼具有良好的低溫沖擊韌性，難以在晃蕩荷載作用下直接發(fā)生破壞，但是殷瓦鋼僅0.7 mm厚，不能單獨承擔(dān)晃蕩荷載，因此判斷殷瓦鋼失效應(yīng)該以木箱頂板發(fā)生破壞為準(zhǔn)。

采用準(zhǔn)靜態(tài)均布壓力加載到整體木箱有限元模型。當(dāng)加載到1.49 MPa時，首層木箱的頂板會發(fā)生強度破壞，此時可能導(dǎo)致首層木箱的殷瓦鋼破壞并造成低溫LNG直接泄漏。

4 結(jié) 論

采用有限元數(shù)值分析方法對GTT NO96薄膜型液艙圍護系統(tǒng)木箱失效模式進行研究，得到如下結(jié)論：

(1) 采用準(zhǔn)靜力均布壓力作用于液艙截面，隨著壓力增加依次出現(xiàn)的失效模式為強度失效、剪切失效、屈曲失效以及殷瓦鋼薄膜失效，所對應(yīng)的木箱位置分別為次層木箱頂板兩層木箱擋板交叉點、次層木箱頂板、次層木箱外圍封板以及首層木箱殷瓦鋼。

(2) 強度失效是GTT NO96木箱結(jié)構(gòu)的首要失效模式，其中次層木箱頂板是整個結(jié)構(gòu)強度最薄弱環(huán)節(jié)，需要在FLNG圍護系統(tǒng)設(shè)計中特別注意。如不滿足FLNG晃蕩荷載要求，則可以采取對頂板加強、改變擋板形式及重新設(shè)計荷載傳遞路徑等措施。

(3) 本文對液艙圍護系統(tǒng)失效模式的研究結(jié)論主要基于準(zhǔn)靜力分析，實際情況中，在晃蕩荷載沖擊荷載作用下，圍護系統(tǒng)的結(jié)構(gòu)響應(yīng)還會產(chǎn)生動力放大效應(yīng)，因此還需要繼續(xù)深入進行動態(tài)沖擊荷載下的結(jié)構(gòu)失效模式研究。

[1] Det Norske Veritas. Classification notes No. 30.9. Sloshing analysis of LNG membrane tanks[S]. Oslo: Det Norske Veritas, 2006.

[2] American Bureau of Shipping. Guidance notes on strength assessment of membrane-type LNG containment systems under sloshing loads[S]. Houston: American Bureau of Shipping, 2009.

[3] Det Norske Veritas. Strength assessment of LNG membrane tanks under sloshing loads[S]. Oslo: Det Norske Veritas, 2011.

[4] Lloyd’s Register. Sloshing assessment guidance document for membrane tank LNG operations[S]. London: Lloyd’s Register, 2009.

[5] Dobashi H, Usami A. Structural response of the insulation box on NO96 membrane containment system [C]. 20th International Offshore and Polar Engineering Conference, 2010: 26.

[6] Arswendy A, Liasjoe O, Moan T. Comparative study of FE models of LNG containment system NO. 96 [C]. 20th International Offshore and Polar Engineering Conference, 2010: 62.

[7] 莊志鵬，劉 俊，唐文勇. 薄膜型LNG船晃蕩沖擊局部強度分析建模方法研究[J]. 船舶工程, 2011, 33(3): 17.

[8] 陸 曄，滕 蓓，祁恩榮. LNG船液艙圍護系統(tǒng)結(jié)構(gòu)極限承載力研究[J]. 艦船科學(xué)技術(shù), 2012, 34(2): 36.

[9] 滕 蓓，祁恩榮，陸 曄, 等. 大型LNG船液艙晃蕩結(jié)構(gòu)動響應(yīng)研究[J]. 船舶科學(xué)技術(shù)，2012, 34(4): 7.

[10] Grazyk M, Berget K, Allers J. Experimental investigation of invar edge effect in membrane LNG tanks[J]. J Offshore Mechanics and Arctic Engineering, 2012, 134(3): 031801.

FailureModesAnalysisofBoxStructureoftheContainmentSysteminFLNG

LIU Ting-ting1， RUAN Shi-lun1， YIN Jiang-zhou1， YUE Qian-jin1， XIE Bin2

(1.StateKeyLaboratoryofStructuralAnalysisforIndustrialEquipment，DalianUniversityofTechnology，Dalian,Liaoning116024，China；2.CNOOCResearchInstitute，Beijing100027，China)

The global and local finite element models of the mainstream containment system GTT NO96 of the floating liquid natural gas unit (FLNG) are established to evaluate the structural safety of tanks under the sloshing load. The failure modes of the models under quasi-static loading are systematically simulated by finite element method. The failure zones and extreme loads for different failure modes are obtained, which could provide a reference for the design of FLNG containment system. The results suggest that the main failure mode of the current design of FLNG container is the strength failure of the cover plate of the secondary insulation box at the intersection between the bulkheads of the primary and secondary boxes, and the extreme load is 0.38 MPa.

containment system； failure mode； floating liquid natural gas unit； sloshing load

TE88

A

2095-7297(2014)01-0050-05

2014-02-21

國家科技重大專項(2011ZX05026-006-06)、工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室自主研究課題(S13207)

劉婷婷(1981—)，女，博士后，主要從事海洋工程結(jié)構(gòu)分析和流體晃蕩沖擊方面的研究。

*通信作者