低压塑壳断路器中电动斥力的三维有限元非线性分析与实验研究

1、低压塑壳断路器中电动斥力的三维有限元非线性分析低压塑壳断路器中电动斥力的三维有限元非线性分析与实验研究与实验研究 低压塑壳断路器中电动斥力的低压塑壳断路器中电动斥力的 三维有限元非线性分析与实验研究三维有限元非线性分析与实验研究 李兴文，陈德桂，向洪岗，李志鹏，刘洪武李兴文，陈德桂，向洪岗，李志鹏，刘洪武 （西安交通大学电气工程学院，陕西（西安交通大学电气工程学院，陕西 西安西安 710049710049） 3-D3-D FINITEFINITE ELEMENTELEMENT NONLINEARNONLINEAR ANALYSISANALYSIS ANDAND EXPERIMENTALEXPE

2、RIMENTAL INVESTIGATIONINVESTIGATION OFOF ELECTRO-ELECTRO-DYNAMICDYNAMIC REPULSIONREPULSION FORCEFORCE ININ MoldedMolded CaseCase CircuitCircuit BreakersBreakersLI Xing-wen, CHEN De-gui, XIANG Hong-gang, LI Zhi-peng, LIU Hong-wu （School of Electrical Engineering, Xian Jiaotong University, Xian 710049

3、, China） A ABSTRACTBSTRACT: : For the optimization design of molded case circuit breakers (MCCBs), it is necessary and important to calculate the electro-dynamic repulsion force acting on the movable conductor. With 3-D finite element nonlinear analysis, according to the equations among current-magn

4、etic field-repulsion force and taking into account the ferromagnet, we introduce contact bridge model to simulate the current constriction between contacts, so Lorentz and Holm force acting on the movable conductor and contact, respectively, can be integrated to calculate. Coupled with circuit equat

5、ions, the opening time of movable contact also can be obtained using iteration with the restriction of contact force. Simulation and experiment for repulsion force and opening time of five different configuration models have been investigated. The results indicate that the proposed method is effecti

6、ve and is capable of evaluating new design of contact systems in MCCBs. K KEYWORDSEYWORDS: : Low-voltage apparatus; Finite element analy-sis; Repulsion force 摘摘要：要：电动斥力的计算对低压塑壳断路器（MCCB）的设计有重要意义。基于电流-磁场-电动斥力之间的方程，并考虑铁磁物质的影响，应用三维有限元非线性分析，引入圆柱导电桥模型作为接触点模拟触头间的电流收缩，统一计算触头间的 Holm 力和动导电杆上的 Lorentz 力。耦合电路方程

7、，并以动触头上的预压力为约束通过迭代运算，从而确定动触头的打开时间。分析了五种不同结构的 MCCB，并对电动斥力和触头打开时间进行了实验研究，结果表明，该方法是有效的，可用于 MCCB 新产品触头系统的设计。 关键词：关键词：低压电器；有限元；电动斥力 1 1 引言引言 当短路电流产生时，低压塑壳断路器（MCCB）中动、静触头在电动斥力的作用下分离，在操作机构的带动下，电弧将在触头之间产生，拉长，并在吹弧力的作用下进入栅片被分割成若干个短弧，而后熄灭。在这个过程的起始阶段，作用在动触头上的电动斥力F F（包括导电回路产生的 Lorentz 力F FL及触头间由于电流收缩产生的 Holm 力F

8、FH）和预压力决定了触头的斥开时间和打开速度，从而对 MCCB 的限流性能产生重要的影响。而且，从式(1)所示的 Holm 公式可以看出，接触点半径r与预压力F FK、触头材料的布氏硬度H、触头表面接触情况（用x描述，其范围一般为 0.30.6，通常取 0.45)有关。而 Holm 力F FH与r，触头半径R、以及电流大小i有关1，这样F FK对F F就有一定的影响；另一方面，F FK的选择也由于发热容许的要求，而受到F F的限制。因此，在进行 MCCB 的触头导电回路的设计时必须进行电动斥力的计算。近年来，随着计算机技术发展和数值计算方法的不断进步，计算机辅助工程在 MCCB 新产品的设计中

9、得到了越来越广泛的应用2-3，大大缩短了产品开发周期，降低了开发成本，同时产品的性能也得到了很大提高。对于电动斥力的计算，国内外的学者做了很多有益的工作。在文献4-5中，由于 Holm 公式已经被证明在计算孤立的触头间电动斥力时是正确的6，因而用式(1)来计算F FH ，而忽略触头区域的电流收缩。这种简化不仅使载流导体中的电流分布和实际相差甚远，而且会进一步影响了磁场的分布，从而使F FL的计算结果有一定的误差。 本文应用三维有限元分析，引入圆柱体导电桥模型作为接触点，以五种不同结构的 MCCB 产品为对象，着重研究了电动斥力和触头斥开时间的计算方法，分析了不同结构的触头系统对于电动斥力的影响

10、。并和实验结果进行了比较，表明仿真所用的方法是正确的，可以应用于 MCCB 新产品的触头系统的开发。同时，计算结果也可作为进一步应用虚拟样机技术进行 MCCB 开断过程仿真时的输入数据。2 2 计算方法计算方法 2.12.1 导电桥模型导电桥模型 应用数值计算方法综合考虑触头间的电流收缩和导电回路对于作用在动导电杆上电动斥力的影响时，必须首先引入一个合理的计算模型来描述触头间的电接触情况。R.holm 在推导式(1)中电动斥力FH解析式时，为了分析上的需要，假定接触导体为超导小球1。本文为了和实际情况更加接近，用位于触头中心的圆柱体导电桥模型来模拟导电斑点，其材料性质也和触头材料相同，半径r可

11、由式(1)所示的 Holm 公式计算。为了确定导电桥高度参数h，对一对圆柱体触头进行了电动斥力仿真，图 1 所示为其截面的示意图。结果发现，在相同的r下，其高度h在 0.10.25mm 范围内对电动斥力的影响不大。表 1 和表 2 分别是仿真条件和结果。鉴于此，在下面的仿真中，高度参数均选为 0.2mm。2.22.2 计算原理计算原理 对于 MCCB 而言，动触头是在一个对转动轴的力矩M M的作用下打开，如图 2所示。 对任何一个单元i，其对于转轴 O 的转矩M Mi为d di和力密度F Fi的向量积，那么在整个动导电杆区域对M Mi进行体积分运算，则可得到作用在其上相对于 O的力矩。从而作用

12、在其上的等效电动力也可以得到，如式(2)所示。而力密度F Fi可通过式(3)计算，其中J Ji和B Bi分别为单元i上的电流密度和磁通密度。文献7通过理论上的计算，指出 MCCB 中涡流几乎不影响电动斥力的数值和相位,这样可以采用恒定场的方程来计算电流密度和磁通密度的分布。在导电体区域，即触头导电回路，电流密度J J满足式(4)和式(5)所示的边界条件。其中s为导体的电导率，在本文中，导电杆和触头分别为铜和银材料；T T为矢量电位，I I为流过导体的电流。 得到了电流密度J J的分布后，在整个场域中，根据磁通密度B B和J J之间的关系式(6)，其中A A为矢量磁位，m 为磁导率，即可得到B

13、B的分布。基于以上的电流-磁场-电动斥力之间的关系，采用三维有限元分析，可以得到作用在动导电杆上和触头上的电动斥力。由于该力和短路电流之间存在单调增的关系，不考虑铁磁物质的影响时，电动斥力与短路电流的平方成正比。而且在特定的短路条件下，短路电流和时间有一定的关系，本文采用振荡回路作为实验电路，那么在触头斥开之前，电路方程为式(7)，这样就可以通过对上述过程的迭代处理，当电动斥力F F等于触头预压力F FK时，迭代结束，此时对应的时间即为触头斥开的时间。式中 I为电流第一半波的有效值，频率为工频 50Hz。3 3 分析模型分析模型对两种额定电流均为 100A 的塑壳断路器产品 CB1、CB2 进

14、行了仿真。其中CB1 静导电杆采用水平面 U 型，而 CB2 采用下进线 U 型。另外，在不改变产品其他形状及尺寸的条件下，将 CB1 的静导电杆改为平板式；将 CB2 的栅片腿长缩短，并加入 U 型电机槽；去掉 CB1 模型的灭弧栅片。研究电动斥力的变化，分别称其为 CB3、CB4 和 CB5。图 3(a)、(b)、(c)、(d)分别为 CB1、CB2、CB3 和 CB4 在 ANSYS 中的剖分图。由于其均具有对称性，因而在仿真过程中对其一半进行分析。其中(a)、(c)右下角为静触头附近区域的静导电杆的俯视及电流流向示意图。CB1 导电回路、动导电杆和动触头以及 CB3 的导电回路电流分布

15、分别如图 4(a)、(b)、(c)所示。和图 3(a)中所示的电流流向相一致，从图 4(a)所示的 CB1 俯视图可以清晰看到由于其静导电回路的 U 型弯曲导致图中 1 所示部分的电流方向和动导电杆相反。另外本文所有分析都是相对于该图所示的坐标系进行的，其中x、y及z轴正方向符合右手法则，在该图中，y轴正方向垂直纸面向外。从图 4(b)可以看出，触头上由于电流收缩，其电流密度最大值比动导杆上高 5 个数量级。而图 4 (c)所示的 CB3 导电回路上的电流分布和图 4 (a)的不同之处就在图中 1 所示部分的电流方向和动导电杆相同。图 5 所示为模型 CB1 和 CB3 动导电杆上的单元z方向

16、上平均磁通密度Bz和y方向上平均电动斥力密度fy沿x轴的分布情况，其原点位于图 2 中的 A 点。可以发现，靠近触头的区域磁场和力密度均比较大，而且该区域远离转轴，力臂较大，因而可以认为，回路电动斥力主要取决于这部分。同时明显看到 CB1的磁场和力密度较 CB3 大。 在统一计算触头间的 Holm 力和动导电杆上的 Lorentz 力的数值分析过程中，可以认为作用在动导电杆上各单元的力归算到其上 B 点（见图 2）的等效合力为回路力，而作用在动触头上各单元的力归算到 B 点的等效合力为触头力，即Holm 力。从表 3 所示的对 CB1 的计算数据，电动斥力和电流的平方近似成正比例关系。当电流为

17、 10kA 时对 CB1 到 CB5 五个模型的电动斥力进行了详细的计算，结果如表 4 所示，其中的导电桥半径是根据式(1)，按照实际产品 CB1 和 CB2 的设计参数确定的。结合表 3 可以看出，不同模型在不同的电流下由触头间电流收缩产生的电动斥力均占总体的 70%以上。对产品 CB1，改变其静导电回路的形状后电动斥力减小，这就表明水平面内的 U 型回路通过改变电流方向可以加强触头区域的磁场，进而增大电动斥力。而不考虑栅片的作用时，回路上的电动斥力减小了 40%左右，说明铁磁物质可以有效加强动导电杆上的磁场。对产品 CB2 ，改变栅片尺寸并加入 U 型电机槽后，电动斥力也有所增大。同时，这种改进也有利于将产气材料固定在电机槽内侧，利用新的气吹灭弧原理提高 MCCB 的开断性能。4 4 实验方法及结果分析实验方法及结果分析 4.14.1 实验方法实验方法 实验