发明内容
为了克服上述技术缺点,本发明的目的是提供了一种基于最小区域的零件平面度误差评定方法,从而提高几何产品测量仪器的平面度误差检测精度。
为达到上述目的,本发明采取的技术方案包含有如下步骤:
步骤1:将被测平面置于测量平台上,在测量空间直角坐标系中测量并获取平面上的点
,
=1, 2,…
n,
n为测点数目且
n为大于3的正整数,所有测点
形成测点集
;任意选取不在同一直线上的3个测点,计算过3个点的平面的法向矢量
,作为被测平面对应的理想平面法向矢量的迭代初值。
步骤2:依次计算测点集
中各测点到理想平面的距离
,并分别记录测点集
中各测点到理想平面的最大距离与最小距离所对应的测点;所记录的最小距离对应的测点形成误差包容区域的低值接触点集合
,所记录的最大距离对应的测点形成误差包容区域的高值接触点集合
。
步骤3:判断集合
中是否只有1个低值接触点且集合
中有1个高值接触点;
如果只有1个高值接触点且只有1个低值接触点,计算低值接触点与高值接触点连线的方向矢量
,测点的投影矢量
=
,跳转到步骤7;
如果低值接触点的数量大于1或高值接触点的数量大于1,跳转到步骤4。
步骤4:判断集合
中是否只有2个高值接触点且集合
中是否只有1个低值接触点;
如果只有2个高值接触点且只有1个低值接触点,测点的投影矢量
等于2个高值接触点连线的方向矢量,跳转到步骤7;
如果高值接触点的数量不为2或低值接触点的数量不为1,则跳转到步骤5。
步骤5:判断集合
中是否只有2个低值接触点且集合
中是否只有1个高值接触点;
如果只有2个低值接触点且只有1个高值接触点,则测点的投影矢量
等于2个低值接触点连线的方向矢量,跳转到步骤7;
如果低值接触点的数量不为2或高值接触点的数量不为1,跳转到步骤6。
步骤6:以2个高值接触点为1个组合,以其中1个组合中的2个高值接触点为计算对象,计算过这2个高值接触点且法向矢量垂直于
的平面
,并判断低值接触点和其余高值接触点是否分别分布在平面
的两侧;遍历所有的高值接触点组合,如果存在某一组合满足低值接触点与其余高值接触点分别分布在对应平面
的两侧,则测点的投影矢量
等于该组合对应的2个高值接触点连线的方向矢量,跳转到步骤7,否则继续执行该步的下面操作;
以2个低值接触点为1个组合,以其中1个组合中的2个低值接触点为计算对象,计算过这2个低值接触点且法向矢量垂直于
的平面
,并判断高值接触点和其余低值接触点是否分别分布在平面
的两侧;遍历所有的低值接触点组合,如果存在某一组合满足高值接触点与其余低值接触点分别分布在对应平面
的两侧,则测点的投影矢量
等于该组合对应的2个低值接触点连线的方向矢量,跳转到步骤7,否则跳转到步骤16。
步骤7:进行坐标变换,使坐标系的
z轴与
平行;被测平面对应的理想平面在
xoy坐标平面内投影为理想直线
,计算理想直线
的方向矢量
;坐标变换后,设测点
在
xoy坐标平面内投影为
,所有的测点
对应的
组成测点集合
。
步骤8:在测点集合
中,剔除坐标值完全相等的测点元素。
步骤9:依次计算测点集合
中各个测点到理想直线
的距离
,并记录并分别记录测点集
中各测点到理想直线
的最大距离和最小距离所对应的测点,所记录的最小距离对应的测点形成误差包容区域的低值接触点集合
,所记录的最大距离对应的测点形成误差包容区域的高值接触点集合
。
步骤10:判断集合
中是否只有1个低值接触点且集合
中是否只有1个高值接触点;
如果只有1个高值接触点与1个低值接触点,则它们为有效接触点
、
,跳转到步骤13;
如果高值接触点的数量大于1或低值接触点的数量大于1,则跳转到下一步。
步骤11:计算集合
中所有高值接触点的横坐标或纵坐标最大值
和最小值
,计算集合
中所有低值接触点的横坐标或纵坐标最大值
和最小值
;
如果
大于等于
,且
小于等于
,则表明符合直线度判别准则,跳转到步骤15;否则,执行步骤12。
步骤12:如果
小于
,则
对应的高值接触点为有效高值接触点
,
对应的低值接触点为有效低值接触点
;如果
大于
,则
对应的高值接触点为有效高值接触点
,
对应的低值接触点为有效低值接触点
;其余接触点为无效接触点。
步骤13:判断
与
的点积是否小于0,如果小于0,将-
赋值给
;此时包容区域的旋转变动的方向矢量
等于
;以过有效接触点
、
的连线为界将
xoy坐标平面划分为
、
2个区域,设
、
的坐标分别为
、
,根据下式可以判断测点
所在的区域,
步骤14:依次计算
区内各个测点
与
的连线的方向矢量
,调整
使
与
的点积为正,并计算
和
的夹角
,
区内所有测点
对应的
组成集合
;依次计算
区内各个测点
与
的连线方向矢量
,调整
使
与
的点积为正,并计算方向矢量和
的夹角
,
区内所有测点
对应的
组成集合
;计算2个集合
、
中最小的夹角,即为包容区域旋转变动量
S;包容区域绕
旋转
S角度后得到理想直线
的方向矢量
;跳转到步骤9。
步骤15:进行步骤7的逆向坐标变换,还原测点的坐标值,同时根据
计算被测平面对应理想平面的法向矢量
;跳转到步骤2。
步骤16:输出平面度误差和对应理想平面的最优参数。
本发明依据平面自身的特征,在平面度误差优化的方向上,给出了包容区域旋转变动的方式,并设计了具体的计算方法。本发明所公布的方法,可计算出平面度误差的最优值和实际平面对应包容平面参数的最优值。
本发明设计了,从步骤9到步骤14,还可以单独用来计算基于最小区域的给定平面内直线度误差。
本发明设计了,对于面对面垂直度误差以及面对线平行度误差,首先进行坐标旋转变换,使基准面的法向矢量或基准线的方向矢量平行于z正,同时被测面的测点作相应的刚体旋转,则被测特征的包容平面在xoy平面变为两个平行线。然后应用给定平面内直线度误差计算的方法,获得面对面垂直度误差以及面对线平行度误差。
本发明还设计了,对于平行平面构成的内表面或外表面,对本发明所公布的方法进行稍加改动,即可计算平行平面的体外作用尺寸或体内作用尺寸。
与现有技术相比,本发明公布了一种基于最小区域法的零件平面度误差评定方法,该方法首先测量并获取被测平面上的测点坐标;然后给出平面的初始参数;根据点到平面的距离,查询与误差包容区域接触的测点,根据接触点的数量与相对位置,应用坐标系变换的方式,将测点投影指定平面内,转化为给定平面内的直线度误差计算;在投影平面内,通过包容区域的旋转变动,搜索与包容区域接触的、满足直线度误差判别准则的测点,获得旋转后的包容区域;然后重新确定投影平面,再按直线度误差计算,依次迭代计算,直到满足平面度判别准则,输出平面度误差以及对应理想平面参数的最优值。本发明可准确计算出满足最小区域判别准则的平面度误差。
本发明的有益效果在于:本方法查询到的接触点相对位置满足最小区域判别准则,计算得到的实际平面对应理想平面的参数和平面度误差均为最优值,计算稳定性好、计算效率高。
具体实施方式
以下结合附图1对本发明的实施方式进行详细说明。
步骤1:将被测平面置于测量平台上,在测量空间直角坐标系中测量并获取平面上的点
,
=1, 2,…
n,
n为测点数目且
n为大于3的正整数,所有测点
形成测点集
;任意选取不在同一直线上的3个测点,计算过3个点的平面的法向矢量
,作为被测平面对应的理想平面法向矢量的迭代初值。
步骤2:依次计算测点集
中各测点到理想平面的距离
,并分别记录测点集
中各测点到理想平面的最大距离与最小距离所对应的测点;所记录的最小距离对应的测点形成误差包容区域的低值接触点集合
,所记录的最大距离对应的测点形成误差包容区域的高值接触点集合
。
步骤3:判断集合
中是否只有1个低值接触点且集合
中有1个高值接触点;
如果只有1个高值接触点且只有1个低值接触点,计算低值接触点与高值接触点连线的方向矢量
,测点的投影矢量
=
,跳转到步骤7;
如果低值接触点的数量大于1或高值接触点的数量大于1,跳转到步骤4。
步骤4:判断集合
中是否只有2个高值接触点且集合
中是否只有1个低值接触点;
如果只有2个高值接触点且只有1个低值接触点,测点的投影矢量
等于2个高值接触点连线的方向矢量,跳转到步骤7;
如果高值接触点的数量不为2或低值接触点的数量不为1,则跳转到步骤5。
步骤5:判断集合
中是否只有2个低值接触点且集合
中是否只有1个高值接触点;
如果只有2个低值接触点且只有1个高值接触点,则测点的投影矢量
等于2个低值接触点连线的方向矢量,跳转到步骤7;
如果低值接触点的数量不为2或高值接触点的数量不为1,跳转到步骤6。
步骤6:以2个高值接触点为1个组合,以其中1个组合中的2个高值接触点为计算对象,计算过这2个高值接触点且法向矢量垂直于
的平面
,并判断低值接触点和其余高值接触点是否分别分布在平面
的两侧;遍历所有的高值接触点组合,如果存在某一组合满足低值接触点与其余高值接触点分别分布在对应平面
的两侧,则测点的投影矢量
等于该组合对应的2个高值接触点连线的方向矢量,跳转到步骤7,否则继续执行该步的下面操作;
以2个低值接触点为1个组合,以其中1个组合中的2个低值接触点为计算对象,计算过这2个低值接触点且法向矢量垂直于
的平面
,并判断高值接触点和其余低值接触点是否分别分布在平面
的两侧;遍历所有的低值接触点组合,如果存在某一组合满足高值接触点与其余低值接触点分别分布在对应平面
的两侧,则测点的投影矢量
等于该组合对应的2个低值接触点连线的方向矢量,跳转到步骤7,否则跳转到步骤16。
步骤7:进行坐标变换,使坐标系的
z轴与
平行;被测平面对应的理想平面在
xoy坐标平面内投影为理想直线
,计算理想直线
的方向矢量
;坐标变换后,设测点
在
xoy坐标平面内投影为
,所有的测点
对应的
组成测点集合
。
步骤8:在测点集合
中,剔除坐标值完全相等的测点元素。
步骤9:依次计算测点集合
中各个测点到理想直线
的距离
,并记录并分别记录测点集
中各测点到理想直线
的最大距离和最小距离所对应的测点,所记录的最小距离对应的测点形成误差包容区域的低值接触点集合
,所记录的最大距离对应的测点形成误差包容区域的高值接触点集合
。
步骤10:判断集合
中是否只有1个低值接触点且集合
中是否只有1个高值接触点;
如果只有1个高值接触点与1个低值接触点,则它们为有效接触点
、
,跳转到步骤13;
如果高值接触点的数量大于1或低值接触点的数量大于1,则跳转到下一步。
步骤11:计算集合
中所有高值接触点的横坐标或纵坐标最大值
和最小值
,计算集合
中所有低值接触点的横坐标或纵坐标最大值
和最小值
;
如果
大于等于
,且
小于等于
,则表明符合直线度判别准则,跳转到步骤15;否则,执行步骤12。
步骤12:如果
小于
,则
对应的高值接触点为有效高值接触点
,
对应的低值接触点为有效低值接触点
;如果
大于
,则
对应的高值接触点为有效高值接触点
,
对应的低值接触点为有效低值接触点
;其余接触点为无效接触点。
步骤13:判断
与
的点积是否小于0,如果小于0,将-
赋值给
;此时包容区域的旋转变动的方向矢量
等于
;以过有效接触点
、
的连线为界将
xoy坐标平面划分为
、
2个区域,设
、
的坐标分别为
、
,根据下式可以判断测点
所在的区域,
步骤14:依次计算
区内各个测点
与
的连线的方向矢量
,调整
使
与
的点积为正,并计算
和
的夹角
,
区内所有测点
对应的
组成集合
;依次计算
区内各个测点
与
的连线方向矢量
,调整
使
与
的点积为正,并计算方向矢量和
的夹角
,
区内所有测点
对应的
组成集合
;计算2个集合
、
中最小的夹角,即为包容区域旋转变动量
S;包容区域绕
旋转
S角度后得到理想直线
的方向矢量
;跳转到步骤9。
步骤15:进行步骤7的逆向坐标变换,还原测点的坐标值,同时根据
计算被测平面对应理想平面的法向矢量
;跳转到步骤2。
步骤16:输出平面度误差和对应理想平面的最优参数。
以下通过实验实例,说明本发明所公布方法计算的有效性。
实验对象1:基于最小区域的给定平面内直线度误差计算
表1 被测直线上的坐标值
测点编号 |
x
|
y
|
z
|
测点编号 |
x
|
y
|
z
|
1 |
154.777 |
208.938 |
-608.996 |
5 |
154.253 |
160.097 |
-608.995 |
2 |
154.708 |
202.788 |
-608.995 |
6 |
154.094 |
143.383 |
-608.996 |
3 |
154.604 |
192.826 |
-608.995 |
7 |
153.876 |
123.145 |
-608.996 |
4 |
154.463 |
179.178 |
-608.995 |
8 |
153.666 |
102.531 |
-608.996 |
在一平面内测量一条直线,共采集8个离散点,测点坐标如表1。测量完成后,在xoy坐标平面内评定直线度误差。本发明所公布方法的计算结果为0.015mm,计算效果参照附图2所示。
由附图2可知,第5号测点为高值接触点(如图中★所示),第1、6号测点为低值接触点(如图中●所示),评定结果符合最优解判别准则。计算得到的最优参数为:平面内直线的方向矢量(0.0104, 0.9999),直线上定点为(154.1735, 151.7400)。
实验对象2:基于最小区域的平面度误差计算
表2 平面1
测点编号 |
x
|
y
|
z
|
测点编号 |
x
|
y
|
z
|
1 |
166.279 |
357.705 |
-671.135 |
10 |
166.365 |
308.424 |
-651.961 |
2 |
166.335 |
328.207 |
-671.135 |
11 |
166.321 |
332.966 |
-651.961 |
3 |
166.372 |
303.4 |
-671.135 |
12 |
166.274 |
361.014 |
-651.961 |
4 |
166.419 |
283.707 |
-671.134 |
13 |
166.277 |
361.013 |
-624.217 |
5 |
166.46 |
265.444 |
-671.134 |
14 |
166.325 |
337.502 |
-621.291 |
6 |
166.495 |
243.694 |
-671.134 |
15 |
166.372 |
310.84 |
-621.291 |
7 |
166.489 |
243.693 |
-651.96 |
16 |
166.41 |
289.931 |
-621.711 |
8 |
166.449 |
266.452 |
-651.961 |
17 |
166.452 |
263.576 |
-621.711 |
9 |
166.414 |
290.133 |
-651.961 |
18 |
166.504 |
241.446 |
-621.711 |
在一个平面内采集18个点,坐标值如表2所示。海克斯康三坐标测量仪评定的结果为0.016mm,本发明所公布方法的计算结果为0.0152 mm,优于三坐标测量机评定结果。
计算的效果如图3所示,第3、17号测点为高值接触点(如图中●所示),第5、9号测点为低值接触点(如图中★所示),符合最优解判别准则。计算得到的最优参数为:平面上一点(167.2993, -155.3035, -125.4085),平面的法向矢量(0.9999, 0.0019, -0.0001)。