CN110703433A - 计算环形孔径二次曲面最接近比较球面曲率半径的方法 - Google Patents

计算环形孔径二次曲面最接近比较球面曲率半径的方法 Download PDF

Info

Publication number
CN110703433A
CN110703433A CN201910976599.5A CN201910976599A CN110703433A CN 110703433 A CN110703433 A CN 110703433A CN 201910976599 A CN201910976599 A CN 201910976599A CN 110703433 A CN110703433 A CN 110703433A
Authority
CN
China
Prior art keywords
quadric
aperture
closest
spherical
curvature
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
Application number
CN201910976599.5A
Other languages
English (en)
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 CN201910976599.5A priority Critical patent/CN110703433A/zh
Publication of CN110703433A publication Critical patent/CN110703433A/zh
Pending legal-status Critical Current

Links

Images

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)
  • Testing Of Optical Devices Or Fibers (AREA)

Abstract

本发明公开了一种计算环形孔径二次曲面最接近比较球面曲率半径的方法。全口径二次曲面看成内环孔径为0的环形孔径二次曲面。本发明克服了全口径和环形孔径二次曲面计算最接近比较球面方法不统一、表达式不直观的缺点。本方法利用非球面参数便可快速、准确、直接地写出各种孔径二次曲面最接近比较球面曲率半径,适用范围广,能显著提高工作效率。

Description

计算环形孔径二次曲面最接近比较球面曲率半径的方法
技术领域
本发明属于先进光学制造技术领域,具体涉及一种计算环形孔径二次曲面最接近比较球面曲率半径的方法。
背景技术
由于非球面比球面可以提供更大的自由度,所以在光学系统中使用非球面可以缩小系统整体外形尺寸,减轻系统重量,改善成像质量等优点。反射式光学系统在天文观测中的应用已十分广泛,由于镜面材料在光学性能上没有特殊的要求,且没有色差问题,因此,它与折射系统相比,可以使用大口径材料,可以轻量化等优点。在反射式光学系统中,使用非球面能够大大改善像质。但是对反射镜而言一般存在中心遮拦,因此反射式光学系统中的反射镜一般是环形孔径或者离轴非球面镜,离轴非球面可以看做是环形孔径非球面的一部分。
通常在加工非球面前,总是先将材料加工成球面,然后在球面的基础上加工出非球面,这个先加工出的球面就称为比较球面。全口径比较球面的选择一般有3种情况:(1)比较球面与非球面在顶点处相切,比较球面曲率半径等于非球面顶点曲率半径,(2)比较球面与非球面在口径边缘处相切,(3)比较球面与非球面在顶点处相切且在口径边缘处相接。
比较球面的选择要要有利于非球面的加工。一般的判断方法是根据非球面度。非球面度是指比较球面与非球面在沿着光轴或法线的方向的偏差。对于同一非球面,当最大非球面度较小时,有利于非球面的加工。对于全口径非球面,第(3)种情况与(1)和(2)相比,最大非球面度最小,因此,通常称第(3)种情况的比较球面为非球面的最接近比较球面。
依据最接近比较球面的特点可知:环形孔径二次曲面最接近比较球面与环形孔径二次曲面内、外边缘相接触,且比较球面的球心在二次曲面的对称轴上。
一般情况下非球面采用其子午截线的方程表示即可,方程式为:
Figure BDA0002233838700000021
式中,c为非球面顶点曲率,k为二次曲面常数,a1,a2为高次项系数,y为纵坐标,表示非球面的半口径坐标,x为横坐标,表示非球面的矢高。
实际应用中,在求解比较球面时,通常忽略非球面的高次项,仅按二次曲面进行计算。
通常在计算二次曲面最接近比较球面时,根据孔径类型把非球面分为全口径、环形孔径以及离轴非球面等,后再利用各自相应的方法计算出最接近比较球面的曲率半径。
全口径二次曲面最接近比较球面曲率半径的求解方法一般过程如下:
首先利用几何的方法,得出最接近比较球面曲率半径表达式为
LA′=-0.25·D2k/R0 (2)
式中R0为非球面顶点曲率半径,x为全口径时非球面矢高。
然后再根据方程式(1)求取x,进而计算出最接近比较球面的曲率半径。
在专利(CN201410324768.4“一种判断旋转轴对称非球面能否采用直接干涉检测的方法”)中还公布了一种利用球差计算全口径二次曲面最接近比较球面的方法:R=R0+LA’/4。LA’为物体位于二次曲面顶点曲率中心时,物体经二次曲面成像的球差值LA′=-0.25·D2k/R0。式中R为全口径二次曲面最接近比较球面的曲率半径,R0为二次曲面顶点曲率半径,D为二次曲面的口径,k为二次曲面常数。
环形孔径二次曲面面最接近比较球面曲率半径的求解主要还是依靠几何与代数的方法,其计算过程如下:
首先利用代数的方法分别将y1=h1=D1/2,y2=h2=D2/2代入方程式(1),计算出x1与x2的值。假设最接近比较球面曲率半径为R,球心坐标为(a,0),利用距离公式可以计算出:
Figure BDA0002233838700000031
将a,x1,y1或a,x2,y2代入下式,从而计算出环形孔径最接近比较球面的曲率半径R值:
在专利(CN201810038961.X“确定光学系统环形口径二次曲面最接近比较球面曲率半径的方法”)中也公布了一种利用球差计算环形孔径二次曲面最接近比较球面的方法,其计算方法如下:
首先计算物体位于二次曲面顶点曲率中心时,物体分别经二次曲面内环和外环口径成像的球差值LA1',LA2'分别为,
Figure BDA0002233838700000033
Figure BDA0002233838700000034
式中D1,D2为两个非球面的孔径,k为二次曲面常数。
利用式(6)计算环形口径二次曲面最接近比较球面球心距离非球面顶点的位置坐标(a,0)。
a=R0+(LA1'+LA2')/4 (7)
将二次曲面参数及y1=D1/2代入方程式(1),计算出x1,因此(x1,y1)即为二次曲面子午截线上的坐标
距离公式为
Figure BDA0002233838700000041
即可算出环形口径二次曲面最接近比较球面的曲率半径。
以上几种方法确实可以计算出各类二次曲面的最接近比较球面,有的计算较为繁琐,有的计算也比较简单。但是不管那一种方法,对于二次曲面的最接近比较球面曲率半径的计算都没有形成统一的表达式,结果也不够直观。
发明内容
发明目的:为了解决现有技术的不足,本发明提供了一种计算环形孔径二次曲面最接近比较球面曲率半径的方法。
技术方案:一种计算环形孔径二次曲面最接近比较球面曲率半径的方法,包括如下步骤:
步骤一、将全口径二次曲面看成内环孔径为0的环形孔径二次曲面;
步骤二、确定与环形孔径二次曲面内、外边缘相接触,且球心在二次曲面对称轴上的比较球面为最接近比较球面;
步骤三、建立最接近比较球面曲率半径的与二次曲面顶点曲率半径、二次曲面顶点曲率中心与最近接比较球面球心的离焦量以及二次曲面与最接近比较球面沿对称轴方向的偏离量几何关系;
步骤四、基于波像差理论计算最近接比较球面球心与二次曲面顶点曲率中心的离焦量以及最近接比较球面与二次曲面顶点沿光轴方向的偏离之间偏离量;
步骤五、利用几何关系即可得到二次曲面最接近比较球面曲率半径R的通用计算公式。
作为优化:二次曲面最接近比较球面曲率半径R的通用计算公式为:
Figure BDA0002233838700000051
式中R0为二次曲面顶点曲率半径,D1,D2分别为非球面内、外环的孔径,k为二次曲面常数;
当D1=0时,通用计算公式变为
Figure BDA0002233838700000052
此时为全口径二次曲面最接近比较球面曲率半径的计算公式。
有益效果:本发明公布了一种计算环形孔径二次曲面最接近比较球面曲率半径的方法,表达式仅与二次曲面自身参数有关,各项物理意义清晰。利用非球面参数可直观、快速地计算出环形孔径二次曲面最接近比较球面曲率半径,适用范围广,能够提高工作效率。
附图说明
下面结合附图和实例对本发明作进一步说明。
图1环形口径二次曲面示意图。
图2离轴二次曲面示意图。
图3为全口径二次曲面示意图。
图4环形孔径非球面及其最接近比较球面示意图。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,以使本领域的技术人员能够更好的理解本发明的优点和特征,从而对本发明的保护范围做出更为清楚的界定。本发明所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例
二次曲面一般分为全口径二次曲面、环形孔径二次曲面以及离轴二次曲面。环形孔径二次曲面如图1所示;离轴二次曲面可以看作是环形孔径二次曲面的一部分,如图2所示;全口径二次曲面可以看作是内环孔径为0的环形孔径二次曲面,如图3所示。这样就可以将各种孔径的二次曲面统一看作为环形孔径二次曲面。
环形孔径二次曲面最接近比较球面是指比较球面与环形孔径二次曲面内、外边缘相接触,且比较球面的球心在二次曲面的对称轴上,如图4所示。图中曲线EOF为二次曲面的子午截线,曲线MAN为二次曲面最接近比较球面的子午截线,Ox为对称轴。O为二次曲面子午截线的顶点与对称轴的交点,A为最接近比较球面的子午截线与对称轴的交点,B二次曲面子午截线顶点曲率半径的圆心,C为最接近比较球面子午截线的圆心,h1=D1/2,h2=D1/2,D1,D2分别为非球面内、外环的孔径。
因此,从图1中可知:
Figure BDA0002233838700000061
式中为最接近比较球面的曲率半径,
Figure BDA0002233838700000063
为二次曲面顶点曲率半径,
Figure BDA0002233838700000064
为二次曲面与最接近比较球面沿对称轴方向的偏离量,
Figure BDA0002233838700000065
为二次曲面顶点曲率中心与最近接比较球面球心的离焦量。
基于波像差理论,参考专利(CN201410324768.4“一种判断旋转轴对称非球面能否采用直接干涉检测的方法”),对于环形孔径非球面的波像差表达式,可写为,
Figure BDA0002233838700000066
当内边缘和外边缘波像差相等时可知:
Figure BDA0002233838700000071
所以此时的离焦量Δ,即二次曲面顶点曲率中心与最近接比较球面球心的偏移量由方程式(10)可得:
Figure BDA0002233838700000073
此时的波像差方程为,
Figure BDA0002233838700000074
由δ(y)=0.5W(y),δ(y)为非球面度,是非球面与比较球面的偏离量,可知:
Figure BDA0002233838700000075
当y=0时,δ为二次曲面与最接近比较球面沿对称轴方向的偏离量,即
Figure BDA0002233838700000076
由方程式(13)可知:
Figure BDA0002233838700000077
将方程式(11)和(14)带入方程式(8),则可以得到二次曲面最接近比较球面曲率半径R的通用计算公式为:
Figure BDA0002233838700000078
式中R0为二次曲面顶点曲率半径,k为二次曲面常数。
当D1=0时,通用计算公式(15)改写为
Figure BDA0002233838700000079
此时为全口径二次曲面最接近比较球面曲率半径的计算公式,该表达式与专利(CN201410324768.4“一种判断旋转轴对称非球面能否采用直接干涉检测的方法”)的结果一致。
二次曲面既可以是凹面,也可以是凸面;当二次曲面为凹面时R0为负值,当二次曲面为凸面时R0为正值。
应用举例1
某环形凹非球面参数为:R0=-1440毫米,k=-1.00486,内环口径D1=2h1=300毫米,外环口径D2=2h2=600毫米。
利用本发明计算最接近比较球面曲率半径过程如下:
发明专利(CN201810038961.X“确定光学系统环形口径二次曲面最接近比较球面曲率半径的方法”)计算过程如下:
Figure BDA0002233838700000082
Figure BDA0002233838700000083
Figure BDA0002233838700000084
利用传统方法计算如下:
Figure BDA0002233838700000086
Figure BDA0002233838700000087
应用举例2
某凹非球面参数为:R0=-1440毫米,k=-1.00486,D=600毫米。
利用本发明计算最接近比较球面曲率半径过程如下:
Figure BDA0002233838700000091
利用传统方法计算如下:
Figure BDA0002233838700000093
因此,从计算过程看,本发明一种计算环形孔径二次曲面最接近比较球面曲率半径的方法,并包含了全口径二次曲面最接近比较球面的计算方法。表达式中各参数仅与非球面自身参数有关,简化了计算过程,且表达式中各项物理意义清晰,有利于理解最接近比较球面的物理意义;从计算结果看,本发明得到的结果与传统计算得到的结果一致,说明本发明的计算方法是准确的、可靠的。

Claims (2)

1.一种计算环形孔径二次曲面最接近比较球面曲率半径的方法,其特征在于:包括如下步骤:
步骤一、将全口径二次曲面看成内环孔径为0的环形孔径二次曲面;
步骤二、确定与环形孔径二次曲面内、外边缘相接触,且球心在二次曲面对称轴上的比较球面为最接近比较球面;
步骤三、建立最接近比较球面曲率半径的与二次曲面顶点曲率半径、二次曲面顶点曲率中心与最近接比较球面球心的离焦量以及二次曲面与最接近比较球面沿对称轴方向的偏离量几何关系;
步骤四、基于波像差理论计算最近接比较球面球心与二次曲面顶点曲率中心的离焦量以及最近接比较球面与二次曲面顶点沿光轴方向的偏离之间偏离量;
步骤五、利用几何关系即可得到二次曲面最接近比较球面曲率半径R的通用计算公式。
2.根据权利要求1所述的计算环形孔径二次曲面最接近比较球面曲率半径的方法,其特征在于:二次曲面最接近比较球面曲率半径R的通用计算公式为:
Figure FDA0002233838690000011
式中R0为二次曲面顶点曲率半径,D1,D2分别为非球面内、外环的孔径,k为二次曲面常数;
当D1=0时,通用计算公式变为
Figure FDA0002233838690000012
此时为全口径二次曲面最接近比较球面曲率半径的计算公式。
CN201910976599.5A 2019-10-15 2019-10-15 计算环形孔径二次曲面最接近比较球面曲率半径的方法 Pending CN110703433A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910976599.5A CN110703433A (zh) 2019-10-15 2019-10-15 计算环形孔径二次曲面最接近比较球面曲率半径的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910976599.5A CN110703433A (zh) 2019-10-15 2019-10-15 计算环形孔径二次曲面最接近比较球面曲率半径的方法

Publications (1)

Publication Number Publication Date
CN110703433A true CN110703433A (zh) 2020-01-17

Family

ID=69198355

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910976599.5A Pending CN110703433A (zh) 2019-10-15 2019-10-15 计算环形孔径二次曲面最接近比较球面曲率半径的方法

Country Status (1)

Country Link
CN (1) CN110703433A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114791668A (zh) * 2022-04-26 2022-07-26 南通大学杏林学院 一种计算二次非球面最佳拟合球面曲率半径的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102997863A (zh) * 2012-11-05 2013-03-27 北京理工大学 一种全口径光学非球面面形误差直接检测系统
CN104048619A (zh) * 2014-07-09 2014-09-17 南通大学 一种判断旋转轴对称非球面能否采用直接干涉检测的方法
CN104678554A (zh) * 2015-03-12 2015-06-03 南通大学 一种计算全口径二次曲面最接近比较球面曲率半径的方法
CN108227186A (zh) * 2018-01-16 2018-06-29 南通大学 确定光学系统环形口径二次曲面最接近比较球面曲率半径的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102997863A (zh) * 2012-11-05 2013-03-27 北京理工大学 一种全口径光学非球面面形误差直接检测系统
CN104048619A (zh) * 2014-07-09 2014-09-17 南通大学 一种判断旋转轴对称非球面能否采用直接干涉检测的方法
CN104678554A (zh) * 2015-03-12 2015-06-03 南通大学 一种计算全口径二次曲面最接近比较球面曲率半径的方法
CN108227186A (zh) * 2018-01-16 2018-06-29 南通大学 确定光学系统环形口径二次曲面最接近比较球面曲率半径的方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114791668A (zh) * 2022-04-26 2022-07-26 南通大学杏林学院 一种计算二次非球面最佳拟合球面曲率半径的方法

Similar Documents

Publication Publication Date Title
CN109029291B (zh) 结合激光差动共焦定位的非球面参数误差干涉测量方法
CN109739019B (zh) 基于矢量像差理论的共轴自由曲面光学系统优化设计方法
CN103809290B (zh) 一种光学系统面形误差相互补偿优化方法
CN111240010B (zh) 一种用于自由曲面测量的可变形镜面形设计方法及装置
US20220244568A1 (en) Ophthalmic lens
CN103234480A (zh) 一种环形凸非球面的快速面形检测方法
CN110703433A (zh) 计算环形孔径二次曲面最接近比较球面曲率半径的方法
CN108227186B (zh) 确定光学系统环形口径二次曲面最接近比较球面曲率半径的方法
CN110727092A (zh) 一种基于自由曲面的离轴反射式两镜扩束系统
US5235357A (en) Spectacle lenses having an inside-out conicoidal front surface and method for forming same
US11294148B2 (en) Method for designing freeform surface imaging systems
JP4573915B2 (ja) 2つの非球面を有する強化眼鏡レンズ
CN108873305B (zh) 一种大视场两反式Golay3稀疏孔径望远镜的设计方法
CN110596893B (zh) 一种基于像差理论确定二次曲面非球面度的方法
WO2024007370A1 (zh) 一种折反射式不动型反摄远变焦系统设计方法
CN207571372U (zh) 监控镜头
CN110703434B (zh) 一种确定环形孔径二次曲面非球面度梯度的方法
CN115576102A (zh) 一种面向制造性能最优的光学优化方法
US9857550B2 (en) Method for reducing the effects of element misalignment in an optical system
CN109212734B (zh) 一种基于自由曲面的稀疏孔径两反光学成像系统
KR20100111497A (ko) 대물렌즈
CN111880305A (zh) 基于地图投影的大型光学望远镜非球面主镜划分的方法
WO2023088588A1 (en) A method for determining an ophthalmic lens adapted to slow down the progression of a vision impairment
CN103792660A (zh) 一种小视场光学系统面形误差相互补偿的自动优化方法
JP2020003765A (ja) 光学系の設計方法、及び、光学設計ソフトウエア

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20200117

RJ01 Rejection of invention patent application after publication