前阵子给一台平板型永磁同步直线电机做方案验证光是Maxwell仿真就跑了大半个月。绕来绕去最核心的其实就三件事空载反电动势、推力输出、磁阻力端部力和齿槽力。这三关过了电机的电磁设计基本就心里有数了这三关里任何一关出现幺蛾子后面做样机必然返工。今天把这套流程完整捋一遍包括平板型直线电机建模时的注意事项、运动设置和Band怎么处理、电流激励怎么给、以及磁阻力那三个力怎么分离出来准备做直线电机预研或者课程设计的朋友可以直接照着操作。先交代一下我用的仿真环境Maxwell瞬态求解器二维模型为主三维模型只在需要校验端部效应时才跑。整个流程不复杂但每一步都有讲究尤其是运动体设置和电流相位做错了波形和力曲线看起来就是不对而且你还找不到原因。1. 仿真前的整体思路模型准备与求解器选型1.1 为什么必须用瞬态场而不是静磁场很多刚接触直线电机仿真的人第一反应是用静磁场求解器因为反电动势、推力、磁阻力这些名字听起来像是某个状态下的结果。但实际上这三样东西全部是位置和时间的函数必须让动子真的动起来才能算出来。静磁场求解器只能给你某个固定位置下的磁链、受力给不出波形更给不出力波动。所以主求解器直接选Transient瞬态场。瞬态场里有两套东西是静磁场没有的一是运动设置MotionSetup二是绕组激励按时间变化。这两个配合起来才能模拟直线电机在运动过程中磁场分布和绕组电压电流的实时变化。顺带说一句三维模型当然也能做而且从磁路对称性来说三维比二维准但三维瞬态模型的计算量不是一般的大调试一次要几个小时。我做这个方案时先用二维模型把所有参数和波形调对最后才在三维里跑了一个工况验证端部效应带来的差别实际发现二维计算结果已经足够指导方案决策。1.2 平板型直线电机的建模参数与结构选择平板型直线电机说白了就是把旋转永磁电机沿径向切开、展平。展平之后问题就来了旋转电机磁路是闭合的直线电机两端是断开的原本均匀的磁场在两头出现畸变这就是磁阻力和推力波动的根源。建模前要把结构形式定下来我这次用的是动初级、定次级的结构也就是初级铁芯带着三相绕组运动次级是贴在背铁上的永磁体固定不动。这种结构在直线电机里最常见仿真设置也最直观。示例电机的参数如下表后面所有仿真都是基于这个模型展开。参数项数值说明额定线速度1 m/s运动方向沿X轴极距20 mm相邻永磁体中心距永磁体材料N35剩磁1.2T左右永磁体厚度3 mm充磁方向垂直气隙气隙长度1 mm机械气隙初级铁芯9槽8极分数槽集中绕组绕组层数双层每槽上下两层每线圈匝数50串联匝数视接法折算电流幅值5 A推力校核工况模型深度80 mm等效叠厚2D仿真的模型深度必须在求解器中单独设置Maxwell默认按1米深度计算如果你不改成实际叠片厚度反电动势和推力会整体偏大几十倍第一次跑的时候差点被这个坑误导。1.3 MotionSetup与Band区域的处理细节Band是瞬态场里最核心的运动容器。直线电机仿真中所有运动部件必须完全包在一个叫Band的对象内Band内部会自动重新剖分网格动子每走一步网格跟着更新。如果Band设置不对整个仿真直接废掉。Band的几何轮廓可以在建模时直接从初级铁芯轮廓偏移出来四周留出至少2到3毫米的余量。注意Band不能和任何静止物体相交尤其不能让永磁体或背铁穿进Band区域否则运动体定义会报错。我做第一个模型时就是Band做小了与气隙处永磁体的边界几乎碰上结果求解器提示网格自交排查了半天才发现是这个原因。MotionSetup的速度设定也有讲究。空载反电动势和磁阻力仿真需要动子匀速运动直接设恒定速度即可比如1 m/s推力仿真如果只需要电磁推力同样用恒定速度模式不必让机械瞬态参与求解。只有当你想模拟带负载启动过程、速度真实变化时才需要勾选机械瞬态并设置负载力、质量、阻尼等参数。2. 空载反电动势仿真核对接线、磁链与电压波形2.1 空载条件的实现方式空载反电动势的实质是电枢电流为零动子运动时绕组中产生的电压。在Maxwell里怎么实现电流为零很关键直接把绕组激励删掉是不行的因为悬空的绕组没有激励源求解器不知道该怎么处理。正确做法是把三相绕组都设为电流源幅值为0A。为什么要用电流源而不是电压源如果用电压源0V理论上绕组短路了那测出来的是短路电压波形完全不对。用0A电流源表示开路绕组两端电压就是真正的反电动势。这个细节我见过不止一个人搞反然后拿到的波形怎么看怎么怪。速度设1 m/s后电频率可以直接算出来f v / (2 * 极距) 1 / (2 * 0.02) 25 Hz极距是20mm一个N极加一个S极的总长度是40mm磁场变化一个完整周期频率就是25Hz。求解时间建议至少取0.12秒也就是三个完整电周期方便后处理时去掉开头暂态、看稳定波形。时间步长取0.0004秒这样每个周期有100个采样点波形足够光滑。2.2 反电动势波形的读取与理论校核求解完成后在后处理中绘制A、B、C三相绕组的电压波形因为电流为0电压就是反电动势。正常结果应该是三相对称的正弦波幅值一致相位互差120度。我这次跑出的相电压幅值约8.7V也就是有效值约6.1V。拿到仿真波形后一定要和理论值对账否则算完心里没底。直线电机每相空载反电动势有效值的估算公式可以写成E 4.44 * f * N * k_w * Φ其中N是每相串联匝数k_w是绕组系数Φ是每极磁通由永磁体磁密和气隙尺寸决定。把示例电机的参数代进去算出来约6.3V和仿真值接近误差在合理范围内。如果偏差超过10%就要回头检查永磁体充磁方向、材料剩磁设置或者模型深度这些基础参数。这里有一个隐含判断20mm极距配合25Hz频率反电动势一个周期的机械位移是40mm和一个极对长度完全对应。如果波形周期不对说明速度和极距换算有问题。2.3 反电动势波形质量能反映什么问题反电动势波形不只是用来对电压数值的它的形状直接告诉你这台电机的电磁设计合不合理。正弦度越好配合正弦电流驱动时的推力波动越小。如果波形削顶或者明显畸变通常是极弧系数选得不好、永磁体形状不理想或者初级铁芯存在局部饱和。三相幅值不相等又是一个常见问题多半是绕组设置不对称、初始终端位置选得不好或者模型本身不具备完全三相对称条件。对于9槽8极这种分数槽结构三相绕组本身是严格对称的如果波形看到幅值差异优先查激励设置里的绕组相序和线圈匹配关系。波形平顶、过零点偏移也和齿槽有关。齿槽的存在会让反电动势波形带上固定的齿槽谐波成分所以看到反电动势谐波偏大时不要急着改绕组先想想是不是齿槽效应导致的。这些异常在后面的磁阻力仿真里会看得更清楚。3. 推力输出仿真让电流和运动协同起来3.1 三相对称电流激励的设置要点推力仿真的核心变化是把空载时的0A电流改为三相对称正弦电流。对表贴式永磁直线电机来说最简单的控制策略是保持d轴电流为零也就是让电流和反电动势同相位。这样搞电磁推力基本与q轴电流成正比理解起来直观仿真设置也方便。电流激励写成这样IA 5 * sin(2 * pi * 25 * time) IB 5 * sin(2 * pi * 25 * time - 2 * pi / 3) IC 5 * sin(2 * pi * 25 * time 2 * pi / 3)单位是安培频率还是25Hz和运动速度严格对应。这里最容易翻车的是相位对齐电流频率必须与反电动势频率一致初相角最好先通过空载反电动势波形确定。我通常先看A相反电动势在time0处是上升过零还是下降过零然后据此设置电流初相角确保初始位置在q轴上。如果相位反了推力会变成负的如果相位差不是90度而是某个任意角平均推力会明显偏小。3.2 推力波形的获取与结果判定在瞬态场的力求解设置中选中所有运动部件初级铁芯、绕组铜线求解X方向的受力。注意绕组也要一起选因为通电线圈在磁场中同样受力只选铁芯会漏掉一部分电磁力。正常结果里推力曲线经过短暂的暂态波动后进入稳定平均值对应电机的稳态电磁推力。我跑的5A工况稳定后平均推力约27N折算下来推力常数约5.4N/A符合这类小推力直线电机的典型水平。推力不是一条光滑直线一定有波动。波动幅度用峰峰值除以平均值来衡量我这次跑出来约正负3%左右。如果波动太大优先怀疑磁阻力叠加和反电动势谐波两个方面。3.3 推力波动与磁阻力的关联解读负载推力曲线里叠加了两种东西电流和磁场相互作用产生的电磁推力分量以及永磁体存在的状态下本来就有的磁阻力波动。前者由电流和反电动势波形决定后者和空载磁阻力完全相同。判断推力波动来自哪边的办法把同一速度下空载测出的磁阻力曲线叠到推力波动曲线上看趋势如果波动主频率和相位与磁阻力高度吻合说明波动主要是磁阻力贡献的如果磁阻力很小但推力波动依然明显大概率是电流谐波和反电动势谐波相互作用的结果。我在实际项目里通常做一组电流幅值扫描从1A到10A看平均推力和波动比例怎么变。推力波动占比随电流增大会有所下降因为电磁推力基底大了同样的磁阻力波动被稀释了。但如果电流大到铁芯饱和波动又会重新变大那就是饱和导致的力谐波了。4. 磁阻力、端部力、齿槽力的仿真分离4.1 三个力的物理来源与关系磁阻力在直线电机里是没有电流时的定位力动子静止时想把动子从某个位置拉开需要的力或者匀速运行时动子必须额外克服的阻力。它不是一个单一来源的力而是端部力和齿槽力两个分量的合成。齿槽力来自初级铁芯开槽后气隙磁导不均匀永磁体总是趋向于停在磁阻最小的位置于是产生一个随位置周期变化的切向力。端部力则是直线电机特有的东西初级铁芯两端是断开的两端齿部面对的磁路和中间完全不同磁力线进出端部产生的附加轴向力。旋转电机没有这个力所以很多只做过旋转电机仿真的人第一次看到磁阻力波形里带有长周期分量时会一脸懵。4.2 分离磁阻力三个分量的实操思路很多人问我怎么在仿真里单独看齿槽力、单独看端部力。其实不需要什么高级操作关键在于建辅助模型做减法。先跑基准模型真实带齿槽的有限长初级空载、电流为0动子匀速运动测得的X方向力就是总磁阻力。再复制一个项目把初级铁芯表面处理成等长的光面铁芯也就是去掉齿和槽保留同样的铁芯长度其它条件不变再跑一次空载。光面铁芯没有齿槽效应测出来的力就是端部力。用总磁阻力减去端部力剩下的就是齿槽力。这个分离方法在工程上是够用的。要注意精度问题当磁路饱和较强时总磁阻力等于端部力加齿槽力的叠加关系会打折扣因为铁芯饱和会改变磁路的等效结构。但针对平板型直线电机的正常工况这个方法给出的分量趋势和主阶次是可靠的。还有一个细节端部力对初级两端的形状极其敏感所以辅助模型的光面铁芯长度必须和真实模型一致端面形状也要尽量接近只是齿槽特征被抹平。如果光面模型只是简单拉了一个平板铁芯端部力会失真。我当时是先复制初级铁芯然后把槽开口区域填平再删掉齿槽的细小几何特征这样端部结构最接近真实情况。4.3 磁阻力波形处理与削弱手段验证分离后的磁阻力曲线放到同一张图里能看出来不同的波动周期。端部力的变化周期往往较长与极距或极对数的周期相关齿槽力的变化周期较短在一个极距里会出现多个波动周期。各个分量的幅度和阶次用FFT分析一下一目了然。削弱磁阻力的常用手段有几种分数槽集中绕组本身就能有效降低齿槽力9槽8极就是典型调整极弧系数也能改变齿槽力幅值一般扫个0.7到0.9就能看到明显差异端部力则要靠辅助齿、端部倒角或调整永磁体端部形状来削弱。仿真验证削弱效果时有个技巧不要只看总磁阻力曲线把分离出的端部力和齿槽力分别列出来看。有些措施对齿槽力有效但对端部力无效甚至可能放大端部力只看总曲线容易被误导。我当时试了两种辅助齿形状总磁阻力峰峰值相差10%左右但拆开看发现一个方案把端部力压低了另一个方案把齿槽力压低了选择哪种取决于你要优化哪个频率段的波动。5. 平板型直线电机仿真常见问题与排查手册5.1 建模阶段的常见错误建模阶段最常见的是Band问题和永磁体充磁方向问题。永磁体充磁方向设反了反电动势会变成零或者波形完全反向。检查方法是看磁力线分布图正常N极和S极交替指向气隙。Region边界也要花心思处理。直线电机模型外边界通常给气球边界但如果只建半模型想节省计算时间要注意永磁体交替极下对称条件不是普通的奇对称设置错了结果会非常奇怪。建议初学阶段老老实实建全模型算力能承担就不要省。5.2 求解设置与网格剖分的坑气隙网格是瞬态场仿真的生命线。气隙里至少要保证两层以上网格永磁体内部和初级齿尖也要加密否则反电动势波形会出现毛刺推力曲线上会有非物理的高频抖动。气隙加密让人感觉像在等批处理跑完但为了提高力计算的稳定性这个成本必须掏。时间步长过大同样会让波形锯齿化。电周期25Hz步长如果取0.002秒一个周期只有20个点FFT做出来全是假的谐波。一般来说每周期50到100个点是底线再多就是纯增加计算量了。5.3 结果异常速查表异常现象可能原因处理办法反电动势幅值偏小永磁体剩磁设置低模型深度未改匝数不对核实材料参数检查模型深度对照绕组设置反电动势三相不对称绕组相序错误初始位置不当线圈连接关系有误检查三相绕组定义调整初始位置推力平均值接近零电流相位与反电动势相位没对齐电流频率与速度不匹配空载跑一遍再定电流初相角核对频率推力曲线高频抖动气隙网格太粗时间步长过大加密气隙网格减小步长磁阻力曲线不是以零为中心端部结构不对称永磁体长度方向有限导致端部力直流分量检查模型端部用长次级模型仿真中途网格自交Band设置与运动部件相交Band过小调整Band几何确保运动件完全包含且不碰静止件5.4 计算资源与效率平衡建议二维模型单次瞬态计算从几分钟到半小时不等主要取决于网格数量和时间步数。第一次做可以把网格设粗一点跑通流程确认激励、Band、后处理都正常再把网格加密跑最终结果。这样做省下来的调试时间非常可观。如果要做参数扫描比如扫极弧系数或者永磁体厚度用Maxwell的参数化功能一次性提交不要手动改一次跑一次。扫描模型建议用相对稀疏的网格得到趋势后再对关键参数点加密验证。6. 个人实操体会与一点补充做这一类直线电机仿真我最大的体会是波形和力曲线是否合理永远要用理论公式先做预判再拿仿真结果核对。不管是反电动势的有效值估算、频率和速度的换算还是推力常数和反电动势常数的量级先算出个大概范围再去看仿真结果思路会清晰很多也不会被一些无关紧要的数值波动带偏。另外一个小技巧后处理阶段不要只看瞬时波形把FFT分析用起来看各阶谐波分量的大小。直线电机推力波动、噪声振动的问题到最后都要靠谐波分析定位到具体阶次养成这个习惯后面会省很多事。这套平板型直线电机的仿真流程我后续还会用在带负载的动态启动过程、多物理场耦合分析上。目前的这套组合拳已经把电磁设计阶段的验证问题解决得差不多了剩下的更多是细化优化层面的工作。