CN110274920B - 体面剖分弹跳射线法分析金属介质目标瞬态特性的方法 - Google Patents

体面剖分弹跳射线法分析金属介质目标瞬态特性的方法 Download PDF

Info

Publication number
CN110274920B
CN110274920B CN201810208704.6A CN201810208704A CN110274920B CN 110274920 B CN110274920 B CN 110274920B CN 201810208704 A CN201810208704 A CN 201810208704A CN 110274920 B CN110274920 B CN 110274920B
Authority
CN
China
Prior art keywords
ray
medium
vector
target
tetrahedron
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
CN201810208704.6A
Other languages
English (en)
Other versions
CN110274920A (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 University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201810208704.6A priority Critical patent/CN110274920B/zh
Publication of CN110274920A publication Critical patent/CN110274920A/zh
Application granted granted Critical
Publication of CN110274920B publication Critical patent/CN110274920B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N22/00—Investigating or analysing materials by the use of microwaves or radio waves, i.e. electromagnetic waves with a wavelength of one millimetre or more

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Image Generation (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种体面剖分弹跳射线法分析金属介质目标瞬态特性的方法。具体步骤如下：建立金属介质混合目标的几何模型，使用四面体单元对介质体进行体剖分，使用三角形单元对金属表面进行面剖分，用来拟合目标的几何外形；为了达到弹跳射线法的理想精度，对照亮区的表面三角形单元进行细分；用射线管模拟电磁波在介质体内和金属表面的传播，追踪射线的场强和路径；当射线离开目标时，计算每根射线的瞬态散射场并累加得到总的瞬态散射信号。本发明采用体面剖分的弹跳射线法可以用于计算亚毫米波段电大尺寸金属介质混合目标瞬态散射信号，避免了在亚毫米波段未知量大的问题，所需时间少，一次计算就可以获得宽频带的RCS。

Description

体面剖分弹跳射线法分析金属介质目标瞬态特性的方法
技术领域
本发明属于微波测量技术领域，特别是一种体面剖分弹跳射线法分析金属介质目标瞬态特性的方法。
背景技术
在我国国防科学研究领域里，导弹系统的设计与仿真，雷达系统的设计与鉴定以及隐身技术的研究等，均需要建立目标(如飞机、坦克等)与环境的电磁散射模型。在亚毫米波段，以坦克目标而言，工作波长远小于坦克的几何尺寸，在这种波段上的电磁散射为高频散射，往往采用弹跳射线法，而以往的弹跳射线法只能用于计算金属目标，而周围的环境及目标覆盖物往往为介质，并且忽略了目标上的一些介质结构对散射场的影响，从而导致计算结果跟真实结果误差较大。
在对这种金属介质混合目标进行电磁特性分析时，数值方法如时域有限差分法(FDTD)，体面积分方程法(VSIE)等，由于受计算时间和内存需求的限制，往往并不适合电大尺寸的情况。传统的SBR算法一般应用于计算金属或者薄涂覆目标的电磁特性，电磁波只在目标表面发生弹跳现象，当涂覆比较厚时，就不再适用。
发明内容
本发明的目的在于提供一种体面剖分弹跳射线法分析金属介质目标瞬态特性的方法，通过采用体面剖分的弹跳射线法，对电大尺寸的目标与周边环境进行高效分析。
实现本发明目的的技术解决方案为：一种体面剖分弹跳射线法分析金属介质目标瞬态特性的方法，步骤如下：
步骤1，建立金属介质混合目标的几何模型，使用四面体单元对介质体进行体剖分，使用三角形单元对金属表面进行面剖分，剖分单元大小要能够拟合物体的几何外形；
步骤2，对照亮区的表面三角形单元进行细分，细分后小三角形的边长为波长的1/4～1/6；
步骤3、将细分后的小三角形看作是一系列的射线管，追踪射线的场强和路径；
步骤4、当射线离开目标时，根据时域积分方程求出每根射线的时域散射场，并进行累加得到总的瞬态散射信号，从而得到目标的宽频带RCS。
进一步地，步骤3所述将细分后的小三角形看作是一系列的射线管，追踪射线的场强和路径，具体包括：
(1)四面体内的路径追踪：因为每个四面体内的介质为均匀介质，射线管在均匀介质中沿着直线传播，所以射线在四面体中的路径为一条射线；
(2)当射线照射介质分界面时，求出射线在介质分界面上的反射系数和透射系数，并求出反射方向和透射方向，然后继续对反射射线和透射射线进行路径追踪和场强追踪；
(3)当射线照射到金属表面时，利用Snell定律追踪反射射线的射线路径和时域场强信号的变化。
进一步地，步骤(1)所述四面体内的路径追踪，具体如下：
因为每个四面体内的介质为均匀介质，射线管在均匀介质中沿着直线传播，所以射线在四面体中的路径为一条射线；
当射线穿入一个四面体时，设该四面体四个顶点为
Figure BDA0001596592260000021
射线在四面体内的入射点为
Figure BDA0001596592260000022
射线传播的方向矢量为
Figure BDA0001596592260000023
则射线所在的直线方程为：
Figure BDA0001596592260000024
要求射线穿出四面体时的交点，首先要求出射线从四面体穿出的所在面，三个候选面分别为
Figure BDA0001596592260000025
对应的法向量分别为
Figure BDA0001596592260000026
由此得：
Figure BDA0001596592260000027
t表示从
Figure BDA0001596592260000028
出发，沿着射线传播的方向矢量
Figure BDA0001596592260000029
到候选面所走过的距离，
Figure BDA00015965922600000210
为射线传播方向的单位矢量；当t为负数时，表示沿着射线传播的反方向所走过的距离；
取t1,t2,t3中最小的正值t0所对应的面即为射线从四面体穿出的面，同时得到穿出点的坐标为：
Figure BDA00015965922600000211
进一步地，步骤(2)中所述求出反射方向和透射方向，具体如下：
一束相位矢量为
Figure BDA00015965922600000212
幅度矢量为
Figure BDA00015965922600000213
的非均匀电磁波从介质1进入到介质2，设分界面的法向量为
Figure BDA0001596592260000031
和
Figure BDA0001596592260000032
之间的夹角为θ1，
Figure BDA0001596592260000033
和
Figure BDA0001596592260000034
之间的夹角为
Figure BDA0001596592260000035
Figure BDA0001596592260000036
与
Figure BDA0001596592260000037
之间的夹角为ρ1，同理，
Figure BDA0001596592260000038
和
Figure BDA0001596592260000039
分别表示在介质2中透射波的幅度矢量和相位矢量，
Figure BDA00015965922600000310
和
Figure BDA00015965922600000311
之间的夹角为θ2，
Figure BDA00015965922600000312
和
Figure BDA00015965922600000313
之间的夹角为
Figure BDA00015965922600000314
Figure BDA00015965922600000315
与
Figure BDA00015965922600000316
之间的夹角为ρ2；
θ2和
Figure BDA00015965922600000317
表示为
Figure BDA00015965922600000318
Figure BDA00015965922600000319
式中
Figure BDA00015965922600000320
Figure BDA00015965922600000321
Figure BDA00015965922600000322
Figure BDA00015965922600000323
式中，介质i中的本征传播常量，只和介质的本身特性有关，表示为：
Figure BDA00015965922600000324
Figure BDA00015965922600000325
εri,μri为介质i的相对复介电常数和磁导率，定义γ1t为：
Figure BDA00015965922600000326
从而得到反射方向和透射方向。
进一步地，步骤(2)中所述求出射线在介质分界面上的反射系数和透射系数具体为：
当电磁波传播到两种介质的交界面处时，电场分解为水平和垂直极化分量，此时反射电场
Figure BDA0001596592260000041
为：
Figure BDA0001596592260000042
其中，
Figure BDA0001596592260000043
是入射电场，
Figure BDA0001596592260000044
是垂直极化方向单位矢量，
Figure BDA0001596592260000045
是水平极化方向单位矢量；R⊥和R//分别是垂直极化和水平极化反射系数；
透射电场
Figure BDA0001596592260000046
为：
Figure BDA0001596592260000047
其中，T⊥和T//分别表示垂直极化和水平极化透射系数；
垂直极化的反射系数和透射系数表示为：
Figure BDA0001596592260000048
平行极化波的反射系数和透射系数为：
Figure BDA0001596592260000049
式中，ε1,μ1和ε2,μ2分别为介质1中和介质2中的相对介电常数和磁导率。
进一步地，步骤(3)中所述当射线照射到金属表面时，利用Snell定律追踪反射射线的射线路径和时域场强信号的变化，具体如下：
在介质和空气界面上计算等效的表面电磁流，根据Huygen原理，等效电流J和等效磁流M表示为：
Figure BDA00015965922600000410
Figure BDA00015965922600000411
当射线照射到介质表面时，不仅会存在感应电流，而且会存在感应磁流，所以射线管的远场散射场
Figure BDA00015965922600000412
表示为：
Figure BDA0001596592260000051
式中，
Figure BDA0001596592260000052
和
Figure BDA0001596592260000053
分别表示感应电流和感应磁流：
Figure BDA0001596592260000054
Figure BDA0001596592260000055
其中，
Figure BDA0001596592260000056
为观察点位置矢量，
Figure BDA0001596592260000057
为目标表面任一面元的位置矢量，r为面元到观察点的距离，
Figure BDA0001596592260000058
为目标表面面元的单位法矢量，
Figure BDA0001596592260000059
分别为物体表面总电场和总磁场，
Figure BDA00015965922600000510
为时间延迟，
Figure BDA00015965922600000511
为射线第一次与目标相交的位置矢量，
Figure BDA00015965922600000512
为入射方向单位矢量，
Figure BDA00015965922600000513
为观察方向单位矢量，
Figure BDA00015965922600000514
表示在介质中传播时所引起的幅度的衰减，
Figure BDA00015965922600000515
表示在介质中传播时所引起的相位的变化，n表示在传播过程中穿行四面体的个数，li表示在第i个四面体中穿行的距离，αi表示在媒质i的电磁波传播的幅度矢量的模，βi表示在媒质i的电磁波传播的相位矢量的模，ρi中表示在媒质i的电磁波传播的幅度矢量和相位矢量的夹角，*表示卷积，F-1表示逆傅里叶变换。
本发明与现有技术相比，其显著优点为：(1)引出了时域弹跳射线方法对不均匀有耗介质的处理方法；(2)由于介质体是由四面体单元剖分得到，金属是由三角形单元剖分得到，故可以拟合任意形状的金属介质混合目标，同时，对介质体内各个四面体单元设置不同的介电参数，可以拟合非均匀介质的介质体目标，在建模上得到了极大的简化；(3)射线在介质体内传播过程中，利用射线在四面体内的求交技术，加快了射线传播的追踪过程，极大的提高了计算速度；(4)对目标照亮区的三角形单元进行细分，从而大大减少了计算所需的内存。
附图说明
图1是本发明体面剖分弹跳射线法分析金属介质目标瞬态特性的方法的体剖分示意图。
图2是本发明体面剖分弹跳射线法分析金属介质目标瞬态特性的方法的四面体内路径追踪示意图。
图3是本发明体面剖分弹跳射线法分析金属介质目标瞬态特性的方法中有耗介质中电磁波的传播方向和角度示意图。
图4是本发明体面剖分弹跳射线法分析金属介质目标瞬态特性的方法中电场强度追踪示意图。
图5是本发明体面剖分弹跳射线法分析金属介质目标瞬态特性的方法中照亮区三角形单元细分示意图。
图6是本发明体面剖分弹跳射线法分析金属介质目标瞬态特性的方法中实施例1模型图。
图7是本发明体面剖分弹跳射线法分析金属介质目标瞬态特性的方法中实施例1结果图。
图8是本发明体面剖分弹跳射线法分析金属介质目标瞬态特性的方法中实施例1结果图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细描述。
本发明体面剖分弹跳射线法分析金属介质目标瞬态特性的方法，步骤如下：
步骤1，建立金属介质混合目标的几何模型，使用四面体单元对介质体进行体剖分，使用三角形单元对金属表面进行面剖分，剖分单元大小要能够拟合物体的几何外形；
步骤2，为了达到弹跳射线法的精度，对照亮区的表面三角形单元进行细分，细分后小三角形的边长为波长的1/4～1/6；
步骤3、将细分后的小三角形看作是一系列的射线管，追踪射线的场强和路径。当射线照射介质分界面时，由于两侧电磁参数的差异，电磁波会在表面发生反射和透射现象，反射和透射射线分别携带一部分能量继续传播，因此有必要求出射线在介质分界面上的反射系数和透射系数以及反射方向和透射方向，然后继续对反射射线和透射射线进行路径追踪和场强追踪。当射线照射到金属表面时，利用Snell定律追踪反射射线的射线路径和时域场强信号的变化；
步骤4、当射线离开目标时，根据时域积分方程求出每根射线的时域散射场，并进行累加得到总的瞬态散射信号，从而可以得到目标的宽频带RCS。
对本发明上述步骤进一步详细描述如下：
一、目标建模
对于金属介质混合目标，对介质部分采用四面体单元剖分建模，对四面体单元映射相应的介电参数，对金属部分采用表面三角形单元剖分建模。这样，便能拟合任意金属介质体目标的形状和描述介质体目标介电参数的非均匀特性。图1为体剖分示意图。
二、射线的路径追踪
由于本方法在建模时用的了四面体来剖分介质体，用三角面元剖分金属表面，以三角面元的中心点来模拟射线管，射线在介质体内射线传播的路径相当于在一个个四面体内穿行的过程。
步骤3所述将细分后的小三角形看作是一系列的射线管，追踪射线的场强和路径，具体包括：
(1)四面体内的路径追踪：因为每个四面体内的介质为均匀介质，射线管在均匀介质中沿着直线传播，所以射线在四面体中的路径为一条射线；
如图2所示，因为每个四面体内的介质为均匀介质，由几何光学原理可知，射线管在均匀介质中沿着直线传播，所以射线在四面体中的路径为一条射线。
当射线穿入一个四面体时，设该四面体四个顶点为
Figure BDA0001596592260000071
射线在四面体内的入射点为
Figure BDA0001596592260000072
射线传播的方向矢量为
Figure BDA0001596592260000073
则射线所在的直线方程：
Figure BDA0001596592260000074
要求得射线穿出四面体时的交点，首先要求得射线从四面体穿出的所在面，三个候选面分别为
Figure BDA0001596592260000075
它们的法向量分别为
Figure BDA0001596592260000076
由此可得
Figure BDA0001596592260000077
t表示从
Figure BDA0001596592260000078
出发，沿着射线传播的方向矢量
Figure BDA0001596592260000079
到候选面所走过的距离，
Figure BDA00015965922600000710
为射线传播方向的单位矢量；当t为负数时，表示沿着射线传播的反方向所走过的距离。
所以取t1,t2,t3中最小的正值t0所对应的面即为射线从四面体穿出的面，同时可得穿出点的坐标为：
Figure BDA00015965922600000711
绝大部分情况下，射线都会从四面体的某个面穿出，但是也会存在从四面体的棱边上穿出的情况，这时的处理就会比较复杂，首先，需要根据t1,t2,t3得到射线穿出的棱边，即t1,t2,t3中最小的两个正值所对应的面的公共棱边，然后根据几何关系，找出所有包含该棱边的四面体，然后把这些四面体的面当作候选面，按照上述方法找出穿出的面，同时也可以得到穿出点的坐标。图2为四面体内路径追踪示意图。
(2)当射线照射介质分界面时，求出射线在介质分界面上的反射系数和透射系数，并求出反射方向和透射方向，然后继续对反射射线和透射射线进行路径追踪和场强追踪；
①反射方向和透射方向的求解
如图3所示，一束相位矢量为
Figure BDA0001596592260000081
幅度矢量为
Figure BDA0001596592260000082
的非均匀电磁波从介质1进入到介质2，设分界面的法向量为
Figure BDA0001596592260000083
Figure BDA0001596592260000084
和
Figure BDA0001596592260000085
之间的夹角为θ1，
Figure BDA0001596592260000086
和
Figure BDA0001596592260000087
之间的夹角为
Figure BDA0001596592260000088
Figure BDA0001596592260000089
与
Figure BDA00015965922600000810
之间的夹角为ρ1，同理，
Figure BDA00015965922600000811
和
Figure BDA00015965922600000812
分别表示在介质2中透射波的幅度矢量和相位矢量，
Figure BDA00015965922600000813
和
Figure BDA00015965922600000814
之间的夹角为θ2，
Figure BDA00015965922600000815
和
Figure BDA00015965922600000816
之间的夹角为
Figure BDA00015965922600000817
Figure BDA00015965922600000818
与
Figure BDA00015965922600000819
之间的夹角为ρ2。
θ2和
Figure BDA00015965922600000820
表示为
Figure BDA00015965922600000821
Figure BDA00015965922600000822
式中
Figure BDA00015965922600000823
Figure BDA00015965922600000824
Figure BDA00015965922600000825
Figure BDA00015965922600000826
式中，介质i中的本征传播常量，只和介质的本身特性有关，表示为：
Figure BDA0001596592260000091
Figure BDA0001596592260000092
εri,μri为介质i的相对复介电常数和磁导率，定义γ1t为：
Figure BDA0001596592260000093
从而得到反射方向和透射方向。图3为有耗介质中电磁波的传播方向和角度示意图。
②射线的场强追踪
当电磁波传播到两种介质的交界面处时，电场分解为水平和垂直极化分量，此时反射电场
Figure BDA0001596592260000094
为：
Figure BDA0001596592260000095
其中，
Figure BDA0001596592260000096
是入射电场，
Figure BDA0001596592260000097
是垂直极化方向单位矢量，
Figure BDA0001596592260000098
是水平极化方向单位矢量；R⊥和R//分别是垂直极化和水平极化反射系数。
透射电场
Figure BDA0001596592260000099
为：
Figure BDA00015965922600000910
其中，T⊥和T//分别表示垂直极化和水平极化透射系数。
垂直极化的反射系数和透射系数可以表示为：
Figure BDA00015965922600000911
平行极化波的反射系数和透射系数为：
Figure BDA00015965922600000912
式中，ε1,μ1和ε2,μ2分别为介质1中和介质2中的相对介电常数和磁导率。图4为电场强度追踪示意图。
(3)当射线照射到金属表面时，利用Snell定律追踪反射射线的射线路径和时域场强信号的变化。
体面剖分的弹跳射线法凭借射线路径追踪和场强追踪，模拟平面波传播的平行射线从目标表面进入散射目标，根据几何光学原理，在不同介质的分边界面上计算反射和透射射线的传播矢量和场强，最后，在介质和空气界面上计算等效的表面电磁流。根据Huygen原理，等效电流J和等效磁流M可以表示为：
Figure BDA0001596592260000101
Figure BDA0001596592260000102
当射线照射到介质表面时，其表面上不仅会存在感应电流，而且会存在感应磁流。所以，射线管的远场散射场可以表示为：
Figure BDA0001596592260000103
式中，
Figure BDA0001596592260000104
和
Figure BDA0001596592260000105
分别表示感应电流和感应磁流对散射场的贡献，可分别表示为：
Figure BDA0001596592260000106
Figure BDA0001596592260000107
其中，
Figure BDA0001596592260000108
为观察点位置矢量，
Figure BDA0001596592260000109
为目标表面任一面元的位置矢量，r为面元到观察点的距离，
Figure BDA00015965922600001010
为目标表面面元的单位法矢量，
Figure BDA00015965922600001011
分别为物体表面总电场和总磁场，
Figure BDA00015965922600001012
为时间延迟，
Figure BDA00015965922600001013
为射线第一次与目标相交的位置矢量，
Figure BDA00015965922600001014
为入射方向单位矢量，
Figure BDA00015965922600001015
为观察方向单位矢量，
Figure BDA00015965922600001016
表示在介质中传播时所引起的幅度的衰减，
Figure BDA00015965922600001017
表示在介质中传播时所引起的相位的变化，n表示在传播过程中穿行四面体的个数，li表示在第i个四面体中穿行的距离，αi表示在媒质i的电磁波传播的幅度矢量的模，βi表示在媒质i的电磁波传播的相位矢量的模，ρi中表示在媒质i的电磁波传播的幅度矢量和相位矢量的夹角，*表示卷积，F-1表示逆傅里叶变换。
三、区三角形单元细分
体面剖分的弹跳射线法需要剖分的三角形边长大约为波长的1/4～1/6可达到理想精度，这大大提高了对计算机内存的要求，为了减小计算所需的内存，只需要对目标进行粗剖分，用来拟合目标的几何外形，然后对照明区的三角形单元进行细分，细分后的小三角形单元大约为波长的1/4～1/6，然后将小三角形当作初始射线管进行时域场强信号追踪和路径追踪，从而大大减少了计算所需的内存。
对于任意形状的大尺寸三角面元，可以取各边中点的连线逐级划分为四个相似的小三角面元，直到小三角形单元大约为波长的1/4～1/6，如图5所示。
实施例1
图6给出了一个立方体模型，内部金属立方体的尺寸为8λ*6λ*5λ，外部立方体的尺寸为10λ*8λ*6λ，中间部分填充为介质，介电参数为2-j0.2，射脉冲的中心频率为f0＝305GHz，所关心的有效频带为300～310GHz，采样频率为40GHz，脉冲的入射方向沿θ＝45°，Φ＝45°，θθ极化，后向观察。图7为金属介质混合立方体模型的时域响应，图8金属介质混合立方体模型的宽频带RCS。
体面剖分的TDSBR计算时间仅为36s，所需内存仅为5M。
本发明提出的体面剖分SBR方法还考虑了电磁波在介质中的传播和弹跳，不仅适用于薄涂覆金属目标，还适用于处理厚涂覆金属目标以及金属和任意介质混合目标的电磁散射问题，并且一次计算就可以获得宽频带内的RCS。

Claims (5)

1.一种体面剖分弹跳射线法分析金属介质目标瞬态特性的方法，其特征在于，步骤如下：
步骤1，建立金属介质混合目标的几何模型，使用四面体单元对介质体进行体剖分，使用三角形单元对金属表面进行面剖分，剖分单元大小要能够拟合物体的几何外形；
步骤2，对照亮区的表面三角形单元进行细分，细分后小三角形的边长为波长的1/4～1/6；
步骤3、将细分后的小三角形看作是一系列的射线管，追踪射线的场强和路径；
步骤4、当射线离开目标时，根据时域积分方程求出每根射线的时域散射场，并进行累加得到总的瞬态散射信号，从而得到目标的宽频带RCS；
步骤3所述将细分后的小三角形看作是一系列的射线管，追踪射线的场强和路径，具体包括：
(1)四面体内的路径追踪：因为每个四面体内的介质为均匀介质，射线管在均匀介质中沿着直线传播，所以射线在四面体中的路径为一条射线；
(2)当射线照射介质分界面时，求出射线在介质分界面上的反射系数和透射系数，并求出反射方向和透射方向，然后继续对反射射线和透射射线进行路径追踪和场强追踪；
(3)当射线照射到金属表面时，利用Snell定律追踪反射射线的射线路径和时域场强信号的变化。
2.根据权利要求1所述的体面剖分弹跳射线法分析金属介质目标瞬态特性的方法，其特征在于，步骤(1)所述四面体内的路径追踪，具体如下：
因为每个四面体内的介质为均匀介质，射线管在均匀介质中沿着直线传播，所以射线在四面体中的路径为一条射线；
当射线穿入一个四面体时，设该四面体四个顶点为
Figure FDA0003273587950000011
射线在四面体内的入射点为
Figure FDA0003273587950000012
射线传播的方向矢量为
Figure FDA0003273587950000013
则射线所在的直线方程为：
Figure FDA0003273587950000014
要求射线穿出四面体时的交点，首先要求出射线从四面体穿出的所在面，三个候选面分别为
Figure FDA0003273587950000015
对应的法向量分别为
Figure FDA0003273587950000016
由此得：
Figure FDA0003273587950000021
t表示从
Figure FDA0003273587950000022
出发，沿着射线传播的方向矢量
Figure FDA0003273587950000023
到候选面所走过的距离，
Figure FDA0003273587950000024
为射线传播方向的单位矢量；当t为负数时，表示沿着射线传播的反方向所走过的距离；
取t1,t2,t3中最小的正值t0所对应的面即为射线从四面体穿出的面，同时得到穿出点的坐标为：
Figure FDA0003273587950000025
3.根据权利要求1所述的体面剖分弹跳射线法分析金属介质目标瞬态特性的方法，其特征在于，步骤(2)中所述求出反射方向和透射方向，具体如下：
一束相位矢量为
Figure FDA0003273587950000026
幅度矢量为
Figure FDA0003273587950000027
的非均匀电磁波从介质1进入到介质2，设分界面的法向量为
Figure FDA0003273587950000028
Figure FDA0003273587950000029
和
Figure FDA00032735879500000210
之间的夹角为θ1，
Figure FDA00032735879500000211
和
Figure FDA00032735879500000212
之间的夹角为
Figure FDA00032735879500000213
Figure FDA00032735879500000214
与
Figure FDA00032735879500000215
之间的夹角为ρ1，同理，
Figure FDA00032735879500000216
和
Figure FDA00032735879500000217
分别表示在介质2中透射波的幅度矢量和相位矢量，
Figure FDA00032735879500000218
和
Figure FDA00032735879500000219
之间的夹角为θ2，
Figure FDA00032735879500000220
和
Figure FDA00032735879500000221
之间的夹角为
Figure FDA00032735879500000222
Figure FDA00032735879500000223
与
Figure FDA00032735879500000224
之间的夹角为ρ2；
θ2和
Figure FDA00032735879500000225
表示为
Figure FDA00032735879500000226
Figure FDA00032735879500000227
式中
Figure FDA00032735879500000228
Figure FDA00032735879500000229
Figure FDA00032735879500000230
Figure FDA0003273587950000031
式中，介质i中的本征传播常量，只和介质的本身特性有关，表示为：
Figure FDA0003273587950000032
Figure FDA0003273587950000033
εri,μri为介质i的相对复介电常数和磁导率，定义γ1t为：
Figure FDA0003273587950000034
从而得到反射方向和透射方向。
4.根据权利要求1所述的体面剖分弹跳射线法分析金属介质目标瞬态特性的方法，其特征在于，步骤(2)中所述求出射线在介质分界面上的反射系数和透射系数具体为：
当电磁波传播到两种介质的交界面处时，电场分解为水平和垂直极化分量，此时反射电场
Figure FDA0003273587950000035
为：
Figure FDA0003273587950000036
其中，
Figure FDA0003273587950000037
是入射电场，
Figure FDA0003273587950000038
是垂直极化方向单位矢量，
Figure FDA0003273587950000039
是水平极化方向单位矢量；R⊥和R//分别是垂直极化和水平极化反射系数；
透射电场
Figure FDA00032735879500000310
为：
Figure FDA00032735879500000311
其中，
Figure FDA00032735879500000313
和T//分别表示垂直极化和水平极化透射系数；
垂直极化的反射系数和透射系数表示为：
Figure FDA00032735879500000312
平行极化波的反射系数和透射系数为：
Figure FDA0003273587950000041
式中，ε1,μ1和ε2,μ2分别为介质1中和介质2中的相对介电常数和磁导率。
5.根据权利要求1所述的体面剖分弹跳射线法分析金属介质目标瞬态特性的方法，其特征在于，步骤(3)中所述当射线照射到金属表面时，利用Snell定律追踪反射射线的射线路径和时域场强信号的变化，具体如下：
在介质和空气界面上计算等效的表面电磁流，根据Huygen原理，等效电流J和等效磁流M表示为：
Figure FDA0003273587950000042
Figure FDA0003273587950000043
当射线照射到介质表面时，不仅会存在感应电流，而且会存在感应磁流，所以射线管的远场散射场
Figure FDA0003273587950000044
表示为：
Figure FDA0003273587950000045
式中，
Figure FDA0003273587950000046
和
Figure FDA0003273587950000047
分别表示感应电流和感应磁流：
Figure FDA0003273587950000048
Figure FDA0003273587950000049
其中，
Figure FDA00032735879500000410
为观察点位置矢量，
Figure FDA00032735879500000411
为目标表面任一面元的位置矢量，r为面元到观察点的距离，
Figure FDA00032735879500000412
为目标表面面元的单位法矢量，
Figure FDA00032735879500000413
分别为物体表面总电场和总磁场，
Figure FDA00032735879500000414
为时间延迟，
Figure FDA00032735879500000415
为射线第一次与目标相交的位置矢量，
Figure FDA00032735879500000416
为入射方向单位矢量，
Figure FDA00032735879500000417
为观察方向单位矢量，
Figure FDA00032735879500000418
表示在介质中传播时所引起的幅度的衰减，
Figure FDA00032735879500000419
表示在介质中传播时所引起的相位的变化，n表示在传播过程中穿行四面体的个数，li表示在第i个四面体中穿行的距离，αi表示在媒质i的电磁波传播的幅度矢量的模，βi表示在媒质i的电磁波传播的相位矢量的模，ρi中表示在媒质i的电磁波传播的幅度矢量和相位矢量的夹角，*表示卷积，F-1表示逆傅里叶变换。
CN201810208704.6A 2018-03-14 2018-03-14 体面剖分弹跳射线法分析金属介质目标瞬态特性的方法 Active CN110274920B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810208704.6A CN110274920B (zh) 2018-03-14 2018-03-14 体面剖分弹跳射线法分析金属介质目标瞬态特性的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810208704.6A CN110274920B (zh) 2018-03-14 2018-03-14 体面剖分弹跳射线法分析金属介质目标瞬态特性的方法

Publications (2)

Publication Number Publication Date
CN110274920A CN110274920A (zh) 2019-09-24
CN110274920B true CN110274920B (zh) 2021-11-23

Family

ID=67958334

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810208704.6A Active CN110274920B (zh) 2018-03-14 2018-03-14 体面剖分弹跳射线法分析金属介质目标瞬态特性的方法

Country Status (1)

Country Link
CN (1) CN110274920B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110907901B (zh) * 2019-10-29 2023-03-31 西安电子科技大学 一种基于tdsbr的海战无源干扰计算及评估方法
CN111965619B (zh) * 2020-08-25 2024-03-08 西安电子科技大学 一种基于射线追踪的复杂目标时域散射信号仿真方法
CN113092884A (zh) * 2021-04-20 2021-07-09 北京航空航天大学 基于弹跳射线法的天线远场预测方法、存储介质和装置
CN113255190B (zh) * 2021-06-07 2023-06-20 西安电子科技大学 一种不规则非均匀介质体雷达散射截面的计算方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106556833A (zh) * 2016-11-24 2017-04-05 上海无线电设备研究所 基于时域弹跳射线法快速近场计算的isar成像仿真方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7880671B2 (en) * 2008-09-18 2011-02-01 Raytheon Company Electromagnetic (EM) solver using a shooting bouncing ray (SBR) technique

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106556833A (zh) * 2016-11-24 2017-04-05 上海无线电设备研究所 基于时域弹跳射线法快速近场计算的isar成像仿真方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
New Results for the Effective Propagation Constants of Nonuniform Plane Waves at the Planar Interface of Two Lossy Media;Jasmin E 等;《IEEE Transactions on Antennas and Propagation》;20030630;第51卷(第6期);第1206-1215页 *
基于射线追踪算法的电波传播特性预测;钱驹;《中国优秀硕士学位论文全文数据库 信息科技辑》;20160531(第5期);I135-10 *
电大尺寸目标电磁散射特性的时域分析;陈磊;《中国优秀硕士学位论文全文数据库 信息科技辑》;20100715(第7期);I135-4 *

Also Published As

Publication number Publication date
CN110274920A (zh) 2019-09-24

Similar Documents

Publication Publication Date Title
CN110274920B (zh) 体面剖分弹跳射线法分析金属介质目标瞬态特性的方法
CN109100692B (zh) 基于迭代物理光学的粗糙面与多个目标复合散射仿真方法
Zaporozhets et al. Bistatic RCS calculations with the vector parabolic equation method
Sukharevsky et al. Modeling of ultrawideband (UWB) impulse scattering by aerial and subsurface resonant objects based on integral equation solving
Maheshwari et al. Application of emission source microscopy technique to EMI source localization above 5 GHz
Borzov et al. Mathematical modeling and simulation of the input signals of short-range radar systems
Cong et al. Numerical modeling of EM scattering from plasma sheath: A review
Zalevsky et al. Calculation of scattering characteristics of aerial radar objects of resonant sizes based on iterative algorithm
CN112986943B (zh) 一种计算蜂窝复合材料目标电磁散射的方法
Perotoni et al. Numerical evaluation of an air-to-air missile radar cross section signature at X-band
CN115169170A (zh) 一种基于非均匀网格模型的复合目标散射半解析快速计算方法
Gorji et al. Physical optics analysis for RCS computation of a relatively small complex structure
Sukharevsky et al. Simulation of scattering characteristics of aerial resonant-size objects in the VHF band
CN113610988A (zh) 一种基于射线追踪的复杂结构目标电磁散射仿真方法
Xiong et al. Modeling of tapered anechoic chambers
Nechitaylo et al. Electromagnetic wave scattering by a screen of finite thickness and conductivity
Ghanmi et al. Radar cross section of modified target using Gaussian beam methods: experimental validation
Martin et al. Radar cross section measurements and simulations of a model airplane in the X-band
CN112731325B (zh) 一种进气道类腔体目标散射测量方法
Mao et al. The Research of Standard Model's RCS Calculation and its Analytical Efficiency Based on FASTEM-STUDIO
Oezdemir et al. A fast and efficient RCS calculation and ISAR image formation tool: pRediCS
Chen et al. Scattering center modeling for low-detectable targets
Setterdahl et al. Signa: A radar signature prediction tool for electrically large targets
Zhu et al. A Simulation Method of Radar Cross Section based on FEKO
CN116882055A (zh) 一种复杂涂覆目标电磁散射模拟与散射机理分析方法及设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant