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

目录手把手教你学 Simulink—— 无轴承永磁同步电机Bearingless PMSM悬浮力与转矩解耦仿真一、Bearingless PMSM 是什么为什么难1.1 与普通 PMSM 的区别1.2 耦合从哪来二、系统总体架构三、关键参数教学默认四、Simulink 建模 Step‑by‑StepStep ① —— 电机本体耦合磁路模型核心■ 气隙磁密分解■ 电磁转矩虚位移法简化■ 径向悬浮力麦克斯韦应力法教学近似Step ② —— 机械运动方程Step ③ —— 解耦控制策略核心中的核心■ 控制目标■ 耦合补偿思路■ 双闭环结构Simulink 实现Step ④ —— 悬浮力观测与解耦度评估Step ⑤ —— 测试场景五、典型结果判读预期数值量级5.1 解耦效果对比最关键5.2 悬浮环性能5.3 转速环性能5.4 解耦度定量六、参数调优与常见坑6.1 关键增益6.2 常见坑七、工程扩展方向八、结论九、下一步可继续手把手教你学 Simulink—— 无轴承永磁同步电机Bearingless PMSM悬浮力与转矩解耦仿真✅ 本讲目标理解 Bearingless PMSM 为什么“一台既要转矩又要悬浮”且二者强耦合​建立含悬浮力模型的电机本体不依赖商业电机库黑盒实现转速–悬浮双闭环 耦合补偿​ 的全数字控制验证解耦效果转矩变化时悬浮力波动、悬浮扰动下转速纹波一、Bearingless PMSM 是什么为什么难1.1 与普通 PMSM 的区别维度普通 PMSMBearingless PMSM​功能只产生电磁转矩转矩 径向悬浮力​绕组三相转矩绕组转矩绕组 悬浮绕组通常两相附加​气隙磁场旋转磁场转矩叠加 2 极或 4 极悬浮磁场​机械结构需要轴承无机械接触磁悬浮​控制目标转速/位置转速 径向位移x,y双目标​耦合关系—转矩电流 ↔ 悬浮力强耦合​1.2 耦合从哪来气隙合成磁场 转矩磁场P 对极 悬浮磁场通常为 2 极或 4 极二者相互作用产生麦克斯韦应力既产生切向转矩也产生径向力**转矩电流变化 → 改变气隙磁密分布 → 影响悬浮力悬浮电流变化 → 改变磁场分布 → 影响转矩 不解耦 → 转速波动引起悬浮振动悬浮调节引起转矩纹波系统失稳。二、系统总体架构┌──────────────┐ │ 直流母线 │ V_dc └──────┬───────┘ │ ├────────────▶ ┌────────────────────┐ │ │ 转矩逆变器3相 │ │ │ SVPWM │→ i_A,i_B,i_C (转矩) └────────────▶ └────────────────────┘ │ ┌────────────────────┐ │ 悬浮逆变器2相 │ │ 独立 H 桥或 4 相 │→ i_α,i_β (悬浮) └──────────┬─────────┘ │ ┌──────▼───────┐ │ 电机本体模型 │ │ • 转矩绕组 │ │ • 悬浮绕组 │ │ • 耦合磁路 │ │ • 麦克斯韦应力 │→ T_e, F_x, F_y, θ, ω └──────┬───────┘ │ ┌────────────▼───────────┐ │ 传感层 │ │ • 光电编码器(θ,ω) │ │ • 电涡流位移传感器(x,y) │ │ • 电流传感器 │ └────────────┬───────────┘ │ ┌────────────▼───────────┐ │ 解耦控制器 │ │ • 转速环 PI → i_q* │ │ • 悬浮环 PI → i_α*,iβ* │ │ • 耦合补偿前馈 │ └─────────────────────────┘三、关键参数教学默认参数转矩绕组悬浮绕组极对数22悬浮磁场极对数定子电阻0.5 Ω0.5 Ω电感5 mH5 mH永磁磁链0.1 Wb—互感耦合L_mt 0.5 mH教学简化转动惯量0.01 kg·m²粘滞摩擦0.001 N·m·s转子质量2 kg径向刚度机械0无轴承纯磁悬浮阻尼0.5 N·s/m控制周期电流环 10 kHz转速/悬浮环 1 kHz开关频率20 kHz四、Simulink 建模 Step‑by‑StepStep ① —— 电机本体耦合磁路模型核心不依赖黑盒用解析耦合模型教学力与转矩可解释■ 气隙磁密分解B_gap B_PM B_torque B_suspensionB_PM永磁体产生的 2 极磁场随转子位置旋转B_torque转矩绕组电流产生的 P 对极磁场B_suspension悬浮绕组电流产生的 2 极磁场可控■ 电磁转矩虚位移法简化T_e (3/2)·P·[λ_pm·i_q (L_d − L_q)·i_d·i_q]与普通 PMSM 一致转矩绕组 dq 变换。■ 径向悬浮力麦克斯韦应力法教学近似关键悬浮力正比于转矩磁场与悬浮磁场的相互作用F_x K_f·[i_α·f_x(θ) i_β·f_y(θ)] K_couple·i_q·g_x(θ) F_y K_f·[i_α·f_y(θ) − i_β·f_x(θ)] K_couple·i_q·g_y(θ)其中K_f悬浮力系数与磁路设计有关f_x(θ), f_y(θ)悬浮磁场空间分布基函数2 极K_couple耦合系数转矩电流 i_q 对悬浮力的影响g_x(θ), g_y(θ)耦合项空间分布 这就是“耦合”的数学表达i_q 不仅产生转矩还通过通过K_couple项影响 F_x,F_y。Simulink 实现MATLAB Functionfunction [Te, Fx, Fy] bearingless_pmsm(id, iq, ialpha, ibeta, theta, prm) % 永磁磁链 lambda_pm prm.lambda_pm; % 电磁转矩标准 PMSM 公式 Te 1.5 * prm.P * (lambda_pm * iq (prm.Ld - prm.Lq) * id * iq); % 悬浮磁场基函数2 极 fx cos(theta); fy sin(theta); % 悬浮力主项 F_main prm.Kf * (ialpha * fx ibeta * fy); F_main_y prm.Kf * (ialpha * fy - ibeta * fx); % 耦合项转矩电流影响悬浮力 F_couple_x prm.Kcouple * iq * cos(theta prm.couple_phase); F_couple_y prm.Kcouple * iq * sin(theta prm.couple_phase); Fx F_main F_couple_x; Fy F_main_y F_couple_y; endStep ② —— 机械运动方程转子旋转J·dω/dt Te − B·ω − TL dθ/dt ω转子径向平移磁悬浮m·ẍ F_x − c·ẋ − k_mech·x (k_mech0无机械轴承) m·ÿ F_y − c·ẏ − k_mech·y 注意径向运动是二阶系统悬浮控制本质是主动阻尼刚度设计。Step ③ —— 解耦控制策略核心中的核心■ 控制目标转速环ω → ω_ref输出 i_q*悬浮环x → 0y → 0输出 i_α, i_β■ 耦合补偿思路问题i_q 变化 → F_x,F_y 变化 → 转子偏移 → 需要 i_α,i_β 补偿解法前馈补偿——在悬浮控制器输出中**减去耦合项预估% 悬浮电流指令含耦合补偿 i_alpha_cmd i_alpha_pid − K_comp·iq·cos(theta couple_phase); i_beta_cmd i_beta_pid − K_comp·iq·sin(theta couple_phase);其中K_comp ≈ K_couple / K_f模型逆补偿。■ 双闭环结构Simulink 实现转速环外环1 kHzerr_w omega_ref − omega; iq_star Kp_w*err_w Ki_w*integral(err_w); iq_star sat(iq_star, −I_max, I_max);悬浮环外环1 kHzerr_x 0 − x; err_y 0 − y; i_alpha_pid Kp_x*err_x Ki_x*integral(err_x); i_beta_pid Kp_y*err_y Ki_y*integral(err_y); % 耦合补偿 i_alpha_star i_alpha_pid − K_comp*iq*cos(thetacouple_phase); i_beta_star i_beta_pid − K_comp*iq*sin(thetacouple_phase);电流环内环10 kHz转矩绕组id,iq → SVPWM标准 PMSM 电流环悬浮绕组i_α,i_β → 独立 PI → H 桥 PWMStep ④ —— 悬浮力观测与解耦度评估解耦度指标定量Decoupling Index ||ΔF due to Δiq|| / ||ΔF due to Δiα,β||理想解耦 → 该项 ≈ 0iq 变化不引起悬浮力变化。Simulink 中实现在补偿前后分别记录 F_x,F_y 对 i_q 阶跃的响应计算传递函数或脉冲响应面积比Step ⑤ —— 测试场景时间工况0~2 s空载启动至 3000 rpm2~4 s悬浮给定阶跃x: 0→0.1 mm4~6 s负载转矩阶跃 0→0.5 N·m6~8 s关闭耦合补偿对比8~10 s恢复补偿 转速反转五、典型结果判读预期数值量级5.1 解耦效果对比最关键工况无补偿悬浮波动有补偿悬浮波动​转矩 0→0.5 N·m 阶跃±0.05 mm​±0.005 mm​转速 3000→−3000 rpm±0.08 mm​±0.01 mm​悬浮指令阶跃转速纹波 ±15 rpm转速纹波 ±3 rpm​ 耦合补偿使悬浮刚度提升约 10 倍转速纹波降低 5 倍。5.2 悬浮环性能指标数值静态悬浮误差 0.005 mm带宽~50 Hz电流环限制抗扰外力冲击恢复时间 20 ms5.3 转速环性能指标数值稳态误差 0.1%负载阶跃恢复 50 ms反转时间 200 ms5.4 解耦度定量Decoupling Index补偿前 0.35 Decoupling Index补偿后 0.03接近 0 → 说明悬浮力与转矩基本实现动态解耦。六、参数调优与常见坑6.1 关键增益参数太小合适太大Kp_w转速响应慢0.1~0.5超调振荡Ki_w静差5~20积分饱和Kp_x悬浮刚度低50~200高频振荡Ki_x静差500~2000噪声放大K_comp补偿不足K_couple/K_f过补偿→不稳定6.2 常见坑现象原因Fix悬浮振荡发散K_comp 过大或符号错半实物辨识 K_couple转速纹波大悬浮环耦合未补偿检查 θ 同步、补偿相位启动抖动初始悬浮力不足预加压定位 电流斜坡高速失稳反电动势影响悬浮模型高速区修正 K_f(ω)解耦度不降补偿模型失配在线辨识耦合项七、工程扩展方向✅更真实模型有限元提取 K_f,K_couple(θ,id,iq) 非线性表齿槽效应、磁饱和、涡流阻尼五自由度轴向角向悬浮✅算法升级滑模控制悬浮环强鲁棒自抗扰控制ADRC统一处理耦合模型预测控制MPC多目标优化✅应用拓展飞轮储能真空高速人工心脏泵无磨损精密机床主轴纳米级悬浮八、结论✅ 你已完成✅ Bearingless PMSM耦合磁路解析模型转矩悬浮力同源于气隙磁场✅双闭环解耦控制转速环 悬浮环 耦合前馈补偿✅ 定量验证补偿后悬浮波动降低 10 倍解耦度从 0.35→0.03✅ 明确调参顺序、耦合补偿物理意义、工程扩展路径