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

电磁场数值计算上机题报告 第一题计算长直接地金属槽中的电场分布。金属槽横截面如图1所示，其侧壁与底面电位均为零，顶盖电位相对值为10。槽内电位函数满足拉普拉斯方程。计算槽内电位分布。要求：（1）先用正方形网格粗分，每边取4个网格计算，取不同的松弛因子，比较其收敛速度。取计算精度为千分之一。（2）划分网格加倍，计算电位分布，并与上面计算结果比较。（一）建立问题的数学物理模型首先列出方程及其边界条件槽内的电位满足二维的拉普拉斯方程：其中的范围是：边界条件是：图1（二）算法设计及其实现在本题中，因为区域为正方形区域，网格采用正方形网格，每边四个网格，因此，每边要有5个结点，网格数m=n=4，比较少，不能用课本中的公式计算收敛因子，取收敛因子为，计算的程序的流程为：①选取计算的场域，并划分网格，网格划分如右图2所示：右图表示网格的划分，共16个网格，一共划分了25个结点，每个节点用相应的下标（i,j）来表示，对应的电位为。根据题意，边界条件的处理如下：图 2其中n为一行对应的节点数，根据题意，这里n=5 ②用表示节点的电位，设经过第n次迭代之后的结点电位用来表示，则超松弛迭代法的差分格式（记）为： 边界条件前面已经给出③给各个节点赋初值，对于非边界上的点（对于边界点的值前面已经赋过初值），如下④迭代计算，直到已经满足精度条件为止，这里精度为0.0001,最后输出计算的结果，输出结