CN106017306B - 基于最小二乘旋转匹配的干涉标准件误差分离方法 - Google Patents

基于最小二乘旋转匹配的干涉标准件误差分离方法 Download PDF

Info

Publication number
CN106017306B
CN106017306B CN201610300921.9A CN201610300921A CN106017306B CN 106017306 B CN106017306 B CN 106017306B CN 201610300921 A CN201610300921 A CN 201610300921A CN 106017306 B CN106017306 B CN 106017306B
Authority
CN
China
Prior art keywords
standard component
interference
focus lens
long focus
ray tracing
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
CN201610300921.9A
Other languages
English (en)
Other versions
CN106017306A (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201610300921.9A priority Critical patent/CN106017306B/zh
Publication of CN106017306A publication Critical patent/CN106017306A/zh
Application granted granted Critical
Publication of CN106017306B publication Critical patent/CN106017306B/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
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02055Reduction or prevention of errors; Testing; Calibration
    • 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)
  • Testing Of Optical Devices Or Fibers (AREA)

Abstract

本发明公开了一种基于最小二乘旋转匹配的干涉标准件误差分离方法,基于矢量光学理论对长焦透镜零位干涉测试系统构建了三维光线追迹仿真系统,利用仿真系统将干涉标准件面形偏差从长焦透镜测试结果中分离,解决了干涉标准件面形偏差在最终的长焦透镜测量结果中引入误差信息的问题,实现长焦透镜透射波前的GRMS和PSD1高精度测量。在误差分离的过程中,关键是找到三维光线追迹仿真系统中干涉标准件的空间姿态与干涉仪检测时的姿态一致处,即干涉标准件的空间姿态匹配问题。提出利用最小二乘原理,通过不断旋转三维光线追迹仿真系统中干涉标准件的空间姿态,使得仿真系统输出波前与干涉仪测试波前差值最小,完成姿态匹配。

Description

基于最小二乘旋转匹配的干涉标准件误差分离方法
技术领域
本发明属于光学元件面形检测领域,具体涉及一种基于最小二乘旋转匹配的干涉标准件误差分离方法。
背景技术
高功率固体激光驱动器采用多路大口径结构,通光口径大、承载激光能量高、系统中使用了大量大口径、高精度的光学元件。高功率固体激光驱动器中光学元件的中频段空间频率成分将导致光束的高频调制与系统的非线性增长,造成光学元件的丝状破坏和降低光束的可聚焦功率。其关于光学元件的评价指标主要有PSD和GRMS,其中PSD是一种基于频谱分析的评价指标,描述面形误差的频率分布,其主要用来评价中频段的波纹差,可以满足光学元件在中频段的质量要求;GRMS是针对元件低频波前的一个关键参数,最终影响系统的聚焦性能。空间滤波透镜、打靶透镜等长焦(8m-30m)透镜的透射波前的GRMS和PSD1会影响激光焦斑分布、激光过孔等系统指标,需要严格控制和精密测量。
在长焦透镜透射波前的干涉检测过程中,大口径的干涉标准件的面形偏差难以控制,其面形偏差会在最后的测量结果中引入误差信息。由于在长焦透镜透射波前的检测过程中,最后由干涉场得到的波前不仅仅是被测透镜的透射波前,还包括了干涉标准件的波前误差。因此,要利用自准直零位干涉检测法(原理如图1所示)实现高精度的长焦透镜透射波前GRMS和PSD1的测量,必须研究如何从最终的测量结果中分离干涉标准件的GRMS和PSD1的数学方法。
发明内容
本发明的目的在于提供一种基于最小二乘旋转匹配的干涉标准件误差分离方法,解决了干涉标准件面形偏差在最终的长焦透镜测量结果中引入误差信息的问题,实现长焦透镜透射波前的GRMS和PSD1高精度测量。
实现本发明目的的技术解决方案为:一种基于最小二乘旋转匹配的干涉标准件误差分离方法,方法步骤如下:
步骤1,基于矢量光学理论对长焦透镜零位干涉测试系统构建三维光线追迹仿真系统,向上述三维光线追迹仿真系统输入长焦透镜的参数。
步骤2,利用干涉仪对上述长焦透镜和干涉标准件组合进行测量,获得实际透射波前Wtest;同时利用干涉仪对上述干涉标准件进行测量,获得干涉标准件的实际面形WS,将实际透射波前Wtest、干涉标准件的实际面形WS输入上述三维光线追迹仿真系统。
步骤3,三维光线追迹仿真系统对长焦透镜零位干涉测试系统进行光线追迹,并输出带有干涉标准件面形误差的仿真透射波前Wt1
步骤4,不断旋转干涉标准件的空间姿态,即旋转干涉标准件的实际面形WS,利用最小二乘法找到三维光线追迹仿真系统中干涉标准件的姿态与干涉仪中相匹配的姿态,即当三维光线追迹仿真系统输出仿真透射波前Wt1与干涉仪测量的实际透射波前Wtest相减的差值最小时,说明二者姿态相匹配。
步骤5,输出与干涉标准件相匹配时对应的干涉系统仿真透射波前W'S,从干涉仪实际测量结果中分离得到长焦透镜的透射波前WT0=Wtest-W'S,即为干涉标准件面形误差分离后的结果。
上述干涉标准件为球面反射镜。
上述步骤1长焦透镜的参数包括长焦透镜的口径、曲率半径、二次项系数和焦距。
本发明与现有技术相比,其显著优点在于:(1)将传统零位干涉测量与三维光线追迹仿真系统相结合进行干涉标准件误差分离,提高了测量精度。
(2)误差分离过程使用最小二乘法进行空间姿态匹配,同时兼顾干涉标准件误差分离速度与精度,使得分离速度更快与精度更高。
附图说明
图1为本发明的基于最小二乘旋转匹配的干涉标准件误差分离方法的流程图。
图2为本发明的长焦透镜零位干涉测试系统检测原理图。
图3为本发明的实施例1中长焦透镜误差分离后波前WT0结果图。
图4为本发明的实施例1中长焦透镜误差分离后波前WT0与真实波前的残差波前。
图5为本发明的实施例误差分离后波前一维PSD曲线图。
具体实施方式
下面结合附图对本发明作进一步详细描述。
结合图1和图2,一种基于最小二乘旋转匹配的干涉标准件误差分离方法,方法步骤如下:
步骤1,基于矢量光学理论对长焦透镜零位干涉测试系统构建三维光线追迹仿真系统,向上述三维光线追迹仿真系统输入长焦透镜的参数;
步骤2,利用干涉仪对上述长焦透镜和干涉标准件组合进行测量,获得实际透射波前Wtest;同时利用干涉仪对上述干涉标准件进行测量,获得干涉标准件的实际面形WS,将实际透射波前Wtest、干涉标准件的实际面形WS输入上述三维光线追迹仿真系统;
步骤3,三维光线追迹仿真系统对长焦透镜零位干涉测试系统进行光线追迹,并输出带有干涉标准件面形误差的仿真透射波前Wt1
步骤4,不断旋转干涉标准件的空间姿态,即旋转干涉标准件的实际面形WS,利用最小二乘法找到三维光线追迹仿真系统中干涉标准件的姿态与干涉仪中相匹配的姿态,即当三维光线追迹仿真系统输出仿真透射波前Wt1与干涉仪测量的实际透射波前Wtest相减的差值最小时,说明二者姿态相匹配;
步骤5,输出与干涉标准件相匹配时对应的干涉系统仿真透射波前W′S,从干涉仪实际测量结果中分离得到长焦透镜的透射波前WT0=Wtest-W′S,即为干涉标准件面形误差分离后的结果。
上述干涉标准件为球面反射镜。
上述步骤1长焦透镜的参数包括长焦透镜的口径、曲率半径、二次项系数、焦距。
实施例1
结合图1和图2,一种基于最小二乘旋转匹配的干涉标准件误差分离方法,方法步骤如下:
步骤1,基于矢量光学理论对长焦透镜零位干涉测试系统构建三维光线追迹仿真系统,向上述三维光线追迹仿真系统输入长焦透镜的参数;
步骤2,利用干涉仪对上述长焦透镜和干涉标准件组合进行测量,获得实际透射波前Wtest;同时利用干涉仪对上述干涉标准件进行测量,获得干涉标准件的实际透射波前WS,将干涉标准件的实际透射波前WS输入上述三维光线追迹仿真系统;
步骤3,三维光线追迹仿真系统对长焦透镜零位干涉测试系统进行光线追迹,并输出带有干涉标准件面形误差的仿真透射波前Wt1
步骤4,不断旋转干涉标准件的空间姿态,利用最小二乘法找到三维光线追迹仿真系统中干涉标准件的姿态与干涉仪中相匹配的姿态,即当三维光线追迹仿真系统输出仿真透射波前Wt1与干涉仪测量的实际透射波前Wtest相减的差值最小时,说明二者姿态相匹配;
步骤5,输出与干涉标准件相匹配时对应的干涉标准件仿真透射波前W′S,从干涉仪实际测量结果中分离得到长焦透镜的透射波前WT0=Wtest-W′S(如图3),即为干涉标准件面形误差分离后的结果。
将本方法的误差分离后透射波前与长焦透镜实际透射波前相减得到残差波前(如图4),残差波前的PV值为0.0750λ(其中λ=632.8nm)。对最终的误差分离后波前进行评价,中频波前均方根RMS值为1.1561nm,低频波前梯度均方根GRMS值为5.0650nm/cm,一维PSD1曲线(如图5)。

Claims (3)

1.一种基于最小二乘旋转匹配的干涉标准件误差分离方法,其特征在于,方法步骤如下:
步骤1,基于矢量光学理论对长焦透镜零位干涉测试系统构建三维光线追迹仿真系统,向上述三维光线追迹仿真系统输入长焦透镜的参数;
步骤2,利用干涉仪对上述长焦透镜和干涉标准件组合进行测量,获得实际透射波前Wtest;同时利用干涉仪对上述干涉标准件进行测量,获得干涉标准件的实际面形WS,将实际透射波前Wtest、干涉标准件的实际面形WS输入上述三维光线追迹仿真系统;
步骤3,三维光线追迹仿真系统对长焦透镜零位干涉测试系统进行光线追迹,并输出带有干涉标准件面形误差的仿真透射波前Wt1
步骤4,不断旋转干涉标准件的空间姿态,即旋转干涉标准件的实际面形WS,利用最小二乘法找到三维光线追迹仿真系统中干涉标准件的姿态与干涉仪中相匹配的姿态,即当三维光线追迹仿真系统输出仿真透射波前Wt1与干涉仪测量的实际透射波前Wtest相减的差值最小时,说明二者姿态相匹配;
步骤5,输出与干涉标准件相匹配时对应的干涉系统仿真透射波前W'S,从干涉仪实际测量结果中分离得到长焦透镜的透射波前WT0=Wtest-W′S,即为干涉标准件面形误差分离后的结果。
2.根据权利要求1所述的基于最小二乘旋转匹配的干涉标准件误差分离方法,其特征在于:上述干涉标准件为球面反射镜。
3.根据权利要求1所述的基于最小二乘旋转匹配的干涉标准件误差分离方法,其特征在于:上述步骤1长焦透镜的参数包括长焦透镜的口径、曲率半径、二次项系数和焦距。
CN201610300921.9A 2016-05-09 2016-05-09 基于最小二乘旋转匹配的干涉标准件误差分离方法 Active CN106017306B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610300921.9A CN106017306B (zh) 2016-05-09 2016-05-09 基于最小二乘旋转匹配的干涉标准件误差分离方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610300921.9A CN106017306B (zh) 2016-05-09 2016-05-09 基于最小二乘旋转匹配的干涉标准件误差分离方法

Publications (2)

Publication Number Publication Date
CN106017306A CN106017306A (zh) 2016-10-12
CN106017306B true CN106017306B (zh) 2018-07-06

Family

ID=57099565

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610300921.9A Active CN106017306B (zh) 2016-05-09 2016-05-09 基于最小二乘旋转匹配的干涉标准件误差分离方法

Country Status (1)

Country Link
CN (1) CN106017306B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107505684B (zh) * 2017-08-25 2020-04-28 南京理工大学 一种镜组的装调方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100598572B1 (ko) * 2005-08-01 2006-07-07 (주)펨트론 표면형상 측정장치 및 그 제어방법
CN101275826A (zh) * 2007-03-29 2008-10-01 富士能株式会社 非球面透镜的面偏移测量方法及装置
CN101451889A (zh) * 2007-12-04 2009-06-10 富士能株式会社 干涉仪装置的系统误差校正方法
CN101495833A (zh) * 2006-07-28 2009-07-29 卡尔蔡司Smt股份公司 确定光学表面的实际形状偏离理想形状的偏差的方法和装置
CN102175142A (zh) * 2011-02-12 2011-09-07 中国科学院光电技术研究所 透镜中心误差干涉测量系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6312373B1 (en) * 1998-09-22 2001-11-06 Nikon Corporation Method of manufacturing an optical system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100598572B1 (ko) * 2005-08-01 2006-07-07 (주)펨트론 표면형상 측정장치 및 그 제어방법
CN101495833A (zh) * 2006-07-28 2009-07-29 卡尔蔡司Smt股份公司 确定光学表面的实际形状偏离理想形状的偏差的方法和装置
CN101275826A (zh) * 2007-03-29 2008-10-01 富士能株式会社 非球面透镜的面偏移测量方法及装置
CN101451889A (zh) * 2007-12-04 2009-06-10 富士能株式会社 干涉仪装置的系统误差校正方法
CN102175142A (zh) * 2011-02-12 2011-09-07 中国科学院光电技术研究所 透镜中心误差干涉测量系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于旋转法的球面干涉仪系统误差标定方法;王乔方等;《云光技术》;20151231;第47卷(第1期);37-40页 *

Also Published As

Publication number Publication date
CN106017306A (zh) 2016-10-12

Similar Documents

Publication Publication Date Title
CN104567738B (zh) 光轴平行度精确测量系统及方法
CN106197258B (zh) 双通道双波长干涉检测装置
CN105588519A (zh) 利用相位差异相位恢复技术检测大口径望远镜面形的方法
CN105352707B (zh) 星敏感器光学系统倍率色差测试设备及测试方法
CN105651187A (zh) 菲涅尔双棱镜厚度的非接触间接测量方法
CN101980056B (zh) 基于m2因子和光闪烁指数确定大气湍流参数方法及装置
Li et al. Measurement of the pairwise kinematic Sunyaev-Zeldovich effect with Planck and BOSS data
CN104048619A (zh) 一种判断旋转轴对称非球面能否采用直接干涉检测的方法
CN103363927A (zh) 平台光电装备的任意轴距多光轴一致性检测装置及方法
CN103234480A (zh) 一种环形凸非球面的快速面形检测方法
CN103471521A (zh) 快速、准确的光学非球面的实时检测方法
CN104504197A (zh) 阿基米德螺旋线平面螺纹测量模型与其偏心参量修正方法
CN106017306B (zh) 基于最小二乘旋转匹配的干涉标准件误差分离方法
CN112903169A (zh) 爆炸应力波传播全场主应力定量分析实验方法
CN204064251U (zh) 透明材料厚度检测仪
CN103983609A (zh) 基于光谱干涉的透明材料折射率及厚度测量装置和测量方法
CN105352915A (zh) 一种折射率二维分布的动态测量方法
CN104777718A (zh) 一种大数值孔径光刻机投影物镜波像差检测方法
WO2017120730A1 (zh) 用渐变孔径光学积分法测量单模光纤模场直径的新方法
CN104280120B (zh) 一种光谱带宽测量方法和装置
CN205192728U (zh) 星敏感器光学系统倍率色差测试设备
CN103808413A (zh) 基于噪声分离的温度与比辐射率分离的方法及装置
CN103344209A (zh) 一种反射镜零重力面形测试方法
CN103163090B (zh) 一种用于反应堆厂房内部的钋气溶胶浓度检测系统
CN202710289U (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
GR01 Patent grant
GR01 Patent grant