★ FEATURED ARTICLE 简介面向电气工程与高压绝缘领域研究者的COMSOL仿真学习资料围绕绝缘材料电击穿与电树枝生成机理展开适合需要理解多物理场建模思路的电气专业学生和工程技术人员。资料包共12个文件其中7个txt技术文档以机理笔记与技术博客为主2个doc分析文档提供仿真分析说明1个html演示页面便于交互查看2张jpg示意图直观展示电树枝形态与电场集中效应压缩后仅111KB结构紧凑、便于按需查阅。内容从材料内部微小缺陷或杂质引发的电场集中切入逐步剖析电荷积累、局部放电、电树枝沿电场方向扩展直至绝缘失效的完整过程并结合COMSOL多物理场建模要点说明材料种类、结构和工作环境对电击穿阈值的影响为聚合物绝缘选型与耐电击穿设计提供参考。目前已有124人学习下载适合作为电树枝仿真入门、机理梳理及科研预研的轻量资料。1. 绝缘材料电击穿为什么难仿真电树枝是多尺度失效的缩影把一块交联聚乙烯试样埋进针-板电极升压到几十千伏几十分钟后从针尖长出一条分叉的碳化通道像一棵倒长的树。这条通道就是电树枝绝缘失效的前兆也是「COMSOL模拟绝缘材料电击穿现象揭示电树枝生成机理」这个标题里要解决的核心问题。它跨了三个尺度宏观电压、微观局部场强放大、介观通道生长。物理上属于高电压绝缘诊断工程上关乎电缆和干式变压器的寿命评估。仿真这件事难在击穿是个概率事件不是简单用「场强超过阈值就击穿」能描述的树枝形态还带有明显的分形和随机性。做这套仿真的读者通常是绝缘方向的研究生、设备寿命评估工程师或者想转行做电气绝缘仿真的从业者。这篇笔记把路径拆成六步几何、材料、电场、判据、生长算法、验证每步都能照着重做。2. 先把针-板电极几何立住曲率半径、材料与维度的选择2.1 针-板结构的三个关键几何量针尖曲率、间距与电场放大倍数实验里最常用的电树枝引发结构就是针-板电极。一根曲率半径很小的金属针插入绝缘介质针尖对着接地平板中间隔着一段绝缘距离。针尖处电场会剧烈集中这个集中程度由三个几何量决定针尖曲率半径 r0、针尖到平板的距离 d、以及试样厚度 t。r0 的典型范围是 1 μm 到 50 μm。曲率半径越小局部场强放大倍数越高越容易在较低外施电压下引发绝缘材料内部局部放电。常见做法是先按实际实验电极的针尖形状定 r0仿真里直接用理想尖点会带来数值奇异后面第五部分专门讲这个坑。d 通常取 1 mm 到 3 mmd 越小平均场强越高树枝萌生等待时间越短仿真里的稳态电场梯度也更陡。电场放大倍数可以用经典双曲面近似估算E_max ≈ 2V / (r0 · ln(1 4d/r0))。这个公式来自针-板电极的解析解近似用来核对数值结果很方便。举例V 20 kVr0 10 μmd 2 mm估出来的 E_max 在数百 kV/mm 量级远超环氧树脂的宏观击穿场强——这正是电树枝能在远低于材料宏观击穿电压下萌生的原因局部场强放大而不是整体场强击穿。建几何时不用刻意去追求显微镜级精度把 r0 和 d 对电场分布的主导作用先复现出来比堆砌几何细节重要得多。2.2 绝缘材料参数取值介电常数、电导率、击穿场强的工程范围材料参数直接决定场强分布形态。对于电树枝仿真最需要关注的三个参数是相对介电常数 εr、体积电导率 σ、以及临界击穿场强 E_th。εr 决定电容性场分布σ 决定直流下的稳态传导场分布E_th 是后续概率生长判据的门槛值。下表列出我常用来做第一轮仿真的参考范围。注意这是工程参考值不是某个标准的规定值实际试样参数必须用实测数据覆盖。材料相对介电常数 εr体积电导率 σ (S/m)直流/工频击穿场强参考 (kV/mm)环氧树脂浇注料3.5 ~ 4.51e-14 ~ 1e-1320 ~ 40交联聚乙烯2.3 ~ 2.51e-16 ~ 1e-1430 ~ 60硅橡胶2.8 ~ 3.21e-14 ~ 1e-1218 ~ 35εr 越高同样几何下场强分布越依赖几何形状而非材料内部差异σ 如果设成零直流稳态下就是纯电容模型电场按介电常数分配这在高分子绝缘里其实并不完全成立因为体积电导率虽然极低但不是零泄漏电流会在长时间加压下改变空间电荷分布。做直流电树枝仿真时我一般先做一组「εr 主导」和一组「σ 主导」的对照看场强分布差异是否显著再决定是否引入空间电荷模型。2.3 用2D轴对称替代三维省算力但不牺牲物理主线的常见做法针-板电极在理想情况下绕针轴旋转对称电树枝生长初期也近似对称。做第一轮仿真时没必要直接上全三维模型用 2D 轴对称截面建立几何COMSOL 会自动按旋转对称处理算出来的场强分布和三维模型基本一致计算量却能降一个数量级以上。具体做法是在几何组件里选择二维使用轴对称坐标系统画出针-板截面。针体用矩形加倒圆角画出来尖端圆角半径就是 r0板电极画成一条接地边界中间是绝缘介质域。边界上针体边界施加高压电位 V板电极边界设接地外部边界用默认的零电荷边界并把计算域往外延展到足够大避免边界效应污染针尖附近的场分布。整个模型域半径取 d 的 5 到 10 倍这个经验值够用。我一般会在这一步把「二维轴对称上的电场结果」和「三维结果」做一次对比测试确认一致性后才放心往下走。很多翻车案例其实就出在几何阶段域开太小、尖端没倒圆角、对称轴选错后面所有结果都是错的。3. 稳态电场求解的完整流程从物理场接口到网格收敛3.1 静电接口建模方程、求解类型与边界条件设置COMSOL 里做电树枝仿真的第一遍求解通常用 AC/DC 模块的静电Electrostatics接口。控制方程就是拉普拉斯方程-∇·(ε0 εr ∇V) 0在纯绝缘域里没有自由电荷。研究类型选稳态Stationary求解器默认用直接法模型规模不大时完全够快。物理场设置的顺序并不复杂。第一步把针电极边界设为「电位」条件输入外施电压 V0第二步把板电极边界设为「接地」V 0第三步材料节点里把绝缘域指认成对应材料填上 εr 和 σ。如果后续要做有损介质的直流稳态场还要在静电接口里开启「传导电流」选项把 σ 参与进求解这时方程变成 -∇·(σ∇V) 0 与电容效应的耦合形态结果更贴近长时间加压后的真实场分布。这一步最常见的问题是把外部边界设成了接地而不是零电荷这相当于在模型边界上强加了一个等势面会把场强分布拉偏。正确做法是让外部边界自然满足零电荷让场线在空气中闭合。判断边界是否设置到位可以直接看后处理里的等势线在边界处是否自然地延伸出去有没有被「切断」的感觉。3.2 针尖网格加密到什么程度才算收敛电场仿真的收敛性几乎完全由针尖附近的网格密度决定。理想尖点是数学奇点局部场强会随网格细化不断增大只要给针尖倒了一个有限圆角 r0场强最大值就应该是有限值并且随网格加密趋于一个稳定平台。网格尺寸量级的设定我通常直接按 r0 的比例来控制。区域推荐网格尺寸说明针尖表面0.1 r0 ~ 0.5 r0决定最大值精度必须细化针尖周围 0.5d 范围0.5 r0 ~ 2 r0捕捉场强衰减梯度其余介质域d/20 ~ d/10对结果影响小可以粗外部空气域渐进增大避免边界反射网格收敛的检验方法不是「网格越小越好」而是把针尖处场强最大值 E_max 记录到一张表里从粗网格到细网格看它变化率。变化率小于 1% 到 2% 就认为收敛了。这个步骤不能省否则后面概率生长模型里用的每一个 E 值都带着网格误差树枝形态必然失真。有个更省事的技巧直接开启物理场控制网格把针尖边界加一层边界层网格厚度设为 r0 的十分之一左右。边界层能很好地捕捉法向场强的快速变化。用全局细化反而浪费算力因为大部分区域场强变化很平缓。3.3 用解析近似核对数值结果多少误差是可以接受的数值结果出来以后第一步不是急着做树枝生长而是拿第一节里的解析近似公式 E_max ≈ 2V / (r0 · ln(1 4d/r0)) 来对账。对不上的原因通常有两个一是针尖倒角不是理想双曲面形状实际几何和近似公式的前提不一致二是网格还没收敛到稳定平台。我一般接受的偏差范围是 10% 以内。偏差超过这个范围优先检查网格收敛曲线和几何圆角半径是否被替换过。为什么这个验证步骤值得做因为后续所有概率生长结果都建立在这张 E 场上如果基础场分布就偏了后面长得再漂亮的树枝都是空中楼阁。这不是玄学而是把数值误差控制在源头。另外如果你的研究对象是工频交流下的电树枝稳态静电场的 E 峰值只是参考更合理的是用 E_rms 作为后续判据输入。这个问题在交流仿真里非常关键第五部分会具体展开。4. 把「击穿条件」翻译成算法阈值场强与电树枝概率生长模型4.1 先定判据从局部场强到损伤累积的映射电树枝生长的物理过程至今没有完全统一的微观解释碰撞电离、电致机械应力、热分解都参与其中。工程仿真里不需要把这些微观机制全部建出来常见做法是用一个等效判据局部场强超过材料临界阈值 E_th 时该位置发生不可逆损伤损伤累积到一定程度就形成树枝通道。问题在于真实电树枝生长表现出明显的随机性——同一样品同样电压每次生长的分支形态都不一样。如果判据写成「E E_th 必生长」仿真出来的树枝会是一条笔直的主干和实际完全不符。因此判据要写成概率形式最常用的是 Weibull 形式P 1 - exp(-(E/E_th)^m)。E 是局部场强E_th 是临界场强m 是形状参数控制随机性的强弱。m 的取值决定了树枝形态。m 很大比如 20 以上高场强区域几乎以概率 1 生长随机性被压制树枝趋于单调m 较小比如 4 到 8次高场强点也有机会被选中自然形成分叉。推荐新手从 m 6 起步把实验结果照片里树枝的分支密度当作调节靶子。E_th 则参考材料的击穿场强表格但要说明E_th 和宏观击穿场强不是同一个量它是局部劣化的临界值通常比宏观值更低。这个参数需要试算几次才能定下来属于模型标定过程。4.2 概率生长核心循环随机数与Weibull概率的最小实现有了概率判据就可以写生长循环。下面这段 MATLAB 代码是我最常用的最小实现把电场矩阵读进来从种子点开始按概率逐点生长。% 电树枝概率生长最小演示 % E_norm : 二维全场强矩阵单位 kV/mm % E_th : 临界劣化场强单位 kV/mm % m : Weibull 形状参数典型 4~10 % rng : 固定随机种子保证结果可复现 rng(42); % 固定种子方便对比参数 E_th 25; % 某环氧树脂标定值 m 6; % 形状参数调小增加分叉 [ny, nx] size(E_norm); % 矩阵尺寸 tree false(ny, nx); % 树枝通道标记 seed_idx find(E_norm max(E_norm(:)), 1); % 选择针尖最亮点做种子 tree(seed_idx) true; for step 1:2000 % 生长步数上限 % 找出所有与树枝相邻但还不是树枝的候选点 candidate find(imdilate(tree, ones(3)) ~tree); if isempty(candidate), break; end % 对每个候选点计算 Weibull 生长概率 E_c E_norm(candidate); P 1 - exp(-((E_c / E_th) .^ m)); P P / sum(P); % 归一化成概率分布 % 按概率随机选一个生长点 r rand(); % 均匀随机数 cumP cumsum(P); pick candidate(find(cumP r, 1)); tree(pick) true; % 忽略visited顺序直接标记生长 end代码逻辑分四块。第一块是初始化和种子选取把场强最大点作为树枝起点对应实验里针尖最先出现放电点。第二块是候选点搜索用 imdilate 求树枝边界纯粹的形态学操作保证每一步只从当前树枝前沿向外扩展。第三块是概率计算核心就是 Weibull 公式 P 1 - exp(-(E/E_th)^m)归一化后变成候选点之间的相对权重。第四块是随机抽样用 rand 生成一个 0 到 1 的随机数在累积分布上找到对应生长的点。参数里最值得玩味的是 m。m 越小概率分布越平坦低场强点也有不小的生长机会树枝分叉多而密m 越大高场强点的权重压倒性领先树枝趋近于沿着最大梯度方向直线生长。E_th 的选取同样敏感设置过高则生长几乎停滞设置过低则整个区域快速变成一坨导电块失去树枝形态。调试顺序是先固定 E_th扫 m再反过来扫 E_th。4.3 把COMSOL的电场结果喂给生长循环数据对接与顺序上面这段代码假设 E_norm 已经是一个规整的二维矩阵但 COMSOL 导出的网格数据不是矩阵而是散点坐标加场值。这一步的数据对接很关键处理不好会少一行、错位一格生长结果完全失真。常见做法是用 COMSOL 的「导出」功能选网格数据集把表达式 E_norm 连同坐标 x 和 y 一起导出成文本文件。然后在 MATLAB 里用 scatteredInterpolant 或 griddata 把散点插值到规整网格上。我一般会用 200 × 200 的规整网格做插值这个密度既能保留场强分布的细节又不会让生长循环慢到难以接受。插值方法推荐 linear不要用 cubic因为电场在针尖附近梯度很大cubic 容易产生过冲把最大值人为抬高。更顺滑的做法是用 LiveLink for MATLAB 直接把 COMSOL 的场数据读进 MATLAB 工作区省掉文件读写和坐标对齐的麻烦。没有 LiveLink 也能跑通只是每次改完几何重新求解后要手动导出一次。这里有一个容易被忽略的顺序问题先做网格收敛测试再用收敛后的网格导出 E 场。如果你先导出了粗网格的 E 场后面长出来的树枝形态即使再好看本质还是在错误输入上做文章。数据对接完成后还有一个可选但推荐的做法把已经生长的树枝点对应的介电常数调低比如从 3.5 降到 2.0重新求解电场再继续下一轮生长。这种「长一步、重算一次场」的迭代方式更接近真实过程因为树枝通道形成后相当于引入了高电导/低介电区域周围场强分布会重新洗牌。代价是每步都要调 COMSOL 求解计算时间会涨几十倍。我建议先跑固定场强的快速版本把模型逻辑调通再考虑升级成迭代重算版本避免一上来就被求解时间劝退。5. 电树枝仿真避坑清单4个高频翻车点和排查顺序5.1 针尖场强一直发散最大值随网格加密只涨不停现象把网格细化了一轮又一轮E_max 非但不收敛反而越涨越高数值已经远超解析近似公式的计算值。原因几何里用了绝对尖点或曲率半径设成了 1e-6 m 以下的极小值这在数学上就是一个奇点场强理论上趋于无穷。网格只是把这个无穷逼近得越来越精确而已。解决给针尖尖端画一个有限尺寸的球冠或圆弧倒角半径设成真实实验针尖的曲率半径水平通常 5 μm 到 20 μm然后做一节「网格尺寸 vs E_max」的收敛曲线确认平台出现后再往下走。这条不解决后面所有概率生长都是假的。5.2 树枝主干太直太对称完全不像真实分支结构现象同样的概率生长代码跑出来的电树枝就是一根直线偶尔才有一条短分叉形态和实验照片里的密分支差很远。原因大概率出在两个地方一是 Weibull 形状参数 m 设得太大高场强点的概率几乎为 1每次随机抽样都稳定选到最大场强点二是随机数种子固定后没有意识到代码里如果存在多处 rand 调用抽样顺序会影响结果。解决把 m 降到 4 到 8 范围重新扫除此之外可以在候选点概率计算之前给 E_norm 加一个 ±5% 的随机扰动人为引入局部不均匀性这在物理上对应材料内部的微小缺陷和气隙。调参时以「分支密度符合实验照片」为目标而不是以「最大场强点必生长」为目标。5.3 直流电场模型模拟工频交流生长速度快得离谱现象把外施电压设为工频有效值用稳态静电场的峰值场强做判据结果树枝生长速度比实验观察快了几个数量级几百步就贯穿整个试样。原因交流下电场随时间正弦变化真正对材料造成累积损伤的是每个周期的过电压区间而不是持续作用的恒稳场强。稳态直流模型默认全场强全程有效相当于把损伤时间高估了一个量级。解决不要直接用 E 峰值喂给概率判据改成先用 E_rms 把生长概率整体压低再配合实验测得的树枝萌生时间曲线标定阈值更严谨的做法是做瞬态求解取一个工频周期内的场强时间序列按周期积分损伤累积量。对新手来说用 E_rms 是性价比最高的修正方案。5.4 换一次网格密度树枝形态就大变样重复性差现象网格加密一档重跑同一次生长模拟树枝通道的形状完全变了甚至总长度都差了一大截。原因概率生长模型本质是对场强分布的随机抽样网格一变局部场强插值结果发生微小变化而随机抽样会放大这些微小差异导致路径整体漂移。解决固定随机种子先做基准对照组然后单独分析网格变化带来的 E 场差异如果 E_max 变化率在 2% 以内则概率差异可接受反之先回去细化网格到收敛。形态对比必须在同一网格基准和同一随机种子下进行否则你比较的其实是随机噪声而不是参数差异。这个习惯一定要养成否则后面每调一个参数都会怀疑人生。6. 仿真结束不算完用三个量化指标验证电树枝形态是否合理概率生长模型跑出来的树枝肉眼看着像不等于物理上合理。我会再用三个量化指标把关分别是盒计数分形维数、通道面积占比、以及生长长度-时间曲线形态。分形维数用盒计数法算把树枝图像覆盖在不同边长 s 的方格下统计覆盖树枝的格子数 N对 log(s) 和 log(N) 做线性回归斜率绝对值就是分形维数。真实电树枝的分形维数通常在 1.3 到 1.9 之间密分支型偏 1.6 以上稀疏型偏 1.4 以下如果你的仿真值低于 1.2说明树枝过直随机性不足。通道面积占比是树枝面积除以整个绝缘视场面积。这个指标能直观反映劣化损伤的严重程度和实验切片图用图像处理算出来的占比对比能快速发现模型是否把不该击穿的大片区域也长满了树。生长长度-时间曲线形态更偏工程。把树枝尖端到针尖的距离随时间作图再转成 log-log 坐标真实实验数据通常会呈现台阶状加速特征萌生期几乎不增长中期快速延伸后期因为通道降低局部场强而变缓。如果仿真曲线是一条光滑直线说明概率模型里缺了一个关键机制——通道形成后的局部场强屏蔽效应。这时候就要回到第四节提到的「长一步、重算一次场」的迭代重算做法。有一回我调试一套环氧树脂电树枝参数试了五组 m 值形态从一根笔直细线到一团模糊墨迹都有却始终对不上实验照片。最后发现不是 m 的问题而是 E_th 设得太低整个高场强区域都成了等效导体树枝失去通道结构。把 E_th 调回标定值后形态立刻收敛。这个教训让我养成了一个习惯每次调参只动一个变量其他全部锁死并且用固定随机种子保证可复现。做这套仿真最怕的就是同时调三个参数结果出了问题都不知道该怪谁。希望这些路径和踩坑记录能帮你在自己的绝缘击穿仿真方向上少走几趟弯路。本文还有配套的精品资源点击获取
阅读完成 · 觉得有帮助？
