1. 从一台烧掉的样机说起:为什么电机设计绕不开电磁仿真
几年前我参与过一款无刷直流电机的开发,样机在台架上跑到额定转速的八成左右时,绕组温度在十几秒内从七十多度飙到一百三十度,紧接着就是绝缘冒烟。拆开之后发现,问题出在定子齿部的磁密已经接近饱和,铁耗远超设计预期,而当时我们的设计流程是"经验公式算一版、打样、测、改",一轮迭代就是三周。那次之后我才真正把 Ansys Maxwell 这类电磁场仿真工具纳入到日常设计流程里,因为靠手算和二维经验公式,根本抓不住齿槽转矩、局部饱和、涡流分布这些三维空间里的细节。
这篇内容想聊的就是Ansys Maxwell 在电机电磁性能仿真这条完整链路。它解决的核心问题是:在还没有实物样机的时候,把电机的磁场分布、反电动势、转矩特性、铁耗铜耗、电感参数算出来,让设计迭代从"打样-测试"变成"建模-仿真-优化-打样"。适合谁看?如果你是从机械或控制转过来做电机设计的工程师,或者已经会画电机结构但没系统跑过电磁仿真,又或者你正在被"仿真跑得慢、结果对不上实测"折磨,那这篇会比较对路。我会按建模、材料、激励、网格、求解、后处理、参数化这条主线,把每一步的取舍逻辑和踩过的坑讲清楚,而不是只丢一堆菜单路径。
需要先说明一点:Maxwell 不是万能的,它算的是电磁场,热和结构要靠 Workbench 里的其他模块耦合。所以整篇文章的定位是"电磁性能仿真",热-固耦合只在必要的地方点到为止。
2. 建模之前先想清楚:2D 还是 3D,这个决定影响后面所有工作量
2.1 二维仿真的适用边界与代价
Maxwell 2D 是绝大多数电机电磁仿真的起点,原因很直接:计算量小、迭代快、参数化扫描便宜。对于径向磁通的旋转电机,如果轴向结构均匀、端部效应不占主导,2D 的结果和实测的偏差通常在可接受范围内。我做过一台 48 槽 8 极的内转子永磁同步电机,2D 算出来的空载反电动势和实测差了不到 4%,转矩脉动趋势也基本一致。
但 2D 有个绕不开的代价:它默认轴向无限长,端部绕组、端部漏磁、斜槽的轴向效果全部丢失。所以 2D 算出来的铜耗会偏小,因为端部绕组电阻没算进去;漏电感也会偏小。我的做法是,2D 用来做拓扑筛选和参数扫描,把转子结构、极槽配合、绕组方案定下来,然后再上 3D 或者用 2D 加端部修正系数去校核。
提示:2D 里如果电机有斜槽,不要直接忽略。可以用多切片法,把轴向切成若干段,每段给一个转子位置偏移,最后把结果叠加,这是 2D 逼近 3D 斜槽效果的常用手段。
2.2 三维仿真什么时候真的必要
3D 不是"更准"的代名词,它是"更贵"的代名词。一次 3D 瞬态求解动辄几小时到几十小时,网格量是 2D 的几十倍。真正需要上 3D 的场景我总结下来有这么几类:轴向磁通电机、横向磁通电机、端部效应显著的短初级直线电机、需要精确算端部漏抗和端部损耗的场合,以及做电磁-热耦合时要把端部绕组的热源分布算准。
直线电机是个典型例子。定位精度和重复精度这类指标,本质上和法向力波动、推力波动强相关,而这些波动很大一部分来自端部效应和边端效应,2D 根本抓不住。所以做直线电机电磁仿真,3D 几乎是必选项。
2.3 对称性利用:把计算量砍掉一半甚至四分之三
Maxwell 里有个特别容易被新手忽略的功能:对称边界。如果电机是整周期对称的,比如 48 槽 8 极,极槽的最大公约数是 8,那你可以只建 1/8 模型,用主从边界(Master/Slave)把周期性接起来。计算量直接降到八分之一,这对参数化扫描来说是救命级别的优化。
但用对称边界有个前提:激励也要对称。如果绕组接法本身不对称,或者你在研究故障工况(比如匝间短路),那对称边界就不能用了。我见过有人为了省时间强行加对称边界,结果算出来的三相电流不平衡,排查了半天才发现是边界条件用错了。
3. 材料库这一步偷懒,后面所有结果都是错的
3.1 硅钢片的 B-H 曲线不能直接用默认值
Maxwell 自带的材料库里有一些通用硅钢片,但如果你用的是特定牌号,比如某厂的 35WW300 或者 50W470,一定要拿到厂家实测的 B-H 曲线和损耗曲线(铁耗的 BP 曲线,也就是不同频率、不同磁密下的单位质量损耗)。默认材料的饱和磁密和损耗系数和实际牌号差得可能很远。
我踩过一次坑:用默认材料算一台电机,铁耗算出来 18W,实测接近 40W。后来换成厂家给的损耗曲线,铁耗算到 36W,和实测就对上了。原因就是默认材料的损耗系数偏乐观。铁耗在 Maxwell 里通常用 Bertotti 模型分离成磁滞损耗、涡流损耗、异常损耗三项,系数都是从损耗曲线拟合出来的,材料不对,系数全错。
3.2 永磁体的退磁曲线和温度修正
永磁体这块,钕铁硼的退磁曲线在常温下是直线,但温度一高,膝点会上移,工作点如果落到膝点以下就会发生不可逆退磁。Maxwell 里可以给永磁体设置温度相关的退磁曲线,做高温工况校核时一定要把温度设对。
具体操作上,在材料属性里把永磁体的矫顽力 Hc 和剩磁 Br 设成温度的函数,或者直接建几个不同温度的材料副本。我一般会算 20℃、80℃、120℃、150℃ 几个点,看最恶劣工况下永磁体工作点离膝点还有多少余量。这个余量我一般留 15% 以上,因为实际工况里局部温度可能比平均温度高不少。
3.3 绕组材料的填充系数与等效
绕组在 Maxwell 里通常不建每一根漆包线,而是建成一个整体区域,然后给一个填充系数(Fill Factor)。这个系数是铜的实际截面积除以槽内总面积。填充系数设得准不准,直接影响铜耗和槽漏抗。手工绕线的填充系数一般在 0.4 到 0.5,机器绕线能到 0.55 到 0.65,扁线电机能到 0.7 以上。
注意:填充系数不仅影响电阻,还影响槽内导热路径。如果你后面要做热耦合,这个值要和你实际工艺对得上,否则热仿真也会偏。
4. 激励与运动设置:瞬态场里最容易出错的地方
4.1 绕组激励的三种给法及适用场景
Maxwell 瞬态场里给绕组加激励,常见有三种方式:外电路(External Circuit)、电流源、电压源。选哪种取决于你要研究什么。
如果你要算反电动势和空载特性,用电流源给 0 电流就行,让转子转起来,看绕组感应出多少电压。如果你要算负载特性和转矩,用电流源给正弦电流最直接,控制好电流幅值和相位角。如果你要研究驱动器的实际影响,比如 PWM 供电下的谐波损耗,那就得用外电路,把逆变器拓扑搭进去。
我个人的习惯是:先用电流源快速扫一遍工作点,把电流角和转矩的关系摸清楚,再用外电路做精细校核。因为外电路求解慢,而且电路参数一多容易出错,不适合做前期探索。
4.2 电流角与 MTPA 的关系
永磁同步电机的电流角(电流矢量相对于 q 轴的角度)直接决定转矩输出。表贴式永磁电机(SPM)的磁阻转矩可以忽略,MTPA 就是电流角为 0,也就是纯 q 轴电流。但内置式永磁电机(IPM)有凸极效应,Ld 小于 Lq,存在磁阻转矩,MTPA 对应一个负的 d 轴电流,电流角在 20 到 40 度之间。
在 Maxwell 里做这个扫描很简单:固定电流幅值,把电流角从 0 扫到 45 度,看平均转矩的最大值出现在哪个角度。这个扫描用 2D 瞬态做,每个角度算一两个电周期,很快就能出结果。我一般会扫 0、10、20、30、40、45 这几个点,先粗定位,再在峰值附近细化。
4.3 运动设置里的 Band 和初始位置
瞬态场要转起来,必须设一个 Band(运动区域),把转子和气隙的一部分包进去,然后给一个转速。这里有两个坑:一是 Band 不能和静止区域的网格重叠,二是初始位置角要和你的转子初始机械角对齐。
初始位置角设错,会导致反电动势相位整体偏移,进而让电流角和转矩的对应关系全乱。我的做法是:先设一个初始位置,跑一个空载,看 A 相反电动势过零点是不是在预期位置,不对就调整初始角,直到相位对上,再去做负载。
5. 网格:不是越密越好,而是要密在刀刃上
5.1 气隙网格的层数与尺寸
气隙是电磁能量交换的核心区域,网格必须够密。Maxwell 里可以用"气隙分层"功能,在气隙里强制生成若干层网格。我一般设 3 到 5 层,每层厚度均匀。层数太少,气隙磁密波形会失真,齿槽转矩算不准;层数太多,网格量暴涨,收益递减。
气隙网格的尺寸,经验上取气隙长度的三分之一到五分之一比较合适。比如气隙 0.5mm,网格尺寸取 0.1 到 0.15mm。这个尺寸下,气隙磁密的波形能比较光滑地还原出来。
5.2 齿尖和槽口的局部加密
定子齿尖和槽口是磁密集中、磁场变化剧烈的地方,这里的网格如果不够密,铁耗和齿槽转矩都会算不准。Maxwell 里可以对这些面或者体单独设网格尺寸,一般取齿尖最小宽度的三分之一左右。
但要注意,局部加密会让网格过渡变差,如果加密区域和周围网格尺寸差太多,求解精度反而会下降。我的做法是设一个网格生长率(Growth Rate),控制在 1.5 以内,让网格从密到疏平滑过渡。
5.3 自适应网格与手动网格的取舍
Maxwell 2D 有自适应网格功能,它会根据场量的误差自动加密,直到收敛。这个功能对新手很友好,但有个问题:自适应网格是基于某个特定工况收敛的,如果你后面要扫很多工况,自适应出来的网格未必对所有工况都合适。
我的习惯是:先用自适应跑一个代表性工况,把网格收敛下来,然后把这个网格固定住,后面所有工况都用这套网格。这样既保证了精度,又保证了不同工况之间的可比性。如果某个工况结果异常,再单独检查网格。
6. 求解设置与收敛:为什么你的仿真跑得慢还不出结果
6.1 时间步长的选择逻辑
瞬态场的时间步长直接决定求解精度和耗时。步长太大,反电动势波形会失真,铁耗算不准;步长太小,求解时间线性增长。经验上,一个电周期至少取 100 到 200 个时间步。如果电机转速高、频率高,步长要相应减小。
我一般会先按每电周期 120 步跑一遍,看波形是否光滑。如果波形有明显锯齿,就加密到 200 步甚至更多。对于要算铁耗的工况,步长还要更小,因为铁耗对磁密变化率敏感,步长太大会把高频分量滤掉。
6.2 非线性收敛与松弛因子
电机材料是高度非线性的,硅钢片的 B-H 曲线在饱和区几乎水平,这会让牛顿迭代很难收敛。Maxwell 里可以调非线性求解的松弛因子和迭代次数。如果发现残差曲线震荡不收敛,可以适当降低松弛因子,或者增加每步的最大迭代次数。
我遇到过一次怎么都不收敛的情况,排查下来是某个区域的网格质量太差,有畸形单元。重新划分网格之后就正常了。所以收敛问题不一定是求解设置的问题,网格质量要先排查。
6.3 仿真跑得慢的几个真实原因
"Ansys Maxwell 仿真很慢"是个高频问题,我总结下来原因通常有这几个:模型里有多余的几何细节(比如倒角、小孔、螺纹),这些对电磁场几乎没影响但会大幅增加网格量;网格全局设得太密,没有做局部控制;时间步长设得太小;用了 3D 但没利用对称性;求解器选了直接求解器而不是迭代求解器。
优化顺序我建议是:先清理几何,把不影响电磁的细节全部去掉;再做网格局部控制,把网格量压下来;然后检查时间步长是否必要那么小;最后考虑对称性和求解器选择。这一套下来,求解时间砍掉一半很常见。
7. 后处理:怎么从一堆场图里读出有用的信息
7.1 反电动势与转矩曲线的提取
空载反电动势是最好提取的:跑一个空载瞬态,在后处理里对绕组区域做磁链对时间的导数,或者直接用感应电压。Maxwell 里可以直接在绕组上定义反电动势的观测量。转矩的话,用虚功法或者麦克斯韦应力张量法,在气隙中间画一条线或者一个面,对转矩做积分。
我一般会把反电动势做 FFT,看谐波含量。反电动势的 THD 直接反映气隙磁密波形的正弦度,THD 高意味着转矩脉动大、铁耗高。如果 THD 超标,就要回去改极弧系数、槽口宽度或者充磁方向。
7.2 查看某一点的场图与场量随时间的变化
"Maxwell 查看某一点的场图"这个需求很常见。操作上,在后处理里用 Field Overlays 或者 Calculator,可以提取任意点的磁密、磁场强度、电流密度等。如果要看某一点随时间的变化,可以在瞬态求解后,用 Fields Reporter 定义一个点,输出该点场量随时间的曲线。
这个功能在排查局部饱和时特别有用。比如你怀疑某个齿尖饱和了,就在那个位置取一个点,看它的磁密峰值有没有超过材料的饱和点。如果超了,就要考虑加大齿宽或者换更高牌号的硅钢片。
7.3 铁耗与铜耗的分离计算
铁耗在 Maxwell 里是后处理算出来的,基于每个单元的磁密波形和材料的损耗曲线。铜耗则是根据绕组电阻和电流算的,注意要把端部绕组的电阻加进去,2D 里默认没有端部。
我一般会把铁耗按定子齿、定子轭、转子分别统计,看哪部分占比最大。如果定子齿的铁耗占比过高,说明齿部磁密偏高,可以考虑优化齿宽或者槽型。铜耗的话,除了直流铜耗,高频工况下还要考虑集肤效应和邻近效应,Maxwell 里可以开绕组的涡流求解来算交流铜耗。
8. 参数化与优化:让仿真真正服务于设计迭代
8.1 参数化扫描的变量选择
Maxwell 的参数化功能可以让你把几何尺寸、材料属性、激励参数都设成变量,然后批量扫描。但变量不是越多越好,我一般会先做灵敏度分析,找出对目标影响最大的几个变量,再针对这几个做精细扫描。
以转矩脉动优化为例,影响最大的通常是极弧系数、槽口宽度、充磁方向这几项。我会先固定其他参数,单独扫极弧系数,看转矩脉动怎么变,找到趋势之后再联合扫描。
8.2 多目标优化的权衡
电机设计很少是单目标优化,通常是转矩、脉动、铁耗、成本几个目标互相打架。Maxwell 本身可以做参数扫描,但多目标优化一般要配合 optiSLang 或者 Workbench 里的 DesignXplorer。
我的经验是,先把目标排序,哪个是硬指标,哪个是软指标。比如转矩是硬指标,必须达标;脉动和铁耗是软指标,在达标的前提下尽量优化。这样优化的时候就有明确的取舍依据,不会陷入"什么都想要"的困境。
8.3 从仿真到实测的闭环
仿真做得再好,最终还是要和实测对上。我一般会在样机测试后,把实测的反电动势、转矩、效率反哺回模型,看偏差在哪里,然后修正模型里的参数(比如填充系数、损耗系数、接触热阻)。这个闭环做几轮之后,模型的置信度会越来越高,后面新项目的仿真结果就敢直接用来做决策了。
9. 几个高频问题的排查思路
9.1 反电动势波形不对
先查初始位置角,再看绕组连接和匝数,最后查网格和气隙。初始位置角错是最常见的,表现为波形整体平移。绕组连接错表现为幅值或相位不对。网格太粗表现为波形有毛刺。
9.2 转矩算出来偏小
先确认电流角和电流幅值给对了没有,再看材料饱和有没有考虑,最后查端部漏磁有没有忽略。2D 里转矩偏小是常态,因为端部效应没算,一般偏小 5% 到 10% 是正常的。
9.3 铁耗和实测差很多
材料损耗曲线是不是厂家实测的,网格是不是够密,时间步长是不是够小,这三项排查完基本就能定位。如果还差,就要考虑加工应力对材料性能的影响,这个在仿真里很难体现,只能靠经验系数修正。
10. 我个人的一点使用体会
Maxwell 这个工具,入门不难,难的是把结果算准。我见过太多人建完模、跑出个云图就交差了,但那个云图到底可不可信,心里没底。我的建议是,每做一个新拓扑,都想办法找一个已知的解析解或者实测数据去对标,哪怕是一个简化的模型。对标的次数多了,你对什么设置会产生什么偏差就有感觉了。
另外,仿真速度这件事,与其抱怨软件慢,不如花时间把几何清理干净、把网格控制好。我做过统计,一个中等复杂度的 2D 瞬态模型,几何清理和网格优化能省下的时间,比换一台工作站还多。工具是死的,人是活的,把功夫下在建模和网格上,比什么都强。