第三章 退磁仿真前处理:几何模型建立、网格划分技巧与边界条件设置

各位工程师朋友,大家好。今天我们进入退磁仿真的核心环节——前处理。说实话,很多人在退磁仿真上栽跟头,不是算法不对,而是前处理没做到位。我自己带团队时,经常跟新人说:前处理花的时间,决定了仿真结果的靠谱程度。今天我就把这些年积累的经验,尤其是永磁体区域网格划分的窍门,跟大家好好聊聊。

退磁仿真前处理 几何模型建立 网格划分技巧 边界与激励源 永磁体建模 气隙处理 倒角与简化 永磁体加密 气隙层网格 网格质量检查 主从边界 电流激励 退磁电流波形

3.1 几何模型建立——别小看这一步

几何模型是仿真的地基。地基歪了,上面盖的楼再漂亮也没用。我见过不少工程师,从CAD导入模型后直接开算,结果算出来的退磁率跟实测差了20%以上。后来一查,原来是模型里永磁体的倒角没处理,导致局部磁场集中。

永磁体建模的几个要点:

  • 尺寸精度:永磁体的长宽高尺寸必须精确到0.01mm级别。特别是厚度方向,直接影响退磁场的计算。我个人习惯在JMAG里重新建模,而不是直接导入CAD,因为JMAG的参数化建模更方便后续调整。
  • 充磁方向定义:这个太重要了。退磁仿真本质上是看永磁体在反向磁场下的表现,充磁方向错了,结果全错。JMAG里用Magnetization Direction来定义,我建议用局部坐标系来指定,尤其是对于瓦形磁钢。
  • 气隙模型:气隙必须单独建模。为什么?因为退磁时气隙中的磁场分布会剧烈变化,没有气隙层,你根本算不准。我一般会在转子外表面和定子内表面之间留出0.5~1mm的气隙层。
💡 我的经验: 对于内置式永磁电机(IPM),永磁体槽的倒角不要忽略。我曾经有一个项目,忽略了0.2mm的倒角,结果退磁仿真显示安全,但实测退磁率高了8%。后来加上倒角,仿真结果就跟实测对上了。

3.2 网格划分技巧——永磁体区域是重中之重

网格划分是退磁仿真里最考验功力的环节。说白了,网格就是仿真的「分辨率」。分辨率不够,细节就丢了。尤其是永磁体区域,退磁现象往往发生在局部,网格太粗根本捕捉不到。

永磁体区域网格划分的核心原则:

  1. 至少3层网格:在永磁体的厚度方向(充磁方向),我要求至少划分3层网格。如果厚度是5mm,每层网格控制在1.5mm以内。为什么?因为退磁时磁场在厚度方向上有明显的梯度变化,一层网格根本反映不出来。
  2. 局部加密:永磁体的四个角部是退磁的高风险区域。我会在这些角部设置局部网格加密,网格尺寸缩小到整体网格的1/3~1/5。JMAG里用Local Mesh Control来实现。
  3. 气隙层网格:气隙层至少划分2~3层网格。气隙中的磁力线分布对退磁影响很大,网格太粗会导致磁阻计算偏差。
🔧 操作技巧: 在JMAG中,我习惯先用Auto Mesh生成基础网格,然后手动调整永磁体区域的网格参数。具体操作:右键点击永磁体部件 → Mesh Settings → 设置Element Size为整体网格的0.3倍,并勾选Curvature Refinement

网格质量检查:

网格画完了,别急着算。先检查质量。我一般看三个指标:

指标 合格标准 我的建议
长宽比 (Aspect Ratio) < 5 永磁体区域控制在3以内
雅可比 (Jacobian) > 0.7 低于0.5的单元必须重画
最小角度 > 20° 三角形单元最好大于30°
⚠️ 避坑指南: 我曾经有一个项目,永磁体网格长宽比达到了8,仿真结果显示退磁率只有2%,但实测到了7%。后来重新划分网格,把长宽比控制在3以内,结果就跟实测吻合了。所以网格质量真的不能马虎。

3.3 边界条件与激励源设置

边界条件和激励源,是仿真的「剧本」。剧本写错了,演员再努力也白搭。退磁仿真中,边界条件主要影响磁场分布,激励源则决定了退磁的「力度」。

边界条件设置:

  • 主从边界 (Master/Slave):对于周期对称模型,用主从边界可以大幅减少计算量。我一般取一个极距作为仿真模型,然后在两侧设置主从边界。注意:主从边界的网格节点必须一一对应,否则会报错。
  • 自然边界 (Natural Boundary):在模型的外围(定子外径和转子内径),设置磁通平行边界(即磁力线沿边界切线方向)。这个设置比较保守,但稳定。
  • 对称边界 (Symmetry):如果电机是偶数极,可以利用对称性只仿真一半或四分之一模型。但要注意:退磁现象在对称位置可能不完全一致,尤其是存在制造公差时。

激励源设置——退磁电流的加载:

退磁仿真最关键的激励源,就是退磁电流。这个电流通常比额定电流大很多,可能是额定电流的3~5倍,甚至更高。

在JMAG中,我一般这样设置:

// 电流激励设置示例(JMAG脚本)
// 三相电流,峰值电流设为额定电流的4倍
double I_rated = 100;  // 额定电流有效值
double I_peak = I_rated * 4 * sqrt(2);  // 峰值电流

// A相电流
CurrentA = I_peak * sin(2*PI*Frequency*Time + PhaseAngle);
// B相电流
CurrentB = I_peak * sin(2*PI*Frequency*Time + PhaseAngle - 2*PI/3);
// C相电流
CurrentC = I_peak * sin(2*PI*Frequency*Time + PhaseAngle + 2*PI/3);
💡 关键点: 退磁电流的相位角(PhaseAngle)非常关键。不同的相位角会导致退磁磁场的方向不同,从而影响退磁区域。我一般会扫描0°~360°的相位角,找出最恶劣的退磁工况。这个扫描过程虽然耗时,但值得做。

退磁电流波形的选择:

退磁电流可以是直流、交流或脉冲波形。对于电机控制器短路故障导致的退磁,我通常用交流电流,频率与电机转速对应。对于外部磁场干扰导致的退磁,用直流电流更合适。

我个人习惯在JMAG的Circuit Editor里搭建一个简单的电路模型,包含三相绕组和电流源,这样更直观。如果直接用Current Source设置,记得勾选Include Skin Effect,因为高频退磁电流的集肤效应会影响绕组电阻。

🔧 实用建议: 在设置激励源之前,先跑一个空载仿真,确认永磁体的空载磁场分布是否正确。如果空载磁场都不对,那退磁仿真就更不用说了。这一步虽然简单,但能帮你排除很多低级错误。

好了,关于退磁仿真的前处理,几何模型、网格划分、边界条件和激励源,我就讲到这里。这些内容看起来琐碎,但每一个细节都可能决定仿真的成败。希望大家在实际操作中,多留个心眼,尤其是永磁体区域的网格,多花点时间绝对值得。

公众号:蓝海资料掘金营,微信deep3321