★ FEATURED ARTICLE 做变压器仿真这几年我最大的感受是电磁场和温度场常常被分成两拨人来做电磁工程师算完损耗交给结构工程师结构工程师再凭经验往热源上一堆最后热点温度差个十几度都不知道该怪谁。后来我把COMSOL的单相变压器电磁场和温度场计算模型整理成一套固定流程从几何导入到热点温度出结果基本两三天能跑完一轮而且数据能和实测温升对得上。这篇就把这套模型的完整思路、参数细节和踩坑过程写出来给正在折腾同类仿真的朋友一个参考。这套模型适合谁如果你手头有一台具体型号的单相变压器比如几kVA到几十kVA的干式变压器想回答绕组最热的地方在哪里铁芯损耗到底能不能把表面温度推到限值以上增加风道后温升能降多少这类问题那这篇就是照着做的。如果你是刚开始接触多物理场耦合想弄清楚损耗怎么样均匀地变成热源也可以从里面的技术细节里捞到不少能直接用的经验。1. 建模前先想清楚单相变压器里热到底从哪来很多新手拿到COMSOL第一件事就是画几何。我反而建议先坐下来把热源盘点一遍。单相变压器的热量来源其实就两类绕组里的铜耗和铁芯里的铁耗。铜耗对应一次侧和二次侧电流在导线电阻上的损耗铁耗则来自交变磁通在硅钢片里产生的磁滞损耗和涡流损耗。这两类损耗的量级、分布位置和随温度的变化规律完全不同如果只是随手找一个总损耗灌进传热模块等于把问题粗暴降维了。1.1 什么情况下才需要做电磁—温度完整耦合不是所有变压器问题都需要同时算电磁场和温度场。如果你的目标只是看某个负载率下的稳定温升那可以不做强耦合按电磁算损耗→温度算分布的单向流程走就行。但如果要研究短路工况、频繁过载或者超大电流密度下的绕组热累积温度对铜电阻率的影响就藏不住了这时候电阻率随温度升高导致铜耗进一步上升形成正反馈必须做双向迭代。我最近处理过一台5kVA单相干式变压器项目额定电压220V/110V额定电流大约22.7A和45.5A自然冷却。这种常规容量变压器的绕组温升范围通常在50K到65K之间铜电阻率随温度的变化约0.4%/K整体铜耗变化也就几个百分点所以单次单向耦合已经足够精确。真正会出问题的场景是温度逼近绝缘耐热等级的极限比如B级绝缘允许130°C绕组热点超过120°C之后每升高1°C带来的风险不是线性的这时仿真里那百分之几的损耗误差就可能转化成是否要改设计的判断误差。1.2 单相变压器天然适合二维建模但要先分清心式和壳式单相变压器主要有心式和壳式两种结构。心式结构通常两个铁芯柱一次、二次绕组套在同一柱上磁场分布相对对称壳式结构则像把铁芯包在绕组外围磁路更短漏磁路径复杂。如果是心式结构沿铁芯柱截面的磁场基本保持一致用二维平面模型就能抓住主要物理过程计算量小、网格调试快。壳式结构因为磁路存在三维转折局部漏磁和铁轭处的涡流分布很难在二维里还原除非做大量工程近似否则建议直接上三维模型哪怕只建四分之一对称模型。我做心式单相变压器仿真时使用的就是二维平面模型但有一个前提绕组沿铁芯柱高度方向的端部效应要靠等效长度修正。COMSOL默认的二维模型把变压器当作无限长处理实际铁芯柱是有高度的绕组的端部漏磁场和散热面积都与无限长假设存在偏差。我的处理办法是给铁芯柱和绕组域一个等效深度把端部漏感折算到窗口区域的磁场里这样二维结果的三维修正就有了着落。1.3 目标定了再决定单向耦合还是双向耦合在项目开始前我会把仿真目标写清楚是要得到磁密分布云图来评估铁芯饱和程度还是要得到绕组热点温度和温升数值来评估绝缘寿命或者两者都要。目标直接决定耦合复杂度。如果只关心磁路设计比如铁芯截面选多大、匝数怎么分配那么电磁场单独求解就够了。如果只关心散热方案对比比如加风道、改风扇、加散热片的效果那么可以手工估算损耗作为热源不必每次都跑电磁场。但如果要出完整的温升报告而且是给第三方看的那就老老实实把损耗从电磁场里解析出来再加载到温度场因为报告里专业审核人员会追问你这个热源的分布依据是什么靠估算很难经得起推敲。2. 电磁场计算几何、材料、线圈激励三件套电磁场部分是整个模型的基础基础设施没打牢后面温度场算得再精细也是空中楼阁。这一节我按实际建模顺序讲先建几何再定义材料最后设线圈激励和损耗提取每一步都有值得注意的细节。2.1 几何层级的取舍我最终保留的部件清单在COMSOL的AC/DC模块中选择磁场mf接口求解变量是磁矢势A。几何方面我最终保留的部件很少铁芯包括心柱和铁轭、两个绕组域、绕组之间的绝缘筒或气道、以及外围的空气域。空气域范围取变压器外轮廓的3到5倍宽度因为磁场在空气中衰减较慢截断得过近会让漏磁场失真影响漏感进而影响绕组电流分配。绕组内部的每匝导线我一般不会单独建模原因有二一是50Hz工频下铜导线的集肤深度大约9mm工程上常见线径大多不超过3mm涡流效应不显著二是每匝导线单独建模的网格量在二维下还可以接受但三维里就是灾难。用多匝线圈域整体代替即可铜耗按整体电流密度均匀分布计算这一近似在工频范围内造成的误差很小。空气域我一般设成矩形并给外边界一个磁绝缘条件。如果变压器外围有金属油箱或安装支架要单独建模并给材料属性因为金属结构件里的涡流损耗会贡献一部分发热。我踩过一次坑就是把支架忽略结果实测温升比仿真高出一截后来发现是支架在漏磁场里产生了额外的涡流热。2.2 硅钢片材料处理非线性B-H曲线和各向异性电导率铁芯材料是最容易让仿真翻车的地方。默认的材料库里往往只有相对磁导率为常数比如4000的线性材料这在高磁密时会出大问题。单相变压器铁芯正常工作磁密通常在1.5T到1.7T之间这个区间恰好是硅钢片B-H曲线的膝点附近磁导率已经明显下降如果按常数磁导率算铁芯磁阻会严重偏低励磁电流偏小饱和趋势完全体现不出来。正确的做法是把硅钢片的B-H曲线以表格形式导入材料属性曲线来源可以是硅钢片厂家的出厂数据也可以在实测磁化曲线上取点。COMSOL里可以做成插值函数然后材料属性里填B-H非线性关系。还有一个细节硅钢片的电导率不是标量而是各向异性的。顺轧制方向和垂直于轧制方向的电导率不同叠片之间又有绝缘涂层所以垂直于叠片方向的电导率应该设得很小约等于零否则涡流会沿着厚度方向横穿叠片流动算出大得离谱的涡流损耗。2.3 原副边线圈的激励方式为什么我推荐电流源起步线圈激励是电磁场计算中最容易漏水的地方。COMSOL的线圈Coil特征支持电流或电压激励。我强烈建议第一步先用电流激励给原边线圈设定为额定电流副边按负载情况折算到原边电流或者直接在副边加对应负载电流这样最容易收敛也方便先确认磁路是否正常。等电流激励的结果验证无误后再改成电压激励或外部电路耦合。电压激励时COMSOL需要把线圈内部电阻和电感自洽计算出来这实际上是求解一个包含感应电动势的电路方程设置不当很容易出现电流突跳。我以前试过在副边直接加开路电压结果副边感应出来的电压和原边完全不成比例排查了很久才发现是线圈的匝数方向定义反了。原副边绕向相反时线圈特征中的绕组方向必须分别定义否则磁动势互相抵消励磁电流会戏剧性增大铁芯磁密跌到几乎为零。2.4 频域计算的损耗输出铜耗、铁耗分别怎么取铁芯损耗在COMSOL磁场接口里可以通过损耗计算特征来实现原理是铁芯在交变磁场下的磁滞和涡流效应共同产生热损耗。如果手头有厂家提供的铁损曲线可以直接用斯坦梅茨公式拟合P k * f^α * B^β。其中k、α、β是材料系数B是磁通密度幅值f是频率。不同牌号硅钢片的系数差异很大直接从网上抄一套系数用往往会差20%以上最好按厂家的单位损耗曲线拟合。绕组铜耗的提取则简单很多在频域求解完成后对绕组的电流密度积分或者直接读取线圈特征中的电阻损耗变量。COMSOL里这个变量已经帮你算好了铜导体的总焦耳热包括电流的直流电阻损耗和涡流效应引起的附加损耗。如果绕组域用的导体简化模型这个值就是按均匀电流密度算出的直流铜耗如果设了单导体模型它会额外考虑绕组内部涡流的影响。实际项目中我还喜欢把损耗结果单独做一个损耗汇总表列出铁芯总损耗、原边铜耗、副边铜耗、总损耗这几项。计算完先核对一下理论值比如5kVA变压器满载时总损耗大约应该在100W到150W范围如果仿真给出来的总损耗和手算经验值差了一个数量级先不要往下走回头查材料属性和激励设置。3. 损耗密度怎样变成温度场的热源从电磁场拿到损耗后不能直接把总损耗除以变压器体积再均分到整个域上。变压器内部温度分布是否准确取决于热源的空间分布是否合理。铁芯损耗主要集中在磁通最密集的柱体和轭部绕组损耗则集中在电流密度最高的部分这两个区域不一定重合而热点往往出现在损耗密度高且散热条件差的交点附近。3.1 总损耗均摊和按分布加载的差别有多大我做过一个对比测试把同样的总损耗分别按均匀体积热源和按电磁场计算出的分布热源加载到同一台变压器上结果绕组热点温度相差了8K到10K。原因很直接绕组内侧靠近铁芯柱散热条件差如果按均匀热源处理会把这部分本应更热的区域稀释掉热点位置也会漂移。对温升评估这种热极限问题来说8K的偏差已经足以影响设计结论了。因此每当我看到有人为了省事直接把电磁场损耗结果用平均值算完再输进热场的我都会建议改成按分布加载。好在COMSOL支持在传热接口里直接调用电磁场算出的损耗变量用变量耦合的方式把热源密度逐点映射过去不需要手工导出再导入数据。省掉的步骤不多但精度提升明显。3.2 铁芯损耗做体积热源从磁密表达式映射到传热域铁芯损耗的空间分布可以直接用磁密幅值B的幂函数来描述。具体做法是在磁场的频域求解结束后为传热接口添加一个体积热源表达式写成P_Fe(B)其中B由mf.normB这个内置变量给出。斯坦梅茨公式在这种映射下显得特别自然每个网格点的热源密度由该点的磁密决定磁密高的地方发热就大。需要注意频域求解出来的B是一个复矢量而铁芯损耗使用的是磁密的幅值实际传给热源表达式的应该是B的模COMSOL里的mf.normB就是这个变量。使用中还要确认损耗公式中的系数单位斯坦梅茨公式通常给的是W/kg需要乘以硅钢片密度换算成W/m³才能作为体积热源使用。这个单位换算错误是很常见的问题我在早期项目里就因为漏乘密度导致铁芯温升低得离谱。如果手头没有斯坦梅茨系数也可以采用COMSOL电磁接口自带的铁损计算功能直接输出单位体积损耗密度变量。这种情况下同样可以把这个变量直接引用到传热接口原理上没有区别只是数据来源从经验公式变成了求解器内部的损耗模型。3.3 绕组损耗的分层密度处理与气隙等效导热绕组损耗的分布同样有讲究。两个绕组在铁芯柱上通常是一内一外排列内层绕组离铁芯近散热路径被铁芯挡住一半外层绕组虽然离铁芯远但朝外壳方向还有空气对流可以散热。所以两个绕组的铜耗虽然数值上可能接近但温度场里的表现完全不同。在把铜耗映射到传热接口时要分别给两个绕组域各自的热源表达式这样才可以观察内层绕组是否比外层绕组温度高、热点是否出现在内层绕组中段。绕组内部如果有多层结构和气道实际工程中还有气隙内的自然对流换热。二维传热模型如果不建立空气流体域就要把气隙处理成等效导热系数。密闭空气在狭窄空间里的自然对流等效导热系数通常在2到8W/(m·K)之间视气隙宽度和温差而定。我一般建议先按4W/(m·K)左右做初始计算再看气隙两侧温差确定是否修正。强迫风冷条件下这个值可以更高但这时不如直接建流体域用共轭传热来得准。4. 温度场求解与耦合的落地流程电磁场给温度场喂了热源温度场反过来又要影响电磁场如何处理这条反馈回路是整个模型的关键。实际工程中我习惯的策略是先单向、后判定也就是默认一次单向映射够用但保留双向迭代的通道当温度跨度较大时再启动第二轮。4.1 固体传热接口的参数清单与边界条件取值温度场部分使用固体传热ht接口。材料热参数里铜的导热系数约400W/(m·K)硅钢片导热系数存在各向异性叠片平面方向约25W/(m·K)法向约4W/(m·K)这个差异在模型里最好用各向异性导热系数体现否则铁芯沿叠片方向的散热会被算得偏快。空气域在这个模型中不参与流体计算而是作为导热介质出现。空气导热系数只有0.026W/(m·K)基本属于隔热体但变压器外壳和铁芯之间的空气层确实承担着散热的作用只是这种散热主要靠自然对流单纯靠导热算会偏低。工程做法是给外壳外表面设置对流热通量边界条件自然冷却条件下对流换热系数取值5到15W/(m²·K)。温度升高时空气自然对流加强等效换热系数也会上升更精细的做法是用温度相关的经验关联式去更新h值。辐射散热在自然冷却变压器中占的比例不低壳体表面温度60°C时向环境辐射的热量占总散热量的两成以上。所以外壳表面一定要设置表面对环境辐射边界发射率取0.9左右环境温度按最严酷工况40°C设置。忽略辐射是温度场仿真偏低的常见原因在温差较大的场合尤为明显。4.2 研究序列怎么排先频域电磁再稳态传热COMSOL里同一个模型文件可以包含多个研究。我常用的顺序是先建一个频域研究求解50Hz下的磁场在此基础上再建一个稳态研究求解温度场。稳态研究所用的热源来自频域结果通过变量引用实现跨研究的数据传递这样在参数扫描时可以一次性把不同电流下的损耗和温升一起算完。这里有个容易犯的错误在稳态研究里COMSOL默认会重新计算所有物理场如果磁场接口也在稳态研究的物理场列表中就会把频域求解的磁场再用稳态方式算一遍结果自然不对。解决方法是把稳态研究里磁场接口的求解开关关闭只保留传热接口参与计算。这个开关一般在研究的物理场接口部分设置把不需要的接口从研究列表中取消勾选即可。4.3 温度对电阻率和损耗的反作用要不要做第二轮迭代启动第二轮双向耦合的条件很简单如果第一轮温度结果里铜温和初始假设温度的差值超过了20K铜电阻变化超过8%我建议做第二轮电磁计算更新电阻率后再重新算损耗和温度场。实际项目中从第一轮温度结果把绕组平均温度提取出来更新铜电导率和斯坦梅茨公式里的温度修正项重新求解电磁场和温度场第二轮迭代通常就能收敛第三轮变化已经很小。有一种情况必须做双向迭代那就是绕组电流密度高、振动或过载工况。比如短路阻抗试验或短时过流工况绕组温度在几十秒内飙升铜耗随温度急剧上升此时如果还按常温电阻率算铜耗计算结果会显著低估热点温度。反过来正常长期运行工况由于温度最终稳定在60到100°C之间单向耦合已经够用没必要为了严谨而牺牲效率。5. 仿真结果校验从磁密云图到实测温升一套模型跑通之后最怕的就是拿着颜色漂亮的云图当成算完了。我的习惯是先做几个快速判据检查再做详细后处理最后再和实测数据对照。很多问题在早期判据阶段就能暴露出来省掉后面的大量返工。5.1 电磁场是否合理的快速判据拿到频域结果后我第一件事看铁芯磁密云图。正常工作状态下铁芯柱和铁轭的磁密应在1.5T到1.7T之间如果云图显示大片区域超过1.9T甚至接近2.0T说明铁芯截面选择偏小磁路已经饱和励磁电流会非常大温升也一定不达标。第二个判据是励磁电流大小。电流激励模式下如果原边电流设定为额定值那么求解后看铁芯损耗和铜耗估算出总功率再与变压器的额定容量相比对。一台正常设计变压器的空载损耗大概只有额定容量的0.5%到1%如果仿真中铁耗占比达到百分之几就要检查铁芯材料属性或损耗系数的单位是否出了问题。第三个判据是磁力线分布。单相变压器的磁力线应沿铁芯柱一路穿过铁轭形成完整闭合回路不会大量泄漏到外壳或空气中。打开磁通密度线的后处理如果看到窗口区域到处是杂乱无章的磁力线多半是绕组激励方向设置或几何简化出了问题这时温度场算得再漂亮也没意义。5.2 热点温度位置和温升数值的解读温度场求解完成后我的关注点通常不在最高温度是多少而在最高温度在哪里。单相变压器绕组热点通常出现在内层绕组沿高度方向的中间偏上位置因为绕组底部通过固定支架向外散热顶部虽然热空气上升但是端部暴露面积大中段被内外层绝缘材料夹住散热条件最差。铁芯的表面温度通常不是最高温度最高温度在铁芯内部。如果云图显示热源密度大的铁心柱中部出现局部高温区而后处理出的热点位于绕组而非铁芯这是正常现象。比较工艺上的热点是导线绝缘最薄弱处如果仿真热点温度超过了绝缘耐热等级允许值就要评估换绝缘材料等级或改进散热结构。对照温升限值时我习惯导出最高温度减去环境温度后的温升值而不是直接报绝对温度。标准环境温度40°C下如果是B级绝缘允许130°C那么对应的温升限值就是90K减去绕组温度测量中的仪表误差实际设计通常还会留5到8K的裕量。5.3 与实测对照时的误差拆分模型跑完想验证可信度最直接的办法是和温升试验曲线对比。实测时把热电偶埋进绕组中段和铁芯表面在额定负载下连续通电直到温升稳定记录温度随时间的变化。仿真则对应地跑一个瞬态温度场研究把负载损耗和铁耗按实测工况加载进去两条温升曲线放在一起对比。误差来源一般分三类。第一类是材料参数误差硅钢片实际铁损比出厂标称值高3%到5%铜电阻率受纯度影响也有波动。第二类是边界条件误差自然对流换热系数h值的估算误差能达到30%辐射发射率表面处理状态不同也有偏差。第三类是几何简化误差二维模型的端部效应修正不充分会低估绕组端部散热。我见过几次对不上数据的情况最后都在这三个方向里找到了原因。实在对不上的时候还有一个技巧用仿真和实测的温差作为修正量反推等效对流换热系数。把h值作为可调参数让仿真热点温度与实测热电偶温度一致这样校正后的模型在后续结构改进对比中就非常可信了。6. 实际操作中绕不开的那几个坑最后这一节记录几件我在这个模型上真实遇过的怪问题。这些问题单独看都不大但都足够让人卡上一两天写出来供大家参考。6.1 空载试验里感应电压把电流顶上去的怪事有次我把原边线圈设成电流激励I22.7A副边设成开路本来预期副边感应出一个110V左右的电压结果副边电压远低于预期而原边电压高得离谱。查了很久才发现副边线圈虽然设了开路但线圈特征里默认的导体模型仍会保留很小的残余导电通路在数值上相当于短路了几个欧姆求解器把副边按有载状态算了。解决起来很简单激励方式改为电流激励之后把不需要参与电路求解的线圈特征从线圈改成终止状态或者显式设电流为零。这类问题在电压激励和电流激励混合使用时会频繁出现排查时要先检查每个线圈的求解状态和激励数值不要一上来就怀疑几何。6.2 薄层绝缘界面处的温度振荡第一次跑稳态传热时我在绕组铜和绝缘纸的分界面处发现温度场出现了明显的锯齿状振荡而不是平滑过渡。原因很典型绝缘纸厚度只有零点几毫米网格尺寸却是毫米级温度梯度在窄域里变化剧烈线性单元精度不够网格又来不及加密数值解在界面上产生了伪振荡。处理办法有两个。一是给薄层域添加边界层网格在绝缘纸厚度方向至少划分3到5层单元二是把绝缘纸用薄层特征替代不用在几何里单独建立薄层域直接给界面赋予热阻。后一种方法能大幅减少网格量在二维模型里尤其高效。不过在薄层热阻的表达式中要仔细确认热导率除以厚度的单位换算漏乘厚度是出现过很多次的低级错误。6.3 铁损模型参数标定钢厂数据到软件的转换硅钢片厂家提供的铁损数据通常是在标准正弦激励、50Hz、不同磁密下的单位质量损耗曲线单位是W/kg。而COMSOL的铁损模型需要斯坦梅茨系数或者归一化的损耗数据这两者的换算并不是简单的单位缩放。磁滞损耗系数和涡流损耗系数在斯坦梅茨公式里的频率指数和磁密指数不同直接用曲线拟合时很容易得到一组数值上正确但物理上没意义的系数。我的经验是如果手头只有几条离散的损耗曲线点先在计算表中按斯坦梅茨公式做回归注意固定频率为50Hz对磁密指数单独拟合。如果拟合结果里磁密指数偏离1.6到2.0太多就要怀疑数据本身是否包含谐波磁通或直流偏磁成分。这类数据来源问题不好查但一旦查清损耗精度的提升会立刻反映在温度场对比上。6.4 多工况参数化扫描的组织方式做负载率扫描时我常把原边电流设成参数I_load取值范围从50%到150%额定电流。在COMSOL的参数化扫描研究中可以直接把电流激励里的电流值绑定到扫描参数上。但如果扫描参数跨越的电流范围很大不同电流下的铁芯饱和程度差异明显非线性求解器有时会在低电流和高电流之间切换时难以收敛。这时候我会把扫描的初始值策略设成来自上一个参数并限制每一步的载荷增量让求解器沿着电流从小到大的路径走比直接在大电流工况下冷启动稳定得多。另外扫描结束后生成的结果文件会很多我习惯在每组扫描结果后面加上自动生成的汇总报告把总损耗、磁密最大值、热点温度三个核心指标直接打印出来。这样几十个工况的扫描结果一眼看过去就能挑出异常工况不用逐个点开云图翻找。整套模型最近一次交付时评审那边拿着仿真热点温度数据和第三方温升试验报告对比误差控制在了5K以内这个精度对工程判断来说已经足够用了。如果你正在为变压器温升问题头疼照着这个思路把电磁场和温度场串起来跑一遍大概率会在第一次或第二次求解时就得到能指导设计的结果。
阅读完成 · 觉得有帮助？
