电磁场有限元法.pptxVIP

电磁场有限元法第1页/共90页 2背景知识有限元法是近似求解数理边值问题的一种数值技术；1968年开始用于求解电磁场问题；有限元法的本质是将微分方程的求解转化为代数方程的求解； 里兹有限元法、伽辽金有限元法最大特点：以适当的形式将解域划分为有限个单元，在每个单元中构造子域基函数，利用里兹变分法或伽辽金法构造代数形式的有限元方程。优点1：具有灵活的离散单元，可以精确地模拟各种复杂的几何结构，求解包含各种复杂形状、复杂媒质的电磁场问题。优点2：所形成的方程组的系数矩阵为稀疏对称阵，利于求解。缺点1：比积分方程法多一维，增加了未知量的数目。缺点1：对于开放问题必须使用边界吸收条件截断计算空间， 增加了一定的计算复杂度。在电磁场计算中，矢量基函数已基本取代了标量基函数；一般情况下，分为频域有限元法和时域有限元法。第2页/共90页 3有限元的基本思路将计算空间离散，划分为有限个小单元，小单元形式简单，数量有限；根据小单元的不同形状，定义单元内的基函数，要求各基函数之间线性无关；基函数是坐标的函数，每个基函数在单元内与各自特定的点或线相关。在这个特定的点或线上，定义在其上的基函数等于1，其它基函数等于0；求解的目标就是单元内这些特定的点或线上的电场值。一旦已知，则单元内任一点的电场值都可以表示为单元内所有基函数的一个线性组合。第3页/共90页 4区域离散的概念为了模拟复杂的区域形状，需要针对不同的问题采用不同的剖分单元形式，通常对于二维问题，我们采用三角形单元剖分；对于三维问题，采用的是四面体单元：第4页/共90页 5有限元边值问题典型的边值问题可用区域内的控制微分方程和包围区域边界上的边界条件来定义：LФ=f 其中L为微分算符，f为激励或者强加函数，Ф是未知量。在电磁学中，控制微分方程包括简单的泊松方程以及复杂的标量波动方程，甚至也有更复杂的矢量波动方程。边界条件有简单的狄利克雷(Dirichlet)条件和诺曼(Neumann)条件，也有复杂的阻抗和辐射边界条件，甚至还有更复杂的高阶条件。第5页/共90页 6求解边值问题两种经典方里兹(Ritz)变分方法 用变分表达式(也称为泛函)表示边值问题，泛函的极小值对应于给定边界条件下的控制微分方程。通过求泛函相对于其变量的极小值，可得到近似解。伽辽金(Galerkin)方法 残数加权方法类型，它通过对微分方程的残数求加权方法来得到方程的解。第6页/共90页 7里兹(Ritz)变分方法LФ=f 的解等于下式泛函对 的解泛函：vj是定义在全域上的展开函数cj是待定的展开系数第7页/共90页 8里兹(Ritz)变分方法将试探函数代入泛函：令其对ci的偏导数为零，从而得到线性代数方程组第8页/共90页 9里兹(Ritz)变分方法其中：（应用算符L的自伴性质）求解该方程组可以得到 LФ=f 的近似解第9页/共90页 10伽辽金(Galerkin)方法假设 是 LФ=f 的近似解，则得到非零的残数为：——使用残数加权法求解微分方程残数加权方法要求wi是所选择的加权函数第10页/共90页 11伽辽金(Galerkin)方法在伽辽金方法中，加权函数与近似解展开中所用的函数相同，这样可得到最精确的解。假设近似解为：则取加权函数选为：因此：得到：在算符L为自伴算符的情况下，伽辽金方法与里兹方法得到相同的方程组。第11页/共90页 12二维标量场有限元分析过程二维边值问题Dirichlet边界条件：混合边界条件：Neumann边界条件： 第12页/共90页 13空间离散1235641234这是二维区域离散的示意图。黑色数字表示节点的全局编码，红色数字表示三角形单元的全局编号。组成每个三角形单元的节点在三角形内有一组局部编码。显然，该局部编码与节点的全局编码有一一对应关系。第13页/共90页 14选择插值基函数使用线性三角形单元，在第e个单元内， 可以近似为：节点坐标带入：解得：其中， 为插值基函数第14页/共90页 15插值基函数其中：第15页/共90页 16当观察点(x,y)位于第i个结点的对边上时：二维插值基函数的性质性质1：性质2：一个单元边的 值与其相对结点处的 值无关，而由该边两端点处的 值确定。从而保证了单元两侧解的连续性结论：第16页/共90页 17用伽辽金法建立公式①②③其中：第17页/共90页 18组合成方程组组合：其中：用矩阵表示为：第18页/共90页 19K矩阵的形成1235641234图中箭头所指为相应三角形单元的起始结点1，并且规定结点1、2、3按顺时针排列。第19页/共90页 20列向量[b]的形成1