电弧研究这件事在电力系统里可以说是既基础又难啃的硬骨头。开关分合闸、故障开断、绝缘配合哪一样都绕不开电弧。可是电弧本身看不见摸不着就算用高速摄像机拍下了燃弧过程录到了弧压弧流波形拿到手的也只是一堆现象数据。想定量描述电弧动态行为尤其是给断路器仿真模型提供可信的弧隙电导变化规律就必须借助数学模型的支撑。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); end3.4 仿真电路与求解器设置电路部分我直接用了Simscape Electrical里的交流电压源、电阻电感负载、以及受控电流源来模拟弧隙。受控电流源的值等于u×g把弧隙等效成一个非线性电导。求解器方面有两点比较关键。第一不能只用定步长求解器因为电弧模型的动态时间常数在微秒级而工频周期是20毫秒两者尺度差距太大。我用了变步长ode15tb求解器它能比较好地处理这种刚性问题。第二仿真最大步长要限制在20微秒以内否则时间常数只有几微秒的动态过程根本算不准确。4. 仿真实例与结果对比同一工况下两个模型的差异4.1 参数设置说明为了让两个模型的对比有意义我用了同一个开断工况分别独立运行。一次典型的仿真参数Cassie模型Uc 800VτC 5μsMayr模型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控制连接又有一堆值得说道的细节。
Mayr与Cassie电弧模型原理及Simulink断路器仿真实战
电弧研究这件事在电力系统里可以说是既基础又难啃的硬骨头。开关分合闸、故障开断、绝缘配合哪一样都绕不开电弧。可是电弧本身看不见摸不着就算用高速摄像机拍下了燃弧过程录到了弧压弧流波形拿到手的也只是一堆现象数据。想定量描述电弧动态行为尤其是给断路器仿真模型提供可信的弧隙电导变化规律就必须借助数学模型的支撑。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); end3.4 仿真电路与求解器设置电路部分我直接用了Simscape Electrical里的交流电压源、电阻电感负载、以及受控电流源来模拟弧隙。受控电流源的值等于u×g把弧隙等效成一个非线性电导。求解器方面有两点比较关键。第一不能只用定步长求解器因为电弧模型的动态时间常数在微秒级而工频周期是20毫秒两者尺度差距太大。我用了变步长ode15tb求解器它能比较好地处理这种刚性问题。第二仿真最大步长要限制在20微秒以内否则时间常数只有几微秒的动态过程根本算不准确。4. 仿真实例与结果对比同一工况下两个模型的差异4.1 参数设置说明为了让两个模型的对比有意义我用了同一个开断工况分别独立运行。一次典型的仿真参数Cassie模型Uc 800VτC 5μsMayr模型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控制连接又有一堆值得说道的细节。
