CN104376211A - 一种坐标测量机测量不确定度评定方法 - Google Patents
一种坐标测量机测量不确定度评定方法 Download PDFInfo
- Publication number
- CN104376211A CN104376211A CN201410654418.4A CN201410654418A CN104376211A CN 104376211 A CN104376211 A CN 104376211A CN 201410654418 A CN201410654418 A CN 201410654418A CN 104376211 A CN104376211 A CN 104376211A
- Authority
- CN
- China
- Prior art keywords
- measurement
- uncertainty
- measuring machine
- coordinate measuring
- temperature
- 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.)
- Pending
Links
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种坐标测量机的测量不确定度标定方法,采取回避坐标测量机的各种误差源之间的复杂关系以及相互间未知的传递规律的方法,通过对坐标测量机进行测量***特性指标分析的方法,说明验证或获取该测量***相应的***特性指标值的方法和途径,并设计了获得坐标测量机的稳定性指标和复现性指标完整的标定方法以及详细的实验方案,考虑温度补偿中的实际温度以及标准光栅尺、被测工件的线膨胀系数对测量结果的影响程度,最终通过方和根法合成获得相对完整的坐标测量机的测量不确定度。
Description
技术领域
本发明涉及坐标测量机不确定度评定领域,具体是一种坐标测量机测量不确定度评定方法。
背景技术
坐标测量机(Coordinate Measuring Machine,CMM)是集机械学、光学、电学、计算机技术以及控制理论于一体的高效、多功能的精密型测量仪器,可以检测尺寸、角度、形状和位置等几何量。与常规的几何量测量仪器相比,坐标测量机所具备的优越性能体现在:测量对象多种多样,即通用性强;测量精度较高,数据准确可靠;可以进行复杂测量,自动化程度高,测量效率很高。目前的坐标测量技术已发展得较为成熟,坐标测量机广泛地应用于电子加工、汽车行业和航天科技等众多机械加工制造领域,成为产品几何量检测中最重要的手段。随着坐标测量技术的发展,以及其应用越来越广泛的同时,坐标测量机测量结果的不确定度评价技术也逐渐受到业内专家的关注和重视。国际标准化组织ISO发布的《测量不确定度表示指南》(Guide to the Expression of Uncertainty inMeasurement,GUM)指出,测量不确定度是评定测量***的重要指标,表示了测量结果的可靠程度,是测量结果中包含的一个重要参数,包含了测量不确定度的测量结果才是完整、可靠的、具有实用价值的。
坐标测量机是极其复杂的几何量测量***,与一般的单一对象测量的比较型测量仪器相比,想要进行坐标测量机面向任务的测量不确定度评价是很困难的。其一,坐标测量机是万能的几何量测量仪器,测量条件和方法也具有很多的选择余地,这就使得多种测量任务具有多种测量不确定度。此外,影响其测量结果的不确定度的误差来源诸多,且这些误差源与测量结果的传递关系一般很难确定,这就使得这些误差源对测量结果的影响难以量化表达。正是由于坐标测量机的误差复杂性和其多功能性,对于没有受过专业培训的测量者而言,想要评定坐标测量机的测量结果的不确定度并不现实。目前,不管是企业还是专业的计量技术机构在使用坐标测量机时,通常只能给出一个被测量的估计值,而没有该估计值的测量不确定度。科学、准确、便捷地评定坐标测量机的测量结果不确定度,成为精密测量与计量技术领域亟待解决的重要难题。
近年来,随着对坐标测量技术的研究日趋完善,评价三坐标测量机不确定度的理论和方法不断被提出,并逐渐形成了产品几何技术规范(Geometrical ProductSpecifications,GPS)ISO15530系列的坐标测量机面向任务的测量不确定度评价体系,该标准提供了为评估坐标测量机测量不确定度的名词、技术和指南,推荐了四种面向任务的评价不确定度的方法。国外的学者和计量研究机构对此非常重视,德国PTB、英国NPL都在按照ISO15530标准体系进行坐标测量机的测量不确定度评定的相关研究。
线性和偏移在坐标测量机的检定规范中有明确的定义和规定,一般以示值误差ΔE的形式来反映。示值误差一般来源于仪器自身,综合体现在坐标测量机的机械结构的探测误差、几何误差和软件误差等在测量结果中的影响程度。ISO10360标准规定示值误差在坐标测量机的验收检测或复检时进行评估,常规以坐标测量机对标准量块(或步距规)检测时得出的最大允许示值误差来进行量化;测量***的分辨力是指测量***识别并反映被测量最微小变化的能力。由于读数最末尾值的舍去与进位,将引入不确定度。且不确定度属于矩形分布(均匀分布)。测量***的重复性是指在测量条件不变的情况下,对被测参量进行连续多次的测量所得结果之间的一致性。坐标测量机的测量结果的重复性实质上反映的是示值的随机误差,其测量结果的重复性所引起的不确定度分量一般用以一组等精度测量结果的标准差定量反映。
测量复现性是指改变测量条件的情况下对同一被测参量进行测量,可改变的测量条件有测量原理、测量方法、测量者、参考标准以及测量时间等,只要有任一种或多种条件发生改变导致的分散性就称为复现性,简而言之,就是条件变化引起的测量结果的不一致性。在产品几何技术规范ISO15530系列中尽管考虑到了测量复现性是面向任务测量不确定度的一个分量,但是由于没有给出其具体的评价方法,因此操作起来难度非常大,从而导致测量结果的质量得不到保障。
坐标测量机的测量不确定度评定常见的方法是灵敏系数解析法。解析方法通常先计算原始误差源的数值,再考虑误差源到测量结果的传递关系,根据不同误差源的相互关系进行合成,最终得到测量结果的不确定度。然而坐标测量机的误差来源较多,且不同误差源之间的关系和传递规律较为复杂,因而,只是通过解析法去对坐标测量机的测量不确定度进行评定异常困难,实用性不强。
发明内容
本发明的目的是提供一种坐标测量机测量不确定度评定方法,以解决现有技术坐标测量机测量不确定度评定存在的问题。
为了达到上述目的,本发明所采用的技术方案为:
一种坐标测量机测量不确定度评定方法,其特征在于:包括坐标测量机测量稳定性标定及测量复现性标定两个部分,其中:
坐标测量机的测量稳定性标定过程如下:
(1)、由同一名测量人员对工件的同一被测参数进行单次测量;
(2)、保持测量环境不变,将坐标测量机保持工作状态,被测工件保持同样的位置及装夹,保持时间≥1小时后,由同一测量人员采用相同的测量方法再次对步骤(1)中的工件的同一被测参数进行单次测量;
(3)、保持相同的时间间隔,且保持条件不变依次进行多组上述的独立测量
(4)、各组测量的平均值按大小顺序排列成顺序统计量x(i),即有:
x(1)≤x(2)≤...≤x(a);
(5)、各组测量为等精度测量,服从均匀分布,选择其中的最大值x(a)减去最小值x(1),所得到的极差即反映了***测量的稳定性;
(6)、根据极差的分布函数可求得***稳定性所引起的不确定度分量为:
坐标测量机的测量复现性的标定过程如下:
(a)、由m名有一定的实际测量经验的且在一定程度上具备相近的测量水准测量人员,进行独立的b次连续测量,m≥5,b≥10;
(b)、每位测量人员测得一组数据后,间隔较短的时间(≤5min)后,另一位测量人员重新开始测量测量过程,不同测量人员之间的操作时间不宜过长,一般不超过5min,主要是考虑减少时间因素在测量复现性中影响程度,因为时间因素在***稳定性中已作考虑;
(c)、各测量人员因测量习惯、实际测量经验的不同其所选择的测量方案也不尽完全相同,但每个测量人员在进行单次重复测量时的测量条件不准随意改变,并应尽快完成单组的重复性测量;
(d)、计算各测量列的平均值,第i位测量人员的测量列的算术平均值记为
(e)、测量人员和测量方案的改变,容易导致上述的m组测量列的平均值之间出现异常值,发现异常值应及时剔除,采用t检验准则剔除可疑的异常值,具体步骤如下:
(e.1)、若认为为测量组中的异常值,将其剔除后计算m组测量列的平均值,即:
(e.2)、计算不包括在内的测量组的标准差:
(e.3)、根据测量次数m以及选定的显著度α,通过查询t分布的检验系数表得到K值;
(e.4)、若则认定为异常值,否则予以保留;
(f)、各组测量之间的测量条件稍有不同,即非等精度测量,因此不能直接用贝塞尔公式计算实验标准差,而是通过计算其合并样本标准差sp(x)来确定测量复现性引起的不确定度分量,过程如下:
(f.1)、分别计算出m组测量结果的各组实验标准差s(xi):
(f.2)、m组测量所包含的测量次数相同,合并样本标准差sp(x)为:
(f.3)、测量复现性是不同测量条件下所得测量值的实验标准差,因此复现性所引起的不确定度分量为:
所述的一种坐标测量机测量不确定度评定方法,其特征在于:运用测量***的特性指标对坐标测量机的测量不确定度进行评定,具体操作步骤如下:
(1)、进行测量***的稳定性标定实验,计算稳定性所引起的不确定度分量uST:
(2)、进行测量***的复现性标定实验,计算复现性误差所引起的不确定度分量uRD:
(3)、进行重复性实验,多次测量后通过贝塞尔公式计算单次测量的标准差后计算出***重复性引入的测量不确定度分量uRT:
(4)、测量***的线性和偏移特性所引起的不确定度分量uLB,通过示值误差ΔE来反映,取均匀分布,即测量***的线性和偏移特性所引起的不确定度分量uLB为:
(5)、测量***的分辨率特性所引起的不确定度分量uRS,分辨率特性所引起的不确定度分量uRS是由于读数的末尾数值的舍入引入的,定义坐标测量机的分辨率为e,取均匀分布,则测量***的分辨率特性所引起的不确定度分量uRS为:
(6)、依据坐标测量机的温度补偿的误差模型,计算环境温度所引起的不确定度分量uT,过程如下:
(6.1)、建立温度补偿的误差数学模型:
δx=x·(αs-αx)·△T+x·αs·△ts-x·αx·△tx,
式中,x表示被测参数的数值大小,αs表示坐标测量机光栅尺的线膨胀系数,αx表示工件的线膨胀系数,ΔT表示测量环境的实际环境温度与标准温度20℃的差值,△ts表示坐标测量机光栅尺的温度相对于实际环境温度的偏差,△tx表示被测工件的温度相对于实际环境温度的偏差,根据该误差数学模型,分析温度修正的各项不确定度分量;
(6.2)、将温度补偿的误差数学模型展开后得到:
δx=x·αs·△T-x·αx·△T+x·αs·△ts-x·αx·△tx,
令
δx=δ1-δ2+δ3-δ4,
其中,δ1=x·αs·△T;δ2=x·αx·△T;δ3=x·αs·△ts;δ4=x·αx·△tx;
针对定义的δ1、δ2、δ3、δ4四个分量分别进行测量不确定度评定;
(6.3)、δ1=x·αs·△T中包含了数学期望为0的变量ΔT,所以其不确定度u(δ1)的数学表达式为:
u(δ1)=x·αs·u(△T);
(6.4)、δ2=x·αx·△T中也包含了数学期望为0的变量ΔT,所以其不确定度u(δ2)的数学表达式为:
u(δ2)=x·αx·u(△T);
在步骤(6.3)和(6.4)中,ΔT为实验室环境温度与标准参考温度20℃的偏差,假设其服从反正弦分布(U形分布),则:
(6.5)、δ3=x·αs·△ts中包含了数学期望为0的变量△ts,所以其不确定度u(δ3)的数学表达式为:
u(δ3)=x·αs·u(△ts),
式中,△ts表示坐标测量机光栅尺的温度相对于实际环境温度的偏差,一般情况下坐标测量机光栅尺温度与实验室环境温度的偏差较小,即u(Δts)≈0℃,因此该分量可忽略不计;
(6.6)、δ4=x·αx·△tx中包含了数学期望为0的变量△tx,所以其不确定度u(δ4)的数学表达式为:
u(δ4)=x·αx·u(△tx),
式中,△tx表示被测工件的温度相对于实际环境温度的偏差,假设其服从反正弦分布(U形分布),则:
(6.7)、合成环境温度所引起的不确定度:
(7)、上述的***特性指标和环境温度的标准不确定度分量均要考虑,计算各分量数值后按照方和根法进行标准不确定度合成:
(8)、按置信概率P=95%,取包含因子k=2,得到扩展不确定度:U95=k·uC,k=2。
本发明回避了坐标测量机的各不同来源的误差源之间的相互关系和传递规律的复杂性以至无法描述这些误差因素在测量过程中引起的不确定度分量的问题,通过基于测量稳定性和复现性标定的基础上评定坐标测量机的测量不确定度,设计相关的实验方案和要求,从而获取相应的测量***的特性指标,并以坐标测量机的六个***特性指标入手进行测量不确定度评定,考虑环境温度的影响,获得准确可靠的坐标测量机测量结果的不确定度,最终得到完整的测量结果的不确定度。
附图说明
图1是测量***的重复性与复现性比较图。
图2是本发明实例的实现流程图。
具体实施方式
一种坐标测量机测量不确定度评定方法,包括坐标测量机测量稳定性标定及测量复现性标定两个部分,其中:
坐标测量机的测量稳定性标定过程如下:
(1)、由同一名测量人员对工件的同一被测参数进行单次测量;
(2)、保持测量环境不变,将坐标测量机保持工作状态,被测工件保持同样的位置及装夹,保持时间≥1小时后,由同一测量人员采用相同的测量方法再次对步骤(1)中的工件的同一被测参数进行单次测量;
(3)、保持相同的时间间隔,且保持条件不变依次进行多组上述的独立测量;
(4)各组测量的平均值按大小顺序排列成顺序统计量x(i),即有:
x(1)≤x(2)≤...≤x(a);
(5)、各组测量为等精度测量,服从均匀分布,选择其中的最大值x(a)减去最小值x(1),所得到的极差即反映了***测量的稳定性;
(6)、根据极差的分布函数可求得***稳定性所引起的不确定度分量为:
坐标测量机的测量复现性的标定过程如下:
(a)、由m名有一定的实际测量经验的且在一定程度上具备相近的测量水准测量人员,进行独立的b次连续测量,m≥5,b≥10;
(b)、每位测量人员测得一组数据后,间隔较短的时间(≤5min)后,另一位测量人员重新开始测量测量过程,不同测量人员之间的操作时间不宜过长,一般不超过5min,主要是考虑减少时间因素在测量复现性中影响程度,因为时间因素在***稳定性中已作考虑;
(c)、各测量人员因测量习惯、实际测量经验的不同其所选择的测量方案也不尽完全相同,但每个测量人员在进行单次重复测量时的测量条件不准随意改变,并应尽快完成单组的重复性测量;
(d)、计算各测量列的平均值,第i位测量人员的测量列的算术平均值记为
(e)、测量人员和测量方案的改变,容易导致上述的m组测量列的平均值之间出现异常值,发现异常值应及时剔除,采用t检验准则剔除可疑的异常值,具体步骤如下:
(e.1)、若认为为测量组中的异常值,将其剔除后计算m组测量列的平均值,即:
(e.2)、计算不包括在内的测量组的标准差:
(e.3)、根据测量次数m以及选定的显著度α,通过查询t分布的检验系数表得到K值;
(e.4)、若则认定为异常值,否则予以保留;
(f)、各组测量之间的测量条件稍有不同,即非等精度测量,因此不能直接用贝塞尔公式计算实验标准差,而是通过计算其合并样本标准差sp(x)来确定测量复现性引起的不确定度分量,过程如下:
(f.1)、分别计算出m组测量结果的各组实验标准差s(xi):
(f.2)、m组测量所包含的测量次数相同,合并样本标准差sp(x)为:
(f.3)、测量复现性是不同测量条件下所得测量值的实验标准差,因此复现性误差所引起的不确定度分量为:
运用测量***的特性指标对坐标测量机的测量不确定度进行评定,具体操作步骤如下:
(1)、进行测量***的稳定性标定实验,计算稳定性所引起的不确定度分量uST:
(2)、进行测量***的复现性标定实验,计算复现性误差所引起的不确定度分量uRD:
(3)、进行重复性实验,多次测量后通过贝塞尔公式计算单次测量的标准差后计算出***重复性引入的测量不确定度分量uRT:
(4)、测量***的线性和偏移特性所引起的不确定度分量uLB,通过示值误差ΔE来反映,取均匀分布,即测量***的线性和偏移特性所引起的不确定度分量uLB为:
(5)、测量***的分辨率特性所引起的不确定度分量uRS,分辨率特性所引起的不确定度分量uRS是由于读数的末尾数值的舍入引入的,定义坐标测量机的分辨率为e,取均匀分布,则测量***的分辨率特性所引起的不确定度分量uRS为:
(6)、依据坐标测量机的温度补偿的误差模型,计算环境温度所引起的不确定度分量uT,过程如下:
(6.1)、建立温度补偿的误差数学模型:
δx=x·(αs-αx)·△T+x·αs·△ts-x·αx·△tx,
式中,x表示被测参数的数值大小,αs表示坐标测量机光栅尺的线膨胀系数,αx表示工件的线膨胀系数,ΔT表示测量环境的实际环境温度与标准温度20℃的差值,△ts表示坐标测量机光栅尺的温度相对于实际环境温度的偏差,△tx表示被测工件的温度相对于实际环境温度的偏差,根据该误差数学模型,分析温度修正的各项不确定度分量;
(6.2)、将温度补偿的误差数学模型展开后得到:
δx=x·αs·△T-x·αx·△T+x·αs·△ts-x·αx·△tx,
令
δx=δ1-δ2+δ3-δ4,
其中,δ1=x·αs·△T;δ2=x·αx·△T;δ3=x·αs·△ts;δ4=x·αx·△tx;
针对定义的δ1、δ2、δ3、δ4四个分量分别进行测量不确定度评定;
(6.3)、δ1=x·αs·△T中包含了数学期望为0的变量ΔT,所以其不确定度u(δ1)的数学表达式为:
u(δ1)=x·αs·u(△T);
(6.4)、δ2=x·αx·△T中也包含了数学期望为0的变量ΔT,所以其不确定度u(δ2)的数学表达式为:
u(δ2)=x·αx·u(△T);
在步骤(6.3)和(6.4)中,ΔT为实验室环境温度与标准参考温度20℃的偏差,假设其服从反正弦分布(U形分布),则:
(6.5)、δ3=x·αs·△ts中包含了数学期望为0的变量△ts,所以其不确定度u(δ3)的数学表达式为:
u(δ3)=x·αs·u(△ts),
式中,△ts表示坐标测量机光栅尺的温度相对于实际环境温度的偏差,一般情况下坐标测量机光栅尺温度与实验室环境温度的偏差较小,即u(Δts)≈0℃,因此该分量可忽略不计;
(6.6)、δ4=x·αx·△tx中包含了数学期望为0的变量△tx,所以其不确定度u(δ4)的数学表达式为:
u(δ4)=x·αx·u(△tx),
式中,△tx表示被测工件的温度相对于实际环境温度的偏差,假设其服从反正弦分布(U形分布),则:
(6.7)、合成环境温度所引起的不确定度:
(7)、上述的***特性指标和环境温度的标准不确定度分量均要考虑,计算各分量数值后按照方和根法进行标准不确定度合成:
(8)、按置信概率P=95%,取包含因子k=2,得到扩展不确定度:U95=k·uC。
本发明提供的基于测量稳定性和复现性标定的坐标测量机测量不确定度评定方法,该方法具体步骤如下:
1)对坐标测量机的***稳定性指标进行标定,得出uST;
2)对坐标测量机的测量复现性指标进行标定,得出uRD;
3)进行重复性实验,通过贝塞尔公式计算单次测量的标准差后计算出uRT;
4)计算通过示值误差ΔE来反映的测量***的线性和偏移特性所引起的不确定度分量uLB;
5)计算由于读数的末尾数值的舍入引入的测量***的分辨率特性所引起的不确定度分量uRS;
6)依据坐标测量机的温度补偿的误差模型,计算环境温度所引起的不确定度分量uT。
7)标准不确定度合成
8)扩展不确定度计算U=k·uC,k=2;
步骤1)中***稳定性所引起的不确定度分量uST表达的是坐标测量机在一段连续时间内对同一工件或标准量进行测量的变化或不一致性,以多组测量数据标准差定量地进行表示。
步骤2)中测量复现性所引起的不确定度分量uRD反映坐标测量机测量结果的分散性,以多组测量数据标准差定量地进行表示。测量复现性是测量***分析重要的量值特性指标,体现了测量中不确定性的一些基本特征,是测量不确定度的重要来源。测量复现性需要在符合测量规范的前提下主动改变某一个或多个会对测量结果产生影响、但影响程度难以量化分析的因素,这些因素在实际测量中一般很难控制或者并没有统一的标准,但确实对测量结果的不确定性造成影响。实质上,测量复现性是不同条件下,主要是不同的测量者、测量方法之间存在的***效应而导致的多次测量结果的变动性。图1表示反映了测量***复现性与重复性的区别。
步骤3)中测量重复性所引起的不确定度分量uRT反映坐标测量机示值的随机误差分量,以一组测量数据标准差定量地进行表示。
步骤4)中线性和偏移所引起的不确定度分量uLB以示值误差ΔE的形式来反映,示值误差以坐标测量机对标准量块检测得到的最大允许示值误差ΔE来进行量化。
步骤5)中***分辨率所引起的不确定度分量uRS是由于读数的末尾数值的舍入引入的。
步骤6)中环境温度所引起的不确定度分量uT由坐标测量机测量过程中采用的温度修正补偿误差引入的,其误差模型为:
δx=x·(αs-αx)·△T+x·αs·△ts-x·αx·△tx,
式中,x表示被测工件的尺寸,αs表示坐标测量机光栅尺的线膨胀系数,αx表示工件的线膨胀系数,ΔT表示测量环境的实际环境温度与标准温度20℃的差值,△ts表示坐标测量机光栅尺的温度相对于实际环境温度的偏差,△tx表示被测工件的温度相对于实际环境温度的偏差。根据该误差数学模型,分析温度修正的各项不确定度分量。将温度修正的数学模型展开后得到:
δx=x·αs·△T-x·αx·△T+x·αs·△ts-x·αx·△tx,
令
δx=δ1-δ2+δ3-δ4,
其中,δ1=x·αs·△T;δ2=x·αx·△T;δ3=x·αs·△ts;δ4=x·αx·△tx。
即环境温度所引起的不确定度分量uT包括四个不确定度分量u(δ1)、u(δ2)、u(δ3)、u(δ4)。
其中,δ1=x·αs·△T中包含了数学期望为0的变量ΔT,所以其不确定度的数学表达式为:
u2(δ1)=[x·△T·u(αs)]2+[x·αs·u(△T)]2=[x·αs·u(△T)]2,
即
u(δ1)=x·αs·u(△T),
式中,ΔT为实验室环境温度与标准参考温度20℃的偏差,假设其服从反正弦分布(U形分布),则
同理,δ2=x·αx·△T中也包含了数学期望为0的变量ΔT,所以其不确定度的数学表达式为:
u(δ2)=x·αx·u(△T),
δ3=x·αs·△ts中包含了数学期望为0的变量△ts,所以其不确定度的数学表达式为:
u(δ3)=x·αs·u(△ts),
式中,△ts表示坐标测量机光栅尺的温度相对于实际环境温度的偏差,一般情况下坐标测量机光栅尺温度与实验室环境温度的偏差较小,即u(Δts)≈0℃,因此该分量可忽略不计。
δ4=x·αx·△tx中包含了数学期望为0的变量△tx,所以其不确定度的数学表达式为:
u(δ4)=x·αx·u(△tx),
式中,△tx表示被测工件的温度相对于实际环境温度的偏差,假设其服从反正弦分布(U形分布),则
合成环境温度所引起的不确定度有:
具体实施例:
下面结合图2及实验室实际测量状况对本发明的具体实施方式作进一步的说明。需要说明的是,对于下列实施方式的说明用于帮助解释和理解本发明,并不构成对本发明的限定。具体实施方式包括如下具体步骤:
1)将被测工件提前放置于实验室恒温(不低于8h);开启三坐标测量机;记录实验室环境温度;清洁被测工件、工作台、校准球和测头;校准测头;装夹工件。
2)坐标测量机测量稳定性的标定
由同一名测量人员对被测参量进行一次测量。保持测量环境不变,将坐标测量机保持工作状态,工件保持同样的位置及装夹。若干小时后(一般不低于1h),由同一测量人员采用相同的测量方法再次对被测参数进行一次测量。相同的时间间隔后,保持条件不变进行a(a≥10)次独立测量。
将a次测量的数据按大小顺序排列成顺序统计量x(i),即有:
x(1)≤x(2)≤...≤x(a),
上述测量为等精度测量,服从均匀分布,选择其中的最大值x(a)减去最小值x(1)得到的极差即反映了***稳定性的影响程度。根据极差的分布函数可求得***稳定性所引起的不确定度分量为:
3)坐标测量机测量复现性的标定
由m名测量人员(m≥5)在三坐标测量机上进行独立的b(b≥10)次连续测量。上述不同测量人员均需有一定的实际测量经验,且在一定程度上具备相近的测量水准。
每位测量人员测得一组数据后,间隔较短的时间后,另一位测量人员应重新选择测头型号、校准测头、选择工件摆放位置和装夹工件、建立坐标系、决定测量方法。需要说明的是不同测量人员之间的操作时间不宜过长,一般不超过5min,主要是考虑减少时间因素在测量复现性中影响程度,因为时间因素在***稳定性中已作考虑。
各测量人员因测量习惯、实际测量经验的不同其所选择的测量方案也不尽完全相同,但每个测量人员在进行单次重复测量时的测量条件不准随意改变,并应尽快完成单组的重复性测量。
实验数据统计并记录于表1中。
表1 ***复现性的测量数据记录表
由于测量人员和测量方案的改变,即各组测量结果之间的测量条件已发生了改变,容易导致上述m组测量列的平均值之间出现异常值,发现异常值须及时剔除。本发明使用t检验准则剔除可疑的异常值,具体步骤如下:
若认为为测量组中的异常值,将其剔除后计算平均值,即:
计算不包括在内的测量组的标准差:
根据测量次数m以及选定的显著度α,通过表2查询t分布的经验系数K(m,α)来判定。若则认定为异常值。否则,予以保留。
表2 t分布的检验系数
由于各组测量之间的测量条件稍有不同,即非等精度测量,因此不能简单地直接用贝塞尔公式计算实验标准差,而是计算其合并样本标准差sp(x)。
首先,分别计算出m组测量结果的实验室标准差s(xi)。
再计算合并样本标准差sp(x):
测量复现性是不同测量条件下所得测量值的实验标准差,因此复现性误差所引起的不确定度分量为:
4)坐标测量机测量重复性的标定
建立合适的坐标系,相同的条件下,对工件的待测参数进行n次(n≥10)的重复连续测量,实验数据记录在表3中。
表3 重复性测量的数据记录表
通过贝塞尔公式计算单次测量的标准差,取表3中的N(N≥3)次测量数据的平均值作为测量结果的估计值,有
5)计算线性和偏移所引起的不确定度分量uLB。
线性和偏移所引起的不确定度分量uLB以示值误差ΔE的形式来反映,示值误差以三坐标测量机对标准量块检测得到的最大允许示值误差ΔE来进行量化,取均匀分布。
6)计算***分辨率所引起的不确定度分量uRS。
***分辨率所引起的不确定度分量uRS是由于读数的末尾数值的舍入引入的,定义三坐标测量机的分辨率为e,取均匀分布。
7)计算环境温度所引起的不确定度分量uT。
测量过程中采用坐标测量机的温度修正补偿方法,误差模型为:
δx=x·(αs-αx)·△T+x·αs·△ts-x·αx·△tx,
式中,x表示被测工件的尺寸,αs表示坐标测量机光栅尺的线膨胀系数,αx表示工件的线膨胀系数,ΔT表示测量环境的实际环境温度与标准温度20℃的差值,△ts表示坐标测量机光栅尺的温度相对于实际环境温度的偏差,△tx表示被测工件的温度相对于实际环境温度的偏差。根据该误差数学模型,分析温度修正的各项不确定度分量。
将该数学模型展开后得到:
δx=x·αs·△T-x·αx·△T+x·αs·△ts-x·αx·△tx,
将四个乘积项按照如下表达式定义:
δ1=x·αs·△T
δ2=x·αx·△T
,
δ3=x·αs·△ts
δ4=x·αx·△tx
则
δx=δ1-δ2+δ3-δ4。
7.1)不确定度分量u(δ1):
δ1=x·αs·△T中包含了数学期望为0的变量ΔT,所以其不确定度的数学表达式为:
u2(δ1)=[x·△T·u(αs)]2+[x·αs·u(△T)]2=[x·αs·u(△T)]2,
即
u(δ1)=x·αs·u(△T),
式中,ΔT为实验室环境温度与标准参考温度20℃的偏差,一般情况下最大偏差为±1℃,假设其服从反正弦分布(U形分布),则
7.2)不确定度分量u(δ2):
同理,δ2=x·αx·△T中也包含了数学期望为0的变量ΔT,所以其不确定度的数学表达式为:
u(δ2)=x·αx·u(△T),
式中u(ΔT)=0.707℃。
7.3)不确定度分量u(δ3):
δ3=x·αs·△ts中包含了数学期望为0的变量△ts,所以其不确定度的数学表达式为:
u(δ3)=x·αs·u(△ts),
式中,△ts表示坐标测量机光栅尺的温度相对于实际环境温度的偏差,一般情况下坐标测量机光栅尺温度与实验室环境温度的偏差较小,即u(Δts)≈0℃,因此该分量可忽略不计。
7.4)不确定度分量u(δ4):
δ4=x·αx·△tx中包含了数学期望为0的变量△tx,所以其不确定度的数学表达式为:
u(δ4)=x·αx·u(△tx),
式中,△tx表示被测工件的温度相对于实际环境温度的偏差,一般情况下该温度偏差通过预先将工件放置实验室恒温处理可控制在±0.2℃的范围内,假设其服从反正弦分布(U形分布),则
7.5)合成环境温度所引起的不确定度
将环境温度引入的各不确定度分量汇总到表4中:
表4 温度引入的不确定度分量预算表
表中,x表示被测参量的数值,αs表示坐标测量机光栅尺的线膨胀系数,αx表示工件的线膨胀系数。
8)测量不确定度合成与报告
根据上述的计算与分析,将各项不确定度分量的评价结果汇总,按表5式样汇总。
表5 坐标测量机的测量不确定度分量汇总表
假设各不确定度分量相互独立,则合成标准不确定度为:
按置信概率P=95%,取包含因子k=2,得扩展不确定度
U95=k·uC,k=2,
以上所述为本发明的较佳且常规的实例而已,但本发明不应只局限该实施例和附图所公开的内容。凡不脱离本发明所公开的精神下完成的等效或修改,都落入本发明的保护范围。
Claims (2)
1.一种坐标测量机测量不确定度评定方法,其特征在于:包括坐标测量机测量稳定性标定及测量复现性标定两个部分,其中:
坐标测量机的测量稳定性标定过程如下:
(1)、由同一名测量人员对工件的同一被测参数进行单次测量;
(2)、保持测量环境不变,将坐标测量机保持工作状态,被测工件保持同样的位置及装夹,保持时间≥1小时后,由同一测量人员采用相同的测量方法再次对步骤(1)中的工件的同一被测参数进行单次测量;
(3)、保持相同的时间间隔,且保持条件不变依次进行多组上述的独立测量;
(4)、各组测量的平均值按大小顺序排列成顺序统计量x(i),即有:
x(1)≤x(2)≤···≤x(a);
(5)、各组测量为等精度测量,服从均匀分布,选择其中的最大值x(a)减去最小值x(1),所得到的极差即反映了***测量的稳定性;
(6)、根据极差的分布函数可求得***稳定性所引起的不确定度分量为:
坐标测量机的测量复现性的标定过程如下:
(a)、由m名有一定的实际测量经验的且在一定程度上具备相近的测量水准测量人员,进行独立的b次连续测量,m≥5,b≥10;
(b)、每位测量人员测得一组数据后,间隔较短的时间(≤5min)后,另一位测量人员重新开始测量测量过程,不同测量人员之间的操作时间不宜过长,一般不超过5min,主要是考虑减少时间因素在测量复现性中影响程度,因为时间因素在***稳定性中已作考虑;
(c)、各测量人员因测量习惯、实际测量经验的不同其所选择的测量方案也不尽完全相同,但每个测量人员在进行单次重复测量时的测量条件不准随意改变,并应尽快完成单组的重复性测量;
(d)、计算各测量列的平均值,第i位测量人员的测量列的算术平均值记为
(e)、测量人员和测量方案的改变,容易导致上述的m组测量列的平均值之间出现异常值,发现异常值应及时剔除,采用t检验准则剔除可疑的异常值,具体步骤如下:
(e.1)、若认为为测量组中的异常值,将其剔除后计算m组测量列的平均值,即:
(e.2)、计算不包括在内的测量组的标准差:
(e.3)、根据测量次数m以及选定的显著度α,通过查询t分布的检验系数表得到K值;
(e.4)、若 则认定为异常值,否则予以保留;
(f)、各组测量之间的测量条件稍有不同,即非等精度测量,因此不能直接用贝塞尔公式计算实验标准差,而是通过计算其合并样本标准差sp(x)来确定测量复现性引起的不确定度分量,过程如下:
(f.1)、分别计算出m组测量结果的各组实验标准差s(xi):
(f.2)、m组测量所包含的测量次数相同,合并样本标准差sp(x)为:
(f.3)、测量复现性是不同测量条件下所得测量值的实验标准差,因此复现性所引起的不确定度分量为:
2.根据权利要求1所述的一种坐标测量机测量不确定度评定方法,其特征在于:运用测量***的特性指标对坐标测量机的测量不确定度进行评定,具体操作步骤如下:
(1)、进行测量***的稳定性标定实验,计算稳定性所引起的不确定度分量uST:
(2)、进行测量***的复现性标定实验,计算复现性误差所引起的不确定度分量uRD:
(3)、进行重复性实验,多次测量后通过贝塞尔公式计算单次测量的标准差后计算出***重复性引入的测量不确定度分量uRT:
(4)、测量***的线性和偏移特性所引起的不确定度分量uLB,通过示值误差ΔE来反映,取均匀分布,即测量***的线性和偏移特性所引起的不确定度分量uLB为:
(5)、测量***的分辨率特性所引起的不确定度分量uRS,分辨率特性所引起的不确定度分量uRS是由于读数的末尾数值的舍入引入的,定义坐标测量机的分辨率为e,取均匀分布,则测量***的分辨率特性所引起的不确定度分量uRS为:
(6)、依据坐标测量机的温度补偿的误差模型,计算环境温度所引起的不确定度分量uT,过程如下:
(6.1)、建立温度补偿的误差数学模型:
δx=x·(αs-αx)·ΔT+x·αs·Δts-x·αx·Δtx,
式中,x表示被测参数的数值大小,αs表示坐标测量机光栅尺的线膨胀系数,αx表示工件的线膨胀系数,ΔT表示测量环境的实际环境温度与标准温度20℃的差值,Δts表示坐标测量机光栅尺的温度相对于实际环境温度的偏差,Δtx表示被测工件的温度相对于实际环境温度的偏差,根据该误差数学模型,分析温度修正的各项不确定度分量;
(6.2)、将温度补偿的误差数学模型展开后得到:
δx=x·αs·ΔT-x·αx·ΔT+x·αs·Δts-x·αx·Δtx,
令
δx=δ1-δ2+δ3-δ4,
其中,δ1=x·αs·ΔT;δ2=x·αx·ΔT;δ3=x·αs·Δts;δ4=x·αx·Δtx;
针对定义的δ1、δ2、δ3、δ4四个分量分别进行测量不确定度评定;
(6.3)、δ1=x·αs·ΔT中包含了数学期望为0的变量ΔT,所以其不确定度u(δ1)的数学表达式为:
u(δ1)=x·αs·u(ΔT);
(6.4)、δ2=x·αx·ΔT中也包含了数学期望为0的变量ΔT,所以其不确定度u(δ2)的数学表达式为:
u(δ2)=x·αx·u(ΔT);
在步骤(6.3)和(6.4)中,ΔT为实验室环境温度与标准参考温度20℃的偏差,假设其服从反正弦分布(U形分布),则:
(6.5)、δ3=x·αs·Δts中包含了数学期望为0的变量Δts,所以其不确定度u(δ3)的数学表达式为:
u(δ3)=x·αs·u(Δts),
式中,Δts表示坐标测量机光栅尺的温度相对于实际环境温度的偏差,一般情况下坐标测量机光栅尺温度与实验室环境温度的偏差较小,即u(Δts)≈0℃,因此该分量可忽略不计;
(6.6)、δ4=x·αx·Δtx中包含了数学期望为0的变量Δtx,所以其不确定度u(δ4)的数学表达式为:
u(δ4)=x·αx·u(Δtx),
式中,Δtx表示被测工件的温度相对于实际环境温度的偏差,假设其服从反正弦分布(U形分布),则:
(6.7)、合成环境温度所引起的不确定度:
(7)、上述的***特性指标和环境温度的标准不确定度分量均要考虑,计算各分量数值后按照方和根法进行标准不确定度合成:
(8)、按置信概率P=95%,取包含因子k=2,得到扩展不确定度:
U95=k·uC,k=2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410654418.4A CN104376211A (zh) | 2014-11-14 | 2014-11-14 | 一种坐标测量机测量不确定度评定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410654418.4A CN104376211A (zh) | 2014-11-14 | 2014-11-14 | 一种坐标测量机测量不确定度评定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104376211A true CN104376211A (zh) | 2015-02-25 |
Family
ID=52555115
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410654418.4A Pending CN104376211A (zh) | 2014-11-14 | 2014-11-14 | 一种坐标测量机测量不确定度评定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104376211A (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104965211A (zh) * | 2015-04-20 | 2015-10-07 | 中国电子科技集团公司第二十研究所 | 一种地面误差不确定度生成方法和电子设备 |
CN105987676A (zh) * | 2016-03-17 | 2016-10-05 | 合肥工业大学 | 一种三坐标测量机机构误差引入的测量不确定度评定方法 |
CN106595960A (zh) * | 2016-12-29 | 2017-04-26 | 广东鸿图科技股份有限公司 | 一种差压式检漏设备准确性及稳定性的检验方法 |
CN108427837A (zh) * | 2018-03-06 | 2018-08-21 | 南昌航空大学 | 一种装配粗差判定方法及*** |
CN109855583A (zh) * | 2018-11-16 | 2019-06-07 | 中国航发西安动力控制科技有限公司 | 三坐标测量机测量不确定度的确定方法 |
CN109990712A (zh) * | 2019-03-28 | 2019-07-09 | 唐山钢铁集团有限责任公司 | 一种测宽仪在线校准方法 |
CN110967184A (zh) * | 2019-12-03 | 2020-04-07 | 合肥工业大学 | 基于振动信号分布特征识别的变速箱故障检测方法和*** |
CN111060049A (zh) * | 2019-12-31 | 2020-04-24 | 广西玉柴机器股份有限公司 | 一种验证三坐标测量机工件坐标系转换稳定性的方法 |
CN111582390A (zh) * | 2020-05-11 | 2020-08-25 | 江苏移动信息***集成有限公司 | 一种大数据预分析的方法及装置 |
CN111623922A (zh) * | 2020-06-24 | 2020-09-04 | 安东仪器仪表检测有限公司 | 压力变送器零位修正值不确定度评定方法 |
CN111678477A (zh) * | 2020-06-20 | 2020-09-18 | 贵阳航发精密铸造有限公司 | 一种涡轮工作叶片最终检验自动检测测量方法 |
CN112197723A (zh) * | 2020-09-29 | 2021-01-08 | 中国航发动力股份有限公司 | 一种坐标测量机的校核标准件及校核方法 |
CN113031514A (zh) * | 2021-03-18 | 2021-06-25 | 西南交通大学 | 一种基于计量学的R-test标定不确定度评定方法 |
CN113405452A (zh) * | 2021-05-10 | 2021-09-17 | 中航西安飞机工业集团股份有限公司 | 一种基于温度补偿的数字化工装飞机坐标系标定方法 |
CN113657819A (zh) * | 2021-10-20 | 2021-11-16 | 北京关键科技股份有限公司 | 基于武器装备质量属性的分析方法、***和介质 |
CN116276407A (zh) * | 2023-05-23 | 2023-06-23 | 北京特思迪半导体设备有限公司 | 减薄机晶圆厚度测量不确定度的确定方法及设备 |
-
2014
- 2014-11-14 CN CN201410654418.4A patent/CN104376211A/zh active Pending
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104965211A (zh) * | 2015-04-20 | 2015-10-07 | 中国电子科技集团公司第二十研究所 | 一种地面误差不确定度生成方法和电子设备 |
CN105987676A (zh) * | 2016-03-17 | 2016-10-05 | 合肥工业大学 | 一种三坐标测量机机构误差引入的测量不确定度评定方法 |
CN106595960A (zh) * | 2016-12-29 | 2017-04-26 | 广东鸿图科技股份有限公司 | 一种差压式检漏设备准确性及稳定性的检验方法 |
CN108427837A (zh) * | 2018-03-06 | 2018-08-21 | 南昌航空大学 | 一种装配粗差判定方法及*** |
CN109855583A (zh) * | 2018-11-16 | 2019-06-07 | 中国航发西安动力控制科技有限公司 | 三坐标测量机测量不确定度的确定方法 |
CN109990712A (zh) * | 2019-03-28 | 2019-07-09 | 唐山钢铁集团有限责任公司 | 一种测宽仪在线校准方法 |
CN110967184B (zh) * | 2019-12-03 | 2021-06-11 | 合肥工业大学 | 基于振动信号分布特征识别的变速箱故障检测方法和*** |
CN110967184A (zh) * | 2019-12-03 | 2020-04-07 | 合肥工业大学 | 基于振动信号分布特征识别的变速箱故障检测方法和*** |
CN111060049A (zh) * | 2019-12-31 | 2020-04-24 | 广西玉柴机器股份有限公司 | 一种验证三坐标测量机工件坐标系转换稳定性的方法 |
CN111582390A (zh) * | 2020-05-11 | 2020-08-25 | 江苏移动信息***集成有限公司 | 一种大数据预分析的方法及装置 |
CN111678477A (zh) * | 2020-06-20 | 2020-09-18 | 贵阳航发精密铸造有限公司 | 一种涡轮工作叶片最终检验自动检测测量方法 |
CN111623922A (zh) * | 2020-06-24 | 2020-09-04 | 安东仪器仪表检测有限公司 | 压力变送器零位修正值不确定度评定方法 |
CN112197723A (zh) * | 2020-09-29 | 2021-01-08 | 中国航发动力股份有限公司 | 一种坐标测量机的校核标准件及校核方法 |
CN112197723B (zh) * | 2020-09-29 | 2022-03-08 | 中国航发动力股份有限公司 | 一种坐标测量机的校核标准件及校核方法 |
CN113031514A (zh) * | 2021-03-18 | 2021-06-25 | 西南交通大学 | 一种基于计量学的R-test标定不确定度评定方法 |
CN113031514B (zh) * | 2021-03-18 | 2022-03-18 | 西南交通大学 | 一种基于计量学的R-test标定不确定度评定方法 |
CN113405452A (zh) * | 2021-05-10 | 2021-09-17 | 中航西安飞机工业集团股份有限公司 | 一种基于温度补偿的数字化工装飞机坐标系标定方法 |
CN113657819A (zh) * | 2021-10-20 | 2021-11-16 | 北京关键科技股份有限公司 | 基于武器装备质量属性的分析方法、***和介质 |
CN116276407A (zh) * | 2023-05-23 | 2023-06-23 | 北京特思迪半导体设备有限公司 | 减薄机晶圆厚度测量不确定度的确定方法及设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104376211A (zh) | 一种坐标测量机测量不确定度评定方法 | |
da Silva Hack et al. | Measurement uncertainty: Literature review and research trends | |
Raghavendra et al. | Engineering metrology and measurements | |
CN109990712B (zh) | 一种测宽仪在线校准方法 | |
Barbato et al. | Measurements for decision making | |
Hall et al. | An introduction to measurement uncertainty | |
Phillips et al. | A careful consideration of the calibration concept | |
Płowucha | Point-plane distance as model for uncertainty evaluation of coordinate measurement | |
Vasilevskyi et al. | Indicators of reproducibility and suitability for assessing the quality of production services | |
CN102269581B (zh) | 一种建筑构件加工精度测评方法 | |
CN108955837A (zh) | 一种质量流量计在线***误差的确定方法及其应用 | |
Salsbury et al. | Measurement uncertainty in the performance verification of indicating measuring instruments | |
Toteva et al. | Methods for selection of measuring instruments | |
CN103837484B (zh) | 一种用于消除光谱乘性随机误差的角度化多变量分析方法 | |
Lamichhane | Investigating Cost-Effective Calibration Methods for Handheld Measurement Devices: A Study on Vernier & Digital Calipers and Micrometers | |
Chunovkina et al. | On the quantitative expression of the accuracy of laboratory and industrial measurements | |
Meškuotienė et al. | Analysis of Fuel’s Amount Measurement Process and Its Reliability in Commercial Transaction | |
Çelebioğlu | Developing a computer program for evaluating uncertainty of some typical dimensional measuring and gauging devices | |
Kondrat’eva et al. | Recommendations on the choice of gas analysis equipment for systems of continuous monitoring and accounting of emissions from thermal power plants | |
Weißensee et al. | Risk of erroneously deciding conformity of measuring instruments | |
Masharipovich et al. | Normalization of characteristics of measurement means and classes of accuracy of measurement means | |
Diesing | New Lessons in Gage Management: Don't Overlook GR&R Testing | |
Ehrlich | Traceability considerations for the characterization and use of measuring systems | |
KARIMOV et al. | Technical science and Innovation | |
Sotirov et al. | STATIC METHOD FOR TAXIMETER VERIFICATION |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150225 |
|
RJ01 | Rejection of invention patent application after publication |