2、永磁同步电机数学模型：dq坐标系下的电压方程、磁链方程、转矩方程
好，咱们进入正题。上一章聊了聊MTPA到底是个啥，这一章得把底子打牢。你想想看，要控制好一台永磁同步电机，你总得知道它“长什么样”吧？数学模型就是它的“身份证”。
我个人习惯，一上来先讲dq坐标系。为什么非得搞个dq坐标系？说白了，就是为了把交流电机那套复杂的三相耦合系统，简化成我们熟悉的直流电机控制方式。你想想，三相交流电在空间里转来转去，看着就头疼。但如果我们站在一个和转子一起旋转的坐标系里看，哎，那些交变的量就变成了直流量。这就是Clark变换和Park变换的功劳，这里不展开，我们直接看变换后的结果。
2.1 dq坐标系下的电压方程
这是最核心的公式之一，我建议你把它刻在脑子里。在dq旋转坐标系下，永磁同步电机的定子电压方程长这样：
ud = Rs * id + d(ψd)/dt - ωe * ψq
uq = Rs * iq + d(ψq)/dt + ωe * ψd这里：
- ud, uq：分别是d轴和q轴的定子电压分量。这就是我们要施加的控制量。
- id, iq：d轴和q轴的定子电流分量。这是我们能测量和控制的。
- Rs：定子相电阻。这个值会随温度变化，我在项目中吃过亏，后面会讲。
- ψd, ψq：d轴和q轴的定子磁链。
- ωe：电角速度。注意是电角度，不是机械角度，差了个极对数。
嗯，这里要注意。公式里有两项：一项是电阻压降 Rs * i，另一项是磁链变化产生的感应电动势。其中 -ωe * ψq 和 +ωe * ψd 这两项，就是我们常说的“反电动势”或者“运动电动势”。它跟转速直接相关，转速越高，这个量越大。
核心要点：电压方程告诉我们，要控制电流，本质上是在和反电动势“打架”。尤其是在高速区，反电动势会很大，如果母线电压不够，你就没法把电流“压”进去，这就是弱磁控制的由来。
2.2 dq坐标系下的磁链方程
有了电压方程，还得知道磁链是怎么来的。磁链方程把电流和磁链联系了起来：
ψd = Ld * id + ψf
ψq = Lq * iq这里：
- Ld, Lq：d轴和q轴的电感。对于内置式永磁同步电机（IPMSM），通常
Lq > Ld。这个差异，就是磁阻转矩的来源，也是MTPA能起作用的关键。 - ψf：永磁体产生的磁链。这是一个常数，由电机设计决定。
你看，d轴磁链由两部分组成：一部分是电流 id 产生的，另一部分是永磁体本身贡献的。而q轴磁链，完全由q轴电流产生。为什么会这样？因为我们在建模时，把永磁体的方向定义为了d轴方向。
个人经验：我在调试一个高速主轴电机时，发现id电流稍微给负一点，转速就能往上冲一截。当时不理解，后来才明白，负的id会产生一个与永磁体相反的磁链，相当于“削弱”了气隙磁场，让反电动势降下来，从而允许更高的转速。这就是弱磁控制，其根源就在这个磁链方程里。
2.3 dq坐标系下的转矩方程
终于到了最激动人心的部分——转矩方程。它告诉我们，我们费了半天劲控制的电流，到底能产生多大的力：
Te = 1.5 * pn * (ψd * iq - ψq * id)把磁链方程代入，可以得到更常用的形式：
Te = 1.5 * pn * [ψf * iq + (Ld - Lq) * id * iq]这里 pn 是电机极对数。
这个公式太重要了，我们拆开来看：
- 第一项：
1.5 * pn * ψf * iq—— 这是永磁转矩。它跟q轴电流成正比，跟直流电机的转矩特性很像。你给iq，它就出力。 - 第二项：
1.5 * pn * (Ld - Lq) * id * iq—— 这是磁阻转矩。因为Ld - Lq是负的（对于IPMSM），所以要想让这一项为正，id必须为负。这就是为什么MTPA控制中，我们需要给一个负的id电流。
避坑指南：我曾经遇到过一台电机，怎么调MTPA都感觉效率不对。后来一查，是电机参数 Ld 和 Lq 因为磁路饱和，跟标称值差了很多。你想想看，公式里的 (Ld - Lq) 变了，你算出来的最优电流角肯定就不对了。所以，参数辨识在实际工程中非常重要，千万别迷信数据手册。
2.4 三个方程的内在联系
这三个方程不是孤立的，它们是一个整体。我习惯这样理解：
- 磁链方程是“桥梁”，把电流和磁链联系起来。
- 转矩方程是“目标”，告诉我们最终要什么。
- 电压方程是“约束”，告诉我们能做到什么程度。
举个例子，你想让电机输出大转矩。从转矩方程看，你需要合适的 id 和 iq。但这两个电流一给出去，通过磁链方程，就会产生磁链。磁链再通过电压方程，就会产生反电动势。如果反电动势超过了母线电压能提供的极限，那你的电流就失控了。
所以，MTPA控制，本质上就是在电压极限圆和电流极限圆的约束下，去寻找那个能让转矩/电流比最大的 id 和 iq 组合。这个我们下一章会详细讲。
| 方程类型 | 表达式 | 物理意义 |
|---|---|---|
| 电压方程 | ud = Rs*id + dψd/dt - ωe*ψq uq = Rs*iq + dψq/dt + ωe*ψd |
电压与电流、磁链变化、反电动势的关系 |
| 磁链方程 | ψd = Ld*id + ψf ψq = Lq*iq |
磁链由电流和永磁体共同产生 |
| 转矩方程 | Te = 1.5*pn*[ψf*iq + (Ld-Lq)*id*iq] | 转矩由永磁转矩和磁阻转矩组成 |
重要提醒：以上所有方程，都是建立在“理想电机”假设上的。实际电机存在铁损、磁路饱和、交叉耦合、温度变化等非线性因素。尤其是 Ld 和 Lq，它们不是常数，会随着电流大小和角度变化。在做高精度控制时，一定要考虑这些非理想特性，否则你的MTPA查表法可能完全不准。
好了，这一章的内容就到这里。dq坐标系下的这三个方程，是后续所有控制策略的基石。你花再多时间去理解它们都不过分。下一章，我们就基于这些方程，正式进入MTPA的核心——如何找到那个最优的电流角。