-
目前,美国联合世界多国计划建造30 m望远镜(TMT),它最终将成为地球上最先进和最强大的光学望远镜之一[1]。当TMT完成后,天文学家将能够研究太阳系和整个银河系及其邻近星系的恒星,并可观测宇宙的边缘,即时间开始的地方形成的星系。
TMT的光学系统由30 m口径的主镜(由492块口径为1.44 m,厚度为45 mm的六边形镜面拼接组成)、一个完全主动的副镜和一个铰接的第三镜组成。第三镜是一个口径为3.5 m×2.5 m的椭圆外形平面反射镜,其为目前世界上最大口径的平面反射镜。大口径反射镜的检测十分困难,而且它的外轮廓是非圆形的,更增加了其制造和测试的难度。
检测大口径平面镜的传统方法是利用三坐标测量仪(CMM)和瑞奇康芒法(Ritchey-Common)。在铣磨和研磨阶段,可以通过三坐标检测获得被测反射镜镜面上离散点的数据点值,通过面形重构获得反射镜的面形误差,但超大口径三坐标测量仪价格非常昂贵,且大型三坐标测量仪精度难以保证,其检测精度仅能达到微米量级[2-3]。为此,可以使用另一种常见的坐标测量机器——摆动臂轮廓仪来进行面形的检测,与传统三坐标相比,其可以将轮廓测量的精度提高一个数量级以上,同时采用连续测量,测量点的密度及测量速度有大幅度的提高,可以很好地满足大口径轮廓测量的需要[4],但是其精度约为0.1 μm RMS。大口径反射镜在抛光过程时,其表面面形可以利用瑞奇康芒法进行测试[5-7],但该方法的瑞奇角产生检测图像压缩,使得被测平面的面形误差与球面本身的像散混合在一起,影响检测结果的准确性;同时大口径反射镜在检测系统的光瞳面上发生变形,使镜面面形误差与系统出瞳光程差之间的对应关系变得复杂,随着瑞奇角的不同而不同,导致检测结果的数据处理变得十分复杂。此外,瑞奇康芒法需要大口径的辅助球面镜,辅助球面镜的口径一般为待测大口径反射镜镜面尺寸的1.2~1.3倍,大口径、高精度球面反射镜的制造相当困难,因此,大口径平面反射镜的瑞奇康芒法检测很难实现。
为了解决大口径反射镜面形检测的困难,尤其是针对30 m望远镜第三反射镜高精度测量的瓶颈,提出了五棱镜扫描和子孔径拼接相结合的方法。五棱镜扫描检测获取低阶面形精度,配合子孔径进行全口径测试,获得全频段面形信息。
-
当大口径反射镜经过研磨和抛光后,可以用五棱镜扫描对其进行测试[8],该方法的原理示意图如图1所示,不管入射光束的角度如何,五棱镜都能使光线偏转90°。入射光束与出射光束之间的夹角仅取决于测试表面的倾斜程度,因此可以通过五棱镜扫描测量光学表面的斜率。
五棱镜扫描测试流程图如图2所示。经五棱镜扫描后可以获得反射镜表面的斜率信息,对斜率数据进行积分能够得到其面形误差,经过像差拟合后,如果低阶像差不明显,反射镜将进行精抛光并结合子孔径拼接测试技术开展后续加工;如果低阶像差明显,该反射镜将重新进行研磨和粗抛光加工,并通过五棱镜扫描再次测量。
大口径反射镜子孔径拼接检测流程如图3所示,其方法原理如图4所示。首先进行子孔径规划,大口径反射镜将被划分成几个小的区域(子孔径),利用小口径干涉仪逐次对大口径反射镜的各区域进行相位测量,利用三角剖分算法和综合优化子孔径拼接算法可以完成对大口径反射镜全口径的面形检测 [9-13]。
首先,根据被测镜和干涉仪的孔径确定子孔径的大小和个数。通过拼接系数可以拼接整个面形,为了正确计算拼接系数,相邻子孔径会有重叠区域,且重叠区域大于子孔径的1/4。其次,通过调整干涉仪与大口径反射镜,使干涉仪对准和遍历各子孔径,利用干涉测量法可以检测各子孔径的相位分布,通过拼接获得全口径面形。在测量每个子孔径面形相位的过程中,由于受测试环境的影响、镜面粗糙度不均匀、测试仪器精度等因素的影响,使得子孔径数据必然存在缺失,将采用三角剖分算法求解和填补子孔径相位数据[14]。
假设将每个子孔径的相位数据转换为对应的全局三维坐标,如图5 (a)所示,首先在X-Y平面上定义均匀网格,将第i和第j个子孔径的相位数据投影到X-Y平面上,三角剖分示意图如图5 (b)所示。对子孔径i的测量数据进行三角剖分,从而得到预先定义的数据点(x, y)的面形值,即利用三点定义一个平面,平面方程为ax+by+z+c=0 (z的系数一定不为0),利用公式(1)求解出系数a,b,c:
$$ \left( {\begin{array}{*{20}{c}} {{x_1}}&{{y_1}}&1 \\ {{x_2}}&{{y_2}}&1 \\ {{x_3}}&{{y_3}}&1 \end{array}} \right)\left( {\begin{array}{*{20}{c}} a \\ b \\ c \end{array}} \right) = \left( {\begin{array}{*{20}{c}} { - {\textit{z}_1}} \\ { - {\textit{z}_2}} \\ { - {\textit{z}_3}} \end{array}} \right) $$ (1) 式中:(xi, yi, zi)( i=1~3)为构成基三角形的三个点的坐标。得到这三个点所确定的平面方程为公式(2):
$$ \textit{z} = - ax - by - c $$ (2) 将点(x, y)的坐标代入到该平面方程中,即可得到对于第i个孔径,该插值点的面形值z。
通过逐次拼接相邻子孔径可以得到整个反射镜的面形误差,但往往会引起拼接误差的累积。文中提出了一种综合优化拼接方法,可以同时计算所有拼接系数。由于拼接平面,各子孔径的未对准误差仅有相对平移量和倾斜量的组合,可以描述为:
$$ w_i^{'} = {w_i} + {p_i} + {a_i}{x_i} + {b_i}{y_i} $$ (3) 式中:wi为干涉仪测量的第i个子孔径的相位数据;wi'为第i个子孔径的实际相位数据;pi为相对的平移系数;ai和bi分别为相对X方向倾斜和Y方向倾斜的系数。
假设重叠区域数为N,使得所有重叠区域相位差的平方和值为最小,用最小二乘法计算最佳拼接系数:
$$ \begin{split} \mathrm{min}={\displaystyle \sum _{i=1}^{N}{\displaystyle \sum _{\begin{array}{l}_{i=1}\\ _{j\ne i}\end{array}}^{N}{\left[\left({w}_{i}+{p}_{i}+{a}_{i}{x}_{i}+{b}_{i}{y}_{i}\right)-\left({w}_{j}+{p}_{j}+{a}_{j}{x}_{j}+{b}_{j}{y}_{j}\right)\right]}^{2}}} \\[-10pt] \end{split}$$ (4) wi和wj为随机相邻子孔径的相位数据,通过三角剖分插值,(xi,yi)和(xj,yj)的坐标相同。因此,公式(4)可以简化为:
$$ \min = \displaystyle\sum\limits_{i = 1}^N {\displaystyle\sum\limits_{\scriptstyle j = 1\atop \scriptstyle j \ne i} ^N {{{\left[\Delta w + \Delta p + \Delta ax + \Delta by\right]}^2}} } $$ (5) 通过最小二乘拟合可以计算出所有的拼接系数,从而精确地拼接出完整的面形信息。
-
当30 m望远镜第三反射镜经过研磨和粗抛光后,镜面的面形误差较小,需要对其面形进行光学检测,为此制定了专门的检测方案规划。
-
TMT项目要求的低阶面形误差检测流程如图6所示。首先扫描出被检面域的低阶信息,然后进行坐标变化处理,将待检面域处理为圆形,拟合圆形面域下的离焦和像散,代入计算公式计算P值,低阶评价指标要求P≤PA=18.8 mas,需要注意的是:mas为TMT项目组自定义的量纲,代表镜面的平滑程度,P值的大小与离焦和像散的RMS值有关,具体计算为:
$$ \begin{split} &\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;P = {C_f}\left| {{Z_4}} \right| + {C_a}\sqrt {{Z_5}^2 + {Z_6}^2} \\ &{C_f} = 1.97\;{\rm{mas}}/(RMS@{\text{μm})},{C_a} = 63.7\;{\rm{mas}}/{(RMS@\text{μm})} \end{split} $$ (6) 式中:
${Z_4}$ 为低阶离焦项;${Z_5}$ 和${Z_6}$ 分别表示0°像散和45°像散。首先,通过五棱镜扫描进行面形测试,如图7所示,其有两个五棱镜,一个棱镜在被测镜上扫描提供测量A,另一个保持固定提供参考测量B,因此大口径反射镜表面的斜率信息可通过(A−B)/2获得[15-16],通过扫描多条母线,可以计算出大反射镜的斜率误差,并通过积分计算出表面的相位分布W(X, Y)。因为面形误差的主要分布是低阶像差,W可以描述为:
$$ \begin{split} {W} = &{A_1} + {A_2}X + {A_3}Y + {A_4}({X^2} + {Y^2}) + {A_5}XY + {A_6}({X^2} - {Y^2}) + \\ &{A_7}X({X^2} + {Y^2}) + {A_8}Y({X^2} + {Y^2}) + {A_9}{({X^2} + {Y^2})^2} + \delta \\[-10pt] \end{split}$$ (7) 式中:Ai (i=1~9)为每种初级像差的系数,它分别表示平移、X倾斜、Y倾斜、离焦、像散(A5、A6)、彗差(A7、A8)和球差系数,δ为面形误差的残余像差,通过最小二乘法,可利用公式(8)拟合最佳系数,因此,可以获得面形误差的低阶像差并循环加工,直到低阶像差的分布满足TMT项目需求。
$$ \begin{split} &\sum {{{\left\{ {W - \left[{A_1} + {A_2}X + {A_3}Y + {A_4}({X^2} + {Y^2}) + {A_5}XY + {A_6}({X^2} - {Y^2}) +\right.}\right.}}} \\ &{{{\left.{\left.{A_7}X({X^2} + {Y^2}) + {A_8}Y({X^2} + {Y^2}) + {A_9}{{({X^2} + {Y^2})}^2}\right]} \right\}}^2} = } {\text{min}}\\[-13pt] \end{split} $$ (8) -
当面形误差的分布主要包含中、高频误差时将通过子孔径拼接进行面形检测。因为30 m望远镜第三反射镜是一个3.5 m×2.5 m的大口径反射镜,如果用小口径干涉仪进行拼接,子孔径数达到上万个,测试和计算时间很长,拼接精度很难保证[17-19]。
为了解决这些问题,30 m望远镜第三反射镜将采用Fizeau干涉仪进行拼接,设备示意图如图8所示,它包含一个小口径干涉仪、一个离轴抛物镜(Off Axis Paraboloid, OAP)和一个口径为1.5 m的参考镜,参考波面和经反射镜返回的测试波面将发生干涉,子孔径的数量大大减少,仅需九个子孔径即可覆盖整个反射镜,子孔径的分布如图9所示。通过开发的拼接算法可以精确地得到全口径的面形误差。
图 8 利用1.5 m参考镜进行拼接干涉检测30 m望远镜第三反射镜示意图
Figure 8. Sketch of testing TMT-M3 by Fizeau interferometry with 1.5 m reference mirror
-
为了解决大口径反射镜,特别是30 m望远镜第三反射镜测试的困难,提出了一种五棱镜扫描与子孔径拼接相结合的新方法。介绍和研究了五棱镜扫描和子孔径拼接计测的基本理论和原理,制定了30 m望远镜第三反射镜的光学测试方案,并研制了30 m望远镜第三反射镜的原理镜进行实验验证,其最终面形误差低阶值为6.61 mas,slopeRMS为0.97 μrad,满足设计指标要求,从而提供了一种高精度光学检测大口径反射镜面形的方法。
Optical testing of the super-large plane mirror (Invited)
-
摘要: 在简要总结了各种检测大口径反射镜难点的基础上,为了实现30 m望远镜(TMT)超大口径第三反射镜的高精度检测,提出了一种融合五棱镜扫描技术和子孔径拼接测试技术的新方法。大口径反射镜分阶段依次进行了五棱镜扫描测试和子孔径拼接检测,对该技术的基本原理和基础理论进行了分析和研究,制定了检测30 m望远镜第三反射镜(口径为3.5 m×2.5 m)的方案,对其测试流程、五棱镜设计、五棱镜扫描像差拟合、拼接最优化算法等进行了详细分析,并对30 m望远镜第三反射镜的原理镜进行了实验验证,其最终拼接检测面形的均方根值(RMS)和斜率均方根值(slopeRMS)分别为28.676 nm和0.97 μrad。Abstract: On the basis of summing up the difficulties of testing large flat mirror, for the sake of testing super-large surface especially for the tertiary mirror of Thirty Meter Telescope (TMT-M3) with precise method, a new technology syncretized pentaprism scanning method and Subaperture Stitching Interferometry (SSI) was proposed. The large mirror was tested by pentaprism scanning and SSI successively, the basic principle of the metrology was introduced and researched, and a reasonable optical testing project of TMT-M3 with the aperture of 3.5 m×2.5 m was made, the flow chart, pentaprism design, aberration fitting of pentaprism scanning, and the optimum stitching algorithm were analyzed particularly. And a prototype of TMT-M3 was developed for experimental verification. The Root Mean Square (RMS) and slope Root Mean Square (slopeRMS) values are 28.676 nm and 0.97 μrad, respectively.
-
Key words:
- optical testing /
- super-large plane mirror /
- TMT /
- SSI /
- pentaprism scan
-
-
[1] TMT International Observatory. Thirty Meter Telescope Astronomy's Next-Generation Observatory [EB/OL]. [2022-01-12]. https://www.tmt.org/. [2] Pan Junhua. Design, Fabrication and Testing of Optical Aspheres [M]. Beijing: Science Press, 1994. (in Chinese) [3] Wang Quandou. Novel profilometer with dual digital length gauge for large aspherics measurement [C]//SPIE, 2000, 4231: 39-46. [4] Su P, Oh C J, Parks R E. et al. Swing arm optical CMM for aspherics [C]//SPIE, 2009, 7426: 74260J. [5] Han S, Novak E, Schuring M. Application of Ritchey-Common test in large flat measurements [C]//SPIE, 2001, 4399: 131-136. [6] Han S, Novak E, Schurig M. Ritchey-common test used for measurement of astronomical optical [C]//SPIE, 2003, 4842: 270-273. [7] Shu K L. Ray-trace analysis and data reduction methods for the Ritchey-Common test [J]. Applied Optics, 1983, 22(12): 1879-1886. doi: 10.1364/AO.22.001879 [8] Mallik P, Zhao C, Burge J H. Measurement of a 2-m flat using a pentaprism scanning system [J]. Optical Engineering, 2007, 46(2): 023602. doi: 10.1117/1.2700386 [9] Murphy P, Fleig J, Forbes G. Subaperture stitching interferometry for testing mild aspheres [C]//SPIE, 2006, 6293: 62930J. [10] Yan Lisong, Zhang Binzhi, Wang Xiaokun, et al. Subaperture stitching testing to flat mirror based on weighting algorithm [J]. Infrared and Laser Engineering, 2021, 50(11): 20210520. (in Chinese) [11] Cai Zhihua, Wang Xiaokun, Hu Haixiang, et al. Non-null stitching test convex aspheric metal mirror [J]. Infrared and Laser Engineering, 2021, 50(11): 20210061. (in Chinese) [12] Kulawiec A, Murphy P, Marco M. Measurement of high-departure aspheres using subaperture stitching with the variable optical null (VON) [C]//SPIE, 2010, 7655: 765512. [13] Zhao C Y, Burge J H. Stitching of off-axis sub-aperture null measurements of an aspheric surface [C]//SPIE, 2008, 7063: 706316. [14] Yan Lisong, Wang Xiaokun, Zheng Ligong, et al. Experimental study on subaperture testing with iterative triangulation algorithm [J]. Opt Express, 2013, 21(19): 22628-22644. doi: 10.1364/OE.21.022628 [15] Qi Erhui, Luo Xiao, Li Ming, et al. Error analysis of scanning pentaprism system in optical testing on large aperture flat mirror [J]. Infrared and Laser Engineering, 2015, 44(2): 639-646. (in Chinese) doi: 10.3969/j.issn.1007-2276.2015.02.044 [16] Geckeler R D. Optimal use of pentaprisms in highly accurate deflectometric scanning [J]. Measurement Science and Technology, 2007, 18(1): 115-125. doi: 10.1088/0957-0233/18/1/014 [17] Su Peng. Absolute measurements of large mirrors [D]. Tucson: The University of Arizona, 2008. [18] Wang X K, Wang L H. Measurement of large aspheric surfaces by annular subaperture stitching interferometry [J]. Chinese Optics Letters, 2007, 11(5): 645-647. [19] Zhang P, Zhao H, Zhou X, et al. Subaperture stitching interferometry using stereovision positioning technique [J]. Opt Express, 2010, 18(14): 15216-15222. doi: 10.1364/OE.18.015216 -