第二章：电机数学模型基础——永磁同步电机(PMSM)的数学模型、电阻电感反电动势系数、转动惯量与阻尼系数
大家好，欢迎来到《Simulink电机参数辨识与标定实战》的第二讲。
今天咱们聊聊永磁同步电机（PMSM）的数学模型。说实话，很多做电机控制的工程师，一上来就调PI参数，模型这块能跳过就跳过。我个人觉得，这其实是个大坑。模型是控制算法的基础，你连电机怎么转的都不清楚，怎么调得好参数？
好，咱们正式开始。
2.1 永磁同步电机的物理模型
先看个最简单的结构。PMSM的定子上有三相绕组，转子上有永磁体。通电后，定子产生旋转磁场，拉着转子跟着转。就这么简单。
但你要做控制，就得用数学语言描述这个过程。嗯，这里要注意，我们通常做两个假设：
- 忽略铁芯饱和——说白了就是认为磁路是线性的
- 忽略涡流和磁滞损耗——这些在高频下才明显，基频控制可以不管
我在项目中遇到过一个问题：有个同事死活调不好电流环，后来发现是模型里没考虑饱和效应。当然，那是大功率电机的情况。一般小功率电机，这两个假设完全够用。
2.2 三相静止坐标系下的数学模型
三相坐标系下的方程长这样：
电压方程：
[ua] [R 0 0 ] [ia] [Laa Lab Lac] d [ia] [ea]
[ub] = [0 R 0 ] [ib] + [Lba Lbb Lbc] — [ib] + [eb]
[uc] [0 0 R ] [ic] [Lca Lcb Lcc] dt[ic] [ec]
其中：
ea = -ψf * ωe * sin(θe)
eb = -ψf * ωe * sin(θe - 2π/3)
ec = -ψf * ωe * sin(θe + 2π/3)看着是不是有点晕？别急，这玩意儿在实际控制中基本不用。为什么？因为三相坐标系下，电感矩阵是时变的，跟转子位置有关。你想想看，每算一步都要更新电感矩阵，计算量太大了。
核心要点：三相模型主要用于仿真验证，实际控制算法都用dq坐标系下的模型。
2.3 dq坐标系下的数学模型——这才是主角
通过Clark变换和Park变换，我们把三相静止坐标系转换到两相旋转坐标系。说白了，就是把交流量变成直流量，这样控制起来就方便多了。
dq坐标系下的电压方程：
ud = R*id + Ld*(did/dt) - ωe*Lq*iq
uq = R*iq + Lq*(diq/dt) + ωe*(Ld*id + ψf)这里：
- ud, uq —— d轴和q轴电压
- id, iq —— d轴和q轴电流
- Ld, Lq —— d轴和q轴电感
- R —— 定子电阻
- ψf —— 永磁体磁链
- ωe —— 电角速度
转矩方程：
Te = 1.5 * p * [ψf*iq + (Ld - Lq)*id*iq]其中p是极对数。注意看，转矩由两部分组成：
- 永磁转矩：1.5 * p * ψf * iq —— 这是主要的
- 磁阻转矩：1.5 * p * (Ld - Lq) * id * iq —— 只有Ld≠Lq时才存在
我记得刚开始做表贴式PMSM（SPMSM）时，Ld≈Lq，磁阻转矩基本为零。后来做内置式（IPMSM），才发现磁阻转矩能贡献不少扭矩。这个在MTPA控制里特别重要。
个人经验：做参数辨识时，我建议先辨识电阻R，再辨识电感Ld和Lq。因为电阻受温度影响大，而电感相对稳定。分开辨识，收敛更快。
2.4 反电动势系数
反电动势系数，也叫Ke，单位是V/(rad/s)或者V/krpm。它跟永磁体磁链ψf的关系是：
Ke = ψf * p或者用线反电动势表示：
E_line = Ke * ωm其中ωm是机械角速度。
怎么测？很简单。把电机拖到某个转速，用示波器看线反电动势的峰值，然后算一下就行。我曾经在项目里遇到过一个问题：测出来的Ke跟标称值差了15%。后来发现是温度没控制好——永磁体的磁通随温度升高会下降，大约每升高10°C，磁通下降0.5%~1%。
避坑指南：我曾经在高温环境下做标定，Ke值偏低，导致电流环增益计算错误。后来我养成了习惯：每次标定前先让电机冷却到室温，或者记录当前温度做补偿。
2.5 转动惯量与阻尼系数
这两个参数主要影响机械动态。运动方程：
Te - TL = J * (dωm/dt) + B * ωm其中：
- J —— 转动惯量，单位kg·m²
- B —— 阻尼系数，单位N·m·s/rad
- TL —— 负载转矩
转动惯量J决定了电机加速的快慢。J越大，加速越慢。阻尼系数B主要是摩擦和风阻，一般很小，但在低速时影响明显。
怎么辨识？我常用的方法是：
- 给电机一个阶跃转矩指令
- 记录转速响应曲线
- 用最小二乘法拟合出J和B
具体来说，如果电机空载（TL=0），转速响应是一个一阶惯性环节：
ωm(s) / Te(s) = 1 / (J*s + B)从阶跃响应中，你可以读出稳态转速和时间常数，然后反推J和B。
关键参数表：
| 参数 | 符号 | 单位 | 典型范围（小功率） |
|---|---|---|---|
| 定子电阻 | R | Ω | 0.1 ~ 10 |
| d轴电感 | Ld | mH | 0.1 ~ 10 |
| q轴电感 | Lq | mH | 0.1 ~ 10 |
| 永磁体磁链 | ψf | Wb | 0.01 ~ 0.5 |
| 反电动势系数 | Ke | V/(rad/s) | 0.01 ~ 0.5 |
| 转动惯量 | J | kg·m² | 1e-5 ~ 1e-3 |
| 阻尼系数 | B | N·m·s/rad | 1e-5 ~ 1e-3 |
2.6 小结与思考
好了，这一章的内容就这些。总结一下：
- PMSM的数学模型，核心就是dq坐标系下的电压方程和转矩方程
- 反电动势系数Ke跟磁链ψf直接相关，受温度影响大
- 转动惯量J和阻尼系数B，决定了电机的机械动态特性
下一章，我会带大家用Simulink搭建PMSM的仿真模型，然后做参数辨识。到时候咱们手把手操作，看看这些参数在仿真里是怎么跑出来的。
最后留个思考题：为什么dq坐标系下的电感Ld和Lq通常不相等？这个问题想明白了，你对IPMSM的理解就上了一个台阶。
咱们下章见。