COMSOL三相电力变压器电磁场和电路耦合计算这个方向我断断续续摸索了大半年。刚开始做的时候真是抓瞎参考手册翻烂了案例库里的模型要么太简单要么和实际工程场景对不上。后来总算把一套能用的仿真流程跑通了能稳定输出三相变压器高低压绕组的电压电流分布、铁心磁通密度分布以及励磁涌流波形。今天把整个思路、关键设置和踩过的坑一次性整理出来希望能给正在做类似问题的朋友省点时间。先说明一下这套方法我是在COMSOL Multiphysics 6.x版本下验证的物理场接口用的是AC/DC模块里的磁场magnetic fieldsmf和电路electrical circuitcir三维瞬态求解。如果是二维模型或者只做频域分析流程会有差异后面我会在相应位置单独说明。1. 整体思路为什么一定要做电磁场和电路的双向耦合1.1 变压器仿真的两种典型做法搞变压器仿真的朋友应该都知道传统的路算法磁路法是把变压器简化成等效电路和磁路用集中参数来算电压电流。好处是快几秒钟出结果但代价是丢掉了空间分布信息——你想看绕组某个位置电流密度是不是过大、铁心某个角落是不是饱和了路算法完全无能为力。另一种是纯有限元磁场仿真。这种方法可以精确计算磁场分布、铁心损耗但有一个很尴尬的问题绕组激励怎么给如果是恒压源或者恒流源直接加边界条件就行。可实际变压器是接到电网上的电压源通过绕组产生磁场磁场的感应又会反过来决定绕组电流——这是一个强耦合的物理过程。如果人为把激励设成理想电压源励磁涌流、绕组间的电磁耦合、漏感影响这些关键现象都算不准。1.2 耦合的本质磁场方程和电路方程互相喂数据COMSOL里做电磁场和电路耦合本质上就是把麦克斯韦方程组算出来的感应电压、电感参数跟基尔霍夫电压电流定律算出来的电路状态量互相传递。磁场侧核心方程是磁矢势形式的安培定律[ \sigma \frac{\partial \mathbf{A}}{\partial t} abla \times \left( \mu^{-1} abla \times \mathbf{A} \right) \mathbf{J}_e ]其中(\mathbf{A})是磁矢势(\mathbf{J}_e)是外加电流密度。绕组区域通过“线圈”特征Coil来计算电流密度公式是[ \mathbf{J}_e \frac{N I(t)}{S} \mathbf{e}_w ]这里(N)是线圈匝数(S)是绕组截面积(\mathbf{e}_w)是绕组方向单位向量(I(t))就是实时电流。电路侧绕组两端电压的方程是[ U(t) R \cdot I(t) \frac{d\psi(t)}{dt} ]其中(\psi(t))是绕组所交链的磁链通过磁矢势(\mathbf{A})沿绕组路径的线积分得到。这两个方程在每一个时间步长内反复迭代电路算出电流(I(t))电流给磁场当激励源磁场算出磁链变化率再把感应电压返回给电路。这种耦合我在实际操作中用的是COMSOL的“线圈电路”Coil with Circuit特征把绕组直接连接到电路接口的元件端口上不需要自己手动写耦合方程但理解背后的交互逻辑非常重要不然收敛出了问题根本不知道从哪排查。1.3 解决方案能拿到什么结果做完整套耦合计算你能拿到的结果包括但不限于高、低压绕组各相的瞬态电流波形和有效值高、低压绕组各相的端电压波形绕组内部任意位置的电流密度分布和线包损耗铁心磁通密度幅值分布判定是否局部饱和系统三相不平衡度、谐波分量空载合闸时的励磁涌流特性这些结果对变压器设计优化、保护整定计算、绕组结构校核都有直接参考价值。2. 建模前的关键决策选二维还是三维铁心怎么简化2.1 模型维度选择第一件事就是要决定做二维还是三维。我见过不少人一上来就建三维完整模型结果网格量几十万、上百万单次瞬态计算跑一两天性价比非常低。如果只看高低压绕组的电压电流分布二维轴对称模型或者二维平面模型往往就够了。二维模型计算速度快网格控制也简单尤其在参数扫描阶段能省下大量时间。但二维模型有一个绕不开的缺陷无法准确考虑铁心的三维结构效应比如铁心拐角处的局部饱和、三相五柱式铁心的旁轭磁路这些都需要三维模型才能准确捕捉。我做这个项目的变压器是三相三柱油浸式铁心结构相对标准最后我选择了三维四分之一模型利用对称性把计算域缩减到原来的四分之一。如果你的模型确实无法利用对称性比如三相绕组排布不对称或有特殊结构那就老老实实建全模型算力成本只能认了。2.2 铁心结构和材料的处理三相电力变压器铁心最常用的是三柱式结构A、B、C三相绕组分别套在三个心柱上上下有铁轭闭合磁路。建模的时候需要注意铁心叠片方向会造成磁导率的各向异性通常轧制方向的磁导率远高于垂直于轧制方向。在COMSOL里可以通过设置各向异性相对磁导率来实现但如果你手头没有实测的各向异性B-H曲线先用各向同性近似也行误差在工程可接受范围内。铁心材料的B-H曲线是决定励磁特性准确性的核心。我用的取向硅钢片直接把厂家提供的直流磁化曲线B-H曲线输入到材料属性中。这里有个关键点COMSOL的低频磁场默认材料属性里B-H曲线是插值函数输入时要把数据点插值方式设为“平滑”不然求导时会出现非物理的磁导率波动直接影响励磁涌流的波形精度。2.3 绕组建模方式绕组在三维模型里有两种建模方式。第一种是把绕组简化成空心圆柱体或多边形柱体通过“线圈”特征赋予匝数和激励这种方式下绕组内部不再区分每匝导线适合关注外部电磁特性的场景。第二种是逐匝建模把每匝导线的几何都画出来这种方式精度最高但几何操作和网格划分成本极高一般只在分析匝间短路、局部放电等内部故障时才用。我推荐常规的电压电流分布计算用第一种简化方式。绕组区域定义成实体材料导电率设为铜的等效值注意考虑填充系数铜占整个绕组区域的体积比例这样算出来的损耗和电阻值才贴近实际。3. 物理场接口配置磁场和电路到底怎么接起来3.1 磁场mf接口的关键设置建立一个三维瞬态研究添加“磁场mf”接口。关键设置有四个地方第一个是域设定。铁心、绕组、空气域都要明确指定材料属性和方程形式。铁心区域选“安培定律”节点导电率设为硅钢片的电导率勾选“磁导率”为各向异性或由B-H曲线定义。绕组区域要点选“线圈”特征选择“线型”为“均匀多匝”然后填入匝数、线圈导线截面面积并把导电率设为铜。空气域包括包围变压器的外部空间采用默认的安培定律磁导率为(\mu_0)导电率为0。第二个是边界条件。模型外边界要设置磁绝缘条件即(\mathbf{n} \times \mathbf{A} 0)模拟磁场不会穿透到计算域外部。如果使用四分之一模型还要在对称面上设置“理想磁导体”或“理想电导体”边界条件具体选哪个要看对称面的类型——电场切向分量为0的用理想电导体磁场法向分量为0的用理想磁导体。第三个是“线圈”特征的激励方式。这里要选“电路”类型把线圈连接到电路接口。COMSOL中线圈节点会生成一个外部端口在电路接口里就能直接跟电路元件连线。这一步就是电磁场和电路耦合的关键枢纽选错激励类型后面全白做。第四个是“电流/电压”约束或“单匝线圈”之类的附加特征用于特殊激励场景。我做的是双绕组变压器仿真高压侧接电源低压侧接负载高压绕组和低压绕组都要用“线圈”特征分别连接到电路的不同回路上。3.2 电路cir接口的搭建电路接口本质上是一个SPICE类比的电路求解器可以放电阻、电感、电容、电压源、电流源等元件。注意这里的阻抗元件是理想集中参数而变压器的绕组阻抗已经隐含在磁场计算中所以不要把绕组的直流电阻重复加进去否则会双算。对于电力变压器仿真我通常这样搭建电路高压侧一次侧接入一个三相交流电压源线电压有效值设380V或根据你实际额定值频率50Hz三相互差120°。电压源串联一个小内阻0.01Ω左右来避免电压源直接并联大电感的数值病态这个小电阻不影响结果但能显著改善收敛性。低压侧二次侧接对称三相负载可以是恒阻抗负载比如纯电阻或者阻感负载也可以是更复杂的非线性负载。如果想把负载电压稳定在某个范围可以在电路里加电压调节逻辑不过这是个闭环控制问题和电磁场计算耦合起来会麻烦不少。我做的工况包括空载、额定负载、短路三种。空载工况低压侧开路额定负载接对称RL负载短路工况则是低压侧直接短接注意这种情况下电流极大需要降低电压源幅值来模拟短路试验的较低电压条件。3.3 耦合接口的接线细节在COMSOL中磁场接口的线圈和电路接口的连接是通过“电路端子”实现的。实际操作中每个“线圈”特征会自动在电路接口中生成一个二端子元件端子编号对应线圈的输入和输出。你在电路接口的“电路”节点下方看到类似“C1”“C2”这样的外部连接元件就是代表磁场接口里的绕组。接线时务必区分线圈的同名端极性。线圈的绕向决定了电压电流的相位关系如果高压侧和低压侧的绕向设置不一致仿真出来的电压波形相位就会出错。一个简单检查方法空载情况下高低压侧电压相位应该基本一致或差180°取决于你定义的参考方向如果发现对不上翻转其中一个线圈方向即可。4. 网格划分与求解策略精度和效率的平衡术4.1 网格策略网格质量直接决定了磁场计算的精度尤其在绕组和铁心交界区域、气隙和拐角附近磁场变化梯度大网格必须加密。我的实践方案是铁心区域采用“扫掠”网格或“自由四面体”网格最大单元尺寸控制在铁心柱直径的1/10左右。铁心叠片结构本身不需要逐层划分等效均质即可但沿硅钢片厚度方向的单元层数不能太少一般至少2层。绕组区域绕组内电流密度分布往往不均匀需要加密。优先选用“扫掠”网格沿绕组厚度方向径向至少分成5~8层以保证电流密度和损耗分布的分辨率。绕组与铁心之间的油道/绝缘区域网格尺寸过渡要平滑用“自由四面体”配合尺寸函数从绕组表面到铁心表面设置一个增长率避免单元尺寸跳跃太大导致计算误差。空气外域远离变压器本体时磁场衰减很快网格可以剧烈放大。外边界最大单元尺寸可以放宽到变压器本体尺寸的2~3倍。整个四分之一模型我通常控制在80万到150万自由度之间单次瞬态仿真0.05s步长0.1ms在32核工作站上大概需要4~8小时。如果你要做参数扫描建议先跑二维模型摸清趋势再用三维模型做几个关键工况的精确验证。4.2 瞬态求解器配置研究类型选择“瞬态”时间范围根据你关注的物理过程来定。如果只看稳态电压电流分布从0开始仿真到0.1s就足够了0.1s是50Hz系统的5个周期足以让瞬态分量衰减到可忽略程度。如果要看空载合闸的励磁涌流至少算到0.3s以上因为励磁涌流的衰减时间常数可能长达数十个周期。求解器我用的是COMSOL默认的“瞬态”求解器但需要手动调整以下几个关键参数非线性方法选择“恒定牛顿”最大迭代次数设为25次。容差因子Tolerance Factor不要用默认的1我实际用0.1这样收敛更严格电流波形的高频分量不容易被数值误差污染。时间步进选择“BDF”最大阶数默认是2对于这种多时间尺度耦合问题够了。最大时间步长设置为周期的1/1000.2ms最小时间步长可以很小但注意别设置太紧导致计算卡死我一般设为1e-6秒。初始条件这是一个非常容易被忽视的坑。空载合闸时如果直接从t0开始施加额定电压相当于一个阶跃激励会在铁心中激起很大的直流衰减分量导致励磁涌流异常高。实际中你要么把电压源的幅值做成斜坡上升函数比如0.01s内从0线性升到额定值要么先算一个“初始条件研究”获得稳态磁场初值再作为瞬态仿真的初值。斜坡上升的方法最简单我强烈推荐。4.3 关于弹塑性变形的一点补充看到搜索热词里有“comsol塑性变形用于查找弹塑性应变变量在迭代未收敛”的信息这里顺带说一句。电磁场耦合仿真一般不涉及结构力学问题但如果你后续做变压器绕组短路时的电动力分析就需要把电磁场计算出的力载荷传递到固体力学接口这时候才会碰到材料非线性和弹塑性问题。如果在这个阶段遇到迭代不收敛常见原因是B-H曲线插值不平滑导致磁场突变或者载荷施加过快。建议把电磁场计算和结构力学计算拆成两步先算稳态或瞬态电磁场导出洛伦兹力密度再作为外加载荷进行结构分析。这种做法虽然需要导入导出但两个物理场的收敛难度都大幅降低。5. 后处理技巧如何精确提取绕组电压电流分布5.1 电流波形和有效值的提取仿真完成后电流数据存在电路接口的“外部连接”节点里。在结果节点下新建“一维绘图组”选择“全局”然后表达式填cir.C1.I代表电路接口中第一个外部连接元件的电流即可画出电流波形。C1、C2这样的命名规则可以在电路接口的定义中找到注意对照好哪个对应高压绕组、哪个对应低压绕组。如果要对电流波形做傅里叶分析谐波分解用“派生值”中的“全局计算”选择表达式cir.C1.I运算符选“FFT”可以得到各次谐波的幅值和相位。我实际用这个功能分析负载电流的畸变率效果很好。电压类似表达式是cir.C1.U或cir.C1.v需要看你是画端口电压还是节点电压。如果你在电路里把绕组的一个端子接地了那么端口电压就是相对地电压如果两个端子都是浮动的端口电压就是绕组两端电压差。5.2 绕组区域电流密度分布要查看绕组内部的电流密度分布需要回到磁场接口后处理。在三维绘图组中选择“体”表达式填mf.Jx、mf.Jy、mf.Jz的矢量和或直接填mf.normJ画电流密度幅值。重点观察径向和轴向的电流密度分布。正常设计下绕组电流密度沿径向应该相对均匀但受漏磁场影响内层和外层导体的电流密度可能有差异高频谐波分量下差异更明显。如果看到某个局部电流密度异常偏大可能意味着该处存在涡流集中绕组结构可能需要优化。5.3 铁心磁通密度和损耗分布铁心区域画mf.normB磁通密度幅值可以看到硅钢片内磁通密度分布是否均匀有没有局部过饱和点。三相变压器正常运行时中性点处铁心的磁通密度往往是三相互补的合成磁通接近零如果看到某个位置磁通密度超过了1.8T甚至2.0T那个位置就在深度饱和区域铁损会显著增加。铁心损耗在磁场接口中自带计算选项。在铁心域上添加“铁损”特征输入铁心材料的损耗参数通常用斯坦梅茨系数(k_h)、(k_e)、(\alpha)、(\beta)等描述求解后直接画mf.FeLoss就能看到瞬态铁损密度。如果只想看平均损耗可以先求解完整周期再在“派生值”里做时间平均。5.4 空间磁场分布和漏磁场变压器外壳、油箱壁和结构件上的杂散损耗跟漏磁场直接相关。在全局或切面绘图组中画mf.normH磁场强度或mf.normB可以直观看到漏磁场从绕组端部扩散路径。端部漏磁场强的区域会产生局部过热风险这是变压器设计中的常见痛点。如果模型包含金属结构件比如油箱夹件、螺栓建议额外检查这些区域的涡流损耗。在AC/DC模块中用电流守恒方程计算结构件上的涡流密度在瞬态计算中可以看到涡流的渗透深度、损耗分布等信息。这一步计算量不小但对大容量变压器的散热设计有重要作用。6. 绕组连接方式对仿真的影响Y接还是Δ接6.1 Y型和Δ型连接的处理三相变压器高低压绕组各有Y接和Δ接两种方式共四种组合Yy、Yd、Dy、Dd。不同连接方式对电路拓扑结构有决定性影响。Y型连接时三相绕组有一端共接中性点。在COMSOL电路接口中中性点用一个节点表示三个绕组的一端都连到这个节点另一端分别接A、B、C三相电源或负载。如果中性点需要接地取决于系统接地方式直接在这个节点上加“地”元件即可。Δ型连接时三相绕组首尾相接成闭合三角形。注意Δ接时绕组内部在一个局部的电压环会在绕组中产生零序循环电流。如果高低压侧是Yd连接低压侧三角形绕组中的环流会直接影响高压侧的电流分配仿真时必须如实搭建。我做的变压器就是典型的Yd11接法高压侧Y接、低压侧Δ接。启动时励磁涌流与这种连接方式直接相关仿真出来的涌流波形和现场实测对得上说明电路拓扑搭建是正确的。6.2 相位关系的验证一个非常有效的验证方法是变压器空载运行低压侧开路仿真得到的低压侧电压与高压侧电压相位差应符合变压器组别标号的相位关系。以Yd11为例高压侧A相电压和低压侧a相电压相位差是330°或者说低压侧滞后高压侧30°。如果你跑出来发现相位差对不上优先检查绕组方向定义和同名端接线。这个错了后面所有分析都没有意义所以仿真初期一定要做一个空载验证。7. 常见问题与排查技巧实录以下是我在实际仿真过程中遇到的最常见的几个坑整理成速查表方便各位对照排查问题现象可能原因解决办法瞬态求解一直不收敛提示“找不到一致的初始值”电压源从0突跳施加到额定值改用斜坡上升激励或先计算稳态初值励磁涌流波形过高高达额定电流20倍以上铁心初始磁链为零合闸相角不合适调整合闸相角电压过零时合闸涌流最大或从静态磁场初始化高压侧和低压侧电压相位差与铭牌组别不符绕组方向设置错误检查“线圈”的绕向必要时反转方向量计算结果有严重的高频振荡电流波形抖动网格太粗或时间步长过大加密绕组区域网格把最大时间步长降为0.1ms三分支系统电流总和不等于零未正确连接中性点或接地检查Y接绕组中性点是否连接是否添加了接地铁心磁通密度严重超出预期B-H曲线插值不平滑或材料参数设置错误检查B-H曲线数据点插值方式改为“平滑三次样条”仿真时间过长一天都没跑完网格量过大或时间范围过长先用等效二维模型或缩短仿真时长调整容差负载电流和手算值偏差超过20%绕组填充系数或导电率设置错误核实绕组等效导电率包含填充系数折算FFT谐波分析结果出现非物理谐波数据时间窗口不是完整周期截取整数倍周期的数据进行FFT分析移除铁心区域后结果几乎不变铁心区域没有激活“安培定律”或磁导率设成了1检查铁心域的材料磁导率确保B-H曲线或相对磁导率值正确这里面最典型、最值得单独说的是“寻找一致的初始状态”这个问题。我第一次跑瞬态时直接撞上这个报错查了很多资料才明白问题不在于模型本身而在于耦合电路系统和磁场系统的初值匹配难度。COMSOL在瞬态计算的第一步就要解一个非线性代数方程组确定所有状态变量的初值。如果电压源突然加一个额定电压相当于电路端强制跳变而磁场端的磁链却来不及跟上这个初值方程组就很可能无解。解决办法就是加一个启动斜坡用平滑的电压上升函数取代阶跃激励。另一个值得注意的点是时间步长的选择。我一开始图计算快把最大时间步长设成1ms结果算出来的电流波形上叠了一层明显的数值毛刺。后来把步长降到0.2ms毛刺消失波形也平滑了。对于50Hz系统0.2ms意味着每个周期有100步对于大部分工程问题这个分辨率已经足够。如果还嫌不够精确可以分段设时间步前两个周期步长设密一些0.05ms确保涌流和暂态过程被精确捕捉之后再放宽到0.2ms。关于收敛性问题还有一个细节如果材料B-H曲线的数据点跨度很大且插值方法使用默认的线性插值可能会在B-H曲线拐点附近产生非物理的微分磁导率即磁导率突变从而引发局部收敛困难。改用平滑三次样条插值通常能解决这个问题。如果手头实测数据点比较稀疏宁可先做一次数据平滑处理再导入。最后分享几个我个人的经验。第一做这类电磁场-电路耦合仿真前处理时间占比可能超过40%几何清理和网格控制一定要舍得花时间后面求解器的麻烦会少很多。第二不要忽略初始条件的作用一个合理的磁场初值能让瞬态收敛快好几倍。第三尽量多和实测数据对标哪怕是用手算的稳态电压电流比对也很有价值能把很多隐性错误提前暴露出来。这个仿真流程后续还可以扩展的方向包括变压器短路电动力与结构耦合分析、温度场与电磁场双向耦合的热分析、GIC地磁感应电流等不平衡直流偏磁下的励磁特性研究。每一步深挖下去都有不少值得折腾的东西后面有机会我再逐个展开聊。