去年做某工业驱动电机方案时我遇到过一件印象很深的事电机在额定工况下所有性能仿真都合格效率和转矩波形也都好看可一旦把短路故障和持续高温叠加上去磁钢工作点直接掉到第二象限深处反电动势掉得肉眼可见。那时候我才真正意识到MAXWELL永磁同步电机退磁仿真不是交付文档里拿来凑数的分析页而是决定这个设计方案敢不敢往下走的关键环节。这篇文章就围绕这个主题展开。我会从磁钢材料模型、仿真模型搭建、退磁判据、最恶劣工况模拟、安全余量设计到实际踩坑的经验把一套完整可复现的退磁仿真思路讲清楚。项目正文里没有提到的细节我会用工程上最常用、最稳妥的做法补上。内容适合三类人看正在做永磁同步电机设计但还没系统做过退磁评估的工程师刚接触Maxwell想用瞬态场做故障工况分析的新手以及需要把退磁仿真结论落到设计改进上的产品开发人员。1. 先从一次返工说开退磁仿真到底在防什么1.1 退磁不是简单的“磁场变弱”永磁体的工作点位于B-H曲线的第二象限。这个象限里的曲线并不完全线性而是存在一个明显的“拐弯”区域。如果磁体受到的退磁磁场强度不大工作点沿着曲线还能回到原来的位置这叫可逆退磁就像橡皮筋拉扯后还能复原。但如果反向磁场强度超过某个临界点把工作点推到膝点拐点以下即使外部磁场撤掉磁体内部的磁畴排列也发生了不可逆的翻转剩磁就回不来了。我见过不少同事把退磁理解成“磁力变小了”其实这个想法很危险。退磁一旦发生是永久性失效。电机表现为反电动势下降、额定转矩下降、转速波动变大严重时还会出现失控和过流。更麻烦的是它不会像温升超标那样通过改善散热就能恢复你把温度降回去磁也不会回来只能换磁钢重新装配。真实的退磁场景主要有三种。一是持续高温尤其是内置式电机转子内部散热差磁钢温度可能比绕组温度还高二是短路冲击无论是对称短路还是单相接地定子绕组会瞬间产生远超额定值的电流形成一个很强的反向电枢磁场直接把磁钢打回膝点以下三是长期老化虽然这个过程缓慢但会逐步消耗磁体的可用工作点余量。三者叠加时退磁风险会成倍放大。1.2 为什么偏偏用Maxwell做这活儿退磁评估有很多手段。最粗糙的是解析公式法通过电枢反应磁场幅值和磁钢矫顽力做一个简单比较判断有没有越过膝点。这种方法速度极快但有个致命缺陷它只能告诉你磁钢整体可能处于什么状态完全看不到局部退磁。实际工程中电机里的退磁恰恰是先从磁钢某个角落开始的比如内置V型磁钢的肩部、表贴式磁钢的端部这些位置磁通集中最容易被反向磁场击穿。Maxwell这种有限元工具的价值就在这里。它把磁钢剖分成大量细小单元每个单元都能计算出磁通密度的大小和方向。你可以直接看到在某个故障时刻磁钢内部哪些单元的工作点落到了膝点以下哪些还在安全区退磁区域的面积占比是多少。这种空间分辨能力解析公式给不了也是设计改进方向判断的依据。再加上Maxwell的瞬态求解器可以模拟短路发生后电流的衰减过程不是简单灌一个恒定大电流进去看稳态结果而是把故障发生的瞬间、电流冲击的峰值、逐渐衰减到稳定的完整过程都算出来。退磁判定在冲击峰值出现的几个毫秒内最关键所以这种瞬态能力几乎是不可替代的。2. 动手之前先搞懂磁钢材料模型和工作点2.1 B-H曲线膝点才是退磁判定的分界线永磁体在电机里工作时承受的磁场方向与其磁化方向相反所以工作点一定落在第二象限。B-H曲线的第二象限通常呈现这样的形状起点是剩磁Br磁密随反向磁场增大而缓慢下降在某一点后突然急剧下降这个急剧下降的转折处就是膝点。退磁仿真最核心的判断逻辑就是看磁钢内部单元的工作点是否越过膝点。注意这里说的“越过”不是只看某一个瞬间。磁场是随转子位置和电流实时变化的工作点在第二象限的轨迹会来回移动。判定退磁的准则是要看最恶劣时刻的工作点位置如果某单元的最小磁通密度Bmin低于对应温度下膝点处的磁通密度Bk就判定该单元发生不可逆退磁。实际工程中我不会只看单点而是统计整个磁钢体积内有多少比例的区域低于Bk。这个比例和电机性能下降幅度直接相关。设计上通常要求正常工况下所有区域都远高于Bk故障工况下允许部分区域短暂越线但面积要有上限。2.2 温度对磁性参数的影响远比你想的大磁钢的Br和Hcj内禀矫顽力都随温度变化。对常用的烧结钕铁硼来说Br的温度系数大约在-0.12%/°C左右Hcj的温度系数更夸张通常在-0.5%/°C到-0.7%/°C之间。这意味着温度从20°C升到120°C剩磁下降约12%而内禀矫顽力可能下降一半以上。内禀矫顽力下降意味着什么它代表磁体抵抗反向磁场的能力变弱了。B-H曲线膝点的位置会随着温度升高明显向右下方移动也就是说允许承受的反向磁场强度变小了。把这一点翻译成工程语言同样一个短路故障在室温下可能只是让磁钢工作点靠近膝点但在高温下就直接越过膝点发生退磁了。这就是为什么退磁仿真不能拿室温材料曲线硬算必须把温度作为关键变量输入。我在实际项目中磁钢温度一般按最恶劣连续工作温度来取。比如某电机定下的转子磁钢最高温度为140°C那我做短路退磁仿真时磁钢材料直接按140°C或更高一档的曲线来设置。具体温度值从热仿真或者实测温度数据里拿但不做折中宁保守勿乐观。2.3 材料库数据导入的常见坑Maxwell自带材料库里有常用磁钢牌号但项目里用的磁钢多半来自指定供应商。供应商给的B-H曲线通常是20°C下的实测数据甚至只有第一象限数据这会出两个问题。第一只有第一象限数据没法做退磁判断。退磁仿真的主战场是第二象限所以一定要找供应商要完整的第二象限B-H数据并且标出各个温度点下的膝点位置。我曾遇到供应商给了一张看起来挺漂亮的曲线仔细一问是典型值不是实测值高温段和外推值差别很大最后换了实测数据后仿真结论直接翻盘。第二很多人搞混了Hcb和Hcj。Hcb是磁化强度矫顽力对应B-H曲线的第二象限Hcj是内禀矫顽力对应M-H曲线。在Maxwell里设置永磁材料时需要给出Br、Hcb更关键的是给完整的内禀退磁曲线。如果只给了Hcb而缺Hcj相关的内禀数据软件在强退磁场下会失真算出来的退磁区域是不准的。这个问题在弱磁调速和故障短路工况下特别突出。有一条很实用的经验拿到供应商的磁钢数据后先在材料设置里把多个温度点的B-H曲线都填进去。Maxwell支持定义材料属性随温度变化填好后就可以在求解参数里设置磁钢温度让软件自动插值取对应曲线。这个方法比手动分别设置不同材料要省事得多也避免在不同温度方案之间切换时改漏材料。3. 搭建退磁仿真模型从几何到求解的全链路设置3.1 模型简化2D还是3D边界条件怎么用绝大多数永磁同步电机的退磁仿真用2D模型就够了尤其是表贴式或常规内置式结构。2D模型计算快方便做多工况扫描对于一个完整的设计迭代周期来说效率优势非常明显。但用2D时要注意端部漏感没算进去短路电流的衰减时间常数会比实际偏快导致退磁结果偏乐观。处理办法是在激励回路里串一个端部漏感修正值这个值可以从3D静磁场计算或者经验公式里取我会把它作为退磁校核模型的一个固定修正项。模型范围上如果电机极对数和绕组具有周期对称性可以用最简极对模型加主/从边界条件比如一对极甚至半个极的结构。这样做的前提是故障工况下的三相电流仍然满足周期性对称。三相短路时三相绕组电流不是完全对称的直流分量会衰减但电机本体结构依然满足旋转周期性所以用周期边界没有问题只是电流激励要按故障回路来设置。如果做单相接地这种非对称故障建议至少用全模型或者半模型否则主从边界条件会把非对称性错误消除掉。3D模型什么时候必须上我通常在面对磁钢分块、斜极、端部磁密特别集中的设计时才会建3D。3D退磁仿真的核心难点在网格数量磁钢内部要做多层剖分气隙网格也要加密算一个瞬态工况往往要数小时到一天。所以我的习惯是先用2D把材料和结构方案定下来最后对最恶劣工况在3D里跑一遍确认结论。3.2 网格剖分磁钢和转角区是最不能省的地方退磁是局部现象网格精度直接影响判断。网格太粗时磁钢表面的磁场梯度被平均掉原本发生在尖角的局部退磁可能完全消失也可能被夸大。我在这上面翻过车某一次网格较粗时算出某工况退磁面积占磁钢的8%加密之后只剩3%而且退磁位置从“整体表面”缩小到“磁钢转角局部”。两种结论对应的设计行动完全不同。网格设置的要点有几条。磁钢区域按层剖分至少分三到四层转角处和倒角处手动加密气隙部分网格要细密至少有若干层单元跨越气隙否则气隙磁场计算失真磁桥、隔磁桥这些磁路关键点要做局部加密因为这些位置的磁通密度极高和退磁区域的磁场分布密切相关。我做退磁仿真前会固定做一次网格无关性验证把磁钢区域的网格尺寸从“粗”到“细”加密一倍对比最关键工况下磁钢单元的最小磁密误差稳定在2%以内才认为这个网格密度可靠。这一步花不了太多时间但能避免后续整个方案评估建立在错误网格上。3.3 激励源与运动设置的正确姿势退磁仿真不能空转得让电机工作在一个真实工况下然后把故障叠加进去。建模时绕组的激励方式有电流源和外电路两种选择。正常工况仿真用电流源比较方便给定三相对称正弦电流幅值和频率按目标工况来。比如做额定转矩下的退磁检查就把电流设定为额定有效值对应的峰值频率按电转速算。如果要做弱磁区仿真还要考虑电流角因为弱磁电流会产生直轴去磁分量对磁钢工作点的影响和短路电流方向类似。外电路一般在模拟故障时用。最典型的做法是搭建一个包含三相绕组、线路电感和短路开关的电路模型在一个稳定运行时刻合上短路开关让电流自己衰减。Maxwell里加外电路需要把绕组终结方式改成外部激励并在Simplorer或者Maxwell Circuit Editor里搭回路。我自己的经验是如果只想快速评估短路退磁直接在激励里把电压源设置成零也能近似等效但电流衰减时间常数会和实际有偏差。要求严一点还是上外电路。运动设置方面做短路仿真时我通常把转速固定不考虑转速骤降。原因是退磁最关键的时刻是短路后最初几个毫秒的电流冲击峰值这期间机械惯性导致转速几乎没变化。固定转速可以简化计算也不会对退磁结论造成实质影响。3.4 求解步长和收敛控制瞬态求解的步长直接决定能不能捕捉到电流冲击峰值。短路发生后的前几个电周期电流变化剧烈如果步长太大峰值会被平滑掉退磁判断偏乐观。我的做法是正常工作段用较粗步长比如每个电周期32到64步在短路发生时刻附近把步长压缩到电角度的0.5°到1°对应的时间这样才能捕捉到峰值的真实高度。非线性求解收敛也需要留意。磁钢材料在高退磁场下工作点接近并越过膝点时磁导率发生剧烈变化迭代容易不收敛。把Maxwell的非线性残差设置得稍低一点比如相对残差设为0.5%同时把每个求解步内的最大迭代次数放宽。如果出现收敛困难优先检查网格而不是盲目放宽残差否则得到的是一个“看起来收敛但物理上错误”的结果。4. 退磁判据与后处理怎样才算“退磁了”4.1 磁通密度云图看什么仿真跑完第一件事是看磁钢内部磁通密度云图。注意看的不是整个电机的磁密云图而是单独把每块磁钢区域提取出来观察磁通密度低点在哪里。退磁区域的位置很有规律。表贴式电机通常是磁钢外表面朝向气隙的边角因为这里的磁路短反向磁场最强内置式电机则通常是V形磁钢靠近磁桥的肩部或尖角位置这些地方磁通集中磁阻突变工作点最容易掉下来。如果云图显示某个局部区域磁密明显低于周围而且低于该温度下膝点对应的磁密这就是退磁的初始点。云图只能给出直观印象需要进一步量化。磁密并不是均匀的真正的判断要看最低值。我一般会在磁钢内部找最小磁密点观察它的位置和数值并把它随电流角、故障时刻的变化趋势提取出来。4.2 退磁比例和最小磁密的量化统计更严格的做法是在后处理里设定一个退磁阈值。定义方式基于该温度下的退磁曲线取膝点对应的磁密Bk低于Bk的单元视为退磁。Maxwell里可以用场计算器写一个表达式把所有磁钢单元里B小于Bk的单元面积求和再除以磁钢总面积得到退磁面积比例。这个比例是我评估方案安全性的核心指标。具体到工程尺度我会同时看两个量磁钢最小磁密Bmin和退磁面积比例。Bmin反映最危险点的状态退磁面积比例反映整体受损害程度。如果Bmin略低于Bk但比例不到百分之几说明只是局部轻微越线可以通过优化磁桥或加厚磁钢来解决如果比例达到10%甚至更高基本宣判这个设计方案在故障工况下不合格需要换材料或改拓扑。场计算器操作对新手有点门槛提供一个替代方案直接导出磁钢区域所有网格单元的磁密数据放到外部脚本里统计。Maxwell可以导出网格单元列表及对应的场量值我用Python脚本做过一版自动化统计输入是仿真结果文件输出每一块磁钢的最小磁密和退磁占比。一开始调试花了一点时间但后面做多方案扫描时效率极高。4.3 用反电动势和转矩做间接验证退磁仿真的结论不能只停留在“磁钢里面磁密低了”这个层面要把它转化为电机外部可测的性能变化。最常见的间接验证指标是反电动势。做法是在退磁仿真前后分别计算相同转速下的绕组反电动势。退磁导致有效磁链下降反电动势幅值按比例下降。比如某电机在短路故障后磁通损失5%反电动势大约也降5%。这个数据和台架实验里的反电动势测量非常吻合可以直接作为仿真判据的交叉验证。转矩也有同样的规律。退磁后电机的转矩常数Km下降同样电流下输出转矩变小。如果要量化评估退磁对整车性能的影响我会计算退磁前后在额定工作点的转矩差以及在高转速区的功率输出能力变化。不过要注意空缺区直轴电感变化会改变电流角分配所以转矩变化不完全等于磁通变化需以电磁仿真结果为准。5. 最恶劣工况模拟短路冲击下的退磁评估5.1 在Maxwell里搭一个能复现的三相短路场景短路工况是退磁仿真里最核心的故障场景。三相突然短路时定子电流会在几个毫秒内冲到额定电流的七八倍甚至十倍以上这个冲击电流产生的电枢磁场直接作用在磁钢上加上温度已经处于最恶劣状态退磁概率最高。搭建步骤先在稳定工况下把电机跑起来稳定几个电周期后合上短路开关。这里有几个细节要处理。短路点要放在绕组出线端短路阻抗要包含实际线路阻抗和接触电阻。很多仿真把短路阻抗设为零认为这是最恶劣情况。实际上零阻抗会让电流衰减时间常数更长冲击峰值确实偏大但真实系统中接触电阻和导线电阻会消耗一部分能量。我一般用相对保守的0到极小阻抗来做确保仿真不会低估退磁。绕组连接方式也有讲究。三相短路分为三相对称短路和单相短路。三相对称短路时三相电流都很大但各相去磁作用和转子位置相关单相短路时非对称电流会产生更大的负序分量对某些位置磁钢的退磁作用反而更严重。所以不能只看三相短路具体工程中要依据电网和驱动器的保护逻辑把可能发生的短路类型都扫一遍。5.2 电压相位角扫描退磁的最恶劣点藏在角度里很多人做短路仿真时随便选了一个转子位置就开始算结果得出退磁风险很低的结论换一个初始角度再算结论完全不同。原因是短路电流的瞬态冲击和初始电压相位角、转子位置角之间存在强耦合。短路瞬间的直流分量大小取决于该时刻三相电流的瞬时值而这个瞬时值和相位角直接相关。我的做法是分两轮扫描。第一轮以30°为步长在0°到360°之间扫描电压相位角跑完所有工况后记录退磁面积比例和磁钢最小磁密。第二轮在最恶劣角度附近用5°甚至2°的细步长加密扫描确定真实的峰值退磁点。某次项目中同样的电机在30°扫描时退磁比例只有2%看着很安全但加密扫描后发现真正的峰点在特定角度下退磁比例达到了6%以上差了至少一倍。这就是角度扫描的价值也是为什么我一定会在故障仿真里预留足够多的工况计算时间。5.3 热-磁耦合为什么这一步不能省退磁仿真的温度设定很关键。磁钢最高工作温度直接决定B-H曲线膝点的位置也就直接决定了短路时会不会越过膝点。很多人拿室温曲线算短路退磁结果算出来一切安全样机实测却退磁了原因就在温度上。正确做法是先通过热仿真或实测拿到磁钢在极限工况下的稳态温度把这个温度作为退磁仿真里磁钢的预设温度。短路过程只有几十到几百毫秒温度在这期间来不及明显变化所以不需要做真正的瞬态热-磁耦合用固定温度即可。但要注意的是这个固定温度必须取极限值而不是平均值。我见过一个项目案例电机额定工况下磁钢温度约90°C但堵转保护动作前有约5秒的过载窗口这个窗口内磁钢表面温度可以冲到130°C以上。如果只按90°C评估膝点位置偏高短路工况看着安全改成130°C后同样工况直接判定退磁。所以温度取值要覆盖控制器最恶劣的保护响应时间不能只看稳态温度。6. 实际项目中的退磁安全余量与设计改进6.1 退磁安全系数怎么取才不会走极端退磁评估最终要落到设计是否安全。我常用的判断逻辑是两套标准正常工作工况下磁钢最小工作点磁密必须高于膝点磁密并留有余量故障工况下允许短时越线但越线面积要控制在小范围内且瞬态过后磁钢不会发生持续性的性能衰减。工程上常把“退磁余量”定义在最恶劣工况下磁钢工作点磁密与膝点磁密的差值上。正常工况建议至少有0.1T到0.2T的余量。这个数来自经验很多磁钢材料在接近膝点时已经出现部分不可逆损伤虽然曲线表象上还在膝点以上但微观上已经有局部磁畴翻转。所以你不能卡着膝点做设计必须留出足够的安全带宽。故障工况的经验法则是退磁面积比例不超过5%且退磁后的电机性能下降不影响安全功能。这个比例怎么确定把退磁后的反电动势降幅和系统堵转保护阈值关联起来。如果反电动势下降超过8%驱动器的某些保护判据可能失效那就必须降低允许的退磁比例。6.2 材料升级和结构优化哪个更划算当退磁仿真判定某个设计方案不安全改进方向有两个换材料或者改结构。换材料是最直接的。NdFeB磁体按耐温等级从低到高有N、M、H、SH、UH、EH等级别耐温越高内禀矫顽力越大膝点在高温下越靠左。将磁钢从N系列换成SH甚至UH系列短路退磁风险能大幅下降。但代价是磁钢成本上升剩磁可能略降需要重新评估电机性能。改结构的方向更多。增加磁钢厚度这是简单粗暴但有效的方法磁钢厚了磁阻增大同样的外部反向磁场在工作点上产生的磁密偏移变小抗退磁能力增强。对内置式结构调整磁桥宽度和形状加强磁钢尖角部位的磁路疏导也能让局部工作点回升。表贴式结构可以加护套或包覆结构但这主要解决机械问题对退磁帮助有限。我个人的经验是优先优化磁桥和磁钢厚度成本增加不大性能影响小。只有当结构优化空间受限时才考虑升牌号。一个项目里通过把磁钢厚度增加0.8mm加上调整磁桥宽度退磁面积比例从9%压到2%以内而换UH牌号只能压到5%还不一定够两相对比结构优化的性价比高很多。6.3 仿真结论要和台架验证互相印证再好的仿真也是模型最终要对上实测才作数。退磁的台架验证主要通过对比退磁前后的反电动势和转矩数据。具体做法是样机在常温下测一组反电动势和转矩基线然后做短路或堵转试验待电机冷却后再测一组两次数据的差异就是不可逆退磁的量化证据。还有一个容易忽略的点每一批磁钢的材料一致性。同一牌号的磁钢不同批次之间的Br、Hcj离散度可能达到3%到5%。仿真用的材料曲线是典型值如果按典型值刚好安全实际批次产出的磁钢可能有一批不合格。量产评估时我会要求供应商提供批次数据的Ppk统计并把最小Hcj作为仿真输入而不是用平均值。设计的最后阶段还建议加一道磁通检验工序逐台检测转子磁钢的磁通量把低于阈值的转子直接拦截在产线外。7. 这类仿真里反复踩过的坑都在这里了7.1 网格越细退磁区域越小这是真的坑有次仿真算出来磁钢退磁面积高达15%按照这个结果整个设计要推翻重来。我当时没有立刻信这个数字做了一次网格加密验证结果面积降到4%。初看很反直觉网格细了退磁区域反而小了。原因是粗网格把局部磁场梯度抹平了平面云图看着整个磁钢表面都低于膝点加密后磁场分布更精确真正的退磁点集中到了磁钢尖角极小的区域里。所以退磁仿真里看到大面积退磁时第一反应应该是检查网格密度尤其是磁钢转角处而不是直接改设计。这个规律反过来也提醒我凡是看到退磁面积在5%上下的结果都值得先做网格敏感性分析再下结论。7.2 不要用静磁场代替瞬态场做退磁判断这个坑我劝过不少人。静磁场求解速度快而且能直接查看额定电流下的磁密分布所以有人拿静磁场灌入大电流看一下磁钢最低磁密就宣布短路退磁安全。问题是短路退磁的关键特征是电流冲击和衰减的瞬态过程。静磁场既无法反映转子运动导致的气隙磁场变化也无法模拟短路发生后电流的直流分量衰减更看不到故障瞬间工作点的完整轨迹。我见过用静磁场算出来“安全”的电机实测试验后反电动势下降明显照样退磁。退磁仿真的底线是用瞬态场尤其是涉及故障时。静磁场最多只能用来快速做方案初筛不能作为交付依据。7.3 磁钢材料温度在软件里很容易设错Maxwell里的材料温度设置有一个很隐蔽的坑如果磁钢材料定义为带温度特性的B-H曲线你需要指定工作温度软件才能插值取到正确的曲线。但如果你在模型里只改了环境温度而没改材料属性里的工作温度软件仍然会使用默认温度下的曲线。我犯过一次这样的错环境温度设置了120°C材料温度忘了改默认25°C的B-H曲线让仿真结果完全失真。这个错误从云图上看不出来因为磁密的分布形态相同只是数值偏乐观。直到和台架数据对比才发现问题。现在我的检查习惯是在求解前逐项核对材料属性里设置的磁钢温度是否和热分析给出的最高工作温度一致。7.4 批量扫描角度脚本自动化值得做退磁仿真涉及大量角度扫描和工况组合手动逐个设置不仅慢还容易漏掉关键工况。我用Maxwell脚本接口写了一个批处理工具能够自动修改激励源的初始相位角、磁钢温度和故障开关时序批量提交求解最后汇总每块磁钢的最小磁密和退磁面积比例。这个工具让我能把一组完整的短路退磁扫描从两三天压缩到半天。脚本思路不复杂把仿真工程文件里的关键参数抽出来循环赋值每个角度提交一次仿真把结果导出到统一表格。如果你手头有多个退磁工况要批量评估这个方法值得花时间搭建起来收益非常可观。另外提醒一句仿真里的磁钢材料数据尤其是高温段的内禀退磁曲线一定以供应商实测数据为准。不同供应商给的同牌号材料高温矫顽力可能差异很大你用软件默认库里的近似曲线做出来的结论可能和量产件完全不同。做这类安全评估材料的准确度比求解精度更容易影响结果别在这上面图省事。
