基于泰勒级数展开的时域积分方程快速算法的制作方法

本发明涉及金属目标瞬态电磁散射特性数值计算技术，属于目标电磁散射特性的快速计算
技术领域：
，具体是一种基于泰勒级数展开的时域积分方程快速算法。
背景技术：
：随着计算电磁学领域研究的深入，传统的频域方法已经不能满足需要。以计算机硬件技术的发展为契机，人们逐步具有了直接在时域对具有宽频带特性的瞬变电磁场的计算分析能力，从而实现了对物理量和物理现象更深刻、更直观的理解。同频域方法相比，在时域求解目标的电磁特性不仅可以直观的揭示目标与电磁波相互作用的机理，而且通过少量的计算就可以获得目标的宽频带信息，这在宽带电磁问题、瞬态电磁问题分析中具有明显的优势。同基于微分方程的时域方法(FDTD、FETD等)相比，基于积分方程的时域方法在求解的未知数数量上具有明显的优势，这是因为积分方程利用格林函数建立源和场的关系，求解区域在边界上，离散后未知数的数量与边界面积成正比。其次，积分方程方法自动满足辐射边界条件，不需要强加吸收边界，而吸收边界是基于微分方程的方法所必需的。由于计算机水平和硬件的限制，传统的时域积分方程方法无法解决大规模的电磁问题。而发展快速算法是解决实际工程问题的必由之路。时域积分方程快速算法的研究是在近20年开展起来的，其中最著名的两种：一是时域平面波算法(PWTD)，二是时域自适应积分方程(TD-AIM)。对于三维目标模型，TD-AIM的计算量和存储量分别为和可见其处理三维目标的效率低下。PWTD算法其原理类似于FMM算法，但是因为时间变量的存在，使得PWTD算法更加复杂，实现起来更加困难。并且该算法的适应性差，当用以求解介质或者有耗媒质等复杂问题时，还需要对算法进行特殊修改。技术实现要素：本发明的目的在于提供一种高效、稳定的分析金属目标的时域电磁散射特性的快速方法。由于对远场部分采用泰勒级数展开重构成聚合、转移、配置的形式进而实现矩阵矢量乘的加速计算，本发明对于求解金属复杂目标散射问题需要更少的计算内存以及计算时间。而且编程相对简单易于实现。实现本发明目的的技术解决方案为：一种基于泰勒级数展开的时域积分方程快速算法，步骤如下：第一步，建立时域电磁场积分方程。金属目标在入射电磁波的作用下，在金属散射体的表面上会产生面感应面电流，感应电流的不断变化，会向外辐射出电磁场，由感应电流产生的电磁场就是散射场，利用入射电磁场和散射电磁场在金属表面满足的边界条件建立时域电磁场积分方程。第二步，将散射体表面上离散得到的子散射体分组。任意两个子散射体间的互耦或自耦根据它们所在组的位置关系而分成近场组对和远场组对。当它们是近场组对时，采用直接数值计算。而当它们为远场组对时，则采用泰勒级数展开成聚合-转移-配置方法计算。第三步，近场阻抗矩阵计算，将金属表面电流密度用空间基函数和时间基函数展开，并在空间域上进行伽辽金测试，时间域上进行点匹配得到矩阵元素值。第四步，远场部分采用泰勒级数展开重构成聚合、转移、配置的形式。第五步，矩阵方程求解以及电磁散射参数的计算。利用时间递推的方式求解每个时刻的电流系数。每个时刻空间某处的散射场贡献来源分为近场组源的贡献和远场组源的贡献。前者直接通过近场矩阵和电流系数相乘得到，后者通过泰勒级数展开成聚合-转移-配置方法快速计算。采用迭代法求解出最终的散射电流系数。并根据互易定理求解出雷达散射截面积。本发明与现有技术相比，其显著优点为：(1)可以快速分析复杂电大目标时域电磁散射特性；(2)远场计算是基于泰勒级数展开的，避免了PWTD算法需要的复杂的时频域转换、谱域积分等操作，编程相对简单易于实现；(3)适用性高，可以只需作较少的改动便可应用于求解介质、有耗、色散等复杂问题的分析中。附图说明图1是本发明的应用于基于泰勒级数展开的时域积分方程快速算法中的多层分组示意图。图2是本发明的应用于基于泰勒级数展开的时域积分方程快速算法中的近远场关系示意图。图3是本发明的应用于基于泰勒级数展开的时域积分方程快速算法中的场源基函数位置向量分布示意图。图4是本发明的应用于基于泰勒级数展开的时域积分方程快速算法中的简易导弹模型不同频点处双站RCS比较图。(a)30MHz；(b)150MHz；(c)270MHz。具体实施方式下面结合附图对本发明作进一步详细描述。本发明一种基于泰勒级数展开的时域积分方程快速算法，步骤如下：第一步，建立时域电磁场积分方程。金属目标在入射电磁波{Einc(r,t),Hinc(r,t)}的作用下，在金属散射体的表面上会产生面感应面电流J(r,t)，感应电流的不断变化，会向外辐射出电磁场，由感应电流产生的电磁场就是散射场{Esca(r,t),Hsca(r,t)}，利用入射电磁场和散射电磁场在金属表面满足的边界条件建立时域电磁场积分方程。αn^(r)×n^(r)×Einc(r,t)+(1-α)ηn^(r)×Hinc(r,t)=αLe{J(r,t)}+η(1-α)Lh{J(r,t)}---(1)]]>Le{J(r,t)}=n^(r)×n^(r)×μ4π∫∫sdS1R∂J(r′,τ)∂t-n(r)×n^(r)×14πϵ∫∫sdS∫0τ▿′·J(r′,τ)Rdt---(2)]]>Lh{J(r,t)}=12J(r,t)-n^(r)×14π∫∫sdS▿×J(r′,τ)R---(3)]]>其中是单位外法向矢量，α取值为0～1之间的实数，μ和ε分别是自由空间的磁导率和介电参数。η是自由空间的波阻抗。R＝|r-r'|,c是自由空间中的光速,τ＝t-R/c是延时。第二步，将散射体表面上离散得到的子散射体分组。我们借鉴MLFMM里面对空间基函数进行分组的思想，将基函数的相互作用转换为组的相互作用。假设空间中有一个刚好可以把整个目标物体包围的立方体，把这个立方体等分成8个子立方体，接着在将每个子立方体等分成8个更小的立方体，如图1(a-d)所示。依次类推，直到达到预先设置的门限值，停止划分。定义第一层为划分最细的一层，第一层中的立方体中如果包含基函数，定义为第一层组，对于其它层组，依次类推即可。任意两个子散射体间的互耦或自耦根据它们所在组的位置关系而分成近场组对和远场组对。对于每一层，我们都会设定一个参数β,4＜β＜6。首先从最高层Nl的各组着手划分，如果两组中心间的距离大于规定的的门限值β，就把其称之为远场组对；然后在其子层中，组对中心间距大于规定的该层的门限值，同时这两个组都没有划分到Nl层的远场组对，那么就规定这样的组对为该层的远场组对；根据以上所述依次类推，对于每一层我们就可以得到其远场组对的信息，如图2所示。当它们是近场组对时，采用直接数值计算。而当它们为远场组对时，则采用泰勒级数展开成聚合-转移-配置方法计算。第三步，近场阻抗矩阵计算。把金属表面电流电流密度J(r,t)用RWG基函数Λn(r)作为空间基函数，三角基函数Tj(t)作时间基函数展开：J(r,t)=Σn=1NSΣj=1NtIn,jΛn(r)Tj(t)---(4)]]>其中Ns是散射体包含的RWG基函数的个数，Nt是时间基函数的个数，In,j是第n个RWG基函数第j个时间步上时间基函数的系数。将式(4)代入(1)，在空间上进行Galerkin测试，在时间上进行点匹配，最终得到如下的矩阵方程：Z0Ii=Vi-Σj=1i-1Zi-jIj---(5)]]>Zi-j=αZEi-j+(1-α)ZMi-j---(6)]]>[ZEi-j]mn=μ04π∫Sm∫SnΛm(r)·Λn(r′)∂τTj(iΔt-R/c)Rds′ds+14πϵ0∫Sm∫Sn▿·Λm(r)▿·Λn(r′)∂τ-1Tj(iΔt-R/c)Rds′ds---(7)]]>[ZMi-j]mn=12∫SmΛm(r)·Λn(r′)Tj(iΔt)ds-14π∫SmΛm(r){n^(r)×∫Sn∂τTj(iΔt-R/c)Λm(r′)×RR2+Tj(iΔt-R/c)Λm(r′)×RR3}ds′ds---(8)]]>第四步，远场部分采用泰勒级数展开重构成聚合、转移、配置的形式。现在通过计算源点r'处的源信号Jn(r',t)对场点r处的辐射贡献来描述该算法的基本原理。假设源点r'所在的空间基函数为Λn(r')，因此r'处的源信号Jn(r',t)可以展开如下：Jn(r′,t)=Σj=1NtIn,jΛn(r′)Tj(t)=Λn(r′)gn(t)---(9)]]>实际操作中我们需要将源信号分解为持续时间更短的分段子信号，保证每段子信号的持续时间满足约束条件。源信号Jn(r',t)可以被分解为L段连续的子信号Jn,l(r',t)，每一段子信号的持续时间为Ts＝(Mt+1)Δt。源信号可以写成如下形式：Jn(r′,t)=Σl=1LJn,l(r′,t)=Σl=1LΛn(r′)gn,l(t)---(10)]]>那么源点r'处第l段子信号在场点r处产生的测试场为：⟨Λm(r),nm×Hn,l(r,t)⟩=14π∫SmΛm(r)·nm×{∫SndS′{[∂τgn,l(τ)Λn(r′)]×RcR2+[gn,l(τ)Λn(r′)]×RR3}}dS≈14πm(∫SmΛm(r)dS·nm×∫SnΛn(r′)dS′×R)[1cR2∂τgn,l(τ)+1R3gn,l(τ)]=14π(mm·nm×mn×R)[1cR2∂τgn,l(τ)+1R3gn,l(τ)]---(11)]]>⟨Λm(r),En,l(r,t)⟩=∫Sm∫SnΛm(r)·Λn(r′)μ0∂τgn,l(τ)4πRds′ds+∫Sm∫Sn▿·Λm(r)▿·Λn(r′)∂τ-1gn,l(τ)4πϵ0Rds′ds≈∫SmΛm(r)ds·∫SnΛn(r′)ds′μ0∂τgn,l(τ)4πR+∫Sm▿·Λm(r)ds∫Sn▿·Λn(r′)ds′∂τ-1gn,l(τ)4πϵ0R=μ04πmm·mn∂τgn,l(τ)R±lmln4πϵ0∂τ-1gn,l(τ)R---(12)]]>其中mm=∫SmΛm(r)ds,mn=∫SnΛn(r′)ds′.]]>现假设源点rn与场点rm分别位于两个组内，如图3所示，两个组分别称为源组和场组，组中心分别为ri和rj，场源基函数之间的矢量可以表示为：R＝rmi+rij-rnj＝Rm-Rn(13)这里，rij＝ri-rj，rmi＝rm-ri，rnj＝rn-rjRm＝rmi+rij/2，Rn＝rnj-rij/2利用泰勒级数展开可以得到如下表达式：Rα=(R·R)α2=[(rmi+rij-rnj)·(rmi+rij-rnj)]α2=rijα[1+(2rmi·r^ijrij+rmi2rij2)+(2rnj·r^jirij+rnj2rij2)-2rmi·r^njrij2]α2≈Rm(α)+Rn(α)---(14)]]>Rm(α)=rijα[12+α(rmi·r^ijrij+rmi2+(α-2)(rmi·r^ij)22rij2)]---(15)]]>Rn(α)=rijα[12+α(rnj·r^jirij+rnj2+(α-2)(rnj·r^ji)22rij2)]---(16)]]>则式(11)(12)可以写成聚合、转移、配置的形式。第五步，矩阵方程求解以及电磁散射参数的计算。利用时间递推的方式求解每个时刻的电流系数。现将任意空间上的散射场按来源分为两部分:一部分由该场点所在组的近场组NFP(α)中的源产生的；另一部分由该场点所在组的远场组FFP(α)中的源产生。基于泰勒级数展开的时域积分方程快速算法的递推公式变为：Z0Ii=Vi-Σα′∈NFP(α)Σj=0i-1Zi-jαα′Ijα′-Σα′∈FFP(α)Σn∈α′(n)⟨Λm(r),αEn,l(r,t)+(1-α)nm×Hn,l(r,t)⟩---(17)]]>近场区域产生的贡献是由经典MOT算法计算得到，主要通过矩阵元素值与电流系数相乘得到。在远场区域，也就是对应组中的远场组对，这些组之间的相互作用，过泰勒级数展开成聚合-转移-配置方法快速计算。由于时刻iΔt以前的Ij在时刻iΔt都是已知的，j＝1,2,3,...,i-1，这样每个时间步求解一次式(17)的矩阵方程，就可以得到每个时刻iΔt上的Ii所以可以递推求出各时刻的电流值。最后可以根据求得的瞬态电流系数计算出我们需要的电磁散射参数。为了验证本发明的正确性与有效性，下面分析了简易导弹模型的的电磁散射特性。算例：简易导弹模型，几何尺寸为长8m，半径0.25m。激励源设置为：调制高斯脉冲，中心频率150MHz，带宽300MHz，入射波的方向θ＝180°,(弹头方向照射)，极化方式VV极化，观察角度0°≤θ≤180°,导弹模型采用0.1λ(λ＝1m)剖分得到8034个三角形，总未知量为12051。分组大小为0.3λ。3个盒子之外为远场区域，近场门限0.3λ*3。使用Δt＝333.3ps，一共计算了600个时间步。CFIE的参数α＝0.5，采用GMRES迭代方法计算，收敛精度为1e-9。图4给出了导弹模型在30MHz、150MHz以及270MHz时的双站RCS曲线并与经典MOT算法计算结果进行对比，表格1统计了本专利提出的方法与传统的时域积分方法消耗的时间和内存对比。表1.计算时间和内存消耗对比从图4中可以看出本发明方法计算的双站RCS数据和MOT计算的结果十分吻合，证明了本发明方法的精确性。从表1看出，采用本发明方法相比经典MOT算法需要更好的计算内存和计算时间。进一步证明了本发明方法的有效性。本发明的实现过程简单，相较于需要的复杂的时频域转换、谱域积分等操作，本发明方法只需在八叉树分组以及多层近远场划分的基础上，对远场组对之间的距离R进行泰勒展开重构成聚合、转移、投射的操作，编程相对简单易于实现。而且本发明方法适用性高，可以只需作较少的改动便可应用于求解介质、有耗、色散等复杂问题的分析中。当前第1页1 2 3