求解电磁场有限元 边界元方程组的有效方法.docVIP

求解电磁场有限元 边界元方程组的有效方法 摘要: 提出了一种求解电磁场有限元 边界元混合法所生成的线性方程组的有效方法 — — — 内观法结合多波前法.由于该线性方程组的系数是一个部分稀疏部分满填充的矩阵,为了加速求解,应用内观法将系数矩阵分为 2块,一块是有限元法形成的稀疏矩阵,另一块是边界元法生成的满阵,然后用多波前法求解稀疏矩阵方程,用高斯约当消去法解满阵方程.采用该方法,计算了二维多层介质柱体的雷达散射截面.计算结果表明,该方法的计算效率远远高于传统的高斯法. 关键词: 有限元法;边界元法;电磁散射;雷达散射截面;内观法;多波前法 　　有限元法是 20世纪 40年代提出的近似求解数理边值问题的一种数值方法.随着计算机技术的迅速发展,它已被广泛地应用于结构分析、 机械制造、 建筑设计、 电磁分析等领域.有限元法对场域的剖分具有很大的灵活性,对不规则、 非均匀的问题具有广泛的适应性,而且形成的是对称、稀疏的带状矩阵,可被高效率地存储和求解. 但是这种方法在分析开放域问题时有困难,因为区域为无限大,剖分时将有无穷多个节点, 形成的矩阵方程将无法求解。. 边界元法把描述场的微分方程通过加权余量法归结为边界上的积分方程,然后对积分方程进行边界分割及插值,从而求得近似解,因此可以使求解的问题降低一维,数据输入简单. 由于边界元法是通过使用适当的G reen函数考虑了Somm erfeld辐射条件,因而这种方法特别适用于分析无限大的开放区域。但是对于不均匀介质问题,由于难以求得 G reen函数,因而应用边界元法求解这类问题时会遇到困难.为了克服有限元与边界元法的缺点,同时保留它们的优点,人们发展了一种新的方法,即有限元与边界元混合法,它在科学和工程计算中有着独特的作用.但用有限元 边界元混合法分析问题时,最后所形成的是一个部分稀疏部分满填充的矩阵,该矩阵不能被高效地存储,同时如果直接求解该矩阵,对于节点数较多的情形,耗费的机时是很多的. 本文提出了一种求解有限元边界元混合法方程组的新方法 — — — 内观法[ 6 ] 结合多波前法[ 7, 8 ],它不直接求解所得矩阵,而是用内观法将原矩阵分为2块,一块是有限元法形成的对称稀疏矩阵,另一块是边界元法生成的满阵,然后用多波前法求解大的对称稀疏矩阵方程,用高斯约当消去法解一个小的满阵方程.多波前法是根据最低阶算法,将稀疏矩阵分解为许多较小的子域,在每个子域内做矩阵分解,消去内部未知量,形成仅与子域边界未知量相关的波前矩阵,再利用降阶算法消去所有相邻子域间的公共未知量. 为了说明该方法的有效性,本文给出了用有限元边界元混合法分析多层介质组成的柱体散射特性的一般方法,并计算其雷达散射截面.计算结果表明,本文方法的效率明显高于高斯 约当消去法. 1、理论分析 首先用有限元边界元混合法分析二维多层介质组成的柱体散射特性,形成矩阵方程,然后用内观法结合多波前法求解该方程组.图 1给出了一任意截面形状的介质柱,它由多层介质组成,假设第 i层介质的相对介电常数和相对磁导率分别为εr , i ,μr , i ,所在区域为Ωi ,最外一层介质所在区域为Ω1 ,其边界为Γ,介质层外的无限大区域为Ω∞. 设有一束 TM 极化的均匀平面电磁波投射到柱体上,入射电场和散射电场都沿 z轴取向,电磁场随时间变化因子为 exp ( j ωt) . 图 1　多层介质组成的二维柱体的横截面 111　有限元公式 对介质层所在区域用有限元法进行分析. 用 表示电场的 z分量 Ez ,则由麦克斯韦方程 可推得, 第 i层介质中的电场 i 满足下面的 He lm ho ltz方程: 0 =ω2 μ0ε 0 ,μ 0 ,ε 0分别为自由空间的磁导率 和介电常数. 两层介质的分界面上 , 　i = i+1 (2) 式中, n为边界的法向方向. 在边界Γ上,令 与式 (1) ～ (3) 等价的变分表达式为 对于最外一层介质 F ( 1 ) = 对于里面的介质层 　　将介质层所在区域Ω剖分为M个小三角形单元,同时边界Γ被分成Ms条线段,假设Ω中共有 N个节点.对每个单元进行插值,并对所有单元求和,然后对每个节点处的场求偏微分,再进行组合,可得 K φ + Cψ = 0 (7) 式中,φ为所有节点处的场组成的列向量;ψ为边界Γ上所有节点处场的法向导数组成的列向量;K, C的计算公式参见文献 [ 9 ]. 显然 K是一个对称、 稀疏矩阵,式 (7)含有 N +Ms个未知量,但只有N个方程. 112　边界元公式 对于介质柱外的区域Ω∞,由于它是一个开放域问题,因此应用边界元法进行分析. 由麦克斯韦方程可推得区域Ω∞中的场满足 下面的 He lm ho ltz方程: 0φ = jk 0 Z0 J z (ρ) (8) 式中, Jz (ρ) 为产生