搞电机设计的人大概都绕不开内置式永磁同步电机IPMSM这个名字。这几年新能源驱动电机、压缩机电机、伺服系统大量采用这个拓扑原因说白了就两条一是永磁体埋在转子铁芯里面结构上比表贴式结实高速运行时不怕离心力把磁钢甩出来二是转子磁路不对称带来的磁阻转矩可以拿来用整个电机的过载能力和弱磁扩速能力都上了一个台阶。但仿真这件事很多刚入门的朋友跟我吐槽过模型倒是不难画但要么算出来的反电动势波形不对要么转矩怎么也调不到目标值更头疼的是电感参数和别人对不上。我之前在好几个项目里反复调试IPMSM模型从空载反电动势、齿槽转矩到额定负载工况把常见的坑基本都踩了一遍。这篇就把我实际用下来的建模流程、参数设置、网格处理和结果判读方法整理出来希望能让正在做IPMSM仿真的人少走点弯路。1. 内容整体设计与思路拆解1.1 IPMSM仿真的核心难点在哪里先聊一个基本问题为什么专门写内置式而不是表贴式表面上看IPMSM只是把磁钢从转子表面挪到内部但这一挪整个电磁场的分析难度就变了。内置式转子结构里磁钢被硅钢片包裹磁通路径变得复杂得多——磁钢产生的磁力线一部分经过气隙进入定子参与机电能量转换另一部分直接在转子内部形成漏磁回路。这块漏磁不仅被磁钢的排列方式影响还和转子隔磁桥的宽度直接挂钩隔磁桥越窄漏磁越小但机械强度越差。仿真里如果不把这个矛盾处理好算出来的空载反电动势和实际测试数据能差到百分之十以上。另一个难点是交直轴电感。IPMSM天然存在凸极效应Ld和Lq不一样这是磁阻转矩的根源。但内部的磁钢、隔磁桥、转子槽形对电感的影响是非线性的随着电流增大磁路饱和程度变化Ld和Lq都在动态变化。这意味着仿真里不能只用一套固定电感参数估算性能得用有限元把不同电流角下的转矩特性扫出来。做联合仿真或者标定的时候这些数据直接决定控制策略里的最大转矩电流比MTPA查表能不能做准。1.2 样机参数与仿真目标设定仿真的第一步不是打开软件而是先把目标定清楚。我这里用一台常见的新能源辅助驱动油泵电机做例子具体参数如下表额定电压320V额定转速4000rpm峰值转速9000rpm额定功率5.5kW。这是个典型的IPMSM设计主要看重的是恒功率弱磁范围够不够宽、高效区覆盖面积够不够大。参数数值备注额定功率5.5 kW连续运行点峰值功率11 kW短时过载工况额定转速4000 rpm额定转矩点峰值转速9000 rpm弱磁运行上限额定电压320 VDC母线电压极数8极转子4对极槽数48槽整数槽绕组磁钢材料N35UH钕铁硼耐温180℃转子结构V型双层磁钢常见内置式布局选V型双层磁钢的布局主要是为了在有限的转子空间里塞进更多磁钢面积同时通过双层结构优化气隙磁场的正弦度。这个结构在量产驱动电机里很常见它的仿真结果对绝大多数IPMSM项目有参考意义。我给自己定的仿真目标是额定点效率不低于94%峰值转矩下磁钢不退磁空载反电动势谐波畸变率控制在10%以内。目标定好之后后面每一个环节的取舍都有依据了。比如网格加密该放在哪、求解精度要到多少、后处理要看哪几条曲线全都可以对着目标来。2. 几何建模与材料参数设置2.1 从图纸到有限元模型的转化拿到电机的二维图纸后第一件事不是急着在有限元软件里画图而是先把图纸上的关键尺寸整理一遍。对于内置式电机最重要的是磁钢槽的精确尺寸和位置因为磁钢槽的形状直接决定隔磁桥宽度而隔磁桥宽度又决定了漏磁系数和转子强度。我在处理这类模型时一般用CAD软件把转子冲片单独导出为DXF格式再导入仿真软件。直接用CAD软件建模的好处是尺寸不容易走样尤其对V型磁钢槽这种带角度的结构手工在仿真软件里画极易出错而且后面调整方案时想参数化也很难做。需要注意的一点是CAD文件导入仿真软件后一定要检查导入的曲面有没有自相交或者重复边的问题。有一次我导入DXF之后没检查直接进入求解结果网格剖分一直报错检查半天才发现转子内圆和转轴外圆有两段重复边删掉之后剖分立刻就正常了。如果只是做电磁仿真三维模型可以简化成二维模型来做。内置式永磁电机沿着轴向的磁场变化相对很小轴向长度远大于气隙长度端部效应对主磁场的影响有限所以二维截面模型既有足够的精度计算速度又快了不止一个数量级。这个简化对前期方案迭代非常友好我个人的习惯是先把二维模型跑通确认真没问题后再根据需要建立三维模型做端部漏感或者斜极分析。2.2 材料属性设置的几个关键细节仿真模型里需要用到的材料主要有三类定转子硅钢片、磁钢和绕组铜线。硅钢片材料的B-H曲线必须设置完整甚至要包含不同频率下的铁损曲线。有人图省事直接选材料库自带的默认硅钢片但那些默认值往往和实际用的牌号差很远。我做这个项目时用的是某国产35JW250牌号硅钢片它的饱和磁密在1.7T左右比材料库默认值略低如果不替换算出来的铁损和电感参数都会偏乐观。磁钢材料方面除了剩磁Br和矫顽力Hc之外最容易被忽略的是磁钢的导磁特性。钕铁硼的回复磁导率通常在1.05左右很多人嫌麻烦直接设为1也就是理想化的不导磁材料。这看似影响不大但在内置式电机里磁钢块之间、磁钢与硅钢片之间的磁场分布都会发生细微变化导致计算出来的磁阻转矩偏大。我对比过同一模型两种设置的差异转矩结果相差约3%对于要卡效率指标的仿真来说这个误差不能忽视。绕组方面定子铜线按75℃工作温度下的电导率来设置比较贴近实际工况。冷态铜的电导率比热态高不少如果按20℃设置铜损会被低估约10%。特别是做效率MAP仿真时温度选取会影响整个高效区的分布走向建议统一按工作温度来。2.3 磁钢充磁方向与转子区域设置这是内置式电机仿真里最容易翻车的地方。表贴式电机的磁钢充磁方向很简单径向或者平行充磁直接选一个就行。但内置式V型布局的磁钢是斜着放的每块磁钢的充磁方向都得单独定义而且必须沿着磁钢厚度的法线方向。我在实际建模中习惯先建立一个局部坐标系让它的一个坐标轴严格垂直于磁钢表面然后给磁钢材料指定这个方向作为充磁方向。V型双层结构里上下两层磁钢的夹角不同需要分别建坐标系。用全局坐标系直接定义充磁方向的后果是磁钢的工作点算不准空载反电动势波形和实测对不上甚至转矩波动都会异常偏大。转子区域的设置为的是让有限元软件能正确识别哪些区域是导体、哪些是铁磁材料。需要在软件里额外指定转子区域的运动属性包括旋转速度、初始位置角。初始位置角的设置直接关系到三相感应电压的相位做空载反电动势分析时角度设错了得到的波形相序就是错的。通常的做法是先跑一个静态场查看A相绕组的磁链过零位置把转子位置对齐到那个角度上再开始瞬态仿真。3. 绕组设置与激励源配置3.1 绕组相序和匝数怎么填一个48槽8极电机是典型的整数槽绕组每极每相槽数为2。绕组的设置要做两件事把不同的槽分给A/B/C三相然后指定每槽导体的匝数和连接方向。在软件里这个工作通常通过绕组编辑器或者槽填充定义来完成。槽号分配必须和实际下线表完全一致这里最容易犯的错误是只看了一个极距范围内的槽号分配剩下的区域直接按规律复制。要注意整数槽绕组虽然具有周期性但并联支路数不同槽分配的结果也不同。我的习惯是先画出槽电势星形图把每相所属槽号验证一遍再填入软件。关于匝数设计上每槽导体数为44采用双层绕组上下层各22根。软件里填写绕组的时候要区分导体数conductors和匝数turns两个概念对单匝线圈来说两者相等但如果是多匝绕制导体数等于匝数乘以并绕根数。我之前接手过别人一个项目文件对方在软件里把匝数填成44实际上设计意图是22匝、每根匝间两根并绕结果转矩直接算大了一倍查了很久才找到问题。这类低级错误在仿真里造成的后果比技术问题更难发现大家填完参数之后一定要先做个小手算验证。3.2 激励源类型选择的学问激励源的设置取决于你关心什么性能。做额定点效率计算时用电流源激励配合转速得到的是该工况下的电磁性能做空载反电动势分析时直接把激励设为零电流或者把绕组Open circuit处理测量的是感应电压。电流源激励下需要根据目标转矩反推电流幅值和电流角。最简单的方法是先用一个初始电流角按照MTPA方向估算电流幅值跑一遍后查看输出转矩再根据偏差调整电流幅值。比如我这台电机目标额定转矩13Nm我先按经典MTPA公式算出d轴电流约-35Aq轴电流约120A。把这些值设置到三相正弦电流源里频率按4000rpm对应的8极频率来设置也就是f n×p/60 4000×4/60 ≈ 267Hz。加上电流源之后要注意相位配合。三相电流Ia、Ib、Ic各自相差120度电角度软件内部会自动配置但电流源的初始相位要和转子初始位置对应起来。这之间如果没有对齐起步瞬间会出现转矩突变虽然稳态结果不受影响但效率计算的积分过程会引入误差。我一般先把转子位置角设置成使A相磁链为零的位置再设置A相电流初始相角为0这样整个加载过程就是平顺的。3.3 端部漏感的等效处理二维模型无法直接模拟绕组端部伸出铁芯的部分但端部漏感对电机性能有明显影响尤其是弱磁工况下的电感和电压约束计算。处理办法有两种一种是在电路方程里额外添加端部漏感参数让外部电路计算时把这一部分补进去另一种更省事的方法是根据经验公式估算端部漏感值在模型参数中通过等效手段加入。对于额定工况仿真端部漏感占比约3~5%影响尚可接受。但做短时过载或者峰值转速以上的弱磁工况分析时端部漏感对电压余量影响就变大了建议务必把它加进去。我通常的做法是按绕组跨距和铁芯长度估算端部电感大约每相0.06mH左右然后写进外电路。这个数值不用特别精确因为它在整个漏感中的占比有限留个大概值对趋势判断足够用。4. 网格剖分与求解设置4.1 气隙和磁钢是剖分的重点网格剖分的质量直接决定有限元解的准确性。对内置式电机来说最关键的剖分区域是气隙、磁钢槽附近的隔磁桥和气隙附近的齿尖。气隙区域我通常单独设置为一个带状区域沿径向剖分至少分成三层网格让磁力线跨越气隙时有足够的离散度。两层和三层的气隙网格计算出的气隙磁密波形可能差5%以上转矩波动差距更大。手动调整气隙部分的网格尺寸目标尺寸设为0.15mm比定子齿部默认网格细了一倍多。磁钢槽的隔磁桥位置因为磁密极高很容易进入深度饱和网格也必须加密。这个部位如果剖分太粗饱和效应算不准漏磁系数就会偏小铁损结果也会失真。我在隔磁桥宽度为0.8mm的区域把网格目标尺寸压到0.1mm级别同时开启自适应剖分让它迭代几次后再看磁密云图是否平滑。4.2 自适应剖分与带宽设置技巧Ansys Maxwell里做瞬态仿真网格剖分如果不加特殊干预算到后期可能出现转矩曲线振荡甚至发散。我的做法是先跑一个静态场开启自适应剖分以储存的磁能为目标量让软件自动迭代加密直到磁能变化小于0.5%为止。用这个剖分结果作为瞬态场的初始网格收敛性好得多。这里还有一个优化技巧把求解带宽band设置成电机旋转区域内部的所有物体随基准坐标系一起运动。band区域必须是凸多边形不能包含倒角或凹槽边界。有些设计师把转子外形做了各种圆弧过渡如果直接复制外轮廓作为band边界很容易出现边界退化导致网格重画失败。我的习惯是做一个小正方形或正多边形包围转子让里面的转子结构全部落入band内部这样网格重画稳定很多。4.3 求解步长与收敛判定瞬态求解的时间步长要匹配基波周期。我的经验是最少一个电周期内采100个点也就是每个电角度3.6度取一个解。对于8极267Hz的工况一个机械周期是0.25ms一个电周期的时间约3.75ms步长取37.5us就够了。如果要做转矩波动分析建议每步对应电角度1度左右再多加密对改善精度作用很有限徒增计算时间。收敛判定不是只看软件报的残差更要看关键物理量是否达到稳态。我一般监视转矩曲线和反电动势波形等到这两个量在一个电周期内的波动模式稳定下来才认定收敛结束。一般需要跑3到5个电周期才能获得稳定的周期解起始的那一个周期数据通常丢弃不用。转速不变的空载工况相对好收敛负载工况因为磁路饱和严重转矩起步过程伴随明显的直流分量需要更多周期才能稳定。遇到这种情形我可以直接把上一周期末的场量作为本周期初值强制软件加载到稳态节省一到两个周期的计算量。5. 典型工况仿真结果与后处理5.1 空载反电动势的波形分析和谐波评价空载反电动势是检验电机设计是否正确的最直观指标。把转速设为4000rpm绕组激励设为零电流待转速稳定后提取A相的反电动势波形。需要注意的细节是软件里感应电压的提取方式是在绕组外部接一个无穷大阻抗表示的电压表或者直接在后处理中通过磁链对时间求导获得。实测下来的波形应该是接近正弦的梯形波反电动势峰值大约在150V左右对应母线电压320V的弱磁能力初步判断是够用的。对波形做FFT分析基波幅值约148V三次谐波占比在4%左右五次谐波在2.5%左右总谐波畸变率约8.5%符合我前面设定的10%以内的指标。反电动势谐波偏大通常意味着转子磁钢布置或者定子绕组分布还有优化空间。比如三次谐波大要考虑是不是V型磁钢的夹角太小导致气隙磁密平顶波形变宽五次谐波大则多和定子齿槽效应有关。这个阶段不要急着改模型先把FFT频谱打出来看看各次谐波的相对大小再对照磁钢夹角和齿槽尺寸关系去调整。5.2 齿槽转矩的提取方法与判断标准齿槽转矩是IPMSM设计里一个很难绕开的指标。它在零电流、低速旋转时表现得最明显表现为转子上一种制动转矩波动。齿槽转矩仿真转速设为1rpm即可这样每个电周期内的点数足够密可以捕捉到完整的波动曲线又不至于因为高速导致数值振荡。齿槽转矩峰峰值做出来大约0.6Nm相对于额定转矩13Nm约占4.6%。这个水平对油泵类应用勉强可接受但对于要求低噪声、低振动的乘用车驱动电机来说偏大一般建议控制在额定转矩的2%以内。降低齿槽转矩的常用手段包括定子斜槽、转子斜极、优化磁钢槽形状和辅助槽位置。在二维模型里斜槽不好直接模拟通常用一个等效的轴向分段偏移来近似或者直接通过调整极弧系数来压低齿槽转矩。还遇到过一个特殊情况同一款转子冲片一套模具压出来的实际齿槽转矩比仿真小不少原因是仿真里定转子为理想同心实际装配存在一定偏心偏心本身会对齿槽转矩起到调制作用使峰峰值略有下降。这种细枝末节不影响到设计方案评估但在对标测试数据时不要强行要求仿真和实测完全一致。5.3 额定负载下的转矩与效率分析额定负载仿真是把电流源加上去驱动转子以4000rpm旋转。这里我用的是电流角约72度电角度的工况外特性上对应最大转矩电流比附近的电流角A相电流峰值约155A。设置完成后查看转矩波形稳定后的平均转矩在13.1Nm左右和设计目标13Nm基本吻合转矩波动峰峰值约1.2Nm占比约9.2%。对于IPMSM9%的转矩波动不算特别突出但还不是最优。我进一步检查了磁场分布发现齿顶局部磁密接近2T出现明显饱和。这部分饱和会导致转矩波动加大我尝试把齿顶形状倒角加大一点让饱和区磁路更平滑结果转矩波动降到7%左右。这说明在IPMSM里转矩波动并不仅仅由齿槽效应决定负载工况下的磁路饱和分布反而更关键。效率计算需要分别提取铜损、铁损和机械损耗。铜损根据电流有效值和相电阻计算相电阻0.04Ω电流有效值110A三相铜损Pcu 3×110²×0.04 1452W。铁损需要由软件根据磁密波形和材料损耗数据直接算包括磁滞损耗和涡流损耗两部分额定点做出来约210W。机械损耗因轴承、风摩等因素估算约100W。那么输出功率为5.5kW总损耗约1762W效率算下来94.8%略高于目标值。这里铁损实测还包含了磁钢中的涡流损耗这部分在IPMSM里不容忽视尤其是高速运行段磁钢涡流损耗占比能超过总铁损的20%以上。5.4 弱磁扩速能力校核IPMSM最大的优势就是弱磁范围宽。这个工况的仿真要考察的是在母线电压320V的限制下电机能否在9000rpm的峰值转速下仍然输出可用的转矩。弱磁分析要采用电压源-电流源混合的求解方式或者直接在特定转速下调整电流的d轴分量把电压幅值压在逆变器能输出的极限之内。额定电压320V对应逆变器最大输出相电压幅值约185V留有10%的裕量后可用电压约166V。我通过扫描不同d轴电流值找到9krpm下满足电压约束的电流组合d轴电流约-186Aq轴电流约55A。这个工作点下输出转矩约5.2Nm。从4000rpm的额定13Nm过渡到9000rpm的5.2Nm功率已经从原来的5.5kW微降到约4.9kW算是预料之中的结果对油泵负载来说功率需求本身也是随转速上升而增加的并不影响使用。弱磁工况仿真要格外注意磁钢的退磁校核。9krpm时d轴电流达到-186A这个去磁磁动势相当大如果磁钢厚度不够或者矫顽力不足工作点容易跌到退磁拐点以下。仿真里要提取磁钢单元的最小磁密确保它保持在磁钢退磁曲线拐点之上。我这台电机用的N35UH在180℃时拐点大约0.1T仿真结果最小磁密为0.32T足有裕量。6. 常见问题排查与仿真提速技巧6.1 转矩曲线发散处理转矩曲线在瞬态仿真里出现锯齿状发散多半是网格退化导致的。具体表现是转矩平均值持续漂移或者上下剧烈跳动跟正常的周期性波动完全两回事。这种时候优先检查band区域是否包住了转子全部运动部件band边缘是否规整其次把气隙网格进一步加密再有就是检查初始位置角和电流相位是否对齐相位不对会在启动瞬间输出异常冲击转矩间接导致网格重画失败。我遇到最多的还是band边界不规整引发的剖分退化。解决方法很简单在几何模型里给旋转区域单独建立一个规则多边形圆形或正方形的band保证转子内孔和转轴等结构都在内部。band不参与电磁计算它只是一个运动边界标记所以形状规整性远比贴合转子外轮廓重要。6.2 电感计算结果对不上的原因电感参数是IPMSM控制策略的基础数据。仿真提取电感有两种方法一种是静态场下施加特定d轴q轴电流依据磁链计算Ld和Lq另一种是瞬态场下的增量电感提取通常用于高频注入法控制算法验证。两种方法做出来的数值不同是正常的。静态法得到的是视在电感反映的是磁链和电流的总比值高频注入法得到的是增量电感反映的是工作点附近的微分斜率。在饱和区域这两个值能差到20%以上。有人拿视在电感去标定控制算法调试时发现转矩响应总是过冲或不足就是因为两条定义下电感含义不一样。另外一个常见的坑是电感矩阵变换时角度搞错。IPMSM的dq轴电感需要通过Park变换把三相磁链变换到旋转坐标系下如果转子位置角没用对变换出来的Ld和Lq会交叉污染Ld里混进q轴分量。建议用软件自带的Clark/Park变换模块去处理磁链不要手动做矩阵运算省得在某个正负号上栽跟头。6.3 计算时间太长的优化经验仿真精度和时间的平衡是门手艺活。我总结了几条经验能让仿真时间大幅缩短而不损失太多精度自适应剖分从静态场开始先把剖分密度固化下来再进瞬态场避免瞬态计算过程中反复重剖分。瞬态场的网格最大尺寸设置合理值不要在齿尖和槽口过度加密加密只集中在气隙、磁钢和隔磁桥三个区域。时间步长不必全程固定不同工况可差异化设置。空载、负载和弱磁工况的步长根据电周期调整不要用一个步长跑到底。开启软件的多核并行计算模式但线程数不宜超过物理核心数超了反而因线程切换变慢。一台8极48槽的IPMSM二维模型在32G内存、8核的普通工作站上做完空载负载弱磁三个工况的瞬态分析大约需要2到3小时。如果还想进一步提速可以尝试在静态场验证方案可行性后再开瞬态细算大部分初步筛选工作用静态场就能完成。6.4 转速范围较大时的重合度处理当转速从4000rpm到9000rpm做全转速扫描时如果每个点都跑一个完整的瞬态过程计算量会成倍上升。我一般先用静态场在多个转速点下计算转矩、磁链和电感对性能趋势做一个初步判断然后用磁链和转矩公式估算各转速下的效率分布。只有表现比较好的两个候选方案才投入瞬态场验证。这个流程把全转速效率MAP的计算时间从几天压到半天以内。某些软件支持参数化扫描把转速、电流幅值、电流角设成变量一次批量跑完几十个工况。做这类批量计算时注意提前检查每个工况的初始位置角是否都自动更新了。位置角不变会导致不同电流角下的仿真转矩都出现系统偏差批量跑完之后还要返工那才是最让人崩溃的。我个人的做法是批量扫描之前先用三个典型工况额定点、峰值转矩点、弱磁极限点手动试算确认设置全部正确后再开放批量任务。磨刀不误砍柴工试算的时间永远比返工的时间短。7. 仿真结果的工程对标与修正7.1 仿真与实测差异的常见来源仿真做完了最后一定要回到实测上做对标。如果仿真结果和实测差别在3%以内说明模型可信5%左右需要排查超过10%基本可以认定某个关键设置出了问题。建模环节里最容易被忽视的误差来源是硅钢片的实际B-H曲线。材料数据手册里给的是标准试样测试值实际冲压、叠压之后由于残余应力和毛刺影响铁芯损耗会增加10%到30%。仿真里用的B-H曲线是从材料库直接调的那这个差异就解释了一部分效率偏差。另一个主要来源是磁钢性能的离散性。仿真里默认所有磁钢材料属性完全一致实际一批磁钢的Br公差通常在±3%左右V型结构里四块磁钢性能不一致会直接影响空载反电动势幅值和波形。装配气隙也是重要的误差源。定转子不同心、轴承游隙造成的实际气隙不均匀会让齿槽转矩的实测值比仿真值小但会引入额外的单边磁拉力导致的振动噪声。测齿槽转矩时最好用专用的低速测试台配合高精度转矩传感器手动盘车的方式测出来的转矩曲线毛刺太多很难和仿真平滑波形做直接对比。7.2 用标定数据反向修订模型模型对标不是一次性工作而是一个迭代过程。我第一次做完对标发现仿真空载反电动势比实测高5%。检查后发现我用的磁钢Br取的是材料手册上限值1.23T而实际到货磁钢的平均Br只有1.18T。把Br改成实测均值后反电动势误差缩到1%以内。这类修正看起来很笨做起来也很耗时但每修正一次后续仿真的可信度就上一个台阶。做过几个项目之后我的经验是要建立一套模型校准记录把每台电机实测得到的磁钢参数、铁芯损耗修正系数、气隙长度偏差做成表格后续类似设计可以直接调用过去的修正系数避免每次都用理论值去碰运气。个人实操体会IPMSM仿真这件事入门门槛不高做深并不容易。很多人一上来就追求复杂的模型、精细的网格却忽略了最基础的材料参数是否准确、绕组设置是否正确。我实际做完这么多项目最大的体会是仿真的价值不在于把结果算得多么漂亮而在于能提前暴露设计里的问题。一套可靠的仿真流程能把样机试制的返工次数降到最低。如果再让我分享一条最想说的实战经验那就是做IPMSM仿真前务必要花时间把磁钢充磁方向和绕组连接关系验证清楚把设计图纸上的尺寸一个一个对应到模型里去。基础的东西做扎实后面所有的高阶分析才有意义。希望这篇分享对你正在做的IPMSM项目有所帮助也欢迎在实际调试中有新发现的朋友多做交流。