李 奇，黄 华，詹杰民，张耀中，郭 霖，许潇楠

(中山大学应用力学与工程系，广东广州 510275)

有关研究表明，在海洋大尺度离岸结构的设计与建造中，对部分结构表面加以布孔可以有效改变绕射波浪的作用机制，结构的可渗透性可以减少波浪对结构散射现象的复杂性，进而导致对波浪的一定的吸收作用，从而降低波浪对结构的总体载荷。另一方面，对于可渗透海底，除了考虑绕射波浪的直接作用外，一般还需考虑海底下孔隙海床中由波浪引起的对固立海工结构底部的渗流作用力［1］。目前，关于结构透空导致波浪减弱现象的研究以及对波浪渗流作用的研究主要限于对Airy微幅波作用问题的分析［2－4］，而透空结构的非线性波作用问题仅在理论上给予了初步探讨，缺乏实算结果说明［5］。考虑部分透空结构主要应用于近岸工程，因而对浅水波作用下透空结构的绕射波作用问题以及波浪渗流作用问题展开研究具有重要的理论和实际意义。

在浅水区，随着波浪特征要素中波高与水深比以及波长与水深比的适度增大，波浪运动的非线性影响将明显增加，继续使用Airy微幅波理论对波浪作用进行预测将产生一定误差，而应用浅水波理论进行波浪作用分析将更趋合理。在一定条件下，椭圆余弦波将成为浅水波理论计算的主要波型之一。Isaacson和苏铭德等［6－7］分别应用椭圆余弦波一阶绕射理论对单一直柱和双柱进行了波浪力计算与分析。邱大洪利用实验结果验证了椭圆余弦波一阶绕射理论用于双柱问题上的准确性［8］。另一方面，Williams［2］及付静、黄华等［9－10］应用 Airy 微幅波理论，分别研究了均匀海和二层海中固立透空柱(群)的水波绕射问题。此外，邱大洪、孙昭晨和黄华等多人应用Biot固结理论和Airy波理论，分别对均匀海与二层海中固立的密实圆柱的波浪渗流作用问题进行了相关分析［11－13］。

本文通过对椭圆余弦波的一阶绕射理论和波浪渗流固结理论加以结合与推广应用，分析了结构表面透空系数、圆环内外径比及水波非线性影响参数等一些新的变化因素对透空圆环墩柱所受波浪作用及波浪渗流作用的可能影响，揭示了浅水波作用下透空结构绕射波浪力及渗流力的各种变化规律。

1 椭圆余弦波对圆环墩柱绕射问题的一阶解

对于固立于浅水区可渗透海床上的大尺度透空圆环墩柱 (如图1所示，外侧表面透空、内柱表面密实)，基于Biot渗流固结理论，可假定墩柱底部与海床间无相对滑动，海床底部不可渗透，海床内土体各向同性、可渗且具有线弹性性质。考虑到波浪在海底处的水质点速度一般远大于海床内的渗流速度，故在浅水波绕射问题分析中可认为波浪场不受海床渗透的影响，从而波浪压力场可按海底不渗透条件而独立求解，进而再考虑波浪压力场对海床的影响而对渗流压力场进行求解［1］。根据浅水波绕射理论，在一定条件下 (例如kd＜0.6，H/d＞0.4 ，(H/d)/ka＜3)［6－7］，更适合应用椭圆余弦波绕射理论。如图1所示，将流域划分为圆环柱内区和圆柱柱外区，设Φ表示椭圆余弦波绕射波总势，对应的定解问题提法为:

图1 直立于浅水区可渗透海床上的透空圆环墩柱Fig.1 Porous vertical circular cylinder resting on porous seabed in shallow water

式中，φi与φs分别为入射与散射波势;φ(j)为对应流域vj(j=2，3)的椭圆余弦波一阶总波势;为对应外流域v2的散射波一阶分量第n个谐波的波势，它满足Sommerfeld辐射条件式 (4)，而kn为对应谐波的波数。式 (5)为椭圆余弦波一阶分量势在圆柱透空外侧表面和内柱密实表面应满足条件，其中ρ、γ及μ分别为海水密度、具有长度量纲的材料系数及流体粘性系数。由椭圆余弦波一阶理论［6］，入射波一阶势解为:

式中，q=K(κ)(kx－ωt)/π;γ=E(κ)/K(κ);K(κ)和E(κ)分别为第一、二类完全椭圆积分;κ′2=1－κ2;cnq为雅克比椭圆函数。另模量κ由方程 κ2K2(κ)=3π2H/4k2d3确定。对上式积分，并将涉及cnq的积分项进行傅里叶级数展开，可得:

进一步应用特征函数法并结合物面条件，可解得椭圆余弦波一阶总绕射波势为:

其中，Jm(x)和(x)分别为m阶第一类Bessel函数和m阶第一类Hankel函数。另有

进一步可给出圆环柱径向表面所受总水平波浪力Fh和力矩Mh分别为:

2 波浪渗流压力解及相应的渗流力解

如图1所示，对海床内渗流场分区，设区域v0中波浪所引起的渗流压力为P(0)=Re(p(0))，区域v1中为P(1)=Re(p(1))，根据Biot的渗流压力场固结理论［1］，对应P(j)的边值问题可表达为:

现求解p(0)和，推广应用特征函数法，可取:

式中，Im(x)为m阶第一类变形Bessel函数，Anml与Bnml为可定系数，需由v0与v1交界面匹配条件并借助傅氏级数展开法加以确定，该条件为:

另外，三个相关参数分别为:

进一步设圆柱底面上所受渗流波浪浮托力和倾覆力矩分别为Fv=Re(fv)及Mv=Re(mv)，依据海工结构最大荷载安全设计原则，可取Cs=0［11－13］，即可设海床不变形及孔隙水不压缩，相应可得:

3 算例与分析

实算中，对所算的最大水平波浪力和力矩分别按因子ρgHad和ρgHad2进行无量纲化。其中引入水波非线性影响参数 λ=(4κ2/3π2)K2(κ)=H/k2d3，它能更有效反映波浪特征与水深的相对变化导致水波非线性特征效应变化而产生的对波浪力的可能影响。此外对所算的最大波浪渗流浮托力和倾覆力矩分别按因子ρgπ(H/2)a2和ρgπ(H/2)a3进行无量纲化。在与水平波浪力和力矩比较时，须将无量纲因子分别转化为ρgHad和ρgHad2，为此可设α=(π/2)(a/d)，它表示浮托力与波浪力两种无量纲力的比例因子，而渗流倾覆力矩与波浪水平力矩两种无量纲力矩的比例因子相应为(2/π)α2。

图2表明，对于浅水长波或半径相对较小的墩柱，ka取值一般相对较小，此时最大无量纲水平波浪力和力矩 (最大无量纲水平力矩为力的0.5倍)随着非线性影响参数λ的增大而增大。当ka取值较大时，情形相反。图3表明，椭圆余弦波一阶绕射理论对最大波浪力的估值明显高于浅水条件下Airy波理论 (λ→0)所算结果，说明在一定浅水条件下，椭圆余弦波绕射理论对透空结构所受波浪作用的估算将更趋合理和安全。

图4表明，与密实圆柱 (G=0)相比，圆环柱表面透空对波浪作用有明显降低效应，但当透空系数增大到一定值时，波浪力的变化较小，说明结构表面应当适度透空。图5表明，随着圆环墩柱内外柱半径比b/a的减小 (透空内域增大)，墩柱的最大水平波浪力 (矩)随之减小，对于透空圆柱(b/a=0)，波浪力减至最小，不过，随着半径比的不断减小，波浪力减小效应明显趋弱，说明半径比可以进行适度取值，以达到墩柱减载和相对稳定的双重效果。

图6和图7表明，对应不同的λ取值，浮托力随ka的增大而单调减小，而倾覆力矩对应的变化趋势为先增后减，说明波浪渗流力与力矩的关系不同于波浪力与力矩的关系，后者在变化趋势保持一致。此外，浮托力随参数λ变化不大，而直接影响墩柱稳定性的倾覆力矩随参数λ的变化较为明显，且变化趋势与水平力的变化趋势 (即水平力矩变化趋势)类似。

图8和图9表明，在一定ka取值范围内圆柱表面透空对浮托力和倾覆力矩均有一定的减弱作用，其中对直接影响墩柱稳定性的倾覆力矩的减弱效应更趋明显，说明立于浅水区可渗透海底上的海工墩柱受浅水波绕射作用时，结构表面透空对波浪渗流力的影响是趋于良性的。

图8 最大无量纲浮托力 (b/a=0，h/a=1，λ=3)Fig.8 Dimensionless seepage lift force

图9 最大无量纲倾覆力矩 (b/a=0，h/a=1，λ=3)Fig.9 Dimensionless seepage overturn moment

图10和图11表明，随半径比增大 (透空内域缩小)浮托力和倾覆力矩在ka的一定取值范围内呈增大趋势，这与半径比对绕射波浪力 (矩)的影响效应相似，只是变化的幅度相对较小一些。结合半径比变化对绕射波浪力 (矩)的影响结果，再次说明通过对墩柱半径比加以调控可达到使墩柱透空减载的同时又确保安全稳定的工程要求。

图12与图13表明，在一定参数条件下，对墩柱稳定性有直接影响的最大渗流倾覆力矩值与波浪力矩值较为接近，说明近岸工程中浅水波对结构的直接作用与波浪引起的渗流作用可能具有相同量级，因此在浅水区海工墩柱的设计上对浅水波诱导的渗流力影响应给予特别重视。

图14和图15表明，尽管由椭圆余弦波理论所算最大无量纲浮托力与浅水条件下Airy波理论所算结果相差不大，但椭圆余弦波理论所算最大无量纲倾覆力矩明显高于Airy波理论的对应结果。考虑到影响海工结构稳定性的主要因素为倾覆力矩，因此，从渗流作用的角度同样可以看到，在一定条件下应用椭圆余弦波理论相对更为合理和安全。

图14 浮托力比较 (b/a=0，G=1，h/a=1，λ=3)Fig.14 Comparison of seepage uplift forces

图15 倾覆力矩比较 (b/a=0，G=1，h/a=1，λ=3)Fig.15 Comparison of seepage overturn moment

4 结论

本文将椭圆余弦波一阶理论有效引入到对浅水区透空墩柱的绕射波浪作用和波浪渗流作用的数学分析与计算之中，取得了浅水波作用下透空圆环单柱绕射波浪场和渗流压力场的解析解式，并据此对墩柱所受波浪力及渗流力进行了实算。结果表明，圆环柱外表面的适度透空对浅水波浪的直接作用和波浪渗流作用均有减弱效应，而前者的减弱效果更为明显，此外，对渗流倾覆力矩的减弱效应较浮托力更为明显。一定条件下，波浪渗流力与波浪直接作用力可能具有相同量级。而圆环墩柱内外柱半径比可适当取值，从而满足墩柱减载和稳定的双重要求。与Airy波浪理论计算结果相比，在一定条件下椭圆余弦波理论因计入了水波非线性因素的可能影响而更为合理和安全。

［1］邱大洪.波浪渗流力学［M］.北京:国防工业出版社，2006:193－203.

［2］WILLIAMS A N，LI W.Water wave interaction with an array of bottom-mounted surface-piercing porous cylinders［J］.Ocean Engineering，2000，27(8):841－866.

［3］曾昭銮，黄华，詹杰民，等.固立于可渗透海床上透空型圆柱的波浪渗流力［J］.中北大学学报:自然科学版，2010(4):327－330.

［4］邹志利，邱大洪.可渗可压缩海床上墩柱底面波浪渗流压力的计算［J］.海洋工程，1993，11(4):47－58.

［5］黄华，詹杰民.可渗透结构的非线性水波绕射作用问题［J］.中山大学学报:自然科学版，2006，45(1):25－28.

［6］ISAACSON M de St Q.Shallow wave diffraction around large cylinder［J］.Journal of Waterway Port Coastal And Ocean Div.，ASCE，1977，103(WW1):69 －82.

［7］苏铭德，潘宇.浅水中群柱受力状况的理论研究［J］.应用数学与力学，1987，8(4):367 －376.

［8］邱大洪，周援衡，贾影.浅水区圆柱墩群上的非线性波浪力［J］.海洋工程，1991，9(1):47－58.

［9］付静，黄华，朱庆勇.二层流体中透空大直径桩柱的波浪特性分析［J］.中山大学学报:自然科学版，2008，47(3):14－18.

［10］付静，黄华，詹杰民，等.二层海洋中波浪对透空圆柱群的绕射作用分析［J］.中山大学学报:自然科学版，2008，47(S2):76 －80.

［11］孙昭晨，邱大洪.作用于可渗可压缩海床上的墩柱底面上的波浪力［J］.海洋学报，1989，11(3):364－371.

［12］邱大洪，陈健.抛石基床上圆柱墩底部的波浪浮托力［J］.海洋学报，1994，16(1):114 －123.

［13］黄华，曾昭銮，詹杰民，等.二层海中立于可渗透海床上的大直径墩柱的波浪渗流力［J］.中山大学学报:自然科学版，2011，50(2):36－41.