高阶有限元结合Chebyshev逼近分析宽带散射特性.pdfVIP

高阶有限元结合Chebyshev 逼近分析宽带散射特性1 刘尧，吴先良 安徽大学计算智能与信号处理教育部重点实验室，安徽合肥（230039 ） 摘 要：将切比雪夫逼近理论应用于目标宽带电磁散射特性分析中，通过求解给定频带的切 比雪夫节点并用高阶有限元方法求得各节点处的目标表面场值，实现整个频带内任一频率点 目标表面场值的快速预测，使得目标电磁散射特性得到快速分析。这一方法避免了在给定频 带内对有限元矩阵的反复计算，大大减少了计算机的内存。 关键词：高阶有限元法，边界积分法，Chebyshev 逼近，宽带雷达散射截面 中图分类号：TN82 1. 引言 [1] 应用有限元-边界积分法 （FE-BI ）可以精确的计算目标的雷达散射截面（RCS ），但 每次计算只能得到一个频率点的 RCS 。而很多电磁应用问题都需要获得宽带频率响应，如 在雷达目标识别中，需要目标的宽带 RCS 以产生一维距离像和合成孔径雷达（SAR）像。 为了获得目标的宽带 RCS，应用有限元法就必须在给定频带内的每个频率点上逐点计算， 当目标的RCS 随着频率变化剧烈时，必须以很小的频率间隔计算才能得到精确的频率响应， 这就意味着在整个频带内矩阵方程求解次数的大幅度增加，占用了大量的CPU 时间。 近年来，应用有限元法求解目标宽带雷达散射截面的方法主要是渐进波形展开（AWE ） 技术[2~5]等，该方法的主要思想是通过求解FE-BI 方程组，得到展开频率点的场值，通过泰 勒级数展开以及padé 近似获得频带内任意频率点的目标场值，但该方法的有效计算频带受 到泰勒级数和padé 近似本身误差的限制，而且计算过程中需要存储稠密的阻抗矩阵及其多 阶导数，占用了大量内存。 基于以上考虑，本文将 Chebyshev 多项式逼近理论[6,7]引入到有限元中分析目标的宽带 电磁散射特性。对于给定的频带，通过坐标变换计算出Chebyshev 节点。应用有限元-边界 积分法求解Helmholtz 方程的变分表达式，计算出这些节点处的场值，从而获得该频带内 场的Chebyshev 多项式逼近，通过对导体和介质柱的分析，验证了本文方法的有效性。 2. 基本理论 2.1 有限元分析 对散射体所在区域用有限元进行分析。对E 和H 极化，散射体所在区域的场分别满足 z z 的Helmholtz 方程可以统一表达为： ∂ 1 ∂Φ ∂ 1 ∂Φ⎛ ⎞ ⎛ ⎞ 2 + +k Φ 0 υ ρ ρ ∈Ω （1） ⎜ ⎟ ⎜ ⎟ 0 ( ) ∂x ⎝u ∂x ⎠ ∂y ⎝u ∂y ⎠ 其中，对E 极化，u µ ，υ ε ；对H 极化，u ε ，υ µ 。Ω表示介质柱所在区 z r r z r r 域。在介质柱的表面上，假设场的法向导数为 1 ∂Φ −ψ （2 ） u ∂n 根据变分原理[8]，与（1）、（2 ）式等价的变分表达式为： 1本课题得到国家自然科学基金项目（No）的资助。 -1- δF Φ 0