电机效率图生成 — 电磁场FEM的损耗分离和等效率线生成方法
电机效率图生成的理论基础
什么是效率图
效率图是什么？我听说过这个名词，但不太理解具体是什么…
横轴是转速，纵轴是转矩，用等高线表示各工作点的电机效率。就像地图上用等高线描绘海拔一样，用"效率95%的线""效率90%的线"来描绘。
那么这有什么重要作用呢？
EV在特定走行模式（如WLTP）中会频繁使用某些工作点，效率图可以评估这些频繁使用区域的平均效率。比如高速路段多的话，想让8000~12000 rpm区域效率最高；城市道路多的话，想优化3000~5000 rpm区域。通过效率图，可以根据实际使用场景来优化设计。
也就是说，根据效率图来调整"岛的位置"，对吧？
完全正确。仅有一个理论最高效率点是不够的。要根据实际使用情况最大化"实用效率"——这是电机设计现场最重要的思想。
损耗分离的物理原理
电机有什么损耗？
电机的全部损耗可分为三大类：
- 铁损 $P_{\mathrm{iron}}$：铁心内磁通密度变化引起的损耗。分为滞后损和涡流损
- 铜损 $P_{\mathrm{cu}}$：绕组电气阻抗产生的热损（焦耳损）
- 机械损 $P_{\mathrm{mech}}$：轴承摩擦与转子风损（风阻）
准确分离和计算这三类损耗是生成高品质效率图的关键。
在低转速和高转速时，哪种损耗比较大？
很好的观察。低转速、高转矩时铜损很大。因为要流很大的电流，$I^2R$效应明显。而高转速时铁损急剧增加。铁损与频率（≈转速）成1~2次方关系。也就是说，随转速增加，损耗的主导因素从"铜损"切换到"铁损"。效率图等高线的形状就反映了这两者的相互作用。
铁损与Steinmetz式
铁损的具体计算方法是怎样的？
最广泛使用的是改进的Steinmetz式（iGSE: improved Generalized Steinmetz Equation）。原始Steinmetz式只适用于正弦波磁通，但电机内的磁通波形含有高次谐波，所以需要iGSE。
经典Steinmetz式用单位体积铁损表示为：
各项的物理意义：
- 第1项 $k_h f B_m^{\alpha}$：滞后损。磁区壁移动的能量散逸。与频率 $f$ 一次方比例。$\alpha$ 是材料常数（通常1.6~2.2）
- 第2项 $k_e (f B_m)^2$：古典涡流损。钢板内诱导涡流产生的焦耳热。与频率的2次方比例
- 第3项 $k_a (f B_m)^{1.5}$：异常涡流损（excess loss）。磁区结构动态变化引起。由Bertotti理论解释
高转速时第2项和第3项效应很大，所以要用薄钢板吗？
完全正确。涡流损与钢板厚度 $d$ 的平方成正比。$k_e \propto d^2 / (12 \rho)$（$\rho$ 是钢板电阻率）。所以高速电机用0.2~0.25 mm厚的电磁钢板。要进一步降低铁损，也可选用非晶铁心，但成本会大幅上升。
铜损（I²R损耗）
铜损是绕组抵抗损，用dq轴电流表示：
其中 $R_s$ 是单相绕组电阻，$i_d$, $i_q$ 是dq轴电流分量。
$R_s$ 随温度变化对吧？绕组变热了阻抗就增大…
敏锐的观察。铜的阻抗温度系数约为0.393%/℃。基于20℃：$R_s(T) = R_{s,20}[1 + 0.00393(T - 20)]$。绕组温度升至120℃时，阻抗增加约39%。所以精确的效率图需要电磁-热耦合分析。另外，高频时会出现表皮效应和邻近效应，使有效阻抗增大，这称为AC损。
机械损
机械损不能通过电磁场分析得到吧？怎样把它加到效率图里？
你说得对，机械损不是FEM直接计算的对象。一般用经验公式估算：
$$ P_{\mathrm{mech}} = P_{\mathrm{bearing}} + P_{\mathrm{windage}} = k_b \omega + k_w \omega^3 $$
轴承损失 $k_b \omega$ 与转速一次方比例，风损 $k_w \omega^3$ 与3次方比例。从实测数据反演出 $k_b$、$k_w$，代入效率图中。
效率定式化
各工作点（转矩 $T$、角速度 $\omega$）的效率为：
公式本身很简单呢。关键在于能否精确求出分母里各项损耗吧？
完全同意。公式是"输出/(输出+损耗)"，中学生都能理解。但要精密求出每个工作点的 $P_{\mathrm{cu}}$、$P_{\mathrm{iron}}$、$P_{\mathrm{mech}}$，需要发挥电磁FEM的全部能力。尤其是 $P_{\mathrm{iron}}$ 的精度直接影响效率图的可信度。
特斯拉Model 3的电机为何采用"IPMSynRM"
特斯拉Model 3后电机采用IPM（嵌入式永磁）与SynRM（同步磁阻）的混合结构。通过同时利用磁石扭矩和磁阻扭矩，可将高效率岛扩展到更宽的转速域。高转速时通过弱磁控制（$i_d < 0$）抑制铁损；低转速时用磁石扭矩获得峰值效率。这种"两全其美"的策略，只有深刻理解效率图等高线物理的工程师才能设计出来——这就是为什么读懂效率图是电机设计者的必备技能。
电机效率图生成的数值计算方法
电磁场FEM定式化
生成效率图时，FEM需要求解什么方程？
电机分析的出发点是用矢量位 $\mathbf{A}$ 表示的扩散方程。二维情况下，截面内的磁通密度只由 $z$ 方向的矢量位 $A_z$ 表示：
$$ \nabla \times \left( \frac{1}{\mu} \nabla \times \mathbf{A} \right) + \sigma \frac{\partial \mathbf{A}}{\partial t} = \mathbf{J}_s + \nabla \times \mathbf{M} $$
左边第1项是磁场扩散，第2项是涡流项。右边 $\mathbf{J}_s$ 是线圈电流源，$\mathbf{M}$ 是永磁体的磁化。
时间方向上求解得到一周的磁通密度时间波形吧？
没错。用时间步进法（Crank-Nicolson法或Backward Euler法）离散化，得到每个时刻的磁通密度分布。FEM离散化后的矩阵方程为：
$$ \left( [S] + \frac{1}{\Delta t} [M] \right) \{A\}^{n+1} = \{F\}^{n+1} + \frac{1}{\Delta t} [M] \{A\}^n $$
$[S]$ 是刚度矩阵（透磁率依赖），$[M]$ 是质量矩阵（导电率依赖），$\{F\}$ 是右边向量（电流源+永磁体）。存在磁饱和时，$[S]$ 对 $B$ 有非线性依赖。需用Newton-Raphson法在每个时间步进行迭代求解。
工作点扫描的实现
生成完整的效率图，需要在转矩-转速平面的很多工作点上做分析吧？这计算量不是很大吗？
这正是效率图生成最大的计算成本难点。比如转速10档×转矩10档=100个点的网格，每个点需要一个电气角周期的瞬态分析。二维分析每点30秒的话，100点共50分钟。三维分析每点30分钟，100点就要50小时。
50小时！有没有加速的办法？
有几种加速技巧：
- DOE（试验设计法）+代理模型：用拉丁超立方或最优LHS采样50~100个点，用克里金法、RBF或高斯过程回归做插值。计算点数可减到1/3~1/10
- 频率响应法：正弦波定常假设成立的工作点，用复数频域分析替代瞬态分析
- 并行计算：各工作点相互独立，可尴尬并行处理。100核可同时运行100个点
- 2D+端部修正：二维截面分析加3D端部效应修正系数。成本是3D的1/100，精度达到80~90%
dq轴电流控制与工作点设置
每个工作点的电流条件怎么确定？转速和转矩确定了，电流也就确定了吧？
不是这样的。同一转矩可以用 $i_d$（弱磁方向）和 $i_q$（转矩生成方向）的无数种组合产生。看IPM电机的转矩公式：
$$ T = \frac{3}{2} p \left[ \psi_m i_q + (L_d - L_q) i_d i_q \right] $$
第1项是磁石扭矩，第2项是磁阻扭矩。$p$ 是极对数，$\psi_m$ 是磁石链交磁通，$L_d, L_q$ 是dq轴电感。生成效率图时，对每个转矩-转速组合，用MTPA（每安培最大转矩）控制或最大效率控制确定最优的 $(i_d, i_q)$。
高转速时电压达到上限怎么办？
切换到弱磁（Field Weakening）控制。向负方向调 $i_d$ 减弱磁通，降低反电动势，实现高转速运行。此时效率会下降，但工作范围扩大。效率图必定包含弱磁区，所以要正确建模基础转速和弱磁区边界。
铁损的后处理计算
FEM分析得到磁通密度的时间波形后，怎样计算铁损？
典型步骤是这样的：
- FEM瞬态分析得到各单元 $B_x(t), B_y(t)$ 的时间波形
- FFT分解为基波和各次谐波，得到幅度 $B_k$ 和频率 $f_k$
- 对各谐波分量应用Steinmetz式，求谐波铁损并求和： $$ P_{\mathrm{iron}} = \sum_{k=1}^{N} \left[ k_h f_k B_k^{\alpha} + k_e (f_k B_k)^2 + k_a (f_k B_k)^{1.5} \right] $$
- 对全部单元进行体积积分，得到整个电机的铁损
JMAG和Maxwell内置了这种后处理自动化，能直接显示铁损云图。
非线性收敛和磁饱和处理
电机铁心会磁饱和吧？FEM怎样处理？
电磁钢板的B-H曲线有强非线性。$B$ 在1.5T附近进入饱和，透磁率急剧下降。FEM在每步逐次更新各单元的透磁率 $\mu$（根据前一步的 $B$ 值），用Newton-Raphson迭代求解。收敛判定通常 $||R||/||R_0|| < 10^{-4}$。
高转矩工作点会让齿部和磁铁桥部分进入深饱和，迭代次数增加，计算时间变长。加速技巧包括：
- B-H曲线用光滑插值（Akima插值比样条稳定）
- 导入松弛系数（under-relaxation factor 0.5~0.8）
- 用前步解作初始值（continuation方法）
电机效率图生成的实务应用
分析工作流程
从头到尾生成效率图，要按什么顺序做？
整体流程分5步：
- 第1步：创建电机形状 — 用CAD或工具内部编辑器建立2D截面模型。利用对称性减小到1个极对（比如8极就用45°）
- 第2步：材料设置 — 设定电磁钢板的B-H曲线（非线性）、Steinmetz系数、绕组导电率·温度系数、永磁体剩余磁感应强度 $B_r$ 和矫顽力 $H_{cJ}$
- 第3步：网格生成+电路耦合 — 在气隙和齿尖部用细网格。将FEM与外部电路（逆变器）耦合，实现dq轴电流控制
- 第4步：工作点扫描运行 — 定义转速×转矩的网格点，每点运行瞬态分析（1~2个电气角周期）。JMAG的效率图研究功能可自动化
- 第5步：后处理 — 汇总各工作点的铁损·铜损·机械损，计算 $\eta$。绘制等效率线云图
网格划分策略
电机分析的网格，哪些地方要特别注意？
电机特有的网格要求有几个：
| 区域 | 推荐单元尺寸 | 原因 |
|---|---|---|
| 气隙 | 3~5层以上 | 转矩精度关键。准确捕捉磁通密度的周向分布 |
| 齿尖·槽口 | 钢板厚的1/2以下 | 影响涡流分布和铁损精度 |
| 磁石-铁心边界 | 磁石厚的1/3以下 | 保证减磁分析精度 |
| 铁心背部 | 相对可粗 | 磁通密度较均匀 |
| 外侧空气域 | 粗网格 | 远处磁通密度快速衰减，对精度影响小 |
旋转的转子和固定的定子网格怎样连接？
在气隙中面设置滑动网格（滑移界面）。转子网格旋转时，通过插值与定子网格的耦合得以维持。JMAG的"MotionMesh"或Maxwell的Band（带形区域）就是这个。在气隙中面节点均匀间距能提高插值精度。
WLTP走行模式与设计最优化
怎样用WLTP模式和效率图来优化设计？
从WLTP速度曲线出发，用车辆行驶方程计算每时刻的电机工作点（转矩、转速），然后在效率图上标出这些点的轨迹。这样就能看出"这个走行模式中哪些区域被频繁使用"。
加权平均效率定义为：
$$ \bar{\eta}_{\mathrm{WLTP}} = \frac{\sum_{i} P_{\mathrm{out},i} \cdot \Delta t_i}{\sum_{i} P_{\mathrm{in},i} \cdot \Delta t_i} = \frac{\sum_{i} T_i \omega_i \Delta t_i}{\sum_{i} \frac{T_i \omega_i}{\eta_i} \Delta t_i} $$
通过优化磁石用量、槽口形状、匝数、钢板材质来最大化 $\bar{\eta}_{\mathrm{WLTP}}$，这是实务设计流程。
原来效率图不只是电机性能图表，而是为整个系统最优化服务的工具啊！
完全同意。所以OEM的EV开发部会让电机设计、电池、控制、车辆动力学的工程师围绕效率图进行讨论。"这个走行模式下要加权平均效率≥95%"的需求，就是驱动设计参数最优化的动力。
与实测的对比验证
FEM生成的效率图与实际电机的效率图能匹配吗？
模型建立合理的话，效率绝对值可达±1~2%的精度。但以下因素容易引起偏差：
- 铁损模型不当：加工应力对B-H曲线的衰减（冲片应力效应）未考虑。实测铁损常比解析值高1.3~2倍
- AC损被忽视：高转速时绕组的表皮效应和邻近效应造成的铜损增加未计入
- 温度依赖性：磁石 $B_r$ 的温度系数 −0.1~−0.12%/℃，高温时扭矩明显下降
- 机械损估算不足：密封摩擦或油冷电机的油搅拌损在经验式中没有充分覆盖
效率图"岛"的读法——专业设计师的关注点
效率图上总有"高效率岛"。岛的形状和位置反映设计本质：岛呈细长形说明工作转速范围狭窄，岛呈圆形宽广说明设计稳健。实务中，设计师把WLTP走行曲线叠加在效率图上，调整岛的位置使频繁使用的工作区域进入高效率区——这叫"效率图最优化"。理论最高效率点不一定重要，关键是根据实际使用工况最大化"实用效率"，这是现场设计的最高智慧。
电机效率图生成的软件对比
主要工具的效率图功能对比
什么软件能生成效率图？推荐哪个？
我把主要工具的效率图功能进行对比：
| 功能 | JMAG-Designer | Ansys Maxwell | Motor-CAD | COMSOL |
|---|---|---|---|---|
| 自动效率图生成 | ◎ 专用功能 | ○ 脚本实现 | ◎ 集成功能 | △ 手工构建 |
| 铁损模型（Steinmetz） | ◎ iGSE + Bertotti | ○ Bertotti 3项 | ○ Steinmetz | ○ 自定义可 |
| AC损计算 | ○ | ○ | ◎（解析法） | ○ |
| 电磁-热耦合 | ○ | ◎（Icepak联动） | ◎（内置热分析） | ◎（多物理场） |
| MTPA/弱磁控制 | ◎ 电路耦合 | ○ 外部电路 | ◎ 内置控制 | △ 手工实现 |
| DOE/优化 | ○（外部耦合） | ◎（optiSLang） | ○ | ○ |
| 计算速度（2D） | ◎ | ○ | ◎（解析并用） | △ |
JMAG-Designer工作流程
JMAG的效率图生成步骤具体怎样？
JMAG有"效率图研究"专用功能，步骤如下：
- 建立基础模型：用一个工作点构建并验证2D瞬态分析模型
- 设置效率图研究：指定转速范围（如500~15000 rpm）和转矩范围（如0~300 Nm），设定网格点数
- 选择控制模式：MTPA + 弱磁控制。设定电压限制·电流限制
- 选择铁损模型：从材料库自动获取Steinmetz系数，或手工输入
- 运行：JMAG自动遍历全部网格点，在探索最优 $(i_d, i_q)$ 的同时做分析
- 查看结果：效率图、损耗分布图、电流矢量轨迹自动生成
JMAG的优势是与日本汽车厂商·零部件商长年共同开发形成的电磁钢板材料库很充实。
Ansys Maxwell工作流程
Maxwell没有自动效率图功能，要怎样生成？
Maxwell单独没有JMAG那样的专用功能，但现在的推荐方法是与Ansys Motor-CAD联动。Motor-CAD先用快速解析法生成粗效率图，然后用Maxwell 2D/3D精密验证关键工作点。这种"混合方法"是现在实务的主流。
如果只用Maxwell的话，可用Python脚本（IronPython）自动化参数扫描。设置Ansys Electronics Desktop的"Optimetrics"（DOE），定义各工作点，并行执行解析。
开源替代方案
商用软件太贵了…有免费生成效率图的办法吗？
有几个开源选择：
- FEMM (Finite Element Method Magnetics)：David Meeker开发的免费软件。可做2D静磁场·涡流分析。用Lua/Python脚本可实现参数扫描自动化。虽然需要脚本编程，但学生和研究人员可接受精度
- Elmer FEM：CSC（芬兰科学计算中心）开发的开源。支持3D电磁分析，但缺电机设计专用GUI和效率图功能
- PyMotor / ADEPT：Python电机设计工具包。用FEMM的Python包装自动生成效率图的研究工具
开源的最大弱点是铁损模型精度和材料库不足。需要自己整理Steinmetz系数，这工作量很大。
效率图"自动生成"的时代——DOE与代理模型的组合
过去生成一张效率图要花几十小时，因为要逐个分析转矩×转速的网格点。现在用DOE（试验设计法）选取采样点，再用高斯过程回归等代理模型进行插值——计算点数从10000降到100~300。Ansys OptiSLang或Altair HyperStudy自动化了这个流程，但采样策略失误会漏掉效率的"谷"，风险存在。
电机效率图生成的故障排除
铁损与实测不符
我生成的效率图高转速时与实测差距很大…铁损好像被低估了。
这个问题很常见。原因和对策整理如下：
| 原因 | 说明 | 对策 |
|---|---|---|
| 加工劣化未考虑 | 冲片或激光加工使钢板边界B-H曲线恶化。铁损实际值增加1.3~2倍 | 应用加工影响系数（Building Factor: 1.3~1.8）。JMAG v23+支持加工劣化模型 |
| Steinmetz系数的温度依存性 | 120℃时滞后损比20℃增加20~30% | 用温度相关的Steinmetz系数。多个温度点拟合数据 |
| 旋转磁场分量忽视 | Steinmetz式适用交变磁场。定子齿部会发生圆形·椭圆旋转磁场 | 采用旋转铁损模型（旋转·交变分离型） |
| 直流磁化偏置 | 永磁体DC偏置导致磁滞回线不对称 | 应用小磁滞环模型 |
高效率岛位置偏移
高效率岛相比实测偏向高转速…什么原因？
高效率岛的位置偏移通常源于电感估算精度不足：
- $L_d, L_q$ 的磁饱和依存性：高负荷时齿部·磁铁桥部进入饱和，$L_q$ 大幅下降，磁阻转矩减少。同样转矩需要更大电流，铜损增加，效率下降。检查FEM是否准确捕捉了这种饱和依存性
- 绕组端部漏磁：2D分析忽视3D端部漏磁通。这部分 $L_d$ 被低估，弱磁控制的效果与实机不符
- 磁石温度衰减：NdFeB磁石 $B_r$ 的温度系数 −0.11%/℃。80℃→120℃升温时 $B_r$ 降约4.4%，磁石转矩下降
非线性收敛失败
高转矩工作点分析时Newton-Raphson迭代发散了。该怎么办？
高转矩工作点铁心深度饱和，透磁率剧变，Newton-Raphson难以收敛。试试这些办法：
- 延伸B-H曲线高磁场段：若实测数据仅到 $B = 2.0$ T，高于此需要恰当外推（向真空透磁率渐近）。外推不当会导致微分不连续
- 松弛系数降至0.3~0.5：限制每步的参数更新幅度来抑制振荡
- 时间步长减半：前一步的解质量好坏直接影响后一步的初值，步长大时易发散
- 检查网格质量：特别是磁铁桥部（磁饱和最严酷）是否有过度扭曲的单元
计算时间爆炸
效率图计算已经3天还没完成…要怎样加速？
冷静。先诊断什么在吃时间：
- 网格点太多？ — 初期用转速5档×转矩5档=25点足够。后用代理模型插值
- 用了3D分析？ — 先用2D掌握全貌，再针对关键点做3D
- 过渡过程计算太长？ — 定常到达后的最后1周期就够。初值用前工作点的解
- 输出数据量过大？ — 保存所有单元·所有时刻的磁通密度会拖累磁盘I/O。仅保存后处理需要的数据
- 没有并行计算？ — 各工作点独立，可按核数同时运行
感谢！我先把网格点降到25个，再用2D重新算。
对。效率图生成的哲学是"粗快全局 → 精化关键点"。最初就追求完美会在计算成本上吃亏。