CN113255190B - 一种不规则非均匀介质体雷达散射截面的计算方法 - Google Patents
一种不规则非均匀介质体雷达散射截面的计算方法 Download PDFInfo
- Publication number
- CN113255190B CN113255190B CN202110632551.XA CN202110632551A CN113255190B CN 113255190 B CN113255190 B CN 113255190B CN 202110632551 A CN202110632551 A CN 202110632551A CN 113255190 B CN113255190 B CN 113255190B
- Authority
- CN
- China
- Prior art keywords
- dielectric constant
- target
- dielectric
- dielectric body
- model
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
一种不规则非均匀介质体雷达散射截面的计算方法,步骤如下:建立目标介质体模型,确定目标介质体模型的介电常数空间分布;基于目标介质体的最大介电常数对目标介质体模型进行网格剖分;生成目标介质体模型的介电常数采样点,由介电常数空间分布得到介电常数采样点对应的介电常数值,提取每个剖分单元的顶点坐标,将每个剖分单元的各顶点和目标介质体模型中每个介电常数采样点进行遍历,确定每个剖分单元内的介电常数采样点,将每个剖分单元的所有介电常数采样点的介电常数值求平均,得到的平均值即为该剖分单元的介电常数值;采用矩量法计算目标介质体的雷达散射截面。本发明可以实现对任意不规则外形、复杂非均匀介质体更精确的散射特性分析。
Description
技术领域
本发明属于电磁散射计算技术领域,特别涉及一种不规则、非均匀介质体的雷达散射截面(RCS)的计算方法,尤其是大气湍流对视距链路及超视距链路的影响效应计算方法。
背景技术
不规则、非均匀介质体的电磁散射特性计算的需求在日益增加,如涂覆非均匀介质的天线散射、脑散射、地壳散射、尾流散射以及大气湍流散射对视距链路、超视距链路效应的应用等技术问题。目前不规则、非均匀介质体的散射特性的计算一般是将介质体近似为简化模型后,再进行电磁散射特性的计算。例如,飞行器运动时引起的大气扰动压强分布是一种不规则、非均匀介质体,现有方法在对其散射特性进行计算时,径向密度梯度模型是将尾流单涡建模为圆柱体,绝热传输模型是将尾流近似为横截面半径与涡距有关的圆,求解湍流散射效应时一般采用估算模型,结构常数和频率决定雷达散射截面的大小。这些计算方法都是将不规则、非均匀介质模型进行近似简化,这些近似降低了电磁散射特性计算的准确性。
发明内容
本发明的目的是提供一种更准确的计算不规则、非均匀介质体的雷达散射截面的计算方法。
为了实现上述目的,本发明采取如下的技术解决方案:
一种不规则、非均匀介质体的雷达散射截面的计算方法,包括以下步骤:
S1、建立目标介质体模型,确定目标介质体模型的介电常数空间分布;
S2、根据目标介质体的介电常数空间分布确定目标介质体的最大介电常数,基于目标介质体的最大介电常数对目标介质体模型进行网格剖分,剖分后得到若干个剖分单元;
S3、根据目标介质体的介电常数空间分布,确定目标介质体模型中每个剖分单元的介电常数值,具体步骤如下:生成目标介质体模型的介电常数采样点,根据目标介质体模型的介电常数空间分布得到目标介质体模型中每个介电常数采样点对应的介电常数值,提取每个剖分单元的顶点坐标,将每个剖分单元的各顶点和目标介质体模型中每个介电常数采样点进行遍历,从而确定每个剖分单元内的介电常数采样点,将每个剖分单元的所有介电常数采样点的介电常数值求平均,得到的平均值即为该剖分单元的介电常数值;
S4、根据目标网格模型和每个剖分单元的介电常数值,采用矩量法计算目标介质体的雷达散射截面。
更具体的,所述体剖分单元为四面体形状。
优选的,采用蒙特卡洛方法确定介电常数空间分布,解决随机场问题。
由以上技术方案可知,本发明方法将不规则、非均匀介质体进行网格剖分,近似为若干小的均匀介质四面体的组合,当小介质体足够多时,各个介电常数均匀的小介质体的耦合散射特性可以近似为非均匀介质体的电磁散射特性,然后对每一个网格剖分单元进行介电常数赋值。通过将不规则、非均匀介质体的介电常数值设为最大值,并且剖分单元满足网格剖分精度要求,使得建模粒度达到最小剖分单元级别,相对于其他近似模型而言提高了结果的准确性,从而实现对任意不规则外形、复杂非均匀介质体的精细化电磁散射建模与更精确的散射特性分析。
附图说明
图1为本发明方法的流程图;
图2为本发明实施例目标介质体模型进行网格剖分后的示意图;
图3为本发明实施例目标介质体的介电常数空间分布情况示意图;
图4为剖分单元的介电常数重新赋值后的模型示意图;
图5为SWG基函数的示意图;
图6为采用本发明方法计算的RCS和采用湍流散射估算模型计算的RCS的结果对比图。
以下结合附图对本发明的具体实施方式作进一步详细地说明。
具体实施方式
下面结合附图对本发明进行详细描述,在详述本发明实施例时,为便于说明,表示器件结构的附图会不依一般比例做局部放大,而且所述示意图只是示例,其在此不应限制本发明保护的范围。需要说明的是,附图采用简化的形式且均使用非精准的比例,仅用以方便、清晰地辅助说明本发明实施例的目的。
如图1所示,本发明方法包括以下步骤:
S1、建立目标介质体模型,确定目标介质体模型的介电常数空间分布。本发明的目标介质体是不规则、非均匀介质体,建模和确定介电常数空间分布可根据不同的情况选择相应的方法实现,例如,如果已知不均匀介质体的物理特性参数分布,如通过ANSYS等软件计算出大气扰动压强起伏空间分布,可基于介质体的外部轮廓建立几何模型后,再根据目标介质体的介电常数和物理特性参数的关系得到模型的介电常数空间分布;如果已知目标介质体的空间谱函数,如湍流空间谱,则可通过蒙特卡洛方法得到目标介质体的介电常数空间分布后,再建立几何模型。
S2、根据目标介质体的介电常数空间分布确定目标介质体的最大介电常数,根据目标介质体的最大介电常数对目标介质体模型进行网格剖分,剖分后得到若干个剖分单元,本实施例的剖分单元为四面体,本发明方法的特点是基于目标介质体的最大介电常数对模型进行网格剖分,但模型的网格剖分方法可采用现有方法,如何进行网格剖分不是本发明的创新点,在此不做赘述。网格剖分与散射计算方法、入射波频率以及介电常数分布相关,对目标介质体模型进行网格剖分时,介电常数越大,剖分尺寸越小,散射计算结果越精确,因此网格剖分时选用目标介质体的介电常数最大值,从而剖分后的剖分单元可以满足最小剖分精度要求,使得每个剖分单元之间的耦合散射特性能表示目标介质体的散射特性,高精度地模拟不均匀介质的散射特性。
S3、根据目标介质体的介电常数空间分布,确定目标介质体模型中每个剖分单元的介电常数值。生成目标介质体模型的介电常数采样点,根据步骤S1中确定的目标介质体模型的介电常数空间分布得到目标介质体模型中每个介电常数采样点(空间采样点)对应的介电常数值,提取每个剖分单元的顶点坐标(V1、V2、V3、V4),将每个剖分单元的各顶点和目标介质体模型中每个介电常数采样点进行遍历,从而确定剖分单元内的介电常数采样点,将每个剖分单元的所有介电常数采样点的介电常数值求平均,得到的平均值即为该剖分单元的介电常数值;
S4、根据目标网格模型和每个剖分单元的介电常数值,采用矩量法计算目标介质体的雷达散射截面。本发明采用矩量法(Method of Moments,MOM)来计算目标介质体的雷达散射截面,矩量法是一种严格的数值方法,精度主要取决于目标几何建模精度和正确的基、权函数的选择以及阻抗元素的计算等,其本身的计算精度可忽略不计,因此可以得到不规则非均匀介质体的RCS精确解。
下面以计算强湍流的后向雷达散射截面(RCS)为例,对本发明方法进行详细说明。
S1、根据强湍流的外部轮廓进行建模,得到目标介质体的几何模型后,根据蒙特卡洛方法,结合湍流的空间谱函数得到强湍流在空间中的介电常数分布情况。
强湍流的介电常数空间起伏分布模型为式中的/>表示大气折射率起伏空间分布,/>大气折射率起伏空间分布认为是由大量谐波叠加而成,可以表示为一高斯随机数用大气湍流折射率起伏的空间谱进行滤波,再通过逆傅里叶变换得到,/>为空间某点(x,y,z)的位置矢量,即源点的位置,/>为厄米复高斯随机数三维数组,/>为三维空间波数,为空间谱,i为虚数单位,e为自然常数,/>K=2π/l,且2π/L0<K<2π/l0,l为目标介质体尺度,L0为湍流运动的外尺度,l0为湍流运动的内尺度。Cn为结构常数,其取值与天气状态以及时间有关,强起伏湍流量级为10-7m-1/3,从而根据得到局部区域强湍流的空间分布。本发明优选使用蒙特卡洛方法来确定强湍流的介电常数空间分布,该方法可以体现出湍流的空间分布和介电常数分布特征,为得到湍流问题的数值解做好了基础,使得湍流的散射特性运算不再只依赖于估算模型。
本实施例建立了“类云朵”形状的不规则、非均匀的目标介质体模型,该几何体在x方向上的最大长度为0.50m,在y方向上的最大长度为0.18m,在z方向上的最大长度为0.25m。该模型的介电常数空间分布服从强湍流情况下介电常数的空间起伏分布模型,根据蒙特卡洛方法生成目标介质体模型的介电常数空间分布,湍流运动内尺度l0设置为0.003m,外尺度为目标体的几何尺寸。
S2、根据目标介质体的介电常数空间分布确定目标介质体的最大介电常数,基于目标介质体的最大介电常数将得到的目标介质体模型进行网格剖分,本实施例得到的强湍流的最大介电常数为1.0002015,实际剖分尺寸满足:当入射频率为1GHz,得到37076个四面体形状的剖分单元,并保存各剖分单元顶点的坐标,剖分结果如图2所示。
S3、本实施例的目标介质体为强湍流,因此目标介质体模型的介电常数分布服从强湍流情况下介电常数的空间起伏分布模型,依据强湍流的介电常数空间起伏分布模型生成13万个介电常数采样点(采样点越多越精确,可根据实际应用场景和应用需求确定采样点数量),得到的介电常数空间分布如图3所示。步骤S2中对目标介质体模型进行了网格剖分,并得到了模型中每个剖分单元(四面体)的顶点坐标,将每个剖分单元的顶点和目标介质体模型中每个介电常数采样点进行遍历,确定由顶点坐标确定的每个剖分单元内的介电常数采样点,将每个剖分单元的所有介电常数采样点的介电常数值求平均,得到的平均值即为该剖分单元的介电常数值。
假设某个剖分单元的顶点坐标分别为V1=(x1,y1,z1)、V2=(x2,y2,z2)、V3=(x3,y3,z3)、V4=(x4,y4,z4)],某一介电常数采样点的坐标为P=(x,y,z),则有:
如果D0和D1、D2、D3、D4同正负,则P点在剖分单元中。
得到剖分单元内介电常数采样点的介电常数值后,将每个剖分单元内的所有介电常数采样点的介电常数值求平均,得到的平均值就是剖分单元的介电常数值,式中的εj表示剖分单元中第j个采样点的介电常数值,ni表示第i个剖分单元的介电常数采样点的个数(包括表面的采样点)。
S4、将计算得到的每个剖分单元的介电常数值赋值给每个剖分单元,重新赋值之前各剖分单元的介电常数值为目标介质体的介电常数最大值,图4为剖分单元的介电常数重新赋值后的模型图,为了便于区分,对剖分单元的介电常数进行重新赋值后,介电常数值有变化的剖分单元采用深浅不同的颜色来表示,从图4可以看出介质体的不均匀性。本发明采用先建模并网格剖分后,再对每个剖分单元进行介电参数赋值,和对每个介电参数不均匀单元分别建模和赋值的现有方法相比,可以通过程序实现自动赋值,大大减小了工作量及时间,降低了建模的复杂程度。
S5、根据目标网格模型和每个剖分单元的介电常数值,采用矩量法计算目标介质体的雷达散射截面。
目标介质体的雷达散射截面式中的k0为空间波数,/>为微分算符,/>为对比率张量,/>为空间Green函数,/>为场点,/>为源点,V表示目标介质体体积,/>为电通密度矢量,ε0为自由空间的介电常数,当整个体域被四面体单元剖分后,采用SWG基函数展开为/>F和H分别表示Full-SWG和Half-SWG(图4),/>为Full-SWG基函数,/>为Half-SWG基函数,/>和分别表示Full-SWG和Half-SWG情况下的未知系数,NF和NH为对应基函数下未知量的个数。
图5为采用湍流散射估算模型得到的RCS和采用本发明方法得到的RCS的结果对比图。根据物理光学法可知,目标的后向RCS随频率的增大呈现一个震荡起伏的趋势,且大气湍流在不同频率下的回波强度变化较大,从图5可以看出来,估算模型计算得到的RCS随频率的变化较小,且呈现稳定的上升趋势,这与理论有所出入,而本发明方法与理论相符,计算结果更为精确。
本发明将离散和局部均匀近似的思想相结合,将不规则、不均匀介质体剖分为若干小的均匀介质四面体,当小介质体足够多时,各个介电常数均匀的小介质体的耦合散射特性可以近似为非均匀介质体的电磁散射特性,并对每个剖分单元自动化赋予介电参数的值,最后使用数值方法计算得到不均匀目标的雷达散射截面。该方法建立了精确的不均匀体介电参数空间分布模型,且可针对空间谱建立随机场的介电参数空间分布特征,为得到湍流问题的数值解做好了基础,从而提高了不规则、不均匀体的电磁散射特性计算的准确性。
以上所述,仅是本发明的较佳实施例而已,并非对本发明做任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (1)
1.一种不规则非均匀介质体雷达散射截面的计算方法,其特征在于,包括以下步骤:
S1、建立目标介质体模型,确定目标介质体模型的介电常数空间分布,步骤如下:基于目标介质体的外部轮廓进行建模,所述目标介质体为强湍流,得到目标介质体的几何模型,强湍流的介电常数空间起伏分布模型为式中的/>表示大气折射率起伏空间分布,根据蒙特卡洛方法生成目标介质体模型的介电常数空间分布;
S3、根据目标介质体的介电常数空间分布,确定目标介质体模型中每个剖分单元的介电常数值,具体步骤如下:生成目标介质体模型的介电常数采样点,根据目标介质体模型的介电常数空间分布得到目标介质体模型中每个介电常数采样点对应的介电常数值,提取每个剖分单元的顶点坐标,将每个剖分单元的各顶点和目标介质体模型中每个介电常数采样点进行遍历,从而确定每个剖分单元内的介电常数采样点,将每个剖分单元的所有介电常数采样点的介电常数值求平均,得到的平均值即为该剖分单元的介电常数值;
S4、根据网格剖分后的目标介质体模型和每个剖分单元的介电常数值,采用矩量法计算目标介质体的雷达散射截面;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110632551.XA CN113255190B (zh) | 2021-06-07 | 2021-06-07 | 一种不规则非均匀介质体雷达散射截面的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110632551.XA CN113255190B (zh) | 2021-06-07 | 2021-06-07 | 一种不规则非均匀介质体雷达散射截面的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113255190A CN113255190A (zh) | 2021-08-13 |
CN113255190B true CN113255190B (zh) | 2023-06-20 |
Family
ID=77186851
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110632551.XA Active CN113255190B (zh) | 2021-06-07 | 2021-06-07 | 一种不规则非均匀介质体雷达散射截面的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113255190B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106446470A (zh) * | 2016-11-17 | 2017-02-22 | 上海无线电设备研究所 | 一种高效并行的非均匀介质频域有限差分方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000046884A (ja) * | 1998-07-28 | 2000-02-18 | Showa Electric Wire & Cable Co Ltd | 空間電荷測定における電荷密度校正方法 |
KR101347629B1 (ko) * | 2012-09-21 | 2014-01-07 | 고려대학교 산학협력단 | 무반사막, 및 무반사막 형성 방법 |
US20160120407A1 (en) * | 2013-06-18 | 2016-05-05 | Northeastern University | Microwave nearfield radar imaging (nri) using digital breast tomosynthesis (dbt) for non-invasive breast cancer detection |
CN103513235B (zh) * | 2013-09-16 | 2016-03-09 | 中国人民解放军国防科学技术大学 | 晴空飞机尾流稳定段雷达散射特性计算方法 |
CN104778293B (zh) * | 2014-01-15 | 2018-06-12 | 南京理工大学 | 非均匀介质目标电磁散射的体积分Nystrom分析方法 |
CN104573240B (zh) * | 2015-01-12 | 2017-09-15 | 西安电子科技大学 | 周期性非均匀介质波导特征模分析的七点频域有限差分方法 |
CN105548971B (zh) * | 2016-01-11 | 2017-12-15 | 西安电子科技大学 | 一种针对海面蒸发波导超视距探测雷达的目标隐身方法 |
US10374316B2 (en) * | 2016-10-21 | 2019-08-06 | At&T Intellectual Property I, L.P. | System and dielectric antenna with non-uniform dielectric |
CN110274920B (zh) * | 2018-03-14 | 2021-11-23 | 南京理工大学 | 体面剖分弹跳射线法分析金属介质目标瞬态特性的方法 |
CN109765538B (zh) * | 2019-01-11 | 2022-12-02 | 西安电子科技大学 | 非均匀介质目标体的散射场确定方法 |
-
2021
- 2021-06-07 CN CN202110632551.XA patent/CN113255190B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106446470A (zh) * | 2016-11-17 | 2017-02-22 | 上海无线电设备研究所 | 一种高效并行的非均匀介质频域有限差分方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113255190A (zh) | 2021-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107576388B (zh) | 一种浅海信道下三维结构声源辐射声场预报方法 | |
Bres et al. | Large-eddy simulations of co-annular turbulent jet using a Voronoi-based mesh generation framework | |
CN112327374B (zh) | Gpu探地雷达复杂介质dgtd正演方法 | |
CN113177356A (zh) | 一种基于深度学习的目标电磁散射特性快速预测方法 | |
CN115169235A (zh) | 一种基于改进生成对抗网络的超表面单元结构逆设计方法 | |
CN106446470A (zh) | 一种高效并行的非均匀介质频域有限差分方法 | |
CN110362877B (zh) | 不确定因素的环境电磁散射特性分析方法 | |
Zhou et al. | Target head direction far field radar cross section reduction based on narrow surface element method | |
CN113255190B (zh) | 一种不规则非均匀介质体雷达散射截面的计算方法 | |
CN116151135B (zh) | 一种电大尺寸目标的电磁仿真方法及系统 | |
CN113758559A (zh) | 基于散射声场分离算法在浅海信道中测量方法 | |
CN104778286B (zh) | 掠海飞行器电磁散射特性快速仿真方法 | |
CN112785052A (zh) | 基于粒子滤波算法的风速风向预测方法 | |
CN109359399B (zh) | 一种计算高速微动类球头柱体目标的电磁散射特性方法 | |
CN109241556B (zh) | 一种随机粗糙体目标建模方法及存储介质 | |
Belibassakis et al. | Three-dimensional Green's function for harmonic water waves over a bottom topography with different depths at infinity | |
CN113177307A (zh) | 一种应用于复杂地形的短波地波仿真方法及系统 | |
CN113821967A (zh) | 基于散射中心模型的大样本训练数据生成方法 | |
CN116484724B (zh) | 基于高阶矩量法电磁建模参数的天线检测方法 | |
Zhou et al. | Offshore surface evaporation duct joint inversion algorithm using measured dual-frequency sea clutter | |
CN113536645A (zh) | 基于高斯过程的带罩反射面天线电性能预测方法 | |
CN109241595A (zh) | 一种基于特征正交分解的天线电性能的快速区间分析方法 | |
Yang et al. | Application and error analysis of narrow-angle parabolic equation method in electromagnetic wave propagation at troposphere | |
Zhao et al. | A new method to measure the dielectric constant in complex electromagnetic environment with black-box modeling regression method | |
CN117332658B (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 |