CN105184033B - 基于阶数步进金属目标的宽频带电磁特性快速预估方法 - Google Patents
基于阶数步进金属目标的宽频带电磁特性快速预估方法 Download PDFInfo
- Publication number
- CN105184033B CN105184033B CN201410277268.XA CN201410277268A CN105184033B CN 105184033 B CN105184033 B CN 105184033B CN 201410277268 A CN201410277268 A CN 201410277268A CN 105184033 B CN105184033 B CN 105184033B
- Authority
- CN
- China
- Prior art keywords
- triangular patch
- threshold value
- impedance matrix
- scattering
- triangle
- 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.)
- Active
Links
Landscapes
- Measuring Magnetic Variables (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于阶数步进金属目标的宽频带电磁特性快速预估方法,步骤如下:将导体目标表面用三角形面元进行剖分离散,从而得到金属目标的时域阶数步进电场积分方程以及磁场积分方程;确定每个三角形面元与其他三角形面元质心之间的距离,根据三角形面元剖分尺寸确定阈值,把任意两个三角形面元的质心距离与该阈值进行比较:当质心距离小于该阈值时,采用传统的矩量法计算阻抗矩阵元素;当质心距离大于阈值时,采用等效偶极子模型法计算阻抗矩阵元素;根据阻抗矩阵,采用迭代法求解出导体表面的散射电流;根据导体表面的散射电磁流,求解出雷达散射截面积。本发明能够加速矩阵的填充,大幅度减少了计算时间,具有很强的实际工程应用价值。
Description
技术领域
本发明属于目标电磁散射特性数值计算技术领域,特别是一种基于阶数步进金属目标的宽频带电磁特性快速预估方法。
背景技术
如何精确、髙效的分析目标的电磁散射问题一直作为计算电磁学的重要使命,被解决的方法也是多种多样。时域阶数步进矩量法是一种精确的数值求解方法,在分析宽频带的电磁辐射和散射问题中得到了广泛的应用。但是,对于求解目标是电大尺寸导体时候,采用传统矩量法进行数值求解,在计算机内存占用量和计算时间两方面都存在难以解决的困难。
发明内容
本发明的目的在于提供一种高效、稳定的基于阶数步进金属目标的宽频带电磁特性快速预估方法,以在计算阻抗矩阵时大幅度减少处理时间。
实现本发明目的的技术解决方案为:一种基于阶数步进金属目标的宽频带电磁特性快速预估方法,步骤如下:
步骤1、用时域平面波模拟脉冲照射导体飞行目标,将导体目标表面用三角形面元进行剖分离散,从而得到金属目标的时域阶数步进电场积分方程以及磁场积分方程;
步骤2、确定步骤1中离散所得的每个三角形面元与其他三角形面元质心之间的距离,根据三角形面元剖分尺寸确定阈值,把任意两个三角形面元的质心距离与该阈值进行比较:当质心距离小于该阈值时,采用传统的矩量法计算阻抗矩阵元素;当质心距离大于阈值时,采用等效偶极子模型法计算阻抗矩阵元素;
步骤3、根据步骤2所得到的阻抗矩阵,采用迭代法求解出导体表面的散射电流;
步骤4、根据导体表面的散射电流,求解出雷达散射截面积。
本发明与现有技术相比,其显著优点为:(1)可解决宽频带问题:时域积分方法可解决宽频带问题的求解;(2)矩阵填充时间大幅度减少:基于等效偶极子模型,求解场源相互作用时候减少计算量,大大缩减了填充矩阵的时间;(3)后时稳定性好:阶数步进时域积分方法是无条件稳定的,克服了时间步进积分方程得后时不稳定的弱点。
附图说明
图1是RWG基函数示意图。
图2是本发明等效偶极子模型。
图3是本发明场源三角形对相互作用时的等效偶极子模型。
图4是本发明实施例1中不同频点下的计算结果示意图,其中(a)频率为50MHz时的双站RCS,(b)频率为100MHz时的双站RCS,(c)频率为150MHz时的双站RCS,(d)频率为200MHz时的双站RCS,(e)频率为250MHz时的双站RCS,(f)频率为300MHz时的双站RCS。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细描述。
步骤1、用时域平面波模拟具有极宽的频谱范围的脉冲照射导体飞行目标,将导体目标表面用三角形面元进行剖分离散,从而得到金属目标的时域阶数步进电场积分方程以及磁场积分方程,具体步骤如下:
(1.1)使用商业软件(ANSYS),将导体目标表面使用一定数量的三角形面元进行剖分离散,对三角形、三角形的节点、三角形的内边分别进行编号,并且记录每个三角形的三个顶点坐标;
(1.2)根据金属表面切向连续的边界条件,建立时域阶数步进电场积分方程以及磁场积分方程:
其中,和分别为观察点处的入射电场和磁场,和分别为观察点处的散射电场和磁场,是金属表面处的单位外法向矢量,是散射目标上处t时刻的感应电流。
入射电磁场在导体目标上产生等效入射电磁流,该等效入射电磁流在目标上生成感应电磁场并产生相应的散射电磁流,由此可推导出等效入射电磁场的表达式:
其中,是散射目标上处τ时刻的感应电流, 为任意场点相对于坐标原点的位置矢量,为任意源点相对于坐标原点的位置矢量,ε0为媒质中的介电常数、μ0为媒质中的磁导率,τ=t-R1/c是滞后的时间,c是电磁波在媒质中的传播速度,表示表示为场点处的梯度算子,为源点处的梯度算子;
(1.3)将感应电流使用RWG空间基函数展开,RWG基函数如图1所示,任意一对基函数是由两个拥有公共边的三角形组成的,对于任意三角形基函数,和分别代表其公共边为m的上下三角形,和分别代表上、下三角形的面积,公共边的边长为lm,由此第m个基函数的定义为:
式中,代表和中任意一点,即为观察点;为基向量;和分别代表和的自由顶点,即与公共边相对的那个点;基函数被用来近似表示表面电流,它具有如下特点:边界边不存在完整的基函数三角形对,因此我们定义边界边上不存在线电流;每一个基函数的公共边处的外法向电流分量是常数且是连续的;对的面微分与面电量密度成正比。
由公式(5)、(6)可以将表面电流展开成如下形式:
其中,是待求的未知电流系数,需要用时间基函数展开;为RWG基函数。再将用时间基函数展开,时间基函数使用加权Laguerre多项式:
加权Laguerre多项的表达式为:
第0阶:L0(t)=e-0.5γt
第1阶:L1(t)=e-0.5γt(1-γt)
第i阶:Li(t)=[(2t-1-γt)Li-1(t)-(t-1)Li-2(t)]/i,i≥2
其中,γ为尺度因子且γ≥0;
对式(8)进行积分以及求导,都有解析表达式:
步骤2、确定步骤1中离散所得的每个三角形面元与其他三角形面元质心之间的距离,根据三角形面元剖分尺寸确定阈值,把任意两个三角形面元的质心距离与该阈值进行比较:当质心距离小于该阈值时,采用传统的矩量法计算阻抗矩阵元素;当质心距离大于阈值时,采用等效偶极子模型法计算阻抗矩阵元素;具体包括以下步骤:
(2.1)确定步骤1中离散所得的每个三角形面元与其他三角形面元质心之间的距离,根据三角形面元剖分尺寸确定阈值,该阈值为三角形面元剖分尺寸的2~4倍,把任意两个三角形面元的质心距离与该阈值进行比较:当质心距离小于该阈值时,采用(2.2)中传统的矩量法计算阻抗矩阵元素;当质心距离大于阈值时,采用(2.3)中的等效偶极子模型法计算阻抗矩阵元素;例如当剖分尺寸为0.1个波长的时候,阀值取0.3个波长得到的结果和全部用矩量法填充矩阵得到的结果相对比,就可以比较吻合。
(2.2)当质心距离小于该阈值时,采用传统的矩量法计算阻抗矩阵元素;
时域电场积分方程以及时域磁场积分方程分别表示为式(3)和式(4),将式(9)、(10)代入式(3)、(4)中得到电场和磁场积分方程最终形式:
对式(11)和(12)进行空间上的Galerkin测试,使用Ns个RWG空间基函数依次对式(11)和(12)进行测试,这样可以得到Ns个方程组,再采用(N0+1)个加权Laguerre多项式进行时间上的Galerkin测试,最终可以得到Ns×Ns的方程组,写成矩阵形式:
其中,表示求解第v阶系数Iv时候的激
励,表示进行空间和时间伽辽金测试后当前阶数产生的电场阻抗矩阵,表示求导
项进行空间和时间伽辽金测试后第v-ω阶产生的电场阻抗矩阵,表示积分项进行空
间和时间伽辽金测试后第v-ω阶产生的电场阻抗矩阵,表示进行空间和时间伽辽金测
试后当前阶数产生的磁场阻抗矩阵,表示求导项进行空间和时间伽辽金测试后第v-
ω阶产生的磁场阻抗矩阵,表示自由项进行空间和时间伽辽金测试后第v-ω阶产生
的磁场阻抗矩阵。
(2.3)当质心距离大于阈值时,采用等效偶极子模型法计算阻抗矩阵元素;
对于RWG基函数,如图2所示,当场源三角形单元质心之间距离大于阀值时候,对基函数的积分可以写成:
其中表示t时刻在r处的感应电流;是偶极子矩量,其值等于有效偶极子长度和有效偶极子电流的乘积,有效偶极子的长度为ln是基函数公共边的长度,和分别是上、下三角形自由顶点为起点到高斯积分点的向量,和分别是上、下三角形观察点到高斯积分点的向量。
结合图3,偶极子模型的等效入射电磁场表示为:
其中,代表空间处t时刻的散射电场值,代表空间处t时刻的散射磁场值, 代表空间任意一点到下三角形质心的距离,代表空间任意一点到上三角形质心的距离,ln代表三角形对公共边的边长;In(t)表示个加权Laguerre多项式,为只与时间基函数相关的电流展开项;η表示波阻抗;τ=t-R2/c是滞后的时间,c是电磁波在媒质中的传播速度;ε为媒质中的介电常数、μ为媒质中的磁导率;标量矢量为观察点到场三角形对公共边中点的向量,为观察点到源三角形对公共边中点的向量。
对等效入射电磁场进行时间以及空间上的伽辽金测试,得到电场积分方程和磁场积分方程的形式如下:
步骤3、根据步骤2所得到的阻抗矩阵,采用迭代法求解出导体表面的散射电流;由步骤2中计算得到的阻抗矩阵,采用迭代法计算出导体表面的感应电流,具体公式如下:
其中,代表散射目标上第v阶上的散射电流,表示进行空间和时间伽辽金测试后第v阶数上产生的阻抗矩阵,这里空间测试依然采用矩量法测试时候用的RWG基函数,时间测试函数依然采用加权Laguerre多项式,是对入射场进行空间和时间伽辽金测试后产生的向量,v、ω代表阶数步进的阶数标号。
步骤4、根据导体表面的散射电流,求解出雷达散射截面积,具体如下:
其中,为空间任意点处的散射电场的θ方向的分量,分别为空间任意点处的散射电场的方向的分量,μ0为自由空间的磁导率,k0为自由空间波数,π为圆周率,r为场点与源点的观察距离,为入射电场的θ方向的分量,为入射电场的方向的分量,为入射磁场的θ方向的分量,为入射磁场的方向的分量。
三维坐标系下,在(θ,φ)方向的双站RCS为:
其中,Es表示散射场的电场分量,Ei分别表示入射场的电场分量,
实施例1
本实例对半径为1m的金属球进行了电磁特性的仿真,仿真的频率为0-300MHz,球心位于(0,0,0),金属球以0.1λ剖分,λ为最高频率对应的自由空间波长,得到2568个三角形(3852条内边),临界值取0.3λ,拉盖尔多项式取80阶,为了验证程序的正确性以及效率,本专利结果与全部用矩量法填充阻抗矩阵的计算方法得到的结果做了比较。图4给出了不同频点处的双站RCS值,其中(a)频率为50MHz时的双站RCS,(b)频率为100MHz时的双站RCS,(c)频率为150MHz时的双站RCS,(d)频率为200MHz时的双站RCS,(e)频率为250MHz时的双站RCS,(f)频率为300MHz时的双站RCS,可以看出本发明的结果与矩量法得到的结果吻合的很好。
本发明利用等效电偶极子模型,快速填充阶数步进时域面积分方程的阻抗矩阵,提高了传统阶数步进时域面积分方程算法的效率,从而求解时域矩阵方程并提取宽频带内的电磁特性。本发明相较于传统的时域阶数步进方法,使用自由空间格林函数,具有很高的计算精度且未知量较小,可以一次计算一个频带内所有频点的电磁特性,能够加速矩阵的填充,在计算阻抗矩阵时大幅度减少了计算时间,并且得到的结果和矩量法填充阻抗矩阵得到的结果吻合的很好,采用时域积分方法可解决宽频带问题的求解,并且阶数步进时域积分方法是无条件稳定的,克服了时间步进积分方程得后时不稳定的弱点,后时稳定性好,具有很强的实际工程应用价值。
Claims (5)
1.一种基于阶数步进金属目标的宽频带电磁特性快速预估方法,其特征在于,步骤如下:
步骤1、用时域平面波模拟脉冲照射导体飞行目标,将导体目标表面用三角形面元进行剖分离散,从而得到金属目标的时域阶数步进电场积分方程以及磁场积分方程;
步骤2、确定步骤1中离散所得的每个三角形面元与其他三角形面元质心之间的距离,根据三角形面元剖分尺寸确定阈值,把任意两个三角形面元的质心距离与该阈值进行比较:当质心距离小于该阈值时,采用传统的矩量法计算阻抗矩阵元素;当质心距离大于阈值时,采用等效偶极子模型法计算阻抗矩阵元素;
步骤3、根据步骤2所得到的阻抗矩阵,采用迭代法求解出导体表面的散射电流;
步骤4、根据导体表面的散射电流,求解出雷达散射截面积。
2.根据权利要求1所述的基于阶数步进金属目标的宽频带电磁特性快速预估方法,其特征在于,步骤1所述将导体目标表面用三角形面元进行剖分离散,从而得到金属目标的时域阶数步进电场积分方程以及磁场积分方程,具体包括以下步骤:
(1.1)将导体目标表面使用三角形面元进行剖分离散,对三角形、三角形的节点、三角形的内边分别进行编号,并且记录每个三角形的三个顶点坐标;
(1.2)根据金属表面切向连续的边界条件,建立时域阶数步进电场积分方程以及磁场积分方程:
其中,和分别为观察点处的入射电场和磁场,和分别为观察点处的散射电场和磁场,是金属表面处的单位外法向矢量,是散射目标上处t时刻的感应电流。
3.根据权利要求1所述的基于阶数步进金属目标的宽频带电磁特性快速预估方法,其特征在于,步骤2中所述计算阻抗矩阵元素,具体包括以下步骤:
(2.1)确定步骤1中离散所得的每个三角形面元与其他三角形面元质心之间的距离,根据三角形面元剖分尺寸确定阈值,该阈值为三角形面元剖分尺寸的2~4倍,把任意两个三角形面元的质心距离与该阈值进行比较:当质心距离小于该阈值时,采用(2.2)中传统的矩量法计算阻抗矩阵元素;当质心距离大于阈值时,采用(2.3)中的等效偶极子模型法计算阻抗矩阵元素;
(2.2)当质心距离小于该阈值时,采用传统的矩量法计算阻抗矩阵元素;
时域电场积分方程以及时域磁场积分方程表示为:
其中,是散射目标上处τ时刻的感应电流, 为任意场点相对于坐标原点的位置矢量,为任意源点相对于坐标原点的位置矢量,ε0为媒质中的介电常数、μ0为媒质中的磁导率,τ=t-R1/c是滞后的时间,c是电磁波在媒质中的传播速度,表示 表示 代表场点处的梯度算子,代表源点处的梯度算子;和分别为观察点处的入射电场和磁场;是金属表面处的单位外法向矢量;
(2.3)当质心距离大于阈值时,采用等效偶极子模型法计算阻抗矩阵元素,其中偶极子模型的等效入射电磁场表示为:
其中,代表空间处t时刻的散射电场值,代表空间处t时刻的散射磁场值, 代表空间任意一点到下三角形质心的距离,代表空间任意一点到上三角形质心的距离,ln代表三角形对公共边的边长,η表示波阻抗,In(τ)代表只与时间基函数相关的电流展开项,τ=t-R2/c是滞后的时间,c是电磁波在媒质中的传播速度,标量矢量 为观察点到场三角形对公共边中点的向量,为观察点到源三角形对公共边中点的向量。
4.根据权利要求1所述的基于阶数步进金属目标的宽频带电磁特性快速预估方法,其特征在于,步骤3中所述采用迭代法求解出导体表面的散射电流,具体公式如下:
其中,代表散射目标上第v阶上的散射电流,表示进行空间和时间伽辽金测试后第v阶数上产生的阻抗矩阵,是对入射场进行空间和时间伽辽金测试后产生的向量,v、ω代表阶数步进的阶数标号。
5.根据权利要求1所述的基于阶数步进金属目标的宽频带电磁特性快速预估方法,其特征在于,步骤4所述雷达散射截面积的公式为:
其中,σ(θ,φ)为在(θ,φ)方向的双站RCS,Es表示散射场的电场分量,Ei分别表示入射场的电场分量,π为圆周率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410277268.XA CN105184033B (zh) | 2014-06-19 | 2014-06-19 | 基于阶数步进金属目标的宽频带电磁特性快速预估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410277268.XA CN105184033B (zh) | 2014-06-19 | 2014-06-19 | 基于阶数步进金属目标的宽频带电磁特性快速预估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105184033A CN105184033A (zh) | 2015-12-23 |
CN105184033B true CN105184033B (zh) | 2019-01-04 |
Family
ID=54906112
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410277268.XA Active CN105184033B (zh) | 2014-06-19 | 2014-06-19 | 基于阶数步进金属目标的宽频带电磁特性快速预估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105184033B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107305536B (zh) * | 2016-04-22 | 2020-11-13 | 南京理工大学 | 混合阶时域不连续伽略金方法 |
CN106529082A (zh) * | 2016-12-02 | 2017-03-22 | 上海无线电设备研究所 | 一种快速计算电大尺寸目标电磁散射特征的方法 |
CN108519622B (zh) * | 2018-04-25 | 2020-01-21 | 中国科学院电子学研究所 | 基于天然场源激励的地下电性目标探测方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101706586A (zh) * | 2009-12-03 | 2010-05-12 | 成都理工大学 | 一种时间域视电阻率的计算方法 |
CN101710187A (zh) * | 2009-12-17 | 2010-05-19 | 成都理工大学 | 一种时间域航空电磁高度校正方法 |
US20130317337A1 (en) * | 2012-05-23 | 2013-11-28 | Chenyu Wu | Continuous Modeling for Dipole Localization from 2D MCG Images with Unknown Depth |
-
2014
- 2014-06-19 CN CN201410277268.XA patent/CN105184033B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101706586A (zh) * | 2009-12-03 | 2010-05-12 | 成都理工大学 | 一种时间域视电阻率的计算方法 |
CN101710187A (zh) * | 2009-12-17 | 2010-05-19 | 成都理工大学 | 一种时间域航空电磁高度校正方法 |
US20130317337A1 (en) * | 2012-05-23 | 2013-11-28 | Chenyu Wu | Continuous Modeling for Dipole Localization from 2D MCG Images with Unknown Depth |
Non-Patent Citations (3)
Title |
---|
Analysis of transient electromagnetic scattering using UV method enhanced time‐domain integral equations with Laguerre polynomials;Wang Q, Shi Y, Li M, et al.;《Microwave and Optical Technology Letters》;20111231;第53卷(第1期);158-163 |
Marching‐on‐in‐order time‐domain volume‐surface integral equations in analysis of transient electromagnetic scattering from composite conducting–dielectric objects;Wang Q Q, Yan C, Shi Y F, et al.;《Microwave and Optical Technology Letters》;20111231;第53卷(第6期);1315-1319 |
基于Laguerre 多项式的边界积分方程时域求解;柴舜连、王生水、陈忠宽、毛钧杰;《电子与信息学报》;20080331;第30卷(第3期);742-745 |
Also Published As
Publication number | Publication date |
---|---|
CN105184033A (zh) | 2015-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103870654A (zh) | 基于并行矩量法与物理光学混合的电磁散射仿真方法 | |
CN105184033B (zh) | 基于阶数步进金属目标的宽频带电磁特性快速预估方法 | |
He et al. | A vector meshless parabolic equation method for three-dimensional electromagnetic scatterings | |
Boag et al. | A fast physical optics (FPO) algorithm for double-bounce scattering | |
Wei et al. | Four dimensional reconstruction using magnetic induction tomography: Experimental study | |
He et al. | A novel parallel parabolic equation method for electromagnetic scatterings | |
CN105630740B (zh) | 基于矩阵泰勒级数展开的电磁分析方法 | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
CN104915465A (zh) | 基于延迟拉盖尔多项式的金属目标瞬态电磁散射分析方法 | |
Voccola et al. | Polarimetric synthetic-aperture inversion for extended targets in clutter | |
CN105277927B (zh) | 飞行器编队瞬态电磁特性时域阶数步进分析方法 | |
CN104778286B (zh) | 掠海飞行器电磁散射特性快速仿真方法 | |
CN107305536A (zh) | 混合阶时域不连续伽略金方法 | |
McCullagh et al. | Nonlinear behavior of Baryon Acoustic Oscillations in redshift space from the Zel'dovich approximation | |
CN104778293B (zh) | 非均匀介质目标电磁散射的体积分Nystrom分析方法 | |
CN106503349B (zh) | 一种类周期结构目标电磁散射特性快速计算方法 | |
CN105182310B (zh) | 一种机动交会下目标近场角闪烁的统计标定方法 | |
CN106294283B (zh) | 基于泰勒级数展开的时域积分方程快速方法 | |
Wiid et al. | Towards electromagnetic characterization of MeerKAT Telescope | |
Ridler et al. | Evaluating and expressing uncertainty in high-frequency electromagnetic measurements: a selective review | |
Chen et al. | Model the electromagnetic scattering from three-dimensional PEC object buried under rough ground by MoM and modified PO hybrid method | |
Anuradha et al. | Discrimination of closely resembling PEC targets based on natural resonant frequencies | |
Kang et al. | Near-to-far-field transformation in FDTD for high resolution radar signal simulation | |
Xiang et al. | Application of barycentric subdivision method for singularity integration in method of moments | |
Xiang et al. | The application of Barycentric subdivision method for numerical integration in method of moments |
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 |