面向指向精度的奇台射电望远镜轮轨误差逆向设计方法

李娜 ,  段宝岩 ,  李向阳 ,  郑彬 ,  吴江

工程(英文) ›› 2022, Vol. 13 ›› Issue (6) : 209 -216.

PDF (1882KB)
工程(英文) ›› 2022, Vol. 13 ›› Issue (6) : 209 -216. DOI: 10.1016/j.eng.2021.11.015

面向指向精度的奇台射电望远镜轮轨误差逆向设计方法

作者信息 +

A Reverse-Design Strategy for the Track Error of the Qi Tai Telescope Based on Pointing Accuracy

Author information +
文章历史 +
PDF (1926K)

摘要

奇台射电望远镜(QTT)的口径为110 m,计划成为世界上最大的全可动望远镜。理想情况下,望远镜的重复指向精度误差应小于2.5 角秒(arcsec);因此,望远镜的结构必须满足超高精度的要求。为此,本文提出了一种轮轨表面的逆向设计方法,用以降低望远镜整体设计与制造的难度。首先,本文使用偏度系数和峰度系数方法验证了轮轨误差测试数据的分布特征。根据这些分布特征,采用了双尺度模型仿真分析了方位轮轨误差。在本文的实验中,长周期、小幅度的误差以大尺度为特征,采用三角函数描述,短周期、大幅度的误差以小尺度为特征,利用分形函数来模拟。本文首先基于双尺度模型,推导了误差对指向精度的影响机理;其次,推导了轮轨误差的均方根(RMS)与望远镜指向精度误差均方根之间的关系;最后,从容许的指向精度误差来推导出轮轨误差容许的RMS值。为了验证所提新方法的有效性,选择了两种典型的射电天文望远镜[绿岸望远镜(GBT)和大型毫米波望远镜(LMT)]作为实验实例。通过对比发现,望远镜指向精度的理论计算值与实测值基本一致,最大误差小于10%。

Abstract

The Qi Tai Telescope (QTT), which has a 110 m aperture, is planned to be the largest scale steerable telescope in the world. Ideally, the telescope's repeated pointing accuracy error should be less than 2.5 arc seconds (arcsec); thus, the telescope structure must satisfy ultra-high precision requirements. In this pursuit, the present research envisages a reverse-design method for the track surface to reduce the difficulty of the telescope's design and manufacture. First, the distribution characteristics of the test data for the track error were verified using the skewness coefficient and kurtosis coefficient methods. According to the distribution characteristics, the azimuth track error was simulated by a two-scale model. The error of the long period and short amplitude was characterized as large-scale and described by a trigonometric function, while the short period and high amplitude error was characterized as small-scale and simulated by a fractal function. Based on the two-scale model, effect of the error on the pointing accuracy was deduced. Subsequently, the relationship between the root mean square (RMS) of the track error and the RMS of the pointing accuracy error of the telescope was deduced. Finally, the allowable RMS value of the track error was derived from the allowable pointing accuracy errors. To validate the effectiveness of the new design method, two typical radio telescopes (the Green Bank Telescope (GBT) and the Large Millimeter Telescope (LMT)) were selected as experimental examples. Through comparison, the theoretical calculated values of the pointing accuracy of the telescope were consistent with the measured values, with a maximum error of less than 10%.Graphical abstractQi Tai Telescope (QTT) with a 110 m caliber, it will become the largest scale steerable telescope in the world. Its repeated pointing accuracy error should be less than 2.5 arc seconds and the telescope structure must satisfy an ultra-high precision requirements. To reduce the difficulty of telescope's design and manufacture, a reverse design method of the track surface was proposed. Firstly, the distribution characteristics of the test data of track error were verified using the skewness coefficient and kurtosis coefficient methods. According to the distribution characteristics, the azimuth track error was simulated by a two scales model. The error of long-period-short-amplitude was defined as the large scale, which was described by the trigonometric function. The short-period-high-amplitude error was recognized as the small scale, and simulated by the fractal function. Based on the two scales model, effect of it on the pointing accuracy was deduced. And then the relationship between the RMS (Root Mean Square) of track error and the RMS value of the telescope's pointing accuracy errors was deduced. Finally, the allowable RMS value of track error was derived from the allowable pointing accuracy errors. To validate the effectiveness of the new design method, two typical radio telescopes (the Green Bank Telescope and the Large Millimeter Telescope) were selected as experimental examples. Through comparison, the theoretical calculated values of the telescope's pointing accuracy are in good agreement with the measured values, with a maximum error less than 10%.

关键词

逆向设计 / 跟踪误差 / 大型射电望远镜 / 双尺度模型 / 指向精度

Key words

Reverse design / Track error / Large radio telescope / Two-scale model / Pointing accuracy

引用本文

引用格式 ▾
李娜,段宝岩,李向阳,郑彬,吴江. 面向指向精度的奇台射电望远镜轮轨误差逆向设计方法[J]. 工程(英文), 2022, 13(6): 209-216 DOI:10.1016/j.eng.2021.11.015

登录浏览全文

4963

注册一个新账户 忘记密码

1、 引言

口径为500 m的球面射电望远镜(FAST)于2016年在中国贵州省建成,是世界上最大、最强的单碟射电望远镜。口径110 m全可动的奇台射电望远镜(QTT)现已开始筹建,建设单位为中国科学院新疆天文台。建成后,它将成为世界上最大口径的全可动射电天文望远镜[1‒2]。虽然FAST在3 GHz以下的工作频段比QTT要灵敏得多,但它的天空覆盖范围有限;相比之下,QTT将在150 MHz至115 GHz的工作范围内运行,并具有全天空覆盖的功能。

QTT的天空覆盖范围和频率覆盖范围将使其成为一种用途非常广泛的仪器,可以观测许多分子线发射以及银河系中心,并大大提高脉冲星引力波测量的准确性。QTT所需的指向精度为2.5角秒(arcsec),这对于质量为6000 t、高度为100 m、反射表面积为26个篮球场面积的望远镜结构来说是极难实现的[3‒4]。以位于墨西哥普埃布拉省Sierra Negra山顶的大型毫米波望远镜(LMT)的指向误差为例[5],在误差补偿之前,LMT的服务环境、机械校准和伺服控制引起的指向误差分别为10.76角秒、5角秒和7.74角秒[6‒7]。但是,传统的设计思想对多个误差源一视同仁;因此,望远镜的最终指向精度非常有限[8]。为了满足QTT设计中极高的指向精度要求,测量中应识别关键误差源;然后,对关键误差源进行优先级控制。

望远镜的结构误差是影响其指向精度的关键因素。误差源主要涉及基座、反射面或轮轨的变形[9]。焊接轮轨是通过焊接制成的,单轨之间采用相对光滑的轮轨接头,可以使整体轮轨经久耐用。以高难度焊接工艺为代价的焊接轮轨适用于高精度和高负载的望远镜[10]。然而,轮轨的加工和焊接过程引入了随机误差,主要包括微观表面误差和宏观轮轨不平度[8,11‒12]。这两个尺度上的误差会导致望远镜的方位轴和俯仰角轴出现误差,从而影响指向精度[13]。

在早期轮轨误差对望远镜影响的研究中,研究人员通过测量仪器获得轮轨误差数据。通过几何转换公式,将不平度直接转化为望远镜方位角和俯仰角的误差[14]。此外,Pisanu等[15]考虑了轮轨误差和温度漂移引起的方位框架变形对指向精度的影响。Kong等[16]研究了轮轨误差测试数据并进行了实验以确定指向精度。然后根据测试数据分析了轮轨误差与指向精度的关系。此外,研究人员还研究了非线性轮轨误差对望远镜指向性能的影响[17]。先前的研究表明,轮轨误差是望远镜设计、建造和运行阶段必须不断考虑的关键误差源。轮轨的结构误差具有明显的多尺度性,但在误差建模中却没有考虑到这一点。因此,指向误差与轮轨误差之间的准确关系尚未建立。这一不足是本文提出的逆向设计策略的基础。

具有有限元建模(FEM)或时域有限差分(FDTD)的高级迭代计算通常用于表面性能预测。但是,这种传统的设计过程受到人为引导的固有误差的影响。本研究基于双尺度误差模型及该误差对望远镜指向的影响机理,提出一种轮轨表面的逆向设计方法,可有效搜索和优化表面结构的核心参数以达到用户要求的表面特性。该模型是一种有效的大型轮轨表面逆向设计技术,有助于减少传统曲面设计中的计算量和人力资源,避免复杂的专业迭代仿真和参数化扫描研究,可为缺乏机械知识的用户服务。

2、 轮轨误差对望远镜指向精度的影响

2.1 轮轨误差测试数据来源与特征分析和验证

根据某些大型射电天文望远镜[如LMT、Robert C.Byrd Green Bank Telescope (GBT)]的大量实验数据分析结果,可以初步假设这些大型望远镜的轮轨误差数据服从高斯分布[18]。在实际工程应用中,这一结论还应通过使用偏度系数法和峰度系数法等测试方法进行验证[19‒20]。这两种方法的详细描述如下:

Sk=i=1nxi-x¯3/ns3(1)

     Uu=i=1n(xi-x¯)4/ns4-3(2)

式中,Sk是偏度系数;Uu是峰度系数;xi是轮轨误差测试数据,i是第i个轮轨误差测试数据;n为测试点数;x¯=1ni=1nxis=1ni=1nxi-x¯2分别是测试数据的平均值和标准差。

如果不均匀测试数据的偏度和峰度都接近于零,则可以确认数据集整体服从高斯分布。对于一组测试数据样本,如果计算的偏度为正,则与数据集遵循标准高斯分布的情况相比,峰度呈现出向较小值侧的偏差。相反,如果计算的偏度为负,则与数据集遵循标准高斯分布的情况相比,峰度呈现出向较大值侧的偏差。如果峰度计算为正,则该分布相对标准高斯分布更锐利。但是,如果计算峰度为负,则该分布相对标准高斯分布更平坦。如果测试数据服从高斯分布,那么偏度和峰度也相应地服从高斯分布,数学期望值为零。偏度和峰度的均方根(RMS)值可以表示如下:

σs=6 (n-2)(n+1) (n+3)            (3)

σu=24n (n-2) (n-3)(n+1)2(n+3) (n+5)(4)

式中,σsσu分别表示偏度和峰度方法的RMS值。

假设轮轨误差的测试数据服从高斯分布,则置信度系数α可以设置为0.05(α=0.05)。当偏度系数Sk > 1.96 σs和峰度系数Uu > 1.96 σu时,假设无效,测试数据集不服从高斯分布。否则,假设成立;也就是说,测试数据集服从高斯分布。

2.2 轨道误差的双尺度模型

轮轨的结构误差主要包含两个方面:单轨加工产生的随机误差和轮轨使用过程中重力引起的变形。这两种类型的误差具有不同的来源和分布特征[21‒22]。加工过程的误差是随机的,呈现高频率和小振幅特征。结果表明,它被确定为小尺度误差,并通过周期函数计算。另一方面,自重变形是系统性的,表现为低频大振幅;因此,它被描述为大尺度误差。基于两种尺度分布的特征,大尺度轮轨误差采用傅里叶级数[23]仿真,小尺度轮轨误差作为拟合残差,采用分形函数[24]建模。最后,建立了轮轨不平度的双尺度误差模型,如式(5)所示。

F(x)=f1(x)+i=1gf2i(Ai,Di,Li,x0i,y0i)=a0+n=1qancos(nω0x)+bnsin(nω0x)+i=1gf2i(Ai,Di,Li,x0i,y0i)

n1i=lg(1/Li) /lg(γ)Mi=lg(Nγn1i/2) /lg(γ)(5)

式中,f1(x)是傅里叶级数函数;f2i(Ai,Di,Li,x0i,y0i)是第i段的Weierstrass-Mandelbrot (W-M)分形函数;a0是常数系数;a1,…,a,b1,…,b是不同的幅度分量;m是傅里叶级数的展开阶;ω0f1(x)的角频率;Ai是W-M分形函数f2(x)的幅度系数;Lif2(x)的采样长度;N为采样点数;Dif2(x)的分形维数;γ是轮轨表面谐波的频率;y0ix0i分别表示f2(x)的纵向和横向位移。

如式(5)所示,轮轨误差的拟合精度由f1(x)f2(x)共同决定。更高的展开级数(m)对应更短的采样长度(L)和更高的函数拟合精度。因此,确定mL的取值是建立双尺度误差模型的关键。优化模型可以表示如下,其中T是转置符号。

F(y)=j=1N(f2i(xj,y)-y(xj))2 /N(6)

f2(x)=A(D-1)n=n1Q1γ(2-D)ncos2πγn(x+Δx)+Δy(7)

findy=(Ai,Di,Li,x0i,y0i)TminF(y)=j=1Nf2i(xj,y)-y(xj)2 /Ns.t.0<x0i<xmaxLmin<Li<Lmaxymin<y0i<ymaxAmin<Ai<AmaxDmin<Di<Dmax(8)

3、 基于指向误差的轮轨误差允许值

3.1 轮轨误差对指向精度误差的影响

为研究轮轨误差对望远镜整体结构和指向精度的影响,本实验建立了4个坐标系,如图1所示。

图1 望远镜坐标系示意图。B:反射体的方位角。

(1)大地坐标系(OXYZ)。该系统的原点位于望远镜方位角轮轨的正中心。Z轴垂直于地球,Y轴正方向指向南方。

(2)方位轴坐标系(OaXaYaZa)。它的原点位于轮轨的中心,而它的Z轴与方位轴重合,并随着方位轴的旋转而偏转。当望远镜沿方位轴无误差且方位角为零时,OaXaYaZa等于OXYZ

(3)俯仰轴坐标系(OeXeYeZe)。它的原点位于方位轴的中心,而它的Z轴与俯仰角轴重合,并随着俯仰角轴的旋转而偏转。当望远镜无轴误差时,俯仰角等于90°,方位角等于0°时,坐标与大地坐标在Z方向只有h的高差。

(4)反射体坐标系(OrXrYrZr)。如果俯仰角等于90°,则它与Z方向俯仰角坐标系的高度差为h1

如图2所示,黑线为方位坐标系中方位框架的初始位置。点1、2、3和4代表4个滚轮,它们的高度值为Z(B+π4)Z(B-π4)Z(B-3π4)Z(B+3π4),其中,Z(x)是轮轨描述函数,B是反射体的方位角。红线代表变形的方位框架。Δx1Δy1分别是方位轴沿X轴和Y轴的旋转,这是由滚轮1和3之间的高度差引起的。类似地,Δx2Δy2是由滚轮2和4之间的高度差引起的旋转。

Δx1=Z(B+π4)-Z(B-3π4)2rcos(π/4)Δx2=Z(B+3π4)-Z(B-π4)2rcos(π/4)   (9)

Δy1=Z(B-3π4)-Z(B+π4)2rcos(π/4)Δy2=Z(B+3π4)-Z(B-π4)2rcos(π/4)   (10)

式中,r是轮轨的半径。

图2 轮轨误差引起的方位坐标系倾斜变形。

此外,轮轨误差导致方位架的变形;这会导致方位架本身的扭转,使俯仰角轴围绕Z轴旋转,如图3所示。由滚轮1和2之间的高度差引起的俯仰轴平移为δ12,由滚轮3和4引起的位移表示为δ34

δ12=he-Z(B+π4)+Z(B-π4)2rcos(π/4)

δ34=heZ(A-3π4)-Z(A+3π4)2rcos(π/4)(11)

式中,he是俯仰轴的长度。

图3 轮轨误差引起的方位架畸变。

轮轨不平度引起的方位坐标系摄动矩阵表示如下:

Rt=1ϕtz-ϕty-ϕtz1ϕtxϕty-ϕtx1

式中,ϕtxϕtyϕtz是坐标系中方位坐标系的误差,可以根据双尺度模型计算得到。

轮轨整体倾斜会导致方位轴倾斜,安装方位轴时初始偏移为零。考虑到轨道误差,坐标变换可表示为:

xr´yr´zr´=RerRaeRtRbaxyz=-ϕtzcosE-ϕtysinEϕtx1=001+-ϕtzcosE-ϕtysinEϕtx0(12)

式中,RerRaeRba分别是俯仰轴坐标系到反射体坐标系、方位轴坐标系到俯仰轴坐标系、大地坐标系到方位轴坐标系的坐标转换矩阵;E是反射器的俯仰角。

由轮轨误差引起的指向误差可以描述如下:

Δt=xr´yr´zr´-xryrzr=-ϕtzcosE-ϕtysinEϕtx0(13)

然后,可以将反射体坐标系中的指向误差转换为大地坐标系中的方位角和俯仰角误差。方位角误差为η=Δθr/cosE,俯仰角误差为β=Δϕrθrϕr是反射体坐标系中的方向角和俯仰角,ΔθrΔϕr是反射体坐标系中指向误差的方向和俯仰角的分量。

ηβ=-(ϕtz+ϕez+ϕaz)-(ϕty+ϕey)tanE-ϕaycosBtanE+   ϕaxsinBtanE-ϕoysecE-ϕtx-ϕaxcosB-ϕaysinB-ϕex-ϕox=(14)

式中,ϕaxϕayϕaz是方位轴坐标系的误差,由轮轨整体倾角和方位轴安装误差引起;ϕexϕeyϕez是俯仰轴坐标系的误差,主要是由方位架变形和俯仰轴安装误差引起;ϕoxϕoy是电轴的误差,主要是由副面、馈源和主面的位置偏差引起。

如果只考虑轮轨误差对指向精度的影响,则ϕax=ϕay=ϕaz=0ϕex=ϕey=0ϕez=ϕox=ϕoy=0。通过将这些参数的取值代入式(14),可以得到指向精度如下:

ηβ=0tanE-1-100122r-122r-122r122r-122r-122r122r122r-h2r2h2r2-h2r2h2r2Z(B+π4)Z(B-π4)Z(B-3π4)Z(B+3π4)(15)

3.2 轮轨误差容许值的逆向求导

目前,望远镜指向精度的误差分布是一个复杂的工程问题,在设计阶段需要综合考虑多个因素[25]。通过结合工程经验、有限元结构分析、控制模型仿真、实测数据、工程估算等多种技术手段,可以确定各种误差源引起的分布式指向误差[26]。由于每个误差源都会影响望远镜的指向精度,因此可以假设由这些误差源引起的指向误差记为Δ1,...,Δi,...,Δn。此外,由于分布指向误差是标量和RMS值,天线的整体指向误差ΔPE可以表示如下[27]:

ΔPE=i=1nΔi2(16)

根据第2节的验证结果,如果轮轨误差的测试数据服从高斯分布,则可以从概率统计的角度进一步推导出轮轨误差与指向精度的关系模型。相应地,在确定了轮轨误差(σ)的RMS后,就可以计算出由其引起的天线方位角和俯仰角的指向误差的RMS值。随后,可以基于获得的误差分布推导出允许的跟踪误差值。望远镜的简化如图4所示。

图4 望远镜的简化示意图。

如图5所示,用4个点来支撑望远镜的方位架和轮轨,分别记为点1、2、3、4。方位架的底部结构是由这4个点组成的正方形。如果只考虑轮轨误差,则指向误差可描述如下:

图5 望远镜方位架和轮轨支撑点的位置。

ηβ=k0tanE-1-100122r-122r-122r122r-122r-122r122r122r-h22r2h22r2-h22r2h22r2Z(B+π4)Z(B-π4)Z(B-3π4)Z(B+3π4)(17)

式中,k是转换系数(k=180π×36001000=206.2648);h2是俯仰轴和轮轨之间的距离;ZB+π4ZB-π4ZB-3π4ZB+3π4分别是4个支撑点的高度,单位为mm。

l1=(-122rtanE+h2r2),l2=(-122rtanE-h2r2),l3=(122rtanE+h2r2),l4=(122rtanE-h2r2)

h1=-122r,h2=122r,h3=122r,h4=-122r

式(17)可简化为以下表达式:

ηβ=k l1Z1+l2Z2+l3Z3+l4Z4h1Z1+h2Z2+h3Z3+h4Z4(18)

假设Z1Z2是两个独立的随机变量,Z1服从正态分布N0,a2Z2服从N0,b2W=k1Z1+k2Z2k1Z1服从正态分布N0,k12a2k2Z2服从N0,k22b2。令随机变量M=kZZ服从正态分布N0,σ2,则z=m/kzm=1/kM的概率密度为:

fM(m)=fz(m)z(m)´=1k12πσe-mk22σ2=12πk2σe-m22k2σ2,-<m<(19)

因此,M=kZ服从N0,k22b2。令X=k1Z1Y=k2Z2。则W=X+Y

 fX(x)=12πk1σe-x22k12σ2,-<x<fY(y)=12πk2σe-y22k22σ2,-<y<(20)

使用卷积定理,可以得到式(21):

fW(w)=-+fX(x)fY(w-x)dxfW(w)=-+fX(y)fY(w-y)dy(21)

式(21)具体表示为:

fW(w)=-+fX(x)fY(w-x)dx=-+12πk1ae-x22k12a212πk2be-(w-x)22k22b2dx=12πk1ak2b-+e-x22k12a2e-(w-x)22k22b2dx=12πk1ak2b-+e-x2(k22b2+k12a2)-2k12a2wx+w2k12a22k12a2k22b2dx(22)

g1=12πk1ak2bg2=k22b2+k12a22k12a2k22b2g3=w2k12a2k22b2+k12a2-k12a2wk22b2+k12a22

则式(22)可进一步表示为:

fWw=g1e-g3-+e-g2x-k12a2wk22b2+k12a22dxlett=x-k12a2w(k22b2+k12a2)=

   g1e-g3-+e-g2t2dt(-<t<+) =

g1e-g3πg2                                           (23)

我们将g1g2g3代入式(19)得到以下等式:

fW(w)=12πk1σk2σ e-w22(k22b2+k12a2) π2k12a2k22b2(k22b2+k12a2)=12π(k22b2+k12a2) e-w22(k22b2+k12a2)(24)

我们假设Z1Z2Z3Z4是相互依赖的,并且服从高斯分布N0,σ2。此外,概 率密度函数可以写成fZiZ=12πσe-z22σ2。则X=l1Z1+l2Z2+l3Z3+l4Z4Y=h1Z1+h2Z2+h3Z3+h4Z4的概率密度可写为:

fXx=12πl12+l22+l32+l42σe-x22l12+l22+l32+l42σ2    (25)

fYy=12πh12+h22+h32+h42σe-x22h12+h22+h32+h42σ2(26)

因此,XY也遵循方差分别为l12+l22+l32+l42σ2h12+h22+h32+h42σ2的高斯分布。根据这些推导,具有高斯分布特征的轮轨误差数据的RMS与望远镜指向误差的RMS值之间的关系可表示为:

σA=k12r2tanE2+1r4σ

σE=k2rσ

σPE=σA2+σE2(27)

式中,σAσEσPE分别是轮轨误差的RMS等于角秒σ时的方位角误差、俯仰角误差和天线指向角误差的RMS。

我们将望远镜总指向误差中由轮轨误差引起的部分表示为θ(RMS)。然后,轮轨处理和安装的不平度要求可以使用式(28)反向推导,这是轮轨误差与指向误差之间的概率统计模型。详细推导过程如下:

σPE=σA2+σE2<θ=kσtanE2+12r2+1r4<θ(28)

因此,

σ<θr2k2rtanE2+2r2+1(29)

4、 实验验证

4.1 GBT的实验验证

GBT是位于美国西弗吉尼亚州格林班克的国家射电天文台(NRAO)的射电望远镜,同时也是世界上几个最大的拥有单碟反射器的望远镜之一。目前,该望远镜由绿岸天文台运营。它高146 m,重7700 t。碟形天线尺寸为100 m×110 m;精度补偿后,最高指向精度为2角秒。GBT是典型的偏置射电望远镜,工作频率范围为0.1~116 GHz。GBT的轮轨由48条轮轨组成,直径为64 m,全长201 m。轮轨误差的RMS为0.0568 mm。

轮轨误差的测试方案描述如下[28‒29]。一块已知厚度和硬度的木板被放置在轮轨上以校准倾角仪。随后,旋转测距仪以测量轮轨。具体测量过程如下。首先将倾角仪安装在测距仪底部的4个滚轮装置上,并在滚轮下方放置一个钢盘。当轮子经过时,倾角仪会倾斜。如果方位架以恒定速度旋转,则通过放大或缩小钢盘厚度与倾角仪读数的比例来获得倾角仪读数与轨迹误差之间的关系。因此,该方法对轮轨的整体轮廓都进行了研究。采用这种基于数据分布特征的验证方案,本实验绘制了GBT轮轨误差测试数据的分布直方图,如图6所示。

图6 GBT轮轨误差测试数据分布直方图。

通过相关系数检验法分析,GBT无标度区间范围为[0.116,0.901],分形维数采用基于小波变换的分形维数法求解,GBT的分形维数为DGBT=1.602,轮轨不平度函数用于建模,ε1=2%,ε2=25%,ε3=20%,GBT天线轮轨不平度的局部描述如图7所示。

图7 GBT轮轨误差对比研究。

为了测量指向精度,本实验首先确定了一个无线电源列表;然后,用望远镜进行了交叉扫描观察。在观测过程中,本实验记录了望远镜在相关时间内的指向精度和辐射功率。对射电点源在俯仰角或方位方向进行扫描观测后,得到了功率在俯仰角或方位方向的变化曲线。指向精度误差的存在导致功率最大值和无线电点源的位置偏差。因此,望远镜在俯仰角方向的指向精度误差是通过在俯仰角方向扫描射电点源得到的,通过在方位方向上扫描得到该方向上的误差。

以GBT为例,本实验根据测试数据和式(17)构建了指向误差模型。然后,对应的轮轨误差允许值由指向误差的RMS值反向推导,根据式(27)给定轨迹误差的允许值,重新计算指向误差并与测量值进行比较(E=45°),具体如表1所示。

表1 计算的指向误差与GBT的测量值之间的比较。GBT的口径为64 m,轮轨的RMS为0.0568 m (°)

RMSPointing error
Calculated values (arcsec)Measured values (arcsec)Relative errors (%)
σA0.25910.23828.77
σE0.25890.24077.56
σPE0.36630.33868.18

4.2 LMT的实验验证

LMT的口径为50 m,副反射面直径为2.5 m。它位于墨西哥普埃布拉州,是世界上最大的毫米波射电望远镜[30‒31]。其观测波段为0.85~4 mm (75~350 GHz),表面精度为0.07 mm,指向精度为1角秒。轮轨直径为39.6 m,由20条轮轨组成。轮轨误差的RMS为0.1697 mm。用于LMT的测试方案与用于GBT的相同。本实验使用基于数据分布特征的验证方案,还绘制了LMT轮轨误差的测试数据分布直方图,如图8和图9所示。

图8 LMT轮轨误差测试数据分布直方图。

图9 LMT轮轨误差对比研究。

表2给出了LMT的计算指向误差与测量值的比较,最大误差为6.5%。

表2 LMT的计算指向误差与测量值的比较。LMT的口径为39.6 m,轮轨的RMS为0.1679 m (°)

RMSPointing error
Calculated values (arcsec)Measured values (arcsec)Relative errors (%)
σA1.23991.30224.78
σE1.23681.32286.50
σPE1.75131.85625.65

4.3 密云实验验证

北京密云50 m望远镜口径为32.5 m,由21条轨迹组成,全长102.1 m,RMS误差为0.3136 mm。根据相关系数检验方法,确定密云望远镜的无标度区间为[0.263,1.156],利用基于小波变换的分形维数法计算出50 m天线的分形维数为D = 1.1521,使用轨迹不平度函数进行建模,ε1=2%,ε2=25%,ε3=20%,密云天线轨迹不平度的局部描述如图10所示。

图10 密云望远镜轮轨误差对比研究。

5、 结论

传统的设计方案无法充分满足QTT等大比例望远镜极高的指向精度要求。为了降低QTT设计与制造的难度,本研究提出了一种轮轨表面的逆向设计方法。基于轮轨表面的误差特性,首先提出了双尺度模型来描述轮轨误差。然后,推导了轮轨误差对指向精度的影响机理。最后,使用这种策略,本研究从轮轨误差引起的指向误差中反向推导出轮轨误差的容许值。在GBT和LMT等望远镜工程上对所提的设计方法进行了实验验证。将所提出的方法与实测数据进行对比分析,结果表明本文所提出的设计方法是高效可行的。

参考文献

[1]

Kärcher HJ. Experience with the design and construction of huge telescope pedestals. In: Proceedings of Second Backaslog Workshop on Extremely Large Telescopes; 2003 Sep 9‒11; Backaskog, Sweden. SPIE; 2003. p. 233‒44.

[2]

Chu TS, Turrin R. Depolarization properties of offset reflector antennas. IEEE Trans Antennas Propaga 1973;21(3):339‒45.

[3]

Bolli P, Orlati A, Stringhetti L, Orfei A, Righini S, Ambrosini R, et al. Sardinia radio telescope: general description, technical commissioning and first light. J Astron Instru 2015;4(3‒4):1550008.

[4]

Antebi J, Kan FW. Precision continuous high-strength azimuth track for large telescopes. In: AngelJRP, GilmozziR, editors. Proceedings of Future Giant Telescopes; 2002 Aug 26‒28; Waikoloa, New Zealand. SPIE; 2002. p. 612‒23.

[5]

Schloerb FP. The large millimeter telescope. In: Proceedings of Ground-Based and Airborne Telescopes II; 2008 Jun 23‒28; Marseille, France. SPIE; 2008 p.626701.

[6]

Kaercher HJ, Baars JWM. Design of the large millimeter telescope/gran telescopio millimetrico (LMT/GTM). In: Butcher HR, editor. Proceedings of Astronomical Telescopes and Instrumentation; 2000 Mar 27‒Apr 1; Munich, Germany. SPIE; 2000. p. 155‒68.

[7]

Zhang XZ, Zhu XY, Kong DQ, Zheng L, Yao C, Zhang HB, et al. Measurements of electronic properties of the Miyun 50 m radio telescope. Res Astron Astrophys 2009;9(3):367‒76.

[8]

Guiar CN, Lansing FL, Riggs R. Antenna pointing systematic error model derivations. Telecommun Acquis Prog Rep 1986;42(88):36‒46.

[9]

Keitzer S, Kimbrell JE, Greenwald D. Deterministic errors in pointing and tracking system I: identification and correction of static errors. In: Proceedings of Acquisition, Tracking, and Pointing V; 1991 Aug 1; Orlando, FL, USA. SPIE; 1991. p. 406‒14.

[10]

Juneja G, Kan FW, Antebi J. Update on slip and wear in multi-layer azimuth track systems. In: Proceedings of Optomechanical Technologies for Astronomy; 2006 May 26‒31; Orlando, FL, USA. SPIE; 2006. p. 6273.

[11]

Deng D, Kiyoshima S. Numerical simulation of welding temperature field, residual stress and deformation induced by electro slag welding. Comp Mater Sci 2012;62(5):23‒34.

[12]

Zong Y, Hu N, Duan B, Yang G, Cao H, Xu W. Manufacturing error sensitivity analysis and optimal design method of cable-network antenna structures. Acta Astronaut 2016;120:182‒91.

[13]

Antebi J, Kan FW. Precision continuous high-strength azimuth track for large telescopes. In: Proceedings of Future Giant Telescopes; 2002 Aug 22‍‒‍28; Waikoloa, HI, USA. SPIE; 2002. p. 612‒23.

[14]

Gawronski W, Baher F, Quintero O. Azimuth-track level compensation to reduce blind-pointing errors of the Deep Space Network antennas. IEEE Antenn Propag Mag 2000;42(2):28‒38.

[15]

Pisanu T, Buffa F, Morsiani M, Pernechele C, Poppi S. Thermal behavior of the Medicina 32-meter radio telescope. In: Proceedings of Modern Technologies in Space- and Ground-based Telescopes and Instrumentation; 2010 Jun 27‒Jul 2; San Diego, CA, USA. SPIE; 2010. p. 35‒44.

[16]

Kong DQ, Wang SG, Wang JQ, Wang M, Zhang HB. A new calibration model for pointing a radio telescope that considers nonlinear errors in the azimuth axis. Res Astron Astrophys 2014;14(6):733‒40.

[17]

Belov Y, Naumov A, Chernikova S, Poperchenko B, Sapozhnikov B. The application of correlation analysis in the process of pointing calibration of 64-m Bear Lakes radio telescope. IEEE Antennas Propag Soc Int Symp 1997;1(1):568‒71.

[18]

Alvarez ML, Torres CT, Rios EH, Santos DC, Cuevas LC, Leon-Huerta A, et al. Metrology and surface adjustment of primary reflector panels on the LMT. In:Proceedings of Advances in Optical and Mechanical Technologies for Telescopes and Instrumentation; 2014 Jun 22‒27; Montréal, QC, Canada. SPIE; 2014. p. 3‒11.

[19]

Doane DP, Seward LE. Measuring skewness: a forgotten statistic? J Stat Educ 2011;19(2):1‒18.

[20]

Joanes DN, Gill CA. Comparing measures of sample skewness and kurtosis. J Roy Stat Soc 1998;47(1):183‒9.

[21]

Deng D, Murakawa H. Prediction of welding distortion and residual stress in a thin plate butt-welded joint. Comp Mater Sci 2008;43(2):353‒65.

[22]

Majumdar A, Bhushan B. Role of fractal geometry in roughness characterization and contact mechanics of surfaces. J Tribol 1990;112(2):205‒16.

[23]

Hu YZ, Tonder K. Simulation of 3-D random rough surface by 2-D digital filter and Fourier analysis. Int J Mach Tools Manufact 1992;32(1‒2):83‒90.

[24]

Qiao B, Liu S. Error assessment in modeling with fractal Brownian motions. Fract 2013;21(3‒4):1‒8.

[25]

Wielebinski R, Junkes N, Grahl BH. The effelsberg 100-m radio telescope: construction and forty years of radio astronomy. J Astron Hist Herit 2011;14(1):3‒21.

[26]

Gawronski W, Baher F, Quintero O. Azimuth-track level compensation to reduce blind-pointing errors of the deep space network antennas. IEEE Antennas Propag Mag 2000;42(2):28‒38.

[27]

Matsuzawa A, Saito M, Iguchi S, Nakanishi K, Saito H. Development of highaccuracy pointing verification for ALMA antenna. In: Proceedings of Groundbased and Airborne Telescopes V; 2014 Jun 22‍‒‍27; Montréal, QC, Canada. SPIE; 2014. p. 1‒7.

[28]

Prestage RM, Constantikes KT, Hunter TR, King LJ, Lacasse RJ, Lockman FJ, et al. The green bank telescope. Proc IEEE 2009;97(8):1382‒90.

[29]

Symmes A, Anderson R, Egan D. Improving the service life of the 100-meter green bank telescope azimuth track. In: Proceedings of Ground-Based and Airborne Telescope II; 2008 Jun 23‒28; Marseille, France. SPIE; 2008. p.1‒12.

[30]

Hughes DH, Correa JJ, Schloerb FP, Erickson N, Romero JG, Heyer M, et al. The large millimeter telescope. In: Proceedings of Ground-Based and Airborne Telescopes III; 2010 Jun 27‒Jul 2; San Diego, CA, USA. SPIE; 2010.p. 1‒13.

[31]

Okubo H, Kawai E, Kondo T, Osaki H, Nakajima J, Ishizuka K. Inclinometer measurement of Kashima 34 m antenna. Technol Dev Cent News 2001;11(19):13‒4.

基金资助

()

AI Summary AI Mindmap
PDF (1882KB)

4434

访问

0

被引

详细

导航
相关文章

AI思维导图

/