《模拟电荷法（2020年8月整理）.pptx》由会员分享，可在线阅读，更多相关《模拟电荷法（2020年8月整理）.pptx（8页珍藏版）》请在金锄头文库上搜索。
1、一 寸 光 阴 不 可 轻 模拟电荷法简介 模拟电荷法是静电场数值计算的主要方法之一。类似于镜像法，模拟电荷法基于静 电场的惟一性定理，将导体电极表面连续分布的自由电荷用位于导体内部的一组离散的 电荷来替代（例如在导体内部设置一组点电荷、线电荷或环电荷等），这些离散的电荷 被称为模拟电荷。然后应用叠加定理，用这些模拟电荷的解析公式计算场域中任意一点 的电位或电场强度。而这些模拟电荷则根据场域的边界条件确定，模拟电荷法的关键在 于寻找和确定模拟电荷。 早在 1950 年，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)是已知的。,1,.,2n n,一 寸 光 阴 不 可 轻 根据叠加定理，对应于所有的 n 个匹配点，可以逐一地列出由设定的 n 个模拟电荷所建 立的电位表达式，即 P11Q1 P12Q2 . P1nQn 01,P Q P Q . P Q , 21 122 2,02 (1),Pn1Q1 Pn2Q2 . PnnQn 0n 式中， P
5、ij 为第 j 个单位模拟电荷在第 i 个匹配点上产生的电位值，它与第 j 个模拟电荷 的类型和位置以及第 i 个匹配点和场域的介电常数有关，而与第 j 个模拟电荷的电荷值 无关，通常称为电位系数。 模拟电荷方程组的求解 针对具体问题列出模拟电荷的线性代数方程组（1）以后，就可以用高斯列主元消 去法求解模拟电荷的电荷值 Qj (j=1, 2, , n)。 模拟电荷的校验与确定 然而，尚不能用这些模拟电荷直接计算场域任意一点的电位或电场强度，必须检验 这些模拟电荷产生的电位是否满足电极表面上非匹配点的边界条件。为此在电极表面上 另外选取 m 个电位校验点，电位校验点一般选在两个相邻的电位匹配点中
6、间。用模拟电 荷对每个电位校验点的电位进行计算，即 k Pk1Q1 Pk 2Q2 . PknQn(k 1, 2, ., m) （2） 将各个电位校验点的电位值k (k=1, 2, , m)与给定的已知电位值0k (k=1, 2, , m)相比， 如果二者的差值满足 k 0k (k 1, 2, ., m) （3） 式中为由计算者预先确定的计算误差，则由式（1）解出的模拟电荷的电荷值有效，即 可以应用这些模拟电荷来计算场域中任意一点的电位或电场强度。如果式（3）不满足， 则应该对第一次假设的模拟电荷的类型、位置、个数做适当调整，重新计算式（1）和 式（2），直到式（3）满足要求为止。 综上所述，模
7、拟电荷法的主要步骤如下： 根据对电极和场域的定性分析和经验，确定一组模拟电荷的类型、位置和数 量； 根据电极表面的几何形状，选定与模拟电荷数量相同的电极表面电位匹配点，,2,3,一 寸 光 阴 不 可 轻 然后建立模拟电荷的线性代数方程组； 解模拟电荷的线性代数方程组，求解模拟电荷的电荷值； 在电极表面另外选定足够数量的电位校验点，校验电极表面的电位计算精度 。如果不符合要求，则重新修正模拟电荷的类型、位置、或数量，再行计算，直到达到 所要求的计算精度为止。一般经过几次修正即能达到要求； 按所得的模拟电荷用解析计算公式计算场域内任意一点的电位或电场强度。 模拟电荷与匹配点的选择 在模拟电荷法中
8、，模拟电荷的选定是相当随意的。这时计算者的直观经验的判断将 起到重要作用。计算者一般按场分布的特点，选择模拟电荷的类型、位置和数量，以求 使所设置的模拟电荷的整体在给定电极表面的总电位有可能满足电极表面的等电位要 求。 模拟电荷与电极表面匹配点的布置对于计算精度有很大的影响。在具体应用时，通 常由于匹配点位于电极表面，所以应首先选定电极表面匹配点的位置，然后再决定相应 模拟电荷的位置。经验表明，遵循以下原则是合宜的： 在电场急剧变化处或在所关心的场域附近，电极表面匹配点和模拟电荷可以 分布密些，但在电极表面的角点处一般不应设置匹配点。 模拟电荷宜正对电极表面匹配点放置，并以落在与电极表面内的垂
9、线上为佳 。此外，若设模拟电荷到电极表面的垂直距离为 a，对应电极表面匹配点左右相邻两个 匹配点之间的距离为 b，则根据经验，二者之间的比值 a/b 取 0.2-1.5，通常取 0.75。 模拟电荷的设定不是越多越好，因为在数值计算中，数值解的误差不仅与离 散误差相关，还与电位系数矩阵的条件数有关。当电极表面匹配点的数量增多时，虽然 离散误差似有所减小，但是同时也导致电位系数矩阵中相邻两行或相邻两列之间的元素 数值相近，因而使电位系数矩阵的行列式的值很小，也就是说，使电位系数矩阵的条件 数增大，即导致所谓“病态”线性代数方程组。这将引起计算时的舍入误差和数字有效 位相减误差的增加，影响计算精度
10、。所以，需要综合考虑。 6常用的模拟电荷及其解析计算公式 常用的模拟电荷类型有点电荷、直线电荷和圆环电荷。这些电荷既可以应用于三维 电场计算也可以应用于二维电场计算。下面，只讨论二维电场计算中应用的典型电荷。 （1）二维平行平面电场中的模拟电荷,一 寸 光 阴 不 可 轻 在二维平行平面电场中，惟一可以采用的模拟电荷是图 1 所示的无限长直线电荷。 若设 y=0 为电位参考点，则该对线电荷在任意点 A 处产生的电位为,22,11,ln,(x x 2 ) ( y y ) 22 2, 2 ,2 ln 4(x x ) ( y y ) 1,（4）,对应的电位系数为,1,1,1(x x )2 ( y y
11、 )2,4(x x )2 ( y y )2,P ln 22（5）,电场强度为,22,1,2,x,2,E ( x x 1 x x 2 ) （6）,22,1,2,y,2,E ( y y 1 y y 2 ) （7）,图 1 一对等值异号的无限长直线电荷 （2）二维轴对称电场中的模拟电荷 在二维轴对称电场中，可以采用的模拟电荷有点电荷、直线电荷和圆环电荷，相关 的电位、电位系数和电场强度解析计算公式如下： 点电荷 图 2 所示的点电荷 q 在任意点 A 处产生的电位为,qq1,4 r4, , 2 (z z )2,（8）,对应的电位系数为,1,4,1 4,P , 2 (z z )2,（9）,电场强度为,
12、一 寸 光 阴 不 可 轻,4,3 2 (z z )2 2,E q（10）,z,q 4,z z ,E ,3 2 (z z )2 2,（11）,图 2 点电荷 直线电荷 图 3 所示长度为 l 的直线电荷在任意点 A 处产生的电位为,4(z zl) 1, ln (z z0 ) 1 （12）,22,2,2,101l,式中， (z z ) ， (z z ) 。,对应的电位系数为,1,4(z zl) 1,P ln (z z0 ) 1 （13）,电场强度为,1,1,1,4,E ,( z z0 z zl ) （14）,11,1,z,4 ,E ( 1 1 ),（15）,图 3 直线电荷 圆环电荷 图 4 所
13、示为半径为 R、电荷量为 q 的圆环电荷，在任意点 A 处产生的电位为,5,一 寸 光 阴 不 可 轻,1,q,2 2, ,K (k)（16）,式中,1 ( R) (z z ) 22,k 2, R 1,K(k)为第一类完全椭圆积分，即,0, 2,d,K (k ) ,1 k 2 sin2 ,（17）,对应的电位系数为,1,1,2 2,P ,K (k)（18）,电场强度为,2,1,1,q,E(k),4 2 , ,E , 1 K (k) ,（19）,22,1,z,2 ,qz z E(k),E ,（20）,式中, ( R)2 (z z )2 R2 2 (z z )2 E(k)为第二类完全椭圆积分，即,
14、2,0,E(k) 1 k 2 sin2 d,（21）,图 4 圆环电荷,6,一 寸 光 阴 不 可 轻 7模拟电荷法的典型算例 （1）球-板电极的电场【参考倪光正工程电磁场原理（第二版）第 98 页例 2-16】 图 5 所示为一个球-板电极结构，导体球的半径为 1m，球心到接地导体板的距离为 5m，球上电位为 100V。应用模拟电荷法计算电场分布，并求球-板之间的电容值。,图 5 球-板电场 应用模拟电荷法求解，使用两个模拟电荷 q1 和 q2，将 q2 放置在球心上，即 y2=5m； q1 放置在导体球的球心对接地导体板镜像位置时对导体球内的反演点上，即 y1=5-1/10=4.9m。同时
15、，为了考虑接地导体板的情况，再考虑两个模拟电荷 q1 和 q2 的镜 像电荷。电位匹配点分别选在导体球的下端点 A 和上端点 B。通过求解模拟电荷的线性 代数方程组，得 q 1.238109C ，q 1.111108C 12,100,7,取 C 点为电位校验点，其电位计算值为 99.9V。导体球表面电位的相对误差为 0.1%。 满足工程计算要求。球-板之间的电容值为 C q1 q2 1.235 nF,8优化模拟电荷法 通过上面的讨论，不难看出模拟电荷法简单、直观，能较有效地满足计算工程中的 电场问题。但是，电场的计算精度主要取决于计算者的经验。为了克服这个不足，后来 在模拟电荷法的基础上，又发
16、展出来了优化模拟电荷法。 优化模拟电荷法是模拟电荷法与优化算法相结合的成果，在确定模拟电荷的类型和 数量后，通过优化算法，自动调整和改变模拟电荷的位置和电荷值来满足电极表面电位 匹配点的边界条件。下面，通过一个棒-板电极的电场计算为例，简述优化模拟电荷法 的基本过程。 图6 所示为棒-板电极系统，圆棒导体的半径为R0.1m，棒的端部为半径为R0.1m,一 寸 光 阴 不 可 轻 的半球导体，距接地导体板的高度为 G0.5m。设导体圆棒的电位为 1V。,2,m,1 m,i1,f x ,0ici,( ) (22),式中，0i 是圆棒导体表面第 i 个电位匹配点上的已知电位值，ci 是圆棒导体表面第 i 个 电位匹配点上用模拟电荷计算的电位值。点电荷和半无限长线电荷的位置坐标满足如下 的约束条件（以确保优化过程中，模拟电荷始终在圆棒导体内部）
