静态场的解法课件.ppt

《静态场的解法课件.ppt》由会员分享，可在线阅读，更多相关《静态场的解法课件.ppt（60页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、静态场的解法第1页，此课件共60页哦3.1 静态场边值问题及唯一性定理静态场的问题大体上可分为两类：（1）分布型问题 （2）边值型问题常遇到的边值问题有三种：（1）全部边界上的位函数是已知的，称为第一类边值问题，又称为狄利克雷（Dirichlet）问题。（2）全部边界上的法线方向的位函数的导数是已知的，称为第二类边值问题，又称为纽曼（Neumann）问题。第2页，此课件共60页哦（3）混合边值问题，又称为第三类边值问题，它是第一类和第二类边值问题的混合型。静态场的边值问题有多种求解方法，大体可分为以下几种：（1）直接求解法：直接积分法、分离变量法、格林函数法等。（2）间接求解法：复变函数法、镜

2、像法等。（3）数值计算法：有限差分法、有限元法、矩量法等。第3页，此课件共60页哦唯一性定理：不管用什么方法，可以如上任何一种方法，也可以依靠判断猜出解答，只要在给定区域内满足所要求解的微分方程，并满足给定的全部边界条件，那么这个解答就是静态场的唯一解答。第4页，此课件共60页哦3.2直接积分法下面举几个例子来介绍这种方法的应用。例3.2.1 空间电荷区如图3.2.1所示，在-l0区域内为负电荷，在0l区域内为正电荷，且电荷分布函数为=Kqx（x范围为-lxl，K为比例常数）。取x=0为电位参考点，在x=l处电场为零。求在-lxl范围内的电位分布和电场分布。第5页，此课件共60页哦解：在-lx

3、l范围内，电位 满足泊松方程 （为材料的介电常数）从图3.2.1可以看出，电位函数 是x的函数，即 =（x）对上式进行积分 第6页，此课件共60页哦第7页，此课件共60页哦在球外：在球体表面根据边界条件可知在r=R处有 第8页，此课件共60页哦对式（3.11）进行积分得 第9页，此课件共60页哦第10页，此课件共60页哦第11页，此课件共60页哦3.3在直角坐标系中的分离变量法 如果边界面的形状适合用直角坐标系表示，那么可以在直角坐标系中求解。在直角坐标系中，位函数的拉普拉斯方程为 位函数是x、y、z的函数，可以表示成三个单变量未知函数的乘积 （3.49）式中f(x)仅为x的函数；g(y)为y

4、的函数，h(z)为第12页，此课件共60页哦z的函数，得上式除以 得 上式第一项仅是x为变量的函数，与y和z无关；而第二项仅随y而变化，第三项仅随z而变化。所以式（3.50）成立的唯一条件是这三项中每一项都是常数。令第一、二、三项分别为常数 ，即 第13页，此课件共60页哦三个常数满足的关系式为（3.54）这样就把偏微分方程（3.48）变成了三个常微分方程，这种方法就是分离变量法。三个常微分方程（3.51）（3.53）可以改写为第14页，此课件共60页哦 下面讨论拉普拉斯方程对二维位场的求解问题。所谓二维位场即是于是有第15页，此课件共60页哦第16页，此课件共60页哦例3.3.1 两块彼此平

5、行的半无限长接地金属板，板间距离为b，两平行板的一端另一块电位为 的极长的金属条，它们之间缝隙极小，但彼此绝缘，如图3.2所示。求两板间的电位分布。解:给定的边界条件为第17页，此课件共60页哦第18页，此课件共60页哦第19页，此课件共60页哦例3.3.2 两块完全相同的T形导体构成导体槽，两块T形导体间有一狭缝，如图3.3.2(a)所示。上板所加的电压为U0，下板接地。求金属槽内的电位分布。图3.3.2 例3.3.2图第20页，此课件共60页哦解：本题所给的场可以分解为两个场的叠加，分解后的两个场如图3.3.2(b)(c)所示。槽内的电位分别为x、y的函数，是一个二维场问题。分解后第一个场

6、是两个距离为d的无穷大的平行板，上板电压为U0，下板接地，其解为(3.3.23)第二个场电位为2，是两个电位为零的无穷大的平行板，并且在x=0处2满足第21页，此课件共60页哦 那么金属槽内的电位分布的解为=1+2，分别求出1和2，也就得出来了，根据唯一性定理，即是要求的唯一解答。1已知，见（3.3.23）式，下面求出2。由例3.3.1的讨论，2可表示为 (3.3.24)式中Cn由x=0处的边界条件求出，即第22页，此课件共60页哦可以确定傅里叶系数Cn为则(3.3.25)金属槽内的电位分布为 (3.3.26)第23页，此课件共60页哦3.4 在圆柱坐标系和球坐标系的分离变量法1.在圆柱坐标系

7、的分离变量法在圆柱坐标系中，拉普拉斯方程为(3.4.1)下面讨论电位不随纵向（z方向）变化的二维场问题，即仅为r和变化的情况。此时拉普拉斯方程变为 (3.4.2)第24页，此课件共60页哦设 （3.4.3）式中R(r)仅为r的函数，把式（3.4.3）代入到式（3.4.2）中得(3.4.4)式（3.4.4）第一项是关于r的函数，第二项是关于的函数，要使上式对于所有的r和都成立，必须每项都等于一个常数，于是有（3.4.5）(3.4.6)第25页，此课件共60页哦对式（3.4.6）求解，其解为(3.4.7)对于所研究的实际问题，位场是单值的，则有所以k必须为整数，令k=n，于是式（3.4.7）变为(

8、3.4.8)用n代替k，并把式（3.4.5）改写为如下形式(3.4.9)它是一个欧拉方程，其解为(3.4.10)第26页，此课件共60页哦（3.4.11）式中的系数由边界条件确定第27页，此课件共60页哦 2.在球坐标系的分离在球坐标系的分离变量法变量法球坐标系中的拉普拉斯方程为在球坐标系中具有轴对称的二维场的解式中常数Am、Bm由边界条件决定。例3.8 无限大介质外加均匀电场，在介质内有一个半径为a的球形空腔，介质的介电常数为，求空腔内、外的电位分布及电场强度。第28页，此课件共60页哦解 本题为球坐标系中具有轴对称性的二维场问题 在空腔内的通解为在介质中的通解为第29页，此课件共60页哦下

9、面利用边界条件确定各个系数。所以B1=0 系数A1、C1、D1可以由r=a时的边界条件求出，当r=a时由1=2 所以可以得出 A1acos=（C1a+D1a-2）cos第30页，此课件共60页哦第31页，此课件共60页哦第32页，此课件共60页哦3.5镜像法镜像法最简单的情况是点电荷对无限大平面的镜像问题。在无限大接地导体平面上方的空气中放置一个电荷q，距平面距离为h，如图3.5.1所示，要求出这个无限大接地导体平面上方任意一点的电位,那么可以设想在平面下方有一个镜像电荷-q，在所研究的区域即平面上方的电位为点电荷q产生的电位和镜像电荷-q产生电位的叠加。第33页，此课件共60页哦图3.5.1

10、 镜像法求解点电荷与无限大接地导体平面形成的位场第34页，此课件共60页哦 建立一个坐标系，如图3.5.2所示，可以很方便地求出平面上方任意一点的电位分布及无限大接地平面上的电荷分布。图3.5.2 在直角坐标系中求解点电荷与无限大接地导体平面问题第35页，此课件共60页哦从图3.5.2可以看出 R=x2+y2+（z-h）21/2 R=x2+y2+（z+h）21/2所以在平面上方空间任意一点P的电位为平面上方空间任意一点P的电场强度为第36页，此课件共60页哦无限大接地导体平面的面电荷分布为例3.9 相交成直角的接地导体平面AOB附近有一个点电荷q，如图3.10所示，中间介质为空气，求空间任意一

11、点P的电位。图3.10直角形导体平面的镜像第37页，此课件共60页哦解:点电荷q位置为（h1,h2,0），在OA面的镜像为-q，位置在（-h1,h2,0）。在OB面的镜像为-q，位置在（-h1,h2,0）。按照平面镜像法则还应该成第三个镜像位置在（-h,-h2，0），第三个镜像为q。这样点电荷q与三个镜像电荷共同作用才能满足原来的边界条件在导体平面AOB上的电位为零。所以本问题可以用三个镜像电荷代替相交成直角的接地导体平面AOB。得到P点的电位为第38页，此课件共60页哦例3.10 一根无限长直导线与地面平行，设导线半径为a，高出地面的高度为h(ha)，求单位长度导线对地的电容。图3.11无限

12、长直导线的镜像法第39页，此课件共60页哦解:本题可以采用镜像法求解，无限长的直导线的线电荷密度为l，则地面这个边界可以用镜像的线电荷密度为-l的直导线代替，如图3.5.4所示。原导线在P点的电场强度为第40页，此课件共60页哦同理镜像电荷的电位为 所以由镜像法可知，地面上方任意一点P的电位为则直导线的电位为所以单位长度导线对地的电容为第41页，此课件共60页哦例3.11 有一个接地导体球半径为a，与球心O相距d1的位置P1点有一个点电荷q1如图3.5.5所示，试求：(1)导体球外的电位函数；(2)球面上感应的面电荷密度；(3)球面上总感应电量。图3.12 点电荷对导体球的镜像第42页，此课件

13、共60页哦解（1）导体球接地时，导体表面电位处处为零。下面在球内区域找到一个镜像电荷q2，用q2来代替球形边界。在球面上任取一点P，则P点的电位为零，即上式对球面上任意一点都是成立的，那么可以在球面上取两个特殊点：一个是OP1与球面的交点，另一个是OP1的延长线与球面的交点，于是有第43页，此课件共60页哦由上面两个方程联立解得这样球外任意一点P的电位为以O为中心建立球坐标系，设OP与OP1夹角为第44页，此课件共60页哦所以由在r=a时的边界条件可知第45页，此课件共60页哦（3）球面上总感应电量为计算结果表明，球面上总感应电量等于镜像电荷的电量。第46页，此课件共60页哦3.6静态场的数值

14、解法1.数值积分法 电磁场分布型问题可以利用数值积分法计算。例如，可以在各向同性电介质中的电场强度E和电位 式中为源区。如为体分布，则需对V体积分；如为面分布，则需对S面积分；如为线分布，则需对l线积分。第47页，此课件共60页哦 再如对于一闭合回路l，电流为I，在场点的磁感应强度B和矢量磁位A也可以通过数值积分法计算 诸如此类的电磁场问题很多，通过矢量计算都可以得到数值积分问题。数值积分的算法很多，如梯形求积算法、辛普生求积算法等。下面通过一个例子介绍梯形求积算法。例3.6.1 求均匀带电直导线的中垂线上一点P的电场强度，设带电直导线的长度为l，带电量为q。第48页，此课件共60页哦解:取棒

15、的中点为坐标系原点，在带电体上取电荷元 图3.6.1 例3.6.1用图（这里加“”表示源点），则dq和P点的电场强度的大小为根据对称性可知，P点的电场强度只有x分量，而y分量Ey=0，即第49页，此课件共60页哦 现在对上式进行数值积分，把被积函数f(x)在积分区间(a,b)取N个小梯形，则第n个小梯形的面积为f(xn)xn，则f(x)在(a,b)上的积分为按照梯形算法，每一个小梯形区间宽度为 ，第n个梯形采样点为第50页，此课件共60页哦则然后编写程序计算数值解。第51页，此课件共60页哦2.有限差分法 在一个闭合边界L所界定的平面域，其定解问题可表述为首先是把求解的场域离散化，即在求解的场

16、域划分成网格，网格的划分有许多种方法，最简单的是正方形网格划分，如图3.6.2所示。然后对偏微分方程进行离散化，对正方形网格可采用五点差分格式。在二维场域中取一点P，则沿x轴方向并通过P点的直线上任意一点的数值ux用泰勒公式展开为第52页，此课件共60页哦 式中h=（x-xP）很小时，4阶以上的高次项可以忽略。P点与其相邻的四个网格点构成了五个网络点，如图3.6.3所示。图3.6.2二维场域的正方形网格剖分 图3.15正方形网格的五点差分格式第53页，此课件共60页哦则1点的值为 3点的值为上两式相加得同理取P点为（xi,yi），则uP=ui,j；u1=ui+1,j；u2=ui,j+1；u3=

17、ui-1,j；u4=ui,j-1。则有ui,j=1/4（ui+1,j+ui,j+1+ui-1,j+ui,j-1-h2F）第54页，此课件共60页哦例3.6.2 如图3.6.4所示为一长直接地金属槽，其侧壁与底面电位均为零，顶盖电位的相对值为10，试求槽中的电位分布。解:在金属槽中题中的边界构成第一类边值问题按照有限差分法，本题的解题过程如下：1）离散化场域。用正方形网格对场域D进行剖分，剖分越细计算精度越高。根据工程需要，沿x、y方向的等分数应取p=30以上，本题中为了形象直观，沿x、y方向的等分数均为p=4，步距h=a/4，如图3.6.5所示。第55页，此课件共60页哦图3.6.4 例3.6

18、.2图图3.6.5场域离散化（2）采用超松弛迭代法，本题中F=0，得到差分方程的迭代运算形式为式中加速收敛因子可以通过一定的方法估算给出，本题中加速收敛因子的最佳取值通过估算取为opt=1.17第56页，此课件共60页哦（3）给出边界条件和初始值。因为本题中给定的是第一类边界条件，则边界条件的差分离散化应采用直接赋值方式，即 给定其他网络点的初始值。可取初始值都为零，但需迭代次数较多，一般可按照其数值规律设定。如本例中可设网格内点的初始值为第57页，此课件共60页哦（4）给出计算流程框图。流程框图中要给定迭代解收敛指标，如本例中可以设DF=10-5，即当各网格内点相邻二次迭代近似值的绝对误差的

19、绝对值小于DF=10-5时，终止迭代循环。（5）编写计算程序。可以采用许多编程语言来编写，如Fortran77、C语言、MATLAB语言等，这里不再给出4.有限元法 以下面的二维拉普拉斯方程的第一类边值问题为例来分析有限元法的原理。设在曲线C为边界的场域S内，电位函数的边值问题为第58页，此课件共60页哦在场域S内的电磁能量为可以看出，电场能量W是电位的函数，称为能量泛函。（3.6.2）可以证明，在边界条件（3.6.2）下，使电场能量W取得极小值的电位函数，必定满足式（3.6.1）。这样，原来的电位函数关于边界条件求解的问题等价为求泛函极限问题（3.6.1）第59页，此课件共60页哦 然后对场域进行剖分与插值，把变分问题离散化。场域的剖分方法很多，二维场域最常用的是三角形单元剖分，如图3.6.6所示，各三角形单元的形状和大小可以是任意的，可以很好地适应边界形状。场域三角形剖分的原则是：任何一个三角形单元的顶点必须同时为其相邻三角形单元的顶点，相邻的单元既不能相互重叠也不能相互分离，各单元节点编号顺序一律按逆时针方向进行。线性插值是最常用的方式。场域中任取一个单元e，如图3.6.7所示，其三个顶点编号顺序标记i、j和m。第60页，此课件共60页哦