3. 电机数学模型：永磁同步电机（PMSM）在 dq 坐标系下的数学模型推导与参数辨识
各位同学，咱们今天来啃一块硬骨头——PMSM的数学模型。说实话，我刚入行那会儿，看到满屏的公式也头疼。但后来我发现，搞懂这个模型，就像拿到了电机控制的“通关密码”。你想想看，没有模型，你连PID调参都只能靠瞎蒙。
我个人习惯，先把物理概念讲清楚，再上公式。这样你心里有底，不会觉得数学是空中楼阁。
3.1 为什么非要搞个 dq 坐标系？
三相电机的物理量（电压、电流、磁链）都是随时间变化的交流量。在ABC三相静止坐标系下，这些量互相耦合，控制起来非常麻烦。说白了，你调A相，B相和C相也跟着变，像个“三头怪”。
我做过一个项目，一开始直接在ABC坐标系下写控制算法，结果调试了两个月，电流环带宽死活上不去。后来改成dq坐标系，一周就搞定了。这就是数学变换的力量。
dq变换的核心思想：把交流量变成直流量。这样，我们就可以像控制直流电机一样控制PMSM了。
- d轴（直轴）：与转子磁极轴线重合，负责产生磁通。
- q轴（交轴）：与d轴垂直，负责产生转矩。
嗯，这里要注意：dq坐标系是旋转坐标系，它跟着转子一起转。所以，在dq轴上看，电机的电压、电流都是恒定值（稳态时）。
3.2 PMSM 在 dq 坐标系下的电压方程
好，咱们直接上干货。PMSM在dq坐标系下的电压方程是这样的：
ud = Rs * id + Ld * (did/dt) - ωe * Lq * iq
uq = Rs * iq + Lq * (diq/dt) + ωe * (Ld * id + ψf)
别急着背，咱们拆开来看：
- ud, uq：d轴和q轴的电压。这就是我们要施加的控制量。
- id, iq：d轴和q轴的电流。这是我们要控制的变量。
- Rs：定子电阻。这个值会随着温度变化，我在项目中吃过它的亏，后面讲参数辨识时会细说。
- Ld, Lq：d轴和q轴的电感。对于表贴式PMSM（SPMSM），Ld ≈ Lq；对于内置式PMSM（IPMSM），Ld < Lq。
- ωe：电角速度。ωe = 极对数 * 机械角速度。
- ψf：永磁体磁链。这是电机的“灵魂参数”，决定了反电动势的大小。
关键点：方程中的 -ωe * Lq * iq 和 +ωe * (Ld * id + ψf) 是交叉耦合项。正是这两项，让dq轴之间产生了耦合。高速时，耦合效应非常明显，如果不做解耦控制，电流环会振荡甚至发散。
3.3 转矩方程与运动方程
有了电压方程，我们还得知道电机能出多大的力。转矩方程如下：
Te = 1.5 * pn * [ψf * iq + (Ld - Lq) * id * iq]
其中 pn 是极对数。这个公式很有意思：
- 第一项：
ψf * iq，这是永磁转矩。说白了，就是永磁体和电流产生的磁场相互作用产生的转矩。 - 第二项：
(Ld - Lq) * id * iq，这是磁阻转矩。对于IPMSM，Ld < Lq，所以这一项是正的。我们可以利用它来“榨取”更多的转矩。
我曾经在一个IPMSM项目中，通过注入负的id电流（弱磁控制），让电机在高速区还能输出不错的转矩。这就是利用了磁阻转矩。
运动方程就简单了：
Te - TL = J * (dωm/dt) + B * ωm
- TL：负载转矩。
- J：转动惯量。
- B：阻尼系数。
- ωm：机械角速度。
3.4 参数辨识——从“纸上谈兵”到“实战”
模型有了，参数从哪来？厂家给的 datasheet 往往不准，或者干脆没有。这时候，就得靠我们自己动手辨识了。
我的经验：参数辨识是仿真和实物的“分水岭”。仿真时参数随便设，实物上参数不准，控制效果一塌糊涂。我建议，拿到一个新电机，先花半天时间把参数摸清楚。
3.4.1 定子电阻 Rs 的辨识
最简单的方法：直流伏安法。
- 给电机施加一个直流电压（注意电流不要超过额定值）。
- 测量稳态时的电压和电流。
- Rs = U / I。
我曾经用这个方法测一个50W的小电机，结果发现电阻值比datasheet大了30%。后来一查，是电机内部接线端子氧化了。所以，实测永远比看资料靠谱。
3.4.2 dq 轴电感 Ld、Lq 的辨识
这里需要一点技巧。我常用的方法是高频注入法：
- 让电机转子固定不动（或者用机械锁住）。
- 在d轴方向注入一个高频电压信号（比如500Hz，幅值5%额定电压）。
- 测量d轴的高频电流响应。
- 根据
Ld = U_high_freq / (2 * π * f_high * I_high_freq)计算。 - 同样的方法，在q轴注入，得到Lq。
注意：高频注入时，电流幅值不能太大，否则会导致电机发热甚至退磁。我一般控制在额定电流的10%以内。另外，注入频率要远高于电机的电气时间常数，但又要低于PWM开关频率的1/10。
3.4.3 永磁体磁链 ψf 的辨识
这个参数可以通过空载反电动势法得到：
- 用外力（或者让电机自由旋转）拖动电机，使其以恒定速度旋转。
- 测量电机的线电压（反电动势）。
- 根据公式
ψf = E_line / (√3 * ωe)计算。
嗯，这里要注意：反电动势的波形应该是正弦波。如果波形畸变严重，说明电机可能有匝间短路或者退磁问题。我遇到过一台电机，反电动势波形像“锯齿”，一查是磁钢碎了。
3.5 参数辨识的实用表格
为了方便大家实操，我整理了一个参数辨识的对照表：
| 参数 | 辨识方法 | 所需设备 | 注意事项 |
|---|---|---|---|
| Rs | 直流伏安法 | 直流电源、万用表 | 注意温度影响，建议在室温下测量 |
| Ld, Lq | 高频注入法 | 示波器、信号发生器 | 注入频率和幅值要合适，避免饱和 |
| ψf | 空载反电动势法 | 示波器、转速计 | 确保转速稳定，波形无畸变 |
| J | 加减速法 | 编码器、电流传感器 | 需要知道负载转矩，或者空载测试 |
3.6 小结与避坑指南
咱们今天把PMSM的数学模型和参数辨识方法捋了一遍。总结一下：
- dq变换是核心，把交流变直流，简化控制。
- 电压方程中的交叉耦合项是高速控制的难点，需要做解耦。
- 转矩方程揭示了永磁转矩和磁阻转矩的贡献，IPMSM可以利用磁阻转矩。
- 参数辨识是连接仿真和实物的桥梁，一定要重视。
避坑指南：我曾经在一个项目中，直接用了datasheet上的电感值做仿真，结果仿真效果完美，一上实物电流环就啸叫。后来实测发现，实际电感值只有datasheet标称值的60%。原因是datasheet是在小信号下测的，而实际工况下电流很大，电感已经饱和了。所以，参数辨识一定要在接近实际工况的条件下进行。
下一章，咱们会基于今天建立的数学模型，开始搭建仿真模型，并手把手调参。到时候你会发现，今天啃的这块硬骨头，会变成你最趁手的工具。