COMSOL变压器电磁场仿真：磁密分布与电路状态联合求解实战

在变压器电磁场设计里，磁密分布和电路运行状态是互为因果的两件事：磁路几何和材料决定了磁链的走法，而磁链又反过来决定绕组的感应电势、电流和阻抗。以前我主要靠路模型快速估算，设计常规油变基本够用；可一旦碰上非标准异形铁心、多绕组或强漏磁场结构，路模型算出的磁密只能看个大概，连铁心拐角会不会局部饱和都判断不了。后来我把COMSOL作为详细验证工具，专门用来做变压器电磁场模型，把磁密云图和电路外特性放在同一个框架里求解。这篇就结合我做过的项目，把COMSOL仿真变压器电磁场的建模思路、磁密分布分析方法和电路状态仿真细节梳理一遍，适合正在入门变压器有限元仿真、或者已经会建简单模型但想看深一层的朋友。

1. 变压器电磁场仿真为什么首选COMSOL：从项目需求反推工具边界

1.1 先想清楚要算到什么程度，再选工具

变压器电磁场仿真不能一上来就开软件。我见过不少人，几何模型拉得很漂亮，最后只输出一张五颜六色的磁密云图，然后就没有然后了——这样的仿真其实没有回答任何工程问题。动手之前，应该先列一个问题清单：是要评估铁心局部饱和？要算短路阻抗？要看负载下的漏磁场分布？还是要给温升、噪声、机械强度计算提供磁密和损耗数据？清单不同，工具的用法完全不同。

如果你只需要某个变压器截面的平均磁密、绕组电感值，路模型或者磁路的解析公式完全够用，用大型有限元软件反而低效。但工程上经常会遇到第二类需求：局部磁密超过1.8T需要判断是否调整铁心截面积；短路时漏磁场不均匀导致绕组电动力集中在端部；或者电容器放电这类瞬态工况下磁路内部的磁密随时间怎么变。这些都要通过场的求解才能得到。COMSOL在这种需求里很合适，因为它把几何建模、材料参数、边界条件、求解控制和后处理放在同一个界面下，模型维护相对省心，不会在几何文件、网格文件和求解器设置之间来回倒腾。

1.2 磁矢势方程与COMSOL“磁场”接口的对应关系

COMSOL的“磁场”（Magnetic Fields）接口，底层求解的是磁矢势A，磁密B = ∇ × A。在频域求解下加正弦激励，方程可以写成：

∇ × (ν ∇ × A) + σ jωA = J

其中ν是磁阻率（磁导率的倒数），σ是电导率，ω是角频率，J是源电流密度。变压器铁心里的硅钢片，实际是叠片结构，每层之间绝缘，涡流路径被切断。三维模型如果一片一片地建叠片，网格量会非常恐怖，也没必要。COMSOL里有“叠片铁心”这种等效材料模型来处理，直接把铁心按照各向异性等效磁导率处理：叠片方向磁导率降低，面内磁导率按硅钢片材料走。这个细节很多入门资料没讲透，实际建模时直接影响磁密分布，特别是叠片方向如果有漏磁路径，等效磁导率设得过高会让结果偏乐观，局部漏磁会被低估。

1.3 COMSOL在磁路-电路联合仿真上的天然优势

变压器电磁场模型很少只看纯磁场，它必然连着外部电路：电压源、负载阻抗、短路阻抗、多绕组串并联。COMSOL的“磁场”接口可以和“电路”（Electrical Circuit）接口强耦合求解，绕组域使用线圈特征，指定匝数、导线截面积，再选择由外部电路驱动，就可以把空间磁场和电路方程同时解。这种耦合的价值在于，你不用手工换算出特定工况的电流再塞进去。比如短路仿真，你真正要算的是：给定一次侧电压，二次侧短路，系统稳定后一次侧电流和二次侧电磁力。这种问题在COMSOL里直接把二次侧用外部短路线圈连起来即可，一次侧接交流电压源，软件同时解磁场和电路，电流会自然收敛到一个合理的短路值，而不是靠人工迭代。

2. 建模前必须抠清楚的几个问题：几何简化、材料B-H曲线与物理场组合

2.1 几何简化，以“磁路闭合路径和漏磁路径”为双重基准

变压器几何看着简单，实际处理好它的简化是仿真成败的关键一步。铁心（心柱、铁轭、拐角）、绕组（高压、低压、可能还有调压绕组）、夹件、油箱、绝缘结构、套管间距……全建出来模型当然完整，但你根本算不动，也没必要。我的简化原则有三条：

第一，铁心主体必须完整。几何上的倒角、螺栓孔如果不影响磁通走向可以去掉，但拐角处的形状尽量保留，因为那里往往是磁密最高、最容易饱和的位置。铁心窗口的高度和宽度、心柱之间的距离直接决定磁路长度和漏磁面积，这些尺寸必须按图纸准确录入。

第二，绕组用环形域代替，不逐匝建模，匝数在“线圈”特征里赋值。这样做的原因是变压器绕组匝数动辄几百上千，逐匝建模的网格量在三维模型里是完全不可行的。绕组区域内部的磁密细节虽然损失了，但外部磁密分布和整个磁路的磁链计算完全够用。如果研究绕组内部的涡流损耗分布，那另说，需要更精细的几何。

第三，漏磁路径相关的部件必须保留，比如铜屏蔽、夹件、油箱壁。这些导磁导体在漏磁作用下会产生涡流和附加损耗，把它们删了，短路工况下的磁密分布和阻抗计算结果会偏得厉害。建模前先问自己一句：哪些金属部件的存在会改变磁通的路径？答案列出来，这些部件基本上都要进模型。

2.2 材料参数：B-H曲线是磁密分布仿真的“命根子”

线性材料模型只能用于很初步的估算。变压器铁心正常工作在接近饱和区，必须用非线性B-H曲线。问题经常出在曲线来源：硅钢片厂家一般给的是50Hz正弦磁化数据，这是交流磁化曲线；如果你做瞬态分析，原则上应该用直流磁化曲线（基本磁化曲线），两者在接近饱和段差别不大，但在低场区差别明显，直接影响励磁电流的基波分量。

COMSOL里用插值函数把厂家数据逐点输入，注意x轴是磁场强度H（单位A/m），y轴是磁通密度B（单位T）。厂家给的数据单位五花八门，可能是A/cm或者Gs，必须先统一量纲再录入。另一个坑是曲线尾部：厂家数据通常只测到H=10000A/m左右，但求解过程中局部单元可能碰到更高的场强，如果你不把曲线尾部适当外推一下，求解器就会报“插值函数外推错误”，或者强行截断导致迭代震荡。我习惯在最后一个数据点之后延伸一条平缓上升的直线，斜率为真空磁导率μ₀，这样既不会出错，也不会人为制造一个虚假的“硬饱和”。

2.3 物理场接口与边界条件的组合

变压器电磁场COMSOL建模最常用的组合是：三维或二维轴对称的“磁场”（mf）接口搭配“电路”（cir）接口。如果变压器结构是旋转对称的，比如单相壳式变压器，可以降维成二维轴对称模型，计算速度快很多，但绕组端部和三相结构就得老老实实用三维。计算域包括铁心域、绕组域和包围它们的空气域或变压器油域。

边界条件上，外边界用磁绝缘（n × A = 0）即可，但空气域要取得足够大。对于需要考察外部漏磁分布的情况，可以在最外层加“无限元域”来模拟开阔空间，避免因空气边界太近而人为压缩漏磁路径。实际做普通变压器的电磁性能验证时，外边界设置在变压器本体尺寸的2倍左右基本够了。如果发现磁力线在空气域边界附近出现不自然的“回流”，说明空气域不够大或者边界条件取错了位置，这时候不要急着调模型参数，先扩大空气域试一次。

3. 磁密分布分析：从云图“看颜色”到定量提取数据

3.1 先看云图，再看箭头图：铁心磁密的判读习惯

仿真跑完，第一件事是看磁密云图。正常变压器铁心在额定励磁下，心柱和铁轭大部分区域的磁密应该在1.6～1.75T之间（取决于硅钢片牌号），云图应该比较均匀。我最关注的是拐角处和T形接缝附近，这些位置因为磁路截面突变，磁通会拥挤，磁密局部可能冲上1.9T甚至更高。如果云图上出现成片的深红色区，说明铁心已进入深度饱和，需要回去改截面或改材料。

除了云图，再叠加一个箭头图，能直观看到磁通走向，特别是有没有不该出现的旁路磁通。例如三相变压器，中心柱磁通和两边柱会在上下铁轭汇流，磁通叠加，铁轭截面往往比心柱大，这种设计细节通过云图一眼就能看出来对不对。箭头图还能帮你发现建模时的低级错误，比如某根铁心柱里磁力线方向与其他柱相反，要么是绕组方向接反，要么是线圈电流设置正负号搞错了。

3.2 把“图”变成“数”：截面、探针和派生值的正确用法

云图只能定性看，工程判断必须定量。我常用的三个手段：

想看某个铁心柱截面上的磁密分布，在该位置定义一个切割面，把面上的磁密值导出，最大值、平均值一目了然。这里要注意切割面的位置要避开拐角圆弧段和T形接缝，否则数值会被局部集中磁密污染，不能代表心柱的真实状态。

想看某一点磁密随电流或频率的变化，用“探针”在几何点或线圈上设置监测量，参数扫描时可以实时看曲线。我经常用探针监测铁心最热点和绕组端部漏磁密。做空载参数扫描时，探针曲线能直接显示饱和趋势，不用跑完一个参数点就去看一次云图。

计算磁链、电感，COMSOL的线圈节点内置了磁链计算，后处理里直接读结果即可，不需要自己在积分里折腾。需要特别提醒的是，电感计算依赖模型里绕组激励设置是否正确，匝数和导线截面积都填对了，算出来的自感和互感才有意义。

3.3 云图异常的判断与模型修正

云图常见异常有几种：铁心内部磁密出现条纹状或者锯齿状分布，十有八九是网格太粗，铁心材料在网格单元边界处产生数值跳跃；磁力线在空气中出现不该有的密集环路，通常意味着计算域外边界取得太近；铁轭磁密明明应该比心柱低，云图上反而高，多半是几何上铁轭和心柱的截面比例没设置对；绕组区域外部出现不正常的强磁场，优先检查绕组电流方向是否一致。遇到这些现象，不要急着调材料参数，先看模型设定，这类问题的根因大多数在网格和几何，不在物理学。把模型文件和云图截图放一起回看，排查效率会高很多。

4. 电路状态仿真的实现路径：绕组激励方式与短路/空载工况模拟

4.1 绕组激励方式的三种选法

在COMSOL的线圈特征里，可以选电流驱动、电压驱动，再或者外部电路。三者的适用场景不一样，选错了轻则算得慢，重则结果完全不对。

做变压器电磁场模型，我几乎不再用纯电流源去算满载情况，因为变压器的性能评估本质上关心的是外电路特性，而不是内禀的磁场。尤其是短路阻抗、电压调整率这类指标，必须通过电压驱动或者外部电路才能自然反映。但有一个例外：在模型调试早期，为了快速验证几何和网格有没有问题，电流源是最稳的选择，因为它绕开了电路方程的非线性。

4.2 空载工况：励磁电流和磁密之间的微妙关系

空载工况下的变压器，一次侧接额定电压，二次侧开路。此时一次侧电流就是励磁电流，它完全由铁心磁路的非线性特性决定：铁心越接近饱和，励磁电流畸变越大。在COMSOL里做频域稳态分析时，你用正弦电压驱动，软件会解出对应电流的数值解。要注意的是，频域解只有一个复数幅值，它实际上是对非线性电流波形做了基波等效。如果你关心电流波形的尖顶程度、谐波含量，必须用瞬态求解器，并且要在一次侧外电路里串联一个很小的电阻，以抑制数值振荡。

空载工况的磁密分布也有其特征：铁心整体磁密偏高，但绕组附近漏磁很少。如果空载励磁电流的计算值和实测值差了超过10%，优先怀疑B-H曲线不是该批次的材料数据。我遇到过几次，把厂家典型曲线换成实测批次曲线后，励磁电流偏差立刻降到3%以内。所以空载工况是验证材料模型最好的试金石。

4.3 短路工况的实现与短路阻抗提取

短路是变压器最严酷的工况之一，算它主要是为两个数据：短路电流和短路阻抗。模型搭建时，把二次侧绕组线圈连接到外部电路的短路支路——直接闭合，或者通过一个很小的电阻模拟绕组直流电阻——一次侧给额定电压。瞬态仿真刚开始会看到一个比较大的冲击电流，然后慢慢稳定，稳定后的电流就是短路电流稳态值。

短路阻抗Zk的常用计算方法是：一次侧施加额定电压，二次侧短路，根据一次侧稳态短路电流Ik计算，Zk = U1 / Ik。在COMSOL里可以用全局表达式直接输出这个比值。漏磁场在这一工况下也是关键输出，因为端部漏磁密越大，绕组电磁力越大。我在处理短路工况时，习惯在后处理中单独画一个绕组端部的局部磁密云图，数值单位换算到mT级别去看，很多机械强度分析需要的数据就来自这里。漏磁场在绕组端部的分布通常呈现“两端强、中间弱”的特征，如果云图显示的规律正好反过来，那就要检查绕组几何的端部形状有没有建错。

4.4 负载工况与电压调整率验证

带上额定负载后，二次侧电流在漏抗上产生压降，输出电压会比空载低。在外部电路里给二次侧接上负载阻抗，然后扫描负载从0到满载，看输出电压和效率曲线，可以筛选出变压器是否满足额定电压调整率承诺。负载工况的磁密分布和空载有很大不同：负载电流产生的磁动势会在绕组周围建立明显漏磁场，铁心磁密反而可能因为负载电流的去磁或助磁作用而有所变化。三相变压器还要注意接线组别对磁路的影响，比如Yd11和Yyn0的铁心磁通分布并不完全相同，仿真时二次绕组的连接方式必须和实际一致，否则负载工况下的磁密分布对不上实测。

5. 求解设置与网格质量控制：那些让仿真“发散”的隐形杀手

5.1 非线性不收敛，第一步查几何和网格，而不是材料

COMSOL电磁仿真发散，尤其是非线性B-H曲线工况下的迭代不收敛，十个里有八个不是物理问题，是几何和网格问题。常见元凶有四个：

第一，铁心和绕组表面存在细小的圆角或尖角，网格在这些位置生成质量很差的狭长单元，导致雅可比矩阵病态。解决方案是删掉直径小于2mm的倒角特征，这些细微几何对宏观电磁场的影响微乎其微。

第二，两个域之间有接触面但几何没有形成完整的共享边界，这会生成重合但不连续的面网格，磁场法向分量的连续性被破坏，求解器反复迭代也无法收敛。排查方法是检查“接触对”设置，或者用“粘合”操作把接触面合并。

第三，空气域和外边界产生退化面，比如非常扁平的拉伸几何，网格密度不够时单元体积为负。这种情况下只需把退化几何重新分割成多个规则形状。

第四，材料曲线在数据尾部斜率突变，导致迭代时试算点落在畸变区域。处理方式是在曲线尾部做光滑外推，保证一阶导数连续。

排查时按这个顺序做：检查网格质量统计（平均质量低于0.4基本是网格坏了）→ 删除细小倒角和尖角特征 → 把B-H曲线尾部数据点去掉几个再插值。大多数情况下，这么处理完，求解器自己就稳定了。

5.2 网格无关性验证：不能跳过的底线检查

很多初学者跳过网格无关性验证，直接看云图下结论。问题是，你的网格稍微加密一点，磁密峰值就变了5%，这个结果敢拿去评审吗？我的做法是固定求解器参数，只细化网格，依次取三档：粗网格、参考网格、细网格，分别提取三个指标——心柱中心区磁密峰值、绕组电感、绕组交流电阻或损耗。如果后两档结果相差在1%以内，说明达到网格无关了；如果还差2%以上，再加密继续验证。

注意，验证指标不能只盯磁密峰值。磁密峰值是局部量，对网格变化比较敏感，而电感和损耗是积分量，对网格变化的敏感度相对低。如果一个模型磁密峰值差3%但电感差只有0.5%，说明网格误差主要集中在局部高磁梯度区域——这个区域恰好是铁心饱和风险最大的地方，不能掉以轻心。网格加密的时间成本会翻倍，但这比后期评审时被质疑数据可靠性再返工划算得多。

5.3 求解器选择的实战经验：频域线性化还是瞬态

变压器电磁场仿真，很多人卡在耗时上。如果可以接受频域稳态近似，优先用频域求解器，速度比瞬态快两个量级。但频域里处理非线性B-H曲线用的是等效磁导率，它把磁滞和非线性做了基波等效。如果关心瞬态电流波形、涌流现象，或者电源里含有谐波，就必须用瞬态求解器。

瞬态求解器要注意时间步长的选择。不要完全依赖自适应步长，我习惯手动设置最大时间步为周期T的1/64或1/128，否则会在磁通过零附近产生蛙跳噪声，体现在电流曲线上就是叠加了高频毛刺。如果模型在瞬态里特别慢，可以先用频域解算出一个稳态初值，再切到瞬态接着跑，能省下大量过渡时间。这个方法在强非线性模型里尤其管用，因为频域解已经接近稳态工作点，瞬态求解器只需要做小幅修正。

6. 仿真结果对标实测：误差根源与一套可复用的修正工作流

6.1 误差最大来源：材料参数与二维简化

把COMSOL仿真结果和变压器出厂试验数据对比，多数偏差来自两个地方。下表是我在实际项目里遇到的主要误差源和量级估算：

材料参数是最大的误差来源。厂家的典型B-H曲线和实际批次硅钢片存在差异，特别是高场段差异可能导致磁密峰值判断偏差10%以上。如果时间允许，最好直接向硅钢片厂家要该批次的实测B-H数据；没有的话，用小规格的环形铁心试样自己测一条磁化曲线，也比直接用典型曲线可靠。第二是二维模型对端部漏磁的忽略。二维模型计算速度快，但绕组端部是三维结构，端部漏磁场计算天然失真，如果研究负载工况的电压调整率，二维结果只能作为趋势参考，不能作为设计结论。三维模型虽然慢，但在漏磁、短路阻抗、端部电动力这些问题上，可靠性明显更高。

6.2 对标流程：先空载、再短路、最后负载

我习惯的验证顺序很固定：先做空载试验对标，再做短路阻抗对标，最后才做负载效率。这个顺序不是随便定的，因为每一步都是下一步的前提。

空载对标看励磁电流和空载损耗。如果励磁电流偏差大，大概率是材料B-H曲线不准；如果励磁电流对但空载损耗对不上，问题出在铁耗模型——频域有限元默认用的等效磁导率不包含磁滞损耗，需要在后处理里叠加 Steinmetz 公式或者磁滞模型才能补上空载损耗。这一步通常迭代两三轮，每一轮只改一个参数（比如换一条B-H曲线），改完重新跑，记录偏差变化趋势。

短路对标看短路阻抗。如果短路阻抗偏差大，更多是漏磁路径几何建模有误，比如夹件位置、油箱壁距离不对，或者绕组之间的主空道尺寸填错。短路阻抗对几何尺寸非常敏感，哪怕主空道差了0.5mm，结果也会体现在百分位上。

负载对标受前两项影响，前两项过了，负载结果一般不会差太远。负载对标时重点看额定负载点的输出电压和效率，如果你在外部电路里设置了正确的负载阻抗，这两个指标能直接和出厂试验报告上的数值逐项核对。

6.3 做一个会“成长”的模型修正档案

我做变压器电磁仿真这几年，最大的体会是：仿真模型的置信度不是一次建出来的，而是长期迭代标定出来的。我会为每个产品系列维护一个文件夹，里面存几何文件、材料曲线、网格设置、求解器参数和对应批次的实测数据。每做一次调整修正就更新材料参数表，慢慢形成一个“该系列专属的材料等效参数库”。下次新项目直接调用这套参数，省去不少试错。

另外，保存COMSOL模型时文件名里带上日期和工况，比如“SFZ11-10000-35_空载_20250112.mph”，避免不同工况的模型互相覆盖。版本管理这种看似朴素的动作，在方案评审或者后期回溯问题时能节省大量时间。我自己就遇到过用错版本模型出图交报告的尴尬，从那以后坚持这套命名规则，再没出过这类问题。

最后分享一个小细节。很多人仿真做完只截图云图贴报告，但我建议把磁密分布的原始数值导出成CSV，配合Python做一次二维热图或者剖面图，这样报告里既能保留全貌又能灵活截取重点区域。我整理结果时，通常会把仿真云图和实测数据放进同一张图里对比，让偏差一目了然。做变压器电磁场仿真，最大的门槛不是软件操作，而是对变压器物理本质的理解——什么时候该相信模型，什么时候该怀疑模型，全靠对磁路、电路和材料特性这三件事的把握。