第1章：波浪理论入门

各位同学好，我是老张。在海洋工程这行摸爬滚打十几年，今天咱们来聊聊波浪理论。说实话，波浪这东西看着简单，但真要把它算清楚，里面的门道可不少。

我个人习惯，做疲劳分析之前，先把波浪搞明白。为什么？因为波浪是海洋结构物最主要的载荷来源。你想想看，一个海上平台天天被波浪拍打，十年二十年下来，焊缝会不会出问题？这就是咱们这门课要解决的核心问题。

1.1 线性波浪理论（Airy波）

线性波浪理论，也叫Airy波理论，是波浪理论里最基础的一个。说白了，它把波浪简化成正弦波来处理。我刚开始做项目时，总觉得这太理想化了，但后来发现，很多工程问题用这个理论就够用了。

Airy波的核心假设有这么几条：

• 流体无黏性——忽略水的黏性影响
• 流动无旋——水质点运动没有旋转
• 波高远小于波长——这是线性化的前提
• 水深均匀——海底是平坦的

嗯，这里要注意，这些假设在实际海洋中并不完全成立。但为什么还要用？因为简单、好用、够用。我在南海某平台的项目中，就是用Airy波做的初步疲劳评估，结果和实测数据误差在15%以内，对于初步设计阶段来说，完全可以接受。

Airy波的速度势函数长这样：

φ = (gH / 2ω) * [cosh(k(z+h)) / cosh(kh)] * sin(kx - ωt)

其中：

• H — 波高
• ω — 角频率
• k — 波数
• h — 水深
• z — 垂向坐标（水面为0，向上为正）

有了速度势，水质点的速度和加速度就能算出来了。这些数据是后面做疲劳分析的基础。

1.2 波浪参数：波高、周期、水深

做波浪分析，三个参数你得烂熟于心：波高、周期、水深。这三者之间的关系，决定了波浪的类型。

参数	符号	典型范围	工程意义
波高	H	0.5~30m	决定波浪力大小
周期	T	3~20s	影响波浪与结构共振
水深	h	10~3000m	决定波浪理论选择

根据水深与波长的比值，波浪可以分为三类：

• 深水波：h/L > 0.5，海底不影响波浪运动
• 浅水波：h/L < 0.05，海底摩擦显著
• 中等水深波：介于两者之间

我曾经在北海做一个导管架平台项目，水深只有30米，但波浪周期长达15秒。你算算看，波长大概350米，h/L只有0.086，属于中等水深。这种情况下用深水波公式算出来的结果，误差能到30%以上。所以选对理论很重要。

1.3 波浪谱：JONSWAP和Pierson-Moskowitz

实际海洋中的波浪不是单一的正弦波，而是由无数不同频率、不同方向的波浪叠加而成。波浪谱就是描述这种能量分布的数学工具。

Pierson-Moskowitz谱（P-M谱）：

这是最早提出的波浪谱之一，适用于充分发展的风浪。它的形式比较简单：

S(ω) = (αg² / ω⁵) * exp[-β(ω₀/ω)⁴]

其中α=0.0081，β=0.74，ω₀=g/U，U是风速。

我个人习惯，在初步设计阶段用P-M谱就够了。但如果你做的是详细设计，特别是对疲劳寿命敏感的结构，我建议用JONSWAP谱。

JONSWAP谱：

JONSWAP谱是在P-M谱基础上加了峰值增强因子γ。它更适用于有限风区的波浪条件。我在墨西哥湾的一个FPSO项目中，实测数据表明JONSWAP谱比P-M谱更吻合，特别是在波浪能量集中的频段。

S_J(ω) = S_PM(ω) * γ^exp[-0.5((ω-ω_p)/(σω_p))²]

γ的典型值在1.0到7.0之间，γ越大，波浪能量越集中。嗯，这里要注意，γ=1.0时JONSWAP谱退化为P-M谱。

1.4 Python波浪模拟

理论讲完了，咱们来点实际的。用Python模拟波浪，其实就三步：生成频谱、离散化、叠加。

下面是我常用的一个代码模板：

import numpy as np
import matplotlib.pyplot as plt

def generate_wave_spectrum(Hs, Tp, gamma=3.3, n_freq=100):
    """
    生成JONSWAP谱
    Hs: 有效波高 (m)
    Tp: 谱峰周期 (s)
    gamma: 峰值增强因子
    """
    g = 9.81
    omega_p = 2 * np.pi / Tp
    omega = np.linspace(0.1, 3.0, n_freq)
    
    # P-M谱
    alpha = 0.0081
    S_pm = alpha * g**2 / omega**5 * np.exp(-1.25 * (omega_p/omega)**4)
    
    # 峰值增强因子
    sigma = np.where(omega <= omega_p, 0.07, 0.09)
    r = np.exp(-0.5 * ((omega - omega_p) / (sigma * omega_p))**2)
    S_j = S_pm * gamma**r
    
    return omega, S_j

# 生成波浪时程
def wave_time_series(omega, S, duration=600, dt=0.5):
    """
    从波浪谱生成时程
    """
    n = len(omega)
    d_omega = omega[1] - omega[0]
    t = np.arange(0, duration, dt)
    eta = np.zeros(len(t))
    
    for i in range(n):
        amplitude = np.sqrt(2 * S[i] * d_omega)
        phase = np.random.uniform(0, 2*np.pi)
        eta += amplitude * np.cos(omega[i]*t + phase)
    
    return t, eta

# 使用示例
omega, S = generate_wave_spectrum(Hs=5.0, Tp=10.0)
t, eta = wave_time_series(omega, S)

plt.figure(figsize=(12, 4))
plt.plot(t[:200], eta[:200])
plt.xlabel('时间 (s)')
plt.ylabel('波面高度 (m)')
plt.title('JONSWAP谱模拟波浪时程')
plt.grid(True)
plt.show()

知识体系总览

下面这张图，是我自己整理的波浪理论知识框架。做项目时，我经常拿出来对照，看看自己卡在哪一步。

这张图把本章的内容串起来了。从左到右，从理论到参数，再到谱分析，最后用Python实现。做疲劳分析时，你每一步都要走扎实。

好了，波浪理论入门就讲到这里。记住，理论是基础，但工程实践才是检验真理的唯一标准。下次咱们聊波浪对结构物的作用力，那才是真正开始算疲劳了。

---

公众号:蓝海资料掘金营，微信deep3321