静态场的边值问题

静态场的边值问题边值问题研究方法解析法数值法分离变量法镜像法复变函数法有限差分法有限元法边界元法矩量法模拟电荷法第2页,共70页，2024年2月25日，星期天§5.1唯一性定理和解的叠加原理

一.唯一性定理1、表述在给定的区域内，泊松方程（或拉普拉斯方程）满足所给定的全部边界条件的解是唯一的。2、边界条件的形式①给定全部边界上的函数值②给出全部边界上函数的法向导数值③给定部分边界上的函数值，而其余边界上给出函数的法向导数值“狄利赫利”边界条件“聂曼”边界条件混合边界条件第3页,共70页，2024年2月25日，星期天3、唯一性定理的证明证明：考虑泊松方程，用反证法在区域τ内存在两个不同的函数和都满足相同的泊松方程，并且在区域边界S上满足同样的边界条件。令则有利用取则对上式两边在区域内作体积分，然后运用散度定理，得第4页,共70页，2024年2月25日，星期天将代入上式得①对于第一类边界条件②对于第二类边界条件③对于第三类边界条件不论对哪类边界条件，面积分都等于零因此有

第5页,共70页，2024年2月25日，星期天由于恒为正值，故上式成立的条件是解之可得讨论：①对于第一类边界条件，因为所以C=0②第二类和第三类边界条件的情况对于求解场函数来说，解是唯一的。4、应当明确①只有在区域的所有边界上给出唯一的边界条件时，边值问题的解才是唯一确定的。②唯一性定理给出了求解电磁场问题的理论依据不论采用什么方法，只要得到的解能够在区域内满足方程而在边界上满足边界条件，这个解就是该边值问题的唯一正确解。第6页,共70页，2024年2月25日，星期天二.解的叠加原理解的可叠加性是方程线性的必然结果。1、表述①对拉普拉斯方程，若、、……都是满足方程的解则其中ai为任意常数。也是方程的解并且、、……都是满足方程的解则其中ai为任意常数。也是方程的解②对泊松方程，若是满足方程的一个任意解第7页,共70页，2024年2月25日，星期天2、证明讨论泊松方程的情况对叠加得到的结果两边作运算，得

因此是方程的解令f=0，即可得到拉普拉斯方程情况的证明3、应用求解边界问题时，可以先将复杂边界条件分解成便于求解的几个边界条件，则总的边界问题解就是这些解的叠加。第8页,共70页，2024年2月25日，星期天例：分解为三个边界问题分解后每个边值问题都只有一个非齐次边界值，求解变得容易。原问题的解应该是三个问题解的叠加第9页,共70页，2024年2月25日，星期天§5.2拉普拉斯方程的分离变量法一.直角坐标系中的分离变量法

直角坐标系中拉普拉斯方程的表达式为令，并代入上式并两边同除以得1、方法介绍①变量的分离第10页,共70页，2024年2月25日，星期天则上式分解成三个独立的全微分方程，即 称为分离常数，分离常数之间满足约束关系②全微分方程解的选取（以为例）第11页,共70页，2024年2月25日，星期天对和也都有与上述相同形式的解。在、、三组可取解中各取其一并相乘，即可得到一个解的表达式。★解的选取并不是任意的，因为存在约束条件如果取，，则必须取a.对于有两个零值边界的方向，其对应的函数取三角函数；b.对于单零值边界方向，对应的函数一般取双曲函数形式；c.而有无限远边界的方向，一般取指数函数形式。d.若位函数与某一坐标变量无关，则该变量对应的函数取成常数，并取作1。★根据边界条件来选择函数：第12页,共70页，2024年2月25日，星期天③拉普拉斯方程的通解本征值：满足齐次边界条件的分离常数可以取一系列特殊值本征值对应的函数称为本征函数或本征解。所有本征解的线性叠加构成满足拉普拉斯方程的通解在许多问题中，单一本征函数不能满足所给的边界条件，而级数形式的通解则可以满足单个解函数所无法满足的边界条件。第13页,共70页，2024年2月25日，星期天2、例题例5.1边界条件如图，求电位分布

o

y

0=U

b

a

0UU=图5－1长方形截面的导体长槽

0=U

0=U

x

解：内部电位与z无关，是二维问题

①根据边界条件写出通解x方向电位有两个零值边界，应取三角函数形式；y方向电位为单零值边界，应取双曲函数形式。第14页,共70页，2024年2月25日，星期天令分离常数则电位通解为②求待定系数将边界条件(a)代入得若对任意y成立，则再利用边界条件(b)得通解变为若对任意y成立，则第15页,共70页，2024年2月25日，星期天再利用边界条件(c)得若对任意x成立，则其中，若用整数n代替i，则上式表示为把边界条件(d)代入上式，得将U0

在（0，a）区间内展开为的傅立叶级数于是其中第16页,共70页，2024年2月25日，星期天对比两边各项系数可得因此，在槽内的电位分布是

0UU=

0=U

0=U

0=U

第17页,共70页，2024年2月25日，星期天例5.2有两块一端弯成直角的导体板相对放置，中间留有一小缝，如图所示。设导体板在x轴和z轴方向的长度远远大于两导体板间的距离b，上导体的电位为U0，下导体接地。求两板间的电位分布。解：电位分布与z无关，这是一个二维拉普拉斯问题

①边界条件第18页,共70页，2024年2月25日，星期天②边界条件的分解

0UUbya==

ybUUxa00==

ybUUUxb000-==

ybUUxb00-==

0=¥=xbU

(A)

(B)

图5－4图5－3的场分解

y

y

x

x

00==yaU

0==bybU

00==ybU

第19页,共70页，2024年2月25日，星期天对于(A)问题，情况与平行板电容器相同，两极间为匀强电场对于(B)问题沿y轴方向存在两个零电位边界，取沿x轴正方向边界无界且电位趋于零，取因此(B)问题的电位通解应为将边界条件(c)代入，得第20页,共70页，2024年2月25日，星期天利用傅立叶展开并比较系数，可求得(B)问题的解③利用叠加原理所以原问题的电位表达式为原问题的解应该是A问题和B问题解的线性叠加第21页,共70页，2024年2月25日，星期天二.圆柱坐标系中的分离变量法1、方法介绍①变量的分离圆柱坐标中拉普拉斯方程为令，代入上式，并在两边同乘以得得上式第二项只与有关，可以先分离出来，令其等于常数将上式代回到式(A)，各项同除以r2，得(A)第22页,共70页，2024年2月25日，星期天令上式左边最后一项等于常数，则上式分离成为两个方程②解的选取

时，

对于对于

时，

时，

时，

时，第23页,共70页，2024年2月25日，星期天对于，是一个n阶贝塞尔方程时

时

时

时

称为n阶第一类贝塞尔函数称为n阶第二类贝塞尔函数，也叫聂曼函数第一类和第二类变形贝塞尔函数其中③本征解的叠加构成通解第24页,共70页，2024年2月25日，星期天例5.3电场强度为的均匀静电场中放入一半径为a的电介质长圆棒，棒的轴线与电场相垂直，棒的电容率为，外部电容率为，求任意点的电位。解：在柱坐标系中求解①根据边界条件，写出通解与z无关，取，则故对于，根据问题的形式，可知所以取n2>0，并且n为正整数只取项第25页,共70页，2024年2月25日，星期天对于，因为所以取于是通解可以表示为②求分区通解r<a时，r=0处电位为有限值，故此区域通解为r>a时，总电位包括外电场和介质棒两部分的贡献外电场电位介质棒电位规定了考虑到当时，介质棒产生的电位是有限值，故第26页,共70页，2024年2月25日，星期天③利用静电场的边界条件求系数将U1和U2的通解表达式代入上面两个边界条件，得比较上面两式两边的系数，得n=1时以上两式联立求解，得第27页,共70页，2024年2月25日，星期天n>1时联立解上面两式求解，得④求出圆柱内外电位⑤求出圆柱内外电场强度第28页,共70页，2024年2月25日，星期天

中的无限长圆柱体是和上面例子完全类似的磁场在均匀外磁场

边值问题。用和分别表示圆柱体内外的磁标位，则它们的边界

条件也和电位的形式相同，即

时,

时,

此时，和的解也与静电场和的解有相同的形式，只需把和解中的和

分别用和代替，用代替，便得到和的解。同样可以得知圆柱体内的磁场强度是均匀的。第29页,共70页，2024年2月25日，星期天

例5.4一导体圆筒的高度为b，半径为a，所给边界条件为

求圆筒内的电位分布函数。解：由问题对称性可知，电位与变量无关，因此应取和;应该取成在零点处为零值；并考虑到电位在r=0处为有限值，故的形式为。因此电位函数的本征解为：的双曲正弦函数，有将边界条件代入上式，得所以应是零阶贝塞尔函数的根第30页,共70页，2024年2月25日，星期天

把所有本征值所对应的本征解相加，得到电位函数的通解将边界条件代入上式，得上式两边同乘以，从0～a对r积分，得

对上式两边分别应用贝塞尔函数积分公式

和正交公式（5－40），可以得到第31页,共70页，2024年2月25日，星期天

所以

例5.5

将例题5.4的边界条件改成求圆筒内的电位分布函数。解：由问题对称性可知，电位与变量无关，因此应取和应该取成在零点处为零值的正弦函数

并考虑到电位在r=0处为的形式为。因此电位函数的通解为有限值,

第32页,共70页，2024年2月25日，星期天

将边界条件代入上式，得

上式两边同乘以，从0～b对z积分，由三角函数的正交性，得代回原式，得到第33页,共70页，2024年2月25日，星期天三.球面坐标系中的分离变量法

球面坐标系中拉普拉斯方程的表达式为

令代入上式，并两边同乘以①变量的分离得引入-m2分离引入n(n+1)分离r分离θ第34页,共70页，2024年2月25日，星期天②解的选取

时对于

时

对于解只有一种形式

时

时

时对于分别称为第一类和第二类连带勒让德函数

分别称为第一类和第二类勒让德函数

③本征解的叠加构成通解第35页,共70页，2024年2月25日，星期天例5.6在均匀电场中放入一半径为a的接地导体球。求任意点的电位和电场强度解：①取球心位于坐标原点，电场方向为极轴方向，建立球坐标系②根据边界条件，求出通解与无关，取m=0，所求场域包括，故只含有第一类勒让德多项式

r

0Ev

x

a

0e

P

z

y

o

q

电位通解可以写成第36页,共70页，2024年2月25日，星期天③根据边界条件确定电位因为导体球接地，所以球内在球外，电位包括均匀电场和导体球两部分的贡献均匀电场的电位因为感应电荷的电位因为时，故总电位利用导体球的边界条件，得不含rn

项第37页,共70页，2024年2月25日，星期天由勒让德多项式的正交公式可得④根据电位求电场因此球外电位是可见，感应电荷对场的贡献相当于一个沿z轴放置的电偶极子，这是由于球面感应电荷的分布恰好是上正下负之故。例5.6略第38页,共70页，2024年2月25日，星期天§5.3镜像法镜像法理论依据：唯一性定理。镜像法基本思路：分界面上的感应源区域外的镜像源无限区域的同种媒质问题分区媒质问题镜像电荷位置选择原则：1、镜像电荷必须位于求解区域以外的空间。

2、镜像电荷的引入不能改变原问题的边界条件。

第39页,共70页，2024年2月25日，星期天一.平面镜像

例5.8真空中一点电荷q位于一无限大接地导体平面的上方，与平面的距离为h。求z>0区域的电位分布。1、导体——介质（电场）用镜像电荷代替导体平面上的感应电荷将区域换成真空解：①给出等效问题z>0区域所满足的边界条件保持不变，即所以，原问题就化作求解电荷在无界真空区域中的问题。区别仅在于我们只取z>0区域的解。第40页,共70页，2024年2月25日，星期天②求解等效问题空间任意点的电位由q和共同产生根据，可知R→∞时，U=0根据，可得

于是③原问题的解第41页,共70页，2024年2月25日，星期天感应电荷密度为

总感应电荷为

④导体表面上的感应电荷总感应电荷恰好等于镜像电荷电量。这一结果是合理的，因为点电荷q所发出的电力线将全部终止于无限大的接地导体平面上。第42页,共70页，2024年2月25日，星期天★点电荷对相交接地平面导体边界的镜像a.两半无限大接地导体平面垂直相交。要满足在导体平面上电位为零，则必须引入3个镜像电荷。如图所示。b.对于非垂直相交的两导体平面构成的边界，若夹角为，则所有镜像电荷数目为2n-1个。第43页,共70页，2024年2月25日，星期天2、介质——介质（电场）

q¢¢

h

R¢¢

),,(zyxP

2e

o

(c)

2e

y

x

q

q¢-

h

R

R¢

),,(zyxP

1e

1e

y

x

o

(b)

-h

例题5.9在1区距离界面h处有个点电荷q

求空间的电位分布q

h

1e

2e

y

x

o

解：①给出等效问题x>0区域，可化作是在充满介质的无界空间中的场问题(b)x<0区域，可化作是在充满介质的无界空间中的场问题(c)第44页,共70页，2024年2月25日，星期天②求解等效问题

两个区域的电位表达式为电位的边界条件是

将电位表达式代入得第45页,共70页，2024年2月25日，星期天联立求解，得

化简得将这两个镜像电荷的表达式代入分区域的电位表达式中即可得到整个空间的电位。③给出原问题的解第46页,共70页，2024年2月25日，星期天3、介质——介质（磁场）例5.10求与分界面平行的无限长线电流I，在空间各点的磁场。解：①原问题的等效1区2区第47页,共70页，2024年2月25日，星期天②求解等效问题设分界面法线方向，切线方向则界面两侧的磁场为利用边界条件将磁场表达式代入，可解得第48页,共70页，2024年2月25日，星期天③给出原问题的解如果区域1是的非磁性材料，区域2是的理想磁导体

则在x=0的分界面上表明：理想磁导体的外侧磁场只有法线分量。这与理想电导体表面只有法线电场分量的情况类似。第49页,共70页，2024年2月25日，星期天二.球面镜像

对于分界面是球面，并且源为点源的静态场问题，可用球面镜像。例5.11在半径为a的接地导体球外M点有一个点电荷q，球心O与M点的距离为d，如图所示。求导体球外的电位分布和球面上的感应电荷。

q'q’Oqr'RRd'dABM解：①设置镜像电荷镜像电荷应在OM连线的球内部分上，设的位置点与O点的距离为

②空间任意点的电位第50页,共70页，2024年2月25日，星期天边界条件对于球面上的任意点都成立考察两点，则有由上面两式解得

求出的镜像电荷及位置是否满足整个球面的零位边界条件？

将

代入表达式，得

第51页,共70页，2024年2月25日，星期天④故这两个点电荷在球外区域产生的位就是原问题的解，即

第52页,共70页，2024年2月25日，星期天④求感应电荷面密度

利用电位可求出导体球面上的感应电荷密度与总感应电荷

a.总感应电荷恰等于镜像电荷b.感应电荷的绝对值小于施感电荷q

表明q发出的电力线一部分终止于导体球面而另一部分则终止于无穷远处。第53页,共70页，2024年2月25日，星期天★讨论：a.导体球不接地且表面上不带过剩电荷b.导体球不接地，并且给出它的电位为U0

c.导体内挖一个球形空腔，空腔内点有一点电荷距球心此时需要在球心处增加一个镜像电荷，并且，新电荷系统由共同组成。此时需要在球心处增加一个镜像电荷，并且，新电荷系统由共同组成。此时它的镜像应该放在腔外的M点上，也就是本例问题的反演，镜像电荷和。腔内的场分布由共同确定。第54页,共70页，2024年2月25日，星期天例5.12假设一个无限大接地导体平面上有一半径为a的半球形导体凸块，在凸块附近有一个点电荷q。求此电荷的镜像。解：①建立坐标系设电荷q和导体平面法线所在的平面为xz平面。o

z

a

z

q¢-

d¢

q¢

q(x,0,z)

-q(x,0,-z)

x

a.电荷q对导体平面xy面的镜像电荷-q，坐标为(x,0,-z)b.电荷q对球面的镜像为②求镜像电荷导体外任意点的场由四个点电荷共同确定。c.镜像电荷对xy平面的镜像为，处在与O的连线上。第55页,共70页，2024年2月25日，星期天三.圆柱面镜像

r

P

o

a

lr¢

R¢

R

lr

j

x

d

o¢

d¢

M

y

b

b

B

例5.13在半径为a的无限长接地导体园柱外有一根与圆柱轴线平行的无限长线电荷，线电荷与圆柱轴线的距离为d，如图所示。求：柱外任意点的电位和柱面上的感应电荷。解：①设置镜像电荷根据场的对称性，可设镜像电荷是一条与圆柱轴平行的线电荷，线密度为，与轴线的距离为。②求解等效问题空间任意点的电位为

第56页,共70页，2024年2月25日，星期天其中

代入边界条件，得上式应对任意都成立，即圆柱面上的电位处处为零。因此应有

由此可以得到上式成立的充分条件是两个方括号部分都等于零第57页,共70页，2024年2月25日，星期天

所以可得到不合理正确解③给出原问题的解利用求得的镜像电荷参数可以得到柱外任意点电位④柱面上的感应电荷面密度和单位长度上的感应电荷分别为

第58页,共70页，2024年2月25日，星期天⑤等量异号平行线电荷的等电位面

lr¢

lr

b

b

y

x

图5—16

平行线电荷的等电位线

假定两线电荷相距为2b，电量分别为，则空间任意点电位为分别表示场点与的距离，可见当取不同值时，就得到不同电位的等位线。若取两线电荷连线为x轴，连线的中垂线为y轴建立直角坐标系，则整理可得第59页,共70页，2024年2月25日，星期天这是一个以常数k为参量的圆族方程，它表示两条平行异号线电荷在二维平面内的等电位线族。等位圆的圆心在半径为k<1时，等位圆在y轴的左侧，电位为负值；k>1时，等位圆在y轴的右侧，电位为正值；k=1时，等位线变为y轴，电位为零；k=0时，等位圆缩为(0,－b)处的点，电位为－∞；k=∞时，等位圆缩为(0,－b)处的点，电位为∞。第60页,共70页，2024年2月25日，星期天例5.14两无限长平行圆柱导体的半径都等于a，轴线之间的距离为2d，如图所示。求：导体柱单位长度的电容。解：用两条平行异号线电荷和作为平行带电圆柱的镜像。首先来确定线电荷的位置b。在右边圆柱边界上选取两个特殊点，设y轴上的任意点为零电位参考点，则空间任意点处的电位为则第61页,共70页，2024年2月25日，星期天即解之得所以空间任意点处的电位为可得带负电圆柱的电位由此将带负电圆柱面的方程代入上式，

第62页,共70页，2024年2月25日，星期天同理，可证明带正电圆柱的电位为两圆柱间的电位差两圆柱单位长度的电容为当时，令，则得第63页,共70页，2024年2月25日，星期天§5.4复变函数法1、复位函数法

2、保角变换法

第64页,共70页，2024年