本文通过ITS-90国际温标宣贯手册中对标准铂电阻温度计在0-419.527℃的参考函数及偏差函数的定义,采用Excel编写计算程序的方法实现了标准销电阻温度计分度值的计算。并通过与国际温标宣贯手册中编号为NO.80136标准铂电阻温度计分度值的数值比对来验证计算结果。
标准铂电阻温度计是根据金属钼的电阻随温度的变化而变化的规律来测量温度的,主要由感温元件、引线、保护管等组成。它属于精密计量器具,具有测量精度高、稳定性好的性能,现用于我国0-961℃的国家基准标准传递,由于其独特的计量特性,又作为温标复现、量值传递和温度测量的重要仪器,在当前各温度检定气程及校准规范中都明确将标准钻电阻温度计作为标准器。
通常所讲的标准铂电阻温度计主要是指测量范围为0-419.527℃的标准铂电阻温度计。作为计量检定机构最主要的温度量传标准,标准铀电阻温度计在送上级部门检定后,证书中只给出了a8、b8、Rtp、WZn、WSn、W(100)、WGa以及在水三相点温度时1mA电流引起的自热效应等参数的数据,由于缺少其它温度点的分度值,使得该仪器的使用范围受到限制。为了获得标准钼电阻温度计任意温度点的分度值,昌晖仪表建议可根据标准铂电阻温度计的参考函数和偏差函数,运用Excel“循环引用”功能来完成“逐次逼近计算”功能,实现标准铂电阻温度计任意分度值的计算。并通过与国际温标宣贯手册中编号为NO.80136#标准铂电阻温度计分度值的数值比对来验证计算结果的准确可靠。
1、标准铂电阻温度计参考函数及偏差函数
标准铂电阻温度计的电阻比W和温度t之间的函数关系Wr(t)作为参考函数,其定义公式为:
-------------------------------------------------------- 公式①
公式①中,Wr(t)为参考函数;Ci为固定的常数;t为温度值,℃。
其他任意一支标准钻电阻温度计的电阻比W(t)与参考函数Wr(t)的差值即为偏差函数△W(t),其定义公式为:
----------------------公式②
公式②中,W(t)为标准铂电阻温度计的电阻比;Wr(t)为参考函数;a8、b8为标准铂电阻温度计的分度常数(由锡和锌凝固点的电阻比计算得到)。
由公式②可以得到:---------------公式③
分度值的计算方法如下:将某个温度点对应的参考函数值Wr(t)带入公式③的右边,得到W1(t),再将W1(t)带入公式③中,得到W2(t):
-----------------------------公式④
如果|W2(t)-W1(t)|<1×10-8,则可认为W2(t)即是该温度点的分度值。如果上述不等式不能满足,则继续用W2(t)代替W1(t),求得W3(t),并重复以上过程直至满足上述不等式为止。
2、Excel函数的循环引用
循环引用就是公式中直接或间接地引用了公式本身所在的单元格。在计算标准铂电阻温度计的电阻比W(t)时,运用Excel“循环引用”功能很容易实现“逐次逼近计算”功能。具体操作如下:
①打开的Excel工作薄,在单元格内输入计算公式,单击【工具】中的【选项】命令。以0℃为例,在“分度表”工作薄的B1单元格内写入公式:B1=参考函数!B1+$C$1(B1-1)+$D$1(B1-1)2----------------------------公式⑤;公式⑤中,“参考上述!B1”是指在“参考函数”工作簿的B1单元格内对应的0℃参考函数值,“$C$1”和“$D$1”分别代表a8、b8,这里用到了Excel函数的绝对值引用。
②点击【重新计算】选项,选中【迭代计算】复选框,设置“最多迭代次数”及“最大误差”。经过实际计算,发现最多迭代次数取3次已满足要求,最大误差中输入1×10-8(IST-1990中要求达到的最大迭代误差),点击【确认】按钮后,即在单元格中显示计算结果,如图1所示。
3、数据比对验证结果
通过上述方法,获得了标准铂电阻温度计在任意温度下的分度值。将计算结果与IST-1990国际温标宣贯手册中的示例进行数据验证,如表1所示。
在Excel计算中,将W(t)、dW(t)/dt的计算结果比IST-1990国际温标宣贯手册中多保留一位。从表1中可以看出,跟IST-1990国际温标宣贯手册的计算结果是一致的,这也证明了采用该计算方法的可靠性。同样,也可以采用该方法计算其他温度段的分度值。
Excel是最常见的办公软件之一,操作简单方便,其数学函数功能也很强大。用户只需了解标准铂电阻温度计的计算公式,就可以依照上述所说的方法自行进行分度值的计算。