三坐标测量机的测头半径补偿与曲面匹配李　春　刘书桂(天津大学精密测试技术与仪器国家重点实验室　天津　300072)

摘要　在非均匀双三次B—样条函数的基础上,导出自由曲面任意点的法矢量通用算法,进而提出自由曲面测头半径补偿公式;为了更好的消除自由曲面测量中的定位误差,提出了应用单纯形法,对测量原始点进行坐标平移和旋转变换,从而较好的解决了曲面匹配问题。关键词　自由曲面　测头补偿　曲面匹配

TheProbeRadiusCompensationofFree-formSurfaceandSurfaceMatchingLi　Chun　LiuShugui(StateKeyLaboratoryofPrecisionMeasuringTechnologyandInstrument,TianjinUniversity,Tianjin300072,China)

Abstract　Basedonnon-uniformB-splines,anewcurrentalgorithmwithnormalvectorofrandomfree-formsurface'spointisdeduced,andmore,aformulawithprobecompensationisproposed.Weofferaarithmeticnamedsimplexmethodinordertoeliminatingorientationerrorintheprocessoffree-formsurfacemeasurement.Itcansettlesurfacematchingwellbyshiftingandrotatingthemeasuringcoordinatesystem.Keywords　Free-formsurface　Probecompensation　Surfacematching

1　引 言三坐标测量机由于其测量精度和智能化程度较高,广泛应用于制造业的CAD/CAM、产品检测和质量控制[1]。用三坐标测量机的球形测头测量自由曲面时,得到的数据是测头中心轨迹,由于测头总有一定的半径r,因此测得的是与被测曲面相距r的包络面。为了得到所需的测量表面,需要求出球心轨迹面所构成的包络面,这个过程被称为测头半径补偿。在实际测量过程中,并不能做到实际曲面和标准曲面完全重合,需要将被测曲面进行旋转、平移等坐标变换,使被测曲面与标准曲面大致重合,从而达到曲面检测的目的,这个过程称之为曲面匹配。2　测头半径补偿方法用球形测头测量曲面时,测头与被测曲面为点接触,测头半径补偿的关键是确定曲面在接触点处的法矢。球测头与被测曲面接触时,球心一定在被测点的法线上,而且被测点一定在球心轨迹面过球心点的法线上。因此不论能否得知被测面的法线方向或是球心面的法线方向,都能对测头半径进行补偿。本文提出了一种新方法,不在测量过程中补偿测头半径,而只是收集测头中心坐标值,然后应用曲面建模理论,计算出球心各点的法矢量值,继而补偿测头半径。(1)自由曲面的偏导数求法首先,根据三坐标测量机所得的原始测量点,我们可以反求出双三次B—样条自由曲面的模型[2]:

S(u,v)=∑ni=0∑mj=0Ni,4(u)Nj,4(v)Pi,j(1)其中Ni,4(u),Nj,4(v)为双三次B—样条基函数,Pi,j为控制预点。先求曲面沿u向的切矢量,即对S(u,v)求偏导:

Su(u,v)=󰀁󰀁uS(u,v)

=∑mj=0Nj,4(v)󰀁󰀁u∑ni=0Ni,4p(u)Pi,j

第24卷第4期增刊 仪　器　仪　表　学　报 2003年8月=∑mj=0Nj,4(v)󰀁󰀁uCj(u)(2)这里,Cj(u)=∑ni=0Ni,4p(u)Pi,j,j=0,…,m为B样条曲线。

而C′(u)=∑n-1i=04Ni,3(u)Pi+1-Piui+4-ui+1,所以:Su(u,v)=∑n-1i=0∑mj=04Ni,3(u)Nj,4(v)Pi+1,j-Pi,jui+4-ui+1(3)同样,我们可以得出:Sv(u,v)=∑ni=0∑m-1j=04Ni,4(u)Nj,3(v)Pi,j+1-Pi,jvi+4-vi+1(4)Su(u,v)和Sv(u,v)分别为曲面上的点沿u向和v向的切矢量。(2)曲面的测头半径补偿公式被测曲面与测头中心轨迹曲面是法向等距面关系。测头中心轨迹曲面上的任意点处的单位矢量可以得出:n󰀁(u,v)=Su(u,v)×Sv(u,v)󰀁Su(u,v)×Sv(u,v)󰀁(5)其中Su(u,v)和Sv(u,v)可以由式(3)(4)得出。根据测头半径值r,可以推出被测实际曲面的补偿公式为:P(u,v)=S(u,v)󰀁rn󰀁(u,v)(6)当测头位于被测曲面法矢量所指的一侧时,取“-”号,反之,取“+”号。(3)计算机仿真结果利用解析曲面进行数字仿真计算,用以考察所述方法的精度。为了方便计算,我们考虑一个椭球面,方程为x23002+y25002+z23002=1。在第一象限内非均匀的取点P′i,j,并用解析方法将其转换为法向等距面上的点Si,j,模拟测量数据。应用本论文所述方法对得到的模拟数据进行曲面拟合与测头半径补偿,得到生成曲面上的点Pi,j,将其与原始数据P′i,j比较,以偏差e=max󰀁P′i,j-Pi,j󰀁为指标。设测头半径r=5mm。表1　最大偏差e的计算结果nu/nvmu=10mu=20mu=300.61.26×10-44.28×10-52.24×10-50.78.25×10-53.57×10-51.75×10-50.82.23×10-51.09×10-56.31×10-60.9——1.58×10-58.29×10-6 设测量点数为mu×mv,令mu=mv=10,20,30,节点控制数nu=nv=(6/10󰀁9/10)mu,取得的偏差如下表,可以看出,偏差e随着测量点数的增大而减小,节点控制数也对偏差有一定的影响。由此可见,本方法具有一定的计算精度。3　曲面匹配曲面匹配是曲面误差评价的基础,在曲面检测和逆向工程中,经常要用到曲面匹配。曲面的匹配实际上是一种受约束的拟合,利用测量点进行拟合的结果是一个其形状与理想轮廓曲面完全相同的曲面。将测量所得的原始点进行适当的平移和旋转后,理论曲面与测量曲面将“贴合”得相当完美。(1)坐标系的“旋合”过程假定被测曲面存在M个原始理论点P′i(x′i,y′i,z′i)(i=1,2,…,M),被测曲面的拟合模型为:z=f(x,y)。将M个原始理论点P′i进行坐标平移和旋转变换,得出一组新的理论点值Pi(xi,yi,zi)(i=1,2,…,M)。xiyizi=A・x′iy′iz′i+x0

y0z0(7)

其中:A=cos󰀁cos󰀁+sin󰀁sin󰀁sin󰀁sin󰀁cos󰀁-cos󰀁sin󰀁+sin󰀁sin󰀁cos󰀁

-sin󰀁cos󰀁+cos󰀁sin󰀁cos󰀁sin󰀁cos󰀁cos󰀁-sin󰀁sin󰀁sin󰀁+cos󰀁sin󰀁cos󰀁cos󰀁cos󰀁x0,y0,z0是三个坐标平移量,󰀁,󰀁,󰀁是三个坐标旋

转量。当x0,y0,z0,󰀁,󰀁,󰀁恰好是理论坐标系与测量坐标系之间的旋转平移关系时,M个新理论点Pi(xi,yi,zi)与拟合模型z=f(x,y)最接近。再定义一个新函数:

F(x0,y0,z0,󰀁,󰀁,󰀁)=∑Mi=1(f(xi,yi)-zi

)2(8)

经过坐标变换后,M个新的理论点Pi与拟合模型z=f(x,y)最接近时,F(x0,y0,z0,󰀁,󰀁,󰀁)为最小,于是,

曲面匹配的问题变成了求:F(x0,y0,z0,󰀁,󰀁,󰀁)→min(9)时的x0,y0,z0,󰀁,󰀁,󰀁的值。(2)单纯形法的求解思路由于式(9)含有非线性运算(含有正弦、余弦),一般情况是先将它进行线性化,将其展开为泰勒级数形

146仪　器　仪　表　学　报 第24卷　式,舍去二次高阶项,这样比较复杂,这里我们应用最优化算法中的单纯形法进行求解。单纯形法是一种不必计算函数梯度的直接搜索法,它能自动适应函数的局部环境变化,在极值附近内将会收缩。算法如下:1)令n维空间的单纯形的n+1个顶点为:x-i(i=1,2,…,n+1)计算函数值Fi=f(x-i)(i=1,2,…,n+1),比较大小,并确定:FR=F(x(R))=max1≤i≤n+1Fi(极大点)FG=F(x(G))=max1≤i≤n+1Fi(次极大点)FL=F(x(L))=max1≤i≤n+1Fi(极小点)2)求出最坏点x(R)的对称点　xN=2xF-x(R)式中:xF=1n∑n+1i=1(x(i)-x(R))。3)若F(xN)≥FR,则将xN缩小为xM,xM=(1-󰀁)xR+󰀁xN　(0<󰀁<1,󰀁≠0.5)如果F(xM)1)如果上述条件满足,并且F(xE)(6x2+y2・󰀁/150)+

cos(x2+y2・󰀁/50)]将其数据进行平移旋转重新得到的数据作为原始理论数据。原始理论数据坐标系及曲面坐标系之间的关系由以上算法计算出(见表2)。原始数据经过坐标平移旋转后,形成了在曲面坐标系中的理论数据。下表是在曲面坐标系下的理论数据与曲面的部分差值(见表3)。因为给定的曲面f(xi,yi)没有误差,所以最终的结果比较好,这也说明了算法具有较高的精确度。

4　结束语应用非均匀双三次B—样条函数对三坐标测量机所得到的测量点进行曲面拟合和测头半径三维补偿,生成和重建被测自由曲面。而单纯形法较好的解决了测量过程中难以准确定位的问题,将测量坐标系与理论坐标系统一起来。

表2　测量坐标系与理论坐标系之间的关系x0y0z0

󰀁󰀂

计算值-20.0-20.020.0-30.0°-10.0°-12.0°

表3　在测量坐标系下的误差xiyif(xi,yi)zif(xi,yi)-zi

-200.0-200.055.98835155.9883480.000003

-200.0-180.0-42.835883-42.835876-0.000007-200.0-160.0-112.730708-112.730707-0.000001-200.0-140.0-125.154689-125.1546900.000001

参考文献1　张国雄.三坐标测量机[M].天津:天津大学出版社,1999.2　朱心雄.自由曲线曲面造型技术.北京:科学出版社,2000.

3　卢红,等.测头半径补偿的方法.组合机床与自动化加工技术,2001,(9).4　高国军.检验自由曲面时精确定位方法的研究.机械科学与技术,1997,(6).5　于源,等.自由曲面测量中曲面匹配的建模及算法分析.机械科学与技术,2001,(5).