CN107560564A - 一种自由曲面检测方法及系统 - Google Patents
一种自由曲面检测方法及系统 Download PDFInfo
- Publication number
- CN107560564A CN107560564A CN201710635739.3A CN201710635739A CN107560564A CN 107560564 A CN107560564 A CN 107560564A CN 201710635739 A CN201710635739 A CN 201710635739A CN 107560564 A CN107560564 A CN 107560564A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- measured
- locations
- parameter
- 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
Links
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
本发明公开了一种自由曲面检测方法，包括：采用三坐标测量设备对包括待测反射球面、投影屏和CCD相机在内的逆向哈特曼检验光路的结构位置参数S进行测量标定；根据测量标定的结构位置参数S，确定待测反射球面的面形误差数据W0；根据测量标定的结构位置参数S和面形误差数据W0，采用泽尼克拟合确定待测球面的面形偏差优化目标；根据确定的面形偏差优化目标，确定初始测量标定的结构位置参数S的各项偏差，并根据所述各项偏差确定待测反射面的实际面形误差ΔW。采用本发明，减小了自由曲面的面形测量误差。
Description
技术领域
本发明涉及一种曲面检测技术领域，尤其涉及一种自由曲面检测方法及系统。
背景技术
由于光学系统的设计要求越来越精密，对其加工精度的要求也相应越来越高，除去一些在日常中用到的凹面反射镜，在试验机生产中对于自由曲面反射镜的需求也越来越多样.要实现对与光学反射面的纳米级精度的检测，现有技术中最为成熟时光学关涉仪检测法，而干涉仪往往存在测量动态范围较小，通用性较差、检测成本高等缺点。一种软件可配置光学测量系统(详见Peng Su等，Software configurable optical test system：acomputerized reverse Hartmann test，Applied Optics，2010，49(23)：4404-4412)采用基于逆向哈特曼检验光路，实现了对反射面面形误差的快速、高测量动态范围和高精度的测量。但该方法需要预知面形，其最为重要的是对光学系统结构参数的预标定工作，进而实现对系统结构误差的校正。预标定的过程是通过繁琐且高要求的三坐标测量装置来实现，对于数值孔径较小的凹反射面的面形测量，微米级的标定精度不会对测量结果造成明显的误差引入。但是对于自由曲面的检测时，尤其当被检曲面类似于凸球面时，系统结构位置参数标定误差会被放大很多，从而产生较为明显的面形测量误差。
发明内容
本发明的目的是提供一种自由曲面检测方法，提高了自由曲面的检测精度。
为解决现有技术存在的问题，本发明提供一种自由曲面检测方法，该方法包括：
s1，采用三坐标测量设备对包括待测反射球面、投影屏和CCD相机在内的逆向哈特曼检验光路的结构位置参数S进行测量标定；
s2，根据测量标定的结构位置参数S，确定待测反射球面的面形误差数据W0；
s3，根据测量标定的结构位置参数S和面形误差数据W0，采用泽尼克拟合，并根据泽尼克拟合系数确定待测球面的面形偏差优化目标；
s4，根据确定的面形偏差优化目标，确定初始测量标定的结构位置参数S的各项偏差，并根据所述各项偏差确定待测反射面的实际面形误差ΔW。
本技术方案首先测量标定结构位置参数实现零位检测，然后通过泽尼克拟合，并根据泽尼克拟合系数确定偏差优化目标，从而实现对待测反射面的实际面形误差的确定，提高了自由曲面面形测量的精度。
另外，所述步骤s1中的结构位置参数为S＝{(xi，yi，zi)，(αi，βi，γi)}i＝1，2，3，其中(xi，yi，zi)表示第i个元件的三维空间坐标位置，(αi，βi，γi)表示第i 个元件关于各个坐标轴的倾角。为后续的零位检测准备了条件。
另外，所述步骤s2包括：
移相并采集条纹：控制投影屏前后移动显示x水平方向和y竖直方向的四步90°移相正弦直条纹，确定条纹间距坐标Dx和Dy相位分布Φx和Φy；
根据所述确定的条纹间距坐标和相位分布，确定待测球面投影在投影屏上的横、纵坐标：
Xm＝Φx·Dx/2π；
Ym＝Φy·Dy/2π；
确定与待测球面面形误差对应的局部斜率(wx，wy)，该局部斜率为：
wx＝(Xm-Xideal)/2dms
wy＝(Ym-Yideal)/2dms
其中，投影坐标数据Xm和Ym的理论值分别为Xideal和Yideal；
对局部斜率(wx，wy)积分得到的面形误差数据W0。
本技术方案实现了对自由曲面面形误差的初始测得值W0的获取。
另外，所述步骤s3包括：
s31，以逆向哈特曼检测光路的结构位置参数S作为变量，不断调整各结构误差参数获取相应的检测结果，对所述检测结果进行泽尼克面形拟合：
其中，W0为初始测得值，Ci，geo和Ci，surf分别为预标定中的结构误差 Wgeo与被测物表面误差Wsurf的第i项泽尼克系数，j表示结构位置参数S中的第j项结构误差因素，Wj为仅改变j因素以后测得的面形数据， N为泽尼克系数的项数，Ci为初始测得值的泽尼克拟合得到的第i项泽尼克系数，Ci，j为改变j因素以后泽尼克拟合得到的第i项泽尼克系数，Zi为第i项泽尼克项；
s32，根据待测球面的面形误差数据W0和仅改变j因素以后测得的面形数据Wj确定泽尼克拟合系数比值。
另外，所述步骤s32包括：
s321，通过差分计算确定：
其中，WΔ，j是j因素的变动造成的测量结果的变动量，ΔCi，j为改变j 因素以后测得的面形数据与初始值各自拟合系数的在第i项上的系数差值；
s322，获取真实面形的泽尼克系数，高阶项(Ci，j-Ci，surf)(i＞3)与离焦项数据C3，j的比值：ri，j＝(Ci，j-Ci，surf)/C3，j＝ΔCi，j/ΔC3，j。
另外，所述步骤s4包括：
s41，确定泽尼克系数的前m项组成的波面Wm为主要结构误差，将所述波面作为优化目标，确定偏差函数：
{fm}＝{Wm}
根据偏差函数{fm}，使收敛时对应的逆向哈特曼检验光路的结构位置参数S在一定误差范围内。
s42采用迭代优化算法，改变结构位置参数S，通过结构误差参数与泽尼克系数的函数方程ψi，j(εj)，求得相应泽尼克系数；εj为位置参数较初始结构参数值的偏差；
s43，考虑所有结构位置参数S，基于真值比值条件ri，j(i＞3)获得初始值中结构位置参数S的各项偏差εj对于测量结果中高阶项的影响：
s44，以面形数据Wsurf为优化目标，确定广义加权偏差函数
{fk}＝{ρkWsurf}，
其中，ρk为加权因子，{fk}为偏差函数；
s45，根据偏差函数{fk}，确定收敛时对应的逆向哈特曼检验光路的结构位置参数Smin；
s46，循环获取机构位置参数Smin，，确定待测反射面的实际面形误差为ΔW。
相应的本发明还提供一种自由曲面检测系统，该系统包括：
三坐标测量设备、控制计算机，以及成像设备，成像用于对待测反射球面进行成像，所述三坐标检测设备用于对；所述成像设备包括：投影仪、CCD相机组，所述CCD组前端装有滤光小孔，所述CCD 组与待测反射面组成逆向哈特曼检验光路；所述待测反射球面面向投影屏和CCD相机，所述待测反射球面与投影屏相对位置使待测反射球面在CCD相机中完整成像。
本发明自由曲面检测方法，基于泽尼克拟合系数比值找到个系统结构位置参数标定误差对检测结构面形的泽尼克拟合系数的影响，提出了涉及泽尼克拟合系数比值的结构测量误差校正方法，主要目的为最小化系统结构误差，从而实现纳米量级的反射球面面形误差的检测。
附图说明
图1是本发明一种自由曲面检测系统中逆向哈特曼检验光路示意图；
图2是本发明一种自由曲面检测方法的第一种实施例的流程示意图；
图3是本发明一种自由曲面检测系统中的立项表面与存在加工误差表面的对比图；
图4是本发明一种自由曲面检测方法中泽尼克拟合得出的表面加工误差真值示意图；
图5是本发明一种自由曲面检测方法中的优化结果示意图。
图中：
1.投影屏，2.CCD相机，3.待测反射球面。
具体实施方式
下面结合附图对本发明进行详细说明。
本发明实施例的基于泽尼克拟合系数比值的自由曲面通用化检测方法，采用逆向哈特曼检验系统，如图1所示，逆向哈特曼检验系统包括投影屏1、CCD相机2、待测反射球面3，投影屏1和CCD相机2并排放置，且CCD相机前端装有滤光小孔4，待测反射球面3 面向投影屏1和CCD相机2，并且待测反射球面3与投影屏1平面近似平行放置，投影屏1和CCD相机2的数据线分别和计算机5相连。
下面说明本发明的另一方面。
参考图2，该图是本发明一种自由曲面检测方法的第一种实施例的流程示意图，该流程包括：
步骤s1，采用三坐标测量设备对包括待测反射球面、投影屏和CCD 相机在内的逆向哈特曼检验光路的结构位置参数S进行测量标定；
步骤s2，根据测量标定的结构位置参数S，确定待测反射球面的面形误差数据W0；
步骤s3，根据测量标定的结构位置参数S和面形误差数据W0，采用泽尼克拟合，并根据泽尼克拟合系数确定待测球面的面形偏差优化目标；
步骤s4，根据确定的面形偏差优化目标，确定初始测量标定的结构位置参数S的各项偏差，并根据所述各项偏差确定待测反射面的实际面形误差ΔW。
下面对具体实现时基于泽尼克拟合系数比值的自由曲面检测方法进行详细说明，其中，具体实现时，以一个口径为50.8mm，曲率半径为250mm的凸球面反射镜为待测自由曲面为例，进行面形误差测量，其中，图3(a)为该凸面镜的理想表面，图3(b)为存在表面误差的同规格凸面镜，图4为图3(b)上的表面误差，包括如下步骤：
(1-1)实验设备的初始化安装与预标定
安置实验装置：由投影屏、待测反射面和前端装有滤光小孔的 CCD相机组成逆向哈特曼检验光路，其中屏与CCD相机并排放置，待测反射球面面向投影屏和CCD相机，并且待测反射球面与投影屏平面近似平行放置，通过投影屏相对地前后平移与旋转微调待测反射球面，使得待测反射球面于CCD相机中完整成像。
预标定工作：利用测量精度可达微米量级的三坐标测量设备对所述的包括待测反射球面、投影屏和CCD相机在内的逆向哈特曼检验光路的结构位置参数S进行测量标定，S＝{(xi，yi，zi)，(αi，βi，γi)}i＝1，2，3，其中(xi，yi，zi)表示第i个元件的三维空间坐标位置，(αi，βi，γi)表示第i个元件关于各个坐标轴的倾角；
(1-2)通过移相求解相位分布
移相并采集条纹：利用计算机控制投影屏前后显示x水平方向和 y竖直方向的四步90°移相正弦直条纹，其对应的条纹间距分别为Dx和Dy，再由CCD相机实时地采集移相正弦条纹经待测球面反射后的像。
求解相位：分别利用四步移相算法求解CCD相机采集到的移相正弦条纹对应的相位分布Φx和Φy，并由此得到其投影在投影屏上的横、纵坐标Xm和Ym为：
Xm＝Φx·Dx/2π，
Ym＝Φy·Dy/2π；
实现“零位”检测：利用预标定中得到的逆向哈特曼检验光路的结构位置参数S，采用光线追迹法得到求解相位步骤中对应投影坐标数据Xm和Ym的理论值分别为Xideal和Yideal；
计算面形误差W0对应的局部斜率(wx，wy)为：
wx＝(Xm-Xideal)/2dms，
wy＝(Ym-Yideal)/2dms，
其中，dms为步骤(1-1)标定测得的待测反射球面到投影屏的距离；对局部斜率(wx，wy)积分得到的面形误差数据W0；
(1-3)泽尼克拟合与结构位置参数优化
步骤一.泽尼克拟合：以所述的逆向哈特曼检验光路的结构位置参数S作为变量，不断调整各结构误差参数，得到相关检测结果，并对其进行泽尼克面形拟合：
其中，W0为初始测得值，Ci，geo和Ci，surf分别为预标定中的结构误差 Wgeo与被测物表面误差Wsurf的第i项泽尼克系数，j表示结构位置参数 S中的第j项结构误差因素，Wj为仅改变j因素以后测得的面形数据， N为泽尼克系数的项数，Ci为初始测得值的泽尼克拟合得到的第i项泽尼克系数，Ci，j为改变j因素以后泽尼克拟合得到的第i项泽尼克系数，Zi为第i项泽尼克项。从而建立各结构位置参数较初始结构参数值的偏差εj与其相应的泽尼克系数的函数方程ψi，j(εj)，即可以通过结构误差数值得到相应的面形对应的泽尼克系数。
步骤二.求解泽尼克拟合系数比值ri，j：通过差分计算
ΔCi，j为改变j因素以后测得的面形数据与初始值各自拟合系数的在第 i项上的系数差值。从而获取相关真实面形的泽尼克系数中，高阶项 (Ci，j-Ci，surf)(i＞3)与离焦项数据C3，j的比值：
ri，j＝(Ci，j-Ci，surf)/C3，j＝ΔCi，j/ΔC3，j；
步骤三.结构参数优化：
1.以泽尼克系数的前m项组成的波面Wm为主要结构误差，将其面形偏差作为优化目标，使整体结构参数误差缩小到一定范围内：
优先校正灵敏度高的角度误差因素，以所述的逆向哈特曼检验光路的结构位置倾角参数(αi，βi，γi)作为优化变量，本步骤中的面形数据 Wm作为优化目标，定义广义加权偏差函数：
以面形数据Wm为优化目标，确定广义加权偏差函数
{fm}＝{Wm}，
根据偏差函数{fm}，使收敛时对应的逆向哈特曼检验光路的结构位置参数S在一定误差范围内。利用迭代优化算法，改变结构位置参数S，通过结构误差参数与泽尼克系数的函数方程ψi，j(εj)，求得相应泽尼克系数，使得偏差函数{fm}迭代减小收敛；其次优化考虑距离误差因素 (xi，yi，zi)，使得相同的偏差函数{fm}迭代减小且收敛。从而将结构位置参数S的参数误差缩小到一定的范围。
2.考虑所有结构位置参数S，基于真值比值条件ri，j(i＞3)获得初始值中结构位置参数S的各项偏差εj对于测量结果中高阶项的影响
定义广义加权偏差函数
{fk}＝{ρkWsurf}，
其中，ρk为加权因子。在一定误差范围内使偏差函数{fk}迭代减小收敛。最终得到偏差函数{fk}收敛时对应的逆向哈特曼检验光路的结构位置参数Smin，重复步骤(1-2)中面形误差求解，进而得到待测反射面的实际面形误差为ΔW。
本发明在逆向哈特曼检验光路基础上，利用计算机辅助测量技术，对实际测量中的系统结构位置参数标定误差进行高精度校正，进而实现可通用化的反射球面高精度检测，其测量精度可达到纳米量级。
再以一个具体的测量实例作进一步说明。
采用本发明的方法对一个口径为50.8mm，曲率半径为250mm的凸球面反射镜进行面形误差测量，包括如下步骤：
(1-1)实验设备的初始化安装与标定
如图1所示，将待测反射凸球面置于逆向哈特曼检验光路中，以 LCD显示器作为投影屏，其像素分辨率为1920×1080，所述的逆向哈特曼检验光路中投影屏和CCD相机并排放置，且CCD相机前端装有滤光小孔，待测反射凸球面面向投影屏和CCD相机，并且待测反射凸球面与投影屏平面近似平行放置，通过相对投影屏前后平移待测反射凸球面，使得CCD相机对待测反射凸球面完整成像；
利用测量精度可达1.9μm的海克斯三坐标测量机对所述的包括待测反射凸球面、投影屏和CCD相机在内的逆向哈特曼检验光路的结构位置参数S进行测量标定；
(1-2)通过移相求解相位分布
利用计算机控制投影屏分别显示x水平方向和y竖直方向的四步 90°移相正弦直条纹，对应的条纹间距都为17.1mm，由CCD相机实时采集移相正弦条纹经待测球面反射后的像。
分别利用四步移相算法求解CCD相机中采集到的移相正弦条纹对应的相位分布Φx和Φy，分别并由此得到其投影在投影屏上的横、纵坐标Xm和Ym为：
Xm＝(Φx/2π)×17.1，Ym＝(Φy/2π)×17.1；
利用预标定中得到的逆向哈特曼检验光路的结构位置参数S，采用光线追迹法得到对应投影坐标数据Xm和Ym的理论值分别为Xideal和 Yideal，得到实际测得的光线投影坐标和光线追迹法得到投影坐标偏差分布。由此可知由于凸球面对于系统结构位置参数标定误差的放大作用，会对最后的面形误差测量引入较大的误差，难以满足高精度测量要求。
计算面形误差W对应的局部斜率(wx，wy)为：
wx＝(Xm-Xideal)/2dms，
wy＝(Ym-Yideal)/2dms，
其中，dms为步骤(1-1)标定测得的待测反射球面到投影屏的距离，为151.981mm。对局部斜率(wx，wy)积分得到的面形误差数据W0；
(1-3)泽尼克拟合与结构位置参数优化
步骤一.泽尼克拟合：以所述的逆向哈特曼检验光路的结构位置参数S作为变量，不断调整各结构误差参数，得到相关检测结果，并对其进行泽尼克面形拟合：
其中，W0为初始测得值，Ci，geo和Ci，surf分别为预标定中的结构误差Wgeo与被测物表面误差Wsurf的第i项泽尼克系数，j表示结构位置参数S 中的第j项结构误差因素，Wi为仅改变j因素以后测得的面形数据，N为泽尼克系数的项数，Ci为初始测得值的泽尼克拟合得到的第i项泽尼克系数，Ci，j为改变j因素以后泽尼克拟合得到的第i项泽尼克系数，Zi为第i项泽尼克项。从而建立各结构位置参数较初始结构参数值的偏差εj与其相应的泽尼克系数的函数方程ψi，j(εj)，即可以通过结构误差数值得到相应的面形对应的泽尼克系数。
步骤二.求解泽尼克拟合系数比值ri，j：通过差分计算
ΔCi，j为改变j因素以后测得的面形数据与初始值各自拟合系数的在第 i项上的系数差值。从而获取相关真实面形的泽尼克系数中，高阶项 (Ci，j-Ci，surf)(i＞3)与离焦项数据C3，j的比值：
ri，j＝(Ci，j-Ci，surf)/C3，j＝ΔCi，j/ΔC3，j。
步骤三.结构参数优化：
1.以泽尼克系数前m项组成的波面Wm作为主要的结构误差，将其面形偏差作为优化目标，使整体结构参数误差缩小到一定范围内：
优先校正灵敏度高的角度误差因素，以所述的逆向哈特曼检验光路的结构位置倾角参数(αi，βi，γi)作为优化变量，本步骤中泽尼克系数的前m项组成的波面Wm为主要结构误差，以面形数据Wm为优化目标，确定广义加权偏差函数
{fm}＝{Wm}，
利用迭代优化算法，改变结构位置参数S，通过结构误差参数与泽尼克系数的函数方程ψi，j(εj)，求得相应泽尼克系数，使得偏差函数{fm}迭代减小收敛；其次优化考虑距离误差因素，使得相同的偏差函数{fm} 迭代减小且收敛。从而将结构位置参数S的参数误差缩小到一定的范围。
2.考虑所有结构位置参数S，基于真值比值条件ri，j(i＞3)获得初始值中结构位置参数S的各项偏差εj对于测量结果中高阶项的影响
定义广义加权偏差函数
{fk}＝{ρkWsurf}，
在限定误差范围内使偏差函数{fk}迭代减小收敛。
得到偏差函数{fk}收敛时对应的逆向哈特曼检验光路的结构位置参数Smin，重复步骤(1-2)中面形误差求解，进而得到待测反射面的实际面形误差为ΔW。
得到偏差函数{fk}收敛时对应的逆向哈特曼检验光路的结构位置参数Smin，重复步骤(1-2)中面形误差求解，进而得到待测反射面的实际面形误差为ΔW。
得到偏差函数{fk}收敛时对应的逆向哈特曼检验光路的结构位置参数Smin，重复步骤中的面形误差求解，进而得到待测反射面的实际面形误差为ΔW，如图5优化结果图所示，其均方根RMS值为 3.4656μm。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。
Claims (7)
1.一种自由曲面检测方法，其特征在于，包括：
s1，采用三坐标测量设备对包括待测反射球面、投影屏和CCD相机在内的逆向哈特曼检验光路的结构位置参数S进行测量标定；
s2，根据测量标定的结构位置参数S，确定待测反射球面的面形误差数据W0；
s3，根据测量标定的结构位置参数S和面形误差数据W0，采用泽尼克拟合，并根据泽尼克拟合系数确定待测球面的面形偏差优化目标；
s4，根据确定的面形偏差优化目标，确定初始测量标定的结构位置参数S的各项偏差，并根据所述各项偏差确定待测反射面的实际面形误差ΔW。
2.根据权利要求1所述的自由曲面检测方法，其特征在于，所述步骤s1中的结构位置参数为S＝{(xi，yi，zi)，(αi，βi，γi)}i＝1，2，3，其中(xi，yi，zi)表示第i个元件的三维空间坐标位置，(αi，βi，γi)表示第i个元件关于各个坐标轴的倾角。
3.根据权利要求1所述的自由曲面检测方法，其特征在于，所述步骤s2包括：
移相并采集条纹：控制投影屏前后移动显示x水平方向和y竖直方向的四步90°移相正弦直条纹，确定条纹间距坐标Dx和Dy，和相位分布Φx和Φy；
根据所述确定的条纹间距坐标和相位分布，确定待测球面投影在投影屏上的横、纵坐标：
Xm＝Φx·Dx/2π；
Ym＝Φy·Dy/2π；
确定与待测球面面形误差对应的局部斜率(wx，wy)，该局部斜率为：
wx＝(Xm-Xideal)/2dms
wy＝(Ym-Yideal)/2dms
其中，投影坐标数据Xm和Ym的理论值分别为Xideal和Yideal；
对局部斜率(wx，wy)积分得到的面形误差数据W0；。
4.根据权利要求1所述的自由曲面检测方法，其特征在于，所述步骤s3包括：
s31，以逆向哈特曼检测光路的结构位置参数S作为变量，不断调整各结构误差参数获取相应的检测结果，对所述检测结果进行泽尼克面形拟合：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>W</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>W</mi>
<mrow>
<mi>g</mi>
<mi>e</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>W</mi>
<mrow>
<mi>s</mi>
<mi>u</mi>
<mi>r</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>g</mi>
<mi>e</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>s</mi>
<mi>u</mi>
<mi>r</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>Z</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<msub>
<mi>Z</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>W</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>Z</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，W0为初始测得值，Ci，geo和Ci，surf分别为预标定中的结构误差Wgeo与被测物表面误差Wsurf的第i项泽尼克系数，j表示结构位置参数S中的第j项结构误差因素，Wj为仅改变j因素以后测得的面形数据，N为泽尼克系数的项数，Ci为初始测得值的泽尼克拟合得到的第i项泽尼克系数，Ci，j为改变j因素以后泽尼克拟合得到的第i项泽尼克系数，Zi为第i项泽尼克项；
s32，根据待测球面的面形误差数据W0和仅改变j因素以后测得的面形数据Wj确定泽尼克拟合系数比值。
5.根据权利要求4所述的自由曲面检测方法，其特征在于，所述步骤s32包括：
s321，通过差分计算确定：
<mrow>
<msub>
<mi>W</mi>
<mrow>
<mi>&Delta;</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>W</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<msub>
<mi>W</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>&Delta;C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>Z</mi>
<mi>i</mi>
</msub>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
其中，WΔ，j是j因素的变动造成的测量结果的变动量，ΔCi，j为改变j因素以后测得的面形数据与初始值各自拟合系数的在第i项上的系数差值；
s322，获取相关真实面形的泽尼克系数，高阶项(Ci，j-Ci，surf)(i＞3)与离焦项数据C3，j的比值：
ri，j＝(Ci，j-Ci，surf)/C3，j＝ΔCi，j/ΔC3，j。
6.根据权利要求1所述的自由曲面检测方法，其特征在于，所述步骤s4包括：
s41，确定泽尼克系数的前m项组成的波面Wm为主要结构误差，将所述波面作为优化目标，确定偏差函数：
{fm}＝{Wm}
根据偏差函数{fm}，使收敛时对应的逆向哈特曼检验光路的结构位置参数S在一定误差范围内。
s42采用迭代优化算法，改变结构位置参数S，通过结构误差参数与泽尼克系数的函数方程ψi，j(εj)，求得相应泽尼克系数；εj为位置参数较初始结构参数值的偏差；
s43，考虑所有结构位置参数S，基于真值比值条件ri，j(i＞3)获得初始值中结构位置参数S的各项偏差εj对于测量结果中高阶项的影响：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>s</mi>
<mi>u</mi>
<mi>r</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>W</mi>
<mrow>
<mi>s</mi>
<mi>u</mi>
<mi>r</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>s</mi>
<mi>u</mi>
<mi>r</mi>
<mi>f</mi>
</mrow>
</msub>
<msub>
<mi>Z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
s44，以面形数据Wsurf为优化目标，确定广义加权偏差函数：
{fk}＝{ρkWsurf}，
其中，ρk为加权因子，{fk}为偏差函数；
s45，根据偏差函数{fk}，确定收敛时对应的逆向哈特曼检验光路的结构位置参数Smin；
s46，循环获取机构位置参数Smin，确定待测反射面的实际面形误差为ΔW。
7.一种自由曲面检测系统，其特征在于，该系统包括，三坐标测量设备、控制计算机，以及成像设备，成像用于对待测反射球面进行成像，所述三坐标检测设备用于对；所述成像设备包括：投影仪、CCD相机组，所述CCD组前端装有滤光小孔，所述CCD组与待测反射面组成逆向哈特曼检验光路；所述待测反射球面面向投影屏和CCD相机，所述待测反射球面与投影屏相对位置使待测反射球面在CCD相机中完整成像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710635739.3A CN107560564B (zh) | 2017-07-28 | 2017-07-28 | 一种自由曲面检测方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710635739.3A CN107560564B (zh) | 2017-07-28 | 2017-07-28 | 一种自由曲面检测方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107560564A true CN107560564A (zh) | 2018-01-09 |
| CN107560564B CN107560564B (zh) | 2019-09-10 |
Family
ID=60974263
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710635739.3A Active CN107560564B (zh) | 2017-07-28 | 2017-07-28 | 一种自由曲面检测方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107560564B (zh) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108436602A (zh) * | 2018-03-30 | 2018-08-24 | 中国石油大学(华东) | 一种利用未修整近圆弧截面砂轮磨削自由曲面的方法 |
| CN108507495A (zh) * | 2018-03-19 | 2018-09-07 | 中国计量大学 | 一种基于逆向哈特曼检测的自由曲面检测方法 |
| CN109307480A (zh) * | 2018-09-30 | 2019-02-05 | 中国计量大学 | 一种透射元件多表面面形检测方法 |
| CN109739019A (zh) * | 2019-03-14 | 2019-05-10 | 南开大学 | 基于矢量像差理论的共轴自由曲面光学系统优化设计方法 |
| CN110245439A (zh) * | 2019-06-20 | 2019-09-17 | 成都飞机工业（集团）有限责任公司 | 一种扩口导管的安装误差建模和补偿方法 |
| CN110645916A (zh) * | 2019-11-14 | 2020-01-03 | 北京航空航天大学 | 基于基准标定板校正位姿的自由曲面测量方法与装置 |
| CN110702009A (zh) * | 2019-11-22 | 2020-01-17 | 中国科学院长春光学精密机械与物理研究所 | 一种基于逆向Hartmann计算机辅助法的三维测量系统 |
| CN110763159A (zh) * | 2019-11-01 | 2020-02-07 | 中国计量大学 | 一种光学偏折显微表面测量装置及方法 |
| CN110793440A (zh) * | 2019-11-01 | 2020-02-14 | 中国计量大学 | 一种光学偏折瞬态测量方法 |
| CN114353699A (zh) * | 2022-01-07 | 2022-04-15 | 中国科学院长春光学精密机械与物理研究所 | 大陡度凸面光学自由曲面高频段像差检测系统及检测方法 |
| CN114593693A (zh) * | 2022-02-16 | 2022-06-07 | 苏州英示测量科技有限公司 | 光学偏折检测方法、电子设备及光学偏折检测系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101650157A (zh) * | 2009-09-18 | 2010-02-17 | 中国科学院长春光学精密机械与物理研究所 | 双曲面凸面反射镜面形误差的检测方法及其装置 |
| CN101943559A (zh) * | 2010-08-12 | 2011-01-12 | 中国科学院光电技术研究所 | 利用三坐标测量机检测大口径非球面光学元件的方法 |
| CN102297663A (zh) * | 2011-04-26 | 2011-12-28 | 北京理工大学 | 基于彩色编码逆向哈特曼式面形测量方法 |
| CN103267495A (zh) * | 2013-04-24 | 2013-08-28 | 中国科学院电工研究所 | 塔式太阳能热发电用单元镜面形检测方法及检测系统 |
| JP2014157106A (ja) * | 2013-02-18 | 2014-08-28 | Kobe Steel Ltd | 形状測定装置 |
| CN106768892A (zh) * | 2016-12-28 | 2017-05-31 | 中国计量大学 | 基于哈特曼‑夏克波前传感器的自由曲面镜片波面拼接方法 |
-
2017
- 2017-07-28 CN CN201710635739.3A patent/CN107560564B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101650157A (zh) * | 2009-09-18 | 2010-02-17 | 中国科学院长春光学精密机械与物理研究所 | 双曲面凸面反射镜面形误差的检测方法及其装置 |
| CN101943559A (zh) * | 2010-08-12 | 2011-01-12 | 中国科学院光电技术研究所 | 利用三坐标测量机检测大口径非球面光学元件的方法 |
| CN102297663A (zh) * | 2011-04-26 | 2011-12-28 | 北京理工大学 | 基于彩色编码逆向哈特曼式面形测量方法 |
| JP2014157106A (ja) * | 2013-02-18 | 2014-08-28 | Kobe Steel Ltd | 形状測定装置 |
| CN103267495A (zh) * | 2013-04-24 | 2013-08-28 | 中国科学院电工研究所 | 塔式太阳能热发电用单元镜面形检测方法及检测系统 |
| CN106768892A (zh) * | 2016-12-28 | 2017-05-31 | 中国计量大学 | 基于哈特曼‑夏克波前传感器的自由曲面镜片波面拼接方法 |
Non-Patent Citations (2)
| Title |
|---|
| DAODANG WANG ETC.: "《Computer-aided high-accuracy testing of reflective surface with reverse Hartmann test》", 《OPTICS EXPRESS》 * |
| 王道档 等: "《基于逆向哈特曼测试光学偏折技术的高精度面形检测方法》", 《第十六届全国光学测试学术交流会摘要集》 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108507495A (zh) * | 2018-03-19 | 2018-09-07 | 中国计量大学 | 一种基于逆向哈特曼检测的自由曲面检测方法 |
| CN108436602A (zh) * | 2018-03-30 | 2018-08-24 | 中国石油大学(华东) | 一种利用未修整近圆弧截面砂轮磨削自由曲面的方法 |
| CN109307480A (zh) * | 2018-09-30 | 2019-02-05 | 中国计量大学 | 一种透射元件多表面面形检测方法 |
| CN109739019A (zh) * | 2019-03-14 | 2019-05-10 | 南开大学 | 基于矢量像差理论的共轴自由曲面光学系统优化设计方法 |
| CN110245439A (zh) * | 2019-06-20 | 2019-09-17 | 成都飞机工业（集团）有限责任公司 | 一种扩口导管的安装误差建模和补偿方法 |
| CN110793440B (zh) * | 2019-11-01 | 2021-05-18 | 中国计量大学 | 一种光学偏折瞬态测量方法 |
| CN110763159A (zh) * | 2019-11-01 | 2020-02-07 | 中国计量大学 | 一种光学偏折显微表面测量装置及方法 |
| CN110793440A (zh) * | 2019-11-01 | 2020-02-14 | 中国计量大学 | 一种光学偏折瞬态测量方法 |
| CN110763159B (zh) * | 2019-11-01 | 2021-08-06 | 中国计量大学 | 一种光学偏折显微表面测量装置及方法 |
| CN110645916A (zh) * | 2019-11-14 | 2020-01-03 | 北京航空航天大学 | 基于基准标定板校正位姿的自由曲面测量方法与装置 |
| CN110702009A (zh) * | 2019-11-22 | 2020-01-17 | 中国科学院长春光学精密机械与物理研究所 | 一种基于逆向Hartmann计算机辅助法的三维测量系统 |
| CN114353699A (zh) * | 2022-01-07 | 2022-04-15 | 中国科学院长春光学精密机械与物理研究所 | 大陡度凸面光学自由曲面高频段像差检测系统及检测方法 |
| CN114353699B (zh) * | 2022-01-07 | 2022-12-16 | 中国科学院长春光学精密机械与物理研究所 | 大陡度凸面光学自由曲面高频段像差检测系统及检测方法 |
| CN114593693A (zh) * | 2022-02-16 | 2022-06-07 | 苏州英示测量科技有限公司 | 光学偏折检测方法、电子设备及光学偏折检测系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107560564B (zh) | 2019-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107560564B (zh) | 一种自由曲面检测方法及系统 | |
| Huang et al. | Review of phase measuring deflectometry | |
| Zhao et al. | Performance analysis and evaluation of direct phase measuring deflectometry | |
| CN109307480B (zh) | 一种透射元件多表面面形检测方法 | |
| CN108507495B (zh) | 一种基于逆向哈特曼检测的自由曲面检测方法 | |
| Zhao et al. | Calibration for stereo vision system based on phase matching and bundle adjustment algorithm | |
| CN102183213B (zh) | 一种基于相位测量偏折术的非球面镜检测方法 | |
| CN108225737A (zh) | 一种新型光学透镜测量方法 | |
| CN106840027A (zh) | 光学自由曲面的像散补偿型干涉检测装置与检测方法 | |
| CN108759798A (zh) | 一种高精度航天器精度测量的实现方法 | |
| CN102494634A (zh) | 一种基于条纹反射的离轴非球面镜检测方法 | |
| CN102200432A (zh) | 非球面物体测量方法和设备 | |
| CN107957251A (zh) | 基于计算机辅助校正的反射球面通用化检测方法 | |
| Wiegmann et al. | Accuracy evaluation for sub-aperture interferometry measurements of a synchrotron mirror using virtual experiments | |
| CN107543683A (zh) | 一种透射元件像差的高精度大动态范围测量系统及测量方法 | |
| CN103353388A (zh) | 一种具摄像功能的双目体式显微成像系统标定方法及装置 | |
| Li et al. | Monocular-vision-based contouring error detection and compensation for CNC machine tools | |
| CN113899321B (zh) | 凹面镜辅助成像的镜面物体三维形貌测量方法及系统 | |
| CN104406770B (zh) | 波像差测量模块的畸变测量装置和畸变校正方法 | |
| Wang et al. | Systematic error control for deflectometry with iterative reconstruction | |
| Kudo et al. | Simulation-based systematic error compensation for nanoprofiler using normal vector tracing method | |
| Zhang et al. | Novel method of positioning optical freeform surfaces based on fringe deflectometry | |
| CN113280755B (zh) | 基于曲面屏相位偏折的大曲率镜面三维形貌测量方法 | |
| Nguyen et al. | Real-time 3D measurement of freeform surfaces by dynamic deflectometry based on diagonal spatial carrier-frequency pattern projection | |
| Chen et al. | A method of sub-aperture slope stitching for testing flat element based on phase measuring deflectometry |
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 |