模拟电荷法.docVIP

模拟电荷法简介 模拟电荷法是静电场数值计算的主要方法之一。类似于镜像法，模拟电荷法基于静电场的惟一性定理，将导体电极表面连续分布的自由电荷用位于导体内部的一组离散的电荷来替代（例如在导体内部设置一组点电荷、线电荷或环电荷等），这些离散的电荷被称为模拟电荷。然后应用叠加定理，用这些模拟电荷的解析公式计算场域中任意一点的电位或电场强度。而这些模拟电荷则根据场域的边界条件确定，模拟电荷法的关键在于寻找和确定模拟电荷。 早在1950年，Loeb在研究棒形电极的电晕放电时，就用了一组虚设在电极内部的电荷来计算棒形电极对地的电场分布，他没有用计算机仅通过手算就完成了计算。到了20世纪60年代末期，M. S. Abou-Seada和E. Nasser借助计算机用模拟电荷法计算了棒形电极和圆柱形电极对地的电场分布。20世纪70年代以后，模拟电荷法在高电压工程的电场计算方面显示了很大的优点。首先H. Singer和H. Steinbigler将模拟电荷法用于二维和三维电场计算，获得了满意的计算结果。之后B. Bachmann将模拟电荷法用于计算具有表面漏电流的电场问题。后来模拟电荷法也被用于高电压电极系统的优化设计，例如通过修正原有电极的外形使电极表面的电场强度均匀分布或使电极表面的最大电场强度最小，以使绝缘材料得到充分利用。 从本质上看，模拟电荷法可以看作是广义的镜像法，但它在数值处理和工程实用方面远优于镜像法。 1．模拟电荷法的基本思想 模拟电荷法的基本思想就是在被求解的场域以外，用一组虚设的模拟电荷来等效代替电极表面的连续分布的电荷，并应用这些模拟电荷的电位或电场强度的解析计算公式来计算电场。模拟电荷的类型（例如点电荷、直线电荷、圆环电荷等）和位置是由计算者事先根据电极的形状和对场分布的定性分析所假设的，而模拟电荷的电荷值则由电极的电位边界条件通过解线性代数组确定。当模拟电荷的电荷值确定后，场中任意一点的电位或电场强度就可以通过叠加定理由这些模拟电荷的电位或电场强度的解析计算公式进行计算。因此，模拟电荷法的理论基础是静电场的惟一性定理。 2．模拟电荷方程组的建立 在场域外（电极内部）设置n个模拟电荷Qj (j=1, 2, …, n)，在给定电位边界条件的电极表面上选定n个电位匹配点，显然，各匹配点上的电位值(0j (j=1, 2, …, n)是已知的。根据叠加定理，对应于所有的n个匹配点，可以逐一地列出由设定的n个模拟电荷所建立的电位表达式，即 (1) 式中， Pij为第j个单位模拟电荷在第i个匹配点上产生的电位值，它与第j个模拟电荷的类型和位置以及第i个匹配点和场域的介电常数有关，而与第j个模拟电荷的电荷值无关，通常称为电位系数。 3．模拟电荷方程组的求解 针对具体问题列出模拟电荷的线性代数方程组（1）以后，就可以用高斯列主元消去法求解模拟电荷的电荷值Qj (j=1, 2, …, n)。 4．模拟电荷的校验与确定 然而，尚不能用这些模拟电荷直接计算场域任意一点的电位或电场强度，必须检验这些模拟电荷产生的电位是否满足电极表面上非匹配点的边界条件。为此在电极表面上另外选取m个电位校验点，电位校验点一般选在两个相邻的电位匹配点中间。用模拟电荷对每个电位校验点的电位进行计算，即 （2） 将各个电位校验点的电位值(k (k=1, 2, …, m)与给定的已知电位值(0k (k=1, 2, …, m)相比，如果二者的差值满足 （3） 式中(为由计算者预先确定的计算误差，则由式（1）解出的模拟电荷的电荷值有效，即可以应用这些模拟电荷来计算场域中任意一点的电位或电场强度。如果式（3）不满足，则应该对第一次假设的模拟电荷的类型、位置、个数做适当调整，重新计算式（1）和式（2），直到式（3）满足要求为止。 综上所述，模拟电荷法的主要步骤如下： （1）根据对电极和场域的定性分析和经验，确定一组模拟电荷的类型、位置和数量； （2）根据电极表面的几何形状，选定与模拟电荷数量相同的电极表面电位匹配点，然后建立模拟电荷的线性代数方程组； （3）解模拟电荷的线性代数方程组，求解模拟电荷的电荷值； （4）在电极表面另外选定足够数量的电位校验点，校验电极表面的电位计算精度。如果不符合要求，则重新修正模拟电荷的类型、位置、或数量，再行计算，直到达到所要求的计算精度为止。一般经过几次修正即能达到要求； （5）按所得的模拟电荷用解析计算公式计算场域内任意一点的电位或电场强度。 5．模拟电荷与匹配点的选择 在模拟电荷法中，模拟电荷的选定是相当随意的。这时计算者的直观经验的判断将起到重要作用。计算者一般按场分布的特点，选择模拟电荷的类型、位置和数量，以求使所设置的模拟电荷的整