3. 波浪理论及其载荷计算
海洋平台所承受的环境载荷中,波浪载荷通常是最主要、也最复杂的控制性载荷。要计算波浪载荷,必须先回答两个问题:一是如何描述波浪本身的运动,二是如何把波浪运动转化为作用在结构上的力。前者属于波浪理论,后者属于载荷模型与适用条件判断。
一、线性波与非线性波理论要点
线性波理论(Airy波)假定波高相对波长很小、水质点运动为简谐运动,波面方程可写为 η = (H/2)·cos(kx − ωt),其中 H 为波高、k = 2π/L 为波数、ω 为圆频率。其色散关系为 ω² = gk·tanh(kd),d 为水深。线性波的最大优点是形式简单、可叠加,适合描述小振幅波浪和做频域分析,但它把波峰与波谷对称处理,无法反映实际海浪“峰陡谷坦”的特征。
当波高与水深、波长之比增大时,需采用非线性波理论。工程上常用的是斯托克斯高阶波(Stokes II、III、V)和椭圆余弦波(Cnoidal)。斯托克斯波适用于相对水深 d/L 较大(约大于 0.1)的深水至中等水深,通过高阶项修正波峰变尖、波谷变平;椭圆余弦波适用于浅水(d/L 较小),能描述波峰陡峭、波谷长而平且存在明显质量输移的浅水波。选择理论时,一般以 d/L 与 H/L 为判据:深水小波用线性波,深水大波用斯托克斯高阶波,浅水用椭圆余弦波,极浅水破碎带则需考虑破碎波模型。
二、Morison方程与绕射理论的适用条件
把波浪运动转换为结构受力,关键在于判断结构尺度与波长的相对关系。当构件直径 D 远小于波长 L(工程判据常取 D/L < 0.2)时,波浪绕过构件时几乎不发生明显散射,此时采用 Morison方程:单位长度水平力 f = 0.5·ρ·C_D·D·|u|·u + ρ·C_M·(πD²/4)·(du/dt),第一项为拖曳力(与速度平方相关,含阻力系数 C_D),第二项为惯性力(与加速度相关,含惯性力系数 C_M)。Morison方程本质上是把构件视为不扰动波场的“细长杆”,适用于导管架平台的弦杆、撑杆等小尺度构件。
当构件直径与波长可比(D/L 较大,通常 D/L > 0.2)时,构件会明显改变周围波场,产生绕射和反射,此时应采用绕射理论,通过求解满足结构物表面边界条件的势流绕射问题得到压力分布,再积分求总力。大型深水平台的立柱、重力式平台的沉箱等大尺度结构属于此类。介于两者之间的过渡区,工程上常做两种方法的对比核算,取较不利结果。
三、波浪载荷计算流程
实际计算通常按以下顺序推进:首先确定设计水位与设计波要素(设计波高、周期、波长),并按 d/L 与 H/L 选定波浪理论,算出设计波面下的水质点速度 u 与加速度 du/dt 沿水深的分布;其次按构件尺度比 D/L 选择 Morison 方程或绕射理论;然后对每个构件沿其长度积分,得到构件水平力与竖向力,再汇总到节点或整体坐标,形成总水平力、总竖向力和倾覆力矩;最后按规范进行载荷组合,与自重、风、流、冰等载荷叠加,用于整体强度与基础校核。
需要强调的是,Morison方程中的 C_D、C_M 并非纯理论常数,而依赖雷诺数、KC 数、表面粗糙度和海生物附着情况。工程中常取 C_D ≈ 0.6~1.2、C_M ≈ 1.5~2.0,海生物附着会显著增大拖曳力,需按规范修正。此外,波浪与流的联合作用、波浪相位对多构件叠加的影响,都是实际计算中必须处理的问题。
四、算例演示
设某导管架平台位于水深 d = 50 m 处,设计波高 H = 8 m、周期 T = 10 s。先由色散关系估算波长:取深水近似 L₀ = gT²/(2π) ≈ 156 m,考虑有限水深修正后 L ≈ 130 m。对直径 D = 1.2 m 的弦杆,D/L ≈ 0.009,远小于 0.2,应采用 Morison 方程;对直径 D = 30 m 的立柱,D/L ≈ 0.23,大于 0.2,应采用绕射理论。
以弦杆某深度处为例,若该点水质点速度 u = 1.5 m/s、加速度 du/dt = 0.8 m/s²,取 ρ = 1025 kg/m³、C_D = 1.0、C_M = 2.0,则单位长度拖曳力约为 0.5×1025×1.0×1.2×1.5² ≈ 1384 N/m,单位长度惯性力约为 1025×2.0×(π×1.2²/4)×0.8 ≈ 1855 N/m,合计约 3.24 kN/m。沿构件分段积分并叠加各构件贡献,即可得到整体波浪载荷。上述数值为教学示例的合理假设,实际工程须以规范系数和详细水动力分析为准。