CN106294898B - 一种加速分析介质目标电磁散射特性的复点源求解方法 - Google Patents
一种加速分析介质目标电磁散射特性的复点源求解方法 Download PDFInfo
- Publication number
- CN106294898B CN106294898B CN201510272526.XA CN201510272526A CN106294898B CN 106294898 B CN106294898 B CN 106294898B CN 201510272526 A CN201510272526 A CN 201510272526A CN 106294898 B CN106294898 B CN 106294898B
- Authority
- CN
- China
- Prior art keywords
- equivalent
- point source
- complex point
- matrix
- group
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种加速分析介质目标电磁散射特性的复点源求解方法。该方法结合了复点源的思想对表面积分的远场区域进行求解。本方法结合八叉树分组的思想对介质目标先进行分组,基函数按组分好,每个组都由一个虚拟的复点源等效球面紧紧包住,近场组和自作用部分仍保留原来的方法加以计算,远场组之间的作用转化为用其包围的复点源来表示,即转化为复点源之间的相互作用。由于复点源具有方向性,可以舍去一些在不影响计算精度条件下的弱相互作用来加速计算远场组之间的作用。本发明在计算均匀介质目标的电磁散射特性分析中能够节省计算时间和降低内存的消耗。具有很强的实际工程应用价值。
Description
技术领域
本发明属于目标电磁散射特性数值计算技术领域,具体是一种基于复点源方法的雷达散射截面快速计算方法。
背景技术
目标的雷达回波特性在军事中具有很重要的意义,提出一种精确而有效的电磁分析模型显得极为重要。如何快速地分析出物体的雷达散射截面,一直是广大学者们致力研究的方向。随着计算机技术的快速发展,与之而产生的电磁场与微波问题的数值分析方法也应运而生。当前计算电磁学也面临着诸多挑战——分析目标越来越复杂精细、目标电尺寸越来越大、材料越来越新颖等。数值计算方法可分为两大类:一类为微分方程类方法,包括有限元法(FEM)、时域有限差分法(FDTD)等。另一类为积分方程类方法,如矩量法(MoM)等。基于积分方程的矩量法具有严格的理论模型,电磁场在无限远处的辐射条件已解析的包含在积分方程类方法中,所以矩量法能够精确地模拟电磁波传播的索末菲辐射条件,因此被广泛运用于各种天线辐射、复杂散射体散射以及静态或准静态等领域。但是矩量法生成的阻抗矩阵是满阵,阻抗矩阵存储量的量级是矩阵阶数的平方O(N2),直接求解的计算复杂度是O(N3),迭代求解的计算复杂度是O(N2)。
最近由Koray Tap等人提出一种复点源方法,该方法利用复点源对矩阵元素进行压缩,先构造虚拟等效球面,并将其包围的基函数用复点源进行展开,互为远场的基函数之间的作用就可以利用相对应复点源之间的作用进行展开和表示。但是该方法仅仅考虑了分析金属目标的情况,对于介质目标则没有考虑,因此计算使用范围不够广泛。而且此方法中,仅仅给出了考虑电场积分的情况,没有考虑磁场积分算子存在的情况,对此提出了扩展和应用。
发明内容
本发明的目的在于提供一种加速分析介质目标电磁散射特性的复点源求解方法。
实现本发明目的的技术解决方案为:一种加速分析介质目标电磁散射特性的复点源求解方法,步骤如下:
步骤1、建立介质目标的离散模型,设置入射波频率为freq,对散射物体进行三角形网格剖分如图1,利用八叉树分组对三角形的边进行分组,得到每个组内的边数和边的编号;包含边信息的组为非空组,不包含边信息的组为空组;
步骤2、对模型剖分得到的三角形离散网格进行基函数的构造,即对边构造RWG基函数用于展开在入射波照射下目标表面产生的感应电流和磁流无论介质内外,总场等于入射场和散射场之和。利用电场和磁场在边界处的切向分量连续特性,建立电场积分方程和磁场积分方程,并加以组合得到表面积分方法的矩阵方程表达形式;
步骤4、将分块矩阵ZEJ,ZHJ,ZEM,ZHM中的远场作用组形成的矩阵用等效面上复点源之间的作用表示,设定角度阈值为θ,当组中心与等效面上等效点形成的矢量与组中心之间的方向矢量之间的夹角大于θ时,舍去该等效点,反之则保留。对矩阵进行组合得到最终的矩阵方程组。
步骤5、利用迭代方法求解得到的矩阵方程组,结合求出的感应电流和感应磁流系数,计算雷达散射截面RCS。
本发明与现有技术相比,其显著优点为:(1)将复点源技术推广到分析介质问题,扩宽了其运用背景;(2)由于复点源具有方向性,可以舍去一部分复点源,从而节省内存和计算时间;(3)与现有的表面积分方法相比,本发明方法具有较低的计算复杂度,内存和时间消耗都大大减少,计算效率有了明显提高。
附图(表)说明
图1是本发明计算目标网格剖分示意图。
图2是本发明基函数示意图。
图3是本发明等效面上复点源的选取及分布示意图。
图4是本发明等效面上复点源舍取示意图。
图5是本发明实施例中某散射体双站RCS曲线图。
表1是本发明与其它方法时间及内存上的比较
具体实施方式
下面结合附图对本发明作进一步详细描述。
步骤1、建立介质目标的离散模型,设置入射波频率为freq,对物体进行三角形网格剖分如图1,利用八叉树分组对三角形的边进行分组,得到每个组内的边数和边的编号;包含边信息的组为非空组,不包含边信息的组为空组;
步骤2、根据权利要求1所述的复点源方法加速分析电磁散射特性的求解方法,其特征在于,所述步骤2中:
在介质外表面和内表面建立电场积分和磁场积分:
其中,j=1,2表示介质目标外表面和内表面,表示的是场点的位置,表示的是源点的位置。表示的是场点到源点之间的距离。是自由空间内的标量格林函数。分别表示外区域和内区域中的波数,P.V.是主值积分项。对(1)-(4)式进行组合,得到矩阵方程表达形式:
其中如图2所示,l表示两三角形公共边的边长,和表示该条边所在的两个三角形和的面积,代表和中任意的一点观察点,基向量和 和分别代表上三角形和下三角形中公共边对应的自由顶点,即除公共边两端点的另外两点。每个基函数都只定义在其对应的两个三角形上,在其他位置其值为零。
将感应电流和感应磁流用RWG基函数进行三维展开,并用伽辽金测试方法对(1)-(4)式进行测试,得到各分块矩阵方程组表达式。其中m和n分别表示第m和n个基函数的编号,L和K分别表示电场和磁场积分,和分别是自由空间和介质中的波阻抗。
2).在等效面上对基函数进行展开,具体步骤如下:
a).计算等效面测试矩阵ZEJ,ZHJ,ZEM,ZHM,对其进行并矢格林函数展开。
d).当第m个和第n个基函数所在的组互为远场组时,结合求得的复点源系数矩阵和如图4所示的复点源角度阈值的取舍,对系数矩阵进行取舍,得到矩阵表达式如下,其中Q表示取舍后的复点源点数,上标1和2分别表示介质表面内和外的情况。
最终就得到了矩阵表达形式中各分块矩阵的具体表达形式,I为待求的电磁流展开系数,V为右边向量。矩阵方程如下:
实施例
本实施例进行了电磁散射的典型仿真,仿真在主频2.8GHz、内存8GB的个人计算机上实现,以介质长方体模型为例,如图4所示,入射波频率为300MHz,入射波的方向θ=0°,长方体沿着x,y,z方向分别为0.8米、0.8米和4米,最细层采用0.2波长分组,共分为3层,第1层、第2层和第三层的非空组尺寸分别为0.25米、0.5米、1米。介电常数εr=2,电导率μr=1。对应的复点源波束宽度b分别为0.12、0.25、0.52,截断角度阈值分别为150°、135°、120°。图5为本发明方法和表面积分方法的RCS曲线图对比结果,从图中的曲线可以看出,本文方法与正确的数值结果吻合,另外与表面积分方法进行时间及内存上的对比(见表1)说明本专利方法能够快速仿真分析介质目标物体的电磁散射特性。
本发明方法 | 表面积分方法 | |
内存(MB) | 745 | 5410 |
计算时间(秒) | 258 | 838 |
表1。
Claims (5)
1.一种加速分析介质目标电磁散射特性的复点源求解方法,其特征在于步骤如下:
步骤1、建立物体的离散模型,设入射波频率为freq,对介质物体进行三角形网格剖分,利用八叉树分组对三角形的边进行分组,得到每个组内的边数和边的编号;包含边信息的组为非空组,不包含边信息的组为空组;
步骤2、对模型剖分得到的三角形离散网格进行基函数的构造,即对边构造RWG基函数用于展开在入射波照射下目标表面产生的感应电流和磁流利用电场和磁场在边界处的切向分量连续特性,建立电场积分方程和磁场积分方程,并加以组合得到表面积分方法的矩阵方程表达形式;
步骤2中,
在介质外表面和内表面建立电场积分和磁场积分:
其中,j=1,2表示介质目标外表面和内表面,表示的是场点的位置,表示的是源点的位置,表示的是场点到源点之间的距离,是自由空间内的标量格林函数,分别表示外区域和内区域中的波数,P.V.是主值积分项;对(1)-(4)式进行组合,得到矩阵方程表达形式:
步骤4、将分块矩阵ZEJ,ZHJ,ZEM,ZHM中的远场作用组形成的矩阵用等效面上复点源之间的作用表示,设定角度阈值为θ,当组中心与等效面上等效点形成的矢量与组中心之间的方向矢量之间的夹角大于θ时,舍去该等效点,反之则保留;对矩阵进行组合得到最终的矩阵方程组;
步骤5、利用迭代方法求解得到的矩阵方程组,结合求出的感应电流和感应磁流系数,计算雷达散射截面RCS。
其中,l表示两三角形公共边的边长,和表示该条边所在的两个三角形和的面积,代表和中任意的一点观察点,基向量和 和分别代表上三角形和下三角形中公共边对应的自由顶点,即除公共边两端点的另外两点;每个基函数都只定义在其对应的两个三角形上,在其他位置其值为零;
将感应电流和感应磁流用RWG基函数进行三维展开,并用伽辽金测试方法对(1)-(4)式进行测试,得到各分块矩阵方程组表达式;其中m和n分别表示第m和n个基函数的编号,L和K分别表示电场和磁场积分,和分别是自由空间和介质中的波阻抗;
4.根据权利要求1所述的加速分析介质目标电磁散射特性的复点源求解方法,其特征在于:步骤4中所述具体步骤为:
2).在等效面上对基函数进行展开,具体步骤如下:
a).计算等效面测试矩阵ZEJ,ZHJ,ZEM,ZHM,对其进行并矢格林函数展开;
d).当第m个和第n个基函数所在的组互为远场组时,结合求得的复点源系数矩阵和复点源角度阈值的取舍,对系数矩阵进行取舍,得到矩阵表达式如下,其中Q表示取舍后的复点源点数,上标1和2分别表示介质表面内和外的情况;
最终得到矩阵表达形式中各分块矩阵的具体表达形式,I为待求的电磁流展开系数,V为右边向量,矩阵方程如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510272526.XA CN106294898B (zh) | 2015-05-25 | 2015-05-25 | 一种加速分析介质目标电磁散射特性的复点源求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510272526.XA CN106294898B (zh) | 2015-05-25 | 2015-05-25 | 一种加速分析介质目标电磁散射特性的复点源求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106294898A CN106294898A (zh) | 2017-01-04 |
CN106294898B true CN106294898B (zh) | 2020-04-10 |
Family
ID=57634522
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510272526.XA Active CN106294898B (zh) | 2015-05-25 | 2015-05-25 | 一种加速分析介质目标电磁散射特性的复点源求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106294898B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110929375B (zh) * | 2019-10-17 | 2021-08-31 | 中国科学院电子学研究所 | 基于二维矩量法和射线追迹法的透镜高效仿真、优化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902841A (zh) * | 2012-08-28 | 2013-01-30 | 同济大学 | 一种计算机载复杂结构雷达横截面积的方法 |
CN103425816A (zh) * | 2013-06-09 | 2013-12-04 | 南京理工大学 | 快速获取金属旋转对称体电磁散射特性的矩阵抽取方法 |
CN103487789A (zh) * | 2013-09-23 | 2014-01-01 | 北京航空航天大学 | 应用于复杂目标多次散射的雷达散射截面快速获取方法 |
CN104200074A (zh) * | 2014-08-18 | 2014-12-10 | 南京理工大学 | 快速获取目标电磁散射特性的多层复波束方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7916067B2 (en) * | 2009-02-11 | 2011-03-29 | The Boeing Company | Removing clutter from radar cross section measurements using spectral tagging |
-
2015
- 2015-05-25 CN CN201510272526.XA patent/CN106294898B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902841A (zh) * | 2012-08-28 | 2013-01-30 | 同济大学 | 一种计算机载复杂结构雷达横截面积的方法 |
CN103425816A (zh) * | 2013-06-09 | 2013-12-04 | 南京理工大学 | 快速获取金属旋转对称体电磁散射特性的矩阵抽取方法 |
CN103487789A (zh) * | 2013-09-23 | 2014-01-01 | 北京航空航天大学 | 应用于复杂目标多次散射的雷达散射截面快速获取方法 |
CN104200074A (zh) * | 2014-08-18 | 2014-12-10 | 南京理工大学 | 快速获取目标电磁散射特性的多层复波束方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106294898A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105930567A (zh) | 一种基于子区域自适应积分的电磁散射特性获取方法 | |
CN104200074B (zh) | 快速获取目标电磁散射特性的多层复波束方法 | |
You et al. | A new method to analyze broadband antenna-radome interactions in time-domain | |
US10302684B2 (en) | Method for determining electric field and related applications | |
CN106294894B (zh) | 快速分析非均匀目标电磁散射特性的有限元边界积分方法 | |
CN112733364B (zh) | 一种基于阻抗矩阵分块的箔条云散射快速计算方法 | |
Li et al. | Fast periodic interpolation method for periodic unit cell problems | |
CN108021533A (zh) | 一种基于广义坐标系求解任意色散材料电磁特性的方法 | |
CN110362877A (zh) | 不确定因素的环境电磁散射特性分析方法 | |
Sevilla et al. | The use of hybrid meshes to improve the efficiency of a discontinuous Galerkin method for the solution of Maxwell’s equations | |
CN111899329A (zh) | 全各向异性媒质的仿真方法、装置、计算机设备及存储介质 | |
CN111144013B (zh) | 高精度介质体目标散射的仿真方法 | |
Ayari et al. | Method of moments versus advanced transverse wave approach for EM validation of complex microwave and RF applications | |
CN113033052B (zh) | 一种蜂窝吸波结构电磁快速数值建模方法 | |
CN104915465A (zh) | 基于延迟拉盖尔多项式的金属目标瞬态电磁散射分析方法 | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
CN106294898B (zh) | 一种加速分析介质目标电磁散射特性的复点源求解方法 | |
Hu et al. | Analysis of scattering from composite conducting dispersive dielectric objects by time-domain volume-surface integral equation | |
CN111931353A (zh) | 一种应用于仿真fss结构的散射场求解方法 | |
Na et al. | A Multiple Huygens Surface Based Ray Tracing Framework with GPU Acceleration | |
CN105205299B (zh) | 电大目标电磁散射特性快速降维分析方法 | |
Huang et al. | Study on Improved CG Algorithm for Reradiation Interference Calculation of Transmission Steel Tower | |
Chobanyan et al. | Volume integral equation-based diakoptic method for electromagnetic modeling | |
Oikawa et al. | Analysis of lightning electromagnetic field on large‐scale terrain model using three‐dimensional MW‐FDTD parallel computation | |
Bilbao et al. | Numerical analysis of electromagnetic fields |
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 |