孫幫碧，張志宏，劉巨斌，鄧 輝

(海軍工程大學(xué)理學(xué)院，湖北武漢 430033）

淺水航道艦船水壓場理論解及其計算

孫幫碧，張志宏，劉巨斌，鄧 輝

(海軍工程大學(xué)理學(xué)院，湖北武漢 430033）

基于淺水波動勢流理論，建立亞臨界航速和線性跨臨界航速淺水波動控制方程?；诒〈俣?，利用傅里葉積分變換法獲得淺水航道艦船水壓場理論解。通過數(shù)值求解，分析艦船航速、水深、航道寬度、色散效應(yīng)對艦船水壓場特性的影響，并與源匯分布法計算結(jié)果進行對比，驗證所建立的淺水航道艦船水壓場數(shù)學(xué)模型和計算方法的正確性與有效性。

淺水;航道;艦船;水壓場;理論解

0 引言

航行艦船引起的水底壓力變化簡稱為艦船水壓場。艦船水壓場是艦船物理場的自身特性，是設(shè)計新型水雷及研制水壓掃雷具以及艦艇自身防護的依據(jù)。Tuck［1］基于漸近匹配展開法，得到淺水薄船水壓場的近似解;Muller［2］利用有限水深興波源勢方法，在船舶縱向剖面上布置源匯，給出薄船、線性自由表面條件下的艦船水壓場計算公式;Sahin和Hyman，Nguyen［3－4］等采用有限水深格林函數(shù)方法，系統(tǒng)開展了潛艇、水面艦艇和氣墊船水壓場的理論建模和數(shù)值計算工作;Jiang［5］和 Gourlay［6－7］基于淺水波動勢流理論對艦船興波和升沉與縱傾作了分析與計算;黎昆［8］、繆濤［9］等利用面元法計算艦船在水底引起的壓力分布。

1 控制方程

假設(shè)流體為理想不可壓流體作無旋運動，艦船沿淺水航道中心線運動，船體兩側(cè)流場左右對稱。如圖1取動坐標(biāo)系固結(jié)于船體，坐標(biāo)原點位于船體水線中心，z軸垂直向上，船頭指向x軸正方向，水深為h，自由表面起伏為ζ，船速恒為V。在動坐標(biāo)系下水流相對于船體以速度V沿x軸負(fù)方向運動。

圖1 坐標(biāo)系Fig.1 Coordinate system

艦船運動引起的流場擾動速度勢Φ應(yīng)滿足下列方程:

如果不考慮式(8）中第三項 (色散項），則式(8）蛻化為亞臨界航速淺水控制方程 (SSW），即

2 理論模型

針對SSW和LTSW方程，結(jié)合船體邊界條件以及岸壁條件即可求出滿足相應(yīng)條件的艦船水壓場理論解。根據(jù)薄船假定，船體邊界條件為

式中:S(x）是船體橫截面面積分布，Sx(x）為S(x）對x求導(dǎo)。

采用數(shù)學(xué)船型，其橫截面面積分布為

2.1 采用SSW方程

2.1.1 矩形航道

1）源匯分布法

為簡便計，省略表示平均的符號 “—”。艦船以亞臨界航速在淺水中航行時，速度勢Φ應(yīng)該滿足控制方程式(9）、船體邊界條件式(10）以及岸壁條件式(6）。通過源匯分布法易于求得速度勢Φ，并根據(jù)拉格朗日積分求得艦船水壓場水底壓力系數(shù)為

根據(jù)流場對稱性，取y＞0。對船體邊界條件式以及岸壁條件式進行傅離葉變換，得

2.2 采用LTSW方程

2.2.1 矩形航道

針對無量綱控制方程式(7），省略平均符號“—”和無量綱符號“*”，采用傅離葉積分變換法，分超臨界航速和亞臨界航速2種情況進行求解。

當(dāng)k在－∞ ～∞間變化時，λ2存在大于和小于2種情況，需要分別討論:

① 當(dāng) λ2＞0時，類似超臨界航速求解方法，有

如果γ=0，不計色散效應(yīng)，則式(33）和式(34）蛻化到SSW方程開闊海域的解。

3 計算結(jié)果與分析

基于SSW方程，采用式(19）和式(20）計算淺水開闊海域亞臨界航速艦船水壓場，2種方法 (傅離葉積分變換法和源匯分布法）計算得到的艦船水壓場縱向通過特性曲線符合良好，如圖2所示。由于未計及色散性影響，故計算得到的艦船水壓場負(fù)壓區(qū)呈現(xiàn)“V”型對稱分布。船首與船尾附近出現(xiàn)正壓峰值，船中下方出現(xiàn)負(fù)壓峰值cpmin(cpmin為水底壓力變化的最小壓力系數(shù)）。在橫距y=0.1L增大到y(tǒng)=0.5L時，負(fù)壓峰值和正壓峰值均減小，但負(fù)壓區(qū)間寬度呈逐漸增大趨勢。

圖2 2種方法計算結(jié)果對比Fig.2 The comparison of two kinds ofmethods

當(dāng)水深傅汝德數(shù)由Fh=0.26逐漸增加到Fh=0.64時，逐漸增大，說明在船中下方產(chǎn)生的水壓信號逐漸明顯，如圖3所示。但在Fh≤0.3時，艦船水壓場縱向通過曲線形狀基本保持不變，說明此時艦船興波對艦船水壓場的影響很小，興波效應(yīng)可以不予考慮。而當(dāng)Fh＞0.3時，與Fh關(guān)系密切，艦船水壓場需要計及興波影響。

當(dāng)水深由0.1L逐漸增大至0.5L時，因為水底距離船體變遠(yuǎn)，因而迅速減小，說明艦船水

圖3 不同水深傅汝德數(shù)下的水壓通過曲線Fig.3 Longitudinal curves of pressure at different Fh

壓場屬于近場，艦船引起的壓力擾動在深度方向衰減較快，如圖4所示。

在矩形航道情況下，為滿足岸壁不可穿透條件，相當(dāng)于在艦船左右兩側(cè)存在一系列艦船鏡像，所有鏡像的貢獻將導(dǎo)致艦船水壓場的負(fù)壓峰值比開闊海域的結(jié)果要大。當(dāng)岸壁寬度逐漸增大時，岸壁效應(yīng)逐漸減弱，逐漸減小，如圖5所示。當(dāng)w＞5L時，艦船水壓場通過曲線與開闊海域計算結(jié)果基本一致，說明此時岸壁影響可以忽略不計。

圖4 不同水深下的水壓通過曲線Fig.4 Longitudinal curves of pressure at different depth

圖5 不同寬度下的水壓通過曲線Fig.5 Longitudinal curves of pressure at differentwidth

基于LTSW方程，考慮色散效應(yīng)，采用式(33）和式(34）計算亞臨界和超臨界航速條件下艦船水壓場的空間分布，如圖6和圖7所示。在h=0.1L及Fh=0.88的亞臨界航速情況下，類似水面興波現(xiàn)象，艦船水壓場分布呈現(xiàn)散波與橫波的特點，艦船水壓場負(fù)壓區(qū)不再呈現(xiàn)對稱的“V”型分布，而具有明顯的后移特征。當(dāng)Fh增大時，散波夾角隨之增加，負(fù)壓峰值進一步后移，甚至有可能移至船尾之后。當(dāng)Fh→1時，散波逐漸消失，在船首、船尾附近出現(xiàn)2道明顯的橫波。在h=0.1L及Fh=1.8的超臨界航速情況下，艦船水壓場的橫波逐漸消失，散波逐漸加強，散波夾角隨Fh增加而減小，且壓力系數(shù)峰值呈現(xiàn)下降趨勢。

圖6 亞臨界航速艦船水壓場Fig.6 Ship hydrodynamic pressure field at subcritical speed

圖7 超臨界航速艦船水壓場Fig.7 Ship hydrodynamic pressure field at supercritical speed

4 結(jié)語

基于勢流理論和薄船假定，在開闊海域和矩形航道中建立了基于SSW和LTSW控制方程的淺水艦船水壓場理論模型，采用傅離葉積分變換法獲取了艦船水壓場理論解，并與源匯分布法計算結(jié)果進行了比較，驗證所建立模型的正確性和有效性。同時分析了艦船航速、水深、航道寬度、色散效應(yīng)等因素對艦船水壓場的影響。結(jié)論如下:

1）在淺水亞臨界航速條件下，基于SSW控制方程，計算得到的艦船水壓場負(fù)壓區(qū)呈現(xiàn)對稱的“V”型分布，船首與船尾附近出現(xiàn)正壓峰值，船中下方出現(xiàn)負(fù)壓峰值。負(fù)壓峰值隨Fh增大而增大，但當(dāng)Fh≤0.3時，艦船水壓場與Fh基本無關(guān)，說明此時艦船興波對水壓場的影響可以忽略。水深對艦船水壓場影響明顯，艦船水壓場隨水深增加衰減較快，說明艦船水壓場屬近場。岸壁的存在相當(dāng)于艦船系列鏡像的作用，并導(dǎo)致艦船水壓場的負(fù)壓峰值增加，但當(dāng)w＞5L時，岸壁影響可以忽略不計。

2）在淺水開闊海域條件下，基于LTSW方程，考慮色散效應(yīng)，計算得到的亞臨界和超臨界航速艦船水壓場空間分布具有表面興波的特征。亞臨界航速艦船水壓場負(fù)壓區(qū)不再是對稱的“V”型分布，而是呈現(xiàn)散波與橫波的特點，并具有明顯的后移特征。在超臨界航速情況下，艦船水壓場的橫波逐漸消失，散波逐漸加強，散波夾角隨Fh增加而減小，且壓力系數(shù)峰值呈現(xiàn)下降趨勢?？紤]色散性效應(yīng)時，由于航道中艦船水壓場理論模型中的被積函數(shù)存在無限多個不等距奇點，計算較為復(fù)雜，留待以后進一步分析和解決。

［1］TUCK E O.Shallow water flows past slender bodies［J］.Fluid Mech，1966(26）:81 －95.

［2］MULLER E.Analysis of the potential flow field and of ship resistance in water of finite depth［J］.I.S.P.，1985(32）:266－277.

［3］SAHIN I，HYMAN M C.Simulation of three-dimensional finite-depth wave phenomenon formoving pressure distributions［J］.Ocean Engineering，2001，28(12）:1621 －1630.

［4］NGUYEN T C.Integral representation of linear free-surface potential flows for bottom pressure calculation［C］.Oceans 2009，MTS/IEEE Biloxi-Marine Technology for Our Future:Global and Local Challenges，2009.

［5］JIANG Tao.Ship waves in shallow Water［R］.Germany:Fortschritt-Berichte VDI，2001.

［6］GOURLAY T P.Flow beneath a ship at small underkeel clearance［J］.Journal of Ship Research，2006，50(3）:250－258.

［7］GOURLAY T P.Slender-body methods for predicting ship squat［J］.Ocean Engineering，2008，35(2）:191 － 200.

［8］黎昆，張志宏，顧建農(nóng)，等.利用面元法計算艦船在水底引起的壓力分布［J］.海軍工程大學(xué)學(xué)報，2011，23(1）:43 －46.

LIKun，ZHANG Zhi-hong，GU Jian-nong，et al.Calculating pressure distribution on water bottom caused by a moving ship with panelmothod［J］.Journal of Naval University of Engineering，2011，23(1）:43 －46.

［9］繆濤，張志宏，顧建農(nóng)，等.面元法求解有限水深船舶興波及水底壓力變化［J］.計算力學(xué)學(xué)報，2012，29(3）:464－469.

MIAO Tao，ZHANG Zhi-hong，GU Jian-nong，et al.The calculation of ship wave and bottom pressure variation in finite depth by panel method［J］.Chinese Journal of Computational Mechanics，2012，29(3）:464－469.

The theoretical solution and calculation of the hydrodynam ic pressure field of ship in the shallow channel

SUN Bang-bi，ZHANG Zhi-hong，LIU Ju-bin，DENG Hui
(College of Science，Navel University of Engineering，Wuhan 430033，China）

Based on the shallow water wave potential flow theory，at low subcritical speed and linear transcritical speed shallow water wave governing equations were established.According to slender ship assumption，The theoretical solution of hydrodynamic pressure field caused by ship moving at shallow water channelwas obtained by using Fourier integral transform method.The influences ofmoving speed of ship，the wide of channel，depth of water and dispersion effect on the ship hydrodynamic pressure field were analyzed by using numerical simulation，and the resultswere compared with the ones of sources distribution method.Themathematicalmodel and the calculation method were validated.

shallow water;channel;ship;hydrodynamic pressure field;theoretical solution

U661.1

A

1672－7649(2014）04－0011－06

10.3404/j.issn.1672－7649.2014.04.002

2013－06－13;

2013－07－12

國防預(yù)研基金資助項目(513030203），海軍工程大學(xué)自然科學(xué)基金資助項目(HGDQNEQJJ12002）

孫幫碧(1990－），男，碩士研究生，研究方向為艦船水壓場。