，

(武漢理工大學 交通學院,武漢 430063)

1 淺水條件下船舶興波計算

船舶興波阻力計算采用的面元法可分為兩類，按照所采用的格林函數(shù)不同，分別為Rankine源方法和Kelvin 源方法(又稱Havelock源方法)。

采用Rankine源格林函數(shù)，在整個流場邊界布置源匯。在船體表面及部分自由面進行源匯分布，設船體表面S上的源強為σs(q)(這里是匯)，其所誘導的速度勢為

(1)

設自由面F上的源強為σf(q)，其所誘導的速度勢為

(2)

則總得速度勢為

φ(p)=φs(p)+φf(p)

(3)

為求解速度勢式(3)，必須求解邊界面上的源強σ。一階Rankine面元法的思想是采用一定數(shù)量的四邊形或三角形小平面(即面元)代替連續(xù)的船體曲面及自由面(以及其它邊界面)。假定每片面元上的源強為常數(shù)σi(q)，則可得到離散形式的速度勢。

將船體曲面S離散為Ns片面元，記第i片面元上Si的源強為σs(i)，則式(1)轉(zhuǎn)化為

(4)

同理，將自由面F共離散為Nf片面元，則式(2)轉(zhuǎn)化為

(5)

將離散后的總速度勢表達式分別代入各個定解條件，并在每片面元上選取一個滿足當?shù)丶s束方程的場點P(稱為配置點)，可得到總共Ns+Nf個線性方程。求解方程組，求得每個配置點所在面元的源強，從而解決問題。

Wigley船船型要素：船長L=8 m，型寬B=0.8 m，吃水T=0.5 m，Wigley單體船在深水中航行時不同F(xiàn)r數(shù)下的興波阻力系數(shù)曲線見圖1。與實驗結(jié)果比較，結(jié)果令人滿意。

圖1 Wigley單體船興波阻力系數(shù)

用面元法計算Wigley船型雙體船和三體船淺水的近域波形。計算的船型：Wigley雙體船，兩片體中縱剖面之間的距離分別為0.2L，0.3L，0.4L和0.5L。對于三體船，不同的側(cè)體排水量，不同的橫向和縱向間距，對其阻力都有極大的影響。本文Wigley三體船，以Wigley單體船為主體，側(cè)體排水量取主體的0.125倍，橫向間距分別為Ys=0.2L，Ys=0.3L,Ys=0.4L，Ys=0.5L，縱向間距Xs=0，得到的波高數(shù)據(jù)見圖2、3。

圖2 Wigley雙體船不同間距下船側(cè)0.3L處波高

圖3 Wigley三體船不同橫向間距下船側(cè)0.3L處波高

由以上波高圖形可以看出，對于Wigley雙體船和三體船，與在深水情況下相同，由于片體之間波系存在相互干擾，當片體與片體或片體與主體之間的間距越小時，其造成的興波越大，在干擾發(fā)生的船艉附近尤為明顯，見表1。

表1 淺水波高計算結(jié)果

2 船行波對直立建筑物的作用

《海港水文規(guī)范》(JTJ213-98)[1]第八章有如下規(guī)定。

當d1≥0.6H時，圖4所示的直墻建筑物上近破波的波浪力按下列規(guī)定確定。

d-建筑物前水深，m；d1-基床上水深，m；Z-靜水面以上高度， m；ps-靜水面處的波浪壓力，kPa；pb-墻底處的波浪壓力，kPa。

圖4近破波的波壓力分布

1)靜水面以上高度Z(m)處的波浪力為零，Z按下式計算。

2)靜水面處的波浪壓力按下列公式計算。

3)墻底處的波浪壓力按下式計算。

pb=0.6ps

4)單位長度墻身上的總波浪力按下列公式計算。

3 計算分析

至今對直墻建筑的破波波浪力還不能以嚴密的理論進行分析計算，在工程上只能采用經(jīng)驗方法，而經(jīng)驗公式中的波浪要素均為粗略的估算，本文則利用船舶興波的數(shù)值計算，以數(shù)值計算得到的較為準確的波形為輸入波，利用經(jīng)驗公式計算波浪力，將船行波的數(shù)值計算與波浪對岸壁的作用相結(jié)合起來計算波浪力。

本文對于近破波的計算，基床底水深d取1.28 m，堤腳水深d1取0.48 m，將在淺水條件下Wigley各種船體船側(cè)0.3L處的近域波高數(shù)據(jù)代入《海港水文規(guī)范》(JTJ213-98)公式計算得到船舶興波對直立建筑的波浪力數(shù)據(jù)，見表2。

表2 波浪力計算結(jié)果

由表2得到的壓力見圖5、6，縱坐標為直墻建筑物各點位置，橫坐標為該點處的波壓力。

可以看出，當間距較小時，雙體船和三體船船波對直立建筑物的近破波波浪力較大；隨著間距的增大，雙體船造成的波浪力逐漸減小?？梢姛o論是對雙體船還是三體船，間距的影響還是很明顯的。由于缺少有關實驗資料，無法比較用船波數(shù)據(jù)計算得到的值與實際值相差多少，這種方法計算的波浪力的準確性尚需進一步檢驗與討論。

圖5 雙體船近破波的波壓力

圖6 三體船近破波的波壓力

4 結(jié)論

采用一階Rankine源面元法對淺水Wigley雙體船和三體船的近域波形進行數(shù)值計算，并以其結(jié)果為輸入波計算了船行波對直立建筑物的近破波波浪力，結(jié)果表明該方法具有一定的可行性[2-4]，其準確性還有待進一步驗證。

[1] 交通部.JTJ213-98海港水文規(guī)范[S].北京：人民交通出版社,1998.

[2] 張 偉.高速排水型船舶興波波形及尾浪的數(shù)值計算[D].武漢：武漢理工大學，2009.

[3] 屈曉婷.破碎波對直立建筑物作用機理的實驗研究 [D].大連：大連理工大學，2008.

[4] 高 高.船舶興波問題面元法計算中自由面離散的穩(wěn)定性分析[J].水動力學研究與進展：A輯,2002(1)：53-60.