留言板

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

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

夏克哈特曼扫描拼接检测平面镜(特邀)

王晶 王孝坤 胡海翔 李凌众 苏航

王晶, 王孝坤, 胡海翔, 李凌众, 苏航. 夏克哈特曼扫描拼接检测平面镜(特邀)[J]. 红外与激光工程, 2021, 50(10): 20210527. doi: 10.3788/IRLA20210527
引用本文: 王晶, 王孝坤, 胡海翔, 李凌众, 苏航. 夏克哈特曼扫描拼接检测平面镜(特邀)[J]. 红外与激光工程, 2021, 50(10): 20210527. doi: 10.3788/IRLA20210527
Wang Jing, Wang Xiaokun, Hu Haixiang, Li Lingzhong, Su Hang. Shack-Hartmann scanning and stiching detection plane mirror (Invited)[J]. Infrared and Laser Engineering, 2021, 50(10): 20210527. doi: 10.3788/IRLA20210527
Citation: Wang Jing, Wang Xiaokun, Hu Haixiang, Li Lingzhong, Su Hang. Shack-Hartmann scanning and stiching detection plane mirror (Invited)[J]. Infrared and Laser Engineering, 2021, 50(10): 20210527. doi: 10.3788/IRLA20210527

夏克哈特曼扫描拼接检测平面镜(特邀)

doi: 10.3788/IRLA20210527
基金项目: 吉林省科技发展计划(20200401065GX);中国科学院青年创新促进会(2019221);国家自然科学基金(61975201);国家自然科学基金面上项目(6207031149)
详细信息
    作者简介:

    王晶,女,博士生,主要从事光学设计与测试技术方面的研究

  • 中图分类号: O439;O436.1

Shack-Hartmann scanning and stiching detection plane mirror (Invited)

  • 摘要: 随着先进光学系统设计与制造的发展,大口径光学系统得到了广泛的应用。然而,大口径平面镜高精度面形的检测手段不足,限制了大口径平面镜的制造与应用。为实现大口径平面反射镜的高精度面形检测,提出一种夏克哈特曼扫描拼接检测平面镜面形的方法。对扫描拼接原理、波前重构算法进行了研究,建立了微透镜阵列成像的数学模型,验证了夏克哈特曼扫描拼接检测原理的可行性。针对一口径为150 mm的平面镜进行了扫描拼接检测实验,拼接得到的全口径面形为0.019λ RMS(λ=635 nm);与干涉检测结果对比,检测精度为0.008λ RMS,结果表明该方法能够实现大口径平面反射镜的高精度检测。
  • 图  1  夏克哈特曼扫描拼接检测平面镜原理示意图

    Figure  1.  Schematic diagram of Shack-Hartmann scanning and splicing plane mirror

    图  2  子孔径中心位置规划示意图

    Figure  2.  Schematic diagram of sub-aperture center location planning

    图  3  单个微透镜成像的傅里叶光学物理过程

    Figure  3.  Fourier optical physics in the imaging of a single microlens

    图  4  单个微透镜成像的仿真结果

    Figure  4.  Simulation results of single microlens imaging

    图  5  微透镜阵列成像的仿真结果

    Figure  5.  Simulation results of microlens array imaging

    图  6  全口径仿真结果。(a) 添加的波像差; (b) 拼接仿真结果; (c) 结果误差

    Figure  6.  Full caliber simulation results. (a) Added wave aberrations; (b) Stitching simulation results; (c) Result error

    图  7  夏克哈特曼扫描拼接检测系统

    Figure  7.  Shack-Hartmann scanning and splicing detection system

    图  8  平面镜边缘位置。(a)上边缘;(b)下边缘;(c) 左边缘;(d) 右边缘

    Figure  8.  Plane mirror edge position. (a) Uupper edge;(b) Lower edge;(c) Left edge; (d) Right edge

    图  9  精调对准过程。(a)未对准;(b)对准

    Figure  9.  Fine alignment process. (a) Misalignment;(b) Alignment

    图  10  全口径斜率检测结果。(a) x方向斜率; (b) y方向斜率

    Figure  10.  Full caliber slope test results. (a) Slope in the x direction; (b) Slope in the y direction

    图  11  (a)全口径面形;(b)干涉仪检测结果

    Figure  11.  (a) Full aperture surface; (b) Interferometer test results

    图  12  残差分布

    Figure  12.  Residual distribution

    表  1  仿真参数

    Table  1.   Simulation parameters

    ParameterValue
    Microlens array diameter/mm 14.61
    Mirror radius /mm 25
    Number of subaperture 37
    Wave length/nm 635
    下载: 导出CSV

    表  2  移动平台量程与精度

    Table  2.   Moving platform range and accuracy

    ParameterValueAccuracy
    x /mm1 0000.01
    y /mm5000.01
    z /mm8000.02
    A90 °4″
    C360 °10″
    下载: 导出CSV

    表  3  实验参数

    Table  3.   Experiment parameters

    ParameterValue
    Detection beam aperture /mm14.6
    Mirror diameter /mm150
    Sampling density17×17
    Scanning step /mm9.1
    Wave length/nm635
    下载: 导出CSV
  • [1] Zhao Xianyu. Research on Hartmann principle based on sub-aperture slope scanning detection and error [J]. Infrared and Laser Engineering, 2019, 48(8): 0813003. (in Chinese) doi:  10.3788/IRLA201948.0813003
    [2] Xiang Quncui, Yong Hengzhao, Yao Quanchu, et al. The large sky area multi-object fiber spectroscopic telescope(LAMOST) [J]. Research in Astron Astrophys, 2012, 12(9): 1197-1242. doi:  10.1088/1674-4527/12/9/003
    [3] Roland GEYL, Jacques PASERI. Optical polishing of the VLT 8.2 m primary mirrors–A report [C]//SPIE, 1996, 2775: 476-479.
    [4] 张林超. 大口径平面反射镜检测方法研究[D]. 长春: 中国科学院长春光学精密机械与物理研究所, 2010.
    [5] Stefan Schulte, Bernd Dörband, Frank Schilike, et al. Interferometric testing of the VLT tertiary mirrors [C]//SPIE, 1999, 3745: 14-24.
    [6] Zhi Ganghan, Lu Lin, Lei Chen, et al. Absolute flatness testing of skip-flat interferometry by matrix analysis in polar coordinates [J]. Applied Optics, 2016, 55(9): 2387-2392. doi:  10.1364/AO.55.002387
    [7] An Qichang, Zhang Jingxu, Fei Yang, et al. Ritchey-Common sparse-aperture testing of the giant steerable science mirror [J]. Applied Optics, 2018, 57(27): 7764-7769. doi:  10.1364/AO.57.007764
    [8] Wang Xiaokun, Zheng Ligong, Zhang Binzhi, et al. Subaperture stitching interferometry for testing a large hyperboloid [J]. Infrared and Laser Engineering, 2009, 38(1): 114-119. (in Chinese) doi:  10.3969/j.issn.1007-2276.2009.01.026
    [9] Guo Fudong, Tang Feng, Lu Yunjun, et al. Rapid adjustment and measurement for subaperture stitching interferometry [J]. Opt. Precision Eng, 2017, 25(10): 2682-07. (in Chinese) doi:  10.3788/OPE.20172510.2682
    [10] Wang Xiaokun. Testing of a large rectangular mirror based on sub-aperture stitching method [J]. Chinese Optics Letter, 2015, 13(13): 40-43.
    [11] 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
    [12] Yan Lisong, Wang Xiaokun, Zheng Ligong, et al. Large-diameter reflector mirror stitching accuracy analysis based on selftest [J]. Infrared and Laser Engineering, 2014, 43(6): 1920-1924. (in Chinese) doi:  10.3969/j.issn.1007-2276.2014.06.040
    [13] Wang Xiaokun, Xue Donglin, Zhang Xuejun. Fabrication and testing of large aspheric system based on common reference [J]. Opt Precision Eng, 2018, 26(04): 743-748. (in Chinese) doi:  10.3788/OPE.20182604.0743
    [14] Mourad Idir, Konstantine Kaznatcheev, Guillaume Dovillaire, et al. A 2 D high accuracy slope measuring system based on a stitching shack hartmann optical head [J]. Optics Express, 2014, 22(3): 2770-2781. doi:  10.1364/OE.22.002770
    [15] Mallik P C V, Chunyu Zhao, Burge J H. Measurement of a 2-meter flat using a pentaprism scanning system [J]. Optical Engineering, 2007, 46(2): 0236021.
    [16] Wei H S, Yan F, Chen X D. Large aperture space optical system testing based on scanning hartmann [J]. Appl Opt, 2017, 56(8): 2078-2083. doi:  10.1364/AO.56.002078
    [17] Wang Xiaokun. Test convex aspheres by subaperture stitching interferometry [J]. Infrared and Laser Engineering, 2013, 42(3): 716-722. (in Chinese) doi:  10.3969/j.issn.1007-2276.2013.03.044
    [18] Wu Jing, Wang Jianli, Lin Xudong. Wave-front test by sub-aperture stitching technique based on Shack-Hartmann wave-front sensor [J]. Opto-Electronic Engineering, 2011, 38(1): 39-43. (in Chinese)
    [19] Burada D R, Pant K K , Bichra M, et al. Experimental investigations on characterization of freeform wavefront using Shack-Hartmann sensor [J]. Optical Engineering, 2017, 56(8): 084107.
    [20] Hongyan Xu, Hao Xian, Yudong Zhang. Comparison of two stitching algorithms for annular subaperture Hartmann Shack testing method[C]//SPIE, 2010, 7849: 78491Q.
    [21] Wang Ruoqiu, Zhang Zhiyu, Xue Donglin, et al. Large-diameter high-efficiency diffractive Fresnel membrane elements for space telescope [J]. Infrared and Laser Engineering, 2017, 46(9): 0920001. (in Chinese) doi:  10.3788/IRLA201746.0920001
  • [1] 高博, 范斌, 王佳, 吴湘, 辛强.  光学元件磁流变加工不确定度误差工艺方法 . 红外与激光工程, 2024, 53(3): 20230595-1-20230595-10. doi: 10.3788/IRLA20230595
    [2] 张誉馨, 黎发志, 闫力松, 严冬, 王超凡, 张宇鹤, 张斌智, 张鑫, 蔡铭宣, 何俊超, 王孝坤.  结合CGH与辅助透镜的长焦距非球面反射镜检测(特邀) . 红外与激光工程, 2022, 51(9): 20220384-1-20220384-6. doi: 10.3788/IRLA20220384
    [3] 苏航, 王孝坤, 程强, 李凌众, 王晶, 李雯研, 吴琼, 唐瓦, 罗霄, 张学军.  子孔径拼接和计算全息混合补偿检测大口径凸非球面(特邀) . 红外与激光工程, 2022, 51(9): 20220576-1-20220576-9. doi: 10.3788/IRLA20220576
    [4] 沈正祥, 王旭, 余俊.  基于离散化路径与拼接的高精度光学检测方法 . 红外与激光工程, 2022, 51(11): 20210144-1-20210144-10. doi: 10.3788/IRLA20210144
    [5] 王孝坤, 戚二辉, 胡海翔, 苏航, 李凌众, 王晶, 罗霄, 张学军.  超大口径平面反射镜的光学检测(特邀) . 红外与激光工程, 2022, 51(1): 20210953-1-20210953-7. doi: 10.3788/IRLA20210953
    [6] 李小明, 王隆铭, 朱国帅.  激光通信一体化SiC/Al摆镜支撑参数优化 . 红外与激光工程, 2021, 50(11): 20210143-1-20210143-8. doi: 10.3788/IRLA20210143
    [7] 闫力松, 张斌智, 王孝坤, 黎发志.  平面镜子孔径加权拼接检测算法(特邀) . 红外与激光工程, 2021, 50(10): 20210520-1-20210520-6. doi: 10.3788/IRLA20210520
    [8] 相超, 王道档, 窦进超, 孔明, 刘璐, 许新科.  光学偏折子孔径拼接面形检测技术 . 红外与激光工程, 2021, 50(11): 20210105-1-20210105-7. doi: 10.3788/IRLA20210105
    [9] 王辰忠, 胡中文, 陈忆, 许明明, 陈力斯.  空间引力波望远镜主反射镜系统的结构设计优化 . 红外与激光工程, 2020, 49(7): 20190469-1-20190469-10. doi: 10.3788/IRLA20190469
    [10] 刘娜, 沈正祥, 马彬, 魏振博, 徐旭东, 王占山.  圆锥近似Wolter-I型X射线望远镜用柱面反射镜面形误差检测方法 . 红外与激光工程, 2018, 47(4): 417001-0417001(7). doi: 10.3788/IRLA201847.0417001
    [11] 谢军, 何锋赟, 王晶, 高阁, 赵天骄, 刘震宇.  经纬仪主镜轴向支撑结构仿真与优化 . 红外与激光工程, 2016, 45(S1): 132-137. doi: 10.3788/IRLA201645.S118001
    [12] 闫公敬, 张宪忠.  子孔径拼接检测凸球面技术研究 . 红外与激光工程, 2016, 45(5): 517002-0517002(4). doi: 10.3788/IRLA201645.0517002
    [13] 戚二辉, 罗霄, 李明, 郑立功, 张学军.  五棱镜扫描技术检测大口径平面镜的误差分析 . 红外与激光工程, 2015, 44(2): 639-646.
    [14] 闫力松, 王孝坤, 郑立功, 王润强, 张学军.  大口径反射镜子孔径拼接自检验精度分析 . 红外与激光工程, 2014, 43(6): 1920-1924.
    [15] 邓万涛, 汪凯巍, 白剑, 张金春.  高精度子孔径拼接中参考面误差的去除方法 . 红外与激光工程, 2014, 43(4): 1194-1199.
    [16] 闫公敬, 张宪忠, 李柱.  子孔径拼接检测光学平面反射镜技术 . 红外与激光工程, 2014, 43(7): 2180-2184.
    [17] 郭玲玲, 张星祥, 张立国, 任建岳.  离轴非球面镜检测中坐标系偏差对参数拟合的影响 . 红外与激光工程, 2013, 42(10): 2800-2804.
    [18] 闫力松, 王孝坤, 罗霄, 郑立功, 张学军.  采用三角剖分算法的子孔径拼接检测技术 . 红外与激光工程, 2013, 42(7): 1793-1797.
    [19] 王孝坤.  大口径凸非球面反射镜子孔径拼接检测 . 红外与激光工程, 2013, 42(3): 716-722.
    [20] 王宁, 贾辛, 邢廷文.  子孔径拼接检测浅度非球面 . 红外与激光工程, 2013, 42(9): 2525-2530.
  • 加载中
图(12) / 表(3)
计量
  • 文章访问数:  237
  • HTML全文浏览量:  94
  • PDF下载量:  62
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-07-30
  • 修回日期:  2021-08-27
  • 刊出日期:  2021-10-20

夏克哈特曼扫描拼接检测平面镜(特邀)

doi: 10.3788/IRLA20210527
    作者简介:

    王晶,女,博士生,主要从事光学设计与测试技术方面的研究

基金项目:  吉林省科技发展计划(20200401065GX);中国科学院青年创新促进会(2019221);国家自然科学基金(61975201);国家自然科学基金面上项目(6207031149)
  • 中图分类号: O439;O436.1

摘要: 随着先进光学系统设计与制造的发展,大口径光学系统得到了广泛的应用。然而,大口径平面镜高精度面形的检测手段不足,限制了大口径平面镜的制造与应用。为实现大口径平面反射镜的高精度面形检测,提出一种夏克哈特曼扫描拼接检测平面镜面形的方法。对扫描拼接原理、波前重构算法进行了研究,建立了微透镜阵列成像的数学模型,验证了夏克哈特曼扫描拼接检测原理的可行性。针对一口径为150 mm的平面镜进行了扫描拼接检测实验,拼接得到的全口径面形为0.019λ RMS(λ=635 nm);与干涉检测结果对比,检测精度为0.008λ RMS,结果表明该方法能够实现大口径平面反射镜的高精度检测。

English Abstract

    • 光学系统中,口径越大,则系统的角分辨能力和能量收集能力越强,意味着观测范围越大,观测物体越清晰,因此大口径光学系统是先进光学系统的必然发展趋势[1]。大口径光学望远镜系统检测时需要大口径平面反射镜作为标准反射镜,平面镜的精度直接决定了望远镜系统的精度。同时,大口径平面镜在很多大口径光学系统中担任着折镜的重要角色,例如LAMOST中反射施密特改正镜MA由24块对角线长为1 100 mm的六边形超薄平面镜组成[2];VLT中M3镜是长短轴分别为1 250 mm和800 mm的椭圆形大口径平面镜[3]。此外,大口径平面镜在大型高精度衍射光栅刻划系统中作为高精度平面导轨使用[4]。然而,大口径平面镜的高精度面形检测手段不足,一直是光学检测领域的难题。

      大口径平面镜的检测方法主要有:大口径干涉仪法[5]、Skip-Flat法[6]、Ritchey-Common法[7]、子孔径拼接法[8-10]、五棱镜扫描法[11]。大口径干涉仪检测方法简单,检测精度最高,但是需要大口径的平面标准镜头,制作难度大,周期长,成本高;Skip-Flat法优点在于通过斜入射减小了干涉仪和参考平面镜的口径,但是,随着入射角的增加,灵敏度会降低;Ritchey-Common法利用口径大于平面镜的球面反射镜替换平面参考镜,球面波斜入射平面镜,由于入射角分布的差异会影响波像差与面形差的转换关系,引入误差,随着入射角的增大,检测灵敏度减小;子孔径拼接法只需要小口径的干涉仪,在干涉仪口径与平面镜口径相差很大时,会引入相当大的拼接误差;五棱镜扫描法采用高精度自准直仪,五棱镜与大行程直线导轨相结合,扫描被测面上若干条轮廓线的法向角度变化,通过积分拼接出全口径面形,受直线导轨的误差影响大,检测精度较低。

      夏克哈特曼扫描拼接检测方法是基于夏克哈特曼传感器探测头(HASO)扫描的斜率检测方法,借鉴子孔径干涉检测[12-14]的拼接思想,将夏克哈特曼传感器采集到的各个子孔径信息拼接出整个口径的斜率信息,利用基于Zernike多项式的模式重构方法[15]恢复检测面形。检测方法简单,采样密度大,可实现高精度的面形检测。文中针对口径150 mm平面镜的检测进行了设计分析和测试,验证了夏克哈特曼扫描拼接检测大口径平面镜的适用性和高效性。

    • 检测光路示意图如图1所示,其主要装置为HASO(夏克哈特曼波前传感器检测头)。激光发出的点光源经过准直镜后产生的平行光被分束镜反射到待测平面镜,从平面镜反射回来的光束透过分束镜被微透镜阵列会聚到位于其焦面的CCD探测器上形成光斑阵列图。

      图  1  夏克哈特曼扫描拼接检测平面镜原理示意图

      Figure 1.  Schematic diagram of Shack-Hartmann scanning and splicing plane mirror

    • 通过质心算法[15]提取出CCD上的离散光斑位置,比较实际光斑与理想光斑的位置偏离,求出检测区域在x方向与y方向的波前斜率SxSy,利用基于Zernike多项式的模式法波前重构进行波前恢复。

      夏克—哈特曼波前传感器采样测量得到的是离散点数据,依据波前斜率与波像差关系,建立如下所示的超定方程组[15]

      $$ \left\{ \begin{matrix} {{S}_{x}}=\displaystyle\sum{_{k=1}^{n}{{C}_{k}}\dfrac{\partial {{Z}_{k}}\left( x,y \right)}{\partial x}} \\ {{S}_{y}}=\displaystyle\sum{_{k=1}^{n}{{C}_{k}}\dfrac{\partial {{Z}_{k}}\left( x,y \right)}{\partial y}} \\ \end{matrix} \right. $$ (1)
      $$ \left[ {\begin{array}{*{20}{c}} {\dfrac{{\partial {Z_{11}}\left( {{{x}},{{y}}} \right)}}{{\partial x}}}&{\dfrac{{\partial {Z_{12}}\left( {{{x}},{{y}}} \right)}}{{\partial x}}}& \cdots &{\dfrac{{\partial {Z_{1n}}\left( {{{x}},{{y}}} \right)}}{{\partial x}}}\\ \vdots &\vdots &\vdots &\vdots \\ {\dfrac{{\partial {Z_{N1}}\left( {{{x}},{{y}}} \right)}}{{\partial x}}}&{\dfrac{{\partial {Z_{N2}}\left( {{{x}},{{y}}} \right)}}{{\partial x}}}& \cdots &{\dfrac{{\partial {Z_{Nn}}\left( {{{x}},{{y}}} \right)}}{{\partial x}}}\\ {\dfrac{{\partial {Z_{11}}\left( {{{x}},{{y}}} \right)}}{{\partial y}}}&{\dfrac{{\partial {Z_{12}}\left( {{{x}},{{y}}} \right)}}{{\partial y}}}& \cdots &{\dfrac{{\partial {Z_{1n}}\left( {{{x}},{{y}}} \right)}}{{\partial y}}}\\ \vdots &\vdots &\vdots &\vdots \\ {\dfrac{{\partial {Z_{N1}}\left( {{{x}},{{y}}} \right)}}{{\partial y}}}&{\dfrac{{\partial {Z_{N2}}\left( {{{x}},{{y}}} \right)}}{{\partial y}}}& \cdots &{\dfrac{{\partial {Z_{Nn}}\left( {{{x}},{{y}}} \right)}}{{\partial y}}} \end{array}} \right] \cdot \left[ {\begin{array}{*{20}{c}} {{C_1}}\\ {{C_2}}\\ \vdots \\ {{C_n}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{S_{x1}}}\\ \vdots \\ {{S_{xN}}}\\ {{S_{y1}}}\\ \vdots \\ {{S_{yN}}} \end{array}} \right] $$ (2)

      式中:k为Zernike多项式的项数;C为Zernike多项式的系数;n为离散数据的个数。

      利用极大范数下最小二乘法[15]求解模式系数,代入Zernike多项式中可以重构出波前。

    • 图2所示,选择待测平面镜的中心子孔径作为基准子孔径,子孔径间重叠区域的波前斜率满足:

      图  2  子孔径中心位置规划示意图

      Figure 2.  Schematic diagram of sub-aperture center location planning

      $$ {W_x}\left( {x,y} \right) = W{1_x}\left( {x,y} \right) + {T_{x1}} + 2{D_1}x = W{2_x}\left( {x,y} \right) + {T_{x2}} + 2{D_2}x $$ (3)
      $$ {W_y}\left( {x,y} \right) = W{1_y}\left( {x,y} \right) + {T_{y1}} + 2{D_1}y = W{2_y}\left( {x,y} \right) + {T_{y2}} + 2{D_2}y $$ (4)

      式中:Txx方向倾斜修正量;Tyy方向倾斜修正量;D为离焦修正量。通过最小二乘算法拟合,使倾斜,离焦误差最小,进行误差解耦[16-20]

      由于夏克哈特曼传感器在扫描检测过程中产生机械移动定位误差,利用微孔径光阑定位装置,在测量前,采用迭代对准方法对偏移子孔径进行迭代对准,提高机械定位精度。

    • 夏克哈特曼波前传感器中的核心元件为微透镜阵列,利用傅里叶光学的方法可以描述微透镜阵列成像的物理过程。其物理过程如图3所示,已知紧靠微透镜的入射波前U0(x,y),经过微透镜调制后的波前U(x,y)传播到焦面CCD上,CCD对焦面上的波面Uf(x,y)的强度分布进行探测,CCD探测到的光斑和强度即为入射光经过微透镜阵列在焦面上的菲涅耳衍射图样及强度分布。

      图  3  单个微透镜成像的傅里叶光学物理过程

      Figure 3.  Fourier optical physics in the imaging of a single microlens

      单个微透镜的复振幅透过率Ts为:

      $$ {{{T}}_{{S}}}\left( {{{x}},{{y}}} \right) = \exp \left( {{{jkn}}\Delta } \right)\exp \left[ { - {{j}}\frac{k}{{2f}}\left( {{x^2} + {y^2}} \right)} \right] $$ (5)

      式中:j为虚数单位;n为微透镜的折射率,为微透镜的中心厚度;f为微透镜的焦距;${{k}} = \dfrac{{2{\rm{\pi }}}}{\lambda }$为波数。

      第一项位相因子仅表示透镜对于入射光波的常量位相延迟,并不影响位相的空间相对分布。即它不会改变光波波面的形状,故可以略去,不予考虑;第二项是透镜的位相调制因子。

      单个微透镜孔径为正方形,设其边长为a,则单个微透镜的孔径透过率函数Tk为:

      $$ {{T_k}} = {\rm{rect}}\left( {\frac{x}{a}} \right){\rm{rect}}\left( {\frac{y}{a}} \right) $$ (6)

      紧靠微透镜后的波前U(x,y)为:

      $$ {{U}}\left( {{{x}},{{y}}} \right) = {u_0}\left( {x,y} \right) \cdot T_s\left( {x,y} \right) \cdot T_k\left( {x,y} \right) $$ (5)

      根据菲涅耳衍射公式,透镜后焦面上的光场分布正比于物体的傅里叶变换[21],则焦面上的强度分布为:

      $$ {{I}} \propto {\left| {F\left( {U\left( {x,y} \right)} \right)} \right|^2} $$ (7)

      利用MATLAB仿真上述过程,如图4所示,为被单个微透镜分割的波前区域会聚在焦面CCD上的光斑成像结果。

      图  4  单个微透镜成像的仿真结果

      Figure 4.  Simulation results of single microlens imaging

      整个微透镜阵列的成像相较于单个微透镜成像,在傅里叶光学上的区别在于孔径透过率函数,微透镜阵列的孔径透过率函数主要由三部分组成:(1)单个微透镜的孔径矩形函数rect(x/a) rect(y/a);(2)微透镜阵列中同形孔径构成的列阵,可视为梳状函数comb(x/a) comb(y/a);(3)将微透镜阵列整体整体示为边长为L的正方形孔径,rect(x/L) rect(y/L),a为单个微透镜孔径,同时也是相邻微透镜中心的间隔。

      则,整个微透镜阵列的孔径透过率函数T(x,y)为:

      $$\begin{split} {{T}}\left( {{{x}},{{y}}} \right) =& \left[ {{\rm{rect}}\left( {\frac{x}{a}} \right){\rm{rect}}\left( {\frac{y}{a}} \right){\rm{comb}}\left( {\frac{x}{a}} \right){\rm{comb}}\left( {\frac{y}{a}} \right)} \right] \cdot\\ &{\rm{rect}}\left( {\frac{x}{L}} \right){\rm{rect}}\left( {\frac{y}{L}} \right) \end{split}$$ (8)

      依据微透镜阵列傅里叶光学成像原理,如图5所示,为平面波经过微透镜阵列会聚在焦面CCD上的光斑成像结果。

      图  5  微透镜阵列成像的仿真结果

      Figure 5.  Simulation results of microlens array imaging

      全口径仿真中各项参数如表1所示;依据子孔径规划,仿真检测37个子孔径,将37个子孔径的光斑阵列图通过质心提取得到波前斜率矩阵,利用Zernike多项式的模式法波前重构和子孔径拼接得到全口径波前。在全口径上添加如图6(a)所示的波像差,RMS=0.182λ,PV=1.365λ,仿真得到的全口径波前与输入像差的误差如图6(c)所示,RMS=1.531e-06λ,PV=1.122e-05λ

      表 1  仿真参数

      Table 1.  Simulation parameters

      ParameterValue
      Microlens array diameter/mm 14.61
      Mirror radius /mm 25
      Number of subaperture 37
      Wave length/nm 635

      图  6  全口径仿真结果。(a) 添加的波像差; (b) 拼接仿真结果; (c) 结果误差

      Figure 6.  Full caliber simulation results. (a) Added wave aberrations; (b) Stitching simulation results; (c) Result error

    • 图7所示为由HASO,五维调节运动平台(xyz,偏转角度A,旋转角度C,其量程与精度如表2所示)和平面镜组成的夏克哈特曼扫描拼接检测系统。

      图  7  夏克哈特曼扫描拼接检测系统

      Figure 7.  Shack-Hartmann scanning and splicing detection system

      表 2  移动平台量程与精度

      Table 2.  Moving platform range and accuracy

      ParameterValueAccuracy
      x /mm1 0000.01
      y /mm5000.01
      z /mm8000.02
      A90 °4″
      C360 °10″
    • 检测过程中需要使出射平行光中心垂直入射子孔径检测区域中心,即满足出射平行光的中心线,载物平台的法线和待测平面镜中心法线“三线合一”的对准调整。

      首先,利用电子水平仪调整平台水平;放置平面镜,找到其上下左右四个边缘位置,如图8所示,并计算平面镜中心位置;安装定心装置(孔径光阑),选择alignment模式,调整偏转角度A,完成精调对准,如图9所示。

      图  8  平面镜边缘位置。(a)上边缘;(b)下边缘;(c) 左边缘;(d) 右边缘

      Figure 8.  Plane mirror edge position. (a) Uupper edge;(b) Lower edge;(c) Left edge; (d) Right edge

      由于HASO内部各元件存在客观的装调误差和制造误差,将检测标准平面镜的光斑位置作为标定结果,利用实际光斑位置对应减去标定光斑位置,可以去除掉系统误差。

      图  9  精调对准过程。(a)未对准;(b)对准

      Figure 9.  Fine alignment process. (a) Misalignment;(b) Alignment

    • 实验中各项参数如表3所示,待检平面镜口径150 mm,子孔径检测口径为14.6 mm,17×17个子孔径可以实现覆盖全口径的检测。利用840D闭环控制系统对各子孔径进行了自动检测。

      表 3  实验参数

      Table 3.  Experiment parameters

      ParameterValue
      Detection beam aperture /mm14.6
      Mirror diameter /mm150
      Sampling density17×17
      Scanning step /mm9.1
      Wave length/nm635

      通过质心提取,求离散斜率矩阵,子孔径拼接,得到了全口径x方向斜率,y方向斜率,如图10所示。

      图  10  全口径斜率检测结果。(a) x方向斜率; (b) y方向斜率

      Figure 10.  Full caliber slope test results. (a) Slope in the x direction; (b) Slope in the y direction

      利用斜率模式法波前重构计算后得到的全口径面形RMS=0.019λ,PV=0.187λ,如图11(a)所示。与干涉仪检测结果图11(b)对比,残差分布如图12所示,PV=0.016λ,RMS=0.008λ。实验结果表明夏克哈特曼扫描拼接检测能够满足平面镜表面面形的高精度检测,验证了夏克哈特曼扫描拼接检测大口径平面镜的可行性。

      图  11  (a)全口径面形;(b)干涉仪检测结果

      Figure 11.  (a) Full aperture surface; (b) Interferometer test results

      图  12  残差分布

      Figure 12.  Residual distribution

    • 为了实现大口径平面镜面形的高精度检测,文中提出了夏克哈特曼扫描拼接检测平面反射镜的方法。阐述了该检测方法原理,利用傅里叶光学建立了微透镜阵列成像的数学模型,规划了子孔径,完成了波前斜率重构算法,进行了原理仿真验证。并对一口径为150 mm的平面镜进行了夏克哈特曼扫描拼接检测的验证实验,检测精度为0.008λRMS,实验结果表明该检测方法在检测平面反射镜时,具有较好的可行性和实用性。另外,若将子孔径倾斜检测,可实现球面、非球面、自由曲面等更复杂曲面面形的检测,从而为先进光学系统的制造提供了依据与保障。

参考文献 (21)

目录

    /

    返回文章
    返回