低压塑壳断路器电动斥力.doc

低压塑壳断路器中电动斥力的低压塑壳断路器中电动斥力的 三维有限元非线性分析与实验研究三维有限元非线性分析与实验研究 李兴文 陈德桂 向洪岗 李志鹏 刘洪武李兴文 陈德桂 向洪岗 李志鹏 刘洪武 西安交通大学电气工程学院 陕西 西安交通大学电气工程学院 陕西 西安西安 710049710049 3 D3 D FINITEFINITE ELEMENTELEMENT NONLINEARNONLINEAR ANALYSISANALYSIS ANDAND EXPERIMENTALEXPERIMENTAL INVESTIGATIONINVESTIGATION OFOF ELECTRO ELECTRO DYNAMICDYNAMIC REPULSIONREPULSION FORCEFORCE ININ MoldedMolded CaseCase CircuitCircuit BreakersBreakers LI Xing wen CHEN De gui XIANG Hong gang LI Zhi peng LIU Hong wu School of Electrical Engineering Xi an Jiaotong University Xi an 710049 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 magnetic 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 equations 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 effective 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 力 耦合电路方程 并以动触头 上的预压力为约束通过迭代运算 从而确定动触头的打开时间 分析了五种不 同结构的 MCCB 并对电动斥力和触头打开时间进行了实验研究 结果表明 该 方法是有效的 可用于 MCCB 新产品触头系统的设计 关键词 关键词 低压电器 有限元 电动斥力 1 1 引言引言 当短路电流产生时 低压塑壳断路器 MCCB 中动 静触头在电动斥力的 作用下分离 在操作机构的带动下 电弧将在触头之间产生 拉长 并在吹弧 力的作用下进入栅片被分割成若干个短弧 而后熄灭 在这个过程的起始阶段 作用在动触头上的电动斥力F F 包括导电回路产生的 Lorentz 力F FL及触头间由 于电流收缩产生的 Holm 力F FH 和预压力决定了触头的斥开时间和打开速度 从 而对 MCCB 的限流性能产生重要的影响 而且 从式 1 所示的 Holm 公式可以看 出 接触点半径r与预压力F FK 触头材料的布氏硬度H 触头表面接触情况 用x描述 其范围一般为 0 3 0 6 通常取 0 45 有关 而 Holm 力F FH与r 触头半径R 以及电流大小i有关 1 这样F FK对F F就有一定的影响 另一方面 F FK的选择也由于发热容许的要求 而受到F F的限制 因此 在进行 MCCB 的触头 导电回路的设计时必须进行电动斥力的计算 近年来 随着计算机技术发展和数值计算方法的不断进步 计算机辅助工 程在 MCCB 新产品的设计中得到了越来越广泛的应用 2 3 大大缩短了产品开发 周期 降低了开发成本 同时产品的性能也得到了很大提高 对于电动斥力的 计算 国内外的学者做了很多有益的工作 在文献 4 5 中 由于 Holm 公式已 经被证明在计算孤立的触头间电动斥力时是正确的 6 因而用式 1 来计算F FH 而忽略触头区域的电流收缩 这种简化不仅使载流导体中的电流分布和实际 相差甚远 而且会进一步影响了磁场的分布 从而使F FL的计算结果有一定的 误差 本文应用三维有限元分析 引入圆柱体导电桥模型作为接触点 以五种不 同结构的 MCCB 产品为对象 着重研究了电动斥力和触头斥开时间的计算方法 分析了不同结构的触头系统对于电动斥力的影响 并和实验结果进行了比较 表明仿真所用的方法是正确的 可以应用于 MCCB 新产品的触头系统的开发 同 时 计算结果也可作为进一步应用虚拟样机技术进行 MCCB 开断过程仿真时的输 入数据 2 2 计算方法计算方法 2 12 1 导电桥模型导电桥模型 应用数值计算方法综合考虑触头间的电流收缩和导电回路对于作用在动导 电杆上电动斥力的影响时 必须首先引入一个合理的计算模型来描述触头间的 电接触情况 R holm 在推导式 1 中电动斥力FH解析式时 为了分析上的需要 假定接 触导体为超导小球 1 本文为了和实际情况更加接近 用位于触头中心的圆柱 体导电桥模型来模拟导电斑点 其材料性质也和触头材料相同 半径r可由式 1 所示的 Holm 公式计算 为了确定导电桥高度参数h 对一对圆柱体触头进 行了电动斥力仿真 图 1 所示为其截面的示意图 结果发现 在相同的r下 其高度h在 0 1 0 25mm 范围内对电动斥力的影响不大 表 1 和表 2 分别是仿真 条件和结果 鉴于此 在下面的仿真中 高度参数均选为 0 2mm 2 22 2 计算原理计算原理 对于 MCCB 而言 动触头是在一个对转动轴的力矩M M的作用下打开 如图 2 所示 对任何一个单元i 其对于转轴 O 的转矩M Mi为d di和力密度F Fi的向量积 那么在整个动导电杆区域对M Mi进行体积分运算 则可得到作用在其上相对于 O 的力矩 从而作用在其上的等效电动力也可以得到 如式 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 B的分布 基于以上的电流 磁场 电动斥力之间的关系 采用三维有限元分析 可以 得到作用在动导电杆上和触头上的电动斥力 由于该力和短路电流之间存在单 调增的关系 不考虑铁磁物质的影响时 电动斥力与短路电流的平方成正比 而且在特定的短路条件下 短路电流和时间有一定的关系 本文采用振荡回路 作为实验电路 那么在触头斥开之前 电路方程为式 7 这样就可以通过对上 述过程的迭代处理 当电动斥力F F等于触头预压力F FK时 迭代结束 此时对应的 时间即为触头斥开的时间 式中 I为电流第一半波的有效值 频率为工频 50Hz 3 3 分析模型分析模型 对两种额定电流均为 100A 的塑壳断路器产品 CB1 CB2 进行了仿真 其中 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 的导电回路电流分布分别如图 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方向上平均磁通密度Bz和 y方向上平均电动斥力密度fy沿x轴的分布情况 其原点位于图 2 中的 A 点 可以发现 靠近触头的区域磁场和力密度均比较大 而且该区域远离转轴 力 臂较大 因而可以认为 回路电动斥力主要取决于这部分 同时明显看到 CB1 的磁场和力密度较 CB3 大 在统一计算触头间的 Holm 力和动导电杆上的 Lorentz 力的数值分析过程中 可以认为作用在动导电杆上各单元的力归算到其上 B 点 见图 2 的等效合力 为回路力 而作用在动触头上各单元的力归算到 B 点的等效合力为触头力 即 Holm 力 从表 3 所示的对 CB1 的计算数据 电动斥力和电流的平方近似成正比例关 系 当电流为 10kA 时对 CB1 到 CB5 五个模型的电动斥力进行了详细的计算 结 果如表 4 所示 其中的导电桥半径是根据式 1 按照实际产品 CB1 和 CB2 的设 计参数确定的 结合表 3 可以看出 不同模型在不同的电流下由触头间电流收 缩产生的电动斥力均占总体的 70 以上 对产品 CB1 改变其静导电回路的形状 后电动斥力减小 这就表明水平面内的 U 型回路通过改变电流方向可以加强触 头区域的磁场 进而增大电动斥力 而不考虑栅片的作用时 回路上的电动斥 力减小了 40 左右 说明铁磁物质可以有效加强动导电杆上的磁场 对产品 CB2 改变栅片尺寸并加入 U 型电机槽后 电动斥力也有所增大 同时 这种改进 也有利于将产气材料固定在电机槽内侧 利用新的气吹灭弧原理提高 MCCB 的开 断性能 4 4 实验方法及结果分析实验方法及结果分析 4 14 1 实验方法实验方法 实验采用单频高压振荡回路模拟短路电流 图 6 为合成回路的原理图 其 中 SP 为实验样机 Ci和Cu分别为电流和电压回路充电电容 Li和Lu分别为电 流和电压回路的振荡电感 H 是一个经过严格标定的无感采样电阻 其阻值为 90mW 误差范围为 0 2 级 K1为主合闸开关 在本实验中 仅需用电流回路就可完成 其振荡频率为 50Hz 第一半波充 电电压有效值和放电电流第一半波有效值比是 46 1 V kA 也就是说 如果要 得到的放电电流第一半波有效值为 10kA 那么电容器组所需的充电电压为 460V 实验开始前 K1打开 通过整流器对电容器组Ci充电 当达到实验所需 的电压后 断开充电电路 即可开始实验 闭合 K1 Ci Li K1 K3和 H 就构成 了一个典型的单频振荡回路 采用示波器和高压探头分别测得电弧电流和电弧 电压波形 采用石英拉压传感器进行电动斥力的测量 当短路电流流过断路器时 在 电动斥力的作用下 动触头试图打开并产生电弧 为了能够测量到电动斥力变 化的波形 且进行多次实验而不发生融焊 那么就必须压紧动导电杆使动触头 不被斥开 在实验中传感器安装在模型上 一个胶木杆一端通过螺纹固定在传 感器上 另一端压紧动导电杆 并可旋转它调节作用在动导电杆上的预压力大 小 4 24 2 电动斥力测量电动斥力测量 为了便于用传感器测量电动斥力 加工了一个将产品 CB1 放大 4 倍的模型 触头材料为铜 x可取 0 45 图 7 为短路电流峰值 9kA 时 作用在动导电杆上 的电动斥力和短路电流波形 其中预压力FK为 25 1N 其具体数据如表 5 所示 从数据中可以看出 峰值电流下的电动斥力计算值和测量所得的力波形峰 值很接近 而且短路电流和电动斥力之间的相位差也很小 考虑到拉压传感器 及胶木杆本身响应时间引起的滞后作用 那么涡流对于电动斥力的影响就很小 了 因此上述引入导电桥 采用静态的分析方法计算电动斥力是合适的 也是 比较准确的 而且 MCCB 一般在几千安电流下作用在动导电杆上的电动斥力就可 以克服触头预压力而使动触头斥开 产品 CB2 的实验也说明了这个问题 这样 也可以采用静态方法计算触头斥开时间而忽略涡流的影响 4 34 3 触头斥开时间实验触头斥开时间实验 以 CB2 断路器产品为研究对象 在预期短路电流有效值为 10kA 条件下进行 开断实验 图 8 为测量到的电弧电流和电弧电压波形 可以明显地看出 当分 流器 H 上的电压上升到 273mV 时 也就是短路电流达到 273mV 90mW 3033A 时 间为 759ms 时 动静触头之间的电压有一个突跳 说明触头此刻在电动斥力的 作用下打开 电弧产生 此时刻前 电流为正弦波 而且由于触头间接触电阻 的存在 电压随着电流上升有一个很小的数值 此时刻后 由于电弧电压的存 在及其限流作用 电流波形不再是正弦波 当电流过零后 电弧熄灭 CB2 的触头预压力为 4 4N 预期短路电流有效值为 10kA 通过迭代 仿真 结果为电流达到 3020A 时间为 674ms 电动斥力达到 4 4N 触头开始分离 和实验结果比较可以发现 由于动触头转动轴摩擦等的影响 实验值略大于计 算值 计算和实验结果表明 上述方法是可行的 根据不同短路情况下电流和时 间的关系 可以计算出触头斥开时间 5 5 结论结论 1 基于电流 磁场 电动斥力之间的关系 应用三维有限元非线性分析 引入导电桥模型来表述触头间的电流收缩 可统一计算作用在动导电杆上和触 头上的电动斥力 再耦合短路电路方程 以触头预压力为约束 经过迭代过程 可计算出触头斥开时间 和实验结果比较表明 仿真方法及结果是正确的 可 以用于 MCCB 新产品的触头导电回路设计 2 电动斥力和电流的平方成正比例关系 3 总的电动斥力主要由动触头及其附近区域的电动斥力提供 而动触头 上由于电流收缩所致的电动斥力占整个电动斥力的 70 以上 4 灭弧栅片对于回路电动斥力有较大影响 5 平面 U 型静导电回路较平板式时的整个电动斥力大 短栅片加入 U 型 电机槽较长栅片时电动斥力大 结构上更容易加入产气材料 从而利用新的气 吹灭弧技术可提高 MCCB 开断性能 6 仿真方法和结果可以进一步用于 MCCB 整个开断过程的仿真 参考文献参考文献 1 Holm R Electric Contacts J 4th Ed Springer Verlag 1967 2 陈旭 陈德桂 Chen Xu Chen Degui 低压限流断路器背后击穿现 象的数值模拟研究 Research on simulation of the back commutation in low voltage current limiting circuit breakers J 中国电机工 程学报 Proceedings of the CSEE 2000 20 3 16 19 3 杨武 荣命哲 王小华 等 Yang Wu Rong Mingzhe Wang Xiaohua et al 考虑电动力效应的高压断路器动力学特性仿真分析 Dynamic simulation of high voltage circuit breaker concerning electrodynamic force J 中国电机工程学报 Proceedings of the CSEE 2003 23 5 103 107 4 Stammberger H Force calculation for the movable contact of circuit break