梁 凱，王亞瓊，馬 超，劉 天

中國(guó)石油天然氣管道工程有限公司，河北廊坊 065000

國(guó)內(nèi)外學(xué)者和工程設(shè)計(jì)人員對(duì)系泊錨鏈的疲勞特性進(jìn)行了大量研究：羅慧明等基于[1]時(shí)域耦合分析，采用雨流計(jì)數(shù)法對(duì)CALM系泊錨鏈的疲勞損傷進(jìn)行了計(jì)算；鄭長(zhǎng)新等[2]以某Spar平臺(tái)錨泊線(xiàn)為研究對(duì)象，基于非線(xiàn)性有限元和雨流計(jì)數(shù)法對(duì)其疲勞性能進(jìn)行了評(píng)估；喬?hào)|生等[3]對(duì)某Spar平臺(tái)的復(fù)合錨泊線(xiàn)的疲勞進(jìn)行了計(jì)算；王文明等[4]采用有限元軟件MSC-MARC對(duì)流花11號(hào)半潛式鉆井平臺(tái)錨鏈在極限環(huán)境載荷下服役10年后的疲勞壽命進(jìn)行了校核；康占賓[5]對(duì)錨鏈在不同平面外彎曲載荷下的累計(jì)損傷疲勞壽命進(jìn)行了分析；杜君峰等[6]對(duì)腐蝕影響下的深水平臺(tái)系泊錨鏈疲勞損傷進(jìn)行了評(píng)估，探討了腐蝕要素對(duì)系泊錨鏈疲勞損傷的影響規(guī)律；宋憲倉(cāng)等[7]研究了二階差頻力對(duì)系泊錨鏈疲勞損傷的影響。

經(jīng)過(guò)對(duì)比研究可以發(fā)現(xiàn)，現(xiàn)階段對(duì)錨鏈疲勞損傷的計(jì)算分析均采用基于雨流計(jì)數(shù)的時(shí)域不規(guī)則波法，而缺乏基于確定性規(guī)則波法的錨鏈疲勞評(píng)估研究和應(yīng)用，主要原因在于基于時(shí)域的不規(guī)則波法最為精確。然而，基于時(shí)域的不規(guī)則波疲勞分析耗時(shí)多，對(duì)計(jì)算機(jī)性能要求高，不能滿(mǎn)足項(xiàng)目初步設(shè)計(jì)時(shí)快速響應(yīng)的要求。因此，本文給出了把基于Longuet-Higgins分布的波浪散布圖離散為規(guī)則波的理論方法，最后以某CLAM系統(tǒng)為例，給出了基于確定性規(guī)則波法進(jìn)行錨鏈疲勞性能評(píng)估的一般流程。

1 疲勞計(jì)算基本理論

1.1 疲勞累計(jì)損傷計(jì)算方法

目前，結(jié)構(gòu)疲勞計(jì)算方法主要有確定性規(guī)則波分析法、確定性不規(guī)則波法和隨機(jī)譜分析法。確定性規(guī)則波分析方法是將波浪散布圖離散為多個(gè)規(guī)則波，分別計(jì)算每個(gè)規(guī)則波作用下的疲勞損傷，而后求和得到總疲勞損傷；確定性不規(guī)則波法是首先對(duì)結(jié)構(gòu)進(jìn)行時(shí)域內(nèi)的動(dòng)力分析，得到結(jié)構(gòu)應(yīng)力時(shí)程，然后采用某種應(yīng)力循環(huán)計(jì)數(shù)方法（如雨流計(jì)數(shù)法）進(jìn)行應(yīng)力循環(huán)計(jì)數(shù)，進(jìn)而計(jì)算結(jié)構(gòu)的疲勞損傷；隨機(jī)性的譜分析方法簡(jiǎn)稱(chēng)譜分析法，該方法也是首先對(duì)結(jié)構(gòu)進(jìn)行動(dòng)力分析，得到結(jié)構(gòu)的應(yīng)力時(shí)程，然后采用快速Fourier變換將應(yīng)力變換到頻域，得到應(yīng)力RAO和功率譜密度函數(shù)，進(jìn)而計(jì)算結(jié)構(gòu)的疲勞損傷。

在上述三種方法中，不規(guī)則波法被認(rèn)為是最為精確的方法，其充分考慮了波浪荷載的隨機(jī)特性及波浪與結(jié)構(gòu)耦合動(dòng)力分析的非線(xiàn)性特性，但該方法耗時(shí)最多，對(duì)計(jì)算機(jī)要求最高。譜分析法計(jì)算速度快，但基于線(xiàn)性假設(shè)，不能考慮荷載及系統(tǒng)非線(xiàn)性因素對(duì)結(jié)構(gòu)疲勞損傷的影響，計(jì)算精度難以保證。規(guī)則波法計(jì)算比不規(guī)則波法快得多，對(duì)計(jì)算機(jī)硬盤(pán)空間的需求也小得多。假如將波浪散布圖離散得較為理想，那么采用規(guī)則波法得到的計(jì)算結(jié)果將和不規(guī)則波法的結(jié)果相差不大。董磊磊等[8]研究表明，若將波浪譜離散的周期區(qū)間離散得更加精細(xì)，則由基于精細(xì)離散規(guī)則波法得到的疲勞壽命將和雨流計(jì)數(shù)法的結(jié)果更加接近。

1.2 Longuet-Higgins法

設(shè)定Zl和Dl分別為波浪散布圖的某短期海況和其相應(yīng)的持續(xù)時(shí)間，其中l(wèi)=1，2，…，k，k為短期海況數(shù)。那么，短期海況Zl出現(xiàn)的概率為：

式中：Dl不是某短期海況實(shí)際發(fā)生的時(shí)間，而是在波浪散布圖中相對(duì)的時(shí)間占比。

根據(jù) Barltrop和 Adams[9]的觀點(diǎn)，若采用Longuet-Higgins[10]分布，則波高為H、周期為T(mén)的規(guī)則波在給定的隨機(jī)海況中出現(xiàn)的聯(lián)合概率密度函數(shù)為：

式中：Z為隨機(jī)海況的波浪譜；Tm為波浪譜平均周期；Tz為波浪譜平均跨零周期；m0為波譜Z的0階譜距；m1為波譜Z的1階譜距；m2為波譜Z的2階譜距；Hs為隨機(jī)海況的有義波高。

用 Hi（i=1，2，…，m） 和 Tj（j=1，2，…，n）分別表示規(guī)則波的波高和周期，將波高和周期分別在區(qū)間內(nèi)的事件稱(chēng)為Bij，則某短期海況Zl中Bij發(fā)生的概率為：

式（3）的積分區(qū)間分別為：

由式（1）、（2）、（3） 可以得到，Bij出現(xiàn)的總概率為：

設(shè)定波浪散布圖發(fā)生的總時(shí)間為Dtotal，則Bij出現(xiàn)的次數(shù)Oij為：

2 算例分析

2.1 數(shù)值模型

本文以某30×104t油輪的CALM系泊系統(tǒng)為例，采用Longuet-Higgins法對(duì)工程區(qū)域的波浪散布圖進(jìn)行處理，得到規(guī)則波散布圖，采用該規(guī)則波散布圖對(duì)系泊系統(tǒng)的錨鏈進(jìn)行疲勞分析。計(jì)算采用OrcaFlex軟件完成。CALM系統(tǒng)的動(dòng)力耦合分析模型見(jiàn)圖1，系泊錨鏈布置見(jiàn)圖2。CALM系統(tǒng)浮筒高5 m，直徑11m，吃水2.5m，在操作工況下油輪通過(guò)兩根φ17in（1in=25.4mm）的尼龍纜連接到浮筒，整個(gè)系統(tǒng)通過(guò)6根等間距布置的懸鏈?zhǔn)藉^鏈固定在海底，錨鏈為無(wú)檔錨鏈，等級(jí)R3S，直徑87mm。

圖1 CALM系統(tǒng)動(dòng)態(tài)耦合分析模型

圖2 系泊錨鏈布置

2.2 張力-壽命(T-N)曲線(xiàn)

錨鏈疲勞累積損傷的計(jì)算是以P-M線(xiàn)性累積損傷模型和T-N曲線(xiàn)為基礎(chǔ)的。T-N曲線(xiàn)一般可表示為：

式中：N為循環(huán)次數(shù)；R為張力范圍與錨鏈參考破斷強(qiáng)度的比值；M為T(mén)-N曲線(xiàn)斜率；K為T(mén)-N曲線(xiàn)截距。本文依照API規(guī)范進(jìn)行錨鏈疲勞計(jì)算校核，對(duì)于無(wú)檔錨鏈，M取值3.0，K取值316。

2.3 環(huán)境條件

在進(jìn)行錨鏈疲勞強(qiáng)度校核時(shí)，本文參照的長(zhǎng)期海況波浪散布圖見(jiàn)表1。

將波浪散布圖中概率小于0.1%的短期海況排除（即概率為0、0.071%的海況），剩余28個(gè)短期海況。根據(jù)CALM系統(tǒng)作業(yè)要求，油輪裝卸作業(yè)時(shí)有義波高不高于3.5 m，因此除發(fā)生概率為0.185%和0.142%兩個(gè)短期海況外，有26個(gè)短期海況適于油輪連接作業(yè)。根據(jù)對(duì)油輪運(yùn)輸頻次統(tǒng)計(jì)分析，在26個(gè)適于油輪作業(yè)的海況下，油輪滿(mǎn)載作業(yè)、油輪壓載作業(yè)和無(wú)系泊油輪作業(yè)三種工況發(fā)生的概率分別為0.3、0.3和0.4。因此可確定三種工況各自包含的短期海況及其發(fā)生概率，每種工況下短期海況的實(shí)際發(fā)生概率為短期海況出現(xiàn)概率與作業(yè)工況發(fā)生概率的乘積。按照E、N、NE、S、SE和SW共6個(gè)方向進(jìn)行統(tǒng)計(jì)，每種工況包含的短期海況數(shù)量見(jiàn)表2。

分別對(duì)各工況各浪向下的短期海況組合采用Longuet-Higgins分布進(jìn)行離散，可得到油輪滿(mǎn)載工況下、油輪壓載工況下以及無(wú)系泊油輪工況下規(guī)則波散布圖各6個(gè)，共18個(gè)規(guī)則波散布圖。為消除分解過(guò)程中產(chǎn)生的不符合實(shí)際的周期較大的規(guī)則波，采用波浪譜平均周期的2倍作為周期閾值進(jìn)行控制，離散水平為波高間隔1 m、周期間隔1 s。限于文章篇幅，本文只列舉油輪壓載工況E方向離散后的規(guī)則波散布圖，見(jiàn)表3。

表1 波浪散布圖

表2 各工況各浪向短期海況數(shù)量

2.4 計(jì)算結(jié)果及校核

對(duì)系統(tǒng)的6根錨鏈進(jìn)行疲勞分析，采用各工況各浪向的規(guī)則波散布圖中的規(guī)則波分別進(jìn)行動(dòng)態(tài)系泊計(jì)算，得到各規(guī)則波下的錨鏈張力；根據(jù)各規(guī)則波出現(xiàn)概率與錨鏈張力數(shù)據(jù)，可得到各規(guī)則波下錨鏈的疲勞損傷；采用線(xiàn)性累積法對(duì)錨鏈疲勞損傷求和，得到錨鏈總的疲勞損傷。按照API規(guī)范要求進(jìn)行錨鏈疲勞壽命校核[11]，計(jì)算及校核結(jié)果見(jiàn)表4。

表3 壓載工況E方向離散后規(guī)則波散布圖

表4 疲勞壽命

由表4可以看出：第一，計(jì)算得到的錨鏈的疲勞壽命安全系數(shù)均大于許用安全系數(shù)3，滿(mǎn)足API規(guī)范要求。第二，根據(jù)波浪散布圖，波浪的主要作用方向在N、NE、E、SE方向，因此位于這些環(huán)境條件相反方向的錨鏈疲勞損傷較大，其他錨鏈疲勞損傷則會(huì)大大減小。

3 結(jié)論

本文以某30×104t油船的CALM系泊系統(tǒng)為研究對(duì)象，基于Longuet-Higgins分布對(duì)工程海域的波浪散布圖進(jìn)行規(guī)則波化離散，最后采用確定性規(guī)則波法對(duì)系泊錨鏈的疲勞累積損傷進(jìn)行了評(píng)估，得到以下結(jié)論：

（1）基于Longuet-Higgins分布對(duì)工程海域的波浪散布圖進(jìn)行了離散，得到了各工況各浪向下的規(guī)則波散布圖，提供了波浪散布圖離散為規(guī)則波的一般方法和流程。

（2）以離散的規(guī)則波波浪散布圖為基礎(chǔ)，采用確定性規(guī)則波法對(duì)CALM系統(tǒng)的錨鏈進(jìn)行了疲勞校核。該方法具有計(jì)算速度快，占用資源少的優(yōu)點(diǎn)，可用于系泊錨鏈的疲勞評(píng)估。

（3）對(duì)于采用多根均布錨鏈系泊的CALM系統(tǒng)，如果作業(yè)區(qū)域環(huán)境載荷方向性明顯，則在主環(huán)境作用相反方向的錨鏈疲勞損傷較大，其他錨鏈可降低錨鏈規(guī)格，從而降低工程造價(jià)。