2、永磁同步电机(PMSM)数学模型:电压方程、磁链方程、转矩方程在dq坐标系下的推导
各位工程师朋友,大家好。今天我们聊一个绕不开的话题——PMSM的数学模型。
说实话,我刚入行那会儿,看到dq变换那一堆公式就头疼。心想:我直接搞三相不行吗?后来在项目里被现实狠狠教育了一回。嗯,今天我就把当年踩过的坑,连同这些方程怎么来的,一次性讲清楚。
2.1 为什么非要搞dq坐标系?
三相静止坐标系下的PMSM方程,你去看一眼就知道——时变系数、强耦合、非线性。说白了,就是一团乱麻。
我举个例子。你在三相坐标系下写电压方程,电感矩阵里全是转子位置θ的函数。每算一步,θ都在变。你想想看,这怎么设计控制器?
dq变换的核心思想,就是跟着转子转。你站在转子上看定子,所有电感都变成了常数。这就是所谓的「去时变、解耦合」。
核心结论:dq坐标系下,PMSM的数学模型从时变系统变成了定常系统。这是矢量控制的理论基石。
2.2 磁链方程:从三相到dq
我们先从磁链说起。三相静止坐标系下,定子磁链长这样:
ψa = Laa·ia + Lab·ib + Lac·ic + ψf·cos(θ)
ψb = Lba·ia + Lbb·ib + Lbc·ic + ψf·cos(θ - 120°)
ψc = Lca·ia + Lcb·ib + Lcc·ic + ψf·cos(θ + 120°)
你看,自感Laa、互感Lab这些,全都是θ的函数。比如Laa = Ls0 + Ls2·cos(2θ)。这玩意儿你直接拿去用,控制器根本没法写。
经过Park变换(也就是dq变换)之后,磁链方程变得极其清爽:
ψd = Ld·id + ψf
ψq = Lq·iq
这里Ld和Lq是常数。ψf是永磁体磁链,也是常数。
我个人习惯:在写代码时,先把ψf标定准确。这个值如果偏了,后面所有转矩估算都会跟着偏。我曾经在一个低速大扭矩项目里,就因为ψf标定差了5%,导致转矩输出一直对不上。
为什么Ld和Lq不一样?因为转子磁路结构不对称。内嵌式PMSM(IPMSM)的Lq通常大于Ld。表贴式(SPMSM)则Ld ≈ Lq。
2.3 电压方程:动态过程的描述
有了磁链,电压方程就好办了。三相静止坐标系下的电压方程我就不写了,太啰嗦。直接看dq坐标系下的结果:
ud = Rs·id + dψd/dt - ωe·ψq
uq = Rs·iq + dψq/dt + ωe·ψd
这里面有几个关键项:
- Rs·id、Rs·iq:电阻压降,低速时占比大
- dψd/dt、dψq/dt:变压器电动势,动态过程才有
- -ωe·ψq、+ωe·ψd:旋转电动势,也叫反电动势。这是速度相关的项
把磁链方程代进去,得到更常用的形式:
ud = Rs·id + Ld·did/dt - ωe·Lq·iq
uq = Rs·iq + Lq·diq/dt + ωe·(Ld·id + ψf)
注意:这里的ωe是电角速度,不是机械角速度。两者差一个极对数:ωe = p·ωm。我见过有新手直接拿机械转速往里代,结果算出来的反电动势全错了。
为什么会这样?因为磁链旋转一圈对应一个电周期,而转子转一圈对应p个电周期。这个关系搞错了,整个模型就废了。
2.4 转矩方程:力从哪来?
转矩方程是大家最关心的。毕竟我们做电机控制,最终要的就是转矩。
从能量角度看,电磁转矩等于磁场储能对机械角度的偏导。推导过程我就不展开了,直接给结论:
Te = 1.5·p·[ψf·iq + (Ld - Lq)·id·iq]
这个公式里有两部分:
| 项 | 含义 | 说明 |
|---|---|---|
| 1.5·p·ψf·iq | 永磁转矩 | 主转矩分量,与iq成正比 |
| 1.5·p·(Ld-Lq)·id·iq | 磁阻转矩 | 只有IPMSM才有,SPMSM此项为0 |
对于SPMSM,Ld = Lq,磁阻转矩项消失。转矩只由iq控制,简单粗暴。
对于IPMSM,磁阻转矩可以利用起来。我做过一个项目,在高速区用负id来产生磁阻转矩,同样的电流下转矩能多出15%-20%。这就是所谓的最大转矩电流比(MTPA)控制。
避坑指南:我曾经在调试一个低速大扭矩无传感器算法时,发现转矩估算值总是偏小。查了两天才发现,是Ld和Lq的参数给反了。IPMSM的Lq > Ld,这个顺序千万别搞错。
2.5 知识体系总览
下面这张图,我把整个dq坐标系下PMSM数学模型的核心逻辑画出来了。你看一眼就能明白各个方程之间的关系:
2.6 工程中的实用要点
理论讲完了,说几个实际干活时要注意的点:
- 参数敏感性:Ld、Lq、ψf这三个参数,在实际运行中会随温度和电流变化。我建议在代码里做在线参数辨识,或者至少做查表补偿。
- 低速时的电压方程:转速很低时,反电动势项ωe·ψq和ωe·ψd几乎为零。这时候电压方程退化成RL电路。无传感器控制在低速区难做,根本原因就在这里——反电动势信号太弱了。
- 转矩估算的精度:如果你用公式Te = 1.5·p·[ψf·iq + (Ld-Lq)·id·iq]来估算转矩,务必确保id和iq的采样同步。我吃过这个亏——电流采样有延时,算出来的转矩一直在抖。
一个小技巧:在调试阶段,我习惯把估算的转矩和实测的转矩(用扭矩传感器)同时录下来,画在同一张图上。偏差超过5%就要查参数了。这个方法帮我快速定位过好几次参数标定错误。
好了,关于PMSM在dq坐标系下的数学模型,今天就讲到这里。这三个方程——磁链方程、电压方程、转矩方程——是后续所有控制算法的基础。你把它吃透了,后面学无传感器控制、MTPA、弱磁控制,都会轻松很多。
公众号:蓝海资料掘金营,微信deep3321