CN104751005A - 一种基于正交实验的平面度误差评定方法 - Google Patents

一种基于正交实验的平面度误差评定方法 Download PDF

Info

Publication number
CN104751005A
CN104751005A CN201510178359.2A CN201510178359A CN104751005A CN 104751005 A CN104751005 A CN 104751005A CN 201510178359 A CN201510178359 A CN 201510178359A CN 104751005 A CN104751005 A CN 104751005A
Authority
CN
China
Prior art keywords
plane
value
flatness
coordinate system
point
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
CN201510178359.2A
Other languages
English (en)
Other versions
CN104751005B (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 Aerospace 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 CN201510178359.2A priority Critical patent/CN104751005B/zh
Publication of CN104751005A publication Critical patent/CN104751005A/zh
Application granted granted Critical
Publication of CN104751005B publication Critical patent/CN104751005B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

一种基于正交实验的平面度误差评定方法，该方法首先测量并获取被测平面上测点坐标；然后根据被测点中三个边缘点确定初始参考平面，并进行坐标系变换，使初始参考平面与xoy平面平行且通过坐标原点，此时理想参考平面z=ax+by中a、b的取值都在[-1,1]之间，创建两个参数三个水平的L9正交表，按照平面度的计算方法，计算九种不同a、b下的平面度值，找出平面度最小点，以该点为中心、a、b取值范围缩小一半创建新的正交表，重复计算比较，直到a、b的搜索范围小于设定要求或多次迭代计算的最小平面度不再变化，则停止搜索，此时的最小平面度误差值即为最优值，输出平面度误差以及相应的平面参数。本发明可准确计算出满足最小区域的平面度误差以及平面参数。

Description

一种基于正交实验的平面度误差评定方法
技术领域
本发明涉及一种基于正交实验的零件平面度误差评定方法，属于精密计量与计算机应用领域，可用于各种情况下几何产品中平面度指标的合格性检测，并为加工过程与加工工艺的改进提供指导。
背景技术
平面是机械零件中最常见的几何要素之一，平面的精度对产品的质量、性能以及装配具有重要的影响。最小包容区域法通过寻找距离最小包容平行平面来获得平面平面度，其评定结果接近理想值，符合国际标准对平面度的定义，目前，国内外学者采用的方式主要有凸包法、计算几何法、遗传算法等。这些方法或多或少存在计算稳定性差、计算效率不高、对采点数量的限制、计算结果精确不高等缺陷，导致最小区域法很难在实际检测中应用。目前市场上一般都采用成熟的最小二乘法近似地计算平面度误差。
发明内容
本发明的主要目的是克服现有平面度误差计算方法中存在的不足，设计了一种基于正交试验法的平面度误差评定方法。本方法不仅提高了平面度误差评定精度，而且算法稳定性好、计算效率高，可以推广应用于其它形状误差评定中。
本发明依据正交试验法以少数分布均匀且整齐可比的点代替多数点、从而简化计算过程的特点，将平面度的评定过程转化为平面矢量方向（a,b,-1）中a、b的快速搜索过程。最终实现平面度误差的快速精密评定。本发明主要包括以下步骤：
步骤1：将被测零件置于测量平台上，在测量空间直角坐标系中获取被测平面的测点坐标，P i （x i ,y i ,z i ），（i=1,2，…，n，n为测量点总数且大于3）；
步骤2：确定初始参考平面
根据被测平面的测点坐标，计算其在X轴、Y轴和Z轴方向的最大坐标差，以其中较大的两个坐标差为搜索方向找出被测平面上的三个边缘点P a （x a ,y a ,z a ）、P b （x b ,y b ,z b ）、P c （x c ,y c ,z c ），通过三边缘点的面为初始参考平面，并求出其方向矢量；
步骤3：坐标系平移旋转
将坐标系移至第一个边缘点，然后进行旋转，使新坐标系的XOY平面与初始参考平面平行，即初始参考平面通过新坐标系的原点、方向矢量为（0,0,1），并计算各测量点在新坐标系中的坐标；
步骤4：在新坐标系中确定实际参考平面z=ax+by中的系数a、b的初始取值范围
此时a、b的搜索区域为以（0,0）为中心，边长为2的正方形，即a、b的取值都在[-1,1]之间；
步骤5：以a、b为参数，创建L9正交表，根据点到面的距离公式计算正交表中9种不同a、b取值情况下各测量点到参考平面z=ax+by的距离                                                和的极值差f(a,b)，并找出极值差f(a,b)最小时对应的a、b值；
步骤6：判断正方形边长是否小于1.0e-7，若是，则跳转到步骤8，否则跳转到步骤7；
步骤7：缩小搜索范围并进行新的计算
以步骤5中找出的f(a,b)最小值时的a、b值为中心，创建边长缩小一半的正方形搜索区域；跳转到步骤5；
步骤8：输出最后一次正交表中平面度误差值f(a,b)的最小值，其值即为最终评定出的平面度误差值T。将此时的平面方向矢量（a,b,-1）进行步骤3中的坐标逆变换即为原坐标系中评定平面的矢量方向。
 本发明的有益效果在于：本方法将正交实验用于平面度误差的评定中，有效的加快了搜索过程，同时在精密要求增加时，计算量增加非常少，能有效、快速实现平面度误差的快速精密评定。
对于本领域技术人员来说，根据和应用本发明公开的构思，能够容易地对本发明方案进行各种变形和改变，应当注意的是，所有这些变形和改变都应当属于本发明的范围。
 具体实施方式
一种计算平面度误差的方法，如附图1所示，包括如下步骤：
步骤1：将被测零件置于测量平台上，在测量空间直角坐标系中获取被测平面的测点坐标，P i （x i ,y i ,z i ），（i=1,2，…，n，n为测量点总数且大于3）；
步骤2：确定初始参考平面
根据被测平面的测点坐标，计算其在X轴、Y轴和Z轴方向的最大坐标差，以其中较大的两个坐标差为搜索方向找出被测平面上的三个边缘点P a （x a ,y a ,z a ）、P b （x b ,y b ,z b ）、P c （x c ,y c ,z c ），通过三边缘点的面为初始参考平面，并求出其方向矢量；
步骤3：坐标系平移旋转
将坐标系移至第一个边缘点，然后进行旋转，使新坐标系的XOY平面与初始参考平面平行，即初始参考平面通过新坐标系的原点、方向矢量为（0,0,1），并计算各测量点在新坐标系中的坐标；
步骤4：在新坐标系中确定实际参考平面z=ax+by中的系数a、b的初始取值范围
此时a、b的搜索区域为以（0,0）为中心，边长为2的正方形，即a、b的取值都在[-1,1]之间；
步骤5：以a、b为参数，创建L9正交表，根据点到面的距离公式计算正交表中9种不同a、b取值情况下各测量点到参考平面z=ax+by的距离和的极值差f(a,b)，并找出极值差f(a,b)最小时对应的a、b值；
步骤6：判断正方形边长是否小于1.0e-7，若是，则跳转到步骤8，否则跳转到步骤7；
步骤7：缩小搜索范围并进行新的计算
以步骤5中找出的f(a,b)最小值时的a、b值为中心，创建边长缩小一半的正方形搜索区域；跳转到步骤5；
步骤8：输出最后一次正交表中平面度误差值f(a,b)的最小值，其值即为最终评定出的平面度误差值T。将此时的平面方向矢量（a,b,-1）进行步骤3中的坐标逆变换即为原坐标系中评定平面的矢量方向。
 附图说明
图1 为本发明的平面度的计算流程图。

Claims (1)

1.一种基于正交试验的平面度误差评定方法，具体步骤如下：
步骤1：将被测零件置于测量平台上，在测量空间直角坐标系中获取被测平面的测点坐标，P i （x i ,y i ,z i ），（i=1,2，…，n，n为测量点总数且大于3）；
步骤2：确定初始参考平面
根据被测平面的测点坐标，计算其在X轴、Y轴和Z轴方向的最大坐标差，以其中较大的两个坐标差为搜索方向找出被测平面上的三个边缘点P a （x a ,y a ,z a ）、P b （x b ,y b ,z b ）、P c （x c ,y c ,z c ），通过三边缘点的面为初始参考平面，并求出其方向矢量；
步骤3：坐标系平移旋转
将坐标系移至第一个边缘点，然后进行旋转，使新坐标系的XOY平面与初始参考平面平行，即初始参考平面通过新坐标系的原点、方向矢量为（0,0,1），并计算各测量点在新坐标系中的坐标；
步骤4：在新坐标系中确定实际参考平面z=ax+by中的系数a、b的初始取值范围
此时a、b的搜索区域为以（0,0）为中心，边长为2的正方形，即a、b的取值都在[-1,1]之间；
步骤5：以a、b为参数，创建L9正交表，根据点到面的距离公式计算正交表中9种不同a、b取值情况下各测量点到参考平面z=ax+by的距离                                                和的极值差f(a,b)，并找出极值差f(a,b)最小时对应的a、b值；
步骤6：判断正方形边长是否小于1.0e-7，若是，则跳转到步骤8，否则跳转到步骤7；
步骤7：缩小搜索范围并进行新的计算
以步骤5中找出的f(a,b)最小值时的a、b值为中心，创建边长缩小一半的正方形搜索区域；跳转到步骤5；
步骤8：输出最后一次正交表中平面度误差值f(a,b)的最小值，其值即为最终评定出的平面度误差值T；将此时的平面方向矢量（a,b,-1）进行步骤3中的坐标逆变换即为原坐标系中评定平面的矢量方向。
CN201510178359.2A 2015-04-16 2015-04-16 一种基于正交实验的平面度误差评定方法 Expired - Fee Related CN104751005B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510178359.2A CN104751005B (zh) 2015-04-16 2015-04-16 一种基于正交实验的平面度误差评定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510178359.2A CN104751005B (zh) 2015-04-16 2015-04-16 一种基于正交实验的平面度误差评定方法

Publications (2)

Publication Number Publication Date
CN104751005A true CN104751005A (zh) 2015-07-01
CN104751005B CN104751005B (zh) 2018-03-23

Family

ID=53590680

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510178359.2A Expired - Fee Related CN104751005B (zh) 2015-04-16 2015-04-16 一种基于正交实验的平面度误差评定方法

Country Status (1)

Country Link
CN (1) CN104751005B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105841640A (zh) * 2016-04-29 2016-08-10 北京航空航天大学 平面度误差评定方法及装置
CN117053687A (zh) * 2023-08-17 2023-11-14 广州市西克传感器有限公司 一种基于激光线扫3d相机的电芯高度段差检测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1710533A1 (en) * 2005-04-06 2006-10-11 Fanuc Ltd Work installation error measuring apparatus
CN102853757A (zh) * 2012-09-24 2013-01-02 大连海事大学 一种平面形状误差在线测量系统和方法
CN103256916A (zh) * 2013-06-10 2013-08-21 陈磊磊 一种基于最小区域的零件平面度误差评定方法
CN103292773A (zh) * 2013-06-18 2013-09-11 陈磊磊 一种基于最小区域的对称度误差评定方法
CN104050372A (zh) * 2014-06-18 2014-09-17 北京工业大学 一种三维几何形状误差自动评定方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1710533A1 (en) * 2005-04-06 2006-10-11 Fanuc Ltd Work installation error measuring apparatus
CN102853757A (zh) * 2012-09-24 2013-01-02 大连海事大学 一种平面形状误差在线测量系统和方法
CN103256916A (zh) * 2013-06-10 2013-08-21 陈磊磊 一种基于最小区域的零件平面度误差评定方法
CN103292773A (zh) * 2013-06-18 2013-09-11 陈磊磊 一种基于最小区域的对称度误差评定方法
CN104050372A (zh) * 2014-06-18 2014-09-17 北京工业大学 一种三维几何形状误差自动评定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王如高 等: "基于正交试验的ＳＬＳ工艺参数多指标优化", 《安徽工程大学学报》 *
雷凯 等: "平面度误差评定算法及程序设计", 《机床与液压》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105841640A (zh) * 2016-04-29 2016-08-10 北京航空航天大学 平面度误差评定方法及装置
CN105841640B (zh) * 2016-04-29 2018-09-18 北京航空航天大学 平面度误差评定方法及装置
CN117053687A (zh) * 2023-08-17 2023-11-14 广州市西克传感器有限公司 一种基于激光线扫3d相机的电芯高度段差检测方法

Also Published As

Publication number Publication date
CN104751005B (zh) 2018-03-23

Similar Documents

Publication Publication Date Title
CN108536648B (zh) 基于多超声波传感器的局部放电非线性模型转换求解与优化方法
Li et al. A modified ICP algorithm based on dynamic adjustment factor for registration of point cloud and CAD model
US10094792B2 (en) Method for analysis of thermal resistance
CN103256916B (zh) 一种基于最小区域的零件平面度误差评定方法
CN105474377B (zh) 代表性目标子集的选择及使用
CN102445174B (zh) 一种基于支持向量回归的多测点平面度评定方法
CN103292773B (zh) 一种基于最小区域的对称度误差评定方法
CN103049596B (zh) 一种提高叶轮零件逆向工程精度的方法
CN103292769B (zh) 一种基于最小区域的平面倾斜度误差检测方法
CN102506805B (zh) 一种基于支持向量分类的多测点平面度评定方法
CN104890131B (zh) 一种基于面形误差斜率的确定性修形加工方法
CN105627935A (zh) 产品料厚三维分析系统及方法
CN104197823A (zh) 射频四极场电极极头三维空间曲面的测量方法
CN104751005A (zh) 一种基于正交实验的平面度误差评定方法
Di Leo et al. Uncertainty evaluation of camera model parameters
Xu et al. Improvement of Monte Carlo method for robot workspace solution and volume calculation
CN105157655A (zh) 一种基于区域搜索的圆度误差快速评定方法
CN104112277A (zh) 一种基于双路径Radon变换的测试曲线的拐点间距计算方法
US8682946B1 (en) Robust peak finder for sampled data
CN102564383B (zh) 基于轴面角确定三轴磁传感器轴间角的方法
CN111177848B (zh) 一种基于有限元模型的应变理论值的获取方法和装置
CN103310106B (zh) 一种零件孔系作用尺寸的计算方法
CN107945279B (zh) 一种评价服装褶裥等级的方法
CN113052918A (zh) 对极式双目相机标定误差评价方法、装置、介质及设备
CN109883896A (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
TA01 Transfer of patent application right

Effective date of registration: 20180208

Address after: 541004 the Guangxi Zhuang Autonomous Region Guilin Qixing District No. 2 Jinji Road

Applicant after: GUILIN University OF AEROSPACE TECHNOLOGY

Address before: No. 74, Jinji Road, Seven Star District, Guilin, the Guangxi Zhuang Autonomous Region

Applicant before: Wang Can

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180323

CF01 Termination of patent right due to non-payment of annual fee