高阶有限元－边界积分法在二维散射中的应用.docVIP

高阶有限元－边界积分法在二维散射中的应用 　　摘 要：以二阶、三阶基函数为例，应用高阶有限元-边界积分法分析了二维散射体电磁散射特性。计算了几种二维方柱(导体和介质)的雷达散射截面，结果与矩量法一致，对三种数值结果进行了误差分析。结果表明：高阶有限元-边界积分法比一阶有限元-边界积分法有着更高的计算精度、收敛速度和计算效率。 　　关键词：高阶有限元;边界积分法;二维散射体;雷达散射截面 　　中图分类号：TN82 文献标识码：B 　　文章编号：1004373X(2008)0500104 　　 　　Application of High-order Finite Element-Boundary Integral 　　Method to 2D Scattering Problems 　　 　　LIU Yao,PANG Yanhong,WU Xianliang 　　(Department of Electronic Science and Technology,Anhui University,Hefei,230039,China) 　　 　　Abstract：Using the second-order and third-order basic functions,the high-order Finite Element Method (FEM) in conjunction with Boundary Integral Method (BIM) is used to analyze electromagnetic scattering characteristics of two-dimensional objects.As the examples,the Radar Cross Sections (RCS) of some two-dimensional square cylinders(conductor and dielectric) are calculated,the results agree well with those obtained by Method of Moments (MOM).The error of result is analyzed.Compared with first-order FEM-BIM,numerical results suggest that the high-order FEM-BIM has higher precision,faster convergence rate and more efficiency. 　　 　　Keywords：high-order finite element;boundary integral method;two-dimensional objects;radar cross section 　　 　　1 引 言 　　 　　有限元法是近似求解数理边值问题的一种数值方法，他于20世纪40年代提出,随着计算机技术的迅速发展，已被广泛地应用在结构分析、机械制造、建筑设计、电磁分析等领域。由于他对场域的剖分具有很大的灵活性，因此对不规则、非均匀的物体具有广泛的适用性，而且形成的矩阵为稀疏矩阵，人们可以利用这一特点快速求解含未知量的稀疏矩阵方程???? [1，2]。但有限元在分析开域问题时遇到困难，因为区域的无限大，有限元剖分时将有无穷多个节点，所形成的矩阵方程将很难求解[3]。多年来对于二维开域问题的研究概括起来主要分为两大类[4-8]：一类是利用吸收边界条件或完全匹配层技术将无限大空间截断，使原散射问题限制在被研究物体周围一个有限的范围内，然后采用有限元进行分析；第二类是应用矩量法或边界积分方法来处理柱体外的无限大空间，柱体所在区域采用有限元法来分析。第一种方法缺点在于扩大了研究区域，产生了较多的节点数，特别是吸收边界的形状以及他与散射体间的距离对于解的精度有很大影响。第二种方法精度较高，产生的节点数较少，但生成的是部分稀疏，部分满秩的矩阵，如果能够快速求解该方程，则该方法有着较好的应用前景，特别对于电尺寸较大的散射体。 　　本文在已有文献工作的基础上，应用高阶有限元-边界积分法分析二维散射体的电磁散射特性。对散射体内、外区域分别应用高阶有限元和边界积分法进行分析。然后通过场的连续性进行耦合，形成待求矩阵方程组，求解该方程组，获得每个结点的场值，进一步求解雷达散射截面。作为算例，分别应用二阶和三阶有限元-边界积分法计算了几种二维方柱(导体和介质)的雷达散射截面，并将结果与矩量法的数值结果作了误差分析。数值结果表明，高阶有限元-边界积分法较之一阶线性有限元-边界积分法有着更高的计算精度、收敛速度和更