基于雷達測波現實方向譜的船舶三自由度搖蕩研究

王作超，石愛國，楊新棟，吳 明，劉 猛

(海軍大連艦艇學院，遼寧 大連116018)

0 引 言

自20世紀60年代，Oudshoom［1］等開始利用雷達獲取海浪信息以來，經過科研人員的不懈努力，這一技術已趨向成熟，也有不少商業(yè)化的產品面世。其中，德國GKSS 開發(fā)的WaMoSⅡ［2］以其產品的穩(wěn)定性、測量的準確性通過了德國船級社(GL)和挪威船級社(DNV)的認證最具代表性，在世界范圍內得到了廣泛應用。我國自2000年，在“遠望”系列測量船首次引進WaMoSⅡ以來，也已陸續(xù)在各型測量艦艇安裝了多套WaMoSⅡ設備。

雷達測波典型的系統(tǒng)組成為:X 波段航海雷達、A/D 轉換器、PC和相應的軟件包。其工作的基本原理是由于海水表面的毛細波與X 波段的雷達波波長相近，產生Bragg 共振現象，使得后向反射的雷達回波得到強化［3］。通過分析雷達回波的時空變化規(guī)律，獲得海浪譜與表面流信息，給出有義波高、主波波向、主波周期、主波周長、流向、流速等信息［4］。與測波浮標、SAR 等其他測波設備相比，WaMoSⅡ最大的特點在于測量的實時性、并可用于船舶搭載，簡便易行。若能將WaMoSⅡ測量的實時海浪信息，結合CFD 技術，解算出此海況下真實的船舶搖蕩信息，這對于提高船只的海上作業(yè)安全、大風浪中的航行策略優(yōu)選、艦載機的起降都大有裨益。

隨著微機電(MEMS)技術的持續(xù)發(fā)展，以其為基礎的運動姿態(tài)測量系統(tǒng)也不斷面世。此類系統(tǒng)一般內置三軸陀螺儀、加速度計，通過其內嵌的微處理器，可輸出校準過的姿態(tài)測量數據。相較其他類型的姿態(tài)測量系統(tǒng)，此類產品具有體積小、低功耗、性能可靠、抗干擾、易于二次開發(fā)等特點。2012年6月在“中國漁政206 號船”(以下簡稱“206 船”)搭載了某型基于微機電的姿態(tài)測量設備(見圖1)，記錄了船舶首搖、縱搖、橫搖及3 個方向加速度時歷。利用“206 船”所載WaMoSⅡ所給出的方向譜數據，基于RNG k-ε 模型［5］，建立了三維數值波浪水池，在水池入口邊界采用邊界造波法生成波浪［6］，在出口處利用阻尼消波法實現消波［7］，模擬了“206 船”的三自由度搖蕩運動，通過與測量數比較，二者吻合較好。

圖1 “206 船”搭載的膜電位測姿儀Fig.1 MEMS motion measurement system of target ship

1 方向譜與短峰不規(guī)則波

1.1 方向譜

海浪的方向譜反映了波浪能量隨頻率和方向的分布情況，實際研究中一般將方向譜看作是統(tǒng)計頻譜和方向擴散函數的乘積［8］

S(ω，θ)=S(ω)D(ω，θ)。

但是統(tǒng)計頻譜是通過長期觀測獲得的海浪場的一般規(guī)律，反映不了實際海浪的特殊性。圖2 給出了WaMoSⅡ記錄的頻譜和對應參數下的ITTC 譜，可以看出，就此海區(qū)而言，WaMoSⅡ實測頻譜的能量更高，也更為集中。本文通過WaMoSⅡ給出的方向譜數據，直接提取方向、頻譜和對應的譜密度值，用于數值水池的構建，可以模擬更為真實的海浪環(huán)境。

圖2 WaMoSⅡⅡ頻譜和ITTC 譜頻對比圖Fig.2 The comparison of spectrum between WaMoSⅡⅡand ITTC

有以下2 點值得注意:1)WaMoSⅡ給出的頻譜和方向譜的譜密度值均為相對值，要結合有義波高將其轉化為絕對值;2)同導航雷達相同，雷達測波系統(tǒng)也會受到降水、雨雪的干擾，影響測量的準確性。為此，WaMoSⅡ在系統(tǒng)中引入了數據質量因子IQ(Quality Index)，只有在IQ 為000 時，才能認為其測量的數據可靠。

1.2 短峰不規(guī)則波

一般可將短峰不規(guī)則波視作無限多頻譜、多方向、振幅變化且相位隨機的微幅波的疊加:

式中:aij為第i 個頻率、第j 個方向成分波的波幅;ωi和ki分別為第i 個頻率成分波的圓頻率和波數;θj是第j 個方向成分波的方向;εij為第i 個頻率、第j 個方向的成分波的隨機相位角，是［0，2π］區(qū)間內均勻分布的隨機數。

根據上文提到的WaMoSⅡ有關數據質量要求，選取了“206 船”于2012年6月13日1000 時刻采集的方向譜信息。其有義波高為2.0 m，平均周期7.34 s，主波向67°，IQ 為000，角度范圍30°-110°，間隔4°，子波頻率0.054 7 ～0.202 3 Hz，間隔0.005 5 Hz。

2 數值模擬

2.1 模擬方案

為了提高數值模擬的計算效率，將原型船按照1:40 尺度進行了縮小。表1 給出了原型船及船模相應尺度的主要參數。根據原型船型線圖繪制的船模如圖3所示。

圖3 目標船船模Fig.3 The model of target ship

2.2 計算域

計算域的設置如圖4所示，入口距船模船首1L，出口距船尾2L，頂部邊界距水線0.7L，底部邊界距水線1L，左、右邊界距船中縱剖面1L。

圖4 數值模擬中的計算域劃分(俯視圖)Fig.4 The calculation domain of simulation(top view)

表1 船模和原型船的主要參數Tab.1 The main parameters of real ship and moel

2.3 網格整體移動

在傳統(tǒng)船舶CFD 計算中，一般采用動網格技術，即船體周圍的網格部分發(fā)生形變，在損失一定網格質量的情況下，完成搖蕩運動的計算。可是采用動網格方法的缺陷也較為明顯:如果船模發(fā)生較大幅度的搖蕩運動，會使得網格質量下降到無法計算，最終使得計算發(fā)散。

本文創(chuàng)新性地提出了網格整體移動方法，求解船舶的三自由搖蕩運動。該方法的基本思路是:在UDF 的編寫中，根據船體所受的力和力矩計算船體的線速度和角速度，使得計算域內的全體網格，隨船體一起運動，這樣就保持了計算域內網格的空間拓撲結構不變，在計算過程中，有效地避免了因網格質量下降導致的數值計算發(fā)散問題。

2.4 邊界條件

邊界條件的具體設置如下:

1)造波區(qū)入口邊界——速度入口，給定波浪沿3 個方向的分速度，根據波高給出水的體積分數及k，ε 值;

2)消波區(qū)入口邊界——速度入口，給定流速沿3 個方向的分速度，根據靜水面位置給出水的體積分數及k，ε 值;

3)上下邊界——速度入口，給定3 個方向流速(u=Uc，v=w=0)及水的體積分數(上邊界為0，下邊界為1);

4)船體——有剪切力無滑移的壁面。

3 數值模擬結果及分析

3.1 數值造波的準確性

為了驗證數值造波的準確性，利用5 點波高儀監(jiān)測數值水池的波面時歷，并利用BMD(貝葉斯方法)法對時歷進行譜分析，結果表明模擬譜與目標譜吻合較好，所采用的數值造波方法可以用作船舶耐波性計算的模擬波浪環(huán)境。

3.2 模擬結果

利用方向譜生成波浪場，待流場發(fā)展較為充分后，加載船體的垂蕩、橫搖和縱搖三自由度運動方程，再待船模運動穩(wěn)定后的某一時刻作為起點，記錄船模的三自由度搖蕩時歷。

由于數值模擬結果與測姿設備所獲數據的初始相位不同，不能進行時歷比較，只能通過對時歷進行譜分析，對數值模擬結果評價。圖7 ～圖9 分別給出了垂蕩、橫搖和縱搖的譜密度值。

圖5 目標船實測橫搖運動時歷Fig.5 The time history of real roll of target ship

圖6 目標船模擬橫搖運動時歷Fig.6 The time history of simulated roll of target ship

圖7 數值模擬與實測數據的橫搖譜密度比較Fig.7 The comparison of roll between simulation and real measurement data

圖8 目標船實測縱搖運動時歷Fig.8 The time history of real pitch of target ship

圖9 目標船模擬縱搖運動時歷Fig.9 The time history of simulated pitch of target ship

圖10 數值模擬與實測數據的縱搖譜密度比較Fig.10 The comparison of pitch between simulation and real measurement data

圖11 目標船實測垂蕩運動時歷Fig.11 The time history of real heave of target ship

圖12 目標船模擬垂蕩運動時歷Fig.12 The time history of simulated heave of target ship

圖13 數值模擬與實測數據的垂蕩譜密度比較Fig.13 The comparison of heave between simulation and real measurement data

從圖7、圖10、圖13 可以看出，本文的數值模擬結果與測姿設備得到的結果是基本吻合的，基本可以體現出運動趨勢，但是譜值中心發(fā)生了少許偏移，原因為:1)所采用的浪級過大，所采數據的實際波高為2 m，船舶在此條件下的運動運動，具有一定的非線性特征;2)數值粘性對于低頻運動影響比較明顯。這也是導致譜誤差較大的原因之一。

4 結 語

本文基于雷達測波當前方向譜，構建了三維數值波浪水池，對船舶的三自由度搖蕩運動進行了模擬，并與實船測姿結果進行了比較。試驗表明:本文給出的船舶三自由搖蕩模擬方法能夠較好地模擬船舶在實際海況中的運動情況。可以預期的是，隨著計算機硬件的不斷發(fā)展，數值模擬算法效率的不斷提升，通過雷達測波系統(tǒng)獲得海上實際情況，再利用CFD 技術進行船舶搖蕩運動預報，對于未來船只的海上作業(yè)必將具有重要且現實的意義。

［1］OUDSHOORN H M.The use of radar in hydro-dynamic surveying［J］.Coastal Engineering，1960(7):59-76.

［2］Reference List［EB/OL］.http://www.oceanwaves.org/start.html.

［3］VALENZUELA G R.Theories for the interaction of electromagnetic and oceanic waves a review［J］.Boundary-Layer Meteorology，1978(19):87-105.

［4］YOUNG I R，ROSENTHAL W，ZIEMER F.A Threedimensional analysis of marine radar images for the determination of ocean wave directionality and surface currents［J］.Journal of Geophysical Research，1985，90:1049-1059.

［5］STEVEN O A，et al.Renormalization group modeling and turbulence stimulations［M］.Near Wall Turbulence Flows.1993:1031-1046.

［6］劉霞，譚國煥，王大國.基于邊界造波法的二階Stokes 波的數值生成［J］.遼寧工程技術大學學報(自然科學版)，2010，29(1):107-111.

LIU Xia，TAN Guo-huan，WANG Da-guo.Numerical simulation of second-order Stokes based on wave-generation method of defining inlet boundary conditions［J］.Journal of Liaoning Technical University(Natural Science)，2010，29(1):107-111.

［7］李凌，林兆偉，尤云祥，等.基于動量源方法的黏性流數值波浪水槽［J］.水動力學研究與進展(A 輯)，2007，22(1):76-82.

LI Ling，LIN Zhao-wei，YOU Yun-xiang，et al.The numerical wave flume of the viscous fluid based on the momentum source method［J］.Journal of Hydrodynamics(Ser.A)，2007，22(1):76-82.

［8］PARK J C，et al.Numerical reproduction of fully nonlinear multi-directional waves by a viscous 3D numerical wave tank［J］.Ocean Engineering，2004，31(11):1549-1565.

［9］吳乘勝，朱德祥，顧民.數值波浪水池及頂浪中船舶水動力計算［J］.船舶力學，2008，12(2):171-179.

WU Cheng-sheng，ZHU De-xiang，GU Min.Computation of hydrodynamic forces for a ship in regular heading waves by a viscous numerical wave tank［J］.Journal of Ship Mechanics，2008，12(2):171-179.