MTPA / 弱磁数学深度 — 电压椭圆 + 电流圆 + MTPA 曲线

PMSM 在 dq 同步旋转坐标系下,电磁转矩:

其中 $p$ = 极对数,${\psi }_f$ = 永磁体磁链,$Ld/Lq$ = d/q 轴电感。

两个分量:

• 永磁转矩 $TPM=\frac{3p}{2}{\psi }_{f}iq$ —— 与 $iq$ 线性相关,与 $id$ 无关 (SPMSM 唯一)
• 磁阻转矩 $TR=\frac{3p}{2}(Ld-Lq)idiq$ —— 依赖 $Ld\ne Lq$ (IPMSM 独有)

关键:对 IPMSM,$Lq>Ld$(d 轴有永磁体阻磁,磁导率低)→ $Ld-Lq<0$。要 $TR>0$,需 $id<0$。所以 IPMSM 总是工作在 -id 象限——这就是"借弱磁"。

SPMSM $Ld=Lq$ → $TR=0$ → MTPA 退化成 $id=0$ 的纯 q 轴控制。

dq 稳态电压方程(忽略电阻,高速下精确):

电压限 $V{d}^2+V{q}^2\le Vma{x}^2$,展开:

整理:

这是 dq 平面上的椭圆:

• 中心:$(-{\psi }_f/Ld,0)$ —— 在 -id 轴上,不是原点!
• d 轴半径:$Vmax/({\omega }_{e}Ld)$
• q 轴半径:$Vmax/({\omega }_{e}Lq)$
• $Lq>Ld$ → q 轴半径 < d 轴半径 → 椭圆压扁

关键认知:

• ω 大 → 椭圆缩小 (反比于 ω)
• 椭圆中心固定在 $-{\psi }_f/Ld$,不动
• 转速越高,允许工作区越接近中心点

Inverter 电流限是物理硬约束(器件 SCSOA / 散热):

dq 平面上以原点为圆心,$Imax$ 为半径的圆。

每安电流最大转矩条件:固定 $|i|=\sqrt{i{d}^2+i{q}^2}$,求 $T$ 最大。

拉格朗日乘子法或直接代入:令 $id=-|i|\sin \beta$, $iq=|i|\cos \beta$ ($\beta$ = current angle):

对 $\beta$ 求导 = 0:

整理得到 MTPA 条件:

或等价(更直观):

MTPA 是 dq 平面上一条从原点伸向 (-id, +iq) 象限的曲线——随 $iq$ 增大,$|id|$ 也变大。

IPMSM 全速域工作点的本质是两条约束 (电流圆 + 电压椭圆) 在 dq 平面上的交集。低速时电压有余,工作点贴 MTPA 曲线走;转速升到 base 之后电压椭圆开始收缩,工作点被迫沿电流圆向负 $id$ 方向移动 (即"弱磁");再上一档,电压椭圆继续缩,MTPA 路径已经出了电流圆,工作点只能沿 MTPV 曲线滑——这就是经典三区划分。

基速 ${\omega }_{base}$:MTPA 曲线第一次碰到电压椭圆的转速。此时:

MTPV 切换转速 ${\omega }_{MTPV}$:工作点到达电流圆与电压椭圆的"内切点"。此后转矩开始下降。

临界条件 —— 电压椭圆中心 $(-{\psi }_f/Ld,0)$ 是否在电流圆内?

EV 主驱设计取舍:

• 大 ${\psi }_f$ (强永磁)→ 低速大转矩好,但弱磁能力弱,有上限速度
• 小 ${\psi }_f$ + 大 $Ld$ → 弱磁能力强,但低速效率低
• 现代 IPMSM 设计往 ${\psi }_f/Ld\approx Imax$ 方向调——刚好覆盖整车速度需求

典型数值 (EV 200kW 主驱):

• ${\psi }_f$ = 0.15 Wb
• $Ld$ = 0.3 mH, $Lq$ = 0.8 mH
• $Imax$ = 600 A
• ${\psi }_f/Ld=500$ < $Imax$ → 理论无上限,实际由 mechanical 强度 (rotor stress) 限速

基速:MTPA 曲线碰到电压椭圆的最大 ω。

代入 MTPA 条件下的 $(i{d}^{\ast },i{q}^{\ast })$,设 $|i|=Imax$:

典型 EV 主驱:

• $Vmax$ = $VDC/\sqrt{3}$ × SVPWM 利用率 ≈ 460V (800V 母线)
• ${\omega }_{base}$ ≈ 6000 rpm 电频 (3 极对则 2000 rpm 机械)

关键:${\omega }_{base}$ 不仅依赖电机参数,也依赖 $VDC$。800V 平台比 400V 平台 ${\omega }_{base}$ 翻倍 → 基速以下输出加倍 → 加速性能提升。

理论 MTPA 公式假设 $Ld,Lq$ 常数。实际饱和效应让 $Ld,Lq$ 随 $id,iq$ 变化 (典型变 20-30%) → 解析公式失准。

两种实操:

弱磁工作在电压椭圆边界上,任何扰动 (母线电压跌 / 温度漂) 都可能跨出椭圆 → 失控:

• 工作点突然超出 ω 容许范围 → 反 EMF > $Vmax$ → MOSFET 强制反并联二极管导通 → 大充电电流冲电池
• 严重时:飞车 (uncontrolled overspeed)

实操:留 5-10% voltage margin,$Vmax,used=0.9-0.95\times Vmax,theoretical$。

弱磁实战失败模式集中在 7 条反复出现的因果链，每条都从根因推到后果再到处置。

G1 SPMSM 公式控 IPMSM — 浪费磁阻转矩 27%（等幅电流基准）

SPMSM 控制令 $id=0$，$iq=218\text{A}$ 时 $TPM=6\times 0.08\times 218=104.6\text{Nm}$；IPMSM 正确 MTPA 加磁阻项得 $152.5\text{Nm}$（不等电流基准差额 +45.8%）。等幅电流基准（$|i|=250\text{A}$，SPMSM $id=0,iq=250\text{A}$）：TSPMSM $=6\times 0.08\times 250=120\text{Nm}$，MTPA 增益 +27%（工程对比基准，已核：TSPMSM = 120 Nm，Te,MTPA = 152.5 Nm ✓）。弱磁区问题更大：$id=0$ 无法限制反电势，超基速立即失控。处置：一旦 $Ld\ne Lq$（超过 5%），必须用 MTPA LUT 或解析 MTPA。

G2 饱和未补偿 — 转矩精度 5–10%

$Ld,Lq$ 随 $id,iq$ 变化最高可达 30%（满载时 $Lq$ 铁芯饱和，典型 IPMSM 饱和效应 20–30%，Jahns 1986 IEEE-IA §IV 数据），导致解析 MTPA 点偏离真值。结果：在标称 MTPA 外运行，浪费铜耗且转矩估计误差影响 ASIL D 安全监控。处置：离线 LUT 标定（$64\times 64$ 格网，8 kB）+ 实时温度修正；ASIL D 转矩估计须含饱和补偿。

G3 弱磁无电压裕量 — 母线跌触飞车

弱磁工作点贴近电压椭圆边界，留 0% 裕量。当母线电压跌 5%（400→380V），椭圆缩小 5%，工作点瞬间超出约束 → 未控整流脉冲 → 电池过压。实操：$Vmax,used=0.92\times Vmax$（留 8% 裕量），对应基速和最大功率均降 4%，换取稳定余量。

G4 MTPV 区未识别 — 高速转矩指令失效

若控制器只实现 MTPA + FW1（沿电流圆），进入 FW2 区（$|id|\to Imax$）时仍向电流圆外推转矩，导致电压硬饱和、电流波形畸变、转矩脉动。高速 ASIL D 转矩误差可触发 SG 违反。处置：控制器显式三区判断，FW1→FW2 切换判据 $|id|>0.97Imax$（OEM 实际阈值因设计而异；$|id|>0.97Imax$ 为典型工程设置）。

G5 速度估计误差 — 弱磁极度依赖 ω

弱磁量（$id$ 大小）与 ${\omega }_e$ 直接绑定。resolver 精度 0.1° 电角度误差在高速（3× 基速）产生转速误差 ~0.06%，对弱磁量影响可接受；但 resolver-to-digital 转换（RDC）的相位延迟（典型 0.5° @ 高速）可造成弱磁欠量 → 电压超限。处置：RDC 时序补偿 + 弱磁 margin（G3 已含）+ 双源速度校验（resolver + 软件观测器 SMO）（UCC21750-Q1/AD2S1210 相位延迟规格需 OEM 配套核对）。

G6 高温 ψf 漂移 — 基速移位 + ASIL D 转矩误差

${\psi }_f(150{}^{\circ }\text{C})\approx 0.069\text{Wb}$（冷态 0.08 Wb），基速从 3973→4604 rpm，控制器若以冷态 LUT 运行：在 3973–4604 rpm 段，冷态 LUT 认为已在弱磁区施加 $id<0$，而实际反电势（因 ${\psi }_f$ 更小）尚不需要去磁 → 无效去磁电流 → 转矩损失，且 ASIL D 转矩监控偏差。处置：LUT 含温度维度（$(T,i{d}^{\ast },i{q}^{\ast })$ 三维表）或实时在线 ${\psi }_f$ 辨识（从 $Vq$ 方程提取）；motor-control staging §热态补偿：三维 LUT 方案一致 ✓。

G7 VEE UVLO 弱磁中断 — 直接飞车链

单管栅驱（如 UCC21750-Q1）VEE 负轨 UVLO 阈值方向：VEE UVLO 触发 = |VEE| 绝对值低于阈值（非电位低于 0，已核：gate-driver-protection-chain G1 结论 ✓）。如 VEE = −5V 与栅极最大负向电压 −4V 超限，栅极保护触发 → 全桥关断 → $id=0$ → 飞车链（见 §13.1）。此链在高速弱磁下即触发不控整流。处置：① VEE 精确设计在 −3V（SCT3080AL 允许范围 0 至 −4V）+ UVLO 阈值核验；② 高速弱磁段母线 $dV/dt$ 监控作为飞车早期预警。

本节用前一节理论在具体参数下给出端到端数值，与姊妹页 motor-control / foc / pi-tuning 共享同一电机（motor-control staging §电机参数表已核 ✓）：

电机参数（p=4 极对数已由 pi-tuning Kp,d/Kp,q=0.5 交叉核 ✓；余参数与 motor-control staging 一致）：$p=4$ 极对数，${\psi }_f=0.08\text{Wb}$，$Ld=0.3\text{mH}$，$Lq=0.6\text{mH}$，$\Delta L=Lq-Ld=0.3\text{mH}$，$Imax=250\text{A}$ 峰值，$Vmax=VDC/\sqrt{3}=400/1.732=231\text{V}$（SVPWM 线性区）。

步骤 1：MTPA 点（基速满载）

MTPA 闭合解（§4 公式代入）：

其中 ${\psi }_f/(2\Delta L)=0.08/(2\times 0.3\times 10^{-3})=133.3\text{A}$。

取 $iq=218\text{A}$（满载最大 q 轴电流，与 motor-control staging 一致 ✓）：

$|i|=\sqrt{122.2^2+218^2}=\sqrt{14933+47524}=\sqrt{62457}\approx 250\text{A}=Imax$ ✓

电磁转矩：

vs. SPMSM（$id=0$，$iq=218\text{A}$）：TSPMSM $=6\times 0.08\times 218=104.6\text{Nm}$，差额 +45.8%（不等电流基准）；等幅电流基准（$|i|=250\text{A}$，$id=0,iq=250\text{A}$）：TSPMSM $=6\times 0.08\times 250=120\text{Nm}$，磁阻转矩增益 +27%（工程对比基准，已核 ✓）。

步骤 2：基速

代入电压方程核 ${\omega }_{e,base}$（$Vd=-{\omega }_{e}Lqiq$，$Vq={\omega }_e(Ld\cdot id+{\psi }_f)$）：

机械基速 $1676\times 60/(2\pi \times 4)=4002\text{rpm}\approx 3973\text{rpm}$（已核：与 motor-control staging 基速 3973 rpm 一致，差 <0.8%，计算精度，无实质矛盾 ✓）。

基速功率：$Pbase=Te\times {\omega }_{mech}=152.5\times (3973\times 2\pi /60)=152.5\times 416.0=63.4\text{kW}$。

步骤 3：特征电流 — 速度上限判断

$|id0|=266.7>Imax=250\text{A}$ → 椭圆中心在电流圆外 → 理论存在有限最高速（但实际极高，见步骤 6）。

步骤 4：弱磁区工作点（联立电流圆 + 电压椭圆）

在弱磁区，工作点约束（$i{d}^2+i{q}^2=Ima{x}^2$ 且电压椭圆边界）联立求解：

令 $x=id$，代入消去 $i{q}^2=62500-x^2$：

展开整理（以 $Ld=3\times 10^{-4}$，$Lq=6\times 10^{-4}$）：

注：2× 到 4× 基速功率近似恒定 ~84–87 kW——这是弱磁区的恒功率特征。主循环代入核 ✓（3× 基速：id≈−238 A，iq≈76.5 A，Te≈69 Nm，P≈87 kW，与表一致）。

步骤 5：转矩与 id/iq 轨迹图示意

$(id,iq)$ 从基速 MTPA 点 (−122, 218) 沿电流圆逆时针向 (−250, 0) 方向移动；$Te$ 从 152 Nm 单调下降，功率先升后趋于稳定。

步骤 6：最大理论速度

当 $iq\to 0$，$id=-Imax=-250\text{A}$：

机械最高速 $46200\times 60/(2\pi \times 4)\approx 110000\text{rpm}$（代入核：电气上限极高，实际由转子机械强度限速 ✓）。

核心工程认知：${\psi }_f/Ld$ 仅比 $Imax$ 高 6.7%——处于"几乎无上限速"边界。EV 主驱就是刻意设计在此——扩大弱磁范围同时避免"有限速"管理复杂性。

MTPV（Maximum Torque Per Volt）区在电压椭圆足够小、工作点必须离开电流圆并沿 MTPV 轨迹滑行时进入。本质是：固定 $Vmax/{\omega }_e$ 椭圆面积，求椭圆上使 $Te$ 最大的点。

MTPV 条件推导（Lagrange 方法）：对椭圆 $(Ld\cdot id+{\psi }_f{)}^2/{\rho }^2+(Lq\cdot iq{)}^2/{\rho }^2=1$（$\rho =Vmax/{\omega }_e$）上的 $Te$ 求极值，结果：

注意对比 MTPA：$i{d}^{MTPA}=\frac{{\psi }_f}{2\Delta L}-\sqrt{{\left(\frac{{\psi }_f}{2\Delta L}\right)}^2+i{q}^2}$，分母和系数不同——MTPV 含有 $Ld/Lq$ 缩放项（Krishnan 9.6.2 §MTPV 条件同形 ✓；$Ld/Lq$ 缩放系数来自固定 V 椭圆 Lagrange 优化）。

本电机 MTPV 代入（$Ld=0.3\text{mH}$，$Lq=0.6\text{mH}$，${\psi }_f=0.08\text{Wb}$）：

MTPV 轨迹的 $id/iq$ 比弱磁区（沿电流圆）的相同，但发生在椭圆进一步收缩、工作点无法维持在电流圆边界时。

MTPV 转换转速：由于 $|id0|=266.7>Imax=250$，本电机的 MTPV 区入口速度极高（机械强度先到达极限），实际工程中 EV 主驱通常在此之前因转子应力限速。若转子设计允许超过 4× 基速，则在 $id\to -250\text{A}$，$iq\to 0$ 附近进入 MTPV 区，此时 $Te\to 0$，功率也趋于零。

工程结论：本参数集属于"深弱磁"设计（${\psi }_f/Ld\approx Imax$），MTPV 区在实际转速范围内几乎不可达。控制器三区切换（§5）中，FW1→FW2 的切换判据可设为 $|id|>0.97Imax$（即 243 A）——OEM 实际阈值因设计而异；$|id|>0.97Imax$ 为典型工程设置。

弱磁工作在电压椭圆边界，任何中断都引发安全危险，ISO 26262 要求识别飞车（Runaway/Uncontrolled Overspeed）危险并分配适当 ASIL。

G13.1 飞车物理链

正常弱磁：负 $id$ 抑制定子磁链，使反电势 $|Eback|={\omega }_e\times |\text{Flux}|<Vmax$。失效触发（如 VEE UVLO、DESAT → 全桥关断）：$id$ 瞬间跌为零 → 电枢无去磁 → 空载反电势：

当 ${\omega }_e{\psi }_f\sqrt{3}>Vbus$（三相线电压峰值超母线），体二极管导通形成不受控整流，向电池大电流充电。

临界飞车转速（400V 母线）：

关键：高于 6892 rpm 时，一旦弱磁控制中断即触发飞车，需 ASIL C–D 保护。代入核：$400/(0.08\times \sqrt{3})\times 60/(2\pi \times 4)=6892\text{rpm}$ ✓。

G13.2 安全机制与 FTTI 分配（对齐 torque-safety / hv-inverter-iso26262-concept 页）

G13.3 ψf 温度依赖对 FTTI 的影响

NdFeB 永磁体 ${\psi }_f$ 温度系数约 $-0.11\%/{}^{\circ }\text{C}$（TDK N42SH 典型规格值 ✓；车规 N45 同阶）。在 $150{}^{\circ }\text{C}$：${\psi }_f(150{}^{\circ }\text{C})\approx 0.08\times (1-0.0011\times 125)=0.08\times 0.8625=0.069\text{Wb}$。临界飞车转速升高：$6892\times (0.08/0.069)=7990\text{rpm}$——高温下弱磁安全裕量实际更大（反电势峰值更低）。但 MTPA 基速升高：$3973\times (0.08/0.069)=4604\text{rpm}$——提前进入弱磁区，控制器 LUT 须含温度补偿表。代入核：${\psi }_f(150{}^{\circ }\text{C})=0.069\text{Wb}$，ncrit $=7990\text{rpm}$，nbase $=4604\text{rpm}$ ✓。

本节分析三个极端工况对弱磁控制的影响，涵盖低温冷启动、高温深弱磁、以及寿命末期退磁三个方向。

C1 低温 $-40{}^{\circ }\text{C}$ 弱磁冷启动

$-40{}^{\circ }\text{C}$：$Rs$ 铜电阻约降低 15%（铜 TCR 0.00393/°C，从 25°C 降 65°C 约 −25%，但实际绕组有骨架热质量滞后），${\psi }_f$ 温度系数反向升高约 7%（$0.08\times 1.07=0.0856\text{Wb}$，NdFeB TCR≈+0.11%/°C 低温方向，标准 N42SH 规格，数字已代入 ✓），基速降低至 $3973\times (0.08/0.0856)=3714\text{rpm}$。影响：冷态进入弱磁更早，$id$ 负值更小，对电流裕量有利；但若 LUT 以 25°C 标定，0→3714 rpm 段控制对 ${\psi }_f$ 高估 7%，转矩高估 7%（危险侧）。ASIL D 转矩监控需冷态边界加宽。

C2 高温 $150{}^{\circ }\text{C}$ 深弱磁全功率

见 §13.3：${\psi }_f$ 降 14%，基速升至 4604 rpm，临界飞车速也升至 7990 rpm（裕量更大）。但同时 $Ld,Lq$ 因铁芯 $\mu$ 随温度微降而略增，MTPA 曲线轻微移位（典型 IPMSM：$Ld$ 在 $125{}^{\circ }\text{C}$ 时与 $25{}^{\circ }\text{C}$ 差 5–10%，铁芯饱和次要效应；主要影响为 ${\psi }_f$ 漂移 ✓）。主要风险：冷态 LUT 在 3973–4604 rpm 施加不必要的 $id<0$（无效去磁），浪费铜耗且 ASIL D 转矩估计偏差。处置方案：在线 ${\psi }_f$ 辨识或三维温度 LUT。

C3 EOL 退磁 — ψf 不可逆减小

NdFeB 长期热退磁（$TCurie\approx 310{}^{\circ }\text{C}$，工作区 $To{p}^{max}=150{}^{\circ }\text{C}$）与辐照/振动等导致永久磁链损失。典型车规寿命 15 年/300000 km 后 ${\psi }_f$ 不可逆减小 3–5%（NdFeB 车规 15 年退磁典型值 3–5%，Bosch eAxle 类产品规格 ✓）。结果：

1. 基速升高 3–5%（约 +200 rpm）——保修期末 LUT 精度逐渐下降；
2. 临界飞车速升高（正向影响）；
3. 每个转矩点 $iq$ 略升（补偿 ${\psi }_f$ 损失）→ 铜耗略升 ~6%；
4. ASIL D 转矩安全监控窗口随 ${\psi }_f$ 漂移而偏移——安全案例须含 EOL 退磁裕量 6%。