基于最小二乘法的平面任意位置椭圆轮廓度误差的精确计算

　　0 引言

　　平面曲线轮廓，如渐开线、椭圆和摆线轮廓等在工程中被广泛应用，因此曲线轮廓度测量、识别和误差评定成为轮廓度测量的重要内容[1]。如为了使发动机活塞能与汽缸良好地贴合，提高发动机热能效应，常将活塞截面加工成中凸变椭圆的几何形状[2]。这种中凸变椭圆的活塞裙部形线较为复杂、轮廓加工精度要求较高。因此，研究椭圆轮廓度误差评价方法对保证活塞及其他椭圆轮廓零件的加工质量有重要的意义。

　　关于椭圆轮廓度误差的评定，国家标准尚未给出明确的定义，也没有给定一种特定的评定算法求解椭圆轮廓度误差。近年来国内外学者开展这方面的研究较少，取得了一些研究成果。比较有代表性的成果有刘书桂等[3]基于最小二乘原理的椭圆误差的评价方法。邹益民等[4]基于几何距离的拟合算法。侯宇等[5]采用有效集法进行数据处理。陈基伟[6]提出的椭圆直接拟合算法。T.S.R.Murthy[7]提出的正交最小二乘法、二项高斯分布法和基于代数距离的最小二乘法来评定椭圆轮廓度误差。这些评定方法对于椭圆误差的评定都有一定的作用和效果。

　　本文通过平面任意位置椭圆方程、椭圆法线方程特点和椭圆本身性质，利用最小二乘原理实现对平面任意位置椭圆轮廓度误差的评定。

　　1 最小二乘椭圆拟合

　　设平面任意位置椭圆(如图1所示)方程为：

　　设为椭圆轮廓上的N(N≥5)个测量点，依据最小二乘原理，所拟合的目标函数为：

　　欲使F为最小，须使

　　由此可得正规方程：

　　解方程(3)可得到A, B ,C , D和E的值，由参考文献[3]可得到拟合出的最小二乘椭圆的五个主参数：位置参数(θ, x , y)及形状参数(a , b )。

　　2 确定过测量点的椭圆法线与最小二乘椭圆的交点

　　如图1所示，设点Mi(Xi、Yi)为过测量点Pi(Xi、Yi)的椭圆法线与最小二乘椭圆的交点，由于交点Mi(Xi、Yi)既在法线上又在最小二乘椭圆[9]上，即点M(Xi、Yi)满足方程组：

　　3 判断各测量点ip的位置

　　设拟合出的最小二乘椭圆的两个焦点为(如图1所示)。依据椭圆的性质可得：

　　测量点pi点到两焦点的距离之和ddi为：

　　圆弧上的速度v=4.995。

　　与求在可行路径1P上的最短时间一样，可求得在可行路径2P上的最短时间是T=96.56322秒。

　　所以从O→A的最短时间T=min{T1,T2}=94.228312秒，最短时间路径如上所述。

　　4 结束语