做过变压器设计或温升校核的人应该都有体会光算电气参数不难难的是把“电-磁-热-流”这四件事放在同一个模型里讲清楚。这个COMSOL仿真模型的核心就是围绕变压器的电磁场、损耗计算和流固耦合温度场展开把绕组发热、铁心发热、油流散热、油箱壁自然对流换热这一整条路径完整还原出来。这篇文章我会从物理本质、建模策略、实操设置到收敛调试把我做这类模型的完整思路和踩过的坑一次讲透适合正在做油浸式变压器温升仿真、或者想入多物理场耦合的新手参考。先说明一点变压器流固耦合温度场本质上是一个“热源在固体里冷却介质在流动两者通过壁面交换热量”的强耦合问题。发热来自电磁损耗损耗大小又和温度有关铜的电阻率随温度升高而增大油的粘度也随温度变化反过来影响对流强度。这一连串互相咬合的物理过程正好是COMSOL最擅长的场景。但“能做”和“做得好”之间隔着大量细节下面逐个展开。1. 先把物理过程拆清楚这个模型到底在解什么1.1 变压器内部的热源与散热路径油浸式变压器的热量来源主要有三块。第一是绕组铜耗也叫负载损耗电流流过导线电阻产生和电流平方成正比第二是铁心损耗包括磁滞损耗和涡流损耗由交变磁通在硅钢片中产生铁心损耗基本只和电压、频率有关和负载关系不大第三是杂散损耗漏磁通在油箱壁、夹件、压板等结构件中感生出涡流造成的附加损耗大容量变压器里这块不可忽略。散热路径则是这样走的绕组和铁心先通过固体导热把热量传到表面再由表面的绝缘油通过对流把热带走。油受热后密度变小自然向上浮形成从底部上升、顶部冷却后下沉的自然对流循环热量最终通过油箱壁传给外部空气。绕组内部还有绝缘纸和油道热量从导线穿过绝缘层再到油道的路径上温度梯度很大热点温度通常在绕组中部偏上的位置而不是最顶端。搞清楚这三条热源路径和一条散热路径模型里需要哪些物理场已经清楚了电磁场算损耗固体传热算绕组和铁心里的温度分布流体场算油的流动传热接口负责把固体和流体两侧的温度和热流耦合起来。1.2 为什么不能只算温度场必须做多物理场耦合有人可能会问我直接把损耗加载到固体上给油侧一个等效对流换热系数不也能算出温度吗对很多工程校核确实是这么干的用经验公式估一个表面换热系数算出来的热点温度在初步设计阶段够用。但问题在于换热系数本身是温度场和流场共同作用的结果不是一个固定值。油的自然对流强度取决于壁面温度和油温之间的温差、油的物性参数以及流道几何。绕组不同位置的温度不一样局部的换热系数也不一样负载变了发热量变了油的流速就跟着变换热能力也在变。用一个常数换热系数去近似等于把最关键的散热机制给钉死了算出来的温度分布自然不可信。要做准确的温升预测就必须把流场解出来让热量传输路径自己“跑”出来这就是流固耦合存在的意义。而且损耗本身也是温度的函数。绕组铜的电阻率随温度升高而增大温度越高损耗越大损耗大了温度再升高形成正反馈。不把电磁场和温度场耦合起来损耗就只能是常温下的定值和实际情况偏差很大。这也是很多新手做变压器温升仿真时忽略的一点。2. 电磁场与损耗计算热源的精度决定温度场的精度2.1 磁场接口选择与几何简化COMSOL里做变压器电磁场AC/DC模块下的“磁场mf”接口是首选频域求解。用频率域的好处是直接拿到稳态下的损耗不用等瞬态波形稳定下来。模型的几何可以简化铁心用硅钢片叠片后的等效实体绕组用多匝线圈域代替不再逐匝画导线。但有一个关键点我想提醒一定要评估“二维轴对称、二维平面、三维全模型”三种方案怎么选。常见的油浸式配电变压器是三相三柱结构但绕组总体呈轴对称分布很多工况下可以用二维轴对称模型计算量小一个数量级。如果你主要关心温度场而电磁场只是用来算损耗的那完全可以用二维电磁场算出损耗分布再映射到三维热流耦合模型里。我自己做单相变压器或者三相变压器单相分析时基本都是这个路子实测下来温度和温升试验结果对得上。三维电磁场模型主要用于需要精确考虑杂散损耗、漏磁在结构件中分布的场景或者油箱壁涡流损耗占比很高的大容量变压器。三维模型的网格量很容易上百万加上流固耦合常规电脑基本跑不动所以一开始就要想清楚取舍。2.2 绕组损耗的准确计算方法绕组损耗分直流损耗和交流损耗。直流损耗好算就是I²RR按工作温度下的电阻率修正。交流损耗则要麻烦一些因为交变漏磁场会在导线内部引起集肤效应和邻近效应让有效电阻比直流电阻大大电流绕组里这个增量可以到百分之几十。COMSOL里的“多匝线圈”功能可以直接算这个。它会自动把线圈域内由正交电流产生的损耗积分出来其中包含了集肤效应和邻近效应的影响。需要注意的是线圈在离散化时每个网格单元里的电流密度是均匀的如果网格太粗集肤效应根本捕捉不到。导线截面上的网格至少要能分辨出一个集肤深度铜在50Hz下的集肤深度大约9.3mm对常见电磁线来说其实不严苛但如果做高频变压器这个约束就会变得非常致命。绕组的绝缘层对传热影响很大但对电磁损耗计算几乎没影响。所以电磁模型里绕组可以做均质化处理也就是把铜和绝缘统一成一个等效域。热模型中再按“铜-绝缘-油道”分层建模。两个模型之间的损耗传递靠的是体积损耗密度分布而不是总损耗数值。2.3 铁心损耗与杂散损耗的处理技巧铁心损耗COMSOL有两种算法。一种是直接用Bertotti模型把磁滞损耗、涡流损耗、异常损耗三项分开用频率和磁通密度幅值表达另一种是直接从频域磁场结算出的磁场中用铁损密度公式计算涡流损耗分量再叠加一个磁滞损耗修正。工程上我更推荐用Bertotti公式因为它的参数可以从硅钢片厂家给的单位损耗曲线上拟合出来。厂家数据表里通常给的是1.0T、1.5T、1.7T下每公斤损耗的瓦数把这些值代入模型校准系数算出来的总铁损和实测值能控制在百分之几以内。注意硅钢片方向性问题如果用的是取向硅钢导磁性和损耗都沿轧制方向最好模型里就要设置各向异性的B-H曲线。杂散损耗是最容易被忽略、也最难算准的一块。它在结构件里的分布极不均匀往往集中在螺栓、夹件尖角等局部位置。网格要足够细才能看到这些局部热点。三维模型如果网格不够杂散损耗算出来普遍偏低。我的建议是如果变压器容量不大、油箱离绕组远直接忽略杂散损耗用1.05~1.1的系数放大铜耗铁耗来包络如果必须算单独建一个包含油箱和结构件的三维模型和主模型分开算。2.4 损耗从电磁场到温度场的传递方式损耗传递有两种做法。第一种最简单粗暴把电磁场算出的总损耗均匀地加载到热模型的绕组域和铁心域上。这个做法误差大因为绕组热点位置往往由局部的损耗密度决定均匀化之后热点温度会偏低。第二种是保留损耗的空间分布。在COMSOL里可以在电磁模型的“派生值”节点中把单位体积损耗密度导出成插值函数或者直接用同一几何模型在传热接口里引用“mf.Qh”变量。如果你的模型是电磁和热一起求解的那直接用mf.Qh作为热源即可这个变量就是电磁损耗密度它会随温度动态变化实现真正的双向耦合。如果两个模型分开做我建议先在电磁模型里把损耗密度导出为数据集再用“体插值”映射到热模型网格上。导出前一定要确认电磁模型的网格密度足够损耗密度在绕组截面上不是均匀的邻近效应会让靠近油道的导线损耗远大于绕组中心。网格太粗会抹平这个梯度温度场自然也就偏了。3. 流固耦合温度场建模这是最花时间、也最容易出问题的一步3.1 流体场怎么选层流还是湍流变压器油在绕组油道里的流动状态直接决定了对流换热强度流动模型的选错会让温度场整个漂移。先判断自然对流的强度。变压器油在绕组油道里的流速通常很低参考瑞利数判断油道特征尺寸在毫米到厘米级别温差只有几十度油的瑞利数在10⁴到10⁶之间处于层流向湍流过渡的边缘。工程上我的经验是配电变压器、油道窄的情况按层流算温升结果偏保守大型变压器、油流速度高的情况要按湍流算。如果选湍流模型比较常用的是k-ε或k-ω。但我提醒一句低雷诺数湍流模型在自然对流中经常不收敛用标准k-ε默认壁函数时边界层网格的y值必须落在合理区间第一层网格厚度要按壁面剪切速度估算。这个要求对流固耦合模型来说很疼——热计算要求边界层网格细湍流壁函数又要求第一层网格不能太细两个条件容易打架。我的建议是先用层流把整体趋势算出来再局部细化边界层网格用k-ω SST模型做对比验证差别在5%以内就不用纠结。3.2 传热接口的核心设置传热接口选择“流体传热ht”它会自动和流体接口耦合在流体域中同时解传热方程和流动方程。固体域中则只有导热不参与流体计算。自然对流的密度变化COMSOL里有两种处理方式完整变密度或Boussinesq近似。Boussinesq近似的思路是除了浮力项其他地方的密度都当作常数浮力项用体积膨胀系数乘以温差来表达。这个近似在温差小几十度以内、膨胀系数小油和水的膨胀系数都满足的情况下精度很高而且数值上稳定得多强烈推荐使用。如果做强迫油循环温差更低Boussinesq也完全够用。油的物性参数务必设置成温度相关的尤其动力粘度。变压器油的粘度随温度变化非常剧烈——从20℃到80℃粘度可以下降一个数量级。如果把粘度设成常数算出来的流速和换热系数都会失真。导热系数和比热容随温度变化小一些但我也建议从资料中把曲线拟合法输进去。温度相关的材料属性会让求解难度上升但这是精度的代价省不得。3.3 边界条件设置的几个关键点流场方面变压器油箱是封闭的所有油箱壁和固体表面都是无滑移壁面条件。绕组和铁心表面也是壁面油只能在它们之间绕行。如果没有外部油泵入口出口都没有流场完全靠浮力驱动这其实是个封闭腔体内的自然对流问题。温度边界方面油箱外壁和空气之间的换热需要用“对流热通量”边界条件换热系数取值很讲究。室内变压器自然冷却时油箱壁表面换热系数大约在5~10W/(m²·K)加上辐射等效换热后可以取到10~15W/(m²·K)。如果做户外变压器还要考虑风速的影响和日照辐射。这个边界条件对顶层油温的影响非常直接建议大家有条件的话用试验数据反推一下换热系数的取值别照搬书上的数。绕组内部结构还要处理一个绝缘层的问题。实际导线外面包着绝缘纸绝缘纸的导热系数只有0.2~0.3W/(m·K)比铜的400W/(m·K)差了几千倍热量从导线到油的主要阻力就在这里。建模时可以在导线外面加一个薄层或者在热传导域中给绕组域设置一个等效的、较低的平均导热系数。用“薄层”边界条件更准确但会显著增加网格量用等效导热系数则更省事。我的做法是电磁模型用等效均质域热模型用分层结构——铜导线定义为实体绝缘纸定义为薄层边界。3.4 耦合方式与求解顺序安排流固耦合温度场的多物理场节点COMSOL里有现成的“非等温流”耦合它会自动把流体传热和流体流动连起来。要确认三件事流体接口是否开启了浮力项、传热接口是否勾选了“流体域”中的对流项、两个接口是否共用了同一个网格。求解顺序上我强烈建议分两步走不要上来就全耦合。第一步固定一个初始温度比如环境温度20℃先把流场算到稳定然后把流场结果作为初值开启非等温流耦合继续求解。这样比直接全耦合稳定得多收敛时间也快。如果还是发散可以把瞬态求解器切出来给定一个很短的时间步让温度场慢慢演化到稳态相当于给非线性迭代一个温暖的起点。4. 网格设计与求解器调试成败都在细节里4.1 边界层网格是温度场精度的生命线流固耦合温度场最关键的网格就是边界层。流体和固体交界处的温度梯度最大如果边界层网格层数不够、第一层厚度太大换热量就会被严重低估。一般的经验是贴近壁面的第一层网格厚度要满足流动边界层的要求同时热边界层需要8~12层网格才能较好捕捉温度梯度。COMSOL的边界层网格功能可以直接在壁面上自动铺。设置时重点关注“厚度调节因子”和“层数”两个参数。厚度调节因子控制第一层网格的厚度与相邻网格的比值取0.1~0.2比较常见层数取8~12层。绕组油道宽度才几毫米油道内部至少要保证4~6层边界层网格否则油道内的热阻会被网格误差吃掉一大块。铁心内部的网格密度反而可以放宽因为铁心是均质发热温度梯度平缓很多。绕组导线内部因为损耗密度不均匀网格适当加密但不必像边界层那样苛刻。4.2 网格单元的细节考虑绕组的几何模型很麻烦导线是圆形的、油道是窄长条、绝缘纸很薄尺寸跨度太大。直接用自由四面体网格经常会出现过度扭曲的单元。我的经验是先对几何进行简化处理把导线截面从圆改成等效正方形油道宽度稍微放大一点只要流动面积大致不变、表面积大致不变温度场的偏差完全可以接受。几何清理工作看着琐碎但能帮你省掉后面跟畸形单元死磕的功夫。很多新手会在同一个域里使用四面体和六面体混合网格这没有不行但要注意交接面的网格必须一致否则耦合变量无法传递。我更倾向于整体用四面体加边界层网格虽然单元数多一些但生成稳定、迭代不闹脾气。只有在绕组结构比较规则的模型中才用扫掠网格生成六面体质量高、收敛快。4.3 求解器配置分离式还是全耦合全耦合求解器把所有物理场的方程放一起解理论上最准确但非线性强时非常容易发散。我的默认做法是电磁场部分用频域直接求解器流固耦合温度场部分用分离式求解器每一步先解流场再解温度场交替迭代。分离式求解器占内存小而且某个物理场发散时更容易定位问题。如果模型规模大单元数超过50万建议打开“几何多重网格”预处理器或者至少使用PARDISO直接求解器内存足够时PARDISO比迭代求解器稳定得多。我遇到过不少情况模型怎么调都不收敛换一个求解器设置就过了所以求解器配置本身就是排查收敛问题的第一站。收敛判据方面默认的容差经常不够。自然对流问题存在多个时间尺度速度场的收敛慢温度场的收敛相对快如果使用统一的容差速度场还没收敛就停了。建议把速度的容差设置得更严格一些同时观察“全局评估”里的油流平均速度它不随迭代变化了才算真正收敛。5. 常见问题排查与工程校验5.1 收敛不了怎么办流固耦合温度场不收敛90%的锅在网络质量和初始值上。先检查网格是否出现负体积或高倾斜度单元再确认初始值是否合理比如把油流初速设为零温度初值设为环境温度一步步往目标工况逼近。如果发散点在油道区域把油道里边界层层数减少、厚度调大先跑通再加密。如果发散点在浮力项先把密度变化关掉让流场和温度场解耦确认网格没毛病后再开启浮力。这个“先解耦、再耦合”的思路是排查耦合发散最快的方法。5.2 结果合理性怎么判断拿到温度场后不要急着出图先做三个校验。第一总损耗和输入功率对得上吗把所有域上的损耗密度积分和电磁模型算出的总损耗比较误差应在几个百分点内。第二油流方向对吗自然对流应该是底部进油、顶部出油绕组周围的流线应该呈上升趋势如果算出来油往下流那肯定是浮力方向设置反了。第三热点位置合理吗绕组热点温度一般在绕组高度70%~80%的位置、靠近油道一侧如果热点在绕组正中央且温度奇高检查一下绝缘层是不是被忽略了。还可以做一个简单的手算校验用总损耗除以油箱表面积得到平均热流密度再估算油箱壁的平均温升和仿真结果对比相差超过20%就要回头查边界条件。5.3 几个容易忽略的工程细节变压器油的发热和冷却过程中油箱壁实际还有辐射换热。虽然油浸式变压器温度不高辐射换热量占总散热量的比例在室内自然通风条件下可以达到20%~30%。在“表面对表面辐射”计算量过大的情况下可以采用等效换热系数把辐射折进去或者用简化的“表面发射率”边界条件。绕组中的绝缘油道如果过小油流可能不是顺畅上升而是出现局部的回流和滞止区。这些位置的散热很差是潜在的热点。仿真的时候看到局部温度异常升高先检查是不是油道几何太窄导致流动受阻再考虑是不是网格问题。油箱内壁和绕组之间的油箱屏蔽、隔板等结构如果只是为了引油流在模型里可以用薄壁面代替不需要画成实体。几何越简单网格越漂亮求解越快。但简化时务必保证油流主要通道的截面形状和尺寸不变这是流场准确的前提。5.4 参数敏感性速查我把自己调试这类模型常用的几个参数和调整方向整理在下面方便大家排查问题时对照。参数对结果的影响调整方向油粘度粘度增大油流减慢换热变弱温升升高确认使用温度相关粘度曲线绝缘纸导热系数直接影响绕组内到油道的温度降分层建模时取0.2~0.3 W/(m·K)油箱外壁换热系数影响顶层油温和油箱壁温按室内/室外环境选5~15 W/(m²·K)边界层网格厚度过厚会低估换热温升偏高厚度调节因子取0.1~0.2损耗密度空间分布分布越均匀热点温度越容易被低估保留电磁损耗密度分布不用均质化浮力项关闭后只有导热温度严重偏高开启Boussinesq近似最后再分享一个我个人的体会这套模型最耗时间的不是求解本身而是电磁场和热流场的网格互相迁就。电磁场想要网格细捕捉集肤效应流场想要边界层网格细捕捉换热两套需求叠在一起网格量很容易失控。我的做法是分成两个模型电磁模型用二维或准三维做把损耗密度精确算出来热流耦合模型用专门优化的三维网格做专注散热路径。两个模型之间用插值函数传递损耗分布。这样既能保证损耗精度又不用牺牲流固耦合网格的质量计算时间也大幅减少。等到需要做参数扫描、找热点随负载变化规律时再考虑把两个模型合并成完整的三维双向耦合模型。建议从分步方案入手跑通了再逐步增加耦合的复杂度这条路走下来最顺。
基于COMSOL的油浸式变压器多物理场流固耦合温升仿真建模
做过变压器设计或温升校核的人应该都有体会光算电气参数不难难的是把“电-磁-热-流”这四件事放在同一个模型里讲清楚。这个COMSOL仿真模型的核心就是围绕变压器的电磁场、损耗计算和流固耦合温度场展开把绕组发热、铁心发热、油流散热、油箱壁自然对流换热这一整条路径完整还原出来。这篇文章我会从物理本质、建模策略、实操设置到收敛调试把我做这类模型的完整思路和踩过的坑一次讲透适合正在做油浸式变压器温升仿真、或者想入多物理场耦合的新手参考。先说明一点变压器流固耦合温度场本质上是一个“热源在固体里冷却介质在流动两者通过壁面交换热量”的强耦合问题。发热来自电磁损耗损耗大小又和温度有关铜的电阻率随温度升高而增大油的粘度也随温度变化反过来影响对流强度。这一连串互相咬合的物理过程正好是COMSOL最擅长的场景。但“能做”和“做得好”之间隔着大量细节下面逐个展开。1. 先把物理过程拆清楚这个模型到底在解什么1.1 变压器内部的热源与散热路径油浸式变压器的热量来源主要有三块。第一是绕组铜耗也叫负载损耗电流流过导线电阻产生和电流平方成正比第二是铁心损耗包括磁滞损耗和涡流损耗由交变磁通在硅钢片中产生铁心损耗基本只和电压、频率有关和负载关系不大第三是杂散损耗漏磁通在油箱壁、夹件、压板等结构件中感生出涡流造成的附加损耗大容量变压器里这块不可忽略。散热路径则是这样走的绕组和铁心先通过固体导热把热量传到表面再由表面的绝缘油通过对流把热带走。油受热后密度变小自然向上浮形成从底部上升、顶部冷却后下沉的自然对流循环热量最终通过油箱壁传给外部空气。绕组内部还有绝缘纸和油道热量从导线穿过绝缘层再到油道的路径上温度梯度很大热点温度通常在绕组中部偏上的位置而不是最顶端。搞清楚这三条热源路径和一条散热路径模型里需要哪些物理场已经清楚了电磁场算损耗固体传热算绕组和铁心里的温度分布流体场算油的流动传热接口负责把固体和流体两侧的温度和热流耦合起来。1.2 为什么不能只算温度场必须做多物理场耦合有人可能会问我直接把损耗加载到固体上给油侧一个等效对流换热系数不也能算出温度吗对很多工程校核确实是这么干的用经验公式估一个表面换热系数算出来的热点温度在初步设计阶段够用。但问题在于换热系数本身是温度场和流场共同作用的结果不是一个固定值。油的自然对流强度取决于壁面温度和油温之间的温差、油的物性参数以及流道几何。绕组不同位置的温度不一样局部的换热系数也不一样负载变了发热量变了油的流速就跟着变换热能力也在变。用一个常数换热系数去近似等于把最关键的散热机制给钉死了算出来的温度分布自然不可信。要做准确的温升预测就必须把流场解出来让热量传输路径自己“跑”出来这就是流固耦合存在的意义。而且损耗本身也是温度的函数。绕组铜的电阻率随温度升高而增大温度越高损耗越大损耗大了温度再升高形成正反馈。不把电磁场和温度场耦合起来损耗就只能是常温下的定值和实际情况偏差很大。这也是很多新手做变压器温升仿真时忽略的一点。2. 电磁场与损耗计算热源的精度决定温度场的精度2.1 磁场接口选择与几何简化COMSOL里做变压器电磁场AC/DC模块下的“磁场mf”接口是首选频域求解。用频率域的好处是直接拿到稳态下的损耗不用等瞬态波形稳定下来。模型的几何可以简化铁心用硅钢片叠片后的等效实体绕组用多匝线圈域代替不再逐匝画导线。但有一个关键点我想提醒一定要评估“二维轴对称、二维平面、三维全模型”三种方案怎么选。常见的油浸式配电变压器是三相三柱结构但绕组总体呈轴对称分布很多工况下可以用二维轴对称模型计算量小一个数量级。如果你主要关心温度场而电磁场只是用来算损耗的那完全可以用二维电磁场算出损耗分布再映射到三维热流耦合模型里。我自己做单相变压器或者三相变压器单相分析时基本都是这个路子实测下来温度和温升试验结果对得上。三维电磁场模型主要用于需要精确考虑杂散损耗、漏磁在结构件中分布的场景或者油箱壁涡流损耗占比很高的大容量变压器。三维模型的网格量很容易上百万加上流固耦合常规电脑基本跑不动所以一开始就要想清楚取舍。2.2 绕组损耗的准确计算方法绕组损耗分直流损耗和交流损耗。直流损耗好算就是I²RR按工作温度下的电阻率修正。交流损耗则要麻烦一些因为交变漏磁场会在导线内部引起集肤效应和邻近效应让有效电阻比直流电阻大大电流绕组里这个增量可以到百分之几十。COMSOL里的“多匝线圈”功能可以直接算这个。它会自动把线圈域内由正交电流产生的损耗积分出来其中包含了集肤效应和邻近效应的影响。需要注意的是线圈在离散化时每个网格单元里的电流密度是均匀的如果网格太粗集肤效应根本捕捉不到。导线截面上的网格至少要能分辨出一个集肤深度铜在50Hz下的集肤深度大约9.3mm对常见电磁线来说其实不严苛但如果做高频变压器这个约束就会变得非常致命。绕组的绝缘层对传热影响很大但对电磁损耗计算几乎没影响。所以电磁模型里绕组可以做均质化处理也就是把铜和绝缘统一成一个等效域。热模型中再按“铜-绝缘-油道”分层建模。两个模型之间的损耗传递靠的是体积损耗密度分布而不是总损耗数值。2.3 铁心损耗与杂散损耗的处理技巧铁心损耗COMSOL有两种算法。一种是直接用Bertotti模型把磁滞损耗、涡流损耗、异常损耗三项分开用频率和磁通密度幅值表达另一种是直接从频域磁场结算出的磁场中用铁损密度公式计算涡流损耗分量再叠加一个磁滞损耗修正。工程上我更推荐用Bertotti公式因为它的参数可以从硅钢片厂家给的单位损耗曲线上拟合出来。厂家数据表里通常给的是1.0T、1.5T、1.7T下每公斤损耗的瓦数把这些值代入模型校准系数算出来的总铁损和实测值能控制在百分之几以内。注意硅钢片方向性问题如果用的是取向硅钢导磁性和损耗都沿轧制方向最好模型里就要设置各向异性的B-H曲线。杂散损耗是最容易被忽略、也最难算准的一块。它在结构件里的分布极不均匀往往集中在螺栓、夹件尖角等局部位置。网格要足够细才能看到这些局部热点。三维模型如果网格不够杂散损耗算出来普遍偏低。我的建议是如果变压器容量不大、油箱离绕组远直接忽略杂散损耗用1.05~1.1的系数放大铜耗铁耗来包络如果必须算单独建一个包含油箱和结构件的三维模型和主模型分开算。2.4 损耗从电磁场到温度场的传递方式损耗传递有两种做法。第一种最简单粗暴把电磁场算出的总损耗均匀地加载到热模型的绕组域和铁心域上。这个做法误差大因为绕组热点位置往往由局部的损耗密度决定均匀化之后热点温度会偏低。第二种是保留损耗的空间分布。在COMSOL里可以在电磁模型的“派生值”节点中把单位体积损耗密度导出成插值函数或者直接用同一几何模型在传热接口里引用“mf.Qh”变量。如果你的模型是电磁和热一起求解的那直接用mf.Qh作为热源即可这个变量就是电磁损耗密度它会随温度动态变化实现真正的双向耦合。如果两个模型分开做我建议先在电磁模型里把损耗密度导出为数据集再用“体插值”映射到热模型网格上。导出前一定要确认电磁模型的网格密度足够损耗密度在绕组截面上不是均匀的邻近效应会让靠近油道的导线损耗远大于绕组中心。网格太粗会抹平这个梯度温度场自然也就偏了。3. 流固耦合温度场建模这是最花时间、也最容易出问题的一步3.1 流体场怎么选层流还是湍流变压器油在绕组油道里的流动状态直接决定了对流换热强度流动模型的选错会让温度场整个漂移。先判断自然对流的强度。变压器油在绕组油道里的流速通常很低参考瑞利数判断油道特征尺寸在毫米到厘米级别温差只有几十度油的瑞利数在10⁴到10⁶之间处于层流向湍流过渡的边缘。工程上我的经验是配电变压器、油道窄的情况按层流算温升结果偏保守大型变压器、油流速度高的情况要按湍流算。如果选湍流模型比较常用的是k-ε或k-ω。但我提醒一句低雷诺数湍流模型在自然对流中经常不收敛用标准k-ε默认壁函数时边界层网格的y值必须落在合理区间第一层网格厚度要按壁面剪切速度估算。这个要求对流固耦合模型来说很疼——热计算要求边界层网格细湍流壁函数又要求第一层网格不能太细两个条件容易打架。我的建议是先用层流把整体趋势算出来再局部细化边界层网格用k-ω SST模型做对比验证差别在5%以内就不用纠结。3.2 传热接口的核心设置传热接口选择“流体传热ht”它会自动和流体接口耦合在流体域中同时解传热方程和流动方程。固体域中则只有导热不参与流体计算。自然对流的密度变化COMSOL里有两种处理方式完整变密度或Boussinesq近似。Boussinesq近似的思路是除了浮力项其他地方的密度都当作常数浮力项用体积膨胀系数乘以温差来表达。这个近似在温差小几十度以内、膨胀系数小油和水的膨胀系数都满足的情况下精度很高而且数值上稳定得多强烈推荐使用。如果做强迫油循环温差更低Boussinesq也完全够用。油的物性参数务必设置成温度相关的尤其动力粘度。变压器油的粘度随温度变化非常剧烈——从20℃到80℃粘度可以下降一个数量级。如果把粘度设成常数算出来的流速和换热系数都会失真。导热系数和比热容随温度变化小一些但我也建议从资料中把曲线拟合法输进去。温度相关的材料属性会让求解难度上升但这是精度的代价省不得。3.3 边界条件设置的几个关键点流场方面变压器油箱是封闭的所有油箱壁和固体表面都是无滑移壁面条件。绕组和铁心表面也是壁面油只能在它们之间绕行。如果没有外部油泵入口出口都没有流场完全靠浮力驱动这其实是个封闭腔体内的自然对流问题。温度边界方面油箱外壁和空气之间的换热需要用“对流热通量”边界条件换热系数取值很讲究。室内变压器自然冷却时油箱壁表面换热系数大约在5~10W/(m²·K)加上辐射等效换热后可以取到10~15W/(m²·K)。如果做户外变压器还要考虑风速的影响和日照辐射。这个边界条件对顶层油温的影响非常直接建议大家有条件的话用试验数据反推一下换热系数的取值别照搬书上的数。绕组内部结构还要处理一个绝缘层的问题。实际导线外面包着绝缘纸绝缘纸的导热系数只有0.2~0.3W/(m·K)比铜的400W/(m·K)差了几千倍热量从导线到油的主要阻力就在这里。建模时可以在导线外面加一个薄层或者在热传导域中给绕组域设置一个等效的、较低的平均导热系数。用“薄层”边界条件更准确但会显著增加网格量用等效导热系数则更省事。我的做法是电磁模型用等效均质域热模型用分层结构——铜导线定义为实体绝缘纸定义为薄层边界。3.4 耦合方式与求解顺序安排流固耦合温度场的多物理场节点COMSOL里有现成的“非等温流”耦合它会自动把流体传热和流体流动连起来。要确认三件事流体接口是否开启了浮力项、传热接口是否勾选了“流体域”中的对流项、两个接口是否共用了同一个网格。求解顺序上我强烈建议分两步走不要上来就全耦合。第一步固定一个初始温度比如环境温度20℃先把流场算到稳定然后把流场结果作为初值开启非等温流耦合继续求解。这样比直接全耦合稳定得多收敛时间也快。如果还是发散可以把瞬态求解器切出来给定一个很短的时间步让温度场慢慢演化到稳态相当于给非线性迭代一个温暖的起点。4. 网格设计与求解器调试成败都在细节里4.1 边界层网格是温度场精度的生命线流固耦合温度场最关键的网格就是边界层。流体和固体交界处的温度梯度最大如果边界层网格层数不够、第一层厚度太大换热量就会被严重低估。一般的经验是贴近壁面的第一层网格厚度要满足流动边界层的要求同时热边界层需要8~12层网格才能较好捕捉温度梯度。COMSOL的边界层网格功能可以直接在壁面上自动铺。设置时重点关注“厚度调节因子”和“层数”两个参数。厚度调节因子控制第一层网格的厚度与相邻网格的比值取0.1~0.2比较常见层数取8~12层。绕组油道宽度才几毫米油道内部至少要保证4~6层边界层网格否则油道内的热阻会被网格误差吃掉一大块。铁心内部的网格密度反而可以放宽因为铁心是均质发热温度梯度平缓很多。绕组导线内部因为损耗密度不均匀网格适当加密但不必像边界层那样苛刻。4.2 网格单元的细节考虑绕组的几何模型很麻烦导线是圆形的、油道是窄长条、绝缘纸很薄尺寸跨度太大。直接用自由四面体网格经常会出现过度扭曲的单元。我的经验是先对几何进行简化处理把导线截面从圆改成等效正方形油道宽度稍微放大一点只要流动面积大致不变、表面积大致不变温度场的偏差完全可以接受。几何清理工作看着琐碎但能帮你省掉后面跟畸形单元死磕的功夫。很多新手会在同一个域里使用四面体和六面体混合网格这没有不行但要注意交接面的网格必须一致否则耦合变量无法传递。我更倾向于整体用四面体加边界层网格虽然单元数多一些但生成稳定、迭代不闹脾气。只有在绕组结构比较规则的模型中才用扫掠网格生成六面体质量高、收敛快。4.3 求解器配置分离式还是全耦合全耦合求解器把所有物理场的方程放一起解理论上最准确但非线性强时非常容易发散。我的默认做法是电磁场部分用频域直接求解器流固耦合温度场部分用分离式求解器每一步先解流场再解温度场交替迭代。分离式求解器占内存小而且某个物理场发散时更容易定位问题。如果模型规模大单元数超过50万建议打开“几何多重网格”预处理器或者至少使用PARDISO直接求解器内存足够时PARDISO比迭代求解器稳定得多。我遇到过不少情况模型怎么调都不收敛换一个求解器设置就过了所以求解器配置本身就是排查收敛问题的第一站。收敛判据方面默认的容差经常不够。自然对流问题存在多个时间尺度速度场的收敛慢温度场的收敛相对快如果使用统一的容差速度场还没收敛就停了。建议把速度的容差设置得更严格一些同时观察“全局评估”里的油流平均速度它不随迭代变化了才算真正收敛。5. 常见问题排查与工程校验5.1 收敛不了怎么办流固耦合温度场不收敛90%的锅在网络质量和初始值上。先检查网格是否出现负体积或高倾斜度单元再确认初始值是否合理比如把油流初速设为零温度初值设为环境温度一步步往目标工况逼近。如果发散点在油道区域把油道里边界层层数减少、厚度调大先跑通再加密。如果发散点在浮力项先把密度变化关掉让流场和温度场解耦确认网格没毛病后再开启浮力。这个“先解耦、再耦合”的思路是排查耦合发散最快的方法。5.2 结果合理性怎么判断拿到温度场后不要急着出图先做三个校验。第一总损耗和输入功率对得上吗把所有域上的损耗密度积分和电磁模型算出的总损耗比较误差应在几个百分点内。第二油流方向对吗自然对流应该是底部进油、顶部出油绕组周围的流线应该呈上升趋势如果算出来油往下流那肯定是浮力方向设置反了。第三热点位置合理吗绕组热点温度一般在绕组高度70%~80%的位置、靠近油道一侧如果热点在绕组正中央且温度奇高检查一下绝缘层是不是被忽略了。还可以做一个简单的手算校验用总损耗除以油箱表面积得到平均热流密度再估算油箱壁的平均温升和仿真结果对比相差超过20%就要回头查边界条件。5.3 几个容易忽略的工程细节变压器油的发热和冷却过程中油箱壁实际还有辐射换热。虽然油浸式变压器温度不高辐射换热量占总散热量的比例在室内自然通风条件下可以达到20%~30%。在“表面对表面辐射”计算量过大的情况下可以采用等效换热系数把辐射折进去或者用简化的“表面发射率”边界条件。绕组中的绝缘油道如果过小油流可能不是顺畅上升而是出现局部的回流和滞止区。这些位置的散热很差是潜在的热点。仿真的时候看到局部温度异常升高先检查是不是油道几何太窄导致流动受阻再考虑是不是网格问题。油箱内壁和绕组之间的油箱屏蔽、隔板等结构如果只是为了引油流在模型里可以用薄壁面代替不需要画成实体。几何越简单网格越漂亮求解越快。但简化时务必保证油流主要通道的截面形状和尺寸不变这是流场准确的前提。5.4 参数敏感性速查我把自己调试这类模型常用的几个参数和调整方向整理在下面方便大家排查问题时对照。参数对结果的影响调整方向油粘度粘度增大油流减慢换热变弱温升升高确认使用温度相关粘度曲线绝缘纸导热系数直接影响绕组内到油道的温度降分层建模时取0.2~0.3 W/(m·K)油箱外壁换热系数影响顶层油温和油箱壁温按室内/室外环境选5~15 W/(m²·K)边界层网格厚度过厚会低估换热温升偏高厚度调节因子取0.1~0.2损耗密度空间分布分布越均匀热点温度越容易被低估保留电磁损耗密度分布不用均质化浮力项关闭后只有导热温度严重偏高开启Boussinesq近似最后再分享一个我个人的体会这套模型最耗时间的不是求解本身而是电磁场和热流场的网格互相迁就。电磁场想要网格细捕捉集肤效应流场想要边界层网格细捕捉换热两套需求叠在一起网格量很容易失控。我的做法是分成两个模型电磁模型用二维或准三维做把损耗密度精确算出来热流耦合模型用专门优化的三维网格做专注散热路径。两个模型之间用插值函数传递损耗分布。这样既能保证损耗精度又不用牺牲流固耦合网格的质量计算时间也大幅减少。等到需要做参数扫描、找热点随负载变化规律时再考虑把两个模型合并成完整的三维双向耦合模型。建议从分步方案入手跑通了再逐步增加耦合的复杂度这条路走下来最顺。