工程电磁场数值分析4(有限元法1D~2D)(免费阅读).pptVIP

工程电磁场数值分析（4）（电磁场有限元法） 华中科技大学电机与控制工程系 陈德智 Email: dzhchen@mail.hust.edu.cn Tel: Office: Room 108, 电机楼 2010.12 第4章 电磁场有限元法（Finite Element Method, FEM） 第4章 电磁场有限元法（FEM） 有限元基本原理与实施步骤：1D FEM 有限元基本原理与实施步骤：2D FEM 有限元方程组的求解 前处理与后处理技术 渐近边界条件 矢量有限元法 求解运动导体涡流问题的迎风有限元法 1. 有限元法基本原理与实施步骤：一维问题 （3）方程离散 （计算系数阵 [K] 和右端项 [b]） 基函数 Ni 只是一阶可导 的，不能严格满足微分方 程，称为“弱解”。 （3）方程离散 （3）方程离散 总体方程 （4）求解方程 一些补充说明：关于有限元的解 一些补充说明： 高斯数值积分 一些Gauss积分点和权值：（关于x=0对称，只给出一半） 一些补充说明： 线性单元与高阶单元 一些补充说明： 二阶单元 一些补充说明： 三阶单元 h方法和p方法的求解精度 By Jianming Jin. The Finite Element Method in Electromagnetics, 2nd Ed., 2002 作业：要独立完成，凡雷同者没分！！ 2. 有限元法基本原理与实施步骤：二维问题 以二维静电场泊松方程的求解为例。 场域离散 二维问题常使用三角形单元离散，便于处理复杂的场域形状，容易实现。 以下把单元e的贡献记为 单元分析：计算单元内积分对系数阵和右端项元素的贡献。 故 ， 写成一般形式，若一个三角形三个顶点编号为i, j, m（逆时针顺序），则 由于单元很小，做单元分析时通常可以取 f (e) 为常数值（可以认为等于三个顶点上的平均值）。因此 上述以节点为序的分析过程对于有限元原理的说明是易于理解的。而在实际编程中，更有效率的是以单元为序，逐个计算单元系数阵[K(e)]，然后合成整体系数阵[K]。单元系数阵[K(e)]定义为 设 i, j, m 是节点的整体编号，元素Kij在整体矩阵中的实际位置是第i行、j列；因此 必须合成到整体矩阵的第i行、j列元素上。 通过上述过程，对于一个“正常”的内部节点就建立起了一个代数方程。“非正常”的节点包括：媒质交界面衔接条件和场域边界条件。 对于静电场问题，媒质分界面衔接条件为 由于有限元方法能够自动满足媒质交界面条件，因此有限元法特别适合于处理多层复杂媒质问题。这是其它方法无可比拟的。 有限元方法的推导过程虽然看起来有些复杂，但是最终结果是非常简单而且优美的。因为边界条件的处理和媒质交界面条件的处理都非常方便，使得有限元方法在处理复杂媒质问题和复杂场域问题时得心应手，获得了广泛的应用，称为最重要的数值分析手段。有人用“功盖四方”来形容有限元，实不为过。 中国人在有限元的发明中有自己独特的贡献。 作业： （1）研究方向为数值计算的同学： 编写一个二维静电场有限元程序，计算右图所示问题，或其它自己找一个问题。 从而 再看边界部分： （1）在节点 i 的对边Gjm上，Ni＝0，故积分贡献为0； 结论：单元边界对积分的贡献为0。所以单元e对系数阵元素的贡献为： （2）在节点 i 的邻边Gij上，由于计算Kij时需要把具有公共邻边的单元的积分累加，此二单元的Ni是连续的；对于单一均匀媒质，要求相邻单元满足 ，故积分的贡献相互抵消。 右端项元素： 公式： 单元矩阵： 整体矩阵合成： 媒质交界面衔接条件 第一个条件是自动满足的（Why？），无须格外处理。 对于第二个条件，前面计算单元边界上积分 时，默认两边j 的法向导数相等，使内边界上的积分结果抵消。因此只要把泊松方程写成 或 满足的条件将是 , 从而也无需另行处理。 媒质交界面衔接条件 第一类边界条件（强加边界条件） 第一类边界节点是指边界上函数值 已知。因此处理方法是，合成整体系数阵之后，将该节点所在行的主元素置1，其它元素均置零，同时将右端项中对应元素设为已知函数值。 要保持对称性；有更简便的做法 第二类边界条件（自然边界条件） 第二类边界节点是指边