超深水鉆井船總強(qiáng)度評(píng)估的應(yīng)力設(shè)計(jì)波法研究

劉仁昌，陳國建，王金峰，趙志堅(jiān)，曹 凱，曹洪建，張利軍，徐海霞，萬家平

(1. 中遠(yuǎn)海運(yùn)重工有限公司，遼寧 大連 116600；2. 必維船級(jí)社（中國）有限公司，上海 200011)

0 引 言

鉆井船是進(jìn)行深遠(yuǎn)海油氣資源開發(fā)的重要裝備，其結(jié)構(gòu)強(qiáng)度是關(guān)鍵安全因子。鉆井船結(jié)構(gòu)不但需要規(guī)范校核，同時(shí)需要設(shè)計(jì)波法直接強(qiáng)度計(jì)算[1-4]。隨著鉆井船創(chuàng)新型設(shè)計(jì)的不斷涌現(xiàn)，船體功能設(shè)計(jì)的拓展，船體變地更加復(fù)雜、獨(dú)特。因此，工程中迫切需求新的設(shè)計(jì)波方法，可針對(duì)獨(dú)特的船體結(jié)構(gòu)更加全面、精準(zhǔn)的進(jìn)行強(qiáng)度極限的校核，以保證結(jié)構(gòu)安全。本文針對(duì)超深水鉆井船總強(qiáng)度評(píng)估的等效設(shè)計(jì)波方法進(jìn)行探索研究，對(duì)比了傳統(tǒng)主控載荷設(shè)計(jì)波和應(yīng)力設(shè)計(jì)波法。在傳統(tǒng)設(shè)計(jì)波強(qiáng)度評(píng)估的同時(shí)，針對(duì)月池角隅探索了以應(yīng)力為目標(biāo)的設(shè)計(jì)波法。通過比較傳統(tǒng)設(shè)計(jì)波和應(yīng)力設(shè)計(jì)波2 套結(jié)果，探索了應(yīng)力設(shè)計(jì)波法在超深水新型鉆井船上的應(yīng)用，進(jìn)一步保證結(jié)構(gòu)安全。

1 基本理論方法與設(shè)計(jì)波方法

1.1 水動(dòng)力理論

本文采用基于三維勢(shì)流理論[5]的水動(dòng)力求解器對(duì)船體波浪載荷進(jìn)行計(jì)算。船舶在微幅（波高波長(zhǎng)比遠(yuǎn)小于1）規(guī)則波中做微幅運(yùn)動(dòng)，流場(chǎng)速度勢(shì)Φ（x，y，z，t）滿足：

速度勢(shì)滿足線性自由面條件：

速度勢(shì)可劃分為入射波和擾動(dòng)勢(shì)兩部分，記

其中：φ0為入射波的一階速度勢(shì)；φp為擾動(dòng)勢(shì)。

在無限水深情況下，入射波速度勢(shì)為：

擾動(dòng)勢(shì)φP在物面SH上適合物面條件：

根據(jù)勢(shì)流理論，將物面條件離散，通過數(shù)值方法求解整個(gè)流場(chǎng)的速度勢(shì)，進(jìn)一步可求解流場(chǎng)所有信息，如附連水質(zhì)量、水動(dòng)壓力及船體運(yùn)動(dòng)[6]等，進(jìn)而計(jì)算船體的水動(dòng)力載荷。

1.2 水動(dòng)力-結(jié)構(gòu)耦合方法

水動(dòng)力載荷包括慣性力和水動(dòng)壓力，需精確傳遞到結(jié)構(gòu)上，并達(dá)到自平衡。慣性力通常可直接加載。為實(shí)現(xiàn)精確的傳遞水動(dòng)力載荷，必維船級(jí)社研發(fā)的水動(dòng)力-結(jié)構(gòu)耦合軟件Homer[7]采用最新研究成果，在有限元網(wǎng)格高斯積分點(diǎn)上直接計(jì)算水動(dòng)壓力[8]，即在船體需要的位置點(diǎn)直接計(jì)算流體勢(shì)函數(shù)：

其中：xs=（xs，ys，zs）表示結(jié)構(gòu)上的點(diǎn)，xh=（xh，yh，zh）表示水動(dòng)力網(wǎng)格上的點(diǎn)。

之后進(jìn)一步求解船體的水動(dòng)壓力載荷：

其中：上標(biāo)“s”為結(jié)構(gòu)計(jì)算對(duì)應(yīng)的物理量

此方法可將水動(dòng)力載荷“完美”傳遞到結(jié)構(gòu)模型。由于直接在結(jié)構(gòu)模型上求解運(yùn)動(dòng)方程，可達(dá)到自平衡。

1.3 應(yīng)力設(shè)計(jì)波方法

主控載荷設(shè)計(jì)波法中，選取若干船體典型主控載荷為目標(biāo)，根據(jù)主控載荷的RAO 峰值、峰值頻率、相位等參數(shù)構(gòu)造等效設(shè)計(jì)波。傳統(tǒng)主控載荷設(shè)計(jì)波法流程簡(jiǎn)單、明確、計(jì)算量小，適合成熟船型。然而，對(duì)結(jié)構(gòu)新穎的船型，某些結(jié)構(gòu)的受力是多重載荷綜合作用結(jié)果，如月池角隅等結(jié)構(gòu)。單一最大主控載荷可能難以覆蓋最危險(xiǎn)情況，這就需要應(yīng)力設(shè)計(jì)波方法進(jìn)行進(jìn)一步分析。

應(yīng)力設(shè)計(jì)波法以關(guān)鍵結(jié)構(gòu)處應(yīng)力為目標(biāo)，是目前船舶與海工領(lǐng)域較為前沿的設(shè)計(jì)波方法。首先，確定船體結(jié)構(gòu)的關(guān)鍵結(jié)構(gòu)節(jié)點(diǎn)，構(gòu)建應(yīng)力RAO 曲線。再通過譜分析獲得應(yīng)力極值，作為設(shè)計(jì)波目標(biāo)值。圖1 為對(duì)比主控載荷設(shè)計(jì)波與應(yīng)力設(shè)計(jì)波流程。與主控載荷設(shè)計(jì)波加載不同，應(yīng)力設(shè)計(jì)波法需在時(shí)域中進(jìn)行弱非線性修正以精確表達(dá)波面起伏（見圖2），求解整個(gè)波浪周期內(nèi)結(jié)構(gòu)響應(yīng)，進(jìn)行強(qiáng)度校核。

圖1 主控載荷設(shè)計(jì)波與應(yīng)力設(shè)計(jì)波方法流程對(duì)比Fig. 1 Comparison of the design wave approach based on main loads and based on stress

應(yīng)力設(shè)計(jì)波法計(jì)算量大，代價(jià)高。由于計(jì)算中針對(duì)特定結(jié)構(gòu)進(jìn)行了全浪向和全頻率范圍的計(jì)算和篩選，并考慮了波面修正，計(jì)算得到的應(yīng)力結(jié)果更符合實(shí)際情況，對(duì)結(jié)構(gòu)的評(píng)判期望更加合理。

圖2 波面若非線性修正前后船表壓力對(duì)比Fig. 2 Comparison of Wave surface pressure before and after correction

2 超深水鉆井船應(yīng)力主控設(shè)計(jì)波研究

本文研究了生存工況下某超深水鉆井船的強(qiáng)度。超深水鉆井船的主尺度參數(shù)如表1 所示。

表1 超深水鉆井船主尺度Tab. 1 Main parameters of ultra-deep-water drilling ship

本船特點(diǎn)在于月池比一般鉆井船大，月池上方井架結(jié)構(gòu)比一般鉆井船重很多，且結(jié)構(gòu)形式更復(fù)雜。因此，針對(duì)月池結(jié)構(gòu)尤其是月池角隅位置，在傳統(tǒng)設(shè)計(jì)波法之外，進(jìn)一步采用應(yīng)力設(shè)計(jì)波方法，全面深入研究角隅處極限強(qiáng)度，以保證船體強(qiáng)度安全。

2.1 船體載荷與結(jié)構(gòu)響應(yīng)計(jì)算

按船體型線和結(jié)構(gòu)，采用有限元軟件FEMAP 建立水動(dòng)力計(jì)算模型和有限元模型，如圖3 所示。

圖3 水動(dòng)力計(jì)算模型和有限元模型Fig. 3 Hydrodynamic calculation model and finite element model

采用三維勢(shì)流理論水動(dòng)力軟件HydroSTAR[9]開展水動(dòng)力性能計(jì)算。計(jì)算得到船體在生存工況下的運(yùn)動(dòng)響應(yīng)RAO 如圖4 所示。90°（橫浪）時(shí)，船體垂蕩、橫搖運(yùn)動(dòng)RAO 峰值最大，但垂蕩和橫搖峰值頻率不同。0°（隨浪）或180°（迎浪）時(shí)，船體縱搖運(yùn)動(dòng)RAO峰值最大。

圖4 船體運(yùn)動(dòng)RAOFig. 4 RAO of ship motion

船體主控載荷RAO，如圖5 所示。包括船中垂向波浪彎矩RAO、1/4 船長(zhǎng)和3/4 船長(zhǎng)處垂直波浪剪力、船中扭矩RAO 和井架重心加速度RAO。對(duì)比運(yùn)動(dòng)RAO 和主控載荷RAO 可看出，兩RAO 峰值緊密相關(guān)。垂直剪力Fz、垂向彎矩VBM 的峰值發(fā)生在與垂蕩、縱搖峰值頻率相近區(qū)域；船舯剖面的扭矩RAO 峰值頻率與橫搖RAO 峰值頻率、井架重心加速度峰值頻率接近。船中井架質(zhì)量巨大，其橫搖對(duì)船中的扭矩產(chǎn)生較大影響。

利用水動(dòng)力-結(jié)構(gòu)耦合軟件HOMER，結(jié)合FEM 模型，將水動(dòng)力載荷加載到有限元模型，求解各個(gè)浪向、各個(gè)頻率下的結(jié)構(gòu)單元的應(yīng)力響應(yīng)RAO。本文選擇月池甲板角隅和外底板角隅處的系列單元作為主控單元，圖6 為兩處單元中面Von Mises 應(yīng)力RAO 曲線。

在月池角隅處單元，其應(yīng)力RAO 同樣呈明顯“雙峰”形態(tài)，即應(yīng)力RAO 峰值和峰值頻率隨浪向不同而改變，第1 個(gè)峰值在波浪頻率0.5 rad/s 附近，此時(shí)入射波浪向?yàn)?80°（迎浪）；第2 個(gè)峰值在波浪頻率0.65 rad/s 附近，入射波浪向?yàn)?35°（頂斜浪）。結(jié)合船體運(yùn)動(dòng)和主控載荷：第1 個(gè)應(yīng)力峰值出現(xiàn)的頻率、浪向與船舯垂向波浪彎矩的頻率、浪向相同。如果以船舯垂向波浪彎矩為主控載荷，獲得的等效設(shè)計(jì)波將使得月池角隅應(yīng)力產(chǎn)生峰值。第2 個(gè)應(yīng)力峰值產(chǎn)生的頻率、浪向與所有主控載荷設(shè)計(jì)波中頻率、浪向都不相同：浪向135°時(shí)，多種載荷綜合作用下的月池角隅應(yīng)力水平更高，月池角隅應(yīng)力對(duì)波浪力的作用更為敏感。

圖5 船體主控載荷RAOFig. 5 RAO of ship main loads

圖6 月池角隅單元應(yīng)力RAOFig. 6 Stress RAO of element at moonpool corner

2.2 主控載荷極值分析

結(jié)合實(shí)際工況對(duì)應(yīng)波浪譜和海況信息進(jìn)行譜分析，預(yù)報(bào)波浪載荷[10]在實(shí)際海況條件下的最大載荷，及最大載荷對(duì)應(yīng)的浪向。生存工況對(duì)應(yīng)的短期海況如表2 所示，對(duì)應(yīng)的波浪譜為JONSWAP 譜。

表2 生存工況短期海況參數(shù)Tab. 2 Short term sea state parameters for survival condition

對(duì)船體波浪載荷RAO 結(jié)果進(jìn)行譜分析，得到百年一遇主控載荷極值及對(duì)應(yīng)的浪向，如表3 所示。表3同時(shí)給出了主控載荷極值對(duì)應(yīng)的RAO 峰值和峰值頻率。

2.3 等效設(shè)計(jì)波參數(shù)確定

根據(jù)譜分析獲得主控載荷極值、浪向角，結(jié)合相應(yīng)的RAO 曲線，可以獲得波浪的頻率、相位。線性波浪幅則根據(jù)載荷極值和RAO 峰值的比率進(jìn)行計(jì)算：

等效規(guī)則設(shè)計(jì)波的最終表達(dá)形式為：

以垂直剪力、垂向彎矩、扭矩、月池角隅應(yīng)力等極值為目標(biāo)的設(shè)計(jì)波如表4 所示。

表3 主控載荷的譜分析結(jié)果Tab. 3 Spectral analysis results of main loads

表4 主控載荷、主控應(yīng)力設(shè)計(jì)波參數(shù)Tab. 4 Design wave parameters for main load and main stress

2.4 船體結(jié)構(gòu)極值響應(yīng)分析

分別將2 套設(shè)計(jì)波加載到船體，進(jìn)行分析，得到結(jié)構(gòu)響應(yīng)。2 套設(shè)計(jì)波在相同結(jié)構(gòu)單元位置如圖7 所示。

表5 給出了應(yīng)力設(shè)計(jì)波法計(jì)算月池角隅應(yīng)力極值及應(yīng)力設(shè)計(jì)波參數(shù)。表6 和表7 給出了2 種設(shè)計(jì)波方法的結(jié)果對(duì)比，可看出利用應(yīng)力設(shè)計(jì)波法得到的單元極限強(qiáng)度所對(duì)應(yīng)的設(shè)計(jì)波浪向、頻率與單元應(yīng)力幅值，都與主控載荷設(shè)計(jì)波法有所不同。從浪向看，由于主控載荷設(shè)計(jì)波中，沒有斜浪設(shè)計(jì)波，故利用主控載荷設(shè)計(jì)波得到的是單元應(yīng)力RAO 的第1 個(gè)峰值所產(chǎn)生的應(yīng)力極值。而利用應(yīng)力設(shè)計(jì)波方法，單元響應(yīng)的極值浪向發(fā)生在頂斜浪，不僅受船體梁載荷影響，更是各種載荷綜合作用的結(jié)果。

應(yīng)力設(shè)計(jì)波得到各個(gè)單元的應(yīng)力極值基本上大于或等于主控載荷設(shè)計(jì)波結(jié)果，如圖8 所示。僅在月池外底板橫向單元上，應(yīng)力水平稍低。從應(yīng)力極值的對(duì)比可看出，利用應(yīng)力設(shè)計(jì)波法得到的應(yīng)力水平通常更高，基本涵蓋主控載荷設(shè)計(jì)波結(jié)果。應(yīng)力隨浪向、頻率的分布，以及得到的對(duì)應(yīng)最大屈服應(yīng)力，對(duì)結(jié)構(gòu)設(shè)計(jì)的加強(qiáng)具有更加明確的針對(duì)性和指導(dǎo)意義。

從計(jì)算量和計(jì)算時(shí)間上看，應(yīng)力設(shè)計(jì)波的計(jì)算量與耗時(shí)遠(yuǎn)大于主控載荷設(shè)計(jì)波。因此，利用應(yīng)力設(shè)計(jì)波法對(duì)創(chuàng)新型、高附加值船型進(jìn)行強(qiáng)度分析和校核，更有效果。

圖7 月池角隅單元位置信息Fig. 7 Location information of elements at moonpool corner

表5 月池角隅單元應(yīng)力設(shè)計(jì)波參數(shù)Tab. 5 Design wave parameters for stress of moonpool corner elements

表6 月池主甲板角隅單元應(yīng)力計(jì)算結(jié)果比較-Von mises stressTab. 6 VM-Stress at the of main deck corner of moonpool

表7 月池外底板角隅單元應(yīng)力計(jì)算結(jié)果比較-Von mises stressTab. 7 VM-Stress of bottom outer shell corner of moonpool

圖8 月池角隅單元總應(yīng)力對(duì)比Fig. 8 Total stress of elements at moonpool corner

3 結(jié) 論

本文針對(duì)超深水鉆井船總強(qiáng)度評(píng)估的等效設(shè)計(jì)波法進(jìn)行探索研究，比較分析主控載荷設(shè)計(jì)波法和應(yīng)力設(shè)計(jì)波新方法的結(jié)果，探索了應(yīng)力設(shè)計(jì)波法在超深水鉆井船上的應(yīng)用。

1）針對(duì)關(guān)鍵結(jié)構(gòu)，應(yīng)力設(shè)計(jì)波法更全面的對(duì)浪向、頻率展開搜索，獲得的應(yīng)力極值對(duì)應(yīng)135°頂斜浪，與傳統(tǒng)方法的應(yīng)力極值浪向180°（迎浪）存在差異。傳統(tǒng)方法以船體梁截面扭矩為主控載荷時(shí)，基于本船型得到的設(shè)計(jì)波為90°橫浪，此時(shí)總扭矩主要由鉆井模塊（重量大、重心高）慣性力引起。然而，此扭矩導(dǎo)致的應(yīng)力更直觀體現(xiàn)在井架支撐結(jié)構(gòu)，由45°斜浪船體梁扭轉(zhuǎn)引發(fā)的月池角隅應(yīng)力卻被忽略，這是傳統(tǒng)設(shè)計(jì)波法在超大型鉆井船（鉆井模塊重量大、重心高）總強(qiáng)度評(píng)估中的缺陷，應(yīng)力設(shè)計(jì)波法直接以目標(biāo)區(qū)域應(yīng)力為關(guān)注點(diǎn)，彌補(bǔ)了此漏洞。

2）從屈服結(jié)果的對(duì)比分析來看，應(yīng)力設(shè)計(jì)波法得到的單元屈服應(yīng)力普遍大于傳統(tǒng)設(shè)計(jì)波法，說明目前工業(yè)界從船體梁整體失效目標(biāo)考慮的設(shè)計(jì)波法在屈服評(píng)估時(shí)的極值涵蓋還不夠全面，關(guān)鍵結(jié)構(gòu)最大屈服應(yīng)力的搜索存在可能遺漏和不足，盡管傳統(tǒng)方法較大程度上可滿足船舶與海工結(jié)構(gòu)設(shè)計(jì)需求。

3）對(duì)比結(jié)果顯示，應(yīng)力設(shè)計(jì)波法得到的應(yīng)力中，存在極少數(shù)不符合規(guī)律的點(diǎn)（屈服應(yīng)力小于傳統(tǒng)方法），有必要針對(duì)該方法進(jìn)一步開展全流程影響因素研究。

4）應(yīng)力設(shè)計(jì)波法目前僅限屈服評(píng)估，對(duì)屈曲和疲勞還無法應(yīng)用。一方面，屈曲和疲勞基于的應(yīng)力分量更多，另一方面，屈曲和疲勞還受應(yīng)力和材料屬性之外的多個(gè)因素影響，更復(fù)雜，以屈曲、疲勞為目標(biāo)的應(yīng)力設(shè)計(jì)波法還有待進(jìn)一步研究，以推廣此方法在更廣范圍的應(yīng)用，發(fā)揮其在結(jié)構(gòu)評(píng)估中的精確性優(yōu)勢(shì)。