环境风载荷下大口径反射面天线波束指向分析
1.
2.
3.
4.
5.
6.
Analysis of beam pointing of large aperture reflector antenna under ambient wind load
1.
2.
3.
4.
5.
6.
通讯作者:
收稿日期: 2023-10-20 修回日期: 2023-11-21
基金资助: |
|
Received: 2023-10-20 Revised: 2023-11-21
作者简介 About authors
栾天(1999—),男,黑龙江大庆人,博士生,从事大口径天线结构变形控制研究,E-mail:
关键词:
Keywords:
本文引用格式
栾天, 薛松, 连培园, 余佳恒, 杜雨轩, 王猛, 赵武林, 卢波, 许谦, 王从思.
LUAN Tian, XUE Song, LIAN Peiyuan, YU Jiaheng, DU Yuxuan, WANG Meng, ZHAO Wulin, LU Bo, XU Qian, WANG Congsi.
现阶段,国内外学者主要采用2种方法来消除风载荷对大口径反射面天线变形的影响。第1种方法是安装天线罩。例如:麻省理工学院的Olmi等[4]为Haystack天文台的37 m口径的射电望远镜安装了天线罩,以减小风载荷对天线变形的影响。第2种方法是采用伺服控制系统来调整天线的指向,从而减小风载荷的影响。Kim等[5]将天线视为理想刚体,通过仿真分析验证了模糊PID(proportional-integral-differential,比例-积分-微分)控制相较于传统控制方法的优越性。Gawronski[6]针对美国国家航空航天局的深空网络天线控制系统,验证了PI(proportional-integral,比例-积分)控制算法的可靠性,但存在对高频段控制效果不佳的问题;此外,还证明了LQG(linear quadratic Gaussian,线性二次高斯)控制和
然而,上述针对天线伺服控制的研究通常将天线的反射面假设为刚性,即不考虑反射面的柔性变形,并将天线视为整体进行控制。但对于大口径反射面天线,其反射面表面的结构变形对波束指向的影响不能忽略,因此须对反射面的变形机理进行研究。现阶段,天线反射面风载特性的研究方法主要包括理论分析、数值模拟、风洞试验以及现场实测等。因受制于成本、普适性和实验难度等,数值模拟是目前最常用的技术手段。数值模拟是指利用CFD(computational fluid dynamics,计算流体动力学)方法,通过计算机对物体的风载特性进行仿真分析,其优势是计算成本低、不受尺寸限制,且可以较方便地获得不同工况下的风载特性[9]。Young等[10]针对风载荷对双子座望远镜的球形保护罩的影响进行了数值模拟,结果表明,不同风场环境下球形保护罩的风载特性不同。Du等[11-12]对相控阵雷达的风载特性进行了数值模拟。王春圆、刘岩等[13-14]对大型射电望远镜的风载特性进行了数值模拟,并通过与风洞试验结果进行对比来验证数值模拟的准确性。Ladd等[15]对GMT(giant Magellan telescope,巨型麦哲伦望远镜)进行了CFD数值模拟,所得结果可为望远镜的结构设计、选址提供参考。综上,国内外学者虽已利用数值模拟及风洞试验等方法对大口径反射面天线的风载特性进行了分析,但对其面形精度和波束指向特性等的分析较少。
为此,笔者以奇台射电望远镜(Qitai radio telescope, QTT)的110 m大口径反射面天线为研究对象,利用CFD数值模拟方法来分析和计算风扰情况下天线反射面表面的风压系数和变形情况,并结合天线结构的变形情况对其波束指向特性的变化规律进行研究,旨在为大口径反射面天线的抗风结构设计与系统控制研究提供理论指导。
1 大口径反射面天线风场的数值模拟
图1
根据图1所示流程,首先采用Space Claim软件对QTT的110 m大口径反射面天线进行风场建模,然后采用Fluent软件对天线反射面表面的风压分布情况进行数值模拟,最后使用Tecplot软件对天线反射面表面的风压系数进行可视化后处理。
1.1 天线计算域模型的构建
由于受到网格生成方法的限制,在天线桁架杆件处无法有效划分出合适的网格,且反射面的背架结构对反射面风压分布的影响较小,故在划分网格时,忽略背架结构,仅建立天线反射面的抛物面模型,以分析其风压分布情况。计算域是指流体流经的区域,对于外部流动,计算域的边界应远离天线模型,以使湍流运动充分发展。若计算域尺寸过大,则会导致网格数量增加,从而造成计算量增大,这会大大降低计算效率,因此须对计算域模型的尺寸进行限制。
反射面天线计算域的尺寸参数包括高度H、上游尺寸
根据阻塞率的计算公式和上述尺寸参数,可得:
式中:γ为阻塞率。
由
图2
图2
大口径反射面天线计算域模型
Fig.2
Computational domain model of large aperture reflector antenna
1.2 天线反射面表面风压系数的计算
在求解前,应对大口径反射面天线计算域模型设置合适的边界条件:入口设为风的速度入口,假设出口处所有物理变量的法向梯度均为零;计算域的上表面及2个侧面为滑动壁,计算域底面和天线反射面为非滑动壁。另外,在计算前还须对整个计算域进行网格划分。鉴于非结构化网格具有较好的自适应性,能够解决复杂流场的网格划分问题,本文采用非结构化网格单元对所构建的计算域模型进行网格划分。
RANS(Reynolds-averaged Navier-Stokes,雷诺平均湍流模型)通过引入合理的湍流模型来对雷诺应力项进行求解,该模型在求平均风压等方面具有明显优势[17]。鉴于
式中:uz 为任一高度处的平均风速;Ti 为任意点处的温度;Cμ 为常数,常取0.09;l为湍流积分尺度。
利用Fluent软件对大口径反射面天线的网格模型进行求解,得到其反射面表面的风压分布情况。通常情况下,采用风压系数来描述风载荷对结构的力作用。结构表面风压系数
式中:p为作用在结构表面的风压,pref为参考点处的风压,
天线反射面表面所受的实际压力为正、背面的压力差,则风压系数
式中:
根据
表1 不同俯仰角下天线反射面表面的风压系数
Table 1
俯仰角 | 反射面正面 | 反射面背面 |
---|---|---|
5° | ||
15° | ||
30° | ||
45° | ||
60° | ||
75° | ||
90° |
由表1所示的天线反射面正、背面的风压系数分布情况可以看出,不同工况下最大变形均出现在反射面边缘,且多为迎风口边缘。这是因为来风在天线反射面的边缘被分开,迎风处反射面正面的风向上产生正压,背面的风向下产生涡流,形成负压,正、背面的作用效果叠加,产生较大风压。此外,在某些俯仰角工况下,风会在天线反射面的背面产生绕流,进而导致反射面边缘处的风压较大。
图3
图3
天线反射面面板表面风压系数随俯仰角的变化趋势
Fig.3
Variation trend of wind pressure coefficient on antenna reflector panel surface with pitch angles
2 大口径反射面天线结构变形及面形精度分析
2.1 天线结构变形的有限元分析
对于所研究的110 m大口径反射面天线,本文只分析其反射面的流固耦合变形,即在有限元模型中将风载荷以风压的形式直接作用于反射面上,从而分析反射面的结构变形。为方便计算,分析时对天线的有限元模型进行简化,只建立反射面及背架结构的模型,将约束点设置在背架结构上,并在反射面上施加压力载荷。
鉴于直接利用天线反射面上网格节点处风压系数计算风载荷的过程非常复杂,为确定用于后文结构变形分析的风载荷值,拟对反射面进行区域划分。反射面表面的风压系数采用所有网格节点处的平均风压系数来近似表示,随后基于该平均风压系数来计算反射面所受的压力,并将压力平均分布至反射面面板的各个顶点上。图4所示为天线反射面表面的平均风压系数的近似计算示意。
图4
图4
天线反射面表面平均风压系数的近似计算
Fig.4
Approximate calculation of average wind pressure coefficient for antenna reflector surface
根据平均风压系数,可以计算出每一块反射面面板上的均布载荷,其计算式如下:
式中:F为反射面面板所受的力,A为单一反射面面板的面积,
在天线有限元模型中,将每一块反射面面板所受的载荷平均分布到4个顶点处,即可将流固耦合问题转换为静力学问题,通过静力学有限元分析可快速计算得到天线反射面的结构变形。
图5
图5
不同俯仰角下天线反射面的变形情况
Fig.5
Deformation of antenna reflector under different pitch angles
由图5可以看出,6 m/s风速下天线反射面变形的变化趋势与风压系数的变化趋势相对应;当俯仰角较小时,天线下端迎风口区域的变形量最大。
图6
图6
不同俯仰角下天线反射面的变形极值
Fig.6
Extreme value of deformation of antenna reflector under different pitch angles
2.2 天线面形精度分析
通常采用结构变形量的均方根误差(root mean square error, RMSE)作为衡量天线反射面变形程度的指标。对不同风速、俯仰角工况下天线反射面变形量的RMSE进行计算,结果如图7所示。
图7
图7
不同风速、俯仰角下天线反射面变形量的RMSE
Fig.7
RMSE of antenna reflector deformation under different wind speeds and pitch angles
由图7可以看到,随着风速的提高,天线反射面变形量的RMSE逐渐增大;在天线俯仰角为60º的工况下,反射面的变形较严重,应作为控制调节的关注重点。
3 大口径反射面天线波束指向分析
结构变形会影响天线的电性能,从而导致增益损失、指向偏差等问题。从天线的方向图中可以看出天线在空间各个方向上接收电磁波的能力。通过绘制不同工况下大口径反射面天线的方向图,即可得到其波束指向的变化情况。
其中:
式中: E 为远区辐射电场,
根据2.1节,可以得到天线反射面上各点的径向位移。由于CFD数值模拟耗时较长,为了快速计算得到不同俯仰角、风速工况下天线电性能的变化,可先提前计算并存储反射面表面的风压系数,随后即可快速计算得到相应风速下反射面的结构变形量,最后根据变形量得到天线的方向图,以分析其波束指向特性。大口径反射面天线波束指向的分析流程如图8所示。
图8
图8
大口径反射面天线波束指向分析流程
Fig.8
Analysis process for beam pointing of large aperture reflector antennas
图9
图10
图11
图11
6 m/s风速下天线第一副瓣电平的变化情况
Fig.11
Variation of the first sublobe level of antenna at 6 m/s wind speed
4 结 论
本文采用CFD数值模拟方法分析了大口径反射面天线表面风压系数的分布情况,从而得到了其风载特性。随后,利用风压系数计算了天线各反射面面板所受的风载荷,并对天线进行有限元仿真分析,得到了天线的结构变形规律。最后,通过天线结构变形量—电磁场的耦合计算,得到了其反射面波束指向特性的变化情况。具体结论和建议如下:
1)天线反射面迎风口处所受的风压较大,且此处的变形较大,在后续的柔性控制中应重点关注。
2)某些工况下天线的结构变形和波束指向偏差较大,故不同工况的控制应具有针对性。
3)由于天线的结构复杂,本文仅考虑了反射面的变形情况,后续研究应综合考虑其他结构的变形。
通过分析大口径反射面天线的风载特性可以快速计算得到其反射面的结构变形以及电性能的变化情况,这可为天线的风扰控制提供参考。
参考文献
大型射电望远镜天线发展动态及机电耦合应用
[J].DOI:10.1360/sspma-2022-0121 [本文引用: 1]
Development of large-aperture radio telescopes and applications of coupled structural-electromagnetic theory
[J].DOI:10.1360/sspma-2022-0121 [本文引用: 1]
大型射电望远镜天线主动面补偿研究进展
[J].DOI:10.1360/sspma2017-00011 [本文引用: 1]
Development of active surface technology of large radio telescope antennas
[J].DOI:10.1360/sspma2017-00011 [本文引用: 1]
风载荷在雷达天线结构设计中的考虑
[J].
Consideration of wind load in radar antenna structure design
[J].
A comparison of radome- and astrodome-enclosed large radio telescopes at millimeter wavelengths: the large millimeter telescope
[J].
Simplified fuzzy-PID controller of data link antenna system for moving vehicles
[C]//
Antenna control systems: from PI to H∞
[J].
Intelligent PD-type fuzzy controller design for mobile satellite antenna tracking system with parameter variations effect
[C]//
A review on control methodologies of disturbance rejections in optical telescope
[J].
AIJ guide for numerical prediction of wind loads on buildings
[J].
Numerical simulations of airflow in telescope enclosures
[J].
Computation of fluctuating wind pressure and wind loads on phased-array antennas
[J].
雷达天线风载特性的数值计算方法及应用研究
[D].
Research and application on numerical method for characteristic computation of wind loads on radar antennas
[D].
巨型射电望远镜风荷载特性的数值模拟研究
[D].
Numerical simulation study on characteristics of wind loads of huge radio telescope
[D].
大型射电望远镜结构风荷载特性研究
[J].DOI:10.3969/j.issn.1007-2276.2015.01.026 [本文引用: 1]
Wind load characteristics of large radio telescope
[J].DOI:10.3969/j.issn.1007-2276.2015.01.026 [本文引用: 1]
Computational fluid dynamics modeling and analysis for the Giant Magellan Telescope (GMT)
[C]//SPIE Astronomical Telescopes + Instrumentation,
CFD模拟技术在大型复杂结构工程中的应用
[D].
Application of CFD simulation technology in large and complex structural engineering
[D].
新疆奇台110米射电望远镜
[J].DOI:10.1360/sspma2014-00039 [本文引用: 1]
Xinjiang Qitai 110 m radio telescope
[J].DOI:10.1360/sspma2014-00039 [本文引用: 1]
A comparison between GO/aperture-field and physical-optics methods of offset reflectors
[J].
大型反射面天线抗风扰控制补偿关键技术研究
[D].
Research on key technologies of control compensation for large antenna under wind disturbance
[D].
A pattern approximation method for distorted reflector antennas using piecewise linear fitting of the exponential error term
[J].
/
〈 |
|
〉 |
