复杂电磁问题的限元、边界积分及混合算法的快速分析技术.pdf

博士论文 快速算法在电磁场有限元仿真中的应用 摘要 由于工程电磁场应用需求的提高，不断追求高效而精确的数值分析方法成为计 算电磁学领域一直以来的研究工作重点。本文针对矢量有限元方法、边界积分方法 及其混合算法开展了一系列研究工作，重点研究了这三种数值方法在电磁散射问题 分析中的应用。 本文首先从矢量有限元基本原理出发，研究了高阶单元以及高阶叠层基函数 的引入对数值精度的改善效果。针对高阶有限元中求解稀疏矩阵的问题，发展了两 网格算法与FGMRES算法相结合的加速技术。 其次将高阶矢量有限元方法结合完全匹配层(PML)技术应用于电磁散射问题分 析。为了有效降低有限元计算区域，发展了一种基于高阶单元的局部共形PML技 术，它能够贴近复杂结构散射体的外廓形状设置截断区域，同时保证对复杂几何边 界的拟合精度。 ‘有限元方法用于分析开域问题时，截断边界的设置往往由于其近似性会影响 到对散射特性的精确分析，因此广泛采用的是边界积分方程类方法分析散射问题。 但是由于该类方法分析开放结构的电磁散射时矩阵性态差，迭代收敛效率很低。本 文研究了两种预条件技术用于加速该类方法的迭代求解。一种是基于多分辨基函数 及其特征谱信息的双步预条件技术，利用多分辨基函数的分层特性，用粗层获得的 特征谱信息构造预条件矩阵，优化基于多分辨基函数的预条件性能。另一种为基于 现特征值聚集在零和无穷大的情况，使得阻抗矩阵不依赖离散密度且有良好的矩阵 性态。 本文最后一部分针对有限元．边界积分混合算法做了一些研究。首先将基于散 度共形和旋度共形的高阶叠层基函数引入到该混合算法中，提高了数值仿真的精 度。根据有限元与边界积分对网格离散要求存在的差异性，以及高阶叠层基函数的 叠层特性，发展了一种高阶FE．低阶BI的混合阶FE．BI算法，可以有效地降低未 知量和提高计算精度。其次利用有限元能够灵活处理复杂媒质(如各向同／异性、双 各向同／异性媒质)，成功地将FE．BI方法用于分析含复杂媒质目标的散射特性。利 用边界积分对格林函数的依赖，结合实镜像快速多极子实现对半空间环境下目标的 散射分析。最后构造了一种由有限元矩阵和边界积分近场部分组成的高效预条件技 术。该预条件矩阵可以很好的近似待求解矩阵，并且易于构造，加速对半空间散射 问题分析时的迭代求解。由于预条件后的FE．BI求解过程中，需要进行预条件矩阵 摘要 博士论文 GMRES在FE．BI内迭代求解中的收敛表现。同时也在计算平台允许的情况下，将 一种基于H．matrix的快速直接求解算法引入到内迭代求解，分析了其内存消耗和 计算复杂度等特性，相比其他求逆的直接算法，其在计算时间和内存消耗上具有一 定的优势。 关键词：有限元方法， 矩量法，快速多极子，Krylov子空间迭代法，预条件 技术，高阶叠层基函数，双各向同／异性媒质 n 博士论文 快速算法在电磁场有限元仿真中的应用 Abstract havebeenseenaconsiderableinresearchon Computationalelectromagnetics surge efficientaccuratenumerical havebeenstimulatedthedemandsfor methods，which by simulations．Theresearchofthisdissertationisfocusonthefinite electromagnetic targets element thefiniteelement— method(FEM)，theboundaryintegralmethod(BI)and boundary