2. 双馈风机数学模型：三相静止坐标系下的DFIG模型
做双馈风机控制，第一步就是要把它的数学模型吃透。我当年刚入行时，总觉得数学推导太枯燥，直接跳过去看控制算法，结果仿真怎么调都不对。后来老老实实回来啃模型，才发现——模型是控制的地基，地基不稳，房子早晚要塌。
这一节，我们就来聊聊三相静止坐标系下的DFIG模型。说白了，就是直接用ABC三相的物理量来描述这台电机。虽然最终控制时我们会把它变换到旋转坐标系，但静止坐标系下的模型是理解一切的起点。
2.1 DFIG的物理结构回顾
双馈风机，全称是双馈异步发电机。它的定子直接挂电网，转子通过背靠背变流器接电网。你想想看，这结构有什么特别之处？
- 定子：三相绕组，直接连接工频电网
- 转子：三相绕组，通过滑环和碳刷引出，接变流器
- 转轴：连接风力机，传递机械功率
嗯，这里要注意：转子绕组是绕线式的，不是鼠笼式。这个区别很关键，因为绕线式转子允许我们从外部注入励磁电流，从而实现变速恒频运行。
核心要点：DFIG的定子和转子都是三相绕组，都可以建立磁链。定子磁链由电网电压决定，转子磁链由变流器控制。两者通过互感耦合在一起。
2.2 三相静止坐标系下的电压方程
我们先写电压方程。我个人习惯把定子和转子分开写，这样逻辑更清晰。
定子电压方程（下标s表示定子）：
u_sa = R_s * i_sa + dψ_sa / dt
u_sb = R_s * i_sb + dψ_sb / dt
u_sc = R_s * i_sc + dψ_sc / dt转子电压方程（下标r表示转子）：
u_ra = R_r * i_ra + dψ_ra / dt
u_rb = R_r * i_rb + dψ_rb / dt
u_rc = R_r * i_rc + dψ_rc / dt这里R_s和R_r分别是定子和转子的每相电阻。ψ是磁链。看起来很简单对吧？但实际用起来，麻烦在后面——磁链方程。
我的经验：在项目中调试DFIG模型时，最容易出错的就是电压方程中的正负号。特别是转子侧，因为转子在旋转，感应电动势的方向容易搞反。我建议你写方程时，先画一个等效电路图，标清楚参考方向，再下笔。
2.3 三相静止坐标系下的磁链方程
磁链方程才是DFIG模型的精髓。为什么？因为定子和转子之间有相对运动，互感是时变的。
定子磁链：
ψ_sa = L_s * i_sa + L_m * [i_ra * cos(θ_r) + i_rb * cos(θ_r + 120°) + i_rc * cos(θ_r - 120°)]
ψ_sb = L_s * i_sb + L_m * [i_ra * cos(θ_r - 120°) + i_rb * cos(θ_r) + i_rc * cos(θ_r + 120°)]
ψ_sc = L_s * i_sc + L_m * [i_ra * cos(θ_r + 120°) + i_rb * cos(θ_r - 120°) + i_rc * cos(θ_r)]转子磁链：
ψ_ra = L_r * i_ra + L_m * [i_sa * cos(θ_r) + i_sb * cos(θ_r - 120°) + i_sc * cos(θ_r + 120°)]
ψ_rb = L_r * i_rb + L_m * [i_sa * cos(θ_r + 120°) + i_sb * cos(θ_r) + i_sc * cos(θ_r - 120°)]
ψ_rc = L_r * i_rc + L_m * [i_sa * cos(θ_r - 120°) + i_sb * cos(θ_r + 120°) + i_sc * cos(θ_r)]其中：
- L_s = L_ls + L_m：定子自感（漏感+互感）
- L_r = L_lr + L_m：转子自感（漏感+互感）
- L_m：定转子之间的互感幅值
- θ_r：转子位置角（电角度）
看到这些cos函数了吗？这就是时变性的来源。转子一转，θ_r就变，互感系数也跟着变。我曾经在仿真中直接用这个模型做控制，结果发现控制器参数根本没法整定——因为系统是时变的，固定参数的PI控制器搞不定。
避坑指南：三相静止坐标系下的DFIG模型，虽然物理意义清晰，但直接用于控制设计非常困难。因为互感矩阵是时变的，导致系统是非线性的。我建议你用它做仿真验证，但控制算法设计一定要用旋转坐标系下的模型。
2.4 转矩方程
电磁转矩是机电能量转换的核心。在静止坐标系下，转矩表达式为：
T_e = -p * L_m * [ (i_sa * i_ra + i_sb * i_rb + i_sc * i_rc) * sin(θ_r)
+ (i_sa * i_rb + i_sb * i_rc + i_sc * i_ra) * sin(θ_r + 120°)
+ (i_sa * i_rc + i_sb * i_ra + i_sc * i_rb) * sin(θ_r - 120°) ]其中p是极对数。这个式子看起来很复杂，但说白了就是定转子电流相互作用产生转矩。你想想看，如果没有转子电流，转矩就是零——所以DFIG的励磁完全由转子侧变流器控制。
2.5 运动方程
机械部分相对简单：
J * dω_m / dt = T_m - T_e - B * ω_m其中：
- J：转动惯量
- ω_m：机械角速度
- T_m：风力机输入的机械转矩
- T_e：电磁转矩
- B：阻尼系数
这个方程是所有电机都通用的，没什么特别。但要注意，DFIG的转速范围通常是0.7~1.3倍同步转速，所以ω_m是变化的。
2.6 模型总结与SVG框架图
好了，我们把三相静止坐标系下的DFIG模型串起来。它由四部分组成：
- 电压方程：描述电压与电流、磁链的关系
- 磁链方程：描述磁链与电流的关系（含时变互感）
- 转矩方程：描述电磁转矩与电流的关系
- 运动方程：描述转速与转矩的关系
这四部分耦合在一起，构成了完整的DFIG数学模型。下面这张图可以帮你理清它们之间的关系：
从这张图可以看出，四个方程是相互耦合的。电压方程需要磁链信息，磁链方程需要电流信息，转矩方程需要电流和位置信息，运动方程需要转矩信息。环环相扣，缺一不可。
实用建议：在实际仿真中，我通常这样搭建DFIG模型：先用磁链方程根据电流计算磁链，再用电压方程根据磁链计算电流的导数，然后积分得到电流。转矩和运动方程则作为辅助输出。这个流程在Simulink或PLECS里都很容易实现。
2.7 小结
三相静止坐标系下的DFIG模型，虽然看起来公式多、耦合复杂，但它是最接近物理实际的描述方式。我建议你花点时间亲手推导一遍这些方程，特别是磁链方程中的互感项。只有真正理解了为什么互感是时变的，你才能明白为什么要做坐标变换。
下一节，我们会把静止坐标系下的模型变换到旋转坐标系。到那时，时变的互感会变成常数，控制设计会变得简单很多。但请记住——旋转坐标系下的模型，本质上就是静止坐标系模型经过数学变换得到的。地基打牢了，上层建筑才能稳固。
公众号:蓝海资料掘金营，微信deep3321