工程电磁场模糊电荷法

1、模拟电荷法石晶洋实践电 1201 班 1121180618模拟电荷法是静电场数值计算的主要方法之一。类似于镜像法，模拟电荷法基于静电场的惟一性定理，将导体电极表面连续分布的自由电荷用位于导体内部的一组离散的电荷来替代（例如在导体内部设置一组点电荷、线电荷或环电荷等），这些离散的电荷被称为模拟电荷。然后应用叠加定理，用这些模拟电荷的解读公式计算场域中任意一点的电位或电场强度。而这些模拟电荷则根据场域的边界条件确定，模拟电荷法的关键在于寻找和确定模拟电荷。早在 1950 年， L. B. Loeb 在研究棒形电极的电晕放电时，就用了一组虚设在电极内部的电荷来计算棒形电极对地的电场分布，他没有用计算

2、机仅通过手算就完成了计算。到了20 世纪 60 年代末期， M. S. Abou-Seada和 E. Nasser借助计算机用模拟电荷法计算了棒形电极和圆柱形电极对地的电场分布。20 世纪70 年代以后，模拟电荷法在高电压工程的电场计算方面显示了很大的优点。首先 H. Singer 和 H. Steinbigler 将模拟电荷法用于二维和三维电场计算，获得了满意的计算结果。之后 B. Bachmann 将模拟电荷法用于计算具有表面漏电流的电场问题。后来模拟电荷法也被用于高电压电极系统的优化设计，例如通过修正原有电极的外形使电极表面的电场强度均匀分布或使电极表面的最大电场强度最小，以使绝缘材料得

3、到充分利用。从本质上看，模拟电荷法可以看作是广义的镜像法，但它在数值处理和工程实用方面远优于镜像法。1模拟电荷法的基本思想模拟电荷法的基本思想就是在被求解的场域以外，用一组虚设的模拟电荷来等效代替电极表面的连续分布的电荷，并应用这些模拟电荷的电位或电场强度的解读计算公式来计算电场。模拟电荷的类型（例如点电荷、直线电荷、圆环电荷等）和位置是由计算者事先根据电极的形状和对场分布的定性分析所假设的，而模拟电荷的电荷值则由电极的电位边界条件通过解线性代数组确定。当模拟电荷的电荷值确定后，场中任意一点的电位或电场强度就可以通过叠加定理由这些模拟电荷的电位或电场强度的解读计算公式进行计算。因此，模拟电荷法

4、的理论基础是静电场的惟一性定理。2模拟电荷方程组的建立在场域外（电极内部）设置n 个模拟电荷Qj (j=1, 2, n)，在给定电位边界条件的电极表面上选定n 个电位匹配点，显然，各匹配点上的电位值0j (j=1,2, , n)是已知的。根据叠加定理，对应于所有的 n 个匹配点，可以逐一地列出由设定的 n 个模拟电荷所建立的电位表达式，即(1)式中， Pij 为第 j 个单位模拟电荷在第 i 个匹配点上产生的电位值，它与第拟电荷的类型和位置以及第 i 个匹配点和场域的介电常数有关，而与第电荷的电荷值无关，通常称为电位系数。j 个模j 个模拟3模拟电荷方程组的求解针对具体问题列出模拟电荷的线性代

5、数方程组（ 1）以后，就可以用高斯列主元消去法求解模拟电荷的电荷值 Qj(j=1, 2, , n)。4模拟电荷的校验与确定然而，尚不能用这些模拟电荷直接计算场域任意一点的电位或电场强度，必须检验这些模拟电荷产生的电位是否满足电极表面上非匹配点的边界条件。为此在电极表面上另外选取 m 个电位校验点，电位校验点一般选在两个相邻的电位匹配点中间。用模拟电荷对每个电位校验点的电位进行计算，即（ 2）将各个电位校验点的电位值k (k=1, 2, m)与给定的已知电位值 0k (k=1, 2, ,m)相比，如果二者的差值满足（ 3）式中 为由计算者预先确定的计算误差，则由式（ 1）解出的模拟电荷的电荷值有

6、效，即可以应用这些模拟电荷来计算场域中任意一点的电位或电场强度。如果式（ 3）不满足，则应该对第一次假设的模拟电荷的类型、位置、个数做适当调整，重新计算式（ 1）和式（ 2），直到式（ 3）满足要求为止。综上所述，模拟电荷法的主要步骤如下：（1）根据对电极和场域的定性分析和经验，确定一组模拟电荷的类型、位置和数量；（2）根据电极表面的几何形状，选定与模拟电荷数量相同的电极表面电位匹配点，然后建立模拟电荷的线性代数方程组；（3）解模拟电荷的线性代数方程组，求解模拟电荷的电荷值；（4）在电极表面另外选定足够数量的电位校验点，校验电极表面的电位计算精度。如果不符合要求，则重新修正模拟电荷的类型、位置

7、、或数量，再行计算，直到达到所要求的计算精度为止。一般经过几次修正即能达到要求；（5）按所得的模拟电荷用解读计算公式计算场域内任意一点的电位或电场强度。5模拟电荷与匹配点的选择在模拟电荷法中，模拟电荷的选定是相当随意的。这时计算者的直观经验的判断将起到重要作用。计算者一般按场分布的特点，选择模拟电荷的类型、位置和数量，以求使所设置的模拟电荷的整体在给定电极表面的总电位有可能满足电极表面的等电位要求。模拟电荷与电极表面匹配点的布置对于计算精度有很大的影响。在具体应用时，通常由于匹配点位于电极表面，所以应首先选定电极表面匹配点的位置，然后再决定相应模拟电荷的位置。经验表明，遵循以下原则是合宜的：（

8、1）在电场急剧变化处或在所关心的场域附近，电极表面匹配点和模拟电荷可以分布密些，但在电极表面的角点处一般不应设置匹配点。（2）模拟电荷宜正对电极表面匹配点放置，并以落在与电极表面内的垂线上为佳。此外，若设模拟电荷到电极表面的垂直距离为 a，对应电极表面匹配点左右相邻两个匹配点之间的距离为 b，则根据经验，二者之间的比值 a/b 取 0.2-1.5 ，通常取 0.75 。（3）模拟电荷的设定不是越多越好，因为在数值计算中，数值解的误差不仅与离散误差相关，还与电位系数矩阵的条件数有关。当电极表面匹配点的数量增多时，虽然离散误差似有所减小，但是同时也导致电位系数矩阵中相邻两行或相邻两列之间的元素数值

9、相近，因而使电位系数矩阵的行列式的值很小，也就是说，使电位系数矩阵的条件数增大，即导致所谓“病态”线性代数方程组。这将引起计算时的舍入误差和数字有效位相减误差的增加，影响计算精度。所以，需要综合考虑。（1）球 - 板电极的电场图5 所示为一个球 - 板电极结构，导体球的半径为h，球心到接地导体板的距离为a，球上电位为100V。应用模拟电荷法计算电场分布图 5 球-板电场应用模拟电荷法求解，使用两个模拟电荷q1 和q2，将q2 放置在球心上，即y2=5m； q1 放置在导体球的球心对接地导体板镜像位置时对导体球内的反演点上，即 y1=h-a*a/(2*h)。同时，为了考虑接地导体板的情况，再考虑

10、两个模拟电荷 q1 和 q2 的镜像电荷。电位匹配点分别选在导体球的下端点A 和上端点 B。通过求解模拟电荷的线性代数方程组取 C 点为电位校验点，得到其电位计算值为U。导体球表面电位的相对误差为 ( 1 U/100)*100% 。看是否满足工程计算要求。#include <stdio.h>#include <math.h>#define PI 3.1416#define E 8.854e-12void main()float a,b,h 。double A2,B2,U1,U2,Ua,m,n,s,c。long double q1,q2 。printf(" 请输

11、入电极距地高度scanf("%f%f",&h,&a)。printf("h=%f,a=%fn",h,a)。:h球半径:an") 。b=a*a/(2*h)。printf("b=%fn",b)。printf("请输入匹配点的电位:U1 U2")。scanf("%lf%lf",&U1,&U2)。printf("U1=%f,U2=%fn",U1,U2)。c=4*PI*E 。c=1/c 。A0=c*(1/a-1/(2*h+a)。A1=c*(1/(a+b)-1/(2*h+a-b)。B0=c*(1/a-1/(2*h-a)。B1=c*(1/(a-b)-1/(2*h-a-b)。q1=(U1*B1-U2*A1)/(A0*B1-A1*B0)。q2=(U1*B0-U2*A0)/(A1*B0-B1*A0)。printf(" 模拟电荷的值为:q1=%e,q2