第9章有限元法.PDFVIP

第９章　 有 限 元 法 与有限差分法类似，有限元法（ＦＥＭ）也是一种将偏微分方程转化为线性代数方程组，进 而求解边值问题的数值方法。 不同的是，有限差分法的基本思想是对微分算子进行近似，而 有限元法是对偏微分方程的解进行近似。 有限元法最早由Ｃｏｕｒａｎｔ于１９４３年提出，以求解势 ［１］ 论中的变分问题 ， 自此之后该方法得到极大的发展，并被广泛应用于结构力学分析及越 来越多的其他领域问题的求解。 在今天，有限元法被认为是一种适用于处理众多工程和数 学问题的通用方法，其中包括微波工程和电磁学中的问题。 ［２］ １９６９年，Ｓｉｌｖｅｓｔｅｒ使用有限元法分析了空心波导中波的传播 ，这是有限元法第一次应用于 微波工程和电磁学中。 自此，研究人员很快认识到这种方法的重要性，并将其成功地应用于求解 各种静电、静磁，以及介质加载波导的问题。 １９７４年，Ｍｅｉ开发了一种将特征函数展开与有限元 ［３］ 法相结合的技术，用于处理如天线和散射分析等开放区域的电磁问题 。 １９８２年，Ｍａｒｉｎ又发展 ［４］ 了一种将有限元法和边界积分方程相结合的方法，用于处理开放区域的散射问题 。 ２０世纪８０年代，随着基于棱边的矢量基函数的提出和发展，有限元法在对矢量电磁场 问题的分析中取得了突破性的进展［５～７］。 矢量有限元可以精确地模拟电场和磁场，并消除 许多由传统的基于节点的标量元引起的问题。 随着矢量有限元的发展，有限元法成为了计 算电磁学中一种非常强大的数值方法。 今天，有限元法是设计天线和微波器件的主要辅助 ［８］ ［９］ 工具，其基本原理及各种应用在诸如Ｓｉｌｖｅｓｔｅｒ ａｎｄ Ｆｅｒｒａｒｉ 及Ｊｉｎ 的许多著作中均有论述。 本章首先通过一个简单的一维算例介绍有限元法的基本原理，然后详细描述用有限元 法求解频域标量和矢量问题的步骤，接着介绍时域有限元法。 对于每一种应用，我们会给出 数值算例来展示有限元法的性能。 本章最后讨论如何使用吸收边界条件（ＡＢＣ）和理想匹配 层（ＰＭＬ）来截断计算区域，从而将有限元法应用到分析无界空间中的电磁问题。 我们还简 单讨论在编程实现有限元法时的一些数值问题。 ９１　 有限元法概述 本节首先简单地介绍加权残差法，然后求解一个简单的一维亥姆霍兹方程，以此来阐明 有限元法的基本原理。 ９１１　 一般原理 使用加权残差法或变分法均可以建立有限元法的公式。 加权残差法直接从边值问题的 偏微分方程出发，而变分法则从边值问题的变分形式出发。 在本章中，我们使用加权残 差法。 为了介绍加权残差法的基本原理，考虑如下的偏微分方程： （９１１） 式中， 为微分算子，φ为待求未知解，ｆ 为源函数。 为了寻求φ 的解，首先用如下一组基函 数将其展开： ３８ 计算电磁学（第二版） （９１２） ＝ 式中，ｖ （ｊ １，２，…，Ｎ）为基函数，其线性组合可表示未知解；ｃ 为相应的未知展开系数。 加ｊ ｊ 权残差法尝试确定ｃ 的方法如下：将式（９１２）代入式（９１１），然后将得到的等式乘以加权ｊ 函数ｗ ，并在整个求解区域Ω 中积分，这样可得ｉ （９１３） 当给定一组加权函数，上式就定义了一个代数方程组，在满足边界条件的要求下求解此方程