电磁场有限元法.ppt

电磁场有限元法 背景知识 有限元法是近似求解数理边值问题的一种数值技术； 1968年开始用于求解电磁场问题； 有限元法的本质是将微分方程的求解转化为代数方程的求解； 里兹有限元法、伽辽金有限元法 最大特点：以适当的形式将解域划分为有限个单元，在每个单元中构造子域基函数，利用里兹变分法或伽辽金法构造代数形式的有限元方程。 优点1：具有灵活的离散单元，可以精确地模拟各种复杂的几何结构，求解包含各种复杂形状、复杂媒质的电磁场问题。 优点2：所形成的方程组的系数矩阵为稀疏对称阵，利于求解。 缺点1：比积分方程法多一维，增加了未知量的数目。 缺点1：对于开放问题必须使用边界吸收条件截断计算空间， 增加了一定的计算复杂度。 在电磁场计算中，矢量基函数已基本取代了标量基函数； 一般情况下，分为频域有限元法和时域有限元法。 有限元的基本思路 将计算空间离散，划分为有限个小单元，小单元形式简单，数量有限； 根据小单元的不同形状，定义单元内的基函数，要求各基函数之间线性无关； 基函数是坐标的函数，每个基函数在单元内与各自特定的点或线相关。在这个特定的点或线上，定义在其上的基函数等于1，其它基函数等于0； 求解的目标就是单元内这些特定的点或线上的电场值。一旦已知，则单元内任一点的电场值都可以表示为单元内所有基函数的一个线性组合。 区域离散的概念 为了模拟复杂的区域形状，需要针对不同的问题采用不同的剖分单元形式，通常对于二维问题，我们采用三角形单元剖分；对于三维问题，采用的是四面体单元： 有限元边值问题 求解边值问题两种经典方 里兹(Ritz)变分方法 用变分表达式(也称为泛函)表示边值问题，泛函的极小值对应于给定边界条件下的控制微分方程。通过求泛函相对于其变量的极小值，可得到近似解。 伽辽金(Galerkin)方法 残数加权方法类型，它通过对微分方程的残数求加权方法来得到方程的解。 里兹(Ritz)变分方法 里兹(Ritz)变分方法 里兹(Ritz)变分方法 伽辽金(Galerkin)方法 伽辽金(Galerkin)方法 二维标量场有限元分析过程 空间离散 选择插值基函数 插值基函数 二维插值基函数的性质 用伽辽金法建立公式 组合成方程组 K矩阵的形成 列向量[b]的形成 列向量[g]的形成 求解方程组 三维有限元分析 空间离散 选择插值基函数 节点插值基函数 三维插值基函数的性质 用伽辽金法建立公式 组合成方程组 K矩阵的形成 通过局部坐标与全局坐标的对应关系，将在每个四面体单元中形成的局部[Ke]矩阵中的所有元素，依次填入全局的[K]矩阵中，最终完成K矩阵的形成。 K矩阵的行与列是节点全局编码的排序，矩阵中每个元素表示行号所对应的节点与列号所对应的节点，在基函数上的相关性。显然，一个节点所对应的基函数只与相邻节点对应的基函数相关。而且两节点基函数之间的相关性是相互的，因此最终形成的K矩阵是一个对称的稀疏阵。 列向量[b]的形成 列向量[g]的形成 时变电磁场中标量有限元的缺点 无论是一维二维还是三维的情况，可以看出，所有的插值基函数都是基于单元节点定义的标量函数。他们被通称为节点有限元。这种有限元常常会伴有非物理的或所谓伪解的出现，而且很难处理材料界面以及导体和介质边缘及角的场的奇异性。克服这种缺陷最好的办法是使用矢量的基函数直接表征待求的场。 齐次矢量波动方程 二维矢量有限元 二维棱边基函数 三角形单元中矢量基函数示意图 矢量基函数的性质 解析法求解S矩阵、T矩阵 高斯积分法求解S矩阵、T矩阵 空间任一四面体单元，到体坐标系下直角四面体单元的映射关系。 高斯积分法求解S矩阵、T矩阵 高斯积分法求解S矩阵、T矩阵 小节 至此，在每个单元内形成了两个6x6的系数矩阵； 矩阵的行表示单元内测试基函数的顺序； 矩阵的列表示单元内棱边基函数的顺序； 矩阵的每个元素表示各相应棱边基函数之间的自作用或互作用关系； 如果求得基函数的系数，则单元内任一点的电场均可通过插值求出。 从单元矩阵到全局矩阵的转换 由于： （1）矢量棱边基函数具有切向连续性； （2）Eie表示在第 i 条棱边上电场的切向分量； （3）电场在分界面上具有切向连续性； 所以，不同单元在相同棱边上的基函数的系数Eie必然相等。 全局矩阵的形成 通过局部编码与全局编码的对应关系，将在每个单元中形成的局部矩阵中的所有元素，依次填入全局矩阵中，最终完成全局矩阵的形成。 全局矩阵的行与列是未知量全局编码的排序，矩阵中每个元素表示行号所对应的未知量与列号所对应的未知量，在基函数上的相关性。； 显然，一个未知量所对应的基函数只与临近未知量所对应的基函数相关。而且两个基函数之间的相关性是相互的，因此最终形成的全局矩阵是一个nxn的对称稀疏阵， n表示未知量的个数，在低阶矢量有限元方