留言板

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

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

消像散的自由曲面棱镜光谱仪光学系统设计

张佳伦 郑玉权 蔺超 纪振华

张佳伦, 郑玉权, 蔺超, 纪振华. 消像散的自由曲面棱镜光谱仪光学系统设计[J]. 中国光学. doi: 10.37188/CO.2019-0049
引用本文: 张佳伦, 郑玉权, 蔺超, 纪振华. 消像散的自由曲面棱镜光谱仪光学系统设计[J]. 中国光学. doi: 10.37188/CO.2019-0049
ZHANG Jia-lun, ZHENG Yu-qun, LIN Chao, JI Zhen-hua. Design of a Freeform Curved Prism Imaging Spectrometer Based on an Anastigmatism[J]. Chinese Optics. doi: 10.37188/CO.2019-0049
Citation: ZHANG Jia-lun, ZHENG Yu-qun, LIN Chao, JI Zhen-hua. Design of a Freeform Curved Prism Imaging Spectrometer Based on an Anastigmatism[J]. Chinese Optics. doi: 10.37188/CO.2019-0049

消像散的自由曲面棱镜光谱仪光学系统设计

doi: 10.37188/CO.2019-0049
基金项目: 国家重点研发计划(2016YFB0500300)
详细信息
    作者简介:

    张佳伦(1994—),男,吉林长春人,硕士研究生,2016年于长春理工大学获得学士学位,主要从事高光谱成像技术、光学系统设计方面的研究。E-mail:jlzhang0807@163.com

    郑玉权(1972—),男,内蒙古通辽人,博士,博士生导师,研究员,主要从事航空航天高光谱成像技术、光学系统设计、光谱辐射测量与定标等方面的研究。E-mail:zhengyq@sklao.ac.cn

    蔺超(1984—),男,内蒙古呼伦贝尔人,硕士,副研究员,主要从事光谱仪器结构设计及光谱定标方面的研究。E-mail:linchao@ciomp.ac.cn

    纪振华(1986—),男,黑龙江双鸭山人,硕士,助理研究员,主要从事空间遥感器设计、光谱辐射定标等方面的研究。E-mail:jizhenhua001@126.com

  • 中图分类号: TP394.1;TH691.9

Design of a Freeform Curved Prism Imaging Spectrometer Based on an Anastigmatism

Funds: National Key Research and Development Program of China (22016YFB0500300)
More Information
  • 摘要: 本文根据像差理论,设计一套计算次镜外反射的Offner型自由曲面棱镜光谱仪初始结构的算法。通过光线追迹获得光线在次镜外反射Offner型光谱仪各光学表面传播的公式,该公式可以确定光学元件的结构参数。应用轴外细光束像散理论分析系统产生的像散,并设定合理阈值作为结构算法的判断依据。在Matlab中迭代优化出符合设计要求的初始结构,获得的初始结构在Zemax中优化。为验证算法效果,本文设计光谱为380−780 nm,数值孔径为0.15,光谱分辨率为6 nm的自由曲面棱镜光谱仪的初始结构。在Zemax中完成优化后系统可达到设计指标且谱线弯曲和色畸变均优于0.1 pixels。设计结果表明采用本文的算法可以快速计算出符合要求的初始结构,大大简化了后续优化的复杂程度。
  • 图  1  次镜外反射型自由曲面棱镜光谱仪结构图

    Figure  1.  Structural diagram of the freeform curved prism imaging spectrometer

    图  2  次镜前自由曲面棱镜光谱仪光路示意图

    Figure  2.  Schematic diagram of the freeform curved prism spectrometer light path before secondary mirror

    图  3  次镜后自由曲面棱镜光谱仪光路示意图

    Figure  3.  Schematic diagram of the freeform curved prism spectrometer light path after secondary mirror

    图  4  轴外细光束像散示意图

    Figure  4.  Schematic diagram of off-axis rays astigmatism

    图  5  本文算法流程图

    Figure  5.  Flow chart of proposed algorithm

    图  6  初始自由曲面棱镜光谱仪结构示意图

    Figure  6.  Schematic diagram of the initial freeform curved prism spectrometer

    图  7  优化后的自由曲面棱镜光谱仪结构示意图

    Figure  7.  Schematic diagram of the optimized freeform curved prism spectrometer

    图  8  光学系统调制传递函数

    Figure  8.  Schematic diagram of MTF

    表  1  设计指标

    Table  1.   Design specifications

    参数设计值
    光谱范围380−780 nm
    数值孔径0.15
    光谱分辨率6 nm
    像元尺寸13 μm × 13 μm
    色散元件材料ZF6
    色散面大小1 mm
    下载: 导出CSV

    表  2  初始结构数据(单位:毫米)

    Table  2.   Initial structural data (Unit: mm)

    表面半径厚度材料偏心量
    1——249.500 000
    2−248.510 000−104.255 000MIRROR−45.990 000
    3−138.305 000−15.990 000ZF6−36.980 000
    4−122.315 000−5.990 000−43.980 000
    5−122.275 0005.990 000MIRROR−45.990 000
    6−122.315 00015.990 000ZF6−43.980 000
    7−138.305 000104.255 000−36.980 000
    8−248.510 000——MIRROR−45.990 000
    9————
    下载: 导出CSV

    表  3  优化后结构数据(单位:毫米)

    Table  3.   Optimized structural data (Unit: mm)

    表面半径厚度材料偏心量
    1——240.045 897
    2−234.113 632−99.964 613MIRROR−58.302 387
    3−138.439 720−18.303 653ZF6−42.950 734
    4−120.219 597−7.986 594−49.982 314
    5−105.485 6317.986 594MIRROR−60.171 600
    6−120.219 59718.303 653ZF6−49.982 314
    7−138.439 72089.943 651−42.950 734
    8−223.969 005——MIRROR−61.713 353
    9————
    下载: 导出CSV
  • [1] 郑玉权. 小型Offner光谱成像系统的设计[J]. 光学精密工程,2005,13(6):650-657. doi:  10.3321/j.issn:1004-924X.2005.06.005

    ZHENG Y Q. Design of compact Offner spectral imaging system[J]. Optics and Precision Engineering, 2005, 13(6): 650-657. (in Chinese) doi:  10.3321/j.issn:1004-924X.2005.06.005
    [2] SANG B, SCHUBERT J, KAISER S, et al. The EnMAP hyperspectral imaging spectrometer: instrument concept, calibration, and technologies[J]. Proceedings of SPIE, 2008, 7086: 708605. doi:  10.1117/12.794870
    [3] FÉRY C. A prism with curved faces, for spectrograph or spectroscope[J]. The Astrophysical Journal, 1911, 34(1): 79.
    [4] MORRISSEY P F. Third-order aberrations of a prism with spherically curved surfaces[J]. Applied Optics, 1994, 33(13): 2539-2543. doi:  10.1364/AO.33.002539
    [5] 刘力. Fery棱镜特性及应用研究[D]. 苏州: 苏州大学, 2013.

    LIU L. Research of the properties and application of Fery prism[D]. Suzhou: Soochow University, 2013. (in Chinese)
    [6] 刘立银. 曲面棱镜Offner高光谱成像系统的设计研究[D]. 南京: 南京理工大学, 2016.

    LIU L Y. Design and research of curved prism Offner hyperspectral imaging system[D]. Nanjing: Nanjing University of Science and Technology, 2016. (in Chinese)
    [7] 石顺祥, 王学恩, 刘劲松. 物理光学与应用光学[M]. 西安: 西安电子科技大学出版社, 1999.

    SHI SH X, WANG X E, LIU J S. Physical Optics and Applied Optics[M]. Xi’an: Xi’an University of Electronic Science and Technology Press, 1999. (in Chinese)
    [8] 潘君骅. 光学非球面的设计、加工与检验[M]. 苏州: 苏州大学出版社, 2004.

    PAN J H. The Design, Manufacture and Test of the Aspherical Optical Surfaces[M]. Suzhou: Soochow University Press, 2004. (in Chinese)
    [9] LOBB D R. Imaging spectrometers using concentric optics[J]. Proceedings of SPIE, 1997, 3118: 339-347. doi:  10.1117/12.283838
    [10] 张云翠, 刘龙, 曹冠英, 等. Fery棱镜光谱仪设计[J]. 红外与激光工程,2009,38(2):287-289. doi:  10.3969/j.issn.1007-2276.2009.02.021

    ZHANG Y C, LIU L, CAO G Y, et al. Design of spectrograph with novel Fery prism[J]. Infrared and Laser Engineering, 2009, 38(2): 287-289. (in Chinese) doi:  10.3969/j.issn.1007-2276.2009.02.021
    [11] 聂云峰, 相里斌, 周锦松, 等. 基于曲面棱镜的宽视场推帚式高光谱成像仪设计[J]. 光谱学与光谱分析,2012,32(6):1708-1711. doi:  10.3964/j.issn.1000-0593(2012)06-1708-04

    NIE Y F, XIANG L B, ZHOU J S, et al. A wide-field push-broom hyperspectral imager based on curved prism[J]. Spectroscopy and Spectral Analysis, 2012, 32(6): 1708-1711. (in Chinese) doi:  10.3964/j.issn.1000-0593(2012)06-1708-04
    [12] 刘光宏, 吴刚, 凌青, 等. Offner成像光谱仪建模及像差分析[J]. 光学学报,2011,31(3):0322001. doi:  10.3788/AOS201131.0322001

    LIU H G, WU G, LING Q, et al. Modelling of Offner imaging spectrometers and aberrations analysis[J]. Acta Optica Sinica, 2011, 31(3): 0322001. (in Chinese) doi:  10.3788/AOS201131.0322001
    [13] FENG L, ZHOU J S, WANG D. A miniature design of near infrared spectrometer with wide field of view[J]. Proceedings of SPIE, 2014, 9298: 92980M.
    [14] 赵渊明. 消像散自由曲面离轴显示系统设计方法研究[D]. 长春: 长春理工大学, 2018.

    ZHAO Y M. Studies on the design of off-axis display systems with anastigmatic free-form surface[D]. Changchun: Changchun University of Science and Technology, 2018. (in Chinese)
    [15] FENG L, ZHOU J S, JING J J, et al. A design of miniature imaging spectrometer with Fery prism[J]. Proceedings of SPIE, 2016, 10000: 100001P.
    [16] DONG W, NIE Y F, ZHOU J S, et al. Design of a concise Féry-prism hyperspectral imaging system based on multi-configuration[J]. Proceedings of SPIE, 2013, 8910: 89100B.
  • [1] 魏雨, 蒋世磊, 孙国斌, 张兴星, 王玉宁.  固态面阵激光雷达接收光学系统设计 . 中国光学, doi: 10.3788/CO.2019-0166
    [2] 吕博, 冯睿, 寇伟, 刘伟奇.  折反射式空间相机光学系统设计与杂散光抑制 . 中国光学, doi: 10.37188/CO.2019-0036
    [3] 黄颂超, 冯云鹏, 程灏波.  非对称轻小型头盔显示器光学系统设计 . 中国光学, doi: 10.37188/CO.2019-0193
    [4] 王琪, 梁静秋, 梁中翥, 吕金光, 王维彪, 秦余欣, 王洪亮.  分孔径红外偏振成像仪光学系统设计 . 中国光学, doi: 10.3788/CO.20181101.0092
    [5] 孟竹, 谭鑫, 巴音贺希格, 王玮, 刘兆武, 赵旭龙.  Seya-Namioka单色仪中光栅曲率半径误差的影响及补偿 . 中国光学, doi: 10.3788/CO.20181104.0623
    [6] 史光辉.  用高斯光学和三级像差理论求变焦距物镜的初始解 . 中国光学, doi: 10.3788/CO.20181106.1047
    [7] 于潇, 姚园, 徐正平.  采用APD阵列的共口径激光成像光学系统设计 . 中国光学, doi: 10.3788/CO.20160903.0349
    [8] 薛金来, 巩岩, 李佃蒙.  N.A.0.75平场复消色差显微物镜光学设计 . 中国光学, doi: 10.3788/CO.20150806.0957
    [9] 杨飞, 金光, 谢金华, 邱振戈, 曲宏松, 贺小军.  基于在轨成像物理机理的立体测绘相机建模与仿真 . 中国光学, doi: 10.3788/CO.20150806.0971
    [10] 曲贺盟, 张新.  高速切换紧凑型双视场无热化红外光学系统设计 . 中国光学, doi: 10.3788/CO.20140704.0622
    [11] 薛珮瑶, 吴耀, 冯茜, 李川.  大视场四象限探测光学系统设计 . 中国光学, doi: 10.3788/CO.20140703.0462
    [12] 王晶, 崔恩坤.  大视场曲面复眼光学系统设计 . 中国光学, doi: 10.3788/CO.20140706.0969
    [13] 李维善, 陈琛, 于国辉, 李臣友, 刘宵婵, 张禹, 王陆军, 刘影.  单机数字立体双通道放映镜头的光学设计 . 中国光学, doi: 10.3788/CO.20130604.0563
    [14] 董科研, 王健, 孙强, 王洪亮, 李全熙, 姜会林.  机载双视场中波红外光学系统优化设计 . 中国光学, doi: 10.3788/CO.20120506.0596
    [15] 李蕊, 刘永基, 王肇圻.  基于个体眼光学结构的波前眼镜设计 . 中国光学, doi: 10.3788/CO.20120505.0512
    [16] 安岩, 孙强, 刘英, 李淳.  交叉型消像散Czerny-Turner结构光谱仪设计 . 中国光学, doi: 10.3788/CO.20120505.0470
    [17] 沈炎龙, 黄珂, 马连英, 郑国鑫, 王大辉.  高功率准分子激光主振荡功率放大系统光学元件的稳定性 . 中国光学,
    [18] 莫仁芸, 朱万彬, 钟俊, 陈璇.  激光位移传感器的光学系统设计 . 中国光学,
    [19] 李萍, 卢振武, 夏利东, 张红鑫, 刘英.  反射式内掩日冕仪的光学设计与杂散光分析 . 中国光学,
    [20] 李宁, 唐勇, 李玉瑶, 耿似玉.  工作距离可调的激光显微操纵器光学系统设计 . 中国光学,
  • 加载中
图(8) / 表ll (3)
计量
  • 文章访问数:  244
  • HTML全文浏览量:  134
  • PDF下载量:  7
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-02-25
  • 修回日期:  2019-04-06
  • 网络出版日期:  2020-06-29

消像散的自由曲面棱镜光谱仪光学系统设计

doi: 10.37188/CO.2019-0049
    基金项目:  国家重点研发计划(2016YFB0500300)
    作者简介:

    张佳伦(1994—),男,吉林长春人,硕士研究生,2016年于长春理工大学获得学士学位,主要从事高光谱成像技术、光学系统设计方面的研究。E-mail:jlzhang0807@163.com

    郑玉权(1972—),男,内蒙古通辽人,博士,博士生导师,研究员,主要从事航空航天高光谱成像技术、光学系统设计、光谱辐射测量与定标等方面的研究。E-mail:zhengyq@sklao.ac.cn

    蔺超(1984—),男,内蒙古呼伦贝尔人,硕士,副研究员,主要从事光谱仪器结构设计及光谱定标方面的研究。E-mail:linchao@ciomp.ac.cn

    纪振华(1986—),男,黑龙江双鸭山人,硕士,助理研究员,主要从事空间遥感器设计、光谱辐射定标等方面的研究。E-mail:jizhenhua001@126.com

  • 中图分类号: TP394.1;TH691.9

摘要: 本文根据像差理论,设计一套计算次镜外反射的Offner型自由曲面棱镜光谱仪初始结构的算法。通过光线追迹获得光线在次镜外反射Offner型光谱仪各光学表面传播的公式,该公式可以确定光学元件的结构参数。应用轴外细光束像散理论分析系统产生的像散,并设定合理阈值作为结构算法的判断依据。在Matlab中迭代优化出符合设计要求的初始结构,获得的初始结构在Zemax中优化。为验证算法效果,本文设计光谱为380−780 nm,数值孔径为0.15,光谱分辨率为6 nm的自由曲面棱镜光谱仪的初始结构。在Zemax中完成优化后系统可达到设计指标且谱线弯曲和色畸变均优于0.1 pixels。设计结果表明采用本文的算法可以快速计算出符合要求的初始结构,大大简化了后续优化的复杂程度。

English Abstract

张佳伦, 郑玉权, 蔺超, 纪振华. 消像散的自由曲面棱镜光谱仪光学系统设计[J]. 中国光学. doi: 10.37188/CO.2019-0049
引用本文: 张佳伦, 郑玉权, 蔺超, 纪振华. 消像散的自由曲面棱镜光谱仪光学系统设计[J]. 中国光学. doi: 10.37188/CO.2019-0049
ZHANG Jia-lun, ZHENG Yu-qun, LIN Chao, JI Zhen-hua. Design of a Freeform Curved Prism Imaging Spectrometer Based on an Anastigmatism[J]. Chinese Optics. doi: 10.37188/CO.2019-0049
Citation: ZHANG Jia-lun, ZHENG Yu-qun, LIN Chao, JI Zhen-hua. Design of a Freeform Curved Prism Imaging Spectrometer Based on an Anastigmatism[J]. Chinese Optics. doi: 10.37188/CO.2019-0049
    • 高光谱技术是一种能准确获得地物二维图像信息和一维光谱信息的遥感手段于20世纪80年代提出,经过30多年的发展已经成为一个光学、机械、电子学等多学科交叉的综合研究方向[1]。自由曲面棱镜是指棱镜前后表面为曲面,两个曲面之间存在偏心或倾斜的棱镜。相比传统棱镜,自由曲面棱镜不需要准直汇聚镜组,有效的减小仪器了的体积与重量。自由曲面棱镜与Offner系统结合可以实现仪器小型化、轻量化设计。在德国ENMAP[2]及我国静止轨道全谱段高光谱探测等诸多项目均有重要应用。

      1911年,CH. Fery提出的Fery棱镜,本质上就是自由曲面棱镜,其通过罗兰圆原理确定元件参数,将单自由曲面棱镜应用于摄谱仪中,获得了较好的结果[3]。20世纪末,P. F. Morrissey基于三级像差理论得到将物点放置于自由曲面棱镜前表面的曲率中心或它的齐明共轭点时球差、彗差、像散具有最小值的结论,这一结果指导了早期自由曲面初始结构的设计[4]。2013年,刘力采用控制像散的方法通过光线追迹得到单自由曲面棱镜光谱仪的初始结构,结果表明能够在一定程度上减少优化复杂性[5]。2016年,刘立银通过光线追迹得到在Offner两臂增加自由曲面棱镜形式的初始结构,大大减小优化环节的工作量[6]。而作为自由曲面棱镜与Offner相结合的另一种形式,在次镜位置增加自由曲面棱镜,具有在相同设计指标下使用棱镜数量更少、容易装调等优点。目前该结构的设计大多基于共轴系统,再用光学设计软件不断优化调整,直至获得一个较好的结构。这就需要设计者具有大量的设计经验,因此设计一种能够快速计算初始结构的算法,尤为重要。

      本文针对这一问题,以次镜外反射型自由曲面棱镜Offner型光谱为基础形式,设计一套可以快速解算出符合要求的初始结构的算法,推导了光线在结构中各光学元件间的传播公式,获得了元件的曲率和偏心量,采用轴外细光束像散理论作为系统的像质评价手段,判断算法得到的结构是否符合设计要求。设计循环结构,经过多次迭代可以快速计算出一组性能较好的初始结构。

    • 图1为次镜外反射自由曲面棱镜结构图。次镜外反射型自由曲面棱镜即将自由曲面棱镜添加在Offner结构的次镜位置。在不破坏Offner结构对称性的前提下实现分光,相比在两臂添加自由曲面棱镜的情况,由于自由曲面棱镜处在系统对称位置,可以明显减少自由曲面棱镜的数量。本文通过分析次镜前的光路发现,在这部分计算中可以确定系统光学元件的曲率半径、偏心量等参数。次镜反射回的光线同样经过这些光学元件,故只需确定光线参数即可。

      图  1  次镜外反射型自由曲面棱镜光谱仪结构图

      Figure 1.  Structural diagram of the freeform curved prism imaging spectrometer

      在分析光学系统时需要规定各参量的正负性,并给出系统的符号确定规则[7]:

      (1)、光线与光轴的夹角:规定由光轴按照锐角旋转到光线方向,如果旋转方向沿顺时针方向,夹角为正,否则为负。

      (2)、光线与法线的夹角:规定由光线按照锐角旋转到法线方向,如果旋转方向沿顺时针方向,则夹角为正,否则为负。

      (3)、法线与光轴的夹角:规定由光轴按照锐角旋转到法线方向,如果旋转方向沿顺时针方向,则夹角为正,否则为负。

      对于Offner系统,当物点与球心处于同一位置时,光经过自由曲面棱镜后色散能力S的大小与棱镜前后两个面的偏心量差值$\Delta h$和棱镜的折射率的变化值$\Delta n$有关[6],即$S = - \Delta n \times \Delta h$。本文设计了光谱范围为380−780 nm,数值孔径为0.15,光谱分辨率为6 nm。使用成都光明的ZF6火石玻璃作为系统的色散元件,该波段材料的折射率变化值$\Delta n = 0.071$,选用的探测器像元尺寸为13 µm × 13 µm,经过计算后∆h = 7 mm。整理后系统参数如表1所示。

      表 1  设计指标

      Table 1.  Design specifications

      参数设计值
      光谱范围380−780 nm
      数值孔径0.15
      光谱分辨率6 nm
      像元尺寸13 μm × 13 μm
      色散元件材料ZF6
      色散面大小1 mm

      图2为该结构的次镜前光路图,通过分析光线在系统各个光学元件间的传播公式可以计算出系统各个元件的参数。根据大量的设计经验确定一些约束条件,这对于存在多个变量的离轴系统是可以接受的。在各光学元件表面分别建立局部坐标系[8],用以计算光线的参数,因此系统涉及大量的坐标变换。

      图  2  次镜前自由曲面棱镜光谱仪光路示意图

      Figure 2.  Schematic diagram of the freeform curved prism spectrometer light path before secondary mirror

      根据Offner的结构特点[9],当物点与主镜在Z方向的位置相同时会获得较好的设计结果[10,11]。因此对于初始结构,给出如下的规定:(1) 初始结构中物点与主镜的球心点位于同一水平位置;(2) 同心自由曲面棱镜存在垂直偏心量,其球心位于同一垂直位置;(3) 为了简化初始结构的计算,Offner结构主镜和三镜的参数相同。设定初始物距l1 = 200 mm,$ - {h_1}$为主镜的偏心量,$-{r_1}$为主镜的曲率半径,$-\Delta {h_{\rm{1}}}$为偏心量的修正量,$-\Delta {r_1}$为曲率半径的修正量。

      $$-{r_1} = -{l_1} - \Delta {r_1},$$ (1)
      $$-{h_1} = - {h_1} - \Delta {h_1}.$$ (2)

      选择平行光轴入射的主光线,根据Snell定律,$n_1$, $n_1'$为主镜物方像方折射率,${i_1}$, $i_1'$为主镜入射角和出射角,由于反射镜表面${n_1} = - n_{\rm{1}}'$,则${i_1} = - i_{\rm{1}}'$。因为系统只存在偏心量,光学元件的光轴相互平行,所以自由曲面棱镜第一个表面入射光线与光轴夹角${u_2}$与主镜出射光线与光轴的夹角$u_1'$相等,即${u_2} = u_1'$。主镜反射光线与光轴的交点为物点经过主镜的像点,系统具有较大偏心量不具备近轴条件,因此像距$l_1'$可以由$l_1' = - {h_1} \div \sin (u_1')$得出。在主镜的曲率中心处建立局部直角坐标系,求出主镜入射点${P_1}$的坐标(${{\rm{z}}_1},{y_1}$)。

      $${\rm{z}}_1^2 + y_1^2 = r_1^2,$$ (3)
      $${{\rm{y}}_1} = {h_1},$$ (4)

      式(3)、(4)联立可解出入射点${P_1}$的坐标,但此坐标是在主镜局部坐标系下求得,如果使用还要进行坐标变换。

      光线到达自由曲面棱镜第一个曲面处,此时光线经过主镜后的像点$A_1'$即为自由曲面棱镜的物点${A_2}$$ - {r_2}$为自由曲面棱镜第一个曲面的曲率半径,以${A_2}$为出发点,根据设计经验和设计指标约束自由曲面棱镜的空间位置${d_2} = {\rm{30\;mm}}$,同时给出修正量$\Delta {d_2}$。第一个曲面的曲率中心位置相对于主镜的曲率中心在水平方向存在偏离量$p{l_1}$。在垂直光轴方向存在偏离量$ - h_{12}'$,自由曲面棱镜第一曲面的偏心量为$-{h_2}$

      $$ - {r_2} = - {r_1} - p{l_1} - ( - l_1' + {d_2}),$$ (5)
      $$-{h_2} = - {h_1} - h_{12}'.$$ (6)

      为了解出光线在自由曲面棱镜第一曲面的入射角、折射角,首先需要知道入射点的位置。主镜上的入射点和第一曲面上的入射点入射光线上,而且光线与光轴的夹角已知。因此可以根据主镜上的入射点坐标与第一曲面的方程联立求出入射点的坐标。以第一曲面的曲率中心建立局部直角坐标系,则主镜的入射点${P_1}$坐标由(${{\rm{z}}_1},{y_1}$)转化为(${\rm{z}}_1',y_1'$)。

      $${\rm{z}}_1' = {{\textit{z}}_1} - p{l_1},$$ (7)
      $$y_1' = {y_1} - ( - h_{{\rm{12}}}'),$$ (8)
      $${\rm{z}}_2^2 + y_2^2 = r_{\rm{2}}^2,$$ (9)
      $${y_2} = \tan ({u_2})({{\textit{z}}_2} - {\textit{z}}_1') + y_1'.$$ (10)

      根据上述方程可以解出入射点${P_2}$坐标(${{\rm{z}}_2},{y_2}$),规定在第一曲面处光轴与法线的夹角为${\alpha _{\rm{2}}}$,连接物点与曲面曲率中心的连线作为辅助轴,辅助轴与光轴的夹角为${\beta _{\rm{2}}}$

      根据入射点可以解出:

      $${\alpha _{\rm{2}}} = \arcsin \left(\frac{{{{\rm{y}}_2}}}{{ - {r_2}}}\right),$$ (11)
      $${\beta _{\rm{2}}} = \frac{{ - h_{12}'}}{{ - {r_1} - p{l_1} - l_1'}}.$$ (12)

      由几何关系可知入射角为法线、光轴、入射光线构成的三角形的外角,存在$\left| {{i_2}} \right| = \left| {{\alpha _2}} \right| + \left| {{u_2}} \right|$。根据符号规则${i_2}$为正值,由Snell定律可以解出$i_2'$,存在出射光线与光轴的夹角$\left| {u_2'} \right| = \left| {i_2'} \right| - \left| {{\alpha _2}} \right|$。由几何关系可知$\angle {P_2}{C_2}A_{\rm{2}}' = {\beta _{\rm{2}}}-{\alpha _{\rm{2}}}$,在$\Delta {P_2}{C_2}A_{\rm{2}}'$中存在$\angle {C_2}A_2'{P_2} = \pi - \angle {P_2}{C_2}A_2' - i_2'$

      根据正弦定理在$\Delta {P_2}{C_2}A_{\rm{2}}'$中存在:

      $$\frac{{{P_2}{C_2}}}{{\angle {C_2}A_2'{P_2}}} = \frac{{{C_2}A_2'}}{{\angle {C_2}{P_2}A_2'}},$$ (13)

      ${C_2}A_2'$是第一曲面的曲率半径,$\angle {C_2}{P_2}A_2' = i_2'$。可以解出折射光线与辅助轴的交点即像点$A_2'$的像高$ - H_2' = {C_2}A_2' \times \sin ({\beta _2})$,像距$l_2' = - {r_2} - H_2' \times\cot ({\beta _2}) $

      光线到达自由曲面棱镜第二个曲面处,我们只需解出曲面入射点${P_3}$的坐标(${{\rm{z}}_3},{y_3}$)就可以采用上述方法获得折射光线的相关参数。$-{r_3}$为曲面的曲率半径,$-{h_3}$为曲面的偏心量,$ - h_{23}'$为自由曲面棱镜第一曲面和第二曲面的偏心量差值,根据设计指标这个值是$-h_{23}' = - \Delta h$。自由曲面棱镜两个曲面沿光轴方向的间距为d3 = 15 mm,修正量为$\Delta {d_3}$

      $$-{r_3} = - {r_2} + {d_3} + \Delta {{\rm{d}}_3},$$ (14)
      $$-{h_3} = - {h_2} - h_{23}'.$$ (15)

      入射点的坐标解算与第一曲面的相同,第一曲面的入射点与第二曲面的入射点都在折射光线上,折射光与光轴的夹角相等,即${u_3} = u_2'$。为了简化计算在第二曲面建立局部直角坐标系,解算入射点的坐标。

      $$z_2' = {z_2},$$ (16)
      $${\rm{y}}_{\rm{2}}^{\rm{'}} = {y_2} - h_{23}',$$ (17)
      $${\rm{z}}_3^2 + y_3^2 = r_3^2,$$ (18)
      $${y_3} = \tan ({u_3})({z_3} - z_2') + y_2',$$ (19)

      根据上述方程可以解出第二曲面入射点${P_3}$的坐标,采用在第一曲面相同的算法可以相继解出。

      $${\alpha _{\rm{3}}} = \arcsin \left(\frac{{{y_3}}}{{ - {r_3}}}\right),$$ (20)
      $${\beta _3} = {\rm{arctan}}\left(\frac{{{\rm{H}}_{\rm{2}}^{\rm{'}} - h_{23}'}}{{{r_3} - l_2'}}\right),$$ (21)
      $$\left| {u_{\rm{3}}'} \right| = \left| {i_{\rm{3}}'} \right| - \left| {{\alpha _{\rm{3}}}} \right|,$$ (22)
      $$ - H_3' = {C_3}A_3' \times \sin ({\beta _3}),$$ (23)
      $$l_3' = - {r_3} - H_3' \times \cot ({\beta _3}).$$ (24)

      光线到达次镜位置,$-{r_4}$为次镜的曲率半径,$-{h_4}$为次镜的偏心量,规定次镜与主镜曲率中心位置相同。入射光线与光轴的夹角相等,即${u_4} = u_3'$。次镜与自由曲面棱镜第二个曲面的间距为d4 = 10 mm,修正量为$\Delta {d_4}$

      $$-{r_4} = - {r_3} + {d_4} + \Delta {d_4}-p{l_1},$$ (25)
      $$-{h_4} = - {h_1}.$$ (26)

      次镜的入射点${P_4}$的坐标。

      $${\rm{z}}_3' = {z_3} + p{l_1} - {d_4} - \Delta {d_4},$$ (27)
      $${\rm{y}}_3' = {y_3} - {h_4} - ( - {h_3}),$$ (28)
      $${\rm{z}}_4^2 + y_4^2 = r_{\rm{4}}^2,$$ (29)
      $${y_4} = \tan ({u_{\rm{4}}})({z_4} - z_{\rm{3}}') + y_{\rm{3}}'.$$ (30)

      根据上述方程可以获得入射点${P_4}$的坐标(${z_{\rm{4}}},{y_{\rm{4}}}$)。

      光线经过次镜完成如图3所示的光路转折,次镜的反射光线还要经过自由曲面棱镜的两个面和三镜才能汇聚于像面,光线入射角、折射角等参数的计算方法相同,这里不予赘述。

      图  3  次镜后自由曲面棱镜光谱仪光路示意图

      Figure 3.  Schematic diagram of the freeform curved prism spectrometer light path after secondary mirror

    • 无限细像散光束的两个主截面方向的光束有各自的焦点,整个细光束形成两条相互垂直的短线,称为焦线。子午光束形成的子午焦线垂直于子午平面,弧矢光束形成的弧矢焦线在子午面内。像散产生的根本原因是,轴外视场光学系统在子午方向和弧矢方向的聚焦能力不同,从而导致弧矢焦面和子午焦面不重合[12,13]。实际上光学系统轴上点发出的细光束与光学系统的光轴重合,所以折射后仍保持同心。但是当细光束的光束轴与投射点法线不重合时,折射后的光束就不再是同心光束,如图4所示。

      图  4  轴外细光束像散示意图

      Figure 4.  Schematic diagram of off-axis rays astigmatism

      图中A为轴外物点,C为折射面的曲率中心,AC为系统的辅轴,系统子午物距与弧矢物距分别为ts,系统子午像距与弧矢像距分别为t′和s′。A发出的轴外细光束经过折射后首先汇聚与点t′处,发散最后汇聚在辅助轴上,两个像点分别对应轴外细光束的子午像点与弧矢像点。子午像点与弧矢像点之间的差值可以作为轴外细光束像散大小的评价标准。根据Snell定律和几何关系,n为物方折射率,n′为像方折射率,分别得到子午方向和弧矢方向的物像关系公式[14]

      $$\frac{{{n'}{{\cos }^2}{i'}}}{{{t'}}} - \frac{{n{{\cos }^2}i}}{t} = \frac{{{n'}\cos {i'} - n\cos i}}{r},$$ (31)
      $$\frac{{{n'}}}{{{s'}}} - \frac{n}{s} = \frac{{{n'}\cos {i'} - n\cos i}}{r},$$ (32)

      整理可以分别得到子午像距与弧矢像距。

      $${t'} = \frac{{{n'}{{\cos }^2}{i'}}}{{\dfrac{{n{{\cos }^2}i}}{t} + \dfrac{{{n'}\cos {i'} - n\cos i}}{r}}},$$ (33)
      $${s'} = \dfrac{{{n'}}}{{\dfrac{n}{s} - \dfrac{{{n'}\cos {i'} - n\cos i}}{r}}},$$ (34)

      子午像距与弧矢像距之间的距离。

      $$xt{s'} = {t'} - {s'},$$ (35)
      $$xt{s'} = \dfrac{{{n'}{{\cos }^2}{i'}\left(\frac{n}{s} - \dfrac{{{n'}\cos {i'} - n\cos i}}{r}\right) - {n'}\left(\dfrac{{n{{\cos }^2}i}}{t} + \dfrac{{{n'}\cos {i'} - n\cos i}}{r}\right)}}{{\left(\dfrac{{n{{\cos }^2}i}}{t} + \dfrac{{{n'}\cos {i'} - n\cos i}}{r}\right) \times \left(\dfrac{n}{s} - \dfrac{{{n'}\cos {i'} - n\cos i}}{r}\right)}}.$$ (36)

      如果系统初始子午物距与弧矢物距相同,即$t = s$。根据Snell定律可知折射角与入射角相关,那么影响像散大小的参数变为入射角、物方折射率、像方折射率、曲率半径。折射率可以由选择的材料确定,因此影响像散变化的因素只与入射角和曲率半径有关。

      在次镜外反射自由曲面棱镜光谱仪光学系统中,初始的子午物距弧矢物距相同,故可以通过追迹光线确定经过每一个光学元件表面后像散值的变化,这样既可以将最后输出的像散值作为系统的评判依据,也可以掌握每一个光学元件表面的像散变化情况。经过主镜得到子午像距$t_1'$弧矢像距$s_1'$$ts$为主镜入射点和自由曲面棱镜第一曲面入射点的间距,得到自由曲面棱镜第一曲面的子午物距${t_2}$与弧矢物距${s_2}$

      $$-ts = -\sqrt {{{(z_1' - {z_1})}^2} + {{(y_1' - {y_1})}^2}},$$ (37)
      $${t_2} = t_1' - \left( {{\rm{ - }}ts} \right),$$ (38)
      $${s_2} = s_1' - \left( {{\rm{ - }}ts} \right).$$ (39)

      将参数带入上述方程就可以解出经过第一曲面后的子午像距弧矢像距,进而得到具体像散值。对于每个光学元件表面,采用相同的方法多次求解可以得到系统的像散值。

    • 根据上述的分析,以光线追迹公式为基础,对次镜外反射结构进行消像散分析。在Matlab中实现程序设计,合理设置修正量,若满足设计要求则可以跳出循环,获得初始结构,如果不符合要求则调整参数继续循环。程序中存在大量的判断结构用于确定光线与光学元件的交点坐标。综合考虑,设计算法的流程,如图5所示。

      图  5  本文算法流程图

      Figure 5.  Flow chart of proposed algorithm

    • 根据具体的循环迭代情况,不断调整修正量的值,使像散值在循环过程中趋于收敛。经过多次循环,系统找到一组符合条件的结构参数,如表2所示。

      表 2  初始结构数据(单位:毫米)

      Table 2.  Initial structural data (Unit: mm)

      表面半径厚度材料偏心量
      1——249.500 000
      2−248.510 000−104.255 000MIRROR−45.990 000
      3−138.305 000−15.990 000ZF6−36.980 000
      4−122.315 000−5.990 000−43.980 000
      5−122.275 0005.990 000MIRROR−45.990 000
      6−122.315 00015.990 000ZF6−43.980 000
      7−138.305 000104.255 000−36.980 000
      8−248.510 000——MIRROR−45.990 000
      9————

      将初始结构输入到Zemax中可以得到光学结构,如图6所示。

      图  6  初始自由曲面棱镜光谱仪结构示意图

      Figure 6.  Schematic diagram of the initial freeform curved prism spectrometer

      对初始结构进行优化,在结构解算中,为了简化系统的复杂程度,设定主镜和三镜参数相同,在具体优化中可以对他们分别进行优化。这样经过几次优化即可得到符合设计要求的光学系统。优化后的参数,如表3所示。

      表 3  优化后结构数据(单位:毫米)

      Table 3.  Optimized structural data (Unit: mm)

      表面半径厚度材料偏心量
      1——240.045 897
      2−234.113 632−99.964 613MIRROR−58.302 387
      3−138.439 720−18.303 653ZF6−42.950 734
      4−120.219 597−7.986 594−49.982 314
      5−105.485 6317.986 594MIRROR−60.171 600
      6−120.219 59718.303 653ZF6−49.982 314
      7−138.439 72089.943 651−42.950 734
      8−223.969 005——MIRROR−61.713 353
      9————

      优化后光学结构,如图7所示。

      图  7  优化后的自由曲面棱镜光谱仪结构示意图

      Figure 7.  Schematic diagram of the optimized freeform curved prism spectrometer

      一个好的光学设计很大程度上取决于初始结构的合理性,使用本文算法解算出的初始结构可以快速优化出符合设计要求的光学结构。同时在追迹光线的过程中,由于计算了光线在每一个表面的像散值,因此可以更容易选择合适的优化路线[15,16]。如图8所示,优化后系统各波段的调制传递函数均接近衍射极限。平均光谱分辨率优于6 nm,谱线弯曲和色畸变均小于0.1 pixel,符合设计要求。

      图  8  光学系统调制传递函数

      Figure 8.  Schematic diagram of MTF

    • 本文根据消像散方法,以轴外细光束像散理论作为判断依据,计算了光线经过每个面后的像散值。通过光线追迹确定光学元件参数,在结构参数解算中,设定一系列修正量,通过循环迭代优化,找到合理的初始结构,将初始结构放到Zemax中进行优化,从而获得符合设计要求光学结构。为了验证算法的正确性,用其设计光谱范围380−780 nm,光谱分辨率为6 nm,数值孔径为0.15的光谱仪。根据上述方法写出Matlab程序,迭代出一个初始结构,在Zemax中快速优化出的光学结构较好。结果证明,采用该算法可以根据设计指标快速计算出合理的初始结构,能够大大减化优化过程。

参考文献 (16)

目录

    /

    返回文章
    返回