CN112685877A - 一种基于离散网格分组的天线罩电性能快速误差分析方法 - Google Patents

一种基于离散网格分组的天线罩电性能快速误差分析方法 Download PDF

Info

Publication number
CN112685877A
CN112685877A CN202011463547.7A CN202011463547A CN112685877A CN 112685877 A CN112685877 A CN 112685877A CN 202011463547 A CN202011463547 A CN 202011463547A CN 112685877 A CN112685877 A CN 112685877A
Authority
CN
China
Prior art keywords
interval
antenna
sector
field
radome
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.)
Granted
Application number
CN202011463547.7A
Other languages
English (en)
Other versions
CN112685877B (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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN202011463547.7A priority Critical patent/CN112685877B/zh
Publication of CN112685877A publication Critical patent/CN112685877A/zh
Application granted granted Critical
Publication of CN112685877B publication Critical patent/CN112685877B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Details Of Aerials (AREA)

Abstract

本发明属于雷达天线技术领域,具体是一种基于离散网格分组的天线罩电性能快速误差分析方法,可用于分析误差对天线罩电性能的影响。本发明通过对天线罩网格进行分组,减少了蒙特卡洛方法中独立随机变量的数目,因而与传统蒙特卡洛方法相比,可减小天线罩误差分析所需的计算量,并提高计算精度。

Description

一种基于离散网格分组的天线罩电性能快速误差分析方法
技术领域
本发明属于雷达天线技术领域,具体是一种基于离散网格分组的天线罩电性能快速误差分析方法,可用于分析误差对天线罩电性能的影响。
技术背景
天线罩是保护天线免受自然环境影响的透波壳,是由天然或人造电介质材料制成的覆盖物,或是由桁架支撑的电介质壳体构成的特殊形状的电磁明窗。设计优良的天线罩,除了具有保护性、传导性、可靠性、隐蔽性和装饰性等功能外,还可以延长整个系统各部分的使用寿命、降低寿命成本和操作成本、简化设计、降低维修成本、保证天线表面和位置的精确度、给天线操作人员创造良好的工作环境。但是天线罩也会对理想天线的电磁辐射产生影响,使理想的天线电性能有所降低。
在加工过程中制造误差因素以及服役过程中外界环境因素的影响下,天线罩的设计参数总是存在一定的误差,从而会对天线罩的电性能产生影响。要保证天线罩的正常、可靠工作,需要有效、准确地预测误差对天线罩电性能的影响程度。蒙特卡洛方法是天线、天线罩误差分析领域中的一种经典方法,然而该方法所需的计算量与随机变量的数目直接相关。
P.Li等人在2017的论文《Power Pattern Tolerance Analysis of Radome Withthe Material Property Error Based on Interval Arithmetic》中利用蒙特卡洛方法分析了一个天线罩厚度误差对电性能的影响,其中用到的离散网格超过20000个,在3000次蒙特卡洛分析下,得到的方向图变化范围较小,难以充分反映特定厚度误差对天线罩电性能的影响程度。
发明内容
本发明的目的在于针对上述现有技术的不足,提供一种可减小天线罩误差分析所需的计算量,提高计算精度的基于离散网格分组的天线罩电性能快速误差分析方法,实现天线罩电性能的快速误差分析。
为实现上述目的,一种基于离散网格分组的天线罩电性能快速误差分析方法,其特征是:包括如下步骤:
(1)根据给定的天线罩模型,对天线罩进行网格划分;
(2)在每个网格中选取一点,按照点的坐标将其归类到10个扇形区域中,实现网格分组;
(3)根据给定的天线罩厚度值变化范围,随机生成1000组厚度误差分布,每一个扇形区域具有相同的厚度误差值;
(4)根据天线罩的结构参数和材料参数,用传输线理论计算1000组厚度误差分布下天线罩的远场方向图,并从中提取增益G1、主波束宽度B1和第一副瓣电平S1这些电性能指标的变化区间;
(5)计算天线罩上各处的透射系数的幅值和相位的变化区间,得到透射系数在复数域形成的扇形区间;
(6)将透射系数形成的扇形区间分割成多边形,扇形圆弧部分的分段数设为10;
(7)利用Minkowski求和将天线罩上各处的透射系数形成的多边形区间相加,得到远场方向图的变化区间,并从中提取增益G2、主波束宽度B2和第一副瓣电平S2这些电性能指标的变化区间;
(8)对比步骤(4)与步骤(8)中得到的天线罩电性能指标变化区间,判断天线罩的电性能指标是否满足预设要求,如果满足,则天线罩电性能分析过程结束,否则,修改步骤(2)中扇形区域的数目以及步骤(6)中扇形圆弧部分的分段数,并重复步骤(2)至步骤(8),直至得到电性能指标满足预设要求的分析结果。
所述的步骤(2)按如下过程进行:
(2a)在划分的网格中选取一点,将该坐标转换为极坐标,极坐标的周向角度记为ψ;
(2b)沿周向角度进行均匀划分,当扇形区域数为M时,从0度开始,间隔一定步长选取M个分隔点,分隔点初始最小值为P1=0,最大值为PM=360,其间任一分隔点为:
Pi=P1+i*360/M,i=1,2,…,M
将M=10代入上式即可得到每一分隔点的值;
(2c)第i个分隔区域的区间范围为[Pi,Pi+1]。
所述的步骤(3)按如下过程进行:
(3a)设天线罩的厚度误差范围为[dL,dU],则以均值为(dL+dU)/2、标准差为(dL+dU)/6生成10个随机数,对于步骤二中10个扇形区域,第一个扇形区域分配一个随机数,即为该区域中所有网格的厚度误差值;
(3b)重复步骤(3a)1000次,得到1000组对应于天线罩网格数目的厚度误差随机分布值,每一组厚度误差的随机分布值均为10个数,与划分的扇形区域数目相等。
本发明通过对天线罩网格进行分组,减少了蒙特卡洛方法中独立随机变量的数目,因而与传统蒙特卡洛方法相比,可减小天线罩误差分析所需的计算量,并提高计算精度。
附图说明
图1是本发明的实现总流程图;
图2是划分10个扇形区域的示意图;
图3是本发明使用的天线与天线罩关系示意图;
图4是本发明中计算透过天线罩后的口径场的子流程图;
图5是本发明仿真使用的某球形天线罩的外形示意图;
图6是不分组时对某球形天线罩进行分析前后所得的远场方向图区间对比图;
图7是采用本发明分组时对某球形天线罩进行分析前后所得的远场方向图区间对比图。
具体实施方式
以下参照附图对本发明作进一步详细描述。
参照图1,一种基于离散网格分组的天线罩电性能快速误差分析方法,其特征是:具体至少包括如下步骤:
步骤一,根据给定的天线罩模型,对天线罩进行网格划分。
在商用模型分析软件中,根据天线罩的结构形式建立天线罩的几何模型,设置网格边长为0.2λ,其中λ为天线的波长,对模型进行网格划分,记网格数目为N。
步骤二,在每个网格中选取一点,按照点的坐标将其归类到10个扇形区域中,实现网格分组。
(2a)在划分的网格中选取一点,将该坐标转换为极坐标,极坐标的周向角度记为ψ;
(2b)沿周向角度进行均匀划分,当扇形区域数为M时,从0度开始,间隔一定步长选取M个分隔点,分隔点初始最小值为P1=0,最大值为PM=360,其间任一分隔点为:
Pi=P1+i*360/M,i=1,2,…,M
将M=10代入上式即可得到每一分隔点的值。
(2c)第i个分隔区域的区间范围为[Pi,Pi+1],划分10个扇形区域的效果如图2所示。
步骤三,根据给定的天线罩厚度值变化范围,随机生成1000组厚度误差分布。
(3a)设天线罩的厚度误差范围为[dL,dU],则以均值为(dL+dU)/2、标准差为(dL+dU)/6生成10个随机数,对于步骤二中10个扇形区域,第一个扇形区域分配一个随机数,即为该区域中所有网格的厚度误差值。
(3b)重复步骤(3a)1000次,得到1000组对应于天线罩网格数目的厚度误差随机分布值,每一组厚度误差的随机分布值均为10个数,与划分的扇形区域数目相等。
步骤四,计算厚度误差分布下的天线罩远场方向图,并提取电性能指标变化区间。
(4a)如图3所示,沿罩体的高度,以天线罩的底面中心为原点、底面为xy平面建立一个坐标系O-xyz,罩体高度沿z方向,其中S为天线口径面,S′为透过天线罩后的口径面。
(4b)根据天线罩的结构参数和材料参数,用传输线理论计算罩体的透射系数TM,并根据已知的天线口径场E(x,y),计算透过天线罩后的口径场:E′(x,y)=E(x,y)·TM,如图4所示,具体包括如下步骤:
(4b1)针对步骤(4a)建立的坐标系,将其x、y、z方向的分量分别用i、j、k表示,已知的天线口径场分布记为E(x,y)。
(4b2)根据天线罩的几何形状与入射的口径场,求出天线罩上各点处的入射角α和极化角β,即将电磁波入射线与入射点处的法线夹角记为入射角α,将电磁波的极化方向与入射平面的夹角记为极化角β,其中入射平面由电磁波入射线与入射点处的法线构成。
(4b3)令天线罩材料的损耗角正切tanδ为0,根据天线罩各处的厚度d、相对介电常数εr,计算蒙皮上各点处的水平极化分量透射系数TH和垂直极化分量透射系数TV
Figure BDA0002832346220000061
Figure BDA0002832346220000062
其中,A=ch(jVd),
Figure BDA0002832346220000063
Figure BDA0002832346220000064
D=ch(jVd),
Figure BDA0002832346220000065
Figure BDA0002832346220000066
ZH=cosα,
Figure BDA0002832346220000067
这些参数均为中间变量;tH、tV分别为TH、TV的模值,ηH、ηV分别为TH、TV的相位。
(4b4)根据水平极化分量透射系数TH和垂直极化分量透射系数TV,得到主极化分量的透射系数:
Figure BDA0002832346220000071
其中,
Figure BDA0002832346220000072
为中间变量。
(4b5)将入射到天线罩上的口径场乘以其对应点处的透射系数,计算透过蒙皮后的口径场:E′(x,y)=E(x,y)·TM
(4c)根据透过天线罩后的口径场E′(x,y),计算带罩天线的远场F(θ,φ),绘制远场方向图,并从该远场方向图中提取增益G1、主波束位置B1和第一副瓣电平S1这些电性能指标,具体包括如下步骤:
(4c1)根据步骤(4b)中求得的透过天线罩后的口径场E′(x,y),按下式计算透过天线罩后的口径场产生的远场F(θ,φ):
Figure BDA0002832346220000073
其中,如图3所示,θ、φ是观察点P在直角坐标系O-xyz中的球坐标角,k0为自由空间传播常数,按公式
Figure BDA0002832346220000074
计算,λ是天线的波长,根据天线工作频率f和光速c,通过公式
Figure BDA0002832346220000075
计算得到。
(4c2)根据透过天线罩后的口径场产生的远场F(θ,φ)绘制加罩后的天线远场方向图,并从方向图中提取增益G1、主波束位置B1和第一副瓣电平S1这些电性能指标的变化区间。
步骤五,根据厚度误差范围计算天线罩上各处的透射系数的幅值和相位的变化区间,得到透射系数在复数域形成的扇形区间。
(5a)天线罩水平极化分量透射系数TH和垂直极化分量透射系数TV可改写为下述表达式:
Figure BDA0002832346220000081
Figure BDA0002832346220000082
其中
Figure BDA0002832346220000083
Figure BDA0002832346220000084
Figure BDA0002832346220000085
均为临时变量,根据厚度d的变化区间[dL,dU],计算临时变量A1的变化区间
Figure BDA0002832346220000086
B1的变化区间
Figure BDA0002832346220000087
以及B2的变化区间
Figure BDA0002832346220000088
(5b)tH的变化区间
Figure BDA0002832346220000089
和ηH的变化区间
Figure BDA00028323462200000810
在复数空间中构成了水平极化分量透射系数TH的扇形区间,按下式计算:
Figure BDA00028323462200000811
Figure BDA00028323462200000812
Figure BDA00028323462200000813
Figure BDA00028323462200000814
其中
Figure BDA00028323462200000815
Figure BDA0002832346220000091
Figure BDA0002832346220000092
Figure BDA0002832346220000093
Figure BDA0002832346220000094
Figure BDA0002832346220000095
Figure BDA0002832346220000096
Figure BDA0002832346220000097
均为临时变量。
(5c)tV的变化区间
Figure BDA0002832346220000098
和ηV的变化区间
Figure BDA0002832346220000099
在复数空间中构成了垂直极化分量透射系数TV的扇形区间,按下式计算:
Figure BDA00028323462200000910
Figure BDA00028323462200000911
Figure BDA00028323462200000912
Figure BDA00028323462200000913
其中
Figure BDA00028323462200000914
Figure BDA00028323462200000915
Figure BDA00028323462200000916
Figure BDA00028323462200000917
Figure BDA00028323462200000918
Figure BDA0002832346220000101
Figure BDA0002832346220000102
Figure BDA0002832346220000103
均为临时变量。
(5d)考虑步骤(4c1)中的远场方向图积分公式,可将透射系数的扇形区间表达进一步改写。水平极化分量透射系数在积分公式中形成的扇形区间的幅值tH′的变化区间记为
Figure BDA0002832346220000104
相位ηH′变化区间记为
Figure BDA0002832346220000105
垂直极化分量透射系数在积分公式中形成的扇形区间的幅值tV′的变化区间记为
Figure BDA0002832346220000106
相位ηV′变化区间记为
Figure BDA0002832346220000107
具体如下:
Figure BDA0002832346220000108
Figure BDA0002832346220000109
Figure BDA00028323462200001010
Figure BDA00028323462200001011
Figure BDA00028323462200001012
Figure BDA00028323462200001013
Figure BDA00028323462200001014
Figure BDA00028323462200001015
其中am=E(x,y)·ΔS,ph=k0sinθ(xcosφ+ysinφ),ΔS为划分的网格单元面积。
步骤六,将透射系数形成的扇形区间分段为多边形区间。
(6a)水平极化分量透射系数形成的扇形区间的四个端点
Figure BDA00028323462200001016
Figure BDA00028323462200001017
的坐标
Figure BDA00028323462200001018
Figure BDA00028323462200001019
Figure BDA00028323462200001020
Figure BDA00028323462200001021
Figure BDA00028323462200001022
如下:
Figure BDA00028323462200001023
Figure BDA0002832346220000111
Figure BDA0002832346220000112
Figure BDA0002832346220000113
其中
Figure BDA0002832346220000114
Figure BDA0002832346220000115
Figure BDA0002832346220000116
之间对应的是扇形的圆弧部分。
(6b)
Figure BDA0002832346220000117
Figure BDA0002832346220000118
之间、
Figure BDA0002832346220000119
Figure BDA00028323462200001110
之间、
Figure BDA00028323462200001111
Figure BDA00028323462200001112
之间用直线连接,
Figure BDA00028323462200001113
Figure BDA00028323462200001114
之间用包围圆弧、且与圆弧相切的切线段连接,将扇形区间的相位范围均分成10份,每一份对应于一段线段,从而将扇形区间用多边形区间表示,该多边形区间完全包络了原扇形区间,到此得到了水平极化分量透射系数形成的多边形区间。
(6c)垂直极化分量透射系数形成的扇形区间的四个端点
Figure BDA00028323462200001115
Figure BDA00028323462200001116
的坐标
Figure BDA00028323462200001117
Figure BDA00028323462200001118
Figure BDA00028323462200001119
Figure BDA00028323462200001120
Figure BDA00028323462200001121
如下:
Figure BDA00028323462200001122
Figure BDA00028323462200001123
Figure BDA00028323462200001124
Figure BDA00028323462200001125
其中
Figure BDA00028323462200001126
Figure BDA00028323462200001127
Figure BDA00028323462200001128
之间对应的是扇形的圆弧部分。
(6d)
Figure BDA00028323462200001129
Figure BDA00028323462200001130
之间、
Figure BDA00028323462200001131
Figure BDA00028323462200001132
之间、
Figure BDA00028323462200001133
Figure BDA00028323462200001134
之间用直线连接,
Figure BDA00028323462200001135
Figure BDA00028323462200001136
之间用包围圆弧、且与圆弧相切的切线段连接,将扇形区间的相位范围均分成10份,每一份对应于一段线段,从而将扇形区间用多边形区间表示,该多边形区间完全包络了原扇形区间,到此得到了垂直极化分量透射系数形成的多边形区间。
(6e)将上述水平极化分量透射系数和垂直极化分量透射系数形成的多边形区间的端点顺序调整,使得起始点在复数域具有最小的虚部,并且所有点按照环绕原扇形的逆时针顺序排列。
步骤七,利用多边形区间求和计算天线罩的远场方向图变化区间,并提取电性能指标变化区间。
(7a)对于步骤中得到的天线罩各处的水平极化分量和垂直极化分量透射系数的多边形区间,逐个进行多边形求和,得到远场方向图中各个离散点处的场值多边形区间。按照如下规则进行:
(7a1)对于任意的两组多边形端点组,分别记为pi,i=1,2,....,m和pj′,j=1,2,...,n,i和j的初始值为1,其合成的多边形端点记为qk,k=1,2,...,t,取s=1。
(7a2)qsx=pix+pjx′,qsy=piy+pjy′;
(7a3)计算从点pi指向点pi+1的矢量与横坐标正方向之间的夹角
Figure BDA0002832346220000121
并计算从点pj′指向pj+1′的矢量与横坐标正方向之间的夹角
Figure BDA0002832346220000122
Figure BDA0002832346220000123
则令i=i+1,否则令j=j+1;
(7a4)若i=m且j=n,则求和过程结束;否则,令s=s+1,转到步骤(5a2)。
(7b)根据步骤(7a)中得到的远场方向图中各个离散点处的场值多边形区间,计算多边形区间各个端点处的幅值,并并中提取其最大值,构成远场方向图的上确界,进而从中提取其最小值,构成远场方向图的下确界,从而得到远场方向图的变化区间,从中提取增益G2、主波束宽度B2和第一副瓣电平S2这些电性能指标的变化区间。
步骤八,判断天线罩的电性能指标是否满足预设要求。
对比步骤四与步骤七中得到的天线罩电性能指标变化区间,判断天线罩的电性能指标是否满足预设要求,如果满足,则天线罩电性能分析过程结束,否则,修改扇形圆弧部分的分段数,并重复步骤二至步骤八,直至得到电性能指标满足预设要求的分析结果。
本发明的优点可通过以下仿真实验进一步说明:
1.仿真参数
某球形天线罩,外形如图5所示,直径为1米,天线罩材料为玻璃钢材料,材料的相对介电常数为4,磁损耗角正切为0.015,罩壁厚度为8mm,壁厚误差为[-0.05mm,0.05mm],罩内天线口径为0.5米,工作频率为9.4GHz,其口径场为等幅同相分布,天线的扫描角取为0°,此时天线仰天,指向罩顶。
2.仿真内容与结果
利用本发明对上述天线罩进行厚度误差作用下的电性能分析,仿真结果如图6和图7所示所示,仿真数据如表1所示。
表1系统的电性能指标区间
Figure BDA0002832346220000131
从上述结果可见,采用不分组的蒙特卡洛方法,归一化增益区间与波束宽度区间上、下限相等(数据保留到小数点后2位),第一副瓣电平区间也较小,而采用本发明提出的分组蒙特卡洛方法后,这些电性能指标的区间均表现出了显著的变化,因而当天线罩厚度存在误差时,采用本方法可以在降低计算量的基础上,提供更为准确的电性能蒙特卡洛分析结果。
上述仿真数据实验证明,本发明可显著改善存在厚度误差时天线罩电性能预测结果的效率与准确度,从而对天线罩设计提供有效的依据。

Claims (7)

1.一种基于离散网格分组的天线罩电性能快速误差分析方法,其特征是:具体至少包括如下步骤:
步骤一,根据天线罩的结构形式建立天线罩的几何模型,设置网格边长为0.2λ,其中λ为天线的波长,对模型进行网格划分,记网格数目为N;
步骤二,在每个网格中选取一点,按照点的坐标将其归类到10个扇形区域中,实现网格分组;
步骤三,根据给定的天线罩厚度值变化范围,随机生成1000组厚度误差分布;
步骤四,计算厚度误差分布下的天线罩远场方向图,并提取电性能指标变化区间;
步骤五,根据厚度误差范围计算天线罩上各处的透射系数的幅值和相位的变化区间,得到透射系数在复数域形成的扇形区间;
步骤六,将透射系数形成的扇形区间分段为多边形区间;
步骤七,利用多边形区间求和计算天线罩的远场方向图变化区间,并提取电性能指标变化区间;
步骤八,对比步骤四与步骤七中得到的天线罩电性能指标变化区间,判断天线罩的电性能指标是否满足预设要求,如果满足,则天线罩电性能分析过程结束,否则,修改扇形圆弧部分的分段数,并重复步骤二至步骤八,直至得到电性能指标满足预设要求的分析结果。
2.根据权利要求1所述的一种基于离散网格分组的天线罩电性能快速误差分析方法,其特征是:所述的步骤二进一步包括如下步骤:
(2a)在划分的网格中选取一点,将该坐标转换为极坐标,极坐标的周向角度记为ψ;
(2b)沿周向角度进行均匀划分,当扇形区域数为M时,从0度开始,间隔一定步长选取M个分隔点,分隔点初始最小值为P1=0,最大值为PM=360,其间任一分隔点为:
Pi=P1+i*360/M,i=1,2,…,M
将M=10代入上式即可得到每一分隔点的值;
(2c)第i个分隔区域的区间范围为[Pi,Pi+1],划分10个扇形区域。
3.根据权利要求1所述的一种基于离散网格分组的天线罩电性能快速误差分析方法,其特征是:所述的步骤三,进一步包括如下步骤:
(3a)设天线罩的厚度误差范围为[dL,dU],则以均值为(dL+dU)/2、标准差为(dL+dU)/6生成10个随机数,对于步骤二中10个扇形区域,第一个扇形区域分配一个随机数,即为该区域中所有网格的厚度误差值。
(3b)重复步骤(3a)1000次,得到1000组对应于天线罩网格数目的厚度误差随机分布值,每一组厚度误差的随机分布值均为10个数,与划分的扇形区域数目相等。
4.根据权利要求1所述的一种基于离散网格分组的天线罩电性能快速误差分析方法,其特征是:所述的步骤四,进一步包括如下步骤:
(4a)沿罩体的高度,以天线罩的底面中心为原点、底面为xy平面建立一个坐标系O-xyz,罩体高度沿z方向,其中S为天线口径面,S′为透过天线罩后的口径面;
(4b)根据天线罩的结构参数和材料参数,用传输线理论计算罩体的透射系数TM,并根据已知的天线口径场E(x,y),计算透过天线罩后的口径场:E′(x,y)=E(x,y)·TM
(4b1)针对步骤(4a)建立的坐标系,将其x、y、z方向的分量分别用i、j、k表示,已知的天线口径场分布记为E(x,y);
(4b2)根据天线罩的几何形状与入射的口径场,求出天线罩上各点处的入射角α和极化角β,即将电磁波入射线与入射点处的法线夹角记为入射角α,将电磁波的极化方向与入射平面的夹角记为极化角β,其中入射平面由电磁波入射线与入射点处的法线构成;
(4b3)令天线罩材料的损耗角正切tanδ为0,根据天线罩各处的厚度d、相对介电常数εr,计算蒙皮上各点处的水平极化分量透射系数TH和垂直极化分量透射系数TV
Figure FDA0002832346210000031
Figure FDA0002832346210000032
其中,A=ch(jVd),
Figure FDA0002832346210000033
Figure FDA0002832346210000034
D=ch(jVd),
Figure FDA0002832346210000035
Figure FDA0002832346210000036
ZH=cosα,
Figure FDA0002832346210000037
这些参数均为中间变量;tH、tV分别为TH、TV的模值,ηH、ηV分别为TH、TV的相位;
(4b4)根据水平极化分量透射系数TH和垂直极化分量透射系数TV,得到主极化分量的透射系数:
Figure FDA0002832346210000041
其中,
Figure FDA0002832346210000042
为中间变量;
(4b5)将入射到天线罩上的口径场乘以其对应点处的透射系数,计算透过蒙皮后的口径场:E′(x,y)=E(x,y)·TM
(4c)根据透过天线罩后的口径场E′(x,y),计算带罩天线的远场F(θ,φ),绘制远场方向图,并从该远场方向图中提取增益G1、主波束位置B1和第一副瓣电平S1这些电性能指标,
(4c1)根据步骤(4b)中求得的透过天线罩后的口径场E′(x,y),按下式计算透过天线罩后的口径场产生的远场F(θ,φ):
Figure FDA0002832346210000043
其中,θ、φ是观察点P在直角坐标系O-xyz中的球坐标角,k0为自由空间传播常数,按公式
Figure FDA0002832346210000044
计算,λ是天线的波长,根据天线工作频率f和光速c,通过公式
Figure FDA0002832346210000045
计算得到;
(4c2)根据透过天线罩后的口径场产生的远场F(θ,φ)绘制加罩后的天线远场方向图,并从方向图中提取增益G1、主波束位置B1和第一副瓣电平S1这些电性能指标的变化区间。
5.根据权利要求1所述的一种基于离散网格分组的天线罩电性能快速误差分析方法,其特征是:所述的步骤五,进一步包括如下步骤:
(5a)天线罩水平极化分量透射系数TH和垂直极化分量透射系数TV可改写为下述表达式:
Figure FDA0002832346210000051
Figure FDA0002832346210000052
其中
Figure FDA0002832346210000053
Figure FDA0002832346210000054
Figure FDA0002832346210000055
均为临时变量,根据厚度d的变化区间[dL,dU],计算临时变量A1的变化区间
Figure FDA0002832346210000056
B1的变化区间
Figure FDA0002832346210000057
以及B2的变化区间
Figure FDA0002832346210000058
(5b)tH的变化区间
Figure FDA0002832346210000059
和ηH的变化区间
Figure FDA00028323462100000510
在复数空间中构成了水平极化分量透射系数TH的扇形区间,按下式计算:
Figure FDA00028323462100000511
Figure FDA00028323462100000512
Figure FDA00028323462100000513
Figure FDA00028323462100000514
其中
Figure FDA00028323462100000515
Figure FDA00028323462100000516
Figure FDA0002832346210000061
Figure FDA0002832346210000062
Figure FDA0002832346210000063
Figure FDA0002832346210000064
Figure FDA0002832346210000065
Figure FDA0002832346210000066
均为临时变量。
(5c)tV的变化区间
Figure FDA0002832346210000067
和ηV的变化区间
Figure FDA0002832346210000068
在复数空间中构成了垂直极化分量透射系数TV的扇形区间,按下式计算:
Figure FDA0002832346210000069
Figure FDA00028323462100000610
Figure FDA00028323462100000611
Figure FDA00028323462100000612
其中
Figure FDA00028323462100000613
Figure FDA00028323462100000614
Figure FDA00028323462100000615
Figure FDA00028323462100000616
Figure FDA00028323462100000617
Figure FDA00028323462100000618
Figure FDA0002832346210000071
Figure FDA0002832346210000072
均为临时变量;
(5d)考虑步骤(4c1)中的远场方向图积分公式,可将透射系数的扇形区间表达进一步改写。水平极化分量透射系数在积分公式中形成的扇形区间的幅值tH′的变化区间记为
Figure FDA0002832346210000073
相位ηH′变化区间记为
Figure FDA0002832346210000074
垂直极化分量透射系数在积分公式中形成的扇形区间的幅值tV′的变化区间记为
Figure FDA0002832346210000075
相位ηV′变化区间记为
Figure FDA0002832346210000076
具体如下:
Figure FDA0002832346210000077
Figure FDA0002832346210000078
Figure FDA0002832346210000079
Figure FDA00028323462100000710
Figure FDA00028323462100000711
Figure FDA00028323462100000712
Figure FDA00028323462100000713
Figure FDA00028323462100000714
其中am=E(x,y)·ΔS,ph=k0sinθ(xcosφ+ysinφ),ΔS为划分的网格单元面积。
6.根据权利要求1所述的一种基于离散网格分组的天线罩电性能快速误差分析方法,其特征是:所述的步骤六,进一步包括如下步骤:
(6a)水平极化分量透射系数形成的扇形区间的四个端点
Figure FDA00028323462100000715
Figure FDA00028323462100000716
的坐标
Figure FDA00028323462100000717
Figure FDA00028323462100000718
Figure FDA00028323462100000719
Figure FDA00028323462100000720
Figure FDA00028323462100000721
如下:
Figure FDA00028323462100000722
Figure FDA0002832346210000081
Figure FDA0002832346210000082
Figure FDA0002832346210000083
其中
Figure FDA0002832346210000084
Figure FDA0002832346210000085
Figure FDA0002832346210000086
之间对应的是扇形的圆弧部分;
(6b)
Figure FDA0002832346210000087
Figure FDA0002832346210000088
之间、
Figure FDA0002832346210000089
Figure FDA00028323462100000810
之间、
Figure FDA00028323462100000811
Figure FDA00028323462100000812
之间用直线连接,
Figure FDA00028323462100000813
Figure FDA00028323462100000814
之间用包围圆弧、且与圆弧相切的切线段连接,将扇形区间的相位范围均分成10份,每一份对应于一段线段,从而将扇形区间用多边形区间表示,该多边形区间完全包络了原扇形区间,到此得到了水平极化分量透射系数形成的多边形区间;
(6c)垂直极化分量透射系数形成的扇形区间的四个端点
Figure FDA00028323462100000815
Figure FDA00028323462100000816
的坐标
Figure FDA00028323462100000817
Figure FDA00028323462100000818
Figure FDA00028323462100000819
Figure FDA00028323462100000820
Figure FDA00028323462100000821
如下:
Figure FDA00028323462100000822
Figure FDA00028323462100000823
Figure FDA00028323462100000824
Figure FDA00028323462100000825
其中
Figure FDA00028323462100000826
Figure FDA00028323462100000827
Figure FDA00028323462100000828
之间对应的是扇形的圆弧部分。
(6d)
Figure FDA00028323462100000829
Figure FDA00028323462100000830
之间、
Figure FDA00028323462100000831
Figure FDA00028323462100000832
之间、
Figure FDA00028323462100000833
Figure FDA00028323462100000834
之间用直线连接,
Figure FDA00028323462100000835
Figure FDA00028323462100000836
之间用包围圆弧、且与圆弧相切的切线段连接,将扇形区间的相位范围均分成10份,每一份对应于一段线段,从而将扇形区间用多边形区间表示,该多边形区间完全包络了原扇形区间,到此得到了垂直极化分量透射系数形成的多边形区间
(6e)将上述水平极化分量透射系数和垂直极化分量透射系数形成的多边形区间的端点顺序调整,使得起始点在复数域具有最小的虚部,并且所有点按照环绕原扇形的逆时针顺序排列。
7.根据权利要求1所述的一种基于离散网格分组的天线罩电性能快速误差分析方法,其特征是:所述的步骤七,进一步包括如下步骤:
(7a)对于步骤中得到的天线罩各处的水平极化分量和垂直极化分量透射系数的多边形区间,逐个进行多边形求和,得到远场方向图中各个离散点处的场值多边形区间。按照如下规则进行:
(7a1)对于任意的两组多边形端点组,分别记为pi,i=1,2,....,m和pj′,j=1,2,...,n,i和j的初始值为1,其合成的多边形端点记为qk,k=1,2,...,t,取s=1;
(7a2)qsx=pix+pjx′,qsy=piy+pjy′;
(7a3)计算从点pi指向点pi+1的矢量与横坐标正方向之间的夹角
Figure FDA0002832346210000091
并计算从点pj′指向pj+1′的矢量与横坐标正方向之间的夹角
Figure FDA0002832346210000092
Figure FDA0002832346210000093
则令i=i+1,否则令j=j+1;
(7a4)若i=m且j=n,则求和过程结束;否则,令s=s+1,转到步骤(5a2);
(7b)根据步骤(7a)中得到的远场方向图中各个离散点处的场值多边形区间,计算多边形区间各个端点处的幅值,并并中提取其最大值,构成远场方向图的上确界,进而从中提取其最小值,构成远场方向图的下确界,从而得到远场方向图的变化区间,从中提取增益G2、主波束宽度B2和第一副瓣电平S2这些电性能指标的变化区间。
CN202011463547.7A 2020-12-11 2020-12-11 一种基于离散网格分组的天线罩电性能快速误差分析方法 Active CN112685877B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011463547.7A CN112685877B (zh) 2020-12-11 2020-12-11 一种基于离散网格分组的天线罩电性能快速误差分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011463547.7A CN112685877B (zh) 2020-12-11 2020-12-11 一种基于离散网格分组的天线罩电性能快速误差分析方法

Publications (2)

Publication Number Publication Date
CN112685877A true CN112685877A (zh) 2021-04-20
CN112685877B CN112685877B (zh) 2023-04-07

Family

ID=75449299

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011463547.7A Active CN112685877B (zh) 2020-12-11 2020-12-11 一种基于离散网格分组的天线罩电性能快速误差分析方法

Country Status (1)

Country Link
CN (1) CN112685877B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113655297A (zh) * 2021-07-30 2021-11-16 中国航空工业集团公司济南特种结构研究所 一种卫通天线罩单程功率传输效率的校正方法
CN114329815A (zh) * 2021-11-29 2022-04-12 西安电子科技大学 考虑函数单调性的天线罩电性能区间分析方法
CN116484522A (zh) * 2023-04-13 2023-07-25 大连理工大学 一种基于级联网络的天线罩电性能补偿加工量确定方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BE701766A (zh) * 1966-07-28 1968-01-24
CN102590656A (zh) * 2012-01-03 2012-07-18 西安电子科技大学 基于远场的天线罩电性能预测方法
CN106654566A (zh) * 2017-01-05 2017-05-10 西安电子科技大学 一种飞行器天线罩的快速厚度设计方法
WO2018121395A1 (zh) * 2016-12-30 2018-07-05 深圳光启高等理工研究院 天线罩的电性能评估方法及装置
CN108920831A (zh) * 2018-06-29 2018-11-30 西安电子科技大学 高速飞行器罩体高温烧蚀对天线电性能影响的快速计算方法
CN111177890A (zh) * 2019-12-11 2020-05-19 西安电子科技大学 基于多边形区间的天线罩电性能不确定性分析方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BE701766A (zh) * 1966-07-28 1968-01-24
CN102590656A (zh) * 2012-01-03 2012-07-18 西安电子科技大学 基于远场的天线罩电性能预测方法
WO2018121395A1 (zh) * 2016-12-30 2018-07-05 深圳光启高等理工研究院 天线罩的电性能评估方法及装置
CN106654566A (zh) * 2017-01-05 2017-05-10 西安电子科技大学 一种飞行器天线罩的快速厚度设计方法
CN108920831A (zh) * 2018-06-29 2018-11-30 西安电子科技大学 高速飞行器罩体高温烧蚀对天线电性能影响的快速计算方法
CN111177890A (zh) * 2019-12-11 2020-05-19 西安电子科技大学 基于多边形区间的天线罩电性能不确定性分析方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PENG LI等: "An Inversion Design Method for the Radome Thickness Based on Interval Arithmetic", 《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》 *
PENG LI等: "Power Pattern Tolerance Analysis of Radome With the Material Property Error Based on Interval Arithmetic", 《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》 *
高文德: "曲面FSS天线罩成型工艺及其误差对电磁特性的影响研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113655297A (zh) * 2021-07-30 2021-11-16 中国航空工业集团公司济南特种结构研究所 一种卫通天线罩单程功率传输效率的校正方法
CN113655297B (zh) * 2021-07-30 2024-04-09 中国航空工业集团公司济南特种结构研究所 一种卫通天线罩单程功率传输效率的校正方法
CN114329815A (zh) * 2021-11-29 2022-04-12 西安电子科技大学 考虑函数单调性的天线罩电性能区间分析方法
CN114329815B (zh) * 2021-11-29 2022-09-16 西安电子科技大学 考虑函数单调性的天线罩电性能区间分析方法
CN116484522A (zh) * 2023-04-13 2023-07-25 大连理工大学 一种基于级联网络的天线罩电性能补偿加工量确定方法
CN116484522B (zh) * 2023-04-13 2024-04-19 大连理工大学 一种基于级联网络的天线罩电性能补偿加工量确定方法

Also Published As

Publication number Publication date
CN112685877B (zh) 2023-04-07

Similar Documents

Publication Publication Date Title
CN112685877B (zh) 一种基于离散网格分组的天线罩电性能快速误差分析方法
CN111177890B (zh) 基于多边形区间的天线罩电性能不确定性分析方法
CN106654566B (zh) 一种飞行器天线罩的快速厚度设计方法
CN106469850B (zh) 一种天线罩的厚度设计方法
CN104036114B (zh) 一种基于机电耦合的六边形有源相控阵天线结构公差的快速确定方法
CN109472066B (zh) 基于单元中心点位移的反射面天线随机误差分析方法
CN102708257A (zh) 一种有源相控阵天线结构公差的快速确定方法
CN105718662B (zh) 基于机电耦合的圆柱共形阵列天线结构公差快速确定方法
CN108984902B (zh) 基于相位调整的飞行器天线罩快速厚度设计方法
CN114329815B (zh) 考虑函数单调性的天线罩电性能区间分析方法
CN106096209B (zh) 一种面向最低成本的反射面天线分块面板尺寸确定方法
CN114169201B (zh) 面向电性能的反射面天线结构加权优化方法
CN111339606A (zh) 一种基于绕射原理的机翼遮挡效应的计算方法
CN114491992B (zh) 一种基于等效偶极矩和物理光学法的高效电磁散射方法
CN116432376A (zh) 基于插值的天线罩电性能区间分析方法
CN111291493B (zh) 一种用于机载预警共形阵俯仰探测波束赋形的设计方法
CN114386271A (zh) 一种考虑互耦效应的随机阵列天线方向图综合方法
CN107622148A (zh) 基于机电耦合的圆柱共形阵列天线结构最佳弹性模量的确定方法
CN117744480A (zh) 飞行器天线罩的稳健变厚度设计方法
CN109472089B (zh) 一种适用于被动毫米波人体安检仪的布局结构及方法
Nayeri et al. Design, full-wave analysis, and near-field diagnostics of reflectarray antennas
Parini et al. Optimizing a CATR quiet zone using an array feed
CN113625062B (zh) 基于泰勒展开法的天线罩电性能估计方法
CN113203897B (zh) 一种任意二维天线阵列的旁瓣抑制与波束形成方法
Vedaee et al. A New Analytical Redesign of a Double-Curvature Reflector Antenna Using Invasive Weed Optimization (IWO) Algorithm

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