基于模拟电荷法的电晕放电后绝缘材料表面电势的计算方法
【专利摘要】本发明是一种基于模拟电荷法的电晕放电后绝缘材料表面电势的计算方法，其特点是，包括建立坐标系、在电晕放电针电极设置模拟电荷及模拟电荷相对应的匹配点、根据叠加原理构建模拟电荷方程组和根据模拟电荷电荷量的大小计算绝缘材料表面的电势分布等步骤。能够对针?板电晕放电后绝缘材料的表面电势进行有效的计算，电晕放电后绝缘材料表面电势计算的准确性高，计算难度低，计算量少，适用于不同形状的电极。具有方法科学，适用性强，效果佳等优点。
【专利说明】
基于模拟电荷法的电晕放电后绝缘材料表面电势的计算方法
技术领域
[0001] 本发明涉及一种表面电势计算方法，尤其涉及一种基于模拟电荷法的电晕放电后 绝缘材料表面电势的计算方法。
【背景技术】
[0002] 随着我国经济持续发展，人民生活水平不断提高，电力需求快速增长，电网规模与 容量不断增大。为提高线路输送能力且保证输电安全，发展更高性能的绝缘材料已成为了 必然趋势。
[0003] 高压线路在输电过程中的安全性会受到绝缘材料性能的影响，电晕放电会引起绝 缘材料的老化甚至失效，绝缘的老化或失效会引发事故，因此采用必要手段对绝缘材料的 性能进行测试是避免事故发生的重要手段之一。目前，采用有限元方法对电晕放电后绝缘 材料表面电势进行计算无法对无限空间进行精确计算，并且在计算针-板电晕放电时，对极 细的针尖附近区域进行剖分会引起起算误差，不仅计算量较大，也无法保证计算的精确程 度。
[0004] 因此，需要采用模拟电荷法对电晕放电后绝缘材料表面电势进行计算，通过调整 模拟电荷的类型及数量，能够有效的进行计算，保证计算精度的同时，还能够降低计算量， 而且适用性广。
【发明内容】
[0005] 本发明的目的是，提供的一种能够对针-板电晕放电后绝缘材料的表面电势进行 有效的计算，电晕放电后绝缘材料表面电势计算的准确性高，计算难度低，计算量少，适用 于不同形状的电极，适用性强的基于模拟电荷法的电晕放电后绝缘材料表面电势的计算方 法。
[0006] 本发明的目的是由以下技术方案来实现的：一种基于模拟电荷法的电晕放电后绝 缘材料表面电势的计算方法，其特征在于，包括如下步骤：
[0007] (a)建立坐标系，并规定X轴为针-板电晕放电的板电极上表面，y轴为针电极轴心；
[0008] (b)在电晕放电针电极设置模拟电荷及模拟电荷相对应的匹配点；
[0009] (c)根据叠加原理，构建模拟电荷方程组，并根据模拟电荷的方程组求解模拟电荷 电荷量的大小；
[0010] (d.)根据模拟电荷电荷量的大小计算绝缘材料表面的电势分布。
[0011] 进一步，步骤(c)和步骤(d)之间还包括步骤(e):根据模拟电荷电荷量的大小计算 边界电势误差，若边界电势误差不满足误差范围要求，再根据电位误差对模拟电荷的位置 及个数进行重置，直到边界电位误差满足误差范围。
[0012] 进一步，步骤(C)中，模拟电荷方程组为：
[0013]
[0014] 其中：约…A为针电极尖端部位电势;％+ι…％,+m为针电极非尖端部位电势;Qi···Q n 为针电极尖端的模拟电荷;Qn+r · · Qn+m为针电极非尖端部位的模拟电荷，当I < i <η，I < j <η +m时，Pij为针电极尖端部位的模拟电荷及镜像电荷产生的电位系数：
，针 电极尖端部分采用的模拟电荷类型为点电荷，
[0015] 其中= I1S模拟电荷与其对应匹配点的距离，I2为镜像电荷与其对应匹配点的距 离；
[0016]当n+Ki<n+m, K j<n+m时，Pij为针电极非尖端部位的模拟电荷及镜像电荷产 生的电位系数
，针电极非尖端部分采用的模拟电 荷类型为线电荷，
[0017] 其中：1:为模拟线电荷上端与其对应匹配点的距离，I2为模拟限电荷下端与其对应 匹配点的距离，In为镜像线电荷上端与其对应匹配点的距离，I 22为镜像线电荷下端与其对 应匹配点的距离，Z1为模拟线电荷上端坐标，Z2为模拟线电荷下端坐标，Z为与模拟电荷对应 的匹配点坐标；
[0018] 对模拟电荷方程组求逆，得出模拟电荷的大小为
[0019] 进一步，步骤(e)中，根据优化模拟电荷法对模拟电荷的位置及数量进行设置:通 过求解电极边界上所有匹配点的已知电位和计算电位之差的平方值之和的最小值来获得 达到误差要求的模拟电荷
讲莫拟电荷j在第i个
匹配点处的电位；％为第i个匹配点的已知电位，若误差值不能满足所需的精确度要求，按 照l<k<2优化模拟电荷的位置及数量，其中k为两个连续的匹配点的距离与匹配点与其对 应模拟电荷之间距离的比值。
[0020] 进一步，步骤(d)中，根据不同介质间的边界条件 ，计算 出材料表面模拟电荷的大小，
[0021]其中：当1彡i彡η，1彡j彡n+m时，fij为针电极尖端部位的模拟电荷及镜像电荷产生 的电场系数
当n+Ki彡n+m, 1彡j彡n+m时，fij为针电极非尖端部位的 模拟电荷及镜像电荷产生的电位系数，
[0022] 进一步，步骤(d)中：根据达到误差允许范围的所有模拟电荷计算电晕放电后绝缘 材料表面电势分布:设任一点为P，根据叠加原理，其电势为
[0023] 其中Q为模拟电荷的电荷量，R为对应模拟电荷与点P之间的距离。
[0024] 本发明的基于模拟电荷法的电晕放电后绝缘材料表面电势的计算方法的有益效 果是:能够对针-板电晕放电后绝缘材料表面的电势分布进行计算，有效地克服了有限元方 法中无法对无限空间进行精确计算以及对极细的针电极尖端进行剖分时引起的计算误差， 提高了电晕放电后绝缘材料表面电势计算的准确性，而且降低了计算难度，减少计算量，适 用于不同形状的电极。具有方法科学，适用性强，效果佳等优点。
【附图说明】
[0025] 下面结合附图和实施例对本发明作进一步描述：
[0026] 图1为本发明的流程图。
[0027] 图2为本发明计算电晕放电后绝缘材料表面电势时模拟电荷以及匹配点设置的示 意图。
[0028] 图3为本发明的计算结果与测试结果对比图。
[0029] 图4为本发明的计算结果与有限元算法计算误差对比图。
【具体实施方式】
[0030] 参照图1和图2,本发明的基于模拟电荷法的电晕放电后绝缘材料表面电势的计算 方法，包括如下步骤：
[0031] (a)建立坐标系，并规定X轴为针-板电晕放电的板电极上表面，y轴为针电极轴心；
[0032] (b)在电晕放电针电极设置模拟电荷及模拟电荷相对应的匹配点；
[0033] (c)根据叠加原理，构建模拟电荷方程组，并根据模拟电荷的方程组求解模拟电荷 电荷量的大小；
[0034] (d)根据模拟电荷电荷量的大小计算绝缘材料表面的电势分布。
[0035] 本实施例中，步骤(C)和步骤(d)之间还包括步骤(e):根据模拟电荷的大小计算边 界电势误差，若边界电势误差不满足误差范围要求，再根据电位误差对模拟电荷的位置及 个数进行重置，然后再执行步骤(C)，直到边界电位误差满足误差范围后执行步骤(d)。
[0036] 本实施例中，步骤(c)中，模拟电荷方程组为：
[0037]
[0038] 其中：約…仍,为针电极尖端部位电势；Ρ,μΓ · · % _为针电极非尖端部位电势;Q1…Qn 为针电极尖端的模拟电荷;Qn+r ·_Qn+m为针电极非尖端部位的模拟电荷，当K i<n，I< j<n +m时，Plj为针电极尖端部位的模拟电荷及镜像电荷产生的电位系数
针 电极尖端部分采用的模拟电荷类型为点电荷，
[0039] 其中= I1S模拟电荷与其对应匹配点的距离，I2为镜像电荷与其对应匹配点的距 离；
[0040] 当n+Ki彡n+m，l彡j彡n+m时，Pij为针电极非尖端部位的模拟电荷及镜像电荷产 生的电位系数，
I，针电极非尖端部分采用的模拟电 荷类型为线电荷，
[0041] 其中：1:为模拟线电荷上端与其对应匹配点的距离，I2为模拟限电荷下端与其对应 匹配点的距离，In为镜像线电荷上端与其对应匹配点的距离，I 22为镜像线电荷下端与其对 应匹配点的距离，Z1为模拟线电荷上端坐标，Z2为模拟线电荷下端坐标，Z为与模拟电荷对应 的匹配点坐标；
[0042] 对模拟电荷方程组求逆，得出模拟电荷的大小为
[0043]本实施例中，步骤（e)中，根据优化模拟电荷法对模拟电荷的位置及数量进行设 置:通过求解电极边界上所有匹配点的已知电位和计算电位之差的平方值之和的最小值来 获得达到误差要求的模拟电荷
为模拟电荷j在 第i个匹配点处的电位;也为第i个匹配点的已知电位。
[0044]若误差值不能满足所需的精确度要求，按照l<k<2优化模拟电荷的位置及数量， 其中k为两个连续的匹配点的距离与匹配点与其对应电荷之间距离的比值。
[0045] 本实施例中，步骤(d)中，根据不同介质间的边界条 计算出材料表面模拟电荷的大小，
[0046] 其中：当1彡i彡n，l彡j彡n+m时，fi j为针电极尖端部位的模拟电荷及镜像电荷产生 的电场系数
当n+Ki彡n+m, 1彡j彡n+m时，fij为针电极非尖端部位的 模拟电荷及镜像电荷产生的电位系数，
[0047] 根据达到误差允许范围的所有模拟电荷计算电晕放电后绝缘材料表面电势分布： 设任一点为P，根据叠加原理，其电契
[0048] 其中Q为模拟电荷的电荷量，R为对应模拟电荷与点P之间的距离。
[0049]以实际计算-IOkV电晕放电后聚酰亚胺薄膜表面电势为例，并与测试结果进行对 比。该实验采用针-板电晕放电方式，针电极为负极，板极为正极，负极电晕电压为-10kV，针 板间距为3cm，在针电极尖端部位选取10个模拟电荷，在针电极非尖端部位选取15个模拟电 荷进行计算，计算结果和测试结果对比如图3所示，图中实线为根据本发明的电势计算方法 得出的绝缘材料表面电势，虚线为实验测得数据。同时，实验测得数据、有限元法计算得到 的绝缘材料表面电势和本发明计算得出的表面电势曲线如图4所示，从图4中可以得出，依 据本发明所提供的方法计算得出的表面电势数值更接近于测量实际值，由此证明本发明的 计算方法有效，且精确度高。
[0050] 综上所述的实施例证明，本发明的基于模拟电荷法的电晕放电后绝缘材料表面电 势的计算方法明显优于有限元法。
[0051] 本发明的实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对 本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行 修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利 要求范围当中。
【主权项】
1. 一种基于模拟电荷法的电晕放电后绝缘材料表面电势的计算方法，其特征在于:包 括如下步骤： (a) 建立坐标系，并规定x轴为针-板电晕放电的板电极上表面，y轴为针电极轴心； (b) 在电晕放电针电极设置模拟电荷及模拟电荷相对应的匹配点； (c) 根据叠加原理，构建模拟电荷方程组，并根据模拟电荷的方程组求解模拟电荷电荷 量的大小； (d) 根据模拟电荷电荷量的大小计算绝缘材料表面的电势分布。2. 根据权利要求1所述基于模拟电荷法的电晕放电后绝缘材料表面电势的计算方法， 其特征在于:步骤(c)和步骤(d)之间还包括步骤(e):根据模拟电荷的电荷量大小计算边界 电势误差，若边界电势误差不满足误差范围要求，再根据电位误差对模拟电荷的位置及个 数进行重置，直到边界电位误差满足误差范围。3. 根据权利要求1所述基于模拟电荷法的电晕放电后绝缘材料表面电势的计算方法， 其特征在于:步骤(c)中，模拟电荷的方程组为：其中：約…为针电极尖端部位电势;…乳+?为针电极非尖端部位电势;Ql…Qn为针 电极尖端的模拟电荷;Qn+1…Qn+m为针电极非尖端部位的模拟电荷，当Ki<n，K j<n+m 时，为针电极尖端部位的模拟电荷及镜像电荷产生的电位系数，，针电 极尖端部分采用的模拟电荷类型为点电荷，其中：1:为模拟电荷与其对应匹配点的距离，12为镜像电荷与其对应匹配点的距离； 当n+Ki<n+m，K j<n+m时，Pij为针电极非尖端部位的模拟电荷及镜像电荷产生的 电位系数：?针电极非尖端部分采用的模拟电荷类 型为线电荷， 其中：1:为模拟线电荷上端与其对应匹配点的距离，h为模拟限电荷下端与其对应匹配 点的距离，In为镜像线电荷上端与其对应匹配点的距离，122为镜像线电荷下端与其对应匹 配点的距离，Zi为模拟线电荷上端坐标，Z 2为模拟线电荷下端坐标，Z为与模拟电荷对应的匹 配点坐标； 对模拟电荷方程组求逆，得出模拟电荷的大小为旧卜!；/5] 1[纠。4. 根据权利要求2所述基于模拟电荷法的电晕放电后绝缘材料表面电势的计算方法， 其特征在于:步骤(e)中，根据优化模拟电荷法对模拟电荷的位置及数量进行设置:通过求 解电极边界上所有匹配点的已知电位和计算电位之差的平方和的最小值来获得达到误差 要求的模拟电荷：为模拟电荷j在第i个匹配点处 的电位;I为第i个匹配点的已知电位， 若误差值不能满足所需的精确度要求，按照l<k<2优化模拟电荷的位置及数量，其中 k为两个连续的匹配点的距离与匹配点与其对应模拟电荷之间距离的比值。5. 根据权利要求1所述基于模拟电荷法的电晕放电后绝缘材料表面电势的计算方法， 其特征在于:步骤(d)中，根据不同介质间的边界条件：|计算出材 料表面模拟电荷的大小， 其中：当1 < n，1 < j < n+m时，f i j为针电极尖端部位的模拟电荷及镜像电荷产生的电 场系数：，当n+Ki彡n+m，1彡j彡n+m时，fij为针电极非尖端部位的模拟 电荷及镜像电荷产生的电位系数，6. 根据权利要求1或5所述基于模拟电荷法的电晕放电后绝缘材料表面电势的计算方 法，其特征在于:步骤(d)中：根据达到误差允许范围的所有模拟电荷计算电晕放电后绝缘 材料表面电势分布:设任一点为P，根据叠加原理，其电势为其中:Q为模拟电荷的电荷量，R为对应模拟电荷与点P之间的距离。
【文档编号】G06F17/50GK106055811SQ201610398213
【公开日】2016年10月26日
【申请日】2016年6月7日
【发明人】张嘉伟, 陈骁, 刘苪彤, 高凤凯
【申请人】东北电力大学