1. 为什么非要算这个“热-流-磁”的耦合账油浸式变压器的热点温度几乎所有做变压器的人都知道它重要但它恰恰又是最难测的一个量。表面温度计贴得到光纤测温布得进的测点也就那么几个但绕组内部最热的那一点藏在绝缘纸和铜导体之间平时根本看不见。我最早接触这个课题就是因为一次温升试验里的困惑顶层油温明明没过限甚至比预计低了好几度可拆开后的绝缘老化痕迹却比预期重得多。后来才意识到那台变压器的绕组热点落在了油流不畅的角落油温不高不代表绕组温度不高。这个现象背后的逻辑要靠多物理场耦合仿真才能讲清楚。过去工程上常用的办法是查标准、查经验公式比如用顶层油温升加一个固定的热点系数去估算绕组热点温度。这种方法对于常规结构、常规冷却方式是够用的可一旦碰上非标设计、强迫油循环、大容量高负载率这些场景经验系数的外推空间就非常有限了。我去做这个仿真项目的初衷就是想把“看不见的热点”变成“算得出来的分布”同时把绝缘油在变压器内部怎么流动、流速大概是多少也一并弄清楚。油流速度分布这件事很多做变压器热计算的人会忽略但它的作用一点都不比温度小。绝缘油在绕组油道里怎么走直接决定热量能不能顺利带走在油箱内部哪些区域形成流动死角哪些地方流速过低都是局部过热的隐患来源。把电磁损耗算出来再把它作为热源丢进温度场让温度分布反过来驱动绝缘油的浮力流动流动又改变温度场——这三者闭环起来才是变压器内部真实发生的过程。这套仿真做下来能解决的问题很具体绕组热点温度具体在哪一匝、偏内还是偏外油流在每层油道里流速大概是多少、是否存在低速区负载率提高以后热点温升大概按什么规律变化。适合的读者大概分三类变压器厂做电磁热设计的工程师做输变电设备状态评估的科研人员还有正在用COMSOL做多物理场耦合课题的研究生。我这里分享的都是自己实际跑过的模型和踩过的坑你可以照着搭也可以结合自己的结构改。2. 三场耦合的能量传递链路从电磁损耗到热油上浮2.1 变压器内部的三个损耗源头先得算清楚做耦合仿真最忌讳的事是几何模型建得漂漂亮亮热源却随便填一个经验值进去。这样算出来的温度场哪怕收敛了也是自欺欺人。电磁损耗是整条耦合链路的起点它算得准不准直接决定了后面所有结果的可信度。油浸式变压器内部的热量来源主要是三个绕组铜损。电流流过绕组导线时的I²R损耗。这里要注意的是铜导体的电导率随温度变化非常明显常温下铜的电导率大约在5.998e7 S/m温度系数约为0.0039/K。也就是说绕组温度从20℃升到90℃时电导率会下降超过20%对应的损耗也会同比例变化。这就是为什么严格的做法必须让电磁场和温度场双向耦合——温度变了电导率变了损耗分布跟着变热点温度又会进一步变化。铁心损耗。铁心中的磁滞损耗和涡流损耗取决于铁心材料的B-H曲线、叠片厚度和工作磁密。在COMSOL里可以用磁场的频域求解来获取硅钢片叠片结构本身很难逐片建模通常用等效的层状电导率来处理——也就是叠片方向电导率极低片内方向电导率正常。杂散损耗。漏磁场在油箱壁、夹件、压板等结构件上感应的涡流损耗。这部分在高电压大容量变压器里特别值得关注因为漏磁场强结构件离得近局部涡流可能相当集中。不过做模型时通常会做取舍关注热点温度的绕组区域必须精细建模结构件可以先简化等主体模型跑通后再加回来。三条损耗路径在COMSOL中的处理方式不太一样我整理了一个对照表损耗类型产生位置COMSOL中获取方式对热点影响绕组铜损高、低压绕组导线磁场接口的Coil特征积分Ohmic损耗密度最直接热源主要来源铁心损耗铁心叠片磁场频域解材料损耗模型影响铁心附近油温和局部油流杂散损耗油箱壁、夹件、压板磁场接口全域求解体积积分离绕组远时影响小近时不能忽略2.2 热量在绕组和铁心里面怎么传递损耗算出来之后接棒的是传热分析。油浸式变压器内部的传热路径很典型铜导体产生的热量先通过绝缘纸和导线绝缘层传导出来再穿过绕组表面的油流边界层交给流过的变压器油带走。绝缘纸的导热系数很低通常只有0.15到0.25 W/(m·K)铜的导热系数是400 W/(m·K)左右两者相差超过两个数量级——这意味着温降的大部分发生在绝缘层和油流边界层上。所以在几何建模阶段如果图省事直接把绕组建成一整块铜不单独划分绝缘纸区域算出来的热点温度会明显偏低。这是我做了好几轮模型后摸索出的一个重要细节热量从铜导线导到绝缘纸表面时绝缘层的温度梯度在总温升中的占比往往能达到20%以上这一层必须建出来。变压器油侧的传热以自然对流为主。对于强迫油循环的变压器油流速度高可说白了还是对流传热——只是从浮力驱动变成了泵驱动。但无论哪种传热系数都强烈依赖油流速度、黏度和油道几何尺寸。这也解释了为什么不能只算温度场不算流场如果把油侧的换热系数作为经验常数直接给定那等于预设了答案完全忽略了绕组结构对油流路径的影响。2.3 绝缘油为什么自己会动起来绝缘油在变压器里的流动主要靠密度差这个“发动机”。靠近绕组和铁心的油被加热后体积膨胀、密度降低在重力作用下往上浮温度低的油密度大往下沉。这种浮力驱动的流动在流体力学里叫自然对流COMSOL里用层流接口配合重力项就能算。这里有个非常关键的选择用可压缩流体直接算密度随温度变化还是用Boussinesq近似。矿物绝缘油的体积热膨胀系数大约在7e-4 1/K的量级温度变化几十度时密度相对变化也就几个百分点。做全可压缩计算容易引起压力-速度耦合的收敛困难计算开销也大。实际项目中我用的做法是Boussinesq近似——密度只在浮力项里随温度线性变化连续性方程和动量方程里的密度当作常数。这样算出来的速度场和温度场在变压器油自然对流的场景下精度足够收敛性却好得多。油的黏度随温度的变化也必须给足。矿物油在20℃时运动黏度可能接近几十厘斯到90℃时降到不到原来的五分之一。黏度直接决定边界层厚度和对流换热强度很多新手容易忽略。我见过有人在论文里把油黏度设成常数算出来的流速和真实情况差了三四倍这就是材料参数给错的典型后果。另外绝缘油黏度做适当截断是必要的仿真中低温区域的油可能黏度大得离谱流动几乎冻结这会带来不必要的数值刚性问题——可以在绝对黏度上加一个下限比如5e-4 Pa·s不会影响精度。流场和温度场的耦合逻辑粗略可以概括成一句话温度差异产生浮力浮力驱动流动流动改变换热换热又反过来重新分配温度。加入电磁场之后整条链路就闭合了电流产生损耗损耗作为热源抬升温度温度影响材料属性材料属性再反向作用于损耗和流场。所以在这个问题上“耦合”不是软件里点一个按钮那么简单它是物理过程中天然存在的仿真只是尽量去还原它。3. 从几何建模到收敛搭建模型的完整实操记录3.1 几何建模阶段的取舍别一上来就建3D整站做变压器耦合仿真最容易犯的错误是第一次建模就想把整台变压器按图纸一比一还原。等你花了两周把几十个零件都画出来之后就会发现光是网格剖分和求解就已经让机器吃不消了。我的经验是先想清楚你关心的物理量是谁再决定哪些几何细节必须保留哪些可以简化。对于热点温度和油流速度分布这两个目标绕组是绝对的主角铁心、油箱、绝缘油是必要的舞台其他的附属结构可以后补。绕组建模时需要注意匝与匝之间的绝缘、层间油道都必须体现因为油道的位置和宽度直接决定油流能带走多少热。具体到几何维度我推荐两套方案方案A快速评估二维轴对称模型把铁心柱、高低压绕组、主油道、油箱抽象成同心圆柱体系。这套方案计算量小单次求解几分钟到十几分钟适合参数扫描比如研究不同负载率下热点温度的变化规律。缺点是没法考虑铁心夹件等结构的局部影响油箱处流动细节也被简化了。方案B精确分析三维周期模型或全模型。取一个对称扇区或一个绕组段作为计算域保留油道的真实几何形状可以看见油流在油道里的真实偏转和低速区分布。网格量会从几十万跳到几百万需要工作站和合理的求解策略。我个人的建议是路线不要跳太快先用二维轴对称把物理过程跑通确认损耗量级合理再考虑扩展到三维。3.2 材料参数表照这个填基本不会翻车材料参数是耦合仿真的地基而且必须区分哪些是常数、哪些必须给温度依赖关系。我踩过的坑基本都发生在参数不准确或忘给温度相关性的时候。这里给一份直接可用的参数表供参考材料参数取值是否随温度变化铜绕组电导率5.998e7 S/m20℃基准温度系数0.0039/K是必须给铜绕组导热系数400 W/(m·K)否绝缘纸导热系数0.18 W/(m·K)近似常数硅钢片铁心叠片方向电导率极小取各向异性否硅钢片铁心导热系数30~40 W/(m·K) 面内5~10 W/(m·K) 叠厚方向否绝缘油密度870 kg/m³参考温度下Boussinesq项中用热膨胀系数绝缘油动力黏度0.008 Pa·s约在80℃需查油品数据是强烈变化绝缘油导热系数0.13 W/(m·K)弱可设常数绝缘油热膨胀系数7e-4 1/K否绝缘油比热容1800~2000 J/(kg·K)近似常数填充材料属性时最容易忽略的是铜的电阻率温度系数。如果我把铜电导率设成常数绕组损耗反馈到温度场的效果就断了热点温度的计算偏差会随着负载率提高越发明显。绝缘油的黏度曲线建议查具体牌号油品数据表或者用经验拟合式不要用单一固定值。3.3 边界条件怎么给才能既贴近物理又能收敛边界条件的设置是仿真“灵魂”所在直接决定你算出来的是物理过程的映射还是一个看起来漂亮的数学解。电磁场边界铁心柱和铁轭的磁导率设置为非线性B-H曲线如果手头没有精确曲线先给定一个常数相对磁导率比如2000~3000也可以对热损耗的分布影响有限。绕组激励建议使用电流或电压激励用COMSOL的Coil特征。我推荐用已知的负载电流作为激励这样可以避免多绕组的电压方程引入额外的收敛难度。计算域外边界通常是磁绝缘取在油箱外足够远的位置或者直接把油箱壳作为边界。温度场和流场边界油箱外壁对空气的自然对流散热表面换热系数取5~10 W/(m²·K)这取决于油箱的实际散热面积与环境风速。绕组内边界和铁心表面默认是流固耦合界面不需要额外设置速度边界层流求解器会自动处理壁面无滑移和热流连续。如果模型包含开口油道段必须在入口和出口设置压力边界条件或循环边界。二维轴对称模型中如果截断了散热器回路通常用“开口边界压力为0”来近似外接油路。重点提醒自然对流必须开重力而且重力方向要和几何模型的坐标轴对应。经常有人搭建二维轴对称模型时忘记加重力或者重力方向反了算出来的油流不是循环而是发散的这就是一个个典型的“跑得通但不说物理”的案例。加上重力后动量方程里要挂一个体积力项——COMSOL里可以直接在流体接口开启“浮力”选择Boussinesq近似并指定热膨胀系数和参考温度。3.4 多物理场耦合节点的设置逻辑COMSOL里做这种耦合有两种常用路径。一种是在“多物理场”节点里直接用预置的耦合——电磁热Electromagnetic Heating把磁场接口和传热接口绑定会自动把焦耳热和磁损耗作为热源送进传热方程非等温流Nonisothermal Flow把流体流动和传热绑定考虑温度对流动的影响以及流动对温度的对流项。然后把温度场的浮力源项、通过流场影响温度分布这块是在多物理场节点中完成的。另一种路径是手动耦合用耦合算子或者通过共享变量来传递数据。这个比较麻烦一般不推荐新手一上来就用。预置的耦合虽然“看不见底层”但由底层公式处理的变量关系更稳健。求解方式上如果要实现双向耦合——也就是温度变化后电导率变化再影响电磁场——就必须让电磁场和传热、流体模块一起求解。这里有个性能问题电磁场的频域求解通常是复数线性方程传热和层流是实数非线性问题放在一个全耦合求解器里会非常吃力。我的做法是使用分离式求解步骤Segregated solver先解电磁场再把损耗密度映射为热源然后解温度和流场更新电导率后再回头解电磁场如此循环直到收敛。这种“外部迭代”的方式在物理上等价于双向耦合数值上却比全耦合稳定得多。4. 热点温度和油流速度到底怎么取出来用4.1 热点定位不是看最大温度云图那么简单的当模型收敛后你会在结果里看到温度分布云图但“热点温度是多少”这个问题的答案远不是看一眼图上的红色最大值就完事。第一步是确认“体域最大值”出现在哪儿。油浸式变压器的热点常见位置是绕组轴向高度的70%~90%靠近上部的位置偏向导线径向的内侧或外侧取决于油流方向。可以在三维结果里直接做“体最大值”探针锁定坐标后手动查看那个位置属于哪一匝、靠近哪个油道。第二步要用截面分析。沿绕组轴向切一个剖面看等温线分布再沿径向做一条水平线输出温度沿径向的曲线这样能判断热点是在绕组的表面还是内部。很多情况下云图显示最热的是贴近油道外侧的一小块区域但那是边界层里的油温不是绕组导体的温度两者的差可能达到5~10K。提取结果时把固体域和流体域分开显示很关键这样才能准确抓到导体内部最热的那一点。第三步是看轴向温度曲线。通常在二维轴对称模型里沿绕组中心线输出温度随轴向坐标的变化你会看到一条从底部到顶部逐渐升高、末端又微微下降的曲线。曲线的最大值位置是热点所在曲线段的斜率变化还能反映油道设计的优劣。如果曲线中间出现突变凸起那通常意味着那个位置的油道流速很低、散热不畅需要检查油道设计。4.2 油流速度分布关注的不只是数字大小油流速度的结果分析直接关系到你对散热能力的判断。变压器自然对流下的油流速度并不快。绕组油道内的典型流速数量级在每秒几厘米左右而不是每秒几米。很多第一次做这个仿真的人看到这个数字会觉得是不是算错了——实际就是这个量级还是那句话自然对流靠浮力驱动驱动压头本来就小。速度场可视化我建议同时看三类图速度流线图可以看到油从底部进入绕组油道、受热上升、到达顶部后从油道出口流出的整体路径。如果流动路径中出现明显的回流或旋涡那个位置就是潜在的热点伴生区。速度矢量截面图重点看油道入口的流速大小。由于底部入口段和顶部出口段存在局部阻力速度往往不是均匀分布而是偏向某一侧。这种偏流现象会造成局部油流速过低是热点形成的重要原因之一。壁面附近的边界层速度剖面油流在壁面处趋于静置边界层厚度决定了换热强度。如果你在COMSOL中开启了层流接口可以在后处理中查看壁面的热通量它会直接反映哪个部位散热好、哪个部位散热差。速度结果还有一个工程层面的用法配合温度场算局部换热系数。在结果导出中把壁面热通量除以当地温升就能得到沿绕组高度方向的换热系数分布。这个数据对后续独立做结构优化、甚至做热网络模型标定都非常值钱。4.3 热量守恒校验仿真到底靠不靠谱先用这个办法自查耦合仿真做出来的温度云图和流线图再漂亮也必须经过校验才能拿去写报告、改设计。我最常用的校验手段是热量平衡核算。方法很简单在后处理中分别对绕组、铁心区域积分体积损耗把总损耗P_total算出来然后在边界上比如油箱外壁积分散热量Q_out。稳态收敛下两者之差应该在百分之几以内。如果偏差超过10%大概率是某些边界条件没加对、部分区域的损耗被重复计入或漏掉了甚至可能是求解并未真正收敛。另一个校验方式是跟标准温升试验公式对表。基于GB/T 1094.2或IEC 60076.2的温升计算方法用仿真得到的平均油温升、绕组平均温升做对比。两者的误差只要在合理范围内我认为10%~15%以内是可以接受的就说明仿真模型的整体物理量落到正确量级了。最后一个手段是网格无关性验证。把网格加密一倍看热点温度变化是否小于1~2K。如果加密后温度变化很大说明旧网格没有解析出边界层的温度梯度热点温度的结果是“伪收敛”。这个步骤虽然在赶进度时总想跳过但一旦要做结论它是不能省略的。5. 折腾两周后我总结的收敛与提速技巧5.1 网格划分策略哪里加密、哪里放松是门手艺耦合仿真的网格剖分最忌讳“雨露均沾”。变压器模型跨度极大从油箱的米级尺度到绕组匝间绝缘的毫米级甚至亚毫米尺度都存在于同一个计算域。如果整体都按细网格画机器跑不动按粗网格画油流边界层又完全丢失。我最终采用的策略是分区控制绕组区域加密铜导体和绝缘纸至少要剖分两层以上确保径向热梯度被解析。近壁区域所有油固交界面必须画边界层网格通常拉5层以上、首层厚度取预估边界层厚度的十分之一。首层厚度可以根据无量纲壁面距离y估算自然对流场景下首层厚度约0.5~1mm网格生长率1.2~1.3。油箱区域远离绕组的地方用粗网格因为那里的温度梯度和流速梯度都很小。铁心内部也可以适当放松热点不会出现在铁心中心。油道入口和出口段网格过渡要平滑不要出现尺寸突变否则会形成数值的伪流动扰动。一个常见的标的网格量二维轴对称模型大约5万到20万单元三维扇区模型可能到200万到500万单元。如果超过这个量级很多先检查是不是简化得不够。5.2 遇到不收敛先从这五个方向排查这个模型我在求解过程中遇到过无数次的“迭代发散”和“残差不降”。总结下来收敛问题大部分出在以下五个方面优先级从高到低初始值给得太激进。一上来就让流场从零速度起步直接算自然对流动量方程很容易发散。解法是先关闭浮力项算一个纯导热的温度场作为存储的初始值或者先不加流体把温度场初步算出来然后再用这个温度场作为流场和温度场耦合计算的初值。这招屡试不爽。黏度过低或流动进入湍流区。虽然自然对流的油流速度很慢但如果油道特别窄、局部流速高局部可能进入过渡流。变压器内部一般仍按层流处理但速度过高时不要强行用层流模型——这时观察局部雷诺数如果超过临界值说明该位置的流动模型需要调整。边界层网格质量差。网格出现大角度倾斜或者相邻单元尺寸比过大会导致动量方程在边界处刚性极强。先修网格而不是硬调求解器。多物理场耦合的反馈太强。温度一个迭代步变化过大会把流场剧烈改动导致发散。我常用的做法是在分离式求解器中把温度方程和流场方程的松弛因子调低比如从默认的1改到0.7左右。Boussinesq近似的参考温度设置不当。如果参考温度不在实际工作温度范围内浮力项会给出很大的虚假驱动力流场会失去稳态。把参考温度设为平均油温或者设为环境温度加预计温升的一半这个“工程直觉”参数会稳很多。5.3 算力不够时的降级路线很多来问我的同行和学生都是在自己的笔记本上跑模型动辄几百万网格的3D模型确实吃不消。这里有一条亲测可行的降级路线按序执行第一步全部用二维轴对称模型跑通物理过程。几何虽然简化了但耦合机制、材料参数、边界条件设置和三维一致。此时你已经能获得热点温度的大致范围和油流速度的数量级。第二步做参数扫描时继续用二维模型比如研究负载率、油道宽度、油泵流量对热点的影响趋势。此时你会得到一组相对变化规律这个规律比绝对数值更能用于设计判断。第三步锁定一个最关心的工作工况再建立三维扇区模型。用二维模型的温度场结果充当三维模型的初始化值收敛难度会显著降低。第四步如果想做精细化热点定位再把几何从扇区模型升级为全模型同时打开电磁场与热流体的双向耦合。此时前期积累的计算经验就派上用场了。6. 最初跑偏的几个方向后来都成了经验这个课题做了快一个月之后回头看最耗时间的不是建模本身而是对物理过程的理解和对结果合理性的判断。这里把踩过的弯路集中记录一下希望后来的人少绕几圈。第一个教训是不要迷信全耦合求解器。最初我图省事建好模型直接开启“电磁-热-流”全耦合觉得一步到位很酷。结果一整个下午都在等迭代不但发散还发散得毫无规律。后来改成电磁场单独求解、热流场耦合求解再把两者放在外层循环稳定性和速度都有了质的提升。这种“分而治之再闭合回路”的策略在复杂的工程仿真里永远比所谓的“全耦合”更务实。第二个教训是结果验证不能只看“形状像不像”。一开始做出来的流线图看起来非常完美但后来算热量平衡时发现总损耗和总散热量差了将近20%。追查后发现我将绕组损耗密度乘以了频域结果的幅值而不是有效值导致热源翻倍。这类隐藏在“看起来正常”结果里的系统性错误如果只靠视觉判断很容易蒙混过关。第三个教训是网格无关性验证不是可有可无的步骤尤其是对于热点温度这种对局部网格非常敏感的物理量。有一次我把绕组绝缘层的网格从2层加密到4层热点温度就变了将近3K。如果不做这个验证我当时就得拿着偏温度的结果去下设计结论。再分享一个实际操作中的小技巧在COMSOL里做流场后处理时速度场常以箭头或流线形式展示但如果想定量分析某个油道截面上的流速分布我建议在那个位置显式定义一个“截面积分算子”直接用“表面积分/面积”求平均速度再用截面上的最小值和最大值判断流动均匀性。这个操作比在云图里用颜色读取数值靠谱得多。最后一件事是关于油流速度分布的实用理解。很多做电力设备的朋友第一次看到仿真里的油流速度只有每秒几厘米时会产生不信任感——“这么慢还能散热吗”实际上自然对流的驱动力本来就弱油流“慢”恰恰是合理的。判断散热能力时不要单看流速绝对大小要看两个指标一是绕组油道内温度边界层的厚度二是油道出口油温与入口油温的差值。只要流量和温升搭配合理慢速油流照样能有效带走热量。理解了这一点你对仿真结果的价值判断就有了真正的抓手。