静态场边值问题

静态场边值问题第1页，共73页，2023年，2月20日，星期四4.1边值问题的分类第一类边值问题：已知位函数在全部边界面上的分布值

边值问题是指存在边界面的电磁问题。根据给定边界条件对边值问题分类：狄里赫利问题（Dirichlet）

第二类边值问题：已知位函数在全部边界面上的法向导数值第三类边值问题：已知一部分边界面上的位函数值,和另一部分边界面上位函数的法向导数值．诺埃曼问题（Neumann）混合边值问题第2页，共73页，2023年，2月20日，星期四边值问题框图一、二类边界条件的线性组合，即已知场域边界上各点电位的法向导数已知场域边界上各点电位值第一类边界条件第二类边界条件第三类边界条件边值问题参考点电位有限值场域边界条件分界面衔接条件自然边界条件微分方程边界条件第3页，共73页，2023年，2月20日，星期四解析法数值法实测法模拟法定性定量边值问题研究方法计算法实验法作图法有限差分法有限元法边界元法矩量法模拟电荷法积分法分离变量法镜像法、电轴法微分方程法保角变换法数学模拟法物理模拟法边值问题研究方法框图第4页，共73页，2023年，2月20日，星期四唯一性定理：在场域V的边界面S上给定位函数或的值，则泊松方程或拉普拉斯方程在场域V内的解唯一。唯一性定理的意义指出了静态场边值问题具有唯一解的条件；为静态场边值问题求解方法提供了理论根据，为结果

正确性提供了判据；唯一性定理是间接法求解拉普拉斯方程(泊松方程)的

理论根据。4.2唯一性定理（UniqunessTheorem）第5页，共73页，2023年，2月20日，星期四4.3直角坐标系中的分离变量法分离变量法是数理方程中应用最广泛的一种方法，它适用于求解具有理想边界条件的典型边值问题。分离变量法是通过偏微分方程求解边值问题。其基本思想是：首先要求给定边界面与坐标面相合，或分段相合；其次要求待求偏微分方程的解可表示为若干个函数的乘积，其中的每个函数分别仅是一个坐标变量的函数。这样，通过分离变量可将一个偏微分方程转化为多个常微分方程来求解。第6页，共73页，2023年，2月20日，星期四分离变量法解题的一般步骤：

根据边界的几何形状和场的分布特征选定坐标系，写出对应的边值问题（微分方程和边界条件）；

分离变量，将一个偏微分方程，分离成几个常微分方程；

解常微分方程，并叠加各特解得到通解；

利用给定的边界条件确定积分常数，最终得函数的解。在直角坐标系中，拉普拉斯方程为：直角坐标系中的分离变量法第7页，共73页，2023年，2月20日，星期四设可以表示为三个函数的乘积，即：当时

代入上式，得其中为分离常数，且※分析与

讨论上式中每项都只是一个变量的函数，其成立的唯一条件是三项中每项都是一个常数，故有第8页，共73页，2023年，2月20日，星期四当时

通解：当时

令其中为实数通解：或者同理可以求得和

利用给定的边界条件确定积分常数，最终得函数的解。双曲函数第9页，共73页，2023年，2月20日，星期四例4.3-1

横截面如图所示的导体长槽，上方有一块与槽相互绝缘的导体盖板，截面尺寸为a×b，槽体的电位为零，盖板的电位为U(x)，求此区域内的电位。在区域0<x<a、0<y<b内边界条件为：①x=0,φ(0,y)=0②x=a,φ(a,y)=0③y=0,φ(x,0)=0④y=b,φ(x,b)=U(x)解：选择直角坐标系第10页，共73页，2023年，2月20日，星期四通解：其中确定常数：当时所以由于在X方向上有重复零点（x＝0和a点），因此函数应为三角函数，即：且令分离变量第11页，共73页，2023年，2月20日，星期四当时故：分离变量法当时因为第12页，共73页，2023年，2月20日，星期四所以故：当时讨论两种情况和分离变量法当第13页，共73页，2023年，2月20日，星期四左右两边同乘以，并在区间(0，a)积分又有因此m=1,3,5,…分离变量法第14页，共73页，2023年，2月20日，星期四对应系数相等当因此分离变量法第15页，共73页，2023年，2月20日，星期四分离变量法的求解步骤

选择适当的坐标系,确定变量的个数；写出方程的通解；利用自然边界条件化简通解；利用电磁边界条件建立待定系数的方程并解方程,求出待定系数。第16页，共73页，2023年，2月20日，星期四4.4圆柱坐标系中的分离变量法圆柱坐标中的拉普拉斯方程（）为仅讨论二维平面场，即与坐标变量无关的情况令，代入上式得化简得

第17页，共73页，2023年，2月20日，星期四令第二项等于（）※分析与

讨论当时

当时

1讨论分离变量法由于是周期性函数，即所以，且为正整数，即。第18页，共73页，2023年，2月20日，星期四当时

当时

2讨论欧拉方程令代入上式，整理得：通解第19页，共73页，2023年，2月20日，星期四例4.4-1一根半径为a、介电常数为的无限长介质圆柱体置于均匀外电场中，且与相垂直。设外电场方向为x轴方向，圆柱轴与z轴重合（如图所示），求圆柱内、外的电位函数。解：选择圆柱坐标系。设园柱内、外的电位分别为、，并假设零电位点在坐标原点。显然、均是与z无关的二维场，都满足拉氏方程：分离变量法第20页，共73页，2023年，2月20日，星期四则，圆柱坐标系中二维场的通解为：令由题意知，场分布对称于x轴，即：故通解中只能包含余弦项即：当时，其中：第21页，共73页，2023年，2月20日，星期四当时，应为有限值，即中不能有项。当时，由介质分界面条件知。切向：第22页，共73页，2023年，2月20日，星期四法向：联立方程（1）、（2）解得：第23页，共73页，2023年，2月20日，星期四介质圆柱体内、外的电场强度为：第24页，共73页，2023年，2月20日，星期四边界条件：1例4.4-1将半径为的无限长导体圆柱置于真空中的均匀电场中，柱轴与垂直，求任意点的电位。2解：第25页，共73页，2023年，2月20日，星期四12分离变量法第26页，共73页，2023年，2月20日，星期四分离变量法第27页，共73页，2023年，2月20日，星期四4.5球坐标系中的分离变量法球坐标中的拉普拉斯方程（）为仅讨论场问题与坐标无关时的情形令，代入上式并整理得令两项分别等于常数和第28页，共73页，2023年，2月20日，星期四引入一个新的自变量则有－－勒让德方程当时，勒让德方程有一个有界解

－－勒让德多项式1讨论第29页，共73页，2023年，2月20日，星期四下面是前几个勒让德多项式当时，；时，当－－勒让德多项式勒让德多项式图形

勒让德多项式具有正交性第30页，共73页，2023年，2月20日，星期四2讨论欧拉方程

通解综上，球坐标中的通解为令第31页，共73页，2023年，2月20日，星期四例4.5-1在均匀电场中，放置一个半径为a的导体球，球心在原点。求球外的电位及电场强度分布。解：边界条件①

通解选择球坐标系，并设导体球为零电位。第32页，共73页，2023年，2月20日，星期四②

第33页，共73页，2023年，2月20日，星期四求解位于接地导体板附近的点电荷产生的电位非均匀感应电荷等效电荷4.6镜像法接地导体球附近有一个点电荷产生的电位等效电荷非均匀感应电荷产生的电位很难求解，可以用等效电荷的电位替代。第34页，共73页，2023年，2月20日，星期四

镜像法的目的：把具有典型边界的静态场的计算问题，转化为无限大均匀媒质空间中的问题求解，以达到简化计算的目的。

镜像法基本思路：在待求解场域外的适当位置，以虚拟电荷替代分界面上导体的感应电荷或媒质的极化电荷。镜像法

镜像法的理论依据：由唯一性定理知：满足同一方程和同样边界条件的位函数的解是相同的。所以引入镜像电荷后，应该：保持原方程和原边界条件不变。

镜像电荷位置选择原则

1镜像电荷必须位于待求解场域以外2镜像电荷的引入不能改变原问题的边界条件第35页，共73页，2023年，2月20日，星期四1.平面边界的镜像法例4.6-1求置于无限大接地平面导体上方，距导体面为处的点电荷的电位，并求出场的分布。一、静电场中的镜像法（除q所在点外的区域）（边界条件）解：当时，

；当时，当时，；确定镜像电荷（位置、数量和大小）上半空间内任意点的电位为平面导体的镜像Pzz第36页，共73页，2023年，2月20日，星期四平面导体上的感应电荷密度为验证边界条件当Z=0时,所以平面导体的镜像z第37页，共73页，2023年，2月20日，星期四对应的场分布图：镜像法第38页，共73页，2023年，2月20日，星期四例4.6-2为无限大接地的导电（）平面（电壁），在出处有一无限长均匀带电的细直导线，导线与轴平行且经过直角坐标（）点，求上半空间（）的电位函数。设细直导线的电荷密度为，则镜像线电荷密度为。解：电壁的作用可以等效为：镜像位置处的镜像线电荷带电体系在空间的电位为zz第39页，共73页，2023年，2月20日，星期四式中不能选为无穷远点，同样式中所以第40页，共73页，2023年，2月20日，星期四例4.6-3

设介电常数分别为和的两种介质，各均匀充满半无限大空间，两者的分界面为平面，在介质1中有一点电荷，距分界面的距离为，如图所示.试求整个空间中任一点的电位函数.设两个区域的电位函数为()和()解：（除q点外的上半空间）（下半空间）当时

当时

当时

第41页，共73页，2023年，2月20日，星期四在介质分界面处，电位函数满足边界条件

解得第42页，共73页，2023年，2月20日，星期四空间的电位分布为：点电荷位于不同介质平面上方的场图第43页，共73页，2023年，2月20日，星期四2.角形区域的镜像法所有相互成角的两块半无限大接地导体平面间的场四都可用镜像法来求解，其像电荷个数为P第44页，共73页，2023年，2月20日，星期四例4.6-4设在点电荷附近有一接地导体球，求导体球外空间的电位及电场分布。边值问题：点电荷对接地导体球面的镜像3.球面边界的镜像法设镜像电荷位于球内，球面上任一点电位为（除q点外的导体球外空间)第45页，共73页，2023年，2月20日，星期四镜像法解得第46页，共73页，2023年，2月20日，星期四由叠加原理，接地导体球外任一点P的电位与电场分别为点电荷位于接地导体球附近的场图镜像电荷不能放在当前求解的场域内导体上总的感应电荷等于镜像电荷接地导体球外的电场计算第47页，共73页，2023年，2月20日，星期四在接地球的基础上判断镜像电荷的个数、大小与位置边值问题：例4.6-5计算不接地金属球附近放置一点电荷时的电场分布。点电荷对不接地金属球的镜像感应电荷分布及球对称性，在球内有两啊个等效电荷

正负镜像电荷绝对值相等

正镜像电荷只能位于球心（除q点外的导体球外空间）放置镜像电荷：+q‘和-q’第48页，共73页，2023年，2月20日，星期四点电荷位于不接地导体球附近的场图任一点电位及电场强度为第49页，共73页，2023年，2月20日，星期四4.柱面边界的镜像法例4.6-6线电荷密度为的无限长带电直线与半径为a的接地无限长导体圆柱的轴线平行，直线到圆柱轴线的距离为，如图所示。求圆柱外空间的电位函数。解：第50页，共73页，2023年，2月20日，星期四两端对求导可得圆柱面上电位为零

第51页，共73页，2023年，2月20日，星期四因为

所以圆柱面上的感应电荷密度为圆柱面上单位长度的感应面电荷为第52页，共73页，2023年，2月20日，星期四例4.6-7设分界面为平面的两个半无限大空间中，分别充满磁导率为和的两种均匀介质，在介质1中存在一平行于分界面的长直线电流I，与分界面的距离为，试求空间的磁场。二、静磁场中的镜像法第53页，共73页，2023年，2月20日，星期四下半空间（）用镜像电流来代替分界面上的磁化电流。上半空间（）用一镜像电流代替分界面上的磁四四四化电流。解：在分界面上第54页，共73页，2023年，2月20日，星期四从而得到相应的磁场为第55页，共73页，2023年，2月20日，星期四镜像法小结

镜像法的理论基础是静态场的唯一性定理；镜像法的实质是用虚设的镜像电荷替代未知电荷的分布，使计算场域变为无限大均匀介质空间；镜像法的关键是确定镜像电荷的个数、大小及位置；应用镜像法解题时，注意：镜像电荷只能放在待求场域以外的区域。叠加时，要注意场的适用区域。第56页，共73页，2023年，2月20日，星期四4.7有限差分法1.二维泊松方程的差分格式通常将场域分成足够小的正方形网格,网格线之间的距离为h,节点0,1,2,3,4上的电位分别用和和表示。（1）（2）二维静电场边值问题：有限差分的网格分割有限差分法（FiniteDifferentialMethod）是基于差分原理的一种数值计算法。其基本思想：将场域离散为许多小网格，应用差分原理，将求解连续函数的泊松方程的问题转换为求解网格节点上的差分方程组的问题。第57页，共73页，2023年，2月20日，星期四将和分别代入式（3），得（4）（5）由（4）–（5）（6）由（4）+（5）（7）设函数在处可微，则沿方向在处的泰勒公式展开为（3）有限差分法第58页，共73页，2023年，2月20日，星期四（9）同理（8）将式（7）和（9）代入式（1），得到泊松方程的五点差分格式即当场域中，得到拉普拉斯方程的五点差分格式有限差分法第59页，共73页，2023年，2月20日，星期四若场域离散为矩形网格，差分格式为：122.边界条件的离散化处理第二类边界条件边界线与网格线相重合的差分格式第一类边界条件给边界离散节点直接赋已知电位值。介质分界面衔接条件的差分格式边界条件的离散化处理其中第60页，共73页，2023年，2月20日，星期四3.差分方程组的求解方法高斯赛德尔迭代法式中

迭代顺序可按先行后列，或先列后行进行。

迭代过程遇到边界节点时，代入边界值或边界差分格式，的直到所有节点电位满足为止。高斯——赛德尔迭代法有限差分法第61页，共73页，2023年，2月20日，星期四松弛迭代法式中——加速收敛因子最佳收敛因子的经验公式：（正方形场域、正方形网格）（矩形场域、正方形网格）欠松弛迭代超松弛迭代迭代发散高斯赛德尔迭代法第62页，共73页，2023年，2月20日，星期四

迭代收敛的速度与电位初始值的给定及网格剖分精细有关

迭代收敛的速度与工程精度要求有程序框图如下：赋边界节点已知电位值赋予场域内各节点电位初始值按超松弛法进行一次迭代，求打印结束N所有内点相邻二次迭代值的最大误差是否小于Y累计迭代次数N=0N=N+1启动迭代解程序框图第63页，共73页，2023年，2月20日，星期四编程题1

横截面如图所示的导体长槽，上方有一块与槽相互绝缘的导体盖板，截面尺寸为a×b，槽体的电位为零，盖板的电位为100V。用高斯赛德尔迭代法编程求解。要求：步长h＝1，x、y方向的网格数为m＝16，n＝10，迭代精度为。计算：迭代次数N与分布。有限差分法第64页，共73页，2023年，2月20日，星期四编程题2

横截面如图所示的导体长槽，上方有一块与槽相互绝缘的导体盖板，截面尺寸为a×a，槽体的电位为零，盖板的电位为100V。用超松弛迭代法编程求解要求：步长h＝1，x、y方向的网格数为m＝10，n＝10，迭代精度为。计算：迭代次数N与分布。有限差分法第65页，共73页，2023年，2月20日，星期四作业：P1274-4、4-8、4-13、4-22、4-29第66页，共73页，2023年，2月20日，星期四双曲正弦双曲余弦a)：其定义域为:(-∞,+∞)；

b)：是奇函数；

c)：在定义域内是单调增a)：其定义域为:(-∞,+∞)；

b)：是偶函数；

c)：其图像过点(0,1)；反双曲正弦函数