1、第二章第二章 静电场静电场Electrostatic field本章研究的主要问题是:在给定的自由电荷分布以及周围空间介质和导体分布的情况下,如何求解电场。注意两点:电荷静止,即:电场不随时间变化,即:本章求解静电场的方法有:分离变量法;镜像法;格林函数法。求解的依据是:唯一性定理。00tE2.1 静电场的标势及其微分方程静电场的标势及其微分方程Scalar potential and differential equation for electrostatic field1.静电场的标势和微分方程静电场的标势和微分方程 静电现象满足以下两个条件:即 电荷静止不动;场量不随时间变化。故把静电
2、条件代入Maxwells equations中去,即得电场满足的方程0;()0jt物理量DE0这两方程连同介质的电磁性质方程 是解决静电问题的基础。根据电场方程 (即 的无旋性),可引入一个标势 。在电磁学中,已知 因为相距为 两点的电势差为由于所以ED0EABldEBA)()(l dldEdl ddzzdyydxxdEE又因为在均匀各向同性的介质中,则有这里 ,故有即此方程称为泊松方程(Poisson equation).若在无源区域内(),上式化为EDEEED)(0)(ED2020此方程称为拉普拉斯方程(Laplace equation)在各种不同条件下求解Poisson equation
3、或Laplace equation是处理静电问题的基本途径。2、静电场的基本问题静电场的基本问题 如果电荷是连续分布的,则观察点 处的标势为这个式子只反映了电荷激发电场这一面,而没有反映电场对电荷的作用另一面。如果空间还有导体存在的活,那么物理机制为x01()()4VxxdVr 考虑到感应情况,诸问题的模拟是:现在,要找出一个电荷对它邻近的电场是怎样作用的,一点上的电场和它邻近的电场又是怎样联)(x导体+-给定电荷分布求空间一点电场分布而场引起导体上感 应电荷分布而感应电荷分布反过来引起系的,即要找出电荷和电场相互作用规律的微分形式,而在导体表面或其他边界上场和电荷的相互作用关系则由边值关系和
4、边界条件反映出来,称之为边值问题。(1)在介质的分界面上,电场满足的边值关系为且为电势所满足的边值关系:2121()0()nEEnDD 在介质分界面附近取两点1和2,而 所以由于 ,故 ,且介质2介质1212112lln1E2E2D1D02l112211221122()()nnE dlE n lEn lElEl 0,21ll021注意:可代替 ,即可代替证:可见而故有即得是连续的电势即在界面上,.21SS 21SSttEE1221()0nEE0 ,02121212112p2p1P1P2lElE222111 ,lElE21ttEE21另外,由方程 可得到:即 也就是说,在两种不同介质的分界面上,
5、电势 满足的关系为)(12DDn22112211()nEEnnSSnn1122SSSSnn112212 (2)在介质与导体的分界面上的情况 由于静电平衡条件,我们知道:导体内部 ;导体表面上的场强与表面导体是等势体;导体内无电荷分布(),电荷只分布在导体的表面上()。因此,在导体与介质的分界面上;00内E0导体1自由电荷介质2常数1即有归纳起来,静电场的基本问题是:求出在每个区域(均匀)内满足泊松方程,在所有SnnE2211 0 ,0 即导体内部内SS常数分界面上满足边值关系和在所研究的整个区域边界上满足边界条件的电势的解。3、利用静电标势来描述静电场的能量利用静电标势来描述静电场的能量 已知
6、在线性介质中静电场的总能量为在静电情形下,能量W可以用电势 和电荷 表出。由 得12WE DdVDE和)()(DDDDDE因此即 若我们考虑的是体系的总能量,则上式的体积分是对全空间进行的。因此上式右边第二项的面积分是对无穷大的面进行的。有限的电荷体系在无穷远处的电势 ,电场 ,而面积r2,故在r时,面积分项的值=0,故有 11()22WdVD dV1122SWdVD dS21rr1讨论:对 的使用注意几点:(1)适用于静电场,线性介质;(2)适用于求总能量(如果求某一部分能量时,面积分项 );(3)不能把 看成是电场能量密度,它只能表示能量与存在着电荷分布的空间有关。真实的静电12WdV12
7、WdV102SD dS21能量是以密度 的形式在空间连续分布,场强大的地方能量也大;(4)中的 是由电荷分布 激发的电势;(5)在静电场中,电场决定于电荷分布。在场内没有独立的运动。因而场的能量就由电荷分布所决定。(6)若全空间充满了介电常数为的介质,且得到电荷分布所激发的电场总能量DEw2112WdV1()()8xxWddVr式中r为 与 点的距离。4、举例讨论、举例讨论例1求均匀电场 的电势。Solution:因为均匀电场中每一点强度 相同,其电力线为平行直线,选空间任一点为原点,并设原点的电势为 。xx0E00Eyoxp0Ex根据 ,得到故得到 这里有个参考点选择问题这里有个参考点选择问
8、题例2均匀带电的无限长直导线的电荷线密度的,求空间的电势。Solution:pppldEpldEpp021)0()()()(12xEldEpp0000)(选取柱坐标:源点的坐标为(0,z),场点的坐标为(R,0),考虑到导线是无限长,电场强度显然与z无关。这里,先求场强 ,后求电势 。场点pRozz电荷源zdqdrE 由于电荷元为 ,因此令(0)(0)rzrzrRez eRez ezdqd330022 3 222 3 2001144114()4()rzzdqdzErrrrRez ez edzdzRzRzdRzdRz2sec ,tg且222222222322222322223222cossec1
9、secsec)tg1(sec)(RdRdRdRRdzRzd而故设p0点与导线的垂直距离为R0,则p点到p0点的电势差为0sinsectg)(22222322322dRRdRRzRzdz200242rrEReReR若选p0为参考点(即 ),则0)(0R00ln2)(RRR0000000ln2ln22)(2)()(000RRRRRdRRdzzRddRRl dERRRRRRRR2.2 唯一性定理唯一性定理Uniqueness theorem 本节内容将回答两个问题:本节内容将回答两个问题:(1)要具备什么条件才能求解静电问题)要具备什么条件才能求解静电问题 (2)所求的解是否唯一)所求的解是否唯一
10、1、静电问题的唯一性定理静电问题的唯一性定理(1)有介质存在的情况 把一个区域V找分为许多小区域Vi,每一个小区域内介电常数为 ,它是各向同性的。每一个区域给定电荷分布iSVkVkiVijsdisdjjV ijSVxx ,)(已知:在每个均匀区域中满足 ,即有几 个区域就是几个泊松方程。在各个均匀区域的交界面上,满足:至此,不知道边界条件,即不知道区域的边界S上的一些条件。这个问题正是唯一性定理所要解决的,下面讨论之。ii2jiijinn)()(,j唯一性定理:唯一性定理:设区域V内给定自由电荷分布 在V的边界S上给定 (i)电势或 (ii)电势的法向导数 ,则V内的电场唯一地被确定。SSn,
11、)(x下面采用的证法:下面采用的证法:证明:设有两组不同的解 和 满足唯一性定理的条件,只要让得 即可。令在均匀区域Vi内有 常数 0 ,22 ii在两均匀区界面上有在整个区域V的边界S上有或者 为了处理边界问题,考虑第i个区域Vi的界面Si上的积分问题,根据格林定理,对已知的任意两个连续 ,nnnnnnjjiijjiijjiijijiji 0 0 SSSSS0 SSSnnn函数 必有:令且 和2()()iiVSdVdSn i222 ()()0 ()iiiiiiiVSiiVSdVdSndVdS 对所有区域求和得到进一步分析:在两个均匀区域Vi和Vj的界面上,由于和 的法向分量相等,又有 ,因此
12、内部分界面的积分为2i()iiiiiVSdVdS ijdSdS ijijjiijjiiiiiijjjjSSSiiiijjjiSSdSdSdSdSdS (这里 )因此 故而在S面上,从而有0ijjijiiiijjiSSdSdSnn nnEEDDjjiijnjinijnin ,iiiiSSdSdS 2()iiiiSVdSd 0 ,0SSn或由于 ,而 ,只有 ,要使 成立,唯一地是在V内各点上都有即在V内任一点上,。由 可见,和 至多只能相差一个常数,但电势的附加常数对电场没有影响,这就是说静电场是唯一的。(2)有导体存在的情况2()0iiidV0)(2i00iiViid2)(0常数 讨论区域是导
13、体外空间V,即V是由导体外表面S1,S2及S包面所围成的空间,当S在无穷远处时,所讨论的区域就是导体外的全空间V。约定:在无穷远处,电场为零,即在S面上 或者表示成 在此基础上,把问题分为两类:A类问题:类问题:已知区域V中电荷分布 ,及所有0S0SVS1S2)(x 导体的形状和排列;每个导体的电势 都给定。B类问题:类问题:已知区域V中电荷分布 ,及所有 导体的形状和排列;每个导体的总电 荷都给定。因为导体面就是边界面,因此上述导体的电势或者总电荷就是边界条件。先用反证法证A类问题。证明:设存在着两个解 和 ,这意味着在区域V内,和 都满足泊松方程:)(x 第 i 个导体的表面为Si 面上,
14、该导体的电势为 。那么,在Si面上,和 都必须等于 。即在S面上,令 则有 应用格林定理:22 ,iiSiSii 0 00222iiSSSiii2()()VSdVdSn i令 ,有式中被积函数 ,要使上式成立,必然在V中每一点上有于是,V中每一点上,。22()iiVSSdVdSdSnn 22 0,0,0 ()0iSSVdV及00)(2常数 但在导体表面上,即得到常数=0,即 ,使得这就说明了对A类问题 有唯一解。再用反证法证再用反证法证B类问题类问题 也设存在两个解 和 ,则有令 代入格林公式中,得00 )(,02 22()iiVSSdVdSdSnn 因为在导体表面Si处,电势并没有给定,但根
15、据电磁学中的知识,导体在静电平衡时为一等势体。虽然 与 不一定相等,但对同一导体而言,故可从积分号内提出来,于是22 0,0,()iSiVSdVdSn即得iS)(应为一确定值iSSii iS2()iiiVSdVdSn 现在分析:因为 中,Si表示电场中第i个导体的表面,导体在静电平衡时,在导体外,紧靠导体表面处的场强方向与导体表面垂直,场强的大小与导体表面对应点的面电荷密度成正比,即从而得到?iSdSn iSdSn nEEnn ,/而)(1)(nnn这样就有式中 和 都表示第i个导体所带的总电荷,又因为它是给定的,即故对每一个导体表面都有此结论。因此得到11()iiiSSSdSdSdSnQQQ
16、Q 0iSdSn 2()0VdVQQ 同理,要使上式成立,必然是即由于 ,此常数对电场无影响,所以此时仍说 是唯一的。唯一性定理(另外一种证明方法)区域V由封闭面S0、S1、S2、等所包围,其中S0是最外包围面。如果V内的电荷密度 分布已知,并且各边界面满足下列条件之一时:0)(20常数 ES0VS1S2 (i)Si面上电势 已知;(ii)Sj面上为等势面。未知常数,并且Sj 面上流出的电通量已知。(iii)Sk面上的电场法线分量En已知。则区域V内电场强度被唯一确定。用反证法证明。证明:设有两上电势 和 ,它们都满足场方程xjcjjjjSNdSn 已知已知kSn 22 ,并满足上述边界条件,
17、则 ,或者 ,和 不必相等,可以相差一个常数,即 要证明场中每一点 成立,只需证明这里因为 ,并 。要使其等于0,则必须 。而由矢量恒等式 EE const 2()0VdV 0)(2 0d2()()()VVdVdVfff)(则有其中因为 所以 即也就是现在考察上式右边的面积分之值。2()()()()()()VVVdVdVdV 22 ,0)(2()()()()VVdVdV2()()()iiViSdVdS a)设Si面满足(i)类边界条件,则故Si面积分为零。b)设Sj面满足(ii)类边界条件,由于 ,故可以将 从积分号内提出来,则有 由于(ii)类边界条件中还包括有给定总通量值,即0)(iS1x
18、Scj未知数故 212)(,xxSxScccjj)()()()()()()jjjjjjjSSSnnjSSdSdSEEdS 从而使得 c)设Sk面满足(iii)类边界条件,则由于在Sk面上En值给定,故 则jjnjnjSSE dSE dS()()0jjSdS()()()()kkkkSSdSEEdS0)(kSnnEE由此可见,满足场方程组和边界条件的 和 必须满足等式即 ,唯一性定理证毕。2、用唯一性定理解决实际问题用唯一性定理解决实际问题例例1有一半径为a的导体球,它的中心恰位于两种均匀无限大介质的分界面上,介质的介质常数分别是 与 。若导体球总电荷为Q,求导体球表面处自由电()()0kkSdS
19、 2()0VdVEE 21荷分布。Solution:设导体球上下两半球各自带电量为q1和q2,则Q=q1+q2又因为导体球是等势体,上下半球电势相等,即Qa2121另外,总电荷Q一定,无限远处电势为0,故满足唯一性定理条件。根据唯一性定理,得到则得21222111 而ararrr222121221122 ,aqaq即故即得到:电荷面密度为:21212112211)(1qQqQqqqQqQq21222111 ,)(22 ,)(2221222222121211aQaqaQaq例例2两同心导体球壳之间充以两种介质,左半球介电常数为 ,右半球介电常数为 。设内球壳半径为a,带电荷为Q,外球壳接地,半径
20、为b,求电场和球壳上的电荷分布。12baS1S221rnSolution:以唯一性定理为依据来解本题。a)写出本题中电势 应满足的方程和边值关系以及边界条件 此区域V为导体球与球壳之间的空间,边界面有两个,即S1和S2,S1是导体球表面,S2是导体球壳内表面,边界条件为:在S1上总电量是Q,在S2上 。在两种介质中,电势都满足Laplace方程,在介质交界面上,电势 连续,电位移矢量的法向分量连续(因为交界面上 )。00f 应满足的定解条件为:现在不论用什么方法,只要求出的点函数 能满足上述条件,那么 就是本题的唯一解。b)根据已知的定解条件,找出电势 的解 由于对称性,选取球坐标,原点在球心
21、,直接积分0 ,)2,1(0212211212已知面上在已知面上在在交界面上SQSnnii)(x)(x可求得解,因为不难看出:在r=b处:0)(1222rrrrii)()(2211中电势右半球中电势左半球DrCBrAbABBbAbr 0 1从而得到同理,在r=b处:即得在两介质的交界面上:)11(1brAbCDDbCbr 0 2)11(2brC21由此得到 A=C又因为在两介质的交界面上,与 ,但 都只与r有关,所以 这样,也满足了Dn连续的条件。到此为止,在条件中,除了在S1面上总电量为Q外,也满足了其它全部条件,而 也只剩下一个待定常数A。现在用 必须满足在S1面上总电量等于Q这个条件来确
22、定A,即 rn021nn21,21,11111111111122222222221212 SrrrrSSSrD dSQADEeerrADEeerrD dSDdSDdSAAedSa 左右1122221222 22 rSSedSaAAaaaaQ左右故从而得到:c)电场和电荷分布情况 根据电势 所得到的结果,有)(221QA)11()(2)11()(2212211brQbrQi相应地,有321222321111)(2)(2rrQerErrQerErr32122223211111)(2)(2rrQEDrrQED由此可见在导体球(r=a)表面上:可见在导体球壳内(r=b)处:|21DD)(2)(2212
23、22222121111aQDDnaQDDnarrarfarrarfff21也可看出:还可进一步求出束缚电荷(极化电荷)分布:已知所以)(2)(221222222121111bQDDnbQDDnbrrbrfbrrbrf壳壳ff壳壳21EDP0321022022321011011)(2)()(2)(rrQEDPrrQEDP而极化电荷体密度:即在两种介质中,极化电荷体密度都为零。在导体球表面上极化电荷面密度分布:0)(10)(12222212211rprpPrrrPPrrrP)(2)()(2)(2120222221201111aQPPnaQPPnarrarparrarp故得到导体球表面上的总电荷 分
24、布:可见在两种介质交界面处:因为 。因而 ,所以注意:在前面计算过程中,难得出导体球面上)(2)(221202222120111aQaQpfpf210prnipifi0nP是常数,但是 或 在每个半球面上虽然都是常数,但 ,即 在球面上不是均匀分布的。现在来说明 不能均匀分布的原因。假定 是均匀分布的,那么由可见,在两个半球面上,因 值不同而不同。导体球内的静电场由 和 共同激发,由于 均匀分布,所以 在球内的电场为零。但 由于非pp21ifipff21ff000()pfn PnEnD ppfffpf均匀分布必将导致它在球内的场不为零,这样导体球就不能达到静电平衡。由此可见,要使导体球达到静电
25、平衡,的分布必须是非均匀的。f2.3 拉普拉斯方程,分离变量法拉普拉斯方程,分离变量法Laplaces equation,method of separate variation 本节内容主要是研讨本节内容主要是研讨Poisson 方程的求解析方法。方程的求解析方法。众所周知,电场是带电导体所决定的。自由电荷只能分布在导体的表面上。因此,在没有电荷分布的区域V里,Poissons equation 就转化为 Laplaces equation,即产生这个电场的电荷都是分布于区域V的边界上,它0 22们的作用通过边界条件反映出来:给定 给定 或导体总电量因此,讨论的问题归结为:怎样求解(通解)L
26、aplaces equation.怎样利用边界条件及边值关系求出积分常数。Laplaces equation可以用分离变量法求通解,其求解条件是:方程是齐次的。边界应该是简单的几何面。SSdSQn n1、用分离变量法求、用分离变量法求Laplaces equation的通解的通解 (1)在直角坐标系中)在直角坐标系中设在数学物理方法中,该方程的通解的 (A、B、C为待定系数)02222222zyx)()()(),(zZyYxXzyx)sincos()sincos()sincos(),(212121zkCzkCykBykBxkAxkAzyxzzyyxx或者写成 (2)在柱坐标系中)在柱坐标系中设
27、该方程的通解为)(;),(222yxzzikyikxikkkkeeezyxzyx01)(1222222zrrrrr)()()(),(zZrRzr)sinh()cosh()sin()cos()()(),(212121kzCkzCnBnBkrNAkrJAzrmm其中,Jm为m阶第一类贝塞尔函数,Nm为m阶第二类贝塞尔函数。如果考虑与z轴无关(k=0)情况,并讨论的区域是 ,故通解为)sin()()()cos()()()1(!)2()1()(02mkrJkrJmkrNnmnkrkrJmmmnnmnm为伽马函数20110)sin()()cos()(ln),(nnnrnnnnnnonrDrCnrBrAr
28、BAr这里A,B,C,D为待定系数。(3)在球坐标系中)在球坐标系中设其通解为0sin1 )(sinsin1)(122222222rrrrrr),()(),(YrRrmnmnnnmnnmmnmnnnmnnmmPrDrCmPrBrAr,1,1)cos()(cos)()cos()(cos)(),(这里 为缔合勒让德(Legendre)函数 对于具有轴对称的问题,m=0(取此轴为极轴)且这里 为勒让德函数,、为待定系数 对于球对称的问题,m=0,n=0。且2、利用边界条件定解利用边界条件定解 说明两点:)(cosmnP01)(cos)(),(nnnnnnPrBrAr)(cosnP为待定系数BArBA
29、r,)(nAnB 第一,如果考虑问题中有i 个区域(均匀分布),必须有i个相应的Laplaces equation.第二,在每个区域的交界面上,应该满足边值关系:边界条件:及导体的总电荷)(面上在ijjjiijiSnnSdSQn SSn或3、举例说明定特解的方法举例说明定特解的方法例例1一个内径和外径分别为R2和R3的导体球壳,带电荷为Q。同心地包围着一个半径为R1的导体球(R10的区域,空间一点的电势为在x0区域电势为:xRo)处有一点电荷Q,求空间的电势分布。Solution:取球心为坐标原点,球心到点电荷Q的方向为x轴,设Q的坐标为(a,0,0)。根据静电平衡条件(现象)。球内的电势为零
30、。故只讨论外空间的电势即可。aQRo 球外空间的电势由Q及球面上感应电荷共同激发的,其电势所满足的定解条件为:用一个象电荷Q来代替球面上的感应电荷,为了不改变原方程,Q必须在球内,并距球心为b,故等效为:(3)0(2)0(1),(1002RRRzyaxQ球外空间一点的电势为Rob QaQxrrP(x,y,z)212222122200)()(41)(41zybxQzyaxQrQrQ 在b R0 ,即近似为两点电荷作用,作用力为排斥力;当Q靠近球面时,此时不论q与Q是否同号,作用力永远为引力,这可由在Q附近的感应电荷与其反号来解释。xeRaaRaQRqaQ22022023020)()2(41xea
31、QqF20410RaooaaQQQQ点电荷Q在球内点电荷Q在球外bb镜象法与光路图比较例例4均匀场中的导体球所产生的电势 由于静电屏蔽,场区域只能在球外。Solution:本题的物理图象是在原有的均匀电场 中放置一中性导体球。此时导体球上的感应电荷也要在空间激发场,故使原来的场空间电场发生了变化,如图所示。由此可见,球外空间任一点的场将是一个均匀场和一个球体感应电荷等效的偶极子的场的迭加。0ER0-+0E 第一步:第一步:用两个点电荷Q激发一均匀场 点电荷Q放在对称轴z=a处,a很大,Q也很大,在坐标原点附近的区域内。第二步:第二步:将一中性导体球放在均匀场中0E+Q-Qz0EaaoEE020
32、0241aQE这样一来,Q相当于两个场源电荷,球面上将出现感应电荷,由象电荷来代替它,即 此时+Q在球面上感应的电量为 ,-Q在球面上感应电量为 ,这仍然保持导体球为电中性(不管导体球接地与否)。根据唯一性定理,导体球外的+Q-QzaaR0bboQaR0QaR0)(0QaR)(0QaR电势就是这四个点电荷分别在某点产生的电势的迭加,即因为aR,则选 略去 和212024020212024020212221220)cos2()cos2()cos2()cos2(41aRRaRRQaRaRRaRRQaRRaaRQRaaRQ外22aRbaRR240224022RaRRb即又因为 皆为小量,应用展开式2
33、12002120021210)cos21()cos21()cos21()cos21(14aRRaRRaRRaRRaRaQaRaQ外aRRRaRaR2020和)(211)1(21为小量xxx则有coscoscos2414cos2cos2cos24)cos1()cos1()cos1()cos1(402300223000230020020000ERRRERaQRaRQaRRaRaQRaRRRRaRaRaRaRaQ外令则 的第一项恰好等于一个原均匀场以o点为参考点电势。第二项恰好等于位于o点的电偶极矩为 的电偶极子的电势。(3)界面为柱面的情况)界面为柱面的情况例例5有一线电荷要密度为的无限长带电直线
34、与半径为R0的接地无限长导体园柱轴线平行。直线与园柱轴zeREP03004)(40300电偶极场即原场外ERRPzEP外线的距离为a(aR0),试求空间的电势分布。Solution:由于导体柱面把整个空间分成柱内、柱外两个区域,而柱内有 ,柱外区域电势满足定解条件:yaaxR0R00内 处于带电直线的电场中的导体园柱,其柱面上要出现感应电荷,空间任一点的电势 就是带电线和感应电荷分别产生的电势的迭加。现在,假定导体园柱面的感应电荷密度为 ,到轴线的距离为b,由于原带电直线不仅带电(均匀)而且是无限长的,导体园柱也是无限长的,故垂直于柱轴的任何平面上的电势分布是完全相同的,即是一个二维场,因此可
35、取一个垂直于柱轴的平面来讨论(2)0(1),(1002RRyax外外外即若取oa连线与圆柱面的交点为电势参考点。则园柱外空间任一点的电势为RobaxrrP(x,y)R0bRrRar0000ln2ln2其中由(2)式得2122212221222122cos2)(cos2)(RbbRybxrRaaRyaxr0)cos2(ln2)cos2(ln202102200021022000aRbRbRRaaRaRRR即要使该等式成立,必有由(4)式,即有bRbRbRRaaRaR02102200210220)cos2(ln)cos2(ln(4)cos2()cos2(3)02102200210220bRbRbRR
36、aaRaR200220200220)(cos2()(cos2(RabRbRbRaRaR比较两边系数,即由(6)式得化简(7)式得到:解这下一元二次方程得到(6)(cos2)(cos2(5)()(2002002022020220RabRbRaRRabRbRaR(7)2020aRabbR0)(202022aRRababaRbab2021 ,其中 b1=a不符合物理要求。故有:因而柱面外任一点的势为相当于球面距离情况相当于平板时电荷情况 20aRb)(lnln2ln2ln20000200000RaRraRaraRRrRar(4)界面为劈形的情况)界面为劈形的情况例例6有两个相交的接地导体平面,其夹角
37、为 ,若在所夹区域内有一电量为Q的点电荷,求下列情况下所夹区域内的电势:rRraRarRaRra000000ln2)(ln2Solution:从上面的例子可以看出,用镜象法处理问题时,只要象电荷都放在所考虑的区域之外,就不会改变电势在该区域内所满足的泊松方程。故检验解是否正确的关键是看它能否满足全部边界条件。4 c);3 )b ;2 )aQ 下面按夹角 不同情况分别讨论其电势分布情况。a、2APB-QQ-QQorr1r2Rr3132 所考虑的区域内,势满足定解条件。为了使A板的电势为零,应在以A板为对称面,将A板上的感应电荷以象电荷-Q放置在与源电荷Q对称的位置“1”处,要使B板的电势为零,应
38、以B板为对称面,将B板上的感应电荷以象电荷-Q放置在与源电荷Q对称的位置“2”处,而且还需在“1”相对于B板的对称位置“3”处放置+Q的象电荷,才能保证 ,不难看出,此时也满足 于是所考虑区域内任一点(2)0(1)(102BARxxQ0A0B的电势为 b、)1111(42130rrrrQ3BAQ+Q+Q-Q-Q-Q12345 要保证上 则必须有5个象电荷,其位置,大小和符号如图示,于是所求区域内电势为 c、0 ,0BA)111111(4543210rrrrrrQ4BAQ+Q+Q-Q-Q-Q123452-Q+Q要保证 则必须有7个象电荷,故电势为 一般说明:只要 满足 偶数的情形,都可用镜象法求
39、解,此时象电荷的个数等于 ,加上原来的电荷总共有 个,这些点电荷都在过原点电荷与两导体面的交线垂直面内。而且都在此垂面与交线的交点为圆心,交点到原点电荷处的距离为半径的圆周上。若 不满足该条件,则象电荷在所求区域内,改变了原方程,否掉。0BA)11111111(476543210rrrrrrrrQ21)2(22.5 格林函数法格林函数法Method of Green function 本节要介绍的是一种用Green定理来求解静电边值问题的方法。即给定区域V内电荷分布 ,和区域V的边界面S上各点的电势 或电势法向导数 ,求区域V内各点的电势值。如果边界条件是给定S上的电势 ,这类边值问题称为第一
40、类边值问题,也称狄利克莱边值问题;如果边值(界)条件是给定S上的 ,这类边值问题称为第三类边值问题,也称诺埃曼边值问题。s)(xSnsSn 在这里,我们将要讨论这些边值问题是怎样借助于有关点电荷的较简单的边值问题而得到解决的。1、点电荷密度的点电荷密度的 函数表示函数表示 因为点电荷分布的特点是在点电荷所在处的电荷密度变为无穷大,而在其他地方电荷密度为零。若在 处有一点电荷Q,则电荷密度可写为显然x(1)0)()(xxxxxxQxVVQdxxQdx)()(对于单位点电荷而言,Q=1,其密度为2、Green函数函数 一个处在 点上的单位点电荷,它所激发的电势方程为假设有一包含 点的某空间区域V,
41、在V的边界S上 有如下边界条件(2)()(xxxx(3)(102xxx(4)1 00SSS或者则把满足边界条件(4)式的(3)式的解称为泊松方程在区域V的第一类或第二类边值问题的Green函数。Green函数一般用 表示,表示单位电荷所在的位置,代表观察点,在(3)式和(4)式中,把 换成G,即Green 函数所满足的方程和边界条件为3、Green公式和边值问题的解公式和边值问题的解x),(xxGx(5)1)(,0)()(1)(002SnxxGxxGxxxxGSS或 在这里,将用Green公式把一般Poisson方程的边值问题的解用Green 函数联系起来。(1)先看)先看Green公式的两种
42、形式公式的两种形式 根据 Gauss 定理,知道当 均为连续,可微的标量点函数,故()nVSSa dVa dSa dS和 ,a2)(a又于是,有式中V为包围面S所围的面积,该式称为Green第一公式。如果上式中的 对调,即 ,同理得到 naa nn 2()()(6)VSdVdSn a和2()()(7)VSdVdSn 将(6)式减去(7)式,得该式称为Green第二公式 Green第一、第二公式是等价的。同时,视方便而选取之。Green公式对解静电问题的意义是:在区域V内找一个待定函数 (为待求),通过这个公式从已知确定未知。(2)边值问题的解)边值问题的解 给定一个区域V,其中给定了22 (8
43、)VSdVdSnn SSnx ,)(且待求的边值问题:相应的Green函数问题是:边界条件:现在,取 满足)(102x)(1)(02xxxxGV给定了S)(xSnGGSS01 0或)(102x 取 满足代入Green第二公式,有因为Green公式中积分,微分都是对变量 进行的,由于Green函数关于源点和场点是对称的,即 ,为方便起见,把变量 换为 ,故有 改为 ,即得)(1),(02xxxxG22(,)()()(,)()()()()VSG x xxxG x xdVxG x xG x xxdSnn )(xxG)(xxG xxx该式左边第二项为得到0011(,)()()()()()()()VSG
44、 x xxxxxdVxG xxG xxxdSnn0011()()()Vxxx dVx0011()()()()()()()VSG xxx dVxxG xxG xxxdSnn故得到这就是用Green函数求解静电问题的一种形式解。讨论几点:a)在区域V中,任一点的势唯一地决定 电荷分布及边界的值0()()()()()()()VSxG xxx dVxG xxG xxxdSnn)(xSSn或 b)如果所取的Green函数属于第一类问题,即 这时则有这实质上就是第一类边值问题的解 c)如果所取的Green函数属于第二类问题,即在这里要说明一点的是:对第二类静电边值问题不能用第二类齐次边界的Green函数,
45、即 ,因0)(SxxG0()()()()()VSG xxxG xxx dVxdSnSnGS010SnG为 Green函数 所代表的物理意义是在 处存在一个单位电荷在空间所激发的电势。因此即代表单位电荷在边界上所激发的电场,由Gauss定理知道由此可见故)(xxGx0SnG01()SG x x dSn()0SG x x dSn0)(SnxxG从而,Green函数在边界上的最简单的形式是取这样且有第二类静电边值问题的Green函数解的形式:式中 为 在边界面S上的平均值。SnxxGS01)(000()1()()()()()()1()()()()VSVSSxxG x xx dVG x xxdSnSx
46、G x xx dVG x xdSx dSnS 1()Sx dSS 在实际问题中,常遇到这类问题:在所考察的区域包含有无穷大的边界面,假如,考察一导体球外的空间电势分布问题,这时所考察的区域是球面和无穷大曲面间包围的区域,所以这时边界面S故有于是故得到01)(0 SnxxGS1()0Sx dSS 0()()()()()VSxxG x xx dVG x xdSn 此式称为外问题的Green函数解的形式。4、Green函数的制作函数的制作 以上的讨论,表面上制裁似乎把静电边值问题的解找到了,其实并作为此,因为只有把问题的Green函数找到了,才能对表达式(第一类边值问题的形式解和第二类边值问题的形式
47、解)作出具体的计算。实际求Green函数本身并不是件容易的事,所以以上解的形式只具有形式解的意义。当然,它把唯一性定理更具体地表达出来了。在这里介绍几种不同区域的Green函数的制作方法。(1)无界空间的Green函数 即在无穷大空间中放一个单位点电荷,求空间某处的电势,也就是Green函数。其中,代表单位电荷的所在位置(源点坐标),代表观察点坐标(场点坐标)。现在,证明上述Green函数是否满足Green函数所满足的微分方程。证明:选电荷所在处为坐标原点,即 ,在球坐2122200)()()(141)(141)(zzyyxxxxxxGxx0 x标系中考虑球对称性,得到而当r=0时,取一小球面
48、S 包围着原点,取 对小球体积V积分,即rxG041)0(rxG141)0(202)0(012rrr12从 函数性质可知,保持小体积V的面积为1,从而有故得到223211114VVSSSdVdVdSrrrrdSr drr 211()14VVxx dVdVr)0()(412rxxr与微分方程比较,即有这里把 与 互换,不变,即有这就说明Green函数具有对称性。(2)上半空间的Green函数 即在接地导体平面的上半空间,由于 ,属于第一类边值问题。)(1141)(0202xxrxxGxx)()(xxGxxG)(xxG0SG 根据镜象法得到:yzor2r1xx21222212220)()()(1)
49、()()(141)(zzyyxxzzyyxxxxG这也可看到(3)球外空间的 Green函数 即在接地导体示外的空间,由 ,属于第一类边值问题。)()(xxGxxG 0SGyzxRRR0 xxrro其中:根据镜象法得212040220212222222222cos21|)(|cos2|)cos(sinsincoscoscos RRRRRRRxRRxrRRRRxxrzyxRzyxR21202021220)cos2)(1)cos2(141)(RRRRRRRRRRxxG 在制作Green函数时,必须注意:求Green函数本身不是很容易的,只有当区域具有简单几何形状时才能得出解析的解,如果 时,Gre
50、en函数法也可以用来解Laplace equation的边值问题。5、Green函数法的应用举例函数法的应用举例例例 在无穷大导体平面上有半径为a的园,园内和园外用极狭窄的绝缘环绝缘,设园内电势为V0,导体板其余部分电势为零,求上半空间的电势。0)(xSolution:静电问题:axyzRP(,z)P(,z)V0此题Green函数满足的形式为相当于无穷大金属平板旁边放置单位电荷求电势问题00 0000002RzRzVz0),()(1),(002zxxGxxxxG其 Green函数为其中:换为柱坐标,且有故Green函数为)11(41)(0rrxxG2222)()()(zzyyxxr222222