很多第一次接触变压器的人都会问变压器运行时为什么总在“嗡嗡”响工程圈里经常用一句话形容这种现象——磁场在“唱歌”。这歌声一点都不浪漫本质是铁芯在交变磁场下发生磁致伸缩加上绕组和结构件受到周期电磁力最后通过空气辐射形成了可听噪声。如果你正被变压器仿真搞得头大又恰好被这种噪声困扰那这篇手记应该能帮你理清思路。我最近做了一台小型干式变压器的 Maxwell 3D 铜片绕组损耗分布仿真过程中把铁芯磁密、铜损热点、电磁力分布以及噪声源头全部过了一遍。这篇文章不打算写成软件说明书而是按我实际的建模思路、参数设置、结果判读顺序来写顺便把踩过的坑都交代清楚。内容面向做变压器、电感、电抗器研发的同行以及正在自学电磁场仿真的学生朋友。你会看到我为什么选 3D、为什么关注铜片损耗云图、怎么从力密度判断“唱歌”位置也会看到一张可以直接抄作业的参数对照表。1. 为什么变压器的磁密分析这么重要从“唱歌”说起1.1 磁致伸缩与噪声产生的物理机理所谓磁致伸缩是指铁磁材料在磁化过程中沿磁化方向发生微小形变。硅钢片的磁致伸缩率通常在 1×10⁻⁶ 到 3×10⁻⁵ 量级这个变形量虽然极微但当它作用在整台变压器的铁芯上而且以 50 Hz 频率反复施加时铁芯就会像喇叭纸盆一样振动起来。这里有个特别容易忽略的细节电磁力与磁密的平方成正比。50 Hz 的交变磁场每秒钟正负翻转 50 次磁密平方之后频率翻倍变成 100 Hz。所以变压器噪声的基频基本是 100 Hz而不是 50 Hz听感上就是那种比较沉闷、连续的“嗡”声。如果电网里还有谐波或者变压器出现了直流偏磁噪声里还会叠加高次谐波和偶次分量声音会变得更尖锐、更复杂。要回答“这台变压器为什么吵”光靠耳朵听不出来只能靠仿真把磁密分布、力密度分布定量化。这也是我在这类项目里坚持先跑电磁场仿真的原因如果没有数据支撑结构降噪就是瞎蒙。1.2 磁场仿真能真正量化哪些工程指标我理解的变压器磁场仿真最终交付的不只是一张漂亮的磁密云图而是几个能直接指导设计的数值指标。第一个是铁芯磁密。磁密越高铁芯体积越小但饱和风险越大。通过仿真可以定位窗口拐角、T 型接缝处的局部高磁密点判断设计是否留足裕量。第二个是损耗分布。铁芯损耗包含磁滞损耗和涡流损耗绕组损耗则主要是铜损。在铜片或铜箔绕组结构中涡流效应会导致电流密度不均匀损耗热点往往集中在端部和边角。这类信息只有三维仿真才能给准。第三个是电磁力分布。力密度高的区域就是振动和噪声的主要激励源后续做结构谐响应、声学分析时都要从这里取载荷。第四个是杂散磁通。漏磁会穿入外壳、夹件、螺栓等金属结构引起局部涡流和发热。仿真能提前发现哪里有“看不见的磁场漏出去”也能为工频磁场抗扰度等测试提供参考。这些指标只要一个没量化设计阶段就容易埋雷。等到样机做出来再测噪声、测温升返工成本往往非常高。1.3 为什么选择 Maxwell 3D 做变压器电磁场仿真低频电磁场仿真的工具不止一种但 Maxwell 3D 在变压器领域用的人确实最多资料也最全。网上搜“maxwell3d 铜片变压器损耗分布仿真图文教程”基本都能找到可参考的案例。我选它的原因主要有三个。第一它有成熟的涡流场和瞬态场求解器。涡流场适合线性频域问题损耗、电感、阻抗都能直接算瞬态场适合非线性、谐波、直流偏磁和磁场动画这类时间相关问题。两种求解器覆盖了变压器仿真的绝大多数场景。第二它和结构、声学仿真工具之间有成熟的耦合接口。电磁力密度可以导到结构软件里做振动响应再把振动结果送到声学模块计算声压级形成完整的“电磁→结构→声学”链条。第三它支持从 2D 到 3D 的渐进式验证。先用 2D 快速看趋势再用 3D 精确定位端部效应工作流很顺。不过 3D 仿真计算量大、建模细节多这也是有些资料只给雏形不给完整参数的原因。下面我把我实际用的建模过程拆开来说重点讲那些不写进官方教程的关键点。2. 仿真模型搭建的完整思路2.1 几何建模让铜片绕组真正“片片可见”我这次仿真对象是一台小型干式变压器铁芯是 EI 型结构绕组采用铜片铜箔绕制。铜片绕组不像普通圆导线那样能简化为一个圆柱体因为铜片在厚度方向上的涡流分布不均匀必须把每一层铜箔建出来。铜在 50 Hz 工频下的趋肤深度大概可以按 δ sqrt(2ρ / (2πfμ₀)) 估算铜的电阻率约 1.7×10⁻⁸ Ω·m算下来大概是 9.3 mm。如果铜片厚度在 0.3~0.5 mm远小于趋肤深度理论上可以不必再细分厚度方向。但实际中铜箔边角处磁场畸变严重涡流回路会在狭窄区域内集中所以我仍然在铜箔的端部和转角处做了更细的网格剖分。几何建模时还有一个容易出问题的点层间绝缘层。铜箔层之间有绝缘材料如果只建铜箔不建绝缘层会让层间接触关系失真。尤其是并联铜箔之间如果有电位差绝缘层对电流分配的影响是不可忽略的。我的建议是在三维软件里把铁芯、铜箔层、绝缘层、骨架、空气域都建出来再导入 Maxwell。空气域要包住变压器整体尺寸至少是模型外轮廓的 1.3 倍否则漏磁会被截断损耗和电感都会失真。2.2 材料与 B-H 曲线硅钢片不是一根直线材料设置是仿真精度最大的影响因素。很多人随便选一个内置钢材就开跑结果磁密和损耗完全对不上。硅钢片材料必须输入实测的直流 B-H 曲线而且分方向。硅钢片在轧制方向的磁导率明显高于垂直方向如果单用一条各向同性曲线铁芯窗口和接缝处的磁密分布会偏。Maxwell 里可以用各向异性材料定义幅度比值按厂家数据给定通常轧制方向和非轧制方向的磁导率差异在 2~10 倍之间。叠片系数也要处理。硅钢片表面有绝缘涂层叠压后有效导磁截面积小于几何截面积。我习惯在材料里直接设叠片系数 0.96并选择叠片方向垂直于铁芯截面。如果不设叠片系数铁芯磁密会被低估损耗会被高估结果看起来“更保守”但实际上是错的。损耗曲线也要给全。铁芯损耗需要输入不同频率下 B-P 曲线Maxwell 才能用插值算出磁滞损耗和涡流损耗的分离值。只给 B-H 曲线不给 B-P 曲线铁耗就是拍脑袋。2.3 激励与求解器电流源还是电压源变压器仿真最常见的激励方式有两种电流源和电压源。我多数情况用电流源因为电流源的收敛性远好于电压源对线圈内部阻抗不敏感设置也简单。如果是涡流场求解器激励频率设为 50 Hz绕组的电流大小按额定工况给。注意线圈必须定义成 Winding并在 Coil 里设置导体匝数、极性以及串联/并联关系。铜片绕组层与层之间如果是并联的每层都要独立建模并指定到同一个 Winding 里否则电流分配失真损耗分布自然也不对。如果是瞬态场求解器时间步长一般取 0.5 ms 以内一个工频周期 20 ms 就有 40 个解。瞬态场建议至少跑两个周期第一个周期用于消除初始暂态第二个周期数据才稳定。当然瞬态场也支持电压源但需要附加外电路里的绕组电阻和漏感参数怎么给都得在电路里补全比电流源繁琐得多。如果只看稳定损耗和谐波涡流场就够如果关心直流偏磁、磁饱和波形、磁场动画那必须上瞬态场。这两种求解器的选择直接决定后面能不能看到你想看的那些量。3. 关键结果解读损耗分布、磁密与“磁场唱歌”位置3.1 铜片损耗分布如何快速定位发热区求解完成后第一步是查看铜片绕组的欧姆损耗体密度也就是 Ohmic Loss。把云图显示切到铜箔层实体上正常情况下可以看到损耗密度在铜片端部和转角位置高于中部这是因为端部磁场的切向分量更大涡流回路集中在局部。我这次仿真里铜片外侧拐角的损耗密度比中间区域高出大约 60%这个结论直接影响了后续设计铜片端部倒角并把绝缘层的厚度稍微加大损耗热点就明显缓和。还有一个经验在调网格之前不要急着分析损耗值。先用默认网格跑一遍记下最大体积损耗密度的位置然后局部加密该区域网格再跑一遍如果最大值变化超过 10%说明前一轮网格不够粗。有些产品需要进行倍频或谐波仿真比如在变频器输出侧使用的变压器。此时高频分量的趋肤深度会小很多铜箔厚度方向的网格必须加密到高频趋肤深度的 1/3 以下否则高频损耗会被严重低估。对 1 kHz 以上的谐波这个要求很苛刻建议先算一次谐波分量再叠加各阶损耗。3.2 铁芯磁密云图饱和不是“全饱和”而是“局部饱和”铁芯磁密看两个东西一个是云图整体分布一个是最大磁密数值。硅钢片饱和点一般看 1.7~1.8 T正常设计工作磁密取 1.3~1.5 T留一定的饱和裕量。我这次仿真额定工况下铁芯最大磁密约 1.46 T出现在窗口内角处。这个位置因为磁路截面积突然变化磁力线拥挤饱和风险最高。外角区域只有 1.1 T说明整体磁路还有余量但窗口内角是明确的薄弱点。如果看到局部磁密超过 1.8 T不用急着改整个铁芯尺寸先看那个高磁密区域是不是几何尖角或者网格奇异点造成的。如果是真实的磁路过窄那就调整窗口尺寸或者提高硅钢片牌号。判断依据很简单把该区域网格加密两轮如果最大磁密变化不大说明是真实饱和如果明显下降说明是网格伪影。瞬态场比涡流场更容易看出饱和问题因为饱和是随瞬时电流变化的。用电流源设置额定励磁电流一个周期内不同相位下磁密峰值位置会移动这种情况在涡流场里看不到必须用瞬态场的动画功能来看。3.3 电磁力分布噪声源到底在哪要说清噪声最终要落到力上。Maxwell 3D 可以输出 Maxwell 应力张量也就是磁曳力密度还能同时算通电导体在磁场中的洛伦兹力密度。两者叠加就是作用在铁芯和绕组上的总电磁力密度分布。“唱歌”的位置多半就是力密度最高的位置。我这次仿真里铁芯接缝区域和铜片绕组端部的力密度明显比其它位置高且呈周期性脉动频率正好是 100 Hz。接缝处磁拉力方向交替变化让铁芯像在被反复拉伸和挤压绕组端部则是轴向力集中会推动绕组整体振动。要看力密度可以在场结果里选择 Force Density直接加载到模型表面或体上观察其大小和方向。云图上颜色越深的地方就是设计需要重点关注的地方。如果要继续做声学预测把力密度数据导到结构有限元软件里做谐响应分析再把振动速度作为边界条件给声学模块就能计算声压级。这条链路虽然长但每一步都非常明确。3.4 让磁场“动”起来瞬态时序后处理磁场做动画是非常直观的后处理适合内部评审也适合给客户演示。用瞬态场求解器跑完一个周期后右键点击场图选择 Animate设定帧数和时间范围就能导出一段磁密矢量随相位变化的动画。动画里能看到磁通沿铁芯循环推进的过程。和网上常见的“同步电机磁场变化动画视频”不同变压器磁路是脉动磁通而不是旋转磁通因此磁场方向反复正反变化。这正好能说明为什么变压器铁芯一般做成闭合回路而同步电机需要气隙旋转磁场两者磁场形态完全不同仿真表现也完全不同。做动画时建议把时间步长控制在 0.5 ms 以内帧数取 40 帧以上否则磁通推进过程看起来卡顿。如果只关心某一相位的瞬时分布也可以不用动画直接选择对应时间点导出云图。4. 实操过程与核心环节实现4.1 从零跑通 Maxwell 3D 变压器损耗仿真的 7 个步骤给第一次用 Maxwell 3D 做变压器仿真的朋友我把操作流程拆成了 7 个步骤按这个顺序走基本不会乱。打开 Maxwell 3D选择求解器。若分析 50 Hz 稳态损耗选 Eddy Current若分析饱和、谐波或瞬态特性选 Transient。导入几何模型。从三维软件导出的格式建议用 STEP 或 SAT注意单位要统一成毫米。模型里包含铁芯、铜片层、绝缘层和空气域。赋材料。铁芯用各向异性硅钢片设置叠片系数铜片用纯铜绝缘层用塑料类材料。材料库缺失的手动输入 B-H 曲线和 B-P 曲线。创建绕组和线圈。每个绕组定义 Winding在 Coil 中设置匝数、电流方向和并联关系。设置激励。给每个绕组加电流或电压激励频率设为 50 Hz。涡流场加正弦稳态电流瞬态场定义随时间变化的电流波形。设置边界和网格。铁芯外表面用默认自然边界外侧空气域设置气球边界铁芯和铜片局部加密网格。求解。先跑一遍默认网格记录最大磁密和最大损耗密度然后局部加密关键区域再跑一遍对比两次结果差异。第 7 步的“加密-对比”是关键千万别省。很多新手跑一遍就出报告结果磁密最大值、铜损最大值全被网格误差带偏后续设计判断全部没有意义。4.2 关键参数设置对照表下面这张表是我这次仿真实际使用的关键参数可以复制过去再按自己产品调整。参数项推荐设置说明求解器Eddy Current / Transient稳态损耗用涡流场瞬态分析用瞬态场频率50 Hz工频变压器实际按产品需求调整铁芯材料各向异性硅钢片 叠片系数 0.96轧制与非轧制方向磁导率分开铜箔材料铜电阻率 1.7×10⁻⁸ Ω·m按工作温度修正电阻率空气域模型外扩 30% 以上避免漏磁被边界截断网格铁芯最大边长 2~4 mm窗口内角和接缝局部加密网格铜箔端部最小边长 0.3~1 mm保证涡流分布分辨率激励电流源额定电流收敛性好推荐首选时间步长瞬态≤0.5 ms一个周期 40 步以上求解迭代最大 15 步误差 1%收敛后检查能量误差这些值不是恒定不变的。比如铁芯更大网格可以放宽铜箔更薄网格可能要更细。核心原则是在关键场变化剧烈的区域加密在均匀区域放松这样才能把计算资源用在刀刃上。4.3 从求解时间到收敛控制3D 仿真别硬跑3D 仿真最现实的障碍是计算时间。我这次模型大约三百万网格单元单台 16 核工作站跑了接近两个小时。如果不做任何对称简化一个复杂三相变压器的 3D 瞬态场算上几天也不奇怪。所以第一件事是找对称面。铁芯和绕组结构如果存在对称性可以切 1/2 或 1/4 模型算完之后镜像显示即可。但要注意绕组激励的对称性三相绕组的三相对称电流在空间上不能简单用同一个对称面截断否则相位关系会出错。第二件事是控制网格总量。没必要全模型都用细网格只在铁芯窗口内角、铜片端部、气隙附近这些磁密梯度大的地方细化。外部空气域网格可以放得很粗因为那里磁场变化平缓。第三件事是设置合理的迭代收敛目标。Maxwell 默认的能量误差要求可能过高反而浪费时间。我一般设置为最大迭代 15 步能量误差低于 1% 就停止。如果连续两步之间目标量变化小于 0.5%也认为收敛。瞬态场还要关注时间步的独立性。用 1 ms 时间步长跑一遍再把时间步长缩到 0.5 ms 跑一遍如果损耗和磁密峰值变化很小说明时间步长足够。如果差异明显继续减小步长直到结果稳定。5. 常见问题与排查技巧实录5.1 最大磁密一下窜到 3 T网格和尖角的锅我第一次仿真时遇到一个特别离谱的结果铁芯最大磁密显示 3.2 T超过了硅钢片物理极限。一开始以为是模型建错了后面一个个排查发现是模型里把铁芯边缘做成了直角又没有在拐角处加密网格导致场奇异点被过分放大。解决办法有三个。第一几何建模时所有尖角都做倒圆角处理哪怕半径只有 0.5 mm也能明显改善磁场奇异。第二在可疑位置附近做局部网格加密对比加密前后最大磁密的变化。第三检查材料 B-H 曲线是否给到了足够高的磁场强度区间如果曲线只给到 H1000 A/m 就断了求解器在插值时会出偏差。经验上网格加密后最大磁密下降超过 10%就要继续加密下降小于 3%基本可以认为收敛了。5.2 损耗数字忽大忽小材料曲线和叠片设置不能省铜损和铁耗的准确性很大程度取决于材料设置。铁耗偏小很可能是因为 B-P 损耗曲线没设全默认材料只计算涡流损耗不计算磁滞损耗。铁耗偏大可能是叠片系数没设置把整块实体铁芯当成了纯硅钢导致涡流回路面积增大。铜损出现异常多半是绕组电流方向或并联关系设置不对。检查方法是看每个 Coil 的电流方向和电流大小是否与产品一致。铜片绕组如果几层并联每层电流分配由求解器自动计算方向错了结果几乎必然不对。还有一种情况是损耗数值“稳定但明显不合理”比如空载损耗比材料手册高出一倍。这时候建议做一次对比仿真把铁芯材料换成线性磁导率看涡流损耗曲线是否线性变化。如果非线性行为不符合理论大概率是网格或材料曲线的问题。5.3 还想预测分贝数从力密度到声场要过三关电磁场仿真得到力密度并不等于得到噪声分贝数。真要预测声压级需要经历“电磁场→结构振动→声辐射”三关。第一关是载荷映射。把 Maxwell 计算出的力密度导出到结构有限元模型注意网格要匹配否则载荷映射会损失空间分辨率。第二关是谐响应分析。施加 100 Hz 及其倍频激励得到变压器表面的振动加速度分布。第三关是声学边界元或有限元计算。以表面振动为边界条件计算远场声压级。这个过程每一步都可能引入误差而且计算量不小。如果项目只是定位噪声源头、确定降噪方向可以只做前两步。我在实际项目里发现力密度云图的高亮区域和样机实测噪声的显著频段位置有很强的相关性用电磁力分布直接判断相对噪声大小在工程上是可行的。5.4 顺带说说工频磁场抗扰度试验和仿真有什么关系变压器成品在验证环节通常要做电磁兼容测试其中有一项是工频磁场抗扰度试验英文全称是 Power Frequency Magnetic Field immunity test行业里常简写为 PFMF对应的基础标准是 IEC 61000-4-8。这项测试会施加一个较强的 50 Hz 磁场到受试设备上验证设备是否误动作或性能下降。变压器产品本身就是一个强磁场源在做 PFMF 试验前可以用 Maxwell 仿真提前看整机外部空间的磁场泄漏范围。我在做这类评估时会把空气域延伸到足够大并在外壳外设置多个测点观察磁场衰减趋势。如果外壳附近磁密仍然很高就要考虑增加屏蔽结构或调整内部元器件布局而不是等到实验室测试时才发现超标。这类仿真不需要建得特别精细重点是磁场泄漏趋势。只要材料和边界条件设置正确趋势判断的准确度足够用于设计早期决策。结尾我做变压器仿真的这几年最深的体会是磁场“唱歌”这件事越是听起来玄越需要用数值把它变成能解释的曲线和云图。第一次在 Maxwell 里调出磁密动画看到磁通沿铁芯循环推进耳边那阵嗡嗡声好像一下子被“看见”了。如果你也遇到一台噪声特别大的变压器我的建议是先别急着加隔音棉、换减震垫花半天把磁场仿真跑一遍找出力密度最高的位置。很多时候问题不在结构本身而是铁芯局部饱和、绕组端部涡流集中、或者接缝处磁拉力脉动太大。改磁路、调整铜片排列、优化接缝结构比事后隔音更治本。对于想深入的朋友下一步可以尝试把 Maxwell 的力密度导到结构软件里做振动响应分析再算声压级这条路走通后你对“磁场唱歌”的理解就不再停留在比喻层面了。 基于Maxwell 3D的变压器铜片绕组损耗与电磁噪声仿真分析
很多第一次接触变压器的人都会问变压器运行时为什么总在“嗡嗡”响工程圈里经常用一句话形容这种现象——磁场在“唱歌”。这歌声一点都不浪漫本质是铁芯在交变磁场下发生磁致伸缩加上绕组和结构件受到周期电磁力最后通过空气辐射形成了可听噪声。如果你正被变压器仿真搞得头大又恰好被这种噪声困扰那这篇手记应该能帮你理清思路。我最近做了一台小型干式变压器的 Maxwell 3D 铜片绕组损耗分布仿真过程中把铁芯磁密、铜损热点、电磁力分布以及噪声源头全部过了一遍。这篇文章不打算写成软件说明书而是按我实际的建模思路、参数设置、结果判读顺序来写顺便把踩过的坑都交代清楚。内容面向做变压器、电感、电抗器研发的同行以及正在自学电磁场仿真的学生朋友。你会看到我为什么选 3D、为什么关注铜片损耗云图、怎么从力密度判断“唱歌”位置也会看到一张可以直接抄作业的参数对照表。1. 为什么变压器的磁密分析这么重要从“唱歌”说起1.1 磁致伸缩与噪声产生的物理机理所谓磁致伸缩是指铁磁材料在磁化过程中沿磁化方向发生微小形变。硅钢片的磁致伸缩率通常在 1×10⁻⁶ 到 3×10⁻⁵ 量级这个变形量虽然极微但当它作用在整台变压器的铁芯上而且以 50 Hz 频率反复施加时铁芯就会像喇叭纸盆一样振动起来。这里有个特别容易忽略的细节电磁力与磁密的平方成正比。50 Hz 的交变磁场每秒钟正负翻转 50 次磁密平方之后频率翻倍变成 100 Hz。所以变压器噪声的基频基本是 100 Hz而不是 50 Hz听感上就是那种比较沉闷、连续的“嗡”声。如果电网里还有谐波或者变压器出现了直流偏磁噪声里还会叠加高次谐波和偶次分量声音会变得更尖锐、更复杂。要回答“这台变压器为什么吵”光靠耳朵听不出来只能靠仿真把磁密分布、力密度分布定量化。这也是我在这类项目里坚持先跑电磁场仿真的原因如果没有数据支撑结构降噪就是瞎蒙。1.2 磁场仿真能真正量化哪些工程指标我理解的变压器磁场仿真最终交付的不只是一张漂亮的磁密云图而是几个能直接指导设计的数值指标。第一个是铁芯磁密。磁密越高铁芯体积越小但饱和风险越大。通过仿真可以定位窗口拐角、T 型接缝处的局部高磁密点判断设计是否留足裕量。第二个是损耗分布。铁芯损耗包含磁滞损耗和涡流损耗绕组损耗则主要是铜损。在铜片或铜箔绕组结构中涡流效应会导致电流密度不均匀损耗热点往往集中在端部和边角。这类信息只有三维仿真才能给准。第三个是电磁力分布。力密度高的区域就是振动和噪声的主要激励源后续做结构谐响应、声学分析时都要从这里取载荷。第四个是杂散磁通。漏磁会穿入外壳、夹件、螺栓等金属结构引起局部涡流和发热。仿真能提前发现哪里有“看不见的磁场漏出去”也能为工频磁场抗扰度等测试提供参考。这些指标只要一个没量化设计阶段就容易埋雷。等到样机做出来再测噪声、测温升返工成本往往非常高。1.3 为什么选择 Maxwell 3D 做变压器电磁场仿真低频电磁场仿真的工具不止一种但 Maxwell 3D 在变压器领域用的人确实最多资料也最全。网上搜“maxwell3d 铜片变压器损耗分布仿真图文教程”基本都能找到可参考的案例。我选它的原因主要有三个。第一它有成熟的涡流场和瞬态场求解器。涡流场适合线性频域问题损耗、电感、阻抗都能直接算瞬态场适合非线性、谐波、直流偏磁和磁场动画这类时间相关问题。两种求解器覆盖了变压器仿真的绝大多数场景。第二它和结构、声学仿真工具之间有成熟的耦合接口。电磁力密度可以导到结构软件里做振动响应再把振动结果送到声学模块计算声压级形成完整的“电磁→结构→声学”链条。第三它支持从 2D 到 3D 的渐进式验证。先用 2D 快速看趋势再用 3D 精确定位端部效应工作流很顺。不过 3D 仿真计算量大、建模细节多这也是有些资料只给雏形不给完整参数的原因。下面我把我实际用的建模过程拆开来说重点讲那些不写进官方教程的关键点。2. 仿真模型搭建的完整思路2.1 几何建模让铜片绕组真正“片片可见”我这次仿真对象是一台小型干式变压器铁芯是 EI 型结构绕组采用铜片铜箔绕制。铜片绕组不像普通圆导线那样能简化为一个圆柱体因为铜片在厚度方向上的涡流分布不均匀必须把每一层铜箔建出来。铜在 50 Hz 工频下的趋肤深度大概可以按 δ sqrt(2ρ / (2πfμ₀)) 估算铜的电阻率约 1.7×10⁻⁸ Ω·m算下来大概是 9.3 mm。如果铜片厚度在 0.3~0.5 mm远小于趋肤深度理论上可以不必再细分厚度方向。但实际中铜箔边角处磁场畸变严重涡流回路会在狭窄区域内集中所以我仍然在铜箔的端部和转角处做了更细的网格剖分。几何建模时还有一个容易出问题的点层间绝缘层。铜箔层之间有绝缘材料如果只建铜箔不建绝缘层会让层间接触关系失真。尤其是并联铜箔之间如果有电位差绝缘层对电流分配的影响是不可忽略的。我的建议是在三维软件里把铁芯、铜箔层、绝缘层、骨架、空气域都建出来再导入 Maxwell。空气域要包住变压器整体尺寸至少是模型外轮廓的 1.3 倍否则漏磁会被截断损耗和电感都会失真。2.2 材料与 B-H 曲线硅钢片不是一根直线材料设置是仿真精度最大的影响因素。很多人随便选一个内置钢材就开跑结果磁密和损耗完全对不上。硅钢片材料必须输入实测的直流 B-H 曲线而且分方向。硅钢片在轧制方向的磁导率明显高于垂直方向如果单用一条各向同性曲线铁芯窗口和接缝处的磁密分布会偏。Maxwell 里可以用各向异性材料定义幅度比值按厂家数据给定通常轧制方向和非轧制方向的磁导率差异在 2~10 倍之间。叠片系数也要处理。硅钢片表面有绝缘涂层叠压后有效导磁截面积小于几何截面积。我习惯在材料里直接设叠片系数 0.96并选择叠片方向垂直于铁芯截面。如果不设叠片系数铁芯磁密会被低估损耗会被高估结果看起来“更保守”但实际上是错的。损耗曲线也要给全。铁芯损耗需要输入不同频率下 B-P 曲线Maxwell 才能用插值算出磁滞损耗和涡流损耗的分离值。只给 B-H 曲线不给 B-P 曲线铁耗就是拍脑袋。2.3 激励与求解器电流源还是电压源变压器仿真最常见的激励方式有两种电流源和电压源。我多数情况用电流源因为电流源的收敛性远好于电压源对线圈内部阻抗不敏感设置也简单。如果是涡流场求解器激励频率设为 50 Hz绕组的电流大小按额定工况给。注意线圈必须定义成 Winding并在 Coil 里设置导体匝数、极性以及串联/并联关系。铜片绕组层与层之间如果是并联的每层都要独立建模并指定到同一个 Winding 里否则电流分配失真损耗分布自然也不对。如果是瞬态场求解器时间步长一般取 0.5 ms 以内一个工频周期 20 ms 就有 40 个解。瞬态场建议至少跑两个周期第一个周期用于消除初始暂态第二个周期数据才稳定。当然瞬态场也支持电压源但需要附加外电路里的绕组电阻和漏感参数怎么给都得在电路里补全比电流源繁琐得多。如果只看稳定损耗和谐波涡流场就够如果关心直流偏磁、磁饱和波形、磁场动画那必须上瞬态场。这两种求解器的选择直接决定后面能不能看到你想看的那些量。3. 关键结果解读损耗分布、磁密与“磁场唱歌”位置3.1 铜片损耗分布如何快速定位发热区求解完成后第一步是查看铜片绕组的欧姆损耗体密度也就是 Ohmic Loss。把云图显示切到铜箔层实体上正常情况下可以看到损耗密度在铜片端部和转角位置高于中部这是因为端部磁场的切向分量更大涡流回路集中在局部。我这次仿真里铜片外侧拐角的损耗密度比中间区域高出大约 60%这个结论直接影响了后续设计铜片端部倒角并把绝缘层的厚度稍微加大损耗热点就明显缓和。还有一个经验在调网格之前不要急着分析损耗值。先用默认网格跑一遍记下最大体积损耗密度的位置然后局部加密该区域网格再跑一遍如果最大值变化超过 10%说明前一轮网格不够粗。有些产品需要进行倍频或谐波仿真比如在变频器输出侧使用的变压器。此时高频分量的趋肤深度会小很多铜箔厚度方向的网格必须加密到高频趋肤深度的 1/3 以下否则高频损耗会被严重低估。对 1 kHz 以上的谐波这个要求很苛刻建议先算一次谐波分量再叠加各阶损耗。3.2 铁芯磁密云图饱和不是“全饱和”而是“局部饱和”铁芯磁密看两个东西一个是云图整体分布一个是最大磁密数值。硅钢片饱和点一般看 1.7~1.8 T正常设计工作磁密取 1.3~1.5 T留一定的饱和裕量。我这次仿真额定工况下铁芯最大磁密约 1.46 T出现在窗口内角处。这个位置因为磁路截面积突然变化磁力线拥挤饱和风险最高。外角区域只有 1.1 T说明整体磁路还有余量但窗口内角是明确的薄弱点。如果看到局部磁密超过 1.8 T不用急着改整个铁芯尺寸先看那个高磁密区域是不是几何尖角或者网格奇异点造成的。如果是真实的磁路过窄那就调整窗口尺寸或者提高硅钢片牌号。判断依据很简单把该区域网格加密两轮如果最大磁密变化不大说明是真实饱和如果明显下降说明是网格伪影。瞬态场比涡流场更容易看出饱和问题因为饱和是随瞬时电流变化的。用电流源设置额定励磁电流一个周期内不同相位下磁密峰值位置会移动这种情况在涡流场里看不到必须用瞬态场的动画功能来看。3.3 电磁力分布噪声源到底在哪要说清噪声最终要落到力上。Maxwell 3D 可以输出 Maxwell 应力张量也就是磁曳力密度还能同时算通电导体在磁场中的洛伦兹力密度。两者叠加就是作用在铁芯和绕组上的总电磁力密度分布。“唱歌”的位置多半就是力密度最高的位置。我这次仿真里铁芯接缝区域和铜片绕组端部的力密度明显比其它位置高且呈周期性脉动频率正好是 100 Hz。接缝处磁拉力方向交替变化让铁芯像在被反复拉伸和挤压绕组端部则是轴向力集中会推动绕组整体振动。要看力密度可以在场结果里选择 Force Density直接加载到模型表面或体上观察其大小和方向。云图上颜色越深的地方就是设计需要重点关注的地方。如果要继续做声学预测把力密度数据导到结构有限元软件里做谐响应分析再把振动速度作为边界条件给声学模块就能计算声压级。这条链路虽然长但每一步都非常明确。3.4 让磁场“动”起来瞬态时序后处理磁场做动画是非常直观的后处理适合内部评审也适合给客户演示。用瞬态场求解器跑完一个周期后右键点击场图选择 Animate设定帧数和时间范围就能导出一段磁密矢量随相位变化的动画。动画里能看到磁通沿铁芯循环推进的过程。和网上常见的“同步电机磁场变化动画视频”不同变压器磁路是脉动磁通而不是旋转磁通因此磁场方向反复正反变化。这正好能说明为什么变压器铁芯一般做成闭合回路而同步电机需要气隙旋转磁场两者磁场形态完全不同仿真表现也完全不同。做动画时建议把时间步长控制在 0.5 ms 以内帧数取 40 帧以上否则磁通推进过程看起来卡顿。如果只关心某一相位的瞬时分布也可以不用动画直接选择对应时间点导出云图。4. 实操过程与核心环节实现4.1 从零跑通 Maxwell 3D 变压器损耗仿真的 7 个步骤给第一次用 Maxwell 3D 做变压器仿真的朋友我把操作流程拆成了 7 个步骤按这个顺序走基本不会乱。打开 Maxwell 3D选择求解器。若分析 50 Hz 稳态损耗选 Eddy Current若分析饱和、谐波或瞬态特性选 Transient。导入几何模型。从三维软件导出的格式建议用 STEP 或 SAT注意单位要统一成毫米。模型里包含铁芯、铜片层、绝缘层和空气域。赋材料。铁芯用各向异性硅钢片设置叠片系数铜片用纯铜绝缘层用塑料类材料。材料库缺失的手动输入 B-H 曲线和 B-P 曲线。创建绕组和线圈。每个绕组定义 Winding在 Coil 中设置匝数、电流方向和并联关系。设置激励。给每个绕组加电流或电压激励频率设为 50 Hz。涡流场加正弦稳态电流瞬态场定义随时间变化的电流波形。设置边界和网格。铁芯外表面用默认自然边界外侧空气域设置气球边界铁芯和铜片局部加密网格。求解。先跑一遍默认网格记录最大磁密和最大损耗密度然后局部加密关键区域再跑一遍对比两次结果差异。第 7 步的“加密-对比”是关键千万别省。很多新手跑一遍就出报告结果磁密最大值、铜损最大值全被网格误差带偏后续设计判断全部没有意义。4.2 关键参数设置对照表下面这张表是我这次仿真实际使用的关键参数可以复制过去再按自己产品调整。参数项推荐设置说明求解器Eddy Current / Transient稳态损耗用涡流场瞬态分析用瞬态场频率50 Hz工频变压器实际按产品需求调整铁芯材料各向异性硅钢片 叠片系数 0.96轧制与非轧制方向磁导率分开铜箔材料铜电阻率 1.7×10⁻⁸ Ω·m按工作温度修正电阻率空气域模型外扩 30% 以上避免漏磁被边界截断网格铁芯最大边长 2~4 mm窗口内角和接缝局部加密网格铜箔端部最小边长 0.3~1 mm保证涡流分布分辨率激励电流源额定电流收敛性好推荐首选时间步长瞬态≤0.5 ms一个周期 40 步以上求解迭代最大 15 步误差 1%收敛后检查能量误差这些值不是恒定不变的。比如铁芯更大网格可以放宽铜箔更薄网格可能要更细。核心原则是在关键场变化剧烈的区域加密在均匀区域放松这样才能把计算资源用在刀刃上。4.3 从求解时间到收敛控制3D 仿真别硬跑3D 仿真最现实的障碍是计算时间。我这次模型大约三百万网格单元单台 16 核工作站跑了接近两个小时。如果不做任何对称简化一个复杂三相变压器的 3D 瞬态场算上几天也不奇怪。所以第一件事是找对称面。铁芯和绕组结构如果存在对称性可以切 1/2 或 1/4 模型算完之后镜像显示即可。但要注意绕组激励的对称性三相绕组的三相对称电流在空间上不能简单用同一个对称面截断否则相位关系会出错。第二件事是控制网格总量。没必要全模型都用细网格只在铁芯窗口内角、铜片端部、气隙附近这些磁密梯度大的地方细化。外部空气域网格可以放得很粗因为那里磁场变化平缓。第三件事是设置合理的迭代收敛目标。Maxwell 默认的能量误差要求可能过高反而浪费时间。我一般设置为最大迭代 15 步能量误差低于 1% 就停止。如果连续两步之间目标量变化小于 0.5%也认为收敛。瞬态场还要关注时间步的独立性。用 1 ms 时间步长跑一遍再把时间步长缩到 0.5 ms 跑一遍如果损耗和磁密峰值变化很小说明时间步长足够。如果差异明显继续减小步长直到结果稳定。5. 常见问题与排查技巧实录5.1 最大磁密一下窜到 3 T网格和尖角的锅我第一次仿真时遇到一个特别离谱的结果铁芯最大磁密显示 3.2 T超过了硅钢片物理极限。一开始以为是模型建错了后面一个个排查发现是模型里把铁芯边缘做成了直角又没有在拐角处加密网格导致场奇异点被过分放大。解决办法有三个。第一几何建模时所有尖角都做倒圆角处理哪怕半径只有 0.5 mm也能明显改善磁场奇异。第二在可疑位置附近做局部网格加密对比加密前后最大磁密的变化。第三检查材料 B-H 曲线是否给到了足够高的磁场强度区间如果曲线只给到 H1000 A/m 就断了求解器在插值时会出偏差。经验上网格加密后最大磁密下降超过 10%就要继续加密下降小于 3%基本可以认为收敛了。5.2 损耗数字忽大忽小材料曲线和叠片设置不能省铜损和铁耗的准确性很大程度取决于材料设置。铁耗偏小很可能是因为 B-P 损耗曲线没设全默认材料只计算涡流损耗不计算磁滞损耗。铁耗偏大可能是叠片系数没设置把整块实体铁芯当成了纯硅钢导致涡流回路面积增大。铜损出现异常多半是绕组电流方向或并联关系设置不对。检查方法是看每个 Coil 的电流方向和电流大小是否与产品一致。铜片绕组如果几层并联每层电流分配由求解器自动计算方向错了结果几乎必然不对。还有一种情况是损耗数值“稳定但明显不合理”比如空载损耗比材料手册高出一倍。这时候建议做一次对比仿真把铁芯材料换成线性磁导率看涡流损耗曲线是否线性变化。如果非线性行为不符合理论大概率是网格或材料曲线的问题。5.3 还想预测分贝数从力密度到声场要过三关电磁场仿真得到力密度并不等于得到噪声分贝数。真要预测声压级需要经历“电磁场→结构振动→声辐射”三关。第一关是载荷映射。把 Maxwell 计算出的力密度导出到结构有限元模型注意网格要匹配否则载荷映射会损失空间分辨率。第二关是谐响应分析。施加 100 Hz 及其倍频激励得到变压器表面的振动加速度分布。第三关是声学边界元或有限元计算。以表面振动为边界条件计算远场声压级。这个过程每一步都可能引入误差而且计算量不小。如果项目只是定位噪声源头、确定降噪方向可以只做前两步。我在实际项目里发现力密度云图的高亮区域和样机实测噪声的显著频段位置有很强的相关性用电磁力分布直接判断相对噪声大小在工程上是可行的。5.4 顺带说说工频磁场抗扰度试验和仿真有什么关系变压器成品在验证环节通常要做电磁兼容测试其中有一项是工频磁场抗扰度试验英文全称是 Power Frequency Magnetic Field immunity test行业里常简写为 PFMF对应的基础标准是 IEC 61000-4-8。这项测试会施加一个较强的 50 Hz 磁场到受试设备上验证设备是否误动作或性能下降。变压器产品本身就是一个强磁场源在做 PFMF 试验前可以用 Maxwell 仿真提前看整机外部空间的磁场泄漏范围。我在做这类评估时会把空气域延伸到足够大并在外壳外设置多个测点观察磁场衰减趋势。如果外壳附近磁密仍然很高就要考虑增加屏蔽结构或调整内部元器件布局而不是等到实验室测试时才发现超标。这类仿真不需要建得特别精细重点是磁场泄漏趋势。只要材料和边界条件设置正确趋势判断的准确度足够用于设计早期决策。结尾我做变压器仿真的这几年最深的体会是磁场“唱歌”这件事越是听起来玄越需要用数值把它变成能解释的曲线和云图。第一次在 Maxwell 里调出磁密动画看到磁通沿铁芯循环推进耳边那阵嗡嗡声好像一下子被“看见”了。如果你也遇到一台噪声特别大的变压器我的建议是先别急着加隔音棉、换减震垫花半天把磁场仿真跑一遍找出力密度最高的位置。很多时候问题不在结构本身而是铁芯局部饱和、绕组端部涡流集中、或者接缝处磁拉力脉动太大。改磁路、调整铜片排列、优化接缝结构比事后隔音更治本。对于想深入的朋友下一步可以尝试把 Maxwell 的力密度导到结构软件里做振动响应分析再算声压级这条路走通后你对“磁场唱歌”的理解就不再停留在比喻层面了。