工程电磁场 第六章电磁场的边值问题.pptVIP

* 如果问题的边界面与直角坐标系的坐标面吻合，则可采用直角坐标系中的分离变量法。 下面通过例子具体说明该方法。 例 求如图所示二维长方形内的电位函数。 解：根据题意，所求区域的电位函数满足的方程及边界条件为 x a y b * 只与x有关 只与y有关 在直角坐标系中方程 可写为 （二维问题，与z无关） ★ 分离变量法的前提即假设待求函数有分离变量形式的解： 上式两端同除以 因此该式成立的条件： 且 ● 为实数 为虚数 为虚数 为实数 ● 为零 为零 ● 为实数 为虚数 为零 或 * 同样的讨论适用于函数 。为满足x=0和x=a的边界条件，应选取 则 因为 将边界条件 将边界条件 于是 称为边值问题 的本征值。它的意义是：在上述边界条件下，分离常数 只有取这些特定值时，方程才有非零解。其解的函 数形式 称为本征函数。 * 对于 因为 将边界条件 于是 得 由于 故 的一般形式 将边界条件 这实际上是将一已知函数展为傅里叶级数。利用傅里叶级数的系数公式得 原问题的解 * 3 有 限 差 分 法 图 3.1 差分网格 3.1 差分表示式 * 二维泊松方程的差分格式 (Difference Form of 2D Poisson’s Equation) （1） 二维静电场边值问题 基本思想：将场域离散为许多网格 ，应用差分原理，将求解连续函数 的微分方程问题转换为求解网格节点上 的代数方程组的问题。 （2） 有限差分的网格分割 * 令 h = x - x0，将 x = x1 和 x3 分别代入式 ( 3 ) （3） 由式（4）+（5） （7） 同理, 沿 x方向在 x0 处的泰勒公式展开为 下 页 上 页 返 回 （4） （5） （6） * 将式(6)、式(7)代入式(1)，得到 当场域中 即 即 五点差分格式 下 页 上 页 返 回 * 上式表明，任一点的电位等于它周围四个点电位的平均值。显然，当h越小，计算越精确。如果待求N个点的电位，就需解含有N个方程的线性方程组。若点的数目较多，用迭代法较为方便。 * 矩形网格剖分 若场域离散为矩形网格，差分格式为 * 3.2 边界条件离散化（Discrete Boundary Condition) 第二类边界条件 第一类边界条件 分界面衔接条件 对称边界条件 其中 介质分界面 对称分界 * 3.3差分方程的数值解法 1. 简单迭代法 图 3.2 节点序号 * 2. 塞德尔(Seidel)迭代法 通常为节约计算时间，对简单迭代法要进行改进，每当算出一个节点的高一次的近似值，就立即用它参与其它节点的差分方程迭代，这种迭代法叫做塞德尔(Seidel)迭代法。塞德尔迭代法的表达式为 此式也称为异步迭代法。由于更新值的提前使用，异步迭代法比简单迭代法收敛速度加快一倍左右，存储量也小。 * 3. 超松驰迭代法 式中α称为松弛因子，其值介于 1 和 2 之间。当其值为1时，超松弛迭代法就蜕变为塞德尔(Seidel)迭代法。 因子α的选取一般只能依经验进行。但是对矩形区域， 当M、 N都很大时，可以由如下公式计算最佳收敛因子α0: * 边界节点赋已知电位值 赋节点电位初始值 累计迭代次数 N=0 N=N+1 按超松弛法进行一次迭代，求 打印 N Y 程序框图 * 例 设如图 所示的矩形截面的长导体槽，宽为 4h，高为3h，顶板与两侧绝缘，顶板的电位为 10V，其余的电位为零，求槽内各点的电位。 * 解：将待求的区域分为 12 个边长为h的正方形网格， 含六个内点，得出差分方程组： * 解以上方程组， 得 * 简 单 迭 代 法 1 2 3 4 5 6 0 0.0 0.0 0.0 0.0 0.0 0.0 1 2.5 0.0 2.5 0.0 2.5 0.0 2 3.125 0.625 3.75 0.625 3.125 0.625 … … … … … … … 10 4.1435 1.5363 5.0698 2.0242 4.1435 1.5363 11 4.1515 1.5419 5.0778 0.0356 4.1515 1.5419 * 超松弛迭代法(α=1.2) * 松弛因子的影响 * * * * * * * * 本章结束 * * 第六章 电磁场的边值问题 * 一、麦克斯韦方程组 * 1、①四个方程