工业热电阻自动测量系统结果不确定度评定实例用于检定工业热电阻的自动测量系统,根据国家计量检定规程(JJG 229—1998)对不确定度分析时可以在0℃点,100℃点,现在A 级铂热电阻的测量为例. B1 冰点(0℃)B1.1 数学模型,方差与传播系数 根据规定,被检的R(0℃)植计算公式为R(0℃)=R i 0=⎪⎭⎫⎝⎛t dt dR t i = Ri0=⎪⎭⎫⎝⎛t dt dR ***0=⎪⎭⎫⎝⎛-t I dt dR R R ℃)( = R i - 0.00391R *(0℃)×)℃(0 0.00391R 0**℃)(R R I -= R i - 0.391×1.00**℃)(R R I -= R i - 0.39 []℃)(0**R R I - 式中: R(0℃)—被检热电阻在0℃的电阻值,Ω;R i —被检热电阻在0℃附近的测得值,Ω;R *(0℃)—标准器在0℃的电阻值,通常从实测的水三点值计算,Ω;R *i —标准器在0℃附近测的值,Ω。
上式两边除以被检热电阻在0℃的变化率并做全微分变为dt 0R =d ()391.0R i+d ⎪⎪⎭⎫⎝⎛⨯-2500399.0**0i R R=dt Ri +dt *0R +dt *iR将微小变量用不确定度来代替,合成后可得方差u 20R t =u 2iR t +u 2t *0R +u 2t *iR (B-2)此时灵敏系数C 1=1,C 2=1,C 3=–1。
B1.2 标准不确定分量的分析计算 B1.2.1 u 2iR t 项分量该项分量是检热电阻在0℃点温度t i 上测量值的不确定度。
包括有:a) 冰点器温场均匀性,不应大于0. 01℃,则半区间为0.005℃。
均匀分布,故 u 1.1=3005.0=0.003℃其估计的相对不确定度为20﹪,即自由度1.1ν=12,属B 类分量。
b) 由电测仪表测量被检热电阻所带入的分量。
本系统配用电测仪表多为6位数字表(K2000,HP34401等),在对100Ω左右测量时仍用100Ω挡,此时数字表准确度为100×106×读数+40×106×量程对工业铂热电阻Pt100来说,电测仪表带入的误差限(半宽)为 被δ=±(100×100×106-+100×40×106-=±0.014Ω化为温度:391.0014.0±=±0.036℃该误差分布从均匀分布,即 u 2.1=3036.0=0.021℃估计的相对不确定度为10﹪,即1.1ν=50,属B 累类分量。
c) 对被检做多次检定时的重复性本规范规定在校准自动测量系统时以一稳定的A 级被检铂热电阻作试样检3次,用极差考核其重复性,经实验最大差为4m Ω以内。
通道间偏差以阻值计时应不大于2m Ω,故连同通道间差异同向叠计在内时,重复性为6m Ω,约0.015℃,则 u 3.1=69.1015.0=0.009℃ 3.1ν=1.8,属A 类分量。
d) 被检热电阻自然效应的影响。
以半区间估计为2m Ω计约5mK 。
这种影响普遍存在,可视为两点分布,故u 4.1=15=5mK估计的相对不确定度为30﹪,即4.1ν=5,属B 类分量。
B1.2.2u *0R t项分量标准器以实测的R *tp 值进行计算,故该误差分量以二等铂电阻温度计检定规程中规定在检定过程中的R tp 的允许变化不超过5mK 来计入。
半区间为2.5mK 。
呈正态分布,即 u 2=58.25.2=0.97mK 自由度为2ν=∞,属B 类分量。
B1.2.3 u *0R t 项分量该分量是使用标准器测量过程中引入的,包括:a) 电测仪表测量标准器R *i 标时引入的标准铂电阻与被检热电阻用同一电测仪表,使用的是100Ω挡,此时数字表的准确度为 100×106读数+40×106×量程而标准器为了排除不同电测带入的系统误差和标准器因应力等引起R tp 值的变化,要求用同一电测测量其R tp 值和R t 值,以比值W t =tp tR R 来计算实际温度,此时如以电测仪器的准确度分别计算对R t 、R tp 的项贡献是不对的,这两项值相关。
推导如下:W t =tp tR R 全微分,得dW i =tp R 1dR i -2tpt R R dR tp (B-3) 对微小变量dR t 、R tp 的计算可以电测仪器的指标及dR t 、R tp 值计入式(B-3) 变为tW δ=tpR 1()66104010010100--⨯⨯+⨯⨯t R -()66104010010100--⨯⨯+⨯⨯tp R =W t -⨯⨯⨯+⨯⨯--661040100110100tp R W t 66104010010100--⨯⨯⨯-⨯⨯tpt R W =(1-W t )/R tp 61040100-⨯⨯⨯ (B-4)所以用比值W 计算时,点测仪器对标准电阻测值引入的极差(以标准铂电阻温度计的典型值计,也可套用某支具体标准铂电阻温度计值)为 ()æ0W δ=()61040100259999601.01-⨯⨯⨯-=6.38×109-换算成温度t δ=3910989.31038.6--⨯⨯=1.6610-⨯℃。
此值为半区间,服从均匀分布,即 u 1.3=3106.16℃-⨯=0.001mK估计的相对不确定度为10﹪,即1.3ν=50,属B 类分量。
b)标准铂电阻温度计的自然影响按二等标准铂电阻温度计检定规程,它的自然允许值不应大于4mK 。
按均匀分布,一半区间计入,即u 2.3=324⨯=1.2mK对此估计的相对不确定度为20﹪,即2.3ν=12,属B 类分量。
b) 标准铂电阻温度计计算温度的计算误差根据ITS-1900,其内插公式的非惟一性为1mK,可按两点分布对待且可靠程度很高故 u 3.3=121⨯=0.5mK B1.3 冰点的标准不确定度分量一览表表B-1 冰点的标准不确定度分量一览表B1.4 冰点的合成的标准不确定度,有效自由度与包含因子 u 2冰c =222222225.02.1001.097.059213+++++++=558.6u 冰c =122.150001.0558.1950211236.234444444+++++=40.46 取40=eff ν,t 95.0(40)=2.02B1.5 冰点的检得结果扩展不确定度u 95=2.02×23.6=47.7mK ,取0.05℃B2 100℃B2.1 数学模式、方差与灵敏系数根据规程,被检工业铂电阻R 100的计算公式为R (100℃)=R b-0.379×()()*100*100*dt dR R Rb -= R b-0.379()()*100**100*dt dR R W Rtpb-⨯(B-5)式中:R (100℃)—被检热电阻在100℃的测得值,Ω; R b —被检热电阻在100℃点附近的测得值,Ω; R *b —标准器的测得值,Ω; W *100—标准器证书上的W 100比值;R *tp —检定时用系统所配电测仪器实测得表准器的水三相点值,Ω。
实际上,标准器全为用比值计算温度,分度表也只有比值的变化率(dW/dt )*,式(B-5)还需变形,因为(dW/dt )*=(dW/dt)* R *tp ,所以R (100)= R b -0.379()()⎥⎥⎦⎤⎢⎢⎣⎡-⨯*100**100*100*dt dR R W dt dR R tp b= R b -0.379⎪⎪⎭⎫⎝⎛⨯-⨯⨯***100**00387.000387.0tp tp tp bR R W R R =R b -0.379⎪⎪⎭⎫ ⎝⎛-⨯00387.000387.0*100*W W b= R b -97.93W *b +97.93W *100 (B-6)式中:W *b —标准器在t b 点测得的阻值R *b 与同一电测实测得的水三相点值与R *tp 之比;其余同式(B-5)。
式(B-6)为本分析的数学模式。
对式(B-6)进行全微分得*100*10093.9793.97W dW dR dR b b +-=将变量以分量不确定度代之并合成可得方差计算公式,即2c u (100℃)=()()222*100*93.9793.97W W R u u u bb +-+ (B-7)B2.2 标准不确定度分量的分析计算 B2.2.1 R b 项本项为检定时该系统测量被检热电阻所引入的不确定度来源。
它包括有:a) 温场均匀性影响恒温水槽水平均匀性为0.01℃。
以标准器所在孔位为基点,以半区间等概率分布计入,则有 003.03201.01.1=⨯=u ℃相对不确定度为20﹪,即有121.1=ν,属B 类分量。
b) 温场的变化波动影响本规范规定了热电阻自校系统的恒温特性为0.04℃/10min ,以测量时间为5min 变化计为 0.02℃,以标准器值为基点在被检项上有半区间、均匀分布的不确定,即 006.0301.02..1==u ℃相对不确定队20﹪,122.1=ν,属B 类分量。
c) 电测仪器测被检电阻所带入的分量在100℃附近被检工业铂热电阻Pt100的阻值约为138.15Ω左右。
电测仪表用1K Ω挡带入的极限误差为()66101010001010051.138--⨯⨯+⨯⨯±=b δΩ±=02385.0 即约±0.02385/0.379=±0.063℃ 呈均匀分布,以半区间计入为 036.03063.03.1==u ℃相对不确定度20﹪,即502.1=ν,属B 类分量。
d) 测量的重复性本规范规定系统的重复行实验以三次等精度重复测量结果的最大差不大于12m Ω(单一通道),而各通道间的偏差允许不大于2m Ω,按同向叠加即使用任一通道的重复行为14m Ω即0.037℃。
根据JJF1059-1999,使用三次间极差法计算单次实验标准差时极差系数为1.64,自由度为1.8,服从t 分布,即 022.069.1037.04.1===单s u ℃ 8.14.1=ν,属A 分量。
e) 被检铂热电阻的自然影响在沸点由于温度较高,根据实验观察,外部热效应影响较少,但被检的阻值较大,内部热效应有一定程度的增大,故仍可按0℃点的估算。
5,55.15.1==νmK u 两点分布,属B 类分量。
B2.2.2 *b W 项该项是检定时对标准器在检定点b t 上的测量所包含不确定度来源。
它包括有:a) 电测仪器测量标准器值时带近来的。
对标准器测量仍为100Ω挡,由此分析式(B-4)得 ⎪⎪⎭⎫ ⎝⎛⨯⨯-±=-31041tp t w R W tδ 此处即⎪⎪⎭⎫ ⎝⎛⨯⨯-±=-31041tp t WR W tδ ⎪⎭⎫⎝⎛⨯⨯-±=-31042539265.1151028.6-⨯±= 以半区间计入,均匀分布处理,即 551.210627.331028.6--⨯=⨯u相对不确定度10﹪,501.2=ν,属B 类分量。