有限差分法计算金属槽内电位分布.doc

班级：物理08-2B姓名：胡艳学号：08070201010有限差分法计算金属槽内电位分布一、选题依据求解电位分布问题是物理学中最常见的问题之一，采用有限差分法解决此类问题是十分有效的。差分方程确定之后，一般选用迭代法求解，这是由于方程组系数矩阵中有大量元素为零，并且系数矩阵形成比较简单，有规律和重复。在迭代过程中常常可以一边形成系数矩阵一边计算，以节省内存，因而迭代法比直接法更常用。迭代法中又以超松弛迭代法最常用，以下给出超松弛迭代方法的公式。对于二维场泊松方程等距剖分差分格式公式为 （2.67）式中称为加速收敛因子或超松弛因子，它的数值决定超松弛程度，影响迭代解收敛的速度。取值范围是1 2 （2.68）加速收敛因子取值因问题而异，对于第一类边值问题，若一方形场域由正方形网格划分，每边的节点数为（n+1），则加速收敛因子可按下式计算： （2.69）若一矩形场域由边长为h的正方形网格划分为mh和nh，且m和n都很大（一般都要大于15），那么加速收敛因子为 （2.70）一般情况下，的最佳值只能是凭经验选取。对于其他形状的场域，也可用等效矩形面积的处理方法，即得出等效矩形面积后，再用式（2.70）求出最佳的。下面我们来看一个计算电位分布的实例。二、处理过程实例：有一长接地金属槽，横截面积如图2.15所示，其侧壁与底面电位均为零，顶盖电位为100V，求槽内电位分布。分析：对于此槽中间区段电位分布，可理想化为二维问题。选定直角坐标系，槽内电位函数满足拉普拉斯方程，构成一类边值问题。 按有限差分法计算步骤，解题过程如下。 选取计算区域，进行离散化处理。取槽内全部区域（或根据对称性取一半区域）进行计算。取正方形网格离散，步长h=0.1a，取x方向节点数为11，y方向节点数为7，每个节点坐标用双下标（i,j）表示。 给出采用超松弛迭代法的差分方程形式，即 场域为矩形，网格数m=10，n=6，按式（2.69）计算得到=1.14。 给出边界条件，本题为第一类边值条件，边界条件利用直接赋值方式，即 内部节点赋初值。由分析可知，场内各处电位值必然介于2和1之间。为加快迭代解收敛速度，采用等差递增方式，即 给出检查迭代解收敛的条件，即各节点相邻两次迭代值相对误差。 画出计算流程图如图2.16所示。 编写计算程序。采用FORTRAN语言编写的计算程序如下。变量与数组说明U（I，J）网格节点的电位值；I，J网格节点的坐标；N迭代次数；G不满足指定误差条件的离散点数；A第n次迭代时网格节点电位值；B第（n+1）次迭代时网格节点电位；C引入加速收敛因子后的网格节点电位；W相邻二次迭代值的相对误差值。源程序DIMENSION U (11,7)DATA U/77*0DO 5 I =2,10U (I,7)=1005 CONTINUEDO 10 I=2,10DO 15 J=2,6U(I,J)=100/6*(J-1)15 CONTINUE10 CONTINUEWRITE(*,11)WRITE(*,22) （U(II,KK),KK=1,11）,II=1,7）N=050 G=ON=N+1DO 20 I=2,10DO 25 J=2,6A=U（I,J）B=0.25*(U(I+1，J)+U(I,J+1)+U(I-1，J)+U(I,J-1)C=A+10456*(B-A)W=ABS（C-A）/C）IF(W.LT.1E-3) GOTO 100G=G+1100 U（I,J）=C25 CONTINUE20 CONTINUEIF(G.GT.0.)GOTO 50WRITE(*,*),N=1，NWRITE(*,*),ALPH=1,1.456WRITE(*,33)WRITE(*,22)(U(II,KK),KK=1,11),II=1,7)11 FORMAT (QHO,15X,28HINTIAL VALUE OF CALCULATION)22 FORMAT(1X,11F7.2)33 FORMAT(1HO,20X,21HRESULT OF CALCULATION)STOP END计算结果输出 INITIAL VALUE OF CALCULATION.00 100 .00 100 .00 100 .00 100 .00 100 .00 100 .00 100 .00 100 .00 100 .00 .00.00 83.33 83.33 83.33 83.33 83.33 83.33 83.33 83.33 83.33 .00.00 66.67 66.67 66.67 66.67 66.67 66.67 66.67 66.67 66.67 .00.00 50.00 50.00 50.00 50.00 50.00 50.00 50.00 50.00 50.00 .00.00 33.33 33.33 33.33 33.33 33.33 33.33 33.33 33.33 33.33 .00.00 16.67 16.67 16.67 16.67 16.67 16.67 16.67 16.67 16.67 .00.00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00N= 15ALPH=104560000 RESULT OF CALCULATION .00 100 .00 100 .00 100 .00 100 .00 100 .00 100 .00 100 .00 100 .00 100 .00 .00 .00 48.33 66.42 73.99 77.19 78.09 77.19 13.99 66.43 48.34 .00 .00 26.92 43.37 52.33 56.69 57.99 56.69 52.33 43.38 26.93 .00 .00 15.99 27.82 35.27 39.25 40.49 39.25 35.27 27.82 15.99 .00 .00 9.22 16.63 21.69 24.54 25.46 24.54 21.69 16.63 9.22 .00 .00 4.25 7.80 10.32 11.78 12.25 11.78 10.32 7.80 4.26 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00Stopprogam terminated.三、分析，总结值得注意的是，加速收敛因子的取值，的数值决定超松弛程度，影响跌代解收敛的速度，加速收敛因子的取值因问题而异，一般情况下，的最佳值只能凭经验选取，对于其他形状的场域，也可用等效矩形的面积的处理方法，即得出等效矩形面积后。再用式 求出最佳的。在解此物理问题时，把它理想化为二维问题，槽内电位函数满足拉普拉斯方程，构成了一类边值问题，采用超松弛迭代法的差分方程形式，用Fortran语言实现N此迭代。并可以求出第n此和（n+1）次迭代时网络节点的电位值，然后，引出加速收敛因子后