CN103292773A - 一种基于最小区域的对称度误差评定方法 - Google Patents

一种基于最小区域的对称度误差评定方法 Download PDF

Info

Publication number
CN103292773A
CN103292773A CN201310241194XA CN201310241194A CN103292773A CN 103292773 A CN103292773 A CN 103292773A CN 201310241194X A CN201310241194X A CN 201310241194XA CN 201310241194 A CN201310241194 A CN 201310241194A CN 103292773 A CN103292773 A CN 103292773A
Authority
CN
China
Prior art keywords
point
measuring
points
area
region
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201310241194XA
Other languages
English (en)
Other versions
CN103292773B (zh
Inventor
钟艳茹
陈磊磊
孙永厚
黄美发
范凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guilin University of Electronic Technology
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201310241194.XA priority Critical patent/CN103292773B/zh
Publication of CN103292773A publication Critical patent/CN103292773A/zh
Application granted granted Critical
Publication of CN103292773B publication Critical patent/CN103292773B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

一种基于最小区域的对称度误差评定方法，该方法主要评定以线为基准的对称度误差。首先测量并获取被测特征以及基准线上测点坐标。通过最小区域法拟合基准直线，然后通过坐标系变换，使基准线对应理想特征的方向矢量与z轴平行，将被测特征投影与xoy平面内，将评定区域划分为若干区域，通过分别在这些区域中计算对应的误差值的方法，计算被测特征参数的最优值以及对应的对称度误差值。本发明可准确计算出满足最小区域的对称度误差以及被测特征参数的最优值。

Description

一种基于最小区域的对称度误差评定方法
技术领域
本发明涉及一种基于最小区域的对称度误差评定方法，属于精密计量与计算机应用领域，可用于各种情况下几何产品中以线为基准对称度的合格性检测，并为加工过程与加工工艺的改进提供指导。
背景技术
对称度是机械几何产品中最常见的检测项目之一，也是影响产品的质量、性能以及装配的重要因素之一。根据国际标准ISO/1101的规定，零件误差计算以最小区域法为仲裁方法。在相同的硬件检测条件下，应用最小区域计算零件误差，可提高检测设备的检测精度。
最小区域法，属于不可微复杂最优化问题，目前，国内外学者主要采用传统优化方法、智能算法、计算几何方法等。这些方法由于存在计算稳定性差、计算效率低、对采点数量有限制、计算结果难以达到精确解等缺陷，导致最小区域法很难在实际检测中应用。目前市场上一般都采用成熟的最小二乘法近似地计算对称度误差。
发明内容
为了克服上述技术缺点，本发明的目的是提供一种基于最小区域的零件对称度误差评定方法。本方法不仅提高了测量仪器检测零件对称度的精度，而且算法稳定性好、计算效率高，可以推广应用于其几何体作用尺寸的计算。
为解决上述问题，本发明一种基于最小区域的零件对称度误差评定方法，主要包括以下步骤：
步骤1：将待检测的零件特征置于测量平台上，在测量空间直角坐标系中测量并提取被测特征和基准线上的测点；应用最小区域的方法拟合基准线，获得基准线对应理想直线的方向矢量                                                
Figure 142819DEST_PATH_IMAGE001
，以及理想直线上一点
Figure 201310241194X100002DEST_PATH_IMAGE002
。
步骤2：进行坐标旋转变换，使
Figure 591118DEST_PATH_IMAGE001
与坐标系z轴平行；坐标旋转变换后，设被测特征的测点在xoy坐标平面上的投影为
Figure 370855DEST_PATH_IMAGE003
，
Figure 201310241194X100002DEST_PATH_IMAGE004
=1, 2,…n, n为测点数目且n为正整数，所有测点
Figure 172589DEST_PATH_IMAGE003
形成测点集
Figure 986961DEST_PATH_IMAGE005
；在集合
Figure 809423DEST_PATH_IMAGE005
中，剔除在xoy平面内重合的测点；坐标旋转变换后，计算
Figure 607615DEST_PATH_IMAGE002
在xoy平面内的投影
Figure 201310241194X100002DEST_PATH_IMAGE006
。
步骤3：任意取2个测点，计算2点连线的方向矢量，将其作为包容区域中心线的方向矢量
Figure 210110DEST_PATH_IMAGE007
，设=
Figure 201310241194X100002DEST_PATH_IMAGE008
；与
Figure 403511DEST_PATH_IMAGE007
垂直的矢量
Figure 157841DEST_PATH_IMAGE009
为。
步骤4：以过点且平行于的直线为界，将xoy平面区域划分为2个区域左半区域C 1、右半区域C 2，通过下式，判断所有测点的归属区域，
Figure 785765DEST_PATH_IMAGE013
当
Figure 201310241194X100002DEST_PATH_IMAGE014
大于0时，
Figure 433915DEST_PATH_IMAGE015
位于C 1区，反之，位于C 2区；其中
Figure 201310241194X100002DEST_PATH_IMAGE016
为辅助点，
Figure 883668DEST_PATH_IMAGE017
等于
Figure 201310241194X100002DEST_PATH_IMAGE018
加
Figure 625359DEST_PATH_IMAGE007
的横坐标，等于
Figure 201310241194X100002DEST_PATH_IMAGE020
加
Figure 964254DEST_PATH_IMAGE007
的纵坐标。
步骤5：计算C 1区内各测点到过O点且方向矢量为
Figure 462231DEST_PATH_IMAGE007
直线(即包容区域中心线)的距离，并记录C 1区内各测点到包容区域中心线的最大距离、最小距离对应的测点
Figure 640403DEST_PATH_IMAGE021
、
Figure 201310241194X100002DEST_PATH_IMAGE022
，则重新计算包容区域垂直方向
Figure 591041DEST_PATH_IMAGE009
；
Figure 4705DEST_PATH_IMAGE023
判断
Figure 357189DEST_PATH_IMAGE009
与计算前的
Figure 706262DEST_PATH_IMAGE009
是否平行；如果平行，则计算C 1区内各点相对包容区域中心线的对称度
Figure 201310241194X100002DEST_PATH_IMAGE024
，并跳转到步骤6；如果不平行，跳转到步骤4。
步骤6：根据
Figure 144196DEST_PATH_IMAGE009
计算整个评定区域的对称度误差
Figure 95972DEST_PATH_IMAGE025
；如果
Figure 568541DEST_PATH_IMAGE024
等于
Figure 681991DEST_PATH_IMAGE025
，则跳转到步骤15；如果
Figure 13746DEST_PATH_IMAGE024
不等于，则跳转到步骤7。
步骤7：按照步骤4的方法，根据
Figure 830710DEST_PATH_IMAGE009
判断测点所在的区域。
步骤8：计算C 2区内各测点到包容区域中心线的距离，并记录C 2区内各测点到包容区域中心线的最大距离、最小距离对应的测点
Figure 201310241194X100002DEST_PATH_IMAGE026
、
Figure 52744DEST_PATH_IMAGE027
，设两点的坐标分别为
Figure 201310241194X100002DEST_PATH_IMAGE028
、
Figure 465270DEST_PATH_IMAGE029
，则应用下式计算包容区域垂直方向
Figure 896864DEST_PATH_IMAGE009
；
Figure 201310241194X100002DEST_PATH_IMAGE030
判断
Figure 78447DEST_PATH_IMAGE009
与计算前的
Figure 799278DEST_PATH_IMAGE009
是否平行；如果
Figure 699101DEST_PATH_IMAGE009
与计算前的
Figure 671736DEST_PATH_IMAGE009
平行，则计算C 2区内各点相对中心线的对称度
Figure 973405DEST_PATH_IMAGE031
，并跳转到步骤9；如果
Figure 333979DEST_PATH_IMAGE009
与计算前的
Figure 252256DEST_PATH_IMAGE009
不平行，则跳转到步骤7。
步骤9：根据
Figure 622058DEST_PATH_IMAGE009
计算整个评定区域的对称度误差
Figure 512653DEST_PATH_IMAGE025
；如果
Figure 716233DEST_PATH_IMAGE031
等于
Figure 325069DEST_PATH_IMAGE025
，则跳转到步骤15；如果
Figure 29719DEST_PATH_IMAGE031
不等于
Figure 40401DEST_PATH_IMAGE025
，则跳转到步骤10。
步骤10：计算各测点到包容区域中心线的距离，并记录各测点到包容区域中心线的最大值，最大值对应的测点即为与包容区域接触的测点；
如果接触点不止1个且接触点分布在O点的两侧，则计算各点相对中心线的对称度
Figure 201310241194X100002DEST_PATH_IMAGE032
，并跳转到步骤15；如果分布在O点一侧或只有1个接触点，则在接触点中，查询距离O点最近的接触点
Figure 414881DEST_PATH_IMAGE033
Figure 201310241194X100002DEST_PATH_IMAGE034
，以
Figure 511013DEST_PATH_IMAGE035
为界将xoy平面区域划分为C 3区、C 4区，通过下式可以判断测点的归属区域，
Figure 201310241194X100002DEST_PATH_IMAGE036
其中，
Figure 753776DEST_PATH_IMAGE037
为参考点，等于减
Figure 164346DEST_PATH_IMAGE009
的横坐标，
Figure 278932DEST_PATH_IMAGE039
等于
Figure 528648DEST_PATH_IMAGE020
减
Figure 513921DEST_PATH_IMAGE009
的纵坐标，
Figure 201310241194X100002DEST_PATH_IMAGE040
，
Figure 230205DEST_PATH_IMAGE041
；当
Figure DEST_PATH_IMAGE042
时，表示
Figure 832087DEST_PATH_IMAGE015
位于C 3区，当
Figure 885494DEST_PATH_IMAGE043
时，表示
Figure 131799DEST_PATH_IMAGE015
位于C 4区。
步骤11：分别连接有效接触点
Figure 612459DEST_PATH_IMAGE033
与C 3区内各个测点
Figure DEST_PATH_IMAGE044
，求各个连线与当前包容区域中心线的夹角集合，夹角的最小值即为旋转角度
Figure 436058DEST_PATH_IMAGE045
，
Figure 201310241194X100002DEST_PATH_IMAGE046
并记录C 3区中对应
Figure 699680DEST_PATH_IMAGE045
的测点
Figure 393967DEST_PATH_IMAGE047
，将该与点
Figure 622003DEST_PATH_IMAGE033
坐标代入下式，
Figure 201310241194X100002DEST_PATH_IMAGE048
得到迭代后的包容区域垂直方向矢量
Figure 420807DEST_PATH_IMAGE009
，并判断
Figure 969600DEST_PATH_IMAGE009
迭代前后的夹角是否大于90度，如果大于90度，则
Figure 57642DEST_PATH_IMAGE009
取反向，跳转到步骤10。
步骤12：根据计算整个评定区域的对称度误差
Figure 54734DEST_PATH_IMAGE025
；如果
Figure 864558DEST_PATH_IMAGE032
等于
Figure 123501DEST_PATH_IMAGE025
，则跳转到步骤15；如果不等于
Figure 146001DEST_PATH_IMAGE025
，则跳转到步骤13。
步骤13：计算各测点
Figure 669386DEST_PATH_IMAGE049
到中心线的距离，并记录各测点到中心线的最小距离，各测点到中心线的最小距离对应的测点即为与包容区域接触的测点；
如果接触点不止1个且接触点分布在O点的两侧，则计算整个评定区域的对称度误差
Figure 99230DEST_PATH_IMAGE025
，并跳转到步骤15；如果接触点分布在一侧或只有1个接触点，则在接触点中，查询距离O点最近的接触点
Figure 481801DEST_PATH_IMAGE033
；按照步骤10中的方式判断所有测点的归属区域。
步骤14：按照步骤11的方法，计算
Figure 22504DEST_PATH_IMAGE009
的值；跳转到步骤13。
步骤15：输出对称度误差
Figure 197133DEST_PATH_IMAGE025
，以及被测特征的参数
Figure 532300DEST_PATH_IMAGE007
。
与现有技术相比，本发明依据对称度自身的特征，在对称度误差优化的方向上，对包容区域进行相应分区和旋转变动，最终计算出对称度误差最优值和被测特征的参数最优值。
本发明的有益效果在于：本方法查询到的接触点相对位置满足最小区域判别准则，计算得到的对称度误差为最优值，计算稳定性好、计算效率高。
对于本领域技术人员来说，根据和应用本发明公开的构思，能够容易地对本发明方案进行各种变形和改变，应当注意的是，所有这些变形和改变都应当属于本发明的范围。
附图说明
图1 为本发明的以直线为基准的对称度的计算流程图。
具体实施方式
以下结合附图对本发明的实施方式进行详细说明。本发明是通过如附图1所示方案实现的，具体步骤如下：
步骤1：将待检测的零件特征置于测量平台上，在测量空间直角坐标系中测量并提取被测特征和基准线上的测点；应用最小区域的方法拟合基准线，获得基准线对应理想直线的方向矢量
Figure 261221DEST_PATH_IMAGE001
，以及理想直线上一点
Figure 746560DEST_PATH_IMAGE002
。
步骤2：进行坐标旋转变换，使
Figure 244538DEST_PATH_IMAGE001
与坐标系z轴平行；坐标旋转变换后，设被测特征的测点在xoy坐标平面上的投影为，=1, 2,…n, n为测点数目且n为正整数，所有测点
Figure 380487DEST_PATH_IMAGE003
形成测点集；在集合
Figure 82044DEST_PATH_IMAGE005
中，剔除在xoy平面内重合的测点；坐标旋转变换后，计算
Figure 51137DEST_PATH_IMAGE002
在xoy平面内的投影
Figure 471754DEST_PATH_IMAGE006
。
步骤3：任意取2个测点，计算2点连线的方向矢量，将其作为包容区域中心线的方向矢量，设=；与垂直的矢量
Figure 206492DEST_PATH_IMAGE009
为
Figure 897367DEST_PATH_IMAGE010
。
步骤4：以过点
Figure 309894DEST_PATH_IMAGE011
且平行于
Figure 72314DEST_PATH_IMAGE009
的直线为界，将xoy平面区域划分为2个区域左半区域C 1、右半区域C 2，通过下式，判断所有测点的归属区域，
Figure 240307DEST_PATH_IMAGE013
当大于0时，
Figure 112765DEST_PATH_IMAGE015
位于C 1区，反之，
Figure 680012DEST_PATH_IMAGE015
位于C 2区；其中
Figure 306166DEST_PATH_IMAGE016
为辅助点，等于加
Figure 368996DEST_PATH_IMAGE007
的横坐标，
Figure 697209DEST_PATH_IMAGE019
等于
Figure 306045DEST_PATH_IMAGE020
加的纵坐标。
步骤5：计算C 1区内各测点到过O点且方向矢量为
Figure 896743DEST_PATH_IMAGE007
直线(即包容区域中心线)的距离，并记录C 1区内各测点到包容区域中心线的最大距离、最小距离对应的测点
Figure 864699DEST_PATH_IMAGE021
、
Figure 491989DEST_PATH_IMAGE022
，则重新计算包容区域垂直方向
Figure 938014DEST_PATH_IMAGE009
；
Figure 68781DEST_PATH_IMAGE023
判断
Figure 614163DEST_PATH_IMAGE009
与计算前的
Figure 197591DEST_PATH_IMAGE009
是否平行；如果平行，则计算C 1区内各点相对包容区域中心线的对称度
Figure 978466DEST_PATH_IMAGE024
，并跳转到步骤6；如果不平行，跳转到步骤4。
步骤6：根据
Figure 963739DEST_PATH_IMAGE009
计算整个评定区域的对称度误差
Figure 680022DEST_PATH_IMAGE025
；如果
Figure 750747DEST_PATH_IMAGE024
等于
Figure 804153DEST_PATH_IMAGE025
，则跳转到步骤15；如果不等于
Figure 655752DEST_PATH_IMAGE025
，则跳转到步骤7。
步骤7：按照步骤4的方法，根据判断测点所在的区域。
步骤8：计算C 2区内各测点到包容区域中心线的距离，并记录C 2区内各测点到包容区域中心线的最大距离、最小距离对应的测点
Figure 211815DEST_PATH_IMAGE026
、
Figure 906101DEST_PATH_IMAGE027
，设两点的坐标分别为
Figure 88821DEST_PATH_IMAGE028
、，则应用下式计算包容区域垂直方向
Figure 935871DEST_PATH_IMAGE009
；
Figure 484664DEST_PATH_IMAGE030
判断
Figure 572706DEST_PATH_IMAGE009
与计算前的是否平行；如果
Figure 304219DEST_PATH_IMAGE009
与计算前的
Figure 973098DEST_PATH_IMAGE009
平行，则计算C 2区内各点相对中心线的对称度
Figure 638565DEST_PATH_IMAGE031
，并跳转到步骤9；如果
Figure 127315DEST_PATH_IMAGE009
与计算前的不平行，则跳转到步骤7。
步骤9：根据
Figure 184450DEST_PATH_IMAGE009
计算整个评定区域的对称度误差
Figure 348715DEST_PATH_IMAGE025
；如果
Figure 996865DEST_PATH_IMAGE031
等于
Figure 537568DEST_PATH_IMAGE025
，则跳转到步骤15；如果
Figure 915460DEST_PATH_IMAGE031
不等于
Figure 47364DEST_PATH_IMAGE025
，则跳转到步骤10。
步骤10：计算各测点到包容区域中心线的距离，并记录各测点到包容区域中心线的最大值，最大值对应的测点即为与包容区域接触的测点；
如果接触点不止1个且接触点分布在O点的两侧，则计算各点相对中心线的对称度
Figure 510706DEST_PATH_IMAGE032
，并跳转到步骤15；如果分布在O点一侧或只有1个接触点，则在接触点中，查询距离O点最近的接触点
Figure 258695DEST_PATH_IMAGE033
，以
Figure 793898DEST_PATH_IMAGE035
为界将xoy平面区域划分为C 3区、C 4区，通过下式可以判断测点的归属区域，
Figure 10116DEST_PATH_IMAGE036
其中，
Figure 892622DEST_PATH_IMAGE037
为参考点，等于
Figure 328599DEST_PATH_IMAGE018
减
Figure 766534DEST_PATH_IMAGE009
的横坐标，
Figure 983888DEST_PATH_IMAGE039
等于
Figure 190879DEST_PATH_IMAGE020
减
Figure 710853DEST_PATH_IMAGE009
的纵坐标，
Figure 636084DEST_PATH_IMAGE040
，
Figure 860392DEST_PATH_IMAGE041
；当
Figure 718626DEST_PATH_IMAGE042
时，表示
Figure 2977DEST_PATH_IMAGE015
位于C 3区，当
Figure 149925DEST_PATH_IMAGE043
时，表示
Figure 584448DEST_PATH_IMAGE015
位于C 4区。
步骤11：分别连接有效接触点
Figure 766031DEST_PATH_IMAGE033
与C 3区内各个测点
Figure 752441DEST_PATH_IMAGE044
，求各个连线与当前包容区域中心线的夹角集合，夹角的最小值即为旋转角度
Figure 652264DEST_PATH_IMAGE045
，
Figure 218375DEST_PATH_IMAGE046
并记录C 3区中对应的测点
Figure 287142DEST_PATH_IMAGE047
，将该
Figure 939840DEST_PATH_IMAGE047
与点
Figure 309642DEST_PATH_IMAGE033
坐标代入下式，
Figure 465816DEST_PATH_IMAGE048
得到迭代后的包容区域垂直方向矢量，并判断
Figure 278232DEST_PATH_IMAGE009
迭代前后的夹角是否大于90度，如果大于90度，则
Figure 717303DEST_PATH_IMAGE009
取反向，跳转到步骤10。
步骤12：根据
Figure 727985DEST_PATH_IMAGE009
计算整个评定区域的对称度误差
Figure 695941DEST_PATH_IMAGE025
；如果等于
Figure 910201DEST_PATH_IMAGE025
，则跳转到步骤15；如果
Figure 40968DEST_PATH_IMAGE032
不等于
Figure 710984DEST_PATH_IMAGE025
，则跳转到步骤13。
步骤13：计算各测点
Figure 28833DEST_PATH_IMAGE049
到中心线的距离，并记录各测点到中心线的最小距离，各测点到中心线的最小距离对应的测点即为与包容区域接触的测点；
如果接触点不止1个且接触点分布在O点的两侧，则计算整个评定区域的对称度误差
Figure 950653DEST_PATH_IMAGE025
，并跳转到步骤15；如果接触点分布在一侧或只有1个接触点，则在接触点中，查询距离O点最近的接触点；按照步骤10中的方式判断所有测点的归属区域。
步骤14：按照步骤11的方法，计算
Figure 245685DEST_PATH_IMAGE009
的值；跳转到步骤13。
步骤15：输出对称度误差
Figure 581988DEST_PATH_IMAGE025
，以及被测特征的参数
Figure 635395DEST_PATH_IMAGE007
。

Claims (1)

1.一种基于最小区域的对称度误差评定方法，其特征在于，具体步骤如下：
步骤1：将待检测的零件特征置于测量平台上，在测量空间直角坐标系中测量并提取被测特征和基准线上的测点；应用最小区域的方法拟合基准线，获得基准线对应理想直线的方向矢量                                                ，以及理想直线上一点
Figure 253189DEST_PATH_IMAGE002
；
步骤2：进行坐标旋转变换，使
Figure 341230DEST_PATH_IMAGE001
与坐标系z轴平行；坐标旋转变换后，设被测特征的测点在xoy坐标平面上的投影为
Figure 201310241194X100001DEST_PATH_IMAGE003
，
Figure 201310241194X100001DEST_PATH_IMAGE005
=1, 2,…n, n为测点数目且n为正整数，所有测点
Figure 545947DEST_PATH_IMAGE003
形成测点集
Figure 167683DEST_PATH_IMAGE006
；在集合中，剔除在xoy平面内重合的测点；坐标旋转变换后，计算
Figure 95505DEST_PATH_IMAGE002
在xoy平面内的投影
Figure 201310241194X100001DEST_PATH_IMAGE007
；
步骤3：任意取2个测点，计算2点连线的方向矢量，将其作为包容区域中心线的方向矢量
Figure 521938DEST_PATH_IMAGE008
，设
Figure 524529DEST_PATH_IMAGE008
=
Figure 201310241194X100001DEST_PATH_IMAGE009
；与
Figure 234865DEST_PATH_IMAGE008
垂直的矢量
Figure 336813DEST_PATH_IMAGE010
为
Figure 201310241194X100001DEST_PATH_IMAGE011
；
步骤4：以过点
Figure 578439DEST_PATH_IMAGE012
且平行于
Figure 819276DEST_PATH_IMAGE010
的直线为界，将xoy平面区域划分为2个区域左半区域C 1、右半区域C 2，通过下式，判断所有测点
Figure 201310241194X100001DEST_PATH_IMAGE013
的归属区域，
Figure 197168DEST_PATH_IMAGE014
当
Figure 201310241194X100001DEST_PATH_IMAGE015
大于0时，位于C 1区，反之，
Figure 198939DEST_PATH_IMAGE016
位于C 2区；其中
Figure 201310241194X100001DEST_PATH_IMAGE017
为辅助点，
Figure 730283DEST_PATH_IMAGE018
等于
Figure 201310241194X100001DEST_PATH_IMAGE019
加
Figure 165944DEST_PATH_IMAGE008
的横坐标，
Figure 672011DEST_PATH_IMAGE020
等于加
Figure 576644DEST_PATH_IMAGE008
的纵坐标；
步骤5：计算C 1区内各测点到过O点且方向矢量为直线(即包容区域中心线)的距离，并记录C 1区内各测点到包容区域中心线的最大距离、最小距离对应的测点
Figure 749317DEST_PATH_IMAGE022
、
Figure 201310241194X100001DEST_PATH_IMAGE023
，则重新计算包容区域垂直方向；
Figure 51171DEST_PATH_IMAGE024
判断
Figure 737367DEST_PATH_IMAGE010
与计算前的
Figure 882041DEST_PATH_IMAGE010
是否平行；如果平行，则计算C 1区内各点相对包容区域中心线的对称度
Figure 201310241194X100001DEST_PATH_IMAGE025
，并跳转到步骤6；如果不平行，跳转到步骤4；
步骤6：根据
Figure 995490DEST_PATH_IMAGE010
计算整个评定区域的对称度误差
Figure 609137DEST_PATH_IMAGE026
；如果
Figure 833445DEST_PATH_IMAGE025
等于
Figure 98204DEST_PATH_IMAGE026
，则跳转到步骤15；如果
Figure 382555DEST_PATH_IMAGE025
不等于
Figure 529502DEST_PATH_IMAGE026
，则跳转到步骤7；
步骤7：按照步骤4的方法，根据
Figure 744452DEST_PATH_IMAGE010
判断测点所在的区域；
步骤8：计算C 2区内各测点到包容区域中心线的距离，并记录C 2区内各测点到包容区域中心线的最大距离、最小距离对应的测点、
Figure 863717DEST_PATH_IMAGE028
，设两点的坐标分别为
Figure 201310241194X100001DEST_PATH_IMAGE029
、，则应用下式计算包容区域垂直方向
Figure 907208DEST_PATH_IMAGE010
；
Figure 201310241194X100001DEST_PATH_IMAGE031
判断
Figure 473319DEST_PATH_IMAGE010
与计算前的
Figure 447091DEST_PATH_IMAGE010
是否平行；如果与计算前的
Figure 194784DEST_PATH_IMAGE010
平行，则计算C 2区内各点相对中心线的对称度，并跳转到步骤9；如果
Figure 907711DEST_PATH_IMAGE010
与计算前的
Figure 704766DEST_PATH_IMAGE010
不平行，则跳转到步骤7；
步骤9：根据
Figure 251285DEST_PATH_IMAGE010
计算整个评定区域的对称度误差
Figure 159198DEST_PATH_IMAGE026
；如果等于
Figure 826251DEST_PATH_IMAGE026
，则跳转到步骤15；如果
Figure 922383DEST_PATH_IMAGE032
不等于
Figure 633987DEST_PATH_IMAGE026
，则跳转到步骤10；
步骤10：计算各测点到包容区域中心线的距离，并记录各测点到包容区域中心线的最大值，最大值对应的测点即为与包容区域接触的测点；
如果接触点不止1个且接触点分布在O点的两侧，则计算各点相对中心线的对称度，并跳转到步骤15；如果分布在O点一侧或只有1个接触点，则在接触点中，查询距离O点最近的接触点
Figure 702437DEST_PATH_IMAGE034
，以
Figure 841294DEST_PATH_IMAGE036
为界将xoy平面区域划分为C 3区、C 4区，通过下式可以判断测点的归属区域，
Figure 201310241194X100001DEST_PATH_IMAGE037
其中，
Figure 346094DEST_PATH_IMAGE038
为参考点，等于减
Figure 784345DEST_PATH_IMAGE010
的横坐标，
Figure 94104DEST_PATH_IMAGE040
等于
Figure 587664DEST_PATH_IMAGE021
减的纵坐标，
Figure 201310241194X100001DEST_PATH_IMAGE041
，；当时，表示
Figure 899194DEST_PATH_IMAGE016
位于C 3区，当
Figure 644165DEST_PATH_IMAGE044
时，表示
Figure 235683DEST_PATH_IMAGE016
位于C 4区；
步骤11：分别连接有效接触点
Figure 195549DEST_PATH_IMAGE034
与C 3区内各个测点
Figure 201310241194X100001DEST_PATH_IMAGE045
，求各个连线与当前包容区域中心线的夹角集合，夹角的最小值即为旋转角度
Figure 784793DEST_PATH_IMAGE046
，
Figure 201310241194X100001DEST_PATH_IMAGE047
并记录C 3区中对应
Figure 564531DEST_PATH_IMAGE046
的测点
Figure 648155DEST_PATH_IMAGE048
，将该
Figure 462528DEST_PATH_IMAGE048
与点
Figure 284990DEST_PATH_IMAGE034
坐标代入下式，
Figure DEST_PATH_IMAGE049
得到迭代后的包容区域垂直方向矢量
Figure 489707DEST_PATH_IMAGE010
，并判断
Figure 688607DEST_PATH_IMAGE010
迭代前后的夹角是否大于90度，如果大于90度，则取反向，跳转到步骤10；
步骤12：根据
Figure 537800DEST_PATH_IMAGE010
计算整个评定区域的对称度误差
Figure 292129DEST_PATH_IMAGE026
；如果
Figure 232404DEST_PATH_IMAGE033
等于，则跳转到步骤15；如果
Figure 920054DEST_PATH_IMAGE033
不等于
Figure 850095DEST_PATH_IMAGE026
，则跳转到步骤13；
步骤13：计算各测点
Figure 390798DEST_PATH_IMAGE050
到中心线的距离，并记录各测点到中心线的最小距离，各测点到中心线的最小距离对应的测点即为与包容区域接触的测点；
如果接触点不止1个且接触点分布在O点的两侧，则计算整个评定区域的对称度误差
Figure 768689DEST_PATH_IMAGE026
，并跳转到步骤15；如果接触点分布在一侧或只有1个接触点，则在接触点中，查询距离O点最近的接触点；按照步骤10中的方式判断所有测点的归属区域；
步骤14：按照步骤11的方法，计算
Figure 770460DEST_PATH_IMAGE010
的值；跳转到步骤13；
步骤15：输出对称度误差
Figure 114854DEST_PATH_IMAGE026
，以及被测特征的参数
Figure 534203DEST_PATH_IMAGE008
。
CN201310241194.XA 2013-06-18 2013-06-18 一种基于最小区域的对称度误差评定方法 Expired - Fee Related CN103292773B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310241194.XA CN103292773B (zh) 2013-06-18 2013-06-18 一种基于最小区域的对称度误差评定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310241194.XA CN103292773B (zh) 2013-06-18 2013-06-18 一种基于最小区域的对称度误差评定方法

Publications (2)

Publication Number Publication Date
CN103292773A true CN103292773A (zh) 2013-09-11
CN103292773B CN103292773B (zh) 2015-11-25

Family

ID=49094034

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310241194.XA Expired - Fee Related CN103292773B (zh) 2013-06-18 2013-06-18 一种基于最小区域的对称度误差评定方法

Country Status (1)

Country Link
CN (1) CN103292773B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104751005A (zh) * 2015-04-16 2015-07-01 王灿 一种基于正交实验的平面度误差评定方法
CN107169681A (zh) * 2017-06-29 2017-09-15 南京航空航天大学 误差分区表征方法
CN108267106A (zh) * 2017-12-30 2018-07-10 唐哲敏 一种快稳简的圆柱度误差评定方法
CN108286957A (zh) * 2017-12-30 2018-07-17 唐哲敏 一种快稳简的平面度误差评定方法
CN108562258A (zh) * 2017-12-30 2018-09-21 唐哲敏 一种快稳简的最大内接圆柱直径评定方法
CN108663009A (zh) * 2018-05-16 2018-10-16 大连大学 连杆类零件的给定方向相对基准要素轴线的平行度误差测量方法
CN115157437A (zh) * 2022-06-28 2022-10-11 中电建路桥集团有限公司 一种预制箱梁质量控制的标准化及数据化方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5986760A (en) * 1997-06-12 1999-11-16 Nikon Corporation Shape measurement method and high-precision lens manufacturing process
JP2002213930A (ja) * 2000-11-14 2002-07-31 Nikon Corp 形状測定方法および高精度レンズの製造方法
JP2004101454A (ja) * 2002-09-12 2004-04-02 Nikon Corp 絶対校正法及び光学素子
CN101393004A (zh) * 2006-12-27 2009-03-25 刘其和 轴键槽专用卡尺

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5986760A (en) * 1997-06-12 1999-11-16 Nikon Corporation Shape measurement method and high-precision lens manufacturing process
JP2002213930A (ja) * 2000-11-14 2002-07-31 Nikon Corp 形状測定方法および高精度レンズの製造方法
JP2004101454A (ja) * 2002-09-12 2004-04-02 Nikon Corp 絶対校正法及び光学素子
CN101393004A (zh) * 2006-12-27 2009-03-25 刘其和 轴键槽专用卡尺

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
杨旭平: "分析法体现基准对称度误差的优化计算", 《广西机械》 *
苑国英等: "对称度误差的测量和评定", 《计量技术》 *
黄美发等: "表面形貌评定方法对比分析", 《机械设计》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104751005A (zh) * 2015-04-16 2015-07-01 王灿 一种基于正交实验的平面度误差评定方法
CN104751005B (zh) * 2015-04-16 2018-03-23 桂林航天工业学院 一种基于正交实验的平面度误差评定方法
CN107169681A (zh) * 2017-06-29 2017-09-15 南京航空航天大学 误差分区表征方法
CN107169681B (zh) * 2017-06-29 2018-01-16 南京航空航天大学 误差分区表征方法
CN108267106A (zh) * 2017-12-30 2018-07-10 唐哲敏 一种快稳简的圆柱度误差评定方法
CN108286957A (zh) * 2017-12-30 2018-07-17 唐哲敏 一种快稳简的平面度误差评定方法
CN108562258A (zh) * 2017-12-30 2018-09-21 唐哲敏 一种快稳简的最大内接圆柱直径评定方法
CN108663009A (zh) * 2018-05-16 2018-10-16 大连大学 连杆类零件的给定方向相对基准要素轴线的平行度误差测量方法
CN115157437A (zh) * 2022-06-28 2022-10-11 中电建路桥集团有限公司 一种预制箱梁质量控制的标准化及数据化方法
CN115157437B (zh) * 2022-06-28 2024-01-16 中电建路桥集团有限公司 一种预制箱梁质量控制的标准化及数据化方法

Also Published As

Publication number Publication date
CN103292773B (zh) 2015-11-25

Similar Documents

Publication Publication Date Title
CN103292773B (zh) 一种基于最小区域的对称度误差评定方法
CN103256916B (zh) 一种基于最小区域的零件平面度误差评定方法
CN104316016B (zh) 一种整体叶盘复杂曲面叶片的纵向测量方法
CN103278126B (zh) 一种基于最小区域的零件球度误差评定方法
CN106248035A (zh) 一种基于点云模型的面轮廓度精确检测的方法及系统
CN101413793A (zh) 复合空间型面几何误差评定方法
CN108682043A (zh) 一种基于参数映射的复杂曲面测量规划方法
WO2021128614A1 (zh) 一种基于特征线的弧面凸轮廓面误差测量与评定方法
CN108332685B (zh) 一种编码结构光三维测量方法
CN103292769B (zh) 一种基于最小区域的平面倾斜度误差检测方法
CN107330142A (zh) 一种基于stl模型的在机检测的测点法矢估算方法
CN103471572A (zh) 全站仪组网测量大型结构件的方法
CN102162728B (zh) 变椭圆活塞裙部横截面线轮廓误差最小区域的测量方法
CN104392030B (zh) 一种基于stl三维模型的曲面加工方法
CN103292654A (zh) 一种计算圆柱体零件作用尺寸的方法
CN104680016A (zh) 基于几何优化逼近的抛物线轮廓最小区域拟合方法
CN103177192B (zh) 五孔探针角度特性曲线数据处理方法
US10132623B2 (en) Method for measuring slant wall thickness dimension of hub
CN102359783B (zh) 基于视觉的移动机器人定位方法
CN104180789A (zh) 基于图形匹配算法的叶片检测方法
CN111409109A (zh) 工业机器人稳定性的测量和评估方法
CN104751005A (zh) 一种基于正交实验的平面度误差评定方法
CN102564383B (zh) 基于轴面角确定三轴磁传感器轴间角的方法
CN106092005A (zh) 三坐标机测量自由曲线的方法
CN103292674B (zh) 一种计算球体零件作用尺寸的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C41 Transfer of patent application or patent right or utility model
CB03 Change of inventor or designer information

Inventor after: Zhong Yanru

Inventor after: Chen Leilei

Inventor after: Sun Yonghou

Inventor after: Fan Kai

Inventor after: Gong Wenfeng

Inventor after: Huang Meifa

Inventor after: Tang Zhemin

Inventor before: Zhong Yanru

Inventor before: Chen Leilei

Inventor before: Sun Yonghou

Inventor before: Huang Meifa

Inventor before: Fan Kai

COR Change of bibliographic data
TA01 Transfer of patent application right

Effective date of registration: 20151029

Address after: 541004 East District School of mechanical and electrical engineering, Guilin University of Electronic Technology, 1 Jinji Road, Qixing District, Guilin, the Guangxi Zhuang Autonomous Region

Applicant after: GUILIN University OF ELECTRONIC TECHNOLOGY

Address before: 541004 graduate school, Guilin University of Electronic Technology, 1 Jinji Road, Qixing District, Guilin, the Guangxi Zhuang Autonomous Region

Applicant before: Chen Leilei

C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151125