永磁同步电机基础:数学模型与核心方程

各位同学,欢迎来到第二章。上一章我们聊了弱磁控制到底要解决什么问题,今天咱们得把基本功打牢——永磁同步电机(PMSM)的数学模型。

说实话,我刚入行那会儿,觉得数学模型就是一堆公式,背下来就行了。直到有一次调试一台高速主轴电机,死活跑不到额定转速以上,翻来覆去查了半天,最后发现是 dq 变换里的角度补偿算错了。嗯,从那以后我才真正理解:公式不是用来背的,是用来理解物理过程的。

这一章,我会带着大家从三相物理模型一步步推导到 dq 坐标系下的方程。你不需要一次性记住所有公式,但一定要理解每个变量代表什么、为什么这么变换。

核心要点:PMSM 的数学模型是所有控制算法的基础。弱磁控制、MTPA、电流环参数整定,全都建立在这套方程之上。

2.1 三相静止坐标系下的 PMSM 模型

先看最原始的模型。一台表贴式永磁同步电机,定子有三相绕组,转子上有永磁体。在三相静止坐标系(ABC 坐标系)下,电压方程长这样:

u_a = R_s * i_a + d(ψ_a)/dt
u_b = R_s * i_b + d(ψ_b)/dt
u_c = R_s * i_c + d(ψ_c)/dt

其中 ψ_a、ψ_b、ψ_c 是各相磁链,包含两部分:自感磁链和永磁体耦合磁链。

磁链方程写出来就比较复杂了:

ψ_a = L_aa * i_a + M_ab * i_b + M_ac * i_c + ψ_f * cos(θ_e)
ψ_b = M_ba * i_a + L_bb * i_b + M_bc * i_c + ψ_f * cos(θ_e - 120°)
ψ_c = M_ca * i_a + M_cb * i_b + L_cc * i_c + ψ_f * cos(θ_e + 120°)

这里 L_aa、L_bb、L_cc 是自感,M_ab、M_bc、M_ca 是互感。θ_e 是转子电角度,ψ_f 是永磁体磁链幅值。

你想想看,这个模型里自感和互感都随转子位置变化——因为磁路在转。这就导致方程是时变的、强耦合的。直接拿这个模型做控制?我试过,非常痛苦。每个控制周期都要重新计算一大堆三角函数和耦合项,实时性根本扛不住。

我的经验:三相模型主要用于电机设计阶段的有限元仿真。做控制算法时,几乎没人直接用这个模型。原因很简单——太复杂,算不过来。

2.2 Clark 变换与 Park 变换

为了解决三相模型的时变耦合问题,前辈们想了个办法:坐标变换。

第一步是 Clark 变换(3s/2s),把三相静止坐标系变换到两相静止坐标系(αβ 坐标系):

[i_α]   [1      -1/2    -1/2  ] [i_a]
[i_β] = [0      √3/2   -√3/2 ] [i_b]
[i_0]   [1/2     1/2     1/2  ] [i_c]

Clark 变换把三个变量变成两个(零序分量通常忽略),但 αβ 轴上的电感仍然随转子位置变化。所以还需要第二步。

第二步是 Park 变换(2s/2r),把 αβ 静止坐标系变换到 dq 旋转坐标系:

[i_d]   [cos(θ_e)   sin(θ_e)] [i_α]
[i_q] = [-sin(θ_e)  cos(θ_e)] [i_β]

Park 变换的关键是:让坐标系跟着转子一起转。这样一来,d 轴始终对准转子磁极方向,q 轴超前 90° 电角度。

我个人习惯把 Park 变换理解为「把相机对准转子拍」——你跟着它转,看到的画面就是静止的。这个比喻虽然不太严谨,但初学者理解起来特别快。

注意:Park 变换中的角度 θ_e 必须是转子电角度,不是机械角度。电角度 = 机械角度 × 极对数。搞错这个,电流环会直接发散。我曾经见过一个同事因为这个 bug 烧了三块驱动板……

2.3 dq 坐标系下的电压方程与磁链方程

经过 Clark + Park 变换后,PMSM 的数学模型变得非常简洁。

磁链方程:

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

这里 L_d 和 L_q 分别是 d 轴和 q 轴电感。对于表贴式 PMSM,L_d ≈ L_q;对于内置式 PMSM,L_q > L_d。这个差异正是磁阻转矩的来源,也是弱磁控制能起作用的关键。

电压方程:

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)

其中 ω_e 是电角速度。注意看,方程里出现了 -ω_e * L_q * i_q 和 +ω_e * (L_d * i_d + ψ_f) 这两项——这就是反电动势和交叉耦合项。

为什么说 dq 模型好?因为在这个坐标系下,所有电感都是常数,不随转子位置变化。你只需要处理两个解耦的直流分量,控制起来简单太多了。

物理量 三相坐标系 dq 坐标系
电感 时变、随位置变化 常数(L_d、L_q)
电流 交流量(正弦波) 直流量(稳态时)
耦合程度 强耦合 弱耦合(仅交叉项)
控制难度 极高 较低(可用 PI 控制)

2.4 转矩方程

终于到了最核心的部分——转矩方程。

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

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

其中 p 是极对数。

这个公式包含两项:

  • ψ_f * i_q —— 永磁转矩,由永磁体磁场和 q 轴电流相互作用产生。这是主要转矩成分。
  • (L_d - L_q) * i_d * i_q —— 磁阻转矩,由 dq 轴电感差异产生。对于内置式 PMSM,L_q > L_d,所以这一项是正的。

说白了,i_q 控制转矩大小,i_d 控制磁场强弱。这就是矢量控制的基本思想——把电流分解成转矩分量和励磁分量,分别控制。

我在做电动汽车驱动系统时,遇到过一个问题:高速时电压饱和,转矩输出不够。后来通过注入负的 i_d(也就是弱磁),利用磁阻转矩补回了一部分转矩损失。这就是弱磁控制的核心——用 i_d 换电压裕量,同时利用磁阻转矩维持输出。

关键洞察:从转矩方程可以看出,即使 ψ_f * i_q 项因为电压限制被削弱,只要 (L_d - L_q) 足够大,磁阻转矩仍然可以撑起一部分输出。这就是为什么内置式 PMSM 比表贴式更适合弱磁运行。

2.5 知识体系总览

为了让大家对本章内容有个整体把握,我画了一张结构图。你可以把它当作 PMSM 数学模型的「地图」——从三相模型出发,经过坐标变换,最终落到 dq 方程和转矩方程。

PMSM 数学模型知识体系 三相静止坐标系模型 时变、强耦合、难控制 Clark 变换 (3s/2s) 三相 → 两相静止 Park 变换 (2s/2r) 静止 → 旋转 dq 坐标系模型 常数电感、弱耦合 磁链方程 ψ_d = L_d·i_d + ψ_f 电压方程 含反电动势与交叉耦合 转矩方程 T_e = 1.5p(ψ_f·i_q + ΔL·i_d·i_q) 核心逻辑:坐标变换将时变系统变为定常系统,使 PI 控制成为可能 弱磁控制本质:通过负 i_d 削弱 ψ_f 的影响,在电压极限内提升转速

2.6 稳态运行与电压极限

在稳态运行时,电流变化率为零(di/dt = 0),电压方程可以简化为:

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

忽略电阻压降(高速时电阻压降占比很小),得到:

u_d ≈ -ω_e * L_q * i_q
u_q ≈ ω_e * (L_d * i_d + ψ_f)

定子电压幅值为:

u_s = √(u_d² + u_q²) ≈ ω_e * √[(L_q * i_q)² + (L_d * i_d + ψ_f)²]

这个公式非常重要。它告诉我们:转速越高,需要的电压越高。但逆变器能输出的最大电压是有限的(受直流母线电压限制)。当 u_s 达到电压极限时,再想提速就必须减小括号里的值——也就是减小 i_q 或者让 i_d 为负。

这就是弱磁控制的物理本质:用负的 d 轴电流抵消永磁体磁链,降低反电动势,从而在电压受限的条件下提高转速。

一个小技巧:在实际项目中,我习惯先根据直流母线电压算出电压极限圆,然后在 i_d-i_q 平面上画出电流极限圆和电压极限椭圆。两个圆的交集就是电机的可行工作区。弱磁控制,说白了就是在这个交集里找最优工作点。

2.7 本章小结

这一章我们走完了从三相模型到 dq 模型的完整推导路径。你不需要记住所有中间公式,但一定要理解以下几点:

  • 坐标变换的目的:把时变系统变成定常系统,让控制变得可行
  • dq 方程的结构:磁链方程简单直接,电压方程包含交叉耦合项
  • 转矩方程的两项:永磁转矩 + 磁阻转矩,后者是弱磁控制的重要抓手
  • 电压极限约束:转速越高电压需求越大,弱磁就是主动降低反电动势