3. SynRM数学模型

同步磁阻电机(SynRM)的无位置传感器控制,本质上是对电机内部磁链与电流关系的实时估计与解耦。要设计观测器、电流环和转速环,必须先建立一套能够准确描述“电压—电流—磁链—转矩”映射关系的数学模型。本节从三相静止坐标系出发,逐步过渡到两相旋转坐标系,给出电压方程、磁链方程与转矩方程,并说明它们在无位置传感器控制中的工程意义。

一、三相静止坐标系下的基本方程

在 ABC 三相静止坐标系中,SynRM 的定子电压方程可写为:

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

其中 R_s 为定子相电阻,ψ_A、ψ_B、ψ_C 为三相定子磁链。由于 SynRM 转子无励磁绕组、无永磁体,磁链完全由定子电流与转子位置共同决定。三相磁链方程可写成矩阵形式:

[ψ_A]   [L_AA  L_AB  L_AC] [i_A]
[ψ_B] = [L_BA  L_BB  L_BC] [i_B]
[ψ_C]   [L_CA  L_CB  L_CC] [i_C]

关键点在于:SynRM 的凸极效应使自感与互感均随转子电角度 θ_e 变化,且变化周期为 2θ_e。这意味着三相坐标系下的电感矩阵是时变的,直接用于控制器设计会带来严重的耦合与非线性。因此,工程上几乎不直接在三相系下做闭环,而是通过坐标变换进入两相旋转坐标系。

二、两相旋转坐标系(d-q 轴)下的方程

采用 d 轴对齐转子磁阻最小方向(即 d 轴电感 L_d 小于 q 轴电感 L_q)的约定,经 Park 变换后,SynRM 在 d-q 旋转坐标系下的电压方程为:

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
ψ_q = L_q·i_q

其中 ω_e 为电角速度,L_d、L_q 为 d、q 轴电感。由于 SynRM 无永磁磁链,ψ_d 与 ψ_q 仅由对应轴电流产生,这使得磁链观测器结构相对简单,但也意味着低速时反电动势幅值很小,位置估计难度显著增加。

将磁链方程代入电压方程,可得电流动态形式:

d i_d/dt = (u_d - R_s·i_d + ω_e·L_q·i_q) / L_d
d i_q/dt = (u_q - R_s·i_q - ω_e·L_d·i_d) / L_q

这两个方程是电流环设计与无位置传感器反电动势观测的核心依据。可以看到,d、q 轴之间存在由 ω_e·L_q·i_q 和 ω_e·L_d·i_d 构成的交叉耦合项,转速越高耦合越强,因此高速区通常需要前馈解耦。

三、转矩方程与磁阻转矩本质

SynRM 的电磁转矩来源于电感凸极效应,其 d-q 轴转矩方程为:

T_e = (3/2)·p·(L_d - L_q)·i_d·i_q

其中 p 为极对数。该式表明:转矩正比于电感差 (L_d - L_q) 与 i_d·i_q 的乘积。若 i_d 或 i_q 为零,则转矩为零。因此 SynRM 必须同时注入 d 轴与 q 轴电流才能产生有效转矩,这与永磁同步电机“i_d=0 控制”有本质区别。

核心结论:SynRM 的转矩由凸极率 ΔL = L_d - L_q 决定,电流角 β(电流矢量与 d 轴夹角)直接决定单位电流产生的转矩大小。最大转矩每安培(MTPA)条件为 i_d = i_q,即 β = 45°(忽略饱和与交叉饱和时)。

四、坐标变换关系与参数对照

从三相静止坐标系到两相旋转坐标系的变换链为:Clarke 变换(ABC→αβ)再 Park 变换(αβ→dq)。下表汇总了两套坐标系下方程的关键差异,便于工程选型与调试对照。

对比项三相静止坐标系 ABC两相旋转坐标系 d-q
电压方程形式三相独立,含时变电感两轴解耦,含交叉耦合项
磁链方程ψ = L(θ_e)·i,矩阵时变ψ_d = L_d·i_d,ψ_q = L_q·i_q
电感参数L_AA、L_AB 等随 2θ_e 变化L_d、L_q 为常数(忽略饱和)
转矩表达需经变换后计算T_e = 1.5p(L_d-L_q)i_d i_q
控制适用性不适合直接闭环电流环、观测器标准坐标系
无位置传感器意义难以提取位置信息可从 ψ_d、ψ_q 或反电动势估计 θ_e

五、工程案例:从方程到观测器接口

假设一台 SynRM 参数为 R_s = 0.5 Ω,L_d = 8 mH,L_q = 20 mH,p = 2,额定电流 10 A。在 d-q 坐标系下,若给定 i_d = 7 A、i_q = 7 A,则转矩为:

T_e = 1.5 × 2 × (0.008 - 0.020) × 7 × 7
    = 1.5 × 2 × (-0.012) × 49
    = -1.764 N·m

负号说明按 L_d < L_q 约定时,需注意电流方向定义。工程中常取 L_q > L_d,并令 i_d 为负或调整 d 轴定义,使转矩为正。这一细节在无位置传感器算法中直接影响位置误差符号,必须在建模阶段统一约定。

调试提示:在编写磁链观测器时,优先使用 d-q 轴电压方程中的 u_d - R_s·i_d 与 u_q - R_s·i_q 积分得到 ψ_d、ψ_q,再通过 θ_e = atan2(ψ_q, ψ_d) 提取位置。低速时积分漂移严重,需引入高频注入或闭环观测器修正。

六、常见误区与复盘

  • 误区一:把 SynRM 当作永磁同步电机处理。SynRM 没有永磁磁链,ψ_d 与 ψ_q 完全由电流决定,转矩公式中不存在永磁转矩项,控制策略必须围绕凸极效应设计。
  • 误区二:忽略交叉饱和。实际 L_d、L_q 随 i_d、i_q 变化,尤其在大电流区,交叉饱和会使转矩公式偏离线性模型,导致观测器位置误差增大。
  • 误区三:d-q 轴定义不统一。不同文献对 d 轴对齐方向、电流正方向、转矩正方向约定不同,直接混用会导致符号错误,建议在项目初期固定一套约定并写入代码注释。
  • 误区四:低速直接依赖反电动势。SynRM 低速时 ω_e·ψ 很小,反电动势信噪比低,必须结合高频注入或脉振电压法才能稳定估计位置。

复盘要点:本节所有方程均建立在“理想线性电感、无交叉饱和、三相平衡”假设之上。实际工程中,参数漂移、逆变器死区、采样延迟都会使模型失配。无位置传感器控制器的鲁棒性,取决于对这些非理想因素的补偿能力,而补偿的起点正是对 d-q 轴电压方程与磁链方程的准确理解。

SynRM 坐标系变换与方程关系图 三相静止坐标系 ABC u_A = R_s·i_A + dψ_A/dt ψ_ABC = L(θ_e)·i_ABC 电感矩阵随 2θ_e 时变 不适合直接闭环控制 Clarke 3/2 变换 两相静止坐标系 αβ u_α = R_s·i_α + dψ_α/dt u_β = R_s·i_β + dψ_β/dt ψ_α、ψ_β 含 θ_e 信息 高频注入常用坐标系 Park θ_e 旋转 两相旋转坐标系 d-q 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,ψ_q = L_q·i_q 电流环与观测器标准坐标系 转矩方程(d-q 轴) T_e = (3/2)·p·(L_d - L_q)·i_d·i_q i_d、i_q 无位置传感器位置估计接口 θ_e = atan2(ψ_q, ψ_d) 或 反电动势积分 低速需高频注入 / 闭环观测器修正 ω_e 反馈 蓝:ABC 坐标系 绿:αβ 坐标系 橙:d-q 坐标系