平面度误差评定方法及装置的制造方法
【专利摘要】本发明涉及一种平面度误差评定方法及装置，其中，该方法包括：初始化旋转角和当前的最佳平面度误差；将零件平面的被测点集投影至旋转角下的投影面；确定投影面内投影点的二维凸包；确定二维凸包内的最小点线距离；判断是否满足最小区域条件；更新所述当前的最佳平面度误差；判断是否满足收敛条件，以确定最终的平面度误差评定结果。本发明计算量小，通过采用凸包快速生成方法生成投影面上投影点集的二维凸包，易于实施。同时，通过改变投影面的旋转角度得到不同的二维凸包以从中找到符合最小区域条件的解，执行步骤简单且易于编程实现，因此容易在实际工程中应用。
【专利说明】
平面度误差评定方法及装置
技术领域
[0001] 本发明涉及精密计量与计算机应用领域，尤其涉及一种平面度误差评定方法及装 置。
【背景技术】
[0002] 平面度误差是评价机械零件的设计、制造、装配是否满足要求的一项重要指标。国 标GB/T11337-2004中的平面度误差的评定方法有:最小包容区域法、最小二乘法、对角线平 面法和三远点平面法。其中最小包容区域法的评定结果小于或等于其他三种评定方法。目 前应用较多的评定方法为最小二乘法和最小包容区域法。最小二乘法虽然从原理上不满足 最小区域条件，但其编程易实现，且计算速度快，可以近似地计算机械零件的平面度误差。 而最小包容区域法虽然原理上符合公差定义，但满足最小包容区域法的平面度评定问题属 于不可微的最优化问题，计算复杂。平面度误差评定的技术难点是如何快速又准确的计算 平面度误差。
[0003] 现有技术中公开了一种基于改进人工蜂群算法的平面度误差最小包容区域的评 定方法，该方法给出了符合最小区域条件的平面度误差评定数学模型。该方法虽然可以用 于平面度误差的评定，但是其建立的数学模型包含两个平面方程参数，并且将平面度误差 问题转化成了二元函数寻优问题，因此该方法计算效率低，限制了其在实际工程中的应用。 此外，现有技术中还公开了一种基于几何搜索逼近的平面度误差最小区域评定算法，该方 法阐述了利用几何优化搜索算法求解平面度误差的过程和步骤，并给出了基于辅助平面的 数学模型。该方法虽然可以用于平面度误差的评定，但是其给出的辅助点及辅助平面的构 造过程繁琐，并且数学模型复杂。
【发明内容】
[0004] 针对现有平面度误差的评定方法计算效率低，在实际工程中的应用受限，数学模 型复杂等缺陷，本发明提出一种平面度误差评定方法，该方法包括：
[0005] SI:初始化旋转角和当前的最佳平面度误差；
[0006] S2:将零件平面的被测点集投影至所述旋转角下的投影面；
[0007] S3:确定所述投影面内所述投影点的二维凸包；
[0008] S4:确定所述二维凸包内的最小点线距离；
[0009] S5:根据所述最小点线距离判断是否满足最小区域条件；
[0010] S6:在满足最小区域条件时，更新所述当前的最佳平面度误差；
[0011] S7:根据更新的所述当前的最佳平面度误差判断是否满足收敛条件，以确定最终 的平面度误差评定结果。
[0012] 可选地，所述方法还包括：
[0013] 在不满足最小区域条件时，更改所述旋转角，并重复执行步骤S2至S7。
[0014] 可选地，所述步骤Sl中初始化旋转角和当前的最佳平面度误差，包括：
[0015] 在测量空间直角坐标系O-XYZ中测量并获取零件平面的被测点点集P，将所述点集 P内第i个点的坐标表示为P1U1, yi，Zl)，并过Z轴建立旋转平面S，所述平面S绕z轴以Θ作为 旋转角进行旋转，将θ〇作为旋转角的迭代初值，并令当前的最佳平面度误差fbd为0;
[0016] 其中，所述坐标系O-XYZ包括X轴、y轴和z轴，i = l，2,…，η，η为被测点的数目，且η 为大于等于3的正整数，所述平面S与平面XOZ的夹角为Θ，且Θ的取值区间为[0，ji]。
[0017] 可选地，所述S2中将零件平面的被测点集投影至所述旋转角下的投影面，包括：
[0018] 将O-XYZ坐标系绕z轴逆时针旋转Θ角，以构建新坐标系O-X7 V Z7，
[0019] 其中，所述新坐标系O-X7 Y' Z7包括X7轴、y'轴和y轴，所述坐标系O-X7 V Z7的y轴 与V轴位于所述平面S内，平面V Of与所述平面S重合；
[0020] 在所述新坐标系O-XY Z7中，将所述被测点坐标表示为P' i，/ ，并将新 坐标系O-XH^下的被测点Iy i(x\，/i，z\)与原始坐标系O-XYZ下的被测点Pi(XiJiJi) 的转换关系表示为下式：
[0021] (X7ilY7i,z/i,l) = (xi,yi,zi,l) · Rz
[0022] 以所述平面S作为投影面，将新坐标系0-f疒Z'下的被测点点集K i投影至旋转投 影面S，得到平面X7OZ'内的投影点集P〃，并将所述点集P〃内第i个点的坐标表示为PWY Z%);且将所述点集P〃内第i个点的坐标可以表示为PW ;
[0023] 其中，其中Rz为坐标系O-XYZ绕z轴逆时针旋转后各个被测点坐标的变换矩阵，且Rz 的表达忒为.
[0024]
[0025] X' i = xi · cos〇-yi · sm〇,z' i = Zi〇
[0026] 可选地，所述S3中确定所述投影面内所述投影点的二维凸包，包括：
[0027]采用凸包快速生成法生成所述平面fOZ'内的投影点集P〃的一个二维凸包Q，所述 二维凸包Q的边数和顶点数均为Il1为二维凸包Q的一条边，p〃4Pp〃1+1为线段^的两个端 点，其中? = 1，2,3···πι。
[0028] 可选地，所述步骤3中采用凸包快速生成法生成所述平面fOZ'内的投影点集P〃的 一个二维凸包Q，包括：
[0029] 在平面X7 OZ7内的投影点集P〃中取两个极端点Pdr和Pul，其中，Pdr为最右最下的极 端点，Pul为最左最上的极端点；
[0030] 采用经过Pdr和Pul两点的直线PdrPul将整个点集P〃分为左区域和右区域，并分别对 所述左区域和右区域进行递归，以获取整个点集s的凸包。
[0031] 可选地，所述分别对所述左区域和右区域进行递归，以获取整个点集s的凸包，包 括：
[0032] 首先设Sl为所述直线PdrPul右边的点集，在点集Sl中找出距离所述直线PdrP ul最远 的点作为PdrPul右边的一个极端点Pb ;
[0033] 连接Pdr和Pb，以及Pb和Pul，把Pdr右侧的点集记为A，Pul右侧的点集记为B;
[0034] 对边PdrPb和点集A、边PbPui和点集B分别递归调用；
[0035] 递归结束后，依次连接所有的最远点，得到点集s 1的凸包，并将该凸包作为点集s 的右凸包；
[0036] 相应地，采用所述右凸包的相同计算方法计算出点集s的左凸包，以得到整个点集 s的凸包。
[0037]可选地，所述步骤S4中确定所述二维凸包内的最小点线距离，包括：
[0038] 梁用W下公忒i+笪m 一維Λ包O毎个顶占刹毎备仂的距离
[0039]
[0040] 其中，& =叉1+111、匕=5^+1-}^、〇 = 21+1-2：1，;[，]_ = 1，2,3'"111，;[为凸包(>)的边的序号，]· 为凸包Q的顶点的序号;对于不同取值的i都对应着一个最大点线距离Hiaxd ij，j = 1，2，3…m， 在这m个j = I，2，3···ηι中找出最小点线距离d〇(0) =minmaxd(0);最小点线距离对应的顶点 为P〃o,且所述最小点线距离对应的边为线段1〇。
[0041 ]另一方面，本发明还提供了一种平面度误差评定装置，该装置还包括：
[0042] 初始化单元，用于初始化旋转角和当前的最佳平面度误差；
[0043] 投影单元，用于将零件平面的被测点集投影至所述旋转角下的投影面；
[0044] 凸包确定单元，用于确定所述投影面内所述投影点的二维凸包；
[0045]最小点线距离确定单元，用于确定所述二维凸包内的最小点线距离；
[0046] 最小区域条件判断单元，用于根据所述最小点线距离判断是否满足最小区域条 件；
[0047] 平面度误差更新单元，用于在满足最小区域条件时，更新所述当前的最佳平面度 误差；
[0048] 收敛条件判断单元，用于根据更新的所述当前的最佳平面度误差判断是否满足收 敛条件，以确定最终的平面度误差评定结果。
[0049] 可选地，所述平面度误差更新单元，还用于：
[0050] 在不满足最小区域条件时，更改所述旋转角，以通过所述投影单元、所述凸包确定 单元、所述最小点线距离确定单元、所述最小区域条件判断单元、所述平面度误差更新单元 以及所述收敛条件判断单元重复执行平面度误差的评定。
[0051] 本发明的平面度误差评定方法及装置，以旋转投影面的旋转角度作为变量，将平 面度误差求解问题转化成了闭区间上的一元函数寻优问题，因此计算量小，并可以结合现 代优化算法高效求解该一元函数寻优问题。并且，由于采用凸包快速生成方法生成投影面 上投影点集的二维凸包，因此易于实施。同时，通过改变投影面的旋转角度得到不同的二维 凸包以从中找到符合最小区域条件的解，执行步骤简单且易于编程实现，因此该方法容易 在实际工程中应用。
【附图说明】
[0052]为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明 的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据 这些附图获得其他的附图。
[0053] 图1为本发明一个实施例的平面度误差评定方法的流程图；
[0054] 图2为本发明一个实施例的平面度误差评定装置的结构示意图；
[0055] 图3为本发明一个实施例的平面度误差评定方法的效果示意图。
【具体实施方式】
[0056] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚地描述，显然，所描述的实施例是本发明 一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有 做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0057] 图1为本发明一个实施例的平面度误差评定方法的流程图；如图1所示，该方法包 括：
[0058] SI:初始化旋转角和当前的最佳平面度误差；
[0059] S2:将零件平面的被测点集投影至所述旋转角下的投影面；
[0060] S3:确定所述投影面内所述投影点的二维凸包；
[0061] S4:确定所述二维凸包内的最小点线距离；
[0062] S5:根据所述最小点线距离判断是否满足最小区域条件；
[0063] S6:在满足最小区域条件时，更新所述当前的最佳平面度误差；
[0064] S7:根据更新的所述当前的最佳平面度误差判断是否满足收敛条件，以确定最终 的平面度误差评定结果。
[0065] 本实施例的平面度误差评定方法，以旋转投影面的旋转角度作为变量，将平面度 误差求解问题转化成了闭区间上的一元函数寻优问题，因此计算量小，并可以结合现代优 化算法高效求解该一元函数寻优问题。并且，由于采用凸包快速生成方法生成投影面上投 影点集的二维凸包，因此易于实施。同时，通过改变投影面的旋转角度得到不同的二维凸包 以从中找到符合最小区域条件的解，执行步骤简单且易于编程实现，因此该方法容易在实 际工程中应用。
[0066] 作为另一种优选的实施方式，所述方法还包括：
[0067] 在不满足最小区域条件时，更改所述旋转角，并重复执行步骤S2至S7。
[0068] 作为另一种优选的实施方式，所述步骤Sl中初始化旋转角和当前的最佳平面度误 差，包括：
[0069] 在测量空间直角坐标系O-XYZ中测量并获取零件平面的被测点点集P，将所述点集 P内第i个点的坐标表示为P1U, yi，Zl)，并过z轴建立旋转平面S，所述平面S绕z轴以Θ作为 旋转角进行旋转，将θ〇作为旋转角的迭代初值，并令当前的最佳平面度误差fbd为0;
[0070] 其中，所述坐标系O-XYZ包括X轴、y轴和Z轴，i = l，2,…，η，η为被测点的数目，且η 为大于等于3的正整数，所述平面S与平面XOZ的夹角为Θ，且Θ的取值区间为[0，ji]。
[0071] 作为另一种优选的实施方式，所述S2中将零件平面的被测点集投影至所述旋转角 下的投影面，包括：
[0072] 将O-XYZ坐标系绕z轴逆时针旋转Θ角，以构建新坐标系O-X7 V Z7，
[0073]其中，所述新坐标系Ol Y' Z7包括X7轴、y'轴和乂轴，所述坐标系Ol疒Z7的乂轴 与V轴位于所述平面S内，平面V Of与所述平面S重合；
[0074] 在所述新坐标系O-XY Z7中，将所述被测点坐标表示为P' i，/ ，并将新 坐标系O-XH^下的被测点Iy i(x\，/i，z\)与原始坐标系O-XYZ下的被测点Pi(XiJiJi) 的转换关系表示为下式：
[0075] (X7I1Y7 i,z/i,l) = (xi,yi,zi,l) · Rz
[0076] 以所述平面S作为投影面，将新坐标系0-f疒Z'下的被测点点集K i投影至旋转投 影面S，得到平面X7OZ'内的投影点集P〃，并将所述点集P〃内第i个点的坐标表示为PWY Z%);且将所述点集P〃内第i个点的坐标可以表示为PW ;
[0077] 其中，其中Rz为坐标系O-XYZ绕z轴逆时针旋转后各个被测点坐标的变换矩阵，且Rz 的表达式为：
[0078:
[0079] X7I = Xi · cosB-yi · SinBjZ7I = Zio
[0080] 作为另一种优选的实施方式，所述S3中确定所述投影面内所述投影点的二维凸 包，包括：
[0081] 采用凸包快速生成法生成所述平面X7OZ'内的投影点集P〃的一个二维凸包Q，所述 二维凸包Q的边数和顶点数均为Il 1为二维凸包Q的一条边，p〃4Pp〃1+1为线段^的两个端 点，其中? = 1，2,3···πι。
[0082] 作为另一种优选的实施方式，所述步骤3中采用凸包快速生成法生成所述平面X' OZ7内的投影点集P〃的一个二维凸包Q，包括：
[0083] 在平面X7 OZ7内的投影点集P〃中取两个极端点Pdr和Pul，其中，Pdr为最右最下的极 端点，Piil为最左最上的极端点；
[0084] 采用经过Pdr和Pul两点的直线PdrPul将整个点集P〃分为左区域和右区域，并分别对 所述左区域和右区域进行递归，以获取整个点集s的凸包。
[0085] 作为另一种优选的实施方式，所述分别对所述左区域和右区域进行递归，以获取 整个点集s的凸包，包括：
[0086] 首先设Sl为所述直线PdrPul右边的点集，在点集Sl中找出距离所述直线PdrPul最远 的点作为PdrPul右边的一个极端点Pb ;
[0087] 连接Pdr和Pb，以及Pb和Pul，把Pdr右侧的点集记为A，Pul右侧的点集记为B ;
[0088] 对边PdrPb和点集A、边PbPui和点集B分别递归调用；
[0089] 递归结束后，依次连接所有的最远点，得到点集si的凸包，并将该凸包作为点集s 的右凸包；
[0090] 相应地，采用所述右凸包的相同计算方法计算出点集s的左凸包，以得到整个点集 s的凸包。
[0091] 作为另一种优选的实施方式，所述步骤S4中确定所述二维凸包内的最小点线距 离，包括：
[0092]采用以下公式计算出二维凸包Q每个顶点到每条边的距离d1J:
[0093]
[0094] 其中，& =叉1+111、匕=5^+1-}^、〇 = 21+1-2：1，;[，]_ = 1，2,3."111，;[为凸包(>)的边的序号，]· 为凸包Q的顶点的序号;对于不同取值的i都对应着一个最大点线距离Hiaxd ij，j = I，2，3…m， 在这m个j = I，2，3···ηι中找出最小点线距离d〇(0) =minmaxd(0);最小点线距离对应的顶点 为P〃o,且所述最小点线距离对应的边为线段1〇。
[0095] 图2为本发明一个实施例的平面度误差评定装置的结构示意图；如图2所示，该装 置还包括：
[0096] 初始化单元10,用于初始化旋转角和当前的最佳平面度误差；
[0097] 投影单元20,用于将零件平面的被测点集投影至所述旋转角下的投影面；
[0098]凸包确定单元30,用于确定所述投影面内所述投影点的二维凸包；
[0099]最小点线距离确定单元40,用于确定所述二维凸包内的最小点线距离；
[0100]最小区域条件判断单元50,用于根据所述最小点线距离判断是否满足最小区域条 件；
[0101] 平面度误差更新单元60,用于在满足最小区域条件时，更新所述当前的最佳平面 度误差；
[0102] 收敛条件判断单元70,用于根据更新的所述当前的最佳平面度误差判断是否满足 收敛条件，以确定最终的平面度误差评定结果。
[0103] 本实施例所述的装置可以用于执行上述方法实施例，其原理和技术效果类似，此 处不再赘述。
[0104] 作为另一种优选的实施方式，所述平面度误差更新单元，还用于：
[0105] 在不满足最小区域条件时，更改所述旋转角，以通过所述投影单元、所述凸包确定 单元、所述最小点线距离确定单元、所述最小区域条件判断单元、所述平面度误差更新单元 以及所述收敛条件判断单元重复执行平面度误差的评定。
[0106] 本实施例所述的装置可以用于执行上述方法实施例，其原理和技术效果类似，此 处不再赘述。
[0107]需要说明的是，对于装置实施例而言，由于其与方法实施例基本相似，所以描述的 比较简单，相关之处参见方法实施例的部分说明即可。
[0108] 下面以一具体的实施例来说明本发明，但不限定本发明的保护范围。
[0109] 首先，在一个待测机械零件平面上采集20个测试点，测试点的坐标值如表1所示：
[0110] 表1机械零件平面的20个测试点
[0113] 基于此，本实施例的平面度误差评定方法的具体步骤如下：
[0114] 步骤1:在测量空间直角坐标系O-XYZ中测量并获取零件平面的被测点点集P，将所 述点集P内第i个点的坐标表示为P1U 1, yi，Zl)，并过Z轴建立旋转平面S，所述平面S绕z轴以 Θ作为旋转角进行旋转，将θ〇作为旋转角的迭代初值，并令当前的最佳平面度误差fbe3StS〇;
[0115] 其中，所述坐标系O-XYZ包括X轴、y轴和z轴，i = l，2,…，n，n为被测点的数目，且η 为大于等于3的正整数，所述平面S与平面XOZ的夹角为Θ，且Θ的取值区间为[0，ji];
[0116] 步骤2:将O-XYZ坐标系绕z轴逆时针旋转Θ角，以构建新坐标系O-X7 VZ7，
[0117] 其中，所述新坐标系O-X7 Y' Z7包括X7轴、y'轴和y轴，所述坐标系O-X7 V Z7的y轴 与V轴位于所述平面S内，平面V Of与所述平面S重合；
[0118] 在所述新坐标系0-X'V Z7中，将所述被测点坐标表示为P' i，/ 〇，并将新 坐标系O-XH^下的被测点Iy i(x\，/i，z\)与原始坐标系O-XYZ下的被测点Pi(XiJiJi) 的转换关系表示为下式：
[0119] (X7ilY7i,Z7I jD = (XijYijZijI) · Rz
[0120] 以所述平面S作为投影面，将新坐标系0-f疒Z'下的被测点点集K i投影至旋转投 影面S，得到平面X7OZ'内的投影点集P〃，并将所述点集P〃内第i个点的坐标表示为PW、 Z%);且将所述点集P〃内第i个点的坐标可以表示为PW i，" 〇 ;
[0121] 其中，其中Rz为坐标系O-XYZ绕z轴逆时针旋转后各个被测点坐标的变换矩阵， =Xi · cosS-yi · SinSjZ7I = Zi;
[0122] 步骤3:采用凸包快速生成法生成所述平面X7 OZ7内的投影点集P〃的一个二维凸包 Q，所述二维凸包Q的边数和顶点数均为m，h为二维凸包Q的一条边，p〃4Pp〃1+1为线段I1的两 个端点，其中? = 1，2,3···πι;
[0123] 步骤4:采用以下公式计算出二维凸包Q每个顶点到每条边的距离d1J:
[0124]
[0125] 兵 T，a = Xi+i-Xi、b = yi+i-yi、c = Zi+i-Zi，i，j = _L，Z，；i."m，i73 ? "Eiy tfJ 姐 tfJ斤，'，J 为凸包Q的顶点的序号;对于不同取值的i都对应着一个最大点线距离Hiaxdij，j = I，2，3…m， 在这m个j = I，2，3···ηι中找出最小点线距离d〇(0) =minmaxd(0);最小点线距离对应的顶点 为P〃o,且所述最小点线距离对应的边为线段1〇;
[0126] 步骤5:判断点集P"中是否存在满足最小区域条件的点：
[0127] 若有且仅有顶点P〃〇到Io的距离为do,且线段Io的端点上存在重合投影点或除端点 位置外还有其他的共线点，则求出线段Io上的任意三个点在点集P中的对应点A、B、C以及顶 点Ρ〃 ο在点集P中的对应点D，跳转到步骤6;
[0128] 若除顶点P〃〇外还有其他顶点到Io的距离为do,则求出线段Io上任意两个点在点集 P中的对应点E、F以及任意两个到Io的距离为do的顶点在点集P中的对应点G、H，跳转到步骤 7；
[0129] 若有且仅有顶点P〃〇到Io的距离为do,且线段Io上有且仅有两个点，则更改旋转角， 跳转到步骤2;
[0130]步骤6:将D点投影至A、B、C三点所构成的平面内，得到投影点D'，算出线段AV与线 段AB的夹角，算出线段AD'与线段AC的夹角β2;算出线段BD'与线段AB的夹角β3;算出线段 BD7与线段BC的夹角说;算出线段CD'与线段BC的夹角算出线段CD'与线段AC的夹角β6;
[0131] = 贝IJA、B、C、D四点满足最小区域条件的三角形准则，计算点 D到A、B、C三点所构成的平面的距离h，跳转至步骤8;
[0132] 若βι+&+β3+β4+β5+β6矣:π，则更改旋转角，跳转到步骤2;
[0133] 步骤7:计算出一个与线段EF以及线段GH平行的平面，并将E、F、G、H四点投影至该 平面得到线段f K和线段C；算出线段f K与线段C把的夹角β〇;
[0134] 若β〇矣0,则E、F、G、H四点满足最小区域条件的交叉准则，计算线段EF与线段GH的 最短距离h，跳转至步骤8;
[0135] 若β〇 = 0，则更改旋转角，跳转到步骤2;
[0136] 步骤8:记录当前的迭代次数；
[0137] 若迭代次数等于1，则将fbe3St的值变更为h的值，跳转到步骤9;
[0138] 若迭代次数大于1，且h小于fbe3St，则变更fbe3St的值为h的值，跳转到步骤9;
[0139] 若迭代次数大于1，且h大于或等于fbe3St，则跳转到步骤9;
[0140] 步骤9:判断当前迭代次数是否达到最大迭代次数；
[0141]若当前迭代次数达到最大迭代次数，则停止迭代，并输出平面度误差以及对应的 满足最小区域条件的高极点和低极点；
[0142]若当前迭代次数未达到最大迭代次数，则更改旋转角，并且重复步骤2至步骤8直 到迭代终止。
[0143] 其中，步骤5至步骤7的过程为判断最小点线距离对应的测试点是否满足最小区域 条件的过程。
[0144] k -袖從冰的实施方式，上述步骤2中的Rz的表达式如下：
[0145]
[0146] 作为另一种优选的实施方式，步骤3中采用凸包快速生成法生成所述平面fOZ'内 的投影点集P〃的一个二维凸包Q，包括：
[0147] 在平面X7 OZ7内的投影点集P〃中取两个极端点Pdr和Pul，其中，Pdr为最右最下的极 端点，Pul为最左最上的极端点；
[0148] 采用经过Pdr和Pul两点的直线PdrPul将整个点集P〃分为左区域和右区域，并分别对 所述左区域和右区域进行递归，以获取整个点集s的凸包。
[0149] 作为另一种优选的实施方式，所述分别对所述左区域和右区域进行递归，以获取 整个点集s的凸包，包括：
[0150] 首先设Sl为所述直线PdrPul右边的点集，在点集Sl中找出距离所述直线PdrP ul最远 的点作为PdrPul右边的一个极端点Pb ;
[0151] 连接Pdr和Pb，以及Pb和Pui，把Pdr右侧的点集记为A，PU1右侧的点集记为B;
[0152] 对边PdrPb和点集A、边PbPui和点集B分别递归调用；
[0153] 递归结束后，依次连接所有的最远点，得到点集si的凸包，并将该凸包作为点集s 的右凸包；
[0154] 相应地，采用所述右凸包的相同计算方法计算出点集s的左凸包，以得到整个点集 s的凸包。
[0155] 在本实施例中更改旋转角的过程结合了现代优化算法中的遗传算法，父代数目为 100,迭代次数为200。其中，旋转角的"更改"过程即为遗传算法中父代交叉、变异的过程;更 改后的旋转角即为遗传算法中父代交叉、变异所产生的子代;遗传算法中的子代适应度为 不同旋转角对应的最小点线距离的相反数。本实施例的计算结果满足交叉准则，平面度误 差为0.041833mm，优于凸壳法的结果（0.04185mm)和最小二乘法的结果（0.04381mm)，并且 运算时间不超过2s。图3为本发明一个实施例的平面度误差评定方法的效果示意图，如图3 所示，☆代表高极点（5号测试点和17号测试点），?代表低极点（16号测试点和20号测试 点），所得到的最小区域面的法向量为(0.01718091,-0.11059279,3.90528946)。
[0156] 以上实施例仅用于说明本发明的技术方案，而非对其限制；尽管参照前述实施例 对本发明进行了详细的说明，本领域的普通技术人员应当理解:其依然可以对前述各实施 例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换;而这些修改或替 换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
【主权项】
1. 一种平面度误差评定方法，其特征在于，包括： S1:初始化旋转角和当前的最佳平面度误差； S2:将零件平面的被测点集投影至所述旋转角下的投影面； S3:确定所述投影面内所述投影点的二维凸包； S4:确定所述二维凸包内的最小点线距离； S5:根据所述最小点线距离判断是否满足最小区域条件； S6:在满足最小区域条件时，更新所述当前的最佳平面度误差； S7:根据更新的所述当前的最佳平面度误差判断是否满足收敛条件，W确定最终的平 面度误差评定结果。2. 根据权利要求1所述的方法，其特征在于，所述方法还包括： 在不满足最小区域条件时，更改所述旋转角，并重复执行步骤S2至S7。3. 根据权利要求1所述的方法，其特征在于，所述步骤S1中初始化旋转角和当前的最佳 平面度误差，包括： 在测量空间直角坐标系0-XYZ中测量并获取零件平面的被测点点集P，将所述点集P内 第i个点的坐标表示为Pl(Xl，yl，Zl)，并过z轴建立旋转平面S，所述平面S绕z轴Wθ作为旋转 角进行旋转，将θ〇作为旋转角的迭代初值，并令当前的最佳平面度误差fbest为0; 其中，所述坐标系0-XYZ包括X轴、y轴和Z轴，i = l，2,…，n，n为被测点的数目，且η为大 于等于3的正整数，所述平面S与平面ΧΟΖ的夹角为Θ，且Θ的取值区间为[〇，ji]。4. 根据权利要求3所述的方法，其特征在于，所述S2中将零件平面的被测点集投影至所 述旋转角下的投影面，包括： 将0-XYZ坐标系绕Z轴逆时针旋转Θ角，W构建新坐标系0-Χ'Υ'Ζ'， 其中，所述新坐标系0-Χ'Υ'Ζ'包括χ/轴、/轴和轴，所述坐标系0-Χ'Υ'Ζ'的Ζ'轴与X' 轴位于所述平面S内，平面X ' ΟΖ '与所述平面S重合； 在所述新坐标系0-Χ ' Υ ' Ζ '中，将所述被测点坐标表示为Pi '（X ' 1，y ' 1，Ζ ' 1)，并将新坐标 系0-Χ'Υ'Ζ'下的被测点Pl'（x'l，y'l，z'l)与原始坐标系0-XYZ下的被测点Pl(Xl，yl，Zl)的转 换关系表示为下式： (x'i,y'i,z'i,l) = (xi,yi,zi,l) · Rz W所述平面S作为投影面，将新坐标系0-Χ'Υ'Ζ'下的被测点点集Pi'投影至旋转投影面 S，得到平面X'OZ'内的投影点集P"，并将所述点集P"内第i个点的坐标表示为Pi"(xi"，zi"); 且将所述点集P"内第i个点的坐标可W表示为Pi"(xi'，zi'）； 其中，其中Rz为坐标系0-XYZ绕Z轴逆时针旋转后各个被测点坐标的变换矩阵，且Rz的表 达式为：x'i = xi · cos白一yi · sin白，z'i = zi。5. 根据权利要求4所述的方法，其特征在于，所述S3中确定所述投影面内所述投影点的 二维凸包，包括： 采用凸包快速生成法生成所述平面X'oz'内的投影点集P"的一个二维凸包Q，所述二维 凸包Q的边数和顶点数均为m，li为二维凸包Q的一条边，p"i和P" 1+1为线段h的两个端点，其 中1 = 1,2,3...m。6. 根据权利要求5所述的方法，其特征在于，所述步骤3中采用凸包快速生成法生成所 述平面X'OZ'内的投影点集P"的一个二维凸包Q，包括： 在平面X'OZ'内的投影点集P"中取两个极端点P化和Pul,其中，Pdr为最右最下的极端点， Pul为最左最上的极端点； 采用经过Pdr和Pul两点的直线PdrPul将整个点集P"分为左区域和右区域，并分别对所述 左区域和右区域进行递归，W获取整个点集S的凸包。7. 根据权利要求6所述的方法，其特征在于，所述分别对所述左区域和右区域进行递 归，W获取整个点集S的凸包，包括： 首先设si为所述直线PdrPul右边的点集，在点集si中找出距离所述直线PdrPul最远的点 作为PdrPul右边的一个极端点Pb ; 连接P化和Pb，W及Pb和Pul，把P化右侧的点集记为A，Pul右侧的点集记为B; 对边PdrPb和点集A、边PbPul和点集B分别递归调用； 递归结束后，依次连接所有的最远点，得到点集si的凸包，并将该凸包作为点集S的右 凸包； 相应地，采用所述右凸包的相同计算方法计算出点集S的左凸包，W得到整个点集S的 凸包。8. 根据权利要求5所述的方法，其特征在于，所述步骤S4中确定所述二维凸包内的最小 点线距离，包括： 采用W下公式计算出二维凸包Q每个顶点到每条边的距离dij:其中，a = xi+广Xi、b = yi+广yi、c = zi+广Zi，i , j = l ,2,3.''111，;[为凸包〇的边的序号〇为凸 包Q的顶点的序号;对于不同取值的i都对应着一个最大点线距离maxdij，j = 1，2，3· ·'m，在运 m个j = 1，2，3· · ·πι中找出最小点线距离do(目）=minmaxd(目）；最小点线距离对应的顶点为Po"， 且所述最小点线距离对应的边为线段1〇。9. 一种平面度误差评定装置，其特征在于，包括： 初始化单元，用于初始化旋转角和当前的最佳平面度误差； 投影单元，用于将零件平面的被测点集投影至所述旋转角下的投影面； 凸包确定单元，用于确定所述投影面内所述投影点的二维凸包； 最小点线距离确定单元，用于确定所述二维凸包内的最小点线距离； 最小区域条件判断单元，用于根据所述最小点线距离判断是否满足最小区域条件； 平面度误差更新单元，用于在满足最小区域条件时，更新所述当前的最佳平面度误差； 收敛条件判断单元，用于根据更新的所述当前的最佳平面度误差判断是否满足收敛条 件，W确定最终的平面度误差评定结果。10.根据权利要求9所述的装置，其特征在于，所述平面度误差更新单元，还用于： 在不满足最小区域条件时，更改所述旋转角，W通过所述投影单元、所述凸包确定单 元、所述最小点线距离确定单元、所述最小区域条件判断单元、所述平面度误差更新单元W 及所述收敛条件判断单元重复执行平面度误差的评定。
【文档编号】G01B11/30GK105841640SQ201610282605
【公开日】2016年8月10日
【申请日】2016年4月29日
【发明人】陶扬, 赵军鹏, 王春洁
【申请人】北京航空航天大学