领导与竞争优势专题帝国传奇郭台铭的领导特质组别：第五组课件.ppt

1、1第5章 静电场边值问题的解法 5.1 边值问题分类边值问题分类5.2 积分法积分法5.3 唯一性定理唯一性定理5.4 分离变量法分离变量法5.5 镜像法镜像法5.6 有限差分法有限差分法习题习题第5章 静电场边值问题的解法2第5章 静电场边值问题的解法本章提要l边值问题分类l唯一性定理l分离变量法l镜像法l有限差分法3第5章 静电场边值问题的解法前面我们讨论了基于库仑定律与叠加原理或高斯定理计算电场的方法，这些方法只能适用于已知电荷分布十分简单的问题。实际上在电工中经常遇到的是这样一类问题：给定空间某一区域内的电荷分布(可以是零)，同时给定该区域边界上的电位或电场(即边值，或称边界条件)，在

2、这种条件下求解该区域内的电位函数或电场强度分布。从数学上来讲，这些问题都是在给定的定解条件边界条件下，求解泊松方程或拉普拉斯方程的定解问题，称为边值问题。4第5章 静电场边值问题的解法静电场边值问题的解法分为解析法和数值法两大类。用解析法得到的场量表达式是准确值，但是它只能解决规则边界的边值问题。本章主要介绍解析法，包括积分法、分离变量法、镜像法。数值解属于近似计算，但对于不规则边界等复杂的静电场问题是非常有用的方法。随着计算机的广泛使用，数值法已成为边值问题求解的主要方法。由于篇幅所限，因而本章只简单介绍数值法中的一种有限差分法。5第5章 静电场边值问题的解法根据给定求解区域边界条件的不同，

3、边值问题分为以下三类。第一类边界条件又称狄利赫莱(Dirichlet)条件：场域边界S上的电位分布已知，即(51)式中rb为相应边界点的位置矢量。5.1 边值问题分类边值问题分类 1bSfrr6第5章 静电场边值问题的解法第二类边界条件又称纽曼条件：场域边界S上电位的法向导数分布已知，即(52)当f2(rb)取零时，称为第二类齐次边界条件。2bSfnrr7第5章 静电场边值问题的解法第三类边界条件又称混合条件：场域边界S上电位及其法向导数的线性组合已知，即(53)在实际问题中，除了给定边界条件外，有时还需要引入某些补充的物理约束条件，称为自然边界条件。在求解边值问题中，自然边界条件非常重要，但

4、它又不是事先给定的，必须根据问题自行确定，举例如下。34bSffnrrrr8第5章 静电场边值问题的解法无限远边界条件：对于电荷分布在有限域的无边界电场问题，在无限远处有(r)|r0，即电位在无限远处趋于零。介质分界面条件：当场域中存在多种媒质时，还必须引入不同介质分界面上的边界条件，即用电位表示的边界条件和 1 2，此边界条件又称为辅助的边界条件。2122snn9第5章 静电场边值问题的解法对于一些具有对称结构的静电场问题，电位函数仅是一个坐标变量的函数。静电场边值问题可归结为常微分方程的定解问题，这时可以直接积分求解电位函数。【例5-1】求真空中球状分布电荷所产生的空间电场强度和电位分布，

5、设电荷体密度为5.2 积积 分分 法法ararr 00 110第5章 静电场边值问题的解法解解 设球状电荷分布内、外的电位分别为1和 2。显然，1满足泊松方程，2满足拉普拉斯方程。由于电荷分布的球对称性，选取球坐标系，因而有221120022222d1 d1ddd1 d0ddrrrrrrrrr (0ra)(ra)11第5章 静电场边值问题的解法边界条件为可解得1和 2的通解为112012002d0ddd0ddr ar arrr ar arrr，11203242CrCrCCr 12第5章 静电场边值问题的解法代入边界条件，得最终得到电位函数的解为21423000,0,2aaCCCC 100()2

6、rara 13第5章 静电场边值问题的解法利用球坐标系中的梯度表达式，求得可以证明，以上结果与应用高斯定理求得的结果完全一致。111022222012()2rrrrdradrdaradrr EeeEee()14第5章 静电场边值问题的解法直接积分法只能求解比较简单的问题，在比较复杂的情况下，必须用其他方法求解。在讨论这些方法之前，了解静态场解的唯一性定理非常重要。静态场解的唯一性定理：满足给定边界条件的泊松方程或拉普拉斯方程的解是唯一的。5.3 唯一性定理唯一性定理15第5章 静电场边值问题的解法图51 包围含有导体的场域16第5章 静电场边值问题的解法【例5-2】证明：唯一性定理。证明 图5

7、1所示为充满均匀介质和置有n个导体的场域。场域空间V的边界为S1，S2，Sn及外边界面S0。设V中存在两个电位函数1和 2，对于给定第一类或第二类边界条件，均满足泊松方程，即(54)2212,17第5章 静电场边值问题的解法令d 1 2，因此有 2 d0(55)利用格林公式可得(56)令 d，代入上式得(57)2dVSVdSn 2ddddVSVdSn 18第5章 静电场边值问题的解法式(57)中场域V的边界面SS0S1S2Sn。如果所设的这两个不同的电位函数的解1和 2，在全部边界面上都应有相同的第一类边界条件或第二类边界条件，则它们在相应边界面Si上的差值 将其代入式(57)，有(58)dd

8、00iiSSn或，2dd0VV19第5章 静电场边值问题的解法这说明，场域V内d的梯度处处为零，即V内所有场点上的 d值与其在各导体表面S1，S2，Sn上的值是相同的。对于第一类边值问题，由于在导体表面上已知 d0，因此整个场域内必有d0，由此得证 1 2，即解是唯一的。对于第二类边值问题而言，即已知各导体表面上的面电荷分布，此时 dC，即电位 1和 2之间可能相差一个常数，但采用相同的电位参考点将导致C0，所以解仍是唯一的。20第5章 静电场边值问题的解法静电场唯一性定理的重要意义在于，求解静电场问题时，不论采用哪一种解法，只要在场域内满足相同的偏微分方程，在边界上满足相同的给定边界条件，就

9、可确信其解答是正确的。21第5章 静电场边值问题的解法当待求电位函数是二个或三个坐标变量的函数时，分离变量法是直接求解偏微分方程定解问题的一种经典方法。对于拉普拉斯方程对应的边值问题，其求解步骤是：根据问题中场域边界的几何形状，选用适当的坐标系；设待求电位函数由两个或三个各自仅含一个坐标变量的函数乘积组成，并代入拉普拉斯方程，借助于“分离”常数，将拉普拉斯方程转换为两个或三个常微分方程；解这些常微分方程并以给定的定解条件决定其中的待定常数和函数后，即可解得待求的电位函数。5.4 分离变量法分离变量法22第5章 静电场边值问题的解法一般而言，当场域边界和某一正交曲线坐标系的坐标面相吻合时，分离变

10、量法往往是一种简便而有效的方法。下面以直角坐标系和圆柱坐标系中的二维场为例说明分离变量法的应用。23第5章 静电场边值问题的解法5.4.1 直角坐标系中的二维场问题直角坐标系中的二维场问题设电位函数为(x，y)，其满足拉普拉斯方程：(59)设电位函数可写成分离变量形式，即(x，y)X(x)Y(y)(510)代入拉普拉斯方程，整理得(511)22222,0 x yxy22221 d1 d0ddXYXxYy24第5章 静电场边值问题的解法令或写成(512)2222221 d1 dddxyXYkkXxYy ，2222221 d0d1 d0dxyXkXxYkYy25第5章 静电场边值问题的解法kx、k

11、y称为分离常数，且满足(513)由式(513)可知，二维拉普拉斯方程在分离变量时只有一个独立的分离常数，即kx或ky，而且它们不能同为实数或虚数，但二者可以同为零。现以方程为例进行说明。220 xykk2221 ddxXkXx 26第5章 静电场边值问题的解法(1)若，解为X(x)a0 xb0(514a)(2)若，解为(514b)20 xk20 xk 1122chshxxxxk xk xX xak xbk xX xa eb e或27第5章 静电场边值问题的解法(3)若，解为X(x)a3coskxxb3sinkxx(514c)以上各式a0、b0、a1、b1、a2、b2、a3、b3均为待定常数。Y

12、(y)的解形式与X(x)完全相同，只是分离常数必须满足式(513)要求。20 xk28第5章 静电场边值问题的解法当kx、ky取不同值时，上述解的线性组合便构成了拉普拉斯方程的通解，即(515)最后，可根据给定的边界条件，通过傅里叶级数展开方法，确定各个待定常数。00001111222213333,()()chshchshcossincossinynynxnxnnynnynnxnnxnkxkxkxkxnnnnnnxnnxnnynnynx ya xbc ydck xdk xak xbk xa eb ec ed eak xbk xck xdk x29第5章 静电场边值问题的解法【例5-3】长直接地

13、金属槽的横截面如图52所示，其侧壁与底面电位均为零，顶盖电位为U0。求槽内电位分布。解解 依题意，本问题为第一类边值问题，即在区域0 xa，0y0)的区域中每一点的电位。通常的求解方法是在直角坐标中解拉普拉斯方程：(536)此方程适用于除了点电荷以外的整个z0区域。解(x，y，z)应满足下列条件：22222220 xyz48第5章 静电场边值问题的解法(1)在接地平面的所有点上，电位为零，即(x，y，z)0(2)在很靠近q的点上，其电位趋近于单个点电荷的电位，即其中r是该点与q之间的距离。(3)在离q非常远的点(x，y或z)上，电位趋于零。(4)电位函数是x和y坐标的偶函数，即 由此看来，构造

14、一个满足所有这些条件的解，确实是很困难的。04qr,x y zx y zx y zx yz49第5章 静电场边值问题的解法从另一个观点来看，在zd处存在正电荷q，它将在导体平面的表面感应出负电荷，产生面电荷密度S。因此，导体平面上方各点的电位将是其中r为从dS到所考察点的距离，而S为整个导体平面的表面。这里的问题在于，必须首先根据边界条件(x，y，0)0求出S。而且，即使求出了导体平面每一点的S，计算上式的面积分也是困难的。下面我们通过例题说明使用镜像法如何大大地简化这些问题。222001,d44Ssqx y zSrxyzd50第5章 静电场边值问题的解法5.5.1 对无限大接地导电平面的镜像

15、对无限大接地导电平面的镜像1点电荷和导体平面点电荷和导体平面【例5-5】设有一点电荷q位于距无限大接地导电平面上方h处，其周围介质的介电常数为，如图54所示。求导体平面上方区域空间点的电位。解解 根据题意可知，电位函数在场域内满足如下边值问题：2 0(除去点电荷所在点)边界条件为|z0051第5章 静电场边值问题的解法图54 点电荷对无限大接地导电平面的镜像(a)无限大接地导电平面上的点电荷；(b)点电荷的镜像52第5章 静电场边值问题的解法可以设想，在场域边界外引入一个与点电荷q呈镜像对称的点电荷qq，并将原来的导体场域用介电常数为的介质所替换。这样，原场域边界面(z0)上的边界条件0保持不

16、变，而对应的边值问题被简化为同一均匀介质空间内两个点电荷的电场计算问题。根据唯一性定理可知，其解的有效区域仅限于图54(b)所示上半部分介质场域。应用镜像法，得待求电位为53第5章 静电场边值问题的解法(537)用直接代入法，很容易证明式(537)的(x，y，z)满足式(536)的拉普拉斯方程。显然，它也满足列于式(536)后面的所有四个条件，因此式(537)是原问题的解，而且鉴于唯一性定理，它是唯一的解。222222121111,44qqzrrxyzhxyzh 54第5章 静电场边值问题的解法无限大接地导电平面上感应电荷的面密度分布为式中负号表示感应电荷与点电荷q的极性相反。对感应电荷作面积

17、分，得上式表明镜像电荷q确实等效于无限大接地导电平面上的全部感应电荷。22222022snzzqhqhDEzxyhh 322220 0d2sSqhddSqqh 55第5章 静电场边值问题的解法【例5-6】如图55(a)所示，正点电荷Q与两个接地且互相垂直的导体半平面的距离分别为d1和d2，求平面上的感应电荷作用在Q上的力。解 采用镜像法求解，镜像电荷有三个，如图55(b)所示。点电荷Q所受力为F1、F2、F3矢量之和，见图55(c)，即22212222220201012()4(2)4(2)4(44)yxxyQQQddddddFeeee56第5章 静电场边值问题的解法图55 点电荷和互相垂直的导

18、体平面(a)实际布置图；(b)实际镜像电荷布置；(c)作用在电荷Q上的力57第5章 静电场边值问题的解法2线电荷和导体平面线电荷和导体平面不难看到，位于无穷大导体平面上方的线电荷的电场，可以根据和它的镜像求出(导体平面已不存在)。【例5-7】求图56所示线电荷及其镜像电荷，求空间电位函数。解解 由高斯定理得点P的电场强度为21PPPeeEEE201022 58第5章 静电场边值问题的解法图56 线电荷的镜像(a)线电荷对无限大接地平面；(b)线电荷的镜像59第5章 静电场边值问题的解法现任取点Q为电位参考点，则点P电位为根据边界条件，在无限大接地导电平面上，即12时，0，代入上式，得C0。场中

19、任意点电位为21QQ111Q12Q20000201ddlnlnlnln2222ln2PC 1222222010lnln22Pxbyxby60第5章 静电场边值问题的解法5.5.2 对无限大介质平面的镜像对无限大介质平面的镜像对于图57所示无限大介质平面上的点电荷边值问题也可采用镜像法。上、下半无限空间中的电场是由点电荷q及其分界面上的束缚电荷共同产生的。对于介质为1的上半空间的电场计算，其分界面上的束缚电荷可归结为在均匀介质1中的镜像点电荷q；对于介质为2的下半空间的电场计算，其分界面上的束缚电荷可归结为在均匀介质2中的点电荷q。镜像电荷q和q的量值，可以通过分界面上的边界条件确定如下。61第

20、5章 静电场边值问题的解法对于分界面上任意点P，由其上的边界条件E1tE2t和D1nD2n，得 对于线电荷与无限大介质平面系统的电场，可类比推得。cos4cos4cos4222121rqrqrq sin4sin4sin4222rqrqrq q2qqq2122121 ,62第5章 静电场边值问题的解法图57 无限大介质平面镜像(a)无限大介质平面上的点电荷；(b)上半空间电场计算的镜像；(c)下半空间电场计算的镜像63第5章 静电场边值问题的解法前几节讨论了求解拉普拉斯方程的解析法。但是对大多数实际问题，住往边界形状复杂，很难用解析法求解，为此需使用数值方法。目前已发展了许多有效的求解边值问题的

21、数值方法，有限差分法是一种较易使用的数值方法。5.6 有限差分法有限差分法64第5章 静电场边值问题的解法有限差分法的思想是将场域离散为许多网格，应用差分原理差分方程代替各个点的偏微分方程，这样得到的任意一个点的差分方程是将该点的电位与其周围几个点相联系的代数方程，对于全部的待求点，就得到一个线性方程组，从而将求解连续函数的微分方程问题转换为求解网格节点上的代数方程组问题。本节简要说明有限差分法的基本原理，并以二维拉普拉斯方程的第一类边值问题为例，说明有限差分法的求解过程。65第5章 静电场边值问题的解法5.6.1 差分表示式差分表示式在xOy平面把所求解区域划分为若干相同的小正方形格子，每个

22、格子的边长都为h，如图58所示。假设某顶点0上的电位是0，周围四个顶点的电位分别为1、2、3和4，将这几个点的电位用泰勒级数展开，就有(538)(539)23231023000112!3!hhhKxxx23233023000112!3!hhhKxxx66第5章 静电场边值问题的解法图58 划分求解区域67第5章 静电场边值问题的解法当h很小时，忽略四阶以上的高次项，得(540)同理，也可得到(541)将式(540)与式(541)相加，并考虑22130202hx22240202hx68第5章 静电场边值问题的解法(542)式(542)表明，任一点的电位等于它周围四个点电位的平均值。显然，h越小，

23、计算越精确。如果待求N个点的电位，就需解含有N个方程的线性方程组。若点的数目较多，则用迭代法较为方便。22220 xy012341469第5章 静电场边值问题的解法5.6.2 差分方程的数值解法差分方程的数值解法如前所述，平面区域内有多少个节点，就能得到多少个差分方程。当这些节点数目较大时，使用迭代法求解差分方程组比较方便。70第5章 静电场边值问题的解法1简单迭代法简单迭代法用迭代法解二维电位分布时，将包含边界在内的节点均以双下标(i，j)表示，i、j分别表示沿x、y方向的标号。次序是x方向表示从左到右，y方向表示从下到上，如图59所示。我们用上标n表示某点电位的第n次迭代值，由式(542)

24、得出点(i，j)第n1次电位的计算公式为(543)1,1,11,114nnnnni jiji jiji j71第5章 静电场边值问题的解法式(543)也叫简单迭代法，它的收敛速度较慢。计算时，先任意指定各个节点的电位值作为零级近似(注意电位在某无源区域的极大、极小值总是出现在边界上，理由请读者自行思考)，将零级近似值及其边界上的电位值代入式(543)求出一级近似值，再由一级近似值求出二级近似值。依此类推，直到连续两次迭代所得电位的差值在允许范围内时，结束迭代。对于相邻两次迭代解之间的误差，通常有两种取法：一种是取最大绝对误差；另一种是取算术平均误差，其中N是节点总数。1,maxkki ji j

25、i j1,1kki ji ji jN72第5章 静电场边值问题的解法图59 节点序号73第5章 静电场边值问题的解法2塞德尔塞德尔(Seidel)迭代法迭代法通常为节约计算时间，需对简单迭代法进行改进，每当算出一个节点高一次的近似值，就立即用它参与其他节点的差分方程迭代，这种迭代法叫做塞德尔(Seidel)迭代法。塞德尔迭代法的表达式为(544)式(544)也称为异步迭代法。由于更新值的提前使用，异步迭代法比简单迭代法收敛速度加快一倍左右，存储量也小。1,1,11,114nnnnni jiji jiji j74第5章 静电场边值问题的解法3超松驰迭代法超松驰迭代法为了加快收敛速度，常采用超松驰

26、迭代法。计算时，将某点的新老电位值之差乘以一个因子以后，再加到该点的老电位值上，作为这一点的新电位值。超松弛迭代法的表达式为(545)式中称为松弛因子，其值介于1和2之间。当其值为1时，超松弛迭代法就蜕变为塞德尔(Seidel)迭代法。1,1,11,1,144nnnnnnni ji jiji jiji ji j1,ni j75第5章 静电场边值问题的解法因子的选取一般只能依经验进行，但是对矩形区域，当沿x、y两个方向的内节点数M、N都很大时，可以由如下公式计算最佳收敛因子0：(546)022222MN76第5章 静电场边值问题的解法 对于其他形状的实际区域，最佳收敛因子的表达式很复杂。实际计算

27、中，往往应用其近似值，通常采用以下几种方法处理：一是将区域等效为近似的矩形区域，再依照式(546)计算0；二是编制可以自动选择收敛因子的计算程序，在起始迭代时取收敛因子为1.5，然后依迭代过程收敛速度的快慢使计算机按程序自动修正收敛因子；三是起始迭代取收敛因子为1，以后逐渐增大，并注意观察迭代过程的收敛速度，当速度减小时，停止增加收敛因子的值，而在以后的迭代中，用最后一个收敛因子的值作为最佳值。77第5章 静电场边值问题的解法【例5-8】设图510所示矩形截面的长导体槽，宽为4h，高为3h，顶板与两侧绝缘，顶板的电位为10 V，其余的电位为零，求槽内各点的电位。解 将待求的区域分为12个边长为

28、h的正方形网格，含六个内点，得出差分方程组：78第5章 静电场边值问题的解法12311 0421414314511 0442361453611 046451479第5章 静电场边值问题的解法图510 例57用图80第5章 静电场边值问题的解法解以上方程组，得【注】以上结果是差分方程组的精确解，但并不是待求各点电位的精确值，这是因为差分方程组本身是原微分组的近似。以下用迭代法求解，简单迭代法和超松弛迭代法的结果分别列于表51和表52，表53给出了收敛因子的影响。16704.1615161V22501.5528161V38205.0932161V43302.0497161V56704.161516

29、1V62501.5528161V81第5章 静电场边值问题的解法表表51 简单迭代法简单迭代法82第5章 静电场边值问题的解法表表52 超松弛迭代法超松弛迭代法(1.2)83第5章 静电场边值问题的解法表表53 松弛因子的影响松弛因子的影响84第5章 静电场边值问题的解法1静态场的许多问题可归结为给定边界条件下求解位函数的泊松方程或拉普拉斯方程的问题，也称为边值型问题。满足给定边界条件的泊松方程或拉普拉斯方程的解是唯一的。2镜像法是一种等效方法，它的基础是静电场解的唯一性定理。利用镜像法求解诸如无限大导体(或介质)平面附近的点电荷、线电荷产生的场，位于无限长圆柱导体附近的平行线电荷产生的场，以

30、及位于导体球附近的点电荷产生的场等特定问题非常有效。镜像法的主要步骤是确定镜像电荷的位置和大小。本章小结本章小结85第5章 静电场边值问题的解法3分离变量法是将一个多元函数表示成几个单变量函数的乘积，从而将偏微分方程分离为几个带分离常数的常微分方程的方法。用分离变量法求解边值型问题，首先要根据边界形状，选择适当的坐标系；然后将偏微分方程在特定的坐标系下分离为几个常微分方程，并得出位函数的通解；最后由边界条件确定通解中的待定常数。86第5章 静电场边值问题的解法4有限差分法应用差分原理将待求场域的空间离散化，把拉普拉斯方程化为各节点上的有限差分方程，并使用迭代法求解差分方程，从而可以求出节点上的

31、位函数值。5由于静电场的解法具有普遍性，因此本章讨论的求解方法不仅适用于静电场，也适用于恒定电场和恒定磁场，在某些情况下，也适用于时变电磁场。87第5章 静电场边值问题的解法知识结构图知识结构图88第5章 静电场边值问题的解法5.1 求截面为矩形的无限长区域(0 xa，0y0半空间充满介电常数为1的介质，z0半空间充满介电常数为2的介质。(1)电量为q的点电荷放在介质分界面上时，求电场强度。(2)电荷线密度为1的均匀线电荷放在介质分界面上时，求电场强度。5.9 在内、外半径分别为a和b的圆柱形区域内无电荷，在半径分别为a和b的圆柱面上电位分别为U和0，求该圆柱形区域内的电位和电场。97第5章

32、静电场边值问题的解法5.10 在半径分别为a和b的两同轴导电圆筒围成的区域内，电荷分布为A/r，A为常数。若介质介电常数为，内导体电位为U，外导体电位为0，求两导体间的电位分布。5.11 电位分别为0和U的半无限大导电平板构成夹角为的角形区域(见题5.11图)，求该角形区域中的电位分布。98第5章 静电场边值问题的解法题5.11图99第5章 静电场边值问题的解法5.12 无限大的导电平板上方距导电平板h处平行放置无限长的线电荷，电荷线密度为1，求导电平板上方的电场。5.13 由无限大的导电平板折成45的角形区，在该角形区中某一点(x0，y0，z0)有一点电荷q，用镜像法求电位分布。5.14 半

33、径为a、带电量为Q的导体附近距球心f处有一点电荷q，求点电荷q所受的力。100第5章 静电场边值问题的解法5.15 内、外半径分别为a、b的导体球壳内距球心为d(da)处有一点电荷q，当(1)导电球壳电位为0；(2)导电球壳电位为；(3)导电球壳上的总电量为Q；时，分别求导电球壳内外的电位分布。5.16 无限大导电平面上有一导电半球，半径为a，在半球体正上方距球心及导电平面h处有一点电荷q(见题5.16图)，求该点电荷所受的力。101第5章 静电场边值问题的解法题5.16图102第5章 静电场边值问题的解法5.17 如题5.17图所示，无限大导电平面上方平行放置一根半径为a的无限长导电圆柱，该导电圆柱轴线与导电平面的距离为h，求导体圆柱与导电平面之间单位长度的电容。103第5章 静电场边值问题的解法题5.17图104第5章 静电场边值问题的解法5.18 两同心导体球壳半径分别为a、b，两导体之间有两层介质，介电常数为1、2，介质界面半径为c，求两导体球壳之间的电容。5.19 利用有限差分法求题5.19图所示区域中各个节点的电位。105第5章 静电场边值问题的解法题5.19图