CN106294283A - 基于泰勒级数展开的时域积分方程快速算法 - Google Patents
基于泰勒级数展开的时域积分方程快速算法 Download PDFInfo
- Publication number
- CN106294283A CN106294283A CN201510266484.9A CN201510266484A CN106294283A CN 106294283 A CN106294283 A CN 106294283A CN 201510266484 A CN201510266484 A CN 201510266484A CN 106294283 A CN106294283 A CN 106294283A
- Authority
- CN
- China
- Prior art keywords
- tau
- prime
- centerdot
- lambda
- alpha
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于泰勒级数展开的时域积分方程快速算法。传统的时域积分方程方法由于计算时间和内存消耗的限制无法解决大规模的电磁问题。而发展快速算法是解决实际工程问题的必由之路。本发明是利用自由空间中格林函数的泰勒级数展开重构成聚合、转移、配置的形式进而实现矩阵矢量乘的加速计算,本发明对于求解金属复杂目标散射问题需要更少的计算内存以及计算时间。而且编程相对简单易于实现,具有很强的实际工程应用价值。
Description
技术领域
本发明涉及金属目标瞬态电磁散射特性数值计算技术,属于目标电磁散射特性的快速计算技术领域,具体是一种基于泰勒级数展开的时域积分方程快速算法。
背景技术
随着计算电磁学领域研究的深入,传统的频域方法已经不能满足需要。以计算机硬件技术的发展为契机,人们逐步具有了直接在时域对具有宽频带特性的瞬变电磁场的计算分析能力,从而实现了对物理量和物理现象更深刻、更直观的理解。同频域方法相比,在时域求解目标的电磁特性不仅可以直观的揭示目标与电
磁波相互作用的机理,而且通过少量的计算就可以获得目标的宽频带信息,这在宽带电磁问题、瞬态电磁问题分析中具有明显的优势。
同基于微分方程的时域方法(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)},利用入射电磁场和散射电磁场在金属表面满足的边界条件建立时域电磁场积分方程。
其中是单位外法向矢量,α取值为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)作时间基函数展开:
其中Ns是散射体包含的RWG基函数的个数,Nt是时间基函数的个数,In,j是第n个RWG基函数第j个时间步上时间基函数的系数。将式(4)代入(1),在空间上进行Galerkin测试,在时间上进行点匹配,最终得到如下的矩阵方程:
第四步,远场部分采用泰勒级数展开重构成聚合、转移、配置的形式。
现在通过计算源点r'处的源信号Jn(r',t)对场点r处的辐射贡献来描述该算 法的基本原理。假设源点r'所在的空间基函数为Λn(r'),因此r'处的源信号Jn(r',t)可以展开如下:
实际操作中我们需要将源信号分解为持续时间更短的分段子信号,保证每段子信号的持续时间满足约束条件。源信号Jn(r',t)可以被分解为L段连续的子信号Jn,l(r',t),每一段子信号的持续时间为Ts=(Mt+1)Δt。源信号可以写成如下形式:
那么源点r'处第l段子信号在场点r处产生的测试场为:
其中
现假设源点rn与场点rm分别位于两个组内,如图3所示,两个组分别称为源组和场组,组中心分别为ri和rj,场源基函数之间的矢量可以表示为:
R=rmi+rij-rnj=Rm-Rn (13)
这里,rij=ri-rj,rmi=rm-ri,rnj=rn-rj Rm=rmi+rij/2,Rn=rnj-rij/2
利用泰勒级数展开可以得到如下表达式:
则式(11)(12)可以写成聚合、转移、配置的形式。
第五步,矩阵方程求解以及电磁散射参数的计算。
利用时间递推的方式求解每个时刻的电流系数。现将任意空间上的散射场按来源分为两部分:一部分由该场点所在组的近场组NFP(α)中的源产生的;另一部分由该场点所在组的远场组FFP(α)中的源产生。基于泰勒级数展开的时域积分方程快速算法的递推公式变为:
近场区域产生的贡献是由经典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进行泰勒展开重构成聚合、转移、投射的操作,编程相对简单易于实现。而且本发明方法适用性高,可以只需作较少的改动便可应用于求解介质、有耗、色散等复杂问题的分析中。
Claims (2)
1.一种基于泰勒级数展开的时域积分方程快速算法,其特征在于步骤如下:
第一步,建立时域电磁场积分方程;利用入射电磁场和散射电磁场在金属表面满足的边界条件建立时域电磁场积分方程;
第二步,将散射体表面上离散得到的子散射体分组;用一个立方体盒子将整个目标物体包围,把这个立方体等分成8个子立方体,接着将每个子立方体等分成8个更小的立方体,依次类推,直到达到预先设置的门限值,停止划分;任意两个子散射体间的互耦或自耦根据它们的位置关系而分成近场作用对和远场作用对;当它们是近场作用对时,计算近场阻抗矩阵,当它们为远场作用对时,采用泰勒级数展开成聚合-转移-配置方法计算源组基函数在场组基函数处产生的场;
第三步,计算近场阻抗矩阵,将金属表面电流密度用空间基函数和时间基函数展开,并在空间域上进行伽辽金测试,时间域上进行点匹配得到矩阵元素值;
第四步,远场作用对之间采用泰勒级数展开重构成聚合、转移、配置的形式计算源组基函数在场组基函数处产生的场;
第五步,矩阵方程求解以及电磁散射参数的计算;利用时间递推的方式求解每个时刻的电流系数,采用迭代法求解出最终的感应电流系数,最后根据求得的瞬态电流系数计算出需要的电磁散射参数。
2.根据权利要求1所述的基于泰勒级数展开的时域积分方程快速算法,其特征在于:所述步骤四中,在八叉树分组以及多层近远场划分的基础上,对远场作用对之间进行泰勒级数展开重构成聚合、转移、投射的操作来加速矩阵矢量乘;
设源点r'所在的空间基函数为Λn(r'),r'处的源信号Jn(r',t)展开如下:
其中空间基函数Λn(r)为RWG基函数,时间基函数Tj(t)为三角基函数,Nt是时间基函数的个数;
源信号Jn(r',t)被分解为L段连续的子信号Jn,l(r',t),每一段子信号的持续时间为Ts=(Mt+1)Δt,Mt为每段子信号的长度,源信号写成如下形式:
源点r'处第l段子信号在场点r处产生的测试电磁场为:
其中Λm(r)为场点r处的测试基函数,nm为场点r处的单位外法向矢量,En,l(r,t),Hn,l(r,t)为源点r'处第l段子信号在场点r处产生的电磁和磁场,μ和ε分别是自由空间的磁导率和介电常数;R=|r-r'|,c是自由空间中的光速,τ=t-R/c是延时;lm,ln分别为第m和n条边的边长,为基函数在其支撑域内的积分;
源点rn与场点rm分别位于两个组内,两个组分别称为源组和场组,组中心分别为ri和rj,场源基函数之间的矢量表示为:
R=rmi+rij-rnj=Rm-Rn (5)
这里,rij=ri-rj,rmi=rm-ri,rnj=rn-rj Rm=rmi+rij/2,Rn=rnj-rij/2
利用泰勒级数展开得到如下表达式:
则将(5-8)代入式(3)(4)将写成聚合、转移、配置的形式来加速矩阵矢量乘。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510266484.9A CN106294283B (zh) | 2015-05-22 | 2015-05-22 | 基于泰勒级数展开的时域积分方程快速方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510266484.9A CN106294283B (zh) | 2015-05-22 | 2015-05-22 | 基于泰勒级数展开的时域积分方程快速方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106294283A true CN106294283A (zh) | 2017-01-04 |
CN106294283B CN106294283B (zh) | 2019-05-07 |
Family
ID=57633456
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510266484.9A Active CN106294283B (zh) | 2015-05-22 | 2015-05-22 | 基于泰勒级数展开的时域积分方程快速方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106294283B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107679287A (zh) * | 2017-09-11 | 2018-02-09 | 三峡大学 | 基于3步4阶隐式泰勒级数法的电磁暂态数值计算方法 |
CN111159637A (zh) * | 2020-01-02 | 2020-05-15 | 西北工业大学 | 一种应用于磁化等离子体计算的电磁波时域精细积分方法 |
WO2024140686A1 (zh) * | 2022-12-27 | 2024-07-04 | 宁波德图科技有限公司 | 电磁近场计算方法及存储介质 |
-
2015
- 2015-05-22 CN CN201510266484.9A patent/CN106294283B/zh active Active
Non-Patent Citations (6)
Title |
---|
JIADE YUAN等: "Efficient Generation of Method of Moments Matrices Using Equivalent Dipole-Moment Method", 《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》 * |
JUNHO YEO等: "Efficient Generation of Method of Moments Matrices Using the Characteristic Function Method", 《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》 * |
XINLEI CHEN等: "Fast Dipole Method for Electromagnetic Scattering From Perfect Electric Conducting Targets", 《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》 * |
姚雨帆等: "应用EDM法和ACA算法快速计算电大开缝导体RCS", 《合肥工业大学学报(自然科学版)》 * |
李金艳: "时域平面波算法研究", 《中国优秀硕士学位论文全文数据库-信息科技辑》 * |
秦宇晨: "时间步进时域体面积分方程及其快速方法的研究", 《中国优秀硕士学位论文全文数据库-信息科技辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107679287A (zh) * | 2017-09-11 | 2018-02-09 | 三峡大学 | 基于3步4阶隐式泰勒级数法的电磁暂态数值计算方法 |
CN107679287B (zh) * | 2017-09-11 | 2021-07-13 | 三峡大学 | 基于3步4阶隐式泰勒级数法的电磁暂态数值计算方法 |
CN111159637A (zh) * | 2020-01-02 | 2020-05-15 | 西北工业大学 | 一种应用于磁化等离子体计算的电磁波时域精细积分方法 |
CN111159637B (zh) * | 2020-01-02 | 2023-06-20 | 西北工业大学 | 一种应用于磁化等离子体计算的电磁波时域精细积分方法 |
WO2024140686A1 (zh) * | 2022-12-27 | 2024-07-04 | 宁波德图科技有限公司 | 电磁近场计算方法及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106294283B (zh) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102156764B (zh) | 一种分析天线辐射和电磁散射的多分辨预条件方法 | |
CN102129523B (zh) | 基于mda和mlssm的分析复杂目标电磁散射的方法 | |
Fan et al. | Multidomain pseudospectral time-domain simulations of scattering by objects buried in lossy media | |
CN103279601B (zh) | 导体目标宽带电磁散射特性的仿真方法 | |
Guan et al. | Electromagnetic scattering from objects above a rough surface using the method of moments with half-space Green's function | |
CN104317984B (zh) | 基于分域建模的船舶电磁散射预测方法及系统 | |
CN102708229A (zh) | 复杂分层媒质结构的矩阵分解结合新奇异值分解方法 | |
CN102592057B (zh) | 周期结构指定频率的本征分析方法 | |
CN104573376A (zh) | 一种时域有限差分法计算电磁散射的瞬态场远场外推方法 | |
CN106446470A (zh) | 一种高效并行的非均匀介质频域有限差分方法 | |
CN106096267A (zh) | 一种腔体电磁散射特性快速计算方法 | |
CN106294283A (zh) | 基于泰勒级数展开的时域积分方程快速算法 | |
CN102054094A (zh) | 一种平面微带电路的快速多层方向性仿真方法 | |
Shi et al. | Comparison of interpolating functions and interpolating points in full-wave multilevel Green's function interpolation method | |
CN103412988A (zh) | 基于相移降阶模型周期结构的三维电磁场仿真模拟方法 | |
CN104731996A (zh) | 一种快速提取电大尺寸金属腔体目标瞬态散射信号的仿真方法 | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
CN104915465A (zh) | 基于延迟拉盖尔多项式的金属目标瞬态电磁散射分析方法 | |
CN103279612B (zh) | 复杂目标雷达回波快速获取的多重网格预条件方法 | |
CN105205299B (zh) | 电大目标电磁散射特性快速降维分析方法 | |
CN105760343A (zh) | 分析复杂多金属目标电磁散射特性的时域积分方法 | |
CN104915324A (zh) | 腔体含介质目标电磁散射混合分析方法 | |
Edlund | A parallel, iterative method of moments and physical optics hybrid solver for arbitrary surfaces | |
Zhao et al. | Sparse matrix canonical grid method for three-dimension rough surface | |
Lv et al. | Higher order hierarchical Legendre basis functions application to the analysis of scattering by uniaxial anisotropic objects |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |