轴承零件圆轮廓最小二乘圆误差的Newton迭代法修正

　　在轴承套圈沟道或钢球等圆轮廓的接触式测量中,对采样点进行评定时,最小二乘圆法是目前通用的评定方法,此方法计算简便但存在一定的误差。为使评定结果更加精确,更逼近真实值,现采用Newton迭代法对最小二乘圆法评定结果进行修正。

　　1 圆轮廓的最小二乘圆通用算法

　　在轴承零件圆轮廓的实际测量过程中,得到的一系列测量点(xi,yi),按照圆的标称形式进行最小二乘圆的优化计算,可以得到圆轮廓半径和圆心的坐标[1-2]。

　　如图1所示,传感器实际测量出来的一系列测量点组成样本集(xi,yi), i=1,2,,,N。

　　样本集(xi,yi)中点到圆心的距离为di,则:

　　采用最小二乘法直接分析圆的轮廓,在对各未知参数求偏导后,无法求解,为了计算方便,可以取点(xi,yi)到圆心距离与圆半径的平方差:

　　对于所有参与评定的采集数据点,须使得最小,令Q(a,b, c)=,对a,b, c求偏导,并等于零,即:

　　但用最小二乘圆法对圆轮廓的采样点进行评定时,不同的仪器各自存在不同的问题。

　　(1)对于圆度仪、圆柱度仪和沟曲率仪等仪器采用极坐标对圆轮廓进行测量评定时,采用的是采样点到圆心的距离与圆半径的平方差的平方,而不是两者差的平方,因此增加了距离与半径之和的权重,计算结果则必然使半径变小,圆心位置发生偏移。

　　(2)对于轮廓仪等仪器采用直角坐标对圆轮廓进行测量评定时,采用的是采样点到圆边缘的距离;而根据最小二乘法的定义应该是实测值与相应圆上点的理论值之间的误差。因此,以此距离代替误差必然会带来分析结果的误差。

　　2 Newton迭代法对最小二乘圆误差的修正

　　现针对上述两种问题,分别提出不同的修正方法。

　　2.1 问题1的修正

　　根据圆的方程,点(xi,yi)到圆心的距离与半径的差为:

　　对于所有参与评定的采集数据点,须使得最小,令Q(A,B,R)=,使分别对A,B,R求偏导,并等于零,即:

　　按照Newton迭代法[3],则有:

　　取(6)~(8)式中的(A,B,R)为计算的初值(A0,B0,R0),根据测量需要的精度,计算(ΔA,ΔB,ΔR)即可。

　　2.2 问题2的修正

　　根据圆的方程R²=(xi-A)²+(yi-B)²有:

　　在直角坐标系中,误差为:

　　对于所有参与评定的采集数据点,须使得最小,令Q(A,B,R)=,使分别对A,B,R求偏导,并等于零。

　　2.2.1 球体的修正

　　测量球体等凸面向上的圆形轮廓时: