很多第一次接触变压器的人都会问:变压器运行时为什么总在“嗡嗡”响?工程圈里经常用一句话形容这种现象——磁场在“唱歌”。这歌声一点都不浪漫,本质是铁芯在交变磁场下发生磁致伸缩,加上绕组和结构件受到周期电磁力,最后通过空气辐射形成了可听噪声。如果你正被变压器仿真搞得头大,又恰好被这种噪声困扰,那这篇手记应该能帮你理清思路。
我最近做了一台小型干式变压器的 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 曲线是否给到了足够高的磁场强度区间,如果曲线只给到 H=1000 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 的力密度导到结构软件里做振动响应分析,再算声压级,这条路走通后,你对“磁场唱歌”的理解就不再停留在比喻层面了。