CN106871819B - 基于最佳补偿位置的非球面顶点曲率半径误差测量方法 - Google Patents
基于最佳补偿位置的非球面顶点曲率半径误差测量方法 Download PDFInfo
- Publication number
- CN106871819B CN106871819B CN201710022242.4A CN201710022242A CN106871819B CN 106871819 B CN106871819 B CN 106871819B CN 201710022242 A CN201710022242 A CN 201710022242A CN 106871819 B CN106871819 B CN 106871819B
- Authority
- CN
- China
- Prior art keywords
- aspherical
- aspheric surface
- curvature radius
- vertex
- compensation position
- 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
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 11
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 15
- 239000000463 material Substances 0.000 claims description 7
- 230000004075 alteration Effects 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims 1
- 238000012986 modification Methods 0.000 claims 1
- 230000004048 modification Effects 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract 1
- 238000004556 laser interferometry Methods 0.000 abstract 1
- 239000011521 glass Substances 0.000 description 17
- 230000003287 optical effect Effects 0.000 description 6
- 238000005457 optimization Methods 0.000 description 4
- 238000005305 interferometry Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012938 design process Methods 0.000 description 2
- 238000013041 optical simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000208140 Acer Species 0.000 description 1
- 241000530560 Etheostoma sagitta Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005693 optoelectronics Effects 0.000 description 1
- 210000001747 pupil Anatomy 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004441 surface measurement Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/255—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures for measuring radius of curvature
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明涉及基于最佳补偿位置的非球面顶点曲率半径误差测量方法,属于非球面面型参数测量领域。旨在利用基于斜率非球面度的最接近比较球面,实现部分补偿法非球面顶点曲率半径误差的测量,构成一种全新的非球面参数误差测量方法。本发明为基于最佳补偿位置的非球面顶点曲率半径误差测量方法,用于激光干涉法测量非球面的顶点曲率半径误差。应用本方法,可以获得非球面的最佳补偿位置,通过测量理想非球面与实际非球面的最佳补偿位置之间的离焦距离,可以计算非球面的顶点曲率半径误差,最终实现非球面面型参数的高精度检测。
Description
技术领域
本发明涉及基于最佳补偿位置的非球面顶点曲率半径误差测量方法,属于非球面面型参数测量领域。旨在利用基于斜率非球面度的最接近比较球面,实现部分补偿法非球面顶点曲率半径误差的测量,构成一种全新的非球面参数误差测量方法。
背景技术
顶点曲率半径是光学非球面的重要面型参数,决定了非球面的基本性质。精确测量顶点曲率半径误差,对于光学非球面的加工和装调非常重要。通常情况下,利用接触法或非接触法可以获得被测面的面形轮廓,然后对面形轮廓直接进行曲率拟合,可以得到被测面的顶点曲率半径。顶点曲率半径的测量值与标称值的差值,即为该非球面的顶点曲率半径误差。
部分补偿法是一种新发展起来的非球面检测方法,属于干涉测量的范畴。在部分补偿法中,准直光经过补偿透镜后,其波前与非球面并不是完全吻合的,因此,反射光再次经过补偿透镜后,不再是准直光。当非准直反射光与参考准直光干涉时就会得到理想干涉条纹,实际干涉条纹与理想干涉条纹的差异就反映了被测非球面的面形误差。但是,该方法无法直接获得被测非球面的顶点曲率半径误差,这是目前需要解决的一大难题。
发明内容
本发明的目的是为了提供基于最佳补偿位置的非球面顶点曲率半径误差测量方法,利用本方法可以测量非球面的顶点曲率半径误差,无需扫描装置,结构简单,可实现非接触、全口径、速度快、精度高的测量。
本发明的目的是通过下述技术方案实现的。
基于最佳补偿位置的非球面顶点曲率半径误差测量方法,具体步骤如下:
步骤1:获取待测非球面名义参数。
名义参数包括待测非球面的口径、顶点曲率半径和非球面系数。
步骤2:根据待测非球面名义参数,设计部分补偿透镜P,获取其相关参数,用于非球面顶点曲率半径误差的测量。
相关参数包括部分补偿透镜P的第一面曲率半径,厚度,材料,第二面曲率半径和口径。
步骤3:在光学仿真软件中建立包含部分补偿透镜的虚拟干涉仪IR,并确定理想非球面的最佳补偿位置;
根据待测非球面名义参数和部分补偿透镜P的相关参数,确定理想非球面的最佳补偿位置,即补偿镜第二面到理想非球面顶点的轴向距离d1:
d1=d2+d3=d2+R0-K·|X0.86| (1)
其中,d1是部分补偿镜第二面到非球面顶点的距离;d2是部分补偿镜第二面到最接近比较球面球心的距离,通过补偿镜的入射光线坐标和出射光线坐标经代数计算进行确定;d3是最接近比较球面球心到非球面顶点的距离。R0是非球面的顶点曲率半径,K=-e2是二次曲面常数,e为偏心率。|X0.86|是以非球面的顶点为原点建立直角坐标系时,0.86口径处的点的横坐标绝对值,即该点的弓形高。
对于(1)式,-K·|X0.86|是非球面的法线像差,d3=R0-K·|X0.86|是法线像差关系式的变形,因此等式d3=R0-K·|X0.86|可以从二次非球面扩展到所有非球面。
步骤4:加工出部分补偿透镜P的实物,建立包含部分补偿透镜P的实际干涉仪IO,并确定实际非球面的最佳补偿位置,即补偿镜第二面到理想非球面顶点的轴向距离d′1。
根据实际干涉仪IO的干涉图进行定位,当干涉图的条纹最稀疏的时候,能够确定实际非球面的最佳补偿位置d′1。
步骤5:测量实际非球面与理想非球面的最佳补偿位置之间的离焦距离Δd:
Δd=d′1-d1 (2)
步骤6:计算实际非球面与理想非球面的最佳补偿位置之间的离焦距离与非球面顶点曲率半径误差之间的关系:
Δd=d′1-d1=(1-K·k(h,R0,K))·ΔR (3)
其中,非球面的口径为2D,则h=0.86D,k(h,R0,K)是通过非球面方程计算出的系数,具体形式为:
步骤7:计算非球面的顶点曲率半径误差:
ΔR=Δd/(1-K·k(h,R0,K)) (5)
有益效果
基于顶点曲率半径误差与最佳补偿位置之间的关系,通过干涉条纹来确定被测非球面的最佳补偿位置,即补偿镜第二面到被测非球面顶点的轴向距离,无需进行多位置的测量和记录,无需扫描装置,就可以获得非球面的顶点曲率半径误差,设计过程和系统结构简单,可实现非接触、全口径、速度快、精度高的测量。
附图说明
附图1是基于最佳补偿位置的非球面顶点曲率半径误差测量方法流程图;
附图2是非球面最佳补偿位置光路图。
具体实施方式
下面结合实施例对本发明进一步说明。
实施例1
基于最佳补偿位置的非球面顶点曲率半径误差测量方法,按以下方式实现:
建立基于最佳补偿位置的非球面顶点曲率半径误差测量系统的流程如附图1所示,具体实施步骤为:
步骤1:获取待测非球面名义参数。
名义参数包括待测非球面的口径、顶点曲率半径和非球面系数。
口径2D=76.2mm;顶点曲率半径R0=889mm;非球面系数K=-1。
步骤2:根据待测非球面名义参数,设计部分补偿透镜P,获取其相关参数,用于非球面顶点曲率半径误差的测量。
相关参数包括部分补偿透镜P的第一面曲率半径,厚度,材料,第二面曲率半径和口径。
根据待测非球面名义参数,依据专利ZL201010270226.5“一种以斜率为优化目标的部分补偿透镜的设计方法”设计部分补偿透镜P:采用目前光电子领域广泛应用的ZEMAX光学设计仿真软件作为设计平台,已知被测非球面口径为2D=76.2mm,顶点曲率半径为R0=889mm,二次曲面常数为K=-1,设计部分补偿透镜的具体步骤为:
第一步:设定系统参数,根据需要设定入瞳直径2D’=80mm和波长λ=532nm。
第二步:计算部分补偿透镜的光学参数,以部分补偿镜的第一面曲率半径R1,厚度d,第二面曲率半径R2为优化变量由ZEMAX软件自动优化。设置补偿镜第一面曲率半径R1为Infinity,厚度d=20mm,材料K9玻璃,折射率n=1.51630,第二面曲率半径R2=1000mm为初始参数。
第三步:设置反射光路及部分补偿系统初始结构参数,设置补偿镜第二面到非球面顶点的距离d0=1800mm作为优化变量由ZEMAX软件自动优化。
第四步:实现以斜率作为优化目标,需要在MeritFunCtion列表中作下述设置:设置MeritFunction菜单ToolS子菜单下DefaultMeritFunction的优化目标为SpotRadius,即弥散圆半径,设置类型为RMS;通过上述设置,ZEMAx软件即可对全口径光线进行分析控制,以优化弥散圆半径,而不是仅控制有限根光线在弥散圆上的位置。
运行Optimization自动优化程序即可对系统进行优化,可以获取分补偿透镜P的相关参数为:
部分补偿透镜P为双凸单透镜,参数为第一面曲率半径578.4mm,厚度22.0mm,材料K9玻璃,折射率n=1.51630,第二面曲率半径3350.0,口径100mm。此时系统弥散圆最大半径为4.269微米。
步骤3:在光学仿真软件中建立包含部分补偿透镜的虚拟干涉仪IR,并确定理想非球面的最佳补偿位置;
根据待测非球面名义参数和部分补偿透镜P的相关参数,确定理想非球面的最佳补偿位置,即补偿镜第二面到理想非球面顶点的轴向距离d1=d2+d3=d2+R0-K·|X0.86|,如附图2所示。
其中,H-部分补偿镜第二面与光轴的交点,O-最接近比较球面的球心,O’-非球面的顶点曲率中心,V-非球面与光轴的交点,A、B-非球面和最接近比较球面的切点。d1-部分补偿镜第二面到非球面顶点的距离;d2-部分补偿镜第二面到最接近比较球面球心的距离,通过补偿镜的入射光线坐标和出射光线坐标经代数计算进行确定;d3-最接近比较球面球心到非球面顶点的距离。R0是非球面的顶点曲率半径,K=-e2是二次曲面常数,e为偏心率。|X0.86|是以非球面的顶点V为原点建立直角坐标系时,点A的横坐标绝对值,即点A的弓形高。-K·|X0.86|是非球面的法线像差,d3=R0-K·|X0.86|是法线像差关系式的变形,因此等式d3=R0-K·|X0.86|可以从二次非球面扩展到所有非球面。
最接近比较球面定义如下,依据谢枫的博士学位论文(北京理工大学,2011)“部分补偿非球面检测法的关键问题研究”:“以非球面上某点处法线与最接近比较球面法线的夹角绝对值表示非球面梯度,称为斜率非球面度,并定义使得最大斜率非球面度最小的比较球面为最接近比较球面”。因此,对于一个非球面,能够唯一确定其最接近比较球面。最接近比较球面与该非球面0.86口径处的点是相切的,此处的斜率非球面度最小,如附图2所示。该口径处的曲率半径等于最接近比较球面的半径,沿径向入射的光线可以原路返回。附图2中使用的是凹非球面,但本方法的建模并不限于这一种类型的非球面,也可以使用凸非球面。
根据待测非球面名义参数和部分补偿透镜P的相关参数,确定理想非球面的最佳补偿位置d1=1828.264097mm。
步骤4:加工出部分补偿透镜P的实物,建立包含部分补偿透镜P的实际干涉仪IO,并确定实际非球面的最佳补偿位置,即补偿镜第二面到理想非球面顶点的轴向距离d′1。
根据实际干涉仪IO的干涉图进行定位,当干涉图的条纹最稀疏的时候,能够确定实际非球面的最佳补偿位置d′1=1830.273mm。
步骤5:测量实际非球面与理想非球面的最佳补偿位置之间的离焦距离Δd=d′1-d1=1830.273-1828.264097=2.008903mm。
步骤6:计算实际非球面与理想非球面的最佳补偿位置之间的离焦距离Δd与非球面顶点曲率半径误差ΔR之间的关系
Δd=d′1-d1=(k-K·k(h,R0,K))·ΔR
其中,非球面的口径为2D,则h=0.86D,k(h,R0,K)是通过非球面方程计算出的系数,具体形式为:
步骤7:计算非球面的顶点曲率半径误差
ΔR=Δd/(1-K·k(h,R0,K))=2.0103mm。
综上,相对测量精度为(ΔR-ΔR0)/R0=(2.0103-2)/889≈0.012‰,其中,ΔR0=2mm是被测非球面的实际顶点曲率半径误差。
本发明的设计过程和系统结构简单,在虚拟干涉仪IR中,可以快速的确定理想非球面的最佳补偿位置。在实际干涉仪Io中,对非球面进行非接触、全口径的测量,基于顶点曲率半径误差与最佳补偿位置之间的关系,通过干涉条纹来确定被测非球面的最佳补偿位置,无需进行多位置的测量和记录,无需扫描装置,就可以获得非球面的顶点曲率半径误差,具有较高的相对测量精度高。
Claims (1)
1.基于最佳补偿位置的非球面顶点曲率半径误差测量方法,其特征在于:具体步骤如下:
步骤1:获取待测非球面名义参数;
名义参数包括待测非球面的口径、顶点曲率半径和非球面系数;
步骤2:根据待测非球面名义参数,设计部分补偿透镜,获取其相关参数,用于非球面顶点曲率半径误差的测量;
相关参数包括部分补偿透镜的第一面曲率半径,厚度,材料,第二面曲率半径和口径;
步骤3:在光学仿真软件中建立包含部分补偿透镜的虚拟干涉仪,并确定理想非球面的最佳补偿位置;
根据待测非球面名义参数和部分补偿透镜的相关参数,确定理想非球面的最佳补偿位置,即部分补偿透镜第二面到理想非球面顶点的轴向距离d1:
d1=d2+d3=d2+R0-K·|X0.86| (1)
其中,d1是部分补偿透镜第二面到非球面顶点的距离;d2是部分补偿透镜第二面到最接近比较球面球心的距离,通过部分补偿透镜的入射光线坐标和出射光线坐标经代数计算进行确定;d3是最接近比较球面球心到非球面顶点的距离;R0是非球面的顶点曲率半径,K=-e2是二次曲面常数,e为偏心率;|X0.86|是以非球面的顶点P为原点建立直角坐标系时,点A的横坐标绝对值,即点A的弓形高;
对于(1)式,-K·|X0.86|是非球面的法线像差,d3=R0-K·|X0.86|是法线像差关系式的变形,因此等式d3=R0-K·|X0.86|可以从二次非球面扩展到所有非球面;
步骤4:加工出部分补偿透镜的实物,建立包含部分补偿透镜的实际干涉仪,并确定实际非球面的最佳补偿位置,即部分补偿透镜第二面到理想非球面顶点的轴向距离d′1;
根据实际干涉仪的干涉图进行定位,当干涉图的条纹最稀疏的时候,能够确定实际非球面的最佳补偿位置d′1;
步骤5:测量实际非球面与理想非球面的最佳补偿位置之间的离焦距离Δd:
Δd=d′1-d1 (2)
步骤6:计算实际非球面与理想非球面的最佳补偿位置之间的离焦距离与非球面顶点曲率半径误差之间的关系:
Δd=d′1-d1=(1-K·k(h,R0,K))·ΔR (3)
其中,非球面的口径为2D,则h=0.86D,k(h,R0,K)是通过非球面方程计算出的系数,具体形式为:
步骤7:计算非球面的顶点曲率半径误差:
ΔR=Δd/(1-K·k(h,R0,K)) (5)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710022242.4A CN106871819B (zh) | 2017-01-12 | 2017-01-12 | 基于最佳补偿位置的非球面顶点曲率半径误差测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710022242.4A CN106871819B (zh) | 2017-01-12 | 2017-01-12 | 基于最佳补偿位置的非球面顶点曲率半径误差测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106871819A CN106871819A (zh) | 2017-06-20 |
CN106871819B true CN106871819B (zh) | 2019-03-26 |
Family
ID=59158313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710022242.4A Active CN106871819B (zh) | 2017-01-12 | 2017-01-12 | 基于最佳补偿位置的非球面顶点曲率半径误差测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106871819B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732185B (zh) * | 2018-05-31 | 2020-08-04 | 哈尔滨工业大学 | 一种非球面光学元件表面紫外预处理轨迹的规划方法 |
CN110889200B (zh) * | 2019-11-04 | 2023-06-02 | 天津大学 | 一种非球面玻璃透镜模压成形预补偿方法 |
CN112066913B (zh) * | 2020-08-28 | 2021-09-21 | 四川大学 | 一种测量球面曲率半径的方法 |
CN113900166A (zh) * | 2021-09-01 | 2022-01-07 | 上海现代先进超精密制造中心有限公司 | 一种用于高次非球面检测的补偿器 |
CN116642432B (zh) * | 2023-04-19 | 2024-07-16 | 西安工业大学 | 一种非球面面形检测方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975990A (zh) * | 2010-09-02 | 2011-02-16 | 北京理工大学 | 一种以斜率为优化目标的部分补偿透镜的设计方法 |
CN102288132A (zh) * | 2011-05-18 | 2011-12-21 | 中国科学院长春光学精密机械与物理研究所 | 采用激光跟踪仪测量非球面顶点曲率半径偏差的方法 |
CN102506761A (zh) * | 2011-10-30 | 2012-06-20 | 中国科学院光电技术研究所 | 一种利用激光跟踪仪测量非球面顶点曲率半径的方法 |
CN103234480A (zh) * | 2013-04-16 | 2013-08-07 | 北京理工大学 | 一种环形凸非球面的快速面形检测方法 |
CN105627947A (zh) * | 2015-12-23 | 2016-06-01 | 中国科学院长春光学精密机械与物理研究所 | 一种旋转对称未知非球面面形误差的测量方法及其测量装置 |
-
2017
- 2017-01-12 CN CN201710022242.4A patent/CN106871819B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975990A (zh) * | 2010-09-02 | 2011-02-16 | 北京理工大学 | 一种以斜率为优化目标的部分补偿透镜的设计方法 |
CN102288132A (zh) * | 2011-05-18 | 2011-12-21 | 中国科学院长春光学精密机械与物理研究所 | 采用激光跟踪仪测量非球面顶点曲率半径偏差的方法 |
CN102506761A (zh) * | 2011-10-30 | 2012-06-20 | 中国科学院光电技术研究所 | 一种利用激光跟踪仪测量非球面顶点曲率半径的方法 |
CN103234480A (zh) * | 2013-04-16 | 2013-08-07 | 北京理工大学 | 一种环形凸非球面的快速面形检测方法 |
CN105627947A (zh) * | 2015-12-23 | 2016-06-01 | 中国科学院长春光学精密机械与物理研究所 | 一种旋转对称未知非球面面形误差的测量方法及其测量装置 |
Non-Patent Citations (2)
Title |
---|
基于Zemax的部分补偿透镜的优化设计;孟晓辰等;《光学学报》;20110610(第06期);0622002(1-7) |
基于斜率非球面度的非球面最接近比较球面定义;谢枫等;《光学学报》;20101115(第11期);3197-3202 |
Also Published As
Publication number | Publication date |
---|---|
CN106871819A (zh) | 2017-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106871819B (zh) | 基于最佳补偿位置的非球面顶点曲率半径误差测量方法 | |
CN105423948B (zh) | 采用变形镜的拼接干涉检测非球面面形的装置 | |
CN109029291A (zh) | 结合激光差动共焦定位的非球面参数误差干涉测量方法 | |
CN107796329B (zh) | 一种凸非球面反射镜面形检测装置及检测方法 | |
CN101858735B (zh) | 一种大口径离轴非球面测量和标定系统 | |
CN103926058B (zh) | 非球面检测中采用自准平面镜测量光轴的方法 | |
CN105627947A (zh) | 一种旋转对称未知非球面面形误差的测量方法及其测量装置 | |
CN110487205B (zh) | 结合色散共焦定位的非球面参数误差干涉测量方法 | |
CN102519397B (zh) | 一种光学球面曲率半径的测量方法 | |
CN102506750A (zh) | 部分补偿非球面反射镜面形检测方法 | |
CN101949691A (zh) | 非零位补偿浅度光学非球面面形检测方法 | |
CN104075668B (zh) | 凸双曲面Hindle检测中的高精度几何量测量方法 | |
CN102155926A (zh) | 一种非球面顶点球曲率半径的测量系统及方法 | |
CN107782254A (zh) | 一种混合补偿式子孔径拼接面形检测方法 | |
CN109855560B (zh) | 一种凸非球面反射镜面形的检测装置及检测方法 | |
CN103217125A (zh) | 一种基于子孔径拼接的高精度平面光学元件面型检测方法 | |
CN108895972A (zh) | 一种基于计算全息的光学元件顶点半径测量的方法和装置 | |
CN106595529B (zh) | 基于虚拟牛顿环的大曲率半径非零位干涉测量方法及装置 | |
CN111929037A (zh) | 光楔补偿器标定系统及其标定方法 | |
JP2013186024A (ja) | 非球面形状計測方法、非球面形状計測装置、光学素子加工装置および光学素子 | |
JP2010281792A (ja) | 非球面体測定方法および装置 | |
CN112902875B (zh) | 一种非球面反射镜曲率半径检测装置及方法 | |
Xi et al. | Intensity response model and measurement error compensation method for chromatic confocal probe considering the incident angle | |
CN105352451B (zh) | 一种基于可变形镜的准万能补偿镜及设计方法 | |
CN110966958B (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 |