工程电磁场--7章--二维泊松方程的有限元法.pptVIP

7.3 二维泊松方程的有限元法 以二维静电场泊松方程的求解为例。 场域离散 二维问题常使用三角形单元离散，便于处理复杂的场域形状，容易实现。 以下把单元e的贡献记为 由于单元很小，做单元分析时通常可以取 f (e) 为常数值（可以认为等于三个顶点上的平均值）。因此 上述以节点为序的分析过程对于有限元原理的说明是易于理解的。而在实际编程中，更有效率的是以单元为序，逐个计算单元系数阵[K(e)]，然后合成整体系数阵[K]。单元系数阵[K(e)]定义为 设 i, j, m 是节点的整体编号，元素Kij在整体矩阵中的实际位置是第i行、j列；因此 必须合成到整体矩阵的第i行、j列元素上。 目标：依据加权余量法，利用分域基，建立离散的代数方程组，即确定系数{Kij} 和{bi}。 单元：互不重叠，覆盖全部场域；每个单元内介质是 单一、均匀的。 节点：网格的交点，待求变量的设置点。 该步骤需要记录的信息： 节点编号、节点坐标 节点属性(激励源、是否边界等) 单元编号 单元节点编号 单元介质 三角形单元内的基函数 设三角形三个顶点处待求函数值分别为u1, u2, u3。如果单元足够小，可以采用线性近似，将单元内任意p点的u(x,y)表示为 代入三个顶点的坐标和函数值，可以解出a、b、c。得到 单元节点的编号按逆时针