3. PMSM在dq坐标系下的数学模型：电压方程、磁链方程、转矩方程与运动方程
好，咱们进入正题。上一章我们把三相静止坐标系下的PMSM模型推了一遍，说实话，那个模型看着就头疼——三个相、互相耦合、还随时间变化。你想想看，要是拿那个模型去写控制算法，光是解耦就能把你绕晕。
所以业界通用的做法，就是做坐标变换。把三相静止的abc坐标系，先变到两相静止的αβ坐标系（Clark变换），再变到两相旋转的dq坐标系（Park变换）。这一章，我们就直接站在dq坐标系下，把PMSM的四个核心方程讲透。
核心观点： dq坐标系下的数学模型，是FOC（磁场定向控制）的基石。没有这个模型，你后面的电流环、速度环全都无从谈起。我个人习惯，在调任何电机之前，先把这四个方程写在白板上，时刻提醒自己：我在控制什么。
3.1 为什么非要用dq坐标系？
说白了，就是为了把交流电机“伪装”成直流电机来控。在abc坐标系下，定子电流产生的磁动势是旋转的，而且三相之间互相影响。但在dq坐标系下，我们把坐标系固定在转子上，跟着转子一起转。
这样一来：
- d轴（直轴）： 与转子永磁体磁场方向对齐。控制d轴电流，就是控制励磁分量。
- q轴（交轴）： 与转子永磁体磁场方向垂直。控制q轴电流，就是控制转矩分量。
你看，是不是清晰多了？d轴管“磁场”，q轴管“转矩”，互不干扰。这就是FOC能实现高性能控制的核心逻辑。
我的经验： 刚接触FOC时，我总搞不清d轴和q轴到底谁管谁。后来在调试一个伺服项目时，发现电机空载电流很大，查了半天才发现是d轴电流给大了，相当于一直在给电机“充磁”却不干活。从那以后，我记住了：d轴是“浪费”的轴，q轴才是“干活”的轴。
3.2 电压方程
先上公式，别怕，我会一句句拆开讲。
ud = Rs * id + Ld * (did/dt) - ωe * Lq * iq
uq = Rs * iq + Lq * (diq/dt) + ωe * (Ld * id + ψf)其中：
- ud, uq： d轴和q轴的定子电压（V）
- id, iq： d轴和q轴的定子电流（A）
- Rs： 定子相电阻（Ω）
- Ld, Lq： d轴和q轴的电感（H）—— 对于表贴式PMSM，Ld ≈ Lq；对于内置式PMSM，Lq > Ld
- ωe： 电角速度（rad/s），ωe = 极对数 * 机械角速度
- ψf： 永磁体磁链（Wb）
你仔细看，这个方程其实可以拆成三部分：
- 电阻压降： Rs * id 和 Rs * iq —— 就是欧姆定律，没啥好说的。
- 电感压降（暂态项）： Ld * (did/dt) 和 Lq * (diq/dt) —— 电流变化时产生的反电动势。
- 运动反电动势（耦合项）： -ωe * Lq * iq 和 +ωe * (Ld * id + ψf) —— 这是dq轴之间的交叉耦合项，也是FOC里需要做前馈补偿的关键。
注意： 我曾经在调试一个高速电机时，忽略了反电动势耦合项，结果电流环怎么调都震荡。后来加上前馈补偿，问题立刻解决。所以，千万别小看那个 -ωe * Lq * iq 和 +ωe * (Ld * id + ψf)，它们在高转速下是主要矛盾。
3.3 磁链方程
磁链方程相对简单，它是电压方程的基础。
ψd = Ld * id + ψf
ψq = Lq * iq其中：
- ψd, ψq： d轴和q轴的定子磁链（Wb）
你看，d轴磁链由两部分组成：一部分是id电流产生的（Ld * id），另一部分是永磁体自带的（ψf）。而q轴磁链完全由iq电流产生（Lq * iq）。
这个方程告诉我们一个道理：
- 如果你想改变d轴磁链，你可以调节id，但永磁体那部分你动不了。
- 如果你想改变q轴磁链，你只能调节iq。
这也是为什么在弱磁控制时，我们给负的id——就是为了抵消一部分永磁体磁链，从而降低反电动势，让电机能在更高转速下运行。
3.4 转矩方程
这是大家最关心的方程，毕竟电机最终是要出力气的。
Te = 1.5 * p * (ψd * iq - ψq * id)把磁链方程代入，可以得到更常用的形式：
Te = 1.5 * p * [ψf * iq + (Ld - Lq) * id * iq]其中 p 是极对数。
这个方程可以拆成两项：
- 永磁转矩： 1.5 * p * ψf * iq —— 这是主转矩，由永磁体和q轴电流相互作用产生。你给iq，它就出力。
- 磁阻转矩： 1.5 * p * (Ld - Lq) * id * iq —— 这是由d轴和q轴电感差异产生的。对于表贴式PMSM（Ld = Lq），这一项为0。对于内置式PMSM（Lq > Ld），这一项不为0，而且可以通过给负的id来利用它，获得更大的总转矩。
关键点： 对于内置式PMSM（IPMSM），磁阻转矩是“免费”的额外转矩。我在做电动汽车驱动项目时，就充分利用了磁阻转矩，在同样的电流下多输出了15%的转矩。这就是所谓的MTPA（最大转矩电流比）控制。
3.5 运动方程
最后一个方程，描述的是电机的机械动态。
Te - TL = J * (dωm/dt) + B * ωm其中：
- Te： 电磁转矩（Nm）
- TL： 负载转矩（Nm）
- J： 转动惯量（kg·m²）
- ωm： 机械角速度（rad/s）
- B： 阻尼系数（Nm·s/rad）
这个方程说白了就是牛顿第二定律的旋转版本：合力矩 = 转动惯量 × 角加速度。Te是驱动力矩，TL是阻力矩，差值决定了电机是加速还是减速。
在实际项目中，转动惯量J是个很头疼的参数。负载变化时，J会变，导致速度环的PI参数需要重新整定。我建议你在做速度环设计时，留出足够的裕量，或者干脆用自适应控制。
3.6 知识体系总览
下面这张图，把四个方程之间的关系画清楚了。你看一遍，应该就能建立起整体认知。
3.7 四个方程之间的逻辑链条
从上图你能看出，这四个方程不是孤立的，它们构成了一个闭环：
- 电压方程 告诉你：给定ud、uq，会产生什么样的id、iq。
- 磁链方程 告诉你：id、iq会产生什么样的ψd、ψq。
- 转矩方程 告诉你：ψd、ψq和id、iq共同作用，产生多大的电磁转矩Te。
- 运动方程 告诉你：Te减去负载TL，决定了电机是加速还是减速，最终输出机械角速度ωm。
而ωm又通过电角速度ωe反馈到电压方程中，形成完整的机电耦合。
我的建议： 在写FOC代码之前，先在纸上把这四个方程的逻辑链条画一遍。我每次接手一个新项目，都会做这件事。它能帮你快速定位问题——比如电机转不起来，是电压不够？还是转矩不够？还是负载太大？
3.8 参数表：表贴式 vs 内置式 PMSM
最后，我用一个表格总结一下两种常见PMSM在dq模型下的参数差异：
| 参数 | 表贴式 PMSM (SPMSM) | 内置式 PMSM (IPMSM) |
|---|---|---|
| Ld 与 Lq 关系 | Ld ≈ Lq | Lq > Ld（通常Lq是Ld的2~3倍） |
| 磁阻转矩 | 几乎为0 | 显著，可利用MTPA |
| 弱磁能力 | 弱（Ld小，弱磁效果有限） | 强（Ld大，弱磁范围宽） |
| 转矩方程简化 | Te = 1.5 * p * ψf * iq | Te = 1.5 * p * [ψf * iq + (Ld-Lq)*id*iq] |
| 典型应用 | 伺服、机器人、低速直驱 | 电动汽车、高速主轴、压缩机 |