《工程电磁场第七章电磁场的边值问题.ppt》由会员分享，可在线阅读，更多相关《工程电磁场第七章电磁场的边值问题.ppt（48页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第六章第六章 电磁场的边值问题电磁场的边值问题 1 静电场和恒定电场的边值问题，可归结为在给定边静电场和恒定电场的边值问题，可归结为在给定边界条件下求解拉普拉斯方程或泊松方程。界条件下求解拉普拉斯方程或泊松方程。常用的方法常用的方法直接法直接法间接法间接法解析法解析法数值法数值法有限差分法有限差分法（FD）有限元方法有限元方法（FEM）矩量法矩量法（MoM）2静电参数(电容及部分电容)静电能量与力有限差分法镜像法,电轴法分离变量法直接积分法数值法解析法边值问题边界条件电位基本方程D 的散度基本物理量 E、D基本实验定律（库仑定律）E 的旋度31.镜像法镜像法 实质实质:是以一个或几个是以一个或
2、几个等效电荷等效电荷代替边界的影响，将原来具代替边界的影响，将原来具有边界的有边界的非均匀非均匀空间变成无限大的空间变成无限大的均匀均匀自由空间，从而使计算过自由空间，从而使计算过程大为简化。程大为简化。依据：依据：惟一性定理。因此，等效电荷的引入必须维持原来的惟一性定理。因此，等效电荷的引入必须维持原来的边界条件不变，从而保证原来区域中静电场没有改变，这是确定边界条件不变，从而保证原来区域中静电场没有改变，这是确定等效电荷的大小及其位置的依据。这些等效电荷通常处于等效电荷的大小及其位置的依据。这些等效电荷通常处于镜像位镜像位置置，因此称为，因此称为镜像电荷镜像电荷，而这种方法称为，而这种方法
3、称为镜像法镜像法。关键：关键：确定镜像电荷的大小及其位置。确定镜像电荷的大小及其位置。局限性：局限性：仅仅对于某些特殊的边界以及特殊分布的电荷才有仅仅对于某些特殊的边界以及特殊分布的电荷才有可能确定其镜像电荷。可能确定其镜像电荷。4 镜镜像像法法是是解解静静电电场场问问题题的的一一种种间间接接方方法法，它它巧巧妙妙地地应应用用唯唯一一性性定定理，理，使某些看来难解的边值问题容易地得到解决。使某些看来难解的边值问题容易地得到解决。使用镜像法时要注意以下三点：使用镜像法时要注意以下三点：（1 1）镜像电荷是虚拟电荷）镜像电荷是虚拟电荷;（2 2）镜像电荷置于所求区域之外的附近区域）镜像电荷置于所求
4、区域之外的附近区域;（3 3）导电体是等位面。）导电体是等位面。5（1）点电荷与无限大的导体平面。）点电荷与无限大的导体平面。介质 导体 q r P 介质 q r P hh 介质 以一个处于镜像位置的点电荷代替边界的影响，使整个空间以一个处于镜像位置的点电荷代替边界的影响，使整个空间变成均匀的介电常数为变成均匀的介电常数为 的空间，则空间任一点的空间，则空间任一点 P 的电位由的电位由 q 及及 q 共同产生，即共同产生，即 考虑到无限大导体平面的电位为零考虑到无限大导体平面的电位为零，求得，求得6 电场线与等位面的分布特性与第二章所述的电偶极子的上半电场线与等位面的分布特性与第二章所述的电偶
5、极子的上半部分完全相同。部分完全相同。由此可见，电场线处处垂直于导体平面，而零电位面与导体由此可见，电场线处处垂直于导体平面，而零电位面与导体表面吻合。表面吻合。电场线等位线 z 7fqo（2）点电荷与导体球。）点电荷与导体球。Padrq 若导体球若导体球接地接地，导体球的电位，导体球的电位为零。为了等效导体球边界的影响，为零。为了等效导体球边界的影响，令镜像点电荷令镜像点电荷q 位于球心与点电荷位于球心与点电荷 q 的连线上。那么，球面上任一点的连线上。那么，球面上任一点电位为电位为 可见，为了保证球面上任一点电位为零，必须选择镜像电荷为可见，为了保证球面上任一点电位为零，必须选择镜像电荷为
6、 8 为为了了使使镜镜像像电电荷荷具具有有一一个个确确定定的的值值，必必须须要要求求比比值值 对对于于球球面面上上任任一一点点均均具具有有同同一一数数值值。由由上上图图可可见见，若若要要求求三三角角形形 OPq 与与 OqP 相似，则相似，则 常数。由此获知镜像电荷应为常数。由此获知镜像电荷应为镜像电荷离球心的距离镜像电荷离球心的距离d 应为应为 这样，根据这样，根据 q 及及 q 即可计算球外空间任一点的电场强度。即可计算球外空间任一点的电场强度。fqOPadrq9l（3）线电荷与带电的导体圆柱。）线电荷与带电的导体圆柱。Pafdr-lO 在圆柱轴线与线电荷之间，离轴线的距离在圆柱轴线与线电
7、荷之间，离轴线的距离d 处，平行放置一根处，平行放置一根镜像电荷镜像电荷 。已知无限长线电荷产生的电场强度为。已知无限长线电荷产生的电场强度为 因此，离线电荷因此，离线电荷r 处，以处，以 为参考点的电位为为参考点的电位为 10 若若令令镜镜像像线线电电荷荷 产产生生的的电电位位也也取取相相同同的的 作作为为参参考考点点，则则 及及 在圆柱面上在圆柱面上 P 点共同产生的电位为点共同产生的电位为 已知导体圆柱是一个等位体，因此，为了满足这个边界条件，已知导体圆柱是一个等位体，因此，为了满足这个边界条件，必须要求比值必须要求比值 为常数。与前同理，可令为常数。与前同理，可令 ，由此得，由此得 1
8、1 （4）点电荷与无限大的介质平面。）点电荷与无限大的介质平面。E 1 1qr0EEtEnq 2 2qE 1 2qeten=+为了求解上半空间的场可用镜像电荷为了求解上半空间的场可用镜像电荷 q 等效边界上束缚电等效边界上束缚电荷的作用，将整个空间变为介电常数为荷的作用，将整个空间变为介电常数为1 的均匀空间。对于下的均匀空间。对于下半空间，可用位于原点电荷处的半空间，可用位于原点电荷处的q 等效原来的点电荷等效原来的点电荷q 与边界与边界上束缚电荷的共同作用，将整个空间变为介电常数为上束缚电荷的共同作用，将整个空间变为介电常数为2 的均匀的均匀空间。空间。12 但是，必须迫使所求得的场符合原
9、先的边界条件，即电场切向但是，必须迫使所求得的场符合原先的边界条件，即电场切向分量保持连续，电位移的法向分量应该相等，即分量保持连续，电位移的法向分量应该相等，即 已知各个点电荷产生的电场强度分别为已知各个点电荷产生的电场强度分别为代入上述边界条件，求得镜像电荷如下：代入上述边界条件，求得镜像电荷如下：13不接地：导体球面电位不为不接地：导体球面电位不为0 0，球面上存在正、负感应电荷球面上存在正、负感应电荷（感应电荷总量为（感应电荷总量为0 0）。）。处理方法：电位叠加原理处理方法：电位叠加原理(5)点电荷对不接地球面导体边界的镜像141 1、先假设导体球面接地，则球面上存在电量为、先假设导
10、体球面接地，则球面上存在电量为 的感应电荷，的感应电荷，镜像电荷可采用前面的方法确定。镜像电荷可采用前面的方法确定。2 2、断开接地。将电量为、断开接地。将电量为 的电荷加到导体球面上，这些电荷必的电荷加到导体球面上，这些电荷必然均匀分布在球面上，以使导体球为等势体。然均匀分布在球面上，以使导体球为等势体。3 3、均匀分布在导体球面上的电荷、均匀分布在导体球面上的电荷 可以用位于球心的等量点可以用位于球心的等量点电荷等效。电荷等效。分析过程分析过程结论：点电荷结论：点电荷q q对对非接地导体球面非接地导体球面的镜像电荷有的镜像电荷有两个两个：15镜像电荷镜像电荷1 1：电量：电量：位置：位置：
11、镜像电荷镜像电荷2 2：电量：电量：位置：位置：位于球心。位于球心。球外空间某点电位为：球外空间某点电位为：点电荷位于不接地导体点电荷位于不接地导体球附近的场图球附近的场图162 2 分离变量法分离变量法 分分离离变变量量法法是是把把一一个个多多变变量量的的函函数数表表示示成成几几个个单单变变量量函函数数乘乘积积的的方方法法。它它首首先先要要求求给给定定边边界界与与一一个个适适当当坐坐标标系系的的坐坐标标面面相相合合；其其次次要要求求在在坐坐标标系系中中，待待求求偏偏微微分分方方程程的的解解可可表表示示为为三三个个函函数数的的乘乘积积，且且其其中中的的每每个个函函数数仅仅是一个坐标的函数。在直
12、角、圆柱、球等坐标系中都可以应用分离变量法。是一个坐标的函数。在直角、圆柱、球等坐标系中都可以应用分离变量法。17直角坐标系中的平行平面场问题直角坐标系中的平行平面场问题 平行平面场中位函数U(x，y)在场域内满足拉普拉斯方程 设定分离变量形式的试探解，即令U(x，y)=X(x)Y(y)，代入方程得 18在x和y取任意值时等式恒成立，这要求两边恒为同一常数。现记该常数为(称为分离常数)：取不同值时，两个常微分方程得到不同形式的解：=0 时，时，时，19位函数U的一般解可记作：20 如果问题的边界面与直角坐标系的坐标面吻合，则可采用直角坐标系中的分离如果问题的边界面与直角坐标系的坐标面吻合，则可
13、采用直角坐标系中的分离变量法。变量法。下面通过例子具体说明该方法。下面通过例子具体说明该方法。例例 求如图所示二维长方形内的电位函数。求如图所示二维长方形内的电位函数。解：根据题意，所求区域的电位函数满足解：根据题意，所求区域的电位函数满足的方程及边界条件为的方程及边界条件为x xa ay yb b21只与只与x x有关有关只与只与y y有关有关在直角坐标系中方程在直角坐标系中方程 可写为可写为（二维问题，与（二维问题，与z z无关）无关）分离变量法的前提即假设分离变量法的前提即假设待求函数有分离变量形式的解：待求函数有分离变量形式的解：上式两端同除以上式两端同除以因此该式成立的条件：因此该式
14、成立的条件：且且 为实数为实数为虚数为虚数为虚数为虚数为实数为实数 为零为零为零为零 为实数为实数为虚数为虚数为零为零或或22同样的讨论适用于函数同样的讨论适用于函数 。为满。为满足足x=0 x=0和和x=ax=a的边界条件，应选取的边界条件，应选取则则因为将边界条件将边界条件将边界条件将边界条件于是于是称为边值问题称为边值问题的的本征值。本征值。它的意义是：在上述边界它的意义是：在上述边界条件下，分离常数条件下，分离常数 只有取这些特只有取这些特定值时，方程才有非零解。其解的函定值时，方程才有非零解。其解的函数形式数形式 称为称为本征函数。本征函数。23对于对于因为将边界条件将边界条件于是于
15、是得得由于由于 故故 的一般形式的一般形式 将边界条件将边界条件 这实际上是将一已知函数展为这实际上是将一已知函数展为傅里叶级傅里叶级数。数。利用傅里叶级数的系数公式得利用傅里叶级数的系数公式得原问题的解原问题的解243 有有 限限 差差 分分 法法 图 3.1 差分网格 3.1 差分表示式差分表示式 25二维泊松方程的差分格式 (Difference Form of 2D Poissons Equation)（1）二维静电场边值问题 基本思想：将场域离散为许多网格，应用差分原理，将求解连续函数 的微分方程问题转换为求解网格节点上 的代数方程组的问题。（2）有限差分的网格分割26令 h=x-x
16、0，将 x=x1 和 x3 分别代入式(3)（3）由式（4）+（5）（7）同理,沿 x方向在 x0 处的泰勒公式展开为下 页上 页返 回（4）（5）（6）27将式(6)、式(7)代入式(1)，得到当场域中即即五点差分格式下 页上 页返 回28上式表明，任一点的电位等于它周围四个点电位的平均值。显然，当h越小，计算越精确。如果待求N个点的电位，就需解含有N个方程的线性方程组。若点的数目较多，用迭代法较为方便。29矩形网格剖分若场域离散为矩形网格，差分格式为303.2 边界条件离散化（Discrete Boundary Condition)第二类边界条件 第一类边界条件 分界面衔接条件 对称边界条
17、件 其中介质分界面对称分界313.3差分方程的数值解法差分方程的数值解法 1.简单迭代法简单迭代法 图 3.2 节点序号 32 2.塞德尔塞德尔(Seidel)迭代法迭代法 通常为节约计算时间，对简单迭代法要进行改进，每当算出一个节点的高一次的近似值，就立即用它参与其它节点的差分方程迭代，这种迭代法叫做塞德尔(Seidel)迭代法。塞德尔迭代法的表达式为 此式也称为异步迭代法。由于更新值的提前使用，异步迭代法比简单迭代法收敛速度加快一倍左右，存储量也小。333.超松驰迭代法超松驰迭代法 式中称为松弛因子，其值介于 1 和 2 之间。当其值为1时，超松弛迭代法就蜕变为塞德尔(Seidel)迭代法
18、。因子的选取一般只能依经验进行。但是对矩形区域，当M、N都很大时，可以由如下公式计算最佳收敛因子0:34边界节点赋已知电位值赋节点电位初始值累计迭代次数 N=0N=N+1按超松弛法进行一次迭代，求 打印 NY程序框图35 例例 设如图 所示的矩形截面的长导体槽，宽为 4h，高为3h，顶板与两侧绝缘，顶板的电位为 10V，其余的电位为零，求槽内各点的电位。36 解解：将待求的区域分为 12 个边长为h的正方形网格，含六个内点，得出差分方程组：37解以上方程组，得 38 简简 单单 迭迭 代代 法法 12345600.00.00.00.00.00.012.50.02.50.02.50.023.125 0.6253.750.625 3.125 0.625104.14351.53635.06982.02424.14351.5363114.15151.54195.07780.03564.15151.541939超松弛迭代法超松弛迭代法(=1.2)40松弛因子的影响松弛因子的影响 414243444546本章结束本章结束4748