圆度误差评定中最小区域法的计算机叠代算法.docVIP

圆度误差评定中最小区域法的计算机叠代算法 回转体零件，其横截面轮廓是否为一正圆，需要与一理想圆进行比较才能得出结论，圆度误差的评定过程就是将被测横截面的实际轮廓与理想圆比较的过程。最小条件法是圆度误差评定中寻找理想圆的基本原则，目前经常采用的方法有最小二乘法、最小外接圆法、最大内切圆法和最小区域法。其中，最小区域法是一种新的优良评定方法，它不仅可以获得最小的误差评定结果，而且对零件的性质有稳定的约束 (通过包边界)，因而是现代测量技术致力研究的评定方法。在传统圆度测量仪中，实现最小区域圆评定的方法是测量仪通过传感器描绘出被测工件的轮廓误差曲线，然后测量人员用同心圆模板来试凑包容轮廓误差曲线，直到符合最小区域圆条件为止。当用计算机对传统圆度仪进行改造时，首先要解决的就是各种评定方法的算法，我们用计算机叠代法推导了最小区域圆的算法。 1 零件坐标系的建立 当被测零件放在测量仪器的工作台上时，零件圆心不能与工作台圆心绝对重合，所以零件坐标系不等于仪器坐标系，在实际测量中，测量点的坐标值是仪器坐标值，而圆度误差测量的基准参考系应是零件基准坐标系。因此，必须将测量点的仪器坐标值转换成零件坐标值。最小区域圆的评定过程就是根据基于仪器坐标的实际测量值，寻找满足最小区域圆条件的零件坐标原点，并将仪器坐标值转换成零件坐标值的过程。 1.1 最小区域圆条件 最小区域圆是指用两同心圆包容基于仪器坐标的实际测量曲线，且两同心圆与测量曲线应至少有内外交替的四点接触，满足此条件的两同心圆半径差为最小，两圆的圆心为满足最小区域圆条件的零件坐标系原点，如图 1 所示。两圆的半径差为被测工件的圆度误差。 1.2 确立零件坐标原点 1.2.1确立初始零件坐标原点 零件坐标系应基于仪器坐标系，在测量点均匀分布的情况下，可采用最小二乘法确立零件初始坐标原点，在测量非均匀分布的情况下，可取前3点求圆确定其圆心为零件初始坐标原 点。设有 n 个测量点 Pi (i=1，2，#8943;，n)，基于仪器坐标值为(Xi，Yi )，所求初始零件坐标原点基于仪器坐标系的坐标值为(X，Y)，则根据最小二乘圆原理可求得X= 2/ nΣXiY= 2/ nΣYi 我们可以将 (X，Y) 作为最初零件坐标原点基于仪器坐标的坐标值，但该值不一定满足最小区域圆条件，它只能作为寻找满足最小区域圆条件的坐标原点的初值。 1.2.2用叠代法求取理想坐标原点 所谓理想坐标原点就是满足最小区域圆条件的两同心圆的圆心，根据最小区域圆的条件，两同心圆的半径差应为最小，所以应使圆心向减少外圆半径和增加内圆半径的方向移动。设与外圆接触的 Pmax 点到圆心的距离矢量为 Rmax，与内圆接触的Pmin 点到圆心的距离矢量为 Rmin，则应使圆心分别向 Rmax 方向和 Rmin 相反的方向移动。设移动步长为 eL，eL 的确定过程如下。 移心方向确定后，移心步长便可由公式确定。如图 2 所示，两同心圆的圆心由 O 移到 O′，被测轮廓上任一被测点 Pi 至两同心圆的中心距离由 ri 变为 ri′，αi 为OPi 与OOcent; (n)的夹角，则 设 PL 为延续外接触点，移心后，外接圆半径 rL 变为 rL′，由式 (1) 得 式中 rL——移心前 PL 点所在外接圆的半径被测圆上任意一点与外接触圆接触的条件为 将式 (1)、(2) 代入式 (3)，便得到使 Pi 点成为外接触点所需移动步长 eLi 的计算公式 对应测量轮廓曲线上的 n 个点，可求出 n 个 eL，其中最小的 eL 即为最佳移动步长。移动后的坐标原点为X ′= X eLcosamax (5) Y′= Y eLsinamax 式中amax——最大距离矢量 Rmax与 X 轴正方向的夹角 2 求取圆度误差 2.1 最小区域圆条件判断算法 根据以上的步长和移动方向求出的圆心还必须满足最小区域圆的另一个条件，即被测轮廓曲线必须有四点分别交替与两同心圆接触。满足此条件的算法为与外圆相接触的两点连线和与内圆相接触的两点连线的交点在轮廓曲线之内。若不满足此条件，则反复叠代，直至找到为止。 2.2 圆度误差的求取 若满足以上条件，则求取两线段的中垂线交点，该点即为符合最小区域圆要求的零件坐标原点。此时可求取零件的圆度误差Fmz=Rmax-Rmin 式中Rmax，Rmin ——分别为轮廓曲线上到零件坐标原点的最大值、最小值，即外圆半径、内圆半径 3 数据处理实例分析 我们用 C 语言完成了以上算法，对石家庄金刚集团所生产的康明斯牌系列发动机上的 6BT 活塞销，在三坐标测量机上的测量数据进行处理。该活塞销的技术参数为 f40±0.0032，圆度公差 f 为 0.0025，圆度公差等级为 5 级，表面