CN110059422B - 频率选择表面曲面天线罩的电磁散射特性仿真方法 - Google Patents
频率选择表面曲面天线罩的电磁散射特性仿真方法 Download PDFInfo
- Publication number
- CN110059422B CN110059422B CN201910330261.2A CN201910330261A CN110059422B CN 110059422 B CN110059422 B CN 110059422B CN 201910330261 A CN201910330261 A CN 201910330261A CN 110059422 B CN110059422 B CN 110059422B
- Authority
- CN
- China
- Prior art keywords
- region
- radome
- scattering
- calculating
- field
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Details Of Aerials (AREA)
Abstract
本发明涉及频率选择表面曲面天线罩的电磁散射特性仿真方法,包括步骤有:(1)对频率选择表面天线罩按照曲率变化情况进行区域分解;(2)对每个区域的电磁特性分别利用子阵策略进行独立计算;(3)将区域叠加,得到频率选择表面天线罩上的初始散射场;(4)用区域分解矩量法计算各个区域间的耦合影响,更新每个区域的激励矢量,通过子阵策略计算出每个区域中新散射场;(5)对各个区域中新的散射场进行叠加,得到所述频率选择表面天线罩的散射总场。本发明解决了复杂外形大型频选天线罩运用全波方法分析时效率低下,甚至无法计算的问题,且降低了频率选择表面天线罩的设计风险和成本、缩短了其设计周期、提高了研究人员工作效率。
Description
技术领域
本发明涉及计算机仿真技术领域,尤其涉及频率选择表面曲面天线罩的电磁散射特性仿真方法。
背景技术
作为天线的保护装置,天线罩广泛应用于军事、民用无线电领域,如飞机的机头雷达罩、地空、空空导引头天线罩、机翼电子战天线罩、机背卫通天线罩等,几乎有天线之处就会有天线罩的存在。频率选择表面(FSS,Frequency Selective Surfaces)是一种通过表面结构设计达到调控电磁波的传输与散射特性的空间滤波结构,将其应用于天线罩,在保护天线的同时可显著提高天线的隐身能力,具有重要价值。
目前在频率选择表面天线罩电磁特性仿真设计中,还未有一个可供设计使用的公式和理论,大部分的研究工作还是出于试验尝试的方法进行设计。对于一些大型频率选择表面天线罩结构,由于复杂外形以及对电磁散射特性的需求破坏了周期边界条件,数值分析方法与全波仿真软件的计算能力受到了极大挑战,在现有的计算机硬件条件下效率低下甚至无法计算。同时随着频率选择表面天线罩中包含的金属/介质层状结构更加复杂,进一步加剧了研究的困难。因此,开发高效高精度计算技术来评估频率选择表面天线罩的电磁特性,对降低频选罩的设计风险和成本、缩短其设计周期、提高研究人员工作效率等工程实践具有重要价值。
发明内容
本发明要解决的技术问题在于,针对现有技术中的缺陷,提供了一种频率选择表面曲面天线罩的电磁散射特性仿真方法,解决了复杂外形大型频选天线罩运用全波方法分析时效率低下,甚至无法计算的问题,在保证计算精度的前提下大大提高了计算效率。
本发明提供的频率选择表面曲面天线罩的电磁散射特性仿真方法,其改进之处在于:所述方法包括如下步骤:
(1)对频率选择表面天线罩按照曲率变化情况进行区域分解;
(2)对每个区域的电磁特性分别利用子阵策略进行独立计算;
(3)将区域叠加,得到频率选择表面天线罩上的初始散射场;
(4)用区域分解矩量法计算各个区域间的耦合影响,更新每个区域的激励矢量,通过子阵策略计算出每个区域中新散射场;
(5)对各个区域中新的散射场进行叠加,得到所述频率选择表面天线罩的散射总场。
其中:步骤(1)包括:
计算曲面模型上各点的高斯曲率和平均曲率,以天线罩模型中部曲率最接近0的一点为基准点,沿天线罩轴向或母线向尖端和边缘两侧寻找,将高斯曲率和平均曲率沿天线罩轴向或母线变化小于1.2倍的区域划分为一个区域,再以划分出的区域边缘为新的基准点,继续沿天线罩轴向或母线向尖端和边缘两侧寻找,继续以变化小于1.2倍的区域划分,直至所述曲面全部划分完成。
其中:步骤(2)基于子阵策略将每个区域均划分为n个更小的子阵;
在每个区域内,采用区域分解矩量法对该区域内的n个子阵进行仿真计算,得到n个子阵中单元的散射场,通过等效计算,得到n个子阵的电流和磁流,利用子阵策略,整合得到该区域的散射场和该区域上的分解耦合的初始电流和磁流。
其中:步骤(5)中,当各个区域中新的散射场进行叠加,得到所述频率选择表面天线罩的散射总场,并与初始散射场的相对误差不满足收敛时,返回步骤(4)重新计算各个区域间的耦合影响,直至所述相对误差满足收敛。
其中:步骤(5)计算所述频率选择表面天线罩的散射总场时,将FSS阵列中所有单元的法向局部坐标系下的散射场转化为局部坐标系下的散射场,再将局部坐标系下的散射场转化为全局坐标系下的散射场,得到FSS阵列的总散射场为:
实施本发明的,具有以下有益效果:
对各个区域的电磁特性,分别利用子阵策略进行独立快速计算:在各个区域内表面可近似看作规则表面,引入子阵策略对其电磁特性进行计算,利用一小规模FSS的电磁特性来外推计算大规模目标FSS的电磁特性,以在保证计算精度的前提下大幅降低计算量,得到各个子区域上的初始电流和磁流;
本发明结合了子阵策略与区域分解矩量法两者的优势,解决了复杂外形大型频选天线罩运用全波方法分析时效率低下,甚至无法计算的问题,在保证计算精度的前提下,降低了频率选择表面天线罩的设计风险和成本、缩短了其设计周期、提高了研究人员工作效率。该方法合理可控,可以通过自由选择区域的个数和大小来调节建模精度与效率。本发明为分区域建模,适用于频选天线罩的电磁散射特性的并行仿真计算。
附图说明
图1是本发明实施例的流程示意图;
图2是本发明实施例的单元全局、法向局部和局部坐标系三种坐标系示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了降低设计风险和成本浪费,需要借助高精度计算方法对复杂外形FSS天线罩的电磁特性进行模拟。目前,受限于计算机硬件资源,运用全波方法分析复杂外形大型FSS天线罩效率低下,甚至无法计算。针对上述问题,本发明提出频率选择表面曲面天线罩的电磁散射特性仿真方法,其流程图如图1所示,具体包括如下步骤:
(1)对频率选择表面天线罩按照曲率变化情况进行区域分解;其原则是在该区域内天线罩表面曲率变化较小,具体由用户设定。
计算曲面模型上各点的高斯曲率和平均曲率,可以以天线罩模型中部曲率最接近0的一点为基准点,沿天线罩轴向或母线向尖端和边缘两侧寻找,将高斯曲率和平均曲率沿天线罩轴向或母线变化小于1.2倍的区域划分为一个区域,然后以划分出的区域边缘为新的基准点,继续沿天线罩轴向或母线向尖端和边缘两侧寻找,并以同样的1.2倍变化准则划分出新的区域。本实施例通过统计曲面模型上各点的高斯曲率和平均曲率,设定区域划分所依据的曲率变化阈值,控制区域分解矩量法的计算精度和效率。
(2)对每个区域的电磁特性分别利用子阵策略进行独立计算;本实施例基于子阵策略将每个区域均划分为n个更小的子阵,在每个区域内,采用区域分解矩量法对该区域内的n个子阵进行仿真计算,得到n个子阵中单元的散射场,通过等效计算,得到n个子阵的电流和磁流,利用子阵策略,整合得到该区域的散射场和该区域上的分解耦合的初始电流和磁流。
(3)将区域叠加,得到频率选择表面天线罩上的初始散射场;
(4)用区域分解矩量法计算各个区域间的耦合影响,更新每个区域的激励矢量,通过子阵策略计算出每个区域中新散射场;
(5)对各个区域中新的散射场进行叠加,得到所述频率选择表面天线罩的散射总场,并与初始散射场的相对误差不满足收敛时,返回步骤(4)重新计算各个区域间的耦合影响,直至所述相对误差满足收敛,得到所述频率选择表面天线罩的散射总场。具体的,本实施例通过对子阵策略进行改进,通过精确控制曲面到平面坐标的变换,从法向局部坐标到局部坐标再到全局坐标的变换,获取单元阵子的散射特性,外推得到整体阵列的散射特性,使目标整体的区域分割和重构过程有法可依。
由于共形FSS中各单元轴向不同,各单元的RCS曲线形状类似,但方向图指向不同,仅在各自的法向局部坐标系下才可等效,因此除了采用区域分解矩量法依次计算得到FSS阵列中每个单元在各自法向局部坐标系下的散射场之外,还需要通过坐标变换来得到全局坐标系下的FSS单元及其阵列的散射场。
为了方便介绍,本实施例取共形FSS阵列中的第n个天线单元,建立如图2所示的坐标示意图。图中on是第n个FSS单元的相位中心,其在全局坐标系下的坐标为(xn,yn,zn)。以on为原点,建立该单元的法向局部坐标系其中轴的正方向为该单元放置的曲面局部外法线方向,和轴的指向以该单元放置在曲面上的具体方位进行定义。为了简化运算,以on为原点,再建立该单元的局部坐标系其中和轴的正方向分别与全局坐标系的x、y和z轴的正方向相同。
对于FSS阵列而言,其总散射场等于各单元散射场的叠加:
式中,En(θ,φ)为第n个FSS单元的散射场,在全局坐标系下,其表示为:
且
将FSS阵列中所有单元的法向局部坐标系下的散射场转化为局部坐标系下的散射场,再将局部坐标系下的散射场转化为全局坐标系下的散射场,代入式(1)即可得FSS阵列的总散射场为:
上式即为共形FSS阵列散射场求解方法的一般表达式,若依次计算FSS阵列中每个单元在各自法向局部坐标系下的散射场,则由式(6)式计算所得的FSS阵列散射总场,完全考虑了单元间的互耦。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (4)
1.一种频率选择表面曲面天线罩的电磁散射特性仿真方法,其特征在于:所述方法包括如下步骤:
(1)对频率选择表面天线罩按照曲率变化情况进行区域分解;
(2)对每个区域的电磁特性分别利用子阵策略进行独立计算;
(3)将区域叠加,得到频率选择表面天线罩上的初始散射场;
(4)用区域分解矩量法计算各个区域间的耦合影响,更新每个区域的激励矢量,通过子阵策略计算出每个区域中新散射场;
(5)对各个区域中新的散射场进行叠加,得到所述频率选择表面天线罩的散射总场;
步骤(1)包括:
计算曲面模型上各点的高斯曲率和平均曲率,以天线罩模型中曲率最接近0的一点为基准点,沿天线罩轴向或母线向尖端和边缘两侧寻找,将高斯曲率和平均曲率沿天线罩轴向或母线变化小于1.2倍的区域划分为一个区域,再以划分出的区域边缘为新的基准点,继续沿天线罩轴向或母线向尖端和边缘两侧寻找,继续以变化小于1.2倍的区域划分,直至所述曲面全部划分完成。
2.根据权利要求1所述的方法,其特征在于:步骤(2)基于子阵策略将每个区域均划分为n个更小的子阵;
在每个区域内,采用区域分解矩量法对该区域内的n个子阵进行仿真计算,得到n个子阵中单元的散射场,通过等效计算,得到n个子阵的电流和磁流,利用子阵策略,整合得到该区域的散射场和该区域上的分解耦合的初始电流和磁流。
3.根据权利要求1所述的方法,其特征在于:步骤(5)中,将各个区域中新的散射场进行叠加,得到所述频率选择表面天线罩的散射总场,并与初始散射场的相对误差不满足收敛时,返回步骤(4)重新计算各个区域间的耦合影响,直至所述相对误差满足收敛。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910330261.2A CN110059422B (zh) | 2019-04-23 | 2019-04-23 | 频率选择表面曲面天线罩的电磁散射特性仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910330261.2A CN110059422B (zh) | 2019-04-23 | 2019-04-23 | 频率选择表面曲面天线罩的电磁散射特性仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110059422A CN110059422A (zh) | 2019-07-26 |
CN110059422B true CN110059422B (zh) | 2023-01-20 |
Family
ID=67320219
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910330261.2A Active CN110059422B (zh) | 2019-04-23 | 2019-04-23 | 频率选择表面曲面天线罩的电磁散射特性仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110059422B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110889216B (zh) * | 2019-11-20 | 2023-11-14 | 上海无线电设备研究所 | 一种曲面频率选择表面天线罩的自适应快速设计方法 |
CN111368398B (zh) * | 2020-02-06 | 2023-12-29 | 北京环境特性研究所 | 一种不确定结构电大目标电磁散射特性分析方法及装置 |
CN112070697B (zh) * | 2020-09-07 | 2024-01-19 | 北京环境特性研究所 | 一种基于逆二维像的局部散射特性评估方法及装置 |
CN115034075B (zh) * | 2022-06-22 | 2023-11-24 | 电子科技大学 | 非规则天线阵列矢量增益方向图的快速精准预测方法 |
CN117057098B (zh) * | 2023-07-13 | 2024-04-26 | 大连理工大学 | 基于物理光学法和等效介质模型的频选天线罩分析方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104992001A (zh) * | 2015-06-19 | 2015-10-21 | 中国人民解放军信息工程大学 | 大规模mimo阵列天线远场辐射场的精确快速计算方法 |
CN106295120A (zh) * | 2016-07-26 | 2017-01-04 | 上海无线电设备研究所 | 一种目标电磁散射特性快速计算方法 |
CN108268696A (zh) * | 2017-12-15 | 2018-07-10 | 西安电子科技大学 | 一种适用于高阶矩量法的fss天线罩建模方法 |
CN108615974A (zh) * | 2018-05-03 | 2018-10-02 | 航天特种材料及工艺技术研究所 | 耗散型宽频带频率选择表面雷达罩 |
CN108664746A (zh) * | 2018-05-17 | 2018-10-16 | 中国人民解放军陆军工程大学 | 一种基于双线传输线模型的频率选择表面近似解析方法 |
-
2019
- 2019-04-23 CN CN201910330261.2A patent/CN110059422B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104992001A (zh) * | 2015-06-19 | 2015-10-21 | 中国人民解放军信息工程大学 | 大规模mimo阵列天线远场辐射场的精确快速计算方法 |
CN106295120A (zh) * | 2016-07-26 | 2017-01-04 | 上海无线电设备研究所 | 一种目标电磁散射特性快速计算方法 |
CN108268696A (zh) * | 2017-12-15 | 2018-07-10 | 西安电子科技大学 | 一种适用于高阶矩量法的fss天线罩建模方法 |
CN108615974A (zh) * | 2018-05-03 | 2018-10-02 | 航天特种材料及工艺技术研究所 | 耗散型宽频带频率选择表面雷达罩 |
CN108664746A (zh) * | 2018-05-17 | 2018-10-16 | 中国人民解放军陆军工程大学 | 一种基于双线传输线模型的频率选择表面近似解析方法 |
Non-Patent Citations (3)
Title |
---|
A Non-overlapping and Non-conformal Domain Decomposition Method with Second Order Transmission Condition for Modelling Large Finite Antenna Arrays;Zhen Peng等;《2010 URSI International Symposium on Electromagnetic Theory》;20101115;第160-163页 * |
基于改进的体面积分方程区域分解方法高效求解有限大频率选择表面结构的电磁散射特性;李先进等;《电波科学学报》;20190215;第34卷(第1期);第76-82页 * |
频率选择表面技术及在天线罩设计方面的应用;王平友;《中国优秀硕士学位论文全文数据库》;20180215(第2期);第13-56页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110059422A (zh) | 2019-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110059422B (zh) | 频率选择表面曲面天线罩的电磁散射特性仿真方法 | |
CN106886656B (zh) | 一种三维阵列天线方向图旁瓣抑制方法 | |
CN106407723B (zh) | 面向低副瓣的稀疏排布阵列天线激励电流幅度的确定方法 | |
CN110047136B (zh) | 一种基于频率选择表面的天线罩三维建模方法 | |
CN111276822B (zh) | 一种天线方向图主瓣可控的天线阵列设计方法 | |
CN104036078B (zh) | 一种基于安装高度的阵列天线辐射和散射性能综合设计方法 | |
CN103246781A (zh) | 基于空间映射的阵列天线雷达散射截面减缩方法 | |
CN104933213A (zh) | 基于空间映射的大规模相控天线阵列宽角扫描优化方法 | |
CN107017468A (zh) | 一种机身共形阵列天线的单元布局与电磁设计方法 | |
CN102122765A (zh) | 直线阵列天线方向图的优化方法 | |
CN107357962A (zh) | 一种基于自适应代理模型的天线罩肋截面尺寸优化方法 | |
CN110276086B (zh) | 基于空间映射的隐身超表面雷达散射截面减缩方法 | |
CN102054094A (zh) | 一种平面微带电路的快速多层方向性仿真方法 | |
Wang et al. | An adjustment method for active reflector of large high-frequency antennas considering gain and boresight | |
CN107729627A (zh) | 机会阵列天线的单元快速选取方法 | |
CN111914364B (zh) | 基于高阶矩量法与投影的频选天线罩建模方法 | |
CN113221370B (zh) | 基于锥面投影的fss天线罩建模方法 | |
CN109212498B (zh) | 一种飞机编队雷达散射截面的快速算法 | |
CN109325306B (zh) | 一种基于局部平面的任意曲面共形阵列建模方法 | |
CN104951617A (zh) | 基于区域分解法的天线优化设计方法 | |
CN113553711B (zh) | 一种任意曲面共形阵列天线设计方法 | |
CN116595626A (zh) | 航天器仪器舱组件多样化布局的双目标优化方法及系统 | |
CN106156431B (zh) | 基于nurbs曲面建模的导体目标电磁散射仿真方法 | |
KR102380267B1 (ko) | 표적 탐지 방법 | |
CN106503349B (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 |