黃 華，曾昭鑾，詹杰民，付 靜，李 奇

(中山大學工學院應用力學與工程系，廣東 廣州 510275)

海床上的建筑物一般會受到孔隙海床中由海洋波浪引起的滲流壓力作用。波浪滲流壓力可對其產生波浪浮托力和傾覆力矩，進而對結構實際承受的波浪載荷產生明顯影響。因此，研究波浪滲流作用具有重要的工程應用上的意義[1]。在假設土介質不變形和孔隙水不可壓，滲流為達西流的條件下，Durand和Monkmeyer[2]對無限厚度海床的情況求得滲流壓力分布的解析解。邱大洪和莊嚴[3]給出了重力式平臺底部所受浮托力及力矩的數值解。孫昭晨和邱大洪[4]基于Biot固結理論,考慮了土介質彈性變形和孔隙流體的壓縮性，給出了圓墩柱下有限深度海床波浪滲流壓力解析解。鄒志利和邱大洪等[5-9]考慮了墩柱直接放置在海床上和放置在埋入海床中拋石基床上兩種情況，給出了安置在海床上大尺度墩樁底面波浪滲流壓力的計算方法和一些數值結果,以及計算了在線性波或橢圓余弦波作用下雙國墩柱波浪滲流力，討論了雙墩柱沿波向串列和并列布置時墩柱上滲流力隨墩距的變化。另外，實際海洋在一定條件下會呈現層化現象，而海洋的層化性將導致海水內波的生成，內波的產生將導致波浪作用機制的改變，對海上建筑物的作用不容忽視。吳建華等[10]計算了層化海中浮筒水波輻射的附加質量和阻尼系數。YEUNG 等[11]對兩層海結構的輻射和繞射問題進行了探討。尤云祥等[12]對兩層海中大直徑圓柱的波浪力進行了計算分析。付靜等[13-14]研究了兩層海中透空圓柱(群)的水動力學特性問題。

本文將Biot固結理論與層化海繞射波理論相結合來研究表面波和內波、結構及海床間相互作用，因而是對均勻海波浪滲流力問題研究的推廣。通過特征函數展開法對二層海中固立直柱的海底滲流壓力場進行求解。在采用分區(qū)法時，本文按海水上下層流域及滲流壓力場內外域，分別使用多種特征級數解，給出了表面波和內波分別所致的滲流壓力數學解，并進而推導了直柱底部所受波浪滲流浮托力和傾覆力矩的解析算式。計算結果反應了結構底部的波浪浮托力和傾覆力矩隨海水上下層密度比等各種參數的變化趨勢，說明內波和表面波所致波浪滲流力與波浪對圓柱的直接水平作用在一定條件下具有相同量級，因而在實際海工結構設計中應對此加以一定重視。

1 二層海中圓柱繞射波場與壓力場的求解

如圖1所示，設上、下層流體的密度分別為ρ1和ρ2(ρ2>ρ1)，厚度分別為h1和h2，水深為d，海床厚度為h。另采用Biot固結理論的相關假設，即假定立于海床上的圓柱結構與海床間無相對滑動，海床底部不透水，土體各向同性、可滲且具有線彈性性質。由于波浪在海底處的水質點速度一般遠大于海床內的滲流速度，故在繞射波浪場的求解中可設海底不可滲透，即海域的波浪場不受海床滲透影響，因而波浪場與滲流場可分別求解。對于波浪場，對應的Stokes一階繞射波勢的邊值問題可表達為：

圖1 二層海中可滲透海底上的直立圓柱

(1)

(2)

(3)

(4)

(5)

(6)

式中，T1=tanhkh1，T2=tanhkh2，λ=ρ1/ρ2。由線性繞射波的相關理論，利用圓柱徑向表面條件可得：

(7)

Z(z)=Z±(z)=

(8)

圓柱側面流場動壓強由伯努利方程可得：

P±=Re(p±),p±=ρgA±e-iω± tZ±(z)·

(9)

(10)

(11)

2 滲流波浪壓力的計算

(12)

(13)

(14)

(15)

(z=-h,r>a)

(16)

(17)

(18)

(19)

(20)

p±(1)=ρ2gA±e-iω± tE±(k)·

(21)

(22)

由式(21)和(22)有:

p±(1)|z = 0=ρ2gA±e-iω± tE±(k)·

(23)

p±(1)=ρ2gA±e-iω± tE±(k)·

(24)

(A±ρ2ga2π)E±(k)e-iω±t·

(25)

(A±ρ2ga3π)e-iω±tE±(k)·

(26)

3 算例與分析

為方便計，實算中設海床不變形及孔隙水不壓縮，即可取Cs=0，另設A+=A-=A。對水平波浪力與垂直浮托力，本文采用ρ1gAπa2進行無量綱化；對水平波浪力矩和傾覆力矩，本文用ρ1gπAa3進行無量綱化；ka為無量綱波數。

圖2表明，在一定的ka范圍，表面波水平波浪力較內波波浪力大。表面波波浪力隨ka的增大而增大，達到極值后減少；內波波浪力則呈單調增的趨勢。在一個ka取值較小范圍內，表面波水平力矩與內波力矩大致相等；隨著ka的增大，兩者都呈現出先增后減再增的變化趨勢，且表面波水平力矩較內波力矩略小，反映了在層化海洋中，內波對結構的直接傾覆作用可能超過表面波作用。

圖3表明，在相同的海況條件下，當ka較小時，表面波水平波浪力較浮托力小,而ka較大時，情形相反。而在一定的ka取值范圍內，內波水平波浪力明顯比內波所致浮托力要小，表明內波引起的浮托力作用不容忽視。

圖4表明，在一定海況條件下，當ka較小時，表面波引起的浮托力及傾覆力矩均超過了內波引起的波浪滲流作用，而隨著ka的增大，情況相反，反映了由內波所致對結構的波浪滲流作用與波浪直接作用一樣可能超過表面波對應的作用，故對內波影響要加以足夠重視。

圖2 h/a=3/4,h1/h2=2/3,λ=0.9，表面波與內波的水平波浪力(矩)

圖3 h/a=3/4,h1/h2=2/3,λ=0.9，表面波及內波波浪力與對應浮托力的比較

圖5表明，表面波、內波引起的水平力矩均較相應的傾覆力矩大，且隨著ka的增大，兩者差距增大。但在ka較小時，兩者相差不大，量級基本相同，因此在結構的實用當中對滲流傾覆力矩應給予充分重視。當ka較大時，表面波及內波引起的傾覆力矩趨弱。

圖6-7表明，在一定海況條件下，隨著上下層流體密度的變化，表面波及內波所致浮托力變化不大，但由其引起的傾覆力矩變化較為明顯；密度比越小，浮托力及傾覆力矩越大，反映出海洋層化性對滲流波浪力具有一定的影響。

4 結 論

本文將兩層海的波浪理論推廣應用于波浪滲流作用場的理論研究，推導了由表面波與內波分別所致可滲透海床內滲流壓力場的解析算式，進而對滲流壓力所致對圓柱底部的浮托力和傾覆力矩進行了計算。結果表明，在一定條件下，內波所致浮托力可能超過內波水平波浪力，而內波波浪滲流作用與其波浪直接作用類似均可能超過表面波作用，因而對內波的影響不容忽視。另一方面，波浪引起的水平力矩和滲流傾覆力矩量級可能接近，因此對用于軟基海底的固立結構，在實際設計和使用中應對滲流傾覆力矩給予關注。此外，上下層流體密度比趨小將可能增強滲流傾覆力矩的作用。

圖4 h/a=3/4,h1/h2=2/3,λ=0.9，表面波與內波所致浮托力及傾覆力矩的比較

圖5 h/a=3/4,h1/h2=2/3,λ=0.9, 水平波浪力矩與傾覆力矩的比較

圖6 d=h=a,h1/h2=2/3，不同λ值下的表面波及內波所致浮托力

圖7 d=h=a,h1/h2=2/3，不同λ值下的表面波及內波所致傾覆力矩

參考文獻：

[1]邱大洪. 波浪滲流力學[M].北京：國防工業(yè)出版社,2006：193-203.

[2]DURAND T J P，MONKMEYER P L. Wave-induced seepage effects on a vertical cylinder [C]. Proceeding of Coastal Engineering，1980:1781-1798.

[3]莊嚴，邱大洪.海上重力式建筑物基底的滲壓荷載[J].海洋學報，1987,9(2)：239-254.

[4]孫昭晨,邱大洪.作用于可滲可壓縮海床上的墩柱底面上的波浪力[J].海洋學報,1989,11(3):364-371.

[5]鄒志利，邱大洪.無限深可壓縮海床上圓柱墩波浪滲流壓力解析解[J].海洋學報,1993,15(3):106-115.

[6]鄒志利，邱大洪.可滲可壓縮海床上墩柱底面波浪滲流壓力的計算[J].海洋工程,1993,11(4):47-58.

[7]邱大洪,鄒志利.雙墩柱波浪滲流力[J].大連理工大學學報,1993,33(2):232-237.

[8]邱大洪,陳健.拋石基床上圓柱墩底部的波浪浮托力[J].海洋學報,1994,16(1):114-123.

[9]邱大洪,臧軍. 群墩上的波浪浮托力[J]. 中國港灣建設,1995(5):8-12.

[10]吳建華, 吳秀恒, 李世漠. 大尺度物體在層化海洋中的波浪繞射和散射理論(Ⅳ)本征函數展開法[J]. 水動力學研究與進展, 1991, 6(3): 46-53.

[11]YEUNG R W, NGUYEN T. Radiation and Diffraction of Waves in a Two-layer Fluid[C]∥Proceedings of the 22th Symposium on Naval Hydrodynamics. Washington, D. C: National Academy Press,1999: 875-887.

[12]尤云祥, 朱偉, 繆國平. 分層海洋中大直徑樁柱的波浪力[J].上海交通大學學報, 2003, 37(8): 1181-1185

[13]付靜，黃華，朱慶勇. 二層流體中透空大直徑樁柱的波浪特性分析[J]. 中山大學學報：自然科學版，2008，47(3):14-18.

[14]付靜,黃華,詹杰民,等. 二層海洋中波浪對透空圓柱群的繞射作用分析[J].中山大學學報:自然科學版,2008,47(2):76-80.