2. 永磁同步电机数学模型:从坐标变换到转矩方程

各位同学好,欢迎来到第二章。这一章我们要啃的,是永磁同步电机(PMSM)控制的理论基石——数学模型。说实话,我刚入行那会儿,觉得数学公式又长又枯燥,心想“我直接调参不就行了?”结果项目一上高速弱磁,电机直接抖得像筛子……从那以后,我再也不敢跳过数学模型了。

这一章我会带着大家,从最直观的三相物理模型出发,一步步推导到我们控制算法里真正用的 dq 轴模型。你跟着走一遍,后面搭建 Simulink 模型时会发现,原来每个模块都是有根有据的。

2.1 为什么需要坐标变换?

先问一个问题:PMSM 定子上有三相绕组,通入三相正弦电流,产生旋转磁场。这个磁场在空间上是旋转的,而转子也在转。如果我们直接在 abc 三相静止坐标系下写方程,你会发现——电感是随转子位置变化的,方程里全是时变系数,求解非常麻烦。

说白了,我们做控制的,希望看到一个“看起来像直流电机”的模型。直流电机的电枢磁场和励磁磁场是解耦的,控制起来多简单!所以,坐标变换的核心思想就是:把交流电机等效成直流电机来控。

核心思路: 三相静止 abc → 两相静止 αβ(Clark 变换)→ 两相旋转 dq(Park 变换)。最终在 dq 坐标系下,电感变成常数,磁链和转矩可以独立控制。

2.2 Clark 变换:从 abc 到 αβ

Clark 变换,也叫 3/2 变换。它的任务是把三相静止坐标系下的变量,投影到两相静止坐标系上。你可以想象成:三个轴(相差 120°)投影到两个垂直的轴(α 轴和 β 轴)上。

数学上,以电流为例,变换公式如下:

i_α = i_a
i_β = (1/√3) * i_a + (2/√3) * i_b

等一下,这里有个细节。上面这个是等幅值变换,也是我个人最常用的形式。还有一种叫等功率变换,系数会差一个 √(2/3)。

变换类型 变换矩阵系数 特点
等幅值变换 2/3 变换前后电流幅值不变,方便调试
等功率变换 √(2/3) 变换前后功率不变,理论推导常用

我的经验: 在 Simulink 里搭建模型时,我习惯用等幅值变换。因为调试时可以直接用示波器看电流波形,幅值对得上,心里踏实。等功率变换在推导转矩公式时更简洁,但实际工程中,等幅值更直观。

Clark 变换的逆变换也很重要,从 αβ 回到 abc:

i_a = i_α
i_b = -0.5 * i_α + (√3/2) * i_β
i_c = -0.5 * i_α - (√3/2) * i_β

嗯,这里要注意,逆变换的系数是 1,不是 2/3。因为正变换已经缩放过一次了。

2.3 Park 变换:从 αβ 到 dq

Clark 变换之后,我们得到了两相静止的 αβ 坐标系。但 αβ 坐标系下的电压、电流仍然是正弦波,控制起来还是麻烦。Park 变换就是要把这些正弦波变成直流量。

怎么做?让坐标系跟着转子转!

d 轴(直轴)对准转子永磁体的 N 极方向,q 轴(交轴)超前 d 轴 90° 电角度。这样,在 dq 旋转坐标系下看,电流和电压都变成了直流分量。

变换公式:

i_d =  i_α * cos(θ) + i_β * sin(θ)
i_q = -i_α * sin(θ) + i_β * cos(θ)

其中 θ 是转子位置的电角度。这个 θ 从哪里来?要么用编码器测,要么用无位置传感器算法估。我们这门课后面会专门讲。

我曾经踩过的坑: Park 变换的角度定义,不同文献可能不一样。有的把 d 轴定义在 α 轴方向上,有的定义在 β 轴方向上。一定要确认你的角度定义和变换矩阵是匹配的,否则电流环会反着转!我建议你在 Simulink 里先给一个直流电压,看电机转不转得起来,再确认角度极性。

2.4 dq 轴数学模型

经过 Clark 和 Park 变换,我们终于得到了 PMSM 在 dq 旋转坐标系下的数学模型。这才是我们写控制算法的真正起点。

电压方程:

u_d = R_s * i_d + L_d * di_d/dt - ω_e * L_q * i_q
u_q = R_s * i_q + L_q * di_q/dt + ω_e * (L_d * i_d + ψ_f)

这里:

  • u_d, u_q:d、q 轴电压
  • i_d, i_q:d、q 轴电流
  • R_s:定子电阻
  • L_d, L_q:d、q 轴电感(对于表贴式 PMSM,L_d = L_q;对于内置式,L_d < L_q)
  • ω_e:电角速度
  • ψ_f:永磁体磁链

你仔细看这个方程,会发现一个有趣的现象:d 轴方程里出现了 q 轴电流的耦合项,q 轴方程里出现了 d 轴电流的耦合项。这就是所谓的交叉耦合。做电流环控制时,如果不做解耦,d 轴和 q 轴会互相影响,动态响应会很差。

磁链方程:

ψ_d = L_d * i_d + ψ_f
ψ_q = L_q * i_q

这个很简单,d 轴磁链由两部分组成:电枢反应磁链(L_d * i_d)和永磁体磁链(ψ_f)。q 轴磁链只有电枢反应。

转矩方程:

T_e = 1.5 * p * [ψ_f * i_q + (L_d - L_q) * i_d * i_q]

这个公式太重要了,我建议你把它贴在工位上。

  • 第一项:ψ_f * i_q,这是永磁转矩,由永磁体和 q 轴电流相互作用产生。
  • 第二项:(L_d - L_q) * i_d * i_q,这是磁阻转矩,由 d、q 轴电感差异产生。

对于表贴式 PMSM(L_d = L_q),磁阻转矩为零,转矩只和 i_q 成正比。对于内置式 PMSM(L_d < L_q),我们可以利用负的 i_d 来产生额外的磁阻转矩,这就是弱磁控制的基础之一。

关键洞察: 从转矩方程可以看出,要控制转矩,本质上就是控制 i_q 和 i_d。i_q 控制主转矩,i_d 控制磁链和磁阻转矩。在基速以下,我们通常让 i_d = 0,只调 i_q;到了高速弱磁区,就需要负的 i_d 来削弱磁场。

2.5 知识体系总览

为了让大家对整个变换过程有个直观的印象,我画了一张图。你可以把它当作本章的“地图”。

PMSM 坐标变换与数学模型知识体系 三相静止 abc u_a, u_b, u_c Clark 两相静止 αβ u_α, u_β Park 两相旋转 dq u_d, u_q 转子位置 θ dq 轴数学模型(核心) 电压方程: u_d = R_s·i_d + L_d·di_d/dt - ω_e·L_q·i_q u_q = R_s·i_q + L_q·di_q/dt + ω_e·(L_d·i_d + ψ_f) 转矩方程: T_e = 1.5·p·[ψ_f·i_q + (L_d - L_q)·i_d·i_q]

2.6 本章小结与工程建议

好了,这一章的内容就到这里。我们来捋一捋重点:

  1. Clark 变换:把三相静止 abc 变到两相静止 αβ,等幅值和等功率两种形式,我推荐你用等幅值。
  2. Park 变换:把两相静止 αβ 变到两相旋转 dq,需要转子位置 θ,注意角度定义的一致性。
  3. dq 电压方程:包含电阻压降、电感变化率和交叉耦合项,是电流环设计的基础。
  4. 转矩方程:永磁转矩 + 磁阻转矩,弱磁控制就是通过调节 i_d 来改变磁链和利用磁阻转矩。

一个小建议: 如果你刚开始学,建议在 Simulink 里搭一个开环的 Clark + Park 变换模块,给一组三相正弦波,看看 αβ 和 dq 的波形。亲眼看到正弦波变成直流量,那种感觉比看十遍公式都管用。