手把手教你学 Simulink—— 磁场调制电机（FMM）的新型控制策略建模与仿真-CSDN博客

✅ 本讲目标：

理解 FMM 的核心：“气隙磁场调制”带来高转矩密度 / 低速大转矩 / 天然减速

建立 FMM 的等效电路 + 调制传递函数 模型（不依赖有限元黑盒）

实现 磁场调制感知的新型控制：调制比优化、低速大转矩 MTPA、弱磁扩速、转矩–悬浮/多端口解耦

定量对比：传统矢量控制 vs FMM‑aware 控制 的转矩密度、效率、调速范围、动态性能

一、为什么 FMM 是“新型控制”的热土？

1.1 FMM 家族与共同特征

类型	代表结构	核心特征
Magnetic Gear (MG)	同轴磁齿轮	无接触传动、调速比 = 极数比对
Flux‑Switching PM (FSPM)	定子永磁型	高转矩密度、转子无永磁
Vernier Machine	磁场调制电机	低速大转矩、转矩密度 ∝ 极数比
Dual‑Stator / Dual‑Rotor	双定子/双转子	多端口功率分流
Modular Spoke‑Type	模块化轮辐式	磁场调制 + 聚磁效应

📌 共同物理本质：气隙中存在“调制磁场”（通常极对数不同于电枢绕组），通过调制实现转矩/速度变换，类似“电磁变速箱”。

1.2 对控制的挑战与机会

挑战	机会（新型控制切入点）
等效极对数动态变化（调制效应）	调制比感知控制：实时估计调制状态优化电流
高凸极率（磁阻转矩占比大）	MTPA 考虑调制项，而非简单 dq MTPA
低速大转矩易饱和	饱和感知弱磁/扩速策略
多端口（双定子/转子）	功率分配 + 解耦控制（类似多电平/多电机协同）
磁场调制带来额外损耗	效率最优电流剖面（损耗模型 + 调制因子）

二、系统总体架构（FMM‑aware 控制）

┌──────────────┐
│ 直流母线       │ V_dc
└──────┬───────┘
       │
       ├────────────▶ ┌────────────────────┐
       │              │ 功率变换器          │
       │              │ （多相/双逆变器）   │→ i_abc (电枢)
       └────────────▶ └────────────────────┘
                      │
                      ┌────────────────────┐
                      │ FMM 本体（调制模型） │
                      │ • 电枢绕组          │
                      │ • 调制磁场（PM/磁阻）│
                      │ • 等效极对数 p_eq    │
                      │ • 调制传递函数 H_mod │→ T_e, ω, θ
                      └──────┬─────────────┘
                             │
                ┌────────────▼───────────┐
                │ 传感与状态估计           │
                │ • 电流/电压             │
                │ • 位置/速度             │
                │ • 调制比估计 p_eq_hat   │
                │ • 磁链/饱和观测         │
                └────────────┬───────────┘
                             │
                ┌────────────▼───────────┐
                │ 新型控制器               │
                │ • 调制感知矢量控制       │
                │ • MTPA+FMM 修正项       │
                │ • 弱磁/扩速（饱和感知）  │
                │ • 多端口功率分配（可选）  │→ i_d*, i_q*, i_x*
                └─────────────────────────┘

三、关键建模：FMM 等效电路 + 调制传递函数

3.1 核心思想（教学简化但物理正确）

把 FMM 看成：

电枢电流 产生电枢磁场（极对数 p_a）
调制磁场（永磁/磁阻结构）具有空间调制函数 M(θ)（周期 = 机械周期 / p_mod）
气隙合成磁场 = 电枢磁场 × M(θ) → 产生等效极对数
```
p_eq = |p_a ± p_mod|
```
电磁转矩正比于 p_eq × 磁链 × 电流（类似 PMSM，但 p_eq 由调制决定）

3.2 调制传递函数（关键创新点）

定义 调制传递函数 H_mod(s, θ)：描述电枢电流 → 有效气隙磁场的传递关系（含空间谐波、饱和、频率依赖）。

教学用简化形式：

B_gap_eff(s) = H_mod(s) · i_dq(s)
H_mod(s) = K_mod / (1 + τ_mod·s) · [1 + Σ_k α_k·cos(k·p_mod·θ)]

其中：

K_mod：调制增益（与结构有关）
τ_mod：调制磁场动态时间常数（磁路弛豫）
α_k：空间谐波系数（磁场调制带来的高次谐波）

3.3 电磁转矩方程（FMM 版）

T_e = (3/2)·p_eq·[λ_pm_eff·i_q + (L_d − L_q)·i_d·i_q] + T_mod_coupling

其中：

λ_pm_eff = λ_pm · H_mod(0)（调制后的有效磁链）
T_mod_coupling：调制磁场与电枢电流的交叉耦合转矩（FMM 特有，传统 PMSM 忽略）

📌 传统矢量控制把 p_eq 当常数、T_mod_coupling 忽略 → 在 FMM 中次优甚至不稳定。

四、Simulink 建模 Step‑by‑Step

Step ① —— FMM 本体模型（调制感知）

参数定义（Vernier‑type FMM 示例）：

prm.p_a = 3;          % 电枢绕组极对数
prm.p_mod = 4;        % 调制磁场极对数（结构决定）
prm.p_eq = prm.p_a + prm.p_mod;  % 等效极对数 = 7
prm.lambda_pm = 0.1;  % 永磁磁链
prm.Ld = 5e-3; prm.Lq = 8e-3;  % 凸极（磁阻转矩大）
prm.K_mod = 1.2;      % 调制增益
prm.tau_mod = 2e-3;   % 调制动态
prm.alpha = [0.1, 0.05];  % 空间谐波

Simulink 实现（MATLAB Function）：

function [Te, lambda_pm_eff] = fmm_torque(id, iq, theta, prm)
    % 调制传递函数（简化：稳态增益 + 空间谐波）
    H_mod = prm.K_mod * (1 + prm.alpha(1)*cos(prm.p_mod*theta) ...
                       + prm.alpha(2)*cos(2*prm.p_mod*theta));
    
    % 有效磁链（调制后）
    lambda_pm_eff = prm.lambda_pm * H_mod;
    
    % 电磁转矩（含调制耦合项）
    T_pm = 1.5 * prm.p_eq * lambda_pm_eff * iq;
    T_rel = 1.5 * prm.p_eq * (prm.Ld - prm.Lq) * id * iq;
    
    % 调制耦合转矩（FMM 特有）
    T_mod_coupling = 0.1 * prm.p_eq * iq * sin(prm.p_mod*theta);
    
    Te = T_pm + T_rel + T_mod_coupling;
end

Step ② —— 机械运动 + 调制动态

% 机械
J·dω/dt = Te − B·ω − TL
dθ/dt = ω

% 调制磁场动态（一阶滞后，可选）
% τ_mod·dH_mod/dt + H_mod = H_mod_ss(θ, i)

Step ③ —— 传统矢量控制（Baseline）

标准 PMSM 矢量控制（忽略调制）：

% 转速环
iq_star = Kp_w*(wr_ref − wr) + Ki_w*∫(wr_ref − wr);
iq_star = sat(iq_star, −Imax, Imax);

% 弱磁（传统）
if wr > wr_base
    id_star = −(wr − wr_base)*K_fw;  % 简单弱磁
else
    id_star = 0;
end

% dq 电流环 → SVPWM

❌ 问题：p_eq 被当成常数（实际由调制决定），T_mod_coupling 未补偿，λ_pm_eff 未随 θ 更新。

Step ④ —— ★新型控制：调制感知矢量控制（FMM‑aware）

■ 改进 1：等效极对数与有效磁链实时估计

% 在线估计（用电流+转矩观测）
p_eq_hat = estimate_p_eq(theta, i_dq, Te_meas);
lambda_pm_eff_hat = Te_meas / (1.5*p_eq_hat*iq + ...) ;  % 逆解

教学用：已知结构参数直接算（p_eq, H_mod(θ)），工程用观测器。

■ 改进 2：MTPA 考虑调制项

传统 MTPA：i_d = (λ_pm − sqrt(λ_pm² + 4(Ld−Lq)²·i_q²)) / (2(Ld−Lq))

FMM‑MTPA：把 λ_pm 替换为 λ_pm_eff(θ)，且考虑 T_mod_coupling 对最优点的偏移：

% FMM‑aware MTPA
lambda_eff = prm.lambda_pm * H_mod(theta);
id_star = (lambda_eff − sqrt(lambda_eff² + 4*(Ld−Lq)²*iq²)) / (2*(Ld−Lq));
% 加调制耦合补偿项
id_star = id_star − K_mod_comp * iq * sin(p_mod*theta);

■ 改进 3：弱磁/扩速（饱和感知）

FMM 低速大转矩易饱和 → 弱磁需更早介入，且考虑调制磁场饱和：

% 调制磁场饱和因子
sat_factor = 1 / (1 + (iq/i_sat)^2);
H_mod_eff = H_mod * sat_factor;

% 弱磁电流（调制感知）
if wr > wr_base
    id_star = −K_fw * (wr − wr_base) * sat_factor;
    % 限制：保证 λ_pm_eff·i_q 不超限
    iq_max = T_max / (1.5*p_eq*lambda_pm_eff_hat);
    iq_star = min(iq_star, iq_max);
end

■ 改进 4：调制耦合前馈补偿

% 在 q 轴电流指令中补偿调制耦合转矩
iq_star = iq_star + K_ff * sin(p_mod*theta) * (TL / (1.5*p_eq*lambda_pm_eff));

Step ⑤ —— 多端口 FMM（双定子）扩展（可选）

若模型是双定子 FMM：

两个逆变器分别控制两套电枢绕组
控制目标：总转矩分配 + 铜损最小 + 热均衡
用 权重分配因子 k ∈ [0,1] 分配 i_q1, i_q2

iq1_star = k * iq_total_star;
iq2_star = (1−k) * iq_total_star;
% k 根据温度/损耗模型动态调整

Step ⑥ —— 测试场景

时间	工况
0~2 s	空载启动至 500 rpm（FMM 低速大转矩区）
2~4 s	负载 0→50 N·m（考验 MTPA）
4~6 s	转速升至 1500 rpm（弱磁区）
6~8 s	切换：传统控制 ↔ FMM‑aware 控制
8~10 s	反转 + 冲击负载

五、典型结果判读（预期数值量级）

5.1 低速大转矩区（500 rpm, 50 N·m）

指标	传统矢量	FMM‑aware
电流幅值	12 A	10.5 A（−12.5%）
铜损	100%	76%
转矩纹波	±3 N·m	±1 N·m
效率	82%	88%

📌 调制感知 MTPA + 耦合补偿 → 同转矩下电流更小、纹波更低。

5.2 弱磁扩速区（1500 rpm）

指标	传统矢量	FMM‑aware
最高稳定转速	1200 rpm	1800 rpm（+50%）
弱磁电流 id	−4 A（过弱磁）	−2.5 A（饱和感知）
高速效率	70%	82%

5.3 动态切换对比（6~8 s）

指标	传统→FMM‑aware 切换后
转矩恢复时间	< 50 ms
转速跌落（负载阶跃）	从 8% → 3%
电流冲击	从 20% → 8%

5.4 调制效应可视化

p_eq_hat 跟踪实际值（误差 < 2%）
λ_pm_eff(θ) 随转角波动（±10%），传统控制忽略 → 转矩纹波来源
T_mod_coupling 占总额定转矩 5~8%，补偿后消失

六、参数调优与常见坑

6.1 关键增益

参数	太小	合适	太大
K_mod_comp	耦合未补偿	0.05~0.2	振荡
K_fw	扩速不足	0.01~0.05	过弱磁→失控
K_ff	动态差	0.1~0.3	超调
MTPA 修正	次优	基于 λ_eff 实时算	数值不稳定

6.2 常见坑

现象	原因	Fix
低速转矩不足	忽略 p_eq 实际值	在线估计或结构参数校准
高速失控	弱磁未考虑调制饱和	加 sat_factor
转矩高频纹波	T_mod_coupling 未补偿	前馈 K_ff·sin(p_mod·θ)
MTPA 点偏移	λ_eff 随 θ 变未更新	每控制周期重算 λ_eff
多端口振荡	分配因子 k 固定	加损耗/热模型动态分配