高压直流电源电-固-热耦合仿真：用COMSOL一次算清发热、温度与热应力

做高压直流电源的工程师应该都撞过这类事样机调试完全正常一到高温满载或者长时间老化发现内部功率电阻的引线端子歪了固定绝缘柱的塑料件变形了甚至灌封的环氧树脂边缘出现细微裂纹。第一次遇到这类故障我第一反应是装配工艺问题或者零件批次质量问题等把温度场和结构场的测试数据摆在一起才反应过来——这其实是电、热、应力三场耦合的结果。这篇文章就讲清楚一件事怎么用COMSOL把高压直流电源内部的发热、温度分布、材料热膨胀和应力形变放进同一个模型里求解。所谓电-固-热耦合本质就是把电流产生的焦耳热、固体内部的热传导、以及温度引起的热应变与应力串成一个完整的因果链在软件里一次算完。这套思路对做高压电源结构设计的工程师、搞热设计的同行、以及刚接触多物理场仿真的学生都有直接的参考价值。下面按我自己做这类项目的习惯来写从问题拆解、模型搭建、参数设置一路讲到求解、后处理和常见坑。大部分内容围绕一个具体案例一台输出5kV级、内部含高压采样电阻与铜汇流排的直流电源模块。这个模型不算复杂但把高压电力电子设备里最典型的发热、散热、形变场景都覆盖到了。1. 先搞明白高压直流电源里的发热与形变是怎么耦合的1.1 一个真实故障现场的复盘几年前我接过一个返修案例一台高压电源模块在老化测试进行到10小时左右操作员发现内部绝缘支撑柱肉眼可见地歪了拆开后看到支撑柱底部的环氧粘接层有放射状裂纹。当时先用热电偶测了表面温度最高点出现在高压采样电阻附近约85℃绝缘柱附近约70℃。说实话70℃对环氧材料来说并不算极端但问题在于这个温度是“局部”的而且整个装配体从室温升到70℃的过程中不同材料的膨胀量差得非常多。那时候我才真正意识到故障根因不是某一个材料不行而是铜母线、陶瓷绝缘件和铝外壳的热膨胀系数不匹配在升温过程中产生了热应力叠加在绝缘柱的安装预紧力上最后把粘接层拉裂了。这个判断如果只靠温度测试是做不出来的必须把电场产热、热量传播和结构变形放在同一个逻辑链条里分析。1.2 三种物理场之间的因果链电-固-热耦合听起来复杂实际上因果关系非常清晰。第一步是电到热电流流过导体和电阻元件产生焦耳热单位体积发热率是qJ·EJ²/σ其中J是电流密度σ是电导率。这个公式就是整个发热分析的核心任何导电域只要有电流流过就会按这个关系产生热源。第二步是热到力温度升高后材料发生热膨胀应变与温差成正比ε_thα·ΔT其中α是热膨胀系数ΔT是相对某个“应力自由参考温度”的温升。这个热应变在约束条件下就转化成了热应力也就是我们常说的“热胀冷缩被卡住”。需要注意第三步是热对电的反向影响几乎所有金属的电导率都随温度升高而下降铜在20℃到70℃的温度范围内电阻率会上升大约19%。这意味着同一个电流在高温下会产生更大的焦耳热形成一个正反馈回路。对多数工程场景这个反馈会造成5%~20%的发热量差异不是小数字。至于结构变形反过来影响电学和热学通常只是微小几何变化工程上可以忽略。所以真正需要“双向耦合”的是电场和温度场结构场直接吃温度场的计算结果就够了。1.3 高压直流电源里到底哪些地方在发热高压直流电源内部的发热源不只是一个。实际产品中主要发热位置大概有这么几类功率半导体器件的导通损耗和开关损耗这个往往是整机最大热源几十瓦到几百瓦很正常高频变压器和电感的铜损加磁芯损耗高压采样电阻、均压电阻和泄放电阻这类元件按U²/R计算损耗例如5kV电压加在250kΩ的电阻上发热就是100W母线排、接线端子和螺栓连接处的接触电阻电流密度大的地方容易形成局部热点。不同发热源在仿真里的处理方式不一样。半导体损耗通常作为体热源直接施加到器件芯片域高压电阻和母线这类导电元件可以直接放进电流场里算焦耳热。建模前先把这个清单列清楚想明白每个发热源是用“真实电场算”还是用“等效热源近似”后面就不会乱。2. 建模方案怎么取舍物理场、耦合方式与几何简化2.1 物理场怎么选才不绕远路COMSOL里跟这个需求相关的物理场有三个基础选项。电流场用Electric Currents这个接口适合处理直流或工频下的传导电流问题它求解的是电势方程能得到电流密度分布和焦耳热分布。如果是几十千赫兹以上的高频开关波形严格来说应该用电磁场或电路耦合但对于热分析而言直接算交流电磁场代价太高工程上更常用的做法是把开关损耗折算成热源电流场只管母线、电阻这类元件的传导和发热。温度场用Heat Transfer in Solids如果机箱内部有风机、需要考虑空气的对流换热细节可以额外加一个流体流动接口做共轭传热代价是网格量和计算时间成倍增加。对于第一版评估模型我建议用等效对流换热系数替代空气域后面需要精确的时候再把空气域加回来。结构场用Solid Mechanics标准线弹性就够用除非温度很高导致材料进入塑性或者要评估蠕变才需要额外启用塑性或蠕变本构。耦合部分直接使用COMSOL的多物理场节点。Joule Heating节点会自动把电流场产生的焦耳热作为热源耦合进传热方程Thermal Expansion节点会自动把温度场作为体载荷施加到固体力学方程。这两个节点是电-固-热耦合的标配不需要手动去写耦合项。2.2 单向耦合够不够什么时候必须双向这个问题我几乎每次都会被问到。先说结论电场和温度场建议做双向耦合温度场到结构场做单向耦合就够了。原因前面提过金属电阻率对温度很敏感铜温升50℃电阻率约增加19%这会让发热量明显偏高。如果不把这个反馈算进去热点温度经常会低估3到10℃对于关注温升裕量的设计来说误差偏大。实际操作中COMSOL允许你把三个物理场全部“全耦合”在一个研究里求解也可以拆成两步。如果模型规模大、内存紧张我习惯先用分离式求解器让软件按电场→温度场→结构场的顺序迭代。计算发现多数情况下迭代十几步就能收敛。结构场之所以可以单向是因为结构变形不会显著改变电流路径和散热路径除非你研究的是大变形柔性结构那另当别论。2.3 几何简化原则不是越细越好很多新手拿到一个模块的CAD图恨不得把每颗螺丝、每个倒角都导进COMSOL结果网格量巨大、求解半天还不收敛。我自己的原则是去掉对电流路径、散热路径和应力路径影响小于1%的特征重点保留母线排截面、电阻器件、绝缘支撑件、散热底板这些热力传导的主干。具体到这个案例几何可以简化成一个组合体一块200mm×120mm×6mm的铝合金散热底板上面通过四根氧化铝陶瓷支撑柱托起一块3mm厚的FR4绝缘板绝缘板上竖直安装一只圆柱形高压陶瓷电阻电阻两侧各接一段铜母线排一进一回。陶瓷支撑柱的螺纹孔直接简化成圆柱螺纹细节不建模。产品如果关于中面对称可以建1/2模型甚至1/4模型切分面上设置对称边界条件计算量能降到原来的四分之一甚至八分之一。3. 材料参数仿真精度的一多半由这里决定3.1 电参数电阻率必须写成温度的函数材料参数是整个仿真里最容易“看起来没问题、实际全错”的地方。最典型的问题就是拿常温电阻率套在一直在发热的元件上。电流场里一定要把电导率定义成温度的函数COMSOL里可以在材料节点用变量表达式直接写例如铜的电阻率可以写成ρ(T)1.75e-8[Ω·m]×(10.0039[1/K]×(T-293.15[K]))。这里0.0039就是铜的电阻温度系数每升高1K电阻率增加0.39%。电阻体的处理要单独说一下。我案例里用的是一个表现损耗的等效体积电阻250kΩ阻值、直径8mm、长度40mm的圆柱体体积电阻率通过公式换算ρeR×A/L250000×(π×0.004²)/0.04算出来大约314Ω·m。这个数值远高于金属但低于绝缘体在COMSOL里直接作为该域的电导率即可它能真实反映电阻体的压降和发热。3.2 热参数与热膨胀系数一张表说明核心差异做热-结构耦合时热膨胀系数CTE的匹配性是决定应力大小的根本原因。常见的几种材料我列了个表方便直接取用作为初值材料电导率σ(S/m)热导率k(W/(m·K))密度(kg/m³)比热容c(J/(kg·K))CTE(ppm/K)弹性模量E(GPa)备注纯铜5.98e7400893038517110屈服强度约70~200MPa视热处理铝合金60632.5e7200270090023.469屈服强度约110~150MPa氧化铝陶瓷96%绝缘2838008507320脆性看第一主应力/抗弯强度FR4绝缘0.3519001200面内13/面外5025各向异性Tg约130~140℃环氧灌封胶绝缘0.61600110035~505温度接近Tg后性能骤降从这张表能直接看出来为什么温度一高界面就容易出问题。铜的CTE是17ppm陶瓷只有7ppm铝是23.4ppm。同样升温50℃铜膨胀0.085%陶瓷膨胀0.035%如果两者刚性连接又没有释放间隙界面上的剪切应力会非常可观。FR4除了面内和面外CTE差异大还有一个玻璃化转变温度Tg的问题温度超过Tg后模量陡降热胀冷缩行为也变这些都必须记在心里。3.3 材料界面的处理方式COMSOL里几何装配有两种模式Form Union和Form Assembly。对于初学者我建议先直接用Form Union把各个零件当成一个连续体里的不同域材料属性自动按域分配界面默认是理想连续电流、热量和位移都连续传递。这样建模最快也不会出现接触设置的问题。如果想要模拟更真实的界面比如功率电阻与FR4板之间的薄层粘接胶、母线排与散热器之间的接触热阻那就用Form Assembly加配对边界再设置接触热阻或者薄层这样能让温差更贴近实测。但代价是引入接触一致性条件和接触热阻参数求解非线性程度上升。第一版模型我一贯的建议是先用Form Union把趋势算出来再在重点界面上细化。4. 实操一步步搭出电-固-热耦合模型4.1 几何建立一个可复现的示例结构下面按一个能直接复现的示例来讲。整体尺寸参考前面说的铝合散热底板200×120×6mm四根直径10mm、高15mm的Al₂O₃陶瓷圆柱支撑柱分布在底板四角距边缘约15mmFR4绝缘板200×120×3mm放在支撑柱上绝缘板中央安装一只圆柱电阻直径8mm、长40mm电阻体两端各连接一条尺寸为150×10×3mm的L型铜母线。建立几何时在COMSOL里直接用长方体、圆柱和拉伸操作就能完成也可以导入三维CAD再布尔运算。建议把每个零件的定义域单独用“选择”分组命名比如“busbar_hv”、“busbar_gnd”、“resistor”、“fr4_board”、“ceramic_post”、“al_base”后面设置材料、边界条件和网格时按组选择会方便很多。4.2 物理场边界条件一步一步配置电流场这步最关键。高压母线排末端的一个端面设置“终端”边界条件类型选“电压”电压值填5000V。低压回流母线排的电路端面设置“接地”。这样电流从高压母线流进电阻体再通过低压母线流回形成完整回路。除这两个边界外其余外边界默认都是电绝缘不需要额外操作。要注意的是COMSOL里电流场必须有一个明确的“参考零势面”否则电势没有定义求解必然失败。传热场的边界条件相对直观。所有暴露在外界空气的表面设置对流热通量对流换热系数按工况填。自然散热时竖直表面取5~8W/(m²·K)水平面朝上取10~15水平面朝下取3~6初次估算统一用8即可。环境温度设293.15K。电阻体表面如果温度预期超过100℃建议额外加“表面辐射”边界发射率填0.85否则温度会高估10~20℃。结构场的边界要小心。把底板下面四个安装孔所在面设置“固定约束”模拟整机安装到机箱上的状态。其他表面不加约束让模型自由热膨胀。这里有个容易犯的错把所有螺栓孔都固定死等于认为结构是完全刚性约束实际螺栓连接有预紧力也有滑移热应力会被人为放大。更稳妥的做法是用“弹簧基础”或“边界约束”里的“指定位移”限定某个方向避免完全静定约束。多物理场节点只需要确保两件事Joule Heating节点把电磁热源加到传热方程Thermal Expansion节点把温度体载荷加到结构方程。Thermal Expansion里有一个“热膨胀参考温度”或者叫应力自由参考温度这个值非常关键。它表示结构在多少温度下是完全没有热应力的。如果你的装配体是在60℃固化粘接然后冷却到25℃那参考温度应该设60℃而不是默认的20℃。这个细节直接影响应力基线和应力符号。4.3 网格划分不同的物理对网格的需求不一样网格是好结果的一半但没必要追求全模型加密。我的策略是分区处理。母线排和电阻体这类长条形结构用扫掠网格沿长度方向划分横截面分个6到10层FR4板厚度方向至少3层面向内用自由三角形后拉伸或者直接扫掠陶瓷支撑柱和底板用自由四面体。关键界面比如电阻体与母线接合处、母线折角处局部细化到0.2~0.5mm。有个细节值得记一下如果算完发现某个位置温度梯度或者应力梯度特别陡先用局部细化网格重算一遍看结果变化大不大。如果细化后应力明显跳变这说明原网格没收敛可能还有几何奇异点而不是材料真的有那么大的应力。整体上这个规模的模型网格数控制在10万到30万单元、自由度在100万以内普通工作站几分钟到十几分钟就能算完。4.4 求解器设置不收敛的时候怎么处理研究类型先选“稳态”。默认的全耦合求解器在电-固-热耦合模型里经常遭遇收敛困难尤其是电阻率随温度剧烈变化的时候。我习惯在求解器设置里切换到“分离”模式把求解步骤按物理场拆成三步先算电流场和温度场再算结构场。分离迭代的次数默认是15次不够就改成50收敛公差保持默认即可。如果还是发散发散一个好用的技巧是启用“辅助扫描”从低电压开始。先让U500V算一遍收敛后再继续扫描U1000V、2000V、5000V相当于给求解器一组递进的初始值。这个做法在工程上特别好用代价只多算几遍但基本能消除大部分非线性收敛问题。如果你要的是瞬态结果比如开机热冲击、断电冷却循环就选“瞬态”研究时间区间设0到3600秒步长自动即可注意所有边界条件要设为随时间变化的函数。4.5 结果后处理从云图到设计结论后处理我习惯按四个维度来看。第一步看温度场找到最高温度的位置确认热点在不在预期的元件上。第二步看电流密度特别是母线折角、螺栓压接点附近的电流拥挤情况电流密度异常高的位置通常就是局部过热的隐患点。第三步看位移把变形用“位移幅值”云图显示同时把变形缩放系数调大比如放大100倍能直观看到哪个零件往哪个方向歪。第四步看应力金属件看von Mises应力陶瓷和环氧这类脆性材料看第一主应力因为它们的失效准则是拉应力主导。从Derived Values里可以直接提取体积最大值、表面最大值和时间平均温度整理成表格后跟材料允许值对照判断。我的参考判据是PCB和环氧件最高温度至少低于T_g 20℃以上陶瓷件第一主应力低于其抗弯强度除以安全系数3铝合金底板von Mises应力低于屈服强度。如果仿真结果和这些判据冲突就要在结构或者散热设计上做改动这就是这个模型的核心价值。5. 常见问题与排查实录5.1 电流算出来了但发热功率对不上这个坑我见过太多次。最常见的原因是电路回路没闭合。电流场必须有一个“终端”和一个“接地”否则仿真软件无从建立电流通路结果经常是电流为零或者发热区在绝缘材料里异常产热。第二个常见原因是把电阻体的电导率设置成了默认的绝缘值或者填数值时单位弄混把电阻率当电导率填。第三个原因是几何单位没换算清楚COMSOL里几何用毫米、材料参数用国际单位时物理场会自动做单位换算但如果你在表达式里手写数值很容易漏掉单位标注。我的排查顺序是先在后处理里看电流密度的流线确认电流路径是不是符合预期再在电阻域上对焦耳热做体积积分把这个功率和手算的U²/R对比偏差超过10%就一定有问题。这个步骤不要省因为热源功率一旦错了后面的温度和应力全都不对。5.2 温度场高到离谱或者干脆不收敛温度爆表的原因分两类。一类是物理上确实可能高温比如散热条件太差、热源功率过大这时要回头检查热源估算和h的取值是不是偏保守了。另一类是数值问题比如电阻率随温度正反馈太强温度越高电阻越大、发热越多稳态解可能根本不存在。遇到这种情况就把“辅助扫描”从更低电压开始或者改瞬态求解让温度从室温逐步上升看看系统是否能稳定在一个平衡点。还有一个容易忽略的因素是辐射。当表面温度超过120℃时辐射散热量占比接近自然对流的30%到50%。如果不加辐射边界热平衡温度会偏高。对于金属外壳内部的元器件因为没有风道、对流很弱h可能只有2~3W/(m²·K)这种情况下辐射反而是主散热路径。先别急着怀疑求解器先确认自己对散热路径的假设是不是贴近现实。5.3 应力云图里出现“漂亮但离谱”的应力集中这是我最常提醒同行的一点。热应力分析里任何几何尖角、锐边、直角过渡都会产生数学上的应力奇异网格越细应力值越高仿佛永远算不准。在真实零件里这些奇异点的应力会被材料塑性、圆角等自然缓解。COMSOL里处理这个问题的思路不是无限加密网格而是辨识奇异区域并且在评估时把距离奇异点一定范围以外的应力作为判断依据或者直接在几何里增加小圆角。另一个非常常见的过约束问题。初学者喜欢把底板的整个底面都固定住实际安装只是几个螺栓点两者计算出来的热应力可能相差数倍。正确的做法是用有限的固定区域模拟安装约束其他区域自由膨胀。如果结构装配关系复杂还可以在关键螺栓处设置弹簧边界让结构有一定的滑移自由度。仿真结果是要指导设计的应力如果虚高反而会误导设计加无谓的加强结构。5.4 仿真温度和实测偏差很大从哪里查起模型和实际样机对不上时从三个方向查。第一查热源整机里是不是还有别的发热源没加进模型比如变压器和电感损耗、半导体开关损耗高压电源内部往往是多热源叠加。第二查界面电阻与散热板、器件与绝缘板之间有没有真实的粘接胶层或导热垫片这层东西的热阻常常成为主要温差来源。用一个简单办法验证——做个瞬态温升测试看起温斜率反算整体的等效热阻和热容然后跟仿真对比。第三查边界机箱内部的空气温度其实高于环境温度仿真直接用环境温度作为对流参考温度会高估散热效果。我在实际项目中吃过一次亏仿真预测电阻表面温度88℃实测却到了103℃。查了半天发现我设的h8是机箱外表面散热但电阻位于封闭箱体内部内部空气几乎静止真正有效的散热路径是通过绝缘板和底板传导对流换热系数远低于外面的自然对流。调整边界条件后误差从15℃降到了2℃以内。5.5 陶瓷绝缘件裂纹的排查故事再分享一个完整的排查案例。某次仿真显示陶瓷绝缘柱最大第一主应力约120MPa低于陶瓷抗弯强度300MPa理论上安全但试制样机偏偏在绝缘柱根部出现裂纹。重新检查后发现问题出在热膨胀参考温度上。这批绝缘柱和FR4板是用环氧胶在60℃固化的固化完成后冷却到室温时整个装配体已经带上了残余热应力。我在模型里把热膨胀参考温度从20℃改成60℃后应力分布完全改变第一主应力峰值从120MPa升到了接近200MPa说明裂纹不是运行温度造成的而是装配固化冷却过程叠加运行温度之后的结果。从那以后我在每个热-结构耦合模型里都会专门确认这个“应力自由参考温度”。6. 写在最后的几点实操体会做这类电-固-热耦合仿真做了几年我最大的体会是模型精度往往不由求解器能力决定而由你对真实产品里热源分布、散热路径和约束方式的理解决定。COMSOL只是把物理规律表达出来你给它错误的边界条件它就给你一个精确的错误结果。花半天时间去搞清楚产品里到底哪几个器件发热、热量通过什么路径排出、安装约束是刚性还是弹性比调一星期求解器都管用。如果只是临时评估风险用单向耦合快速出温度趋势就行如果要出设计报告或者做故障分析务必做双向电-热耦合并把电阻率写成温度的函数。最后分享一个小技巧在COMSOL里把变形结果放大30到50倍导出一段动画拿去找结构工程师讨论往往比几十页的数据报告更能说明问题。热应力这种东西光看数值不容易形成直觉一旦看到结构“扭起来”方案取舍很快就清晰了。