3. 双馈电机数学模型（上）：三相静止坐标系下的电压方程、磁链方程、转矩方程
说实话，搞双馈电机这么多年，我始终觉得数学建模是道绕不过去的坎。很多人一上来就盯着坐标变换、矢量控制，结果连最基本的电压方程都写不利索。嗯，这就像盖楼不打地基，早晚要塌。
今天咱们就踏踏实实把三相静止坐标系下的数学模型啃下来。这部分内容虽然看起来公式多，但说白了就是三个方程：电压方程、磁链方程、转矩方程。你只要把这三大块理清楚，后面的坐标变换、控制策略就都好办了。
核心要点：三相静止坐标系下的数学模型是DFIG最原始的物理描述，所有后续的简化、变换都基于此。我个人习惯先把这组方程写熟，再去做坐标变换，这样心里有底。
3.1 双馈电机的物理结构回顾
在列方程之前，咱们先快速过一下DFIG的结构。双馈电机本质上是一台绕线式异步电机，转子侧有滑环和电刷，通过变频器给转子绕组通入交流励磁电流。
定子绕组直接挂在电网上，转子绕组通过背靠背变流器也挂在电网上。这就是「双馈」名字的由来——定子和转子都能向电网馈电。
你想想看，这种结构的好处是什么？转子侧变流器只需要处理转差功率，容量大约是电机额定功率的30%左右，成本一下子就降下来了。我在做1.5MW机组项目时，对这个优势体会特别深。
3.2 三相静止坐标系下的电压方程
好，咱们直接上干货。三相静止坐标系，就是咱们最熟悉的ABC坐标系。定子和转子各有三相绕组，一共六套绕组。
电压方程写起来其实很简单，就是基尔霍夫电压定律：绕组上的电压降等于电阻压降加上感应电动势。
定子侧电压方程：
u_A = R_s * i_A + dψ_A / dt
u_B = R_s * i_B + dψ_B / dt
u_C = R_s * i_C + dψ_C / dt转子侧电压方程：
u_a = R_r * i_a + dψ_a / dt
u_b = R_r * i_b + dψ_b / dt
u_c = R_r * i_c + dψ_c / dt这里有个细节要注意：定子侧的物理量下标用大写字母ABC，转子侧用小写字母abc。我刚开始学的时候经常搞混，后来养成习惯，写方程前先标清楚是哪一侧。
我的经验：实际工程中，定子电阻R_s通常很小，但千万别忽略它。我曾经在调试一台2MW机组时，因为忽略了定子电阻压降，导致空载电压计算偏差了3%，最后花了整整两天才找到原因。
把六个方程写成矩阵形式，看起来更紧凑：
[u_ABC] = [R_s] * [i_ABC] + d/dt [ψ_ABC]
[u_abc] = [R_r] * [i_abc] + d/dt [ψ_abc]其中电阻矩阵是对角阵：
[R_s] = diag(R_s, R_s, R_s)
[R_r] = diag(R_r, R_r, R_r)3.3 三相静止坐标系下的磁链方程
电压方程里出现了磁链ψ，那磁链又是怎么来的？说白了，磁链等于自感磁链加上互感磁链。双馈电机有六套绕组，它们之间都存在互感，所以磁链方程比电压方程复杂得多。
定子磁链方程：
ψ_A = L_s * i_A + M_s * i_B + M_s * i_C + M_sr * cos(θ_r) * i_a + M_sr * cos(θ_r - 120°) * i_b + M_sr * cos(θ_r + 120°) * i_c
ψ_B = M_s * i_A + L_s * i_B + M_s * i_C + M_sr * cos(θ_r + 120°) * i_a + M_sr * cos(θ_r) * i_b + M_sr * cos(θ_r - 120°) * i_c
ψ_C = M_s * i_A + M_s * i_B + L_s * i_C + M_sr * cos(θ_r - 120°) * i_a + M_sr * cos(θ_r + 120°) * i_b + M_sr * cos(θ_r) * i_c转子磁链方程：
ψ_a = L_r * i_a + M_r * i_b + M_r * i_c + M_sr * cos(θ_r) * i_A + M_sr * cos(θ_r + 120°) * i_B + M_sr * cos(θ_r - 120°) * i_C
ψ_b = M_r * i_a + L_r * i_b + M_r * i_c + M_sr * cos(θ_r - 120°) * i_A + M_sr * cos(θ_r) * i_B + M_sr * cos(θ_r + 120°) * i_C
ψ_c = M_r * i_a + M_r * i_b + L_r * i_c + M_sr * cos(θ_r + 120°) * i_A + M_sr * cos(θ_r - 120°) * i_B + M_sr * cos(θ_r) * i_C看着是不是有点晕？别急，我来帮你拆解一下。
这里面有几个关键参数：
- L_s：定子每相绕组的自感
- L_r：转子每相绕组的自感
- M_s：定子两相绕组之间的互感
- M_r：转子两相绕组之间的互感
- M_sr：定子和转子绕组之间的最大互感
- θ_r：转子位置角，即转子a相轴线与定子A相轴线之间的夹角
注意：定转子之间的互感系数包含cos(θ_r)项，这意味着互感是随转子位置变化的。这就是为什么三相静止坐标系下的方程是时变的，处理起来很麻烦。我们后面做坐标变换，目的就是为了消除这个时变性。
写成矩阵形式：
[ψ_ABC] = [L_ss] * [i_ABC] + [L_sr(θ_r)] * [i_abc]
[ψ_abc] = [L_rs(θ_r)] * [i_ABC] + [L_rr] * [i_abc]其中[L_sr(θ_r)]和[L_rs(θ_r)]互为转置，且都包含θ_r的三角函数。
3.4 三相静止坐标系下的转矩方程
转矩方程是电机数学模型中最有意思的部分。它把电磁能量转换成了机械运动。
双馈电机的电磁转矩可以用磁场储能对转子位置角的偏导来求，也可以直接用电流和磁链来表达。我这里给出工程上最常用的形式：
T_e = -n_p * M_sr * { [i_A * i_a + i_B * i_b + i_C * i_c] * sin(θ_r)
+ [i_A * i_b + i_B * i_c + i_C * i_a] * sin(θ_r + 120°)
+ [i_A * i_c + i_B * i_a + i_C * i_b] * sin(θ_r - 120°) }其中n_p是电机的极对数。
这个公式看起来复杂，但它的物理意义很清晰：电磁转矩是定转子电流相互作用产生的，而且与转子位置角θ_r密切相关。
避坑指南：我曾经在仿真时直接用这个公式算转矩，结果发现转矩波形有高频振荡。后来排查了半天，发现是电流采样中混入了噪声，导致sin项计算不准确。建议在实际应用中，先用滤波器处理电流信号，再代入转矩公式。
转矩方程还有另一种更简洁的表达方式：
T_e = n_p * (i_ABC^T * ∂L_sr(θ_r)/∂θ_r * i_abc)这种形式在理论推导中很常用，但实际计算时还是展开成上面的三角函数形式更方便。
3.5 三相静止坐标系模型的局限性
方程都列出来了，但你可能会问：这玩意儿能用吗？
说实话，直接用三相静止坐标系模型做控制，非常困难。原因有三：
- 时变性：定转子互感随θ_r变化，导致系统是时变的，控制器的设计变得极其复杂。
- 强耦合：六套绕组之间互相耦合，A相电流的变化会影响B相、C相，甚至影响转子侧的所有相。
- 非线性：转矩方程中包含电流的乘积项和三角函数，系统呈现强非线性。
所以，在实际工程中，我们几乎不会直接用三相静止坐标系模型来设计控制器。它的主要价值在于：
- 作为理论推导的起点
- 用于电磁场分析和有限元仿真
- 验证坐标变换后模型的正确性
我的建议：三相静止坐标系模型就像是一张「原图」，虽然看起来杂乱无章，但所有信息都在里面。你只有理解了这张原图，才能明白坐标变换到底在做什么——它是在「重新拍照」，让画面变得更清晰、更容易处理。
3.6 知识体系总览
为了帮你理清思路，我画了一张图，把本章的核心逻辑串起来：
从这张图可以看得很清楚：电压方程、磁链方程、转矩方程三者相互关联，共同构成了完整的数学模型。而它们共同的局限性——时变性、强耦合、非线性——正是我们下一章要做坐标变换的根本原因。
总结一下：三相静止坐标系下的数学模型虽然在实际控制中很少直接使用，但它是理解双馈电机本质的必经之路。我个人建议你把这三个方程亲手推导一遍，哪怕花上半天时间也值得。等你把坐标变换学完再回头看，会发现这些「原始方程」其实很美。
公众号:蓝海资料掘金营，微信deep3321