CN108067939B - 一种空间点位测量基准误差补偿方法 - Google Patents
一种空间点位测量基准误差补偿方法 Download PDFInfo
- Publication number
- CN108067939B CN108067939B CN201611022909.2A CN201611022909A CN108067939B CN 108067939 B CN108067939 B CN 108067939B CN 201611022909 A CN201611022909 A CN 201611022909A CN 108067939 B CN108067939 B CN 108067939B
- Authority
- CN
- China
- Prior art keywords
- measurement
- error
- coordinate value
- space point
- point location
- 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
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q15/00—Automatic control or regulation of feed movement, cutting velocity or position of tool or work
- B23Q15/007—Automatic control or regulation of feed movement, cutting velocity or position of tool or work while the tool acts upon the workpiece
- B23Q15/12—Adaptive control, i.e. adjusting itself to have a performance which is optimum according to a preassigned criterion
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种空间点位测量基准误差补偿方法。该方法可有效解决零件实际测量坐标系与理论测量坐标系不符,导致测量结果存在测量基准误差的问题。在进行计算时仅考虑测量法矢方向测量结果的误差,依据最小二乘法构建目标函数,求解实测值与理论值之间的坐标变换矩阵,进而对实测值进行测量基准误差补偿。本发明方法可实现对测量基准误差的有效补偿,提高测量精度,为零件加工精度的评估提供有效数据依据。
Description
技术领域
本发明涉及一种误差补偿方法,尤其是一种用于空间点位测量误差补偿方法,具体地说是一种用于空间点位测量的基准误差补偿方法。
背景技术
零件在数控加工过程中,为获取零件当前加工状态,通常采用在线检测或三坐标测量机测量的方式对零件点位进行测量。在零件测量过程中,由于零件理论测量坐标系与实际测量坐标系间存在偏差,使得零件空间点位测量结果存在测量基准误差。这种误差将会对零件加工精度的评估带来很大的影响。
查阅相关技术和文献发现,蔺小军(2013)在学术期刊《计量学报》2013,34(2),p128-133发表了论文“叶片型面测量编程与误差处理技术”公开了一种消除测量***误差的方法,基于已知叶片CAD模型提取测量点理论坐标值,采用微平面法计算法向矢量,并采用ICP算法进行配准消除***误差,该方法可在保证测量精度的前提下提高叶片零件测量效率。
刘元朋(2005)在学术期刊《中国机械工程》2005,16(12),p1080-1082发表了论文“复杂曲面测量数据最佳匹配问题研究”公开了一种复杂曲面类零件测量数据匹配方法,通过初试匹配和精确匹配来实现测量数据的最佳匹配,其中精确匹配依据最小二乘原理构造目标函数,应用L-BFGS-B算法进行精确匹配,可有效解决复杂曲面类零件测量数据与曲面的匹配问题。
使用探头进行测量时,探头沿测量点法矢方向对测量点误差进行测量。如图1所示,当零件实际测量坐标系与理论坐标系不一致,采用探头测量时,测量误差仅存在于测量法矢方向。上述方法在进行测量数据匹配时,对问题进行了简化,未考虑测量法矢对测量误差的影响,采用上述方法无法对测量基准误差进行有效补偿。
发明内容
本发明的目的是针对零件测量时由于实际测量坐标系与理论测量坐标系之间的偏差导致测量结果存在测量基准误差的问题,发明了一种空间点位测量基准误差补偿方法。
本发明的技术方案是:
一种空间点位测量基准误差补偿方法,包括以下步骤:
步骤1,从测量机或机床获取零件空间测量点的实际测量坐标值;
步骤2,基于零件设计模型及测量文件,获取零件空间测量点理论坐标值及每个测量点的测量法矢;
步骤3,将空间点位实际测量坐标值,理论坐标值,不同点处的测量法矢及测量误差分别定义为pa,pt,v和e;
步骤4,假定将空间点位理论值pt经三次平移及三次旋转后得到p′t,其中平移旋转量为δx,δy,δz,εx,εy,εz;
步骤5,将p′t中每点与平移旋转前对应点在测量方向上的误差同该点实际测量值误差的差值平方和设为最小二乘法目标函数,用于求解步骤2中的平移旋转量δx,δy,δz,εx,εy,εz;
步骤6,当目标函数值最小时,空间点位平移旋转量即为空间点位测量实际坐标系与理论坐标系间的位置关系,为求解空间点位平移旋转量,将目标函数对平移旋转量δx,δy,δz,εx,εy,εz分别求导并赋值为0;
步骤7,将求导方程联立得到方程组,求解空间点位平移旋转量δx,δy,δz,εx,εy,εz;
步骤8,根据空间点位平移旋转量,将空间点位理论值进行平移和旋转得到变换后的空间点位坐标值p′t;
步骤9,根据步骤6求解的变换后空间点位坐标值p′t及空间点位实际测量值pa沿测量方向求解经基准误差补偿后的空间点位误差e′;
步骤10,依据误差e′及空间点位理论值pt求解经基准误差补偿后的空间点位测量值p′a。
步骤11,将经基准误差补偿后的空间点位测量值输入测量机测量***或机床***作为判断零件加工误差的依据。
所述的空间点位理论值经平移旋转后得到的p′t在平移旋转量较小,省略计算过程中的高阶无穷小量后,可用下式表示:
其中δx,δy,δz,εx,εy,εz分别为沿X,Y,Z方向的平移量和旋转量;
所述的最小二乘法目标函数可用下式表示:
所述的目标函数对平移旋转量的求导过程计算如下:
所述的平移旋转量求解过程如下:
将求导方程联立得到方程组如下式所示:
所述的经基准误差补偿后的空间点位测量误差e′可表示为:
所述的经基准误差补偿后的空间点位测量值p′a可表示为:
针对空间点位测量误差仅存在于测量法矢方向的特点,本发明以测量法矢方向误差为最小二乘法目标函数,可准确求解空间点位变换矩阵,实现对空间点位测量基准误差的补偿。
本发明的有益效果如下:
1、通过考虑测量法矢对测量误差的影响,构建最小二乘法目标函数,求解空间点位实测坐标系与理论坐标系间的变换矩阵,进而实现对测量结果的基准误差补偿;
2、通过实现对测量结果的基准误差补偿,可有效提高测量精度,为零件加工精度的评估提供数据依据。
附图说明
图1为零件测量误差示意图。
其中:1.零件理论位置,2.零件实际位置,3.测量探头,4.零件测量点位实际空间位置,5.零件测量点处的测量法矢,6.零件测量点位理论空间位置,7.沿测量法矢方向的零件测量误差。
图2为零件测量点位示意图。
其中:8.标准S试件,9.S试件上的空间测量点位,10.空间点位处的测量法矢。
表1为S试件空间测量点位理论坐标、测量法矢及公差值。
表2为S试件空间点位经误差补偿前后的坐标值与误差值。
具体实施方式
下面结合附图和实施例对本发明提出的空间点位测量基准误差补偿方法进行说明,但本发明专利并不限于本实例。
以专利“综合检测数控铣床精度的"S"形检测试件及其检测方法(ZL200710048269.7)”中提出的S试件为例说明本发明提出的方法,如图2中8所示。在S试件上选取15个空间测量点位,测量点位理论坐标、测量法矢及公差如表1所示。
表1
空间点位 | p<sub>t</sub>/理论坐标值/mm | v/测量法矢 | 公差/mm |
点1 | (-52.975,-70.603,23) | (-0.976,0.071,0.206) | ±0.05 |
点2 | (-51.491,-20.917,23) | (-0.973,-0.006,0.233) | ±0.05 |
点3 | (-50.243,28.747,23) | (-0.979,0.126,0.161) | ±0.05 |
点4 | (-28.083,71.213,23) | (-0.482,0.858,0.178) | ±0.05 |
点5 | (-4.153,76.97,23) | (-0.017,0.982,0.189) | ±0.05 |
点6 | (20.317,73.269,23) | (0.306,0.935,0.179) | ±0.05 |
点7 | (55.853,40.676,23) | (0.895,0.425,0.137) | ±0.05 |
点8 | (70.763,-6.683,23) | (0.967,0.256,-0.010) | ±0.05 |
点9 | (87.823,-53.218,23) | (0.853,0.507,-0.123) | ±0.05 |
点10 | (126.998,-81.098,23) | (0.186,0.972,-0.140) | ±0.05 |
点11 | (151.731,-81.533,23) | (-0.159,0.978,-0.136) | ±0.05 |
点12 | (174.337,-71.932,23) | (-0.638,0.753,-0.160) | ±0.05 |
点13 | (191.549,-26.65,23) | (-0.987,0.074,-0.143) | ±0.05 |
点14 | (191.776,23.054,23) | (-0.981,-0.006,-0.195) | ±0.05 |
点15 | (193.367,72.737,23) | (-0.979,0.070,-0.190) | ±0.05 |
2、假定将空间点位理论值经平移旋转后得到p′t,在平移旋转量较小,省略计算过程中的高阶无穷小量后,可用下式表示:
其中δx,δy,δz,εx,εy,εz分别为沿X,Y,Z方向的平移量和旋转量;
3、将平移旋转后得到的p′t中每点与平移旋转前对应点在测量方向上的误差同该点实际测量值误差的差值平方和设为最小二乘法目标函数,用下式表示:
4、当目标函数值最小时,空间点位平移旋转量即为空间点位测量实际坐标系与理论坐标系间的位置关系,为求解空间点位平移旋转量,将目标函数对平移旋转量δx,δy,δz,εx,εy,εz分别求导并将值赋为0;
5、将求导方程联立得到方程组,求解空间点位平移旋转量,用下式表示:
经上式计算得到空间点位平移旋转量结果如下:
(δx,δy,δz,εx,εy,εz)=(0.0168,-0.037,0.4834,-0.0008,-0.0007,0.0004)
6、根据空间点位平移旋转量,得到变换后的空间点位坐标值p′t,进而求解p′t与空间点位实测值pa沿测量方向的误差,即经基准误差补偿后的空间点位误差e′,如表2所示,补偿前平均误差为0.066mm,经误差补偿后平均误差为0.017mm;
7、依据误差e′及空间点位理论值pt求解经基准误差补偿后的空间点位测量值p′a,如表2所示。
表2
8、将经基准误差补偿后的空间点位测量值输入测量机测量***作为S试件加工误差判断依据。
本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。
Claims (4)
1.一种空间点位测量基准误差补偿方法,包括以下步骤:
步骤1,从测量机或机床获取零件空间点位实际测量坐标值;
步骤2,基于零件设计模型及测量文件,获取零件空间点位理论坐标值及每个测量点的测量法矢;
步骤3,将空间点位实际测量坐标值,空间点位理论坐标值,不同点处的测量法矢及测量误差分别定义为pa,pt,v和e;
步骤4,假定将空间点位理论坐标值pt经三次平移及三次旋转后得到p′t,其中平移旋转量为δx,δy,δz,εx,εy,εz;
步骤5,将p′t中每点与平移旋转前对应点在测量方向上的误差同该点实际测量值误差的差值平方和设为最小二乘法目标函数,用于求解步骤4中的平移旋转量δx,δy,δz,εx,εy,εz;
步骤6,当目标函数值最小时,平移旋转量即为空间点位测量实际坐标系与理论坐标系间的位置关系,为求解平移旋转量,将目标函数对平移旋转量δx,δy,δz,εx,εy,εz分别求导并赋值为0;
步骤7,将求导方程联立得到方程组,求解平移旋转量δx,δy,δz,εx,εy,εz;
步骤8,根据平移旋转量,将空间点位理论坐标值进行平移和旋转得到变换后的空间点位坐标值p′t;
步骤9,根据步骤8求解的变换后空间点位坐标值p′t及空间点位实际测量坐标值pa沿测量方向求解经基准误差补偿后的空间点位测量误差e′;
步骤10,依据误差e′及空间点位理论坐标值pt求解经基准误差补偿后的空间点位实际测量坐标值p′a;
步骤11,将经基准误差补偿后的空间点位测量值输入测量机测量***或机床***作为判断零件加工误差的依据;
所述的空间点位理论坐标值经平移旋转后得到的p′t,省略计算过程中的高阶无穷小量后,可用下式表示:
其中δx,δy,δz,εx,εy,εz分别为沿X,Y,Z方向的平移量和旋转量;
所述的目标函数对平移旋转量的求导过程计算如下:
所述的最小二乘法目标函数可用下式表示:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611022909.2A CN108067939B (zh) | 2016-11-18 | 2016-11-18 | 一种空间点位测量基准误差补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611022909.2A CN108067939B (zh) | 2016-11-18 | 2016-11-18 | 一种空间点位测量基准误差补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108067939A CN108067939A (zh) | 2018-05-25 |
CN108067939B true CN108067939B (zh) | 2021-08-03 |
Family
ID=62160375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611022909.2A Active CN108067939B (zh) | 2016-11-18 | 2016-11-18 | 一种空间点位测量基准误差补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108067939B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109443273B (zh) * | 2018-09-28 | 2020-08-07 | 易思维(杭州)科技有限公司 | 利用三维测量***对待测工件进行精确定位的方法 |
CN110421406B (zh) * | 2019-07-14 | 2021-04-20 | 深圳市烨嘉为技术有限公司 | 基于偏心差控制的刀具动态自适应补偿方法 |
CN111561924B (zh) * | 2020-05-21 | 2022-08-30 | 哈尔滨工业大学 | 一种磁信标的校正方法及基于旋转磁偶极子的定位方法 |
CN112084599B (zh) * | 2020-09-09 | 2022-02-11 | 南京航空航天大学 | 一种面向航空航天箱体部件结构变形的补偿方法 |
CN112729214B (zh) * | 2020-11-27 | 2022-06-14 | 成都飞机工业(集团)有限责任公司 | 一种基于试验数据修正基准点坐标的测量方法 |
CN116931507B (zh) * | 2023-09-18 | 2024-01-12 | 成都飞机工业(集团)有限责任公司 | 一种群孔穿孔控制方法、装置、存储介质及电子设备 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0063466A2 (en) * | 1981-04-17 | 1982-10-27 | Toyota Jidosha Kabushiki Kaisha | Method for detecting a position of a moving body |
US4705016A (en) * | 1985-05-17 | 1987-11-10 | Disco Abrasive Systems, Ltd. | Precision device for reducing errors attributed to temperature change reduced |
CN101000285A (zh) * | 2007-01-16 | 2007-07-18 | 成都飞机工业(集团)有限责任公司 | 综合检测数控铣床精度的“s”形检测试件及其检测方法 |
CN101545775A (zh) * | 2009-05-05 | 2009-09-30 | 东南大学 | 利用数字地图计算像片方位元素和建筑物高度的方法 |
CN102168991A (zh) * | 2011-01-29 | 2011-08-31 | 中北大学 | 三轴矢量传感器与安装载体间安装误差的标定补偿方法 |
CN102840870A (zh) * | 2012-09-19 | 2012-12-26 | 中天启明石油技术有限公司 | 三维正交方向传感器几何误差的修正方法 |
CN105458372A (zh) * | 2015-12-29 | 2016-04-06 | 北京理工大学 | 基于非可展直纹面的侧铣误差补偿装置及其刀位规划方法 |
CN105538038A (zh) * | 2016-01-27 | 2016-05-04 | 清华大学 | 机床平动轴几何误差辨识方法 |
-
2016
- 2016-11-18 CN CN201611022909.2A patent/CN108067939B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0063466A2 (en) * | 1981-04-17 | 1982-10-27 | Toyota Jidosha Kabushiki Kaisha | Method for detecting a position of a moving body |
US4705016A (en) * | 1985-05-17 | 1987-11-10 | Disco Abrasive Systems, Ltd. | Precision device for reducing errors attributed to temperature change reduced |
CN101000285A (zh) * | 2007-01-16 | 2007-07-18 | 成都飞机工业(集团)有限责任公司 | 综合检测数控铣床精度的“s”形检测试件及其检测方法 |
CN101545775A (zh) * | 2009-05-05 | 2009-09-30 | 东南大学 | 利用数字地图计算像片方位元素和建筑物高度的方法 |
CN102168991A (zh) * | 2011-01-29 | 2011-08-31 | 中北大学 | 三轴矢量传感器与安装载体间安装误差的标定补偿方法 |
CN102840870A (zh) * | 2012-09-19 | 2012-12-26 | 中天启明石油技术有限公司 | 三维正交方向传感器几何误差的修正方法 |
CN105458372A (zh) * | 2015-12-29 | 2016-04-06 | 北京理工大学 | 基于非可展直纹面的侧铣误差补偿装置及其刀位规划方法 |
CN105538038A (zh) * | 2016-01-27 | 2016-05-04 | 清华大学 | 机床平动轴几何误差辨识方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108067939A (zh) | 2018-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108067939B (zh) | 一种空间点位测量基准误差补偿方法 | |
CN108563186B (zh) | 一种五轴球头铣削几何误差补偿方法 | |
JP5244786B2 (ja) | 差分較正 | |
Tsutsumi et al. | Identification of angular and positional deviations inherent to 5-axis machining centers with a tilting-rotary table by simultaneous four-axis control movements | |
Nojedeh et al. | Tool path accuracy enhancement through geometrical error compensation | |
CN109032069B (zh) | 一种采用电涡流位移传感器的非接触式R-test测量仪球心坐标计算方法 | |
Wan et al. | Calibration and compensation of machine tool volumetric error using a laser tracker | |
CN110125455B (zh) | 一种用于机器人钻孔中优化钻头位姿的方法 | |
Linares et al. | Impact of measurement procedure when error mapping and compensating a small CNC machine using a multilateration laser interferometer | |
JP6955991B2 (ja) | 空間精度補正方法、及び空間精度補正装置 | |
JP5968749B2 (ja) | 幾何誤差同定方法と当該幾何誤差同定方法を使用した数値制御方法、数値制御装置及びマシニングセンタ | |
KR101593330B1 (ko) | 볼바와 진직도 데이터를 이용한 다축 정밀 제어 기계의 직각도 측정 방법 | |
Sortino et al. | Compensation of geometrical errors of CAM/CNC machined parts by means of 3D workpiece model adaptation | |
Camboulives et al. | Calibration of a 3D working space by multilateration | |
Aguado et al. | Improving a real milling machine accuracy through an indirect measurement of its geometric errors | |
Zimmermann et al. | Self-calibration of rotary axis and linear axes error motions by an automated on-machine probing test cycle | |
US11982529B2 (en) | Method for registering a total station in the reference system of a CAD model | |
CN105094051B (zh) | 运动平台***的平面定位补偿方法 | |
Acosta et al. | Verification of machine tools using multilateration and a geometrical approach | |
KR102035334B1 (ko) | 4축 공작기계의 기하학적 오차 측정 방법 | |
Barari | Automotive body inspection uncertainty associated with computational processes | |
CN107421476A (zh) | 一种空间孔位测量基准误差补偿方法 | |
Khameneifar et al. | A new methodology for evaluating position and orientation errors of airfoil sections | |
Gąska et al. | Analysis of changes in coordinate measuring machines accuracy made by different nodes density in geometrical errors correction matrix | |
Zargarbashi et al. | A model based method for centering double ball bar test results preventing fictitious ovalization effects |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |