计算机控制光学表面成形技术的驻留时间算法

吴清飞+任志英+高诚辉+林春生

摘要：

为了提高镜片的加工精度与效率，利用计算机控制光学表面成形技术（CCOS）的抛光方法对光学镜片进行抛光全过程动态仿真。根据Preston方程建立材料去除函数模型，对抛光过程中压力、转速以及工件与抛光磨头相对半径比对抛光去除速率的影响进行分析。为建立球面镜片的动态全过程仿真，结合卷积原理，推导加工残余误差与去除函数和驻留时间三者间的线性关系，根据镜片的对称性，将元素个数从2m+1点简化为m+1点，以提高运算效率。最后为获得仿真最小残余误差，采用非负最小二乘法求解驻留时间。结果表明，材料去除速率函数类似于高斯分布，抛光后能使镜片面形误差收敛，对模拟表面进行仿真，半径为100 mm的镜片其初始表面形貌粗糙度的均方根值从0.467 μm收敛到0.028 μm，轮廓最大高度从6.12 μm收敛到1.48 μm。对实测表面进行加工仿真同样令其表面形貌粗糙度的均方根值从3.007 μm收敛到0.107 μm，轮廓最大高度从160.73 μm收敛到13.76 μm，因此提出的驻留时间求解方法对于球面镜片抛光全过程动态仿真有一定的可行性。

关键词：

光学镜片； 球形磨头； 抛光去除； 动态仿真； 快速迭代法

中图分类号： TH 161文献标志码： Adoi： 10.3969/j.issn.10055630.2017.04.007

Abstract：

In order to improve the precision and efficiency of the lens manufacture，this paper describes that the optical lenses using computer controlled optical surfacing technology（CCOS） performs the whole process of dynamic simulation.The material removal function model is established according to Preston equation to analyze the influence of pressure，rotating speed and the ratio of radius of the workpiece and the spherical polishing tool over the material removal rate.To establish a dynamic simulation of the whole process of the lens，considering the principle of convolution to derive a linear algebraic equation of the residual error，removal function and dwell time and its symmetry，we can reduce the number of points from to that improves operational efficiency.Finally，a nonnegative least squares method gives a solution to minimize residual error.The results showed that the material removal rate was similar to a Gaussian distribution function.It can make that the lens surface shape error converges after polishing.Structure surfaces radius is 100mm，and its initial rootmeansquare is 0.467 μm converging to 0.028 μm.Peak value converges from 6.12 μm to 1.48 μm.The simulation of the measured surface also make its initial RMS from 3.007 μm converge to 0.107 μm and PV converge from 160.73 μm to 13.76 μm.Thus，the proposed dwell time solving method is feasible for whole process of dynamic simulation of spherical polishing.

Keywords：

optical lens； spherical polishing tool； polishing removal； dynamic simulation； fast iterative algorithm

引言

光學镜片具有高分辨率及良好的光学性能，其广泛应用在民用市场和军用市场中且需求量仍在不断增加。光学镜片的质量评定都是在最后一道工序抛光工艺结束后进行的，因此抛光工艺会直接影响到光学镜片的表面质量[1]。目前光学镜片的抛光方法主要有传统抛光法、化学抛光法、气囊抛光法[2]和磁流变抛光法[34]等。随着近半个世纪计算机技术的迅猛发展，诞生了以美国为首的一些发达国家率先发展起来的计算机控制光学表面成形技术（computer controlled optical surfacing，CCOS），其工作原理为采用一个比被抛光元件尺寸小得多的抛光磨头（抛光磨头的直径一般小于工件直径的1/4），而后通过计算机控制抛光磨头相对于工件的压力、速度以及驻留时间从而实现对工件表面的抛光加工[5]。endprint

Aspden于20世纪80年代率先对材料去除机理与数学模型进行研究，而后Jones又在Aspden的数学模型基础上，创新性地提出了采用卷积迭代法的驻留时间计算模型，为该技术奠定了重要理论基础[67]。Cheung等[8]利用Preston公式基于平面建立仿真模型，通过研究材料去除速率和驻留时间函数来预测构造的结构表面抛光后的表面形貌特征，并通过实验论证了预测结果的合理性。周旭升[9]基于平面工件建立了双旋转抛光磨头的去除函数模型，通过研究面形误差收敛的优化控制方法，提出了基于加工时间的脉冲迭代法和基于加工精度的脉冲迭代法，并分析了抛光盘的偏心率和转速比等对去除函数修形能力的影响。李全胜等[10]基于对光学非球面元件模型的研究，提出了求解CCOS中驻留时间的三种算法：基于加工仿真的驻留时间算法，基于傅里叶变换（FFT）的驻留时间算法，基于FFT与滑动平均的算法，并且应用滑动平均的算法对非球面模型进行加工仿真，运算效率比较高且获得的驻留时间分布平缓，加工后残余误差较小。罗丽丽等[11]对大口径的光学元件驻留时间的求解提出了基于大规模稀疏矩阵的非负最小二乘法，并以此方法对典型的320 mm×320 mm的微浮雕面形进行仿真加工，结果表明，该算法不仅可以提高运算效率而且能够使面形误差收敛至0.067λ（λ=632.8 nm）。

综上所述，目前国内外很多学者采用CCOS抛光技术对平面或非球面进行了大量的研究工作，但据多家企业调查发现，中小型非球面的主要加工方式是在研抛成球面镜片表面上镀上一层特殊的光学树脂，因此球面镜片需求量依旧很大，研究球面镜片的抛光过程具有重要的意义。本文基于球面光学元件，采用CCOS抛光方法依据Preston方程推导材料去除函数模型，最后进行模拟加工，验证了模型的合理性。

1.1CCOS工作原理

在抛光过程中有众多因素会对光学镜片的表面材料去除造成影响，例如相对压力、相对速度、环境温度和抛光液及其浓度等。现在描述CCOS表面材料去除最成功的模型仍是Preston在1927年提出的著名的Preston方程。即

dzdt=KPV

（1）

式中：dz/dt为单位时间内光学镜片表面的去除量；K为比例常数，由抛光过程中除压力与速度外的其他影响因素决定，如抛光液的材料、浓度、环境温度、湿度等；P为表面接触区域的瞬时接触压力；V为某一点的瞬时相对抛光速度。

因此，当加工时间为t时，工件的表面材料去除量

1.2去除函数推导

如图1所示为半径为Rp的球形抛光磨头对半径为r的球形工件抛光的几何示意图。其中球形工件的转速为ω1，球形抛光磨头的转速为ω2。点O为球形工件的中心，O1为球形抛光磨头的中心。

根据赫兹接触理论[12]可知兩弹性球体间的接触区域为圆形区域。基于赫兹假设，不考虑两弹性球体间的接触摩擦力，可计算在两球体间的接触压力为

P（x，y）=P01-x2+y2a232

（4）

式中：P0为抛光磨头中心点的压强；a为抛光接触区域圆的半径。

如图2所示为球形磨头与工件表面在抛光接触时某一点的速度分解示意图。其中v1为工件自转产生的线速度，v2为球形抛光磨头转动产生的线速度，θ为两者之间的夹角。因为接触区域的半径a很小，且接触区域近似看作平面，所以在接触区域上各点相对于抛光磨头转动中心的距离相等。

2影响去除函数的因素

2.1压力对去除速率的影响

根据Preston方程，当其他抛光条件都不变时，抛光过程中材料去除量与压力成正比，即材料去除量随着压力的增大而增大。在抛光过程中，当抛光磨头施加在工件表面的力保持不变时，接触区域上的压力分布最大值与接触区域的大小有关，根据赫兹接触理论，压力分布服从式（4）。其中中心点处压力P0为

P0=5F2πa2

（8）

式中：a=3π（k1+k2）rRp4（r+Rp）1/3，k1=1-μ21πE1，k2=1-μ22πE2，μ1、μ2、E1、E2分别为工件与抛光磨头的泊松比和弹性模量，r、Rp分别为工件与抛光磨头的半径；F为抛光磨头作用于工件的法向量作用力。

当保持工件与磨头的材料不变且抛光磨头施加在工件表面上的力F也不变时，由于球面半径与球形磨头半径比值不同而引起压力的变化曲线如图3所示，压力的变化表现为中心大并向四周递减的趋势。从图4可以看出，当半径比从0.1逐渐增加到0.2时，压力P的峰值急剧减小；而随着半径比的继续增大，压力P峰值的减小也越来越缓慢。因此在相同的抛光条件下，随着工件半径与磨头半径比值的增大，会导致材料去除率的峰值随之减小。但在实际的抛光过程中，随着半径比的减小，会导致接触区域的减小，因此应合理的选择工件与球形磨头的半径比，才能获得更好的抛光去除效率。

2.2工件与磨头转速比对去除速率的影响

根据Preston方程，抛光时材料去除量和工件与球形磨头的相对速度成正比，即材料去除量随着相对速度的增大而增大。

图5为不同转速比下去除速率的变化曲线，从图中可以看出当保持其他抛光条件不变时，令球形磨头转速不变，改变工件的转速，可以看出工件材料的去除速率随着转速比ω1/ω2的降低而减小。当转速比小于1/10时，材料去除速率的减小也越来越慢，说明当转速比小于1/10时，工件的转速对材料去除的影响比较小。

图6所示为在以球面圆心为坐标原点建立的直角坐标系y=a/2处，工件材料去除速率曲线变化图。因为合速度是由工件转速与球形磨头转速矢量合成的，当工件转速与磨头转速的作用相互抵消时，会造成在该加工点的材料去除速率为零。在该转速比下，越往外合速度越大，而压力分布正好相反，因此产生的去除函数速率曲线呈倒W型，如图6（b）。endprint

3驻留时间的求解

由CCOS抛光技术原理可知，工件表面的材料去除量是由去除函数和驻留时间来控制的，因此通过控制抛光磨头在各驻留点的驻留时间可以精确控制工件的表面材料去除，从而达到对工件表面的面形误差进行修正。因此，求解驻留时间函数是CCOS抛光技术的关键步骤之一。对于三维表面形貌的评定，其基准面的提取是表面评定的关键技术，目前国内外基准提取常用的主要方法有最小二乘多项式拟合法、高斯滤波法和小波滤波法[1314]。本文采用高斯滤波法来提取加工前后的表面形貌基准。

3.1材料去除向量

当磨头处于ri时，抛光区域的大小为半径为a的圆，设r和θ为建立在该抛光区域内的极坐标值。其轨迹示意图如图7所示。

3.3仿真结果分析

3.3.1

构造表面仿真

应用分形几何中的WM函数来构造待抛光的粗糙表面，获得1 024×1 024采樣点，球面半径为100 mm的初始表面如图9所示。计算后可得其初始面形误差均方根（RMS）值为0.467 μm，最大峰值（PV）为6.12 μm。采用本文研究的算法，运用MATLAB对工件进行仿真加工来验证算法的正确性。得到抛光后的球面工件如图10所示，从图中可看出，仿真加工后的面形误差得到收敛，计算得最终的面形误差RMS值为0.028 μm，PV为1.48 μm。

3.3.2

实测表面仿真

如图11所示为利用美国KLATencor公司的MicroXam100型光学轮廓仪进行测量所得的实际光学元件表面，该仪器的最小采样间距为0.078 μm，最大量程为1 000 mm，分辨率为0.01 nm。运用MATLAB采用高斯滤波法计算得其初始面形误差RMS值为3.007 μm，PV为160.7 μm。通过MATLAB软件对其进行模拟仿真加工后得到表面如图12所示，经计算其RMS值为0.107 μm，PV为13.36 μm。综上所述，不管是对仿真构造表面还是实际表面，本文所述的仿真迭代法都能使面形误差得到收敛，效果较理想，从而证明本模型的有效性与合理性。

3.4运算效率比较

本文基于其对称性，计算时将元素个数从2m+1减为m+1后采用非负最小二乘法求解驻留时间。对同等大小的工件，保证其他条件都相同，取不同的采样点数可获得其运算时间随采样点数变化的曲线如图13所示。从图中可以看出，当采样点数小于2 048时，其运算时间都较短，但当采样点数大于2 048后，运算时间近似于指数式增长。如当点数为8 192时，运算耗费的时间达到1 643 s，而在点数减半为4 096时所花费的时间仅为210 s，因此采用对称性将运算时的元素个数减半能大大缩短运算时间。

4结论

本文基于旋转对称的球形工件采用CCOS抛光技术对其进行模拟仿真加工，通过分析材料去除模型，求解驻留时间分布函数来进行研究。综合以上分析，本文结论可归纳如下：

（1） 抛光磨头只需径向进给，且其去除函数形状类似于高斯分布，能较好的使面形误差得到收敛，具有较高的加工效率，最终能够有效的减小工件的表面残余误差；

（2） 加工过程中，球形磨头与工件的接触区域面积较小，能够减小磨头移动到工件边缘时因压力的突变而引起的边缘效应；

（3） 根据工件的对称性，采用合并运算的方式，有效的提高了运算效率。最终在得到材料去除函数和残余误差后，应用非负最小二乘法求得驻留时间分布。模拟加工结果表明，该算法能使面形误差得到收敛，从而证明了模型的正确性与合理性。

参考文献：

[1]杨力.先进光学制造技术[M].北京：科学出版社，2001.

[2]SONG J F，YAO Y X，XIE D G，et al.Effects of polishing parameters on material removal for curved optical glasses in bonnet polishing[J].Chinese Journal of Mechanical Engineering，2008，21（5）：2933.

[3]YANG Z M，HU H P，GAO Y B.The research in algorithm for MRF forming dwell time of optical components[J].Advanced Materials Research，2010，102104：750753.

[4]SINGH A K，JHA S，PANDEY P M.Mechanism of material removal in ball end magnetorheological finishing process[J].Wear，2013，302（1/2）：11801191.

[5]袁巨龙，张飞虎，戴一帆，等.超精密加工领域科学技术发展研究[J].机械工程学报，2010，46（15）：161177.

[6]ASPDEN R，MCDONOUGH R，NITCHIE F R.Computer assisted optical surfacing[J].Applied Optics，1972，11（12）：27392747.

[7]JONES R A.Computer controlled optical surfacing with orbital tool motion[J].Optical Engineering，1986，25（6）：256785.

[8]CHEUNG C F，KONG L B，HO L T，et al.Modelling and simulation of structure surface generation using computer controlled ultraprecision polishing[J].Precision Engineering，2011，35（4）：574590.

[9]周旭升.大中型非球面计算机控制研抛工艺方法研究[D].长沙：国防科学技术大学，2007.

[10]李全胜，成晔，蔡复之，等.计算机控制光学表面成形驻留时间算法研究[J].光学技术，1999（3）：5659.

[11]罗丽丽，何建国，王亚军，等.计算机控制光学表面成形中大规模驻留时间求解[J].强激光与粒子束，2011，23（12）：32073212.

[12]YANG M Y，LEE H C.Local material removal mechanism considering curvature effect in the polishing process of the small aspherical lens die[J].Journal of Materials Processing Technology，2001，116（2/3）：298304.

[13]REN Z Y，GAO C H，HAN G Q，et al.DTCWT robust filtering algorithm for the extraction of reference and waviness from 3D Nano scalar surfaces[J].Measurement Science Review，2014，14（2）：8793.

[14]REN Z Y，GAO C H，SHEN D.3D profile filter algorithm based on parallel generalized Bspline approximating Gaussian[J].Chinese Journal of Mechanical Engineering，2015，28（1）：148154.

（编辑：张磊）endprint