原位检测系统中触发式测头的误差分析与补偿

王立成暋黄信达暋丁暋汉华中科技大学数字制造装备与技术国家重点实验室,武汉,430074摘要:在分析了原位检测系统中触发式测头误差来源的基础上,提出了自动调整测头偏心的方法,并用标定球对测头的实际作用半径进行了校准.
同时,面向实际测量任务,提出了基于在线标定和双线性插值技术的逐点测头半径补偿方法,实验验证,所提出的半径补偿方法比商用软件提供的方法更准确,提高了测量结果的准确性.
关键词:原位检测;触发式测头;在线标定;半径补偿中图分类号:TH161;TB92暋暋暋暋暋暋DOI:10.
3969/j.
issn.
1004-132X.
2012.
15.
003ErrorAnalysisandCompensationforTouchTriggerProbeofOn-machineMeasurementSystemWangLicheng暋HuangXinda暋DingHanStateKeyLaboratoryofDigitalManufacturingEquipmentandTechnology,HuazhongUniversityofScienceandTechnology,Wuhan,430074Abstract:OnthebasisofanalysisontheerrorsfromtouchtriggerprobeofOMM,amethodtocorrecttheoffsetofprobewasgiven,andtheactualeffectiveradiusofprobewascalibratedbymeansofacalibrationball.
Orientedtopracticalmeasurementtask,anewproberadiuscompensationmeth灢odbasedonon-linecalibrationandbilinearinterpolationtechniqueswaspresented.
Theverificationexperimentsofradiuscompensationshowthattheaccuracybythismethod,comparedtothemethodprovidedbythecommercialsoftware,canbemuchimproved.
Keywords:on-machinemeasurement(OMM);touchtriggerprobe;on-linecalibration;radiuscompensation收稿日期:2012—03—06基金项目:国家重点基础研究发展计划(973计划)资助项目(2011CB706804);国家自然科学基金资助重点项目(50835004)0暋引言原位检测系统通常是指基于数控机床本体的多轴坐标测量系统,它在不重复装夹情况下可对被加工零件进行在机测量.
采用原位检测技术可以及时发现加工误差,并据此修正加工工艺参数,对零件进行补加工;同时可避免零件加工和测量在不同设备上进行而产生的重复定位误差和工时耗费,因此受到了研究者的重视,并得到越来越多的应用[1飊4].
该系统在大型工件、复杂曲面或多表面工件制造中的应用需求更为迫切.
由于原位检测时的环境较为恶劣,且影响测量精度的因素众多,所以原位检测系统的测量精度控制比坐标测量机的测量精度控制更困难.
目前,原位检测系统一般采用触发式测头,针对触发式测头精度的研究主要集中在坐标测量机检测环境下的测头预行程误差分析[5飊8]和测头半径补偿[9飊11]两个方面,但考虑到原位检测所处环境的复杂性,对其中测头精度的分析和误差补偿技术的研究有待进一步深入.
本文通过分析和实验,对原位检测系统中的触发式测头进行了校准,并提出了一种基于在线标定技术的测头半径补偿方法.
1暋原位检测系统的构成与工作流程本文涉及的零件原位检测实验平台由MIRKON五轴加工中心、ReinshawOMP40测头系统及DelCam/PowerInspect自动原位检测软件构成.
利用该平台,可以对已加工零件进行原位测量,获取被测对象的几何信息,在此基础上,自主开发用于测量误差补偿、形状误差评估和加工误差修正等工艺过程的算法,来解决复杂零件制造过程中出现的问题,实现复杂零件的加工-测量一体化.
系统的工作流程如图1所示.
2暋测量误差分析与测头校准原位检测的误差来源为硬件误差、测量规划误差、软件算法误差和环境引起的误差等几个方面.
其中,硬件中的测头的误差对测量结果准确性影响不容忽视,在测量前,需要对测头精度进行·4771·中国机械工程第23卷第15期2012年8月上半月图1暋基于原位检测的加工测量一体化方法流程图校准.
2.
1暋测头误差测头误差主要为测头安装偏心误差、测头预行程误差和测头半径补偿方法误差.
测头安装偏心误差指测头中心轴线和机床主轴中心轴线的偏差.
该误差通常通过千分表进行手工调整,也可先利用原位检测技术计算出测头偏心量,再通过软件进行补偿.
触发式测头本身存在各向异性,导致不同探测距离、触测速度、测杆长度以及测头直径对测头的触发时机产生影响,因而存在一个有规律变化的系统误差———测头预行程误差.
可以采用多次触发的方式,减小触测速度变化过大带来的影响,也可以通过标定实验的方法对预行程误差进行检测和补偿.
接触式测量获取的数据是测头中心点的坐标,其与被测对象表面接触点之间相差一个测头半径的距离,为了得到待测对象几何尺寸的真实数据,需要进行半径补偿.
但受到测头制造精度和被测对象表面特性的影响,半径补偿值不可避免地存在着误差,故本文提出了一种改进的逐点半径补偿方法.
2.
2暋坐标系定位原位检测的一个特点是待测零件的加工坐标系、测量坐标系和模型坐标在系理论上完全重合,省去了测量数据之间的坐标变换,方便进行误差计算.
但在使用标定球校准之前,需要调整机床坐标系原点,使其位于标定球的理想球心,保证测量规划的位置和实际测量的位置一致.
对准的流程如图2所示,表1给出了机床坐标系原点对标定球球心的对准结果.
图2暋坐标系定位流程图表1暋机床坐标系原点对标定球心自动对准实验数据毺m暋(x,y,z)初始原点坐标(-30.
1,28.
8,-8.
0)一次调整后的原点坐标(-2.
0,1.
0,11.
0)二次调整后的原点坐标(-2.
0,1.
0,0.
5)2.
3暋测头半径校准因为测头存在预行程误差,所以测头探针头标称直径不能作为测量时的实际作用直径使用,需要进行测头校准.
具体的步骤如下:(1)采用ISO推荐的25点测球法[12]对标称半径为Rcal的标定球进行自动测量,测量结果记录在表2中,测量得到的原始数据是对应测点位置上测头探针头球心的坐标.
暋表2暋实际测量得到的测头探针头中心点坐标mm暋点序(x,y,z)点序(x,y,z)1(10.
926,0.
022,-0.
014)14(5.
443,5.
494,7.
723)2(7.
702,7.
753,-0.
014)15(-0.
028,7.
748,7.
711)3(-0.
028,10.
930,-0.
015)16(-5.
486,5.
480,7.
703)4(-7.
734,7.
728,-0.
015)17(-7.
751,0.
022,7.
708)5(-10.
933,0.
022,-0.
015)18(-5.
497,-5.
446,7.
718)6(-7.
754,-7.
703,-0.
014)19(-0.
028,-7.
723,7.
731)7(-0.
028,-10.
928,-0.
014)20(5.
452,-5.
458,7.
736)8(7.
722,-7.
728,-0.
014)21(3.
845,-1.
582,10.
109)9(9.
318,3.
894,4.
176)22(1.
574,3.
890,10.
092)10(-3.
884,9.
330,4.
158)23(-3.
894,1.
624,10.
086)11(-9.
348,-3.
838,4.
164)24(-1.
632,-3.
848,10.
099)12(3.
847,-9.
333,4.
18)25(-0.
028,0.
022,10.
932)13(7.
719,0.
022,7.
733)暋暋(2)对这25点用非线性最小二乘拟合算法拟合球面[13],得到拟合球面的球心坐标为(-2.
8,1.
7,1.
1)毺m,拟合球面的半径Rfit=·5771·原位检测系统中触发式测头的误差分析与补偿———王立成暋黄信达暋丁暋汉10灡9286mm.
(3)得到的拟合半径是标定球半径和探针头实际作用半径的和,在假设标定球精度比测头探测精度高很多的前提下,可以得到测头实际作用的半径:Rprobe=Rfit-Rcal=10.
9286-15.
8758/2=2.
9907mm经过上述步骤完成测头半径校准后,将校准过的测头半径输入到数控系统,作为实际测量中的测头半径补偿值进行半径补偿.
以上的校准过程简单实用,是触发式坐标测量系统进行测头校准和半径补偿的必备过程.
但依据这种校准方法进行半径补偿,尚未考虑测头在实际测量点处由于触发方向、触发条件不同所导致的补偿半径的误差,而这种误差在原位测量中的效应会更趋明显,成为测量不确定度的一个重要来源.
因此,需要有一种面向测量任务的原位标定技术来进一步完善测头半径的补偿方法.
3暋测头误差补偿3.
1暋考虑测头偏心的误差修正由于探针或测头安装导致的测头偏心误差会使实际测量位置偏离基于CAD模型规划的测量点,从而发生测量点与名义值之间的误对准.
修正测头偏心误差就是根据实际测量点坐标寻找与之对应的名义点坐标,重新计算二者之间的误差.
测头偏心修正方法如图3所示.
图3暋测头偏心修正示意图图3中,PN是不考虑偏心时CAD模型上的规划测点,P曚O是存在偏心时得到的测头实际中心点,P曚T是完成半径补偿后的被测表面接触点,当存在测头偏心时,实际测点P曚T所对应的CAD模型上的名义点为P曚N点,通过计算|P曚TP曚N|,才能获得准确的测量结果.
通常情况下,偏心量和测头探针头半径相比很小,所以假设被测曲面名义点变化所引起的曲面曲率变化可以忽略,即规划的测点PN处法矢n和实测后名义点P曚N处的法矢相同,那么可以利用已知的PN、P曚O和n,通过内积公式求出P曚O和P曚N两点之间的距离:Dcorrect=|P曚OP曚N|=(P曚O-P曚N)·n又因为测头探针头补偿半径Rprobe已知,所以计算出实际测点P曚T与修正后的名义点P曚N之间的误差e(e=Rprobe-Dcorrect)后,可完成测头偏心误差的单点修正.
3.
2暋基于在线标定技术的半径补偿在基于标定球的测头校准基础上,对触发式测头原位检测的半径补偿方法改进.
改进方法的基本思想是:将被测零件和标准球同时放置在工作台上,零件测量前,利用标定技术建立测头半径补偿的误差列表,并使用双线性插值算法得到曲面任意位置的法矢信息及其对应的测头半径补偿值,再根据实际测量点的法矢信息,建立实际测量位置和误差列表之间的映射关系,通过搜索误差补偿列表,找到对应的测头半径补偿值,实现逐点半径补偿.
详细的算法流程如图4所示.
图4暋半径补偿流程图3.
3暋实验验证为了评价本文中提到的补偿算法,实验中将一个标定球面作为标准曲面进行了验证实验(图5).
首先在标定球的上半球面上采集181个点,并以采集点作为节点在半个球面上划分网格,根据每个节点处的测头探针头的实际作用半径建立半径补偿列表(表3).
然后通过双线性差值算法得到半个球面上不同法矢处的半径补偿值,图6给出了半个球面不同法矢方向上的测头半径补偿值的变化.
·6771·中国机械工程第23卷第15期2012年8月上半月图5暋原位检测在线标定与测量实验表3暋基于标准球面的网格节点半径补偿误差列表mm暋U向V向0曘15曘30曘45曘60曘75曘90曘0曘2.
99422.
99152.
98962.
98982.
98762.
98882.
988112曘2.
99422.
99282.
99212.
98952.
98942.
99002.
989724曘2.
99422.
99422.
99142.
99062.
99102.
98872.
993136曘2.
99422.
99302.
99202.
99172.
98982.
99182.
992348曘2.
99422.
99502.
99182.
99372.
98952.
99262.
992560曘2.
99422.
99292.
99202.
99182.
99122.
99192.
992672曘2.
99422.
99482.
99162.
99292.
99102.
99232.
992584曘2.
99422.
99332.
99162.
99142.
99002.
99152.
992896曘2.
99422.
99492.
99102.
99202.
99162.
98992.
9917108曘2.
99422.
99202.
99012.
99162.
98882.
98792.
9892120曘2.
99422.
99292.
99142.
99052.
99112.
99162.
9888132曘2.
99422.
99352.
99282.
99252.
98992.
99052.
9917144曘2.
99422.
99362.
99292.
99252.
99042.
99032.
9948156曘2.
99422.
99342.
99242.
99202.
99292.
99352.
9950168曘2.
99422.
99272.
99482.
99032.
99022.
98982.
9923180曘2.
99422.
99582.
99262.
99052.
99082.
99002.
9951192曘2.
99422.
99442.
99332.
99082.
98952.
98832.
9933204曘2.
99422.
99272.
99392.
98942.
99152.
98992.
9919216曘2.
99422.
99412.
99032.
98802.
98872.
98972.
9884228曘2.
99422.
99202.
99092.
98592.
98602.
98612.
9866240曘2.
99422.
99422.
99042.
98802.
98632.
98582.
9857252曘2.
99422.
99242.
99042.
98952.
98732.
98582.
9895264曘2.
99422.
99372.
99062.
98872.
98882.
98622.
9893276曘2.
99422.
99222.
99122.
98802.
99042.
98782.
9873288曘2.
99422.
99422.
99202.
98812.
99022.
99042.
9898300曘2.
99422.
99322.
98992.
98922.
99022.
98982.
9894312曘2.
99422.
99252.
98942.
98702.
98832.
98802.
9909324曘2.
99422.
99252.
98902.
99122.
98852.
98812.
9894336曘2.
99422.
99272.
98932.
98792.
98902.
98592.
9904348曘2.
99422.
99362.
99112.
98702.
98802.
98752.
9876暋暋在同一个环境下再测量一个直径经过标定的球面.
将本文的补偿方法和商用测量软件所带的补偿方法进行对比,结果见图7和表4.
从表4数据可见,所提出插值补偿方法得到的测量结果中,极差、均值误差和标准差都比商用软件所使用的定值补偿方法要小,从而验证了本文提出方法的有效性.
表4暋不同半径补偿方法的统计结果比较毺m暋补偿方法最大值最小值极差均值误差标准差软件方法5.
996-3.
2559.
2511.
0262.
622本文方法2.
921-3.
9686.
8890.
3661.
752图6暋测头半径补偿误差分布示意图图7暋不同半径补偿方法单点误差比较结果3.
4暋适用性分析本文提出的半径补偿是面向现场标定和具体测量任务的误差补偿方法,假设被测表面的法矢信息没有显著的变化,测头探针可以比较准确地获得规划时被测点处的法矢.
这在测量大型平坦曲面时,或者在所使用的原位检测测头存在较大的系统误差的情况下,该方法具有更明显的优势.
另外,与传统的测头探针半径补偿取一个固定的值相比,该方法更准确反应了实际测量对象的特点和实际测量的过程,可以减小测量的不确定度.
4暋结论(1)给出了原位检测工作流程,实现了复杂零件的原位检测.
(2)分析了原位检测系统中触发式测头的误差来源,对测头偏心和预行程导致的误差进行了校准.
(3)提出了基于在线标定技术的测头半径补偿方法,实验验证该方法比商用软件中采用的半径补偿方法精度更高.
参考文献:[1]暋MearsL,RothJT,DjurdjanovicD,etal.
QualityandInspectionofMachiningOperations:CMMInte灢grationtotheMachineTool[J].
JournalofManu灢facturingScienceandEngineering,2009,131(5):051006.
[2]暋ChoMW,KimGH,SeoT.
IntegratedMachiningErrorCompensationMethodUsingOMMDataand·7771·原位检测系统中触发式测头的误差分析与补偿———王立成暋黄信达暋丁暋汉ModifiedPNNAlgorithm[J].
InternationalJournalofMachineToolsandManufacture,2006,46(12/13):1417飊1427.
[3]暋ZelenyJ,JandaM.
AutomaticOn-machineMeas灢urementofComplexParts[J].
ModernMachineryScienceJournal,2009,3:92飊95.
[4]暋李斌,孙忠林,彭芳瑜,等.
叶片零件在线检测测点规划问题研究[J].
中国机械工程,2008,19(5):540飊542.
[5]暋WeckenmannA,EstlerT,PeggsG,etal.
ProbingSys灢temsinDimensionalMetrology[J].
CIRPAnnals-ManufacturingTechnology,2004,53(2):657飊684.
[6]暋EstlerWT,PhillipsSD,BorchardtB,etal.
ErrorCompensationforCMMTouchTriggerProbes[J].
PrecisionEngineering,1996,19(2/3):85飊97.
[7]暋WozniakA,DoboszM.
InfluenceofMeasuredOb灢jectsParametersonCMMTouchTriggerProbeAc灢curacyofProbing[J].
PrecisionEngineering,2005,29(3):290飊297.
[8]暋刘丽冰,刘又午,赵小松,等.
加工中心在线检测误差补偿技术研究[J].
中国机械工程,2000,11(3):313飊316.
[9]暋XiongZhenhua,LiZexiang.
ProbeRadiusCompen灢sationofWorkpieceLocalization[J].
JournalofManufacturingScienceandEngineering,2003,125(1):100飊104.
[10]暋SrinivasanS,KovvurY,AnandS.
ProbeRadiusCompensationforOn-machineMeasurementofSculpturedSurfaces[C]//ProceedingsofASMEInternationalMechanicalEngineeringCongressandExposition.
Anaheim,CA,USA,2004:913飊920.
[11]暋ErkanT,MayerR,WozniakA.
SurfaceProbingSimulatorfortheEvaluationofCMMProbeRadi灢usCorrectionSoftware[J].
InternationalJournalofAdvancedManufacturingTechnology,2011,55(1/4):307飊315.
[12]暋海克斯康测量技术(青岛)有限公司.
实用坐标测量技术[M].
北京:化学工业出版社,2007.
[13]暋ShakarjiCM.
Least-squaresFittingAlgorithmsoftheNISTAlgorithmTestingSystem[J].
Jour灢nalofResearchoftheNationalInstituteofStand灢ardsandTechnology,1998,103(6):633飊641.
(编辑暋张暋洋)作者简介:王立成,男,1976年生.
华中科技大学数字制造装备与技术国家重点实验室博士研究生.
主要研究方向为曲面加工技术、精密测量技术与应用.
黄信达,男,1991年生.
华中科技大学数字制造装备与技术国家重点实验室博士研究生.
丁暋汉,男,1963年生.
华中科技大学数字制造装备与技术国家重点实验室教授、博士研究生导师.
(上接第1773页)[2]暋郑岗,赵亮.
金属厚度的脉冲涡流无损检测研究[J].
传感器与微系统,2006,25(4):35飊40.
[3]暋SophianAli,TianGuiyun,TaylorD,etal.
AFea灢tureExtractionTechniqueBasedonPrincipalCom灢ponentAnalysisforPulsedEddyCurrentNDT[J].
NDT&EInternational,2003,36(1):37飊41.
[4]暋TianGY,SophianA,TaylorD,etal.
MultipleSen灢sorsonPulsedEddy-currentDetectionfor3-DSubsurfaceCrackAssessment[J].
SensorsJournal,2005,5(1):90飊96.
[5]暋LebrunB,JayetY,BabouxJC.
PulsedEddyCurrentApplicationtotheDetectionofDeepCracks[J].
Mate灢rialsEvaluation,1995,53(11):1296飊1300.
[6]暋BieberJA,ShaligramSK,RoseJH,etal.
Time-gatingofPulsedEddyCurrentSignalsforDefectCharacterizationandDiscriminationinAircraftLap-joints[J].
ReviewofProgressinQNDE,1997,16B:1915飊1921.
[7]暋SmithRA,HugoGR.
DeepCorrosionandCrackDetectioninAgingAircraftUsingTransientEddyCurrentNDE[J].
ReviewofProgressinQNDE,1999,19:1401飊1408.
[8]暋杨宾峰,罗飞路,张玉华,等.
飞机多层结构中裂纹的定量检测及分类识别[J].
机械工程学报,2006,42(2):63飊67.
[9]暋HeYunze,PanMengchun,LuoFeilu,etal.
PulsedEddyCurrentImagingandFrequencySpectrumAnalysisforHiddenDefectNondestructiveTestingandEvaluation[J].
NDT&EInternational,2011,44(4):344飊352.
[10]暋周德强,田贵云,尤丽华,等.
基于频谱分析的脉冲涡流缺陷检测研究[J].
仪器仪表学报,2011,32(9):1948飊1953.
[11]暋ParkDG,AnganiCS,KimGD.
EvaluationofPulsedEddyCurrentResponseandDetectionoftheThicknessVariationintheStainlessSteel[J].
IEEETransactionsonMagnetics,2009,45(10):3893飊3896.
[12]暋PlotnikovYA,NathSC,RoseCW.
DefectChar灢acterizationinMulti-layeredConductiveCompo灢nentswithPulsedEddyCurrent[J].
ReviewofProgressinQNDE,2002,21A:1976飊1983.
(编辑暋张暋洋)作者简介:周德强,男,1979年生.
江南大学机械工程学院讲师、博士,无锡国盛精密模具有限公司工程师.
主要研究方向为电磁无损检测与评估.
发表论文10余篇.
李暋勇(通讯作者),男,1978年生.
西安交通大学航天航空学院副教授、博士.
张秋菊,女,1963年生.
江南大学机械工程学院教授、博士研究生导师.
尤丽华,女,1955年生.
江南大学机械工程学院副教授.
·8771·中国机械工程第23卷第15期201