2. SynRM 数学模型

同步磁阻电机(SynRM)的转子没有永磁体,也没有励磁绕组,转矩完全来自交直轴磁路磁阻差异所产生的磁阻转矩。要完成参数辨识与磁链建模,第一步必须把物理量之间的数学关系写清楚:在 abc 三相静止坐标系中,电压、磁链、电流都是随时间变化的交流量,且电感矩阵随转子位置周期变化;通过 Park 变换进入 dq 旋转坐标系后,交流量变为直流量,电感矩阵对角化,模型才便于控制器实现和参数提取。本节围绕这两套坐标系,梳理电压方程、磁链方程与转矩方程,并说明它们与后续辨识工作的接口关系。

一、abc 坐标系:物理直观但强耦合

在 abc 三相静止坐标系下,SynRM 的定子电压方程由电阻压降与磁链变化率构成。由于转子凸极效应,定子三相自感与互感均随转子电角度 θr 呈 2θr 周期变化,因此三相之间不是解耦的,磁链方程中会出现随位置变化的电感矩阵。

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(θr)  L_ab(θr)  L_ac(θr)] [i_a]   [ψ_fa]
[ψ_b] = [L_ba(θr)  L_bb(θr)  L_bc(θr)] [i_b] + [ψ_fb]
[ψ_c]   [L_ca(θr)  L_cb(θr)  L_cc(θr)] [i_c]   [ψ_fc]

其中 Rs 为定子相电阻,ψfa、ψfb、ψfc 为永磁磁链分量。对 SynRM 而言,转子无永磁体,该分量理论上为零,但实际中可能因材料剩磁或测量偏置而存在微小偏置,辨识时通常作为可选项处理。abc 坐标系的优点是物理意义直接、便于与实测三相电压电流对应;缺点是电感矩阵时变且耦合,直接用于控制器设计会带来大量三角函数运算,因此工程上几乎不直接使用。

二、dq 坐标系:解耦后的核心模型

通过 Park 变换,将三相量投影到随转子同步旋转的 dq 轴上。d 轴取转子磁阻最小的方向(直轴),q 轴取磁阻最大的方向(交轴)。变换后,SynRM 的电压方程与磁链方程显著简化:

u_d = R_s·i_d + dψ_d/dt − ω_e·ψ_q
u_q = R_s·i_q + dψ_q/dt + ω_e·ψ_d

ψ_d = L_d(i_d, i_q)·i_d
ψ_q = L_q(i_d, i_q)·i_q

T_e = (3/2)·p·(ψ_d·i_q − ψ_q·i_d)
    = (3/2)·p·(L_d − L_q)·i_d·i_q   (线性磁路近似)

式中 ωe 为电角速度,p 为极对数。电压方程中的 ωeψq 与 ωeψd 是旋转反电势交叉耦合项,它们体现了 dq 轴之间的能量交换。磁链方程表明,d 轴磁链主要由 id 建立,q 轴磁链主要由 iq 建立,但实际铁芯存在交叉饱和,Ld 与 Lq 都是 id、iq 的二元函数,这正是参数辨识需要重点处理的对象。

关键结论:SynRM 的转矩只包含磁阻转矩项,没有永磁转矩项。当 Ld > Lq 时,正 iq 与正 id 配合可产生正向转矩;若 id 为零,则转矩为零。因此 SynRM 必须施加一定的 d 轴电流来建立磁链,这与永磁同步电机“id=0 控制”有本质区别。

三、从方程到辨识:参数关系与数据流

参数辨识的目标,是从可测的电压、电流、转速数据中反推出 Rs、Ld(id,iq)、Lq(id,iq) 以及磁链曲线 ψd(id,iq)、ψq(id,iq)。dq 电压方程提供了最直接的观测通道:在稳态下,dψ/dt 项为零,方程退化为代数关系,便于用最小二乘拟合;在动态下,则需借助积分或状态观测器处理微分项。下图给出从 abc 采样到 dq 参数提取的整体数据流与方程接口。

SynRM 数学模型与参数辨识数据流 abc 坐标系 u_a u_b u_c i_a i_b i_c L(θr) 时变耦合 Park 变换 θr → dq 3/2 变换 dq 坐标系 u_d = Rs·i_d + dψ_d/dt − ωe·ψ_q u_q = Rs·i_q + dψ_q/dt + ωe·ψ_d ψ_d = L_d·i_d,ψ_q = L_q·i_q T_e = 1.5p(L_d−L_q)i_d i_q 参数辨识 Rs, Ld, Lq ψd(id,iq) 实测数据 电压 / 电流 / 转速 温度 / 位置角 磁链建模 查表 / 多项式 / 神经网络 交叉饱和补偿 转矩校验 T_e 与实测转矩对比 残差分析 迭代修正 辨识用核心方程接口 稳态:u_d = Rs·i_d − ωe·Lq·i_q,u_q = Rs·i_q + ωe·Ld·i_d 动态:dψ_d/dt = u_d − Rs·i_d + ωe·ψ_q,dψ_q/dt = u_q − Rs·i_q − ωe·ψ_d 转矩:T_e = 1.5p(ψ_d·i_q − ψ_q·i_d),用于验证辨识结果的一致性

四、典型场景与常见误区

在工程实践中,dq 模型最常见的用途有三类:一是电流环设计,需要 Ld、Lq 作为前馈解耦参数;二是无位置传感器控制,需要磁链模型估算转子位置;三是效率优化,需要根据磁链曲线寻找最小铜耗工作点。以一台 5.5 kW、4 极 SynRM 为例,若在额定工况下测得 id=8 A、iq=12 A、Ld=0.12 H、Lq=0.035 H,则按线性公式估算转矩约为 1.5×2×(0.12−0.035)×8×12 ≈ 24.5 N·m。但实际铁芯饱和会使 Ld 在大电流下明显下降,若仍用空载电感值计算,转矩误差可能超过 15%。

常见误区:把 Ld、Lq 当作常数处理。SynRM 的凸极比本身就不高,交叉饱和又进一步压缩有效凸极差,若忽略 id、iq 对电感的影响,会导致转矩估算偏差、电流环解耦失效,甚至无位置传感器控制失稳。正确做法是建立 Ld(id,iq)、Lq(id,iq) 二维查表或拟合模型。

方程/参数 abc 坐标系 dq 坐标系 辨识用途
电压方程 三相独立,含时变电感 含 ωe·ψ 交叉耦合项 稳态拟合 Rs、Ld、Lq
磁链方程 L(θr) 矩阵耦合 ψd=Ld·id,ψq=Lq·iq 构建 ψ(i) 曲线
转矩方程 需先求三相磁链 T_e=1.5p(ψd·iq−ψq·id) 校验辨识一致性
电感特性 随 θr 周期变化 随 id、iq 饱和变化 二维查表建模