《用有限差分法解静电场边值问题.docx》由会员分享，可在线阅读，更多相关《用有限差分法解静电场边值问题.docx（7页珍藏版）》请在三一文库上搜索。
1、用有限差分法解静电场边值问题一、目的1. 掌握有限差分法的原理与计算步骤；2. 理解并掌握求解差分方程组的超松弛迭代法，分析加速收敛因子a的作用；3. 学会用有限差分法解简单的二维静电场边值问题，并编制计算程序。二、方法原理有限差分法是数值计算中应用得最早而又相当简单、直观的一种方法。应用有限差分法通常所采取的步骤是：采用一定的网格分割方式离散化场域。进行差分离散化处理。用离散的、只含有限个未知数的差分方程组，来近似代替场域内具有连续变量的偏微分方程以及边界上的边界条件(也包括场域内不同媒质分界面上的衔接条件)。 结合选定的代数方程组的解法，编制计算机程序，求解由上面所得对应于待求边值 问题的
2、差分方程组，所得解答即为该边值问题的数值解。现在，以静电场边值问题.：2 . -:2 ?一 +一 =0在 D 中(1)；:x2；:y2" = f(s)(2)为例，说明有限差分法的应用。f(s)为边界点s的点函数，二位场域D和边界L示于图5.1-1 中。y.x图5.1-1有限差分的网格分割1. 离散化场域应用有限差分法时，首先需从网格划分着手决定离散点的分布方式。通常采用完全有规律的方式，这样在每个离散点上可得出相同形式的差分方程，有效地提高解题速度。如图 5.1-1所示，现采用分别与 x, y轴平行的等距(步距为 h)网格线把场域 D分割成足够多的 正方形网格。各个正方形的顶点(也即
3、网格线的交点)称为网格的结点。这样，对于场域内 典型的内结点0,它与周围相邻的结点 1、2、3和4构成一个所谓对称的星形。2. 差分格式造好网格后，需把上述静电场边值问题中的拉普拉斯方程(1)式离散化。设结点 0上的电位值为。结点1、2、3和4上的电位值相应为 勺、%，%和甲4,则基于差分原理的应 用，拉普拉斯方程(1)式在结点0处可近似表达为少+饱+物-4甲1=0（3）这就是规则正方形网格内某点的电位所满足的拉普拉斯方程的差分格式，或差分方程。对于场域内的每一个结点，关系式（3）式都成立，都可以列出一个相同形式的差分方程。但是，对于近邻边界的结点，其边界不一定正好落在正方形网格的结点上，而可
4、能如图5.1-2所示。其中1、2为边界线上的结点，p、q为小于1的正数。仿上所述，可推得对这 些近邻边界结点的拉普拉斯方程的差分格式为12- 3411、°(4)+ + +（）平0 =0p（1 p） q（1 q） 1 p 1 q p q式中：CP和甲2分别是给定边界条件函数 f在对应边界点处的值，是已知的。qh0h °ph4图5.1-2近邻边界的结点3. 边界条件的近似处理为了求解给定的边值问题，还必须对边界条件，以及具体问题中可能存在的分界面上的衔接条件，进行差分离散化处理，以构成相应的差分边值问题。这里，我们只考虑正方形网 格分割下的边界条件的近似处理。第一类边界条件如果
5、网格结点正好落在边界L上，因此对应于边界条件（2）式的离散化处理，就是把点函数f （s）的值直接赋予对应的边界结点。如果边界L不通过网格分割时所引进的结点（例如图5.1-2中的1、2结点是边界线L与网格线的交点，并不是网格分割时所引进的网格结 点），那末在紧邻边界的结点的差分格式应选用（4）式，这时，把点函数 f的值直接赋予边界线L与网格线的交点1和2。第二类边界条件应当指出，从实际电场问题的分析出发，如图 5.1-3所示，以电力线为边界的第二类齐 次边界条件是常见的一种情况。=0cn L(5)这时，可沿着场域边界外侧安置一排虚设的网格结点，显然，对于边界结点0,由于该处竺=0，故必有 = Q
6、，因此相应于边界条件（5）式的差分计算格式为.n29+谿+气-4 %=0（6）边界sh图5.1-3第二类齐次边界的一种情况图5.1-4对称线上结点的差分格式同样，在许多工程问题中，常常能够判定待求电场具有某些对称性质，这样只需要计算某一对称部分的场就能完全决定整个场的分布。为此，还必须导出位于场的对称线上的结点所满足的差分计算格式。以对称线与网格结点相重合为例（见图5.1-4 ）,设AA'线为一对称线，对于位于对称线上的任一结点0,由拉普拉斯方程（因对称性，必有 ？=饱）可得相应的差分计算格式是2Q +平2+Q-4 %=0（7）媒质分界面上的衔接条件在此选取两种情况进行差分离散化的处理
7、。分界面与网格线相重合的情况；设分界面L与网格线相重合，如图5.1-5所示，在两种媒质£a和£b中电位都满足拉普拉斯方程。容易导得，两种媒质分界面上衔接条件在结点0的差分格式为其中；b二-1：2 冬1 K1 KO-3-4 -40 =：0(8)图5.1-5分界面与网格线相重合图5.1-6 分界面L对网格呈对角线形态分界面对于网格呈对角线形态的情况:如图5.1-6所示，分界面L对于网格呈对角线形态，在两种媒质 端和&b中电位 滞E满足拉普拉斯方程。容易导得，两种媒质分界面上衔接条件在结点0的差分格式为赤(12)莒 J "I 0 =0(9)其中K =。；b总之，
8、类似以上的分析处理方法，可以逐个导得各种类型的边界条件和衔接条件差分离 散化的计算格式。限于篇幅，在此不再展开。4. 差分方程组的求解在对场域D内各个结点（包括所有场域内点和有关的边界结点）逐一列出对应的差分方程，组成差分方程组后，就可选择一定的代数解法，以算出各离散结点上待求的电位值。注 意到差分方程组的系数一般是有规律的，且各个方程都很简单，包含的项数不多（最多不超过5项），因此，对于有限差分法，通常都采用逐次近似的迭代方法求解。在迭代法的应用中，为加速迭代解收敛速度，一般采用的是超松弛迭代法。由于编写计算机程序的需要，每一网格结点的位置由双下标（i, j）予以识别，如图 5.1-7所示。
9、对于差分方程（3）式，采用超松弛迭代法（规定迭代的运算顺序是：从左下角开始做起，即小的先做；对固定的如下迭代公式算得i, j小的先做。)，则关于结点0迭代到第(n+1)次时的近似值，应由矶汁=*(in)j)(n) :"(n) (n 1) :"(n 1) ( n)( (i 1,j). (i,j 1)- (i 1 j) . (i,j)- 4 .(i,j)4(10)3(i-1,j)0(i,j)4(i,j-1)2(i,j+1)1(i+1,j)图5.1-7 结点的双下标(i, j)标号式中：口称为加速收敛因子，其取值范围是1V 口2,当aA 2时，迭代过程将不收敛。加速收敛因子口有一
10、个最佳取值问题，但随具体问题而异。对于第一类边值问题，若一正方形场域由正方形网格分割(每边结点数为m+1)，则最佳收敛因子ct0可按下式计算2：-0 =TT1 sin-m在更一般的情况下，O0只能凭借经验取值。应当指出，为加速迭代解收敛速度，在迭代运算前，恰当地给定各内点的初始值(即所 谓第0次近似值)也是一个有效的途径。5.迭代解收敛程度的检验在超松弛迭代法的应用中，还必须涉及迭代解收敛程度的检验问题。对此，通常的处理 方法是：迭代一直进行到所有内结点上相邻两次迭代解的近似值满足修正条件(11)(12)时，终止迭代。即将式(12)作为检查迭代解收敛程度的依据。其中：W是指定的最大允许误差。6
11、.有限差分法的程序框图启动图5.1-8程序框图三、上机作业1 .设有一个长直接地金属矩形槽，（a=2b），如题5.1-1图所示，其侧壁与底面电位均为零，顶盖电位为 100V （相对值），求槽内的电位分布。Jhy中=100V1 bF=0中=中=0 (o aIB-题5.1-1图 矩形接地金属槽具体要求：编写一个计算机程序（用你熟悉的程序语言）求相邻两次迭代值的指定的最大允许误差小于10-5的迭代收敛解。 以步距h=旦的正方形网格离散化场域，然后应用有限差分法求电位中的数值解。40 根据场分布的对称性，试以半场域为计算对象，并以步距将该半场域由正方40形网格予以分割，然后用有限差分法求电位华的数值解
12、。 分别取r为n个不同的值和最佳值o（0,求电位甲的数值解，以此分析加速收敛因子的作用。从迭代收敛时的迭代次数和最终数值解这两方面总结自已的看法。用计算机描绘等位线分布。a b取中心点P（-,-）处电位的精确解（解析解）与数值解进行比较，说明误差范围。2.2 2在导电纸模拟的研究中，制备了如题5.1-2图所示的二维恒定电流场模型，其内两种导电媒质的电导率分别为和：2 （新=2技），它们在场域的对角线 L上接合。电极间外施电压10V。用有限差分法求此由两种导电媒质构成的二维恒定电流场内的电位分布。具体要求：写出该问题的数学模型，即恒定电流场的边值问题。按正方形网格离散化场域，步距 h=1cm的，然后用有限差分法求两种导电媒质 中电位的数值解。W=10-5 迭代收敛的指标：当网格内点相邻二次迭代近似值的绝对误差均小于 时，终止迭代循环。用计算机描绘等位线分布。