CN110220477A - 红外干涉仪测量光学非球面反射镜面形的检测装置及方法 - Google Patents

红外干涉仪测量光学非球面反射镜面形的检测装置及方法 Download PDF

Info

Publication number
CN110220477A
CN110220477A CN201910534152.2A CN201910534152A CN110220477A CN 110220477 A CN110220477 A CN 110220477A CN 201910534152 A CN201910534152 A CN 201910534152A CN 110220477 A CN110220477 A CN 110220477A
Authority
CN
China
Prior art keywords
aspherical
cos
sin
measured
standard
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910534152.2A
Other languages
English (en)
Other versions
CN110220477B (zh
Inventor
乐丽珠
朱光宇
王贵全
李晓斌
杨静
邵毅
姜兴平
张伟涛
于振龙
王春兴
王乔方
李在丙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
KUNMING NORTH INFRARED TECHNOLOGY Co Ltd
Original Assignee
KUNMING NORTH INFRARED TECHNOLOGY Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by KUNMING NORTH INFRARED TECHNOLOGY Co Ltd filed Critical KUNMING NORTH INFRARED TECHNOLOGY Co Ltd
Priority to CN201910534152.2A priority Critical patent/CN110220477B/zh
Publication of CN110220477A publication Critical patent/CN110220477A/zh
Application granted granted Critical
Publication of CN110220477B publication Critical patent/CN110220477B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/2441—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)

Abstract

红外干涉仪测量光学非球面反射镜面形的检测装置及方法，涉及光学检测技术领域，尤其是一利用红外干涉仪进行检测的红外干涉仪测量光学非球面反射镜面形的检测装置及方法。本发明的红外干涉仪测量光学非球面反射镜面形的检测装置及方法，该检测装置包括红外激光干涉仪、空间滤波器、第一准直镜、参考反射镜、分光镜、标准镜头、标准镜头夹持器、待测镜头夹持器、第二准直镜、红外探测器、移相器、导轨和计算机；本发明为生产、科研部门提供了一种非球面ΡV、RMS等波面参数的激光干涉法测量方法，为非球面波面参数评价提供了有效可靠的依据。

Description

红外干涉仪测量光学非球面反射镜面形的检测装置及方法
技术领域
本发明涉及光学检测技术领域，尤其是一利用红外干涉仪进行检测，适合批量生产的在线检测的红外干涉仪测量光学非球面反射镜面形的检测装置及方法。
背景技术
在光学***中采用非球面光学元件，不仅可以减少光学元件的数量和减轻重量，缩小***尺寸，简化仪器结构，而且有利于校正像差，改善成像质量，增加光学设计的自由度，降低成本。
近年来，由于锗Ge等红外材料的特殊加工特性及单点金刚石车削技术的发展，红外光学非球面元件加工技术及其应用得到了快速的发展，与此同时，对光学非球面元件的面形参数的检测提出了更高的要求。尽管目前发展出多种非球面的检测技术，但非球面面形检测仍然是光学检测中的难点和热点。目前，非球面元件面形的测量方法主要分为接触式测量法、阴影法及激光干涉法。
接触式测量主要借助轮廓仪或三坐标测量仪等设备对非球面元件进行多个离散点的测量，然后经过数据处理，拟合得到面形误差。由于离散点间的空白区域在扫描过程中信息缺失，导致非球面元件的表面特征无法仅通过多个离散点来直观地表征，而且即使光学元件具有相似二维的表面，其三维形貌也会千差万别。而且设备测量头与非球面元件的接触可能会给表面带来一定的划痕。
阴影法测量主要观察阴影分布的图形和阴影图的明暗对比，对非球面元件的面形继续定性检测，该方法存在主观、定量困难、灵敏度低等缺点。
传统的激光干涉法检测光学非球面反射面形要借助补偿装置，有的补偿器本身可能还有非球面元件，而要实现非球面的高精度检测，还要为其设计补偿器，因此又给加工和装调带来很多困难。另外检测不同非球面反射镜所用的补偿装置各不相同、通用性差，且补偿装置加工成本高，而且并不能单独测量补偿器的补偿效果。
发明内容
本发明所要解决的就是现有非球面反射镜面形测量困难的问题，提供一种用红外干涉仪进行检测，适合批量生产的在线检测的红外干涉仪测量光学非球面反射镜面形的检测装置及方法。
本发明的红外干涉仪测量光学非球面反射镜面形的检测装置及方法，其特征在于该检测装置包括红外激光干涉仪、空间滤波器、第一准直镜、参考反射镜、分光镜、标准镜头、标准镜头夹持器、待测镜头夹持器、第二准直镜、红外探测器、移相器、导轨和计算机，红外激光干涉仪、空间滤波器、第一准直镜、分光镜、标准镜头夹持器和导轨从左至右设置，标准镜头放置在标准镜头夹持器上，待测镜头夹持器放置在导轨上沿导轨滑动，待测的非球面反射镜放置在待测镜头夹持器上，红外激光干涉仪、空间滤波器、第一准直镜、标准镜头以及待测的非球面反射镜处于同一光轴上，分光镜的中点位于该光轴上；参考反射镜设置在分光镜上方，第二准直镜设置在分光镜下方，红外探测器设置在第二准直镜下方；移相器与标准镜头夹持器连接，带动标准镜头夹持器移动，计算机分别与红外探测器、移相器连接；
具体检测步骤为：
步骤S1：测量待测非球面反射镜与标准球面反射镜的面形误差
将相应的标准球面镜安装在标准镜头夹持器上，待测非球面反射镜安装在待测镜头夹持器上，调整红外激光干涉仪、待测非球面反射镜的位置和角度，使红外激光干涉仪射出的标准球面波入射到待测非球面反射镜表面上，经过非球面反射镜反射后按原路返回，与红外激光干涉仪出射的标准球面波生产干涉形成干涉条纹，通过红外激光干涉仪测试，得到非球面反射镜与标准球面的波相差W1(x，y)；
步骤S2：导出待测非球面反射镜与标准球面镜面形误差的Zernike系数矩阵
通过测试***自带的红外激光干涉仪测试软件导出非球面反射镜与标准球面的波相差W1(x，y)的36项Zernike系数矩阵Z1,Z2,Z3……Z36；
步骤S3：求解非球面最接近球参数(R,a)
在步骤S1的波相差W1(x，y)坐标系中，非球面的矢高方程为：
式中，x为离非球面回转轴的径向距离；y为相应的垂直距离；c为顶点曲率，c＝1/R0，R0为顶点曲率半径；k为二次常数；an为非球面系数；
非球面度A计算：
式中(x,y)表示非球面上的点，a为非球面最接近球面的球心位置，R为非球面最接近球面半径，最大非球面度Amax为非球面度A在从中心到边缘所有位置中的最大值，非球面的最接近球面是指满足最大非球面度最小的球参数(R，a),即求解
求解Amax(R,a)＝max{A(R,a,x)}＝min的球参数(R，a)，输入被测非球面口径及非球面参数(包括顶点曲率半径R0、二次常数k、非球面系数an,如(1)式所示)，以非球面顶点R0作为起始点，搜索求解得到非球面方程满足最大非球面度A最小的球参数(R，a)；
步骤S4：计算待测非球面最接近球面偏差曲线方程
将计算得到的非球面最接近球面的参数(R，a)代入(2)式，求解得到非球面最接近球面偏差曲线；
步骤S5：利用多项式方程拟合非球面最接近球面偏差曲线方程
非球面最接近球面偏差曲线方程为偶次方多项式，采用下式所示的偶次方多项式方程拟和非球面最接近球面偏差曲线
y1(x)＝a0+a1x2+a2x4+a3x6+a4x8+a5x10 (3)
对X轴进行归一化处理，用(3)式拟合步骤S4得到的非球面最接近球面偏差曲线，求解偏差曲线的拟合系数a0、a1、a2、a3、a4、a5；
步骤S6：计算待测非球面最接近球面偏差曲线的泽尼克系数矩阵Q
36项泽尼克系数矩阵表达式如下所示:
Z0＝1；
Z1＝ρ*cos(a)；
Z2＝ρ*sin(a)；
Z3＝2*ρ^2-1；
Z4＝ρ^2*cos(2*a)；
Z5＝ρ^2*sin(2*a)；
Z6＝(3*ρ^2-2)*ρ*cos(a)；
Z7＝(3*ρ^2-2)*ρ*sin(a)；
Z8＝6*ρ^4-6*ρ^2+1；
Z9＝ρ^3*cos(3*a)；
Z10＝ρ^3*sin(3*a)；
Z11＝(4*ρ^2-3)*ρ^2*cos(2*a)；
Z12＝(4*ρ^2-3)*ρ^2*sin(2*a)；
Z13＝(10*ρ^4-12*ρ^2+3)*ρ*cos(a)；
Z14＝(10*ρ^4-12*ρ^2+3)*ρ*sin(a)；
Z15＝20*ρ^6-30*ρ^4+12*ρ^2-1；
Z16＝ρ^4*cos(4*a)；
Z17＝ρ^4*sin(4*a)；
Z18＝(5*ρ^2-4)*ρ^3*cos(3*a)；
Z19＝(5*ρ^2-4)*ρ^3*sin(3*a)；
Z20＝(15*ρ^4-20*ρ^2+6)*ρ^2*cos(2*a)；
Z21＝(15*ρ^4-20*ρ^2+6)*ρ^2*sin(2*a)；
Z22＝(35*ρ^6-60*ρ^4+30*ρ^2-4)*ρ*cos(a)；
Z23＝(35*ρ^6-60*ρ^4+30*ρ^2-4)*ρ*sin(a)；
Z24＝70*ρ^8-140*ρ^6+90*ρ^4-20*ρ^2+1；
Z25＝ρ^5*cos(5*a)；
Z26＝ρ^5*sin(5*a)；
Z27＝(6*ρ^2-5)*ρ^4*cos(4*a)；
Z28＝(6*ρ^2-5)*ρ^4*sin(4*a)；
Z29＝(21*ρ^4-30*ρ^2+10)*ρ^3*cos(3*a)；
Z30＝(21*ρ^4-30*ρ^2+10)*ρ^3*sin(3*a)；
Z31＝(56*ρ^6-105*ρ^4+60*ρ^2-10)*ρ^2*cos(2*a)；
Z32＝(56*ρ^6-105*ρ^4+60*ρ^2-10)*ρ^2*sin(2*a)；
Z33＝(126*ρ^8-280*ρ^6+210*ρ^4-60*ρ^2+5)*ρ*cos(a)；
Z34＝(126*ρ^8-280*ρ^6+210*ρ^4-60*ρ^2+5)*ρ*sin(a)；
Z35＝252*ρ^10-630*ρ^8+560*ρ^6-210*ρ^4+30*ρ^2-1；
由于理论非球面与最接近球面的偏差曲面是三维旋转对称曲面，由旋转对称性可得到：36项Zernike系数项中只有第1项Z0、第4项Z3、第9项Z8、第16项Z15、第25项Z24和第36项Z35系数不为零，其余项的系数全部为零，a0与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有常数项相对应,a1与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有二次项相对应,a2与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有四次项相对应通过，a3与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有六次项相对应，a4与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有八次项相对应，a5与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有十次项相对应，通过求解a0、a1、a2、a3、a4、a5与36项Zernike系数项中第1项、第4项、第9项、第16项、第25项和第36项系数的对应关系，可得到理论非球面与最接近球面偏差的Zernike系数项Q1、Q4、Q9、Q16、Q25、Q36，其余项的系数全部为零；
步骤S7：系数矩阵Z与系数矩阵Q之差求解
将干涉仪测试软件导出的波相差W1(x，y)的36项Zernike系数矩阵Z与计算得到的非球面与最接近球面偏差的Zernike系数矩阵Q的对应项相减，得到新的Zernike系数矩阵T；
步骤S8：被测非球面的三维形貌
对步骤S7得到的Zernike系数矩阵T进行三维拟合，得到被测非球面的面形参数。
本发明提供的红外干涉法测量光学非球面反射镜面形的检测方法，为生产、科研部门提供了一种非球面ΡV、RMS等波面参数的激光干涉法测量方法，为非球面波面参数评价提供了有效可靠的依据，测量过程中还具有以下有益效果：
1)采用基于红外干涉仪测量光学非球面反射镜面形的检测方法，实现了红外干涉法测量光学非球面反射镜面形的定量检测，数据自动采集、分析、计算、存储和打印，***运行稳定，操作使用方便；
2)本发明给出的测试方法具有操作简单、测值可靠、快速，适合批量生产的在线检测。
附图说明
图1为本发明检测装置原理图。
图2为本发明激光干涉仪测试结果。
图3为本发明非球面与最接近球面图。
图4球参数(R，a)求解程序流程图。
图5非球面最接近球面偏差曲线。
图6最接近球面偏差曲线的拟合系数求解程序流程图。
图7为本发明非球面最接近球面偏差三维形貌。
图8本发明非球面最接近球面偏差三维形貌求解程序流程图。
图9为本发明被测非球面三维形貌图。
图10为本发明被测非球面三维形貌求解程序流程图。
其中，红外激光干涉仪1，空间滤波器2，第一准直镜3，参考反射镜4，分光镜5，标准镜头6，标准镜头夹持器7，非球面反射镜8，待测镜头夹持器9，第二准直镜10，红外探测器11，移相器12，导轨13，计算机14。
具体实施方式
实施例1：一种红外干涉仪测量光学非球面反射镜面形的检测装置及方法，其特征在于该检测装置包括红外激光干涉仪、空间滤波器、第一准直镜、参考反射镜、分光镜、标准镜头、标准镜头夹持器、待测镜头夹持器、第二准直镜、红外探测器、移相器、导轨和计算机，红外激光干涉仪、空间滤波器、第一准直镜、分光镜、标准镜头夹持器和导轨从左至右设置，标准镜头放置在标准镜头夹持器上，待测镜头夹持器放置在导轨上沿导轨滑动，待测的非球面反射镜放置在待测镜头夹持器上，红外激光干涉仪、空间滤波器、第一准直镜、标准镜头以及待测的非球面反射镜处于同一光轴上，分光镜的中点位于该光轴上；参考反射镜设置在分光镜上方，第二准直镜设置在分光镜下方，红外探测器设置在第二准直镜下方；移相器与标准镜头夹持器连接，带动标准镜头夹持器移动，计算机分别与红外探测器、移相器连接；
具体检测步骤为：
步骤S1：测量待测非球面反射镜与标准球面反射镜的面形误差
将相应的标准球面镜安装在标准镜头夹持器上，待测非球面反射镜安装在待测镜头夹持器上，调整红外激光干涉仪、待测非球面反射镜的位置和角度，使红外激光干涉仪射出的标准球面波入射到待测非球面反射镜表面上，经过非球面反射镜反射后按原路返回，与红外激光干涉仪出射的标准球面波生产干涉形成干涉条纹，通过红外激光干涉仪测试，得到非球面反射镜与标准球面的波相差W1(x，y)；
步骤S2：导出待测非球面反射镜与标准球面镜面形误差的Zernike系数矩阵
通过测试***自带的红外激光干涉仪测试软件导出非球面反射镜与标准球面的波相差W1(x，y)的36项Zernike系数矩阵Z1,Z2,Z3……Z36；
步骤S3：求解非球面最接近球参数(R,a)
在步骤S1的波相差W1(x，y)坐标系中，非球面的矢高方程为：
式中，x为离非球面回转轴的径向距离；y为相应的垂直距离；c为顶点曲率，c＝1/R0，R0为顶点曲率半径；k为二次常数；an为非球面系数；
非球面度A计算：
式中(x,y)表示非球面上的点，a为非球面最接近球面的球心位置，R为非球面最接近球面半径，最大非球面度Amax为非球面度A在从中心到边缘所有位置中的最大值，非球面的最接近球面是指满足最大非球面度最小的球参数(R，a),即求解
求解Amax(R,a)＝max{A(R,a,x)}＝min的球参数(R，a)，输入被测非球面口径及非球面参数(包括顶点曲率半径R0、二次常数k、非球面系数an,如(1)式所示)，以非球面顶点R0作为起始点，搜索求解得到非球面方程满足最大非球面度A最小的球参数(R，a)；非球面方程满足最大非球面度A最小的球参数(R，a)
对于口径为D＝110mm，R0＝252.9，K＝0.434的非球面，通过步骤S3求解，得到最大非球面度A最小的球参数R＝251.5795，a＝251.5795；
步骤S4：计算待测非球面最接近球面偏差曲线方程
将计算得到的非球面最接近球面的参数(R，a)代入(2)式，求解得到非球面最接近球面偏差曲线；
步骤S5：利用多项式方程拟合非球面最接近球面偏差曲线方程
非球面最接近球面偏差曲线方程为偶次方多项式，采用下式所示的偶次方多项式方程拟和非球面最接近球面偏差曲线
y1(x)＝a0+a1x2+a2x4+a3x6+a4x8+a5x10 (3)
对X轴进行归一化处理，用(3)式拟合步骤S4得到的非球面最接近球面偏差曲线，求解偏差曲线的拟合系数a0、a1、a2、a3、a4、a5；
步骤S6：计算待测非球面最接近球面偏差曲线的泽尼克系数矩阵Q
36项泽尼克系数矩阵表达式如下所示:
Z0＝1；
Z1＝ρ*cos(a)；
Z2＝ρ*sin(a)；
Z3＝2*ρ^2-1；
Z4＝ρ^2*cos(2*a)；
Z5＝ρ^2*sin(2*a)；
Z6＝(3*ρ^2-2)*ρ*cos(a)；
Z7＝(3*ρ^2-2)*ρ*sin(a)；
Z8＝6*ρ^4-6*ρ^2+1；
Z9＝ρ^3*cos(3*a)；
Z10＝ρ^3*sin(3*a)；
Z11＝(4*ρ^2-3)*ρ^2*cos(2*a)；
Z12＝(4*ρ^2-3)*ρ^2*sin(2*a)；
Z13＝(10*ρ^4-12*ρ^2+3)*ρ*cos(a)；
Z14＝(10*ρ^4-12*ρ^2+3)*ρ*sin(a)；
Z15＝20*ρ^6-30*ρ^4+12*ρ^2-1；
Z16＝ρ^4*cos(4*a)；
Z17＝ρ^4*sin(4*a)；
Z18＝(5*ρ^2-4)*ρ^3*cos(3*a)；
Z19＝(5*ρ^2-4)*ρ^3*sin(3*a)；
Z20＝(15*ρ^4-20*ρ^2+6)*ρ^2*cos(2*a)；
Z21＝(15*ρ^4-20*ρ^2+6)*ρ^2*sin(2*a)；
Z22＝(35*ρ^6-60*ρ^4+30*ρ^2-4)*ρ*cos(a)；
Z23＝(35*ρ^6-60*ρ^4+30*ρ^2-4)*ρ*sin(a)；
Z24＝70*ρ^8-140*ρ^6+90*ρ^4-20*ρ^2+1；
Z25＝ρ^5*cos(5*a)；
Z26＝ρ^5*sin(5*a)；
Z27＝(6*ρ^2-5)*ρ^4*cos(4*a)；
Z28＝(6*ρ^2-5)*ρ^4*sin(4*a)；
Z29＝(21*ρ^4-30*ρ^2+10)*ρ^3*cos(3*a)；
Z30＝(21*ρ^4-30*ρ^2+10)*ρ^3*sin(3*a)；
Z31＝(56*ρ^6-105*ρ^4+60*ρ^2-10)*ρ^2*cos(2*a)；
Z32＝(56*ρ^6-105*ρ^4+60*ρ^2-10)*ρ^2*sin(2*a)；
Z33＝(126*ρ^8-280*ρ^6+210*ρ^4-60*ρ^2+5)*ρ*cos(a)；
Z34＝(126*ρ^8-280*ρ^6+210*ρ^4-60*ρ^2+5)*ρ*sin(a)；
Z35＝252*ρ^10-630*ρ^8+560*ρ^6-210*ρ^4+30*ρ^2-1；
由于理论非球面与最接近球面的偏差曲面是三维旋转对称曲面，由旋转对称性可得到：36项Zernike系数项中只有第1项Z0、第4项Z3、第9项Z8、第16项Z15、第25项Z24和第36项Z35系数不为零，其余项的系数全部为零，a0与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有常数项相对应,a1与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有二次项相对应,a2与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有四次项相对应通过，a3与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有六次项相对应，a4与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有八次项相对应，a5与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有十次项相对应，通过求解a0、a1、a2、a3、a4、a5与36项Zernike系数项中第1项、第4项、第9项、第16项、第25项和第36项系数的对应关系，可得到理论非球面与最接近球面偏差的Zernike系数项Q1、Q4、Q9、Q16、Q25、Q36，其余项的系数全部为零；
对于口径为D＝110mm，R0＝252.9，K＝0.434的非球面，通过步骤S6计算得到，Q1＝0.00532367，Q1＝0，Q3＝0，Q4＝0.0000561，Q5＝0，Q5＝0，Q5＝0，Q5＝0，Q9＝-0.00532295，Q10＝0，Q11＝0，Q12＝0，Q13＝0，Q14＝0，Q15＝0，Q16＝-0.00005607，Q17＝0，Q18＝0，Q19＝0，Q20＝0，Q21＝0，Q22＝0，Q23＝0，Q24＝0，Q25＝-0.00000073，Q26＝0，Q27＝0，Q28＝0，Q29＝0，Q30＝0，Q31＝0，Q32＝0，Q33＝0，Q34＝0，Q35＝0，Q36＝-0.00000001；
步骤S7：系数矩阵Z与系数矩阵Q之差求解
将干涉仪测试软件导出的波相差W1(x，y)的36项Zernike系数矩阵Z与计算得到的非球面与最接近球面偏差的Zernike系数矩阵Q的对应项相减，得到新的Zernike系数矩阵T；
步骤S8：被测非球面的三维形貌
对步骤S7得到的Zernike系数矩阵T进行三维拟合，得到被测非球面的面形参数；
通过步骤S7计算得到，Zernike系数矩阵T：T1＝0.0756，T2＝0.05，T3＝0，，T4＝0.048，T5＝0.001，T6＝-0.002，T7＝0.008，T8＝-0.004，T9＝-0.039，T10＝0，T11＝-0.001，T123＝0.001，T13＝0，T14＝-0.002，T15＝0.002，T16＝-0.002，T17＝0.002，T18＝0.001，T19＝0，T20＝0，T21＝-0.001，T22＝0，T23＝-0.002，T24＝0.002，T25＝0.005，T26＝0，T27＝0，T28＝0.001，T29＝0.001，T30＝0，T31＝0.001，T32＝-0.001，T33＝-0.001，T34＝-0，T35＝0.001，T36＝-0.004。通过步骤S8计算得到被测非球面的面形参数。

Claims (1)

1.一种红外干涉仪测量光学非球面反射镜面形的检测装置及方法，其特征在于该检测装置包括红外激光干涉仪、空间滤波器、第一准直镜、参考反射镜、分光镜、标准镜头、标准镜头夹持器、待测镜头夹持器、第二准直镜、红外探测器、移相器、导轨和计算机，红外激光干涉仪、空间滤波器、第一准直镜、分光镜、标准镜头夹持器和导轨从左至右设置，标准镜头放置在标准镜头夹持器上，待测镜头夹持器放置在导轨上沿导轨滑动，待测的非球面反射镜放置在待测镜头夹持器上，红外激光干涉仪、空间滤波器、第一准直镜、标准镜头以及待测的非球面反射镜处于同一光轴上，分光镜的中点位于该光轴上；参考反射镜设置在分光镜上方，第二准直镜设置在分光镜下方，红外探测器设置在第二准直镜下方；移相器与标准镜头夹持器连接，带动标准镜头夹持器移动，计算机分别与红外探测器、移相器连接；
具体检测步骤为：
步骤S1：测量待测非球面反射镜与标准球面反射镜的面形误差
将相应的标准球面镜安装在标准镜头夹持器上，待测非球面反射镜安装在待测镜头夹持器上，调整红外激光干涉仪、待测非球面反射镜的位置和角度，使红外激光干涉仪射出的标准球面波入射到待测非球面反射镜表面上，经过非球面反射镜反射后按原路返回，与红外激光干涉仪出射的标准球面波生产干涉形成干涉条纹，通过红外激光干涉仪测试，得到非球面反射镜与标准球面的波相差W1(x，y)；
步骤S2：导出待测非球面反射镜与标准球面镜面形误差的Zernike系数矩阵
通过测试***自带的红外激光干涉仪测试软件导出非球面反射镜与标准球面的波相差W1(x，y)的36项Zernike系数矩阵Z1,Z2,Z3……Z36；
步骤S3：求解非球面最接近球参数(R,a)
在步骤S1的波相差W1(x，y)坐标系中，非球面的矢高方程为：
式中，x为离非球面回转轴的径向距离；y为相应的垂直距离；c为顶点曲率，c＝1/R0，R0为顶点曲率半径；k为二次常数；an为非球面系数；
非球面度A计算：
式中(x,y)表示非球面上的点，a为非球面最接近球面的球心位置，R为非球面最接近球面半径，最大非球面度Amax为非球面度A在从中心到边缘所有位置中的最大值，非球面的最接近球面是指满足最大非球面度最小的球参数(R，a),即求解
求解Amax(R,a)＝max{A(R,a,x)}＝min的球参数(R，a)，输入被测非球面口径及非球面参数(包括顶点曲率半径R0、二次常数k、非球面系数an,如(1)式所示)，以非球面顶点R0作为起始点，搜索求解得到非球面方程满足最大非球面度A最小的球参数(R，a)；
步骤S4：计算待测非球面最接近球面偏差曲线方程
将计算得到的非球面最接近球面的参数(R，a)代入(2)式，求解得到非球面最接近球面偏差曲线；
步骤S5：利用多项式方程拟合非球面最接近球面偏差曲线方程
非球面最接近球面偏差曲线方程为偶次方多项式，采用下式所示的偶次方多项式方程拟和非球面最接近球面偏差曲线
y1(x)＝a0+a1x2+a2x4+a3x6+a4x8+a5x10 (3)
对X轴进行归一化处理，用(3)式拟合步骤S4得到的非球面最接近球面偏差曲线，求解偏差曲线的拟合系数a0、a1、a2、a3、a4、a5；
步骤S6：计算待测非球面最接近球面偏差曲线的泽尼克系数矩阵Q
36项泽尼克系数矩阵表达式如下所示:
Z0＝1；
Z1＝ρ*cos(a)；
Z2＝ρ*sin(a)；
Z3＝2*ρ^2-1；
Z4＝ρ^2*cos(2*a)；
Z5＝ρ^2*sin(2*a)；
Z6＝(3*ρ^2-2)*ρ*cos(a)；
Z7＝(3*ρ^2-2)*ρ*sin(a)；
Z8＝6*ρ^4-6*ρ^2+1；
Z9＝ρ^3*cos(3*a)；
Z10＝ρ^3*sin(3*a)；
Z11＝(4*ρ^2-3)*ρ^2*cos(2*a)；
Z12＝(4*ρ^2-3)*ρ^2*sin(2*a)；
Z13＝(10*ρ^4-12*ρ^2+3)*ρ*cos(a)；
Z14＝(10*ρ^4-12*ρ^2+3)*ρ*sin(a)；
Z15＝20*ρ^6-30*ρ^4+12*ρ^2-1；
Z16＝ρ^4*cos(4*a)；
Z17＝ρ^4*sin(4*a)；
Z18＝(5*ρ^2-4)*ρ^3*cos(3*a)；
Z19＝(5*ρ^2-4)*ρ^3*sin(3*a)；
Z20＝(15*ρ^4-20*ρ^2+6)*ρ^2*cos(2*a)；
Z21＝(15*ρ^4-20*ρ^2+6)*ρ^2*sin(2*a)；
Z22＝(35*ρ^6-60*ρ^4+30*ρ^2-4)*ρ*cos(a)；
Z23＝(35*ρ^6-60*ρ^4+30*ρ^2-4)*ρ*sin(a)；
Z24＝70*ρ^8-140*ρ^6+90*ρ^4-20*ρ^2+1；
Z25＝ρ^5*cos(5*a)；
Z26＝ρ^5*sin(5*a)；
Z27＝(6*ρ^2-5)*ρ^4*cos(4*a)；
Z28＝(6*ρ^2-5)*ρ^4*sin(4*a)；
Z29＝(21*ρ^4-30*ρ^2+10)*ρ^3*cos(3*a)；
Z30＝(21*ρ^4-30*ρ^2+10)*ρ^3*sin(3*a)；
Z31＝(56*ρ^6-105*ρ^4+60*ρ^2-10)*ρ^2*cos(2*a)；
Z32＝(56*ρ^6-105*ρ^4+60*ρ^2-10)*ρ^2*sin(2*a)；
Z33＝(126*ρ^8-280*ρ^6+210*ρ^4-60*ρ^2+5)*ρ*cos(a)；
Z34＝(126*ρ^8-280*ρ^6+210*ρ^4-60*ρ^2+5)*ρ*sin(a)；
Z35＝252*ρ^10-630*ρ^8+560*ρ^6-210*ρ^4+30*ρ^2-1；
由于理论非球面与最接近球面的偏差曲面是三维旋转对称曲面，由旋转对称性可得到：36项Zernike系数项中只有第1项Z0、第4项Z3、第9项Z8、第16项Z15、第25项Z24和第36项Z35系数不为零，其余项的系数全部为零，a0与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有常数项相对应,a1与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有二次项相对应,a2与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有四次项相对应通过，a3与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有六次项相对应，a4与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有八次项相对应，a5与Z0、Z3、Z8、Z15、Z24、Z35共6项的所有十次项相对应，通过求解a0、a1、a2、a3、a4、a5与36项Zernike系数项中第1项、第4项、第9项、第16项、第25项和第36项系数的对应关系，可得到理论非球面与最接近球面偏差的Zernike系数项Q1、Q4、Q9、Q16、Q25、Q36，其余项的系数全部为零；
步骤S7：系数矩阵Z与系数矩阵Q之差求解
将干涉仪测试软件导出的波相差W1(x，y)的36项Zernike系数矩阵Z与计算得到的非球面与最接近球面偏差的Zernike系数矩阵Q的对应项相减，得到新的Zernike系数矩阵T；
步骤S8：被测非球面的三维形貌
对步骤S7得到的Zernike系数矩阵T进行三维拟合，得到被测非球面的面形参数。
CN201910534152.2A 2019-06-19 2019-06-19 红外干涉仪测量光学非球面反射镜面形的检测装置及方法 Active CN110220477B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910534152.2A CN110220477B (zh) 2019-06-19 2019-06-19 红外干涉仪测量光学非球面反射镜面形的检测装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910534152.2A CN110220477B (zh) 2019-06-19 2019-06-19 红外干涉仪测量光学非球面反射镜面形的检测装置及方法

Publications (2)

Publication Number Publication Date
CN110220477A true CN110220477A (zh) 2019-09-10
CN110220477B CN110220477B (zh) 2020-09-29

Family

ID=67814247

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910534152.2A Active CN110220477B (zh) 2019-06-19 2019-06-19 红外干涉仪测量光学非球面反射镜面形的检测装置及方法

Country Status (1)

Country Link
CN (1) CN110220477B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1995943A (zh) * 2007-01-04 2007-07-11 四川大学 一种大口径非球面镜全场检测方法
CN101709955A (zh) * 2009-11-24 2010-05-19 中国科学院长春光学精密机械与物理研究所 子孔径拼接干涉检测光学非球面面形的装置
CN101949691A (zh) * 2010-09-07 2011-01-19 中国科学院长春光学精密机械与物理研究所 非零位补偿浅度光学非球面面形检测方法
JP2011122857A (ja) * 2009-12-08 2011-06-23 Fujifilm Corp 非球面体測定方法および装置
CN102997863A (zh) * 2012-11-05 2013-03-27 北京理工大学 一种全口径光学非球面面形误差直接检测***
CN103471521A (zh) * 2012-10-29 2013-12-25 南通大学 快速、准确的光学非球面的实时检测方法
CN104236856A (zh) * 2014-09-10 2014-12-24 中国科学院上海光学精密机械研究所 物镜成像***的波像差检测装置及其***误差校正方法
CN105444693A (zh) * 2015-11-25 2016-03-30 中国科学院长春光学精密机械与物理研究所 一种浅度非球面的面形误差测量方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1995943A (zh) * 2007-01-04 2007-07-11 四川大学 一种大口径非球面镜全场检测方法
CN101709955A (zh) * 2009-11-24 2010-05-19 中国科学院长春光学精密机械与物理研究所 子孔径拼接干涉检测光学非球面面形的装置
JP2011122857A (ja) * 2009-12-08 2011-06-23 Fujifilm Corp 非球面体測定方法および装置
CN101949691A (zh) * 2010-09-07 2011-01-19 中国科学院长春光学精密机械与物理研究所 非零位补偿浅度光学非球面面形检测方法
CN103471521A (zh) * 2012-10-29 2013-12-25 南通大学 快速、准确的光学非球面的实时检测方法
CN102997863A (zh) * 2012-11-05 2013-03-27 北京理工大学 一种全口径光学非球面面形误差直接检测***
CN104236856A (zh) * 2014-09-10 2014-12-24 中国科学院上海光学精密机械研究所 物镜成像***的波像差检测装置及其***误差校正方法
CN105444693A (zh) * 2015-11-25 2016-03-30 中国科学院长春光学精密机械与物理研究所 一种浅度非球面的面形误差测量方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
CHRIS J. EVANS等: "Visualization of surface figure by the use of Zernike polynomials", 《APPLIED OPTICS》 *
PAUL E. MURPHY等: "Interference imaging for aspheric surface testing", 《APPLIED OPTICS》 *
李庆国等: "一种快速检测非球面面形的新方法", 《光学仪器》 *
王孝坤等: "干涉法实时测量浅度非球面技术", 《光学精密工程》 *
莫卫东等: "确定非球面检测中最佳参考球面波的方法", 《中国激光》 *
袁吕军等: "非球面数字波面检测技术", 《应用光学》 *
贺俊等: "使用红外干涉仪测量非球面面形", 《光学精密工程》 *

Also Published As

Publication number Publication date
CN110220477B (zh) 2020-09-29

Similar Documents

Publication Publication Date Title
CN105627947B (zh) 一种旋转对称未知非球面面形误差的测量方法及其测量装置
CN101858735B (zh) 一种大口径离轴非球面测量和标定***
EP2369319A2 (en) Aspheric object measuring method and apparatus
CN110531531B (zh) 卡塞格林光学***主次反射镜的装调方法
CN107782254A (zh) 一种混合补偿式子孔径拼接面形检测方法
CN107990833A (zh) 光学镀膜透镜中心厚度测量装置及方法
CN210773918U (zh) 非接触式透镜中心厚度测量装置
CN114252023B (zh) 用于非球面计算全息检测的计算机辅助调整装置及方法
CN110220477A (zh) 红外干涉仪测量光学非球面反射镜面形的检测装置及方法
CN104048596B (zh) 补偿器及补偿器与干涉仪共轴的调节方法
CN104634275A (zh) 一种基于牛顿环的非球面实时干涉测量装置及方法
IL309531A (en) Optically based verification of parallelism between internal planes
Herffurth et al. Assessing surface imperfections of freeforms using a robotic light scattering sensor
IL309527A (en) Methods and systems for verifying parallelism between internal planes
JP4768904B2 (ja) 光学素子又は光学系の物理量測定方法
Gebhardt et al. Alignment turning and assembly of the Sentinel 4 optical modules
CN110686618B (zh) 结合全反射角定位的非球面参数误差干涉测量方法及***
CN114137736B (zh) 大口径光学元件装调装置及其装调方法
Riediger et al. OCT measurement of aspheric polymer lenses for adaptive assembly of micro optical imaging objectives
Wang et al. Convex aspherical surface testing using catadioptric partial compensating system
CN108507488A (zh) 基于轴向扫描的锥镜面形检测***及检测方法
JP5620289B2 (ja) 被検面形状測定方法および被検面形状測定装置と被検面形状測定プログラム
Peng et al. Optimization of sub-cells orientation for assembly of a high-quality transmission sphere
Greenleaf Self-calibrating surface measuring machine
KR102028699B1 (ko) 렌즈의 3차원 형상 측정 방법 및 이를 위한 시스템

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant