李 琴 王賀龍 趙燕滔 張 帷 王立寧

天津市產(chǎn)品質(zhì)量監(jiān)督檢測(cè)技術(shù)研究院，天津 300384

1 概述

CNAS-CL01：2006《檢測(cè)和校準(zhǔn)實(shí)驗(yàn)室能力認(rèn)可準(zhǔn)則》要求，檢測(cè)實(shí)驗(yàn)室應(yīng)具有并應(yīng)用測(cè)量不確定度評(píng)定程序[1]。實(shí)驗(yàn)室在進(jìn)行檢測(cè)參數(shù)不確定度評(píng)定時(shí)，計(jì)算工作量大且較繁瑣。使用Delphi 7[2]開(kāi)發(fā)測(cè)量不確定度分析軟件，可實(shí)現(xiàn)輸入檢測(cè)結(jié)果及不確定度各分量的值后直接報(bào)告測(cè)量結(jié)果的擴(kuò)展不確定度。

2 測(cè)量不確定度分析原理

測(cè)量不確定度是表征賦予被測(cè)量值分散性的非負(fù)參數(shù)[3]。進(jìn)行測(cè)量不確定度評(píng)定，首先分析檢測(cè)參數(shù)的不確定度來(lái)源，建立測(cè)量過(guò)程的模型。然后，進(jìn)行標(biāo)準(zhǔn)不確定度分量的評(píng)估和計(jì)算，包括A類不確定度分量的評(píng)估和B類不確定度分量的評(píng)估。A類評(píng)估對(duì)于常規(guī)檢測(cè)，測(cè)量系統(tǒng)穩(wěn)定、又在B類評(píng)估中考慮了儀器的漂移和環(huán)境條件的影響，可以采取預(yù)先評(píng)估的方式，再依據(jù)實(shí)際測(cè)量值得到A類不確定度。對(duì)于大部分檢測(cè)工作，不確定度各分量不相關(guān)[4]。依據(jù)各分量的不確定度值，可得到B類不確定度。最后計(jì)算合成標(biāo)準(zhǔn)不確定度、擴(kuò)展不確定度，報(bào)告結(jié)果。

3 軟件開(kāi)發(fā)過(guò)程

3.1 界面設(shè)計(jì)

為了方便使用人員錄入數(shù)據(jù)，在軟件的界面設(shè)計(jì)上采用一次可以錄入全部數(shù)據(jù)的方式，同時(shí)也方便使用人員核對(duì)已經(jīng)錄入的數(shù)據(jù)。例如輸入A類不確定度的實(shí)驗(yàn)值如下圖。

圖1 軟件界面設(shè)計(jì)

3.2 計(jì)算方法設(shè)計(jì)

根據(jù)項(xiàng)目目前輸入數(shù)據(jù)的需求：A類不確定度實(shí)驗(yàn)值至多20個(gè)數(shù)據(jù)，A類不確定度實(shí)測(cè)值至多10個(gè)，B類不確定度值至多20個(gè)，同時(shí)考慮軟件擴(kuò)展性采用靜態(tài)生成數(shù)組的方式，以下以A類不確定度的實(shí)驗(yàn)值說(shuō)明其擴(kuò)展性。

X: array[1..20]of TsuiEdit：//聲明20個(gè)A類不確定度的輸入控件

當(dāng)A類不確定度的值可能增加到25個(gè)時(shí)，首先在界面上新增加5個(gè)輸入控件，修改

X: array[1..20]of TsuiEdit：

為

X: array[1..25]of TsuiEdit：，

再將新增加的5個(gè)控件放入數(shù)組中如下：

X[21]:= suiEdit21：X[22]:= suiEdit22：X[23]:= suiEdit23：X[24]:= suiEdit24：X[25]:= suiEdit25：

而計(jì)算步驟不需要修改，重新編譯生成新的軟件即可完成。

軟件中另外一個(gè)可能會(huì)變化的需求是分布情況以及計(jì)算方法，為此特意編寫(xiě)了根據(jù)分布名稱求計(jì)算分布值的函數(shù)getM，在新增分布情況的計(jì)算方法是只需要修改getM函數(shù)即可，計(jì)算如下：

function getM（s: String）: Double：

begin

if s = '兩點(diǎn)分布' then Result := 1

else if s = '正態(tài)分布' then Result := 3

else if s = '矩形分布 ' then Result := Sqrt（3）

else if s = '三角形分布' then Result := Sqrt（6）

else if s = '反正弦分布' then Result := Sqrt（2）

else if s = '梯形分布' then Result := 2

else

begin

Result := -1：

ShowMessage（'未知的分布'）：

end：

end：

3.3 計(jì)算結(jié)果設(shè)計(jì)

為了方便用戶使用，在計(jì)算結(jié)果上采用使用用戶自行選擇計(jì)算精度的方式，如下圖。

用戶根據(jù)檢測(cè)要求，自行錄入精度，軟件再根據(jù)錄入的精度顯示用戶的計(jì)算結(jié)果，在一定程度上方便了用戶的使用。

4 使用軟件分析并報(bào)告測(cè)量不確定度

利用已開(kāi)發(fā)的測(cè)量不確定度分析軟件進(jìn)行測(cè)量不確定度分析，將預(yù)評(píng)估的測(cè)量值、實(shí)際測(cè)量值、各分量的不確定度值、各分量的分布情況輸入后，可獲得不確定度報(bào)告。

4.1 評(píng)定煙氣中二氧化硫測(cè)定的測(cè)量不確定度

檢測(cè)依據(jù)方法HJ/T57-2000[5]，使用KN-900SN煙氣分析儀，最大允許誤差D=±5%。事先對(duì)某鍋爐煙氣二氧化硫濃度測(cè)量進(jìn)行20次重復(fù)獨(dú)立測(cè)量，測(cè)量結(jié)果見(jiàn)下表。

表1 進(jìn)行n=20次獨(dú)立重復(fù)測(cè)量的測(cè)量值

HJ/T 57-2000規(guī)定對(duì)同一工況進(jìn)行連續(xù)m=3次讀數(shù)，測(cè)量結(jié)果分別為145mg/m3、141mg/m3、140mg/m3，取其平均值142mg/m3作為測(cè)量結(jié)果。

數(shù)值輸入如下。

圖3 煙氣中二氧化硫測(cè)定測(cè)量不確定度分析

結(jié)果報(bào)告如下。

圖4 煙氣中二氧化硫測(cè)定測(cè)量不確定度報(bào)告

4.2 評(píng)定金屬材料抗拉強(qiáng)度的測(cè)定測(cè)量不確定度

檢測(cè)依據(jù)方法GB/T228-2002，使用的儀器設(shè)備：電子拉伸試驗(yàn)機(jī)、0～25mm千分尺、0～150mm游標(biāo)卡尺。加工10塊標(biāo)稱厚度8. 00mm和標(biāo)稱寬度15. 00mm試樣。對(duì)這10個(gè)試樣進(jìn)行抗拉強(qiáng)度測(cè)量，測(cè)量結(jié)果見(jiàn)表2。

數(shù)值輸入如下。

圖5 金屬材料抗拉強(qiáng)度的測(cè)定測(cè)量不確定度分析

表2 金屬材料抗拉強(qiáng)度s重復(fù)性測(cè)量數(shù)據(jù)

結(jié)果報(bào)告如下。

圖6 金屬材料抗拉強(qiáng)度的測(cè)定測(cè)量不確定度報(bào)告

5 結(jié)語(yǔ)

上述分析軟件的建立，實(shí)現(xiàn)了按照新修訂的JJF1059.1-2012要求進(jìn)行A類不確定度的預(yù)評(píng)估，實(shí)現(xiàn)了分量不相關(guān)時(shí)B類不確定度的評(píng)定，實(shí)現(xiàn)了檢測(cè)實(shí)驗(yàn)室擴(kuò)展不確定度的評(píng)定與報(bào)告，大大簡(jiǎn)化了實(shí)驗(yàn)室測(cè)量不確定度的評(píng)估工作，提高了工作效率，切實(shí)可行。

[1]中國(guó)合格評(píng)定國(guó)家認(rèn)可委員會(huì). CNAS-CL01:2006，檢測(cè)和校準(zhǔn)實(shí)驗(yàn)室能力認(rèn)可準(zhǔn)則

[2]梁水. Delphi開(kāi)發(fā)典型模塊大全. 人民郵電出版社

[3]國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局. JJF1059. 1-2012. 測(cè)量不確定度評(píng)定與表示

[4]中國(guó)合格評(píng)定國(guó)家認(rèn)可委員會(huì). CNAS-GL05:2011. 測(cè)量不確定度要求的實(shí)施指南

[5]HJ/T57-2000 固定污染源排氣中二氧化硫的測(cè)定定電位電解法. 中國(guó)環(huán)境科學(xué)出版社

[6]GB/T228. 1-2010 金屬材料 拉伸試驗(yàn) 第1部分：室溫試驗(yàn)方法. 中國(guó)標(biāo)準(zhǔn)出版社