FPSO立管平臺簡化疲勞分析及軟件開發(fā)

嚴鋒 陳正豪 黃渙青

摘? ? 要：FPSO在服役期間，其立管平臺會受到交變應(yīng)力而引起疲勞破壞。常用的疲勞分析方法，分為簡化疲勞方法、譜分析方法及斷裂力學(xué)分析方法。本文借助Sesam軟件，對某FPSO立管平臺進行有限元分析，并根據(jù)ABS規(guī)范和復(fù)合辛普森數(shù)值算法，編寫了一款簡化疲勞計算軟件并應(yīng)用于該實際工程，對有限元計算結(jié)果進行簡化疲勞評估，后期驗證了計算的正確性與精確性。

關(guān)鍵詞：FPSO立管平臺;簡化疲勞分析;復(fù)合辛普森數(shù)值算法

中圖分類號：U661.4? ??? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標識碼：A

Simplified Fatigue Assessment of FPSO Riser Platform

and Software Development

YAN Feng， CHEN Zhenghao， HUANG Huanqing

（ Guangzhou Marine Engineering Corporation， Guangzhou 510250 ）

Abstract： The Fatigue failure will occur because of alternating stress that are produced by the variable loads imposed on the riser platform of FPSO. Commonly used fatigue assessment methods are divided into simplified fatigue assessment method， spectral assessment method and the fracture mechanics assessment method. The finite element analysis （FEA） of a FPSO riser platform is carried out with Sesam software to obtain the fatigue stress at each critical location. According to ABSs rules and composite Simpson formula， a simplified fatigue calculation software is written and applied to the fatigue assessment of the FEA results， and the correctness and accuracy of the calculation is verified in the later stage.

Key words： FPSO Riser Platform; Simplified fatigue assessment; Composite Simpson formula

1? ? ?前言

海洋結(jié)構(gòu)物受到交變循環(huán)載荷時，會產(chǎn)生交變的應(yīng)力，進而引起疲勞破壞。FPSO在位工作時，受到波浪、流或設(shè)備周期性載荷的作用，因此需要對FPSO立管平臺進行疲勞評估。

疲勞評估主要有三種方式：簡化疲勞分析法：譜分析法;基于斷裂力學(xué)的疲勞分析法。

簡化疲勞分析法作為一種間接疲勞估計方法，常用于結(jié)構(gòu)節(jié)點疲勞強度的篩選。當(dāng)通過簡化疲勞分析法計算得到的結(jié)果滿足要求時，無需進行更深入一步的分析;如果計算結(jié)果不滿足要求時，則需要采用其它的方法進行評估。

2? ? ?簡化疲勞法的基本理論

海洋工程疲勞損傷，通常使用線性累積理論來進行預(yù)報，常用的有Palmgren-Miner累積損傷理論[2]，其認為每一次循環(huán)載荷損傷為獨立事件，總的損傷可以寫為：

（1）

式中：D為累積損失度;J為不同范圍的應(yīng)力數(shù);ni為第i個應(yīng)力范圍Si的循環(huán)次數(shù);Ni為在應(yīng)力范圍Si下直到失效時的應(yīng)力循環(huán)次數(shù)，可以由S-N曲線確定。

2.1? ?S-N曲線

許多船級社如ABS、DNV、API通過對標準構(gòu)件進行疲勞試驗，得到某種材料構(gòu)件應(yīng)力幅值與循環(huán)次數(shù)的關(guān)系曲線，稱為S-N曲線。

ABS規(guī)范根據(jù)構(gòu)件的樣式、連接型式和所處的環(huán)境，給定了相應(yīng)的S-N曲線，并根據(jù)節(jié)點焊接形式和連接方式的不同，指定了S-N曲線適用范圍。

應(yīng)力范圍S與疲勞破壞次數(shù)N的關(guān)系，如下式所示：

（2）

式中：A、m為疲勞強度系數(shù)，由疲勞試驗獲得。

一組典型的雙直線型S-N曲線，如圖1所示。

2.2? ?簡化疲勞損傷計算

ABS規(guī)范定義了相應(yīng)累積損傷公式，對于單直線S-N曲線，累積疲勞表達式為：

（3）

對于雙直線型S-N曲線，累積疲勞表達式為：

（4）

式中：D為損傷值;A、C、m、r為S-N曲線參數(shù);γ為Weibull分布形狀參數(shù);δ為Weibull尺度分布;NT為設(shè)計壽命內(nèi)應(yīng)力循環(huán)次數(shù)。一般設(shè)計壽命為25年時，取NT=1E8;Γ、Γ0為不完全伽馬函數(shù)。

3? ? 不完全Gamma函數(shù)求法

對雙直線型S-N曲線，不完全Gamma函數(shù)可以簡化為：

（5）

（6）

不完全Gamma函數(shù)通常是無法直接積分求得的，對于如何求解不完全Gamma函數(shù)，文獻[3]提出了眾多的方法：切比雪夫法、辛普森法、泰勒級數(shù)近似、遞推法、非線性回歸擬合參數(shù)化方法。本文采用復(fù)合辛普森法計算不完全Gamma函數(shù)積分值。

3.1? ?復(fù)合辛普森法

復(fù)合辛普森法[4]將定積分區(qū)間[a，b]分成n等分，每一份稱為一子區(qū)間，其步長為? ? ? ? ? ? ，分點為XK=a+kh，k=0，1，2，…n。

根據(jù)定積分的性質(zhì)，區(qū)間[a，b]上的積分值等于每一個子區(qū)間積分值之和，即 （7）

對每一子區(qū)間[XK，XK+1]使用辛普森法求得積分值，記區(qū)間中點為XK+1/2，區(qū)間[a，b]的積分值可以用下式來計算。

（8）

當(dāng)f（x）在區(qū)間[a，b]上有連續(xù)四階導(dǎo)數(shù)時，[a，b]區(qū)間的復(fù)合辛普森公式余項為：

（9）

復(fù)合辛普森法可以在每個小區(qū)間上采用數(shù)值積分求積公式進行計算，當(dāng)子間區(qū)足夠多時，求積公式誤差非常小，數(shù)值精度可以由余項RS來控制。

3.2? ?不完全Gamma函數(shù)的數(shù)值解

本文利用Python語言，在軟件層面實現(xiàn)了不完全Gamma函數(shù)的數(shù)值解法?！禖CS-船體結(jié)構(gòu)疲勞強度指南（2007）》[5]附錄1 Gamma函數(shù)表給出了一些不完全Gamma函數(shù)的近似解，對比其函數(shù)表數(shù)據(jù)與自編不完全Gamma函數(shù)計算程序結(jié)果，可以驗證軟件的精確度。

對比結(jié)果表明：最大誤差不超過0.065%。誤差原因主要是因為數(shù)據(jù)截斷造成的誤差，由此可證明自編軟件計算不完全Gamma函數(shù)是可靠而準確的。

4? ? 簡化疲勞分析程序

根據(jù)前面所述簡化疲勞法的基本原理，基于ABS規(guī)范，利用Python語言開發(fā)了一款簡化疲勞分析軟件。該軟件內(nèi)置了ABS規(guī)范的S-N曲線數(shù)據(jù)，只需要輸入幾個基本參數(shù)，就可以快速的算出熱點疲勞損傷度。

該軟件可以利用拉格郎日插值法快速計算熱點應(yīng)力，只需提取有限元分析的4個單元中心應(yīng)力結(jié)果和位置，就可以計算出熱點應(yīng)力。

在Design Parameters中輸入設(shè)計年限、修正板厚、設(shè)計安全因子等參數(shù);在Fatigue Assessment Parameters中選用S-N曲線、Weibul形狀參數(shù)及參考循環(huán)次數(shù)，點擊計算即可完成單直線損傷、雙直線損傷及許用應(yīng)力的計算。

4.1? ?算法比較下

（1）傳統(tǒng)手工計算

提取有限元結(jié)果→插值計算熱點應(yīng)力→查ABS規(guī)范得出S-N曲線參數(shù)→進行板厚修正→查表

函數(shù)插值點附近值→進行雙線型插值計算

函數(shù)近似值→計算損傷值。

傳統(tǒng)手工計算，每節(jié)點計算一次都需要查一次S-N曲線表、查4次Gamma函數(shù)表并進行二次雙線性插值計算。對許用應(yīng)力計算因涉及到循環(huán)算法，因此無法求出，只能計算疲勞損傷。

（2）自編軟件計算

提取有限元結(jié)果→自動計算熱點應(yīng)力→輸入計算參數(shù)/選用S-N曲線→計算損傷值→采用循環(huán)算法計算許用應(yīng)力值[1]。

使用軟件計算，無需查表和進行插值計算，效率大幅度提高，特別是節(jié)點越多工況越多時，效率提高越明顯，并且可以利用循環(huán)算法計算許用應(yīng)力。

4.2? 工程驗證

對某FPSO直升機甲板平臺，分別使用傳統(tǒng)手工計算和軟件計算方法，對其中某些節(jié)點進行了計算并作對比，表明兩種結(jié)果是一致的。

5? ? 立管平臺簡化疲勞分析

基于以上理論，采用DNV的Sesam軟件，對立管受力最大工況進行有限元分析，并根據(jù)ABS指南對可能發(fā)生破壞的節(jié)點進行疲勞分析。

5.1? ?計算模型

模型采用直角坐標系：坐標原點為FPSO尾垂線與基線交點;X方向為船長方向，從船尾指向船首;Y方向為船寬方向，從右舷指向左舷;Z方向為垂直方向，從基線指向主甲板。

該計算模型如圖2所示，包括：立管平臺上部、下部和部分FPSO船體。主船體模型縱向范圍自FR67-FR92，橫向范圍自距中25.48（m）縱艙壁至船體外板部份，垂向范圍從船底板至主甲板。

對于疲勞評估點，ABS規(guī)范規(guī)定必須細化為txt（t為板厚）大小的網(wǎng)格。模型采用板梁模型，板采用4節(jié)點單元模型;除細化區(qū)、加強筋、縱骨等采用梁單元模擬，細化區(qū)內(nèi)縱骨采用板單元模擬。

主船體部分采用一個縱骨間距劃分網(wǎng)格，即0.975 m;上下立管采用0.5 m劃分網(wǎng)格;對疲勞評估區(qū)，采用細化網(wǎng)格模型。在疲勞熱點附近采用txt大小劃分網(wǎng)格，并均勻過渡。典型疲勞評估模型，如圖3所示。

5.2? ?立管載荷

立管載荷采用OrcaFlex軟件進行計算。

根據(jù)立管系統(tǒng)設(shè)計方案，結(jié)合立管系統(tǒng)受到的環(huán)境載荷，對立管系統(tǒng)與FPSO之間的相互作用進行動態(tài)耦合分析，可得到立管頂端節(jié)點的張力時歷曲線。

整個模型由FPSO船舶模型和55根立管組成，其中立管采用深海鋼質(zhì)懶波立管SLWR。計算模型如圖4所示。

水動力計算時，波浪方向取0～360 ?，間隔為30 ?。

計算結(jié)果表明：在橫浪時FPSO的響應(yīng)幅度最大。對橫浪1年一遇、100年一遇的海況進行分析，得出各個立管頂端節(jié)點張力時歷曲線。圖5～圖8為橫浪時立管41的載荷。

立管最大載荷，如表1所示。

5.3? ?船舶運動加速度

FPSO立管平臺在位工作時，主要受到船體運動產(chǎn)生的加速度載荷和立管受到的交變載荷。

加速度由AQWA水動力計算得出，僅計算橫浪耦合情況，此時加速度幅值最大，如表2所列。

5.4? ?評估結(jié)果

對圖3所示節(jié)點進行疲勞評估。因篇幅所限，這里僅列出滿載生存計算結(jié)果：

對熱點處進行應(yīng)力外推，使用自編軟件進行拉格郎日插值，插值應(yīng)力為97.7 Mpa;使用自編軟件進行損傷計算，節(jié)點損傷為0.168，小于1，因此可認為疲勞滿足使用要求;許用應(yīng)力值為320 MPa。

6? ? 結(jié)論

設(shè)計時必須考慮FPSO立管平臺的疲勞性能。由于結(jié)構(gòu)復(fù)雜，節(jié)點眾多，疲勞計算非常繁雜。本文根據(jù)了ABS海工規(guī)范，編寫了相應(yīng)的簡化疲勞分析軟件，并應(yīng)用在實際工程項目中，取得了良好的結(jié)果。后期對軟件的正確性與精確性進行了論證，證明本疲勞計算軟件計算結(jié)果是精確的。

參考文獻

[1] ABS. Guide for the fatigue assessment of offshore structures[M]. American?Bureau of Shipping.2010.

[2] 曾春華，鄒十踐編譯.疲勞分析方法及應(yīng)用[M].北京：國防工業(yè)出版社，1991.

[3] 劉穎.不完全伽瑪函數(shù)計算算法改進與應(yīng)用[D].電子科技大學(xué)， 2015.

[4] 馬東升，董寧.數(shù)值計算方法[M].北京：機械工業(yè)出版社， 2015.

[5] CCS.船體結(jié)構(gòu)疲勞強度指南[M].上海：上海規(guī)范研究所， 2007.