手把手教你学 Simulink—— 无轴承永磁同步电机（Bearingless PMSM）悬浮力与转矩解耦仿真

目录手把手教你学 Simulink—— 无轴承永磁同步电机Bearingless PMSM悬浮力与转矩解耦仿真一、为什么 Bearingless PMSM 的解耦是命门1.1 普通 PMSM vs Bearingless PMSM1.2 两种主流结构二、系统总体架构三、关键参数教学默认径向两自由度四、Simulink 建模 Step-by-StepStep ① —— ★气隙磁场与力-转矩耦合模型核心■ 气隙磁动势叠加■ 径向悬浮力解析表达式关键结论■ 转矩方程含悬浮电流耦合■ Simulink 实现MATLAB FunctionStep ② —— 转子悬浮动力学Step ③ —— 转矩环 FOC标准但考虑耦合Step ④ —— ★悬浮力控制αβ 系核心解耦■ 位置环外环■ 力-电流逆解解耦关键■ 悬浮电流环2 相独立Step ⑤ —— ★完整交叉耦合补偿框图Step ⑥ —— 偏心观测与主动阻尼进阶Step ⑦ —— 测试场景五、典型结果判读预期数值量级5.1 悬浮精度稳态 1~5 s5.2 转矩-悬浮耦合残差5~5.5 s 转矩阶跃5.3 突加扰动恢复3~3.1 s5.4 零速悬浮5.5~7 s5.5 解耦前后刚度对比六、参数调优与常见坑6.1 关键参数6.2 常见坑七、工程扩展方向八、结论手把手教你学 Simulink—— 无轴承永磁同步电机Bearingless PMSM悬浮力与转矩解耦仿真✅ 本讲目标理解 Bearingless PMSM转矩绕组 悬浮绕组或单绕组复用​ 的两种主流拓扑建立气隙磁场叠加 → 转矩 径向力​ 的解析模型实现转矩-悬浮双闭环解耦控制FOC 转矩环 悬浮力前馈解耦验证μ级悬浮精度、转矩稳态误差1%、偏心扰动抑制一、为什么 Bearingless PMSM 的解耦是命门1.1 普通 PMSM vs Bearingless PMSM维度普通 PMSMBearingless PMSM​功能只产生转矩同时产生转矩 径向悬浮力​轴承机械轴承/磁轴承无机械轴承气隙悬浮​绕组三相转矩绕组三相或多余相复用/附加悬浮绕组​控制变量i_d, i_qi_d, i_q转矩 i_α, i_β悬浮​耦合关系无转矩电流→影响悬浮力悬浮电流→影响转矩​气隙大mm 级极小0.3~0.5mm偏心即擦膛​ 核心矛盾电流既产生转矩又产生径向力二者在数学上耦合物理上必须解耦控制。1.2 两种主流结构结构绕组配置解耦难度单绕组复用型4 相/3 相特殊连接​三相电流中注入特定谐波分量产生悬浮力高强耦合双绕组型转矩绕组悬浮绕组​定子两套独立绕组如 3 相转矩 2 相悬浮中磁路耦合仍存 本讲以双绕组型3 相转矩 2 相悬浮径向两自由度为教学默认——物理清晰、解耦可解析。二、系统总体架构┌──────────────┐ │ 直流母线 V_dc │ └──────┬───────┘ │ ├────────────▶ ┌────────────────────┐ │ │ 转矩逆变器3相 │→ i_a,i_b,i_c └────────────▶ └────────────────────┘ │ ┌──────────────┐ │ 悬浮逆变器2相│→ i_x, i_y径向 └──────┬─────────┘ │ ┌─────────▼──────────┐ │ ★Bearingless PMSM │ │ • 转子永磁体 │ │ • 转矩绕组3相 │ │ • 悬浮绕组2相 │ │ • 气隙磁场叠加 │→ Te, F_x, F_y, ω, x, y └─────────┬──────────┘ │ ┌─────────▼──────────┐ │ 转子动力学悬浮 │ │ m·ẍ F_x - k_x·x - c·ẋ d(t) │ │ m·ÿ F_y - k_y·y - c·ẏ d(t) │ │ J·ω̇ Te - Bω - T_load │ └─────────┬──────────┘ │ ┌─────────▼──────────┐ │ ★解耦控制器 │ │ • 转矩环 FOCdq │ │ • 悬浮环αβ 力控制 │ │ • 交叉耦合前馈补偿 │ │ • 偏心观测主动阻尼 │ └────────────────────────┘三、关键参数教学默认径向两自由度参数值极对数转矩p_t2悬浮力极对数 p_s1与转矩极对数差 1满足悬浮条件转矩绕组 R_t0.5 Ω转矩绕组 L_d/L_q4/5 mH永磁磁链 λ_pm0.08 Wb悬浮绕组 R_s0.3 Ω悬浮绕组 L_x/L_y2 mH互感 L_ts转矩-悬浮0.5 mH磁路耦合转子质量 m1 kg转子惯量 J0.002 kg·m²气隙 g00.4 mm额定转速3000 rpm悬浮刚度目标 100 N/mm控制周期电流环 10 kHz悬浮环 5 kHz开关频率20 kHz悬浮条件p_s p_t ± 1——这是 Bearingless 电机能产生径向力的数学前提麦克斯韦应力法推导。四、Simulink 建模 Step-by-StepStep ① —— ★气隙磁场与力-转矩耦合模型核心这是本讲最硬核的部分。用麦克斯韦应力法推导径向力与转矩。■ 气隙磁动势叠加转矩绕组产生以p_t对极的旋转磁场F_t(θ,t) F_tm·cos(p_t·θ - ω_e·t) ...悬浮绕组产生以p_s p_t - 1对极的旋转磁场F_s(θ,t) F_sx·cos(p_s·θ - ω_e·t) F_sy·sin(p_s·θ - ω_e·t)气隙合成磁场​ → 用麦克斯韦径向力公式F_r(θ) (B²(θ) / (2μ0))·g0 径向应力积分得合力■ 径向悬浮力解析表达式关键结论经推导标准文献结果当p_s p_t - 1时% 径向力x, y 方向 F_x K_f · [ i_q_t · i_x_s (λ_pm/L_d)·i_x_s ] % 近似线性化 ≈ K_fx · i_x_s K_couple_x · i_q_t; F_y K_f · [ i_q_t · i_y_s (λ_pm/L_d)·i_y_s ] ≈ K_fy · i_y_s K_couple_y · i_q_t;其中K_fx ≈ K_fy悬浮力电流刚度N/AK_couple_x K_fx · (λ_pm / (L_d·i_nom))转矩电流对悬浮力的耦合系数反向耦合悬浮电流也通过互感影响转矩在电压方程中体现■ 转矩方程含悬浮电流耦合% 转矩仍主要由转矩绕组 i_q_t 产生 Te 1.5·p_t·[ λ_pm·i_q_t (L_d-L_q)·i_d_t·i_q_t ] 1.5·p_t·L_ts·(i_x_s·di_q_t/dt 项近似) ≈ 1.5·p_t·λ_pm·i_q_t 悬浮电流影响小但存在■ Simulink 实现MATLAB Functionfunction [Te, Fx, Fy] bearingless_pmsm_model(... idt, iqt, ixs, iys, didt, diqt, dixs, diys, omega_e, prm) % 参数 p_t prm.p_t; p_s prm.p_s; lam prm.lambda_pm; Ld prm.Ld; Lq prm.Lq; Lts prm.Lts; % 转矩主要项 悬浮耦合小项 Te 1.5*p_t*(lam*iqt (Ld-Lq)*idt*iqt) ... 0.5*p_t*Lts*(ixs*diqt - iqt*dixs iys*diqt - iqt*diys); % 耦合修正 % 悬浮力刚度 K_f prm.K_f; % 由气隙/极数决定如 5 N/A K_couple K_f * lam / (Ld * prm.I_nom); % 转矩-悬浮耦合 % 径向力线性化模型 Fx K_f * ixs K_couple * iqt; Fy K_f * iys K_couple * iqt; % 注精确模型含 cos(p_t θ) 调制教学用线性化已够 end 工程做法用线性化耦合模型做控制器设计用非线性查表做高精度验证。Step ② —— 转子悬浮动力学% 径向两自由度 m 1; c 0.5; k_x 0; % 被动刚度≈0无机械轴承 % x 方向 m * ddx/dt² Fx_cmd - c*dx/dt d_x(t); % d_x: 偏心/扰动 % y 方向 m * ddy/dt² Fy_cmd - c*dy/dt d_y(t); % 偏心扰动重力不平衡 d_x m*g m*e*omega²*cos(omega*t); d_y m*e*omega²*sin(omega*t);Step ③ —— 转矩环 FOC标准但考虑耦合% 转速环 err_w w_cmd - w_fb; iqt_star Kp_w*err_w Ki_w*∫err_w; iqt_star sat(iqt_star, -Iqmax, Iqmax); % 弱磁区 if w w_base idt_star -weak_mag(w); else idt_star 0; end % dq 电流环转矩绕组 vq_t R_t*iqt Lq*diqt/dt w_e*(Ld*idtlam) w_e*Lts*iys; % 耦合项 vd_t R_t*idt Ld*didt/dt - w_e*Lq*iqt w_e*Lts*ixs; % PI 解耦前馈 vq_cmd Kp_i*(iqt_star-iqt) w_e*(Ld*idtlam) - w_e*Lts*iys; vd_cmd Kp_i*(idt_star-idt) - w_e*Lq*iqt - w_e*Lts*ixs; 注意w_e*Lts*iys项——这是悬浮电流对转矩绕组反电动势的耦合必须前馈减去。Step ④ —— ★悬浮力控制αβ 系核心解耦悬浮力目标让转子居中x→0, y→0。■ 位置环外环% x 方向 Fx_star Kp_x*(0 - x_fb) Kd_x*(0 - dx_fb) M_x*ddx_cmd; % 前馈质量加速度 % y 方向 Fy_star Kp_y*(0 - y_fb) Kd_y*(0 - dy_fb) M_y*ddy_cmd;■ 力-电流逆解解耦关键由Fx K_f·ixs K_couple·iqt反解% 消除转矩耦合项 ixs_star (Fx_star - K_couple * iqt_fb) / K_f; iys_star (Fy_star - K_couple * iqt_fb) / K_f; 这一步就是解耦核心把iqt对悬浮力的影响前馈减掉悬浮电流只听位置环的。■ 悬浮电流环2 相独立% x 相悬浮绕组 vx_s R_s*ixs Lx*dixs/dt w_e*Lts*iqt; % 转矩耦合项 vx_cmd Kp_s*(ixs_star - ixs) - w_e*Lts*iqt; % 前馈解耦 % y 相 vy_s R_s*iys Ly*diys/dt w_e*Lts*iqt; vy_cmd Kp_s*(iys_star - iys) - w_e*Lts*iqt;Step ⑤ —— ★完整交叉耦合补偿框图把所有耦合路径画清楚Simulink 中以前馈实现耦合路径补偿方式iqt → Fx/Fy转矩电流产生额外径向力悬浮力逆解中减K_couple·iqtixs/iys → Te悬浮电流经互感影响转矩转矩电压方程中减ω·Lts·iysiqt → 悬浮反电动势悬浮电压方程中减ω·Lts·iqt转子偏心 → 气隙不均 → 力刚度变化在线观测偏心角调K_f(θ)Step ⑥ —— 偏心观测与主动阻尼进阶% 用 x,y 位置估计偏心角 ecc_angle atan2(y_fb, x_fb); ecc_mag sqrt(x_fb² y_fb²); % 主动阻尼在力指令中加与速度同相的阻尼略大于机械阻尼 Fx_star ... (c_active - c_mech)*dx_fb;Step ⑦ —— 测试场景时间工况验证目标0~1 s启动悬浮x,y 从 0.1mm→0悬浮建力1~3 s加速至 3000 rpm转矩-悬浮并行3~3.1 s突加径向扰动 5 Nx 向​悬浮恢复3.1~5 s稳态悬浮恒速解耦残差5~5.5 s转矩阶跃 50%→100%​转矩变化对悬浮的耦合影响5.5~7 s转速降至 0不落膛零速悬浮稳定性7~8 s偏心质量不平衡激励主动阻尼效果五、典型结果判读预期数值量级5.1 悬浮精度稳态 1~5 s指标无解耦普通 FOC 套用解耦控制​位置 RMS 误差0.05 mm0.002 mm2 μm​最大偏移含不平衡0.15 mm0.01 mm​气隙利用率临界擦膛安全裕度 75%​5.2 转矩-悬浮耦合残差5~5.5 s 转矩阶跃指标无解耦前馈全解耦​转矩阶跃引起 x 偏移0.03 mm​0.001 mm​恢复时间20 ms2 ms​悬浮力电流畸变明显无​5.3 突加扰动恢复3~3.1 s指标数值最大偏移0.008 mm恢复时间±2μm3 ms​力电流峰值2.5 A限幅 3 A5.4 零速悬浮5.5~7 s指标数值位置漂移1 μm是否需要特殊低频算法否解耦模型在 ω0 仍成立因悬浮力与 ω 无关5.5 解耦前后刚度对比控制器悬浮刚度 (N/mm)普通 PI无解耦40且随转速变化解耦前馈​150恒定​六、参数调优与常见坑6.1 关键参数参数太小合适太大K_f力刚度悬浮无力3~8 N/A—K_couple补偿不足精确计算值过补偿振荡Kp_x 位置环刚度低200~500 N/m位置超调Kd_x 微分阻尼不足10~30 N·s/m噪声放大悬浮电流环 Kp力跟踪慢1~3高频振Lts 互感解耦漏项实测值过补偿6.2 常见坑现象原因Fix转子低速偏心振荡解耦未含 ω0 时的互感项检查 Lts 在 v 方程中的前馈转矩阶跃后悬浮跳变K_couple 符号错由Fx Kf·ixs Kc·iqt定符号高速悬浮力下降气隙磁密饱和未建模加 λ_pm(ixs) 饱和修正零速时悬浮发散位置微分噪声用状态观测器出速度两相悬浮电流不平衡绕组阻抗不对称在线辨识 R_sx/R_sy擦膛气隙接触扰动超过最大悬浮力限幅紧急降速七、工程扩展方向✅更真实模型精确麦克斯韦应力积分有限元磁密分布导入转子动力学含陀螺效应高速时 x-y 耦合悬浮绕组谐波与转矩绕组空间谐波交互温升→电阻/磁链漂移→解耦参数自适应✅算法升级滑模控制悬浮强鲁棒自抗扰 ADRC统一处理耦合扰动MPC 多变量转矩两向悬浮协同神经网络在线辨 K_couple✅结构拓展单绕组复用型 Bearingless4 相 CLA 连接轴向-径向联合悬浮3 自由度多转子/叠片转子拓扑✅应用高速泵无密封、无磨损飞轮储能真空高速人工心脏泵无接触、低溶血半导体真空传输洁净无油八、结论✅ 你已完成✅ Bearingless PMSM一套气隙、两重功能​ 的耦合机理麦克斯韦应力法✅线性化耦合模型FxKf·ixsKc·iqtTe1.5pλ·iqt耦合小项✅三层解耦转矩 FOC 悬浮力逆解减耦 电压方程互感前馈✅ 定量验证悬浮精度 2μm、转矩阶跃耦合偏移 0.03mm→0.001mm、刚度 40→150 N/mm✅ 明确p_s p_t ± 1是悬浮力产生的数学必要条件