[发明专利]一种计算球体零件作用尺寸的方法

说明书

技术领域

本发明涉及一种计算球体零件作用尺寸的方法，属于精密计量与计算机应用领域，可用于各种情况下球形几何产品的合格性检测，并为加工过程及加工工艺的改进提供指导。

背景技术

球形零件在航空航天设备以及精密仪器中具有广泛应用，球度误差具有重要的研究意义。最大内接球法，是计算实际球面轮廓的半径最大的内接球面，最大内接球体现了被测内球面所能包容的最大配合外球面，由此获得的球度误差可视为被测内球面与最大配合外球面之间的最大间隙；最小外接球法，是计算实际球面轮廓的半径最小的外接球面，最小外接球体现了被测外球面所能通过的最小配合内球面，由此获得的球度误差可视为被测外球面与最小配合内球面之间的最大间隙。最大内接球和最小外接球统称为球体的作用表面，最大内接球的半径和最小外接球的半径统称为球体的作用尺寸。

球体的作用尺寸计算，属于不可微复杂最优化问题，目前，国内外学者主要采用传统优化方法、智能算法、计算几何方法等。这些方法由于存在计算稳定性差、计算效率低、对采点数量有限制、计算结果难以达到精确解等缺陷，导致最大内接球法与最小外接球法很难在实际检测中应用。目前市场上一般都采用成熟的最小二乘法近似地计算零件球度误差以及球体的作用尺寸。

发明内容

为了克服上述技术缺点，本发明的目的是提供一种计算球体零件作用尺寸的方法。本方法不仅提高了测量仪器测量球体作用尺寸的精度，而且算法稳定性好、计算效率高，可以推广应用于其几何体作用尺寸的计算。

本发明一种计算球体零件作用尺寸的方法，用于计算球体的最大内接尺寸时，主要包括以下步骤：

步骤1：将被测球置于测量平台上，在测量空间直角坐标系中测量并获取球表面上的点                                                ，=1, 2,…, n代表测点数目且n为大于4的正整数，所有测点形成测点集；取不在同一平面上的4个测点，计算4点组成的四面体外接球的球心，作为计算的迭代初始值。

步骤2：依次计算测点集中各测点到球心的距离，并记录各测点距离球心的最小距离，最小距离为被测球体的半径R，所记录的最小距离对应的测点组成误差包容区域内边界的接触点集合。

步骤3：判断包容区域内边界的接触点的数量是否等于1；

如果接触点的数量等于1，则该接触点为有效接触点，此时包容区域平移的方向矢量等于球心坐标减去接触点坐标；

如果接触点的数量不等于1，则跳转到步骤4。

步骤4：判断包容区域内边界的接触点的数量是否等于2；

如果接触点的数量等于2，则2个接触点为有效接触点，此时包容区域平移的方向矢量等于球心坐标减去2个接触点连线的中点坐标；

如果接触点的数量不等于2，则跳转到步骤5。

步骤5：判断包容区域内边界的接触点的数量是否等于3；

如果接触点的数量等于3，则3个接触点为有效接触点，此时包容区域平移的方向矢量等于球心坐标减去3个接触点的外接圆心坐标；

如果接触点的数量不等于3，则跳转到步骤6。

步骤6：判断包容区域接触点内边界的数量是否大于等于4；

如果接触点的数量大于等于4，以3个接触点为1个组合，以其中1个组合中的3个接触点为计算对象，计算过3个接触点的平面，判断球心与其余接触点是否分别分布在该平面的两侧；遍历所有的组合，如果所有组合都不满足球心与其余接触点分别分布在对应平面的两侧，则表示满足判别准则，跳转到步骤8；如果存在某一组合满足球心与其余接触点分别分布在对应平面的两侧，则该组合对应的3个测点为有效接触点，则此时包容区域平移的方向矢量等于球心坐标减去3个有效接触点外接圆心坐标；

如果接触点的数量不等于4，则跳转到步骤2。

步骤7：依次计算包容区域变动到内边界与各个非接触测点接触时，球心虚拟位置；

首先计算测点和其中1个有效接触点的垂直中分面，然后计算垂直中分面与过球心且方向矢量为的直线的交点，即为包容区域变动到内边界与测点接触时球心的位置，遍历所有的非接触测点，计算对应的，并计算到球心的距离，得到集合；在集合中，剔除与异向的对应元素，然后查询集合中的最小值，即为包容区域的平移变动量；根据包容区域的平移变动量以及移动方向矢量，计算球心的坐标；跳转到步骤2。

步骤8：输出球度误差、球的球心坐标以及球的半径R。