漂浮式海上风力发电机：波浪载荷的主动陷波阻尼控制

漂浮式海上风力发电机（Floating Offshore Wind Turbines, FOWT）是深海可再生能源开发的战略制高点。不同于近浅海（< 50 m）直接打入海床的传统单桩固定式风机，以 NREL OC3-Hywind 单立柱浮体（Spar-buoy）为代表的深海浮体平台常年服役于 100 至 300 米以上的深远海域。

深远海风能资源虽强劲平稳，却带来了严峻的流体结构耦合动力学挑战。深吃水圆柱浮体在俯仰回转方向（俯仰角 $\theta$ ）物理上呈现为巨大的水下单摆系统。由于周围海水对深海浮体的辐射阻尼极为微弱（阻尼比 $\zeta \approx 0.03$ ），当海浪主波频接近结构固有频率时，极易激发复合型谐振摆动，严重威胁塔架的结构疲劳寿命。

本案例基于简化的单自由度（1-DOF）Spar 浮体俯仰动力学模型，深入剖析传统 PD 控制为何在海浪激励下必然失效，并展示如何利用二阶陷波滤波器（Notch Filter）化解波频扰动与结构谐振之间的频域冲突。

物理建模与参数

浮体平台参数基于经典的 NREL OC3-Hywind 基准模型。在俯仰方向（船舶力学标准第 5 自由度），动力学遵循牛顿-欧拉转动平衡方程：

I_{\text{tot}} \ddot{\theta} + B_{55} \dot{\theta} + C_{55} \theta = M_{\text{wave}} - M_{\text{ctrl}}

参数	物理意义	基准值
$I_{55}$	钢结构自身俯仰转动惯量	$5.00 \times 10^{10}\ \text{kg}\cdot\text{m}^2$
$A_{55}$	流体附加质量转动惯量（随动水体）	$1.70 \times 10^{10}\ \text{kg}\cdot\text{m}^2$
$I_{\text{tot}}$	总体虚拟俯仰惯量 ( $I_{55} + A_{55}$ )	$6.70 \times 10^{10}\ \text{kg}\cdot\text{m}^2$
$B_{55}$	水动力辐射阻尼系数	$8.65 \times 10^8\ \text{N}\cdot\text{m}\cdot\text{s/rad}$
$C_{55}$	静水恢复刚度（稳心与重力力矩）	$3.111 \times 10^9\ \text{N}\cdot\text{m/rad}$
$f_n$	浮体俯仰固有共振频率	$0.0343\ \text{Hz}\ (\omega_n = 0.2155\ \text{rad/s},\ T_n \approx 29.2\ \text{s})$
$\zeta$	开环天然阻尼比	$\approx 0.03$ ( $Q \approx 16.7$ ，共振尖峰高达 $+24.5\ \text{dB}$ )

自由度定义（船舶工程 "55" 命名法）

在船舶工程和海洋流体动力学中，空间浮体具有 6 个标准自由度：

Surge（纵荡 / 沿 X 轴平动）
Sway（横荡 / 沿 Y 轴平动）
Heave（垂荡 / 沿 Z 轴升沉）
Roll（横摇 / 绕 X 轴转动）
Pitch（纵摇俯仰 / 绕 Y 轴转动）
Yaw（首摇 / 绕 Z 轴偏航）

下标 $55$ 表示第 5 自由度（俯仰运动）在第 5 自由度上引起的力矩反力。

传递函数与频域分离特征

输入总力矩到浮体俯仰角 $\theta$ 的开环传递函数为：

G(s) = \frac{\Theta(s)}{M(s)} = \frac{1}{I_{\text{tot}} s^2 + B_{55} s + C_{55}}

频域特征鲜明地呈现出两组迥异的频段分布：

浮体自身结构俯仰固有谐振频率位于 $\omega_n = 0.2155\ \text{rad/s}$ ( $f_n = 0.0343\ \text{Hz}$ )。
大洋主导海浪能谱峰值频率集中在 $f_p \approx 0.10\ \text{Hz}$ 附近（ $\omega_p = 0.628\ \text{rad/s}$ ，涌浪周期 $T_p \approx 10\ \text{s}$ ）。

由于海浪主频大约是浮体共振频率的 3 倍（ $\omega_p / \omega_n \approx 2.91$ ），系统天然存在频域分离断层：

结构固有谐振为超低频（ $0.034\ \text{Hz}$ ）。
海浪涌浪激励为相对高频（ $0.10\ \text{Hz}$ ）。

为何常规 PD 反馈控制必然失效

常规 PD 控制器依赖微分项（ $K_d s$ ）来提供人工阻尼。微分环节的幅频增益随频率以 $+20\ \text{dB/decade}$ 快速放大：

|K_d \cdot j\omega| = K_d \cdot \omega

在海浪频率（ $0.10\ \text{Hz}$ ）处，微分项对每弧度的波浪致倾斜给出 $\sqrt{K_p^2 + (K_d\omega)^2} \approx 561\ \text{MN}\cdot\text{m}/\text{rad}$ 的反力矩，尽管波浪并不是它需要阻尼的对象。在本课的波浪（ $20\ \text{MN}\cdot\text{m}$ ）下，峰值力矩仅约 $1.25\ \text{MN}\cdot\text{m}$ ——为 $\pm 50\ \text{MN}\cdot\text{m}$ 限值的 $2.5\%$ ，只有当波频倾角达到约 $5.1^\circ$ 时才会触及。这些力矩几乎全部花在波浪频率上，对共振毫无帮助；不过 PD 确实把共振阻尼从 $\zeta = 0.030$ 提高到约 $0.054$ 。

接入二阶陷波滤波器的智能主动阻尼

解决方案是在反馈控制回路中串联一个针对海浪优势频率（ $f_0 = 0.10\ \text{Hz}$ ）的二阶陷波滤波器（Notch Filter） $N(s)$ ：

N(s) = \frac{s^2 + \omega_0^2}{s^2 + 2\zeta_n \omega_0 s + \omega_0^2} \quad \text{其中 } \omega_0 = 2\pi f_0

在海浪主频处（ $s = j\omega_0$ ），分子严格为零：

|N(j\omega_0)| = 0 \quad (\text{理论衰减 } -\infty\ \text{dB})

该滤波器使控制器对一阶波浪运动“失明”，让浮体顺从地随波峰起伏：在本课中，峰值力矩和耗散能量下降约 70%。这并非没有代价：在结构共振点（ $0.034\ \text{Hz}$ ）处 $|N| \approx 0.93$ ，但滤波器使反馈滞后 $21^\circ$ ，因而 PD 所增加的阻尼缩减约 70%（闭环 $\zeta \approx 0.037$ 带陷波， $0.054$ 不带， $0.030$ 开环）。

降阶边界与教学可视化

本模型聚焦于解耦的单自由度俯仰姿态。旋转风轮、系泊缆绳与波浪起伏构成逼真的物理背景，在保证直观教学体验的同时免除高阶流体结构多体仿真的庞大开销。

参阅相关组件

StateSpace：Spar 浮体水动力学的紧凑二状态状态空间形式。
TransferFunction：二阶陷波滤波器传递函数实现。
Integrator：角加速度至俯仰角的两阶积分链路。
Saturation：执行机构最大力矩物理硬饱和。
BodePlot：伯德图频域陷波凹坑分析。
uPlotDisplay：角位移、角速度与控制力矩实时多通道曲线。