CN106294920B - 分析介质目标电磁散射特性的时域准显式方法 - Google Patents

分析介质目标电磁散射特性的时域准显式方法 Download PDF

Info

Publication number
CN106294920B
CN106294920B CN201510326023.6A CN201510326023A CN106294920B CN 106294920 B CN106294920 B CN 106294920B CN 201510326023 A CN201510326023 A CN 201510326023A CN 106294920 B CN106294920 B CN 106294920B
Authority
CN
China
Prior art keywords
time
medium
equation
quasi
transient
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
Application number
CN201510326023.6A
Other languages
English (en)
Other versions
CN106294920A (zh
Inventor
陈如山
丁大志
樊振宏
呼延龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Tech University
Original Assignee
Nanjing Tech University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Tech University filed Critical Nanjing Tech University
Priority to CN201510326023.6A priority Critical patent/CN106294920B/zh
Publication of CN106294920A publication Critical patent/CN106294920A/zh
Application granted granted Critical
Publication of CN106294920B publication Critical patent/CN106294920B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Resistance Or Impedance (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种分析介质目标电磁散射特性的时域准显式方法。步骤如下：建立时域体积分方程；将待求未知量进行时间上用时间基函数进行离散，空间上用非共形的四面体网格进行离散；形成待求解的矩阵方程，未知电流为介质瞬态体电流；利用准显式方法高效求解矩阵方程，得到介质的瞬态体电流系数，再根据体电流系数确定瞬态电磁散射参量。本发明中利用准显式方法求解的非共形时域体积分方程方法可以灵活的处理待求目标的网格离散，形成的相对稠密的迭代矩阵可以利用准显式方法的优势加速求解。

Description

分析介质目标电磁散射特性的时域准显式方法
一技术领域
本发明属于电磁仿真技术领域，特别是一种分析介质目标电磁散射特性的时域准显式方法。
二背景技术
介质目标的电磁散射特性分析是电磁问题中的非常重要研究领域，雷达目标的形状和体积等物理量都是通过对雷达散射截面等参数进行计算得出的，且雷达散射截面积是雷达系统对目标“可观测性”的一个重要指标。因此，对于各种目标散射特性的研究在这些应用领域具有特别重要的现实意义。
近年，随着宽频带电磁散射系统的快速发展，瞬态电磁散射特性的分析越来越引起科研学者和工程人员的关注。相比于其它方法，时域体积分方程方法非常适合于介质目标瞬态电磁散射特性的分析,尤其适合非均匀介质目标瞬态电磁散射特性的分析(NoelT.Gres,Arif A.Ergin and Eric Michielssen,“Volume-integral-equation-basedanalysis of transient electromagnetic scattering from three-dimensionalinhomogeneous dielectric objects,”Radio Science,vol.36,no.3,pp.379–386,2001.)。但是当分析的介质目标存在高度不均匀介电常数，或者是存在多尺度的时候，网格的处理成为了普通时域体积分方程面临的难题，虽然可以利用非共形网格的时域积分方程方法进行分析，但是由于迭代矩阵相对稠密，求解效率较低。
三发明内容
本发明的目的在于提供一种更加高效和准确地分析非均匀或存在多尺度的介质目标电磁散射特性的时域准显式方法。
实现本发明目的的技术解决方案为：一种分析介质目标电磁散射特性的时域准显式方法，步骤如下：
步骤1，建立介质时域体积分方程；
步骤2，将待求未知量采用四阶拉格朗日基函数进行时间上的离散，并采用非共形的四面体单元进行空间上的离散；
步骤3，形成准显式方法待求解的矩阵方程，未知电流为介质瞬态体电流；
步骤4，利用预测-校正方法求解矩阵方程，得到介质的瞬态体电流系数，再根据体电流系数确定瞬态电磁散射参量。
本发明与现有技术相比，其显著优点是：(1)可以更加灵活和准确地分析非均匀或存在多尺度的介质目标瞬态电磁散射特性，并且对离散网格具有鲁棒性；(2)可以高效地求解矩阵方程。
四附图说明
图1是本发明中算例模型示意图。
图2是本发明实施例中介质目标在不同频率点处的双站雷达散射截面结果图。
图3是本发明实施例中介质目标上某一点处的电流随时间的变化曲线图。
五具体实施方式
下面结合附图对本发明作进一步详细描述。
本发明分析介质目标电磁散射特性的时域准显式方法，步骤如下：
步骤1，建立介质时域体积分方程；
令电磁波照射到介质目标上，在介质体内产生感应体电流J，根据介质的电场边界条件，即总电场等于入射电场与散射电场之和，得到介质时域体积分方程TD-VIE，如下：
Einc(r,t)+Esca(r,t)＝Etot(r,t) (1)
其中，Einc表示照射在介质目标上的电磁波的入射电场，Etot表示总电场，Esca表示介质目标在电磁波照射后产生的散射电场，t为当前时刻，瞬态散射电场的表达形式为：
Figure BDA0000736980290000021
将式(2)代入式(1)，则式(1)重新改写为：
Figure BDA0000736980290000022
其中，V表示四面体单元，μ0表示自由空间的磁导率，ε表示介电常数，ε0表示自由空间的介电常数，εr为介质体的相对介电常数，r为场的位置坐标，r′为源的位置坐标，R为场源点之间的距离，c表示真空中的光速，
Figure BDA0000736980290000023
表示对时间函数的积分，
Figure BDA0000736980290000024
表示对时间函数的求导，▽为梯度算子。
步骤2，将待求未知量采用四阶拉格朗日基函数进行时间上的离散，并采用非共形的四面体单元进行空间上的离散，具体如下：
介质目标的瞬态感应体电流可离散表示如下：
Figure BDA0000736980290000031
其中：
Figure BDA0000736980290000032
式中，fn(r)为半个SWG基函数，Tl(t)为四阶拉格朗日时间基函数，
Figure BDA0000736980290000035
为第n个未知量在第l时刻的待求瞬态电流系数，n为空间未知量编号、l为时间步数编号，NV为空间未知量个数、Nt为时间步数。
步骤3，形成准显式方法待求解的矩阵方程，未知电流为介质瞬态体电流，具体如下：
为了得到一阶常微分方程的形式，将式(3)在时间上进行两次求导：
Figure BDA0000736980290000033
将式(6)空间上采用伽辽金测试，时间上采用点测试，得准显式的时域体积分的矩阵方程形式：
其中
Figure BDA0000736980290000041
式中，
Figure BDA0000736980290000044
其中，
Figure BDA0000736980290000045
为建立场源之间联系的时域阻抗矩阵，G为稀疏的Gram矩阵，表示积分主值项上场源基函数作用形成的矩阵，表示第i个时间步的激励，Δt表示每个时间步长，Ii和Ij分别是第i和第j个时间步的待求未知量的系数，i、j为时间步的编号，fv(r)为半个SWG基函数，v表示介质，
Figure BDA0000736980290000047
和
Figure BDA0000736980290000048
分别表示第m个和第n个空间基函数，S为积分区域，和
Figure BDA00007369802900000410
分别表示第m个和第n个介质三角形的外法向量，gj(iΔt,R)为时域格林函数、Tj为时间基函数。
步骤4，利用预测-校正方法求解矩阵方程，得到介质的瞬态体电流系数，再根据体电流系数确定瞬态电磁散射参量，采用预测-修正型线性k步法，具体如下：
在假设式(7)矩阵方程中的Ii(i＝0:j)已知的情况下可以利用式(7)来计算
Figure BDA00007369802900000412
的值。为了得到Ij的值，应该对式(7)进行时间积分。为此，要用到预测-修正型线性k步法，并且用逐次超松弛(SOR)来提高预测-修正的精确性和稳定性。生成的时间步进如下：
假定在第j时间步，j＝k,…,Nt-1，设置初始值：
Figure BDA0000736980290000051
(4.1)计算式(7)中等号右边在一个时间步长内不会发生变化的固定部分
Figure BDA0000736980290000052
Figure BDA0000736980290000053
(4.2)通过预测系数p、Ij与的过去值预测Ij：
Figure BDA0000736980290000055
上式中，l、k表示时刻的变量；
(4.3)将预测的Ij值及式(7)右边的固定部分
Figure BDA0000736980290000056
代入式(7)计算
Figure BDA0000736980290000057
Figure BDA0000736980290000058
其中，
Figure BDA0000736980290000059
表示第j时刻电流导数的系数在预测后的初始值，
Figure BDA00007369802900000510
表示场源基函数距离在一个时间步内的基函数对形成的矩阵；
(4.4)重复操作以下步骤，直到结果Ij收敛，具体如下：
(a)用修正系数c修正Ij：
Figure BDA00007369802900000511
(b)将修正的
Figure BDA00007369802900000512
值与式(7)右边的固定部分
Figure BDA00007369802900000513
代入式(7)计算
Figure BDA00007369802900000514
Figure BDA00007369802900000515
其中，表示第j时刻电流导数的系数在第q次修正后的值；
(c)对应用逐次超松弛SOR：
(d)检查收敛性，如果结果收敛，即χPECE为预设的修正收 敛精度，那么
Figure BDA0000736980290000061
在步骤(c)中，α∈[0,1]是SOR参数，选择α＝1就是不使用SOR；选择α＜0.5将对
Figure BDA0000736980290000062
更依赖于前面时刻的系数；选择α＞0.5将对
Figure BDA0000736980290000063
更依赖于当前时刻的系数。
通过上述方法得到介质的瞬态体电流系数，再根据体电流系数确定瞬态电磁散射参量。
实施例1
为了验证本发明方法的准确性与高效性，下面给出了一个介质球壳的介质目标的瞬态电磁特性的分析，其中，球壳的内半径为0.25m，厚度为0.05m,相对介电常数为2，如图1。瞬态电磁散射的双站RCS的计算结果与解析值Mie计算的结果相比较吻合较好，图2所示，并验证了在(-0.239,-0.091,0.067)位置处的电流稳定性，如图3。
本算例中，入射电场采用调制高斯平面波，其表达式如下：
Figure BDA0000736980290000064
其中，极化方向
Figure BDA0000736980290000065
传播方向tc＝7σ,
Figure BDA0000736980290000067
Einc(r,t)的频谱的中心频率为f0＝50MHz，最高频率为100MHz，fbw为频带宽 度。时间步长Δt＝0.3lm，总时间步Nt＝300，lm是光米(light meter)，即光在自由空间中 传播1m距离所花的时间。表1给出了求解时间的数据对比，体现出本方法的高效性。
表1计算效率对比
求解时间(分钟)
隐式方法 146.87
准显式方法 28.88
综上所述，本发明与传统的时域体积分方程方法相比，基于非共形网格方法求解的时域体积分方程方法可以更加灵活的处理待求目标的网格离散，特别是对于不均匀的介质体或者是存在多尺度模型的情况。并针对稠密矩阵迭代较慢的情况，引入了准显式的求解方法，加速了系统的求解速度。

Claims (5)

1.一种分析介质目标电磁散射特性的时域准显式方法，其特征在于，步骤如下：
步骤1，建立介质时域体积分方程；
步骤2，将待求未知量采用四阶拉格朗日基函数进行时间上的离散，并采用非共形的四面体单元进行空间上的离散；
步骤3，形成准显式方法待求解的矩阵方程，未知电流为介质瞬态体电流；
步骤4，利用预测-校正方法求解矩阵方程，得到介质的瞬态体电流系数，再根据体电流系数确定瞬态电磁散射参量。
2.根据权利要求1所述的分析介质目标电磁散射特性的时域准显式方法，其特征在于，步骤1中所述建立介质时域体积分方程，具体如下：
令电磁波照射到介质目标上，在介质体内产生感应体电流J，根据介质的电场边界条件，即总电场等于入射电场与散射电场之和，得到介质时域体积分方程TD-VIE，如下：
Einc(r,t)+Esca(r,t)＝Etot(r,t) (1)
其中，Einc表示照射在介质目标上的电磁波的入射电场，Etot表示总电场，Esca表示介质目标在电磁波照射后产生的散射电场，t为当前时刻，瞬态散射电场的表达形式为：
Figure FDA0002276159400000011
将式(2)代入式(1)，则式(1)重新改写为：
Figure FDA0002276159400000012
其中，μ0表示自由空间的磁导率，ε表示介电常数，ε0表示自由空间的介电常数，εr为介质体的相对介电常数，r为场的位置坐标，r′为源的位置坐标，R为场源点之间的距离，c表示真空中的光速，
Figure FDA0002276159400000013
表示对时间函数的积分，表示对时间函数的求导，
Figure FDA0002276159400000015
为梯度算子。
3.根据权利要求1所述的分析介质目标电磁散射特性的时域准显式方法，其特征在于，步骤2中所述将待求未知量采用四阶拉格朗日基函数进行时间上的离散，并采用非共形的四面体单元进行空间上的离散，具体如下：
介质目标的瞬态感应体电流离散表示如下：
Figure FDA0002276159400000016
其中：
Figure FDA0002276159400000021
式中，fn(r)为半个SWG基函数，Tl(t)为四阶拉格朗日时间基函数，为第n个未知量在第l时刻的待求瞬态电流系数，n为空间未知量编号、l为时间步数编号，NV为空间未知量个数、Nt为时间步数。
4.根据权利要求1所述的分析介质目标电磁散射特性的时域准显式方法，其特征在于，步骤3中所述形成准显式方法待求解的矩阵方程，未知电流为介质瞬态体电流，具体如下：
为了得到一阶常微分方程的形式，将式(3)在时间上进行两次求导：
将式(6)空间上采用伽辽金测试，时间上采用点测试，得准显式的时域体积分的矩阵方程形式：
Figure FDA0002276159400000024
式中，
Figure FDA0002276159400000025
Figure FDA0002276159400000026
Figure FDA0002276159400000031
式中，
其中，
Figure FDA0002276159400000033
为建立场源之间联系的时域阻抗矩阵，G为稀疏的Gram矩阵，Vi exp表示第i个时间步的激励，Δt表示每个时间步长，Ii和Ij分别是第i和第j个时间步的待求未知量的系数，i、j为时间步的编号，fv(r)为半个SWG基函数，v表示介质，
Figure FDA0002276159400000034
和
Figure FDA0002276159400000035
分别表示第m个和第n个空间基函数，S为积分区域，
Figure FDA0002276159400000036
和
Figure FDA0002276159400000037
分别表示第m个和第n个介质三角形的外法向量，gj(iΔt,R)为时域格林函数、Tj为时间基函数。
5.根据权利要求1所述的分析介质目标电磁散射特性的时域准显式方法，其特征在于，步骤4中所述利用预测-校正方法求解矩阵方程，得到介质的瞬态体电流系数，再根据体电流系数确定瞬态电磁散射参量，采用预测-修正型线性k步法，具体如下：
假定在第j时间步，j＝k,…,Nt-1，设置初始值：
Figure FDA0002276159400000038
(4.1)计算式(7)中等号右边在一个时间步长内不会发生变化的固定部分
Figure FDA0002276159400000039
Figure FDA00022761594000000310
(4.2)通过预测系数p、Ij与
Figure FDA00022761594000000311
的过去值预测Ij：
上式中，l、k表示时刻的变量；
(4.3)将预测的Ij值及式(7)右边的固定部分
Figure FDA00022761594000000313
代入式(7)计算
Figure FDA00022761594000000314
Figure FDA0002276159400000041
其中，表示第j时刻电流导数的系数在预测后的初始值，
Figure FDA0002276159400000043
表示场源基函数距离在一个时间步内的基函数对形成的矩阵；
(4.4)重复操作以下步骤，直到结果Ij收敛，具体如下：
(a)用修正系数c修正Ij：
Figure FDA0002276159400000044
(b)将修正的值与式(7)右边的固定部分代入式(7)计算
Figure FDA0002276159400000047
Figure FDA0002276159400000048
其中，
Figure FDA0002276159400000049
表示第j时刻电流导数的系数在第q次修正后的值；
(c)对
Figure FDA00022761594000000410
应用逐次超松弛SOR：
Figure FDA00022761594000000411
(d)检查收敛性，如果结果
Figure FDA00022761594000000412
收敛，即
Figure FDA00022761594000000413
χPECE为预设的修正收敛精度，那么
Figure FDA00022761594000000414
在步骤(c)中，α∈[0,1]是SOR参数，选择α＝1就是不使用SOR；选择α＜0.5将对
Figure FDA00022761594000000415
更依赖于前面时刻的系数；选择α＞0.5将对
Figure FDA00022761594000000416
更依赖于当前时刻的系数。
CN201510326023.6A 2015-06-12 2015-06-12 分析介质目标电磁散射特性的时域准显式方法 Active CN106294920B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510326023.6A CN106294920B (zh) 2015-06-12 2015-06-12 分析介质目标电磁散射特性的时域准显式方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510326023.6A CN106294920B (zh) 2015-06-12 2015-06-12 分析介质目标电磁散射特性的时域准显式方法

Publications (2)

Publication Number Publication Date
CN106294920A CN106294920A (zh) 2017-01-04
CN106294920B true CN106294920B (zh) 2020-02-11

Family

ID=57650201

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510326023.6A Active CN106294920B (zh) 2015-06-12 2015-06-12 分析介质目标电磁散射特性的时域准显式方法

Country Status (1)

Country Link
CN (1) CN106294920B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115062465A (zh) * 2022-06-13 2022-09-16 南京理工大学 非均匀介质旋转对称体的多维不确定性问题快速分析方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103279601A (zh) * 2013-05-17 2013-09-04 南京理工大学 导体目标宽带电磁散射特性的仿真方法
CN104699870A (zh) * 2013-12-06 2015-06-10 南京理工大学 电大复杂有耗介质目标电磁散射抛物线快速仿真方法
CN104778286A (zh) * 2014-01-09 2015-07-15 南京理工大学 掠海飞行器电磁散射特性快速仿真方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050234686A1 (en) * 2004-02-06 2005-10-20 Hongwei Cheng Analysis method and system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103279601A (zh) * 2013-05-17 2013-09-04 南京理工大学 导体目标宽带电磁散射特性的仿真方法
CN104699870A (zh) * 2013-12-06 2015-06-10 南京理工大学 电大复杂有耗介质目标电磁散射抛物线快速仿真方法
CN104778286A (zh) * 2014-01-09 2015-07-15 南京理工大学 掠海飞行器电磁散射特性快速仿真方法

Also Published As

Publication number Publication date
CN106294920A (zh) 2017-01-04

Similar Documents

Publication Publication Date Title
CN110276109B (zh) 一种高超声速飞行器等离子体鞘套电磁特性的仿真方法
CN107167807B (zh) 三维天线阵列综合孔径辐射计分段图像反演方法
CN113567943B (zh) 基于saim与cat获取载体平台宽带rcs的方法
CN116306890A (zh) 一种反演蜂窝等效介电常数的神经网络的训练方法
CN111931353A (zh) 一种应用于仿真fss结构的散射场求解方法
CN106294920B (zh) 分析介质目标电磁散射特性的时域准显式方法
CN108090311A (zh) 基于路径特征集与支持向量回归的空气击穿电压预测方法
Huang et al. A parallel FDTD/ADI-PE method for ultralarge-scale propagation modeling of ILS signal analysis
CN104778293B (zh) 非均匀介质目标电磁散射的体积分Nystrom分析方法
CN104915326A (zh) 基于等效原理的区域分解阶数步进时域积分方法
CN105630740B (zh) 基于矩阵泰勒级数展开的电磁分析方法
CN105184033B (zh) 基于阶数步进金属目标的宽频带电磁特性快速预估方法
CN106202595B (zh) 分析介质目标瞬态电磁散射特性的时域非共形网格方法
Gan et al. Improved T-matrix method for simultaneous reconstruction of dielectric and perfectly conducting scatterers
CN109239771A (zh) 一种基于非均匀背景介质的弹性波成像方法
Xiong et al. Low frequency modeling for electromagnetic analysis of arbitrary anechoic chambers
CN105095546A (zh) 分析多尺度导体目标电磁散射特性的混合阶Nystrom方法
CN109657196B (zh) 一种基于空域强耦合路径加速的mom-ka混合方法
CN114491992A (zh) 一种基于等效偶极矩和物理光学法的高效电磁散射方法
CN112016225A (zh) 面积分自洽耦合分析方法、装置、设备及存储介质
Johansson et al. Comparison between two phase-retrieval methods for electromagnetic source modeling
Essiben et al. Electromagnetic compatibility of aperture antennas using electromagnetic band gap structure
CN105095154A (zh) 分析超高速飞行目标电磁散射的高阶体面积分方程方法
Du et al. Three-dimensional completion method with uniformity decision mechanism for spectrum sensing data
RodriguezVarela et al. Multi-level spherical wave expansion for fast near-field to far-field transformation

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