基于经验模型的熔断器燃弧计算仿真

来源:用户上传      作者:陈永洁

　　摘要：熔断器的燃弧特性决定了其分断过电压与分断能力，所以，为了研究其燃弧特性，建立了经验模型。本文首先利用ANSYS有限元软件计算了熔断器的弧前时间;分析了电弧燃烧的物理过程，并进行了一定简化，基于此基础上构建了熔体烧蚀、石英砂烧蚀、电弧电压和外部电路模型以及电弧与石英砂之间的热交换等;然后通过Matlab计算，将获得的电流、电压曲线和实验曲线经过拟合确定与电弧相关的物理参数，引入燃弧过程中电导率温度变化的曲线;最后，基于该模型，提出实践研究对象，并分析了电弧燃弧过程和熔体结构尺寸，并介绍了电弧过程中电导率随温度的变化曲线;最后，基于该模型分析了电弧燃烧过程和熔体结构尺寸对弧压峰值、过零时间的影响。而基于仿真实践研究发现，随着熔体弧压峰值的上升，减少了过零时间;保持熔体总长度的状态，更改断口数或断口长度不会影响电弧特性。
　　关键词：半经验数学模型;熔断器;燃弧;电弧
　　熔断器的短路可以分为两个阶段，分别是弧前和燃弧。弧前在燃弧阶段，会产生较大的燃弧能量和过电压。而过电压会给设备绝缘和系统产生安全威胁，导致设备发生损坏或是导致短路电流重燃，出现燃弧灾害。燃弧能量过高会导致高压线路熔断器爆炸，从而发生分断失败。为了提高熔断器的可靠性，需要对于其电弧特性进行科学合理的控制，并计算电弧模型[1]。本文首先分析了熔断器的分段机理，并计算了熔断器弧前时间。并以此为基础来考虑燃弧物理过程，构建数学模型，同时考虑电弧温度、电导率关系。最后，采用Matlab 来求解模型，能够获得众多的变量，包括电弧温度、电弧形态的扩展、电流弧压曲线等。通过实践研究来持续优化模型，并分析熔w结构参数改变对燃弧特性的影响规律。
　　一、熔断器弧前计算
　　短路故障电流经过熔断器时，因为狭颈的电流密度最大，所以这一部分会聚集大量热量，最先出现熔断，并出现电弧。因此，这一时间可以认为是弧前时间。而弧前时间计算可以使用的方法非常多，包括有限元法、有限差分法等众多计算方法。J.G.leach 、A.Wright等人最早将有限差分法应用到弧前时间计算领域，但随着各种仿真软件的快速出现，如Ansys、Ansoft 等，所以弧前时间计算开始使用有限元法，并且这一方法也综合考虑到了熔化潜热、熔化吸热以及石英砂导热等物理过程的其他影响因素[2]。
　　例如，有一1500Vdc 600A 额定熔断器，利用有限元软件进行APDL编程，可以得到其对称熔体结构的弧前时间，熔体材料是银，各种数据为常规参数。边界条件：除了上边界恒定温度20K之外，其余的所有界面绝热;由电流上边界流入，下边界流出，大小为上升率10A/us的斜坡电流。基于这一数据，从而得到计算之后熔体温度情况。
　　二、熔断器燃弧过程数学模型
　　构建已知弧前时间的燃弧过程数学模型，然后构建如下图1所示的单个狭颈熔体结构的研究对象。在该模型中，石英砂处于十分致密的状态下，将熔体包裹在中间。弧形区域是图1中的灰色部分，从狭颈中心开始电弧起弧，这一阶段的温度回逐渐上升，并且开始向周围扩展，侵蚀周围的石英砂。与此同时，石英砂会吸收能量，在大量能量被吸收之后，开始出现温度下降，且限制电流过零[3]。
　　（一）模型假设
　　电弧分为阴极压降区域、弧柱区域、阳极压降区域三个区域。电弧燃烧时的总功率表现为电弧的温升、石英砂熔化与汽化、熔体的熔化和汽化三种形式。由于阳极和阴极的长度只有10 -3 mm，且弧压值仅为10V 左右，电弧柱占据了两个电极之间的间隙长度。在模型中，忽视了电弧温度对石英砂的热传导作用，仅考虑到了电功率对熔体的烧蚀。此时，电弧温度保持在18000 K左右。当电弧电压缓慢降低时，弧柱区域内的电弧将会和周围的石英砂之间产生热交换，从而导致电弧温度缓慢下降[4]。
　　基于此，可以假设燃弧过程如下：1.电流均匀流过每个狭径，所有狭径同时起弧;2.弧宽为熔体的总宽度，即电弧发展方向为长宽;3.熔化温度到汽化点部分为电弧区;4.电弧温度从18000k开始，在电弧电压缓慢下降阶段开始下降;5.电极附近电弧的电功率降用于熔化烧蚀，即弧长的延伸方向。
　　（二）建立电弧模型
　　根据上述假设构建本文的电弧数学模型，具体包含电弧石英砂热交换、电弧电压、石英砂烧蚀、熔体烧蚀、外电路等。
　　1.熔体烧蚀模型
　　基于电热场守恒原则来，构建熔体烧蚀温度场模型，以此来对于电弧长度方向的扩展情况求解。首先，基于熔体内部的内热源项、热传导作用，引入物质的焓，计算焓与温度的对应关系，从而能够得到温度场分布情况，更加了解温度和银熔体汽化点温度两者之间的大小。焓的单位为KJ/KG，主要表示为单位质量物质的全部内能。银熔体汽化点温度低于温度的话，则表示该部分才是电弧扩展部分[5]。
　　2.数学方程
　　（1）
　　 （2）
　　其中，Cp 为定压比热容;H 为焓;为电流密度;a为熔化流体分数;Lr 为银的熔化潜热;λ是热导率ρy表示密度;γ为电导率;k为温度系数;ΔH 为潜热项;G 为银的汽化潜热;β为汽化气体分数;h 为不考虑潜热物质的焓。
　　3.熔体潜热处理
　　对于银片熔化潜热、汽化潜热进行计算，引入焓来表示模型，得到银片温度场分布情况。在潜热过程中，逐渐注入能量，焓增加但温度不变。要想使温度、焓相对应，则假设在潜热时发生温度改变，但是改变程度轻微，函数实现温度与焓之间的相互转化。
　　4.边界条件
　　电弧不仅会对银片烧蚀，使其向长延伸，还会对周围石英砂烧蚀，从而实现厚度延伸。所以，根据能量平衡原则建立模型，石英砂烧蚀的能量来源于弧柱电功率，即石英砂汽化吸热和弧柱电功率彼此之间的相应功率一致。根据模型可以获得每一个时间步长内电弧不同切片上的电弧截面的扩展，熔体的初始总宽和电弧宽度相同，表示电弧厚度在扩展[7]。

nlc202205301612

转载注明来源:https://www.xzbu.com/1/view-15432449.htm