前两年做电力设备绝缘结构优化时,遇到一个很蹊跷的问题:某环氧浇注绝缘子在短时耐压测试中,电压还没达到材料标称击穿场强,局部位置就出现了明显的碳化通道。当时第一反应是材料有内部缺陷,但反复切片和电镜观察后,发现通道起始位置更像热源而非缺陷。后来用Comsol Multiphysics搭了一个热击穿模型,才彻底想明白——这根本不是传统意义上的电击穿,而是焦耳热在绝缘材料内部不断累积,最终把材料的绝缘性能热失控掉的过程。这篇就围绕Comsol热击穿模型展开:热击穿是怎么发生的、为什么高压场景下很容易触发、模型里需要哪些关键设置、稳态和瞬态怎么配合着判断临界电压,以及我踩过的一些坑。适合做高压绝缘设计、电缆附件开发、功率模块封装的朋友参考,也适合刚接触热电多物理场仿真的研究生用来理解整个建模思路。1. 先搞清楚:热击穿盯上的不是电场极限,而是热平衡1.1 三种击穿里,热击穿是最容易被误判的绝缘材料的失效,工程上通常分三类。第一类是纯电击穿,电子在强电场下被直接加速,发生碰撞电离,击穿时间在纳秒到微秒级。这种和材料的本征击穿场强直接相关,做实验的时候加压速度要非常快,不然测量到的其实已经不是纯电击穿了。第二类是局部放电引发的侵蚀击穿。绝缘内部或表面存在气隙、气泡,在电场下反复微放电,慢慢把材料蚀穿。这个时间尺度可以到月甚至年,常见于运行中的电缆、电机绕组。第三类就是热击穿,本质不是场强超高,而是泄漏电流产生的焦耳热和散热之间失去平衡。热击穿的时间跨度很大,可能是几秒,也可能在长期运行后缓慢发展。它的判断标准不是电场超过了多少MV/m,而是单位时间内产生的热能是否大于散掉的热能。也就是说,这是一个热平衡问题,不只是单纯的场强问题。这种差别在实验里很容易骗人。耐压测试里看到材料失效,第一反应都是怀疑场强设计不够、局部几何结构造成电场集中。但有些情况下,加压过程比较慢,泄漏电流一直在给材料加热,等温度升高到某个程度,电导率大幅上升,电流再次上升,温度就收不住了——最后材料烧穿、碳化,外观上看起来和电击穿留下的通道几乎没有区别。我第一次遇到的时候就被它骗过一次。1.2 为什么高压场景特别容易触发热击穿热击穿的驱动源是焦耳热,而焦耳热功率密度在电场下等于电导率乘以电场强度的平方,也就是 q σE²。这里有个容易被忽视的点:虽然绝缘材料电导率极低,但在高场强下E非常大,平方之后产生的体积热不可低估。举个例子,某环氧材料常温电导率约1e-14 S/m,如果在局部位置场强达到10 MV/m,体热功率密度大约 q 1e-14 × (1e7)² 1 W/m³。单独看这个数不大,但考虑到绝缘层厚度小、温升引起的电导率呈指数上升,外加散热路径差,温度就会慢慢爬升。另外一个重要原因是高压绝缘结构里总存在电场不均匀的地方,比如电极毛刺、绝缘层内部气隙、电缆附件半导电层切口。这些位置的场强远高于平均场强,热源高度集中。在COMSOL模型里最常见的情况,最后温度最高点就是这些高场强区域。本质上,高压设备绝缘设计不只要控制最大场强,还要保证在最苛刻工况下,焦耳热有足够的散热通道,让热平衡点存在。这就是热击穿模型要解决的工程问题。注意:热击穿这个词在某些行业叫热老化或热不稳定性,模型物理完全一致,读文献时留意一下术语就行。2. 失稳的数学本质:电流-热耦合的正反馈怎么走向失控2.1 电导率对温度的指数依赖热击穿模型的第一个核心输入,就是电导率随温度的变化关系。绝缘材料的电导率不是常数,而是强烈依赖温度。物理上,这来自离子电导的主导机制——温度升高,离子迁移率变大,电导率上升。工程上常用两种表达式。Arrhenius形式:σ(T) σ0 × exp(-Ea / (kB × T))其中Ea是活化能,kB是玻尔兹曼常数,T是绝对温度。这种形式物理意义清晰,适合与文献值比对。另一种是简化指数形式:σ(T) σ0 × exp(α × (T - T0))在温度变化范围不大的情况下,两者基本等价,参数α、T0可以通过拟合Arrhenius表达式得到。对初学者,第二种形式更好调参,因为它直接反映每升高1K,电导率涨多少。这里有个工程经验:很多仿真文献直接给常温电导率,不给你活化能。如果找不到Ea,可以用典型高分子绝缘材料0.5~1.0 eV这个区间来估计。Ea取0.8 eV时,温度从20℃升到80℃,电导率大约上升两个数量级以上——这个放大效应足够触发热击穿。所以选错活化能,仿真结果能差出一个量级,这不是危言耸听。2.2 热平衡方程与失稳判据在COMSOL里,热击穿模型实际上是两个物理场的耦合。电流场(电流守恒加欧姆定律):∇·J 0, J -σ(T)∇V固体传热:ρCp ∂T/∂t - ∇·(k∇T) σ|∇V|²右边那一项就是焦耳热源。稳态情况下,温度对时间的导数为零,方程变成一个带源项的非线性热传导方程。为什么非线性?因为源项里的σ(T)又依赖T,而且是强烈指数依赖。数学上,这类非线性稳态方程不是对所有电压都有解——当电压超过某个临界值Vc,别说稳定的稳态温度分布不存在,甚至任何稳态解都不存在。可以用一个集中的热平衡示意来理解。把绝缘材料看作一个热容很小的集中热节点,发热功率P(T)正比于σ(T)×V²,散热功率Q(T)主要指对流散热和向周围的热传导。随着T升高,P(T)会指数上升,而Q(T)最多线性上升。两条曲线的关系有三种情况:电压较低时,曲线P(T)和Q(T)在低温处有个交点,这个交点是稳定平衡点。加扰动后温度会自动回到这个点。电压提高后,曲线往上移,交点逐渐靠近切点。达到临界电压时,两条曲线相切;超过临界电压,二者没有交点,任何初始温度都会持续升温直到材料破坏。这个切点就是热击穿临界点,对应的电压就是临界电压。这就是整个模型的物理判据。3. 从零搭一个热击穿模型:几何、材料与耦合3.1 几何简化与拓扑仿真几何越简单越好,尤其是第一版模型。我常用二维轴对称平板或圆柱结构模拟绝缘试样:厚度d 2 mm(沿轴向)半径R 10 mm上下各留一层电极区域,厚度0.5 mm,材料设为铜电极要不要建出来?主流做法是直接在绝缘层上下端面加终端边界条件,不建电极几何。如果建电极,好处是热量可以通过电极散走,更接近实际;坏处是多了一个高导电高导热域,网格和求解压力会大一些。我的建议是:第一版不建电极,用边界条件代替;后续做热传导校验时再把电极加上。几何拓扑上特别注意一点:COMSOL转CAD内核时对某些导入几何的狭长面、微小圆角支持并不好,在网格生成阶段容易报不支持的拓扑错误。从CAD导入复杂结构前,建议用删除小特征、修复几何等操作把微小圆角去掉,或者干脆在COMSOL内部用基本几何体搭一个简化模型。绝缘击穿模型精度对几何细节不敏感,对材料参数和边界更敏感。3.2 材料参数:热击穿模型的三件套热击穿模型需要三组核心参数,缺一不可:电参数:相对介电常数、电导率σ(T)热参数:导热系数k、密度ρ、比热容Cp失效判据:材料临界温度,比如环氧的玻璃化转变温度或分解温度一组典型环氧绝缘料的参数如下:参数数值说明常温电导率 σ01e-14 S/m20℃下实测值活化能 Ea0.8 eV决定电导率温度灵敏度导热系数 k0.2 W/(m·K)聚合物典型值密度 ρ1200 kg/m³常规环氧树脂密度比热容 Cp1500 J/(kg·K)典型聚合物比热相对介电常数 εr4.0直流建模时影响瞬态响应这些参数未必需要特别精确,因为热击穿模型的工程价值更多在于判断临界电压在哪个量级、温度最高点出现在哪,而不是给出某个极小点处的精确温度值。前提是电导率随温度的规律要选对,绝对值可以后面用实测数据校准。3.3 物理场选择与多物理场耦合COMSOL里搭热击穿模型,推荐用以下组合:电流场:Electric Currents(ec)固体传热:Solid Heat Transfer(ht)耦合方式:用多物理场节点下的Electromagnetic Heating,或直接在传热方程里手动添加热源项 Q ec.Qrh手动加源项和用内置耦合节点的差别:内置节点会自动处理焦耳热和电导率温度依赖的双向耦合,省事且不容易漏项;手动加源项目前仅适合想完全控制方程的用户。我给新手的第一建议永远是直接用内置节点。边界条件设置:上端面:电势V 0,接地下端面:电势V V0,在全局定义里设一个参数,方便扫描侧面:电绝缘热边界:上下端面设置为固定温度T0 293.15 K,模拟与电极、外部导体紧密接触时散热能力强;侧面设置对流热通量,对流换热系数h取5~10 W/(m²·K),外部温度T0这里有个容易出错的地方:如果把上下端面也设成绝热,那绝缘层几乎是闷烧状态,临界电压会低得离谱。实际高压绝缘结构大多通过金属附件散热,不能一律绝热。至于对流换热系数,我之前见过有人默认填100 W/(m²·K),那相当于强风冷条件,会明显高估散热能力,让热击穿临界电压偏高。自然对流在室内环境一般就5~10,除非你明确知道自己要模拟风冷或油冷工况。4. 稳态扫描的取舍:算不出解不一定是模型错4.1 辅助扫描:从低压到高压爬上去热击穿模型最典型的非线性行为,是稳态解随电压变化的曲线:低压时,温度随电压平缓上升;接近临界电压时,温度增长越来越快;超过临界电压后,稳态解直接消失。因此,不要上来就在目标高电压下求解。正确做法是用COMSOL的辅助扫描,也就是参数化扫描加继续法:定义参数V0,初始值设为1 kV。使用参数化扫描,扫描范围从1 kV到某个最大值,比如40 kV,步长1 kV。求解器设置为辅助扫描,它会在每个参数步以上一步解为初值继续迭代。这个以上一步解为初值非常关键。高非线性问题对初值极其敏感,直接让求解器从零初值去算40 kV下的稳态,即使存在稳态解,Newton迭代也很容易发散到看似无解。一步步扫描时,你会看到温度—电压曲线在后段明显上弯。此时就需要加密电压步长,把临界点附近的曲线形状算出来。比如在35 kV附近算不出时,改为34.0、34.2、34.4这样的小步长继续逼近,能更准地锁定临界电压区间。4.2 不收敛意味着什么:物理无解还是数值问题这是整个热击穿建模里最容易吵起来的问题:同一模型,有人说稳态能收敛,有人说不能。原因多半不在模型本身,而在你从哪个方向逼近这个电压。稳态求解器报错时,先做三个排查:检查是否从低电压逐步扫描上来的。直接高压起算产生无解是很常见的。检查最后一步解的残差。如果残差在发散,大概率是物理上已无稳定稳态;如果残差持续下降但下降到一半卡住,可能是网格或数值问题。尝试更严格的相对容差(1e-6)和更小的步长。若仍然不收敛,再判断为物理临界。还有一类报错来自初值计算失败,COMSOL会在日志里提示找不到一致的初始值或迭代初始残差过大。这种情况先做一步无电流纯热传导的求解,把温度场算个大概,再用它作为电流耦合求解的初始条件,往往能绕过去。提示:热击穿模型的稳态不收敛,有时是好消息——它说明你已经逼近或超过了临界电压。不要急着调网格、换求解器,先复核电压-温度曲线是否已经出现明显上弯。在COMSOL中,如果稳态求解器不收敛,还可以尝试把非线性方法从自动改成恒定牛顿,并打开阻尼因子,或者使用分离式求解步骤:先求解电流场再求解热场,交替迭代。分离式在部分强耦合问题上比全耦合稳健,但收敛速度慢一些。我的习惯是:全耦合优先,分离式作为备选。5. 瞬态分析补位:温度失控的时间与击穿路径5.1 什么时候必须转瞬态稳态模型只能告诉你临界电压在哪和临界时温度场怎么分布,但它回答不了两个工程问题。第一,如果实际工作电压略低于临界电压,材料会不会在某个长时间尺度上仍然击穿?这是因为稳态解的稳定性还取决于扰动大小和热积累速度。某些情况下,稳态解存在,但稳定裕度很低,一个外部热扰动或负载波动就可能把系统推到失稳区。这时需要瞬态分析看响应过程。第二,若电压高于临界电压,温度从初始状态爬升到破坏温度需要多久?这个时间决定了保护装置能否及时动作、耐压测试该按什么时间窗口判定结果。瞬态设置的要点:研究类型选时域显式(Time Dependent)。初始值用低压稳态解,比如1 kV下的稳态温度场,然后突然把电压阶跃到目标电压。时间序列设为对数增长,例如 1e-4、1e-3、1e-2、...、100 s,便于捕捉早期快速变化和后期缓慢增长。若算到后期出现温度急剧上翘,说明热击穿正在发生,记录到达临界温度的时刻。这里有个容易犯的错误:瞬态时间步长过大,导致温度越过临界温度后继续上涨很小,曲线看起来平缓。问题在于自适应步长把一段快速上升过程漏掉了。解决方法是显式输出热源功率密度或最大温度随时间的变化曲线,用最大温度—时间曲线的斜率判断是否出现指数式上升,而不是只看最终状态。5.2 从瞬态结果判断击穿时间与击穿位置瞬态求解完成后,后处理不要只盯温度云图。我通常同时看三个量:最大温度随时间变化曲线(T_max vs t)最大电流密度随时间变化曲线某个监测点(比如电极边缘高场强处)的电导率随时间变化曲线三个量如果同时出现指数式攀升,而且曲线尾部斜率越来越大,基本可以确定热击穿进程已不可逆。取材料临界温度(比如环氧170℃开始明显分解)对应的时间,就是工程意义上的击穿时间。击穿位置由于高电场与热源正反馈的共同作用,通常不会分布在全域,而是集中在场强最高且散热最差的点。后处理时可以做一个温度高于临界温度的等值面提取,直接看出击穿通道的起始点。很多时候,这个点不会是几何正中心,而在电极边缘附近——因为电场集中效应在边缘最明显。这也是为什么高压绝缘设计里,边缘的圆角半径、毛刺打磨永远是质量控制重点。6. 结果提炼:从温度云图到临界电压预测6.1 电流-电压特性曲线:读取临界电压的最直接方法做完参数化扫描后,用全局计算提取每个电压稳态解下的总泄漏电流。方法是派生值—全局计算,表达式用边界上的法向电流密度积分,或者直接在终端边界上读取电流值。然后以电压为横轴、电流为纵轴作图。温度稳定的热击穿模型,低压段I-V曲线基本是直线,因为电导率几乎没变;电压继续升高,曲线开始明显弯曲;到临界附近,电流急剧升高,甚至呈垂直线上升。通常让电压步长足够密后,曲线的近似垂直段就是临界电压区间。判断临界电压还有一个辅助手段:把每个稳态解下的最大温度也提取出来,画出电压-最大温度曲线。临界电压附近这条曲线同样会急速上升。两条曲线交叉验证,比单独看I-V曲线可靠,因为I-V曲线有时受边界条件设置影响较大。6.2 与理论解析解互验对于平板几何和一维等温假设,热击穿临界电压有一个近似解析式可以用于量级校验:Vc ≈ sqrt(8 × k_th × T0² / (σ0 × Ea/kB))其中k_th是导热系数,Ea是活化能,kB是玻尔兹曼常数,T0是环境温度。代入前面示例参数,k_th 0.2 W/(m·K),T0 300 K,σ0 1e-14 S/m,Ea/kB约9280 K,得到的Vc在几十千伏量级。这个解析解虽然粗糙,它假设材料内温度、场强均匀,实际会有梯度,但用来校核仿真量级非常有用——如果COMSOL算出来和解析值差了几个数量级,那多半是边界条件、单位或参数输入有问题,而不是物理差异。我在实际项目里一般把仿真临界电压与解析量级控制在2倍以内,作为模型可信度的一条快速检验线。超过2倍就要回头查参数和边界了。6.3 后处理里值得注意的细节电导率用对数坐标显示:绝缘材料电导率跨越几个数量级,线性色标会把低电导率区域显示成一片底色,无法体现局部分布。观察热通量矢量:如果散热设计合理,热通量应该从高温热点指向低温边界;如果发现热通量方向异常,比如指向高场区域,检查网格或边界设置。用数据留样显示失效区域时,把阈值设成材料临界温度,而不是某个相对温度。电力工程师看到局部温度超过200℃比温度升高了150K更能直接对应失效判据。7. 建模中反复折磨人的细节:踩坑实录7.1 电导率表达式引起的单位与数量级问题热击穿模型里电导率跨度的巨大是数值问题的主要来源。常温下环氧可能只有1e-14 S/m,温度升高到160℃时可能到1e-7 S/m,跨越了7个数量级。COMSOL默认求解器对这类尺度的系统会有一定压力。几个实际经验:电导率不要直接写1e-14 S/m就完了。建议把σ(T)用解析函数定义,方便后续校准和批量修改。COMSOL里的解析函数可以直接引用状态变量T。如果电导率下限过低,比如低于1e-18 S/m,矩阵条件数会恶化,甚至报奇异矩阵。可以给电导率设置一个下限,比如 max(σ(T), 1e-16 S/m),避免数值发散。物理上,极低温或极低电场时电荷输运机制可能已经转变,这个下限也有物理合理性。有些用户从文献里拿到的是μS/cm,直接填进COMSOL就错得离谱。1 μS/cm 1e-4 S/m,换算时务必仔细。7.2 残留特征与CAD导入的坑热击穿本身不需要移动网格,但很多人是从结构仿真、流体仿真转来做电气仿真的,旧文件里残留着移动网格、变形几何等特征,转过来后求解器总会报出现非物理变形之类的错误。热击穿模型用不上这些。如果是从旧文件改造,建议新建研究,删掉不必要的特征。CAD导入时,转换为CAD内核时不支持的拓扑这类报错多发生在带圆角、过渡面的3D几何上。解决办法有两个:在CAD软件里导出前就做几何清理;或者在COMSOL里用删除细小特征和修复节点处理。网格划分前最好先用对象编辑检查一下最小尺寸和最大边比例,不然到求解阶段才发现网格质量差,返工成本高。7.3 排查迭代未收敛的通用套路热击穿模型求解报错,总是非线性不收敛,但原因可以有十几种。我的排查顺序固定如下:先看日志,定位哪一步、哪个物理场报错。检查该物理场的残差是否在某个数量级反复震荡——震荡多半是网格太粗或时间步长过大。缩小电压步长,或瞬态时间步长,重新扫。换求解器:全耦合改分离式;直接求解器PARDISO改MUMPS或GMRES。把初始条件设为上一个低压稳态解,或零电流纯热解。如果以上全试过仍不收敛,物理上可能真的无解了,承认它,把该点视为临界点。这个套路在COMSOL的很多非线性问题上都能用。很多做电流热耦合的人卡在No convergence几天,其实只是高估了全耦合求解器在强非线性上的鲁棒性,分开迭代、缩小步长两步就解决了。7.4 快速收敛的小技巧:把热分析先跑通我的个人习惯是先跑一个不带电流耦合的纯热源模型,直接指定热源为一个常数值,看温度场分布和量级是否合理;然后再开启电导率耦合,加高电压扫描。这样把传热模型正确性和电流耦合正确性两步分开验证,可以避免两个误差源混在一起时根本没法定位问题。如果你发现纯热源模型的温度就不对,后面再怎么调电导率都白搭。个人这几年做绝缘类仿真的体会是,热击穿模型的难点从来不在物理场设置,而在如何把真实材料的行为用可靠的参数表达出来。尤其电导率的温度依赖曲线,直接决定临界电压数值,有条件的话建议实测,或者至少从同材料文献里引用可靠数据,不要拍脑袋给Ea。模型算完,最好再用不同电压下的I-V特性做一次对比验证,能大幅提高结果可信度。另外一个很实用的小技巧:给客户或评审展示热击穿结果时,不要只贴一张热云图,把电压-最大温度曲线、I-V曲线和热通量矢量图放在一起,会让人很快理解你做了完整的热平衡分析,而不是只跑了个云图。