工程电磁场数值分析5(有限元法).pptVIP

工程电磁场数值分析（有限元法） 华中科技大学电机与控制工程系 陈德智 2010.11 第4章 电磁场有限元法（Finite Element Method, FEM） 第4章 电磁场有限元法（FEM） 有限元的基本原理与实施步骤 有限元方程组的求解 前处理与后处理技术 渐近边界条件 矢量有限元法 求解运动导体涡流问题的迎风有限元法 1. 有限元法的基本原理与实施步骤 场域离散 以二维静电场泊松方程的求解为例。二维问题常使用三角形单元离散，便于处理复杂的场域形状，容易实现。 以下把单元e的贡献记为 单元分析：计算单元内积分对系数阵和右端项元素的贡献。 故 ， 写成一般形式，若一个三角形三个顶点编号为i, j, m（逆时针顺序），则 由于单元很小，做单元分析时通常可以取 f (e) 为常数值（可以认为等于三个顶点上的平均值）。因此 通过上述过程，对于一个“正常”的内部节点就建立起了一个代数方程。“非正常”的节点包括：媒质交界面衔接条件和场域边界条件，稍后再讨论。 上述以节点为序的分析过程对于有限元原理的说明是易于理解的。而在实际编程中，更有效率的是以单元为序，逐个计算单元系数阵[K(e)]，然后合成整体系数阵[K]。单元系数阵[K(e)]定义为 设i, j, k是节点的整体编号，