电磁场数值计算上机题报告 .pdfVIP

电磁场数值计算上机题报告 第一题 计算长直接地金属槽中的电场分布。金属槽横截面如图1所示，其侧壁与底 面电位均为零，顶盖电位相对值为 10。槽内电位函数满足拉普拉斯方程。计算 槽内电位分布。 要求：（1）先用正方形网格粗分，每边取4个网格计算，取不同的松弛因子， 比较其收敛速度。取计算精度为千分之一。（2）划分网格加倍，计算电位分布， 并与上面计算结果比较。 （一）建立问题的数学物理模型 首先列出方程及其边界条件 槽内的电位满足二维的拉普拉 斯方程： 2 2 2    0 2 2 x y 其中 的范围是： x,y 0xa,0ya 边界条件是：   0 x0 xa  0; 10 y0 ya 图1 （二）算法设计及其实现 在本题中，因为区域为正方形区域，网格采用正方形网格，每边四个网格，因 此，每边要有5个结点， 网格数m=n=4，比较少，不能用课本中的 公式计算收敛因子，取收敛因子为 1.4， 计算的程序的流程为： ①选取计算的场域，并划分网格，网格划 分如右图2所示： 右图表示网格的划分，共 16 个网格，一 共划分了25 个结点，每个节点用相应的下 标（i,j）来表示，对应的 电位为u(i,j) 根据题意，边界条件的 。 处理如下： 图 2 u(1,j ) u (i ,1)u n(,1)  u(k, j) 10(2k n 1) 其中n 为一行对应的节点数，根据题意，这里n=5 u(n) ② 用u(i,j)表示节点的电位，设经过第n 次迭代之后的结点电位用 来表 i,j 示，则超松弛迭代法的差分格式（记u u(i,j)）为： i,j  (n 1 ) n( ) n ( ) n ( )  n ( 1 ) n ( 1 ) n ( ) u  u  ( u  u  u  u 4 u) i, j ,i j 4  1i , j , i j1  1i , j  , i j1 , i j 边界条件前面已经给出 ③给各个节点赋初值，对于非边界上的点 （对于边界点的值前面已经赋过初 值），如下 10 u(i,j) (j 1) n1 ④迭代计算，直到已经满足精度条件为止，这里精度为0.0001,最后输出计算 的结果，输出结果见生成的数据文件shuju.txt 中 ⑤计算框图如图3所示 图 3 计算框图 ⑥用fortran90语言编写计算机程序，计算各点的电位，程序清单如下： program homework1 implicit none integer i,j,k,n real,a ocatable :: u(:,:) !定义动态数组 real a,b,c,e