CN102519388B - 大口径离轴凸非球面反射镜面形检测方法 - Google Patents

大口径离轴凸非球面反射镜面形检测方法 Download PDF

Info

Publication number
CN102519388B
CN102519388B CN2011103348520A CN201110334852A CN102519388B CN 102519388 B CN102519388 B CN 102519388B CN 2011103348520 A CN2011103348520 A CN 2011103348520A CN 201110334852 A CN201110334852 A CN 201110334852A CN 102519388 B CN102519388 B CN 102519388B
Authority
CN
China
Prior art keywords
aperture
aspheric surface
convex aspheric
sub
measured
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.)
Active
Application number
CN2011103348520A
Other languages
English (en)
Other versions
CN102519388A (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.)
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Original Assignee
Changchun Institute of Optics Fine Mechanics and Physics of CAS
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 Changchun Institute of Optics Fine Mechanics and Physics of CAS filed Critical Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority to CN2011103348520A priority Critical patent/CN102519388B/zh
Publication of CN102519388A publication Critical patent/CN102519388A/zh
Application granted granted Critical
Publication of CN102519388B publication Critical patent/CN102519388B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Instruments For Measurement Of Length By Optical Means (AREA)

Abstract

本发明涉及一种大口径离轴凸非球面反射镜面形检测方法,该方法包括如下步骤:利用与待测非球面反射镜基准子孔径区域面形相应的计算全息元件对基准子孔径区域进行零位补偿干涉测量;利用与待测非球面反射镜离轴子孔径区域面形相应的计算全息元件对离轴子孔径区域进行零位补偿干涉测量:求解待测凸非球面反射镜全口径面形分布;将各离轴子孔径在xyz坐标系内的坐标换算为XYZ坐标系下的坐标,通过计算得到待测凸非球面反射镜全孔径相位分布;从全口径相位分布中剔去系统调整误差获得待测凸非球面反射镜全口径面形误差分布。本发明数据处理和数学运算简单,实验操作简单易行,检测成本很低,测试时间短。

Description

大口径离轴凸非球面反射镜面形检测方法
技术领域
本发明涉及一种检测光学非球面面形的方法,特别涉及一种检测大口径离轴凸非球面反射镜面形的方法。
背景技术
在光学系统中使用非球面元件,能矫正像差,改善像质,而且可以减小光学系统的尺寸和重量,因此非球面光学元件正越来越多的被用于各个领域。在众多应用领域中,凸非球面元件是较为常用的核心部件。
但是凸非球面面形的测量一直是光学检测中的难点,以往我们一般采用无像差点法对其进行测量,其实质是:若待测表面是理想的,而点光源精确置于其中的一个几何焦点上,则由表面反射的光线形成球面波前,其球心与另一几何焦点重合。通过求解波前的变形量及其方向,并据此确定引起波前形变的非球面形状偏差。无像差点法是一种经典、可靠的检验二次曲面的方法,其优点是设计理论误差为零;但缺点是当检验大口径非球面时,就需要更大口径的Hindle球或者大口径的辅助平面,尤其是检测凸非球面,辅助镜面的口径往往是待测镜面的几倍,这就造成了材料和工艺上的难度。此外利用该方法在检测反射系统时,非球面的中心部位往往无法检测。
目前,设计补偿透镜利用零位补偿法对凸非球面进行测量,仍是检测凸非球面镜最常用的方法之一。但是凸非球面补偿器的设计要比凹非球面补偿器的设计要复杂得多,凸非球面的补偿器至少要有两片透镜组成,有时补偿器本身可能还含有非球面元件,而要实现该非球面的高精度检测还要为其再设计一套补偿器,因此又给加工和装调带来了很多的困难。此外,通过计算全息(CGH)进行零位补偿可以较好的实现对非球面的检测,但是对于高陡度、大偏离量的非球面,计算全息的线频密度非常高,从而使其刻划制作成为不可能。此外,由于检测凸非球面需要汇聚波前作为参考波面,这就需要所使用的干涉仪的口径和补偿透镜或CGH的口径大于待测凸非球面的口径,大口径的补偿透镜和CGH的制作和装调在目前仍然存在很多的困难。
由此可得,要想完成对大口径离轴凸非球面的检测,最常规的方法之一是制作零位补偿透镜,但是补偿器将由几块透镜组成,结构复杂,且补偿透镜本身可能还含有非球面,此外补偿器的口径将非常大,因此该方法存在很多瓶颈;另一种常规的方法是制作衍射元件计算全息完成对其测量,但是对于大偏离量大口径的凸非球面,需制作很大口径且线频密度很高的计算全息元件,目前计算全息制作受口径和刻划线频密度的限制,从而使这种方法很难实现。此外,利用子孔径拼接法测定多个子孔径的相位分布,通过子孔径拼接算法可以完成对浅度凸非球面的拼接测量,但是对于大口径、大偏离量的凸非球面的拼接检测,子孔径数目很多,运算非常复杂,且不可避免的带来拼接误差传递和累积。
发明内容
本发明要解决的技术问题是提供一种将计算全息补偿与子孔径拼接技术联用,补偿元件结构简单、子孔径数目少、数据处理和数学运算简单、能够实现大口径离轴凸非球面面形精确测量的大口径离轴凸非球面反射镜面形检测方法。
为了解决上述技术问题,本发明的大口径离轴凸非球面反射镜面形检测方法包括如下步骤:
1)将干涉仪和待测凸非球面反射镜分别安装在第一和第三调整机构上;设干涉仪的光轴为Z轴,待测凸非球面反射镜的顶点为坐标系XYZ的原点,则过待测凸非球面反射镜的顶点且垂直于干涉仪光轴的面为XOY面;
2)将待测凸非球面反射镜中心区域的子孔径作为基准子孔径;利用与待测凸非球面反射镜基准子孔径区域面形相应的计算全息元件对基准子孔径区域进行零位补偿干涉测量:
将计算全息元件安装在第二调整机构上;利用第一、第二、第三调整机构调整干涉仪、计算全息元件和待测凸非球面反射镜之间的位置关系,使得干涉仪出射的光束经计算全息元件零位补偿后能够沿法向入射到待测凸非球面反射镜的基准子孔径区域,该光束经待测凸非球面反射镜反射后与干涉仪的参考波面形成干涉条纹,得到待测凸非球面反射镜基准子孔径的相位分布;
3)利用与待测凸非球面反射镜离轴子孔径区域面形相应的计算全息元件对离轴子孔径区域进行零位补偿干涉测量:
依次将与离轴子孔径区域面形相应的计算全息元件安装在第二调整机构上,并利用第一、第二、第三调整机构调整待测凸非球面反射镜、计算全息元件与干涉仪的相对位置关系,使得干涉仪的出射光束经各计算全息元件零位补偿后能够沿法向入射到待测凸非球面反射镜上对应的离轴子孔径区域,该光束经待测凸非球面反射镜反射后与干涉仪的参考波面形成干涉条纹,得到待测凸非球面反射镜的各离轴子孔径区域相位分布;
4)求解待测凸非球面反射镜全口径面形分布:
选择测量各子孔径相位分布时各子孔径的顶点o为子孔径坐标系xyz的原点,过各子孔径顶点o且垂直于干涉仪光轴的平面为xoy面;
利用式(1)求出各子孔径重叠区域相位差的平方和值S并令其为最小值;
S = Σ j 1 ≠ 0 N 1 Σ i 1 ⋐ W 0 , W j 1 n { W 0 ( x 0 - j 1 i 1 , y 0 - j 1 i 1 ) - [ W j 1 ( x j 1 - j 1 i 1 , y j 1 - j 1 i 1 ) + p j 1 x j 1 - j 1 i 1
+ a j 1 x j 1 - j 1 i 1 + b j 1 y j 1 - j 1 i 1 ] } 2 + Σ j 2 ∩ j 3 ≠ 0 N 2 Σ i 2 ⋐ W j 2 , W j 3 n { [ W j 2 ( x j 2 - j 2 i 2 , y j 2 - j 2 i 2 )
+ p j 2 x j 2 - j 2 i 2 + a j 2 x j 2 - j 2 i 2 + b j 2 y j 2 - j 2 i 2 ]
- [ W j 3 ( x j 3 - j 3 i 2 , y j 3 - j 3 i 2 ) + p j 3 x j 3 - j 3 i 2 + a j 3 x j 3 - j 3 i 2 + b j 3 y j 3 - j 3 i 2 ] } 2
(1)
式(1)中,N1是离轴子孔径与基准子孔径的重叠区域数,N2是离轴子孔径间的重叠区域数,n是重叠区域内的采样点数;
Figure GDA00003368248000045
为测得的基准子孔径相位分布中基准子孔径与第j1个离轴子孔径的重叠区域中i1点的相位,为测得的第j1个离轴子孔径相位分布中基准子孔径与第j1个离轴子孔径的重叠区域中i1点的相位;
Figure GDA00003368248000047
Figure GDA00003368248000048
分别是第j1个离轴子孔径相对基准子孔径沿X方向的倾斜系数、沿Y方向的倾斜系数和Z方向相对平移系数;
Figure GDA000033682480000410
为测得的第j2个离轴子孔径相位分布中第j2个离轴子孔径与第j3个离轴子孔径的重叠区域中i2点的相位,
Figure GDA000033682480000411
为测得的第j3个离轴子孔径相位分布中第j2个离轴子孔径与第j3个离轴子孔径的重叠区域中i2点的相位;
Figure GDA000033682480000413
Figure GDA000033682480000414
分别是第j2个离轴子孔径相对基准子孔径沿X方向的倾斜系数、沿Y方向的倾斜系数和Z方向相对平移系数;
Figure GDA000033682480000416
Figure GDA000033682480000417
分别是第j3个离轴子孔径相对基准子孔径1沿X方向的倾斜系数、沿Y方向的倾斜系数和Z方向相对平移系数;
利用最小二乘拟合,对各个系数分别求偏导并令其值为零,得到离轴子孔径相对基准子孔径的最佳拼接因子aj、bj和pj
∂ S ∂ p j = 0 ∂ S ∂ a j = 0 ∂ S ∂ b j = 0 - - - ( 2 )
式(2)中,1≤j≤M-1,M为待测凸非球面反射镜子孔径的数量;
将各离轴子孔径在xyz坐标系内的坐标换算为XYZ坐标系下的坐标,利用式(3)得到拼接后的各离轴子孔径的相位分布wj(X,Y);并结合公式:w0(X,Y)=w0(x,y)获得待测凸非球面反射镜全口径相位分布W(X,Y);
wj(X,Y)=w0(X,Y)-(pj+ajX+bjY)    (3)
5)从全口径相位分布W(X,Y)中剔去系统调整误差ε(X,Y)获得待测凸非球面反射镜全口径面形误差分布e(X,Y):
e(X,Y)=W(X,Y)-ε(X,Y)    (4)
ε(X,Y)=AX+BY+C(X2+Y2)+D     (5)
A B C D = ΣXX ΣXY ΣX ( X 2 + Y 2 ) ΣX ΣYX ΣYY ΣY ( X 2 + Y 2 ) ΣY Σ ( X 2 + Y 2 ) X Σ ( X 2 + Y 2 ) Y Σ ( X 2 + Y 2 ) 2 Σ ( X 2 + Y 2 ) ΣX ΣY Σ ( X 2 + Y 2 ) N - 1 ΣXW ΣYW Σ ( X 2 + Y 2 ) W W - - - ( 6 )
(12)其中,A、B、C、D分别为Z方向平移、X方向倾斜、Y方向倾斜和离焦项,N为求解系数A、B、C、D所需的数据点数目,N大于等于4;W=W(X,Y);
6)误差与精度分析:
判断e(X,Y)的均方根值与待测凸非球面反射镜面形精度要求的设计值之间的相对偏差是否大于20%;若大于20%,则重复步骤1)~5),对待测凸非球面反射镜的面形进行检测,直至得到的e(X,Y)的均方根值与待测凸非球面反射镜面形精度要求的设计值之间的相对偏差小于20%,此时得到的e(X,Y)即为待测凸非球面反射镜的面形误差分布。
本发明对检测结果进行精度分析和判别以其获得准确的待测凸非球面反射镜面形误差分布,检测精度可达λ/50以上,λ=632.8nm。
本发明根据待测凸非球面反射镜的数学方程划分子孔径区域,设计与各子孔径对应的单块计算全息元件,利用单块计算全息元件补偿干涉仪出射的光束,测得各子孔径的相位分布,然后通过子孔径拼接就能够准确的实现对非球面反射镜面形的测量。本发明降低了补偿元件的设计难度、简化了补偿元件的结构、避开了大口径计算全息元件不可制作和大偏离量非球面计算全息线频密度太高不可刻划的难点,有效的减少了子孔径拼接数目。本发明物理概念明确,数据处理和数学运算简单,实验操作简单易行,检测成本很低,测试时间短。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1是实现本发明大口径离轴凸非球面反射镜面形检测方法的装置结构示意图。
图2是本发明的大口径离轴凸非球面反射镜面形检测方法流程图。
图3是本发明的测量大口径离轴凸非球面反射镜的子孔径分布示意图。
图4a、4b分别为基准子孔径和离轴子孔径在测量相位分布时的位置示意图,图4c为子孔径侧视图。
图5a、5b分别为两个有重叠区域的离轴子孔径在测量相位分布时的位置示意图。
具体实施方式
如图1所示,实现本发明的大口径离轴凸非球面反射镜面形检测方法的装置包括激光干涉仪11、计算全息元件12、第一调整机构13、第二调整机构14、第三调整机构15和待测大口径离轴凸非球面16。
所述第一调整机构13采用目前公知的三维精密调整架,可以精确控制干涉仪11沿X方向、Y方向和沿Z方向(光轴方向)的平动。第二调整机构14采用目前公知的六维精密调整架,可以精确调整计算全息元件12的俯仰、扭摆和旋转以及沿X方向、Y方向和沿Z方向(光轴方向)的平动。第三调整机构15采用目前公知的六维精密调整架,可以精确调整被测凸非球面16的俯仰、扭摆和旋转以及沿X方向、Y方向和沿Z方向(光轴方向)的平动。
如图2所示,本发明的大口径离轴凸非球面反射镜面形检测方法流程如下:
1)首先,根据待测凸非球面反射镜的数学方程划分子孔径区域并设计与各子孔径对应的单块计算全息元件12;所述计算全息元件为振幅型(CGH)。
子孔径划分和单块计算全息元件12设计的原则是:一是各子孔径间有重叠区域;二是各子孔径区域所对应的计算全息元件的线频密度不太高,在目前的工艺刻划水平能够实现。
2)根据光学设计结果,加工和制作计算全息元件12;
首先,根据设计结果选购计算全息元件基板,然后对玻璃基板进行加工和检测,并将最终检测结果代入光学设计中进行复算,复算通过后,我们设计了刻蚀设备接口文件,将其以激光直写或电子束直写方式制作光刻掩膜版,用光刻方法将图形刻划到预先制备并通过高精度检测的玻璃基板上,形成振幅型CGH,最后进行光学镀膜以实现增加光强的目的。
3)利用计算全息元件对各子孔径区域进行零位补偿干涉测量;
子孔径划分如图3所示;首先在第二调整机构14上安装好与待测凸非球面反射镜中心区域(基准子孔径)对应的计算全息元件12,利用第一、第二、第三调整机构13、14、15调整好干涉仪11、计算全息元件和待测凸非球面反射镜16之间的位置关系,使得干涉仪11出射的光束经计算全息元件零位补偿后能够沿法向入射到待测凸非球面反射镜16上,并精确的位于待测凸非球面反射镜16的中心区域,该光束经待测凸非球面反射镜16反射后与干涉仪11的参考波面形成干涉条纹,从而可以测定得到待测凸非球面反射镜16中心区域(基准子孔径)的相位分布。
在第二调整机构14上依次更换与离轴子孔径对应的计算全息元件12,并利用第一、第二、第三调整机构13、14、15调整待测凸非球面反射镜16、计算全息元件12与干涉仪11的相对位置关系,使得干涉仪11的出射光束经各计算全息元件12零位补偿后能够沿法向精确的入射到待测凸非球面反射镜16各对应离轴子孔径区域,该光束经待测凸非球面反射镜16反射后与干涉仪11的参考波面形成干涉条纹,从而可以测定得到待测凸非球面反射镜的各离轴子孔径区域相位分布。
通过全局优化子孔径拼接算法分析求解待测凸非球面反射镜全口径面形分布。
为了精确获得全口径面形误差分布,本发明将采用误差均化的全局优化拼接算法,具体分析和求解过程如下:
通过重叠区域的数据求解相邻两个子孔径的相对调整误差,可以实现子孔径两两拼接。多次利用两个子孔径的拼接原理就可以实现多个子孔径的拼接。但是这样往往会造成误差传递和累积,从而降低了整个孔径的检测精度。因此,在子孔径拼接过程中存在综合优化的问题。当拼接区域大于两个时,假设共有M个子孔径拼接,可以先选定其中任一个子孔径作为基准,为了便于定位和测量,一般选择待测凸非球面反射镜中心区域的子孔径(基准子孔径)作为参考标准。
由于是对大口径非球面反射镜各区域进行零位补偿测量,因此各子孔径间位置的相对失调量仅带来相对平移和倾斜。
设干涉仪的光轴为Z轴,待测凸非球面反射镜的顶点为坐标系XYZ的原点,则过待测凸非球面反射镜的顶点且垂直于干涉仪光轴的面为XOY面;选择测量各子孔径相位分布时各子孔径的顶点o为子孔径坐标系xyz的原点,过各子孔径顶点o且垂直于干涉仪光轴的平面为xoy面。
如图3所示,假使基准子孔径1的相位分布为w0(X,Y),则拼接后离轴子孔径相位分布wj(X,Y)与基准子孔径1的相位分布关系为:
wj(X,Y)=w0(X,Y)-(pj+ajX+bjY)    (3)
w0(X,Y)=w0(x,y)
(X,Y)是待测凸非球面反射镜内任意一点在坐标系XYZ内的坐标,M为子孔径的数量(包括基准子孔径和其它子孔径);aj、bj和pj分别是第j个离轴子孔径相对基准子孔径沿X方向的倾斜系数、沿Y方向的倾斜系数和Z方向相对平移系数。利用最小二乘法,使得所有重叠区域相位差的平方和值S为最小,可得下式:
S = Σ j 1 ≠ 0 N 1 Σ i 1 ⋐ W 0 , W j 1 n { W 0 ( x 0 - j 1 i 1 , y 0 - j 1 i 1 ) - [ W j 1 ( x j 1 - j 1 i 1 , y j 1 - j 1 i 1 ) + p j 1 x j 1 - j 1 i 1
+ a j 1 x j 1 - j 1 i 1 + b j 1 y j 1 - j 1 i 1 ] } 2 + Σ j 2 ∩ j 3 ≠ 0 N 2 Σ i 2 ⋐ W j 2 , W j 3 n { [ W j 2 ( x j 2 - j 2 i 2 , y j 2 - j 2 i 2 )
+ p j 2 x j 2 - j 2 i 2 + a j 2 x j 2 - j 2 i 2 + b j 2 y j 2 - j 2 i 2 ]
- [ W j 3 ( x j 3 - j 3 i 2 , y j 3 - j 3 i 2 ) + p j 3 x j 3 - j 3 i 2 + a j 3 x j 3 - j 3 i 2 + b j 3 y j 3 - j 3 i 2 ] } 2 = min
(1)
式中N1是离轴子孔径与基准子孔径的重叠区域数,N2是离轴子孔径间的重叠区域数,n是重叠区域内的采样点数;如图4a、4b、4c所示,
Figure GDA00003368248000095
为测得的基准子孔径1相位分布中基准子孔径1与第j1个离轴子孔径2的重叠区域3中第i1点的相位,
Figure GDA00003368248000096
为测得的第j1个离轴子孔径2相位分布中基准子孔径1与第j1个离轴子孔径2的重叠区域3中第i1点的相位;
Figure GDA00003368248000101
Figure GDA00003368248000102
Figure GDA00003368248000103
分别是第j1个离轴子孔径2相对基准子孔径1沿X方向的倾斜系数、沿Y方向的倾斜系数和Z方向相对平移系数;如图5a、5b所示,
Figure GDA00003368248000104
为测得的第j2个离轴子孔径4相位分布中第j2个离轴子孔径4与第j3个离轴子孔径5的重叠区域6中第i2点的相位,为测得的第j3个离轴子孔径5相位分布中第j2个离轴子孔径4与第j3个离轴子孔径5的重叠区域6中第i2点的相位;
Figure GDA00003368248000107
Figure GDA00003368248000108
分别是第j2个离轴子孔径5相对基准子孔径1沿X方向的倾斜系数、沿Y方向的倾斜系数和Z方向相对平移系数;
Figure GDA00003368248000109
Figure GDA000033682480001010
Figure GDA000033682480001011
分别是第j3个离轴子孔径5相对基准子孔径1沿X方向的倾斜系数、沿Y方向的倾斜系数和Z方向相对平移系数。
利用最小二乘拟合,对各个系数分别求偏导并令其值为零可得:
∂ S ∂ p j = 0 ∂ S ∂ a j = 0 ∂ S ∂ b j = 0 - - - ( 2 )
式中1≤j≤M-1,M为待测凸非球面反射镜子孔径的数量。利用(1)、(2)式就可以得到各离轴子孔径相对基准子孔径的最佳拼接因子ai、bi和pi(aj、bj和pj);然后,利用式(3)得到拼接后的各离轴子孔径的相位分布wj(x,y),再将各离轴子孔径在xyz坐标系内的坐标换算为待测凸非球面反射镜坐标系XYZ下的坐标,从而把所有离轴子孔径的位相数据校正统一到相同的基准上。
以上是一种全局拼接的方法,可以将拼接累积误差在整个拼接区域“扩散”开,起到了误差均化的作用,从而可以大大地减小拼接累积误差,很好地实现多个子孔径的拼接。
实验中一般以待测元件中心区域的子孔径为基准,通过以上的综合优化拼接可以获得待测凸非球面反射镜的整个面形的相位分布,但是所得的相位分布中必将存在一定的系统调整和定位误差,必须将其消去。假定系统高阶误差可以忽略,则经过综合优化拼接获得的全孔径波前误差的相位测量值可表示为:
W(X,Y)=e(X,Y)+ε(X,Y)     (7)
式中:(X,Y)是待测凸非球面反射镜上任意一点在XYZ坐标系内的坐标;e(X,Y)表示(X,Y)点面形误差;ε(X,Y)为(X,Y)点系统调整误差。W(X,Y)是全孔径拼接得到的(X,Y)点相位。通常整个系统的装调定位误差为四项,分别对应为Z方向平移、X方向倾斜、Y方向倾斜和离焦项,即:
ε(X,Y)=AX+BY+C(X2+Y2)+D     (5)
利用最小二乘拟合,可求得系统调整误差系数的最优解:
A B C D = ΣXX ΣXY ΣX ( X 2 + Y 2 ) ΣX ΣYX ΣYY ΣY ( X 2 + Y 2 ) ΣY Σ ( X 2 + Y 2 ) X Σ ( X 2 + Y 2 ) Y Σ ( X 2 + Y 2 ) 2 Σ ( X 2 + Y 2 ) ΣX ΣY Σ ( X 2 + Y 2 ) N - 1 ΣXW ΣYW Σ ( X 2 + Y 2 ) W W - - - ( 6 )
(12)式(6)中,N为求解系数A、B、C、D所需的数据点数目,N大于等于4;W=W(X,Y)。
因此,从拼接所得的全口径相位分布中剔去系统调整误差即可获得准确的待测凸非球面反射镜全口径面形误差分布,从而完成对大口径待测凸非球面反射镜面形的高精度检测。
5)误差与精度分析;
计算全息补偿与子孔径拼接联用检测大口径离轴凸非球面的方法受很多因素的影响:如计算全息元件制作精度、子孔径拼接算法精度、机构调整精度以及测试环境等影响。我们将对检测结果进行精度分析和判别,判断e(X,Y)的均方根值与待测凸非球面反射镜面形精度要求的设计值之间的相对偏差是否大于20%;若大于20%,则重复以上步骤,对待测凸非球面反射镜的面形进行检测,直至得到的e(X,Y)的均方根值与待测凸非球面反射镜面形精度要求的设计值之间的相对偏差小于20%,此时得到的e(X,Y)即为待测凸非球面反射镜的面形误差分布。

Claims (1)

1.一种大口径离轴凸非球面反射镜面形检测方法,其特征在于包括如下步骤:
1)将干涉仪(11)和待测凸非球面反射镜(16)分别安装在第一和第三调整机构(13、15)上;设干涉仪的光轴为Z轴,待测凸非球面反射镜的顶点为坐标系XYZ的原点,则过待测凸非球面反射镜的顶点且垂直于干涉仪光轴的面为XOY面;
2)将待测凸非球面反射镜中心区域的子孔径作为基准子孔径;利用与待测凸非球面反射镜基准子孔径区域面形相应的计算全息元件对基准子孔径区域进行零位补偿干涉测量:
将计算全息元件(12)安装在第二调整机构(14)上;利用第一、第二、第三调整机构(13、14、15)调整干涉仪(11)、计算全息元件(12)和待测凸非球面反射镜(16)之间的位置关系,使得干涉仪(11)出射的光束经计算全息元件(12)零位补偿后能够沿法向入射到待测凸非球面反射镜(16)的基准子孔径区域,该光束经待测凸非球面反射镜(16)反射后与干涉仪(11)的参考波面形成干涉条纹,得到待测凸非球面反射镜(16)基准子孔径的相位分布;
3)利用与待测凸非球面反射镜离轴子孔径区域面形相应的计算全息元件对离轴子孔径区域进行零位补偿干涉测量:
依次将与离轴子孔径区域面形相应的计算全息元件安装在第二调整机构(14)上,并利用第一、第二、第三调整机构(13、14、15)调整待测凸非球面反射镜(16)、计算全息元件与干涉仪(11)的相对位置关系,使得干涉仪(11)的出射光束经各计算全息元件零位补偿后能够沿法向入射到待测凸非球面反射镜(16)上对应的离轴子孔径区域,该光束经待测凸非球面反射镜(16)反射后与干涉仪(11)的参考波面形成干涉条纹,得到待测凸非球面反射镜的各离轴子孔径区域相位分布;
4)求解待测凸非球面反射镜全口径面形分布:
选择测量各子孔径相位分布时各子孔径的顶点o为子孔径坐标系xyz的原点,过各子孔径顶点o且垂直于干涉仪光轴的平面为xoy面;
利用式(1)求出各子孔径重叠区域相位差的平方和值S并令其为最小值;
S = Σ j 1 ≠ 0 N 1 Σ i 1 ⋐ W 0 , W j 1 n { W 0 ( x 0 - j 1 i 1 , y 0 - j 1 i 1 ) - [ W j 1 ( x j 1 - j 1 i 1 , y j 1 - j 1 i 1 ) + p j 1 x j 1 - j 1 i 1
+ a j 1 x j 1 - j 1 i 1 + b j 1 y j 1 - j 1 i 1 ] } 2 + Σ j 2 ∩ j 3 ≠ 0 N 2 Σ i 2 ⋐ W j 2 , W j 3 n { [ W j 2 ( x j 2 - j 2 i 2 , y j 2 - j 2 i 2 )
+ p j 2 x j 2 - j 2 i 2 + a j 2 x j 2 - j 2 i 2 + b j 2 y j 2 - j 2 i 2 ]
- [ W j 3 ( x j 3 - j 3 i 2 , y j 3 - j 3 i 2 ) + p j 3 x j 3 - j 3 i 2 + a j 3 x j 3 - j 3 i 2 + b j 3 y j 3 - j 3 i 2 ] } 2
(1)
式(1)中,N1是离轴子孔径与基准子孔径的重叠区域数,N2是离轴子孔径间的重叠区域数,n是重叠区域内的采样点数;为测得的基准子孔径(1)相位分布中基准子孔径(1)与第j1个离轴子孔径(2)的重叠区域(3)中i1点的相位,
Figure FDA00003368247900025
为测得的第j1个离轴子孔径(2)相位分布中基准子孔径(1)与第j1个离轴子孔径(2)的重叠区域(3)中i1点的相位;
Figure FDA00003368247900026
Figure FDA00003368247900027
Figure FDA00003368247900028
分别是第j1个离轴子孔径(2)相对基准子孔径(1)沿X方向的倾斜系数、沿Y方向的倾斜系数和Z方向相对平移系数;
Figure FDA00003368247900029
为测得的第j2个离轴子孔径(4)相位分布中第j2个离轴子孔径(4)与第j3个离轴子孔径(5)的重叠区域(6)中i2点的相位,
Figure FDA000033682479000210
为测得的第j3个离轴子孔径(5)相位分布中第j2个离轴子孔径(4)与第j3个离轴子孔径(5)的重叠区域(6)中i2点的相位;
Figure FDA000033682479000211
Figure FDA000033682479000212
分别是第j2个离轴子孔径(5)相对基准子孔径(1)沿X方向的倾斜系数、沿Y方向的倾斜系数和Z方向相对平移系数;
Figure FDA000033682479000215
Figure FDA00003368247900031
分别是第j3个离轴子孔径(5)相对基准子孔径(1)沿X方向的倾斜系数沿Y方向的倾斜系数和Z方向相对平移系数;
利用最小二乘拟合,对各个系数分别求偏导并令其值为零,得到离轴子孔径相对基准子孔径的最佳拼接因子aj、bj和pj
∂ S ∂ p j = 0 ∂ S ∂ a j = 0 ∂ S ∂ b j = 0 - - - ( 2 )
式(2)中,1≤j≤M-1,M为待测凸非球面反射镜子孔径的数量;
将各离轴子孔径在xyz坐标系内的坐标换算为XYZ坐标系下的坐标,利用式(3)得到拼接后的各离轴子孔径的相位分布wj(X,Y);并结合公式:w0(X,Y)=w0(x,y)获得待测凸非球面反射镜全口径相位分布W(X,Y);
wj(X,Y)=w0(X,Y)-(pj+ajX+bjY)    (3)
5)从全口径相位分布W(X,Y)中剔去系统调整误差ε(X,Y)获得待测凸非球面反射镜全口径面形误差分布e(X,Y):
e(X,Y)=W(X,Y)-ε(X,Y)    (4)
ε(X,Y)=AX+BY+C(X2+Y2)+D    (5)
A B C D = ΣXX ΣXY ΣX ( X 2 + Y 2 ) ΣX ΣYX ΣYY ΣY ( X 2 + Y 2 ) ΣY Σ ( X 2 + Y 2 ) X Σ ( X 2 + Y 2 ) Y Σ ( X 2 + Y 2 ) 2 Σ ( X 2 + Y 2 ) ΣX ΣY Σ ( X 2 + Y 2 ) N - 1 ΣXW ΣYW Σ ( X 2 + Y 2 ) W W - - - ( 6 )
(12)其中,A、B、C、D分别为Z方向平移、X方向倾斜、Y方向倾斜和离焦项,N为求解系数A、B、C、D所需的数据点数目,N大于等于4;W=W(X,Y);
6)误差与精度分析:
判断e(X,Y)的均方根值与待测凸非球面反射镜面形精度要求的设计值之间的相对偏差是否大于20%;若大于20%,则重复步骤1)~5),对待测凸非球面反射镜的面形进行检测,直至得到的e(X,Y)的均方根值与待测凸非球面反射镜面形精度要求的设计值之间的相对偏差小于20%,此时得到的e(X,Y)即为待测凸非球面反射镜的面形误差分布。
CN2011103348520A 2011-10-28 2011-10-28 大口径离轴凸非球面反射镜面形检测方法 Active CN102519388B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011103348520A CN102519388B (zh) 2011-10-28 2011-10-28 大口径离轴凸非球面反射镜面形检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011103348520A CN102519388B (zh) 2011-10-28 2011-10-28 大口径离轴凸非球面反射镜面形检测方法

Publications (2)

Publication Number Publication Date
CN102519388A CN102519388A (zh) 2012-06-27
CN102519388B true CN102519388B (zh) 2013-11-27

Family

ID=46290392

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011103348520A Active CN102519388B (zh) 2011-10-28 2011-10-28 大口径离轴凸非球面反射镜面形检测方法

Country Status (1)

Country Link
CN (1) CN102519388B (zh)

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103196389B (zh) * 2013-04-02 2015-09-30 中国科学院光电技术研究所 检测计算全息基片面形和材料不均匀性误差的装置和方法
CN103335607B (zh) * 2013-06-28 2015-12-02 中国科学院西安光学精密机械研究所 大口径金属离轴椭球镜或抛物镜面形检验方法
CN103386640B (zh) * 2013-07-29 2016-01-20 中国科学院长春光学精密机械与物理研究所 大口径反射镜加工机床精确对准方法
CN104006759A (zh) * 2014-05-05 2014-08-27 中国科学院长春光学精密机械与物理研究所 大口径大偏离量非球面反射镜抛光过程中复合检测方法
CN104142129B (zh) * 2014-07-18 2016-11-23 中国科学院长春光学精密机械与物理研究所 离轴三反非球面系统凸非球面次镜面形拼接检测方法
CN104776811B (zh) * 2015-04-10 2017-11-28 中国科学院国家天文台南京天文光学技术研究所 大口径、批量化离轴子镜的面形检测方法
CN106091978B (zh) * 2016-06-01 2019-02-19 西安工程大学 斜入式激光干涉法测量中干涉条纹图像的拼接方法
CN106225712A (zh) * 2016-08-01 2016-12-14 中国科学院长春光学精密机械与物理研究所 一种离轴三反非球面光学系统共基准检测与加工方法
CN106225713A (zh) * 2016-08-02 2016-12-14 中国科学院长春光学精密机械与物理研究所 一种离轴三反非球面光学系统共基准检测与加工方法
CN107782254B (zh) * 2017-09-29 2019-05-21 中国科学院长春光学精密机械与物理研究所 一种混合补偿式子孔径拼接面形检测方法
CN108267094B (zh) * 2018-01-12 2020-04-14 暨南大学 一种基于旋转cgh的非圆柱面干涉拼接测量系统及方法
CN109855560B (zh) * 2018-12-28 2020-05-19 中国科学院长春光学精密机械与物理研究所 一种凸非球面反射镜面形的检测装置及检测方法
CN110108228A (zh) * 2019-04-26 2019-08-09 中国科学院国家天文台南京天文光学技术研究所 计算全息和球面反射镜检测不同离轴量非球面的方法
CN110625629B (zh) * 2019-10-29 2023-05-02 中国科学院长春光学精密机械与物理研究所 一种大口径光学复杂曲面的立式加工检测集成方法
CN110966954A (zh) * 2019-10-31 2020-04-07 中国科学院长春光学精密机械与物理研究所 大口径光学元件面形拼接检测方法及设备
CN110986823A (zh) * 2019-12-19 2020-04-10 华中科技大学 一种大口径改正镜面形测量系统和方法
CN111536896B (zh) * 2020-04-09 2022-06-07 中国科学院长春光学精密机械与物理研究所 一种激光干涉面形检测自动检测装置与方法
CN111623957B (zh) * 2020-05-11 2022-03-25 中国科学院光电技术研究所 一种用于x射线聚焦镜拼接干涉检测的点云配准拼接方法
CN112461156A (zh) * 2020-11-12 2021-03-09 中国科学院长春光学精密机械与物理研究所 大口径凸非球面光学元件的面形精度检测装置及检测方法
CN112517940B (zh) * 2020-11-19 2022-08-09 中国人民解放军国防科技大学 一种超精密车削加工在位检测装置及其应用方法
CN113175897B (zh) * 2021-05-08 2023-03-14 中国科学院国家天文台南京天文光学技术研究所 离轴高次类椭球面反射镜的组合补偿面形检测系统及方法
CN114252244B (zh) * 2021-12-17 2022-09-20 华中科技大学 一种二次曲面反射镜拼接干涉检测子孔径规划方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101013027A (zh) * 2007-01-31 2007-08-08 中国人民解放军国防科学技术大学 大口径大相对孔径非球面镜中高频误差检测装置与方法
CN101709955A (zh) * 2009-11-24 2010-05-19 中国科学院长春光学精密机械与物理研究所 子孔径拼接干涉检测光学非球面面形的装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5245402A (en) * 1992-06-15 1993-09-14 The United States Of America As Represented By The Secretary Of The Army General aspherical surface optical testing device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101013027A (zh) * 2007-01-31 2007-08-08 中国人民解放军国防科学技术大学 大口径大相对孔径非球面镜中高频误差检测装置与方法
CN101709955A (zh) * 2009-11-24 2010-05-19 中国科学院长春光学精密机械与物理研究所 子孔径拼接干涉检测光学非球面面形的装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Measurement of large aspheric surfaces by annular subaperture stitching interferometry;王孝坤等;《Chinese Optics Letters》;20071110;第5卷(第11期);第645-647页 *
子孔径拼接干涉检测凸非球面的研究;王孝坤等;《光学学报》;20100731;第30卷(第7期);第2022-2026页 *
子孔径拼接干涉检测离轴非球面研究;王孝坤等;《光子学报》;20110131;第40卷(第1期);第92-97页 *
子孔径拼接干涉测量一大口径双曲面;王孝坤等;《红外与激光工程》;20090225;第38卷(第1期);第114-119页 *

Also Published As

Publication number Publication date
CN102519388A (zh) 2012-06-27

Similar Documents

Publication Publication Date Title
CN102519388B (zh) 大口径离轴凸非球面反射镜面形检测方法
CN101709955B (zh) 子孔径拼接干涉检测光学非球面面形的装置
WO2017101557A1 (zh) 面形检测装置及方法
CN101687302B (zh) 非球面、其幅值的布置及包含其的投射物镜
CN101949691A (zh) 非零位补偿浅度光学非球面面形检测方法
US7948638B2 (en) Scanning interferometric methods and apparatus for measuring aspheric surfaces and wavefronts
CN102506750A (zh) 部分补偿非球面反射镜面形检测方法
CN107782254B (zh) 一种混合补偿式子孔径拼接面形检测方法
CN101915556B (zh) 可用于低反射率光学球面面形检测的偏振点衍射干涉系统
CN101419062B (zh) 一种大口径双曲面次镜检测系统
CN104006759A (zh) 大口径大偏离量非球面反射镜抛光过程中复合检测方法
CN103969787A (zh) 一种离轴四反镜头的初装定位方法
CN106441816B (zh) 计算全息法测量长焦距透镜透射波前的检测装置及检测方法
CN103776389A (zh) 一种高精度非球面组合干涉检测装置与方法
CN102155926A (zh) 一种非球面顶点球曲率半径的测量系统及方法
CN103196389B (zh) 检测计算全息基片面形和材料不均匀性误差的装置和方法
CN102636130B (zh) 一种大动态范围测量非球面光学元件表面轮廓装置
CN104142129A (zh) 离轴三反非球面系统凸非球面次镜面形拼接检测方法
CN100585362C (zh) 一种大口径非球面镜全场检测方法
US20050225774A1 (en) Method for measuring and manufacturing an optical element and optical apparatus
CN105444693A (zh) 一种浅度非球面的面形误差测量方法
CN102927930B (zh) 采用平行光管拼接检测超大口径反射镜面形误差的方法
CN102393565B (zh) 一种折反式逆补偿器
CN103175481B (zh) 一种离轴光学非球面镜顶点半径的测量方法和装置
CN105627945A (zh) 非球面元件中心与外圆中心偏离量的测量装置及测量方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant