3. 变桨距控制基础模型：空气动力学基础、桨叶受力分析、变桨距控制系统的数学模型

各位工程师朋友，欢迎来到第三章。这一章，咱们要啃一块硬骨头——变桨距控制的基础模型。

说实话，我刚开始接触风电那会儿，觉得变桨不就是转个角度嘛，有什么难的？后来在项目里栽了跟头才明白，不懂空气动力学和受力分析，你连PID参数都调不明白。嗯，咱们一步步来。

3.1 空气动力学基础：风是怎么“推”动叶片的？

先问大家一个问题：风真的在“推”叶片吗？

其实不完全是。风经过叶片时，主要产生两种力——升力和阻力。这和飞机机翼的原理很像，只不过叶片是竖着转的。

核心公式就一个，贝茨（Betz）定律：

P = 0.5 * ρ * A * V³ * Cp

其中：

• P：风轮捕获的功率（W）
• ρ：空气密度（kg/m³）——海拔高了，密度小了，出力就下降
• A：风轮扫掠面积（m²）——叶片越长，捕风能力越强
• V：风速（m/s）——注意是三次方！风速翻倍，功率变8倍
• Cp：风能利用系数——理论上限是0.593（贝茨极限）

我在项目里遇到过一件事：有个同事把Cp当成固定值0.48来算，结果仿真出来的功率曲线和实测差了20%。后来一查，他用的Cp是额定风速附近的，低风速段根本没校准。所以啊，Cp曲线一定要分段标定。

3.2 桨叶受力分析：叶片到底在扛什么？

叶片转起来，受力很复杂。我把它拆成三个主要方向：

1. 气动力：升力+阻力，作用在翼型上
2. 离心力：叶片旋转产生的径向拉力
3. 重力：叶片自重，方向随位置变化

你想想看，叶片每转一圈，重力方向就变一次。这会导致交变载荷，疲劳问题就是这么来的。

具体到变桨控制，我们最关心的是桨叶根部弯矩。公式长这样：

M_b = ∫(0→R) r * dF_aero(r) dr

其中：

• M_b：根部弯矩（Nm）
• r：到轮毂中心的距离（m）
• dF_aero(r)：r处的气动力微元（N）

说白了，就是每个叶片段产生的力，乘以力臂，再积分。这个弯矩直接决定了变桨执行器需要多大的扭矩。

3.3 变桨距控制系统的数学模型

好，有了空气动力学和受力分析的基础，咱们可以搭数学模型了。

变桨距控制系统，我习惯把它分成三个环节：

环节	输入	输出	典型模型
执行器（变桨电机+减速器）	角度指令 β_ref	实际桨距角 β	一阶惯性 + 延迟
气动环节	桨距角 β、风速 V、转速 ω	气动转矩 T_aero	非线性查表（Cp曲线）
传动链	气动转矩 T_aero、发电机转矩 T_gen	转速 ω	二阶质量-弹簧-阻尼

执行器模型，我常用这个：

G(s) = K / (τ*s + 1) * e^(-Td*s)

其中：

• K：增益（通常≈1）
• τ：时间常数（典型值0.1~0.3秒）
• Td：延迟时间（通信+计算，约20~50ms）

气动环节是非线性的，没法用简单传递函数描述。我一般用查表法：

T_aero = 0.5 * ρ * A * R * V² * Ct(λ, β)

其中Ct是转矩系数，和Cp的关系是：Ct = Cp / λ。

传动链模型，我建议用两质量块模型：

J_r * dω_r/dt = T_aero - K_s * (θ_r - θ_g) - D_s * (ω_r - ω_g)
J_g * dω_g/dt = K_s * (θ_r - θ_g) + D_s * (ω_r - ω_g) - T_gen

这里：

• J_r, J_g：风轮和发电机的转动惯量
• K_s, D_s：传动轴的刚度和阻尼
• θ_r, θ_g：风轮和发电机的角位置

为什么要用两质量块？因为单质量块模型忽略了传动轴的扭振。而扭振频率通常在1~3Hz，正好和变桨控制的带宽重叠，搞不好就会共振。

下面这张图，是我自己总结的变桨控制模型结构，你看一眼就明白了：

你看，整个模型就是：风速和桨距角进气动环节，算出转矩；转矩进传动链，算出转速；转速又反馈回去影响叶尖速比λ，从而影响气动环节。这是一个闭环。

我个人习惯，在搭建这个模型时，先用Simulink搭个简单的版本，跑通后再加细节。比如一开始可以忽略延迟，用一阶惯性代替执行器。等基本功能对了，再慢慢加延迟、加非线性、加限幅。

最后说一句：模型永远是对现实的近似。我在现场见过太多“仿真完美、现场崩盘”的例子。所以，模型建好后，一定要用实测数据校验。哪怕只校验几个工况点，也比纯仿真强百倍。

---