目录
- 气动力到载荷——BEM理论
-
- 一、为什么要研究 BEM 理论？
- 二、BEM理论的核心思想
- 三、第一部分：动量理论（Momentum Theory）
-
- 1 控制体模型
- 2 轴向诱导因子
- 3 Betz 极限
- 四、第二部分：叶素理论（Blade Element Theory）
-
- 1 叶片截面受力
- 2 相对风速
- 五、BEM理论的耦合
-
- BEM计算流程
-
- Step1 叶片离散
- Step2 初始假设
- Step3 计算速度三角形
- Step4 查翼型表
- Step5 计算气动力
- Step6 更新诱导因子
- Step7 迭代
- 六、气动力如何变成载荷？
-
- 叶素载荷分解
- 最终得到的载荷
- 七、为什么BEM在风电行业如此重要？
-
- 常见BEM修正
- 八、一个工程总结
- 九、BEM在风电软件中的应用
- 十、给风电工程师的一句话
气动力到载荷——BEM理论
——风电机组气动分析的核心方法（深入浅出版）
一、为什么要研究 BEM 理论？
在风电机组设计中，有一个核心问题：
风吹到叶片上 → 叶片产生气动力 → 气动力变成转矩 → 转矩变成结构载荷
工程师需要回答三个问题：
直接用 CFD（计算流体力学） 可以算，但计算量巨大。 因此工程界普遍采用一种 高效且准确的工程模型：BEM 理论。
BEM = Blade Element Momentum
它把两个理论结合起来：
- 动量理论（Momentum Theory）
- 叶素理论（Blade Element Theory）
从而实现 快速计算风机气动力和载荷。
二、BEM理论的核心思想
BEM 的基本思想其实很直观：
把一整片叶片切成很多“小叶片”，分别计算每一段的气动力，然后累加。
如下图所示：
每个叶片截面都有：
- 相对风速
- 攻角
- 升力
- 阻力
通过这些量，可以计算该截面的 气动力贡献。
最后把所有叶素相加，就得到：
- 转矩
- 推力
- 功率
- 结构载荷
三、第一部分：动量理论（Momentum Theory）
在研究叶片之前，我们先看一个更简单的问题：
一股空气流过一个“能量提取盘”会发生什么？
这就是 风机动量理论模型。
1 控制体模型
风机可以看成一个 作用盘（Actuator Disk）：
气流变化：
| 上游 | V∞ |
| 叶轮处 | V |
| 下游 | Vw |
空气通过风机后：
- 速度下降
- 动能被风机吸收
2 轴向诱导因子
BEM 中有一个非常重要的变量：
轴向诱导因子
a
=
V
∞
−
V
V
∞
a = \\frac{V_\\infty – V}{V_\\infty}
a=V∞V∞−V
意义：
- 表示风机让风速降低了多少
典型值：
| 无风机 | 0 |
| 正常风机 | 0.2~0.4 |
| 极限状态 | 0.5 |
3 Betz 极限
根据动量理论可以推导：
风机最多只能利用 59.3% 的风能
这个著名结论叫：
Betz Limit
C
p
m
a
x
=
16
27
=
0.593
C_p^{max} = \\frac{16}{27} = 0.593
Cpmax=2716=0.593
现实风机通常：
| 小风机 | 0.35 |
| 商用机组 | 0.45~0.5 |
四、第二部分：叶素理论（Blade Element Theory）
动量理论只考虑 整个风轮。
但叶片其实是 翼型结构。
因此需要空气动力学。
1 叶片截面受力
叶片截面气动力：
翼型会产生两种力：
| 升力 L | 垂直气流 |
| 阻力 D | 平行气流 |
公式：
升力：
L
=
1
2
ρ
V
r
e
l
2
c
C
L
L = \\frac{1}{2} \\rho V_{rel}^2 c C_L
L=21ρVrel2cCL
阻力：
D
=
1
2
ρ
V
r
e
l
2
c
C
D
D = \\frac{1}{2} \\rho V_{rel}^2 c C_D
D=21ρVrel2cCD
其中：
| ρ | 空气密度 |
| Vrel | 相对风速 |
| c | 叶片弦长 |
| CL | 升力系数 |
| CD | 阻力系数 |
2 相对风速
叶片不是静止的，而是在旋转。
因此叶片看到的风速是：
来流 + 旋转速度
速度三角形：
V
r
e
l
=
(
V
(
1
−
a
)
)
2
+
(
Ω
r
(
1
+
a
′
)
)
2
V_{rel} = \\sqrt{(V(1-a))^2 + (\\Omega r (1+a'))^2}
Vrel=(V(1−a))2+(Ωr(1+a′))2
参数：
| Ω | 角速度 |
| r | 半径 |
| a | 轴向诱导因子 |
| a’ | 切向诱导因子 |
五、BEM理论的耦合
现在问题来了：
动量理论得到：
- a
- a’
叶素理论需要：
- a
- a’
所以必须 迭代求解。
BEM计算流程
基本流程：
Step1 叶片离散
将叶片分为：
r1, r2, r3 … rn
Step2 初始假设
假设：
a = 0.3
a' = 0.01
Step3 计算速度三角形
得到：
- 相对风速
- 攻角
Step4 查翼型表
获得：
CL
CD
Step5 计算气动力
计算：
- 推力
- 转矩
Step6 更新诱导因子
根据动量理论：
更新
a
a'
Step7 迭代
直到：
|a_new – a_old| < tolerance
六、气动力如何变成载荷？
BEM不仅计算功率，还计算 结构载荷。
从叶素力可以分解：
| 法向力 | 推力 |
| 切向力 | 转矩 |
叶素载荷分解
计算：
法向力：
F
n
=
L
cos
ϕ
+
D
sin
ϕ
F_n = L \\cos \\phi + D \\sin \\phi
Fn=Lcosϕ+Dsinϕ
切向力：
F
t
=
L
sin
ϕ
−
D
cos
ϕ
F_t = L \\sin \\phi – D \\cos \\phi
Ft=Lsinϕ−Dcosϕ
最终得到的载荷
积分整个叶片：
T
=
∫
F
t
r
d
r
T = \\int F_t r dr
T=∫Ftrdr
得到：
转矩
T
h
r
u
s
t
=
∫
F
n
d
r
Thrust = \\int F_n dr
Thrust=∫Fndr
得到：
推力
这些载荷最终作用到：
| 叶根 | Blade Root Moment |
| 轮毂 | Hub Load |
| 塔架 | Tower Load |
也就是你在 Bladed / GT-Sim 中看到的：
- blade_root_Mxy
- tower_bottom_Mxy
- yaw_bearing_load
七、为什么BEM在风电行业如此重要？
BEM有三个巨大优势：
| 快 | 比CFD快1000倍 |
| 准 | 工程精度足够 |
| 易扩展 | 可加入修正模型 |
常见BEM修正
工程软件都会加入：
| Prandtl Tip Loss | 叶尖损失 |
| Glauert correction | 高诱导修正 |
| Dynamic stall | 动态失速 |
| Yaw correction | 偏航修正 |
这些模型保证：
BEM仍然是风机气动计算的核心方法。
八、一个工程总结
可以用一句话总结 BEM：
BEM = 用动量理论描述风轮整体 + 用叶素理论描述叶片局部
最终实现：
风速 → 气动力 → 转矩 → 功率 → 结构载荷
九、BEM在风电软件中的应用
所有主流软件都基于 BEM：
| Bladed | 整机载荷 |
| OpenFAST | 气动仿真 |
| QBlade | 叶片设计 |
| WT-Perf | 功率预测 |
甚至很多 CFD计算 也会用：
BEM作为初始解。
十、给风电工程师的一句话
理解 BEM 的意义在于：
你会真正理解 风 → 气动力 → 结构载荷 的传递链。
这对于理解：
- 极限载荷
- 疲劳载荷
- 风模型影响
至关重要。
