CN108362202A - 倾斜波面干涉测量非球面中的参数确定方法 - Google Patents

倾斜波面干涉测量非球面中的参数确定方法 Download PDF

Info

Publication number
CN108362202A
CN108362202A CN201810108589.5A CN201810108589A CN108362202A CN 108362202 A CN108362202 A CN 108362202A CN 201810108589 A CN201810108589 A CN 201810108589A CN 108362202 A CN108362202 A CN 108362202A
Authority
CN
China
Prior art keywords
aspherical
node
calculated
center
point
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
CN201810108589.5A
Other languages
English (en)
Other versions
CN108362202B (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.)
Weihai Jiarui Optoelectronic Technology Co ltd
Original Assignee
Suzhou 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 Suzhou University of Science and Technology filed Critical Suzhou University of Science and Technology
Priority to CN201810108589.5A priority Critical patent/CN108362202B/zh
Publication of CN108362202A publication Critical patent/CN108362202A/zh
Application granted granted Critical
Publication of CN108362202B publication Critical patent/CN108362202B/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

Landscapes

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

Abstract

本发明属于光学精密测试领域,具体涉及一种倾斜波面干涉测量非球面中的参数确定方法。包括如下步骤:首先根据给定的非球面方程和口径,按角度等间隔划分非球面;然后设定入射球面波的中心,利用公式计算球面波的中心对于每个节点的镜像点;然后以镜像点向对应的节点坐标位置发射光线构成反射波面,根据计算出的反射波面与参考球面波之间的相位差的分布,判断出可测量干涉条纹的分布范围,移动球面波中心的位置,重复前述过程,直到整个非球面的范围都能够被测量,即确定了所有球面波源的位置,完成了点源阵列的设计。本方法的步骤清晰、计算准确、适用范围广,使用计算机程序可以快速完成设计,计算的过程中无须采用近似估计,计算结果准确。

Description

倾斜波面干涉测量非球面中的参数确定方法
技术领域
本发明属于光学精密测试领域,具体涉及一种倾斜波面干涉测量非球面中的参数确定方法。
背景技术
在旋转对称的光学系统中,使用非球面元件,可以让系统中的元件数量变少,同时能够达到更高的性能。然而,对于非球面的测量则比对球面的测量要困难许多。干涉法是测量光学元件表面面形的一种常见的方法,在球面的测试中,可以加工更高精度的球面,从而实现零位测量,而在非球面的测量中,必须根据待测的非球面方程,加工一块面形相同但精度却更高的非球面才能实现零位测量,显然这种方法不可能实现。因此,需要采取其他测量手段实现检测,例如加工一块计算全息面进行零位检测,使用轮廓仪面形进行扫描、使用结构光投影等方法。近年来,德国斯图加特大学的Osten等人提出了一种多重倾斜波面测量的方法(Eugenio Garbusi, Goran Baer, Wolfgang Osten, Advanced studies on themeasurement of aspheres and freeform surfaces with the Tilted-waveInterferometer),该方法在干涉光源中引入了多个轴外的点源,使其产生多束球面波来补偿被测件不同局部面形。中国发明专利CN103528539A、CN103575229B和CN103759668A也对使用该方法进行自由曲面测量进行讨论。在这一系统的设计中,产生球面波的点源阵列的位置计算是其中的一个关键环节,在非球面的测量中,针对其旋转对称的特点,设计一种简单而精确的点源位置计算方法,充分挖掘该技术的潜力,对于发展非球面的干涉测量技术是大有裨益的。
发明内容
本发明的目的在于提供一种在基于点源阵列非球面干涉仪中,根据仪器的设计指标以及待测的非球面参数,计算点源排列的一种方法。使用该方法可以准确的计算出点源阵列中所有点源的位置,进而顺利完成整个仪器的设计工作。由于非球面都是旋转对称的,因此,只需要以任意一条直径为基础,确定点源的分布情况,再进行旋转对称,就可以设计出二维平面上分布的点源阵列。
本发明的技术方案如下:
倾斜波面干涉测量非球面中的参数确定方法,包括如下步骤:
(1)根据给定的非球面方程和口径,按角度等间隔划分非球面,计算每个节点对应位置的单位切向矢量和单位法向矢量;
(2)设定入射球面波的中心,根据步骤(1)中计算得到的每个节点处的单位切向矢量和单位法向矢量,利用公式计算球面波的中心对于每个节点的镜像点;
(3)以步骤(2)得到的镜像点坐标分别作为点光源,向对应的节点坐标位置发射光线构成反射波面,根据设计指标要求及计算出的反射波面与参考球面波之间的相位差的分布,判断出可测量干涉条纹的分布范围;
(4)移动球面波中心的位置,重复步骤(2)至步骤(3)的过程,直到整个非球面的范围都能够被测量,即确定了所有球面波源的位置,完成了点源阵列的设计。
所述步骤(1)中每个节点对应位置的单位切向矢量和单位法向矢量计算过程如下:
(1a)建立直角坐标系,根据非球面的参数,以非球面上的任意一条直径为截线,给定非球面的参数方程:
,(1);
其中,为非球面的矢量方程;为参数,在此选为轴的夹角;分别为该参数方程在直角坐标系中的坐标;
(1b)设待测的非球面口径为,根据定义的坐标系,的变换范围为,据此可计算得到参数的变化范围是,其中,
,(2);
(1c)根据步骤(1a)中非球面的参数方程,可以计算出非球面的单位法向矢量和单位切向矢量的表达式:
,(3);
,(4);
,(5);
其中,为弧长参数,为单位切向矢量,为单位法向矢量;
(1d)根据步骤(1b)中计算出的,按照其弧度(角度)划分出个等间隔的区段,则相邻两个节点之间的间隔为:
,(6);
那么,对应每个节点的值就是,…,,…,
(1e)将步骤(1d)中计算得到的一系列参数代入步骤(1c)中计算出的单位切矢和单位法矢中,即可得到每个节点对应的切向矢量及法向矢量
设定入射球面波的中心为,所述步骤(2)中该球面波中心对于步骤(1)中得到的每个节点的镜像点的计算过程如下:
(2a)根据公式(1)、(3)、(5),将每个节点的切矢和法矢移动至该节点的坐标处,则有:
,(7);
(2b)定义公式(7)中的,根据下式(8)可以计算出球面波中心的镜像点:
,(8);
其中,表示镜像点的坐标;
(2c)将每个节点的值代入式(7)和(8),就可以得到一组镜像点的坐标组
所述步骤(3)中可测量干涉条纹分布范围的判断过程包括如下步骤:
(3a)以步骤(2c)得到的镜像点坐标组中的每一对坐标分别作为点源,向对应的节点坐标位置发射光线,则每束光线的方程以及入射光线与反射光线的夹角可以表示为:
,(9);
,(10);
(3b)根据设计指标要求的CCD像元尺寸以及干涉条纹分辨率的条件,可以计算出入射光线与出射光线之间允许的最大夹角
,(11);
其中,为激光的波长,表示一根条纹最少需要用几个像素表示;
(3c)将步骤(3a)中计算出的每个节点对应的夹角与步骤(3b)中计算出的最大夹角进行比较,找出以为圆心发出的球面波能够测量的非球面角度范围,再代入非球面方程,确定能够测量的孔径范围
所述公式(11)中的取值为2或3。
本发明与现有技术相比,有益效果在于:本发明提供了一种在基于点源阵列非球面干涉仪中,根据仪器的设计指标以及待测的非球面参数,计算点源排列的一种方法,由于在这种干涉仪中,点源阵列的参数确定是仪器设计的关键环节,直接决定了仪器的整体性能;本发明提供的方法步骤清晰、计算准确、适用范围广,使用计算机程序可以快速完成设计,计算的过程中无须采用近似估计,计算结果准确。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1为本发明的计算流程图。
图2为本发明实施例中,上中下三个波源位置入射和反射以后的光线图,其中1、2、3分别为上中下三个位置的被测椭球面的局部面形;4、5、6分别为上中下三个位置的波源。
图3、图4、图5分别表示本发明实施例中,上中下三个波源位置发出光线经过反射以后,满足干涉条纹采集要求的角度变化范围曲线图;其中,横坐标为非球面方程中的参数,纵坐标为入射光线和反射光线之间的夹角
具体实施方式
下面结合附图及实施例对本发明技术方案作进一步描述。
实施例:
已知待测的非球面的方程为的一个椭圆形,直径为,以圆心为原点建立坐标系,CCD的像元尺寸为,这是一个比较典型的参数,很多工业相机的像元尺寸都能够达到这个值,分辨率较高的相机像元尺寸是上述值的几分之一;测量波长选为,这是常见的氦氖激光器的发射波长;为了保证设计有一定的冗余,每2个像素包含一根条纹,即,根据式(11)计算出的最大允许偏差角。再根据非球面方程计算出椭圆线切向与法向向量,以角度间隔为对曲线进行分割,并计算每个镜像点的出射光线与入射光线之间的夹角,在原点处出射的光线,能够形成干涉的区域为,相应的y轴向范围为,长度为,图4中反映了在此范围内角变化的曲线,图中的最大值,满足上述最大允许偏差的要求。再将点光源的位置沿y轴上移,重复上述计算,得到能够形成干涉的区域为,对应的y轴向范围为,长度为,图3中反映了在此范围内角变化的曲线,图中的最大值,也满足上述最大允许偏差的要求。因为选取的椭球面是沿光轴对称的,所以将点光源的位置沿y轴下移,得到能够形成干涉的区域为,对应的y轴向范围为,长度为,图5中反映了在此范围内角变化的曲线,图中的最大值。因此使用上述三个位置的点光源,可以测量给定椭球面的口径为。至此,就完成了点光源的入射位置计算。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明,本文所定义一般原理可以在不脱离本发明的精神或范围的情况下,在其他实施例中实现。本发明的范围由权利要求及其等同物限定。

Claims (5)

1.倾斜波面干涉测量非球面中的参数确定方法,其特征在于:包括如下步骤:
(1)根据待测的非球面方程和口径,按角度等间隔划分非球面,计算每个节点对应位置的单位切向矢量和单位法向矢量;
(2)设定入射球面波的中心,根据步骤(1)中计算得到的每个节点处的单位切向矢量和单位法向矢量,利用公式计算球面波的中心对于每个节点的镜像点;
(3)以步骤(2)得到的镜像点坐标分别作为点光源,向对应的节点坐标位置发射光线构成反射波面,根据设计指标要求及计算出的反射波面与参考球面波之间的相位差的分布,判断出可测量干涉条纹的分布范围;
(4)移动球面波中心的位置,重复步骤(2)至步骤(3)的过程,直到整个非球面的范围都能够被测量,即确定了所有球面波源的位置,完成了点源阵列的设计。
2.根据权利要求1所述的倾斜波面干涉测量非球面中参数的确定方法,其特征在于:所述步骤(1)中每个节点对应位置的单位切向矢量和单位法向矢量计算过程如下:
(1a)建立直角坐标系,根据非球面的参数,以非球面上的任意一条直径为截线,给定非球面的参数方程:
,(1);
其中,为非球面的矢量方程;为参数,在此选为轴的夹角;分别为该参数方程在直角坐标系中的坐标;
(1b)设待测的非球面口径为,根据定义的坐标系,的变换范围为,据此可计算得到参数的变化范围是,其中,
,(2);
(1c)根据步骤(1a)中非球面的参数方程,可以计算出非球面的单位法向矢量和单位切向矢量的表达式:
,(3);
,(4);
,(5);
其中,为弧长参数,为单位切向矢量,为单位法向矢量;
(1d)根据步骤(1b)中计算出的,按照其弧度(角度)划分出个等间隔的区段,则相邻两个节点之间的间隔为:
,(6);
那么,对应每个节点的值就是,…,,…,
(1e)将步骤(1d)中计算得到的一系列参数代入步骤(1c)中计算出的单位切矢和单位法矢中,即可得到每个节点对应的切向矢量及法向矢量
3.根据权利要求1所述的倾斜波面干涉测量非球面中参数的确定方法,其特征在于:设定入射球面波的中心为,所述步骤(2)中该球面波中心对于步骤(1)中得到的每个节点的镜像点的计算过程如下:
(2a)根据公式(1)、(3)、(5),将每个节点的切矢和法矢移动至该节点的坐标处,则有:
,(7);
(2b)定义公式(7)中的,根据下式(8)可以计算出球面波中心的镜像点:
,(8);
其中,表示镜像点的坐标;
(2c)将每个节点的值代入式(7)和(8),就可以得到一组镜像点的坐标组
4.根据权利要求1所述的倾斜波面干涉测量非球面中参数的确定方法,其特征在于:所述步骤(3)中可测量干涉条纹分布范围的判断过程包括如下步骤:
(3a)以步骤(2c)得到的镜像点坐标组中的每一对坐标分别作为点源,向对应的节点坐标位置发射光线,则每束光线的方程以及入射光线与反射光线的夹角可以表示为:
,(9);
,(10);
(3b)根据设计指标要求的CCD像元尺寸以及干涉条纹分辨率的条件,可以计算出入射光线与出射光线之间允许的最大夹角
,(11);
其中,为激光的波长,表示一根条纹最少需要用几个像素表示;
(3c)将步骤(3a)中计算出的每个节点对应的夹角与步骤(3b)中计算出的最大夹角进行比较,找出以为圆心发出的球面波能够测量的非球面角度范围,再代入非球面方程,确定能够测量的孔径范围
5.根据权利要求4所述的倾斜波面干涉测量非球面中参数的确定方法,其特征在于:所述公式(11)中的取值为2或3。
CN201810108589.5A 2018-02-02 2018-02-02 倾斜波面干涉测量非球面中的参数确定方法 Active CN108362202B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810108589.5A CN108362202B (zh) 2018-02-02 2018-02-02 倾斜波面干涉测量非球面中的参数确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810108589.5A CN108362202B (zh) 2018-02-02 2018-02-02 倾斜波面干涉测量非球面中的参数确定方法

Publications (2)

Publication Number Publication Date
CN108362202A true CN108362202A (zh) 2018-08-03
CN108362202B CN108362202B (zh) 2020-09-25

Family

ID=63004606

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810108589.5A Active CN108362202B (zh) 2018-02-02 2018-02-02 倾斜波面干涉测量非球面中的参数确定方法

Country Status (1)

Country Link
CN (1) CN108362202B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103528539A (zh) * 2013-09-26 2014-01-22 南京理工大学 基于点源阵列的非零位干涉系统
CN103852030A (zh) * 2014-03-17 2014-06-11 南京理工大学 用于倾斜波面非零位干涉系统的自由曲面面形重构方法
US20150185452A1 (en) * 2012-10-25 2015-07-02 Young Optics Inc. Wide-angle projection system
WO2016027874A1 (ja) * 2014-08-21 2016-02-25 公立大学法人大阪市立大学 応力可視化装置および力学物性値可視化装置
CN106197314A (zh) * 2016-07-19 2016-12-07 南京理工大学 一种获得倾斜波面干涉系统中点源阵列分布的规划方法
CN107560565A (zh) * 2017-08-24 2018-01-09 南京理工大学 基于动态分时倾斜载频干涉的面形检测装置及检测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150185452A1 (en) * 2012-10-25 2015-07-02 Young Optics Inc. Wide-angle projection system
CN103528539A (zh) * 2013-09-26 2014-01-22 南京理工大学 基于点源阵列的非零位干涉系统
CN103852030A (zh) * 2014-03-17 2014-06-11 南京理工大学 用于倾斜波面非零位干涉系统的自由曲面面形重构方法
WO2016027874A1 (ja) * 2014-08-21 2016-02-25 公立大学法人大阪市立大学 応力可視化装置および力学物性値可視化装置
CN106197314A (zh) * 2016-07-19 2016-12-07 南京理工大学 一种获得倾斜波面干涉系统中点源阵列分布的规划方法
CN107560565A (zh) * 2017-08-24 2018-01-09 南京理工大学 基于动态分时倾斜载频干涉的面形检测装置及检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
沈华: "基于多重倾斜波面的光学自由曲面非零位干涉测量关键技术研究", 《中国优秀硕士学位论文全文数据库》 *

Also Published As

Publication number Publication date
CN108362202B (zh) 2020-09-25

Similar Documents

Publication Publication Date Title
CN108050933B (zh) 角锥棱镜回射光斑定位精度检测装置与方法
JP4290733B2 (ja) 3次元形状計測方法及びその装置
CN102168955B (zh) 一种光学球面曲率半径的检测方法
US7064817B1 (en) Method to determine and adjust the alignment of the transmitter and receiver fields of view of a LIDAR system
WO2017101557A1 (zh) 面形检测装置及方法
CN109307480B (zh) 一种透射元件多表面面形检测方法
CN109373897B (zh) 一种基于激光虚拟标尺的测量方法
KR101596290B1 (ko) 두께 측정 장치 및 두께 측정 방법
CN108489421A (zh) 一种条纹投影检测平面元件面形方法及装置
CN108036740A (zh) 一种基于多视角的高精度实时三维彩色测量系统及其方法
CN103925890B (zh) 一种基于光束畸变的三维角度测量系统
CN103471561A (zh) 一种三维小角度测量装置及方法
CN103134443B (zh) 一种大口径大径厚比反射镜面形自准直检测装置及方法
CN111412850B (zh) 一种基于单摄像头的高温三维数字图像相关测量系统及方法
JP2016038300A (ja) 光波面センサにおける光スポット重心位置取得方法、光波面計測方法およびプログラム
CN106352985A (zh) 一种非对称空间外差光谱仪结构
CN108362202A (zh) 倾斜波面干涉测量非球面中的参数确定方法
CN107817088B (zh) 离轴抛物镜光轴方向的标定方法及系统
JP2001235317A (ja) 光学球面曲率半径測定装置
CN205940927U (zh) 摆镜特性参数测试装置
CN107941469A (zh) 一种等边三棱镜顶角偏差测量方法
CN108278977B (zh) 测量仪及测量方法
CN113358324A (zh) 一种基于空间相移的散斑干涉烧蚀测量系统及方法
CN112985780A (zh) 一种光学系统倍率色差的测量方法
CN108507488B (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20221111

Address after: 264200 workshop of zaobu Road South and Jinnuo Road West, Gushan Town, economic and Technological Development Zone, Weihai City, Shandong Province

Patentee after: WEIHAI JIARUI OPTOELECTRONIC TECHNOLOGY CO.,LTD.

Address before: 215009 CREE Road, Suzhou hi tech Zone, Suzhou, Jiangsu Province, No. 1

Patentee before: SUZHOU University OF SCIENCE AND TECHNOLOGY