CN103292769B - 一种基于最小区域的平面倾斜度误差检测方法 - Google Patents
一种基于最小区域的平面倾斜度误差检测方法 Download PDFInfo
- Publication number
- CN103292769B CN103292769B CN201310241235.5A CN201310241235A CN103292769B CN 103292769 B CN103292769 B CN 103292769B CN 201310241235 A CN201310241235 A CN 201310241235A CN 103292769 B CN103292769 B CN 103292769B
- Authority
- CN
- China
- Prior art keywords
- plane
- contact point
- point
- calculate
- containment 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.)
- Active
Links
Landscapes
- Image Analysis (AREA)
Abstract
一种基于最小区域的平面倾斜度误差检测方法，该方法首先测量并获取被测平面上测点坐标；然后给出被测平面的初始参数，查询被测要素与误差包容区域接触的测点，根据接触点的坐标确定包容区域旋转的方向，并通过作辅助圆锥的方法计算变动量，获得变动后的包容区域参数；然后在查询接触点，计算变动量，依次迭代计算，直到满足判别准则，输出平面倾斜度误差以及平面参数的最优值。本发明可准确计算出满足最小区域的平面倾斜度误差以及平面参数最优值。
Description
技术领域
本发明涉及一种基于最小区域的平面倾斜度误差检测方法，属于精密计量与计算机应用领域，可用于各种情况下几何产品中平面倾斜度指标的合格性检测，并为加工过程与加工工艺的改进提供指导。
背景技术
平面是机械零件中最常见的几何要素之一，平面相对基准的倾斜精度对产品的质量、性能以及装配具有重要的影响。根据国际标准ISO/1101的规定，在计算零件误差时应符合最小条件，即最小区域。在相同的硬件检测条件下，应用最小区域计算零件误差，可提高检测设备的检测精度。
最小区域法，属于不可微复杂最优化问题，目前，国内外学者主要采用传统优化方法、智能算法、计算几何方法等。这些方法由于存在计算稳定性差、计算效率低、对采点数量有限制、计算结果难以达到精确解等缺陷，导致最小区域法很难在实际检测中应用。目前市场上一般都采用成熟的最小二乘法近似地计算平面倾斜度误差。
发明内容
本发明的主要目的是克服现有平面倾斜度误差计算方法中存在的不足，设计了一种基于最小区域法的平面倾斜度误差检测方法。本方法不仅提高了平面倾斜度误差评定精度，而且算法稳定性好、计算效率高，可以推广应用于其它形状误差评定中。
本发明依据平面倾斜度自身的特征，在优化的方向上，对包容区域进行相应旋转变动，最终计算出平面参数最优值和平面倾斜度误差最优值。本发明主要包括以下步骤：
步骤1：将被测零件置于测量平台上，在测量空间直角坐标系中分别获取基准特征与被测平面的测点坐标。
步骤2：判断基准特征是平面还是直线；
如果基准特征为平面，则应用最小区域的方法拟合基准平面，得到基准平面的法向矢量；如果基准特征为直线，则应用最小区域的方法拟合基准直线，得到基准直线的方向矢量；将基准平面的法向矢量与基准直线的方向矢量统一记作基准矢量。
步骤3：随机给出被测平面的初始法向矢量，并保证与成图纸规定的理想角度；进行坐标变换，使坐标系z轴与矢量平行，同时被测平面的测点与作相同的变换；然后沿z轴正向进行坐标平移，使所有的测点均位于xoy平面之上；设坐标变换后测点的坐标为，其中=1,2,…,n，n代表测点数目且为大于3的正整数，所有测点形成测点集；变换后判断的第3个分量是否小于0，如果小于0，则将赋值给，设调整后的为；以每个测点为锥顶分别作辅助圆锥，辅助圆锥的半锥角等于理想角度的余角，辅助圆锥轴线的方向矢量等于(0,0,1)；包容区域与xoy平面相交为平行直线，计算xoy平面内包容直线的方向矢量，计算在xoy坐标平面上的投影。
步骤4：各个测点对应的辅助圆锥与xoy平面相交为圆，计算各个相交圆的半径；应用下式，依次计算各个测点到包容区域的距离，
并分别记录各个测点到包容区域的最大距离和最小距离所对应的测点，所记录的最小距离对应的测点形成包容区域的低值接触点集合，所记录的最大距离对应的测点形成包容区域的高值接触点集合；最大距离与最小距离之差为平面的倾斜度误差。
步骤5：判断低值接触点数量、高值接触点数量是否均为1；
如果低值接触点、高值接触点数量不均为1，则跳转到步骤6；
如果低值接触点、高值接触点数量均为1，则低值接触点、高值接触点均为有效接触点，跳转到步骤7。
步骤6：在xoy平面内，计算接触点对应的圆与平行直线的切点，进行坐标变换，使与坐标系的纵坐标平行，然后计算所有高值接触点的横坐标最大值和最小值，计算所有低值接触点的横坐标最大值和最小值；
如果大于等于，且小于等于，则表明符合判别准则，则跳转到步骤9；如果小于，则分别查询、所对应的高值接触点与低值接触点，所查询到的接触点为有效接触点；如果大于，则分别查询、所对应的高值接触点与低值接触点，所查询到的接触点为有效接触点。
步骤7：计算有效高值接触点与有效低值接触点对应圆的外公切线；公切线将xoy平面划分为两个区域、；设有效高值接触点为，计算对应圆与包容区域的切点为；设有效低值接触点为，计算对应圆与包容区域的切点为；此时包容区域的旋转变动的方向矢量等于。
步骤8：依次计算的各测点对应的圆与对应的圆的外公切线；然后分别计算的各测点对应的圆与对应的圆的外公切线；然后依次计算所有公切线与夹角，并查询最小的夹角以及对应的测点；最小夹角即为包容区域的变动量；对包容区域旋转角度，并重新计算；跳转到步骤4。
步骤9：通过与理论夹角，计算被测平面的包容平面的法向矢量，并计算被测平面的倾斜度误差。
本发明的有益效果在于：本方法查询到的接触点相对位置满足最小区域判别准则，计算得到的实际平面对应理想平面的参数和平面倾斜度误差均为最优值，计算稳定性好、计算效率高。
对于本领域技术人员来说，根据和应用本发明公开的构思，能够容易地对本发明方案进行各种变形和改变，应当注意的是，所有这些变形和改变都应当属于本发明的范围。
附图说明
图1为本发明的平面倾斜度的计算流程图。
具体实施方式
一种基于最小区域的平面倾斜度误差检测方法，如附图1所示，包括如下步骤：
步骤1：将被测零件置于测量平台上，在测量空间直角坐标系中分别获取基准特征与被测平面的测点坐标。
步骤2：判断基准特征是平面还是直线；
如果基准特征为平面，则应用最小区域的方法拟合基准平面，得到基准平面的法向矢量；如果基准特征为直线，则应用最小区域的方法拟合基准直线，得到基准直线的方向矢量；将基准平面的法向矢量与基准直线的方向矢量统一记作基准矢量。
步骤3：随机给出被测平面的初始法向矢量，并保证与成图纸规定的理想角度；进行坐标变换，使坐标系z轴与矢量平行，同时被测平面的测点与作相同的变换；然后沿z轴正向进行坐标平移，使所有的测点均位于xoy平面之上；设坐标变换后测点的坐标为，其中=1,2,…,n，n代表测点数目且为大于3的正整数，所有测点形成测点集；变换后判断的第3个分量是否小于0，如果小于0，则将赋值给，设调整后的为；以每个测点为锥顶分别作辅助圆锥，辅助圆锥的半锥角等于理想角度的余角，辅助圆锥轴线的方向矢量等于(0,0,1)；包容区域与xoy平面相交为平行直线，计算xoy平面内包容直线的方向矢量，计算在xoy坐标平面上的投影。
步骤4：各个测点对应的辅助圆锥与xoy平面相交为圆，计算各个相交圆的半径；应用下式，依次计算各个测点到包容区域的距离，
并分别记录各个测点到包容区域的最大距离和最小距离所对应的测点，所记录的最小距离对应的测点形成包容区域的低值接触点集合，所记录的最大距离对应的测点形成包容区域的高值接触点集合；最大距离与最小距离之差为平面的倾斜度误差。
步骤5：判断低值接触点数量、高值接触点数量是否均为1；
如果低值接触点、高值接触点数量不均为1，则跳转到步骤6；
如果低值接触点、高值接触点数量均为1，则低值接触点、高值接触点均为有效接触点，跳转到步骤7。
步骤6：在xoy平面内，计算接触点对应的圆与平行直线的切点，进行坐标变换，使与坐标系的纵坐标平行，然后计算所有高值接触点的横坐标最大值和最小值，计算所有低值接触点的横坐标最大值和最小值；
如果大于等于，且小于等于，则表明符合判别准则，则跳转到步骤9；如果小于，则分别查询、所对应的高值接触点与低值接触点，所查询到的接触点为有效接触点；如果大于，则分别查询、所对应的高值接触点与低值接触点，所查询到的接触点为有效接触点。
步骤7：计算有效高值接触点与有效低值接触点对应圆的外公切线；公切线将xoy平面划分为两个区域、；设有效高值接触点为，计算对应圆与包容区域的切点为；设有效低值接触点为，计算对应圆与包容区域的切点为；此时包容区域的旋转变动的方向矢量等于。
步骤8：依次计算的各测点对应的圆与对应的圆的外公切线；然后分别计算的各测点对应的圆与对应的圆的外公切线；然后依次计算所有公切线与夹角，并查询最小的夹角以及对应的测点；最小夹角即为包容区域的变动量；对包容区域旋转角度，并重新计算；跳转到步骤4。
步骤9：通过与理论夹角，计算被测平面的包容平面的法向矢量，并计算被测平面的倾斜度误差。
Claims (1)
1.一种基于最小区域的平面倾斜度误差检测方法，其特征在于，具体步骤如下：
步骤1：将被测零件置于测量平台上，在测量空间直角坐标系中分别获取基准特征与被测平面的测点坐标；
步骤2：判断基准特征是平面还是直线；
如果基准特征为平面，则应用最小区域的方法拟合基准平面，得到基准平面的法向矢量；如果基准特征为直线，则应用最小区域的方法拟合基准直线，得到基准直线的方向矢量；将基准平面的法向矢量与基准直线的方向矢量统一记作基准矢量；
步骤3：随机给出被测平面的初始法向矢量，并保证与成图纸规定的理想角度；进行坐标变换，使坐标系z轴与矢量平行，同时被测平面的测点与作相同的变换；然后沿z轴正向进行坐标平移，使所有的测点均位于xoy平面之上；设坐标变换后测点的坐标为，其中=1,2,…,n，n代表测点数目且为大于3的正整数，所有测点形成测点集；变换后判断的第3个分量是否小于0，如果小于0，则将赋值给，设调整后的为；以每个测点为锥顶分别作辅助圆锥，辅助圆锥的半锥角等于理想角度的余角，辅助圆锥轴线的方向矢量等于(0,0,1)；包容区域与xoy平面相交为平行直线，计算xoy平面内包容直线的方向矢量，计算在xoy坐标平面上的投影；
步骤4：各个测点对应的辅助圆锥与xoy平面相交为圆，计算各个相交圆的半径；应用下式，依次计算各个测点到包容区域的距离，
并分别记录各个测点到包容区域的最大距离和最小距离所对应的测点，所记录的最小距离对应的测点形成包容区域的低值接触点集合，所记录的最大距离对应的测点形成包容区域的高值接触点集合；最大距离与最小距离之差为平面的倾斜度误差；
步骤5：判断低值接触点数量、高值接触点数量是否均为1；
如果低值接触点、高值接触点数量不均为1，则跳转到步骤6；
如果低值接触点、高值接触点数量均为1，则低值接触点、高值接触点均为有效接触点，跳转到步骤7；
步骤6：在xoy平面内，计算接触点对应的圆与平行直线的切点，进行坐标变换，使与坐标系的纵坐标平行，然后计算所有高值接触点的横坐标最大值和最小值，计算所有低值接触点的横坐标最大值和最小值；
如果大于等于，且小于等于，则表明符合判别准则，则跳转到步骤9；如果小于，则分别查询、所对应的高值接触点与低值接触点，所查询到的接触点为有效接触点；如果大于，则分别查询、所对应的高值接触点与低值接触点，所查询到的接触点为有效接触点；
步骤7：计算有效高值接触点与有效低值接触点对应圆的外公切线；公切线将xoy平面划分为两个区域、；设有效高值接触点为，计算对应圆与包容区域的切点为；设有效低值接触点为，计算对应圆与包容区域的切点为；此时包容区域的旋转变动的方向矢量等于；
步骤8：依次计算的各测点对应的圆与对应的圆的外公切线；然后分别计算的各测点对应的圆与对应的圆的外公切线；然后依次计算所有公切线与夹角，并查询最小的夹角以及对应的测点；最小夹角即为包容区域的变动量；对包容区域旋转角度，并重新计算；跳转到步骤4；
步骤9：通过与理论夹角，计算被测平面的包容平面的法向矢量，并计算被测平面的倾斜度误差。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310241235.5A CN103292769B (zh) | 2013-06-19 | 2013-06-19 | 一种基于最小区域的平面倾斜度误差检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310241235.5A CN103292769B (zh) | 2013-06-19 | 2013-06-19 | 一种基于最小区域的平面倾斜度误差检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103292769A CN103292769A (zh) | 2013-09-11 |
| CN103292769B true CN103292769B (zh) | 2015-11-25 |
Family
ID=49094031
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310241235.5A Active CN103292769B (zh) | 2013-06-19 | 2013-06-19 | 一种基于最小区域的平面倾斜度误差检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103292769B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107169681B (zh) * | 2017-06-29 | 2018-01-16 | 南京航空航天大学 | 误差分区表征方法 |
| CN108286957A (zh) * | 2017-12-30 | 2018-07-17 | 唐哲敏 | 一种快稳简的平面度误差评定方法 |
| CN108267106B (zh) * | 2017-12-30 | 2020-11-27 | 桂林电子科技大学 | 一种快稳简的圆柱度误差评定方法 |
| CN108562258A (zh) * | 2017-12-30 | 2018-09-21 | 唐哲敏 | 一种快稳简的最大内接圆柱直径评定方法 |
| CN110285781B (zh) * | 2019-07-10 | 2023-09-29 | 桂林电子科技大学 | 一种相对于基准面的平面平行度快速评定方法 |
| CN112781518B (zh) * | 2021-01-20 | 2022-07-19 | 广东百年基业水利水电建设有限公司 | 一种房屋形变监测方法及系统 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1467476A (zh) * | 2002-06-07 | 2004-01-14 | 富士写真光机株式会社 | 样品的倾斜误差调整的辅助方法 |
| CN101010644A (zh) * | 2004-09-03 | 2007-08-01 | 日矿金属株式会社 | 板状材料加工面的确定方法、加工方法及其装置 |
| JP2010117260A (ja) * | 2008-11-13 | 2010-05-27 | Epson Toyocom Corp | 姿勢検出装置の補正パラメーター作成方法、姿勢検出装置の補正パラメーター作成用装置及び姿勢検出装置 |
| CN101750044A (zh) * | 2008-11-28 | 2010-06-23 | 红塔烟草(集团)有限责任公司 | 一种使用测针中心坐标和矢量补偿精密测量3d轮廓面的方法 |
| CN102162728A (zh) * | 2010-12-17 | 2011-08-24 | 南京工程学院 | 变椭圆活塞裙部横截面线轮廓误差最小区域评定方法 |
| CN103115601A (zh) * | 2013-02-19 | 2013-05-22 | 南京工程学院 | 轴类零件圆柱度的超差测定方法 |
| KR101266394B1 (ko) * | 2011-12-16 | 2013-05-22 | 한국과학기술연구원 | 거리측정 센서의 입사각 영향에 따른 오차 보정방법 |
-
2013
- 2013-06-19 CN CN201310241235.5A patent/CN103292769B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1467476A (zh) * | 2002-06-07 | 2004-01-14 | 富士写真光机株式会社 | 样品的倾斜误差调整的辅助方法 |
| CN101010644A (zh) * | 2004-09-03 | 2007-08-01 | 日矿金属株式会社 | 板状材料加工面的确定方法、加工方法及其装置 |
| JP2010117260A (ja) * | 2008-11-13 | 2010-05-27 | Epson Toyocom Corp | 姿勢検出装置の補正パラメーター作成方法、姿勢検出装置の補正パラメーター作成用装置及び姿勢検出装置 |
| CN101750044A (zh) * | 2008-11-28 | 2010-06-23 | 红塔烟草(集团)有限责任公司 | 一种使用测针中心坐标和矢量补偿精密测量3d轮廓面的方法 |
| CN102162728A (zh) * | 2010-12-17 | 2011-08-24 | 南京工程学院 | 变椭圆活塞裙部横截面线轮廓误差最小区域评定方法 |
| KR101266394B1 (ko) * | 2011-12-16 | 2013-05-22 | 한국과학기술연구원 | 거리측정 센서의 입사각 영향에 따른 오차 보정방법 |
| CN103115601A (zh) * | 2013-02-19 | 2013-05-22 | 南京工程学院 | 轴类零件圆柱度的超差测定方法 |
Non-Patent Citations (1)
| Title |
|---|
| 表面形貌评定方法对比分析;黄美发,等;《机械设计》;20120531;第29卷(第5期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103292769A (zh) | 2013-09-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103292769B (zh) | 一种基于最小区域的平面倾斜度误差检测方法 | |
| CN108536648B (zh) | 基于多超声波传感器的局部放电非线性模型转换求解与优化方法 | |
| CN103256916B (zh) | 一种基于最小区域的零件平面度误差评定方法 | |
| CN107167169B (zh) | 基于机器视觉系统的指针式仪表读数识别测量方法 | |
| CN103292773B (zh) | 一种基于最小区域的对称度误差评定方法 | |
| CN102663767B (zh) | 视觉测量系统的相机参数标定优化方法 | |
| CN103278126B (zh) | 一种基于最小区域的零件球度误差评定方法 | |
| CN102622747B (zh) | 一种用于视觉测量的摄像机参数优化方法 | |
| WO2020168685A1 (zh) | 一种三维扫描视点规划方法、装置及计算机可读存储介质 | |
| CN106248035A (zh) | 一种基于点云模型的面轮廓度精确检测的方法及系统 | |
| CN104344795A (zh) | 点云剖面量测系统及方法 | |
| CN104869641A (zh) | 基于AP优化的Wi-Fi室内定位方法 | |
| CN105631939B (zh) | 一种基于曲率滤波的三维点云畸变校正方法及其系统 | |
| CN103471519A (zh) | 一种应用无棱镜光电全站仪测量输变电塔架变形的方法 | |
| CN103278095B (zh) | 一维可伸缩式靶标装置及其标定特征点的产生方法 | |
| CN102679926B (zh) | 用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法 | |
| CN105157655A (zh) | 一种基于区域搜索的圆度误差快速评定方法 | |
| CN201576308U (zh) | 一种复杂曲面的配准系统 | |
| WO2013121196A1 (en) | Measurement apparatus and method | |
| CN104680016A (zh) | 基于几何优化逼近的抛物线轮廓最小区域拟合方法 | |
| CN113510748B (zh) | 基于双目视觉的机器人空间位置重复性和准确度测量方法 | |
| CN108021037B (zh) | 三元组目标阵列近场效应修正方法 | |
| CN103292655B (zh) | 一种有基准约束的圆柱体的作用尺寸的计算方法 | |
| CN104680586A (zh) | 空间任意位置椭球面最小区域拟合方法 | |
| CN104751005A (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: Bao Jiading Inventor after: Tang Zhemin Inventor after: Sun Yingchao Inventor after: Huang Meifa Inventor after: Fan Kai Inventor after: Zhan Fan Inventor after: Tang Yirui Inventor before: Zhong Yanru Inventor before: Chen Leilei Inventor before: Gong Wenfeng Inventor before: Sun Ying Inventor before: Sun Yingchao |
|
| 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 |