CN104697478A - 基于恒定接触力的直纹面自适应测量方法 - Google Patents

基于恒定接触力的直纹面自适应测量方法 Download PDF

Info

Publication number
CN104697478A
CN104697478A CN201510114726.2A CN201510114726A CN104697478A CN 104697478 A CN104697478 A CN 104697478A CN 201510114726 A CN201510114726 A CN 201510114726A CN 104697478 A CN104697478 A CN 104697478A
Authority
CN
China
Prior art keywords
measurement
contact force
probe
along
measured
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
CN201510114726.2A
Other languages
English (en)
Other versions
CN104697478B (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.)
Shaanxi Hanji Precision Machinery Co ltd
Xiamen University
Original Assignee
Xiamen University
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 Xiamen University filed Critical Xiamen University
Priority to CN201510114726.2A priority Critical patent/CN104697478B/zh
Publication of CN104697478A publication Critical patent/CN104697478A/zh
Application granted granted Critical
Publication of CN104697478B publication Critical patent/CN104697478B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/20Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring contours or curvatures, e.g. determining profile

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • A Measuring Device Byusing Mechanical Method (AREA)

Abstract

基于恒定接触力的直纹面自适应测量方法,涉及曲面零件测量。首先对被测表面进行网格划分,对被测复杂曲面单元节点进行接触力和测头坐标的获取,通过次梯度的大小和方向变化对其恒定接触力的探针位移变化进行预测;在进行表面测量时重复利用平行测量路径上单元节点的位置坐标,有效提高测量效率;无需更换测量硬件设备即可自适应完成由直纹面构成的复杂曲面测量,实用性强,可有效提高测量效率减少操作流程。能够实现基于恒定接触力的直纹面自适应的测量方法,同时,使用单元法对测点进行直接半径补偿,使用次梯度法预测第二个接触点的调节位移,在路径规划中,对单元交叉节点测量数据进行重复使用,有效提高了测量的效率。

Description

基于恒定接触力的直纹面自适应测量方法
技术领域
本发明涉及曲面零件测量,尤其是涉及一种基于恒定接触力的直纹面自适应测量方法。
背景技术
坐标测量仪是集机械、控制、计算机技术为一体的大型的精密测量仪器被广泛应用于工业领域中,特别适用于测量曲面等比较难以评定加工形位质量的零件上。测量***作为测量仪的最重要组成部分对测量精度有着重要的影响。测头部分一般由接触式测头11和探针12固连组成(见图1所示,图中标记P为被测曲面),利用导轨实现沿相应方向上的运动并使用光栅尺对其进行移动位置测量,通过测头对被测件进行探测,实现测量任务。
国内坐标测量仪多应用微平面法修正半径,然而这种方法采用“3+2”进行修正,测点太多,测量速度比较慢。深圳大学的全荣等(全荣等.三坐标测量机测头半径对测量数据的影响及其消除方法[J].计量学报,1993,14(4):251一255)提出了一种用函数修正测头半径的方法,这种方法可以更好的补偿精度。武汉理工大学的卢红等(卢红等.测头半径补偿的方法[J].组合机床与自动化加工技术,2001,9:39一41)提出了一种类似于求测头球心等距曲面的方法实现测头半径的补偿。自适应测量方面,王平江(王平江.曲面测量、建模及数控加工集成研究[博士学位论文].武汉:华中理工大学,1996)提出了等弧长均匀网格划分技术然而需要大量的人机交互方法对其进行规划;Y.F.Zhang(YF.Zhang,A.Y.C.Nee,J.YH.Fuh.et al.A Neural Network Approach to Determining Optimal Inspection Sampling Size for CMM[J].Computer Integrated Manufaeturing System,1996,3(9))提出了采用神经网络的方法对其采样点计算,然而需要进行大量加工实验来形成样本。
发明内容
本发明的目的在于提供一种基于恒定接触力的直纹面自适应测量方法。
本发明包括以下步骤:
1)对被测表面进行x、z方向等间隔网格划分;
2)移动测头至被测表面零件起始位置;
3)沿着被测表面法向方向施加一个微小位移,通过探针刚度和位移量获取接触力大小,具体算法为:
Fy=(y1-y0)×ky
其中,Fy为探针和被测表面接触力在坐标轴y轴方向上的分量;ky为探针在测量坐标系y方向上的刚度;y0为探头初始位置的y轴方向坐标值;y1为施加微小位移后的探头y轴方向坐标值;
4)根据零件装夹方式和其几何特征指定其沿着测量坐标系的大致运动路径方向;
5)沿着指定路径方向移动探针到第二测量点附近区域,其位置为起始位置沿测量方向上位移长度为ε,调节接触力与起始位置处Fy相等,沿着指定路径垂直方向移动探针到第三测量点附近区域,其位置为起始位置沿测量方向垂直方向上位移长度为ε,调节接触力与起始位置处Fy相等;
6)计算有限单元法向矢量:
n x = ( y 1 - y 0 ) ( z 2 - z 0 ) - ( y 2 - y 0 ) ( z 1 - z 0 ) n y = ( x 1 - x 0 ) ( z 2 - z 0 ) - ( x 2 - x 0 ) ( z 1 - z 0 ) n z = ( x 1 - x 0 ) ( y 2 - y 0 ) - ( x 2 - x 0 ) ( y 1 - y 0 )
7)计算补偿后测量值:
x 0 ′ = Q x - n x · x 0 y 0 ′ = Q y - n y · y 0 z 0 ′ = Q z - n z · z 0
其中,(x′0,y′0,z′0)为补偿后接触点测量值,(Qx,Qy,Qz)为探针实测点坐标值;
8)计算沿测量方向上次梯度向量变化:
y ( k + 1 ) - y ( k ) α y k = g y ( k )
其中,为x、y、z方向上的步长;
9)根据次梯度变化预测下一测量点x、y、z方向上调节位移量;
Δy k + 1 > Δy k , ( g y ( k + 1 ) · g y ( k ) > 0 , g y ( k + 1 ) > g y ( k ) ) 时,较上一测量点沿y轴向上倾斜;
&Delta;y k + 1 < &Delta;y k , ( g y ( k + 1 ) &CenterDot; g y ( k ) > 0 , g y ( k + 1 ) < g y ( k ) ) 时,较上一测量点沿y轴向下倾斜;
&Delta;y k + 1 = &Delta;y k , ( g y ( k + 1 ) = g y ( k ) ) 时,与上一测量点等高;
时,形成凸点,减小测量步长重新测量;
时,形成凹点,减小测量步长重新测量;
10)测量下一单元至测量结束。
在步骤10)中,测量时,可重复利用上一单元重复单元节点,测量单元2时,可读取单元1的B1坐标值,避免重复测量;测量下一平行路径时,可读取上一条路径上重复的单元结点;在测量单元6~10时,可读取单元1~5的B2、C2、D2、E2、F2点,避免重复测量。
本发明首先由操作人员将探头移动至测量起始位置,探头将测量初始接触力,并保持恒定接触力不变沿着规定测量方向上进行网格化测量,通过对其次梯度的大小和方向的变化达到预测位移调整量的目标。
本发明属于曲面类零件测量技术领域,特别涉及一种通过直接接触式测头来进行接触力的获取并通过次梯度法对其进行探针位移控制以达到准确测量其形位表面的方法。该测量方法首先对被测表面进行网格划分,对被测复杂曲面单元节点进行接触力和测头坐标的获取,通过次梯度的大小和方向变化对其恒定接触力的探针位移变化进行预测。在进行表面测量时重复利用平行测量路径上单元节点的位置坐标,有效提高测量效率。该方法无需更换测量硬件设备即可自适应完成由直纹面构成的复杂曲面测量,实用性强,可有效提高测量效率减少操作流程。
本发明能够实现基于恒定接触力的直纹面自适应的测量方法,同时,使用单元法对测点进行直接半径补偿,使用次梯度法预测第二个接触点的调节位移,在路径规划中,对单元交叉节点测量数据进行重复使用,有效提高了测量的效率。
附图说明
图1为现有的测头组成示意图。
图2为被测直纹面表面。
图3为被测表面网格划分。
图4为单元法矢。
具体实施方式
以下实施例将结合附图对本发明作进一步的说明。
参见图2~4,本发明实施例包括以下步骤:
1)对被测表面进行x、z方向进行等间隔网格划分;
2)移动测头至被测表面零件起始位置(见图3A1点);
3)沿着被测表面法向方向施加一个微小位移,通过探针刚度和位移量获取接触力大小,具体算法为:
Fy=(y1-y0)×ky
其中,Fy为探针和被测表面接触力在坐标轴y轴方向上的分量。ky为探针在测量坐标系y方向上的刚度;y0为探头初始位置的y轴方向坐标值;y1为施加微小位移后的探头y轴方向坐标值;
4)根据零件装夹方式和其几何特征指定其沿着测量坐标系的大致运动路径方向,如图3箭头方向位置所示;
5)沿着指定路径方向移动探针到第二测量点附近区域,其位置为起始位置沿测量方向上位移长度为ε,调节接触力与起始位置处Fy相等。沿着指定路径垂直方向移动探针到第三测量点附近区域,其位置为起始位置沿测量方向垂直方向上位移长度为ε,调节接触力与起始位置处Fy相等;
6)计算有限单元法向矢量(如图4所示):
n x = ( y 1 - y 0 ) ( z 2 - z 0 ) - ( y 2 - y 0 ) ( z 1 - z 0 ) n y = ( x 1 - x 0 ) ( z 2 - z 0 ) - ( x 2 - x 0 ) ( z 1 - z 0 ) n z = ( x 1 - x 0 ) ( y 2 - y 0 ) - ( x 2 - x 0 ) ( y 1 - y 0 )
7)计算补偿后测量值:
x 0 &prime; = Q x - n x &CenterDot; x 0 y 0 &prime; = Q y - n y &CenterDot; y 0 z 0 &prime; = Q z - n z &CenterDot; z 0
其中,(x′0,y′0,z′0)为补偿后接触点测量值,(Qx,Qy,Qz)为探针实测点坐标值。
8)计算沿测量方向上次梯度向量变化:
y ( k + 1 ) - y ( k ) &alpha; y k = g y ( k )
其中,为x、y、z方向上的步长。
9)根据次梯度变化预测下一测量点x、y、z方向上调节位移量。 &Delta;y k + 1 > &Delta;y k , ( g y ( k + 1 ) &CenterDot; g y ( k ) > 0 , g y ( k + 1 ) > g y ( k ) ) 时较上一测量点沿y轴向上倾斜; &Delta;y k + 1 < &Delta;y k , ( g y ( k + 1 ) &CenterDot; g y ( k ) > 0 , g y ( k + 1 ) < g y ( k ) ) 时较上一测量点沿y轴向下倾斜; &Delta;y k + 1 = &Delta;y k , ( g y ( k + 1 ) = g y ( k ) ) 时与上一测量点等高; g y ( k + 1 ) > 0 , g y ( k ) < 0 时,形成凸点,减小测量步长重新测量;时,形成凹点,减小测量步长重新测量。
10)测量下一单元至测量结束。测量时可重复利用上一单元重复单元节点,如图3所示,测量单元2时可读取单元1的B1坐标值,避免重复测量;测量下一平行路径时可读取上一条路径上重复的单元结点;在测量单元6~10时,可读取单元1~5的B2、C2、D2、E2、F2点,避免重复测量。

Claims (2)

1.基于恒定接触力的直纹面自适应测量方法,其特征在于包括以下步骤:
1)对被测表面进行x、z方向等间隔网格划分;
2)移动测头至被测表面零件起始位置;
3)沿着被测表面法向方向施加一个微小位移,通过探针刚度和位移量获取接触力大小,具体算法为:
Fy=(y1-y0)×ky
其中,Fy为探针和被测表面接触力在坐标轴y轴方向上的分量;ky为探针在测量坐标系y方向上的刚度;y0为探头初始位置的y轴方向坐标值;y1为施加微小位移后的探头y轴方向坐标值;
4)根据零件装夹方式和其几何特征指定其沿着测量坐标系的大致运动路径方向;
5)沿着指定路径方向移动探针到第二测量点附近区域,其位置为起始位置沿测量方向上位移长度为ε,调节接触力与起始位置处Fy相等,沿着指定路径垂直方向移动探针到第三测量点附近区域,其位置为起始位置沿测量方向垂直方向上位移长度为ε,调节接触力与起始位置处Fy相等;
6)计算有限单元法向矢量:
n x = ( y 1 - y 0 ) ( z 2 - z 0 ) - ( y 2 - y 0 ) ( z 1 - z 0 ) n y = ( x 1 - x 0 ) ( z 2 - z 0 ) - ( x 2 - x 0 ) ( z 1 - z 0 ) n z = ( x 1 - x 0 ) ( y 2 - y 0 ) - ( x 2 - x 0 ) ( y 1 - y 0 )
7)计算补偿后测量值:
x 0 &prime; = Q x - n x &CenterDot; x 0 y 0 &prime; = Q y - n y &CenterDot; y 0 z 0 &prime; = Q z - n z &CenterDot; z 0
其中,(x′0,y′0,z′0)为补偿后接触点测量值,(Qx,Qy,Qz)为探针实测点坐标值;
8)计算沿测量方向上次梯度向量变化:
y ( k + 1 ) - y ( k ) &alpha; y k = g y ( k )
其中,为x、y、z方向上的步长;
9)根据次梯度变化预测下一测量点x、y、z方向上调节位移量;
&Delta; y k + 1 > &Delta; y k , ( g y ( k + 1 ) &CenterDot; g y ( k ) > 0 , g y ( k + 1 ) > g y ( k ) ) 时,较上一测量点沿y轴向上倾斜;
&Delta; y k + 1 < &Delta; y k , ( g y ( k + 1 ) &CenterDot; g y ( k ) > 0 , g y ( k + 1 ) < g y ( k ) ) 时,较上一测量点沿y轴向下倾斜;
&Delta; y k + 1 = &Delta; y k , ( g y ( k + 1 ) = g k ( k ) ) 时,与上一测量点等高;
时,形成凸点,减小测量步长重新测量;
时,形成凹点,减小测量步长重新测量;
10)测量下一单元至测量结束。
2.如权利要求1所述基于恒定接触力的直纹面自适应测量方法,其特征在于在步骤10)中,测量时,重复利用上一单元重复单元节点,测量单元2时,读取单元1的B1坐标值,避免重复测量;测量下一平行路径时,读取上一条路径上重复的单元结点;在测量单元6~10时,读取单元1~5的B2、C2、D2、E2、F2点,避免重复测量。
CN201510114726.2A 2015-03-16 2015-03-16 基于恒定接触力的直纹面自适应测量方法 Expired - Fee Related CN104697478B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510114726.2A CN104697478B (zh) 2015-03-16 2015-03-16 基于恒定接触力的直纹面自适应测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510114726.2A CN104697478B (zh) 2015-03-16 2015-03-16 基于恒定接触力的直纹面自适应测量方法

Publications (2)

Publication Number Publication Date
CN104697478A true CN104697478A (zh) 2015-06-10
CN104697478B CN104697478B (zh) 2016-12-28

Family

ID=53344850

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510114726.2A Expired - Fee Related CN104697478B (zh) 2015-03-16 2015-03-16 基于恒定接触力的直纹面自适应测量方法

Country Status (1)

Country Link
CN (1) CN104697478B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109084722A (zh) * 2018-06-20 2018-12-25 华中科技大学 一种自适应采样的复杂曲面接触式测量方法
CN109446470A (zh) * 2018-11-16 2019-03-08 厦门大学 基于非接触式检测的弧齿锥齿轮加工用模具磨损预测方法
CN111692995A (zh) * 2020-06-12 2020-09-22 中国航发哈尔滨轴承有限公司 试拉件坐标系的垂直差速测法
CN113009324A (zh) * 2021-05-24 2021-06-22 中国电子科技集团公司第二十九研究所 一种曲面多探针测试治具

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0389062A1 (en) * 1989-03-22 1990-09-26 Koninklijke KPN N.V. Communication system, comprising a central station and a number of decentralised stations connected to that central station via an optical transmission network
EP1378378A1 (fr) * 2002-07-04 2004-01-07 Société de Technologie Michelin Estimation de l'usure d'un pneu.
US6773551B1 (en) * 1999-03-15 2004-08-10 Voith Paper Patent Gmbh Device for determining the characteristics of a running material web
CN102607484A (zh) * 2012-03-27 2012-07-25 西北工业大学 一种变间距自适应搜寻测量点的飞机壁板法矢确定方法
CN102798364A (zh) * 2012-08-17 2012-11-28 河南工业职业技术学院 反求工程中异常测量数据点的剔除方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0389062A1 (en) * 1989-03-22 1990-09-26 Koninklijke KPN N.V. Communication system, comprising a central station and a number of decentralised stations connected to that central station via an optical transmission network
US6773551B1 (en) * 1999-03-15 2004-08-10 Voith Paper Patent Gmbh Device for determining the characteristics of a running material web
EP1378378A1 (fr) * 2002-07-04 2004-01-07 Société de Technologie Michelin Estimation de l'usure d'un pneu.
CN102607484A (zh) * 2012-03-27 2012-07-25 西北工业大学 一种变间距自适应搜寻测量点的飞机壁板法矢确定方法
CN102798364A (zh) * 2012-08-17 2012-11-28 河南工业职业技术学院 反求工程中异常测量数据点的剔除方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张征等: ""接触力学自适应无网格计算***设计"", 《计算机仿真》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109084722A (zh) * 2018-06-20 2018-12-25 华中科技大学 一种自适应采样的复杂曲面接触式测量方法
CN109084722B (zh) * 2018-06-20 2019-08-13 华中科技大学 一种自适应采样的复杂曲面接触式测量方法
CN109446470A (zh) * 2018-11-16 2019-03-08 厦门大学 基于非接触式检测的弧齿锥齿轮加工用模具磨损预测方法
CN111692995A (zh) * 2020-06-12 2020-09-22 中国航发哈尔滨轴承有限公司 试拉件坐标系的垂直差速测法
CN113009324A (zh) * 2021-05-24 2021-06-22 中国电子科技集团公司第二十九研究所 一种曲面多探针测试治具

Also Published As

Publication number Publication date
CN104697478B (zh) 2016-12-28

Similar Documents

Publication Publication Date Title
CN104697478A (zh) 基于恒定接触力的直纹面自适应测量方法
Rezchikov et al. Mathematical models for estimating the degree of influence of major factors on performance and accuracy of coordinate measuring machines
CN102221341B (zh) 基于随机并行梯度下降优化技术的快速数字图像相关测量方法
CN109918614B (zh) 一种基于模态学习的全局动应变测量方法
CN102735386A (zh) 考虑弯曲刚度的斜拉索索力数值计算方法
CN103115593A (zh) 扫描测头标定方法
CN104358194A (zh) 基于激光跟踪仪的轨道静态平顺性测量与分析方法
CN102001024A (zh) 基于加工机床的自由曲面原位测量测量方法
CN103542820B (zh) 一种检测风洞内表面平整度的方法
CN102426170A (zh) 脆性材料亚表面损伤层微裂纹全息反演检测方法
CN103862327A (zh) 一种球头球心位置检测方法
Yang et al. An accurate and fast method to inspect rail wear based on revised global registration
CN103557980B (zh) 体外预应力筋张拉力精确测试方法
Franco et al. Theoretical analysis of straightness errors in coordinate measuring machines (CMM) with three linear axes
CN105937886A (zh) 形状测量装置、加工装置及形状测量装置的校正方法
CN109202539A (zh) 一种复合材料弱刚度异形结构在线检测方法
Fung et al. An approach to on-machine motion error measurement of a linear slide
Sheng et al. A novel bridge curve mode measurement technique based on FOG
CN104597513A (zh) 一种地球物理磁场大数据预处理值的获取方法
Gan et al. Bridge continuous deformation measurement technology based on fiber optic gyro
CN106871784A (zh) 一种激光测量方法、装置以及***
CN102889866A (zh) 以石墨烯键长作为计量基准的长度计量溯源方法
Yuan et al. Calibration and compensation of dynamic Abbe errors of a coordinate measuring machine
Gu et al. Theoretical analysis of crack propagation measurement for brittle materials based on virtual principal strain field
CN108375337B (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: Yao Bin

Inventor after: Sun Weifang

Inventor after: Xiang Mingxin

Inventor after: Zhou Bin

Inventor after: Feng Wei

Inventor after: Shen Zhihuang

Inventor after: Zhang Xianglei

Inventor before: Yao Bin

Inventor before: Sun Weifang

Inventor before: Feng Wei

Inventor before: Shen Zhihuang

Inventor before: Zhang Xianglei

COR Change of bibliographic data
TA01 Transfer of patent application right

Effective date of registration: 20160524

Address after: Xiamen City, Fujian Province, 361005 South Siming Road No. 422

Applicant after: XIAMEN University

Applicant after: SHAANXI HANJI PRECISION MACHINERY CO.,LTD.

Address before: Xiamen City, Fujian Province, 361005 South Siming Road No. 422

Applicant before: Xiamen University

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

Granted publication date: 20161228

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