中国测试  2019, Vol. 45 Issue (9): 21-26

文章信息

王斌
WANG Bin
U型截面薄壁件测量精度优化方法研究
Measurement accuracy optimization method study of the thin wall piece with U shape section
中国测试, 2019, 45(9): 21-26
CHINA MEASUREMENT & TEST, 2019, 45(9): 21-26
http://dx.doi.org/10.11857/j.issn.1674-5124.2019030043

文章历史

收稿日期: 2019-03-11
收到修改稿日期: 2019-05-08
U型截面薄壁件测量精度优化方法研究
王斌     
核工业理化工程研究院,天津 300180
摘要:为提高三坐标测量机在测量一种U型截面薄壁件时的测量精度,分别从温度和测量力两方面研究精度优化方法。基于物体热膨胀效应的原子热振动原理分析铝合金材料的热膨胀影响因素,将温度影响引入测量能力评估算法,得到可实施温度补偿的边界条件,利用Minitab对变温测量实验数据进行分析并拟合出该薄壁件在径向方向上的热膨胀函数,实现在温度边界范围内对测量结果进行精度优化;针对RDS XXT测头系统的内部构造,通过对测头结构进行力学分析并建立采点触发力学模型,确定该测头系统的采点触发模式为力矩触发,利用材料力学理论分析力矩触发模式下的被测件形变机理。测力及测量截面位置交叉测量实验结果表明:测杆长度越长,测力越小,且越靠近工件底部,测量误差越小,测量精度越高。
关键词三坐标测量机    薄壁件    温度补偿    测量变形    测量精度    
Measurement accuracy optimization method study of the thin wall piece with U shape section
WANG Bin     
Research Institute of Physical and Chemical Engineering of Nuclear Industry, Tianjin 300180, China
Abstract: In order to improve the measurement accuracy of thin wall piece with U shape section by CMM, the precision optimization method is studied from two aspects of temperature and measuring force respectively. Based on the principle of atomic thermal vibration, the influence factors of thermal expansion of aluminum alloy materials are analyzed, the temperature influence is introduced into the measurement capability evaluation algorithm, the boundary conditions of temperature compensation can be implemented, the experimental data of variable temperature measurement are analyzed by Minitab, and the thermal expansion function of the thin-walled parts in radial direction is fitted. The temperature compensation of the measurement results can be realized within the temperature boundary range. Aiming at the internal structure of the RDS XXT probe system, through the mechanical analysis of the probe structure and the establishment of the mining point trigger mechanical model, it is determined that the picking point trigger mode of the probe system is torque trigger. Based on the theory of material mechanics, the deformation mechanism of the measured parts under the torque trigger mode is analyzed, and the experimental results show that the measurement accuracy can be improved effectively by selecting the length of the measuring rod and measuring the height of the cross section reasonably.
Key words: CMM     thin wall piece     temperature compensation     measuring deformation     measurement accuracy    
0 引 言

工业坐标测量机是国内外几何量测量领域普遍使用的精密测量设备,其测量精度的影响因素很多,包括静态因素和动态因素两大类[1-2]。静态因素主要包括温度、振动、量值溯源等;动态因素主要包括测量策略、测量力等。通过研究和分析国内外计量测试行业最新技术发展动态发现,为了有效减小测量误差,普遍采取制定国家或行业标准对坐标测量进行规范,如德国的VDA标准,严格规范了坐标机在使用环境、量值溯源等方面的技术要求,并强制贯标执行,使静态误差因素得到有效控制;欧洲的AUKOM坐标测量行业认证机构则通过实行技术资格认证制度,特别是对于形位公差的测量,建立了AUKOM GD&T专项认证标准,进而对坐标机的测量策略进行统一和规范,有效减少了动态因素所导致的误差[3-4]

通过上述分析可以看出,三坐标测量机的测量多数误差来源可通过制定行业标准或加强操作的一致性加以消除,但是对于工件自身温度差异以及结构特点在测量时产生的测量误差,就需要结合材料的热膨胀系数和测力作用机理加以分析并优化,从而提高测量精度。国内外学者在这一领域进行了广泛研究。Valkov等[5]研究了航空铝合金板材的热膨胀系数,认为材料厚度对热膨胀变形有重要影响,并介绍了适于做热膨胀系数测定的制样方法,并采用光学法测量热膨胀系数。Kompan等[6]分析了铝合金产品的不同设计结构与热膨胀系数的关系。上述研究主要解决铝合金材料热膨胀系数的测定精度、测定条件等问题[7],缺乏对采用坐标测量机测量铝合金产品时工件温度与测量精度的关系分析,特别是对于薄壁件,其温度对测量精度的影响和补偿方法,尚未做相关深入研究。Laaouina等[8]研究了接触式测头采点触发时测力在工件表面的作用机理,以及工件由测力所导致的形变规律。Pan等[9]通过分析坐标机广泛采用的压电式触发测头的采点触发模式,建立力学模型,通过计算得到了测力对测量精度的影响关系式。上述研究虽然是针对压电式坐标机测头进行的分析研究,但对于本文研究具有重要的参考和借鉴意义[10]

本文将对一种中心对称剖面呈U型的铝合金薄壁回转零件进行测量精度的温度补偿方法和测力影响消除方法研究。通过建立一种带入温度影响因素的测量能力评估算法,确定温度误差补偿的边界条件并拟合补偿算式;提出一种针对环形触发测头系统测力分析的力学模型,研究测力导致的形变机理,减小测力影响,提高测量精度。

1 工件温度对测量精度的影响分析及补偿方法 1.1 金属热膨胀模型分析及精度优化方法

热膨胀特性是物体的基本属性,由原子的热振动效应产生。根据原子理论模型可知,原子间势能u(r)与原子间的相互距离r存在一定的函数关系,因此,当原子间发生相互振动时,其间距由r0变为r0+x,而势能变为:

$u(r) = u({r_0} + x)$ (1)

由热力学统计原理及热膨胀系数概念,可得到物体的热膨胀系数 $ {\alpha _{{\rm{CTE}},T}}$

$ \begin{split} {\alpha _{{\rm{CTE}},T}} =&\dfrac{1}{r_0} \dfrac{{{\rm{d}}\overline x }}{{{\rm{d}}T}} = \dfrac{1}{{{r_0}}} \left[ \dfrac{{\beta K}}{{{\alpha ^2}}} + \left(\dfrac{{2\delta {K^2}}}{{{\alpha ^3}}} - \dfrac{{16\gamma {\beta ^2}}}{{{\alpha ^4}}}\right)T -\right.\\& \left.\left(\dfrac{{18{\gamma ^2}\beta }}{{{\alpha ^6}}} + \dfrac{{{21{\beta ^2}\delta }}}{{{\alpha ^6}}} + \dfrac{{558\gamma \delta }}{{4{\alpha ^5}}}\right){K^3}{T^2} + \cdots \right] \end{split}$ (2)

上式表征的是不同材料的热膨胀系数随着温度的不同而发生变化,由于铝合金的热膨胀系数受组分比例、金相组织、相间缺陷、相间差异性及加工工艺的热效应等因素影响[11],其热效应很难用经典热膨胀系数计算方法得到。那么,在实际测量时,为研究铝合金的热膨胀特性,通常采用铝合金试件温度变化引起特定尺寸方向伸长量进行测试,绘制关系曲线,进而对不同温度时材料的热膨胀效应进行分析。

1.2 精度优化边界的温度条件

根据测量能力评估方法,测量设备的测量误差至少要小于被测参数设计公差的 $\dfrac{1}{4}$ ,而被测工件的热形变所导致的偏差应与设备自身测量误差合并考虑,累计误差与被测参数设计因满足如下关系[12]

$ \left( {{\delta _{{\text{热形变}}}} + {\delta _{{\text{测量误差}}}}} \right) < T/3 $ (3)

其中, ${\delta _{{\text{热形变}}}} = L \times \Delta T \times C$ ,由公式(3)可以得到

$ \left( {L \times \Delta T \times C + {\delta _{{\text{测量误差}}}}} \right) < T/3 $ (4)

式中:L——零件尺寸,mm;

ΔT——工件温度与测量温度差值,℃;

C——膨胀系数。

那么:

$ \Delta T < (T/3 - {\delta _{{\text{测量误差}}}})/(C \times L) $ (5)

将U型截面薄壁件的直径及设计公差、坐标机的标称精度以及铝合金热膨胀系数带入公式(5),可得到:

$ \Delta T < {7.1^ \circ }{\rm{C}} $

即:当U型截面薄壁件的温度与标准环境温度差小于7.1 ℃时,可以应用温度补偿手段对测量值进行修正。温度差大于该值,则需进行初步恒温再进行测量。考虑到三坐标测量机的工作环境温度为(20±3)℃,因此,在该温度差范围内开展了温度影响实验,以确定最佳的温度补偿算式。

1.3 温度补偿算式拟合

选取100件U型截面薄壁件,分别在17~23 ℃的温度范围内,测量各试件的直径、圆度、高度等设计参数,同时,测量过程中保持测量环境温度恒定。测量实验结果表明,温度影响主要体现在径向尺寸上。为此,以100个试件在不同温度条件下的测量平均值作为拟合元素,利用Minitab数据分析软件进行分析,生成温度系数函数模型如下:

$\phi D = 0.005\;229 \times T + C$ (6)

其中,C为U型薄壁件的名义直径。温度系数函数模型可通过计算得到不同温度条件下的补偿值,可用于测量室在不具备恒温条件下的实时测量结果修正,保证了测量结果的准确性。

1.4 温度补偿结果验证实验

通过选取5个实验件,将其分别恒温至17~23 ℃范围内的各温度点,测量各温度点下的测量值并作为比对数据,利用比对数据与标准温度下(20 ℃)时的测量值(参考数据)进行对比,验证补偿效果,考核测量精度。实验结果如表1所示。

表 1 温度补偿验证实验数据
序号 状态 差值/mm
17℃ 18 ℃ 19 ℃ 20 ℃ 21 ℃ 22 ℃ 23 ℃
1 补偿前 −0.017 −0.011 −0.003 0.002 0.005 0.012 0.016
补偿后 −0.003 0.002 −0.003 0.001 0.002 0.003 −0.002
2 补偿前 −0.016 −0.010 −0.004 0.001 0.006 0.014 0.017
补偿后 0.001 −0.004 0.003 0.002 0.002 0.001 0.003
3 补偿前 −0.020 −0.012 −0.006 0.003 0.008 0.013 0.018
补偿后 0.003 −0.002 0.004 0.002 0.005 0.003 -0.002
4 补偿前 −0.019 −0.013 −0.004 0.002 0.007 0.015 0.021
补偿后 0.001 0.003 −0.002 0.003 0.004 −0.002 0.001
5 补偿前 −0.018 −0.012 −0.005 0.002 0.006 0.014 0.018
补偿后 −0.002 0.003 0.002 0.001 0.003 0.004 −0.002

由表可知,5个测试件中补偿前测量误差最大达到0.021 mm,而补偿后最大误差仅为0.005 mm,测量精度有了提高了4倍多。通过温度补偿,在工件未完全恒温的状态下,有效提高了坐标机的测量精度。

2 测力对精度的影响机理及优化方法 2.1 测头触发模式与测力的关系

三坐标测量机用测头可分为动态测头与静态测头两大类。Zeiss G2坐标机的RDS XXT测头系统是可实现动态(扫描式)和静态(触发式)测量的双模式信号采集系统,采用触发采点测量时,可以看作为静态测量,此时,测量力对工件造成的影响与测头系统的结构及触发模式有关。RDS XXT测头内部结构如图1(a)所示,1为测杆,2为环形保持架。测杆安装在环形保持架上,环形保持架的外壁布有均匀的触点,当测杆没有受到测量力时,环形保持架的外壁与触点均保持接触;当测杆的球状端部与工件接触时,至少会引起一个触点脱离接触,从而引起电路的断开,产生阶跃信号,完成测量触发采点[13]

图 1 测头结构及触发受力示意图

接触测量时,由于测头和工件存在相对运动,当被测表面有几何形状误差或对基准要素有位置误差,有可能使测头瞬时离开被测工件表面而产生测量误差。为了分析方便,将测针和测杆简化为等截面的测杆单元,质量为M,长度为L,作用于工件上的静测量力为F0。当测针相对工件表面测量时,测杆单元绕支点做左右摆动,摆角为θ,则测杆单元的力学方程为

$J\ddot \theta = - \frac{M}{2}\sin \theta + (F - {F_0})L$ (7)

由式(8)可以看出,测头系统的触发力学方程与测杆摆角关系较大,测头系统的结构决定了其摆动距离不大于图1中触点球径,而当θ较小时,该力学模型可简化为力矩方程,即测量力与测杆长度成反比。说明RDS XXT测头系统是由力矩触发采点的,可通过改变测杆长度控制测量力的大小。

图1(b)所示,测头位于坐标系的Z轴的原点,测杆长度为P,环形保持架半径OC的投影与X轴夹角 ${\varphi _0}$ ,测头触碰零件产生的触发力为F(θ,φ),摆角θ和方位角φ直接影响触发力的方向。那么,如果当测头采点方向与零件柱面法线方向相同,极角θ和方位角φ不超过触点半径,则测力 $F$ 可使测杆围绕AB轴旋转,当环形保持架C离开触点时信号触发。因此,为保证采点触发而加载的测力将会使工件产生接触形变。此外,该测量力还会对测针造成弯曲,产生弯曲位移,进而导致工件弯曲变形。

2.2 测力所致形变对测量精度的影响

假设零件和测头不变,那么形变量只与测量力有关[14],RDS XXT测头系统配置有3种不同的采点触发测量力模式,分别为敏感性触发、标准型触发和加力触发,分别对应不同的表面粗糙度状况的工件,敏感性触发模式主要用于测量表面光洁度较高的塑料、玻璃等工件的测量;标准型用于大多数金属材料的机加产品;加力触发用于表面粗糙的复合材料检测。因此,对于铝合金薄壁工件的检测,应采用标准型触发模式,其测力与测杆长度关系曲线如图2所示。

图 2 测力与测杆长度曲线图

由图可以看出,随着测杆长度增加,3种测力模式的测力均呈减小趋势,当测杆长度为30 mm时,3种测力模式的测力差值为0.05 N,而当测杆长度为150 mm时,测力差值为小于0.01 N。

2.2.1 接触形变量计算

使用接触式量仪进行测量时,测头和被测工件的测量面直接接触,为了保证仪器与工件表面的可靠接触,需要一定大小的力作用于工件上,但测量力的存在也会使工件和测头产生接触变形。接触变形量的大小与工件的材料、测量力的大小有关[15]

$ {\delta _{{\text{接触}}}} = 1.231\sqrt[3]{{{{\left( {\frac{F}{E}} \right)}^2}\frac{{{R_1} + {R_2}}}{{{R_1}{R_2}}}}} $ (8)

式中:F——测量力,N;

E——工件的弹性模量,GPa;

R1——工件半径,mm;

R2——测头半径,mm。

将铝合金材料弹性模量,测头球半径为1.5 mm,代入公式(9),得到不同测力条件下的变形量如表2所示。

表 2 测力所致接触形变量计算结果
测杆长度/mm 测力/mN 接触变形量/μm
30 150 1.8
40 110 1.5
50 90 1.3
75 55 0.9
100 40 0.7
125 30 0.6
150 25 0.5

根据以上分析计算可以看出,当测杆长为30 mm时,对应测量力为150 mN,工件的接触变形量为1.8 μm,同时由于测量过程中工件接触的均匀性,考虑到其接触变形量与被测参数公差比值较小以及测量机自身测量误差(1.9 μm),其变形影响可忽略。

2.2.2 弯曲形变量计算

工件在测量力的作用下也会产生弯曲变形,弯曲变形量主要受测量位置以及测量力大小的影响。根据材料力学理论,工件在测量力的作用下产生的弯曲变形量[16]

$ {\delta _{{\text{弯曲}}}} = \frac{{F{L^3}}}{{3EI}} $ (9)

式中:F——测量力,N;

E——工件弹性模量,GPa;

I——工件截面惯性矩,N·m。

工件的弯曲变形不仅和工件的材料有关,其装夹位置的影响也较大。U型截面薄壁件的剖面尺寸结构如图3所示。

图 3 惯性矩简化模型图

这样,可根据下式计算出U型截面薄壁件的惯性矩:

$I = \frac{{Be_1^3 - b{h^3} + ae_2^3}}{3}$ (10)

式中:B——工件外径,mm;

e1——工件外底面至装卡截面距离,mm;

b——工件内径,mm;

h——工件内底面至装卡截面距离,mm;

a/2——工件壁厚,mm;

e2——工件端面至装卡截面距离,mm。

弯曲变形量与装卡位置的关系,用工件装夹高度H与工件直径的比值i作为装夹位置,考虑到该薄壁件的轴向长度,装卡位置差ΔH设置2 mm。选取测力为55 mN,测杆长度为75 mm,将各参数代入公式(10),计算出不同测量截面测量力与弯曲变形量的关系如表3所示。

表 3 测力所致弯曲形变量计算结果
位置ΔH/mm 25 mN 40 mN 55 mN 90 mN 150 mN
2 0.5 0.6 0.8 0.9 1.2
4 0.6 0.9 1.1 1.4 1.7
6 1.2 1.3 1.6 1.9 2.2
8 2.2 2.4 2.5 2.8 3.1
10 4.0 4.2 4.5 5.2 6.2

由表可以看出,随着测量力的增加,其弯曲变形量也在增加,并且随着测量位置的升高,其弯曲变形量最大达到了6.2 μm,可见弯曲变形的影响是造成测量误差的主要原因[17]。所以,在测量速度、测量长度一定的条件下,测量采点位置越靠近工件底部,误差越小,精度越高。

2.3 测力及测量位置与精度关系验证实验

2.3.1 测力与精度关系验证实验

选取5个U型薄壁件,编号记为1#~5#,在同一测量截面高度(ΔH=2 mm)以不同的测量力连续测量10次,记录并计算10次测量结果的极差值,考核其测量精度,测量实验结果如表4所示。

表 4 测力与精度关系验证实验数据
编号 25 mN 40 mN 55 mN 90 mN 150 mN
1# 2.5 3.6 4.8 5.2 6.1
2# 2.1 3.8 4.5 5.8 7.2
3# 1.9 3.2 5.1 4.9 6.7
4# 2.2 3.5 4.4 5.3 6.4
5# 2.3 3.6 4.7 5.1 6.9

由表测量结果可以看出,随着测量力的不断增加,测量误差越来越大,精度逐渐下降。可见,测量力导致的弯曲形变对精度的影响较为显著,与理论计算结果一致。

2.3.2 测量位置与精度关系验证实验

仍选取2.3.1选用的5个试验件,以相同的测量力(40 mN)在不同的截面高度连续测量10次,记录并计算10次测量结果的极差值,考核其测量精度,测量实验结果如表5所示。

表 5 测量位置与精度关系验证实验数据
编号 2 mm 4 mm 6 mm 8 mm 10 mm
1# 1.8 2.5 2.9 3.7 4.9
2# 2.0 2.8 3.2 4.1 5.3
3# 1.7 2.2 2.6 3.4 5.2
4# 2.2 2.8 3.5 4.2 5.5
5# 1.5 2.3 3.0 4.0 5.3

由表可以看出,随着测量截面高度的不断增加,测量误差越来越大,精度逐渐下降。可见,测量位置距离定位位置越远,测力导致的弯曲形变效应越显著,同样与理论计算结果一致。

3 结束语

本文针对U型截面薄壁件在加工过程中会由热效应或环境温度热传导而产生变形的特点,通过变温实验建立了尺寸温度补偿算式,保证测量结果的准确性;通过分析三坐标测量机RDS XXT测头系统的结构及采点触发原理,确定了该测头系统的静态测量力所导致工件形变主要为弯曲变形,对于U型截面薄壁易变形件,测量精度实验结果表明,通过合理配置测杆长度及采点高度,可有效提高此工件的测量精度。

参考文献
[1]
张国雄. 三坐标测量机[M]. 天津: 天津大学出版社, 1999.
[2]
罗哉. 关节臂式坐标测量机测量力误差分析及补偿[J]. 仪器仪表学报, 2017(5): 1159-1166. DOI:10.3969/j.issn.0254-3087.2017.05.014
[3]
YOUNES E, ABDERRAZAK E O, AHMED C. Experimental investigation of dynamic errors in coordinate measuring machines for high speed measurement[J]. International Journal of Precision Engineering and Manufacturing, 2018, 19(8): 1115-1124. DOI:10.1007/s12541-018-0132-x
[4]
HUSSAIN M, GU J, ENGEL R, et al. Designed experiment to find the optimal combination of the factors for the coordinate measuring machine (CMM) to measure cylindricity of engine cylinder bore[C]// International Conference on Industrial Engineering & Operations Management, 2015.
[5]
VALKOV S, NEOV D, BEZDUSHNYI R, et al. Study of the microstructure, crystallographic structure and thermal stability of Al-Ti-Nb alloys produced by selective electron beam alloying[J]. Journal of Surface Investigation: X-ray, Synchrotron and Neutron Techniques, 2018, 12(3): 436-441. DOI:10.1134/S1027451018030187
[6]
KOMPAN T A, KONDRATIEV S V, KORENEV A S, et al. Measurement of the thermal expansion coefficient for ultra-high temperatures up to 3000 K[J]. International Journal of Thermophysics, 2018, 39(3): 40. DOI:10.1007/s10765-017-2353-0
[7]
程湘红. 三坐标测量机测量误差分析及处理[J]. 工程技术研究, 2018(1): 113-114, 160. DOI:10.3969/j.issn.1671-3818.2018.01.053
[8]
LAAOUINA L, NAFI A, MOUCHTACHI A, et al. Application of CMM separation method for identifying absolute values of probe errors and machine errors[C]// International Conference on Engineering & Mis, 2016.
[9]
PAN F, NIE L, BAI Y, et al. Geometric errors measurement for coordinate measuring machines[J]. IOP Conference Series Earth and Environmental Science, 2017, 81(1): 012117.
[10]
KALUS J, JORGENSEN J K. Measuring deformation and mechanical properties of weld lines in short fibre reinforced thermoplastics using digital image correlation[J]. Polymer Testing, 2014, 36: 44-53. DOI:10.1016/j.polymertesting.2014.03.017
[11]
丁鸿章. 工程金属材料热膨胀系数的计算[J]. 浙江工业大学学报, 2000, 28(4): 358-366. DOI:10.3969/j.issn.1006-4303.2000.04.018
[12]
黄晓明, 孙杰, 李剑锋, 等. 7050-T7451铝合金预拉伸板材热膨胀系数变化规律及机理[J]. 中国有色金属学报, 2013(12): 3282-3288.
[13]
刘育楠, 董学军, 刘红刚, 等. 三坐标测量机测头系统精度验证方法研究[J]. 工业计量, 2016, 26(6): 41-45.
[14]
邵竑文, 孟亮, 刘斯宇. 论测量力对几何尺寸测量结果的影响[J]. 计量与测试技术, 2015, 42(12): 52-54.
[15]
李志渤, 黄强先, 赵剑, 等. 接触测量中的微探球力变形研究[J]. 中国机械工程, 2014(14): 1126-1129.
[16]
熊伟, 葛志华, 唐威, 等. 轮毂轴承接触应力及变形量计算与分析[J]. 机械工程师, 2018(10): 119-121. DOI:10.3969/j.issn.1002-2333.2018.10.041
[17]
邱江月, 夏金利. 探讨测量力对测试结果的影响[J]. 工业计量, 2010(2): 365-366.