劉 旭, 陳曉明, 蔡連財(cái)

(1.中遠(yuǎn)海運(yùn)特種運(yùn)輸股份有限公司， 廣州 510623； 2.中國(guó)船舶重工集團(tuán)國(guó)際工程有限公司， 北京 100024 )

半潛船主要通過(guò)一系列下潛、上浮操作或碼頭滾裝的方式來(lái)裝卸和運(yùn)輸不可分割的超大型設(shè)備[1-2]，是重要的新型海洋工程特種運(yùn)輸船型，目前最大載重噸已發(fā)展到11萬(wàn)t，承運(yùn)的單件貨物最大重量也達(dá)到6.4萬(wàn)t。半潛船船寬與船長(zhǎng)、型深、吃水的比值明顯大于常規(guī)運(yùn)輸船。半潛船運(yùn)輸?shù)奶卮笮驮O(shè)備和海洋平臺(tái)類別繁雜，每次航行過(guò)程中貨物的綁扎狀態(tài)都相差較大，船舶運(yùn)動(dòng)響應(yīng)變化也較大，稍有不慎，很容易引起貨物移動(dòng)、傾倒，甚至引發(fā)全船舶傾覆等安全事故。[3-4]在航運(yùn)過(guò)程中，為確保航行安全，半潛船每一次運(yùn)載特大件貨物時(shí)均需要針對(duì)具體的運(yùn)載狀態(tài)進(jìn)行耐波性分析，并依據(jù)計(jì)算結(jié)果進(jìn)行牢固綁扎。因此，全面掌握半潛船在不同運(yùn)載航行狀態(tài)下的水動(dòng)力特性和耐波性規(guī)律是非常有必要的。

作為一種高技術(shù)高附加值的特種船型，半潛船設(shè)計(jì)難度很大，2006年之前僅有荷蘭具備完整的建造設(shè)計(jì)能力，近些年，我國(guó)通過(guò)引進(jìn)吸收的方式也逐漸具備一定的設(shè)計(jì)能力，但仍與世界先進(jìn)水平存在巨大差距。迄今為止半潛船還沒(méi)有統(tǒng)一的技術(shù)要求，國(guó)內(nèi)外半潛船的關(guān)鍵技術(shù)資料、研究報(bào)告和文獻(xiàn)資料極少，針對(duì)半潛船耐波性的研究數(shù)量更少，白鑫等[5]、楊索賢等[6]分別利用HydroStar軟件預(yù)報(bào)半潛船的波浪載荷，得出規(guī)范公式無(wú)法對(duì)波浪載荷做出合理的評(píng)估而采用直接算法較為合理的結(jié)論。莫瑞芳等[7]利用Sesam/PatranPre軟件進(jìn)行零航速耐波性分析，研究浪向?qū)Υ故?、橫搖的響應(yīng)規(guī)律。目前，公開的文獻(xiàn)缺少半潛船耐波性試驗(yàn)分析及數(shù)值模擬與試驗(yàn)的對(duì)比分析等方面的研究，不利于半潛船技術(shù)的發(fā)展。

基于此，本文以5萬(wàn)噸級(jí)半潛船為研究對(duì)象，開展模型試驗(yàn)和實(shí)船試驗(yàn)、數(shù)值模擬方法在半潛船耐波性分析領(lǐng)域的應(yīng)用研究，且通過(guò)對(duì)比分析探討半潛船航行狀態(tài)下的水動(dòng)力規(guī)律，為同類型半潛船耐波性分析提供有益參考。

1 耐波性試驗(yàn)

半潛船耐波性試驗(yàn)方法包括模型試驗(yàn)和實(shí)船試驗(yàn)。耐波性模型試驗(yàn)主要包括靜水試驗(yàn)、規(guī)則波試驗(yàn)、不規(guī)則波試驗(yàn)和瞬態(tài)波試驗(yàn)。[8]其中,不規(guī)則波模型試驗(yàn)主要通過(guò)船舶模型在不規(guī)則波中的運(yùn)動(dòng)響應(yīng)試驗(yàn)，以確定船體在不同狀態(tài)下各種頻率響應(yīng)函數(shù)及加速度等。耐波性實(shí)船試驗(yàn)項(xiàng)目包括船舶運(yùn)動(dòng)響應(yīng)及垂向加速度等測(cè)試。[9]不規(guī)則波模型試驗(yàn)和實(shí)船試驗(yàn)?zāi)茌^真實(shí)地反映天然海浪作用下的船舶運(yùn)動(dòng)規(guī)律，前者是船舶設(shè)計(jì)期間確定船舶耐波性的重要依據(jù)，后者是評(píng)價(jià)和檢驗(yàn)船舶耐波性優(yōu)劣的最終依據(jù)。本文分別考慮了半潛船的這兩種試驗(yàn)情況。

1.1 試驗(yàn)參數(shù)和試驗(yàn)條件

為得到半潛船準(zhǔn)確的水動(dòng)力特性，給后續(xù)綁扎設(shè)計(jì)提供依據(jù)，需要開展半潛船的耐波性試驗(yàn)研究。目標(biāo)船型為我國(guó)開發(fā)的某5萬(wàn)噸級(jí)半潛船，目標(biāo)船型的實(shí)際尺寸參數(shù)和試驗(yàn)條件見表1。

耐波性模型試驗(yàn)在荷蘭MARIN水池進(jìn)行，采用先進(jìn)的搖板式多單元造波系統(tǒng)及消波系統(tǒng)。根據(jù)船舶與海洋工程模型試驗(yàn)原理和相關(guān)規(guī)范[8,10]，模型試驗(yàn)?zāi)Ｐ涂s尺比取1∶49.5，見圖1。

圖1 5萬(wàn)噸級(jí)半潛船模型

表1 半潛船模型參數(shù)

描述輕載工況滿載工況實(shí)船工況(空載工況)排水量/t54 461.0869 191.4450 275.50型長(zhǎng)/m216.70216.70216.70型寬/m434343艉吃水/m8.009.817.50艏吃水/m8.009.817.00初穩(wěn)心高值/m9.254.9218.28橫搖轉(zhuǎn)動(dòng)半徑/m17.2021.5012.97

1.2 不規(guī)則波試驗(yàn)方法分析

(1)

不規(guī)則波試驗(yàn)的RAO主要通過(guò)不規(guī)則波試驗(yàn)獲得半潛船在白噪聲波譜中的運(yùn)動(dòng)響應(yīng)，應(yīng)用譜分析技術(shù)將相應(yīng)的時(shí)域信號(hào)進(jìn)行傅里葉變換，得到頻域響應(yīng)，再進(jìn)行平滑處理，最終計(jì)算得出船舶在不同浪向角下各自由度的RAO。

2 數(shù)值模擬

2.1 水動(dòng)力分析數(shù)值方法

水動(dòng)力分析數(shù)值方法是船舶耐波性分析的基本技術(shù)手段。與之前的研究?jī)H數(shù)值模擬不同[5-7]，本文將三維繞輻射勢(shì)流理論的OCTOPUS計(jì)算程序?qū)Π霛摯哪筒ㄐ灶A(yù)報(bào)結(jié)果與模型試驗(yàn)和實(shí)船試驗(yàn)結(jié)果進(jìn)行對(duì)比分析，以此驗(yàn)證數(shù)值模擬方法的合理性。在計(jì)算過(guò)程中，流場(chǎng)的求解基于頻域格林函數(shù)邊界元法進(jìn)行，一階問(wèn)題采用基于Korvin源勢(shì)法，二階采用Molin法。OCTOPUS水動(dòng)力數(shù)值方法具體的實(shí)現(xiàn)原理見圖2。

實(shí)際海況是不規(guī)則且隨機(jī)的，譜分析法可給出特定不規(guī)則波海況下船體的運(yùn)動(dòng)響應(yīng)。[11]本文計(jì)算采用JONSWAP 波浪譜，其定義為

(2)

圖2 OCTOPUS水動(dòng)力數(shù)值方法具體的實(shí)現(xiàn)原理

航行船舶的響應(yīng)譜的譜矩可表達(dá)為

(3)

式(3)中：Sζ(ω)為海浪譜密度函數(shù);ω為波浪圓頻率。

2.2 計(jì)算模型

水動(dòng)力計(jì)算模型包括濕表面模型和質(zhì)量模型。濕表面模型是基于半潛船的肋骨型線圖中的50個(gè)典型肋位，其中：艉部1/4船長(zhǎng)包含15個(gè)剖面；艏部1/4船長(zhǎng)包含15個(gè)剖面；舯部包含20個(gè)剖面。為保證計(jì)算精度，計(jì)算區(qū)域長(zhǎng)度和寬度方向均實(shí)取船寬的5倍。船體的網(wǎng)格沿船長(zhǎng)方向約55個(gè)，船寬方向約20個(gè)，在型線變化比較大的地方網(wǎng)格密度適當(dāng)加大，四周自由液面的網(wǎng)格與船體網(wǎng)格相匹配。船體水動(dòng)力計(jì)算模型見圖3。

圖3 半潛船水動(dòng)力計(jì)算模型

3 模型與實(shí)船試驗(yàn)工況及結(jié)果分析

3.1 不規(guī)則波模型試驗(yàn)與實(shí)船試驗(yàn)工況

為完整考察半潛船的航行耐波性特性，分別測(cè)試以下典型工況。其中，定義從艉部方向傳來(lái)波浪入射角為 0°，艏部方向傳來(lái)波浪入射角為180°。

1) 輕載工況1(吃水=8 m)，Hs=6.50 m,Tp=9.00 s，浪向90°/180°，航速10 kn。

2) 輕載工況2(吃水=8 m)，Hs=6.50 m,Tp=12.50 s，浪向?yàn)?0°/180°，航速為10 kn。

3) 滿載工況1(吃水=9.81 m)，Hs=6.50 m,Tp=9.00 s，浪向?yàn)?0°/180°，航速為10 kn。

4) 滿載工況2(吃水=9.81 m)，Hs=6.50 m,Tp=12.50 s，浪向?yàn)?0°/180°，航速為10 kn。

5) 實(shí)船工況(吃水=9.0 m)，Hs=1.5～4.0 m,Tp=10～12 s，浪向?yàn)?35°/-45°/-60°，航速為9.0～11.8 kn。

3.2 模型試驗(yàn)橫搖和垂蕩運(yùn)動(dòng)分析

半潛船模型試驗(yàn)的橫搖、垂蕩運(yùn)動(dòng)情況見圖4、圖5、圖6和圖7。

圖4 RMS橫搖運(yùn)動(dòng)幅值圖5 橫搖RAO

圖6 RMS垂蕩運(yùn)動(dòng)幅值圖7 垂蕩RAO

由圖4與圖5可知：

1) 在尾斜浪和橫浪之間區(qū)域的海況下，由于波浪的遭遇頻率接近于橫搖運(yùn)動(dòng)自然周期，船舶處于共振區(qū)，因此半潛船的橫搖激勵(lì)最大，橫搖響應(yīng)最為劇烈。這與之前的結(jié)論是一致的。[9]

2) 橫搖運(yùn)動(dòng)對(duì)裝載狀態(tài)非常敏感，在輕載工況下船舶的最大橫搖角幅值遠(yuǎn)大于滿載工況，主要原因是該裝載工況下船舶固有周期13.6 s更加接近波浪遭遇周期。

由圖6和圖7可知：

1) 在橫浪的海況下，半潛船的垂蕩運(yùn)動(dòng)響應(yīng)最大，其中滿載工況垂蕩RAO達(dá)到1.18，較為危險(xiǎn)，其他浪向垂蕩性能相對(duì)較好。

2) 在滿載工況下垂蕩運(yùn)動(dòng)大于輕載工況。顯然裝載狀態(tài)對(duì)垂蕩運(yùn)動(dòng)存在一定影響。

3) 在小于0.7 rad/s波浪頻率范圍內(nèi)，垂蕩受波浪高度影響較大，對(duì)波浪遭遇周期不太敏感。

3.3 試驗(yàn)與數(shù)值模擬的結(jié)果對(duì)比

3.3.1自然橫搖周期對(duì)比

模型試驗(yàn)與數(shù)值模擬的自然橫搖周期差異非常小，自然橫搖周期對(duì)比見表2。在3%以內(nèi)。該差異是船舶的尺度效應(yīng)所致，因此可忽略。

表2 自然橫搖周期對(duì)比 s

3.3.2運(yùn)動(dòng)響應(yīng)傳遞函數(shù)對(duì)比

半潛船不規(guī)則波縱搖RAO的模型試驗(yàn)結(jié)果和OCTOPUS預(yù)報(bào)結(jié)果比較見圖8。由圖8可知：

1) 半潛船縱搖對(duì)裝載狀況不太敏感。

2) OCTOPUS預(yù)報(bào)結(jié)果非常接近于在模型試驗(yàn)中測(cè)量到的縱搖RAO，故本文采用的數(shù)值方法可行，計(jì)算結(jié)果的精確度較好。

3) 半潛船模型試驗(yàn)成本較高，且試驗(yàn)數(shù)據(jù)非常有限，難于完整描述半潛船耐波性規(guī)律，尤其在高頻和低頻區(qū)間。OCTOPUS則可給出相對(duì)準(zhǔn)確的縱搖RAO高頻和低頻計(jì)算結(jié)果，有助于實(shí)際工程應(yīng)用。

圖8 縱搖RAO模型試驗(yàn)結(jié)果和OCTOPUS預(yù)報(bào)結(jié)果

3.3.3加速度預(yù)報(bào)對(duì)比

為進(jìn)一步驗(yàn)證在不同工況下船舶運(yùn)動(dòng)響應(yīng)的預(yù)報(bào)精度，選擇半潛船駕駛臺(tái)和主甲板上2個(gè)參考點(diǎn)進(jìn)行觀測(cè)，坐標(biāo)分別為駕駛臺(tái)：X=198.40 m，Y=0.00 m，Z=42 m；主甲板：X=0.00 m，Y=21.50 m，Z=13.00 m；圖9～圖12分別給出了輕載工況(吃水=8 m)，Hs=6.50 m,Tp=12.50 s在不同浪向情況下的駕駛臺(tái)和主甲板參考點(diǎn)數(shù)值模擬和模型試驗(yàn)的垂向加速度和橫向加速度值對(duì)比結(jié)果,見圖9～圖12。

圖9 駕駛臺(tái)處垂向加速度值對(duì)比結(jié)果圖10 駕駛臺(tái)處橫向加速度值對(duì)比結(jié)果

圖11 主甲板處垂向加速度值對(duì)比結(jié)果圖12 主甲板處橫向加速度值對(duì)比結(jié)果

由圖9～圖12可知:駕駛臺(tái)和主甲板參考點(diǎn)的垂向加速度和橫向加速度值數(shù)值計(jì)算曲線與模型試驗(yàn)曲線的趨勢(shì)基本吻合，數(shù)值模擬結(jié)果略為保守:在頂浪和艏斜浪的海況下，駕駛臺(tái)和主甲板處垂向加速度下相對(duì)較大，表明縱搖運(yùn)動(dòng)起著重要的作用，垂向加速度和浪向有關(guān)；在艉斜浪和橫浪的海況下，駕駛臺(tái)處橫向加速度最大，橫向加速度同時(shí)受浪向和所處位置的影響。半潛船在實(shí)際綁扎分析時(shí)需要全面考慮加速度處于最不利的情況。

3.3.4實(shí)船試驗(yàn)與數(shù)值模擬橫搖運(yùn)動(dòng)對(duì)比分析

對(duì)在某海域?qū)嶋H營(yíng)運(yùn)中的半潛船進(jìn)行實(shí)船耐波性試驗(yàn), 實(shí)測(cè)該船航行過(guò)程中的橫搖運(yùn)動(dòng)特性，并與OCTOPUS數(shù)值模擬結(jié)果進(jìn)行對(duì)比，從而檢驗(yàn)數(shù)值模擬方法的可靠性和實(shí)用性。兩者對(duì)比結(jié)果見圖13。

由圖13可知：實(shí)船試驗(yàn)的橫搖RAO值與數(shù)值模擬在同等海況下的計(jì)算值是比較接近和吻合的，數(shù)值模擬結(jié)果具有較高的準(zhǔn)確性，可用于實(shí)船耐波性預(yù)報(bào)。但兩者也存在一定差別，主要原因是試驗(yàn)時(shí)的海洋環(huán)境條件存在一定的不確定因素，而數(shù)值模擬的數(shù)學(xué)模型基于一定假設(shè)條件，無(wú)法完全真實(shí)模擬實(shí)際海況。在實(shí)際工程應(yīng)用中，數(shù)值模擬應(yīng)考慮一定安全裕度。

4 結(jié)束語(yǔ)

本文結(jié)合實(shí)例分析，利用試驗(yàn)方法和數(shù)值模擬方法研究某5萬(wàn)噸級(jí)半潛船的耐波性，得出以下結(jié)論，可對(duì)同類型半潛船的開發(fā)設(shè)計(jì)提供一定的參考：

1) 半潛船橫向加速度與浪向、所處半潛船垂向位置相關(guān)。半潛船垂向加速度則與浪向相關(guān)。

2) 半潛船橫搖、垂向和縱搖運(yùn)動(dòng)響應(yīng)與裝載狀態(tài)和浪向相關(guān)。在輕載工況下的船舶橫搖自然周期接近于一般航行海域的波浪遭遇周期，半潛船最大橫搖角發(fā)生在輕載艉斜浪工況；半潛船最大垂蕩運(yùn)動(dòng)發(fā)生在滿載橫浪工況。這證明在大風(fēng)浪運(yùn)載航行的情況下，半潛船應(yīng)該盡可能避免頂浪和橫浪航行，

根據(jù)實(shí)際營(yíng)運(yùn)經(jīng)驗(yàn)，盡可能保持船舶與浪向在20°夾角的情況下航行。

3) OCTOPUS數(shù)值模擬結(jié)果的分布規(guī)律與模型試驗(yàn)、實(shí)船試驗(yàn)相符合，且能夠提供更為完整的預(yù)報(bào)結(jié)果，說(shuō)明數(shù)值模擬方法對(duì)于計(jì)算半潛船耐波性是一種有效的工具，可作為半潛船的貨物綁扎分析主要技術(shù)手段。