CN110703434B - 一种确定环形孔径二次曲面非球面度梯度的方法 - Google Patents

一种确定环形孔径二次曲面非球面度梯度的方法 Download PDF

Info

Publication number
CN110703434B
CN110703434B CN201910976600.4A CN201910976600A CN110703434B CN 110703434 B CN110703434 B CN 110703434B CN 201910976600 A CN201910976600 A CN 201910976600A CN 110703434 B CN110703434 B CN 110703434B
Authority
CN
China
Prior art keywords
annular aperture
quadric
asphericity
aperture
quadric surface
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
CN201910976600.4A
Other languages
English (en)
Other versions
CN110703434A (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.)
Nantong University
Original Assignee
Nantong University
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 Nantong University filed Critical Nantong University
Priority to CN201910976600.4A priority Critical patent/CN110703434B/zh
Publication of CN110703434A publication Critical patent/CN110703434A/zh
Application granted granted Critical
Publication of CN110703434B publication Critical patent/CN110703434B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/0012Optical design, e.g. procedures, algorithms, optimisation routines

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Lenses (AREA)
  • Grinding And Polishing Of Tertiary Curved Surfaces And Surfaces With Complex Shapes (AREA)

Abstract

本发明公开了一种基于三级像差理论确定环形孔径二次曲面非球面度梯度的方法,全口径二次曲面可以看作是内环孔径为0的环形孔径二次曲面。本发明克服了全口径和环形孔径二次曲面计算非球面度梯度没有函数形式,表达式不统一、不直观的缺点,理论可靠,结果简单。利用非球面参数便可快速、准确、直接地写出各种孔径二次曲面非球面度梯度的表达式,适用范围广,能显著提高工作效率,为非球面的加工和检测方案的制定,以及为工程时间和造价评估提供依据。

Description

一种确定环形孔径二次曲面非球面度梯度的方法
技术领域
本发明属于先进光学制造技术领域,具体涉及一种确定环形孔径二次曲面非球面度梯度的方法。
背景技术
非球面的加工一般都是从最接近比较球面开始的,无论是基于确定性数控技术的先进光学制造还是传统的光学冷加工都需要计算非球面的非球面度、非球面度梯度等参数。数值的大小和分布直接影响到所采用的加工工艺、加工余量大小以及工程造价等。最大非球面度一定程度上可以反映加工的难度,但不能看其绝对值,还与非球面的孔径、非球面系数等参数有关。真正能够反映加工难度的是非球面度的变化值,即非球面度梯度。因此,最接近比较球面、非球面度、非球面度梯度等参数的计算成为光学非球面加工的一个重要环节,为非球面加工工艺和检测方法的确定提供重要依据,也是评估工程时间和造价的重要因素。
环形孔径二次曲面最接近比较球面是指比较球面与环形孔径二次曲面内、外边缘相接触,且比较球面的球心在二次曲面的对称轴上。全口径二次曲面看成内环孔径为0的环形孔径二次曲面。全口径二次曲面的中性带(非球面度取极大值的位置)的位置在。但环形孔径中性带的位置较难确定,目前还没有一个简便的方法可以直接确定。
由于环形孔径(或离轴)非球面的非球面度计算基本限于代数的方法,因此,对环形孔径非球面度梯度的计算更是困难,目前也没发现相关报道。
发明内容
发明目的:为了解决现有技术的不足,本发明提供了一种确定环形孔径二次曲面非球面度梯度的方法。
技术方案:一种确定环形孔径二次曲面非球面度梯度的方法,包括如下步骤:
步骤一、依据物点位于二次曲面顶点曲率中心时的球差表达式建立环形孔径二次曲面的波像差方程;
步骤二、利用环形孔径非球面最接近比较球面特征,确定最接近比较球面曲率中心与非球面顶点曲率中心的偏移量并代入环形孔径二次曲面的波像差表达式,得到仅含有环形孔径二次曲面参数的波像差表达式W(y)、非球面度δ(y),非球面度与波像差之间的关系为:δ(y)=0.5W(y);
步骤三、对非球面度方程求一阶导数,即可得到二次曲面非球面度梯度方程;当dW/dy=0时即可确定环形口径二次曲面在最接近比较球面时中性带的位置,将此位置参数带入非球面度方程式即可得到最大非球面度,全口径非球面可以看做内环孔径为零的环形孔径非球面,环形孔径非球面的最接近比较球面是指比较球面与环形孔径非球面内、外边缘相接触,且比较球面的球心在二次曲面的对称轴上。
环形孔径二次曲面在最接近比较球面时非球面度梯度方程为:
式中k为二次曲面常数,R0二次曲面顶点曲率半径,ɑ为遮拦比,即环形孔径二次曲面内、外环孔径之比;
当遮拦比ɑ=0时,环形孔径二次曲面演变为全口径二次曲面,环形孔径二次曲面的结果适用于全口径非球面,其非球面度梯度表达式为
环形口径二次曲面在最接近比较球面时,
非球面度取极值的位置为
最大非球面度为
有益效果:本发明提供了一种确定环形孔径二次非球面在最接近比较球面时最大非球面度的方法,可直观、快速地确定环形孔径二次非球面在最接近比较球面时非球面度以及最大非球面度,理论可靠,判断准确,适用范围广,节省时间,能够提高工作效率。
附图说明
下面结合附图和实例对本发明作进一步说明。
图1为环形孔径二次曲面及其最接近比较球面示意图。
图2为环形孔径二次曲面非球面度梯度曲线。
图3为全口径二次曲面非球面度梯度曲线。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,以使本领域的技术人员能够更好的理解本发明的优点和特征,从而对本发明的保护范围做出更为清楚的界定。本发明所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例
环形孔径二次曲面最接近比较球面是指比较球面与环形孔径二次曲面内、外边缘相接触,且比较球面的球心在二次曲面的对称轴上,如图1所示。图中曲线EOF为二次曲面的子午截线,曲线MAN为二次曲面最接近比较球面的子午截线,Ox为对称轴。O为二次曲面子午截线的顶点与对称轴的交点,A为最接近比较球面的子午截线与对称轴的交点,B二次曲面子午截线顶点曲率半径的圆心,C为最接近比较球面子午截线的圆心,h1=D1/2,h=D/2,D1,D分别为非球面内、外环的孔径。
一般情况下非球面采用其子午截线的方程表示即可,方程式为:
式中,c为非球面顶点曲率,k为二次曲面常数,a1,a2为高次项系数,y为纵坐标,表示非球面的半口径坐标,x为横坐标,表示非球面的矢高。
实际应用中,在求解非球面度时,通常忽略非球面的高次项,仅按二次曲面进行计算。
基于波像差理论,参考专利(CN201410324768.4“一种判断旋转轴对称非球面能否采用直接干涉检测的方法”),对于环形孔径非球面的波像差表达式,可写为,
当内边缘和外边缘波像差相等时可知:
所以此时的离焦量Δ为:
此时的波像差方程为,
为了书写方便,令h1=αh,α称为遮拦比。方程式(5)改写为:
由δ(y)=0.5W(y)可知:
所以环形孔径二次曲面非球面度梯度表达式为:
当dδ/dy=0时,非球面度取极值,求解后可得:
方程式(9)即为非球面度取得极值的位置,将其带入方程式(7)可知非球面度的极值δmax为,
当ɑ=0时,即为全口径情况下在最接近比较球面时,非球面度取得极大值的位置以及最大非球面度,
这个结果与全口径非球面得出的结果是一致的。
二次曲面既可以是凹面,也可以是凸面;当二次曲面为凹面时R0为负值,当二次非球面为凸面时R0为正值。
具体实施例1
某环形凹非球面参数为:R0=-1440毫米,k=-1.00486,内环口径D1=2h1=300毫米,外环口径D2=2h2=600毫米。将参数带入方程式(8),即可快速得到非球面度梯度表达式为:,据此可以画出非球面度梯度曲线如图2所示。非球面度取极值的位置为:/>,即距非球面中心237.17mm的圆上。最大非球面度为:
具体实施例2
某凹非球面参数为:R0=-1440毫米,k=-1.00486,D=600毫米。将参数带入方程式(11),即可快速得到非球面度梯度表达式为:,据此可以画出非球面度梯度曲线如图2所示。非球面度取极值的位置为:,即距非球面中心212.1mm的圆上。最大非球面度为:
因此,从计算过程看,本发明不仅首次给出了二次曲面在最接近比较球面时非球面度梯度的函数表达式,并且统一了全口径和环形口径非球面非球面度梯度计算方法。表达式中各参数仅与非球面自身参数有关;从计算结果看,本发明得到的结果与传统计算得到的结果一致,说明本发明的计算方法是准确的、可靠的。

Claims (3)

1.一种确定环形孔径二次曲面非球面度梯度的方法,其特征在于:包括如下步骤:
步骤一、依据物点位于二次曲面顶点曲率中心时的球差表达式建立环形孔径二次曲面的波像差方程;
步骤二、利用环形孔径非球面最接近比较球面特征,确定最接近比较球面曲率中心与非球面顶点曲率中心的偏移量并代入环形孔径二次曲面的波像差表达式,得到仅含有环形孔径二次曲面参数的波像差表达式W(y)、非球面度δ(y),非球面度与波像差之间的关系为:δ(y)=0.5W(y);
步骤三、对非球面度方程求一阶导数,即可得到二次曲面非球面度梯度方程;当dW/dy=0时即可确定环形口径二次曲面在最接近比较球面时中性带的位置,将此位置参数带入非球面度方程式即可得到最大非球面度,全口径非球面可以看做内环孔径为零的环形孔径非球面,环形孔径非球面的最接近比较球面是指比较球面与环形孔径非球面内、外边缘相接触,且比较球面的球心在二次曲面的对称轴上。
2.根据权利要求1所述的确定环形孔径二次曲面非球面度梯度的方法,其特征在于:环形孔径二次曲面在最接近比较球面时非球面度梯度方程为:
式中k为二次曲面常数,R0二次曲面顶点曲率半径,ɑ为遮拦比,即环形孔径二次曲面内、外环孔径之比;
当遮拦比ɑ=0时,环形孔径二次曲面演变为全口径二次曲面,环形孔径二次曲面的结果适用于全口径非球面,其非球面度梯度表达式为
3.根据权利要求1所述的确定环形孔径二次曲面非球面度梯度的方法,其特征在于:环形口径二次曲面在最接近比较球面时,
非球面度取极值的位置为
最大非球面度为
CN201910976600.4A 2019-10-15 2019-10-15 一种确定环形孔径二次曲面非球面度梯度的方法 Active CN110703434B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910976600.4A CN110703434B (zh) 2019-10-15 2019-10-15 一种确定环形孔径二次曲面非球面度梯度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910976600.4A CN110703434B (zh) 2019-10-15 2019-10-15 一种确定环形孔径二次曲面非球面度梯度的方法

Publications (2)

Publication Number Publication Date
CN110703434A CN110703434A (zh) 2020-01-17
CN110703434B true CN110703434B (zh) 2024-04-12

Family

ID=69198607

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910976600.4A Active CN110703434B (zh) 2019-10-15 2019-10-15 一种确定环形孔径二次曲面非球面度梯度的方法

Country Status (1)

Country Link
CN (1) CN110703434B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102937421A (zh) * 2012-10-29 2013-02-20 南通大学 旋转轴对称光学非球面的实时检测方法
CN103017681A (zh) * 2012-12-07 2013-04-03 南通大学 近抛物面的旋转轴对称凹非球面的实时检测方法
CN103196391A (zh) * 2013-04-16 2013-07-10 北京理工大学 一种近抛物面的环形凹非球面的快速面形检测方法
CN104048619A (zh) * 2014-07-09 2014-09-17 南通大学 一种判断旋转轴对称非球面能否采用直接干涉检测的方法
CN106289107A (zh) * 2016-09-29 2017-01-04 西安交通大学 一种lc‑slm误差补偿方法及其非球面面型检测方法
CN108919486A (zh) * 2018-08-20 2018-11-30 南京理工大学 低温真空环境下光学窗口的热光学特性表征与优化方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102937421A (zh) * 2012-10-29 2013-02-20 南通大学 旋转轴对称光学非球面的实时检测方法
CN103017681A (zh) * 2012-12-07 2013-04-03 南通大学 近抛物面的旋转轴对称凹非球面的实时检测方法
CN103196391A (zh) * 2013-04-16 2013-07-10 北京理工大学 一种近抛物面的环形凹非球面的快速面形检测方法
CN104048619A (zh) * 2014-07-09 2014-09-17 南通大学 一种判断旋转轴对称非球面能否采用直接干涉检测的方法
CN106289107A (zh) * 2016-09-29 2017-01-04 西安交通大学 一种lc‑slm误差补偿方法及其非球面面型检测方法
CN108919486A (zh) * 2018-08-20 2018-11-30 南京理工大学 低温真空环境下光学窗口的热光学特性表征与优化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
任志君 等.非球面光学零件工艺.《光学零件制作工艺学》.2019, *

Also Published As

Publication number Publication date
CN110703434A (zh) 2020-01-17

Similar Documents

Publication Publication Date Title
CA2467170A1 (en) Method and apparatus relating to the optical zone of an optical element
CN103234480A (zh) 一种环形凸非球面的快速面形检测方法
CN110703434B (zh) 一种确定环形孔径二次曲面非球面度梯度的方法
CN111240010B (zh) 一种用于自由曲面测量的可变形镜面形设计方法及装置
CN103809290A (zh) 一种光学系统面形误差相互补偿优化方法
CN108446423B (zh) 一种光学元件面形加工的工艺与参数选择及其应用方法
CN108227186B (zh) 确定光学系统环形口径二次曲面最接近比较球面曲率半径的方法
CN110703433A (zh) 计算环形孔径二次曲面最接近比较球面曲率半径的方法
Moroni et al. Performance improvement for optimization of the non-linear geometric fitting problem in manufacturing metrology
CN111307102B (zh) 一种光学加工中获得离轴非球面面型特征的方法
CN110596893B (zh) 一种基于像差理论确定二次曲面非球面度的方法
CN112163294A (zh) 一种圆柱度误差快速评定方法
CN107300750A (zh) 光学成像镜头
CN106468817A (zh) 光学成像镜头
CN108873305A (zh) 一种大视场两反式Golay3稀疏孔径望远镜的设计方法
CN111880305B (zh) 基于地图投影的大型光学望远镜非球面主镜划分的方法
CN204556980U (zh) 矫正像差眼镜片
CN103792660B (zh) 一种小视场光学系统面形误差相互补偿的自动优化方法
CN104375281A (zh) 一种眼镜片矫正像差的方法及其眼镜片
CN207408657U (zh) 光学无热化镜头
CN207408658U (zh) 大视角光学无热化镜头
CN221037155U (zh) 一种法兰rj环连接槽加工专用试块
CN103760684A (zh) 自由曲面镜片质量评价的方法
CN109471244A (zh) 一种基于机器视觉的大视场双远心镜头
Fang Research on Evaluation Algorithm of Roundness Measurement Data Based on Adaptive Genetic Algorithm

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