CN103017681B - 近抛物面的旋转轴对称凹非球面的实时检测方法 - Google Patents
近抛物面的旋转轴对称凹非球面的实时检测方法 Download PDFInfo
- Publication number
- CN103017681B CN103017681B CN 201210521650 CN201210521650A CN103017681B CN 103017681 B CN103017681 B CN 103017681B CN 201210521650 CN201210521650 CN 201210521650 CN 201210521650 A CN201210521650 A CN 201210521650A CN 103017681 B CN103017681 B CN 103017681B
- Authority
- CN
- China
- Prior art keywords
- wave aberration
- wave
- aberration
- actual
- theoretical
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Testing Of Optical Devices Or Fibers (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种近抛物面的旋转轴对称凹非球面的实时检测方法,用平面自准直法仿真出凹非球面相对于最接近抛物面的波像差,在极坐标下利用zernike多项式进行拟合,将极坐标下的zernike方程转化为直角坐标下的形式;使用数字波面干涉仪测量出非球面相对于抛物面的波像差,将实际波像差的矩阵和理论波像差的矩阵统一到同一坐标系下,让两个波像差的像素一一对应,然后将两个波像差的矢高做差法运算,即,即可得到非球面实际面形与理论面形的残差分布。本发明具有快速、准确、检测范围广等优点,具有广阔的市场前景。
Description
技术领域
本发明属于先进光学制造和检测技术领域。
技术背景
近抛物面的旋转轴对称凹非球面主要是指接近抛物面的凹椭球面或凹双曲面以及以凹椭球面或凹双曲面为基础的高次曲面。高精度光学非球面元件面形的检测主要采用干涉检测技术。在该项技术中,无像差点检测、零位补偿干涉检测技术广泛应用于非球面抛光阶段的面形检测。
所谓的无像差点检测是指根据费马原理,光线从一点传到另外一点,经过任意多次折射或反射,其光程为极大值或极小值,也就是说光程是定值,光学上把这样的点成为无像差点,利用无像差点检测非球面的方法称为无像差检测。
此类无像差点检测方法具有一定的缺点,具体表现在无像差点检测主要用于检测旋转轴对称二次曲面,不能检测旋转轴对称高次曲面。
零位补偿干涉检测技术是指利用光学设计软件,如ZEMAX,CODE V等,设计一种带有特定波像差的光学系统,称之为零位补偿器,其中的零位补偿器的设计是基于理想非球面的,检验光束经由数字波面干涉仪出射至补偿器,光束经过补偿器再经被检非球面反射,再次经过补偿器后回到干涉仪,从而实现待检非球面元件面形的检测。
此类零位补偿检测不但能够检测旋转轴对称二次非球面也可以检测旋转轴对称高次非球面。但是这种检测方法也有一定的缺点,具体表现在针对不同面形的非球面元件,需要设计不同的补偿器,同时为了获得高精度的测量结果,要求在设计补偿器时,一方面使之能够很好地校正非球面波前差,另一方面要求补偿器各元件的厚度,曲率半径,空气间隔、同心度等公差分配合理。这样补偿器的误差极易产生鬼像,而导致衍射环的出现,并由于补偿器其中某些元件的反射光与参考光发生相互干涉,从而在像面上出现一些伪干涉条纹,由于这些伪干涉条纹与检测光同时发生相位移动,因此对检测结果影响很大。补偿器的精度不但受设计结果的影响,还会受装调的影响,补偿器自身精度的检测也是个难题。补偿检测光路调整复杂,耗时。
发明内容
本发明的目的在于提供一种方便、准确的近抛物面的旋转轴对称凹非球面的实时检测方法。
本发明的技术解决方案是:
一种近抛物面的旋转轴对称凹非球面的实时检测方法,其特征是:利用光学设计软件,如ZEMAX,CODE V等,用平面自准直法仿真出凹非球面(二次曲面或者高次曲面)相对于最接近抛物面的波像差,称为非球面相对于抛物面的理论波像差,将此波像差,在极坐标下利用zernike多项式(取前36项或前37项均可以)进行拟合,令x=rcosθ,y=rsinθ,将极坐标下的zernike方程转化为直角坐标下的形式;使用数字波面干涉仪(如zygo,wyko,fisba,esdi等)利用球面镜头搭建平面自准直光路测量出非球面相对于抛物面的波像差,称为非球面相对于抛物面的实际波像差。这个实际波像差用离散的三维矩阵(x,y,z)表示,x,y表示像素的位置,z表示对应像素位置波像差的矢高。根据实际波像差的三维矩阵,确定实际波面的有效像素,以此为依据对理论波像差在直角坐标系下的zernike多项式进行像素划分,将zernike多项式表示的理论波像差转化为矩阵(x',y',z')的形式,保证与实际波像差有效像素的分布相同,将实际波像差的矩阵和理论波像差的矩阵统一到同一坐标系下,让两个波像差的像素一一对应,然后将两个波像差的矢高做差法运算,即Δz=z′-z,即可得到非球面实际面形与理论面形的残差分布,从而实现对近抛物面的旋转轴对称凹非球面的实时检测。
标准平面镜口径不小于待测非球面镜口径,标准平面镜中心开中心孔,中心孔的大小不大于凹非球面中心盲区的大小。
数字波面干涉仪能够用自准直法测量出待测光学非球面相对于最接近抛物面的全口径波像差,且需要去除平移(piston),倾斜(tilt)等误差。
本发明不仅克服了无像差点不能检测旋转轴对称高次非球面,也克服了传统零位补偿检验中补偿器专用性,装调复杂、耗时等缺点,只需要大于等于待检非球面镜口径的标准平面镜,具有快速、准确、检测范围广等优点,具有广阔的市场前景。本检测方法检测非球面的最大非球面度和非球面度梯度取决于数字波面干涉仪内CCD阵列像元的大小和数目。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是本发明的工作原理图。
图2是本发明双曲面相对于最接近抛物面的理论波像差图。
图3是本发明双曲面检测光路图。
图4是本发明双曲面相对于最接近抛物面的实际波像差图。
图5是本发明实际双曲面面形与理论面形的残差分布图。
具体实施方式
一种近抛物面的旋转轴对称凹非球面的实时检测方法,其特征是:利用光学设计软件,如ZEMAX,CODE V等,用平面自准直法仿真出凹非球面(二次曲面或者高次曲面)相对于最接近抛物面的波像差,称为非球面相对于抛物面的理论波像差,将此波像差,在极坐标下利用zernike多项式(取前36项或37项均可以)进行拟合,令x=rcosθ,y=rsinθ,将极坐标下的zernike方程转化为直角坐标下的形式;使用数字波面干涉仪(如zygo,wyko,fisba,esdi等)利用球面镜头搭建平面自准直光路测量出非球面相对于抛物面的波像差,称为非球面相对于抛物面的实际波像差。这个实际波像差用离散的三维矩阵(x,y,z)表示,x,y表示像素的位置,z表示对应像素位置波像差的矢高。根据实际波像差的三维矩阵,确定实际波面的有效像素,以此为依据对理论波像差在直角坐标系下的zernike多项式进行像素划分,将zernike多项式表示的理论波像差转化为矩阵(x',y',z')的形式,保证与实际波像差有效像素的分布相同,将实际波像差的矩阵和理论波像差的矩阵统一到同一坐标系下,让两个波像差的像素一一对应,然后将两个波像差的矢高做差法运算,即Δz=z′-z,即可得到非球面实际面形与理论面形的残差分布,从而实现对近抛物面的旋转轴对称凹非球面的实时检测。
标准平面镜口径不小于待测非球面镜口径,标准平面镜中心开中心孔,中心孔的大小不大于凹非球面中心盲区的大小。
数字波面干涉仪能够用自准直法测量出待测光学非球面相对于最接近抛物面的全口径波像差,且需要去除平移(piston),倾斜(tilt)等误差。
利用光学设计软件,如ZEMAX,CODE V等,对近抛物面的旋转轴对称凹非球面进行平面自准直检测进行仿真,检测光路如图2所示,得到非球面相对于最接抛物面的理论波像差如图3所示。
利用数字波面干涉仪1检测待检非球面,2为干涉仪镜头,3为待检非球面镜,4为标准平面镜。标准平面镜置于干涉仪镜头与待检非球面镜之间,到干涉仪镜头焦点的距离小于rd/(2D)的范围均可,其中r为待检非球面镜的顶点曲率半径,d为标准平面镜中心孔的直径,D为待检非球面镜口径,为了调整的方便通常将标准平面镜置于干涉仪镜头焦点前,如图2所示。从而得到待测非球面3相对于最接近抛物面的实际波像差,去除平移(piston),倾斜(tilt)等误差,如图4所示。
将实际波像差的矩阵和理论波像差的矩阵统一到同一坐标系下,让两个波像差的像素一一对应,然后将两个波像差的矢高做差法运算,即Δz=z′-z,即可得到非球面实际面形与理论面形的残差分布,如图5所示。
Claims (1)
1.一种近抛物面的旋转轴对称凹非球面的实时检测方法,其特征是:用平面自准直法仿真出凹非球面相对于最接近抛物面的波像差,即非球面相对于抛物面的理论波像差,将此波像差,在极坐标下利用zernike多项式,取前36项或前37项进行拟合,令x=rcosθ,y=rsinθ,将极坐标下的zernike方程转化为直角坐标下的形式;使用数字波面干涉仪利用球面镜头搭建平面自准直光路测量出非球面相对于抛物面的波像差,即非球面相对于抛物面的实际波像差,这个实际波像差用离散的三维矩阵(x,y,z)表示,x,y表示像素的位置,z表示对应像素位置波像差的矢高;根据实际波像差的三维矩阵,确定实际波面的有效像素,以此为依据对理论波像差在直角坐标系下的zernike多项式进行像素划分,将zernike多项式表示的理论波像差转化为矩阵(x',y',z')的形式,保证与实际波像差有效像素的分布相同,将实际波像差的矩阵和理论波像差的矩阵统一到同一坐标系下,让两个波像差的像素一一对应,然后将两个波像差的矢高做差法运算,即Δz=z'-z,即可得到非球面实际面形与理论面形的残差分布,从而实现对近抛物面的旋转轴对称凹非球面的实时检测;所述使用数字波面干涉仪利用球面镜头搭建平面自准直光路是将标准平面镜置于干涉仪球面镜头与待检非球面镜之间,标准平面镜到干涉仪镜头焦点的距离小于rd/(2D),其中r为待检非球面镜的顶点曲率半径,d为标准平面镜中心孔的直径,D为待检非球面镜口径。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310413646.8A CN103471522B (zh) | 2012-12-07 | 2012-12-07 | 检测范围广的凹非球面的实时检测方法 |
CN 201210521650 CN103017681B (zh) | 2012-12-07 | 2012-12-07 | 近抛物面的旋转轴对称凹非球面的实时检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210521650 CN103017681B (zh) | 2012-12-07 | 2012-12-07 | 近抛物面的旋转轴对称凹非球面的实时检测方法 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310413646.8A Division CN103471522B (zh) | 2012-12-07 | 2012-12-07 | 检测范围广的凹非球面的实时检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103017681A CN103017681A (zh) | 2013-04-03 |
CN103017681B true CN103017681B (zh) | 2013-12-25 |
Family
ID=47966560
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310413646.8A Expired - Fee Related CN103471522B (zh) | 2012-12-07 | 2012-12-07 | 检测范围广的凹非球面的实时检测方法 |
CN 201210521650 Expired - Fee Related CN103017681B (zh) | 2012-12-07 | 2012-12-07 | 近抛物面的旋转轴对称凹非球面的实时检测方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310413646.8A Expired - Fee Related CN103471522B (zh) | 2012-12-07 | 2012-12-07 | 检测范围广的凹非球面的实时检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN103471522B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103234480A (zh) * | 2013-04-16 | 2013-08-07 | 北京理工大学 | 一种环形凸非球面的快速面形检测方法 |
CN103196391A (zh) * | 2013-04-16 | 2013-07-10 | 北京理工大学 | 一种近抛物面的环形凹非球面的快速面形检测方法 |
CN103791854A (zh) * | 2014-01-23 | 2014-05-14 | 中国科学院长春光学精密机械与物理研究所 | 高空间分辨率子孔径拼接方法 |
DE102015119274B4 (de) * | 2015-11-09 | 2018-07-12 | Björn Habrich | Verfahren und Vorrichtung zur Bestimmung der räumlichen Position eines Gegenstandes mittels interferometrischer Längenmessung |
CN106907991B (zh) * | 2017-02-24 | 2019-06-07 | 湖北航天技术研究院总体设计所 | 一种基于补偿器的离轴非球面镜零位检测对准方法 |
CN108227186B (zh) * | 2018-01-16 | 2020-03-20 | 南通大学 | 确定光学系统环形口径二次曲面最接近比较球面曲率半径的方法 |
CN110703434B (zh) * | 2019-10-15 | 2024-04-12 | 南通大学 | 一种确定环形孔径二次曲面非球面度梯度的方法 |
CN112697398B (zh) * | 2020-12-10 | 2023-09-19 | 中国科学院光电技术研究所 | 一种空间位置变化前后两次检测波像差残差的计算方法 |
CN117091532A (zh) * | 2023-08-25 | 2023-11-21 | 同济大学 | 一种非球面高精度干涉仪绝对测量装置及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1731232A (zh) * | 2005-09-05 | 2006-02-08 | 长春理工大学 | 一种光学非球面检测准万能补偿镜 |
US7276717B2 (en) * | 2004-11-05 | 2007-10-02 | Canon Kabushiki Kaisha | Measuring apparatus, exposure apparatus and device manufacturing method |
CN101650157A (zh) * | 2009-09-18 | 2010-02-17 | 中国科学院长春光学精密机械与物理研究所 | 双曲面凸面反射镜面形误差的检测方法及其装置 |
CN102155926A (zh) * | 2011-03-09 | 2011-08-17 | 浙江大学 | 一种非球面顶点球曲率半径的测量系统及方法 |
CN102175150A (zh) * | 2011-01-27 | 2011-09-07 | 南京理工大学 | 具有对点、测试双探测器的红外干涉检测装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008256517A (ja) * | 2007-04-04 | 2008-10-23 | Hoya Corp | 収差測定方法 |
-
2012
- 2012-12-07 CN CN201310413646.8A patent/CN103471522B/zh not_active Expired - Fee Related
- 2012-12-07 CN CN 201210521650 patent/CN103017681B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7276717B2 (en) * | 2004-11-05 | 2007-10-02 | Canon Kabushiki Kaisha | Measuring apparatus, exposure apparatus and device manufacturing method |
CN1731232A (zh) * | 2005-09-05 | 2006-02-08 | 长春理工大学 | 一种光学非球面检测准万能补偿镜 |
CN101650157A (zh) * | 2009-09-18 | 2010-02-17 | 中国科学院长春光学精密机械与物理研究所 | 双曲面凸面反射镜面形误差的检测方法及其装置 |
CN102175150A (zh) * | 2011-01-27 | 2011-09-07 | 南京理工大学 | 具有对点、测试双探测器的红外干涉检测装置 |
CN102155926A (zh) * | 2011-03-09 | 2011-08-17 | 浙江大学 | 一种非球面顶点球曲率半径的测量系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103471522B (zh) | 2016-03-02 |
CN103017681A (zh) | 2013-04-03 |
CN103471522A (zh) | 2013-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103017681B (zh) | 近抛物面的旋转轴对称凹非球面的实时检测方法 | |
CN102937421B (zh) | 旋转轴对称光学非球面的实时检测方法 | |
CN101876540B (zh) | 基于多波前透镜补偿器的非球面绝对测量系统 | |
CN101240999B (zh) | 一种大型光学平面的干涉测量装置与方法 | |
CN100507488C (zh) | 非球面镜衍射图像检测方法与装置 | |
CN101650157B (zh) | 双曲面凸面反射镜面形误差的检测方法及其装置 | |
CN102997864B (zh) | 一种大口径光学非球面镜检测系统 | |
CN108871733B (zh) | 大口径光学系统近场检测装置及其测量方法 | |
CN101285732B (zh) | 一种大口径抛物面镜检测系统 | |
CN103234480A (zh) | 一种环形凸非球面的快速面形检测方法 | |
CN104048619B (zh) | 一种判断旋转轴对称非球面能否采用直接干涉检测的方法 | |
CN103557791B (zh) | 一种大口径非球面主镜二次常数测量装置与方法 | |
CN101241000B (zh) | 高陡度凸二次非球面的无像差点法子孔径拼接测量方法 | |
CN104006759A (zh) | 大口径大偏离量非球面反射镜抛光过程中复合检测方法 | |
CN110082073B (zh) | 一种在子孔径拼接检测光学系统透射波前中调整平面反射镜倾斜的装置和方法 | |
CN108061514B (zh) | 一种利用轴向扫描光干涉法检测非球面的动态建模方法 | |
CN100567932C (zh) | 扇形离轴非球面镜拼接测量系统 | |
CN103196391A (zh) | 一种近抛物面的环形凹非球面的快速面形检测方法 | |
CN105352451B (zh) | 一种基于可变形镜的准万能补偿镜及设计方法 | |
CN107131846A (zh) | 一种用于超大口径凸扁球面检测的光学系统 | |
CN102519611A (zh) | 共光路径向剪切数字波面干涉仪 | |
CN104034262B (zh) | 一种非球面镜二次常数测量方法 | |
CN206803957U (zh) | 用于超大口径凸扁球面检测的光学系统 | |
CN104776811B (zh) | 大口径、批量化离轴子镜的面形检测方法 | |
CN109579739A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131225 Termination date: 20151207 |
|
EXPY | Termination of patent right or utility model |