CN102980532A - 采用三坐标测量仪拼接测量大口径非球面面形的方法 - Google Patents

采用三坐标测量仪拼接测量大口径非球面面形的方法 Download PDF

Info

Publication number
CN102980532A
CN102980532A CN2012105706445A CN201210570644A CN102980532A CN 102980532 A CN102980532 A CN 102980532A CN 2012105706445 A CN2012105706445 A CN 2012105706445A CN 201210570644 A CN201210570644 A CN 201210570644A CN 102980532 A CN102980532 A CN 102980532A
Authority
CN
China
Prior art keywords
sub
aperture
measuring machine
sigma
coordinates measuring
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
CN2012105706445A
Other languages
English (en)
Other versions
CN102980532B (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.)
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Original Assignee
Changchun Institute of Optics Fine Mechanics and Physics of CAS
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 Changchun Institute of Optics Fine Mechanics and Physics of CAS filed Critical Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority to CN201210570644.5A priority Critical patent/CN102980532B/zh
Publication of CN102980532A publication Critical patent/CN102980532A/zh
Application granted granted Critical
Publication of CN102980532B publication Critical patent/CN102980532B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

采用三坐标测量仪拼接测量大口径非球面面形的方法,涉及一种测量大口径非球面面形的方法,解决现有对大口径非球面的检测无法采用三坐标测量仪实现检测的问题;对被检测大口径非球面的尺寸进行划分,相邻子孔径有四分之一的区域为重合区域;采用三坐标测量仪测量第一子孔径区域的数据,并在重叠区域贴上三个靶标,采用三坐标测量仪测定靶标的数据;调整被检测大口径非球面的位置,采用三坐标测量仪再次测定三个靶标的数据;移去三个靶标,获得第二子孔径区域的面形数据;采用迭代算法进行两孔径的基准统一,求解拼接因子,如果有多个子孔径,采用两两拼接的方法实现对大口径非球面面形的检测。本方法操作简单,检测成本低。

Description

采用三坐标测量仪拼接测量大口径非球面面形的方法
技术领域
本发明涉及一种测量非球面面形尤其是大口径非球面面形的方法。
背景技术
对于中低精度的非球面元件和高精度光学非球面加工过程的检验,利用三坐标测量仪对其进行轮廓测量,一直是一种最常用和实用的方法。 该方法通过测量非球面表面相对某一测量基准的绝对矢高,然后通过计算机软件与理论值相比较,获得绝对矢高相对于理论矢高的偏差,通过计算机软件进行误差分析、数据拟合等计算可以获得非球面面形误差分布。
但是该方法受三坐标测量仪测试范围的限制,对于中小口径的非球面,其外形尺寸一般在三坐标测量仪的检测范围以内,能够直接完成对其面形的测量,但对于尺寸超过三坐标测量仪的大口径非球面的测量,三坐标测量仪却无法实现。
发明内容
本发明为解决现有对大口径非球面的检测无法采用三坐标测量仪实现检测的问题,提供一种采用坐标测量仪拼接测量大口径非球面面形的方法。
采用三坐标测量仪拼接测量大口径非球面面形的方法,该方法由以下步骤实现:
步骤一、对被检测大口径非球面的尺寸进行划分,获得多个子孔径;所述多个子孔径在三坐标测量仪的检测范围内,相邻子孔径有四分之一的区域为重合区域;
步骤二、采用三坐标测量仪检测并采集子孔径区域的面形数据,并在该子孔径区域与相邻的另一个子孔径区域的重叠区域贴上三个靶标,采用三坐标测量仪测定三个靶标的数据;
步骤三、调整被检测大口径非球面的位置,使三坐标测量仪能够检测到另一个子孔径区域的面形数据;采用三坐标测量仪测定三个靶标的数据;移去三个靶标,三坐标测量仪获得另一个子孔径区域的面形数据;
步骤四、判断是否是最后一个子孔径,如果是,则执行步骤五;如果否,返回执行步骤二;
步骤五、采用迭代算法将步骤二和步骤三获得的子孔径区域的面形数据进行坐标变换,获得在相同基准面上的两个子孔径数据;
步骤六、判断步骤五获得的两个子孔径的重叠区域是否大于四分之一,如果否,则返回执行步骤三;如果是,则采用最小二乘法获得两个子孔径的拼接系数,最终实现对大口径非球面面形的检测。
本发明的有益效果:本发明扩大了三坐标测量仪的测试范围,能够准确的实现对大口径非球面面形的测量。该方法物理概念明确、数据处理简单、测试效率很高。利用本发明方法能够精确实现对大口径非球面面形的测量,该方法操作简便、数据分析和运算简单、测量时间短、检测成本低。
附图说明
图1为本发明所述的三坐标测量仪拼接测量大口径非球面面形方法的流程图;
图2为本发明所述的三坐标测量仪拼接测量大口径非球面面形方法中子孔径拼接示意图。
具体实施方式
具体实施方式一、结合图1和图2说明本实施方式,本实施方式所述的三坐标测量仪拼接测量大口径非球面面形方法,该方法应用的装置主要包括三坐标测量仪、被检测大口径非球面、ViScan光学探头和镜面靶标等;通过三坐标测量仪对非球面各子孔径区域进行面形测量,对各子孔径数据进行分析和求解即可得到大口径非球面全口径的面形信息。
具体步骤为:
一、建立全局坐标系;建立三坐标测量仪的全局坐标系,作为大口径非球面镜体拼接检测的统一基准。以三坐标测量基准平面上任意相互垂直的两条直线分别定义为x轴和y轴,以与基准平面垂直的直线定义为z轴,坐标系如图2中所示。
二、子孔径规划:根据三坐标测量仪的行程和大口径非球面的尺寸划分子孔径的大小并计算子孔径的个数。划分的原则为:各子孔径在三坐标测量仪的检测范围内,子孔径的个数为最小,但各子孔径区域有1/4以上的重叠区域4。
三、子孔径测试;假定利用两次拼接测量即可完成全口径的测量,子孔径拼接示意图如图2所示,阴影部分代表第一子孔径1和第二子孔径2的重叠区域4。利用三坐标测量仪检测并采集到子第一孔径1区域的面形分布Z1(x1,y1,z1)。
四、测量靶标坐标;在第一子孔径1和第二子孔径2的重叠区域4贴上三个靶标3,靶标3为中心带十字线的圆形粘贴,结合图2,利用三坐标测量仪的ViScan光学探头能够精确测定到三个靶标3各自中心十字叉点的坐标Z11(x11,y11,z11),Z12(x12,y12,z12),Z13(x13,y13,z13)。
五、镜体位置调整;将大口径非球面镜体的位置进行调整,使三坐标测量仪行程能够完全测定第二子孔径2区域。
六、靶标坐标再次测定;在第2个测试位置,利用三坐标测量仪的ViScan光学探头再次测定到三个靶标各自中心十字叉点的坐标Z21(x21,y21,z21),Z22(x22,y22,z22),Z23(x23,y23,z23)。
七、其它子孔径测试;移去靶标3,利用三坐标测量仪测定得到第二子孔径2区域的面形分布Z2(x2,y2,z2)。
如有多个子孔径才能完成拼接测量,则重复步骤三至步骤七;
八、基准统一;通过坐标变换,可以求解得到两个子孔径之间的平移量和旋转量,从而使坐标系统一。
具体步骤算法如下:
由于不同子孔径测量之间的位置调整都在相同的基面上,因此位置调整仅会带来x方向和y方向的平移以及绕z方向的转动。假定两次测量两个子孔径在x方向、y方向上的相对平移量分别为dx、dy,绕z轴的转动角度量为θ,则根据坐标变换理论,相邻两个子孔径上三个靶标中心点位置坐标有如下关系:
( x 11 , y 11 , 1 ) = ( x 21 , y 21 , 1 ) . V ( x 12 , y 12 , 1 ) = ( x 22 , y 22 , 1 ) . V ( x 13 , y 13 , 1 ) = ( x 23 , y 23 , 1 ) . V - - - ( 1 )
其中V为两个坐标系之间的变换矩阵,推导可得V的表述如下式:
V = cos θ sin θ 0 - sin θ cos θ 0 d x d y 0 - - - ( 2 )
因为三个特征点两次测量的坐标值是已知的,通过公式(1)-(2),利用迭代算法即可得到第二子孔径2相对第一子孔径1在x方向、y方向上的平移量dx、dy和绕z轴的转动量θ,从而可以将两个子孔径统一到相同的基准上。
九、求解重叠区域;将两个子孔径统一基准后,分析测量区域,在某一区域两次测量均有采样数据,则该区域即为重叠区域4,判定重叠区域4是否大于子孔径面积的1/4,如大于1/4则进行下一步处理,否则重复步骤五至步骤九。
十、求解拼接因子;通过最小二乘拟合求解相邻两个子孔径间的调整系数,可以将两个子孔径统一基准,但不可能完全校准实际测量时的位置调整量,对于离轴非球面的检测,相邻两个子孔径间的调整误差将会引入五种初级像差,即为平移、倾斜、离焦、像散和彗差。
假使以第一子孔径1为基准,则第二子孔径2的面形分布Z2与基准第一子孔径1面形分布Z1的关系为:
Z = Z 2 + a 1 + a 2 x 1 + a 3 y 1 + a 4 ( x 1 2 + y 1 2 ) + a 5 x 1 y 1 + a 6 ( x 1 2 - y 1 2 ) + a 7 x 1 ( x 1 2 + y 1 2 ) + a 8 y 1 ( x 1 2 + y 1 2 ) - - - ( 3 )
其中a1、a2、a3和a4分别是第二子孔径2相对基准第一子孔径1的相对平移系数、沿x方向的倾斜系数、沿y方向的倾斜系数和相对离焦系数,a5和a6是相对象散系数,a7和a8是相对彗差系数。
利用最小二乘法,根据公式(4)可求解拼接系数:
A=B-1C         (4)
其中:
A = a 1 a 2 a 3 a 4 a 5 a 6 a 7 a 8 - - - ( 5 )
B = Σ x 1 2 Σ x 1 y 1 Σ x 1 ( x 1 2 + y 1 2 ) Σ x 1 2 y 1 Σ x 1 ( x 1 2 - y 1 2 ) Σ x 1 2 ( x 1 2 + y 1 2 ) Σ x 1 y 1 ( x 1 2 + y 1 2 ) Σ x 1 Σ x 1 y 1 Σ y 1 2 Σ y 1 ( x 1 2 + y 1 2 ) Σ x 1 y 1 2 Σ y 1 ( x 1 2 - y 1 2 ) Σ x 1 y 1 ( x 1 2 + y 1 2 ) Σ y 1 2 ( x 1 2 + y 1 2 ) Σy 1 Σ x 1 ( x 1 2 + y 1 2 ) Σ y 1 ( x 1 2 + y 1 2 ) Σ ( x 1 2 + y 1 2 ) 2 Σ x 1 y 1 ( x 1 2 + y 1 2 ) Σ ( x 1 4 - y 1 4 ) Σ x 1 ( x 1 2 + y 1 2 ) 2 Σ y 1 ( x 1 2 + y 1 2 ) 2 Σ ( x 1 2 + y 1 2 ) Σ x 1 2 y 1 Σ x 1 y 1 2 Σ x 1 y 1 ( x 1 2 + y 1 2 ) Σ x 1 2 y 1 2 Σ x 1 y 1 ( x 1 2 - y 1 2 ) Σ x 1 2 y 1 ( x 1 2 + y 1 2 ) Σ x 1 y 1 2 ( x 1 2 + y 1 2 ) Σ x 1 y 1 Σ x 1 ( x 1 2 - y 1 2 ) Σ y 1 ( x 1 2 - y 1 2 ) Σ ( x 1 4 - y 1 4 ) Σ x 1 y 1 ( x 1 2 - y 1 2 ) Σ ( x 1 2 - y 1 2 ) 2 Σ x 1 ( x 1 4 - y 1 4 ) Σ y 1 ( x 1 4 - y 1 4 ) Σ ( x 1 2 - y 1 2 ) Σ x 1 2 ( x 1 2 + y 1 2 ) Σ x 1 y 1 ( x 1 2 + y 1 2 ) Σ x 1 ( x 1 2 + y 1 2 ) 2 Σ x 1 2 y 1 ( x 1 2 + y 1 2 ) Σ x 1 y 1 2 ( x 1 2 + y 1 2 ) Σ x 1 2 ( x 1 2 + y 1 2 ) 2 Σ x 1 y 1 ( x 1 2 + y 1 2 ) 2 Σ x 1 ( x 1 2 + y 1 2 ) Σ x 1 y 1 ( x 1 2 + y 1 2 ) Σ y 1 2 ( x 1 2 + y 1 2 ) Σ y 1 ( x 1 2 + y 1 2 ) 2 Σ x 1 y 1 2 ( x 1 2 + y 1 2 ) Σ y 1 ( x 1 4 - y 1 4 ) Σ x 1 y 1 ( x 1 2 + x 1 2 ) 2 Σ y 1 2 ( x 1 2 + y 1 2 ) 2 Σ y 1 ( x 1 2 + y 1 2 ) Σ x 1 Σ y 1 Σ ( x 1 2 + y 1 2 ) Σ x 1 y 1 Σ ( x 1 2 - y 1 2 ) Σ x 1 ( x 1 2 + y 1 2 ) Σ y 1 ( x 1 2 + y 1 2 ) Σ n ij - - - ( 6 )
C = Σ x 1 Δz Σy 1 Δz Σ ( x 1 2 + y 1 2 ) Δz Σ x 1 y 1 Δz Σ ( x 1 2 - y 1 2 ) Δz Σ x 1 ( x 1 2 + y 1 2 ) Δz Σ y 1 ( x 1 2 + y 1 2 ) Δz ΣΔz - - - ( 7 )
Δz=Z1-Z2 (8)
对于同轴非球面的检测,相邻两个子孔径间的调整误差仅会引入平移、倾斜和彗差三种初级像差,即仅需要求解拼接系数a1、a2、a3、a7和a8
十一、全口径拼接,根据公式(4)-(8)即可求解出第二子孔径2相对于第一子孔径1由于调整引入的误差,从而把两个子孔径拼接起来。即使有多个子孔径拼接,利用两两拼接原理,就可完成对大口径非球面面形的准确测量。

Claims (2)

1.采用三坐标测量仪拼接测量大口径非球面面形的方法,其特征是,该方法由以下步骤实现:
步骤一、对被检测大口径非球面的尺寸进行划分,获得多个子孔径;所述多个子孔径在三坐标测量仪的检测范围内,相邻子孔径有四分之一的区域为重叠区域;
步骤二、采用三坐标测量仪检测并采集子孔径区域的面形数据,并在该子孔径区域与相邻的另一个子孔径区域的重叠区域贴上三个靶标,采用三坐标测量仪测定三个靶标的数据;
步骤三、调整被检测大口径非球面的位置,使三坐标测量仪能够检测到另一个子孔径区域的面形数据;采用三坐标测量仪测定三个靶标的数据;移去三个靶标,三坐标测量仪获得另一个子孔径区域的面形数据;
步骤四、判断是否是最后一个子孔径,如果是,则执行步骤五;如果否,返回执行步骤二;
步骤五、采用迭代算法将步骤二和步骤三获得的子孔径区域的面形数据进行坐标变换,获得在相同基准面上的两个子孔径数据;
步骤六、判断步骤五获得的两个子孔径的重叠区域是否大于四分之一,如果否,则返回执行步骤三;如果是,则采用最小二乘法获得两个子孔径的拼接系数,最终实现对大口径非球面面形的检测。
2.根据权利要求1所述的采用三坐标测量仪拼接测量大口径非球面面形的方法,其特征在于,所述三个靶标为中心带十字线的圆形粘贴。
CN201210570644.5A 2012-12-25 2012-12-25 采用三坐标测量仪拼接测量大口径非球面面形的方法 Active CN102980532B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210570644.5A CN102980532B (zh) 2012-12-25 2012-12-25 采用三坐标测量仪拼接测量大口径非球面面形的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210570644.5A CN102980532B (zh) 2012-12-25 2012-12-25 采用三坐标测量仪拼接测量大口径非球面面形的方法

Publications (2)

Publication Number Publication Date
CN102980532A true CN102980532A (zh) 2013-03-20
CN102980532B CN102980532B (zh) 2015-06-10

Family

ID=47854774

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210570644.5A Active CN102980532B (zh) 2012-12-25 2012-12-25 采用三坐标测量仪拼接测量大口径非球面面形的方法

Country Status (1)

Country Link
CN (1) CN102980532B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103257033A (zh) * 2013-05-15 2013-08-21 中国科学院光电技术研究所 基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法
CN103776391A (zh) * 2014-02-19 2014-05-07 厦门大学 用于磨削阶段大口径非球面光学元件轮廓测量方法
CN103822605A (zh) * 2014-03-18 2014-05-28 厦门大学 大口径光学元件轮廓的一次拼接测量装置
CN106248352A (zh) * 2016-09-13 2016-12-21 湖北航天技术研究院总体设计所 一种大口径平面镜拼接检测对准方法
CN109341535A (zh) * 2018-12-13 2019-02-15 中国航空工业集团公司北京长城计量测试技术研究所 一种高效高精度大型卫星天线装配测量系统
CN111308593A (zh) * 2020-03-30 2020-06-19 南京理工大学 用于手机镜头的超薄型反射式自由曲面设计方法
CN114252244A (zh) * 2021-12-17 2022-03-29 华中科技大学 一种二次曲面反射镜拼接干涉检测子孔径规划方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3309743B2 (ja) * 1996-11-27 2002-07-29 富士ゼロックス株式会社 形状測定方法及び装置
JP2005043249A (ja) * 2003-07-23 2005-02-17 Brother Ind Ltd 3次元形状検出装置、撮像装置、及び、3次元形状検出方法
JP2010133751A (ja) * 2008-12-02 2010-06-17 Topcon Corp 形状測定装置およびプログラム
CN101946559A (zh) * 2008-02-19 2011-01-12 松下电器产业株式会社 感应加热烹调器

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3309743B2 (ja) * 1996-11-27 2002-07-29 富士ゼロックス株式会社 形状測定方法及び装置
JP2005043249A (ja) * 2003-07-23 2005-02-17 Brother Ind Ltd 3次元形状検出装置、撮像装置、及び、3次元形状検出方法
CN101946559A (zh) * 2008-02-19 2011-01-12 松下电器产业株式会社 感应加热烹调器
JP2010133751A (ja) * 2008-12-02 2010-06-17 Topcon Corp 形状測定装置およびプログラム

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王孝坤等: "子孔径拼接干涉检测非球面技术的研究", 《光学测试新理论、新技术》 *
王孝坤等: "子孔径拼接技术在大口径高陡度非球面检测中的应用", 《强激光与粒子束》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103257033A (zh) * 2013-05-15 2013-08-21 中国科学院光电技术研究所 基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法
CN103776391A (zh) * 2014-02-19 2014-05-07 厦门大学 用于磨削阶段大口径非球面光学元件轮廓测量方法
CN103822605A (zh) * 2014-03-18 2014-05-28 厦门大学 大口径光学元件轮廓的一次拼接测量装置
CN106248352A (zh) * 2016-09-13 2016-12-21 湖北航天技术研究院总体设计所 一种大口径平面镜拼接检测对准方法
CN109341535A (zh) * 2018-12-13 2019-02-15 中国航空工业集团公司北京长城计量测试技术研究所 一种高效高精度大型卫星天线装配测量系统
CN109341535B (zh) * 2018-12-13 2020-05-29 中国航空工业集团公司北京长城计量测试技术研究所 一种高效高精度大型卫星天线装配测量系统
CN111308593A (zh) * 2020-03-30 2020-06-19 南京理工大学 用于手机镜头的超薄型反射式自由曲面设计方法
CN114252244A (zh) * 2021-12-17 2022-03-29 华中科技大学 一种二次曲面反射镜拼接干涉检测子孔径规划方法及系统
CN114252244B (zh) * 2021-12-17 2022-09-20 华中科技大学 一种二次曲面反射镜拼接干涉检测子孔径规划方法及系统

Also Published As

Publication number Publication date
CN102980532B (zh) 2015-06-10

Similar Documents

Publication Publication Date Title
CN105318891B (zh) 一种星敏感器基准立方镜安装误差的标定装置
CN102980532B (zh) 采用三坐标测量仪拼接测量大口径非球面面形的方法
CN101709955B (zh) 子孔径拼接干涉检测光学非球面面形的装置
CN102430959A (zh) 数控机床转台运动误差的快速检测方法
CN102288132B (zh) 采用激光跟踪仪测量非球面顶点曲率半径偏差的方法
CN103175486B (zh) 一种圆柱度误差的拼接干涉测量装置及方法
Chao et al. Calibration of laser beam direction for optical coordinate measuring system
CN102128599B (zh) 接触式非球面面形检验装置
CN101216299A (zh) 六自由度关节式柔性坐标测量机高精度标定方法
CN105571527A (zh) 一种转台摆角精密测量方法
CN109141225A (zh) 基于圆光栅的轴系五、六自由度误差测量方法及测量系统
CN104458124A (zh) 一种质心测量方法
CN103542791A (zh) 法兰多参数集成测量方法及其测量装置
CN103591888A (zh) 大口径离轴非球面光学元件几何参数的测算方法
CN103791868A (zh) 一种基于虚拟球的空间标定体及其标定方法
CN102944204A (zh) 利用二维测量机检测截面轮廓度的方法
CN102873586B (zh) 数控加工工件曲率半径快速在线测量装置
CN110806571A (zh) 一种多结构光传感器空间姿态标定件及其标定方法
CN107990856A (zh) 一种超量程工件的空间位置误差检测方法
CN106989670A (zh) 一种机器人协同的非接触式高精度大型工件跟踪测量方法
CN102506761B (zh) 一种利用激光跟踪仪测量非球面顶点曲率半径的方法
Liu et al. A coaxiality measurement model based on multi-systematic errors separation for turbine shaft
CN103294920A (zh) 一种基于最小区域的零件圆柱度误差评定方法
CN102927930B (zh) 采用平行光管拼接检测超大口径反射镜面形误差的方法
CN105627945A (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
C14 Grant of patent or utility model
GR01 Patent grant