-
非球面相比球面有更多的自由度,可以有效地降低光学系统的像差、体积和质量 [1-3] 。因此,非球面被越来越广泛地用于先进的光学系统中,例如空间望远镜[4]、光刻机[5]、星载遥感相机[6]、星载光谱仪[7]和激光武器[8]等领域。
非球面参数是表征一个非球面的整体形状的基本变量。对于回转对称非球面,非球面的表达式[3]可以表示为:
$$ z = \frac{{{r^2}}}{{R + \sqrt {{R^2} - (1 + K){r^2}} }} + \sum\limits_{i = 2}^n {{A_{2i}}{r^{2i}}} $$ (1) 式中:z为非球面在对称轴z轴方向的矢高;
$r = \sqrt {{x^2} + {y^2}} $ 为非球面上的点(x, y, z)距离z轴的径向距离;R为非球面的顶点曲率半径;K为非球面的圆锥系数;A2i为高阶非球面系数。其中R、K和A2i即为非球面参数,决定了非球面的面型,表征了非球面的基本光学性质。随着光学设计和加工水平的提高,包含高阶项的高阶非球面逐步得到应用,对上述非球面参数的测量成为光学系统加工和装调过程中的迫切需求。首先,非球面的顶点曲率半径R是最重要的参数,学者们基于各种原理对R的测量方法进行了研究。清华大学的王浩等人[9]采用光线追迹法测量了一个抛物面的R;浙江大学的杨甬英等人[10]利用部分零补偿干涉法,通过多次等步长移动部分补偿镜与被测镜之间的间距求解出补偿距离,进而求解出R;亚拉巴马大学亨维尔分校的Ying Pi等人[11]利用二次非球面离轴三共焦位置的干涉图特征测量R;北京理工大学的杨佳苗等人[12]采用差动共焦方法测量了抛物面的R。这些方法假设非球面的圆锥系数K不存在误差,当K存在误差时,求解存在偏差。
为了消除K的影响,一些学者研究了可以同时测量R和K的方法。墨西哥国立自治大学的Rufino Díaz-Uribe等人[13]采用节点工作台法,对一个椭球面的R和K进行了测量;哥伦比亚国立大学的Yobani Mejía等人[14]采用一系列Hartmann图像测量了一个近似抛物面的顶点曲率半径和K;中国科学院长春光学精密机械与物理研究所的陈新东等人[15]采用零补偿干涉法定位被测面,采用激光跟踪仪测量光路间隔,以求解二次非球面的顶点曲率半径和K;中国科学院长春光学精密机械与物理研究所的熊玲等人[16]采用摆臂式轮廓测量方法扫描被测面获取被测面表面轮廓点云分布,再通过拟合计算非球面参数,实现了顶点曲率半径和K的同时测量,不过使用的计算模型暂不适用于高阶非球面。
对R、K和A2i的同时测量是高阶非球面在将来广泛应用时不可避免的课题。已有一些方法对其进行了探索。采用三坐标测量的点云数据进行拟合求解非球面参数是最常用的方法,因为三坐标测量机具有测量精度高和自由度高的两大优势,在测量领域中得到了广泛的使用。在拟合方法中,拟合算法因为耦合问题往往存在一定的求解偏差,需要特殊的拟合算法。例如新加坡南洋理工大学的Zhaomin Wang等人[17]采用更高阶的分量对R和K进行求解,再求解需要的A2i,但该方法没有考虑更高阶分量的非球面参数误差可能对R和K造成的影响,会导致R和K的求解偏差,进而导致高阶非球面参数A2i的求解偏差。
笔者在前期的工作中[18-21]关注了部分补偿干涉法中非球面参数误差与剩余波前,以及补偿距离即补偿镜与被测非球面间距之间的关系,构建了剩余波前、补偿距离和被测非球面参数之间的数学模型,并提出虚实结合迭代算法实现了高阶非球面参数误差的求解。其中,补偿距离的测量精度对被测非球面的测量精度有着直接的影响。以下两种方式可获得补偿距离:第一种[18-20]是针对部分补偿镜设计消球差透镜组,将补偿镜从猫眼位置向实际被测位置移动,通过导轨记录这段位移,进而求解补偿距离。该方式操作较为复杂:部分补偿镜与消球差透镜组相接触,在安装并取下消球差透镜组过程中易导致补偿镜姿态发生变化;被测镜在测量过程中也需要进行移动,对导轨的精度有着较高的要求。第二种[21]是采用差动共焦的方案实现补偿距离的测量,但是系统非常复杂,在部分补偿干涉系统之外还需要复杂的差动共焦测距系统。
为了实现对高阶非球面参数误差的高精度测量,同时简化测量过程和测量装置,文中提出了一种基于虚实结合干涉仪的非球面参数误差测量系统,该系统在部分补偿干涉法测量剩余波前的基础上,复用干涉仪利用轴向扫描透镜的猫眼波前定位方法测量补偿距离,结合虚实结合迭代算法求解出被测非球面的参数误差。该系统在保证测量精度的同时,系统复杂程度相比之前采用差动共焦的方案得到了有效的降低。
-
实验系统如图3所示。测量系统包括部分补偿干涉系统(partial compensated interferometry system, PCIS)和补偿距离测量系统,而补偿距离测量系统包含猫眼波前定位系统(cat’s-eye wavefront positioning system)和测量CL轴向位移的迈克尔逊激光外差干涉仪。其中猫眼波前定位系统复用了部分补偿干涉系统,仅多出了一个CL,且该CL是可安装与拆除的,因此系统结构简单,易于装调且适用性广。
当未安装CL时,PC、SUT与干涉仪构成部分补偿干涉系统,用于测量系统最佳补偿波前。实验所用的干涉仪是Zygo Dynafiz 4 in (1 in=2.54 cm)动态干涉仪。调节PC与SUT之间的间距,当干涉仪采集到的剩余波前的PV值最小时,固定PC与SUT,并测量获得此时的剩余波前,即最佳补偿波前。
之后在平移台上安装CL,置于干涉仪出射的准直光束中,CL会将准直光会聚,当会聚光会聚在PC后表面或SUT表面附近时,干涉仪中可以采集到反射回干涉仪的猫眼波前形成的干涉图。根据猫眼波前PV值最小特性定位位置1和位置2。同时平移台上安装有一个角锥棱镜,该棱镜为迈克尔逊激光外差干涉测距系统的一部分,可记录CL移动的位移。通过迈克尔逊双频激光干涉测距系统,可以高精度地测量出位置1到位置2之间的位移。激光外差测距系统中使用的双频激光器为PT-1105 C,数字采集卡的型号为PT-1313 B。
SUT的名义参数,即参数的设计值如表1所示,是一个凸的四阶非球面。
表 1 被测非球面的名义参数
Table 1. Nominal parameters of SUT
Parameter R/mm K A4 Designed value 35.8256 −0.6291 1.4398×10−7 PC的设计参数如表2所示,采用一个双胶合透镜,最后一面没有镀增透膜以保证采用猫眼波前测量补偿距离时测量光返回到干涉仪。
表 2 PC的设计参数
Table 2. Designed parameters of PC
Surface 1 2 3 Radius/mm 58.86 −41.54 −1074.58 Thickness/mm 15.3 20 — Aperture/mm 46 50.8 — Glass ZF4 ZF7 — CL采用一个双胶合透镜,透镜的标定参数如表3所示。
表 3 会聚透镜的参数
Table 3. Parameters of CL
Surface 1 2 3 Radius/mm 92.233 −66.841 −196.175 Thickness/mm 4 2.5 — Aperture/mm 25.4 25.4 — Glass H-K9L H-ZF2 — -
在实际部分补偿干涉系统中测量实际的剩余波前,当剩余波前PV值最小时,空间坐标归一化的剩余波前如图4(a)所示。在虚拟干涉仪中仿真获得虚拟的剩余波前,归一化的虚拟剩余波前如图4(b)所示。实际剩余波前与虚拟剩余波前的偏差如图4(c)所示。因为回转对称非球面的参数只会产生回转对称的像差,非回转对称像差由随机面形误差造成。对实际剩余波前与虚拟剩余波前的偏差进行泽尼克分解,只保留回转对称的面形误差,以消除非回转对称的面形误差,可以获得如图4(d)所示的回转对称的剩余波前偏差。
图 4 剩余波前与偏差。(a) 实际剩余波前;(b) 初始虚拟剩余波前;(c) 实际与初始虚拟剩余波前的偏差;(d) 偏差中的回转对称部分
Figure 4. Residual wavefront and deviation. (a) Real compensation wavefront; (b) Initial virtual compensation wavefront; (c) Deviation between the real and initial virtual compensation wavefront; (d) Rotationally symmetric component of the deviation
在猫眼波前测量系统中,会聚光在PC后表面与SUT表面附近反射回干涉仪的猫眼波前可被动态干涉仪实时监测,同时迈克尔逊双频激光干涉仪可以检测PC的位置。图5为猫眼波前在位置1与位置2附近时PV值与会聚透镜位置之间的关系。位置1为会聚光会聚在SUT表面附近时的情况,位置2为会聚光会聚在PC后表面附近时的情况。
通过位置1与位置2之间的差值,可以获得采用猫眼波前探测时会聚光焦点从PC后表面移动至PC表面所需要的位移l。再通过虚拟干涉系统,利用这个位移求解出实际最佳补偿距离dr。进行了10次测量的测量结果如表4所示。
表 4 补偿距离测量结果
Table 4. Measured result of compensation distances
No. l/mm dr/mm 1 45.1666 16.1898 2 45.1688 16.1904 3 45.1714 16.1910 4 45.1639 16.1892 5 45.1626 16.1888 6 45.1622 16.1887 7 45.1588 16.1879 8 45.1545 16.1868 9 45.1568 16.1874 10 45.1604 16.1883 10次测量的平均值为16.1888 mm,当置信度系数
$k = 3$ 时,平均值的A类不确定度为0.0013 mm。利用参考文献[21]中的求解模型求解,迭代过程结果如图6所示。特别注意的是,CL和PC与设计值均存在制造误差,需要对CL和PC进行标定,标定的内容有CL和PC的各面曲率、透镜厚度以及折射率。实际求解时使用的是标定结果。
迭代计算获得的非球面参数误差如表5所示。表5中将松下公司三坐标测量仪UA3P测量的表面轮廓采用高斯-牛顿法[22]进行拟合的结果作为实际值,用以比较求解的非球面参数。通过虚实结合迭代算法求得的顶点曲率半径为35.7482 mm,相对误差为2.321‱;圆锥系数为−0.6286,相对误差为7.94 8‱;四阶非球面系数为−4.883×10−7,相对误差为2.995%。
表 5 非球面参数求解结果
Table 5. Calculated result of aspheric surface parameters
Parameter R/mm K A4 Real value by UA3P 35.7569 −0.6291 −4.741×10−7 Proposed method 35.7482 −0.6286 −4.883×10−7 Error −0.0083 0.0005 −1.42×10−8 Relative error 2.321‱ 7.948‱ 2.995% 最终获得的虚拟波前如图7(b)所示,虚拟波前与图7(a)中的实际波前的偏差如图7(c)所示,消除非回转对称的面形误差、只保留回转对称的面形误差如图7(d)所示。其中图7(c)所示的波前偏差相比初始的图4(c)中的波前偏差有所降低,且波前偏差中主要是非回转对称的偏差。
图 7 虚实结合干涉仪剩余残差。(a) 实际剩余波前;(b) 优化后虚拟剩余波前;(c) 实际与优化后虚拟剩余波前的偏差;(d) 偏差中的回转对称部分
Figure 7. Residual error in virtual-real combination interferometer. (a) real compensation wavefront; (b) optimized virtual compensation wavefront; (c) deviation between the real and optimized virtual compensation wavefront; (d) rotationally symmetric component of the deviation
-
在求解非球面参数误差的过程中,补偿距离是一个重要参数,直接影响非球面参数的测量精度。在求解的过程中采用10次补偿距离测量的平均值16.1888 mm,由表4中的数据可以计算出置信度为3时的A类不确定度为0.0013 mm。为分析其影响,将补偿距离设置为16.1875~16.1901 mm,间距为0.0001 mm,分别求解非球面参数,得到的结果如图8所示。
从图8可以看出,R与补偿距离成近似的线性关系,补偿距离的测量精度直接决定了R的测量精度,当补偿距离的不确定度为0.0013 mm时,R从35.7495 mm线性减小至35.7469 mm,最大偏差相较测量结果35.7482 mm的偏差为0.0013 mm,可认为此次测量结果R的不确定度为0.0013 mm。K和A4与补偿距离没有明显的线性关系,但是K和A4在此次测量中呈现出明显的负相关关系。公式(2)为含有四次非球面参数误差的公式[21]:
$$ \Delta {D_4} = \frac{{(K + \Delta K + 1)}}{{8{{\left( {R + \Delta R} \right)}^3}}} - \frac{{(K + 1)}}{{8{R^3}}} + \Delta {A_4} $$ (4) 式中:ΔD4为非球面面形偏差的幂级数拟合的四阶项;ΔR为顶点曲率半径误差;ΔK为圆锥系数误差;ΔA4为非球面的四阶参数误差。由此可见,ΔK和ΔA4确实成负相关。
经分析可知,K的不确定度为0.0002,A4的不确定度为6×10−10。
除去上述补偿距离测量误差的影响,实际测量结果与作为真值的UA3P的测量结果仍存在一定偏差。造成该偏差的可能原因有两个:一是测距的迈克尔逊干涉仪的阿贝误差和余弦误差造成的测量补偿间距时的系统误差,需要对导轨的直线度、导轨移动方向与光轴之间的夹角进行标定;二是非回转对称的装调误差与非回转对称的像差造成的系统误差,可利用曲率半径接近的标准球面镜进行标定,为了消除PC非旋转对称面形误差的影响,需要对其面形进行标定。在高精度测量应用中,上述误差可通过标定进一步减小,拟在后续的研究中开展。
Virtual-real combination interferometry for aspheric surface parameter error measurement with cat-eye-wavefront positioning (invited)
-
摘要: 非球面广泛应用于光学系统中。表征非球面的基本参数包括顶点曲率半径、圆锥系数以及高阶非球面系数,它们贯穿了非球面的设计、制造、检测、装调过程。对非球面参数的高精度测量是加工和装调的前提。提出了一种基于虚实结合干涉仪的部分补偿非球面参数误差测量系统。在该系统中,采用部分补偿干涉法测量剩余波前,复用干涉仪利用轴向扫描透镜的猫眼波前定位法测量补偿镜与被测镜之间的间距,采用虚实结合迭代算法进行参数误差求解。该系统仅需在部分补偿干涉光路中引入会聚透镜,装调简单,测量精度高。通过一个四次非球面参数误差测量实验验证了该方法的有效性和精度。
-
关键词:
- 仪器科学与技术 /
- 非球面参数误差测量系统 /
- 虚实结合干涉法 /
- 猫眼波前定位
Abstract: Aspheric surfaces are widely used in optical systems. Aspheric surface parameters, including vertex radius of curvature, conic constant, and high-order aspheric coefficients participate in optical design, manufacture, measurement, alignment and assembly. Precision measurement for aspheric surface parameters is the basis for manufacture, alignment and assembly. A partial compensation measurement system is proposed based on virtual-real combination interferometer for aspheric surface parameter error measurement. In this measurement system, residual wavefronts are measured by partial compensation interferometry. The interferometer is multiplexed in cat-eye-wavefront positioning method to measure the compensation distance which is the distance between compensator and aspheric surface under test. Aspheric surface parameter errors are calculated by virtual-real combination iterative algorithm. The system only needs to introduce a converging lens into the optical path of partial compensation interferometry, which is easy to align and assemble, and has high measurement accuracy. The effectiveness and accuracy are verified through the measurement of a 4th-order aspheric surface. -
图 4 剩余波前与偏差。(a) 实际剩余波前;(b) 初始虚拟剩余波前;(c) 实际与初始虚拟剩余波前的偏差;(d) 偏差中的回转对称部分
Figure 4. Residual wavefront and deviation. (a) Real compensation wavefront; (b) Initial virtual compensation wavefront; (c) Deviation between the real and initial virtual compensation wavefront; (d) Rotationally symmetric component of the deviation
图 7 虚实结合干涉仪剩余残差。(a) 实际剩余波前;(b) 优化后虚拟剩余波前;(c) 实际与优化后虚拟剩余波前的偏差;(d) 偏差中的回转对称部分
Figure 7. Residual error in virtual-real combination interferometer. (a) real compensation wavefront; (b) optimized virtual compensation wavefront; (c) deviation between the real and optimized virtual compensation wavefront; (d) rotationally symmetric component of the deviation
表 1 被测非球面的名义参数
Table 1. Nominal parameters of SUT
Parameter R/mm K A4 Designed value 35.8256 −0.6291 1.4398×10−7 表 2 PC的设计参数
Table 2. Designed parameters of PC
Surface 1 2 3 Radius/mm 58.86 −41.54 −1074.58 Thickness/mm 15.3 20 — Aperture/mm 46 50.8 — Glass ZF4 ZF7 — 表 3 会聚透镜的参数
Table 3. Parameters of CL
Surface 1 2 3 Radius/mm 92.233 −66.841 −196.175 Thickness/mm 4 2.5 — Aperture/mm 25.4 25.4 — Glass H-K9L H-ZF2 — 表 4 补偿距离测量结果
Table 4. Measured result of compensation distances
No. l/mm dr/mm 1 45.1666 16.1898 2 45.1688 16.1904 3 45.1714 16.1910 4 45.1639 16.1892 5 45.1626 16.1888 6 45.1622 16.1887 7 45.1588 16.1879 8 45.1545 16.1868 9 45.1568 16.1874 10 45.1604 16.1883 表 5 非球面参数求解结果
Table 5. Calculated result of aspheric surface parameters
Parameter R/mm K A4 Real value by UA3P 35.7569 −0.6291 −4.741×10−7 Proposed method 35.7482 −0.6286 −4.883×10−7 Error −0.0083 0.0005 −1.42×10−8 Relative error 2.321‱ 7.948‱ 2.995% -
[1] Liang Zijian, Yang Yongying, Zhao Hongyang, et al. Advances in research and applications of optical aspheric surface metrology [J]. Chinese Optics, 2022, 15(2): 161-186. (in Chinese) [2] Chen Baohua, Wu Quanying, Tang Yunhai, et al. Processing method of large off-axis aspheric mirror [J]. Optics and Precision Engineering, 2021, 29(5): 1095-1102. (in Chinese) doi: 10.37188/OPE.20212905.1095 [3] Malacara D. Optical Shop Testing[M]. 3rd ed. New York: Wiley-Interscience, 2006. [4] Wang Fengpu, Li Xinnan, Xu Chen, et al. Optical testing path design for LOT aspheric segmented mirrors with reflective-diffractive compensation [J]. Chinese Optics, 2021, 14(5): 1184-1193. (in Chinese) doi: 10.37188/CO.2020-0218 [5] Li Meixuan, Wang Li, Dong Lianhe, et al. Application of aspheric technique in immersion lithography lighting system [J]. Infrared and Laser Engineering, 2017, 46(11): 1118003. (in Chinese) [6] Zhang Chao, Xing Hui, Song Junru, et al. Measurement of optical axis eccentricity of a large-aperture concave ellipsoid mirror [J]. Infrared and Laser Engineering, 2021, 50(12): 20210196. (in Chinese) doi: 10.3788/IRLA20210196 [7] Wang Xin, Liu Qiang, Shu Rong. Application of Schmidt optical system with wide-field of view and fast focal ratio to aerospace imaging spectrometer [J]. Optics and Precision Engineering, 2019, 27(3): 533-541. (in Chinese) doi: 10.3788/OPE.20192703.0533 [8] Lv Chenyang, Zhan Renjun, Cui Ying. Collimating lens design and aspheric optimization method of laser cauterization gun [J]. Infrared and Laser Engineering, 2021, 50(3): 20200219. (in Chinese) doi: 10.3788/IRLA20200219 [9] Wang H, Li Y, Zeng L, et al. A simple ray tracing method for measuring the vertex radius of curvature of an aspheric mirror [J]. Optics Communications, 2004, 232(1): 61-68. [10] Yang Y, Zhang L, Shi T, et al. Determination of aspheric vertex radius of curvature in non-null interferometry [J]. Applied Optics, 2015, 54(10): 2838-2844. doi: 10.1364/AO.54.002838 [11] Pi Y, Reardon P J. Determining parent radius and conic of an off-axis segment interferometrically with a spherical reference wave [J]. Optics Letters, 2007, 32(9): 1063-1065. doi: 10.1364/OL.32.001063 [12] Yang J, Qiu L, Zhao W, et al. Laser differential confocal paraboloidal vertex radius measurement [J]. Optics Letters, 2014, 39(4): 830-833. doi: 10.1364/OL.39.000830 [13] Díaz-Uribe R, Cornejo-Rodriguez A. Conic constant and paraxial radius of curvature measurements for conic surfaces [J]. Applied Optics, 1986, 25(20): 3731-3734. doi: 10.1364/AO.25.003731 [14] Mejía Y, Díaz-Uribe R, Pacheco A L, et al. Measuring conic constant and vertex radius of fast convex conic surfaces from a set of Hartmann patterns [J]. Optics Communications, 2016, 363: 166-175. doi: 10.1016/j.optcom.2015.11.011 [15] Chen Xindong, Li Ruigang. Research on geometric parameter measurement method using laser tracker in null lens asphere testing [J]. Chinese Journal of Lasers, 2015, 42(5): 0508007. (in Chinese) [16] Xiong L, Qi E, Luo X, et al. Stitching swing arm profilometer test for large aperture aspherics [J]. Chinese Optics Letters, 2019, 17(11): 112201. doi: 10.3788/COL201917.112201 [17] Wang Z, Qu W, Asundi A. A simplified expression for aspheric surface fitting [J]. Optik, 2017, 140: 291-298. [18] Hao Q, Li T, Hu Y, et al. Vertex radius of curvature error measurement of aspheric surface based on slope asphericity in partial compensation interferometry [J]. Optics Express, 2017, 25(15): 18107-18121. doi: 10.1364/OE.25.018107 [19] Hao Q, Li T, Hu Y. Partial compensation interferometry for measurement of surface parameter error of high-order aspheric surfaces [C]//2017 International Conference on Optical Instruments and Technology: Optical Systems and Modern Optoelectronic Instruments, 2018, 10616: 106160L. [20] Hao Q, Li T, Hu Y, et al. Partial compensation interferometry measurement system for parameter errors of conicoid surface [J]. Review of Scientific Instruments, 2018, 89(6): 063102. doi: 10.1063/1.5027146 [21] Hao Q, Tao X, Hu Y, et al. Interferometric measurement of high-order aspheric surface parameter errors based on virtual-real combination iterative algorithm [J]. Optics Express, 2021, 29(17): 27014-27030. doi: 10.1364/OE.435252 [22] Sha D, Quan S, Zhu Q, et al. Research on optical aspherical surface’s fitting method [J]. Acta Photonica Sinica, 1995, 24(2): 152-157. (in Chinese)