3. 异步电机数学模型：从电压方程到dq坐标系
各位同学，今天我们来啃一块硬骨头——异步电机的数学模型。
说实话，我刚入行那会儿，看到这一堆方程就头大。心想：不就是个电机嘛，转就完了，搞这么复杂干嘛？直到我在项目里被电流震荡、转矩脉动折磨得死去活来，才明白——不懂数学模型，你连问题出在哪都不知道。
好，咱们一步步来。
3.1 三相静止坐标系下的电机方程
先看最原始的模型。三相异步电机，定子有三相绕组，转子也有三相绕组（绕线式或笼型等效）。
嗯，这里要记住一个核心思想：电机本质上就是一个电磁耦合系统。电压加在绕组上，产生电流，电流产生磁链，磁链变化又感应出反电动势。说白了，就是电生磁、磁生电的循环。
3.1.1 电压方程
定子三相电压方程，写成矩阵形式：
[u_A] [R_s 0 0 ] [i_A] d [ψ_A]
[u_B] = [0 R_s 0 ] [i_B] + — [ψ_B]
[u_C] [0 0 R_s] [i_C] dt[ψ_C]转子三相电压方程类似，只是转子绕组是短路的（笼型电机），所以转子电压为零：
[0 ] [R_r 0 0 ] [i_a] d [ψ_a]
[0 ] = [0 R_r 0 ] [i_b] + — [ψ_b]
[0 ] [0 0 R_r] [i_c] dt[ψ_c]这里R_s是定子电阻，R_r是转子电阻。注意转子侧的量都是折算到定子侧的，这个折算系数我在后面会讲。
3.1.2 磁链方程
磁链方程就复杂了。因为定转子之间、各相之间都有互感，而且互感还随转子位置变化。
[ψ_s] [L_s L_sr] [i_s]
[ψ_r] = [L_rs L_r ] [i_r]展开写就是6×6的矩阵：
ψ_A = L_AA·i_A + L_AB·i_B + L_AC·i_C + L_Aa·i_a + L_Ab·i_b + L_Ac·i_c
ψ_B = L_BA·i_A + L_BB·i_B + L_BC·i_C + L_Ba·i_a + L_Bb·i_b + L_Bc·i_c
...（共6个方程）你看，光是写出来就够头疼的。而且L_Aa这些互感系数是转子位置角θ的函数：
L_Aa = L_m·cos(θ)
L_Ab = L_m·cos(θ + 120°)
L_Ac = L_m·cos(θ - 120°)为什么会这样？因为转子在转，定子A相和转子a相之间的相对位置一直在变，互感自然跟着变。
3.1.3 转矩方程
电磁转矩可以从能量角度推导：
T_e = n_p · L_m · [(i_A·i_a + i_B·i_b + i_C·i_c)·sin(θ) +
(i_A·i_b + i_B·i_c + i_C·i_a)·sin(θ+120°) +
(i_A·i_c + i_B·i_a + i_C·i_b)·sin(θ-120°)]看着就复杂，对吧？而且同样包含sin(θ)项。
3.1.4 运动方程
这个相对简单：
T_e - T_L = J · (dω_m/dt) + B · ω_m其中T_e是电磁转矩，T_L是负载转矩，J是转动惯量，B是阻尼系数，ω_m是机械角速度。
运动方程描述的是电机轴上的力学平衡。说白了就是：电磁转矩减去负载转矩，剩下的用来加速。
3.2 为什么要变换到dq坐标系？
刚才那一堆方程，你看出问题了吗？
三个字：太复杂。
三相静止坐标系下，电压方程是6阶的，磁链方程里全是时变系数，转矩方程里三角函数满天飞。这样的模型，别说做实时控制，就是离线仿真都费劲。
那怎么办？
聪明的前辈们想了个办法：坐标变换。
把三相静止的ABC坐标系，变换到两相旋转的dq坐标系。这样一来，时变的互感系数变成了常数，交流量变成了直流量。
你想想看，控制直流量是不是比控制交流量容易多了？PID调参都简单不少。
3.3 dq坐标系下的数学模型
经过Clark变换（3s/2s）和Park变换（2s/2r），我们得到dq坐标系下的方程。
3.3.1 dq坐标系下的电压方程
u_sd = R_s·i_sd + dψ_sd/dt - ω_e·ψ_sq
u_sq = R_s·i_sq + dψ_sq/dt + ω_e·ψ_sd
0 = R_r·i_rd + dψ_rd/dt - (ω_e - ω_r)·ψ_rq
0 = R_r·i_rq + dψ_rq/dt + (ω_e - ω_r)·ψ_rd注意看，这里出现了ω_e·ψ_sq和ω_e·ψ_sd项，这叫旋转反电动势。它代表的是坐标系旋转带来的耦合项。
我个人习惯把这两项叫做「交叉耦合项」。在高速时，这两项的影响很大，如果不做解耦补偿，电流环会很难调。
3.3.2 dq坐标系下的磁链方程
这个就清爽多了：
ψ_sd = L_s·i_sd + L_m·i_rd
ψ_sq = L_s·i_sq + L_m·i_rq
ψ_rd = L_m·i_sd + L_r·i_rd
ψ_rq = L_m·i_sq + L_r·i_rq所有电感系数都是常数！L_s是定子自感，L_r是转子自感，L_m是互感。没有三角函数，没有时变参数。
嗯，这里要注意：L_s = L_ls + L_m，L_r = L_lr + L_m。L_ls和L_lr分别是定转子漏感。
3.3.3 dq坐标系下的转矩方程
转矩方程也变得极其简洁：
T_e = n_p · L_m · (i_sq·i_rd - i_sd·i_rq)或者用磁链表示：
T_e = n_p · (ψ_sd·i_sq - ψ_sq·i_sd)更常用的形式是：
T_e = n_p · (ψ_rd·i_sq - ψ_rq·i_sd) · (L_m/L_r)如果采用转子磁场定向（把d轴对准转子磁链方向），那么ψ_rq = 0，转矩方程进一步简化为：
T_e = n_p · (L_m/L_r) · ψ_rd · i_sq你看，转矩只和转子磁链ψ_rd、定子q轴电流i_sq成正比。这不就和直流电机一样了吗？
3.3.4 dq坐标系下的运动方程
运动方程形式不变：
T_e - T_L = J · (dω_m/dt) + B · ω_m只是T_e现在用dq量来计算了。
3.4 知识体系总览
下面这张图，是我画的本章节的知识结构。你可以把它当作一个地图，随时回来看看自己走到哪了。
3.5 小结与工程建议
好了，总结一下本章的核心要点：
- 三相静止模型：物理意义清晰，但方程复杂，含时变参数，不适合控制
- 坐标变换：Clark + Park，把交流量变直流量，时变参数变常数
- dq模型：电压方程有交叉耦合项，磁链方程简洁，转矩方程线性化
- 转子磁场定向：进一步简化，转矩只由ψ_rd和i_sq决定
1. 参数不准（尤其是转子时间常数）
2. 忽略了逆变器非线性（死区效应、管压降）
3. 采样延时和PWM更新不同步
这些工程细节，我们后面会一一讲到。
数学模型的建立，是矢量控制的第一步，也是最关键的一步。模型建对了，后面的控制器设计才有根基。模型建错了，后面再怎么调参也是白费力气。
下一节，我们会基于这个dq模型，开始设计电流环和速度环的控制器。到时候你会发现，有了这个数学模型，控制器的设计思路会变得非常清晰。
公众号:蓝海资料掘金营，微信deep3321