二维静电场有限元分析软件包SEF2D2.docVIP

电磁场数值计算实验指导书 电信教研室 华北电力大学 二〇〇一年八月 目录 1 有限差分法……………………………………………………………1 2 矩量法……………………………………………………………6 3 有限单元法……………………………………………………………9 4 附录……………………………………………………………10 1 有限差分法 有限差分法（Finite Differential Method）是基于差分原理的一种数值计算法。其基本思想：将场域离散为许多小网格，应用差分原理，将求解连续函数的泊松方程的问题转换为求解网格节点上的差分方程组的问题。 二维泊松方程的差分格式 图1-1 有限差分的网格分割 二维静电场边值问题： (1-1) (1-2) 通常将场域分成足够小的正方形网格,网格线之间的距离为,节点上的电位分别用和表示。 设函数在处可微，则沿方向在处的泰勒公式展开为 （1-3） 将和分别代入式（1－3）,得 （1-4） （1-5） 由（1-4）-（1-5）得 （1-6） (1-4)+(1-5)得 （1-7） 同理 （1-8） （1-9） 当场域中得到拉普拉斯方程的五点差分格式 1.2 边界条件的离散化处理 图1-2 边界条件的离散化处理 若场域离散为矩形网格(如图2-2示)，差分格式为： (1-10) （1）第一类边界条件:给边界离散节点直接赋已知电位值 （2）对称边界条件:合理减小计算场域，差分格式为： (1-11) 图1-3 边界条件的离散化处理 （3）第二类边界条件:边界线与网格线相重合的差分格式： (1-12) (4)介质分界面衔接条件 的差分格式 (1-13) 其中 1.3 差分方程组的求解方法 (1) 高斯——赛德尔迭代法 (1-14) 式中： ( 迭代顺序可按先行后列，或先列后行进行。 ( 迭代过程遇到边界节点时，代入边界值或边界差分 格式，直到所有节点电位满足为止。 （2）超松弛迭代法 (1-15) 式中：——加速收敛因子 迭代收敛的速度与有明显关系 表1-1 迭代收敛的速度与的关系 收敛因子（） 1.0 1.7 1.8 1.83 1.85 1.87 1.90 2.0 迭代次数（） 1000 269 174 143 122 133 171 发散 最佳收敛因子的经验公式： （正方形场域、正方形网格） （矩形场域、正方形网格） ( 迭代收敛的速度与电位初始值的给定及网格剖分精细有关 ( 迭代收敛的速度与工程精度要求有关 借助计算机进行计算时，其程序框图1－5所示 图1-5 迭代解程序框图 练习题 试用超松弛迭代法求解接地金属槽内电位的分布 已知： 给定边值 如图1－6示 给定初值 误差范围 选取 计算：迭代次数N=? 分布。 图1-6接地金属槽的网格剖分 (2) 按对称场差分格式求解电位的分布 已知： 给定边值：如图1-7示 给定初值 误差范围: 图1-7 接地金属槽内半场域的网格剖分 计算：1）迭代次数，，将计算结果保存到文件中； 2）按电位差画出槽中等位线分布图。（选做，可以不做） 2 矩量法 矩量法（method of moment）在电磁场分析中有着广泛的应用。其概念相当简单，基本上是用未知场的积分方程去计算给定媒质中场的分布。 在静电学中，在由点的电荷分布在点产生的电位分布可以表示为 (2-1) 这里实质上是电位分布的源，是点和点间的距离。然而一般情况下是未知的，而源区电位的分布是给定的。因此，为了求出空间每个地方的电位分布，我们必须估计源区的电荷分布。 设的一个解是 (2-2) 这里是源区一些离散位置上预先选定的电荷分布，是待定未知系数，以式（2-2）代入式（2-1）得 (2-3) (2-4) 这里。所以考虑在位置的电荷，可以表示为下述电位的线性组合，即： (2-5) 所以 (2-6) 由于在源区是已知的，所以未知系数可以由