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

1、低压塑壳断路器中电动斥力的三维有限元非线性分析与实验研究            低压塑壳断路器中电动斥力的 三维有限元非线性分析与实验研究 李兴文，陈德桂，向洪岗，李志鹏，刘洪武 （西安交通大学电气工程学院，陕西 西安710049） 3-D Finite Element Nonlinear Analysis and Experimental Investigation of Electro-dynamic Repulsion Force in Molded Case Circ

2、uit BreakersLi Xing-wen, Chen De-gui, Xiang Hong-gang, Li Zhi-peng, Liu Hong-wu （School of Electrical Engineering, Xian Jiaotong University, Xian 710049, China）Abstract: For the optimization design of molded case circuit breakers (MCCBs), it is necessary and important to calculate the electro-dynami

3、c repulsion force acting on the movable conductor. With 3-D finite element nonlinear analysis, according to the equations among current-magnetic field-repulsion force and taking into account the ferromagnet, we introduce contact bridge model to simulate the current constriction between contacts, so

4、Lorentz and Holm force acting on the movable conductor and contact, respectively, can be integrated to calculate. Coupled with circuit equations, the opening time of movable contact also can be obtained using iteration with the restriction of contact force. Simulation and experiment for repulsion fo

5、rce and opening time of five different configuration models have been investigated. The results indicate that the proposed method is effective and is capable of evaluating new design of contact systems in MCCBs.           Keywords: Low-voltage apparatus; finit

6、e element analy-sis; repulsion force 摘要：电动斥力的计算对低压塑壳断路器（MCCB）的设计有重要意义。基于电流-磁场-电动斥力之间的方程，并考虑铁磁物质的影响，应用三维有限元非线性分析，引入圆柱导电桥模型作为接触点模拟触头间的电流收缩，统一计算触头间的Holm力和动导电杆上的Lorentz力。耦合电路方程，并以动触头上的预压力为约束通过迭代运算，从而确定动触头的打开时间。分析了五种不同结构的MCCB，并对电动斥力和触头打开时间进行了实验研究，结果表明，该方法是有效的，可用于MCCB新产品触头系统的设计。 关键词：低压电器；有限元；电动斥力 1 引言 当短路

7、电流产生时，低压塑壳断路器（MCCB）中动、静触头在电动斥力的作用下分离，在操作机构的带动下，电弧将在触头之间产生，拉长，并在吹弧力的作用下进入栅片被分割成若干个短弧，而后熄灭。在这个过程的起始阶段，作用在动触头上的电动斥力F（包括导电回路产生的Lorentz力FL及触头间由于电流收缩产生的Holm力FH）和预压力决定了触头的斥开时间和打开速度，从而对MCCB的限流性能产生重要的影响。而且，从式(1)所示的Holm公式可以看出，接触点半径r与预压力FK、触头材料的布氏硬度H、触头表面接触情况（用x描述，其范围一般为0.30.6，通常取0.45)有关。而Holm力FH与r，触头半径R、以及电流大

8、小i有关1，这样FK对F就有一定的影响；另一方面，FK的选择也由于发热容许的要求，而受到F的限制。因此，在进行MCCB的触头导电回路的设计时必须进行电动斥力的计算。近年来，随着计算机技术发展和数值计算方法的不断进步，计算机辅助工程在MCCB新产品的设计中得到了越来越广泛的应用2-3，大大缩短了产品开发周期，降低了开发成本，同时产品的性能也得到了很大提高。对于电动斥力的计算，国内外的学者做了很多有益的工作。在文献4-5中，由于Holm公式已经被证明在计算孤立的触头间电动斥力时是正确的6，因而用式(1)来计算FH ，而忽略触头区域的电流收缩。这种简化不仅使载流导体中的电流分布和实际相差甚远，而且会

9、进一步影响了磁场的分布，从而使FL的计算结果有一定的误差。 本文应用三维有限元分析，引入圆柱体导电桥模型作为接触点，以五种不同结构的MCCB产品为对象，着重研究了电动斥力和触头斥开时间的计算方法，分析了不同结构的触头系统对于电动斥力的影响。并和实验结果进行了比较，表明仿真所用的方法是正确的，可以应用于MCCB新产品的触头系统的开发。同时，计算结果也可作为进一步应用虚拟样机技术进行MCCB开断过程仿真时的输入数据。2 计算方法      论文低压塑壳断路器中电动斥2.1 导电桥模型 应用数值计算方法综合考虑触头间的电流收缩和导电回路对于作用在动导电杆上电动

10、斥力的影响时，必须首先引入一个合理的计算模型来描述触头间的电接触情况。R.holm在推导式(1)中电动斥力FH解析式时，为了分析上的需要，假定接触导体为超导小球1。本文为了和实际情况更加接近，用位于触头中心的圆柱体导电桥模型来模拟导电斑点，其材料性质也和触头材料相同，半径r可由式(1)所示的Holm公式计算。为了确定导电桥高度参数h，对一对圆柱体触头进行了电动斥力仿真，图1所示为其截面的示意图。结果发现，在相同的r下，其高度h在0.10.25mm范围内对电动斥力的影响不大。表1和表2分别是仿真条件和结果。鉴于此，在下面的仿真中，高度参数均选为0.2mm。2.2 计算原理 对于MCCB而言，动触

11、头是在一个对转动轴的力矩M的作用下打开，如图2所示。 对任何一个单元i，其对于转轴O的转矩Mi为di和力密度Fi的向量积，那么在整个动导电杆区域对Mi进行体积分运算，则可得到作用在其上相对于O的力矩。从而作用在其上的等效电动力也可以得到，如式(2)所示。而力密度Fi可通过式(3)计算，其中Ji和Bi分别为单元i上的电流密度和磁通密度。文献7通过理论上的计算，指出MCCB中涡流几乎不影响电动斥力的数值和相位, 这样可以采用恒定场的方程来计算电流密度和磁通密度的分布。在导电体区域，即触头导电回路，电流密度J满足式(4)和式(5)所示的边界条件。其中s为导体的电导率，在本文中，导电杆和触头

12、分别为铜和银材料；T为矢量电位，I为流过导体的电流。  得到了电流密度J的分布后，在整个场域中，根据磁通密度B和J之间的关系式(6)，其中A为矢量磁位，m 为磁导率，即可得到B的分布。基于以上的电流-磁场-电动斥力之间的关系，采用三维有限元分析，可以得到作用在动导电杆上和触头上的电动斥力。由于该力和短路电流之间存在单调增的关系，不考虑铁磁物质的影响时，电动斥力与短路电流的平方成正比。而且在特定的短路条件下，短路电流和时间有一定的关系，本文采用振荡回路作为实验电路，那么在触头斥开之前，电路方程为式(7)，这样就可以通过对上述过程的迭代处理，当电动斥力F等于触头预压力FK时，迭代结束，此

13、时对应的时间即为触头斥开的时间。式中 I为电流第一半波的有效值，频率为工频50Hz。3 分析模型对两种额定电流均为100A的塑壳断路器产品CB1、CB2进行了仿真。其中CB1静导电杆采用水平面U型，而CB2采用下进线U型。另外，在不改变产品其他形状及尺寸的条件下，将CB1的静导电杆改为平板式；将CB2的栅片腿长缩短，并加入U型电机槽；去掉CB1模型的灭弧栅片。研究电动斥力的变化，分别称其为CB3、CB4和CB4在ANSYS中的剖分图。由于其均具有对称性，因而在仿真过程中对其一半进行分析。其中(a)、(c)右下角为静触头附近区域的静导电杆的俯视及电流流向示意图。CB1导电回路、动导电杆和动触头以

14、及CB3的导电回路电流分布分别如图4(a)、(b)、(c)所示。和图3(a)中所示的电流流向相一致，从图4(a)所示的CB1俯视图可以清晰看到由于其静导电回路的U型弯曲导致图中1所示部分的电流方向和动导电杆相反。另外本文所有分析都是相对于该图所示的坐标系进行的，其中x、y及z轴正方向符合右手法则，在该图中，y轴正方向垂直纸面向外。从图4(b)可以看出，触头上由于电流收缩，其电流密度最大值比动导杆上高5个数量级。而图4 (c)所示的CB3导电回路上的电流分布和图4 (a)的不同之处就在图中1所示部分的电流方向和动导电杆相当电流为10kA时对CB1到CB5五个模型的电动斥力进行了详细的计算，结果如表4所示，其中的导电桥半径是根据式(1)，按照实际产品CB1和CB2的设计参数确定的。结合表3可以看出，不同模型在不同的电流下由触头间电流收缩产生的电动斥力均占总体的70%以上。对产品CB1，改变其静导电回路的形状后电动斥力减小，这就表明水平