3、永磁同步电机数学模型:坐标变换(Clark/Park)、电压方程、磁链方程、转矩方程,为弱磁控制奠定基础

永磁同步电机(PMSM)的数学模型是分析其稳态与动态性能、设计高性能控制算法(尤其是弱磁控制)的理论基石。在三相静止坐标系下,PMSM 的电压方程和磁链方程呈现为强耦合、非线性、时变系数的微分方程组,直接求解与控制器设计极为困难。通过坐标变换,将三相交流量等效为两相直流分量,不仅实现了磁链与转矩的解耦,更将时变系数转化为常数系数,使得基于 PI 调节器的矢量控制成为可能。弱磁控制的本质,是在逆变器电压极限圆和电流极限圆约束下,通过调节直轴(d轴)去磁电流分量来削弱永磁体产生的气隙磁场,从而在基速以上扩展转速范围。因此,深入理解坐标变换的数学本质、电压/磁链/转矩方程的物理意义,是掌握弱磁控制策略的前提。

PMSM 坐标变换与数学模型关系图 三相静止坐标系 ABC / U,V,W 电压方程(时变耦合) 磁链方程(非线性) 两相静止坐标系 α-β 轴 等幅值/等功率变换 消除三相耦合 两相旋转坐标系 d-q 轴(随转子旋转) 交流→直流 实现磁链/转矩解耦 Clark T = 2/3 Park θ = ωt d-q 轴数学模型(核心方程) 电压方程 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) 磁链方程 ψ_d = L_d·i_d + ψ_f ψ_q = L_q·i_q 转矩方程 T_e = 1.5·p_n·[ψ_f·i_q + (L_d - L_q)·i_d·i_q] 磁阻转矩(凸极机) 弱磁控制约束条件 电压极限椭圆 (L_d·i_d + ψ_f)² + (L_q·i_q)² ≤ (u_max/ω_e)² 高速时电压约束收紧 电流极限圆 i_d² + i_q² ≤ i_max² 逆变器容量限制

3.1 坐标变换:从三相静止到两相旋转

3.1.1 核心原理:Clark 变换与 Park 变换的数学本质

坐标变换的核心思想是利用数学上的线性变换,将三相绕组中对称分布的正弦交流量映射到等效的两相绕组中。Clark 变换(3s/2s)将三相静止坐标系(A-B-C)变换为两相静止坐标系(α-β),其变换矩阵为:

| i_α |   | 1   -1/2   -1/2 |   | i_A |
| i_β | = | 0   √3/2  -√3/2 | × | i_B |
                (等幅值变换,系数 2/3 已省略)

Park 变换(2s/2r)则将 α-β 轴系旋转 θ 角(θ 为转子 d 轴与 α 轴的电角度),得到随转子同步旋转的 d-q 轴系。变换矩阵为:

| i_d |   | cosθ   sinθ |   | i_α |
| i_q | = | -sinθ  cosθ | × | i_β |

经过 Clark+Park 级联变换后,基波电流在 d-q 轴系中表现为直流量,PI 调节器可实现对直流量的无静差跟踪。需要特别指出,等幅值变换与等功率变换的系数差异(2/3 与 √(2/3))会影响转矩方程中的系数,工程应用中必须保持变换与方程系数的一致性。

3.1.2 操作流程:坐标变换的工程实现步骤

  1. 采样与重构:采集两相电流 i_A、i_B,通过 i_C = -(i_A + i_B) 重构第三相电流。
  2. Clark 变换:将 i_A、i_B、i_C 变换为 i_α、i_β,注意根据应用选择等幅值或等功率变换系数。
  3. 转子位置获取:通过旋变/编码器读取转子机械角度,乘以极对数 p_n 得到电角度 θ。
  4. Park 变换:利用 sinθ、cosθ 将 i_α、i_β 变换为 i_d、i_q。
  5. 反变换:控制器输出的 u_d、u_q 经 Park 逆变换和 Clark 逆变换得到三相调制波,用于 SVPWM 生成。
经验提醒:坐标变换中的角度 θ 必须与转子磁极位置严格对齐。若编码器安装存在零位偏差,会导致 d-q 轴电流耦合,表现为转矩输出波动和效率下降。建议在系统调试初期进行编码器零位标定,通常通过注入 d 轴电流并观察 q 轴响应来校准。

3.1.3 案例拆解:表贴式电机坐标变换的数值验证

某表贴式 PMSM(L_d = L_q = 0.5 mH,ψ_f = 0.1 Wb,p_n = 4),额定电流 10 A,额定转速 3000 rpm。在某一时刻,三相电流采样值为 i_A = 8.66 A,i_B = -4.33 A,i_C = -4.33 A,转子电角度 θ = 30°。经 Clark 变换得 i_α = 8.66 A,i_β = 7.5 A;再经 Park 变换得 i_d = 8.66×cos30° + 7.5×sin30° = 11.25 A,i_q = -8.66×sin30° + 7.5×cos30° = 2.17 A。该结果表明,若直接采用 i_d = 0 控制,则需将电流矢量旋转至 q 轴方向,实际 q 轴电流仅 2.17 A,远未达到额定值。此案例说明坐标变换的正确性直接影响电流环的给定精度。

3.1.4 常见误区与总结

  • 误区一:忽略变换系数。等幅值变换与等功率变换的系数不同,若混用会导致转矩计算偏差达 1.5 倍。
  • 误区二:角度方向错误。Park 变换的旋转方向必须与电机旋转方向一致,否则 d-q 轴定义颠倒,控制失效。
  • 误区三:忽略零序分量。对于星形连接无中线系统,零序电流为零,但三角形连接或故障工况下需额外处理。

3.2 电压方程与磁链方程:d-q 轴解耦模型

3.2.1 核心原理:电压方程的物理意义

在 d-q 旋转坐标系下,PMSM 的定子电压方程为:

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)

式中,R_s 为定子电阻,L_d、L_q 为 d-q 轴电感,ω_e 为电角速度,ψ_f 为永磁体磁链。方程右侧第二项为变压器电动势(电流变化引起),第三项为运动电动势(旋转引起)。d 轴方程中的 -ω_e·L_q·i_q 和 q 轴方程中的 +ω_e·(L_d·i_d + ψ_f) 是交叉耦合项,也是弱磁控制中需要重点处理的项。磁链方程为:

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

磁链方程表明,d 轴磁链由电枢反应磁链 L_d·i_d 与永磁体磁链 ψ_f 叠加而成;q 轴磁链仅由 q 轴电流产生。弱磁控制的本质就是通过负的 i_d 削弱 ψ_d,从而降低反电动势,使电机在高速区仍能满足电压约束。

3.2.2 判断方法:耦合项对控制性能的影响分析

在基速以下,通常采用 i_d = 0 控制,此时 d-q 轴耦合项 -ω_e·L_q·i_q 和 ω_e·ψ_f 随转速升高而增大。当转速较高时,若不进行解耦补偿,d-q 轴电流环会相互干扰,导致电流响应振荡。工程上常用的解耦方法有两种:

  • 前馈解耦:在电流环 PI 输出上叠加解耦项 u_d_ff = -ω_e·L_q·i_q,u_q_ff = ω_e·(L_d·i_d + ψ_f)。该方法依赖精确的电机参数。
  • 反馈解耦:通过状态观测器估计耦合项并补偿,对参数鲁棒性更强,但实现复杂。
参数敏感性分析:当 L_d、L_q 因磁路饱和而发生变化时(尤其在弱磁深度工况下),前馈解耦的补偿精度下降。建议在弱磁控制中采用查表法或在线参数辨识修正电感值,以保证解耦效果。

3.2.3 案例拆解:弱磁临界点的电压约束计算

某内置式 PMSM(IPM)参数:L_d = 0.3 mH,L_q = 0.8 mH,ψ_f = 0.12 Wb,R_s = 0.05 Ω,p_n = 4,直流母线电压 U_dc = 300 V,最大相电压幅值 u_max = U_dc/√3 ≈ 173 V(SVPWM 线性调制区)。额定转速 2000 rpm(ω_e = 2π×2000/60×4 ≈ 837.8 rad/s)。在额定工况下,忽略电阻压降,电压幅值近似为:

u ≈ ω_e·√[(L_d·i_d + ψ_f)² + (L_q·i_q)²]

若采用 i_d = 0 控制,i_q = 10 A,则 u ≈ 837.8×√[(0.12)² + (0.008)²] ≈ 100.7 V,未超过 173 V。当转速升至 3500 rpm(ω_e ≈ 1466 rad/s)时,同样电流下 u ≈ 176 V,已超过电压极限。此时必须施加负的 i_d(如 -20 A),使 L_d·i_d + ψ_f = 0.3m×(-20) + 0.12 = 0.114 Wb,电压降至 u ≈ 1466×√[(0.114)² + (0.008)²] ≈ 167 V,满足约束。该案例清晰展示了弱磁控制的必要性。

3.2.4 常见误区与总结

  • 误区一:忽略电阻压降。在低速大扭矩工况下,电阻压降占比不可忽略,若在弱磁计算中忽略 R_s,会导致电压裕量估计偏大。
  • 误区二:电感参数恒定假设。IPM 电机的 L_d、L_q 随电流增大而显著下降(磁路饱和),弱磁深度越大,电感越小,需使用增量电感而非表观电感。
  • 误区三:解耦项符号错误。d 轴解耦项为负,q 轴为正,符号错误会加剧耦合而非消除耦合。

3.3 转矩方程:从电磁转矩到弱磁控制策略

3.3.1 核心原理:转矩方程的构成与物理含义

PMSM 的电磁转矩方程为:

T_e = 1.5·p_n·[ψ_f·i_q + (L_d - L_q)·i_d·i_q]

式中第一项 ψ_f·i_q 为永磁转矩,由永磁体磁场与 q 轴电流相互作用产生;第二项 (L_d - L_q)·i_d·i_q 为磁阻转矩,源于 d-q 轴电感不对称(凸极效应)。对于表贴式电机(L_d = L_q),磁阻转矩为零;对于内置式电机(L_d < L_q),当 i_d 为负时,磁阻转矩为正,可辅助提升转矩输出。弱磁控制中,负的 i_d 在削弱磁链的同时,也引入了磁阻转矩,因此 IPM 电机在弱磁区仍能保持较好的转矩输出能力。

3.3.2 操作流程:弱磁控制策略的设计路径

  1. 确定约束边界:根据逆变器容量确定电流极限圆 i_max,根据直流母线电压确定电压极限椭圆 u_max/ω_e。
  2. 划分运行区域:基速以下采用最大转矩电流比(MTPA)控制;基速以上进入弱磁区,沿电压极限椭圆与电流极限圆的交点轨迹运行。
  3. 计算弱磁电流指令:常用方法包括查表法(离线优化)、公式法(基于电压闭环的 PI 调节)和模型预测控制(在线优化)。
  4. 电压反馈弱磁:通过比较电压指令幅值与 u_max 的差值,经 PI 调节器输出负的 i_d 修正量,实现自动弱磁。
  5. 过调制与六步法:在深度弱磁区,可切换至过调制或六步方波模式,进一步扩展电压利用率。
风险提示:弱磁控制中 i_d 为负值,会产生去磁效应。若 i_d 过大(超过退磁电流阈值),可能永久性退磁永磁体,导致电机性能不可逆下降。设计时必须校核最恶劣工况下的去磁安全裕量,并设置 i_d 下限保护。

3.3.3 案例拆解:IPM 电机弱磁区转矩输出对比

某 IPM 电机参数:L_d = 0.2 mH,L_q = 0.5 mH,ψ_f = 0.08 Wb,p_n = 3,i_max = 30 A,u_max = 150 V。在转速 6000 rpm(ω_e = 1885 rad/s)时,电压极限椭圆短半轴为 u_max/ω_e = 0.0796 Wb。若采用纯 i_d = 0 控制,则需 ψ_f = 0.08 Wb > 0.0796 Wb,无法满足电压约束,电机无法运行。采用弱磁控制,令 i_d = -20 A,则 L_d·i_d + ψ_f = 0.2m×(-20) + 0.08 = 0.076 Wb,满足约束。此时 i_q = √(30² - 20²) = 22.36 A。转矩 T_e = 1.5×3×[0.08×22.36 + (0.2m - 0.5m)×(-20)×22.36] = 1.5×3×[1.789 + 0.134] ≈ 8.65 N·m。若采用表贴式电机(L_d = L_q = 0.2 mH),同样条件下转矩 T_e = 1.5×3×0.08×22.36 ≈ 8.05 N·m。可见 IPM 电机在弱磁区因磁阻转矩贡献,转矩输出略高。

3.3.4 常见误区与总结

  • 误区一:认为弱磁必然导致转矩大幅下降。实际上,IPM 电机的磁阻转矩可在弱磁区提供额外转矩,合理设计可延缓转矩衰减。
  • 误区二:忽略电流极限约束。弱磁时 i_d 增大,i_q 必须相应减小,否则总电流超过逆变器容量,导致过流故障。
  • 误区三:电压极限椭圆计算错误。u_max 应取逆变器实际最大输出相电压,而非直流母线电压。SVPWM 线性区为 U_dc/√3,过调制区可提升至 2U_dc/π。

3.4 知识结构总结与工程要点

数学模型 关键表达式 物理意义 对弱磁控制的作用
Clark 变换 i_α = (2/3)(i_A - i_B/2 - i_C/2) 三相→两相静止,消除零序 为 Park 变换提供输入
Park 变换 i_d = i_αcosθ + i_βsinθ 两相静止→旋转,交流变直流 实现 d-q 轴解耦控制
电压方程 u_d = R_s·i_d + L_d·di_d/dt - ω_e·L_q·i_q d 轴回路电压平衡 耦合项影响高速稳定性
磁链方程 ψ_d = L_d·i_d + ψ_f d 轴磁链由电枢反应+永磁体构成 负 i_d 削弱 ψ_d 实现弱磁
转矩方程 T_e = 1.5·p_n·[ψ_f·i_q + (L_d-L_q)·i_d·i_q] 永磁转矩+磁阻转矩 弱磁区磁阻转矩补偿
工程要点复盘:(1)坐标变换是 PMSM 矢量控制的基础,必须保证变换系数、角度方向、零位标定正确;(2)电压方程中的耦合项在高速时不可忽略,需采用前馈或反馈解耦;(3)磁链方程揭示了弱磁控制的本质——通过负 i_d 削弱永磁磁链,以电压约束换取转速扩展;(4)转矩方程表明 IPM 电机在弱磁区具有磁阻转矩优势,应充分利用;(5)弱磁控制设计必须同时满足电压极限椭圆和电流极限圆双重约束,并考虑参数变化和退磁风险。