做变压器仿真的同行应该都有这种感觉稳态场算起来一切岁月静好磁场分布均匀、磁密云图漂亮可一旦切换到三相短路工况模型就开始各种闹脾气。我最近刚把一个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曲线没有它就不要谈磁密变压器铁心正常工作时是工作在饱和点附近的尤其是短路时铁心会进入深度饱和所以材料参数绝不能给一个固定的相对磁导率。我见过太多人直接在材料节点里填了一个μr2000结果算出来的磁密分布跟实际严重不符因为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 磁密到底是看峰值还是看分布云图后处理时我们习惯抓到某一帧看磁密云图比如t20ms时的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曲线、时间步长背后都有物理依据。不要盲目追求网格多细、模型多精美先理解每一步在算什么、为什么这样算才能得到真正可信的结果。希望这篇内容能帮你在自己的仿真项目中少走几条弯路。