2. 电机数学模型与坐标变换仿真基础

要在Simulink中搭建一个可用于无位置传感器控制验证的永磁同步电机(PMSM)本体模型,第一步不是拖模块,而是把电机的数学描述写清楚。三相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

其中 R_s 为定子相电阻,ψ_a、ψ_b、ψ_c 为三相绕组磁链。磁链由永磁体磁链与定子电流产生的磁链叠加而成,并随转子电角度 θ_e 变化。由于三相绕组在空间上互差120°电角度,磁链表达式包含 cosθ_e、cos(θ_e−2π/3)、cos(θ_e+2π/3) 等项,导致电感矩阵随转子位置周期性变化。这意味着abc坐标系下的电压方程是时变系数微分方程,直接用于仿真虽然可行,但计算量大、耦合强,不利于控制器设计。

为解决这一问题,需要引入坐标变换。Clarke变换将三相静止坐标系下的量投影到两相静止αβ坐标系,其等幅值形式为:

i_α = (2/3)(i_a − 0.5·i_b − 0.5·i_c)
i_β = (2/3)((√3/2)·i_b − (√3/2)·i_c)

若三相电流满足 i_a + i_b + i_c = 0,则 i_α = i_a,i_β = (i_a + 2i_b)/√3。Clarke变换把三个耦合量压缩为两个正交量,但αβ轴上的电感仍随转子位置变化,因为转子永磁体方向尚未与坐标轴对齐。

Park变换进一步把αβ坐标系旋转到与转子磁链方向一致的dq坐标系。定义d轴与转子永磁体磁链方向重合,q轴超前d轴90°电角度,则:

i_d =  i_α·cosθ_e + i_β·sinθ_e
i_q = −i_α·sinθ_e + i_β·cosθ_e

经过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)

其中 ω_e 为电角速度,ψ_f 为永磁体磁链,L_d、L_q 为dq轴电感。这两个方程是后续电流环设计、反电动势观测器和滑模观测器的基础。转矩方程则为 T_e = 1.5·p·[ψ_f·i_q + (L_d−L_q)·i_d·i_q],p 为极对数。对于表贴式电机 L_d≈L_q,磁阻转矩项可忽略,转矩近似与 i_q 成正比。

PMSM 坐标变换与 dq 轴动态模型信号链 abc 三相静止坐标系 u_a, u_b, u_c i_a, i_b, i_c ψ_a, ψ_b, ψ_c 时变电感矩阵 L(θ_e) Clarke 变换 3/2 变换 abc → αβ αβ 两相静止坐标系 i_α, i_β u_α, u_β ψ_α, ψ_β 电感仍随 θ_e 变化 Park 变换 旋转变换 αβ → dq 转子电角度 θ_e 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] 耦合项:−ω_e·L_q·i_q 与 +ω_e·(L_d·i_d + ψ_f) 关键参数 R_s = 2.875 Ω L_d = L_q = 8.5 mH ψ_f = 0.175 Wb p = 4 J = 0.003 kg·m²

在Simulink中搭建电机本体模型时,推荐采用“dq轴方程 + 坐标反变换”的结构,而不是直接搭abc三相时变方程。具体做法是:用两个积分器分别求解 i_d 和 i_q,输入为 u_d、u_q 和 ω_e,输出 i_d、i_q 再经Park反变换和Clarke反变换得到三相电流。这样模型内部是常系数,数值稳定性好,也便于后续接入观测器。

搭建步骤可以按以下顺序进行。第一步,在Simulink中新建空白模型,从Simscape Electrical或Simulink基础库中选取增益、求和、积分、三角函数、乘法等模块。第二步,建立dq轴电流积分回路:将 u_d 减去 R_s·i_d 再加上 ω_e·L_q·i_q,除以 L_d 后送入积分器得到 i_d;同理,将 u_q 减去 R_s·i_q 再减去 ω_e·(L_d·i_d + ψ_f),除以 L_q 后积分得到 i_q。第三步,用 i_d、i_q 和 θ_e 计算 T_e,再经机械方程 dω_m/dt = (T_e − T_L − B·ω_m)/J 得到机械角速度,乘以极对数得到 ω_e,积分得到 θ_e,反馈回Park变换。第四步,用θ_e对i_d、i_q做Park反变换得到i_α、i_β,再用Clarke反变换得到i_a、i_b、i_c,作为电机模型的三相电流输出。

参数设置方面,以一台典型表贴式PMSM为例:定子电阻 R_s = 2.875 Ω,dq轴电感 L_d = L_q = 8.5 mH,永磁体磁链 ψ_f = 0.175 Wb,极对数 p = 4,转动惯量 J = 0.003 kg·m²,阻尼系数 B = 0.0005 N·m·s。这些数值属于教学示例,实际使用时需替换为所研究电机的铭牌或参数辨识结果。仿真步长建议取 1e-5 s 或更小,求解器选用 ode4 或 ode23tb,以保证电流环动态不被数值误差淹没。

坐标变换的等幅值约定与等功率约定会影响系数。本节采用等幅值Clarke变换,即 i_α = i_a 在平衡条件下成立,转矩公式中的1.5系数即来源于此。若改用等功率变换,系数需相应调整,全文与图示必须保持一致。

模型搭好后,先做开环验证:给定 u_d = 0、u_q = 常数,观察 i_q 是否跟随、转速是否上升、三相电流是否对称。若 i_d 出现明显非零稳态,说明Park变换角度或反变换符号有误;若转速不上升,检查转矩方向与负载符号。只有开环波形正确后,才适合接入无位置传感器观测器进行闭环验证。