基于CFD的艦船遭遇畸形波數(shù)值模擬的驗證研究

蔡 烽，張本輝，吳 明，楊 波，王 驍，石愛國

（海軍大連艦艇學(xué)院 航海系，遼寧 大連 116018）

0 引 言

畸形波（freak wave）是海洋中高且陡的大波，其持續(xù)時間很短,但出現(xiàn)的偶然性和巨大的破壞性，對船舶航運和海洋工程結(jié)構(gòu)物等極具威脅，已引發(fā)多起海上事故[1]，因此畸形波越來越引起人們的關(guān)注，它的發(fā)生機理及工程應(yīng)用問題已成為當(dāng)前海洋物理學(xué)界和船舶水動力學(xué)界的一個研究熱點問題[2-5]。海軍艦艇需要在各種海況下航行，海浪環(huán)境深刻地影響艦艇的作戰(zhàn)效能，因此如何更加準(zhǔn)確地把握遠(yuǎn)洋海區(qū)風(fēng)浪環(huán)境的非線性規(guī)律和特點，對畸形波進行監(jiān)測和預(yù)報，進一步計算艦船遭遇畸形波時的耐波性能、載荷響應(yīng)以及預(yù)報其是否存在傾覆等風(fēng)險，對完善現(xiàn)有的船舶設(shè)計和操縱規(guī)范、保障艦艇航行安全以及提升海軍戰(zhàn)斗力具有顯著作用。

關(guān)于畸形波的數(shù)值模擬方法，從考慮波浪調(diào)制不穩(wěn)定性[6-7]的非線性波動方程（例如三階、四階非線性薛定諤方程）出發(fā)，可以研究畸形波的發(fā)生機理，但是據(jù)此很難把握畸形波發(fā)生的時空條件，因此不宜用于開展艦船遭遇畸形波相關(guān)的耐波性試驗?；贚onguet-Higgins模型是實驗室模擬產(chǎn)生畸形波的有效手段和常用方法，簡單易用，可以實現(xiàn)畸形波在實驗室的定時定點生成[8]。黃國興[9]采用人工干預(yù)組成波隨機初相位的方法，使部分組成波初相位相同，可以得到包含畸形波的波列，但模擬的效率比較低，且不能控制畸形波的生成時間和生成地點。Kriebel[10]采用一個基本隨機波列和一個瞬態(tài)波列線性疊加的雙波列疊加模型模擬了畸形波；裴玉國[11]采用三波列疊加模型優(yōu)化了畸形波的模擬。這兩種模擬方法都可以實現(xiàn)畸形波的定時定點生成，但基于瞬態(tài)波列與隨機波列的疊加，瞬態(tài)波列的能量所占的比例會影響整個模擬波列的譜的結(jié)構(gòu)。對于大尺度的畸形波，采用瞬態(tài)波列方法將影響波浪序列的統(tǒng)計特性。劉贊強[12]采用改進的相位調(diào)制方法來模擬畸形波，既滿足波浪序列的統(tǒng)計特性又可保持目標(biāo)譜的結(jié)構(gòu)，且模擬效率較高，是一種有效的模擬畸形波的方法。因此，本文提出了計及航速航向的相位調(diào)制方法來模擬畸形波（頂浪航行時僅需要考慮航速的影響），進而利用CFD（計算流體力學(xué)）實現(xiàn)了艦船頂浪狀態(tài)下定時定點遭遇畸形波，并與Bennett等人的水池試驗數(shù)據(jù)對比進行驗證，吻合較好。

1 英國南安普頓大學(xué)艦船遭遇畸形波水池實驗

對于畸形波的數(shù)值模擬、演化規(guī)律及其對近岸結(jié)構(gòu)物響應(yīng)的研究，國內(nèi)外學(xué)者做了大量的工作，但對于艦船遭遇畸形波的研究還相對較少。國外的研究艦船遭遇畸形波的團隊主要有兩個：德國工業(yè)大學(xué)的Clauss[13]和英國南安普頓大學(xué)的Bennett[14]，國內(nèi)則偏重于理論的研究[15]，本文利用Bennett等人的水池實驗結(jié)果對數(shù)值模擬結(jié)果進行驗證，下面對水池試驗相關(guān)情況進行簡要介紹。

1.1 水池試驗條件

拖曳水池的長、寬、深分別為60 m×3.7 m×1.86 m，最大可用拖曳速度為4.5 m/s?；尾ň劢沟奈恢迷诰嚯x造波機25 m處，可生成兩維非規(guī)則海浪的周期范圍為T≥0.7 s。造波機的槳板是由用戶自定義輸入技術(shù)所控制的；獲得的波浪測量結(jié)果顯示消波區(qū)的反射低于10%。

1.2 船模相關(guān)參數(shù)

圖1 英國利爾德級護衛(wèi)艦的型線圖和船模示意圖Fig.1 The body plan and the experimental set-up of leander class frigate hull

水池試驗采用的艦船為英國利安德級護衛(wèi)艦，曾是20世紀(jì)60～80年代英國海軍的主力戰(zhàn)艦。Andrew和Lloyd在1981年時曾用此船模研究過搖蕩運動和甲板上浪。船?？s尺比為1：43.62，總長為2.6 m，型線圖和船模示意如圖1所示，主要的船型參數(shù)如表1所示。

表1 利安德級護衛(wèi)艦的主要型值數(shù)據(jù)Tab.1 Principal particulars of Leander class frigate

1.3 基于相位優(yōu)化技術(shù)生成畸形波

Bennett等人采用了相位優(yōu)化技術(shù)生成畸形波，相位優(yōu)化技術(shù)由Clauss等人[16]提出，波浪的計算機控制優(yōu)化過程如圖2所示。

給定目標(biāo)畸形波的聚焦時間、上跨零點周期、最大波高以及峰前波陡，構(gòu)造約束條件和目標(biāo)函數(shù)；初始的相位是隨機生成的，利用序列二次規(guī)劃算法對相位進行修正，波浪序列自動生成、測量、評估和修正，直到滿足目標(biāo)值。該相位優(yōu)化技術(shù)的優(yōu)勢是利用物理的水池可以將波浪非線性效果考慮在內(nèi)，實際上同時完成自我驗證的過程。

圖2 計算機控制的波浪生成優(yōu)化過程Fig.2 Computer controlled optimization of waves

1.4 Bennett對畸形波的定義

Bennett等人采用的畸形波定義主要由兩個條件組成：

式中：HR為最大波高，ηR為峰值，Hs為有義波高，AI被稱為畸形指數(shù)，CI被稱為波峰指數(shù)。

2 畸形波海浪環(huán)境生成的數(shù)值模擬研究

采用CFD方法模擬艦船在波浪中的運動時，假設(shè)在固定坐標(biāo)系中，航速為U0，浪舷角為χ；則長峰非規(guī)則波的波高方程為[17]：

設(shè)在位置 x=xc、y=yc，時刻 t=tc時生成畸形波，調(diào)制初相 θi使部分（或者全部）組成波在 x=xc、y=yc，t=tc時 ηi（xc，yc，tc）為正，則在此疊加的波高會增大。 令組成波數(shù)M=M1+M2，則（3）式可以寫為：

在此令后M2個組成波的合成波波面η2（x,y, t）在預(yù)定位置處聚焦出現(xiàn)大波，需要調(diào)制后M2個組成波的初相位 θi，使 ηi（xc，yc，tc）。

（1） 當(dāng) ki（xccos χ+ycsin χ-U0cos χtc）-ωitc＜0 時，令整數(shù) N=int［（ki（xccos χ+ycsin χ-U0cos χtc）-ωitc）/2 π ］，此時 N＜0，（6）式可以寫為：

調(diào)制 θi，使-π/2＜ki（xccos χ+ycsin χ-U0cos χtc）-ωitc-2Nπ+θi＜π/2，這樣 cos（ ki（xccos χ+ycsin χ-U0cos χtc）-ωitc-2Nπ+θi）＞0，此時 ηi（xc，yc，tc）＞0，η2（xc，yc，tc）＞0，由于-2π＜ki（xccos χ+ycsin χ-U0cos χtc）-ωitc-2Nπ＜0，θi在下述區(qū)間隨機取值：

（2） 當(dāng) ki（xccos χ+ycsin χ-U0cos χtc）-ωitc≥0 時，令整數(shù) N=int［ （ki（xccos χ+ycsin χ-U0cos χtc）-ωitc）/2 π ］，此時 N≥0，（6）式可以寫為

調(diào)制 θi，使-π/2＜ki（xccos χ+ycsin χ-U0cos χtc）-ωitc-（2N+1 ） π+θi＜π/2，這樣 cos（ xccos χ+ycsin χ-U0cos χtc）-ωitc-（2N+ 1） π+θi＞0，此時 ηi（xc，yc，tc）＞0，η2（xc，yc，tc）＞0，θi的確定方法與情況（1）中所述的相同，在此不再贅述。

艦船頂浪航行時，浪舷角為χ=0；則方程（3）可簡化為：

基于此數(shù)學(xué)模型，可以對艦船定時定點遭遇畸形波的波浪環(huán)境進行數(shù)值模擬，為利用CFD方法進行相關(guān)研究奠定基礎(chǔ)。

對于CFD而言，采用Bennett等人的相位優(yōu)化技術(shù)還存在一定的技術(shù)困難。即使初始條件相同，相位優(yōu)化解并不依賴于初始隨機相位，優(yōu)化之后的相位選擇并不唯一，海浪的隨機特性并沒有完全丟失，數(shù)值模擬試驗的相位無法與水池實驗的相位完全一致，導(dǎo)致生成的畸形波海浪環(huán)境必然存在著差異。因此，盡可能地對相位進行篩選，使各項波浪指數(shù)與水池實驗值接近，以期與水池試驗具有可比性。另外，船體周圍的波浪場實際上是由入射波（即畸形波）、艦船反射波和輻射波構(gòu)成，本節(jié)主要考慮畸形波的定時定點生成，在加載船模進行耐波性相關(guān)數(shù)值計算時，F(xiàn)luent軟件可以自動將艦船與流場的流固耦合運動考慮在內(nèi)。

3 低速頂浪航行時數(shù)值模擬和水池試驗對比

3.1 基于CFD進行船模水池實驗的方案

利用Fluent前處理軟件Gambit生成的利安德級護衛(wèi)艦船體如圖3所示。

圖3 利安德級護衛(wèi)艦船體Fig.3 Leander class frigate hull

Bennett等人通過控制造波機的波浪生成以及拖車的運動，確保船模的中部剛好遭遇畸形波波峰，本文在利用CFD進行相關(guān)的數(shù)值模擬時，在預(yù)定時刻令聚焦位置剛好位于船模中部，從而實現(xiàn)了艦船定時定點遭遇畸形波，并在此記錄波形。

3.1.1 計算域

計算域設(shè)置為長方體，如圖4所示，計算域尺寸為：入口距船首1L，出口距船尾2L，頂部邊界距水線0.5L，底部邊界距水線1.5L，左、右邊界距船中縱剖面1L。

圖4 頂浪流場的計算域Fig.4 The computational domain of head-sea filed

3.1.2 網(wǎng)格劃分

為充分發(fā)揮結(jié)構(gòu)及非結(jié)構(gòu)網(wǎng)格優(yōu)勢，對計算域分為船體近流場區(qū)及遠(yuǎn)流場區(qū)進行處理，近流區(qū)域采用非結(jié)構(gòu)網(wǎng)格進行網(wǎng)格生成，船體面網(wǎng)格采用三角形網(wǎng)格，從船體表面網(wǎng)格以一定比例外推而生成近密外疏的非結(jié)構(gòu)四面體網(wǎng)格。遠(yuǎn)流場區(qū)完全布設(shè)結(jié)構(gòu)網(wǎng)格，由自由面向上下兩底邊漸疏，在保證網(wǎng)格質(zhì)量的同時保持合適的網(wǎng)格數(shù)量，整個計算域網(wǎng)格劃分效果如圖5所示。

網(wǎng)格總數(shù)為1 381 643，網(wǎng)格質(zhì)量分布如表2所示。

3.1.3 邊界條件設(shè)置

計算域及邊界條件名稱如圖6所示。

具體設(shè)置如下：

入口邊界（in）：速度入口（velocity-inlet）條件，設(shè)置入口處流體的速度、入口處水的體積分?jǐn)?shù)、湍動能k，耗散率ε；

圖5 網(wǎng)格劃分效果Fig.5 Meshing effect

表2 網(wǎng)格質(zhì)量評價Tab.2 The evaluation of mesh quality

出口邊界（out）：壓力出口（pressure-outlet）條件，設(shè)定靜壓力、底面位置和自由面高度，回流的湍動能k及耗散率ε則采用Fluent的推薦值；

圖6 計算域及邊界名稱Fig.6 Computational domain and name of boundary conditions

上下邊界（top、bot）—速度入口，給定三個方向流速（u=U0、v=w=0）及水的體積分?jǐn)?shù)（0）；

左右邊界（port、stab）—滑移的壁面；

船體（ship）—有剪切力無滑移的壁面。

3.1.4 數(shù)值模擬的初始條件

數(shù)值波浪水池采用的靶譜為Jonswap譜，其水池試驗參數(shù)如表3所示。

表3 水池試驗的主要參數(shù)Tab.3 The main parameters in experiment

聚焦時刻為tc=5 s，聚焦位置位于船舶的重心處即xc=0 m。

3.2 數(shù)值模擬結(jié)果及分析

3.2.1 利安德護衛(wèi)艦頂浪航行流場分析

Bennett等人利用相位優(yōu)化技術(shù)在水池中定時定點生成了畸形波，為檢驗數(shù)據(jù)采集系統(tǒng)的不確定度，特進行了三次重復(fù)實驗，其波面的測量結(jié)果如圖7所示。

圖7 測量過程中波面的不確定度Fig.7 Wave profiles showing uncertainty in experimental measurements

由圖 7可知，畸形波的聚焦時間為tc=6.2 s，AI接近2.36，CI接近1.37，說明利用相位優(yōu)化技術(shù)生成可以定時定點生成畸形波，另外，測量過程的相對誤差較小。

為了分析艦船反射波、輻射波對生成的畸形波波浪場的影響，不加載艦船的波浪場中，在船模中心位置設(shè)置浪高儀進行時歷監(jiān)測，并作為參考值。在圖5的三維計算域中，加載了船模之后，分別在距離船模中心橫向距離y=0.5 m、1.5 m、2.5 m處設(shè)置浪高儀，進行波浪時歷監(jiān)測，并將測量結(jié)果同參考值進行比對，對圖8所包含的極值大波進行特征統(tǒng)計，其畸形特征參數(shù)如表4所示。

由圖 8和表4可知，加載了船舶搖蕩之后，對入射波浪場（畸形波）產(chǎn)生了干擾影響，與參考值對比可知，這種干擾影響相對較小。干擾主要是由船舶在行進過程中產(chǎn)生的反射波和繞射波造成的，隨著遠(yuǎn)離船模，干擾作用在不斷減小。

圖8 不同橫向位置處測得的波浪時歷Fig.8 Wave profiles measured in different transverse locations

表4 不同橫向位置處測得波浪時歷的畸形參數(shù)比對Tab.4 The freak-wave parameter of wave profiles measured in different transverse locations

3.2.2 利安德護衛(wèi)艦頂浪航行搖蕩運動數(shù)值模擬

當(dāng)船舶流場發(fā)展比較充分的時候，取船模搖蕩運動比較穩(wěn)定后的某一時刻（t=1.5 s）作為記錄起點，數(shù)值模擬船模在包含畸形波的非規(guī)則波浪中做縱搖和垂蕩運動，當(dāng)艦船遭遇畸形波時的時刻（t=4.5～6 s，每隔0.25 s取一幅圖片）的瞬時態(tài)勢如圖 9所示。

圖9 U0=0.6 m/s時艦船頂浪遭遇畸形波的過程Fig.9 The process of ship encountering with freak waves at forward speed U0=0.6 m/s

由圖 9可知，艦船在包含畸形波的非規(guī)則波浪中頂浪航行時，當(dāng)遭遇的波浪比較大時，會出現(xiàn)砰擊、上浪等強非線性現(xiàn)象，會對船舶產(chǎn)生一系列不利于航行的搖蕩響應(yīng)，甚至危及船舶的航行安全。

圖9中利安德級護衛(wèi)艦船模在包含畸形波的非規(guī)則波浪航行時，測得的縱搖時歷和垂蕩時歷分別如圖10和圖11所示。

由圖10和圖11可知，當(dāng)艦船遭遇比較大的波浪時，垂蕩和縱搖值相對較大，特別遭遇畸形波時，更加明顯，與圖9對應(yīng)可知，產(chǎn)生了砰擊和上浪現(xiàn)象。

采用下跨零點法對U0=0.6 m/s工況下所采集的縱搖時歷和垂蕩時歷進行處理，求得其最大縱搖值和最大垂蕩值（與水池實驗一致，此處采用的是谷—峰最大值）如表5所示。

圖10 U0=0.6 m/s工況下所采集到的縱搖時歷Fig.10 The pitch profile obtained at forward speed U0=0.6 m/s

圖11 U0=0.6 m/s工況下所采集到的垂蕩時歷Fig.11 The heave profile obtained at forward speed U0=0.6 m/s

表5 低速頂浪航態(tài)下數(shù)值模擬值與水池實驗值的比對Tab.5 The comparison between CFD numerical simulation values and wave-tank experimental values at forward low-speed

由表 5可知，CFD數(shù)值模擬的初始條件同水池試驗基本一致，則獲取的ξmax（m）和θmax（deg）比較接近，誤差值可以接受，證明了低速頂浪航行情況下本文數(shù)值模擬的有效性。

4 中速頂浪航行數(shù)值模擬和水池試驗對比

當(dāng)艦船以不同的航速遭遇畸形波時，其運動響應(yīng)又會有何不同，結(jié)合Bennett等人的水池實驗情況，本節(jié)設(shè)計了另一航速的試驗，試驗的主要參數(shù)如表6所示。

表6 試驗的主要參數(shù)Tab.6 The main parameters in experiment

在圖5的三維計算域中，距離船模中心橫向距離y=2.5 m處設(shè)置浪高儀進行波浪時歷監(jiān)測，如圖12所示。

圖12 橫向位置y=2.5 m處測得的波浪時歷Fig.12 Wave profiles measured in transverse location y=2.5 m

利用下跨零點法對圖12中的波浪時歷所包含的極值大波進行特征統(tǒng)計，將畸形參數(shù)與水池實驗給定值進行比對，如所表7示。

表7 水池試驗和CFD數(shù)值模擬參數(shù)的比對Tab.7 The comparison between CFD numerical simulation and wave-tank experimental values

在U0=1.4 m/s的工況下，艦船遭遇包含畸形波的非規(guī)則海浪所對應(yīng)的縱搖和垂蕩時歷分別如圖13和圖14所示。

圖13 U0=1.4 m/s工況下所采集到的縱搖時歷Fig.13 The pitch profile obtained at forward speed U0=1.4 m/s

圖14 U0=1.4 m/s工況下所采集到的垂蕩時歷Fig.14 The heave profile obtained at forward speed U0=1.4 m/s

采用下跨零點法對U0=1.4 m/s工況下所采集的縱搖時歷和垂蕩時歷進行畸形參數(shù)統(tǒng)計，求得其最大縱搖值和最大垂蕩值（與水池實驗一致，此處采用的是谷-峰最大值）如表8所示。

表8 中速頂浪航態(tài)下CFD數(shù)值模擬值與水池實驗值的比對Tab.8 The comparison between CFD numerical simulation and wave-tank experimental values at forward middle-speed

由表 8可知，CFD數(shù)值模擬的初始條件同水池試驗一致，則獲取的ξmax（m）和θmax（deg）比較接近，誤差值可以接受，證明了中速頂浪航行情況下本文數(shù)值模擬的有效性。

將表5和表8對比可知，隨著航速的增大，垂蕩值增大，縱搖值減少，CFD數(shù)值模擬得出的趨勢與水池實驗的趨勢相同，得出這樣的結(jié)論與艦船在畸形波中的搖蕩運動有關(guān)。當(dāng)艦船遭遇畸形波時，船頭有可能被掩埋在波高較大的海浪中，當(dāng)船從波峰中穿過時會出現(xiàn)明顯的上浪現(xiàn)象。當(dāng)高速遭遇到海浪時，船頭掩埋和甲板上浪會更加嚴(yán)重；在一定程度上使縱搖運動減弱；然而，垂蕩運動并不因為上浪嚴(yán)重而受到影響，而是隨著航速的增大而增大。

5 試驗結(jié)果與船舶設(shè)計規(guī)范的對比

由于畸形波的突發(fā)特性，可以在毫無預(yù)兆的情況下出現(xiàn)，意味著艦船有隨時遭遇畸形波的可能，特別是在海況比較惡劣的情況下，導(dǎo)致海難事故時有發(fā)生，因此在船舶設(shè)計階段，就需要將最為惡劣的情況考慮在內(nèi)，以防悲劇的發(fā)生。

Bennett等人對水池試驗的搖蕩運動時歷使用中心差分方法得到垂蕩加速度和縱搖加速度，并與英國勞安德船級社的規(guī)范值進行比對，發(fā)現(xiàn)縱搖加速度遠(yuǎn)小于規(guī)范值，而垂蕩加速度則有可能大于船級社的規(guī)范值（特別是航速較高的情況下）。在利用CFD進行數(shù)值模擬時，U0=0.6 m/s和U0=1.4 m/s兩種工況下的垂蕩加速度如圖15所示。U=1.4m/s U=0.6m/s

圖15 兩種工況下的垂蕩加速度對比Fig.15 The comparison of the heave acceleration under two different work conditions

求取圖15中垂蕩加速度曲線的最大值，并將其與重力加速度進行歸一化，進而與英國勞安德船級社的規(guī)范值（2010）對比，如表9所示。

由表 9 可知，當(dāng) U0=0.6 m/s時，amax（g）＜0.39，并沒有超越船級社的規(guī)范值，但是當(dāng)U0=1.4 m/s時，amax（g）＞0.39，超越了船級社的規(guī)范值，這在船舶設(shè)計的過程中需要予以考慮。

表9 試驗得到的垂蕩加速度與船級社規(guī)范值的比對Tab.9 The comparison between heave acceleration value with Lioyd’s Register Rules

6 結(jié) 論

本文基于隨機波浪的Longuet-Higgins模型，在相位調(diào)制方法的基礎(chǔ)上考慮了航速航向的影響，調(diào)制部分組成波的初相位，可以實現(xiàn)艦船定時定點遭遇畸形波。主要可以得出以下結(jié)論：

（1）利用CFD方法對頂浪情況下艦船定時定點遭遇畸形波進行了數(shù)值模擬，并將數(shù)值模擬結(jié)果與英國南安普頓水池實驗值進行比對，吻合較好，證明了本文數(shù)值模擬的有效性。

（2）當(dāng)艦船遭遇相同的畸形波時，隨著航速的增大，垂蕩值增大，縱搖值減小，CFD數(shù)值模擬得出的趨勢與水池試驗的結(jié)論一致。

（3）當(dāng)艦船高速遭遇畸形波時，某些動力響應(yīng)值（比如垂蕩加速度）超出了船級社的規(guī)范值，有可能存在一些安全隱患，在船舶設(shè)計過程中是應(yīng)該考慮在內(nèi)的。

參 考 文 獻：

[1]董艷秋,紀(jì) 凱,黃衍順.波浪中船舶橫搖穩(wěn)性的研究[J].船舶力學(xué),1999,3(2)：1-6.Dong Yanqiu,Ji Kai,Huang Yanshun.Study on the ship stability in rolling in waves[J].Journal of Ship Mechanics,1999,3(2)：1-6.

[2]劉首華,牟 林,劉克修,等.畸形波研究的進展及存在問題[J].地球科學(xué)進展,2013,28(6)：665-673.Liu Shouhua,Mu Lin,Liu Kexiu,et al.Research advance and problems in freak waves[J].Advances in Earth Science,2013,28(6)：665-673.

[3]高志一,于福江,徐福祥,等.畸形波生成條件預(yù)報方法研究進展[J].海洋通報,2011,30(3)：351-356.Gao Zhiyi,Yü Fujiang,Xü Fuxiang,et al.Progress in the method of freak wave condition forecasting[J].Martine Science Bulletin,2011,30(3)：351-356.

[4]Cui Cheng,Zhang Ningchuan,Pei Yuguo,Liu Qianlie.Numerical study on generation and evolution of freak waves[J].Journal of Ship Mechanics,2012,12,16(12)：51-56.

[5]沈玉稿.畸形波的數(shù)值模擬及其與海洋結(jié)構(gòu)物相互作用研究[D].上海：上海交通大學(xué),2013.Shen Yügao.Numerical simulation of freak wave and the interactions between freak wave and offshore structures[D].Shanghai：Shanghai Jiao Tong University,2013.

[6]Zhang Yunqiu,Zhang Ningchuan.Numerical simulation and mechanism analysis of freak waves[J].Acta Oceanologica Sinica,2007,26(5)：116-124.

[7]Lo E,Mei C C.A numerical study of water-wave modulation based on a higher-order nonlinear Schroedinger equation[J].Journal of Fluid Mechanics,1985,150：395-416.

[8]劉贊強.畸形波模擬與其在核電取水結(jié)構(gòu)物作用探討[D].大連：大連理工大學(xué),2011.Liu Zanqiang.Freak wave simulation and its interaction with nuclear power station ingarage structure[D].Dalian：Dalian University of Technology,2011.

[9]黃國興.畸形波的模擬方法及基本特性研究[D].大連：大連理工大學(xué),2002.Huang Guoxing.Research on the simulation method and basic characteristics of freak waves[D].Dalian：Dalian University of Technology,2002.

[10]Kriebel D L.Efficient simulation of extreme waves in a random sea[C]//Abstract for Rogue Waves 2000 Workshops.Brest,France,2000：1-2.

[11]裴玉國.畸形波的生成及基本特性研究[D].大連：大連理工大學(xué),2007.Pei Yüguo.The generation of freak waves and its behaviors[D].Dalian：Dalian University of Technology,2007.

[12]劉贊強,張寧川,俞聿修.改進的相位調(diào)制法模擬畸形波：I-理論模型與驗證[J].水動力研究與進展,A輯,2010,25(3)：383-390.Liu Zanqiang,Zhang Ningchuan,Yü Yüxiu.The generation of freak waves based on a modified phase modulation method：I-Theory and validation[J].Chinses Journal of Hydrodynamics,Series A,2010,25(3)：383-390.

[13]Clauss G F,M Klein,Kauffeldt A.Limiting loads and motions of ships in extreme sea states[C]//13th Congress of Int1.Maritime Assoc.of Mediterranean IMAM 2009,12-15 Oct.2009.Istanbul,Turkey,2009.

[14]Bennett S S,Hudson D A,Temarel P.The influence of forward speed on ship motions in abnormal waves：Experimental measurements and numerical predictions[J].Journal of Fluids and Structures,2013,154-172.

[15]董艷秋,楊冠聲,陳學(xué)闖.關(guān)于一種危害船舶安全的波浪—畸形波的探討[J].船舶力學(xué),2003,7(2)：33-38.Dong Yanqiu,Yang Guansheng,Chen Xuechuang.Study on the great damage caused by freak wave[J].Journal of Ship Mechanics,2003,7(2)：33-38.

[16]Gunther F G,Christian E S,Marco Klein.Generation of rogue waves with predefined steepness[C]//Proceedings of OMAE 2006,June 4-9,2006.Hamburg,Germany,2006.

[17]吳 明.不規(guī)則波中艦船搖蕩運動的數(shù)值模擬及預(yù)報研究[D].大連：海軍大連艦艇學(xué)院,2013 Wu ming.The research on prediction and simulation of ship motions in irregular waves[D].Dalian：Dalian Naval Academy,2013.