工程电磁场数值计算FEM程序（EMF2D）.ppt

 LASTND 节点总数 N0 LASTEL 三角元总数 E0  X, Y（N0） 节点 N 的 x、y 坐标 NODEI, NODEJ, NODEM（E0） 三角元 e 的三顶点编号 SD 激励源密度 p  CPM 媒质特征参数 ? ● 变量与哑元说明 前处理：三角元剖分数据文件 U（N0） 节点位函数值，开始时还表示待解线性代数方程的右端项 给定问题的物理参数 LBN 强制边界条件的节点总数 BNODE（LBN） 强制边界节点编号 BPTAL（LBN） 强制边界节上给定的位函数值 NX1, NY1 取决于特定问题，为构造强制 边界条件而引入的工作单元 待求的数值解 U（x、y ） L（N0） 系数矩阵主对角线元素 Kii 的地址向量 SL 系数矩阵下三角变带宽一维压缩存贮的总长度 K（SL） 系数矩阵下三角变带宽一维压缩存贮数组 强制边界条件对应的 变量与数组 NUMBER1 源区单元从 1 起始编号所得最后的编号 构造系数矩阵 源区单元总数  PARAMETER （LND=700, LEL =1200, LNZ =30100） PARAMETER （LB=100） REAL K，K0 INTEGER E， SL， BNODE(LB) COMMON /LKK/ L(LND)，K0(LND)，K(LNZ) /NM/ NN1，N1，N COMMON /XY/ LASTND，X(LND)，Y(LND)* /ELEM/ LASTEL，NUMBER1，NODEI(LEL)，NODEJ(LEL)，NODEM(LEL) COMMON /VA/ U(LND) DIMENSION BPTAL(LB) OPEN (6，FILE=’EMFOUT.DAT’，FORM= ’FORMATTED’) OPEN (5，FILE= ’AGEN.DAT’，FORM= ’FORMATTED’STATUS= ’OLD’) ● 程序（EMF2D） ( 1－9 行) 说明语句 ( 10－11 行) 打开数据文件 C IN