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数学模型的核心逻辑画出来了。你看一眼就能明白各个方程之间的关系:

PMSM dq坐标系数学模型核心逻辑 三相静止坐标系 时变、耦合、非线性 Park变换 dq旋转坐标系 定常、解耦、线性 磁链方程 ψd = Ld·id + ψf ψq = Lq·iq 电压方程 ud = Rs·id + dψd/dt - ωe·ψq uq = Rs·iq + dψq/dt + ωe·ψd 转矩方程 Te = 1.5·p·[ψf·iq + (Ld-Lq)·id·iq] 应用:矢量控制 / 无传感器控制 / MTPA 基于上述三个方程进行控制器设计和观测器设计 注:SPMSM中Ld=Lq,磁阻转矩项为0;IPMSM中Lq>Ld,可利用磁阻转矩

2.6 工程中的实用要点

理论讲完了,说几个实际干活时要注意的点:

  1. 参数敏感性:Ld、Lq、ψf这三个参数,在实际运行中会随温度和电流变化。我建议在代码里做在线参数辨识,或者至少做查表补偿。
  2. 低速时的电压方程:转速很低时,反电动势项ωe·ψq和ωe·ψd几乎为零。这时候电压方程退化成RL电路。无传感器控制在低速区难做,根本原因就在这里——反电动势信号太弱了。
  3. 转矩估算的精度:如果你用公式Te = 1.5·p·[ψf·iq + (Ld-Lq)·id·iq]来估算转矩,务必确保id和iq的采样同步。我吃过这个亏——电流采样有延时,算出来的转矩一直在抖。

一个小技巧:在调试阶段,我习惯把估算的转矩和实测的转矩(用扭矩传感器)同时录下来,画在同一张图上。偏差超过5%就要查参数了。这个方法帮我快速定位过好几次参数标定错误。

好了,关于PMSM在dq坐标系下的数学模型,今天就讲到这里。这三个方程——磁链方程、电压方程、转矩方程——是后续所有控制算法的基础。你把它吃透了,后面学无传感器控制、MTPA、弱磁控制,都会轻松很多。


公众号:蓝海资料掘金营,微信deep3321