CN1885097A - 利用计算机对非球面检验光学系统进行调整的方法 - Google Patents
利用计算机对非球面检验光学系统进行调整的方法 Download PDFInfo
- Publication number
- CN1885097A CN1885097A CN 200510016907 CN200510016907A CN1885097A CN 1885097 A CN1885097 A CN 1885097A CN 200510016907 CN200510016907 CN 200510016907 CN 200510016907 A CN200510016907 A CN 200510016907A CN 1885097 A CN1885097 A CN 1885097A
- Authority
- CN
- China
- Prior art keywords
- aspheric mirror
- wave aberration
- tested aspheric
- sigma
- optical system
- 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.)
- Pending
Links
Images
Landscapes
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
本发明针对利用计算机对非球面检验光学系统进行调整的方法。首先完成系统光路的粗装调,然后实时采集一幅干涉图,获得被检非球面镜的波像差,利用赛德尔像差多项式作为基底拟合获得被检非球面镜失调引入波像差,根据拟合波像差与失调引入波像差的相互对应关系得到失调量的数值,利用计算机控制调整机构实现系统精确定位。再采集一幅干涉图,当干涉条纹为“零”条纹时,则调整结束;否则,对新采集的干涉图拟合计算失调量并且重新调整,直至调整结果满足要求。本发明结合了光学系统的检验、检验结果的数据分析处理以及优化算法、实现了微小平移量和转动量的控制与精确调整,使调整过程在理论指导下可视、定量、有序地完成,提高了干涉检验的效率。
Description
技术领域:本发明属于非球面检验技术领域,特别涉及利用计算机对非球面检验光学系统进行调整的方法。
背景技术:光学系统的计算机辅助装调技术是通过计算机对实时检测结果和理论结果进行对比分析,经过必要的优化处理,给出光学系统中各零件失调的方位和量值。在装调过程中,可以有目的地使各零件达到最佳的位置状态,提高装调的效率。从八十年代开始,国内外光学工作者(美国Santa Barbara研究中心、美国Arizona大学光学中心、美国Ball Aerospace技术公司、日本Nikon公司、中国科学院长春光机所、哈尔滨工业大学)就已经开始研究计算机辅助装调技术,但它主要是针对加工完成后的近衍射极限的成像光学系统进行装调。在非球面干涉检验过程中,由于传统的手工调整,实现光路的准确定位非常困难,因此如何高效、准确地实现非球面检验光学系统精确定位,提高干涉检验的效率,则是现有技术中有待解决的问题。
发明内容:为了解决上述背景技术非球面干涉检验过程中,采用传统的手工调整周期长、实现光路的准确定位非常困难,本发明的目的是从实时采集的干涉图中分离出非球面检验光学系统失调引入波像差,计算出调整机构的调整量数值,采用计算机控制调整调整机构实现非球面检验光学系统的精确定位。为此本发明将提供一种利用计算机对非球面检验光学系统进行调整的方法,提高干涉检验的效率,为加工出合格的非球面面形提供保障。
本发明的技术方案:
1.首先将非球面检验光学系统中的数字波面干涉仪、光学补偿器和被检非球面镜置于气浮防震平台上,调整数字波面干涉仪、光学补偿器和被检非球面镜之间的相对位置,直至在数字波面干涉仪观察到非球面检验光学系统出瞳波面的干涉条纹,则完成非球面检验光学系统光路的粗装调;
2.在步骤1的基础上,再利用数字波面干涉仪实时采集一幅干涉图,获得被检非球面镜的波像差为W(x,y),波像差W(x,y)包含被检非球面镜的实际加工面形残留波像差为Wm(x,y)以及光学补偿器和被检非球面镜的失调引入的波像差为Wa(x,y);
3.根据步骤2再利用刚体位移理论得到非球面检验系统中的光学补偿器和被检非球面镜失调引入的波像差Wa(x,y);
利用赛德尔(seidel)像差多项式作为基底拟合获得被检非球面镜失调引入波像差Wa(x,y),然后根据拟合波像差与失调引入波像差的相互对应关系得到失调量Dx、Dy、Dz以及θx、θy的数值;
4.按照实际要求,利用第一调整机构精密调整光学补偿器的位置,利用第二调整机构精密调整被检非球面镜的位置,第一调整机构实现被检非球面镜的平动调整(即X、Y、Z三个方向的平动),第二调整机构实现被检非球面镜的转动调整(即俯仰和扭摆);根据步骤3所获得的失调量数据,利用计算机控制第一调整机构和第二调整机构调整光学补偿器和被检非球面镜位置的精确定位;
5.根据步骤4采集一幅干涉图,当数字波面干涉仪所获得的干涉条纹为“零”条纹时,则非球面检验光学系统中的光学补偿器和被检非球面镜调整结束;否则,返回步骤3,直至调整光学补偿器和被检非球面镜的结果满足要求,则计算机完成对非球面检验光学系统的调整装配。
优点:由于本发明在非球面检验光学系统调整过程中,从实时采集的干涉图中,分离了系统失调引入波像差,结合了光学系统的检验、检验结果的数据分析处理以及优化算法、实现了微小平移量和转动量的控制与精确调整。从干涉图中分离非球面检验光学系统失调引入波像差,给出了同轴非球面分离失调引入波像差的数学模型,对于离轴非球面则将其看作同轴母镜的一部分,获取各采样点在同轴母镜下的坐标值,在同轴条件下对离轴非球面检验光学系统进行拟合计算,获得失调量的数值结果,采用计算机辅调整技术,使整个非球面检验光学系统的调整过程在理论指导下可视、定量、有序地完成,从而提高干涉检验的效率。实践证明,该方案具有很高的调整精度,收敛速度很快,实验表明:经过四次调整后,系统像差的均方根值(RMS值)为0.028λ(λ=632.8nm),缩短了干涉检验周期。
附图说明:
图1为本发明的流程图。
具体实施方式:
本发明的技术方案的流程图如图1所示。
非球面检验光学系统由Zygo干涉仪、光学补偿器和被检非球面镜三部分组成。由零位补偿检验原理可知,在零位补偿检验结果中,既包含有被检非球面的实际加工残留的波像差,也包含光路调整过程中失调量引入的波像差。因此,要获得被检非球面镜镜面的实际加工残留的波像差数据,必须从干涉检验结果中分离系统失调引入的波像差,对非球面检验光学系统进行精确调整。
本发明可以直接用于同轴非球面检验光学系统,也可以用于离轴非球面检验光学系统。
当用于离轴非球面检验光学系统时,将被检非球面镜看作是同轴母镜的一部分,根据被检离轴非球面镜设计要求,获得离轴非球面和其同轴母镜之间的关系,将被检离轴非球面镜上各点的坐标转换到同轴母镜坐标系中,得到被检离轴非球面镜上各采样点在同轴母镜下的坐标,使用上述同轴非球面的拟合条件对离轴非球面的失调量进行拟合;
失调量的确定过程是数据处理过程也是优化过程。算法上即是建立失调量和失调量引入的波像差之间的关系。波像差W(x,y)包括被检非球面镜镜面的实际加工残留误差为Wm(x,y)和由第一调整机构和第二调整机构的失调引入的误差为Wa(x,y)(其中(x,y)代表镜面上某一点的位置)。因此干涉图波像差W(x,y)可用下式表示:
W(x,y)=Wm(x,y)+Wa(x,y) (1)
所以要获得合格的非球面面形,必须从测量结果中去除失调引入波像差Wa(x,y)。将被检非球面镜看作是一个六自由度的刚体。分别为沿xyz轴的平动Dx,Dy,Dz和绕x,y,z轴的转动θx,θy,θz,其中由于z为光轴,所以绕z轴的旋转θz因其具有回转对称性可以忽略。根据刚体位移理论我们可以获得被检非球面镜引入波像差:
Wa(x,y)=-2Dz+2x(cDx+θy)+2y(cDy-θx)+p2c2Dz
在上述式(2)中:p2=x2+y2(x,y为镜面坐标),c=1/r,r为被检非球面镜的曲率半径,k为被检非球面镜的二次曲面常数;
-2Dz:平移; 2x(cDx+θy):X-倾斜;
2y(cDy-θx):Y-倾斜; p2c2Dz:离焦;
xp2c3kDx:X-彗差; yp2c3kDy:Y-彗差; 球差
考虑到离轴非球面可以看作同轴母镜的一部分,因此(2)式也适用于对离轴非球面的调整量误差拟合,在获得采样点在同轴母镜下的坐标后,可以使用同轴非球面的拟合条件对离轴非球面的失调量误差进行拟合。
采用赛德尔(seidel)像差多项式作为基底拟合获得被检非球面镜失调引入波像差Wa(x,y),设:
Wa(x,y)=c0+c1x+c2y+c3(x2+y2)+c4x(x2+y2)+c5y(x2+y2)+c6(x4+y4)(3)
采用最小二乘法,即F(x,y)=[w(x,y)-wa(x,y)]2=min,求解系数c0、c1、c2、c3、c4、c5、c6的最小二乘解。
令:
由于F(x,y)=min的充要条件为:
由(4)式,可以获得一个包含7个未知量的线性方程组:
设:
则:C=A-1B
在求得C的最小二乘解以后,然后根据拟合波像差与失调引入波像差的相互对应关系得到失调量Dx、Dy、Dz以及θx、θy的数值。在实践中应根据获得的近似解重新调整光路,然后再求解、再调整,直到光学调整达到满意的精度为止。在光学调整精确的情况下,测得的结果必然是一个不包含失调量误差的面形误差分布,使用这种方法可以较为精确地实现失调量误差与加工误差的分离。
按照实际要求,设计两套调整机构:将光学补偿器与第一调整机构固定连接,用于精密调整光学补偿器的位置;将被检非球面镜第二调整机构固定连接,用于精密调整被检非球面镜的位置,并且确保调整机构具有合适的行程范围。第一调整机构主要实现被检非球面镜的平动调整(即X、Y、Z三个方向的平动),第二调整机构主要实现被检非球面镜的转动调整(即俯仰和扭摆);将所获得的失调量数据存储于计算机的存储器中,利用计算机控制第一调整机构和第二调整机构调整光学补偿器和被检非球面镜的精确位置;
第一调整机构和第二调整机构的控制系统由计算机、运动控制卡、驱动器及步进电机构成,利用步进电机对调整机构进行控制。
根据计算获得的失调量数据,采用计算机控制调整机构,精确控制与调整光学补偿器和被检非球面镜的微小平移量和转动量,实现非球面检验光学系统精确定位。
我们结合实际加工任务,利用计算机对一离轴非球面检验光学系统进行光路调整,检验系统粗装调后,利用Zygo数字波面干涉仪实时采集一幅干涉图,波像差为1.070λ(PV值)、0.164λ(RMS值)(λ=632.8nm),经过四次调整后系统像差的RMS值为0.028λ。据此看出,该方法收敛速度很快,具有很高的精度,缩短了系统检验周期。
Claims (1)
1、利用计算机对非球面检验光学系统进行调整的方法,其特征在于调整步骤如下:
(1)首先将非球面检验光学系统中的数字波面干涉仪、光学补偿器和被检非球面镜置于气浮防震平台上,调整数字波面干涉仪、光学补偿器和被检非球面镜之间的相对位置,直至在数字波面干涉仪观察到非球面检验光学系统的出瞳波面的干涉条纹,则完成非球面检验光学系统光路的粗装调;
(2)在步骤1的基础上,再利用数字波面干涉仪实时采集一幅干涉图,获得被检非球面镜的波像差为W(x,y),波像差W(x,y)包含被检非球面镜的实际加工面形残留波像差为Wm(x,y)以及光学补偿器和被检非球面镜的失调引入的波像差为Wa(x,y);
(3)根据步骤2再利用刚体位移理论得到非球面检验系统中的光学补偿器和被检非球面镜的失调引入的波像差Wa(x,y);
利用赛德尔(seidel)像差多项式作为基底拟合获得被检非球面镜失调引入波像差Wa(x,y),然后根据拟合波像差与失调引入波像差的相互对应关系得到失调量Dx、Dy、Dz以及θx、θy的数值;
(4)按照实际要求,利用第一调整机构精密调整光学补偿器的位置,利用第二调整机构精密调整被检非球面镜的位置,第一调整机构实现被检非球面镜的平动调整(即X、Y、Z三个方向的平动),第二调整机构实现被检非球面镜的转动调整(即俯仰和扭摆);根据步骤3所获得的失调量数据,利用计算机控制第一调整机构和第二调整机构调整光学补偿器和被检非球面镜的位置精确定位;
(5)根据步骤4采集一幅干涉图,当数字波面干涉仪所获得的干涉条纹为“零”条纹时,则非球面检验光学系统中的光学补偿器和被检非球面镜调整结束;否则,返回步骤3,直至调整光学补偿器和被检非球面镜的结果满足要求,则计算机完成对非球面检验光学系统的调整装配。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200510016907 CN1885097A (zh) | 2005-06-22 | 2005-06-22 | 利用计算机对非球面检验光学系统进行调整的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200510016907 CN1885097A (zh) | 2005-06-22 | 2005-06-22 | 利用计算机对非球面检验光学系统进行调整的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1885097A true CN1885097A (zh) | 2006-12-27 |
Family
ID=37583344
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200510016907 Pending CN1885097A (zh) | 2005-06-22 | 2005-06-22 | 利用计算机对非球面检验光学系统进行调整的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1885097A (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100432749C (zh) * | 2007-07-25 | 2008-11-12 | 中国科学院上海技术物理研究所 | 离轴三反射光学系统装调中光轴变换方法 |
CN102322817A (zh) * | 2011-06-13 | 2012-01-18 | 中国科学院长春光学精密机械与物理研究所 | 一种使用菲索干涉仪检测被测元件的快速自动调整的方法 |
CN102368114A (zh) * | 2011-11-15 | 2012-03-07 | 中国科学院西安光学精密机械研究所 | 基于波像差检测的光学系统面形补偿装调方法 |
CN102478386A (zh) * | 2010-11-30 | 2012-05-30 | 苏州大学 | 三次位相板精度的检测方法及系统 |
CN103322942A (zh) * | 2013-07-15 | 2013-09-25 | 中国科学院光电技术研究所 | 基于波像差的光学系统各光学件面形的检测方法 |
CN103471521A (zh) * | 2012-10-29 | 2013-12-25 | 南通大学 | 快速、准确的光学非球面的实时检测方法 |
CN105547183A (zh) * | 2015-12-21 | 2016-05-04 | 中国科学院长春光学精密机械与物理研究所 | 一种复位被测非球面空间位置的调整方法 |
CN106225713A (zh) * | 2016-08-02 | 2016-12-14 | 中国科学院长春光学精密机械与物理研究所 | 一种离轴三反非球面光学系统共基准检测与加工方法 |
CN106907991A (zh) * | 2017-02-24 | 2017-06-30 | 湖北航天技术研究院总体设计所 | 一种基于补偿器的离轴非球面镜零位检测对准方法 |
CN107271150A (zh) * | 2017-06-27 | 2017-10-20 | 南京理工大学 | 一种基于干涉测量及灵敏度矩阵的计算机辅助装调方法 |
CN109839086A (zh) * | 2017-11-27 | 2019-06-04 | 长光华大基因测序设备(长春)有限公司 | 一种变形镜刚体位移误差测量的方法及装置 |
CN110634372A (zh) * | 2019-09-29 | 2019-12-31 | 中国科学院长春光学精密机械与物理研究所 | 一种光学系统装调策略验证系统 |
CN112697398A (zh) * | 2020-12-10 | 2021-04-23 | 中国科学院光电技术研究所 | 一种空间位置变化前后两次检测波像差残差的计算方法 |
CN115586640A (zh) * | 2022-10-21 | 2023-01-10 | 中国人民解放军国防科技大学 | 一种同轴两反系统分步计算的装调方法及装置 |
-
2005
- 2005-06-22 CN CN 200510016907 patent/CN1885097A/zh active Pending
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100432749C (zh) * | 2007-07-25 | 2008-11-12 | 中国科学院上海技术物理研究所 | 离轴三反射光学系统装调中光轴变换方法 |
CN102478386A (zh) * | 2010-11-30 | 2012-05-30 | 苏州大学 | 三次位相板精度的检测方法及系统 |
CN102478386B (zh) * | 2010-11-30 | 2013-07-03 | 苏州大学 | 三次位相板精度的检测方法及系统 |
CN102322817A (zh) * | 2011-06-13 | 2012-01-18 | 中国科学院长春光学精密机械与物理研究所 | 一种使用菲索干涉仪检测被测元件的快速自动调整的方法 |
CN102368114A (zh) * | 2011-11-15 | 2012-03-07 | 中国科学院西安光学精密机械研究所 | 基于波像差检测的光学系统面形补偿装调方法 |
CN103471521A (zh) * | 2012-10-29 | 2013-12-25 | 南通大学 | 快速、准确的光学非球面的实时检测方法 |
CN103471521B (zh) * | 2012-10-29 | 2015-10-21 | 南通大学 | 快速、准确的光学非球面的实时检测方法 |
CN103322942A (zh) * | 2013-07-15 | 2013-09-25 | 中国科学院光电技术研究所 | 基于波像差的光学系统各光学件面形的检测方法 |
CN105547183B (zh) * | 2015-12-21 | 2017-11-28 | 中国科学院长春光学精密机械与物理研究所 | 一种复位被测非球面空间位置的调整方法 |
CN105547183A (zh) * | 2015-12-21 | 2016-05-04 | 中国科学院长春光学精密机械与物理研究所 | 一种复位被测非球面空间位置的调整方法 |
CN106225713A (zh) * | 2016-08-02 | 2016-12-14 | 中国科学院长春光学精密机械与物理研究所 | 一种离轴三反非球面光学系统共基准检测与加工方法 |
CN106907991A (zh) * | 2017-02-24 | 2017-06-30 | 湖北航天技术研究院总体设计所 | 一种基于补偿器的离轴非球面镜零位检测对准方法 |
CN107271150A (zh) * | 2017-06-27 | 2017-10-20 | 南京理工大学 | 一种基于干涉测量及灵敏度矩阵的计算机辅助装调方法 |
CN107271150B (zh) * | 2017-06-27 | 2019-05-07 | 南京理工大学 | 一种基于干涉测量及灵敏度矩阵的计算机辅助装调方法 |
CN109839086A (zh) * | 2017-11-27 | 2019-06-04 | 长光华大基因测序设备(长春)有限公司 | 一种变形镜刚体位移误差测量的方法及装置 |
CN109839086B (zh) * | 2017-11-27 | 2021-01-15 | 长春长光华大智造测序设备有限公司 | 一种变形镜刚体位移误差测量的方法及装置 |
CN110634372A (zh) * | 2019-09-29 | 2019-12-31 | 中国科学院长春光学精密机械与物理研究所 | 一种光学系统装调策略验证系统 |
CN112697398A (zh) * | 2020-12-10 | 2021-04-23 | 中国科学院光电技术研究所 | 一种空间位置变化前后两次检测波像差残差的计算方法 |
CN112697398B (zh) * | 2020-12-10 | 2023-09-19 | 中国科学院光电技术研究所 | 一种空间位置变化前后两次检测波像差残差的计算方法 |
CN115586640A (zh) * | 2022-10-21 | 2023-01-10 | 中国人民解放军国防科技大学 | 一种同轴两反系统分步计算的装调方法及装置 |
CN115586640B (zh) * | 2022-10-21 | 2024-07-09 | 中国人民解放军国防科技大学 | 一种同轴两反系统分步计算的装调方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1885097A (zh) | 利用计算机对非球面检验光学系统进行调整的方法 | |
CN102661719B (zh) | 用于非球面子孔径拼接测量的近零位补偿器及面形测量仪和测量方法 | |
WO2021203707A1 (zh) | 一种激光干涉面形检测自动检测装置与方法 | |
Fang et al. | Manufacturing and measurement of freeform optics | |
CN100462673C (zh) | 大口径大相对孔径非球面镜中高频误差检测装置与方法 | |
CN102997863B (zh) | 一种全口径光学非球面面形误差直接检测系统 | |
CN106441153B (zh) | 一种大口径非球面元件轮廓高精度检测方法及装置 | |
CN105423948A (zh) | 采用变形镜的拼接干涉检测非球面面形的装置及方法 | |
CN101949691A (zh) | 非零位补偿浅度光学非球面面形检测方法 | |
CN1746648A (zh) | 一种大口径深型非球面镜检测系统 | |
CN102620683B (zh) | 子孔径拼接检测非球面调整误差补偿方法 | |
CN104006759A (zh) | 大口径大偏离量非球面反射镜抛光过程中复合检测方法 | |
CN102095385A (zh) | 新型球面绝对测量系统及方法 | |
CN112504177B (zh) | 多功能立式零位重叠扫描干涉测量装置 | |
CN106918301A (zh) | 平面面形子孔径拼接干涉测量装置和测量方法 | |
CN102717305A (zh) | 一种光学自由曲面原位测量方法 | |
CN106595473A (zh) | 非球面模具的在位测量系统及其测量方法和测量检验方法 | |
Semenov et al. | M1 and M2 mirror manufacturing for ARIES project: current status | |
CN114252023A (zh) | 用于非球面计算全息检测的计算机辅助调整装置及方法 | |
West et al. | Alignment and use of the optical test for the 8.4-m off-axis primary mirrors of the Giant Magellan Telescope | |
CN108592820A (zh) | 基于动态波前调制结合计算全息片的干涉面形检测方法 | |
Liu et al. | Performance-controllable manufacture of optical surfaces by ultra-precision machining | |
CN112504176B (zh) | 光学面形局部大误差的补偿干涉测量装置与方法 | |
CN110887638B (zh) | 一种光学系统像平面绘制装置与方法 | |
CN110109230A (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 | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |