永久磁石电动机的不可逆减磁分析
永久磁石电动机减磁的理论基础
不可逆减磁的定义
老师，"减磁"一旦发生就无法恢复了吗？
好问题。减磁分为两种："可逆减磁"和"不可逆减磁"。可逆减磁在温度降低后会恢复。问题在于不可逆减磁。当B-H曲线的工作点落到knee point（拐点）以下时，永久磁石的磁力会永久性地降低。即使消除外部原因也无法恢复。
永久性！这是怎样的机制呢？
简单说，永久磁石内部有"磁区"。在正常状态下，这些磁区方向一致，产生强磁力。但在强反磁场或高温条件下，部分磁区会反向翻转。当压力超过knee point时，翻转的磁区无法回到原来的方向——这就是不可逆减磁。
实际的电动机中什么时候会发生呢？
最危险的是高温和大电流的组合。比如，电动车爬坡（大扭矩=大电流）或逆变器短路故障时。线圈铜耗会使磁石温度上升到180°C或200°C，矫顽力大幅下降，工作点迅速跌到knee point以下。
B-H曲线与knee point
knee point在B-H曲线上具体是什么位置？能给个公式吗？
永久磁石的减磁曲线（第二象限的B-H曲线）通常直线下降，但在某个磁场强度处会急剧弯曲，磁通密度迅速下降。这个"转折点"就是knee point。数学上：
其中 $B_r(T)$ 是温度 $T$ 下的剩余磁通密度，$\mu_{\text{rec}}$ 是回线透磁率，$H_k(T)$ 是knee point磁场强度。不可逆减磁的判定条件是：
其中 $H_{\text{op}}$ 是磁石工作点的磁场强度（包含反磁场）。
工作点跌到knee point以下就不行了。但knee point本身也随温度变化吧？
完全正确。温度升高时knee point会向右上移动——也就是说，高温下用更弱的反磁场就能引发不可逆减磁。磁石在高温下变得"脆弱"。这就是为什么减磁分析中温度依赖的B-H曲线至关重要。
温度系数和磁石等级
温度变化对磁石特性影响有多大？想看数字。
NdFeB（钕铁硼）磁石的温度依赖关系由以下公式表示：
代表性的温度系数是这样的：
| 磁石种类 | $\alpha_B$ (%/°C) | $\alpha_H$ (%/°C) | 最高使用温度 |
|---|---|---|---|
| NdFeB (N级) | -0.10 〜 -0.12 | -0.55 〜 -0.65 | 80°C |
| NdFeB (SH级) | -0.10 〜 -0.12 | -0.50 〜 -0.58 | 150°C |
| NdFeB (UH/EH级) | -0.10 〜 -0.12 | -0.45 〜 -0.55 | 180〜200°C |
| SmCo | -0.03 〜 -0.04 | -0.15 〜 -0.30 | 300°C |
| 铁氧体 | -0.18 〜 -0.20 | +0.30 〜 +0.40 | 250°C |
NdFeB的 $\alpha_H$ 是-0.5〜-0.6%/°C，这意味着升温100°C矫顽力就降低50〜60%！？
是的。这就是为什么电动汽车电动机的设计中，如果磁石温度会达到150°C，必须选用SH级或UH级的高矫顽力磁石。用便宜的N级磁石，在夏天爬坡时就会减磁，性能大幅下降。
铁氧体的 $\alpha_H$ 是正的？温度升高矫顽力反而增加？
有意思吧。铁氧体在高温下变强，但低温下变弱。所以铁氧体电动机的减磁风险在冬天的冷启动。NdFeB的情况相反。这是材料学的本质差异。
减磁裕度的定义
knee point到底在哪里确定了，该怎样设计来避免减磁呢？
这就是减磁裕度（Demagnetization Margin: DM）的定义。用它来定量化：
$H_{\text{op}}$ 是磁石内部最小磁场强度（最严酷的位置），$H_k(T)$ 是该温度下的knee point磁场。DM为正时安全，为负时发生减磁。
工程实践中的标准是这样的：
| DM值 | 评价 | 应用领域 |
|---|---|---|
| DM > 30% | 充足余裕 | 消费品，可靠性要求低 |
| DM 20〜30% | 标准余裕 | 汽车（常规工况） |
| DM 10〜20% | 余裕不足 | 需要注意，必须灵敏度分析 |
| DM < 10% | 危险区域 | 应改变磁石等级 |
汽车一般20%以上。但同一个条件设置下，结果变化很大吧…
正是这点难度所在。温度设为多少、电流多大、磁石公差怎么考虑，结果就天差地别。所以减磁分析的核心是"最坏工况的积累"。
实际运行中引发减磁的条件
电动机的哪些运行工况最危险？
危险度从高到低排列：
- 逆变器短路故障：最大电流的5〜10倍过电流流动，d轴方向产生巨大反磁场。最严酷工况。
- 低速大扭矩连续运行：电动汽车爬坡、牵引。铜耗大，磁石温度急速上升。
- 弱磁控制（Field Weakening）：高速转速时在d轴注入负方向电流来削弱磁通。反磁场直接增大。
- 铁氧体低温环境：-40°C冷启动。
弱磁控制也危险…高速巡航时磁石温度也高…
对。高温+弱磁的双重打击。实际设计中，要在最高速度×连续运行×最高冷媒温度的条件下验证减磁裕度。
磁石的"记忆"和"遗忘"——磁区壁的钉扎效应
永久磁石的磁化是通过内部磁区（域）的方向一致性来维持的。当外加反磁场时，磁区壁开始移动，部分磁区反向。当压力超过knee point时，磁区壁会突破结晶界面或析出物的钉扎位置，无法回到原点。这是不可逆减磁的微观机制。高保磁力NdFeB磁石（SH/UH/EH级）通过添加Dy（镝）或Tb（铽）来增强磁区壁钉扎力。稀土元素价格波动直接影响电动机设计成本。
永久磁石电动机减磁的数值计算方法
电磁场FEM的公式推导
减磁分析具体是用什么方程在计算机上求解呢？
基础是磁矢量势 $\mathbf{A}$ 的麦克斯韦方程。二维电动机截面分析的支配方程是：
其中 $\nu = 1/\mu$ 是磁气阻系数，$\mathbf{J}_s$ 是施加电流密度，$\sigma \frac{\partial \mathbf{A}}{\partial t}$ 是涡电流项，$\mathbf{M}_r$ 是永久磁石的剩余磁化矢量。
$\mathbf{M}_r$ 项代表永久磁石。减磁发生时这个值会变化？
正确。减磁后 $\mathbf{M}_r$ 会减小。但重要的是，磁石内部不是均匀减磁，而是发生"部分减磁"。磁石的边角处反磁场集中，最先从那里开始减磁。
FEM计算时，按单元检查磁通密度B和磁场强度H的工作点是否跌到knee point以下。
非线性B-H曲线的处理
非线性B-H曲线在数值上怎么处理？
用牛顿-拉夫逊法（NR法）迭代计算。每次迭代更新透磁率 $\mu$。磁石区域的本构关系为：
其中 $\mu_{\text{rec}}$ 是回线透磁率（通常1.02〜1.10）。knee point以下时：
$\mu_{\text{demag}}$ 是knee point以下的非线性透磁率，从B-H曲线的斜率通过插值求得。
knee point前后模型转换，收敛肯定很困难…
是的。knee point处B-H曲线斜率变化剧烈，NR法容易收敛不好。对策是：
- 在knee point附近用样条插值使B-H曲线光滑连接
- 使用阻尼（低松弛因子：0.3〜0.7）
- 收敛判定设置为 $\|R\|/\|R_0\| < 10^{-4}$
- 分小步加载磁场
热-电磁耦合分析
温度特别重要，只做电磁分析够吗？必须加热分析？
准确的减磁评估需要热-电磁耦合分析。流程是这样的：
- 电磁场分析计算损耗分布（铁损+铜损+磁石涡电流损）
- 损耗作为热分析的热源，计算温度分布
- 根据温度分布在各单元更新磁石B-H曲线
- 用更新的B-H曲线重新做电磁分析 → 回到第1步
反复迭代到温度收敛。通常2〜5次迭代就收敛了。
计算时间想必很长…
确实。实务上经常用"单向耦合"的偷懒做法。先定常热分析预测磁石最高温度，用那个温度的B-H曲线做一次电磁分析。严格的双向耦合只在最后确认阶段做。
部分减磁建模
你之前说的"部分减磁"怎样在模型中表现？
部分减磁建模主要有两种方法：
- 单元剩余磁化更新法：对每个有限单元判定工作点，跌到knee point以下的单元，$B_r$ 值减小。JMAG的"减磁分析"功能就是这样。
- 回线追踪法：减磁后的磁石沿着回线（斜率 $\mu_0 \mu_{\text{rec}}$）回缩。新回线的截距就是减磁后的 $B_r'$。
减磁后的电动机特性要正确评估，必须拿更新的 $B_r'$ 再跑一遍电磁分析。减磁前后的扭矩纹波和齿槽扭矩波形会改变，影响NVH。
永久磁石电动机减磁的工程应用
分析流程
怎样从头开始做减磁分析？
典型的减磁分析步骤是这样的：
- 建立电动机模型：定子和转子的二维截面模型。利用旋转周期对称性，只建1个极距（比如8极的话45°）。
- 输入材料数据：磁石B-H曲线从20°C到200°C，每隔20°C一条。钢板的B-H曲线和铁损曲线。
- 设置工作条件：转速、电流幅值、相位角（d-q轴电流）、通电方式。
- 设置温度：磁石温度。单向耦合的话直接输入预计的最高温度；双向耦合的话连接热分析。
- 运行电磁分析：过渡分析至少计算1个电气周期（比如60°）。
- 判定减磁：提取磁石各单元的工作点，与knee point对比。
- 部分减磁重算（如需）：更新减磁单元的 $B_r$，重新分析。评估扭矩降低率。
"每个时刻"都要看？不是只看最坏时刻就行？
好问题。转子位置不同，反磁场大小变化，所以"最严酷时刻"需要找。通常d轴电流最大时或3相短路初期最危险。但磁石内部不同位置的最坏时刻可能错开，所以全时步判定最安全。
网格策略
减磁分析的网格有什么要点吗？
减磁分析的网格至关重要。网格质量决定部分减磁的分布能否正确捕捉：
| 区域 | 推荐单元尺寸 | 原因 |
|---|---|---|
| 磁石（整体） | 0.3〜0.5 mm | 准确捕捉部分减磁分布 |
| 磁石边角 | 0.1〜0.2 mm | 反磁场集中，减磁最先发生 |
| 气隙 | 间隙长度的1/3〜1/5 | 直接影响扭矩精度 |
| 定子齿尖 | 0.3〜0.5 mm | 抓磁饱和 |
| 回铁 | 1.0〜2.0 mm | 可粗些（磁通相对均匀） |
磁石网格跟结构分析应力集中部分一个感觉。
比喻很好。减磁的"磁石角部反磁场集中"就像结构的"应力集中"。网格太粗会漏掉减磁。至少磁石厚度方向要5层以上单元。
最坏工况条件设置
"最坏工况"具体怎样定义？
汽车OEM的典型最坏工况设置是这样的：
- 磁石温度：最高使用温度（比如180°C）。部分OEM还加10〜20°C冗余。
- 电流：最大电流（定格的1.5〜3倍）。考虑短路故障的话5〜10倍。
- 磁石公差：$B_r$ 下限公差（-3〜-5%），$H_{cj}$ 下限公差（-5〜-10%）。
- 磁石厚度：尺寸公差下限（厚度越小反磁场越大）。
全部条件同时最坏会不会过度设计？
概率上讲全最坏同时出现的机率很低。现在有些公司用"Six Sigma统计设计法"，考虑各参数的分布。但安全性相关功能，保守做法（全积累）的OEM还是多数。
结果评价和判定标准
分析结果出来了，怎样判定合格？
判定标准分三层：
- 减磁率：磁石整体平均 $B_r$ 降低比例。通常≤5%可以接受。
- 减磁面积比：跌到knee point以下的单元占面积的比例。≤10%是目标。
- 扭矩降低率：减磁前后扭矩的降低。≤2〜3%目标。
最终判定用扭矩降低率最直观。磁石角部局部减磁，如果对扭矩影响小，就不算问题。
电动汽车"坡道试验"——减磁分析救了一款车
某日本OEM开发电动汽车电动机样机，试验中发现"长坡低速登坡后，平地加速变肉"。原因是磁石部分减磁。分析团队用减磁模拟重现了工况，发现磁石角部2mm范围的区域跌到knee point以下。解决方案：从N42SH升级到N42UH，并在磁石角部加C0.5倒角。磁石成本增加约15%，但避免了召回——这是减磁分析的价值体现。
永久磁石电动机减磁的软件比较
JMAG-Designer
减磁分析最强的软件是什么？先讲JMAG。
JMAG是日本JSOL公司开发，专攻电机设计。减磁分析方面：
- 减磁专用功能：GUI里直接选"减磁分析"，自动执行单元knee point判定和部分减磁重算。
- 磁石材料库：与信越化学、TDK、日立金属等直接合作，内置温度依赖B-H曲线。
- 减磁分布图：彩色等高线展示减磁率分布，一目了然。
- 参数扫描：温度、电流变量自动扫描，生成减磁裕度地图。
Ansys Maxwell
Ansys Maxwell怎样？
Ansys Maxwell是全球最大CAE公司的电磁求解器。减磁分析方面：
- Icepak热耦合：同一Workbench平台内Maxwell和Icepak无缝耦合，包括冷却液流场的详细温度。
- Demagnetization Module：温度依赖B-H曲线设置，自动部分减磁判定更新。
- 脚本自动化：IronPython批量参数扫描。
- 3D分析：支持端部效应、磁石斜切等复杂几何。
Altair Flux
Flux（前Cedrat）是否在减磁上有特色？
Altair Flux来自法国，在学术严谨性上口碑好：
- 部分减磁模型精度：各单元独立追踪回线，算法精细。
- FEM-BEM耦合：用边界元法处理开放域，空气网格无需细分。
- FluxMotor：电动机专用前处理器，模板快速建模。
功能比较矩阵
一览表对比一下？
| 功能 | JMAG | Maxwell | Flux | COMSOL |
|---|---|---|---|---|
| 减磁专用GUI | ◎ | ○ | ○ | △（手工设置） |
| 部分减磁模型 | ◎ | ○ | ◎ | △ |
| 磁石厂商DB | ◎（国内充实） | ○ | ○ | △ |
| 热耦合 | ○ | ◎（Icepak） | ○ | ◎（内置） |
| 3D支持 | ○ | ◎ | ○ | ◎ |
| 脚本自动化 | ○ | ◎ | ○ | ◎（Java API） |
| 中文支持 | ◎ | ○ | △ | ○ |
| 价格带 | 中 | 高 | 中 | 中〜高 |
日本电机厂商JMAG多，海外Maxwell多的感觉？
完全同意。日本汽车OEM/Tier1选JMAG压倒性多，磁石厂商的材料数据联动是关键。欧美用Ansys Maxwell、MotorCAD（Speed/Motor-CAD）。中国JMAG和Maxwell齐驾并驱。最后的选型标准往往是"我想用的磁石材料这个软件里有没有"。
减磁分析的"材料数据"才是关键
减磁分析精度最左右的因素其实不是求解算法，而是磁石材料数据品质。完整的温度依赖B-H曲线（特别是knee point附近）决定了精度。通用手册只有20°C和100°C两个点，中间温度都靠线性插值，这个插值在knee point处精度不够。JMAG与信越化学等磁石厂的数据合作，提供从-10°C到200°C的详细曲线，这是它在日本电机设计中被广泛采用的原因。
永久磁石电动机减磁的故障排除
减磁判定过度
老师，分析说磁石一半以上减磁，但实机没这么严重…怎么回事？
减磁过度评估的典型原因有：
| 原因 | 症状 | 对策 |
|---|---|---|
| 网格太粗 | 磁石角部反磁场模糊 | 磁石细分到0.2mm以下 |
| B-H曲线外推 | knee point以下数据缺失，线性外推 | 获取磁石厂实测完整数据 |
| 温度过高 | 设定温度超过实际 | 热耦合确认实际温度 |
| 磁化方向不对 | 配向偏差导致局部反磁场过大 | 确认磁石厂配向数据 |
非线性收敛失败
减磁分析NR法不收敛。怎么办？
减磁分析特有的收敛问题，knee point曲率突变是主因。试试这些：
- 阻尼系数：NR更新步长乘以0.3〜0.5（JMAG的"松弛系数"，Ansys的"Damping Factor"）
- 电流分步：不要一步到最大，分10%→50%→100%
- B-H曲线光滑化：knee point附近补充数据点，用样条插值
- 初值改善：前一转子角度的解作为下一步初值
与实测数据偏离
模拟说不减磁，但实机减磁了。反过来也有吗？
有的。模拟和实机不符的主要原因：
- 磁石温度被低估：只测表面温度，内部更热。IPM电动机的磁石被铁包围散热难。
- 过渡性温度峰值：定常分析看不到短时间内的温度尖峰（短路故障的几秒）。
- 磁石个体差异：同批磁石保磁力差异±10%。设计用典型值（typ.），实际最坏要用最小值（min.）。
- 装配初期减磁：磁石插入转子时的冲击或着磁不良引起初期减磁。
模拟不只是算，还要考虑材料数据质量、制造偏差…综合工程呀。
正是。减磁分析是"跑个软件"而已，还要靠工程师的经验判断才能做好。
设计检查清单
最后，减磁分析做完提交前要检查什么？
设计检查清单。全部检查过再出结果：
- □ 磁石B-H曲线用的是温度依赖数据（不是只有20°C）
- □ knee point以下的B-H曲线有实测数据（没有的话注意外推精度）
- □ 磁石温度是最坏工况值
- □ 电流条件包括最大电流或短路故障
- □ 考虑了磁石公差（$B_r$ 下限、$H_{cj}$ 下限）
- □ 网格收敛性验证过（磁石网格加倍细化，结果变化≤10%）
- □ 减磁裕度（DM）符合公司标准
- □ 部分减磁后的扭矩降低率评估过
- □ 如用铁氧体，低温条件也检查过
具体怎样用算法求解电动机不可逆减磁？
老师解释得太清楚了！电动机减磁的困惑解除了。
价值
详细
错误