基于雨流計(jì)數(shù)法的錨泊線疲勞強(qiáng)度評(píng)估

劉 娟，孫 政，劉敬喜

(華中科技大學(xué) 船舶與海洋工程學(xué)院，湖北 武漢 430074)

基于雨流計(jì)數(shù)法的錨泊線疲勞強(qiáng)度評(píng)估

劉 娟，孫 政，劉敬喜

(華中科技大學(xué) 船舶與海洋工程學(xué)院，湖北 武漢 430074)

深海 FPSO 一般由多根錨泊線組成的系泊系統(tǒng)進(jìn)行定位，為保證結(jié)構(gòu)的安全性，開(kāi)展錨泊線疲勞強(qiáng)度評(píng)估分析顯得至關(guān)重要。本文采用改進(jìn)的雨流計(jì)數(shù)法結(jié)合 T-N 曲線和 Miner 線性疲勞累計(jì)損傷理論對(duì)某 FPSO 錨泊線進(jìn)行疲勞性能分析。針對(duì)錨泊線的疲勞損傷，分析其疲勞損傷的不利位置，從而為 FPSO 系泊系統(tǒng)的維護(hù)提供重要依據(jù)。

雨流計(jì)數(shù)法；錨泊線；疲勞強(qiáng)度；張力時(shí)間歷程

0 引 言

深海 FPSO 長(zhǎng)期系泊于海上，在風(fēng)浪流等海洋環(huán)境載荷的作用下，不可避免地會(huì)產(chǎn)生往復(fù)運(yùn)動(dòng)，從而造成了其系泊系統(tǒng)長(zhǎng)期遭受交變載荷作用而產(chǎn)生疲勞損傷。一旦其系泊系統(tǒng)發(fā)生疲勞斷裂，將會(huì)造成嚴(yán)重的安全事故，因此對(duì) FPSO 的系泊系統(tǒng)進(jìn)行疲勞損傷評(píng)估顯得極為重要。

目前，國(guó)內(nèi)外針對(duì)深海裝備系泊系統(tǒng)的疲勞性能開(kāi)展了相關(guān)研究工作。胡毅對(duì) FPSO 系泊系統(tǒng)進(jìn)行研究分析了它的時(shí)域受力得到了錨泊線的張力時(shí)歷曲線。Luo 分析了時(shí)域范圍內(nèi)錨泊線動(dòng)力響應(yīng)的非高斯特性，并使用張力循環(huán)計(jì)數(shù)計(jì)算錨泊線的疲勞損傷。Larsen 和 Mathisen 基于 T-N 曲線對(duì)系泊系統(tǒng)的疲勞性能進(jìn)行研究，并推廣到半潛式平臺(tái)和轉(zhuǎn)塔式系泊系統(tǒng)的疲勞損傷研究。Gilberto Bruno等針對(duì) FPSO 系泊系統(tǒng)的疲勞性能分析采用了 Dirlik 方法、雨流計(jì)數(shù)法、窄帶譜分析法和寬帶譜修正法進(jìn)行計(jì)數(shù)，最后對(duì)其疲勞壽命進(jìn)行計(jì)算。

本文以作業(yè)于南海 S4 海域 FPSO 系泊系統(tǒng)的錨泊線為對(duì)象開(kāi)展疲勞性能分析。計(jì)算 64 個(gè)短期海況條件下 FPSO 系泊系統(tǒng)的動(dòng)力響應(yīng)，并采用雨流計(jì)數(shù)法對(duì)錨泊線張力時(shí)間歷程進(jìn)行統(tǒng)計(jì)獲得其疲勞載荷譜。結(jié)合 T-N 曲線并應(yīng)用 Miner 線性累積損傷理論對(duì)其 64 個(gè)海況條件下的疲勞損傷進(jìn)行累積，分別獲得錨泊線中鋼鏈與鋼纜在長(zhǎng)期海況下的疲勞損傷。針對(duì)錨泊線的疲勞損傷，分析其疲勞損傷的不利位置，為 FPSO 系泊系統(tǒng)的維護(hù)提供重要依據(jù)。

1 疲勞強(qiáng)度評(píng)估步驟及原理

1.1 疲勞強(qiáng)度評(píng)估步驟

為了精確評(píng)估錨泊線的疲勞損傷，采用雨流計(jì)數(shù)法，對(duì)錨泊線張力時(shí)間歷程進(jìn)行統(tǒng)計(jì)分析獲得其疲勞載荷譜。然后結(jié)合 T-N 曲線并應(yīng)用 Miner 線性累積損傷理論進(jìn)行疲勞損傷評(píng)估[1–7]，整個(gè)評(píng)估過(guò)程如圖 1所示。

圖 1 錨泊線疲勞損傷評(píng)估過(guò)程Fig. 1 The evaluation process of fatigue strength

1.2 雨流計(jì)數(shù)法程序的實(shí)現(xiàn)方法

本文采用雨流計(jì)數(shù)法[8]對(duì)纜線張力時(shí)歷曲線進(jìn)行統(tǒng)計(jì)，提取張力范圍大小及循環(huán)次數(shù)。雨流計(jì)數(shù)法的計(jì)算機(jī)程序采用四峰谷值比較法來(lái)實(shí)現(xiàn)，主要步驟如下：

1）峰谷值點(diǎn)的提取

雨流計(jì)數(shù)法處理的是峰谷值依次出現(xiàn)的數(shù)據(jù)，因此得到錨泊線張力時(shí)歷數(shù)據(jù)結(jié)果后要提取其峰谷值。對(duì)張力時(shí)間歷程的峰谷值點(diǎn)提取分為 2 步：第 1 步是將數(shù)據(jù)中連續(xù)等值點(diǎn)進(jìn)行去除，保留其中一個(gè)點(diǎn)；第2 步是將剩下數(shù)據(jù)點(diǎn)的非峰谷值點(diǎn)進(jìn)行去除。對(duì)于某點(diǎn)是否屬于峰谷值點(diǎn)的判定，可以采用作差乘積法來(lái)實(shí)現(xiàn)。即對(duì)去除連續(xù)等值點(diǎn)之后的非首位點(diǎn) Pi，若有

上式成立，則認(rèn)為此點(diǎn)為非峰谷值點(diǎn)，應(yīng)將其去除。對(duì)于首位點(diǎn)，直接認(rèn)為其為峰谷值點(diǎn)，不再進(jìn)行判定。

2）循環(huán)數(shù)的提取

針對(duì)峰谷值序列的典型段，采用四點(diǎn)法雨流計(jì)數(shù)原則進(jìn)行統(tǒng)計(jì)：記下雨滴流過(guò)的最大峰谷值，作為一個(gè)循環(huán)；從載荷歷程中刪除雨滴流過(guò)的部分，對(duì)剩余歷程段，重復(fù)雨流計(jì)數(shù)法，直至再無(wú)剩余歷程為止。

其條件如下：

① 升式。如果 a ＜ b，b ≤ d，c 〉 a，則記錄一個(gè)全循環(huán) bcb′，如圖 2 所示。得到的幅值為 Sa=|b–c|/2，平均值為

圖 2 波形 1Fig. 2 Up waveform

② 降式。如果 a 〉 b，b ≥ d，c ≤ a，則記錄一個(gè)全循環(huán) bcb′，如圖 3 所示。得到的幅值為 Sa=|b–c|/2，平均值為

圖 3 波形 2Fig. 3 Down waveform

重復(fù)上述方法計(jì)數(shù)后，剩下的是發(fā)散收斂波。這個(gè)峰谷值序列已經(jīng)無(wú)法按照以上計(jì)數(shù)原則進(jìn)行循環(huán)提取。根據(jù) ASTM[9]標(biāo)準(zhǔn)中的介紹，可按照半循環(huán)計(jì)數(shù)法，直接從發(fā)散波中提取半循環(huán)。

基于以上流程，本文采用 Matlab 軟件平臺(tái)編制了雨流計(jì)數(shù)法程序，用于提取錨泊線張力時(shí)歷曲線中張力范圍大小及循環(huán)次數(shù)，以計(jì)算纜線的疲勞損傷。該Matlab 程序基于無(wú)量綱隨機(jī)載荷數(shù)據(jù)進(jìn)行程序運(yùn)算，并與手工處理結(jié)果對(duì)比驗(yàn)證，結(jié)果表明程序處理結(jié)果與手工處理結(jié)果一致，驗(yàn)證了雨流計(jì)數(shù)法程序的正確性。

1.3 錨泊線的 T-N 曲線

Miner 理論中使用的 S-N 曲線[10]是表示一定循環(huán)特征下標(biāo)準(zhǔn)試件的疲勞強(qiáng)度與疲勞壽命之間關(guān)系的曲線，通常它是由大量實(shí)驗(yàn)結(jié)果來(lái)確定的。在錨泊線的疲勞分析中一般使用的是 T-N 曲線，即只考慮張力作用產(chǎn)生的疲勞，而不考慮彎曲、扭轉(zhuǎn)等別的作用因素，其計(jì)算公式為：

式中：N 為材料達(dá)到破壞時(shí)的循環(huán)次數(shù)；R 為某海況下纜線張力范圍與最小斷裂強(qiáng)度的比值；M 和 K 為材料 T-N 曲線的參數(shù)，分別代表斜率與截距。鋼索材料的 T-N 曲線和相應(yīng)的 M，K 值一般采用美國(guó)石油協(xié)會(huì)標(biāo)準(zhǔn) API-2SK 規(guī)范推薦[11]。

1.4 Miner 線性累積損傷理論

將錨泊線所處的長(zhǎng)期海況離散成一系列短期海況i = 1，2，…，n。每種海況所產(chǎn)生的疲勞損傷為 Di，則產(chǎn)生的累積疲勞損傷為：

每種短期海況由系泊系統(tǒng)所處海洋環(huán)境的風(fēng)、浪、流參數(shù)及方向確定，此外還要描述每種海況的發(fā)生概率 pi。系泊系統(tǒng)的疲勞壽命為：

單個(gè)海況下的疲勞損傷可由下式計(jì)算：

式中：ni為短期海況 i 下的應(yīng)力循環(huán)次數(shù)；Ri為 i 海況下標(biāo)準(zhǔn)化處理后的應(yīng)力范圍；為海況的期望值。

2 錨泊線疲勞損傷評(píng)估算例

2.1 算例分析

本文計(jì)算模型為 1 艘外轉(zhuǎn)塔式 FPSO，工作海域水深 520 m?？傞L(zhǎng) 200 m，船寬 40 m，型深 16.15 m，吃水 8 m。

錨泊纜的布置形式影響系統(tǒng)的剛度。本文的系泊系統(tǒng)由 9 根系泊線組成，每 3 根系泊線 1 組，線組之間的夾角為 120°，每組 3 根系泊線之間的夾角為 5°。錨泊纜的布置情況可參考圖 4。系泊系統(tǒng)采用懸鏈線式系泊方式，懸鏈線選用鋼鏈-鋼纜-鋼鏈的組合形式（與錨、轉(zhuǎn)塔連接部分為鋼鏈，中間部分為鋼纜）。這種組合形式的優(yōu)點(diǎn)是：兩端的鋼鏈耐磨，能夠延長(zhǎng)系泊線的工作壽命；中間的鋼纜質(zhì)量小，可減輕系泊系統(tǒng)總體的重量，以適用更深的工作水域。

表 1 為系泊纜線的材料參數(shù)。

在本文的計(jì)算中，鋼鏈為普通的無(wú)檔鏈節(jié)，即取M = 3.0，K = 316；鋼索為螺旋狀絞繩，取 M = 5.05，，其中 Lm為鋼索的平均載荷與最小斷裂強(qiáng)度的比值，故在不同海況下 K 值是變化的。本文研究疲勞分析對(duì)象為錨泊線中間的鋼鏈及鋼纜，并認(rèn)為鋼鏈與鋼纜的的張力響應(yīng)相同。

圖 4 FPSO 系泊系統(tǒng)錨泊線布置Fig. 4 The arrangement of mooring lines

表 1 錨泊線組成及材料參數(shù)Tab. 1 The material parameter of mooring lines

2.2 計(jì)算結(jié)果

2.2.1 錨泊線張力時(shí)間歷程及疲勞載荷譜

基于多體水動(dòng)力學(xué)軟件 AQWA 中的纜索動(dòng)力學(xué)功能計(jì)算 FPSO 系泊系統(tǒng)的時(shí)域響應(yīng)[7]，可以得到 FPSO南海海域極限海況 180° 浪向角的錨泊纜張力時(shí)間歷程，從中提取每根系泊纜的張力最大值、最小值、平均值和標(biāo)準(zhǔn)差（見(jiàn)表 2）。在模擬中計(jì)算時(shí)長(zhǎng) 10 000 s，步長(zhǎng) 0.1 s。一年內(nèi)有 3 156 個(gè) 10 000 s，所以錨泊纜一年的疲勞損傷為 10 000 s 模擬時(shí)間內(nèi)產(chǎn)生的疲勞損傷的 3 156 倍。

從表 2 可看出，2 號(hào)錨泊纜 L2的張力最大。圖 5給出了其張力響應(yīng)時(shí)間歷程曲線。本文以受力最大的2 號(hào)錨泊纜為例，采用編制的雨流計(jì)數(shù)法程序?qū)ζ溥M(jìn)行統(tǒng)計(jì)，獲取該錨泊線的張力范圍概率分布直方圖—疲勞載荷譜。

疲勞分析中波浪的長(zhǎng)期狀態(tài)可看作由許多短期海況序列組成。每個(gè)短期海況可以用波高 HS、平均跨零周期 TZ、波譜、流速、風(fēng)速等參數(shù)來(lái)描述，且每個(gè)海況都有其確定的發(fā)生概率 P。表 3 為南海 S4 海域一年的波浪散布圖。

依據(jù)上表波浪散布圖資料，可將長(zhǎng)期海況離散成64 個(gè)短期海況。根據(jù)波高 HS、平均跨零周期 TZ、流速、風(fēng)速，發(fā)生概率等參數(shù)特點(diǎn)，選取其中具有代表性的 6 個(gè)海況條件：10，21，28，38，50，58，分別得到其疲勞載荷譜如圖 6 所示。

表 2 南海海域極限海況 180 域浪向角系泊系統(tǒng)纜線動(dòng)力響應(yīng)統(tǒng)計(jì)值Tab. 2 The dynamic response of mooring lines in extreme ocean condition

圖 5 纜線 L2 張力響應(yīng)時(shí)歷曲線Fig. 5 The tension response of Mooring line L2

表 3 南海 S4 海域波浪雙變量散布圖Tab. 3 The Scatter diagrams of Nanhai S4sea waters

圖 6 典型海況下的疲勞載荷譜Fig. 6 The fatigue load spectrum of typical ocean condition

從圖 6 可以看出：1）短期海況下的纜線張力疲勞載荷譜并不服從常見(jiàn)的理論概率分布，而僅僅反映各短期海況下的錨泊線張力范圍；2）較小張力范圍出現(xiàn)的概率最大，張力范圍越大，其出現(xiàn)概率越小。

2.2.2 長(zhǎng)期海況下的累積疲勞損傷計(jì)算

根據(jù) 2.4 Miner 線性累積損傷理論，F(xiàn)PSO 系泊系統(tǒng)在長(zhǎng)期海況作用下產(chǎn)生的疲勞損傷 D 可以通過(guò) n 個(gè)短期海況下的疲勞損傷 Di累積得到，表 4 及表 5 給出了鋼鏈和鋼纜在不同海況下的疲勞損傷計(jì)算結(jié)果。

表 4 鋼鏈疲勞損傷計(jì)算表Tab. 4 The fatigue damage of steel chain

表 5 鋼纜疲勞損傷計(jì)算表Tab. 5 The fatigue damage of wire rope

3 結(jié) 語(yǔ)

本文采用雨流計(jì)數(shù)法對(duì) FPSO 系泊系統(tǒng)進(jìn)行疲勞損傷評(píng)估，得到以下結(jié)論：

1）通過(guò)計(jì)算 64 個(gè)短期海況下 FPSO 系泊系統(tǒng)的疲勞損傷并累積起來(lái)，可以發(fā)現(xiàn)錨泊線的疲勞損傷主要發(fā)生在鋼鏈部分，鋼鏈的疲勞損傷遠(yuǎn)大于鋼纜的疲勞損傷。在 1 年的長(zhǎng)期海況條件下，鋼鏈的疲勞損傷約為鋼纜疲勞損傷的 1 643 倍。這表明鋼纜的疲勞性能要顯著優(yōu)于鋼鏈的疲勞性能。在實(shí)際工程應(yīng)用中，鋼鏈的疲勞情況需要引起特別注意。在它的日常維護(hù)中，應(yīng)重點(diǎn)注意鋼鏈的維護(hù)保養(yǎng)。

2）在計(jì)算每個(gè)短期海況條件下的疲勞損傷時(shí)可以發(fā)現(xiàn)，較小的張力范圍引起的疲勞損傷幾乎可以忽略不計(jì)。在累積計(jì)算長(zhǎng)期海況條件下的疲勞損傷時(shí)可以發(fā)現(xiàn)，若干波高 Hs較小的短期海況所引起疲勞損傷也可以忽略不計(jì)。

3）在不考慮各短期海況分布概率的影響，鋼鏈與鋼纜在各個(gè)短期海況條件下的疲勞損傷如圖 7 所示。可以發(fā)現(xiàn)，在同一短期海況條件下，鋼鏈的疲勞損傷也要遠(yuǎn)大于鋼纜的疲勞損傷。在跨零周期不變 Tz的情況下，隨著波高 Hs的增大，F(xiàn)PSO 的運(yùn)動(dòng)變得劇烈會(huì)導(dǎo)致鋼鏈或者鋼纜的疲勞損傷的增大。

圖 7 疲勞損傷隨海況變化圖Fig. 7 The fatigue damage in different ocean condition

[1]LUO Y. Tension response statistics and fatigue analysis of catenary mooring lines[C]//Proceeding of first 90 European offshore mechanics symposium. Trondheim, Norway: International Society of Offshore & Polar Engineers, 1990: 213–220.

[2]LARSEN K, MATHISEN J. Reliability-based fatigue analysis of mooring lines[C]//Proceedings of the 15th international conference on offshore mechanics and arctic engineering. Florence, Italy: American Society of Mechanical Engineers, 1996.

[3]MATHISEN J, H?RTE T, MOE V, et al. DEEPMOOR-design methods for deep water mooring systems. Calibration of a fatigue limit state[R]. DNV Report, No. 98–3110, 1999.

[4]MATHISEN J, LARSEN K. Risk-based inspection planning for mooring chain[J]. Journal of Offshore Mechanics and Arctic Engineering, 2004, 126(3): 250–257.

[5]SIGURDSSON G, LOTSBERG I, LANDET E. Risk based inspection of FPSOs[C]//International conference on offshore mechanics and arctic engineering, OMAE 2000. New Orleans: OMAE, 2000.

[6]VAZQUEZ-HERNANDEZ A O, ELLWANGER G B, SAQRILO L V S. Fatigue analysis and reliability of floating production systems mooring lines in deepwater[C]//Proceedings of the 17th international conference on offshore mechanics and arctic engineering. Lisbon, Portugal, 1998.

[7]劉元丹, 劉敬喜, 譚安全. 單點(diǎn)系泊FPSO風(fēng)浪流載荷下運(yùn)動(dòng)及其系泊力研究[J]. 船海工程, 2011, 40(6): 146–149.

[8]董樂(lè)義, 羅俊, 程禮. 雨流計(jì)數(shù)法及其在程序中的具體實(shí)現(xiàn)[J]. 航空計(jì)測(cè)技術(shù), 2004, 24(3): 38–40.

[9]ASTM. Standard practices for cycle counting in fatigue analysis: ASTME1049-85(2005)[S]. ASTM International, 2005.

[10]DNV. Offshore standard-position mooring[R]. DNV OS-E301, 2004.

[11]API. Recommended practice for design and analysis of stationkeeping systems for floating structures: API RP 2SK[S]. American Petroleum Institute, 2005.

Evaluation on fatigue strength of mooring line based on rain-flow count method

LIU Juan, SUN Zheng, LIU Jing-xi
(College of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan 430074, China)

The mooring system of FPSO composed of multiple mooring lines. It is critical to evaluate the fatigue strength of the mooring line to guarantee the safety of the structure. In this paper improved rain-flow count method, T- N curve and Miner linear fatigue cumulative damage theory were used to evaluate the fatigue performance of mooring line. From the analyses of the fatigue damage of mooring line, we can get the information of the disadvantageous position, which provide important basis for the maintenance of FPSO mooring system.

rain-flow count method；mooring line；fatigue strength；stress-time curve

TP301.6

A

1672–7619(2016)12–0059–06

10.3404/j.issn.1672–7619.2016.12.012

2016–03–29；

2016–04–20

國(guó)家自然科學(xué)基金資助項(xiàng)目(51579110)；國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)資助項(xiàng)目(2012AA112601)

劉娟(1989–)， 女，碩士研究生，主要研究方向?yàn)榇芭c海洋結(jié)構(gòu)物設(shè)計(jì)與制造。