第4章静态场边值问题的解法

第四章静态场边值问题的解法,4.1分离变量法4.1.1直角坐标系中的分离变量法4.1.2圆域内的二维场问题4.1.3球坐标中的分离变量法4.2有限差分法,静电场和恒定电场的分析归结为求解相应的泊松方程或拉普拉斯方程。给定边值的泊松方程和拉普拉斯方程有唯一点解。三类给定边值：,当媒质不均匀时，作为定解条件还需加入辅助边界条件,如果场域扩展为无界区域，还需提出无限远处的边界条件。微分方程与边界条件一起构成边值问题。,4.1.1直角坐标系中的分离变量法定解问题一长直接地金属槽截面如图。其侧壁与底面的电位均为零，而顶盖电位40。求槽内电位分布。,解：设金属槽的长度远大于截面尺度，忽略边缘效应，将问题简化为二维场来分析。,4.1分离变量法,问题转化为如下定解问题,令(x,y)=X(x)Y(y)，代入方程，得,只有当左右两边都是常数时，上式才成立，令此常数为。将偏微分方程转化为两个常微分方程。,边界条件转化为：,1需先解下列边值问题。,是待定常数，要解出使方程有非零解的值和此非零解X(x)。该边值问题称为常微分方程在此边值条件下的固有值（特征值）问题。称为该问题的固有值（特征值），X(x)称为该问题的固有（特征）函数。,分三种情况讨论。,（1）设0，令2，常微分方程的通解为,代入边界条件A0，Bsina0，B不能为零，否则只有零解。sina0，,所以固有值和固有函数分别为,2再解Y（y）,3满足方程和部分边界条件的一组特解是：,这组解满足方程和部分边界条件，但不一定满足所有边界条件。,由于方程是线性的，可应用叠加原理，将所有特解叠加,只要无穷级数收敛，且能关于x和y逐项微分两次，则（x,y）与n(x,y)一样满足方程和部分边界条件。适当选择Dn和En，可使（x,y）满足方程和所有边界条件。,在上两式两边分别乘并作积分。,无穷级数的系数应满足,无穷级数解满足原定解条件。通常取4项就能得到足够精确的结果。,分离变量法的主要步骤：（1）分离变量。将偏微分方程的定解问题化为常微分方程的定解问题（线性齐次偏微分方程）。（2）确定固有值和固有函数。当边界条件是齐次的时，利用其求固有值，并求出满足零边界条件的非零解。（3）求解其他常微分方程。得到满足齐次边界条件的偏微分方程的特解Un（x，y）。（4）将所有Un（x，y）叠加，利用其中的常数使其满足偏微分方程其余的定解条件。,教材例4.1.1求如图长方体积中的电位函数。边界条件为除zc面电位不为零外，其他各表面的电位都为零。Zc表面上给定的电位函数为U（x，y）。,解：,分离变量，令(x,y,z)=X(x)Y(y)Z(z),有二个独立的本征值。边界条件可分解为：,X(0)=X(a)=0Y(0)=Y(b)=0,利用齐次边界条件求出本征值和本征函数。,与上述方法一样，可求出,其中m1，2，。n=1,2,。,满足部分齐次边界条件的偏微分方程的一组特解为,为使解满足所有的边界条件，将所有特解叠加。,从已知边界条件,根据线性齐次方程的叠加原理，通过调整系数可使(x,y,z)满足所有边界条件。,将U（x，y）展开成双重傅立叶级数,比较系数可得,原定解问题的解是,4.1.2圆域内的二维场问题设在一圆形场域D内电位函数满足拉普拉斯方程，场域边界L上给定为第一类边值，求此二维场中位函数的分布。解：选取柱坐标，待求的边值问题为,区域内,边界L上,满足自然周期条件：（r,）=（r,+2）,令（r，）R(r)()代入方程,偏微分方程边值问题转化为常微分方程边值问题。,固有值问题,1先解固有值问题,（1）当0，令2，,a和b中至少有一不为零。,若a不为零，,若b不为零，,都有22n，,即：,固有值固有函数,2再解R(r),欧拉方程，可令r=et，代入方程得,当0，,当n2，n1，2，.,3得到一组满足方程和自然周期条件的特解。,其中,满足方程的一般解可表示为,再利用边界条件=f(s)及其他自然边界条件，即可定出其中的各常数。,选定坐标系。通常给定的边界面应与相应当坐标面相重合，以得到问题的最简洁定解条件。,例4.2.1一横截面半径为r0，介电常数为1的长直介质圆柱体放置在均匀的外电场中（场强的数值为E0，它的方向与介质圆柱的轴线垂直），均匀场中介质的介电常数为2，求圆柱体放入后，场中的电位和电场分布。,解：当长圆柱轴向长度远大于横截面半径时，其中间段场的分析可简化为二维场问题。采用柱坐标，柱内和柱外电位函数1和2分别满足泊松方程。,由于无限远电位不是零，可令原点电位为零。10，(r=0),外电场可用电位函数表示：0E0 x=-E0rcos在无限远处介质圆柱体产生的影响应当消失。2E0rcos（r）边界处,利用上述解题方法，得柱内外电位的通解分别是,利用边界条件可定出各系数。,以下利用此题的具体特征解题。（1）场的分布对称于x轴，,(r,)=(r,-)，,(2)因场的对称性，y轴是等位线,(r,/2)=0,n=1,3,5.为奇数，代入,可得满足,时，,n只能取1。否则，可导致另外与x轴斜交的零等位线。,所以：,圆柱内外电位1(r,)和2(r,)的一般形式为,10，(r=0),得b10。,2E0rcos（r）,得a2E0。,因而有,从,可求出,圆柱内场强均匀分布，且与外场方向一致。,圆柱内电场强度,4.1.3球坐标中的分离变量法,如果求解球空间或场有球面边界，应用球坐标更方便。球坐标中的拉普拉斯方程为,令n(n+1),（任一实数总可写成这种形式）。,即在球坐标下，电位的解为,其中,Pn(x)可用罗德利克公式表示,勒让德多项式具有正交性。,在一定的条件下函数f(x)可按固有函数Pn(x)展开。,教材例4.3.1在均匀外场E0中放置一半径为a的介质球，求的介电常数为，求外为空气，计算球内外的电位函数。解：采用球坐标，电位与无关。设球外电位为1，球内电位为2。外电场E0引起的电位为0-E0rcos。r趋于无穷时，10-E0rcos,r=0时，2有界。,在整个区域内，电位满足方程,用分离变量法，可解出电位（r,），只有n为整数时，方程在区间内才有有界的解。,r=0时，2有界,Rn（r）有界，,r趋于无穷时，10-E0rcos,右边可看成将cos展开成勒让德函数的级数。比较两边系数，左边只有n1的项不为零。,，且A1E0。,Cn中只有C1不为零。,可解出,例4.2.2已知半径为r1的球形域内的电位u满足拉普拉斯方程，求u的分布，其中在边界球面上,解：采用球坐标。从方程和边界条件分析，u与无关，定解问题是：,分离变量,只有当n整数，方程在内才有有界的解,欧拉方程的解,要使u有界，Rn有界，C20。,用叠加原理,利用边界条件,比较系数：C0=1/3,C2=2/3,Cn=0(n2),原定解问题的解：,当场域边界的几何形状比较复杂时，很难用解析法进行分析。应采用数值计算法。有限差分法将连续场域内的问题转化为离散系统的问题，通过离散化模型上各离散点的数值解来逼近连续场域内的真实解。1差分原理设有一函数f(x)，当独立变量x有一微小增量xh，相应f(x)的增量为：f(x)=f(x+h)-f(x),称为函数f(x)的差分。不同于增量为无限小的微分，差分被称为有限差分。当h很小时，f(x)df(x).,4.2有限差分法,中心差分f(x)=f(x+h/2)-f(x-h/2),一阶差商：,二阶差商,偏导数也可用差商近似表示。因而偏微分方程可表示为差分方程（代数方程）。,2以二维场为例，将边值问题转化为一组差分方程组（代数方程组）。设边值问题是,(1)决定离散点的分布方式。按正方网格划分，网格边长（步长）h，网格线的交点称结点。设结点O上的电位为(xo,yo)=o,结点1，2，3，4上的电位为1，2，3，4。,任一点x的电位,考虑1，3两点x1=xo+h,x3=xo-h,边界条件也可进行离散化处理，对第一类边值，可直接把点函数f(s)的值赋予各边界结点。,3差分方程的解法设将场域划分如图.边界上的值分别为f1,f16。在各内点上作出差分，泊松方程变成下列差分方程组,解出关于1，2.9的代数联立方程组，即可求出各点的函数值。,算法简单迭代法，以解拉普拉斯方程为例。（1）设定内点初值，用计算机解题时，可都取零值。,（2）按一固定顺序(从左到右，从下到上)依次利用,计算内点o点的新值。即o点的新值就是围绕该点的4个点的电位的平均值。,如（j，k）点在第n1次迭代时按下式计算：,当所有的内点都计算完后，用他们的新值代替旧值，完成一次迭代。再进行下一次迭代。直到每一点计算得到的新值与旧值之差小于指定的范围。这种方法的特点是用前一次迭代的得到的结点电位作为下一次迭代时的初值。,超松弛法简单迭代法收敛慢。超松弛法的改进：,（1）,即计算（j，k）点时，左边点（j-1,k）和下面点（j，k-1）用的是新值。这种迭代方法称为高斯赛德尔迭代法。,（2）将上式写成增量的形式，,引进加速收敛因子，在12之间。,加速解的收敛。2时，迭代过程将发散。,最佳收敛因子0的取值随问题而异。对第一类边值问题，正方形场域，网格按正方形划分，每边结点数p1，则,一长直接地金属槽截面如图。其侧壁与底面的电位均为零，而顶盖电位4100。求槽内电位分布。,例4.2.3,解：二维场第一类边值问题。将二维场域划分成正方形网格，步距ha/4。场域内任一点电位应满足二维拉普拉斯方程的差分计算格