三坐标测量机测量平面度的样本点优化选取方法

摘要

本发明提供一种三坐标测量机测量平面度的样本点优化选取方法，首先，测量同一生产系统依次间隔抽样的零件确定系统误差模式。然后，将零件表面系统误差模式与禁忌搜索算法整合，改进传统禁忌搜索算法过分依赖初始解的缺点。最后，通过搜索生成三坐标测量机测量零件表面平面度时的测量样本点。本发明生成的测量样本点能够通过测试较少的样本点而快速的覆盖整个平面的极限点。利用选取出的样本点形成的测量样本点集来规划测量路径，可以有效的提高三坐标测量机的测量精度以及测量效率，更加精确地反映零件表面的平面度形貌和加工纹理，改进了传统测量样本点选取方法不能兼顾精度与效率的缺点以及禁忌搜索算法搜索样本点需要人为修正初始样本点的缺陷。

权利要求

1.一种三坐标测量机测量平面度的样本点优化选取方法，其特征在于,包括以下步骤：第一步：预检测

采用三维高分辨率表面形貌测量技术对多件同一生产系统依次生产的零件表面进行测量作为预检测，得到三维高密度点云数据，以及以彩色表示的表面形貌图像；

第二步：确定系统误差模式

根据点云数据与表面形貌图像，记录表面高度最高与最低的极限点出现的位置以及出现频率，统计结果后，出现最高与最低极限点次数较多的位置即是该零件的系统误差所在位置；在系统误差所在位置分别选取若干个样本点作为样本点集，计算样本点集内每一样本点的平面度偏差，将拥有最大平面度偏差的样本点作为之后搜索过程的初始样本点；

第三步：设定搜索过程参数

根据被测零件表面积设定搜索过程参数，包括搜索移动步长△S，禁忌列表容量C，最大妥协移动数CM，以及总移动数M；从初始样本点(x0,y0)点开始移动，则搜索可移动的方向有(x0+△s,y0),(x0-△s,y0),(x0,y0+△s),(x0,y0-△s)四个方向，其中，△S为搜索移动步长；将以上四步移动作为候选移动列表，分别对其进行平面度偏差计算，从中选取平面度偏差最大的位置点作为下一位置，同时当前点(x0,y0)写入禁忌列表，在一定禁忌表容量内不得出现；

第四步：妥协移动

当所有候选位置点平面度偏差都小于现有位置点的偏差时，依旧选取备选位置点中偏差较大的点作为下一步的位置点，将这一次移动称为妥协移动；第三步中设定的最大妥协移动次数为CM，当妥协移动次数大于CM时，表明搜索过程无法找到更大的平面偏差，因此最优样本点就是之前搜索过程中找到的最大偏差位置点，搜索结束；

第五步：搜索结束条件

当搜索移动次数达到第三步中设定的总移动数M时，搜索过程结束，最优样本点为搜索过程中平面度偏差最大的样本点。

2.根据权利要求1所述的三坐标测量机测量平面度的样本点优化选取方法，其特征在于，第一步中，三维高分辨率表面形貌测量技术横向与纵向最小测量间隔为0.15mm，垂直方向高度测量精度达1μm，测量范围达到500mm*300mm，测点密度达到每平方毫米40个测量点。

3.根据权利要求1所述的三坐标测量机测量平面度的样本点优化选取方法，其特征在于，第二步中，应用MATLAB软件记录平面度最大偏差样本点出现频率与位置，以及读入所需点云数据并且计算平面度偏差，确定制造系统的系统误差模式。

4.根据权利要求1所述的三坐标测量机测量平面度的样本点优化选取方法，其特征在于，第三步中，搜索移动步长△S小于被测平面边长的1/20且大于1mm；

禁忌列表容量C为3-5；

妥协移动数CM为3-10；

总移动数M根据被测平面表面积设定。

5.根据权利要求4所述的三坐标测量机测量平面度的样本点优化选取方法，其特征在于，应用MATLAB软件进行搜索过程模拟。

6.根据权利要求1所述的三坐标测量机测量平面度的样本点优化选取方法，其特征在于，第一步中，采用三维高分辨率表面形貌测量技术对八件以上同一生产系统依次生产的零件表面进行测量。

说明书

三坐标测量机测量平面度的样本点优化选取方法

技术领域

[0001] 本发明涉及零件表面平面度测量技术领域，具体地，涉及一种三坐标测量机测量平面度的样本点优化选取方法。

背景技术

[0002] 在精密制造过程中，对于零件的尺寸、形位精度要求越来越高。为了能够对零件尺寸、形位精度进行有效的测量与控制，提高测量结果的准确性与测量过程的效率显得尤为重要。作为零件重要的形位精度，表面平面度的制造精度要求越来越高，这对平面度的精确测量提出了更高的要求。其中，三坐标测量机是目前工业领域内最为常用的平面度测量设备。影响三坐标测量机测量平面度准确性的主要因素有两个：一是三坐标测量机本身的测量精度；二是测点样本选取是否能充分体现平面的实际特征。要得到精确的平面度测量结果，不仅需要测量精度高的测量仪器，也对测量方法尤其是样本点选取方法提出了更高的要求。因此，作为影响三坐标测量机测量平面度精度与测量效率的主要因素，三坐标测量机测量点选取优化方法研究至关重要。在测量过程中，测量样本点选取越多，覆盖范围越大，所得到的测量结果就越准确，但是考虑到增加样本点的数量会导致测量时间延长，影响测量效率，所以需要优化测量样本点的选取来得到较高的测量精度。

[0003] 已有的技术中，Badar在论文“Intelligent search-based selection of sample points for straightness and flatness estimation”(《ASME制造科学与工程杂志》2003年第125卷第2期，263-271页)中提到规划测量路径点时应当以最少的测量点保证最优的测量点数据分布，即效率与精度的平衡决策问题，提出针对测量直线度与平面度的三坐标测量机样本点缩减策略。该方法应用了传统的禁忌搜索算法和结合了传统禁忌搜索算法的混合搜索算法对最优测量样本点进行选取，由于传统禁忌搜索算法自身对初始解的要求过高，因此在无法提供优秀初始解时不能有效的优化测量样本点，必须通过人为修正初始解改善样本点搜索过程精度与效率。

[0004] 专 利 号 为US20130030773 名 称 为“COORDINATE MEASURING SYSTEM DATA REDUCTION”的美国专利，提出的基于零件几何特征信息与尺寸公差对零件进行分块分析缩减测量样本容量的方法，通过对零件CAD模型的预处理与分割，能够根据几何特征对零件的测量样本点进行缩减，对不同的分块采用不同的样本缩减方法进行处理。该技术虽然能够很好的适应不同几何形状与尺寸的零件，但是对于不同几何特征的分块进行样本缩减时没有将测量精度作为优化目标，因此不能保证缩减后的测量样本达到原始精度，特别是针对同一零件存在多个不同几何特征分块的情况时，不能有效地保证测量精度，违背了三坐标测量机样本点优化选取的目标。

发明内容

[0005] 针对现有技术中的缺陷，为了克服三坐标测量机测量零件表面平面度测量点选取方法不能同时满足测量精度与测量效率的不足，本发明提出了一种三坐标测量机测量平面度的样本点优化选取方法，该方法是基于改进禁忌搜索算法的三坐标测量机样本点优化选取方法，以测量精度与样本容量缩减作为优化目标，首先利用三维高分辨率表面形貌测量技术(水平分辨率可达150μm，垂直分辨率可达1μm)产生的三维高密度点云数据(每平方毫米可达40个测量点)寻找系统误差模式，选取出更加准确的初始样本点，解决了一般禁忌搜索算法搜索结果过分依赖于初始解的问题；之后将平面度偏差作为搜索目标，利用改进型禁忌智能搜索算法进行其余的样本点搜索，选取可以完全表征表面平面度变化的极限样本点。该方法可以对零件表面进行全局搜索，能够通过测试少量的样本点得到高精度的测量结果，实现保证测量精度的前提下缩减样本容量的目的，此外针对不同生产条件的零件该方法都具有优秀的适应性。

[0006] 本发明是按照下述技术方案实现的。本发明方法包括下述步骤：

[0007] 第一步：预检测

[0008] 采用三维高分辨率表面形貌测量技术对多件同一生产系统依次生产的零件表面进行测量作为预检测，得到三维高密度点云数据，以及以彩色表示的表面形貌图像；

[0009] 第二步：确定系统误差模式

[0010] 根据点云数据与表面形貌图像，记录表面高度最高与最低的极限点出现的位置以及出现频率，统计结果后，出现最高与最低极限点次数较多的位置即是该零件的系统误差所在位置；在系统误差所在位置分别选取若干个样本点作为样本点集，计算样本点集内每一样本点的平面度偏差，将拥有最大平面度偏差的样本点作为之后搜索过程的初始样本点；

[0011] 第三步：设定搜索过程参数

[0012] 根据被测零件表面积设定搜索过程参数，包括搜索移动步长△S，禁忌列表容量C，妥协移动数CM，以及总移动数M；从初始样本点(x0,y0)点开始移动，则搜索可移动的方向有(x0+△s,y0),(x0-△s,y0),(x0,y0+△s),(x0,y0-△s)四个方向，其中，△S为搜索移动步长；将以上四步移动作为候选移动列表，分别对其进行平面度偏差计算，从中选取平面度偏差最大的位置点作为下一位置，同时当前点(x0,y0)写入禁忌列表，在一定禁忌表容量内不得出现；

[0013] 第四步：妥协移动

[0014] 当所有候选位置点平面度偏差都小于现有位置点的偏差时，依旧选取备选位置点中偏差较大的点作为下一步的位置点，将这一次移动称为妥协移动；第三步中设定的最大妥协移动次数为CM，当妥协移动次数大于CM时，表明搜索过程无法找到更大的平面偏差，因此最优样本点就是之前搜索过程中找到的最大偏差位置点，搜索结束；

[0015] 第五步：搜索结束条件

[0016] 当搜索移动次数达到第三步中设定的总移动数M时，搜索过程结束，最优样本点为搜索过程中平面度偏差最大的样本点。

[0017] 优选地，第一步中，三维高分辨率表面形貌测量技术横向与纵向最小测量间隔为0.15mm，垂直方向高度测量精度达1μm，测量范围达到500mm*300mm，测点密度达到每平方毫米40个测量点。

[0018] 优选地，第二步中，应用MATLAB软件记录平面度最大偏差样本点出现频率与位置，以及读入所需点云数据并且计算平面度偏差，确定制造系统的系统误差模式。

[0019] 优选地，第三步中:

[0020] 搜索移动步长△S小于被测平面边长的1/20且大于1mm；

[0021] 禁忌列表容量C为3-5；

[0022] 妥协移动数CM为3-10；

[0023] 总移动数M根据被测平面表面积设定。

[0024] 优选地，应用MATLAB软件进行搜索过程模拟。

[0025] 优选地，第一步中，采用三维高分辨率表面形貌测量技术对八件以上同一生产系统依次生产的零件表面进行测量。

[0026] 本发明提供一种三坐标测量机测量平面度的样本点优化选取方法，本发明的有益效果具体为：本发明生成的测量样本点，能够通过测试较少的样本点而快速的覆盖整个平面的极限点，利用选取出的样本点形成的测量样本点集来规划测量路径，可以有效的提高三坐标测量机的测量精度，更加精确地反映出零件表面的平面度形貌和加工纹理，改进了传统测量样本点选取方法不能兼顾精度与效率的缺点，以及禁忌搜索算法搜索样本点需要人为修正初始样本点的缺陷。同时，由于利用三维高分辨率表面形貌测量技术，保证了零件表面系统误差所在位置的精确定位，为后续解决系统误差模式的质量控制与改进加工方法提供了先决条件，具有重要的工程实用价值。

附图说明

[0027] 通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：

[0028] 图1为本发明提供的三坐标测量机测量平面度的样本点优化选取方法的方法流程图；

[0029] 图2为本发明根据三维高分辨率表面形貌测量技术测得的发动机缸体顶面表面形貌图像确定系统误差模式示意图；

[0030] 图3为本发明三维高密度点云数据经过按区块简化预处理后的发动机缸体顶面高度坐标图像示意图。

具体实施方式

[0031] 下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明的保护范围。

[0032] 请参阅图1，一种三坐标测量机测量平面度的样本点优化选取方法，包括以下步骤：

[0033] 第一步：预检测

[0034] 采用三维高分辨率表面形貌测量技术对多件同一生产系统依次生产的零件表面进行测量作为预检测，得到三维高密度点云数据，以及以彩色表示的表面形貌图像。

[0035] 第二步：确定系统误差模式

[0036] 根据点云数据与表面形貌图像，记录表面高度最高与最低的极限点出现的位置以及出现频率，统计结果后，出现最高与最低极限点次数较多的位置即是该零件的系统误差所在位置；在系统误差所在位置分别选取若干个样本点作为样本点集，计算样本点集内每一样本点的平面度偏差，将拥有最大平面度偏差的样本点作为之后搜索过程的初始样本点。

[0037] 第三步：设定搜索过程参数

[0038] 根据被测零件表面积设定搜索过程参数，包括搜索移动步长△S，禁忌列表容量C，妥协移动数CM，以及总移动数M；从初始样本点(x0,y0)点开始移动，则搜索可移动的方向有(x0+△s,y0),(x0-△s,y0),(x0,y0+△s),(x0,y0-△s)四个方向，其中，△S为搜索移动步长；将以上四步移动作为候选移动列表，分别对其进行平面度偏差计算，从中选取平面度偏差最大的位置点作为下一位置，同时当前点(x0,y0)写入禁忌列表，在一定禁忌表容量内不得出现。

[0039] 第四步：妥协移动

[0040] 当所有候选位置点平面度偏差都小于现有位置点的偏差时，依旧选取备选位置点中偏差较大的点作为下一步的位置点，将这一次移动称为妥协移动；第三步中设定的最大妥协移动次数为CM，当妥协移动次数大于CM时，表明搜索过程无法找到更大的平面偏差，因此最优样本点就是之前搜索过程中找到的最大偏差位置点，搜索结束。

[0041] 第五步：搜索结束条件

[0042] 当搜索移动次数达到第三步中设定的总移动数M时，搜索过程结束，最优样本点为搜索过程中平面度偏差最大的样本点。

[0043] 优选地，第一步中，三维高分辨率表面形貌测量技术横向与纵向最小测量间隔为0.15mm，垂直方向高度测量精度达1μm，测量范围达到500mm*300mm，测点密度达到每平方毫米40个测量点。

[0044] 优选地，第二步中，应用MATLAB软件记录平面度最大偏差样本点出现频率与位置，以及读入所需点云数据并且计算平面度偏差，确定制造系统的系统误差模式。

[0045] 优选地，第三步中，

[0046] 搜索移动步长△S小于被测平面边长的1/20且大于1mm；

[0047] 禁忌列表容量C为3-5；

[0048] 妥协移动数CM为3-10；

[0049] 总移动数M根据被测平面表面积设定。

[0050] 优选地，应用MATLAB软件进行搜索过程模拟。

[0051] 优选地，第一步中，采用三维高分辨率表面形貌测量技术对八件以上同一生产系统依次生产的零件表面进行测量。

[0052] 具体地，下面以某汽车发动机厂生产的某型号直列四缸发动机缸体顶面的测量样本点选取过程为例，结合附图对本发明的具体实施作进一步描述。

[0053] 本实施例中以最大平面度负偏差为例说明搜索过程，最大平面度正偏差操作过程完全相同。为简化搜索过程方便演示，对示例所测量的缸体顶面点云数据进行了分区块简化预处理，缩减了数据量，不影响方法展示。

[0054] 如图1、图2、图3所示，一种三坐标测量机测量平面度的样本点优化选取方法，包括下述步骤：

[0055] 第一步：如图1所示，按照流程图进行预检测。采用三维高分辨率表面形貌测量技术对八件同一生产系统生产、依次间隔抽样额零件表面进行测量作为预检测，得到三维高密度点云数据，以及用彩色表示的表面形貌图像。其中，发动机顶面横向与纵向最小测量间隔为0.15mm，垂直方向高度测量精度可达1μm，测量点总数可达80万。

[0056] 第二步：确定系统误差模式。根据点云数据与表面形貌图像，记录表面高度最高与最低的极限点出现的位置以及出现频率，统计结果后，出现最高极限点次数最多的位置如图2所示第一列和第三列线圈标注的位置，出现最低极限点次数最多的位置为第二列和第四列线圈标注的位置。其他出现极限点较多的位置包括3号与4号缸孔间的右侧部位，以及1号缸孔左侧上方。在上述位置共选取5个样本点作为样本点集，如图3所示，分别计算每一样本点的平面度偏差，将拥有最大平面度负偏差的样本点(-125，-240)作为搜索过程的初始样本点。

[0057] 第三步：设定搜索步长为10mm，禁忌列表容量为3，妥协移动数3以及总移动数10。从初始样本点(-125,-240)点开始移动，则搜索可移动的方向有(-115,-240),(-135,-240),(-125,-230)，(-125，-250)四个方向。以上四步移动作为候选移动列表，分别对其进行平面度偏差计算，从中选取平面度偏差最大的位置点(-135，-240)作为下一位置，同时当前点(-125,-240)写入禁忌列表，在3次移动内不得出现。

[0058] 第四步：重复迭代进行第三步搜索过程，经过5次搜索移动，其中存在3次所有候选样本点的平面度偏差都不优于当前位置平面度偏差的情况，所以妥协移动达到3次，触发搜索结束条件，搜索过程停止。

[0059] 第五步：最终所得的平面度最大负偏差样本点为(-135，-250，-84.9187)。同样条件下以最大正偏差作为搜索目标进行搜索，从初始解(-125，-300)开始搜索，经过4次搜索后妥协移动达到3次，终止搜索，所得的平面度最大正偏差样本点为(-125，-290，-84.8716)。利用所选取出的样本点集进行测量精确度达到91.67％。

[0060] 以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影响本发明的实质内容。