4. 电机模型基础(三):PMSM在dq坐标系下的数学模型

好,咱们继续往下走。前面两章我们把Clark变换和Park变换聊透了,也搞清楚了为什么要把三相静止的玩意儿搬到旋转的坐标系里去看。说白了,就是为了把那些随时间变化的交流量,变成直流量来控制。这一章,我们就来正式推导PMSM在dq坐标系下的数学模型。

我个人觉得,这是整个FOC电流环建模里最核心的一步。你后面写的所有代码,调的所有参数,追根溯源,都跑不出这几个方程。所以,咱们得把它吃透。

核心要点: dq坐标系下的数学模型,是FOC控制的“宪法”。所有控制策略都必须遵守它。

4.1 为什么要用dq坐标系?

你想想看,在三相静止坐标系(ABC)下,电机的电压方程是一组时变系数的微分方程。电感矩阵里全是跟转子位置θ有关的项,解起来非常头疼。控制上更是麻烦,因为你要跟踪的是正弦变化的量。

但在dq坐标系下,我们相当于坐在转子上看世界。对于表贴式PMSM(SPMSM)来说,电感变成了常数;对于内置式PMSM(IPMSM),电感虽然d、q轴不一样,但也是常数。这样一来,时变系统就变成了定常系统,控制难度直接降了一个维度。

我的经验: 我在做第一个FOC项目时,一开始没太在意这个“定常化”的意义,直接在ABC坐标系下硬调PI参数,结果怎么调都调不好,电流波形乱七八糟。后来老老实实回到dq模型,才发现问题出在耦合项上。所以,别嫌推导麻烦,这是值得的。

4.2 dq坐标系下的磁链方程

我们先从磁链说起。磁链是电压和转矩的源头。

在dq坐标系下,定子磁链由两部分组成:一部分是定子电流自己产生的,另一部分是永磁体贡献的。方程如下:

ψd = Ld * id + ψf
ψq = Lq * iq

其中:

  • ψd, ψq:d轴和q轴的定子磁链
  • Ld, Lq:d轴和q轴的电感。对于SPMSM,Ld = Lq;对于IPMSM,Ld < Lq
  • id, iq:d轴和q轴的定子电流
  • ψf:永磁体磁链,这是一个常数(在温度不变的前提下)

这里有个关键点:永磁体只贡献d轴的磁链。为什么?因为我们在做Park变换时,把d轴对准了转子永磁体的N极方向。所以,永磁体的磁场天然就在d轴上。q轴上是没有永磁体磁链的。

注意: 这里的Ld和Lq是“视在电感”或者叫“同步电感”,不是我们测三相绕组时得到的线电感或相电感。它们是通过坐标变换推导出来的等效电感。在Matlab/Simulink里建模时,直接填这两个值就行。

4.3 dq坐标系下的电压方程

有了磁链,我们就可以写电压方程了。dq坐标系下的定子电压方程,是FOC电流环设计的直接依据。

ud = Rs * id + d(ψd)/dt - ωe * ψq
uq = Rs * iq + d(ψq)/dt + ωe * ψd

把磁链方程代入,得到更常用的形式:

ud = Rs * id + Ld * d(id)/dt - ωe * Lq * iq
uq = Rs * iq + Lq * d(iq)/dt + ωe * (Ld * id + ψf)

我们来拆解一下这个方程,看看每一项的物理意义:

项 物理意义 备注
Rs * id / Rs * iq 电阻压降 跟直流电机一样,就是欧姆定律
Ld * d(id)/dt / Lq * d(iq)/dt 电感压降(变压器电动势) 电流变化引起的,是我们要控制的对象
-ωe * Lq * iq d轴反电动势(运动电动势) 耦合项! 由q轴电流和转速引起
ωe * Ld * id q轴反电动势(运动电动势) 耦合项! 由d轴电流和转速引起
ωe * ψf q轴永磁体反电动势 这就是我们常说的“反电动势”,跟转速成正比

看到没?d轴方程里出现了q轴电流,q轴方程里出现了d轴电流。这就是交叉耦合。转速越高,耦合越强。如果不做解耦控制,你调d轴的PI,q轴会跟着抖;调q轴的PI,d轴也会受影响。这就是为什么很多初学者发现电流环调不稳的原因。

我曾经在一个高速主轴电机项目上吃过这个亏。电机额定转速12000rpm,我直接用基本的PI控制,结果电流环带宽死活上不去,一加速就震荡。后来加了前馈解耦,把 -ωe * Lq * iq 和 ωe * (Ld * id + ψf) 这两项补偿掉,问题才解决。解耦的具体方法,我们后面章节会详细讲。

4.4 dq坐标系下的转矩方程

转矩方程是电机控制的目标。我们费这么大劲控制电流,最终是为了控制转矩。

在dq坐标系下,电磁转矩的通用表达式为:

Te = 1.5 * p * [ ψf * iq + (Ld - Lq) * id * iq ]

其中 p 是电机极对数。

这个方程包含两项:

  • ψf * iq:永磁转矩。这是主要的转矩成分,由q轴电流和永磁体相互作用产生。对于所有PMSM,这一项都存在。
  • (Ld - Lq) * id * iq:磁阻转矩。这一项只有在Ld ≠ Lq时才存在,也就是内置式PMSM(IPMSM)的专利。因为IPMSM的转子磁路不对称,d轴磁阻大(电感小),q轴磁阻小(电感大),所以Ld < Lq。利用这个差异,我们可以通过注入负的id(弱磁)来产生额外的转矩。

重要结论:

  • 对于表贴式PMSM (SPMSM):Ld = Lq,磁阻转矩为0。转矩只与iq成正比。控制策略就是让id=0,只控制iq。
  • 对于内置式PMSM (IPMSM):Ld < Lq,磁阻转矩不为0。我们可以利用这个特性,通过控制id和iq的分配,在相同电流幅值下获得更大的转矩,这就是最大转矩电流比(MTPA)控制。

4.5 运动方程

最后,我们还需要一个方程把电磁转矩和电机的机械运动联系起来。这就是运动方程:

Te - Tl = J * d(ωm)/dt + B * ωm

其中:

  • Te:电磁转矩(我们算出来的)
  • Tl:负载转矩(外部施加的)
  • J:转动惯量(电机转子+负载折算过来的)
  • ωm:机械角速度(rad/s)。注意,ωe = p * ωm
  • B:粘滞摩擦系数

这个方程告诉我们:电磁转矩减去负载转矩,剩下的部分用来加速(或减速)电机,以及克服摩擦。在速度环和位置环的设计中,这个方程是核心。

4.6 知识体系总览

说了这么多,我们来画一张图,把这一章的知识点串起来。这样你脑子里能有个整体的框架。

PMSM dq坐标系数学模型知识体系 PMSM dq坐标系数学模型 坐标变换 (Clark + Park) 磁链方程 ψd = Ld·id + ψf ψq = Lq·iq 电压方程 ud = Rs·id + Ld·d(id)/dt - ωe·Lq·iq uq = Rs·iq + Lq·d(iq)/dt + ωe·(Ld·id+ψf) 转矩方程 Te = 1.5·p·[ψf·iq + (Ld-Lq)·id·iq] 永磁转矩 + 磁阻转矩 运动方程 Te - Tl = J·d(ωm)/dt + B·ωm 电磁转矩 - 负载转矩 = 加速 + 摩擦 应用:FOC电流环设计基础 电压方程 → 电流环PI参数设计 | 转矩方程 → 控制策略选择 (id=0 / MTPA / 弱磁)

这张图把整个逻辑串起来了:通过坐标变换,我们得到了dq坐标系下的四个核心方程。磁链方程是基础,电压方程是控制对象,转矩方程是控制目标,运动方程连接了电机和机械负载。后面的电流环设计,就是围绕着电压方程和转矩方程展开的。

4.7 本章小结

嗯,这一章内容不少,但都是硬核干货。我们来快速回顾一下:

  • 磁链方程:定义了d、q轴磁链与电流、永磁体的关系。永磁体只贡献d轴磁链。
  • 电压方程:包含了电阻压降、电感压降和反电动势项。注意交叉耦合项的存在,这是后续解耦控制的关键。
  • 转矩方程:包含永磁转矩和磁阻转矩。SPMSM只有永磁转矩,IPMSM两者都有。
  • 运动方程:连接了电磁转矩和机械运动,是速度环设计的基础。

这些方程,你最好能手推一遍。别光看,拿笔在纸上写一写。我当年学的时候,光这个推导就写了满满三页纸。写完之后,你对FOC的理解会上一个台阶。