有限差分法在静态电磁场计算中的应用.pdfVIP

维普资讯 2000年 第 18卷 第 4期 长 春 邮 电 学 院 学 报 2000 Vo1．18 No．4 JOURNAL OF CHANGCHUN TJ0N INSTITUTE 文 章编 号 ：1000—1794 (2000)04—0027—04 有 限差分法在静态 电磁场 计算 中的应用 二) or／, 摘耍 ：利用有限差分法求解二维静 电场 问题 ，所得 的数值结果与解析结果相一致 ，分析讨论 了有限差分法 的精度和收敛性问题 。分析结果表明：当 取最佳值时收敛最快，步长越小精 度越高 。 关t词 塑里苎坌鎏 ；竖丝生 精度 存 中圈分类号 ：0 |41．4 文献标识码 ：A 引 言 在工程上碰到的电磁场 问题有许多是属于静态 电磁场计算 ，这就要求求得满足实际 条件的麦克斯韦方程的精确解答，然而只有一些典型几何形状和结构相对简单的问题才 有可能求得严格的解析解 。对于实际的工程 电磁场 问题，由于结构、边界复杂，要求它 的解析解并非易事 ，而且在绝大多数情况下都是十分 困难的。为此 ，发展了电磁场数值 计算方法，以弥补解析法的不足，有限差分法便是其 中的一种 。 l 有 限差分法原理 如图 1所示 ，在 由一边界 C界定的二维场域 D 内，电位 函数 满足如下边值 问题 一 耄+雾=一 =F (1) c= ，(户) (2) 其 中 P为体 电荷密度 ，￡为该点处 的介 电常数 。 对于给定 的边值问题 ，应用差分法 ，就是把所研究的区域 ，离散化为很多网格和节 点 ，如图 1所示。对任一节点用差商代替微商，从而把所研究 区域的偏微分方程变为各 个节点上的差分方程组，则求场的微分方程的解就变为求联立差分方程组解的问题 。利 用差商代替微商 ，二维泊松方程 (1)可近似表达为 鸭 + + + 一 4％ = h F (3) 兰墓品翟信zoo息o-产08业1部9中青年科学基垒项目(97。1Bl9) 作者筒 舟：何晓东 (1963-- )，女 ．吉林长春人 ．长春邮电学院剐教授 ，理学博士 ，主要从事光 电技术研究 维普资讯 28 长 春 邮 电 学 院 学 报 第 ls卷 式 中^为网格步长。当用网格将区域划分后，对场域 内每一网格节点可写出类似的式子。 为 了求解给定 的边值问题 ，除了对场域 内泊松方程进行差分离散化外 ，还必须对边 界条件进行差分离散 化处 理 ，以构成相应的差分边值 ／jlD ( )l13 问题 。若网格节点恰好落在 {lJf 30j2t／ 3 边界 C上 ，则只需直接把位 lfh I 函数 I ( 表征边 界节 l一 点)的值赋给对应 的边界节 O 图l 场域 的 网橹划分 固2 不等距 罔格 点 点；若 网格节点不落在边界 c上 ，如图 2所示 ，其 中 l、2为边界线上 的点，P、q为小于 1的正数 ，采用泰勒公式进 行差分离散化处理 ，可推得这些节点的泊松方程 的差分形式为 + 一．佟+南 传+南 一(古+吉j一