馮乾棟，胡嘉駿，張 凡，耿彥超

(中國船舶科學(xué)研究中心，江蘇 無錫 214082)

破損艦船在波浪中的響應(yīng)預(yù)報(bào)與試驗(yàn)研究

馮乾棟，胡嘉駿，張 凡，耿彥超

(中國船舶科學(xué)研究中心，江蘇 無錫 214082)

為研究艦船破損后在波浪中的響應(yīng)，開展艦船在完整正浮狀態(tài)以及破損后傾斜狀態(tài)的模型試驗(yàn)?；谌S線性勢流理論，使用三維頻域方法對模型破損前后的響應(yīng)做數(shù)值預(yù)報(bào)。與試驗(yàn)結(jié)果相比，數(shù)值預(yù)報(bào)結(jié)果符合較好，驗(yàn)證了預(yù)報(bào)方法的合理性。最后，通過與該艦船模型破損前的響應(yīng)結(jié)果對比，研究分析破損對艦船在波浪中響應(yīng)的影響。

破損艦船；模型試驗(yàn)；三維頻域方法

0 引 言

雖然航海技術(shù)不斷進(jìn)步，由于擱淺、碰撞等意外情況，艦船仍然難免面臨著出現(xiàn)結(jié)構(gòu)破損導(dǎo)致海水淹艙的危險(xiǎn)。船舶破損艙室進(jìn)水后，船身不僅吃水增加，并且會失去正浮產(chǎn)生傾斜，如果再遭遇惡劣海況，船身結(jié)構(gòu)極可能會進(jìn)一步破損，從而造成海難。因此，對于艦船破損后在波浪中的響應(yīng)研究很有現(xiàn)實(shí)意義。

自從“愛沙尼亞”號滾裝船海難后，不少國外學(xué)者對艦船破損后的響應(yīng)預(yù)報(bào)以及試驗(yàn)展開研究。Emin Korkut等[1-2]分別制作了完整以及破損的滾裝船模型，通過水池試驗(yàn)比較了滾裝船破損前后的運(yùn)動響應(yīng)與波浪載荷。Chan等[3-4]使用二維切片理論，在時(shí)域內(nèi)考慮了船體大幅值運(yùn)動時(shí)劇烈變化的濕表面影響，預(yù)報(bào)了破損滾裝船在規(guī)則斜浪中的非線性運(yùn)動與波浪載荷。Lee等[5]考慮船體多種破損情況，使用了源匯分布法求解了破損后非對稱剖面的二維輻射和繞射流場，利用線性切片理論預(yù)報(bào)了船舶破損后的波浪載荷，并展開了模型試驗(yàn)與數(shù)值預(yù)報(bào)做了對比。

在國內(nèi)，吳明遠(yuǎn)、青兆熹[6-7]采用三維方法在頻域內(nèi)計(jì)算了船體周圍的速度勢分布，預(yù)報(bào)了船舶橫傾后的6自由度運(yùn)動與波浪載荷。這些研究都獲得了有益的結(jié)果，但是都尚未開展模型試驗(yàn)，其理論預(yù)報(bào)結(jié)果缺乏模型試驗(yàn)的驗(yàn)證。

本文通過開展模型試驗(yàn)與數(shù)值預(yù)報(bào)，研究艦船發(fā)生第一類艙室破損，由艙室淹水導(dǎo)致的船身傾斜與吃水增加對船體運(yùn)動與波浪載荷的影響。在模型試驗(yàn)方面，本文分別開展完整正浮與破損后橫傾的艦船模型對比試驗(yàn)，測量該模型在2種狀態(tài)下在頂浪規(guī)則波中的響應(yīng)。另外，本文基于三維勢流理論，使用三維奇點(diǎn)分布法求解非定常流場，進(jìn)而對模型的運(yùn)動響應(yīng)與波浪載荷進(jìn)行數(shù)值預(yù)報(bào)，并將數(shù)值預(yù)報(bào)同模型試驗(yàn)結(jié)果進(jìn)行比較驗(yàn)證。

1 破損艦船三維線性頻域預(yù)報(bào)方法

1.1 三維速度勢求解

(1)

式中：ωe為遭遇頻率；n為船體表面單位法矢量，在細(xì)長船體假設(shè)下，本文取簡化形式的mj，即(m1,m2,m3)=(0,0,0)， (m4,m5,m6)=(0,n3,-n2)。

(2)

通過簡化，帶有航速的自由面條件與零航速情況形式相同，只需以遭遇頻率ωE代替自然頻率ω，從而使得有航速速度勢的求解可轉(zhuǎn)換為零航速問題。

簡化后的船體三維速度勢可以通過三維奇點(diǎn)分布法來求解：采用面元法將連續(xù)物面(船體平均濕表面)離散，使用零航速的三維自由面Green函數(shù)建立分布源方程組，求解得到面元中心處的源強(qiáng)與Green函數(shù)后，使用分布源積分方法即可求解出船體周圍的三維速度勢。

1.2 破損艦船在規(guī)則波中的運(yùn)動與波浪載荷計(jì)算

假設(shè)艦船發(fā)生第一類艙室破損，即艙的頂部位于水線以下，船體破損后海水灌滿整個(gè)艙室，因此本文不考慮淹艙水晃蕩對于船體運(yùn)動與載荷的影響。

(3)

式中，η為船體的穩(wěn)態(tài)運(yùn)動響應(yīng)；M為質(zhì)量矩陣；C為靜水回復(fù)力矩陣；A和B分別為附加質(zhì)量矩陣以及附加阻尼矩陣；fw為波浪激勵(lì)力。

由于船舶破損后具有一定橫傾，船體質(zhì)量分布、靜水浮力不再關(guān)于中縱剖面對稱，所以破損船運(yùn)動方程中的質(zhì)量矩陣M與靜水回復(fù)力矩陣C與破損前有所不同。

(4)

完好的船舶由于水線面左右對稱，靜水恢復(fù)矩陣中有5項(xiàng)不為0，而非正浮破損船的水線面不再對稱，其靜水回復(fù)力矩陣C中有9項(xiàng)不為0：

(5)

式中：Awp為水線面面積;Sx為水面線對x軸的靜矩;Sy為水面線對y軸的靜矩；Jxy為水面線對原點(diǎn)的靜矩;hx為橫穩(wěn)性高;hy為縱穩(wěn)性高。

通過式(3)求解得到六自由度運(yùn)動后，船中剖面處的波浪載荷Q可以通過下式計(jì)算：

(6)

式中p為各面元處的脈動壓力。

2 船模試驗(yàn)

本文分別對完整正浮與破損后橫傾的艦船模型進(jìn)行了試驗(yàn)，測量了該模型2種狀態(tài)下在頂浪規(guī)則波中的響應(yīng)。船模試驗(yàn)在中國船舶科學(xué)研究中心05耐波性水池中進(jìn)行，該水池長69 m,寬46 m，水深4 m，在長、短兩邊布置有搖板式造波機(jī)。試驗(yàn)采用某深V型艦船模型，其主要參數(shù)如表1所示。模型為分段龍骨梁型玻璃鋼自航船模，共分4段，用測量梁連接。

表1 船模主要參數(shù)Tab.1 Principal dimensions of ship model

為了開展破損后橫傾的船體模型試驗(yàn)，本次試驗(yàn)通過使用壓載來模擬艦船2處破損后的淹艙水。在重心前后60 cm處分別放置7 kg壓載，其在橫剖面內(nèi)的位置如圖1所示，測得船體產(chǎn)生橫傾7.8°(左傾)。圖2為水池中的試驗(yàn)?zāi)Ｐ汀?/p>

試驗(yàn)都在規(guī)則波中進(jìn)行，波高取82.5 mm，使用的波長范圍：λ/L=0.4～2.5， 船模航速V=0，0.813 m/s，對應(yīng)的Fr=0，0.128。在船模重心處放置陀螺儀，用于測量船體垂蕩與縱搖響應(yīng)。在水池測量小車上，位于模型前約1.5 m處安裝1臺伺服式浪高儀，用來測量試驗(yàn)時(shí)模型遭遇的波浪。

圖1 壓載在橫剖面內(nèi)的位置Fig.1 Position of load in crosssection

圖2 水池中的試驗(yàn)?zāi)Ｐ虵ig.2 Test model in the pool

3 破損前后響應(yīng)對比分析

根據(jù)三維線性頻域理論對該船模破損前后在規(guī)則波中的響應(yīng)進(jìn)行計(jì)算，并與試驗(yàn)結(jié)果比較分析。圖3和圖4分別為破損前后的縱搖、垂蕩響應(yīng)，圖5為破損前后的垂向波浪彎矩。

圖3 縱搖RAO比較Fig.3 Comparison of pitch′s RAO

圖4 垂蕩RAO比較Fig.4 Comparison of heave′s RAO

圖5 船中垂向波浪彎矩RAO比較Fig.5 Comparison of V.B.M′s RAO in mid-ship

由圖3～圖5可看出:

1)試驗(yàn)結(jié)果與理論預(yù)報(bào)結(jié)果均顯示：當(dāng)船體發(fā)生破損后，船體縱搖響應(yīng)減小，在V=0.813 m/s時(shí)更為明顯；船體垂蕩運(yùn)動基本不變。

2)與完整正浮時(shí)相比，破損后的船體垂向波浪彎矩在峰值附近有較明顯的增大，試驗(yàn)結(jié)果在V=0.813 m/s時(shí)峰值增加8.32%。

3)理論預(yù)報(bào)能夠正確反映出船體破損后在波浪中的響應(yīng)變化，且與試驗(yàn)結(jié)果符合較好，但理論預(yù)報(bào)的縱搖響應(yīng)與波浪彎矩峰值比試驗(yàn)值偏大。

4 結(jié) 語

為了研究艦船破損前后在波浪中響應(yīng)的變化，本文開展了船模水池試驗(yàn)并做了理論預(yù)報(bào)。通過對比分析，理論預(yù)報(bào)結(jié)果與試驗(yàn)結(jié)果符合較好，驗(yàn)證了理論計(jì)算方法。理論預(yù)報(bào)與試驗(yàn)結(jié)果均顯示：該艦船發(fā)生破損后，縱搖響應(yīng)減小，垂蕩響應(yīng)基本不變，垂向波浪彎矩在峰值附近有較明顯的增大。垂向波浪彎矩的增加對破損后的艦船結(jié)構(gòu)十分不利，有可能會使得艦船出現(xiàn)進(jìn)一步的破損。通過開展對艦船破損前后在波浪中響應(yīng)的變化的研究，可為對艦船破損后生命力的研究提供依據(jù)。

[1] KORKUT E,ATLAR M,INCECIK A.An experimental study of global loads acting on an intact and damaged Ro-Ro ship model[J].Ocean Engineering,2005,32:1370-1403.

[2] KORKUT E,ATLAR M,INCECIK A.An experimental study of motion behavior with an intact and damaged Ro-Ro ship model[J].Ocean Engineering, 2004,31:483-512.

[3] CHAN HS,ATLAR M,INCECIK A.Large-amplitude motion responses of a Ro-Ro ship to regular oblique waves in intact and damaged conditions[J].Marine Science and Technology,2002,7:91-99.

[4] CHAN HS,ATLAR M,INCECIK A.Global wave loads on intact and damaged Ro-Ro ships in regular oblique waves[J].Marine Structures,2003,16:323-344.

[5] LEE YW,CHAN HS,PU YC,et al.Global wave loads on a damaged ship[R].Lloyd′s Register,2012.

[6] 吳明遠(yuǎn).破損船舶運(yùn)動與波浪載荷的三維計(jì)算方法研究[D].哈爾濱：哈爾濱工程大學(xué),2007.

WU Ming-yuan.Study of 3D calculation method of damaged ships' motion and wave load[D].Harbin:Harbin Engineering University.2007.

[7] 青兆熹.破損船舶波浪運(yùn)動載荷計(jì)算方法研究[D].哈爾濱：哈爾濱工程大學(xué),2009.

QING Zhao-xi.Study of calculation method of damaged ships′ motion and wave load[D].Harbin: Harbin Engineering University.2007.

[8] 戴仰山,沈進(jìn)威.船舶波浪載荷[M].北京：國防工業(yè)出版社,2007.

DAI Yang-shan,SHEN Jin-wei.Ship wave loads[M].Beijing:National Defense Industry Press,2007.

[9] 孫伯起.低航速物體在波浪上的運(yùn)動理論和波浪阻尼計(jì)算[R].無錫：CSSRC研究報(bào)告,1990.

SUN Bo-qi.Ship′s motion theory and wave resistance calculation with low speed[R].Wuxi:CSSRC Research Report,1990.

[10] 鄒元杰,段文洋,任慧龍，等.船舶在波浪中脈動壓力預(yù)報(bào)的三維方法[J].哈爾濱工程大學(xué)學(xué)報(bào),2002,23(1):20-25.

ZOU Yuan-jie, DUAN Wen-yang,REN Hui-long,et al.Three-dimensional method for prediction of oscillation pressure on ship waves[J].Journal of Harbin Engineering University,2002,23(1):20-25.

Prediction and model test of responses of damaged vessels in waves

FENG Qian-dong,HU Jia-jun,ZHANG Fan,GENG Yan-chao

(China Ship Scientific Research Center,Wuxi 214082,China)

Aimed at studying the responses of damaged vessels in waves, intact and damaged vessel models′ tests were carried out. Numerical predictions to the responses of intact and damaged vessels in waves were made by using 3D linear potential theory in frequency domain. Numerical predictions match the results of models tests well, verifying the prediction theory. Finally, the influence of damage to the vessel′s responses in waves was analyzed by comparing the responses with the intact vessel′s.

damaged vessels;model test;three-dimensional method in frequency domain

2013-12-27;

2014-02-19

馮乾棟(1989-)，男，碩士研究生，研究方向?yàn)榇芭c海洋工程環(huán)境載荷。

U661.32

A

1672-7649(2014)11-0029-04

10.3404/j.issn.1672-7649.2014.11.006