电磁场数值计算作业报告.docxVIP

《电磁场数值计算》—有限元法报告第一题（一）问题描述及数学物理模型的建立①①9 ④ 5 ⑤ 6 （12） 128 ③⑥（11） 11② 3 ⑦ 4 ⑩7 ① 1 ⑧ 2 ⑨ 10有一矩形场区，尺寸为69，如图1所示，在区域内部J=0，的左边界A=0,右边界A=100，上下边界满足：，媒质均为的磁导率为，利用有限元法求A的分布首先，此问题的微分方程及其定解条件为： 图1 求解场域及网格划分 如果利用有限元来求解，那么对应的泛函为（注：分别为上边界和下边界）： 第一类边界条件为：（二）场域的离散化及单元节点的编号处理 网格的划分如图1所示，先求出每个单元的刚度矩阵，然后再进行整体合成，得到系数矩阵，具体的做法可参考教材，注意，这里左边界和右边界属于第一类边界，因此，利用下式进行计算： 这里，因为J=0,所以=0， 节点编号如图1所示，这里，由于网格的划分和节点的编号都不是很规则，但节点数目和网格数目不是很多，所以，在程序中，通过穷举的办法，找到相应单元的（i，j，m）的值，以及相应的节点坐标。 各单元的（i,j,m）（按逆时针方向）如下： 1(3,7,1) 2(3,8,7) 3(9,8,3)