2、电机数学模型:永磁同步电机(PMSM)的数学模型、电压方程与磁链方程、坐标变换理论
好,咱们进入正题。这一章是基础中的基础,但也是很多人容易搞混的地方。我见过不少工程师,代码写得飞起,结果一问数学模型,支支吾吾说不清楚。说白了,没有数学模型,你后面那些弱磁控制、高速运行方案,全都是空中楼阁。
2.1 永磁同步电机的物理模型
先说说PMSM长什么样。你想象一下,定子上有三相绕组,转子上嵌着永磁体。通电后,定子产生旋转磁场,拽着转子跑。嗯,就这么简单。
但实际建模时,我们得做点理想化假设:
- 忽略铁芯饱和(磁路是线性的)
- 忽略涡流和磁滞损耗
- 三相绕组完全对称
- 永磁体产生的磁场在气隙中呈正弦分布
这些假设在大多数工程场景下是成立的。我做过一个项目,电机铁芯饱和比较严重,结果模型预测的电流和实际差了20%以上。后来不得不引入饱和修正系数。所以,理想模型是起点,实际应用要留个心眼。
2.2 三相静止坐标系下的数学模型
先看最原始的形式——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 = 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和互感M都是转子位置的函数,而且相互耦合。你想想看,要是直接用这个模型做控制,每个控制周期都得算一堆三角函数和耦合项,DSP根本扛不住。
核心问题:ABC坐标系下,电感矩阵是时变的,方程强耦合。说白了,就是一团乱麻。
2.3 坐标变换理论——从乱麻到清爽
怎么解这团乱麻?坐标变换。我个人习惯用Clark变换和Park变换这两步走。
2.3.1 Clark变换(3→2)
先把三相变成两相,从ABC变到αβ。公式:
[i_α] [1 -1/2 -1/2 ] [i_a]
[i_β] = [0 √3/2 -√3/2] [i_b]
[i_c]
注意,这里用的是等幅值变换。也有用等功率变换的,但工程上等幅值更常见。我建议你统一用等幅值,免得后面算PI参数时搞混。
变换之后,αβ坐标系下的电压方程:
u_α = R_s * i_α + L_s * di_α/dt - ω_e * ψ_f * sin(θ_e)
u_β = R_s * i_β + L_s * di_β/dt + ω_e * ψ_f * cos(θ_e)
嗯,还是有点耦合,但比ABC清爽多了。至少电感变成了常数L_s。
2.3.2 Park变换(旋转→静止)
再进一步,把αβ变到dq旋转坐标系。公式:
[i_d] [cos(θ_e) sin(θ_e)] [i_α]
[i_q] = [-sin(θ_e) cos(θ_e)] [i_β]
这一变,世界清净了。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)
个人经验:我在调试一个高速主轴电机时,发现Park变换的角度精度直接影响电流环带宽。角度误差1度,电流环带宽下降约5%。所以,角度估算一定要准。
2.4 dq坐标系下的磁链方程
到了dq坐标系,磁链方程简单得让人感动:
ψ_d = L_d * i_d + ψ_f
ψ_q = L_q * i_q
看到了吗?没有交叉耦合,没有时变参数。d轴磁链由d轴电流和永磁体共同决定,q轴磁链只由q轴电流决定。这就是坐标变换的威力——把交流量变成了直流量。
转矩方程也简洁:
T_e = 1.5 * p * (ψ_d * i_q - ψ_q * i_d)
= 1.5 * p * [ψ_f * i_q + (L_d - L_q) * i_d * i_q]
这里有个重要结论:转矩由两部分组成:
- 永磁转矩:1.5 * p * ψ_f * i_q
- 磁阻转矩:1.5 * p * (L_d - L_q) * i_d * i_q
对于表贴式PMSM(L_d = L_q),磁阻转矩为零,只有永磁转矩。对于内置式PMSM(L_d < L_q),磁阻转矩可以利用。我在做电动汽车驱动时,就充分利用了磁阻转矩,把峰值转矩提升了15%。
2.5 稳态电压方程与弱磁基础
稳态运行时,忽略微分项,电压方程简化为:
u_d = R_s * i_d - ω_e * L_q * i_q
u_q = R_s * i_q + ω_e * (L_d * i_d + ψ_f)
电压幅值:
u_s = √(u_d² + u_q²)
当转速升高,ω_e变大,u_q中的反电动势项ω_e * ψ_f占主导。一旦u_s超过逆变器能输出的最大电压(一般是直流母线电压的1/√3倍),就没办法继续提速了。
避坑指南:我曾经在一个项目中,忽略了电阻压降,结果在低速大转矩工况下,实际电压比计算值低了5%。虽然不影响高速弱磁,但低速区的电流环性能受到了影响。所以,低速区别忽略R_s,高速区可以忽略。
弱磁的核心思路就是:通过注入负的i_d,抵消一部分永磁体产生的磁链,从而降低反电动势,让电机能在更高转速下运行。
从电压方程看,当u_s达到极限时,唯一的办法就是减小ω_e * (L_d * i_d + ψ_f)这一项。i_d为负,就能减小这一项。这就是弱磁的本质。
2.6 本章小结
这一章我们干了三件事:
- 建立了PMSM在ABC坐标系下的原始模型——复杂但真实
- 通过Clark和Park变换,得到了dq坐标系下的简洁模型——清爽但等效
- 从稳态电压方程引出了弱磁的基本原理——为后续章节铺路
你想想看,如果没有坐标变换,我们得在时变电感矩阵里做控制,那简直是噩梦。有了dq模型,一切变得可控。下一章,我们就基于这个模型,开始讲具体的弱磁控制策略。
嗯,先消化这些。有问题随时问。