留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

猫眼波前定位虚实结合干涉非球面参数误差测量 (特邀)

胡摇 陶鑫 郝群

胡摇, 陶鑫, 郝群. 猫眼波前定位虚实结合干涉非球面参数误差测量 (特邀)[J]. 红外与激光工程, 2022, 51(9): 20220497. doi: 10.3788/IRLA20220497
引用本文: 胡摇, 陶鑫, 郝群. 猫眼波前定位虚实结合干涉非球面参数误差测量 (特邀)[J]. 红外与激光工程, 2022, 51(9): 20220497. doi: 10.3788/IRLA20220497
Hu Yao, Tao Xin, Hao Qun. Virtual-real combination interferometry for aspheric surface parameter error measurement with cat-eye-wavefront positioning (invited)[J]. Infrared and Laser Engineering, 2022, 51(9): 20220497. doi: 10.3788/IRLA20220497
Citation: Hu Yao, Tao Xin, Hao Qun. Virtual-real combination interferometry for aspheric surface parameter error measurement with cat-eye-wavefront positioning (invited)[J]. Infrared and Laser Engineering, 2022, 51(9): 20220497. doi: 10.3788/IRLA20220497

猫眼波前定位虚实结合干涉非球面参数误差测量 (特邀)

doi: 10.3788/IRLA20220497
基金项目: 国家自然科学基金(51735002);中国科学院战略性先导科技专项(XDA25020317)
详细信息
    作者简介:

    胡摇,女,副教授,博士,主要从事智能干涉仪、计算成像方面的研究

  • 中图分类号: TH74

Virtual-real combination interferometry for aspheric surface parameter error measurement with cat-eye-wavefront positioning (invited)

Funds: National Natural Science Foundation of China (51735002);Strategic Priority Program of Chinese Academy of Science (XDA25020317)
  • 摘要: 非球面广泛应用于光学系统中。表征非球面的基本参数包括顶点曲率半径、圆锥系数以及高阶非球面系数,它们贯穿了非球面的设计、制造、检测、装调过程。对非球面参数的高精度测量是加工和装调的前提。提出了一种基于虚实结合干涉仪的部分补偿非球面参数误差测量系统。在该系统中,采用部分补偿干涉法测量剩余波前,复用干涉仪利用轴向扫描透镜的猫眼波前定位法测量补偿镜与被测镜之间的间距,采用虚实结合迭代算法进行参数误差求解。该系统仅需在部分补偿干涉光路中引入会聚透镜,装调简单,测量精度高。通过一个四次非球面参数误差测量实验验证了该方法的有效性和精度。
  • 图  1  基于猫眼波前定位的虚实结合部分补偿干涉系统。PC:部分补偿镜;SUT:被测面;CL:会聚透镜

    Figure  1.  Virtual-Real partial compensation interferometry system based on cat-eye-wavefront positioning. PC: partial compensator; SUT: surface under test; CL: convergent lens

    图  2  猫眼波前定位原理图

    Figure  2.  Schematic of cat-eye-wavefront positioning

    图  3  非球面参数误差测量系统

    Figure  3.  Aspheric surface parameter errors measurement system

    图  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

    图  5  猫眼波前随位置变化

    Figure  5.  Deviation of cat’s-eye wavefront with the position

    图  6  参数误差求解迭代结果

    Figure  6.  Iteration results of parameter error calculation

    图  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

    图  8  误差分析

    Figure  8.  Error analysis

    表  1  被测非球面的名义参数

    Table  1.   Nominal parameters of SUT

    ParameterR/mmKA4
    Designed value35.8256−0.62911.4398×10−7
    下载: 导出CSV

    表  2  PC的设计参数

    Table  2.   Designed parameters of PC

    Surface123
    Radius/mm58.86−41.54−1074.58
    Thickness/mm15.320
    Aperture/mm4650.8
    GlassZF4ZF7
    下载: 导出CSV

    表  3  会聚透镜的参数

    Table  3.   Parameters of CL

    Surface123
    Radius/mm92.233−66.841−196.175
    Thickness/mm42.5
    Aperture/mm25.425.4
    GlassH-K9LH-ZF2
    下载: 导出CSV

    表  4  补偿距离测量结果

    Table  4.   Measured result of compensation distances

    No.l/mmdr/mm
    145.166616.1898
    245.168816.1904
    345.171416.1910
    445.163916.1892
    545.162616.1888
    645.162216.1887
    745.158816.1879
    845.154516.1868
    945.156816.1874
    1045.160416.1883
    下载: 导出CSV

    表  5  非球面参数求解结果

    Table  5.   Calculated result of aspheric surface parameters

    ParameterR/mmKA4
    Real value by UA3P35.7569−0.6291−4.741×10−7
    Proposed method35.7482−0.6286−4.883×10−7
    Error−0.00830.0005−1.42×10−8
    Relative error2.321‱7.948‱2.995%
    下载: 导出CSV
  • [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)
  • [1] 李明泽, 侯溪, 赵文川, 王洪, 李梦凡, 胡小川, 赵远程, 周杨.  非球面光学表面缺陷检测技术现状和发展趋势(特邀) . 红外与激光工程, 2022, 51(9): 20220457-1-20220457-20. doi: 10.3788/IRLA20220457
    [2] 程乙轮, 谭逢富, 何枫, 侯再红, 秦来安, 王浩, 黄志刚, 吴德成.  漫反射成像法的激光参数测量系统设计 . 红外与激光工程, 2022, 51(9): 20210921-1-20210921-8. doi: 10.3788/IRLA20210921
    [3] 张誉馨, 黎发志, 闫力松, 严冬, 王超凡, 张宇鹤, 张斌智, 张鑫, 蔡铭宣, 何俊超, 王孝坤.  结合CGH与辅助透镜的长焦距非球面反射镜检测(特邀) . 红外与激光工程, 2022, 51(9): 20220384-1-20220384-6. doi: 10.3788/IRLA20220384
    [4] 吕晨阳, 战仁军, 崔莹.  激光热灼枪准直镜头设计与非球面优化方法 . 红外与激光工程, 2021, 50(3): 20200219-1-20200219-8. doi: 10.3788/IRLA20200219
    [5] 刘岩, 雷柏平, 范斌, 边疆.  基于视觉测量的目标定位技术及结构参数优化 . 红外与激光工程, 2020, 49(S2): 20200191-20200191. doi: 10.3788/IRLA20200191
    [6] 赵冬娥, 王思育, 马亚云, 张斌, 李诺伦, 李沅, 褚文博.  基于涡旋光与球面波干涉的微位移测量研究 . 红外与激光工程, 2020, 49(4): 0413005-0413005-6. doi: 10.3788/IRLA202049.0413005
    [7] 王聪, 陈佳夷, 栗孟娟, 王海超, 李斌.  基于干涉测量的Ф1.3 m非球面反射镜定心 . 红外与激光工程, 2020, 49(1): 0113001-0113001(6). doi: 10.3788/IRLA202049.0113001
    [8] 俞家勇, 卢秀山, 田茂义, 贺岩, 吕德亮, 胡善江, 王延存, 曹岳飞, 黄田橙.  机载激光雷达测深系统定位模型与视准轴误差影响分析 . 红外与激光工程, 2019, 48(6): 606005-0606005(9). doi: 10.3788/IRLA201948.0606005
    [9] 李博, 杨军, 黄楠, 石玉松.  基于激光干涉法测水介质正弦压力的动态校准技术 . 红外与激光工程, 2019, 48(8): 805003-0805003(7). doi: 10.3788/IRLA201948.0805003
    [10] 李家琨, 冯其波, 包传辰, 杨婧, 赵斌陶.  转轴转角定位误差测量方法与误差分析 . 红外与激光工程, 2019, 48(2): 217001-0217001(6). doi: 10.3788/IRLA201948.0217001
    [11] 钟波, 陈贤华, 王健, 周炼, 石琦凯, 邓文辉.  高精度离轴非球面透镜的制造与检测 . 红外与激光工程, 2018, 47(7): 718003-0718003(8). doi: 10.3788/IRLA201847.0718003
    [12] 刘楚, 钟凯, 史杰, 靳硕, 葛萌, 李吉宁, 徐德刚, 姚建铨.  迈克尔逊干涉法精确测量太赫兹频谱及目标速度 . 红外与激光工程, 2018, 47(11): 1117006-1117006(7). doi: 10.3788/IRLA201847.1117006
    [13] 张芳, 高教波, 王楠, 赵宇洁, 吴江辉, 郑雅卫.  变间隙法布里-珀罗干涉式长波红外光谱成像系统 . 红外与激光工程, 2017, 46(3): 318001-0318001(7). doi: 10.3788/IRLA201746.0318001
    [14] 李美萱, 王丽, 董连和, 赵迎.  非球面技术在浸没式光刻照明系统中的应用 . 红外与激光工程, 2017, 46(11): 1118003-1118003(6). doi: 10.3788/IRLA201746.1118003
    [15] 詹翔空, 李政勇, 张伊, 王明宇, 王海洋, 吴重庆.  基于Stocks矢量测量与干涉法的径向偏振光束重建 . 红外与激光工程, 2017, 46(4): 427002-0427002(6). doi: 10.3788/IRLA201746.0427002
    [16] 李树旺, 邵士勇, 梅海平, 饶瑞中.  大气吸收性成分的光热干涉法测量 . 红外与激光工程, 2016, 45(S1): 163-168. doi: 10.3788/IRLA201645.S111002
    [17] 王孝坤.  非零位补偿检验非球面技术 . 红外与激光工程, 2015, 44(9): 2702-2709.
    [18] 王孝坤.  异形口径离轴非球面光学加工与测试技术 . 红外与激光工程, 2014, 43(9): 2959-2963.
    [19] 邓勇, 刘宁, 曹红蓓, 郭龙秋.  基于Nd:YAG激光回馈干涉效应的PZT精密测量技术与系统 . 红外与激光工程, 2014, 43(10): 3434-3438.
    [20] 代雷, 隋永新, 吴迪.  基于激光干涉法的曲率半径精密检测系统 . 红外与激光工程, 2013, 42(8): 2221-2225.
  • 加载中
图(8) / 表(5)
计量
  • 文章访问数:  244
  • HTML全文浏览量:  57
  • PDF下载量:  42
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-07-20
  • 修回日期:  2022-08-30
  • 刊出日期:  2022-09-28

猫眼波前定位虚实结合干涉非球面参数误差测量 (特邀)

doi: 10.3788/IRLA20220497
    作者简介:

    胡摇,女,副教授,博士,主要从事智能干涉仪、计算成像方面的研究

基金项目:  国家自然科学基金(51735002);中国科学院战略性先导科技专项(XDA25020317)
  • 中图分类号: TH74

摘要: 非球面广泛应用于光学系统中。表征非球面的基本参数包括顶点曲率半径、圆锥系数以及高阶非球面系数,它们贯穿了非球面的设计、制造、检测、装调过程。对非球面参数的高精度测量是加工和装调的前提。提出了一种基于虚实结合干涉仪的部分补偿非球面参数误差测量系统。在该系统中,采用部分补偿干涉法测量剩余波前,复用干涉仪利用轴向扫描透镜的猫眼波前定位法测量补偿镜与被测镜之间的间距,采用虚实结合迭代算法进行参数误差求解。该系统仅需在部分补偿干涉光路中引入会聚透镜,装调简单,测量精度高。通过一个四次非球面参数误差测量实验验证了该方法的有效性和精度。

English Abstract

    • 非球面相比球面有更多的自由度,可以有效地降低光学系统的像差、体积和质量 [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为高阶非球面系数。其中RKA2i即为非球面参数,决定了非球面的面型,表征了非球面的基本光学性质。随着光学设计和加工水平的提高,包含高阶项的高阶非球面逐步得到应用,对上述非球面参数的测量成为光学系统加工和装调过程中的迫切需求。

      首先,非球面的顶点曲率半径R是最重要的参数,学者们基于各种原理对R的测量方法进行了研究。清华大学的王浩等人[9]采用光线追迹法测量了一个抛物面的R;浙江大学的杨甬英等人[10]利用部分零补偿干涉法,通过多次等步长移动部分补偿镜与被测镜之间的间距求解出补偿距离,进而求解出R;亚拉巴马大学亨维尔分校的Ying Pi等人[11]利用二次非球面离轴三共焦位置的干涉图特征测量R;北京理工大学的杨佳苗等人[12]采用差动共焦方法测量了抛物面的R。这些方法假设非球面的圆锥系数K不存在误差,当K存在误差时,求解存在偏差。

      为了消除K的影响,一些学者研究了可以同时测量RK的方法。墨西哥国立自治大学的Rufino Díaz-Uribe等人[13]采用节点工作台法,对一个椭球面的RK进行了测量;哥伦比亚国立大学的Yobani Mejía等人[14]采用一系列Hartmann图像测量了一个近似抛物面的顶点曲率半径和K;中国科学院长春光学精密机械与物理研究所的陈新东等人[15]采用零补偿干涉法定位被测面,采用激光跟踪仪测量光路间隔,以求解二次非球面的顶点曲率半径和K;中国科学院长春光学精密机械与物理研究所的熊玲等人[16]采用摆臂式轮廓测量方法扫描被测面获取被测面表面轮廓点云分布,再通过拟合计算非球面参数,实现了顶点曲率半径和K的同时测量,不过使用的计算模型暂不适用于高阶非球面。

      RKA2i的同时测量是高阶非球面在将来广泛应用时不可避免的课题。已有一些方法对其进行了探索。采用三坐标测量的点云数据进行拟合求解非球面参数是最常用的方法,因为三坐标测量机具有测量精度高和自由度高的两大优势,在测量领域中得到了广泛的使用。在拟合方法中,拟合算法因为耦合问题往往存在一定的求解偏差,需要特殊的拟合算法。例如新加坡南洋理工大学的Zhaomin Wang等人[17]采用更高阶的分量对RK进行求解,再求解需要的A2i,但该方法没有考虑更高阶分量的非球面参数误差可能对RK造成的影响,会导致RK的求解偏差,进而导致高阶非球面参数A2i的求解偏差。

      笔者在前期的工作中[18-21]关注了部分补偿干涉法中非球面参数误差与剩余波前,以及补偿距离即补偿镜与被测非球面间距之间的关系,构建了剩余波前、补偿距离和被测非球面参数之间的数学模型,并提出虚实结合迭代算法实现了高阶非球面参数误差的求解。其中,补偿距离的测量精度对被测非球面的测量精度有着直接的影响。以下两种方式可获得补偿距离:第一种[18-20]是针对部分补偿镜设计消球差透镜组,将补偿镜从猫眼位置向实际被测位置移动,通过导轨记录这段位移,进而求解补偿距离。该方式操作较为复杂:部分补偿镜与消球差透镜组相接触,在安装并取下消球差透镜组过程中易导致补偿镜姿态发生变化;被测镜在测量过程中也需要进行移动,对导轨的精度有着较高的要求。第二种[21]是采用差动共焦的方案实现补偿距离的测量,但是系统非常复杂,在部分补偿干涉系统之外还需要复杂的差动共焦测距系统。

      为了实现对高阶非球面参数误差的高精度测量,同时简化测量过程和测量装置,文中提出了一种基于虚实结合干涉仪的非球面参数误差测量系统,该系统在部分补偿干涉法测量剩余波前的基础上,复用干涉仪利用轴向扫描透镜的猫眼波前定位方法测量补偿距离,结合虚实结合迭代算法求解出被测非球面的参数误差。该系统在保证测量精度的同时,系统复杂程度相比之前采用差动共焦的方案得到了有效的降低。

    • 首先介绍虚实结合干涉仪,以及如何利用它实现非球面的面型参数误差的测量。

      虚实结合部分补偿非球面测量系统如图1所示,由实际的部分补偿干涉系统(图1(a)),虚拟的部分补偿干涉系统(图1(b))以及数据处理部分(图1(g))组成。在测量过程中,由实际干涉仪部分补偿镜PC与被测非球面镜SUT构成实际的部分补偿干涉系统,测量获得实际的剩余波前。按照实际干涉仪、PC的标定结果,SUT的设计参数以及补偿距离的设计值构成虚拟的部分补偿干涉系统,仿真获得虚拟的剩余波前。当实际的SUT存在面形误差时,实际补偿距离可能发生变化,实际的剩余波前与虚拟的剩余波前也存在偏差。调整虚拟干涉系统中的补偿距离,并优化SUT的面形,当实际剩余波前与虚拟剩余波前相同时,实际干涉系统中的SUT与虚拟干涉系统中的SUT相同。此时虚拟干涉系统中加载的面形误差即为实际SUT的面形误差。考虑到面型参数误差造成的面形误差对剩余波前的影响是回转对称的,若提取出剩余波前的回转对称误差,并通过修改虚拟干涉仪的非球面参数使得剩余波前的回转误差最小,此时虚拟干涉仪中的非球面参数误差即为实际SUT的非球面参数误差。

      图  1  基于猫眼波前定位的虚实结合部分补偿干涉系统。PC:部分补偿镜;SUT:被测面;CL:会聚透镜

      Figure 1.  Virtual-Real partial compensation interferometry system based on cat-eye-wavefront positioning. PC: partial compensator; SUT: surface under test; CL: convergent lens

      之前的研究[21]表明,上述非球面参数误差求解的过程中需要两个关键测量数据:一是实际部分补偿干涉系统测得的剩余波前;二是某种方法测得的补偿距离。为了提高补偿距离的测量精度和测量稳定性,降低系统复杂性,文中提出猫眼波前定位法,并在实际干涉系统中加入轴向扫描用的会聚透镜CL,具体原理如下。

    • 猫眼波前定位法的原理如图2所示。干涉仪发出的准直激光束经由CL形成会聚光,入射之后的任意表面后发生反射,再次经过CL后返回干涉仪,此时干涉仪测量获得的波前定义为猫眼波前(cat’s-eye wavefront)。如图2(b)所示,以PC的后表面定位为例,当会聚光焦点恰好在该表面时,光束上下翻转近似原路返回,得到的猫眼波前峰谷(peak-to-valley,PV)值最小。如图2(a)或图2(c)所示,当CL发生轴向偏移时,由于焦点偏离待定位面,引入了离焦为主的像差,干涉仪检测到的猫眼波前也会发生变化,猫眼波前的PV值会随之增大。

      图  2  猫眼波前定位原理图

      Figure 2.  Schematic of cat-eye-wavefront positioning

      利用猫眼位置的波前PV达到最小值的特性进行定位,结合测位移装置可以对补偿距离进行测量,其原理如图1 (c)~(f)所示。对于每一个光学表面,猫眼波前都会有相似的表现。利用这一点,可以通过检测猫眼波前PV值最小的位置对CL的后续表面进行精确定位。

      在实际系统中可以通过实际干涉仪测量猫眼波前,调整CL使光会聚在补偿镜最后一面,此时猫眼波前PV值最小,对应CL的位置记为位置1,如图1(c)所示。继续调整CL使光会聚在被测非球面镜,此时猫眼波前PV值再次达到最小,对应CL的位置记为位置2,如图1(e)所示。位置1与位置2之间的间距即为实际干涉系统中CL的移动距离,可以采用测位移装置获取,记为l

      为保证实际干涉仪补偿距离dr与虚拟干涉仪中补偿距离dv一致,进行如下优化:将标定得到的CL准确参数加载进虚拟干涉系统,优化CL与PC之间的间距,将光会聚在PC的最后一面,使得猫眼波前的PV值最小,此时会聚透镜的位置即为虚拟系统中的位置1,如图1(d)所示。之后调整CL的轴向位置向PC靠近距离l,即将CL定位至位置2,如图1(f)所示。此时在虚拟系统中加入SUT,优化SUT轴向位置,使得猫眼波前PV值最小,此时PC与SUT之间的间距即为虚拟系统中的补偿距离dv。且上述过程完全模拟实际测量过程,因此实际干涉系统中的补偿距离dr=dv

    • 结合上述原理,参数误差的测量过程主要分两步,第一步使用干涉仪、PC、SUT组成实际部分补偿干涉系统,调整PC与SUT的间距,使得剩余波前的PV值最小,此时可以认为补偿距离是实际最佳补偿距离,测量此时部分补偿法的剩余波前。第二步将CL置于PC与干涉仪之间,按照1.2节流程测得猫眼位置1和2之间的位移l

      计算的流程主要分两步,第一步是利用位置1和位置2之间的位移l,通过1.2节末尾所述方法求解实际干涉系统中的最佳补偿距离dr。第二步是利用参考文献[21]所述的虚实结合迭代算法求解出被测非球面的面形参数误差,此处不再赘述。

    • 实验系统如图3所示。测量系统包括部分补偿干涉系统(partial compensated interferometry system, PCIS)和补偿距离测量系统,而补偿距离测量系统包含猫眼波前定位系统(cat’s-eye wavefront positioning system)和测量CL轴向位移的迈克尔逊激光外差干涉仪。其中猫眼波前定位系统复用了部分补偿干涉系统,仅多出了一个CL,且该CL是可安装与拆除的,因此系统结构简单,易于装调且适用性广。

      图  3  非球面参数误差测量系统

      Figure 3.  Aspheric surface parameter errors measurement system

      当未安装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

      ParameterR/mmKA4
      Designed value35.8256−0.62911.4398×10−7

      PC的设计参数如表2所示,采用一个双胶合透镜,最后一面没有镀增透膜以保证采用猫眼波前测量补偿距离时测量光返回到干涉仪。

      表 2  PC的设计参数

      Table 2.  Designed parameters of PC

      Surface123
      Radius/mm58.86−41.54−1074.58
      Thickness/mm15.320
      Aperture/mm4650.8
      GlassZF4ZF7

      CL采用一个双胶合透镜,透镜的标定参数如表3所示。

      表 3  会聚透镜的参数

      Table 3.  Parameters of CL

      Surface123
      Radius/mm92.233−66.841−196.175
      Thickness/mm42.5
      Aperture/mm25.425.4
      GlassH-K9LH-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后表面附近时的情况。

      图  5  猫眼波前随位置变化

      Figure 5.  Deviation of cat’s-eye wavefront with the position

      通过位置1与位置2之间的差值,可以获得采用猫眼波前探测时会聚光焦点从PC后表面移动至PC表面所需要的位移l。再通过虚拟干涉系统,利用这个位移求解出实际最佳补偿距离dr。进行了10次测量的测量结果如表4所示。

      表 4  补偿距离测量结果

      Table 4.  Measured result of compensation distances

      No.l/mmdr/mm
      145.166616.1898
      245.168816.1904
      345.171416.1910
      445.163916.1892
      545.162616.1888
      645.162216.1887
      745.158816.1879
      845.154516.1868
      945.156816.1874
      1045.160416.1883

      10次测量的平均值为16.1888 mm,当置信度系数$k = 3$ 时,平均值的A类不确定度为0.0013 mm。

      利用参考文献[21]中的求解模型求解,迭代过程结果如图6所示。特别注意的是,CL和PC与设计值均存在制造误差,需要对CL和PC进行标定,标定的内容有CL和PC的各面曲率、透镜厚度以及折射率。实际求解时使用的是标定结果。

      图  6  参数误差求解迭代结果

      Figure 6.  Iteration results of parameter error calculation

      迭代计算获得的非球面参数误差如表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

      ParameterR/mmKA4
      Real value by UA3P35.7569−0.6291−4.741×10−7
      Proposed method35.7482−0.6286−4.883×10−7
      Error−0.00830.0005−1.42×10−8
      Relative error2.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  误差分析

      Figure 8.  Error analysis

      图8可以看出,R与补偿距离成近似的线性关系,补偿距离的测量精度直接决定了R的测量精度,当补偿距离的不确定度为0.0013 mm时,R从35.7495 mm线性减小至35.7469 mm,最大偏差相较测量结果35.7482 mm的偏差为0.0013 mm,可认为此次测量结果R的不确定度为0.0013 mm。KA4与补偿距离没有明显的线性关系,但是KA4在此次测量中呈现出明显的负相关关系。公式(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非旋转对称面形误差的影响,需要对其面形进行标定。在高精度测量应用中,上述误差可通过标定进一步减小,拟在后续的研究中开展。

    • 文中提出了一个基于虚实结合干涉仪的部分补偿非球面参数误差测量系统。该系统采用部分补偿法测量部分补偿后的剩余波前,采用猫眼波前定位方法测量补偿距离,即补偿镜与非球面被测镜之间的间距。利用剩余波前和补偿距离迭代优化计算出非球面被测镜的参数误差。

      通过实验测量了一个凸的非球面,其顶点曲率半径为35.7482,相对误差为2.321‱;圆锥系数为−0.6286,相对误差为7.948‱;四阶非球面系数为−4.883×10−7,相对误差为2.995%。

      上述原理和系统研究为高阶非球面的面形误差和参数误差一体化干涉测量奠定了技术基础,装置简单,无需移动被测面,有望应用于大口径非球面的参数误差测量。

参考文献 (22)

目录

    /

    返回文章
    返回