波浪荷载载各种工程中的确定

波浪荷载在各种工程中的确定在海洋工程中，无论是在石油钻井平台还是跨海工程，波浪荷载对结构的破坏都是不容忽视的因素。

在海上大跨度桥梁的建设中，无论是施工过程还是整体设计，波浪荷载的研究都有重大工程意义，特别是对于诸如斜拉桥、悬索桥桥塔等大型墩式结构，更是如此。

波浪力的计算需要两方面理论的支持：波浪运动理论及波浪荷载计算理论。

前者研究波浪的运动，后者在已知波浪运动的前提下计算波浪对水中物体的作用。

对于规则波，常采用的波浪运动理论有Airy 理论、Stokes 理论、椭圆余弦波以及孤立波理论。

Airy 理论以静水面代替波面，适用于振幅较小、水深较大的情况；Stokes 理论可以考虑波高的2阶以及更高阶项，Airy 理论可认为是Stokes 的1阶形式；椭圆余弦波计算较为繁琐，工程运用仍较少；孤立波理论用于考虑孤立波，即水质点相对水体移动的非振动波。

关于波浪荷载计算理论，不同的结构形式是不同的。

而小直径桩的波浪荷载计算主要采用试验测量及经验分析的方法。

其中，使用最广泛的是Morrison 于1952年提出的莫里森公式，这一公式本身以及有关的试验测量理论和测量资料，都有了很大的进展，已被许多国家的设计规范所采纳。

下面我将对波浪荷载理论及其在近海结构、跨海结构、钻井平台结构中的运用作简要叙述。

1 常用的波浪运动理论1.1 微幅波理论微幅波理论是应用势函数来研究波浪运动的一种线性波浪理论。

（1）水深无限时推进波的势函数：sin()2kz gH e kx t φωω=- H 为波高，ω为波浪圆频率，2T πω=， k 为波数，2k L π=。

在无限水深的推进波中波周期T 与波长L 0不是独立的，他们之间具有一定的关系：200022gT L L gT c T ππ====0c 为波速。

（2）水深有限时推进波的势函数：()sin()2gH chk d z kx t chkdφωω+=⋅- 在有限水深的推进波中波周期T 与波长L 的关系为：222gT L thkd L gT c thkd T ππ====假定波浪在浅水中推进时，其波周期T 保持不变，则：00L c thkd L c == 它说明了在微幅波理论适用的范围内，波浪由深水向浅水推进时的波长与波速变化规律。

1.2 有限振幅波理论有限振幅波理论仍假定波浪运动为势运动，可以用速度势来描述波浪运动的状态。

基本方程：222222220x x z z φψφψ∂∂∂∂+=+=∂∂∂∂ φ和ψ为波浪运动的速度势和流函数。

z 坐标垂直于波得速方向，x 坐标沿流速方向。

由基本方程和边界条件可求出波速为：22414221121126[1()2sh ch ch g c k J o J k ch sh ++=⋅++除了上述两种最常用的理论之外，还有孤立波理论、椭圆余弦波理论等理论方法值得研究，详述请参考邱大洪编著的《波浪理论及其在工程上的应用》一书，书中2至4章有详细叙述。

2 波浪荷载的计算理论2.1 直墙上的波浪力（1）直墙上的立波浪力二阶浅水立波的波压力公式：222()1()sin 2[1(32)cos 2]41(csc 22cos 2)2ch z kd ch z kd p z A A cth kd cth kd chkd chkdA h kd cth kd ωτωτωτ++=-+-+-++令 kd →∞，得二阶深水立波的波压力公式：22211sin cos 2(1cos 2)24z z p z Ae A A e ωτωτωτ=-++-+ （2）直墙上的破波浪力1）远破波的波压力计算在静水面以上高度为H （推进波的波高）处得波压力为0，静水面处得波压力P s 为12s p K K H γ=其中 γ——海水容重1K ——与海底坡度i 有关的系数2K ——与坡坦有关的系数在静水面以上的波浪力的分布按直线变化。

在静水面以下z=H/2处得波浪附加应力P s 为0.7z s p p =在墙底处的波浪附加压力P d 为当d/H ≤1.7时，0.6d s p p =当d/H>1.7时，0.5d s p p =墙底面上的波浪浮托力P u 为2d u bp p u = （b 为墙底宽度，u 为波浪浮托力分布图的折减系数，可取0.7） 2）近破波的波压力计算静水面处得波压力P s 为：当2/3≥d 1>1/3时，111.25(1.80.16)(10.13)s H H p H d d γ=--； 当1/3≥d 1≥1/4时，1111.25[(13.936.4)(0.67) 1.03](10.13)s d H H p H d d d γ=--+-。

在墙底处的波浪附加压力P d 为0.6d s p p =在计算单位长度堤身上的近破波的总波浪力P 时运用下面的公式：当2/3≥d 1>1/3时，11.25(1.90.17)s H p H d γ=-； 当1/3≥d 1≥1/4时，1111.25[(14.838.8)(0.67) 1.1]s d H p Hd d d γ=--+。

2.2 墩柱上的波浪力（1）Morison 法这个方法的基本假定是认为当墩柱尺度与波长相比较小时，墩柱的存在并不影响波动场，故作用在墩柱上的波浪力，除了与墩柱尺度相关外，取决于未被墩柱扰动的波动场内的墩柱轴线处得水质点运动速度和加速度，计算规则波对小直径桩柱及柱群的作用力时，通常采用Morison 提出的波浪力方程。

Morison 方程是一种带有经验性质的半理论公式，它包含两项，即惯性力和速度力，惯性力项的形式与无粘性流体的波动理论的解相同，而速度力项的形式则与稳定流中的物体上产生的阻力相仿，这个方程的运用，要求墩柱直径D 与波长L 之比较小，在一般情况下，当D/L 小于等于0.15时才较适用。

（2）MacCamy 和Fuchs 的绕射理论该理论假定流体是无粘性的，运动时有势的，并利用了线性化得自由水面边界条件，故只在波动幅度响度较小的情况下才能适用，当然当D/L 大于0.25时，非线性对波浪力的影响一般小于5%，这样的差别通常在工程设计中是可以允许的。

对于无粘性假设，实验表明，当波高H 与墩柱直径之比H/D 小于等于1.0时，由于流体粘性所引起的阻力对波浪力的影响一般不超过5%。

因此可以认为，线性化的绕射理论的适用范围是H/D 小于等于1.0。

综上，两种方法的适用范围概括如下：1） 当H/D>1.0且H/D>0.15时，两方法都不适用；2） 当H/D<1.0且H/D<0.15时，运用Morison 方程，而且只记其中的惯性力项；3） 当H/D<1.0且H/D>0.15时，可采用不考虑流体粘性效应的波浪绕射理论；4）当H/D>1.0且H/D<0.15时，采用Morison方程。

各种墩柱形式下的波浪力计算详见《波浪理论及其在工程上的应用》书中的第八章。

3 具体实例3.1 东海大桥70+m跨低墩基础波浪力（1）计算基础波浪力东海大桥非通航孔桥主要采用60m、70m跨度，基础数量庞大。

跨海大桥桥墩基础除承受上部结构传来的荷载以外，还承受巨大的海洋波浪荷载。

在下部结构设计中，波浪产生的水平力为主要设计荷载。

在设计中有必要对作用于基础的波浪荷载做深入细致的研究。

对于跨海大桥桥墩基础，波浪荷载主要为作用于墩身、承台和桩的水平力；其次还有向上的垂直力，即浮托力。

平均低水位：-1.34m；平均高水位：1.86m；平均水位：0.23m；设计低水位：2.48m；重现期(50年一遇水位：3.60m；重现期100 年一遇水位：3.70m。

此设计采用重现期(50 年一遇水位，校核采用重现期100年一遇水位。

东海大桥桥址区沿波浪方向的设计波浪要素见下表。

桥址区主要浪向设计波浪要素桩基设计波高按累积频率1%波高计。

波浪力按Morison半经验公式计算，本桥墩基础各构件尺度较小，按孤立建筑物计算其波浪力，当迎波面宽度与波长比D/L<0.2时，采用流体的绕流理论计算；当D/L>0.2 时，考虑流体不连续而采用流体的绕射理论计算，计算中认为构件全宽上流体压力梯度不变。

波浪产生的与波浪平行的正向力俗称波浪力，其由速度分力和惯性分力组成，其计算公式如下：速度分力212D DrDHP C K =⋅⋅力矩2max32D DrDH LM C Kπ=⋅⋅惯性分力2122MrAHP C K =⋅⋅力矩2Im44ax MrAH LM C Kπ=⋅⋅其中，C D 为速度力系数；C M 为惯性力系数；r 为水容重；H 为计算波高；D 为构件迎波面宽度；A 为构件的截面积；L 为波长。

作用于构件的正向水平总波浪力P 为P D 和P 1的合力，在D/L<0.2时按下式计算：合力 max Im cos cos sin D ax P P t t P t =∂∂-∂合力矩 22max max Im m (10.25/)D ax D ax M M M M =+⋅在D/L>0.2时，圆形柱体主要受惯性力的作用，最大水平总波浪力为：合力 max Im ax P P =合力矩 max max I M M =（2）计算水流力2/(2)P krAV g =其中k 为物体形状系数；r 为水容重；A 为构件的截面积；V 为水流速度；g 为常数9.8。

跨海大桥工程基础设计必须考虑波浪荷载的作用，桥墩波浪物理模型试验表明，《海港水文规范》中波浪力计算方法对桥墩基础的小尺度圆形承台而言，基本上反映了作用于结构上波浪力的实际情况，但对单桩来说，计算值偏小，所以在桥梁基础设计中，单桩波浪力和群桩系数需认真考虑。

3.2半潜式平台波浪载荷计算3.2.1半潜式平台波浪载荷计算(1)计算条件平台主尺度和工作环境是浮体长度90. 0 m ；浮体宽度14. 0 m ；浮体高度6. 0 m ；两浮体间距50. 00 m ；立柱直径<9. 0 m ；立柱中心纵向间距31. 0 m ；立柱中心横向间距50. 00 m ；上部平台长度72. 0 m ；上部平台宽度63. 7 m ；上部平台高度5. 2 m ；主甲板高( 距基线) 30. 0 m; 上甲板高( 距基线) 35. 2 m ；最大工作水深200 m ；最大风速51. 6 m/ s ；最大波高18 m ；潮流流速3. 00 kN ；拖航排水量16 080 t ；拖航吃水6 m 。

（2）计算模型平台结构模型的总坐标系原点位于2 个浮体底部的对称中心, X 轴沿着平台纵向指向首部, Y 轴从右舷指向左舷, Z 轴则以垂直向上为正。

海洋平台总体结构十分复杂, 若完全按照其实际情况进行分析是无法实现的, 当然也没有必要, 因此需建立简化力学模型。