CN107764203A - 基于部分补偿法的双波长相移干涉非球面测量方法及装置 - Google Patents

基于部分补偿法的双波长相移干涉非球面测量方法及装置 Download PDF

Info

Publication number
CN107764203A
CN107764203A CN201711019962.1A CN201711019962A CN107764203A CN 107764203 A CN107764203 A CN 107764203A CN 201711019962 A CN201711019962 A CN 201711019962A CN 107764203 A CN107764203 A CN 107764203A
Authority
CN
China
Prior art keywords
wavefront
tested
wavelength
dual wavelength
aspherical
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
CN201711019962.1A
Other languages
English (en)
Other versions
CN107764203B (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201711019962.1A priority Critical patent/CN107764203B/zh
Publication of CN107764203A publication Critical patent/CN107764203A/zh
Application granted granted Critical
Publication of CN107764203B publication Critical patent/CN107764203B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

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

Landscapes

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

Abstract

本发明属于光学精密测试技术领域,涉及一种基于部分补偿法的双波长相移干涉方法及实现装置。该方法构建部分补偿法双波长相移干涉仪,获取两个单波长的被测波前包裹相位;建模部分补偿法双波长理想干涉仪,获得两个单波长的剩余波前和包裹相位;采用误差分离和消除算法消除被测波前包裹相位中的已知及未知波前变化量,最后采用逆向迭代优化重构被测非球面的面形误差。装置包括第一激光器和第二激光器、第一狭缝和第二狭缝、第一平面镜和第二平面镜、第一分束镜和第二分束镜、扩束镜、准直镜、标准平面镜、移相器、部分补偿镜、被测非球面、成像镜头和含稀疏阵列传感器的干涉图采集组件。本发明特别适用于小面形误差的陡度非球面,以及大面形误差的模压非球面和自由曲面的加工质量测量。

Description

基于部分补偿法的双波长相移干涉非球面测量方法及装置
技术领域
本发明属于光学精密测试技术领域,特别涉及双波长相移干涉检测方法的改进,用于实现大非球面度或大面形误差的光学非球面面形误差的高精度测量。
背景技术
利用非球面取代传统球面可用更少的光学元件改善系统成像性能,同时具有增加设计自由度、减小体积和减轻系统重量等优点。在航天相机和天文望远镜等精密光学系统中,某些非球面特别是陡度光学非球面的制造和测试甚至对光学系统的更新换代起决定性的作用。然而非球面波前特别是陡度非球面波前的高精度检测一直是非球面光学元件设计和制造中面临的测量难题。
如果被测波前与参考球面存在较大的偏离,产生的密集干涉条纹不满足奈奎斯特条件,即每条干涉条纹少于两个像素,从而超过了相移干涉测量技术(PSI)可检测到的波前斜率最大极限。目前,常用的扩展PSI量程的非球面检测方法有环形或圆形子孔径拼接法,长波长法,高密度探测器法,双波长相移干涉法(TWPSI)和亚奈奎斯特法(SNI)。
子孔径拼接技术是将全口径分割成许多个满足奈奎斯特条件的小区域,以求解每一个子孔径,子孔径拼接法需要精密运动的机械扫描装置和复杂的拼接算法去实现子孔径的自动切换和拼接,检测时间也相对较长。长波长法和高密度探测器法分别采用昂贵的长波长光源和高密度探测器扩大PSI可测动态范围,大大增加了仪器系统的成本。为了测量高陡度波前,亚利桑那大学光学研究中心先后提出了双波长相移干涉法(TWPSI)和亚奈奎斯特法(SNI)。SNI是基于被测波前一阶导数或斜率的连续性假设重建被测波前,采用稀疏阵列传感器记录被测波前相位,可实现相位变化超过π的高斜率非球面波前检测。TWPSI不需要对被测波前的预知,TWPSI采用等效波长的解包裹波前作为参考修正单波长2π模糊性解包裹波前,可实现相位变化超过π的陡度非球面波前检测。
专利号为US4832489的专利首先提出双波长相移干涉装置和方法,用于精确重构陡度非球面的轮廓,采用了泰曼-格林干涉光路结构,运用等效合成波长解包裹波前修正单波长解包裹波前的2π模糊性。对与参考表面偏离达数百个波长的被测表面实现了单波长精度测量,但是需要针对双波长干涉仪以及被测非球面特别设计消色差零位补偿镜,增加了系统设计难度和实现的通用性。
专利号为CN201310571673.8和CN201410342492.2的专利先后提出了一种双波长同时相移干涉测量方法,对采集到的双波长混合相移干涉条纹图,分别提取两个单波长的包裹相位,然后计算出合成波长相位。该方法主要适用于阶梯、凹槽等具有突变的微结构表面的三维形貌测量。
专利号为CN201620294905.9的专利提出了一种斐索式双波长干涉测试装置及其合成波长相位提取方法,需要针对部分光路进行消色差设计,主要解决合成波长相位提取速度和不同波长移相误差对合成波长的影响问题。
由于具体应用领域不同,上述专利各自实现了相关领域的关键问题解决,因此有各自的优缺点和适用范围。对于TWPSI应用于陡度非球面测量,补偿镜和干涉仪的色差一直是TWPSI实际应用需要解决的关键问题,如果对补偿镜和干涉仪进行消色差设计必然大大增加仪器设计难度和周期。
部分补偿法是近年来逐渐发展起来的一种非球面非零位检测方法,它一般仅采用单透镜部分补偿被测非球面的法线像差以减少被测非球面斜率,降低了补偿镜的设计难度,减少了仪器系统设计周期,是一种实现非球面通用化检测的有效技术途径。针对目前非球面特别是陡度非球面的测量需求,本发明提出一种将部分补偿法和TWPSI相结合的方法,用于大非球面度或大面形误差的非球面测量,该方法中采用算法消除了干涉仪和补偿镜的色差及其误差放大效应,因此不需要特别针对干涉仪和补偿镜进行消色差设计;同时针对由于部分补偿法与TWPSI结合违背了零位条件存在的在回程误差,提出了校正方法,可实现被测非球面面形误差的高精度测量。
发明内容
为了突破双波长相移干涉用于陡度波前测量的限制,本发明的目的在于提供一种基于部分补偿法的双波长相移干涉方法及实现装置,特别适用于小面形误差的陡度非球面,以及大面形误差的模压非球面和自由曲面的加工质量测量。
为了实现大非球面度或大面形误差非球面面形误差的高精度测量,本发明提出一种基于部分补偿法的双波长相移干涉方法,方法的具体步骤如下:
步骤一、构建双波长相移干涉仪,称为双波长实际干涉仪,获得被测非球面波前与两个单波长λ1和λ2对应的包裹相位为当标准平面镜由移相器驱动,完成对不同波长相应干涉条纹的相移时,由探测器采集相应波长的相移干涉图,并利用相移公式进一步可求得两种波长的实际包裹相位。
步骤二、建模基于部分补偿法的双波长理想干涉仪,求解双波长理想干涉仪剩余波前对应的包裹相位根据双波长实际干涉仪光学元件位置和参数建模双波长理想干涉仪,获得双波长理想干涉仪像面的剩余波前,则对应于波λ1和λ2的剩余波前的包裹相位的求解分别表示为:
式(1a)和(1b)中的符号“\”表示取余。
步骤三、求解只与被测非球面面形误差相关的包裹相位消除已知波前变化量ΔW0。步骤一获得的包裹相位分别减去步骤二获得的包裹相位则与被测非球面面形误差相关的包裹相位分布表示为:
步骤四、求解对应于合成波长λeq只与被测非球面面形误差相关的像面波前Weq,E。根据求解Weq,E,对应于合成波长λeq的、只与被测非球面面形误差相关的像面波前表示为:
式(3)中合成波长λeq为λ1﹒λ2/(λ21),为合成波长λeq的相邻像素包裹相位之差,ΔMeq为EWPUA的干涉级次修正因子。
步骤五、求解只与被测非球面面形误差相关的像面波前W1E或W2E。这里以波长λ1为例,获得对应于波长λ1的、只与被测非球面面形误差相关的像面波前分别表示为:
式(4)中,为单波长λ1的相邻像素包裹相位之差,ΔM1E为求解波长λ1的干涉级次修正因子。
步骤六、采用误差分离和消除算法,消除未知波前变化量ΔWE,求解修正后只与面形误差相关的单波长像面波前W1T或W2T。采用像面波前Weq,E修正单波长的像面波前W1E或W2E,获得修正后的单波长像面波前W1rev或W2rev。在相位解包裹路径上将Weq,E和W1E差值的相邻像素对应波前相减后表示为:
△Weq,dif=Weq,E-W1E-[Weq,E'-W1E'] (5)
式(5)中,Weq,E′和W1E′分别为在解包裹路径上与Weq,E和W1E相邻像素的波前。修正后的解包裹波前表示为:
式(6)中,<>表示对其中的取值进行四舍五入取整。
步骤七、求解被测波前W1T或W2T。运用W10或W20和修正后的单波长像面波前W1rev或W2rev,获得W1T或W2T。则对应于波长为λ1的实际干涉仪像面波前表示为:
W1T=W1rev+W10 (7)
同理,W2T也可通过步骤五至七求解。
步骤八、在已建模的双波长理想干涉仪中,在波长λ1或λ2的理想干涉仪的标准非球面上附加泽尼克边缘系数,用来描述被测非球面的面形误差,然后以Zernike边缘系数为变量,以步骤六中求解的像面波前W1T或W2T为优化目标,在光学设计软件中利用光线追迹和逆向迭代优化,精确重构出被测非球面的面形误差Wasp_fig
为了实现基于部分补偿法的双波长相移干涉方法的技术方案,本发明还提出该方法的实现装置,采用如下技术方案:
一种基于部分补偿法的双波长相移干涉方法的实现装置,包括:第一激光器、第二激光器、第一狭缝、第二狭缝、第一平面镜、第一分束镜、扩束镜、准直镜、第二分束镜、标准平面镜、移相器、部分补偿镜、被测非球面、第二平面镜、成像镜头、含探测器的干涉图采集组件;
第一激光器和第二激光器输出不同工作波长的激光束,两束激光光束分别经第一平面镜反射及第一分束镜透射后合光,合光后的双波长激光光束经扩束镜和准直镜后变为宽口径平行光束,宽口径平行光束透过第二分束镜后入射到标准平面镜分为透射检测光束和反射参考光束。透射检测光束入射到部分补偿检测光路:先经部分补偿镜透射、再经被测非球面反射返回,再次透过部分补偿镜和标准平面镜。反射参考光束和返回的透射检测光束发生干涉,干涉光束经第二分束镜反射后再经第二平面镜反射入射到成像镜头,最后成像于探测器。
所述的第一激光器和第二激光器为单波长激光器;所述的部分补偿检测光路包括部分补偿镜和被测非球面,所述的部分补偿镜一般为单透镜;在第二分束镜的反射光方向配置有第二平面镜,在第二平面镜后面配置成像镜头;所述的探测器为稀疏阵列传感器,避免干涉条纹频率超过奈奎斯特频率时,探测器记录的干涉条纹图对比度下降,相位信息丢失或失真。
本发明的特点及有益效果在于:
(1)提出一种基于部分补偿法的双波长相移干涉方法,其特点在于,部分补偿法和双波长相移干涉的结合,不需要针对被测非球面和干涉仪设计特定的消色差补偿器,部分补偿镜多为单透镜,从而降低了系统设计难度和成本,缩短了仪器设计周期,利于实现非球面的通用化测量。
(2)提出一种基于部分补偿法的双波长相移干涉方法,其特点在于,通过方法的步骤一至三消除干涉仪组件和补偿镜引入的色差及其导致的已知被测波前变化量,从而可获得只与被测波前面形误差相关的像面波前。
(3)提出一种基于部分补偿法的双波长相移干涉方法,其特点在于,通过方法的步骤四到七消除了非球面面形误差导致的未知被测波前变化量,从而可获得被测波前像面的高精度求解。
(4)提出一种基于部分补偿法的双波长相移干涉方法,其特点在于,通过方法的步骤八中的逆向迭代优化重构,可实现被测非球面面形误差的高精度重构。
(5)提出实现基于部分补偿法的双波长相移干涉方法的实现装置,其特点在于,采用改进型斐索干涉光路结构,通过共光路结构保证双波长相移干涉仪的同轴度,降低仪器调试难度。
(6)本发明的方法及实验装置可以解决小面形误差的陡度非球面,以及大面形误差的模压非球面和自由曲面的加工质量测量难题。
附图说明
图1为部分补偿法双波长相移干涉方法流程图;
图2为基于部分补偿法的双波长相移干涉非球面测量方法的实现装置;
图3为与椭球面相关的W1E、WeqE和W1rev与对应真值的偏差;
图4为求解的椭球面Wasp_fig及其与真值的偏差;
图中标记:1-第一固体激光器、2-第二固体激光器、3-第一狭缝、4-第二狭缝、5-第一平面镜、6-第一分束镜、7-扩束镜、8-准直镜、9-分束镜、10-标准平面镜、11-移相器、12-部分补偿镜、13-被测非球面、14-第二平面镜、15-成像镜头、16-稀疏阵列探测器。
具体实施方式
对于双波长相移干涉测量,如果在相位测量过程中被测波前WT始终不变,则光强分布可表示为:
然而,实际分时采集λ1和λ2的两组数据时被测波前会发生改变,如果在用λ2进行测量时,被测波前与用λ1测量时改变了△WT,若定义Meq=λeq/[(λ12)/2]为误差放大倍率,则式(1)变为:
式(2)表明了双波长相移干涉的误差放大效应,在传统的双波长相移干涉仪中,需要针对干涉仪和被测非球面设计消色差的零位补偿镜,即使如此,也只能测量面形误差较小的非球面,才能使双波长相移干涉在零位测量状态,否则大面形误差的非球面导致从非球面返回检测光和参考不共路,违背零位条件的双波长相移干涉的误差放大效应将对这些波前变化量进行放大,从而导致双波长相移干涉无法工作。
本发明提出的方法主要针对干涉仪及其补偿透镜的色差和被测非球面的面形误差导致的被测波前改变量进行消除,也是双波长相移干涉实际应用面临的最大问题。
在双波长理想干涉仪中,非球面设置为不包含面形误差的标准非球面,通过光线追迹出对应的两个单波长λ1和λ2的剩余波前W10和W20,该剩余波前之差为被测波前已知变化量,即为干涉仪及其补偿透镜的色差。在双波长λ1和λ2的实际干涉仪中,被测非球面为标准非球面附加面形误差后的非球面,与面形误差相关的被测波前为W1E和W2E,则被测非球面波前表示为:
W1T=W10+W1E (3a)
W2T=W20+W2E (3b)
式(3a)和(3b)等式两边分别相减,则波长λ2的被测波前相对于波长λ1的被测波前的变化量ΔWT表示为:
△WT=△W0+△WE (4)
式(4)中被测波前变化量分ΔWT为ΔW0和ΔWE,ΔW0为与被测非球面面形误差无关的色差导致的已知波前变化量,ΔWE为与被测非球面面形误差相关的未知波前变化量。
下面结合附图和实施例对本发明提出的方法和实现装置做详细说明,阐释如何利用基于部分补偿法的双波长相移干涉方法中的步骤,在所提出的实现装置中,先后消除与色差相关的已知波前变化量和与面形误差相关的未知波前变化量,并运用逆向迭代优化重构技术,最终实现被测非球面面形误差的高精度测量。
实施例1:大非球面度椭球面测量。
本实施例的实现装置如图2所示,第一激光器(1)和第二激光器(2)输出的工作波长λ1和λ2分别为532nm和556nm,产生的等效合成波长λeq为12.325μm,可测量的最大相邻光程差之差为6.162μm。
本实施例中,被测椭球面(13)的通光口径D为580mm,非球面顶点曲率半径R0为-1179.447mm,非球面相对口径为D/R0为1/2,被测椭球面相对于顶点球和最佳参考球的最大非球面度分别为281.85μm和67.80μm,是一个非球面度相当大的非球面,本实施例在大非球面度椭球面附加小面形误差产生陡度波前。
设计的用于补偿椭球面的部分补偿镜(12)的玻璃材料为K9,前后表面曲率半径r1和r2分别为240mm和-153mm,部分补偿镜(12)的厚度为20mm,是一个简单的双凸透镜,很容易实现高精度的加工和检测,部分补偿镜(12)和被测椭球面(13)的距离为1349mm。
利用基于部分补偿法的双波长相移干涉方法,在所提出的实现装置中,实现椭球面面形误差的高精度测量具体实施如下:
1)通过双波长相移干涉仪获得被测波前全口径范围的包裹相位分布在标准椭球面上附加面形误差,预设面形误差的PV值和RMS值分别为0.7623λ1和0.2199λ1。此时,波长λ1和λ2的实际干涉仪像面上理论波前的PV值分别为366.8741λ1和354.0456λ2
2)建模双波长的理想干涉仪,通过光线追迹获得的剩余波前W10和W20,其PV值分别为4.6227λ1和33.1537λ2,则全口径范围内波长λ1和λ2的剩余波前点对点相减之差ΔW0值的PV值为17.2004μm,剩余波前变化量很大,如果不进行消除,将严重影响后续包裹相位的解相精度。
3)通过本发明提出方法的步骤一到三,首先从被测波前包裹相位消除了△W0这部分与色差相关的已知波前变化量,获得只与被测波前的面形误差相关包裹相位接着通过本发明提出方法的步骤四,采用合成波长解包裹算法对解包裹,求解出对应于合成波长λeq的解包裹波前Weq,E,该波前没有2π模糊性数据,但存在未知波前变化分量的误差放大效应;通过步骤五,采用单波长解包裹算法分别对解包裹获得与λ1和λ2对应的解包裹波前W1E和W2E,该波前只与被测非球面的面形误差相关,但均含有未能正确解包裹的2π模糊性数据。
4)通过本发明提出方法的步骤六,对Weq,E和W1E、W2E应用合成波长修正单波长算法,即用Weq,E修正W1E或W2E中有2π模糊性的解包裹波前获得与被测椭球面面形误差相关的像面波前W1rev、W2rev
5)通过本发明提出方法的步骤七,运用剩余波前W10、W20,像面波前W1rev、W2rev,求解出对应于波长λ1和λ2的被测波前W1T和W2T,其PV值分别为366.8741λ1和354.0456λ2
图3为与椭球面相关的W1E、WeqE和W1rev与对应真值的偏差,(a)、(b)和(c)分别为单波长解包裹算法、合成波长解包裹算法以及本发明提出的方法求解的与面形误差相关像面波前与对应真值的偏差。图3表明,通过本发明提出的方法中的步骤一到七,获得与预设真值一致的被测波前测量值,仿真偏差值仅为10-13λ1
从本实施例上述数据可以看出,全口径范围内波长λ1和λ2的被测波前点对点相减之差△WT值的PV值达17.8352μm,即测量波长λ2时的被测波前PV值相对于λ1变化了33.5248λ1。如果对被测波前变化量不进行误差分离和消除,双波长相移干涉仪中的误差放大效应,使其在像面波前被放大后将产生近百个波长的波前变化量,这将导致根本不能获得实际干涉仪对应的像面波前。
5)通过本发明提出方法的步骤八,在建模的双波长理想干涉仪中,将求解的被测波前W1T和W2T,分别作为单波长λ1和λ2理想干涉仪的目标波前。标准椭球面上附加37项泽尼克边缘系数作为变量,利用光线追迹和逆向迭代优化,重构出被测非球面的面形误差Wasp_fig
逆向迭代优化法重构求解的椭球面面形误差PV值和RMS值分别为0.7672λ1和0.2213λ1。图4所示为本实施例中求解的椭球面面形误差Wasp_fig及其与真值的偏差,与预设真值的PV值偏差和RMS值偏差分别优于λ1/200和λ1/500。如果采用传统干涉仪中面形误差和干涉仪像面波前的二分之一近似函数关系求解将引入非常大的误差。
至此,在本实施例中,通过本发明提出的方法和实验装置获得了大非球面度被测椭球面面形误差的高精度测量。
实施例2:大面形误差抛物面测量。
本实施例中,被测非球面(13)为大面形误差的抛物面,其通光口径D为108mm,顶点曲率半径R0为-1727.2mm,抛物面相对口径为D/R0为1/16,被测抛物面相对于顶点球和最佳参考球的最大非球面度分别为0.206μm和0.052μm,是一个浅度非球面。
设计的用于补偿抛物面的部分补偿镜(12)的玻璃材料为K9,前后表面曲率半径r1和r2分别为760mm和-4965mm,也是一个简单的双凸透镜,部分补偿镜(12)的厚度为20mm,部分补偿镜(12)和被测非球面(13)的距离为2985mm。
运用本发明提出的方法,在所提出的实现装置中获得抛物面的面形误差测量的具体实施步骤,与实施例1中相同,不予重复。实施例2与实施1的区别在于,在抛物面上预设的面形误差的PV值和RMS值分别为152.4616λ1和39.6496λ1,因此,被测非球面(13)是一个大面形误差的非球面,与波长λ1和λ2对应的被测波前W1T和W2T的PV值分别为483.5274λ1和464.2669λ2,求解的被测抛物面面形误差的PV值和RMS值分别为152.4616λ1和39.6496λ1,与预设真值的偏差的PV值和RMS值偏差分别为0.7830λ1和0.2281λ1
至此,在大面形误差的抛物面实施例中,通过本发明提出的方法和实验装置获得了其面形误差的高精度测量。

Claims (5)

1.一种基于部分补偿法的双波长相移干涉非球面测量方法,其特征在于包括以下步骤:
步骤一、构建部分补偿法双波长相移干涉仪,称为双波长实际干涉仪,分时采集λ1和λ2的相移干涉图,获取与两个单波长λ1和λ2对应的被测非球面波前的包裹相位
步骤二、建模基于部分补偿法的双波长理想干涉仪,其中非球面设置为不包含面形误差的标准非球面,通过光线追迹获得双波长理想干涉仪像面的剩余波前W10和W20,求解与两个单波长对应的剩余波前的包裹相位求解公式如下式(1a)(1b)所示,符号“\”表示取余:
步骤三、采用误差分离和消除算法,首先消除已知波前变化量ΔW0。将步骤一的结果,即被测非球面波前的包裹相位分别减去步骤二的结果,即对应的剩余波前包裹相位获得只与被测非球面面形误差相关的包裹相位分布求解式如下式(2a)(2b)所示:
步骤四、根据求解对应于合成波长λeq的,只与被测非球面面形误差相关的像面波前Weq,E,实现如下式(3)所示:
式(3)中的为合成波长λeq的相邻像素包裹相位之差,ΔMeq为干涉级次修正因子;
步骤五、根据求解被测波前中对应于λ1和λ2的,只与被测非球面面形误差相关的波前W1E或W2E,以波长λ1为例,实现如下式(4)所示:
ΔM1E为求解时的干涉级次修正因子;
步骤六、采用误差分离和消除算法,消除未知波前变化量ΔWE。根据Weq,E,W1E或W2E,用Weq,E修正具有2π模糊性的W1E或W2E数据,获得对应于波长λ1或λ2的,只与被测非球面面形误差相关的像面波前W1rev或W2rev。实现如下式(5)所示:
ΔWeq,dif=Weq,E-W1E-[Weq,E'-W1E'] (5a)
式(5a)中,Weq,E′和W1E′分别为在解包裹路径上与Weq,E和W1E相邻像素的波前,式(5b)中,<>表示对其中的取值进行四舍五入取整;
步骤七、根据W10、W1rev,或W20、W2rev,求解双波长实际干涉仪的解包裹波前W1T或W2T,求解式如下式(6a)和(6b)所示:
W1T=W1rev+W10 (6a)
W1T=W1rev+W10 (6b)
步骤八、采用逆向迭代优化法,重构被测非球面的面形误差Wasp_fig。在建模的双波长理想干涉仪中,在波长λ1或λ2的理想干涉仪的标准非球面上附加泽尼克边缘系数,用来描述被测非球面的面形误差,然后以Zernike边缘系数为变量,以W1T或W2T为优化目标,在光学设计软件中利用光线追迹和逆向迭代优化,从而精确重构出被测非球面的面形误差Wasp_fig
2.一种基于部分补偿法的双波长相移干涉的实现装置,其特征在于,包括:第一激光器(1)、第二激光器(2)、第一狭缝(3)、第二狭缝(4)、平面镜(5)、第一分束镜(6)、扩束镜(7)、准直镜(8)、第二分束镜(9)、标准平面镜(10)、移相器(11)、部分补偿镜(12)、被测非球面(13)、第二平面镜(14)、成像镜头(15)、含稀疏阵列传感器的干涉图采集组件(16);
在第一激光器(1)的输出光束方向配置有第一狭缝(3)和第一平面镜(5),在第二激光器(2)的输出光束方向配置有第二狭缝(4)和第一分束镜(6),在第一平面镜(5)的反射光方向配置有第一分束镜(6),在第一分束镜(6)的透射光方向配置有扩束镜(7)和准直镜(8),在准直镜(8)输出的宽口径平行光束方向配置有第二分束镜(9),在第二分束镜(9)的透射光方向配置标准平面镜(10),在标准平面镜(10)位置配置有移相器(11),在标准平面镜(10)的透射光方向配置有部分补偿镜(12),在部分补偿镜(12)的后面配置被测非球面(13),在第二分束镜(9)的反射光方向配置有第二平面镜(14),在第二平面镜(14)的反射光方向配置有成像镜头(15),在成像镜头(15)的输出光方向配置包含探测器的干涉图采集组件(16)。
3.根据权利要求2所述的一种基于部分补偿法的双波长相移干涉的实现装置,其特征在于:第一固体激光器(1)和第二固体激光(2)输出不同工作波长的激光器,产生等效合成波长λeq,测量的最大相邻光程差之差为λeq/2。
4.根据权利要求2所述的一种基于部分补偿法的双波长相移干涉的实现装置,其特征在于:所述的部分补偿镜(12)包括单透镜,用以部分补偿被测非球面(13)的法线像差和减少被测波前斜率。
5.根据权利要求2所述的一种基于部分补偿法的双波长相移干涉的实现装置,其特征在于:所述干涉图采集组件(16)包括稀疏阵列传感器。
CN201711019962.1A 2017-10-27 2017-10-27 基于部分补偿法的双波长相移干涉非球面测量方法及装置 Active CN107764203B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711019962.1A CN107764203B (zh) 2017-10-27 2017-10-27 基于部分补偿法的双波长相移干涉非球面测量方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711019962.1A CN107764203B (zh) 2017-10-27 2017-10-27 基于部分补偿法的双波长相移干涉非球面测量方法及装置

Publications (2)

Publication Number Publication Date
CN107764203A true CN107764203A (zh) 2018-03-06
CN107764203B CN107764203B (zh) 2019-08-23

Family

ID=61270681

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711019962.1A Active CN107764203B (zh) 2017-10-27 2017-10-27 基于部分补偿法的双波长相移干涉非球面测量方法及装置

Country Status (1)

Country Link
CN (1) CN107764203B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108627084A (zh) * 2018-04-27 2018-10-09 华中科技大学 一种基于静止的迈克尔逊干涉仪的激光器波长校准系统
CN108955569A (zh) * 2018-09-27 2018-12-07 成都太科光电技术有限责任公司 大口径长焦距菲索型球面干涉测试装置
CN109029291A (zh) * 2018-08-16 2018-12-18 北京理工大学 结合激光差动共焦定位的非球面参数误差干涉测量方法
CN109737892A (zh) * 2019-01-08 2019-05-10 北京卫星环境工程研究所 基于区域定位拟合算法的数字莫尔移相干涉面形测量方法
CN110307962A (zh) * 2019-06-28 2019-10-08 苏州维纳仪器有限责任公司 检测光学系统任意波长透射波前的方法
CN110487205A (zh) * 2019-07-31 2019-11-22 北京理工大学 结合色散共焦定位的非球面参数误差干涉测量方法
CN111209689A (zh) * 2020-02-14 2020-05-29 北京理工大学 非零干涉非球面测量回程误差去除方法及装置
CN111240010A (zh) * 2020-01-08 2020-06-05 北京理工大学 一种用于自由曲面测量的可变形镜面形设计方法及装置
CN111947592A (zh) * 2020-06-29 2020-11-17 中国科学院上海光学精密机械研究所 一种动态双波长移相干涉测量装置和测量方法
CN113587844A (zh) * 2021-07-27 2021-11-02 中国科学院长春光学精密机械与物理研究所 移相干涉测量系统及测量方法
CN114964033A (zh) * 2022-02-14 2022-08-30 上海大学 一种使用波长移相法测量超表面形貌分布的方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002071327A (ja) * 2000-09-01 2002-03-08 Nippon Telegr & Teleph Corp <Ntt> 表面形状測定方法及び表面形状測定器
CN1587950A (zh) * 2004-07-08 2005-03-02 北京理工大学 一种用部分补偿透镜实现非球面面形的干涉测量方法
WO2009088120A1 (en) * 2008-01-04 2009-07-16 Nano System Co., Ltd. Self-compensating apparatus of interference-fringe-generating position to compensate optical path difference result from objectives lens working distance variation for environment temperature variation in regard of interference lens
CN101876540A (zh) * 2010-05-07 2010-11-03 中国科学院光电技术研究所 基于多波前透镜补偿器的非球面绝对测量系统
CN102506750A (zh) * 2011-10-28 2012-06-20 中国科学院长春光学精密机械与物理研究所 部分补偿非球面反射镜面形检测方法
CN103630086A (zh) * 2013-11-13 2014-03-12 华南师范大学 一种基于单色ccd的双波长同时相移干涉测量方法
CN103776388A (zh) * 2012-03-02 2014-05-07 浙江大学 非球面非零位干涉检测中部分补偿透镜对准方法
CN105758295A (zh) * 2016-03-22 2016-07-13 华南师范大学 一种时空混合匹配的双波长同时相移干涉测量方法
CN106949853A (zh) * 2017-04-12 2017-07-14 北京理工大学 基于液晶计算全息图的同步移相干涉测量系统及方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002071327A (ja) * 2000-09-01 2002-03-08 Nippon Telegr & Teleph Corp <Ntt> 表面形状測定方法及び表面形状測定器
CN1587950A (zh) * 2004-07-08 2005-03-02 北京理工大学 一种用部分补偿透镜实现非球面面形的干涉测量方法
WO2009088120A1 (en) * 2008-01-04 2009-07-16 Nano System Co., Ltd. Self-compensating apparatus of interference-fringe-generating position to compensate optical path difference result from objectives lens working distance variation for environment temperature variation in regard of interference lens
CN101876540A (zh) * 2010-05-07 2010-11-03 中国科学院光电技术研究所 基于多波前透镜补偿器的非球面绝对测量系统
CN102506750A (zh) * 2011-10-28 2012-06-20 中国科学院长春光学精密机械与物理研究所 部分补偿非球面反射镜面形检测方法
CN103776388A (zh) * 2012-03-02 2014-05-07 浙江大学 非球面非零位干涉检测中部分补偿透镜对准方法
CN103630086A (zh) * 2013-11-13 2014-03-12 华南师范大学 一种基于单色ccd的双波长同时相移干涉测量方法
CN105758295A (zh) * 2016-03-22 2016-07-13 华南师范大学 一种时空混合匹配的双波长同时相移干涉测量方法
CN106949853A (zh) * 2017-04-12 2017-07-14 北京理工大学 基于液晶计算全息图的同步移相干涉测量系统及方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
WANG XIAOKUN: "Testing of aspheric mirror by non-null compensation method", 《INFRARED AND LASER ENGINEERING》 *
刘惠兰等: "利用部分补偿透镜进行非球面面形测量", 《北京理工大学学报》 *
吕宁等: "一种基于部分补偿与稀疏阵列的深度非球面测量新方法", 《四川大学学报(自然科学版)》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108627084A (zh) * 2018-04-27 2018-10-09 华中科技大学 一种基于静止的迈克尔逊干涉仪的激光器波长校准系统
CN109029291A (zh) * 2018-08-16 2018-12-18 北京理工大学 结合激光差动共焦定位的非球面参数误差干涉测量方法
CN108955569A (zh) * 2018-09-27 2018-12-07 成都太科光电技术有限责任公司 大口径长焦距菲索型球面干涉测试装置
CN108955569B (zh) * 2018-09-27 2023-10-27 成都太科光电技术有限责任公司 大口径长焦距菲索型球面干涉测试装置
CN109737892A (zh) * 2019-01-08 2019-05-10 北京卫星环境工程研究所 基于区域定位拟合算法的数字莫尔移相干涉面形测量方法
CN110307962A (zh) * 2019-06-28 2019-10-08 苏州维纳仪器有限责任公司 检测光学系统任意波长透射波前的方法
CN110307962B (zh) * 2019-06-28 2020-10-27 苏州维纳仪器有限责任公司 检测光学系统任意波长透射波前的方法
CN110487205A (zh) * 2019-07-31 2019-11-22 北京理工大学 结合色散共焦定位的非球面参数误差干涉测量方法
CN111240010B (zh) * 2020-01-08 2021-04-20 北京理工大学 一种用于自由曲面测量的可变形镜面形设计方法及装置
CN111240010A (zh) * 2020-01-08 2020-06-05 北京理工大学 一种用于自由曲面测量的可变形镜面形设计方法及装置
CN111209689A (zh) * 2020-02-14 2020-05-29 北京理工大学 非零干涉非球面测量回程误差去除方法及装置
CN111947592B (zh) * 2020-06-29 2021-10-08 中国科学院上海光学精密机械研究所 一种动态双波长移相干涉测量装置和测量方法
CN111947592A (zh) * 2020-06-29 2020-11-17 中国科学院上海光学精密机械研究所 一种动态双波长移相干涉测量装置和测量方法
CN113587844A (zh) * 2021-07-27 2021-11-02 中国科学院长春光学精密机械与物理研究所 移相干涉测量系统及测量方法
CN114964033A (zh) * 2022-02-14 2022-08-30 上海大学 一种使用波长移相法测量超表面形貌分布的方法

Also Published As

Publication number Publication date
CN107764203B (zh) 2019-08-23

Similar Documents

Publication Publication Date Title
CN107764203B (zh) 基于部分补偿法的双波长相移干涉非球面测量方法及装置
CN102589416B (zh) 用于非球面测量的波长扫描干涉仪及方法
CN104215176B (zh) 高精度光学间隔测量装置和测量方法
CN106768280B (zh) 一种基于多波长无透镜傅里叶变换数字全息的振动检测装置
US12018930B2 (en) Transient digital moire phase-shifting interferometric measuring device and method for the surface shape of an optical element
CN105318847A (zh) 基于系统建模的非球面非零位环形子孔径拼接方法
CN102589440A (zh) 连续变角度数字全息干涉测量的方法和装置
CN100405003C (zh) 条纹计测装置的变换系数校正方法和装置及条纹计测装置
CN107064067B (zh) 一种双干涉仪的空气折射率廓线的测量系统及方法
CN110082071B (zh) 一种直角棱镜光学平行差的测量装置及方法
Zhang et al. High-precision calibration method for shear ratio based on the shearing wavefront feature extraction of a phase plate
CN114322848A (zh) 一种球面波前曲率半径测量装置及测量方法
CN104792269B (zh) 一种对线性相移误差不敏感的光纤端面高度值的解算方法
CN110631510B (zh) 一种基于迈克尔逊结构的高精度测角装置及测角方法
CN105674875A (zh) 一种全视场低频外差点衍射干涉仪
CN106352985B (zh) 一种非对称空间外差光谱仪结构
Baer et al. Measurement of aspheres and free-form surfaces in a non-null test interferometer: reconstruction of high-frequency errors
US20070019210A1 (en) Time-delayed source and interferometric measurement of windows and domes
Stašík et al. Aspheric surface measurement by absolute wavelength scanning interferometry with model-based retrace error correction
US5777738A (en) Interferometric measurement of absolute dimensions of cylindrical surfaces at grazing incidence
CN103575229A (zh) 共光路型多重倾斜波面补偿非零位干涉测量装置
CN112781727A (zh) 一种基于棱镜的横向剪切干涉光谱成像仪及成像方法
Liu et al. A novel aspheric surface testing method using part-compensating lens
CN205809297U (zh) 一种激光测距系统
CN113029526B (zh) 一种综合孔径共相误差估计方法及装置

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