企业官网建设流程全解析

电弧研究这件事，在电力系统里可以说是既基础又难啃的硬骨头。开关分合闸、故障开断、绝缘配合，哪一样都绕不开电弧。可是电弧本身看不见摸不着，就算用高速摄像机拍下了燃弧过程，录到了弧压弧流波形，拿到手的也只是一堆现象数据。想定量描述电弧动态行为，尤其是给断路器仿真模型提供可信的弧隙电导变化规律，就必须借助数学模型的支撑。Mayr电弧模型和Cassie电弧模型，正是在这个需求下被反复验证的两个经典工具。

我先说结论：这两个模型都不是"哪个更先进、哪个被淘汰"的关系，它们各自抓住电弧物理过程的不同侧面。Cassie模型适合描述大电流燃弧阶段的弧隙动态，Mayr模型擅长刻画电流过零前后的电导骤变。实际工程里，更多时候是把它们配合着用，甚至做成组合模型。这篇文章我会从两个模型的物理假设和数学推导讲起，再给出一套在MATLAB/Simulink里的建模仿真步骤，用同一组断路器参数对比两者的仿真波形差异，最后聊聊参数整定里那些文档不会写、但实战中非常关键的坑。想搭电弧模型做开断仿真，或者在做高压电器研发、电力系统暂态分析的同行，这篇应该对你有用。

1. 两个模型的前世今生：电弧的"数学模型化"到底解决了什么问题

1.1 电弧为什么难建模

先聊点背景。电弧本质上是一个等离子体通道，温度可以到上万开尔文，弧柱直径随电流变化剧烈，而且整个过程发生在毫秒甚至微秒级时间内。电弧间隙的电阻不是固定值，它随弧道电离程度、散热条件、气流吹拂状态实时变化。想用简单的欧姆定律去描述它，根本行不通。

但电力系统仿真又特别需要电弧的动态电阻模型。比如断路器开断短路电流时，弧隙电阻从接近零的燃弧状态，到电流过零瞬间骤增到近乎绝缘状态，这个电导变化过程直接决定了断路器能否成功开断。如果弧隙电阻变化建模不准，电弧熄灭后的恢复电压计算、截流过电压评估、开断失败的重击穿分析全都失真。所以从上世纪三十年代开始，研究者们就想方设法用数学公式描述电弧电导的动态过程，Mayr和Cassie就是这条探索路上最著名的两个成果。

1.2 两个模型的诞生背景与定位

Cassie模型是1939年提出的。当时主要为了研究空气开关和高压断路器在大电流燃弧阶段的电弧特性。Cassie把电弧看作一个直径可变、但内部温度和能量密度基本恒定的气体通道，能量主要通过气流对流和辐射耗散。这个假设对应的是大电流、强气流吹拂的燃弧场景。

Mayr模型是1943年提出的，恰恰相反，它把电弧当成一个直径基本固定、但温度和电导率剧烈变化的等离子体柱，能量损耗主要是热传导和热游离过程控制的径向散热。这个假设对应的是小电流、电弧即将过零熄灭时的场景。

两个模型相差四年提出，但出发点完全互补。Cassie抓的是"电弧通道的几何尺寸变化"，Mayr抓的是"电弧通道的热状态变化"。

1.3 从工程需求理解模型的价值

现代工程里为什么这两个几十年前的老模型还在用？因为它在宏观层面抓住了电弧电导变化的本质特征，计算量又小，非常适合嵌入电力系统电磁暂态仿真程序。比它更精细的模型当然有，比如基于CFD的磁流体动力学电弧模型，能算出电弧内部温度场、流场、电磁场的分布，精度高但计算开销惊人，一个开断过程可能要跑几小时甚至几天。对于系统级仿真来说，几百个节点同时计算，显然不能每个断路器都上CFD。所以Mayr和Cassie这类集总参数模型依然是主流选择。

2. 模型原理解析：从物理假设到数学表达

2.1 Mayr电弧模型的核心思想

Mayr模型的出发点是电弧弧柱的能量平衡。把电弧看作一个圆柱形等离子体，它的电导率取决于弧柱温度，而温度又由注入能量和耗散能量的平衡决定。电弧电流注入焦耳热，弧柱同时向周围介质散热，当两者平衡时，电弧稳定燃烧。如果平衡被打破，弧柱电导就会动态变化。

数学上，Mayr模型用微分方程描述电导的动态变化：

1/g × dg/dt = (1/τM) × (u×i/P0 - 1)

其中：

• g是电弧电导（西门子）
• u是弧隙电压（伏）
• i是电弧电流（安培）
• P0是电弧的散热功率（瓦），在Mayr模型里被假设为常数
• τM是电弧时间常数（秒），表征电导变化对能量失衡的响应速度

这个公式的物理含义非常直观：当注入功率u×i大于散热功率P0时，弧柱能量增加，电导上升；反之电导下降。变化的速度由时间常数τM决定。

2.2 Cassie电弧模型的核心思想

Cassie模型的假设完全不同。它认为电弧通道截面半径可变，但弧柱内部的温度、电场强度保持恒定。能量的耗散不是常数，而是和电弧通道的截面积成正比，也就是和电弧电压的平方相关。这种假设更适合描述电流很大时电弧通道被强烈电离、弧柱整体像导体一样的状态。

Cassie模型的微分方程形式：

1/g × dg/dt = (1/τC) × (u²/Uc² - 1)

其中：

• Uc是电弧电压常数（伏），可以理解为电弧稳态燃弧时的弧压降
• τC是Cassie时间常数

注意这里的区别：Mayr模型的判定条件是u×i与P0比较，也就是用功率判断。Cassie模型用u²与Uc²比较，也就是用电压判断。Cassie认为，电弧电导的变化主要取决于弧压偏离稳态弧压的程度——弧压超过Uc越多，电流通道扩张越快，电导增长越快。

2.3 两方程的统一视角：能量平衡框架

如果把两个方程并排看，其实它们都属于"能量平衡"框架的变体。通用形式可以写成：

1/g × dg/dt = (1/τ) × (输入能量/耗散能量 - 1)

Mayr把耗散能量固定为常数P0，认为输入能量是实时的功率u×i； Cassie把耗散能量表示为与电压平方相关的项Uc²，输入能量用u²来代表。

所以选择哪个模型，本质上是选择哪种能量耗散规律更贴合你要仿真的电弧阶段。

实际测试中，这两者有一个非常直观的差别：Mayr模型的电导对电流过零附近的变化极其敏感，表现为电弧电导在过零前快速下降、过零后几乎归零；而Cassie模型在大电流燃弧阶段的电导变化相对平缓，稳定地维持一个"大电导、低弧压"的状态。这个差异我后面用仿真波形给你们看。

注意：两个模型里的时间常数τ不是随便取的，它跟灭弧介质的冷却速度、气压、触头开距都有关系。不同文献里给的范围很宽，从零点几微秒到几十微秒都有。后面第5节我专门讲参数怎么选。

3. 仿真实操：在Simulink里从零搭建两个电弧模型

3.1 建模思路与准备

我把整套仿真放在MATLAB/Simulink环境里做，因为它的Simscape Electrical和自定义S-Function都比较好用。当然你在ATP-EMTP或PSCAD里也能实现，只是实现方式不同——ATP里可以用TACS控制的时变电阻，PSCAD里可以用自定义元件。我这里用最朴素的方式：用微分方程直接搭自定义模块。

仿真的思路是这样的：电路拓扑用最经典的单相短路开断电路，电源经断路器的弧隙向负载供电，断路器开断时弧隙两端产生电弧，我们的模型就是计算这个弧隙的动态电阻。

需要的参数：

• 系统电压：10kV
• 短路电流：20kA（有效值）
• 工频：50Hz
• 断路器开断时刻：设在对电流过零前3ms开始分闸

3.2 Cassie模型的Simulink实现

先建Cassie模型。它的微分方程是：

dg/dt = (g/τC) × (u²/Uc² - 1)

需要两个输入：弧隙电压u和当前电导g的反馈。弧压u由电路计算得到，g反馈回电路决定电流。

在Simulink里我习惯用MATLAB Function模块写这个方程，代码非常短：

function g = cassie_arc(u, g0, Uc, tauC, dt) persistent g_prev; if isempty(g_prev) g_prev = g0; end dg_dt = (g_prev / tauC) * (u*u / (Uc*Uc) - 1); g = g_prev + dg_dt * dt; g = max(g, 1e-6); % 防止电导为0或负值导致仿真崩溃 g_prev = g; end

这里我给电导加了下限1e-6 S，这是因为电弧电导在理论上可以趋近零，但数值计算中一旦电导为负，电流就会反转，产生完全错误的波形。这是一个非常实用的防错处理。

3.3 Mayr模型的Simulink实现

Mayr模型的实现稍微不一样，因为它的微分方程涉及功率项：

dg/dt = (g/τM) × (u×i/P0 - 1)

注意它需要电流i作为输入，而电流本身又和电导g有关：i = u × g。所以存在一个代数环问题。我实际搭建时先解出电流，再计算电导增量，做一个离散时间的迭代求解。

function g = mayr_arc(u, g_prev, P0, tauM, dt) i = u * g_prev; dg_dt = (g_prev / tauM) * (u * i / P0 - 1); g = g_prev + dg_dt * dt; g = max(g, 1e-6); end

3.4 仿真电路与求解器设置

电路部分我直接用了Simscape Electrical里的交流电压源、电阻电感负载、以及受控电流源来模拟弧隙。受控电流源的值等于u×g，把弧隙等效成一个非线性电导。

求解器方面有两点比较关键。第一，不能只用定步长求解器，因为电弧模型的动态时间常数在微秒级，而工频周期是20毫秒，两者尺度差距太大。我用了变步长ode15tb求解器，它能比较好地处理这种刚性问题。第二，仿真最大步长要限制在20微秒以内，否则时间常数只有几微秒的动态过程根本算不准确。

4. 仿真实例与结果对比：同一工况下两个模型的差异

4.1 参数设置说明

为了让两个模型的对比有意义，我用了同一个开断工况，分别独立运行。

一次典型的仿真参数：

• Cassie模型：Uc = 800V，τC = 5μs
• Mayr模型：P0 = 20kW，τM = 1μs
• 初始电导：g0 = 1e5 S（燃弧初期的强导通状态）

为什么这样取值？后面第5节会详细讲计算依据。这里先看结果差异。

4.2 电弧电压波形对比

电流过零前大约1毫秒这段区间，两个模型的电压波形差异非常明显。

Cassie模型的电弧电压基本被"锁定"在Uc附近。电流从几十安培到几千安培变化时，弧压始终维持在700到900V这个窄范围内波动。这是因为Cassie假设电弧电压只取决于通道截面，而截面变化滞后于电流变化，所以弧压表现出"钳位"效应。

Mayr模型的电弧电压则完全不是这个走势。随着电流向过零推进，Mayr的电弧电压缓慢上升，在接近过零时迅速拉高到几千伏，形成一个明显的"尖峰"。这是因为Mayr的电导在电流变小时快速下降，为了继续传导同样的电流，弧压被迫升高。这个现象在真实的断路器开断过程中也存在——电流过零前电弧电压急剧上升，通常被称为"弧前电压尖峰"。

4.3 电弧电导变化趋势对比

电导变化更能看出两者本质差异。

Cassie电导在燃弧区间保持一个相对稳定的高值，比如10⁴到10⁵ S量级，说明弧道处于强导通状态。到电流过零时才看到电导开始下降，但下降曲线比较平滑。

Mayr电导则是全程跟随电流变化。电流峰值时电导达到最大值，电流一减小，电导就开始快速衰减。在过零附近，电导可以从几千西门子骤降到零点几西门子，跨度接近四个数量级。

这直接印证了我前面说的：Cassie模型擅长描述"导通期间"的弧道，Mayr模型擅长描述"开断瞬间"的电导骤变。

4.4 对开断能力评估的影响差异

这个差异落到工程应用上，就是评估断路器开断能力时用错模型会得到完全不同的结论。

举个例子，我用同一个断路器模型分别计算它开断20kA短路电流时的弧后介电恢复强度。用Cassie模型得到的弧后弧隙电阻恢复速度偏慢，因为Cassie认为电弧电导变化受限于通道截面收缩，而截面在过零后仍然保留一定尺寸。用Mayr模型得到的弧后电阻恢复速度则极快，因为热游离在小电流下迅速衰减，电导呈指数式下降。如果只用一个模型去校核"电流过零后弧隙能否承受恢复电压"，很可能得到截然相反的判断。

所以不少商业软件在断路器模块里直接内置了Mayr-Cassie组合模型：大电流区间走Cassie，接近过零时切换到Mayr。这种切换方式我在第6节的组合模型应用里再展开说。

5. 参数整定与灵敏度分析：仿真结果准不准就在这几个数

5.1 参数初值怎么估

很多初学者上来就问"Mayr和Cassie参数用默认值行不行"。答案很干脆：不行。这两个模型的参数直接决定仿真结果的量级，必须根据你的工程对象来整定。

先说Mayr的P0。它代表电弧稳态燃弧时的散热功率。工程估算上，可以用稳态电弧电压U_arc与额定电流I_arc的乘积来近似。比如某型号真空断路器在小电流燃弧阶段稳定弧压约200V，对应电流100A，那么P0 ≈ 200×100 = 20kW。注意这个值不是常数，它在燃弧不同阶段可能不同，但作为参数整定的起点已经足够。

再说法定的τ。时间常数没法直接测量，但可以用经验公式估算：τ ≈ k × (燃弧时间常数参考值)。对于空气开关，τ通常取数微秒；SF6断路器相对快一些，0.5到2微秒；真空断路器在小电流阶段可以到5微秒左右。这个值强烈影响电导变化速度，建议大家做一次参数扫描来确定，我下面会说具体方法。

Cassie模型的两个参数里，Uc相对好确定，它就是电弧稳态燃弧时的弧压。你可以用短路试验波形里燃弧平台期的弧压平均值来标定。τC同样用参数扫描拟合。

5.2 参数校准的实操方法

有条件的团队，建议用真实断路器开断实验的录波数据做参数拟合。

方法不复杂：对录波的弧压弧流波形做平滑处理后，先用公式g(t) = i(t) / u(t)算出实测电导曲线，然后把仿真电导曲线和实测电导曲线对齐，在一个包含至少一个完整的燃弧到过零的时间窗口内，用最小二乘或遗传算法优化τ和P0（或Uc），使误差最小。

我常用MATLAB的Global Optimization Toolbox做这个拟合，目标函数就是仿真电导与实测电导的均方根误差。整个流程跑下来，参数精度可以做到实测波形的90%以上匹配度。

5.3 灵敏度分析的结论

我做过一轮灵敏度扫描，结论如下：

• Mayr模型对τM最敏感。τM从1μs变到5μs，电弧电压峰值可以从3kV降到1.2kV，电导过零恢复速度也差一个数量级。这个参数定了整个模型的"性格"。
• Cassie模型对Uc更敏感，因为Uc直接锁定稳态弧压。Uc变化10%，燃弧期的弧压就跟着变化10%。τC的影响反而平缓一些，主要影响过零附近的过渡速度。
• 初始电导g0对燃弧初期波形影响很大，但对过零之后的动态过程影响不大。所以如果你只关心开断性能，初始值设得不要太离谱就行。

提示：做参数扫描时千万别用默认的幅值很大的阶跃电流作为激励，电弧模型对这种突变非常敏感，容易产生虚假振荡影响拟合效果。用实际开断波形或者一个完整的正弦短路电流作为激励，拟合出来的参数才可信。

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

6.1 仿真跑着跑着就发散了怎么办

这是老问题。电弧模型本质上是强非线性方程，电导变化好几个数量级，数值求解非常容易出现发散。我排查的顺序一般是这样：

第一步，先检查初始电导。g0如果设成很大的值比如1e8 S，电路求解器很容易把电压算到极小值，导致数值精度问题。我一般控制在1e5以下。

第二步，看求解器。电弧模型微分方程带有微秒级的时间常数，和整个工频电路耦合后会形成刚性系统。默认的ode45基本跑不动，换成ode15s或ode23tb会好很多，而且需要把最大步长设置到微秒级。

第三步，检查电导下限值。我前面提过那个max(g, 1e-6)的处理，就是为了防发散。没有这个保护，电导一旦跌到0附近，除法运算直接报错或者输出NaN。

6.2 电弧电压波形出现异常振荡

有段时间我用Cassie模型仿真燃弧过程，波形总是在过零前出现高频振荡。后来定位到原因：电路里的寄生电容和电弧电导组成了一个LC回路，而我的变步长求解器在这个局部回路上没走足够多的步数。

解决办法很直接：把最大仿真步长从20微秒改到5微秒，同时在电弧模块输出端并联一个很小的数值阻尼电容（比如0.1μF级别），振荡立刻消失。这个电容在物理上也有对应——断路器触头之间本来就存在杂散电容。

6.3 Mayr和Cassie怎么组合使用

最后分享一个组合方案。我现在的断路器开断模型就是这么做的：

用一个开关逻辑判断电弧电流绝对值。当电流大于某个阈值（比如几百安培）时，认为电弧处于大电流燃弧阶段，启用Cassie方程。当电流降到阈值以下、接近过零时，切换为Mayr方程。

切换点一定要做平滑过渡处理。最简单的办法是用一个sigmoid权重函数在两个模型输出的电导之间做加权：

g_combo = w × g_cassie + (1 - w) × g_mayr

其中w从1渐变到0，过渡时间取10到20微秒。直接硬切换会出现电导跳变，反映到电压波形上就是一个毛刺，非常难看。

注意：组合模型的难点不只是切换逻辑，还有参数匹配。两个模型的参数要保证在切换点附近，电导值和电导变化率都不能有太大的阶跃，否则数值上虽然不报错，物理上却完全不合理。

6.4 小电流开断场景下模型选择

还有一个容易踩的坑。上面这套组合方案针对的是短路电流开断这种大电流工况。如果你仿真的是空载变压器分闸或者小电感电流开断，电流幅值本身就不大，这时候直接用Mayr模型反而更好。因为小电流下电弧本来就处于热游离主导的状态，Cassie的大电流假设完全不适用。

我做过一组对比，开断2A容性小电流时，用Cassie模型算出的截流值比实测高了三倍，Mayr模型则和实测结果很接近。这个案例说明一个道理：模型选择不是越复杂越好，而是越贴合物理场景越好。

结尾

我在做断路器仿真这条路上，Mayr和Cassie这两个模型前前后后折腾了大半年才真正吃透。回头最深的体会是：这两个模型公式都不长，但背后那两个物理假设——散热恒定还是电压锁定——决定了它们各自的适用范围，也决定了你在工程中能不能用对场景、整定对参数。纸上谈兵看公式觉得两者就是个数学变形，真到仿真波形对不上实测的时候，才会意识到物理图像远比公式重要。如果你正卡在电弧仿真的参数调试上，建议先从第5节的灵敏度分析入手，把时间常数和特征电压参透，很多问题会迎刃而解。后面有时间我再写一篇组合模型在ATP-EMTP里的实现细节，那里面的元件封装和TACS控制连接又有一堆值得说道的细节。