做变压器仿真的同行应该都有这种感觉：稳态场算起来一切岁月静好，磁场分布均匀、磁密云图漂亮，可一旦切换到三相短路工况，模型就开始各种闹脾气。我最近刚把一个110kV油浸式变压器的三相短路电磁场仿真彻底跑通，电流、电压、磁密三个量全部对上了理论估算范围，中间踩坑无数，所以把这套方法和心得完整写出来，希望能帮到同样在搞变压器短路仿真的朋友。
这篇内容基于我用的Comsol Multiphysics，版本不重要，核心思路通用。主要聊三件事：三相短路瞬态过程的物理本质，怎么在Comsol里把一个能算的变压器模型搭起来，以及如何从时间波形和磁密云图里读出可信的数据。适合正在做变压器短路电磁力校核、继电保护定值验证或者想搞懂绕组电动力来源的工程师和研究生。如果你只是会建立一个静磁场模型，这篇文章能帮你跨过从稳态到瞬态那道坎。
1. 变压器三相短路到底在短路什么？先把物理过程捋清楚
1.1 为什么短路仿真比稳态仿真更考验建模功底
我在刚接触这个课题的时候，第一反应是“不就是把三个绕组都接到故障上嘛，稳态计算加上负载不就完了？”真做起来才知道完全不是一回事。稳态仿真只需要关心工作点附近的磁场分布，而三相短路是严重的瞬态电磁过程，必须考虑磁场能量的突变、绕组电流的直流分量衰减、铁心饱和程度的大幅变化，还有绕组之间的电磁耦合。这些因素在稳态模型里根本体现不出来。三相短路瞬间，由于磁链不能突变，某一相绕组的电流会突然叠加一个很大的非周期分量，使得总电流在头几个周波内可能达到稳态短路电流的1.5到2.5倍，这个冲击电流直接决定了绕组能够承受的电磁机械力。而绕组受到的电磁力又与磁密和电流的乘积成正比，所以这三者必须放在同一个瞬态场里求解，才能看到真实的时空演化过程。
用Comsol做这个仿真的核心优势是它支持磁场模块与电路模块的强耦合。你可以把变压器绕组的端部接成三相短路逻辑，也可以在外部电路里设置三相电压源以及短路开关，然后让Comsol在每个时间步上同时更新磁场分布和电路电流。这种耦合关系在只做磁场分析或者只做电路仿真的工具里是看不到的。真正的难点在于，几何模型、材料非线性、电路拓扑、网格剖分、时间步长，任何一环处理不好，瞬态结果就会发散，或者波形出现明显的非物理振荡。
1.2 三相短路时刻的选择：A相过零还是B相峰值？这对结果影响很大
很多初做短路仿真的朋友会忽略“短路时刻”这个概念。三相短路并不是随便在一个时间点把断路器合上就行。故障发生的瞬间，各相电压的相位角不同，导致电流的直流分量大小完全不同。以最常见的三相短路为例，如果A相电压在过零点附近发生短路，A相电流的非周期分量会非常大；如果恰好A相电压处于峰值附近，那么A相电流几乎不含非周期分量，直接过渡到对称短路电流。最恶劣的短路工况是让某一相电压过零时短路，这一相会产生最大的冲击电流。因此，仿真时必须明确指定短路初始相位角，而且最好设置至少两到三个不同的短路相位角来对比。
我在模型里用参数化扫描做了三个短路相位角：0度、30度、90度（以A相电压过零为基准）。结果非常直观，0度时A相电流峰值比90度时高出了大约1.8倍。这个信息对后续的机械力校核至关重要，最大电磁力对应的就是最恶劣的短路相位角。如果只做单工况，很容易低估绕组受力风险。所以建议你拿到项目的第一件事，就是先把短路时刻这个参数定义好，留成一个可以扫描的变量，而不是固化在某个具体时间。
2. 从零搭一个能算的变压器模型：几何、材料与边界条件
2.1 几何建模：2D轴对称还是3D全模型？
先说一个我自己的判断：如果你要分析的是三相变压器，就别偷懒用2D轴对称模型。三相变压器的铁心是三相三柱或三相五柱结构，每相绕组在空间上并不是轴对称的。2D模型只能模拟单相或完全对称的某些切面，无法正确反映三相绕组之间的互感耦合以及短路时三相电流的不平衡。我最终用的是3D模型，但做了一定简化：把铁心的叠片效应用等效各向异性电导率代替，绕组采用多匝线圈域，并用均匀电流密度等效铜导线，不建每匝导线的实体。
几何建模有几个关键尺寸需要谨慎：绕组高度、铁心窗宽度、高低压绕组之间的主空道距离、绕组端部到铁轭的距离。这些尺寸直接决定了漏磁场的路径和大小。漏磁场是产生绕组短路电动力和杂散损耗的主要来源，所以千万不能为了画图方便把绕组和铁心画成紧贴在一起。我在建模时特意保留了实际的结构间隙，虽然网格数量上去了，但结果可信度高很多。
另外，三相变压器的高压绕组通常是Y接或D接，低压绕组也有多种连接组别，几何位置上的相序排列决定了短路时各相电流的空间相位。我建议在模型里保持A、B、C三相绕组的空间分布与实际一致，不要为了对称性而任意调整位置。这样后处理提取每相电流时，才能正确对应到绕组的实际相别。
2.2 材料参数与非线性B-H曲线：没有它就不要谈磁密
变压器铁心正常工作时是工作在饱和点附近的，尤其是短路时铁心会进入深度饱和，所以材料参数绝不能给一个固定的相对磁导率。我见过太多人直接在材料节点里填了一个μr=2000，结果算出来的磁密分布跟实际严重不符，因为2000是线性假设，只要磁密稍微超过1.7T，铁心的磁导率会急剧下降到几十甚至个位数。正确做法是导入真实的铁心B-H曲线，而且曲线要覆盖到过饱和区域，比如至少到2.2T以上。短路时铁心某区域的磁密可以达到2.0T以上，如果你的曲线只到1.8T，插值外推的误差会非常大。
我用的办法是从公开的硅钢片手册里提取一段典型的B-H曲线，再在Comsol材料库里创建插值函数。注意B-H曲线的数据点要足够多，尤其是在饱和区域，斜率变化剧烈，数据点太少会导致非线性迭代不收敛。绕组材料设定也非常关键，铜绕组的电导率在常温下是58MS/m左右，但在瞬态过程中绕组会发热，电导率会变化。如果你只关心前几十个毫秒的短路过程，可以忽略温度变化，假设电导率恒定；如果要做热-结构耦合或者长时间短路，就必须加上温度依赖。
2.3 短路工况的电路耦合设置：外部电路还是内置线圈？
Comsol的磁场模块有两种绕组激励方式：一种是直接把电流密度加载到线圈域上，另一种是通过“线圈”特征连接到外部电路。对于变压器短路仿真，必须用第二种方式，因为短路电流不是预设好的，而是由磁场和电路共同决定的。你需要在模型外部建立一个电路拓扑，包含三相电压源、线路阻抗、断路器（用可变电阻模拟）、以及变压器的绕组。
我最开始用了一个很简单的电路模型：忽略电源内阻和线路阻抗，直接让三个绕组端部短接，电压源加在另一个绕组侧。结果算出来的短路电流是纯感性回路，相位和衰减趋势都还可以，但因为没有考虑系统阻抗，电流幅值偏大。后来我把系统阻抗加了进去，用短路容量反算阻抗值，波形和理论计算就基本吻合了。因此建议你在设置电路参数时，先算清楚变压器本身的短路阻抗是多少，然后分配到高低压侧的系统阻抗。
Comsol的电路接口支持电阻、电感、电容、开关等元件，你要把短路操作模拟成某个时刻开关闭合或者理想变压器短路。我用的是“电路”接口，在模型中加入一个可控开关，在指定时间后闭合，从而实现三相短路。这里有个容易出错的地方：开关的闭合时刻要和前面提到的短路相位角对应上，因为电压源相位是随时变化的，必须用正弦函数配合时间偏移来设定电源的初始相位。
3. 电流、电压、磁密的提取姿势：后处理里最容易忽略的细节
3.1 磁密到底是看峰值还是看分布云图？
后处理时我们习惯抓到某一帧看磁密云图，比如t=20ms时的B分布。但瞬态短路过程里，磁密在空间上不是一个均匀的同步变化，不同位置的磁密达到峰值的时间点并不同步。比如铁心柱中心的磁密峰值可能出现在短路后10ms，而旁轭位置的峰值可能出现在14ms。所以如果你只关注某一个固定时刻，极有可能漏掉真实的最大磁密。正确做法是定义一组“域点探针”或“域最大值探针”，在整个求解时间范围内记录不同区域磁密随时间的变化曲线，然后从曲线上读取每个区域的峰值和对应时刻。
我在模型里设置了四组探针：A相心柱中心、B相心柱中心、上铁轭拐角、绕组主空道中心。这四组探针得到的曲线形态差异很大。铁心柱中心的磁密基本保持正弦形状，但在短路瞬间会有一个大的涌起，幅度比稳态高很多；铁轭拐角处的磁密则出现了明显的畸变，因为短路时大量漏磁通从这里通过，局部饱和导致波形平顶。有了这些时间曲线，你不仅能看到最大磁密是多少，还能看出饱和区域随时间的演化，这就比一张静态云图信息量大得多。
3.2 电压与电流波形的时间对齐问题
从电路角度，三相短路时电压会突然跌落，电流会急剧上升，两者在时间上并不是完全同步的。原因在于变压器绕组是强感性负载，电流滞后于电压接近90度。仿真中如果你把三相电源电压和绕组电流画在同一张图里，需要先注意二者的相位基准。我用Comsol的全局变量分别导出了A相电压源电压和A相绕组电流，画在一起后能看到明显的相位差。这个相位差不是误差，而是物理过程导致的，可以用来验证你的模型是否正确。如果电压和电流的相位差接近0度或180度，那说明你很可能把电路里的参考方向设置反了，或者绕组一致性设置错误。
另外，短路瞬间的电压波形上经常会出现高频衰减分量，这是磁场能量重新分配产生的暂态分量，在数值上表现为电压振荡，会随着时间逐渐衰减。如果你的时间步长太大，这个振荡可能被抹掉，或者出现非物理的持续振荡。我建议在初始阶段把时间步长设小，比如0.01ms，等振荡衰减后再逐渐加大步长，这样既能保证瞬态细节，又能减少总的计算时间。
3.3 从磁密变化反推短路电流冲击
一个很容易忽略的技巧是，利用磁密波形来验证电流计算结果独立于电路求解。因为磁密B和电流I在安培环路定律下应该是自洽的，如果你的电流探针和磁密探针来自同一个模型，它们的一致性是理论上必然的，但如果你的网格太粗，电流和磁密的对应关系就会变差。我在做了几组网格对比后发现，当绕组区域网格尺寸从5mm细化到2mm时，电流峰值变化不到2%，但磁密峰值变化了接近8%。这说明磁密对网格更加敏感，因为磁密是电流的导数型物理量，对空间分辨率要求更高。因此，如果你想重点研究磁密，网格细化的优先级要高于电流。
实际操作中，我通常会在绕组和铁心气隙附近使用边界层网格，并设置最大单元尺寸不超过1mm到3mm。铁心内部可以放宽到10mm，因为那里的磁密变化比较平滑。如果算完发现磁密云图出现锯齿状等值线，那几乎可以肯定是网格不够密，不要急着调求解器。
3.4 绕组电流不平衡度的检查
三相变压器发生三相短路，在理想对称系统里三相电流应该是幅值相等、相位相差120度。但因为实际几何结构很难做到完全对称，尤其是三相五柱结构和各相漏磁路径不同，会导致三相电流存在一定的不平衡。我仿真出来的结果，B相电流峰值通常比A、C相稍微大一点，这和理论分析中B相磁路较短、互感耦合更强的结论一致。如果你发现三相电流差异超过5%，就要检查是不是几何模型里相别搞错了，或者边界条件存在不对称设置。通过记录三相电流的有效值和峰值，还能辅助判断保护装置对故障相别识别的逻辑是否合理。
4. 求解器踩坑实录：网格怎么剖、步长怎么设、不收敛怎么办
4.1 不同短路时刻下网格敏感度对比
我一开始图省事，直接用静磁场模型用的较粗网格，做一个均匀的网格剖分，结果瞬态计算到几十步之后就发散。排查了半天，发现是气隙和绕组内区域的网格太粗，导致磁导率更新时出现振荡。后来我做了网格敏感度对比，分别在绕组区域采用10mm、5mm、2mm尺寸，计算短路后前20ms的电流波形。结果发现10mm的网格在短路时刻出现明显的数值毛刺，电流峰值偏小约6%，而5mm和2mm的结果差异小于1%。这说明对于三相短路这种包含快速瞬变的仿真，绕组区域网格至少要达到漏磁场特征尺寸的1/3以下。我最终选用了2mm，虽然自由度增加了不少，但求解稳定性和结果精度都让人放心。
另外，铁心的叠片结构等效为各向异性电导率后，网格不能画成完全各向同性，需要在叠片方向加密，否则涡流损耗和屏蔽效果会算不准。这一点对绕组涡流损耗评估尤为重要，因为短路过程中绕组中的涡流会引起额外的损耗和发热。
4.2 时间步长对电流峰值的影响：实测对比
Comsol的瞬态求解器默认使用自适应时间步长，但如果你设置的相对容差不够严格，它在电流峰值附近可能会跳过一些关键时刻，导致峰值被低估。我做过一次对比：用固定时间步长0.05ms和自适应步长分别计算，结果自适应步长的峰值比固定步长低了不少。原因是自适应步长在某些快速变化的阶段，以一个较大的步长跨过了极值点。所以建议在短路发生后的一到两个周期内，使用固定时间步长，步长取5e-5秒到1e-4秒，之后可以切换回自适应。我在仿真中设置了时间序列：0到1ms用0.01ms步长，1到10ms用0.05ms步长，10到50ms用0.1ms步长，这样总计算点仍然不多，但峰值捕捉准确。
如果你发现电流波形在峰值附近有明显的三角波形状，那大概率是步长太大，二阶精度退化，导致曲线变成线性插值。加密时间步后，波形会变得光滑，峰值也会更接近理论值。
4.3 不收敛的经典根因：B-H曲线斜率为0与初始值失配
非线性瞬态求解不收敛，最常见的原因是B-H曲线在某个数据点出现斜率突变为0，导致雅可比矩阵奇异。Comsol里默认的插值方法在数据点稀疏时会产生过度震荡，尤其是超出数据范围时会线性外推，而外推的磁导率可能变得非常低，使得迭代过程在铁心区域出现严重的数值振荡。解决方法是把B-H曲线的数据范围设置得足够宽，并且在材料设置里把最大相对磁导率和最小相对磁导率限制一下，比如最大设为1e5、最小设为1.1，这样可以避免某些区域磁导率下降到接近真空的水平。同时要确保铁心区域初始磁通密度为0，但绕组电流初值也要为0。如果初始值里残留了稳态磁场，瞬态一开始就会产生虚假的电磁冲击，导致不收敛。
还有一个容易忽略的问题：绕组线圈的“线圈长度”或“线圈截面积”设置必须与实际几何匹配。如果线圈截面积设小了，等效电阻过大，电流波形会衰减得比真实情况快，严重时会导致求解器因为电流过小而失去耦合。我会在建模时把几何尺寸直接参数化，然后从几何中提取线圈截面积，避免人为手填出错。
4.4 一次真实的“负磁密”排查过程
我调试过程中遇到过非常诡异的现象：铁心某些区域的磁密出现了负值，而且负区域的形状随着时间不断变化，看起来像波动。一开始我以为是材料B-H曲线定义错了，后来一步步排查，发现是网格坐标系导致的。铁心材料我设置了各向异性电导率，需要指定轧制方向。我错误地把轧制方向设置在全局坐标系的X方向，但铁心在不同位置的方向不同，这就导致材料属性在某些区域定义错误，磁密计算出来方向混乱。修正办法是在每个铁心域单独设置材料方向，或者使用材料坐标系来定义轧制方向。改完之后，磁密分布就完全正常了。
如果你也遇到类似的负磁密或磁密非物理分布，先别怀疑物理模型，优先检查材料方向、坐标系和网格的雅可比行列式。另外，铁心叠片方向在3D模型中的等效电导率设置，也会强烈影响涡流路径，一旦方向错，结果里会出现明显的局部磁密“坑”。
4.5 求解器的直接法还是迭代法？
Comsol里磁场模块可以用直接求解器（MUMPS）或迭代求解器（GMRES）。对于三维非线性瞬态模型，我强烈推荐MUMPS。迭代求解器在瞬态非线性问题里很容易因预处理不完善而不收敛，尤其在你的模型包含电路耦合时，电路变量和场变量的尺度差异巨大，迭代求解器很难同时处理好。MUMPS虽然内存占用高，但稳定性好，而且往往能一次收敛。我试过用GMRES跑了4小时还没算完第一个毫秒，换成MUMPS后总计算时间反而缩短了。所以对于有电压源和大量非线性材料的模型，不要犹豫，直接选MUMPS，内存够就上。
如果你的模型真的太大，MUMPS内存不够，可以用分离式求解器，把磁场和电路分开迭代，但需要手动设置容差和最大迭代次数，而且很容易出现不收敛。我的经验是，老老实实买一台内存大的机器，或者把网格稍微放粗，也比折磨求解器强。
5. 几个让结果更有说服力的后处理小技巧
5.1 用“最恶劣相位角”作为最终输出工况
我在做项目总结时，会同时给出三个不同短路相位角下的电流峰值和磁密峰值，并明确标注出最恶劣工况对应的相位角。这样不仅展示了仿真的全面性，也让结构工程师或者保护工程师可以直接取用最不利数据。不要只给一张漂亮的磁密云图，要给出一个包括相位角、最大电流、最大磁密、最大电磁力的表格。你可以用Comsol的参数化扫描功能批量计算，然后导出结果到一个表格里，这样工作会高效很多。
5.2 磁密云图和磁力线叠加显示
后处理时我喜欢把磁密云图作为底图，同时在叠加磁力线场图。在短路初期，你会看到磁力线大量从铁心表面逃逸到空气和绕组区域，这正是漏磁通异常增大的表现。抓住这个时机的磁力线分布，能直观解释为什么绕组的端部往往受到最大的轴向力。如果只看磁密云图，这样的受力分析线索并不明显。你也可以在绕组区域做一个切面，单独显示该面内的电流密度和磁密法向分量，用来评估绕组导体表面受到的应力分布。
5.3 导出数据时别忘了保存所有时间步的探针数据
很多人在后处理中喜欢直接截取云图，但我更建议把探针数据和全局变量全部导出成CSV，然后用Matlab或Python做二次处理。因为Comsol后处理画图功能虽然方便，但要画多条曲线对比时还是外部工具更方便。比如你把A、B、C三相电流叠加判断对称性，或者把磁密峰值曲线包络出来看衰减速度，这些在Comsol里也能做，但可能不如做Python脚本灵活。我一般会在导出时勾选“所有时间步”，再在后处理阶段用程序绘制波形，这样后续调整图例、标注幅值都更方便。
6. 最后，关于这套模型还能怎么扩展
如果你已经成功跑通了上面这套对称三相短路仿真，可以继续把模型往两个方向扩展：一是考虑变压器分接开关位置变化，不同分接档位对应不同的短路阻抗和匝比，短路电流也会不同；二是考虑外部线路参数不对称，比如单相接地短路或两相短路。Comsol的电路接口支持你随意改变故障类型，只需要修改开关连接关系即可。另外，如果你关心短路后的热效应，可以把磁场的损耗计算结果导入固体传热接口做顺序耦合，看看绕组温升是否在允许范围内。
我个人做下来的最大体会是，三相短路仿真的成败并不在于软件功能强大与否，而在于对物理过程理解得是否透彻。每一个参数设置，比如短路相位角、B-H曲线、时间步长，背后都有物理依据。不要盲目追求网格多细、模型多精美，先理解每一步在算什么、为什么这样算，才能得到真正可信的结果。希望这篇内容能帮你在自己的仿真项目中少走几条弯路。
