CN107240780B - 基于面片积分公式的伞状天线结构优化设计方法 - Google Patents

基于面片积分公式的伞状天线结构优化设计方法 Download PDF

Info

Publication number
CN107240780B
CN107240780B CN201710388179.6A CN201710388179A CN107240780B CN 107240780 B CN107240780 B CN 107240780B CN 201710388179 A CN201710388179 A CN 201710388179A CN 107240780 B CN107240780 B CN 107240780B
Authority
CN
China
Prior art keywords
umbrella
rib
antenna
triangle
calculating
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
CN201710388179.6A
Other languages
English (en)
Other versions
CN107240780A (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 CN201710388179.6A priority Critical patent/CN107240780B/zh
Publication of CN107240780A publication Critical patent/CN107240780A/zh
Application granted granted Critical
Publication of CN107240780B publication Critical patent/CN107240780B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01QANTENNAS, i.e. RADIO AERIALS
    • H01Q15/00Devices for reflection, refraction, diffraction or polarisation of waves radiated from an antenna, e.g. quasi-optical devices
    • H01Q15/14Reflecting surfaces; Equivalent structures
    • H01Q15/16Reflecting surfaces; Equivalent structures curved in two dimensions, e.g. paraboloidal
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01QANTENNAS, i.e. RADIO AERIALS
    • H01Q15/00Devices for reflection, refraction, diffraction or polarisation of waves radiated from an antenna, e.g. quasi-optical devices
    • H01Q15/14Reflecting surfaces; Equivalent structures
    • H01Q15/141Apparatus or processes specially adapted for manufacturing reflecting surfaces

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Aerials With Secondary Devices (AREA)

Abstract

本发明公开了一种基于面片积分公式的伞状天线结构优化设计方法,具体步骤包括:输入伞状天线结构参数与电参数;计算伞状天线最优焦距;计算天线肋的分段数;计算肋上点的坐标;计算相邻肋间点的坐标;生成伞状天线所有节点坐标;计算面片平面方程系数;计算三角形边方程系数;计算三角形投影面积;计算三角形面片轴向误差;输出伞状天线轴向精度;判断轴向精度是否满足要求;输出伞状天线结构参数;更新伞状天线结构参数。本发明考虑了伞状天线采用面片拼合而成的特点,基于面片积分公式对伞状天线进行轴向精度分析,可指导伞状天线的机械结构设计与机电集成优化设计。

Description

基于面片积分公式的伞状天线结构优化设计方法
技术领域
本发明属于雷达天线技术领域,具体涉及雷达天线领域中的一种基于面片积分公式的伞状天线结构优化设计方法。
背景技术
作为星载可展开天线的一种,伞状天线是一类较早开始研究并投入应用的星载可展开天线。利用柔性金属丝网与刚性肋的特点设计而成的伞状天线正被逐渐被应用于高增益、轻质量的星载天线设计中。伞状天线采用柔性金属丝网构成天线反射面,不可避免地引入面片拼合误差,即伞状天线概念阶段需要进行考虑的原理误差。与此同时,伞状天线在轨受到热冲击、姿态变轨等外部载荷影响,其反射面形状也会发生变形,导致电性能恶化。如何有效准确计算伞状天线的表面误差是进行伞状天线详细机械结构设计与机电集成优化设计的前提。
李小平和徐德红在论文“网状可展开天线两种网面成形方式分析”(电子机械工程,2010年第26卷第1期,38-40)总结了网状可展开天线两种网面成形方式,给出了伞状天线的表面均方根值误差计算公式;但该公式仅仅考虑了伞状天线的原理误差,对于在任一工况下的表面变形难以给出通用的计算过程。王从思等和冷国俊等分别在论文“天线表面误差的精确计算方法及电性能分析”(电波科学学报,2006年第26卷第3期,403-409)和“天线最佳吻合轴向误差的精确计算方法”(电波科学学报,2009年第24卷第5期,826-831)中提出了针对地基圆抛物面天线的表面误差计算方法与改进方法,其均采用的是节点误差进行叠加的方式获得天线表面误差均方根值,无法适应于具有较大面片拼合而成的伞状天线上。丁波在“星载网状可展开天线结构分析与优化设计”(西安电子科技大学2016年硕士学位论文)中提出了采用面积坐标进行误差计算的思路,但该方法以理想反射面天线为基准,没有考虑到伞状天线焦距发生变化的特性。因此,针对伞状天线的结构特点,考虑伞状天线面片拼合与焦距变化的特殊性,采用基于面片积分公式的方法对伞状天线进行表面误差计算,以期准确地获得伞状天线在任一工况下的表面误差,并以此开展伞状天线机械结构设计与机电集成优化设计。
发明内容
本发明的目的是克服上述现有技术的不足,提供一种基于面片积分公式的伞状天线结构优化设计方法。该方法基于面片积分公式,考虑到伞状天线采用面片拼合而成的特点,以伞状天线变化后的焦距为基准,通过积分运算的方法获得伞状天线轴向精度,可指导伞状天线的机械结构设计与机电集成优化设计。
本发明的技术方案是:基于面片积分公式的伞状天线结构优化设计方法,包括如下步骤:
(1)输入伞状天线结构参数与电参数
输入用户提供的伞状天线的结构参数与电参数;其中结构参数包括伞状天线口径、焦距、偏置距离、肋数和轴向精度设计要求;电参数包括工作波长;
(2)计算伞状天线最优焦距
根据用户提供的天线结构参数,按照下式计算伞状天线最优焦距
Figure GDA0002228122600000031
其中,fs表示伞状天线最优焦距,下标s表示区别于理想天线的伞状天线,f表示用户输入的伞状天线结构参数中的焦距,π表示圆周率,N表示肋数;
(3)根据用户提供的天线结构参数与电参数计算天线肋的分段数;
(4)根据用户提供的天线结构参数与肋的分段数,计算肋上点的坐标;
(5)计算相邻肋间点的坐标
根据相邻肋构成抛物柱面的特性,结合肋上点坐标计算相邻肋间点的坐标;根据伞状天线圆形口径的闭合特性,计算第N根肋与第1根肋构成的肋间点的坐标;
(6)生成伞状天线所有节点坐标
将计算得到的肋上点、相邻肋间点的坐标与原点坐标合并在一起得到伞状天线所有节点坐标;
(7)计算面片平面方程系数
根据伞状天线所有节点信息,按照下式计算面片平面方程系数
Figure GDA0002228122600000032
Figure GDA0002228122600000033
Figure GDA0002228122600000034
其中,A、B、C分别为面片平面方程的三个系数,a、b、c分别为三角形面片顶点,xa、ya、za表示编号为a的三角形顶点直角坐标,xb、yb、zb表示编号为b的三角形顶点直角坐标,xc、yc、zc表示编号为c的三角形顶点直角坐标;
(8)计算三角形边方程系数
根据伞状天线节点坐标,按照下式计算三角形三条边方程的系数
Figure GDA0002228122600000041
Figure GDA0002228122600000042
Figure GDA0002228122600000043
其中,a、b、c分别为三角形面片顶点,Kab、Lab分别为以三角形顶点a、b构成的边方程的两个系数,下标ab表示三角形顶点ab构成的边,Kac、Lac分别为以三角形顶点a、c构成的边方程的两个系数,下标ac表示三角形顶点a、c构成的边,Kbc、Lbc分别为以三角形顶点b、c构成的边方程的两个系数,下标bc表示三角形顶点b、c构成的边,xa、ya、za表示编号为a的三角形顶点直角坐标,xb、yb、zb表示编号为b的三角形顶点直角坐标,xc、yc、zc表示编号为c的三角形顶点直角坐标;
(9)计算三角形投影面积
根据伞状天线节点坐标与三角形边方程系数,按照下式计算三角形投影面积
Figure GDA0002228122600000044
其中,S表示三角形投影面积,xa、xb、xc表示三角形编号为a、b、c的三个顶点在x方向上的直角坐标,Kab、Lab分别为以三角形顶点a、b构成的边方程的两个系数,下标ab表示三角形顶点ab构成的边,Kac、Lac分别为以三角形顶点a、c构成的边方程的两个系数,下标ac表示三角形顶点a、c构成的边,Kbc、Lbc分别为以三角形顶点b、c构成的边方程的两个系数,下标bc表示三角形顶点b、c构成的边,dydx表示在三角形投影面上进行积分运算,积分变量分别为y、x分量;
(10)计算三角形面片轴向误差
10a)根据伞状天线结构参数、面片平面方程系数和最优焦距,计算面片内部点轴向误差
Δ=A·x+B·y+C-(x2+y2)/4f-(f-fs)
其中,Δ表示面片内部点轴向误差,A、B、C分别为面片平面方程的三个系数,f为伞状天线结构参数中的焦距,fs为伞状天线最优焦距,下标s表示区别于理想天线的伞状天线,x、y表示节点直角坐标;
10b)根据面片内部点轴向误差,计算三角形面片的轴向误差
Ω=∫∫Δ2dydx
其中,Ω表示三角形面片的轴向误差平方值,Δ表示面片内部点轴向误差,dydx表示在三角形投影面上进行积分运算,积分变量分别为y、x分量;
(11)输出伞状天线轴向精度
根据三角形面片轴向误差与三角形投影面积,按照下式计算伞状天线轴向精度
Figure GDA0002228122600000051
其中,δ表示伞状天线轴向精度,Ωi表示第i个三角形面片的轴向误差平方值,Si表示第i个三角形投影面积,n表示三角形面片数目,Σ表示加和符号;
(12)判断轴向精度是否满足要求
判断伞状天线轴向精度是否满足轴向精度设计要求,如果满足要求则转至步骤(13),否则转至步骤(14);
(13)输出伞状天线结构参数
当伞状天线轴向精度满足轴向精度设计要求时,输出伞状天线结构参数;
(14)更新伞状天线结构参数
当伞状天线轴向精度不满足轴向精度设计要求时,更新伞状天线结构参数,转至步骤(1)。
步骤(3)中天线肋的分段数按照下式进行选择计算:
Figure GDA0002228122600000061
其中,λ为工作波长,D为伞状天线口径,m为天线肋的分段数,m取满足上式条件的整数。
步骤(4)中根据用户提供的天线结构参数与肋的分段数,按照下式计算肋上点的坐标:
Figure GDA0002228122600000062
Figure GDA0002228122600000063
Figure GDA0002228122600000064
其中,xi,j、yi,j、zi,j分别为肋上点的x向坐标、y向坐标与z向坐标,下标i表示肋编号,下标j表示所在肋上点编号,D表示伞状天线口径,m表示天线肋的分段数,π表示圆周率,N表示肋数,f表示伞状天线焦距,肋编号i的取值范围为从1到N,肋上点编号j的取值范围为从1到m。
步骤(5)中:
5a)根据相邻肋构成抛物柱面的特性,结合肋上点坐标按照下式计算相邻肋间点的坐标:
Figure GDA0002228122600000071
Figure GDA0002228122600000072
Figure GDA0002228122600000073
其中,xi,j,k、yi,j,k、zi,j,k分别为相邻肋间点的x向坐标、y向坐标与z向坐标,下标i表示肋编号,下标j表示所在肋上点编号,下标k表示相邻肋间点在对应肋上点之间的编号,肋编号i的取值范围为从1到N-1,肋上点编号j的取值范围为从2到m,m表示天线肋的分段数,相邻肋间点在对应肋上点之间的编号k的取值范围为从1到j-1,xi,j、yi,j、zi,j分别表示第i根肋上第j个肋上点的x向坐标、y向坐标与z向坐标,xi+1,j、yi+1,j、zi+1,j分别表示与第i根肋相邻的第i+1根肋上第j个肋上点的x向坐标、y向坐标与z向坐标;
5b)根据伞状天线圆形口径的闭合特性,按照下式计算第N根肋与第1根肋构成的肋间点的坐标:
Figure GDA0002228122600000074
Figure GDA0002228122600000075
Figure GDA0002228122600000076
其中,xN,j,k、yN,j,k、zN,j,k分别为第N根肋与第1根肋构成的肋间点的x向坐标、y向坐标与z向坐标,下标N表示第N根肋编号,下标j表示第N根肋上点编号,下标k表示第N根肋与第1根肋构成的肋间点在对应肋上点之间的编号,肋上点编号j的取值范围为从2到m,m表示天线肋的分段数,第N根肋与第1根肋构成的肋间点在对应肋上点之间的编号k的取值范围为从1到j-1,xN,j、yN,j、zN,j分别表示第N根肋上第j个肋上点的x向坐标、y向坐标与z向坐标,x1,j、y1,j、z1,j分别表示第1根肋上第j个肋上点的x向坐标、y向坐标与z向坐标。
本发明的有益效果:本发明首先输入伞状天线结构参数与电参数,根据结构参数与电参数信息计算伞状天线最优焦距与天线肋的分段数;其次,依次计算肋上点的坐标与相邻肋间点的坐标,并生成伞状天线所有节点坐标;再次,根据节点坐标计算面片平面方程系数与边方程系数,以此计算三角形投影面积;然后,结合伞状天线最优焦距计算三角形面片轴向误差,并输出伞状天线轴向精度;最后,判断轴向精度是否满足设计要求,如果满足要求则输出伞状天线结构参数,否则更新天线结构参数,实现伞状天线结构优化设计。
与现有技术相比,本发明具有以下优点:
1.本发明基于面片积分公式,考虑到伞状天线采用面片拼合而成的特点,以伞状天线最优焦距为基准,通过积分运算获得伞状天线轴向精度;
2.本发明与之前进行轴向精度分析的方法相比,既考虑了伞状天线面片拼合的特点,又考虑了伞状天线由于肋带来的焦距变化,还可以对任一工况下的表面变形进行计算,具有较强的通用性。
以下将结合附图对本发明做进一步详细说明。
附图说明
图1为本发明的流程图;
图2为伞状天线结构示意图;
图3为伞状天线投影示意图。
具体实施方式
下面结合附图1,对本发明具体实施方式作进一步的详细描述:
本发明提供了一种基于面片积分公式的伞状天线结构优化设计方法,包括如下步骤:
步骤1,输入用户提供的伞状天线的结构参数与电参数;其中结构参数包括伞状天线口径、焦距、偏置距离、肋数和轴向精度设计要求;电参数包括工作波长;
步骤2,根据用户提供的天线结构参数,按照下式计算伞状天线最优焦距
Figure GDA0002228122600000091
其中,fs表示伞状天线最优焦距,下标s表示区别于理想天线的伞状天线,f表示用户输入的伞状天线结构参数中的焦距,π表示圆周率,N表示肋数;
步骤3,根据用户提供的天线结构参数与电参数计算天线肋的分段数,其中分段数按照下式进行选择计算
Figure GDA0002228122600000092
其中,λ为工作波长,D为伞状天线口径,m为天线肋的分段数,m取满足上式条件的整数;
步骤4,根据用户提供的天线结构参数与肋的分段数,按照下式计算肋上点的坐标
Figure GDA0002228122600000101
Figure GDA0002228122600000102
Figure GDA0002228122600000103
其中,xi,j、yi,j、zi,j分别为肋上点的x向坐标、y向坐标与z向坐标,下标i表示肋编号,下标j表示所在肋上点编号,D表示伞状天线口径,m表示天线肋的分段数,π表示圆周率,N表示肋数,f表示伞状天线焦距,肋编号i的取值范围为从1到N,肋上点编号j的取值范围为从1到m;
步骤5,计算相邻肋间点的坐标
5a)根据相邻肋构成抛物柱面的特性,结合肋上点坐标按照下式计算相邻肋间点的坐标
Figure GDA0002228122600000104
Figure GDA0002228122600000105
Figure GDA0002228122600000106
其中,xi,j,k、yi,j,k、zi,j,k分别为相邻肋间点的x向坐标、y向坐标与z向坐标,下标i表示肋编号,下标j表示所在肋上点编号,下标k表示相邻肋间点在对应肋上点之间的编号,肋编号i的取值范围为从1到N-1,肋上点编号j的取值范围为从2到m,m表示天线肋的分段数,相邻肋间点在对应肋上点之间的编号k的取值范围为从1到j-1,xi,j、yi,j、zi,j分别表示第i根肋上第j个肋上点的x向坐标、y向坐标与z向坐标,xi+1,j、yi+1,j、zi+1,j分别表示与第i根肋相邻的第i+1根肋上第j个肋上点的x向坐标、y向坐标与z向坐标;
5b)根据伞状天线圆形口径的闭合特性,按照下式计算第N根肋与第1根肋构成的肋间点的坐标
Figure GDA0002228122600000111
Figure GDA0002228122600000112
Figure GDA0002228122600000113
其中,xN,j,k、yN,j,k、zN,j,k分别为第N根肋与第1根肋构成的肋间点的x向坐标、y向坐标与z向坐标,下标N表示第N根肋编号,下标j表示第N根肋上点编号,下标k表示第N根肋与第1根肋构成的肋间点在对应肋上点之间的编号,肋上点编号j的取值范围为从2到m,m表示天线肋的分段数,第N根肋与第1根肋构成的肋间点在对应肋上点之间的编号k的取值范围为从1到j-1,xN,j、yN,j、zN,j分别表示第N根肋上第j个肋上点的x向坐标、y向坐标与z向坐标,x1,j、y1,j、z1,j分别表示第1根肋上第j个肋上点的x向坐标、y向坐标与z向坐标;
步骤6,将计算得到的肋上点、相邻肋间点的坐标与原点坐标合并在一起得到伞状天线所有节点坐标;
步骤7,根据伞状天线所有节点信息,按照下式计算面片平面方程系数
Figure GDA0002228122600000114
Figure GDA0002228122600000115
Figure GDA0002228122600000121
其中,A、B、C分别为面片平面方程的三个系数,a、b、c分别为三角形面片顶点,xa、ya、za表示编号为a的三角形顶点直角坐标,xb、yb、zb表示编号为b的三角形顶点直角坐标,xc、yc、zc表示编号为c的三角形顶点直角坐标;
步骤8,根据伞状天线节点坐标,按照下式计算三角形三条边方程的系数
Figure GDA0002228122600000122
Figure GDA0002228122600000123
Figure GDA0002228122600000124
其中,a、b、c分别为三角形面片顶点,Kab、Lab分别为以三角形顶点a、b构成的边方程的两个系数,下标ab表示三角形顶点ab构成的边,Kac、Lac分别为以三角形顶点a、c构成的边方程的两个系数,下标ac表示三角形顶点a、c构成的边,Kbc、Lbc分别为以三角形顶点b、c构成的边方程的两个系数,下标bc表示三角形顶点b、c构成的边,xa、ya、za表示编号为a的三角形顶点直角坐标,xb、yb、zb表示编号为b的三角形顶点直角坐标,xc、yc、zc表示编号为c的三角形顶点直角坐标;
步骤9,根据伞状天线节点坐标与三角形边方程系数,按照下式计算三角形投影面积
Figure GDA0002228122600000125
其中,S表示三角形投影面积,xa、xb、xc表示三角形编号为a、b、c的三个顶点在x方向上的直角坐标,Kab、Lab分别为以三角形顶点a、b构成的边方程的两个系数,下标ab表示三角形顶点ab构成的边,Kac、Lac分别为以三角形顶点a、c构成的边方程的两个系数,下标ac表示三角形顶点a、c构成的边,Kbc、Lbc分别为以三角形顶点b、c构成的边方程的两个系数,下标bc表示三角形顶点b、c构成的边,dydx表示在三角形投影面上进行积分运算,积分变量分别为y、x分量;
步骤10,计算三角形面片轴向误差
10a)根据伞状天线结构参数、面片平面方程系数和最优焦距,计算面片内部点轴向误差
Δ=A·x+B·y+C-(x2+y2)/4f-(f-fs)
其中,Δ表示面片内部点轴向误差,A、B、C分别为面片平面方程的三个系数,f为伞状天线结构参数中的焦距,fs为伞状天线最优焦距,下标s表示区别于理想天线的伞状天线,x、y表示节点直角坐标;
10b)根据面片内部点轴向误差,计算三角形面片的轴向误差
Ω=∫∫Δ2dydx
其中,Ω表示三角形面片的轴向误差平方值,Δ表示面片内部点轴向误差,dydx表示在三角形投影面上进行积分运算,积分变量分别为y、x分量;
步骤11,根据三角形面片轴向误差与三角形投影面积,按照下式计算伞状天线轴向精度
Figure GDA0002228122600000131
其中,δ表示伞状天线轴向精度,Ωi表示第i个三角形面片的轴向误差平方值,Si表示第i个三角形投影面积,n表示三角形面片数目,Σ表示加和符号;
步骤12,判断伞状天线轴向精度是否满足轴向精度设计要求,如果满足要求则转至步骤13,否则转至步骤14;
步骤13,当伞状天线轴向精度满足轴向精度设计要求时,输出伞状天线结构参数;
步骤14,当伞状天线轴向精度不满足轴向精度设计要求时,更新伞状天线结构参数,转至步骤1。
本发明的优点可通过以下仿真实验进一步说明:
1.仿真条件:
伞状天线口径10m,焦距10m,偏置距离0,肋数为18。
伞状天线结构示意图如图2所示,伞状天线投影示意图如图3所示。
2.仿真结果:
采用本发明的方法进行基于面片积分公式的伞状天线轴向精度,并输出伞状天线轴向精度。采用本方法得到的伞状天线轴向精度为6.64mm。
综上所述,本发明首先输入伞状天线结构参数与电参数,根据结构参数与电参数信息计算伞状天线最优焦距与天线肋的分段数;其次,依次计算肋上点的坐标与相邻肋间点的坐标,并生成伞状天线所有节点坐标;再次,根据节点坐标计算面片平面方程系数与边方程系数,以此计算三角形投影面积;然后,结合伞状天线最优焦距计算三角形面片轴向误差,并输出伞状天线轴向精度;最后,判断轴向精度是否满足设计要求,如果满足要求则输出伞状天线结构参数,否则更新天线结构参数,实现伞状天线结构优化设计。
与现有技术相比,本发明具有以下优点:
1.本发明基于面片积分公式,考虑到伞状天线采用面片拼合而成的特点,以伞状天线最优焦距为基准,通过积分运算获得伞状天线轴向精度;
2.本发明与之前进行轴向精度分析的方法相比,既考虑了伞状天线面片拼合的特点,又考虑了伞状天线由于肋带来的焦距变化,还可以对任一工况下的表面变形进行计算,具有较强的通用性。
本实施方式中没有详细叙述的部分属本行业的公知的常用手段,这里不一一叙述。以上例举仅仅是对本发明的举例说明,并不构成对本发明的保护范围的限制,凡是与本发明相同或相似的设计均属于本发明的保护范围之内。

Claims (4)

1.基于面片积分公式的伞状天线结构优化设计方法,其特征在于,包括如下步骤:
(1)输入伞状天线结构参数与电参数
输入用户提供的伞状天线的结构参数与电参数;其中结构参数包括伞状天线口径、焦距、偏置距离、肋数和轴向精度设计要求;电参数包括工作波长;
(2)计算伞状天线最优焦距
根据用户提供的天线结构参数,按照下式计算伞状天线最优焦距
Figure FDA0002228122590000011
其中,fs表示伞状天线最优焦距,下标s表示区别于理想天线的伞状天线,f表示用户输入的伞状天线结构参数中的焦距,π表示圆周率,N表示肋数;
(3)根据用户提供的天线结构参数与电参数计算天线的肋的分段数;
(4)根据用户提供的天线结构参数与肋的分段数,计算肋上点的坐标;
(5)计算相邻肋间点的坐标
根据相邻肋构成抛物柱面的特性,结合肋上点坐标计算相邻肋间点的坐标;根据伞状天线圆形口径的闭合特性,计算第N根肋与第1根肋构成的肋间点的坐标;
(6)生成伞状天线所有节点坐标
将计算得到的肋上点、相邻肋间点的坐标与原点坐标合并在一起得到伞状天线所有节点坐标;
(7)计算面片平面方程系数
根据伞状天线所有节点信息,按照下式计算面片平面方程系数
Figure FDA0002228122590000021
Figure FDA0002228122590000022
Figure FDA0002228122590000023
其中,A、B、C分别为面片平面方程的三个系数,a、b、c分别为三角形面片顶点,xa、ya、za表示编号为a的三角形顶点直角坐标,xb、yb、zb表示编号为b的三角形顶点直角坐标,xc、yc、zc表示编号为c的三角形顶点直角坐标;
(8)计算三角形边方程系数
根据伞状天线节点坐标,按照下式计算三角形三条边方程的系数
Figure FDA0002228122590000024
Figure FDA0002228122590000025
Figure FDA0002228122590000026
其中,a、b、c分别为三角形面片顶点,Kab、Lab分别为以三角形顶点a、b构成的边方程的两个系数,下标ab表示三角形顶点ab构成的边,Kac、Lac分别为以三角形顶点a、c构成的边方程的两个系数,下标ac表示三角形顶点a、c构成的边,Kbc、Lbc分别为以三角形顶点b、c构成的边方程的两个系数,下标bc表示三角形顶点b、c构成的边,xa、ya、za表示编号为a的三角形顶点直角坐标,xb、yb、zb表示编号为b的三角形顶点直角坐标,xc、yc、zc表示编号为c的三角形顶点直角坐标;
(9)计算三角形投影面积
根据伞状天线节点坐标与三角形边方程系数,按照下式计算三角形投影面积
Figure FDA0002228122590000031
其中,S表示三角形投影面积,xa、xb、xc表示三角形编号为a、b、c的三个顶点在x方向上的直角坐标,Kab、Lab分别为以三角形顶点a、b构成的边方程的两个系数,下标ab表示三角形顶点ab构成的边,Kac、Lac分别为以三角形顶点a、c构成的边方程的两个系数,下标ac表示三角形顶点a、c构成的边,Kbc、Lbc分别为以三角形顶点b、c构成的边方程的两个系数,下标bc表示三角形顶点b、c构成的边,dydx表示在三角形投影面上进行积分运算,积分变量分别为y、x分量;
(10)计算三角形面片轴向误差
10a)根据伞状天线结构参数、面片平面方程系数和最优焦距,计算面片内部点轴向误差
Δ=A·x+B·y+C-(x2+y2)/4f-(f-fs)
其中,Δ表示面片内部点轴向误差,A、B、C分别为面片平面方程的三个系数,f为伞状天线结构参数中的焦距,fs为伞状天线最优焦距,下标s表示区别于理想天线的伞状天线,x、y表示节点直角坐标;
10b)根据面片内部点轴向误差,计算三角形面片的轴向误差
Ω=∫∫Δ2dydx
其中,Ω表示三角形面片的轴向误差平方值,Δ表示面片内部点轴向误差,dydx表示在三角形投影面上进行积分运算,积分变量分别为y、x分量;
(11)输出伞状天线轴向精度
根据三角形面片轴向误差与三角形投影面积,按照下式计算伞状天线轴向精度
Figure FDA0002228122590000041
其中,δ表示伞状天线轴向精度,Ωi表示第i个三角形面片的轴向误差平方值,Si表示第i个三角形投影面积,n表示三角形面片数目,Σ表示加和符号;
(12)判断轴向精度是否满足要求
判断伞状天线轴向精度是否满足轴向精度设计要求,如果满足要求则转至步骤(13),否则转至步骤(14);
(13)输出伞状天线结构参数
当伞状天线轴向精度满足轴向精度设计要求时,输出伞状天线结构参数;
(14)更新伞状天线结构参数
当伞状天线轴向精度不满足轴向精度设计要求时,更新伞状天线结构参数,转至步骤(1)。
2.如权利要求1所述的基于面片积分公式的伞状天线结构优化设计方法,其特征在于,步骤(3)中天线的肋的分段数按照下式进行选择计算:
Figure FDA0002228122590000042
其中,λ为工作波长,D为伞状天线口径,m为天线肋的分段数,m取满足上式条件的整数。
3.如权利要求1所述的基于面片积分公式的伞状天线结构优化设计方法,其特征在于,步骤(4)中根据用户提供的天线结构参数与肋的分段数,按照下式计算肋上点的坐标:
Figure FDA0002228122590000051
Figure FDA0002228122590000052
Figure FDA0002228122590000053
其中,xi,j、yi,j、zi,j分别为肋上点的x向坐标、y向坐标与z向坐标,下标i表示肋编号,下标j表示所在肋上点编号,D表示伞状天线口径,m表示天线肋的分段数,π表示圆周率,N表示肋数,f表示伞状天线焦距,肋编号i的取值范围为从1到N,肋上点编号j的取值范围为从1到m。
4.如权利要求1所述的基于面片积分公式的伞状天线结构优化设计方法,其特征在于,步骤(5)中:
5a)根据相邻肋构成抛物柱面的特性,结合肋上点坐标按照下式计算相邻肋间点的坐标:
Figure FDA0002228122590000054
Figure FDA0002228122590000055
Figure FDA0002228122590000056
其中,xi,j,k、yi,j,k、zi,j,k分别为相邻肋间点的x向坐标、y向坐标与z向坐标,下标i表示肋编号,下标j表示所在肋上点编号,下标k表示相邻肋间点在对应肋上点之间的编号,肋编号i的取值范围为从1到N-1,肋上点编号j的取值范围为从2到m,m表示天线肋的分段数,相邻肋间点在对应肋上点之间的编号k的取值范围为从1到j-1,xi,j、yi,j、zi,j分别表示第i根肋上第j个肋上点的x向坐标、y向坐标与z向坐标,xi+1,j、yi+1,j、zi+1,j分别表示与第i根肋相邻的第i+1根肋上第j个肋上点的x向坐标、y向坐标与z向坐标;
5b)根据伞状天线圆形口径的闭合特性,按照下式计算第N根肋与第1根肋构成的肋间点的坐标:
Figure FDA0002228122590000061
Figure FDA0002228122590000062
Figure FDA0002228122590000063
其中,xN,j,k、yN,j,k、zN,j,k分别为第N根肋与第1根肋构成的肋间点的x向坐标、y向坐标与z向坐标,下标N表示第N根肋编号,下标j表示第N根肋上点编号,下标k表示第N根肋与第1根肋构成的肋间点在对应肋上点之间的编号,肋上点编号j的取值范围为从2到m,m表示天线肋的分段数,第N根肋与第1根肋构成的肋间点在对应肋上点之间的编号k的取值范围为从1到j-1,xN,j、yN,j、zN,j分别表示第N根肋上第j个肋上点的x向坐标、y向坐标与z向坐标,x1,j、y1,j、z1,j分别表示第1根肋上第j个肋上点的x向坐标、y向坐标与z向坐标。
CN201710388179.6A 2017-05-27 2017-05-27 基于面片积分公式的伞状天线结构优化设计方法 Active CN107240780B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710388179.6A CN107240780B (zh) 2017-05-27 2017-05-27 基于面片积分公式的伞状天线结构优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710388179.6A CN107240780B (zh) 2017-05-27 2017-05-27 基于面片积分公式的伞状天线结构优化设计方法

Publications (2)

Publication Number Publication Date
CN107240780A CN107240780A (zh) 2017-10-10
CN107240780B true CN107240780B (zh) 2020-04-17

Family

ID=59986017

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710388179.6A Active CN107240780B (zh) 2017-05-27 2017-05-27 基于面片积分公式的伞状天线结构优化设计方法

Country Status (1)

Country Link
CN (1) CN107240780B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107729705B (zh) * 2017-11-29 2020-10-16 中国电子科技集团公司第五十四研究所 一种面天线单块面板精度的测算方法
CN108110432B (zh) * 2017-12-18 2020-03-10 中国电子科技集团公司第五十四研究所 一种反射面天线面板分块方法及反射面天线
CN112685943B (zh) * 2021-01-05 2023-06-02 南阳理工学院 一种伞状柔性张拉膜聚光器聚光性能预测方法和系统
CN113764899B (zh) * 2021-08-04 2022-11-18 同济大学 一种肋网式可展开天线的网面安装方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6560570B1 (en) * 1999-11-22 2003-05-06 Sandia Corporation Method and apparatus for connecting finite element meshes and performing simulations therewith
CN104484537A (zh) * 2014-12-30 2015-04-01 西安电子科技大学 一种金属桁架式天线罩肋截面尺寸的机电集成优化方法
CN106338656A (zh) * 2016-08-23 2017-01-18 哈尔滨哈玻拓普复合材料有限公司 大型地面天线罩空间肋阵的散射云图获取方法
CN106354906A (zh) * 2016-08-19 2017-01-25 西安电子科技大学 一种基于三角形面片误差的网状可展开天线找形方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6560570B1 (en) * 1999-11-22 2003-05-06 Sandia Corporation Method and apparatus for connecting finite element meshes and performing simulations therewith
CN104484537A (zh) * 2014-12-30 2015-04-01 西安电子科技大学 一种金属桁架式天线罩肋截面尺寸的机电集成优化方法
CN106354906A (zh) * 2016-08-19 2017-01-25 西安电子科技大学 一种基于三角形面片误差的网状可展开天线找形方法
CN106338656A (zh) * 2016-08-23 2017-01-18 哈尔滨哈玻拓普复合材料有限公司 大型地面天线罩空间肋阵的散射云图获取方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Analysis of Arbitrarily Shaped Dielectric Radomes Using Adaptive Integral Method Based on Volume Integral Equation;Jing-Li Guo;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20060731;第54卷(第7期);全文 *
星载网状可展开天线结构分析与优化设计;丁波;《西安电子科技大学2016年硕士论文库》;20151201;全文 *

Also Published As

Publication number Publication date
CN107240780A (zh) 2017-10-10

Similar Documents

Publication Publication Date Title
CN107240780B (zh) 基于面片积分公式的伞状天线结构优化设计方法
CN102968532B (zh) 65m口径大型反射面天线结构机电集成设计方法
CN106446469B (zh) 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法
Tian et al. Data-driven modelling and optimization of stiffeners on undevelopable curved surfaces
CN107016161B (zh) 基于Zernike多项式和三角函数的赋形反射面天线型面描述方法
CN112364467B (zh) 一种放松反射面天线远场分析电磁网格尺寸的方法
CN106991210B (zh) 一种基于机电耦合模型的赋形反射面天线电性能预测方法
CN107357962A (zh) 一种基于自适应代理模型的天线罩肋截面尺寸优化方法
CN107391880A (zh) 一种精确高效的天线罩电性能分析方法
CN102073754A (zh) 基于误差因素的反射面天线机电综合分析方法
CN107169174B (zh) 最优焦距下伞状天线机电集成优化设计方法
CN107590309B (zh) 基于近似计算公式的网状天线电性能时域特性分析方法
You et al. A Quick Calculation Method for Radiation Pattern of Submillimeter Telescope with Deformation and Displacement
CN108399274B (zh) 一种共形承载天线的机电综合分析方法
CN116467998A (zh) 一种基于等效电路的极化转换超表面快速计算方法及系统
Garza et al. High-order Chebyshev-based Nyström methods for electromagnetics
CN114169201B (zh) 面向电性能的反射面天线结构加权优化方法
CN106934097B (zh) 面向电性能的空间网状天线关键动力学模态选取方法
CN109408957B (zh) 一种考虑确定性误差的伞状天线随机误差分析方法
CN113468759A (zh) 一种电磁问题确定方法、系统及存储介质
Lei Radiation pattern analysis of reflector antennas using CAD model-based physical optics method
CN113722902A (zh) 基于神经网络的赋形反射面天线最佳吻合参数估计方法
CN109408956B (zh) 基于单元节点位移的反射面天线平均功率方向图建模方法
Chen et al. Topology optimization for the layout design of radar absorbing coatings in cavities
Zhao et al. Robust design optimization of benchmark aerodynamic case based on polynomial chaos expansion

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