工程电磁场计算方法及精度分析

模拟电荷法是静电场数值计算的主要方法之一。类似于镜像法，模拟电荷法基于静电场的唯一性定理，将导体电极表面连续分布的自由电荷用位于导体内部的一组离散的电荷来替代（例如在导体内部设置一组点电荷、线电荷或环电荷等），这些离散的电荷被称为模拟电荷。然后应用叠加定理，用这些模拟电荷的解析公式计算场域中任意一点的电位或电场强度。而这些模拟电荷则根据场域的边界条件确定，模拟电荷法的关键在于寻找和确定模拟电荷。

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

从本质上看，模拟电荷法可以看作是广义的镜像法，但它在数值处理和工程实用方面远优于镜像法。

模拟电荷法的基本思想就是在被求解的场域以外，用一组虚设的模拟电荷来等效代替电极表面的连续分布的电荷，并应用这些模拟电荷的电位或电场强度的解析计算公式来计算电场。模拟电荷的类型（例如点电荷、直线电荷、圆环电荷等）和位置是由计算者事先根据电极的形状和对场分布的定性分析所假设的，而模拟电荷的电荷值则由电极的电位边界条件通过解线性代数组确定。当模拟电荷的电荷值确定后，场中任意一点的电位或电场强度就可以通过叠加定理由这些模拟电荷的电位或电场强度的解析计算公式进行计算。因此，模拟电荷法的理论基础是静电场的唯一性定理。

在场域外（电极内部）设置n 个模拟电荷Qj（j=1，2，…，n），在给定电位边界条件的电极表面上选定n 个电位匹配点，显然，各匹配点上的电位值φ0j（j=1，2，…，n）是已知的。根据叠加定理，对应于所有的n 个匹配点，可以逐一地列出由设定的n 个模拟电荷所建立的电位表达式，即

式中，Pij 为第j 个单位模拟电荷在第i 个匹配点上产生的电位值，它与第j 个模拟电荷的类型和位置以及第i 个匹配点和场域的介电常数有关，而与第j 个模拟电荷的电荷值无关，通常称为电位系数。

针对具体问题列出模拟电荷的线性代数方程组（1-12）以后，就可以求解模拟电荷的电荷值Qj（j=1，2，…，n）。

然而，尚不能用这些模拟电荷直接计算场域任意一点的电位或电场强度，必须检验这些模拟电荷产生的电位是否满足电极表面上非匹配点的边界条件。为此在电极表面上另外选取m 个电位校验点，电位校验点一般选在两个相邻的电位匹配点中间。用模拟电荷对每个电位校验点的电位进行计算，即

将各个电位校验点的电位值φk（k=1，2，…，m）与给定的已知电位值φ0k（k=1，2，…，m）相比，如果二者的差值满足

式中Δ 为由计算者预先确定的计算误差，则由式（1-12）解出的模拟电荷的电荷值有效，即可以应用这些模拟电荷来计算场域中任意一点的电位或电场强度。如果式（1-14）不满足，则应该对第一次假设的模拟电荷的类型、位置、个数做适当调整，重新计算式（1-12）和式（1-13），直到式（1-14）满足要求为止。

综上所述，模拟电荷法的主要步骤如下：(www.xing528.com)

（1）根据对电极和场域的定性分析和经验，确定一组模拟电荷的类型、位置和数量。

（2）根据电极表面的几何形状，选定与模拟电荷数量相同的电极表面电位匹配点，然后建立模拟电荷的线性代数方程组。

（3）解模拟电荷的线性代数方程组，求解模拟电荷的电荷值。

（4）在电极表面另外选定足够数量的电位校验点，校验电极表面的电位计算精度。如果不符合要求，则重新修正模拟电荷的类型、位置或数量，再行计算，直到达到所要求的计算精度为止。一般经过几次修正即能达到要求。

（5）按所得的模拟电荷用解析计算公式计算场域内任意一点的电位或电场强度。

在模拟电荷法中，模拟电荷的选定是相当随意的。这时计算者的直观经验的判断将起到重要作用。计算者一般按场分布的特点，选择模拟电荷的类型、位置和数量，以求使所设置的模拟电荷的整体在给定电极表面的总电位有可能满足电极表面的等电位要求。

模拟电荷与电极表面匹配点的布置对于计算精度有很大的影响。在具体应用时，通常由于匹配点位于电极表面，所以应首先选定电极表面匹配点的位置，然后再决定相应模拟电荷的位置。经验表明，遵循以下原则是合宜的：

（1）在电场急剧变化处或在所关心的场域附近，电极表面匹配点和模拟电荷可以分布密些，但在电极表面的角点处一般不应设置匹配点。

（2）模拟电荷宜正对电极表面匹配点放置，并以落在与电极表面内的垂线上为佳。此外，若设模拟电荷到电极表面的垂直距离为a，对应电极表面匹配点左右相邻两个匹配点之间的距离为b，则根据经验，二者之间的比值a/b 取0.2～1.5，通常取0.75。

（3）模拟电荷的设定不是越多越好，因为在数值计算中，数值解的误差不仅与离散误差相关，还与电位系数矩阵的条件数有关。当电极表面匹配点的数量增多时，虽然离散误差似有所减小，但是同时也导致电位系数矩阵中相邻两行或相邻两列之间的元素数值相近，因而使电位系数矩阵的行列式的值很小，电位系数矩阵的条件数增大，即导致所谓“病态”线性代数方程组。这将引起计算时的舍入误差和数字有效位相减误差的增加，影响计算精度。