企业官网建设流程全解析

这几年做高压器件仿真，被问得最多的问题之一就是：“这个间隙到底会不会打火？” 每次听到这个问题，我第一个想到的工具就是 COMSOL。电击穿仿真在 COMSOL 里并不是一个“一键求解”的黑盒操作，它需要你同时懂物理图像、模块选型、边界条件和求解器脾气。这篇内容就围绕 COMSOL 电击穿仿真，把背后原理、建模思路、实操路线和踩坑记录完整整理出来，给正在做绝缘评估、放电分析或高压结构设计的朋友一份能直接抄作业的参考。

1. 电击穿仿真到底在算什么：从物理图像到软件方程

1.1 击穿不是“电压到某个值就放电”这么简单

电击穿在工程里常见的场景很多：针板电极之间的空气间隙、PCB 上走线间距、绝缘子沿面爬电、变压器油隙、甚至高压开关柜里的母线相间距离。本质上，电击穿是绝缘介质在外加电场下失去绝缘能力的过程，但物理机制并不统一——气体间隙主要是电子碰撞电离和流注发展，固体绝缘则通常是电-热耦合下的热失稳，液体的机理更复杂。

COMSOL 里做电击穿仿真，第一步不是急着点“研究”，而是想清楚你要算的是哪种击穿，因为不同机理对应不同物理场接口：

• 如果只想知道电场分布，看哪里场强超标，那用 AC/DC 模块里的“静电”接口就够了，配合帕邢定律或经验击穿场强做判据。
• 如果要算固体绝缘的热击穿过程，需要“电流”接口（Electric Currents）耦合“固体传热”，让电导率随温度变化，观察焦耳热导致的温度/电流正反馈。
• 如果要模拟气体放电的流注传播，那就得用“等离子体”模块，里面带了电子输运、反应源项和泊松方程，计算代价也大得多。

我见过不少人直接拿等离子体模块去算所有电击穿问题，结果算半天不收敛，其实他需要的只是电场分布。选题之前先把物理图画清楚，能省一半时间。

1.2 COMSOL 里的几个“电击穿接口”怎么选

COMSOL 官方把电磁场仿真的接口挂在 AC/DC 模块下面，和电击穿直接相关的接口和我的习惯选法如下：

当你不确定选哪个时，我的经验是：先从最便宜的模型入手，算出结果后和实测或经验公式对比，心中有数了再升级模型。直接上等离子体模块就像用高射炮打蚊子，又慢又难收敛，没必要。静电模块算出场强分布、再用帕邢曲线查击穿电压，是很多工程评估的常规动作，速度快、结论清晰，也方便汇报。

1.3 击穿判据：帕邢定律、流注判据和限值评估

拿到电场分布之后，怎么判定“会不会击穿”？最粗略的是看最大场强是不是超过了材料/气体的击穿场强。空气在标准大气压下的均匀场击穿场强大约是 3 kV/mm，但这只是经验值——非均匀场（针板电极）下击穿电压会明显低，不能简单拿 3 kV/mm 去套。

更讲究一点的做法是帕邢定律：对于气体间隙，击穿电压是气压 p 和间隙距离 d 乘积 p·d 的函数，曲线呈 U 形，左边是真空高场强区，右边是高压高距离区。在很多 SF6、空气、氮气场景里，我习惯把帕邢曲线数值查出来放进 COMSOL 的全局表达式里，用后处理判断“给定电压下若满足 pd 条件下的击穿电压低于外加电压→击穿风险高”。

流注判据是更物理的方法：沿着电场线积分有效电离系数 α，当积分超过大约 18~20（经典流注判据），认为流注可以自持发展，击穿即将发生。这个判据在科研论文里很常见，做工程仿真时也可以在 COMSOL 里后处理实现。你需要先把电离系数 α 和电场强度的函数关系写出来（用全局常微分/表达式），然后做个路径积分。我第一次实现这个判据的时候走了不少弯路，后来发现用“派生值→线积分”就能做，不需要额外写代码。

2. 建模前的准备工作：几何、材料与网格

2.1 几何结构怎么画：不是越精细越好

电击穿仿真的几何结构通常不复杂——针板间隙就是一根细针加一块平板，球隙就是两个球，沿面爬电就再加一段绝缘子。但几何里有个常见坑：尖端。实际针尖不可能是数学上的无限尖端点，总有曲率半径，仿真里如果不做处理，尖端附近场强会趋近无穷大，网格稍微一变，结果天翻地覆。

我的处理习惯是：给针尖加一个小倒角或把端点削成半径 0.05~0.5 mm 的圆弧，具体取决于你要模拟的真实针尖尺度。这样电场最大值才是有限值，才能和后处理判据对比。一定要记住：几何上的“理想无限尖”不是真实物理，是数学奇点，COMSOL 会如实算出奇异性，但这个结果没有工程意义。

在 COMSOL 里做空心圆柱等轴对称结构时，建议直接选二维轴对称或二维轴对称+旋转，而不是建完整三维模型。针板间隙本质上是轴对称的（如果边界条件也对），用二维轴对称可以把三维网格量减少几万倍，求解速度显著提升。如果产品结构确实不是轴对称的，比如绝缘子有伞裙还不均匀，那才需要建三维。

2.2 材料参数：介电常数、电导率和击穿场强

材料参数直接决定电场分布是否可信。静电接口只需要相对介电常数：空气 1.00059，环氧树脂常用 3.5~5，SF6 约 1.002，陶瓷类按 6~10 选。电流接口还需要电导率，固体绝缘材料电导率往往是温度的函数，我通常会写成σ=σ₀·exp(-Ea/(k_B·T)) 这样的Arrhenius形式，否则热击穿的正反馈机制就体现不出来。

气体介质二次电子发射系数和电离系数只有在等离子体模块里才需要。空气放电常见的参数包括：电子迁移率 µₑ ≈ 0.04 m²/(V·s)（随约化电场变化），中性气体密度 N 按 2.4e25 1/m³ 算。这些参数如果不想自己翻文献，COMSOL 等离子体模块内置材料库就有 air 的默认反应列表，可以直接用，算常见情况够用。

记住一个原则：宁可参数粗一点也别把量纲搞错。COMSOL 里默认是国际单位制，气体密度单位是 1/m³，电导率单位是 S/m，电压单位是 V。我见过有人把密度输入成 cm³ 量级导致结果差六个数量级，回头查半天。开机之前先把单位意识建立起来。

2.3 网格设计与尖端加密

电击穿仿真的网格关键是局部加密，尤其是高场强区域。针尖附近曲率半径小，电场梯度大，网格不够密就无法捕获场强峰值。我常用的做法是：在针尖附近添加一个细小的圆球/圆形域，把“较细化”甚至“极细化”的网格施加在该区域，其余区域用常规粗细；边界处加边界层网格，保证电场法向变化被解析。

判断网格是否够密，最稳的方法是做网格收敛性测试：把针尖局部网格尺寸从 1/10 尖端半径缩到 1/50 尖端半径，看最大电场是否收敛。如果电场值随网格加密还在明显漂移，说明几何奇异点影响仍然存在，需要继续加密或进一步把尖端倒角加大。

顺带提一句，有些情况会用到“移动网格”接口——比如流注发展伴随气体密度变化或者电极变形时，但这个属于进阶玩法，新手不建议一上来就碰。先固定几何，把场算明白，再考虑移动网格。

3. 三条实操路线：静电判据、热击穿和气体放电

3.1 路线一：静电评估模式（最快、最常用）

静电评估模式适合工程初步校核。具体操作是这样：

第一步，新建模型，选择“二维轴对称”空间维度，添加“AC/DC 模块→静电”物理场接口。几何就弄一个针板间隙：针电极半径 0.5 mm，间隙 10 mm，平板电极厚度 2 mm。

第二步，设置边界条件：针电极赋电位 20 kV，平板电极接地，外部边界按默认“零电荷”。如果是开放边界想模拟无限远处，可以添加“无限元域”，或者把求解域边界拉远到离电极 10 倍间隙以上，工程上足够。

第三步，划分网格，针尖局部细化，求解。后处理里画出电场模的流线图，看一眼最大电场出现在哪。再用“派生值→体最大值”查看 Emax，比如算出 Emax=15 kV/mm。

第四步，对比判据：均匀场空气击穿场强按 3 kV/mm，但针板结构在 10 kV 左右就可能起晕，因为不均匀系数很大。这时最好用帕邢定律或经验曲线去查“这个间隙结构对应的实际击穿电压是多少”，不能只看 Emax。你可以把电压做参数化扫描：从 5 kV 扫到 30 kV，把每个电压下的 Emax 和对应判据关联，就能得到一个“阈值电压”，这就是工程里最关心的答案。

这个模式的好处是快、论断清晰，用来回答领导或客户“会不会打火”够了。缺点是忽略了空间电荷对场分布的影响——当放电发展到一定程度，空间电荷会显著削弱或增强局部电场，那就要上等离子体模块了。

3.2 路线二：固体绝缘热击穿（电流-传热耦合）

热击穿是很多变压器、电缆附件的真实失效模式。固体绝缘在常温下电导率很低，但温度升高后电导率指数上升，漏电流增大，焦耳热增多，形成正反馈。COMSOL 里实现这个过程的设置如下：

物理场选择“电流”接口，材料设置电导率随温度变化表达式。比如我用过的环氧配方里，电导率大致是 σ(T)=1e-12·exp(12000/300−12000/T) S/m 这种量级（具体参数得看材料测试数据）。再加“固体传热”接口，两者多物理场耦合，热源就是“电磁热”的焦耳热项。边界条件方面，电极两端一个接高电位、一个接地；外表面设置对流换热系数，典型自然对流大约是 5~10 W/(m²·K)，别设成绝热，否则温度会虚高。

求解时先做一个稳态扫描：电压从低往高扫，观察温度解是否存在。在临界电压以下，温度场收敛到低温值；越过临界电压，可能就没有稳态解了，数据发散或温度冲向极高值——这个临界电压就是热击穿电压。这种“解不存在”的电压点，正是我要找的击穿界限。

这里有个心得：做参数扫描时光线性的传递过程比较麻烦。解决方法是先锁定线性求解器，配合辅助扫描慢慢逼近临界点。如果直接瞬态求解，时间步长会被拉得极其小，容易卡死。还有，启用“电流连续性开关”可以避免击穿后电流不连续导致的数值振荡，但默认设置通常不需要手动改。

3.3 路线三：气体放电流注动力学（等离子体模块）

等离子体模块是电击穿仿真里最“硬核”的一条路。它求解的是电子密度、离子密度的输运方程和电位方程耦合系统，能计算出流注从电极附近产生、发展到对面电极的全过程。设置步骤相对繁琐，我给出一个能跑通的参考配置：

• 物理场：等离子体→直流放电接口，加上静电接口（等离子体接口会自动耦合泊松方程）。
• 反应列表：使用材料库内嵌的空气反应——包含电子碰撞电离（e+O₂→O₂⁺+2e，e+N₂→N₂⁺+2e）、附着过程、复合过程。默认列表足以复现空气击穿的趋势，不用自己写反应，除非你对定量精度要求特别高。
• 边界条件：负电极（针尖）设电压 -20 kV 或 +20 kV，阳极接地；电极表面二次电子发射系数设典型值 0.02~0.05。
• 初始条件：电子密度和离子密度给一个很小的背景值，比如 1e9 1/m³，不要从 0 开始，否则数值上会有除零问题。
• 网格：尖端局部加密强于其他一切，流注通道后处理时自然能看到细丝状结构。

瞬态求解时，这个模型是极度刚性的。我的建议是选择 BDF 求解器（默认），最大阶数设 2，绝对容差开小一点比如 1e-8 或 1e-9，时间范围从皮秒量级开始，因为电场建立和电子响应在皮秒到纳秒尺度就完成了。如果你要看到流注在几十纳秒内穿透间隙，那时间步长差不多是 0.1~1 ns，步数会非常多，一般需要算数小时甚至更久。

后处理的重点是电子密度云图和电场时序图。流注头部电场会被空间电荷增强，形成局部场增强，你可以观察到电场峰值和电子密度峰值同处移动。做这类仿真，耐心比技巧更重要——跑之前把网格、参数、初始条件全检查一遍，别指望一次成功。

4. 求解器与瞬态参数：最容易翻车的环节

4.1 稳态还是瞬态：这是个选择问题

静电评估肯定是稳态求解，一步到位，没啥好说。热击穿最好是稳态扫描，除非你要看击穿时温度随时间的演化过程——那是瞬态需求。气体放电的流注传播只能瞬态，因为它本来就是一个随时间发展的非平衡过程。

有个常见误区：稳态求解不收敛就认为“击穿了”，这不准确。稳态不收敛可能只是数值原因，比如初始值太差、非线性太强、网格太差。要先排查数值问题，再做物理结论。我通常的做法是：先缩小载荷（电压）到很小，确认问题能收敛，再逐步增大电压，每一步以上一步的解作为初始值，这样能大幅增加收敛性。

4.2 容差、时间步长与求解器配置

COMSOL 默认求解器在多数情况下表现尚可，但电击穿问题常常需要手动干预。

对于静电模块的稳态问题，默认的直接求解器（PARDISO）足够快。对于电流-传热耦合，建议打开“非线性方法→自动牛顿”，并把最大迭代次数调高到 25 以上，阻尼因子设为 0.9，这样不容易因为迭代中途震荡而失败。

瞬态气体放电问题，我用过的最稳配置是：

• 求解器：BDF，最大阶数 2，最小阶数 1
• 时间步进：设置初始步长 1e-13 s，最大步长 1e-10 s
• 容差：绝对容差 1e-8，相对容差 0.001
• 辅助扫描：先把电压固定为击穿电压以下跑一遍预解，有了合理初值再加大电压

这些参数不是拍脑袋定的，而是来自一个核心矛盾：放电中电子输运的时间尺度在皮秒量级（电子迁移率高、速度快），而流注整体发展需要几十纳秒，如果时间步长太大，电子密度就跳变剧烈；步长又会被刚性系统拖得很小。BDF 二阶可以在中高频段保持稳定，但阶数太高误差反而大。这套配置我试过很多次，能稳定跑完流注过程。

4.3 初始值：决定成败的隐形因素

电击穿仿真的初始值设置是很多人忽视的地方。尤其是等离子体模块，如果初始电子密度设得太低（比如 1e6 1/m³），数值上近似“无”的状态会让伽辽金有限元里的非线性迭代找不到方向；设得过高又相当于提前注入了大量电离种子，击穿会过早发生。我常用 1e9~1e12 1/m³ 作为背景密度，这和气体中实际存在的天然电离水平在同一数量级。

电势的初始值也有讲究。静电稳态还好，电流-传热和等离子体问题最好用“先前解”或“辅助扫描扫出来的解”作为初值，不要从零电压零温度开始，否则非线性迭代很容易被初始猜测带崩。COMSOL 里可以设置“依赖其他研究步骤的解”，这一步操作能让很多高难度模型从“不收敛”变成“顺利跑通”。

5. 常见问题与排查技巧实录

5.1 不收敛：先查网格密度还是边界条件？

每次遇到“不收敛”的报错，我的排查顺序固定如下：先看网格质量。电击穿仿真的不收敛十有七八是网格局部太粗糙，电场在尖端附近梯度极大，但网格无法解析，导致残差在高场区集中爆炸。可以用“网格→统计”看一下最小单元质量，如果低于 0.1，先把那个区域重新画一遍。

再看边界条件是否矛盾。比如两个电极都赋了电压却有一个接地、或者悬浮电位设定冲突、或者开放边界直接设零势能，都可能导致无解或振荡。还有一个小技巧：把电压降到原来十分之一，如果模型能收敛，基本确认是强非线性问题；如果“收敛性极好”但电场结果显示处处均匀，那大概率是约束加多了，检查边界条件有没有被误用。

5.2 等离子体模块常见负密度与振荡

等离子体模块算到一半，电子密度出现负值，这是经典的“数值病”。原因大多出在时间步长太大或空间离散过粗，导致输运方程在对流占优区域产生非物理振荡。处理办法：细化流注路径上的网格，减小最大时间步长，或者把输运项的对流格式改为“流线扩散稳定”。COMSOL 里在等离子体接口的离散设置中开启稳定项，很多时候能直接消掉负密度。

还有个被忽视的原因：反应列表里存在极快反应，时间尺度远小于步长，导致刚性振荡。这时把对应反应的速率常数条件改为更平滑的函数，或用查表插值，也能改善。总之，负密度是信号不是最终结果，系统在告诉你“你的时空分辨率跟不上物理过程了”。

5.3 网格依赖性：换网格后结果完全不同怎么办

仿真做完一调整网格，最大场强从 10 kV/mm 变成 30 kV/mm，这种情况我经历过太多次。原因基本就两种：一种是尖端几何存在奇异性，网格越细场强越大；另一种是边界层附近网格不足，法向电场解析不准确。

对于几何奇异点问题，解决方案不是无穷加密，而是把尖端改为有实际曲率半径的形状，然后做收敛性测试；对于边界层问题，添加边界层网格并设定第一层厚度与趋肤深度/德拜长度同数量级即可。做电击穿仿真，强烈建议在交付结论前花半天时间做网格收敛性测试，不然后面换了网格结果变了，汇报时就尴尬了。

5.4 电击穿仿真问题排查速查表

6. 个人实操中的几点体会

最后分享几个我在反复跑电击穿仿真后积累的小心得，算不上什么精深理论，但确实能让工作顺很多。

第一，项目时间紧的时候，先跑静电评估永远是对的。它能快速定位危险区域，给我下一步细化方向。等到客户追问“流注发展和击穿时间是多少”时，再启动等离子体模块也不迟。我自己用这个思路完成了好几个开关柜绝缘评估项目，每个都省了一两周时间。

第二，气体放电仿真不要执着于精确复现文献里的每一条曲线。工程上更关键的是“能不能击穿、击穿电压大约是多少”这类定性结论。参数不完美、反应列表简化，只要趋势和数量级对，就有工程价值。科研导向的朋友另说，但工程实践里，快速可靠的趋势判断往往比精确数字更重要。

第三，网格收敛性测试不是可选项，是必选项。哪怕只是把最大单元尺寸减半再算一次，比较两个电场最大值，你就能摸清结果的置信区间。凡是没做网格测试就敢交付的电击穿报告，都在赌运气。每次交付前把这个测试结果放进报告附件，客户和评审都会更信任你的结论。

COMSOL 电击穿仿真这条路，说难确实难在物理机制多样、数值问题多；说容易也很容易，只要把物理图像、模块选型、边界条件和网格密度这几件事理清楚，大部分问题都能在合理时间内搞明白。希望这篇整理能帮你少走几圈弯路，把时间花在真正有价值的分析和设计上。