Примеры расчета неопределенности измерений - Школа метрологии
Главная » Полезные статьи » Примеры расчета неопределенности измерений

Примеры расчета неопределенности измерений

Оценивание неопределенности измерений (количественных величин) – одна из важных задач, стоящих перед каждой лабораторией. Требование к оцениванию неопределенности измерений заложено в межгосударственном стандарте ГОСТ ISO/IEC 17025-2019, а также политике ILAC-G17:2002.

Чтобы получить предварительную информацию по теме Неопределенность измерений, посмотрите, пожалуйста, наше обучающее видео:

Международное метрологическое сообщество давно уже разработало и приняло основные принципы концепции неопределенности, закрепив их в серии международных документов JCGM (Joint Committee for Guides in Metrology), а также документах ISO/IEC Guide 98. Разработано много дополнительных руководств по различных подходам оценивания неопределенности измерений в конкретных областях испытаний/измерений (EA, EURACHEM, Nordtest, EUROLAB и т.д.).

Несмотря на хорошо разработанные и представленные в документах свободного доступа теоретические принципы и подходы оценивания, и даже наличия конкретных примеров, у специалистов лабораторий все еще остаются вопросы по практической реализации этих принципов и подходов для измерений, проводимых в их лаборатории.

С целью наглядного представления основных принципов концепции неопределенности измерений и подхода моделирования предлагаем Вашему вниманию решение нескольких несложных, но часто встречающейся в практике многих испытательных лабораторий, задач:

Пример 4. Оценивание неопределенности измерений жесткости воды комплексонометрическим методом

Измеряемая величина – жесткость воды.

Методика измерений – ГОСТ 31954-2012 «Вода питьевая. Методы определения жесткости» (п. 4).

Метод испытаний – комплексонометрический метод, основан на образовании комплексных соединений трилона Б с ионами щелочноземельных элементов. Определение показателя проводят титрованием пробы раствором трилона Б при рН = 10 в присутствии индикатора.

Объект испытаний – проба питьевой воды.

1 Функция измерений

Жесткость воды Жi, °Ж, для каждого параллельного определения согласно ГОСТ 31954-2012 (п. 4.6.1) вычисляется в соответствии со следующей функцией измерений

где Стр – концентрация раствора трилона Б, ммоль/дм3;

F – множитель разбавления исходной пробы воды при консервировании, F = 1;

K – коэффициент поправки к концентрации раствора трилона Б;

Vтр – объем раствора трилона Б, израсходованный на титрование, см3;

Vпр – объем пробы, взятой для анализа, см3.

Расхождение между результатами двух параллельных определений не должно превышать предела повторяемости r = 10 % для значений жесткости св. 0,4 °Ж согласно ГОСТ 31954 -2012 (п. 4.7, таблица 1).

Измеренное значение жесткости Ж, °Ж, пробы питьевой воды вычисляется как среднее арифметическое двух параллельных определений в соответствии со следующей функцией измерений

где Fr – поправочный коэффициент, учитывающий допускаемое расхождение между результатами параллельных определений жесткости воды, F1 = 1.

Измеренное значение жесткости округляется до первого десятичного знака.

2 Результаты испытаний

При измерении жесткости пробы питьевой воды были получены результаты испытаний, представленные в таблице 1.

Таблица 1 – Исходные данные

3 Анализ входных величин

3.1 Концентрация раствора трилона Б, Стр

Значение оценки величины молярнойконцентрации раствора трилона Б рассчитывается исходя из процедуры приготовления раствора и принимается равным
Стр = 25 ммоль/дм3. Величина рассматривается как постоянная.

Погрешности, связанные с установлением значения молярнойконцентрации раствора трилона Б учитываются через коэффициент поправки K (см. п. 3.3).

3.2 Множитель разбавления исходной пробы воды при консервировании, F

Значение оценки величины принимается равным 1, т.к. разбавление пробы воды не проводится.

Неопределенность величины F не оценивается.

3.3 Коэффициент поправки к концентрации раствора трилона Б, K

Значение оценки величины коэффициента поправки к концентрации раствора трилона Б Ki при единичном определении вычисляется по формуле

где V1 – объем раствора иона магния, см3;

V2 – объем раствора трилона Б, израсходованный на титрование, см3;

F2 – поправочный коэффициент, учитывающий неопределенность измерений концентрации 25 ммоль/дм3 раствора ионов магния, F2 = 1.

На основании данных, приведенных в таблице 1, значения коэффициента поправки к концентрации раствора трилона Б Ki при трех определениях составили: K1 = K2 =K3 = 1,0000.

Расхождение между результатами трех параллельных определений коэффициента поправки к концентрации раствора трилона Б не превышает допускаемого значения, равного 0,001 (согласно ГОСТ 25794.1-83 «Реактивы. Методы приготовления титрованных растворов для кислотно-основного титрования», п. 1.11).

Значение оценки величины K вычисляется как среднее арифметическое трех параллельных определений коэффициента поправки к концентрации раствора трилона Б по формуле и составляет

где F3 – поправочный коэффициент, учитывающий допускаемое расхождение между результатами повторных измерений коэффициента поправки, F3 = 1.

Стандартная неопределенность измерений коэффициента поправки u(K) определяется по закону распространения неопределенностей путем суммирования относительных стандартных неопределенностей всех влияющих величин, входящих в функции измерений (3) и (4), по формуле

3.3.1Объем раствора иона магния V1

Значение оценки величины составляет V1 = 10 см3 и определяется по шкале пипетки с одной меткой 2-го класса точности номинальной вместимостью 10 см3 по ГОСТ 29169-91 «Посуда лабораторная стеклянная. Пипетки с одной отметкой».

Стандартная неопределенность измерений u(V1), см3, состоит из следующих основных вкладов:

– стандартная неопределенность, обусловленная отклонением действительного объема пипетки от номинального, u(VΔ,1), см3;

– стандартная неопределенность, обусловленная влиянием температурных колебаний на изменение объема жидкости при проведении измерений, u(Vt,1), см3.

Стандартная неопределенность, обусловленная отклонением действительного объема пипетки от номинального, u(VΔ,1) рассчитывается на основании информации о пределах допускаемой погрешности номинальной вместимости используемой пипетки с одной отметкой 2-го класса точности ±ΔV1 = ±0,04 см3, установленных ГОСТ 29169-91 (п. 7.1, таблица 1) при температуре 20 °С, в предположении треугольного распределения погрешности в установленных пределах по формуле и составляет

Стандартная неопределенность, обусловленная влиянием температурных колебаний на изменение объема жидкости при проведении измерений, u(Vt,1) рассчитывается в предположении прямоугольного распределения исходя из возможных пределов изменения объема жидкости при условии колебания температуры в лаборатории при измерении в пределах (20 ± 2) °С (±Δt = ±2 °С) и коэффициента объемного расширения жидкости kV, ºС-1, (равного 0,00021 ºС-1 для воды) по формуле и составляет

Стандартная неопределенность измерений u(V1) вычисляется путем суммирования стандартных неопределенностей перечисленных выше вкладов по формуле и составляет

3.3.2 Поправочный коэффициент, учитывающий неопределенность измерений концентрации 25 ммоль/дм3 раствора ионов магния, F2

Значение оценки величины F2 принимается равным 1.

Стандартная неопределенность измерений величины u(F2) вычисляется на основании относительной стандартной неопределенностью измерений концентрации 25 ммоль/дм3 раствора ионов магния u(C), по формуле

где С – концентрация раствора ионов магния, С = 25 ммоль/дм3.

Раствор ионов магния с концентрацией 25 ммоль/дм3 готовится из стандарт-титра магний сернокислый 7-водный с(1/2MgSO4·7H2O) = 0,1 моль/дм3 (ТУ 2642-001-56278322-2008). Значение оценки величины C рассчитывается по формуле и составляет

где c1 – приписанная концентрация стандарт-титра сульфата магния, с1 = 0,1 моль/дм3 (рассматривается как постоянная величина);

FV  – поправочный коэффициент, учитывающий разведение стандарт-титра, FV = 1;

1/4 – коэффициент, учитывающий фактор эквивалентности стандарт-титра сульфата магния и его разведение;

1000 – коэффициент перевода моль в ммоль.

Относительная стандартная неопределенность измерений концентрации 25 ммоль/дм3 раствора ионов магния u(C)/C рассчитывается по формуле

Значение оценки величины FV принимается равным 1.

Стандартная неопределенность измерений поправочного коэффициента, учитывающего разведение стандарт-титра, FV, вычисляется на основании информации о разведении стандарт-титра при приготовлении 25 ммоль/дм3 раствора ионов магния: в мерную колбу вместимостью 2000 см3 выливают содержимое ампулы и доливают до метки дистиллированной водой, по формуле

где V3 – объем приготовленного раствора ионов магния с концентрацией 25 ммоль/дм3, см3.

Значение оценки величины составляет V3 = 2000 см3 и определяется по шкале мерной колбы номинальной вместимостью 2000 см3 2-го класса точности по ГОСТ 1770-74«Посуда мерная лабораторная стеклянная. Цилиндры, мензурки, колбы, пробирки. Общие технические условия».

Стандартная неопределенность измерений u(V3), см3, состоит из следующих основных вкладов:

– стандартная неопределенность, обусловленная отклонением действительного объема мерной колбы от номинального, u(VΔ,3), см3;

– стандартная неопределенность, обусловленная влиянием температурных колебаний на изменение объема жидкости при проведении измерений, u(Vt,3), см3.

Стандартная неопределенность, обусловленная отклонением действительного объема мерной колбы от номинального, u(VΔ,3) рассчитывается на основании информации о пределах допускаемой погрешности номинальной вместимости используемой мерной колбы 2-го класса точности ±ΔV3 = ±1,20 см3, установленных ГОСТ 1770-74 (п. 2.4, таблица 7) при температуре 20 °С, в предположении треугольного распределения погрешности в установленных пределах по формуле и составляет

Стандартная неопределенность, обусловленная влиянием температурных колебаний на изменение объема жидкости при проведении измерений, u(Vt,3) рассчитывается в предположении прямоугольного распределения исходя из возможных пределов изменения объема жидкости при условии колебания температуры в лаборатории при измерении в пределах (20 ± 2) °С (±Δt = ±2 °С) и коэффициента объемного расширения жидкости kV, ºС-1, (равного 0,00021 ºС-1 для воды) по формуле и составляет

Стандартная неопределенность измерений u(V3) вычисляется путем суммирования стандартных неопределенностей перечисленных выше вкладов по формуле и составляет

Таким образом, стандартная неопределенность измерений u(FV), рассчитанная по формуле (12), составляет

Относительная стандартная неопределенность измерений концентрации 25 ммоль/дм3 раствора ионов магния, рассчитанная по формуле 11), составляет

Стандартная неопределенность измерений величиныF2 вычисляется по формуле (9) и в соответствии с приведенными выше расчетами составляет

3.3.3 Объем раствора трилона Б, израсходованный на титрование, V2

Значения оценки величины V2 определяются по шкале бюретки номинальной вместимостью 25 см3 2-го класса точности по ГОСТ 29251-91 «Посуда лабораторная стеклянная. Бюретки. Часть 1. Общие требования». Полученные при проведении испытаний пробы питьевой воды значения объемов представлены в таблице 1.

Стандартная неопределенность измерений u(V2), см3, состоит из следующих основных вкладов:

– стандартная неопределенность, обусловленная погрешностью измерения объема сливаемой жидкости бюретки, u(VΔ,2), см3;

– стандартная неопределенность, обусловленная погрешностью при снятии показаний со шкалы бюретки, u(Vc,2), см3;

– стандартная неопределенность, обусловленная отличием температуры, при которой проводятся измерения, от температуры, при которой нормируется погрешность измерения объема сливаемой жидкости бюретки, u(Vt,2), см3.

Стандартная неопределенность, обусловленная погрешностью измерения объема сливаемой жидкости, u(VΔ,2) вычисляется на основании информации о пределах погрешности измерения сливаемой жидкости используемой бюретки ±ΔV2 = ±0,1 см3, установленных ГОСТ 29251-91 (п. 6.1, таблица 1), в предположении треугольного распределения погрешности в установленных пределах по формуле и составляет

Стандартная неопределенность, обусловленная погрешностью при снятии показаний со шкалы бюретки, u(Vc,2) рассчитывается на основании информации о цене наименьшего деления используемой бюретки dV2 = 0,1 см3, в предположении прямоугольного распределения вероятностей погрешности снятия показаний в пределах ±dV2/2 по формуле и составляет

Стандартная неопределенность, обусловленная отличием температуры, при которой проводятся измерения, от температуры, при которой нормируется погрешность измерения объема сливаемой жидкости бюретки, u(Vt,2) рассчитывается в предположении прямоугольного распределения исходя из возможных пределов изменения объема жидкости при условии колебания температуры в лаборатории при измерении в пределах (20 ± 2) °С (±Δt = ±2 °С) и коэффициента объемного расширения жидкости kV, ºС-1, (равного 0,00021 ºС-1 для воды) по формуле и составляет

Стандартная неопределенность измерений u(V2) вычисляется путем суммирования стандартных неопределенностей перечисленных выше вкладов по формуле и составляет

3.3.4 Поправочный коэффициент, учитывающий допускаемое расхождение между результатами повторных измерений коэффициента поправки,F3

Значение оценки величины F3 принимается равным 1.

Стандартная неопределенность измерений поправочного коэффициента, учитывающего допускаемое расхождение между результатами повторных измерений коэффициента поправки, u(F3) определяется на основании установленного ГОСТ 25794.1-83 (п. 1.11) допускаемого расхождения между результатами трех определений величины K, rk = 0,001, которое рассматривается как критическая разность согласно ГОСТ Р ИСО 5725-6-2002 «Точность (правильность и прецизионность) методов и результатов измерений. Часть 6. Использование значений точности на практике» (п. 5). С учетом того, что за измеренное значение величины K принимается среднее арифметическое трех параллельных определений, стандартная неопределенность u(F3) рассчитывается по формуле и составляет

Стандартная неопределенность измерений коэффициента поправки к концентрации раствора трилона Б u(K) рассчитывается по формуле (5) и составляет

3.4 Объем раствора трилона Б, израсходованный на титрование, Vтр

Значения оценки величины Vтр определяются по шкале бюретки номинальной вместимостью 25 см3 2-го класса точности по ГОСТ 29251-91 «Посуда лабораторная стеклянная. Бюретки. Часть 1. Общие требования». Полученные при проведении испытаний пробы питьевой воды значения объемов представлены в таблице 1.

Стандартная неопределенность измерений u(Vтр), см3, состоит из следующих основных вкладов:

– стандартная неопределенность, обусловленная погрешностью измерения объема сливаемой жидкости бюретки, u(VΔ), см3;

– стандартная неопределенность, обусловленная погрешностью при снятии показаний со шкалы бюретки, u(Vc), см3;

– стандартная неопределенность, обусловленная отличием температуры, при которой проводятся измерения, от температуры, при которой нормируется погрешность измерения объема сливаемой жидкости бюретки, u(Vt), см3.

Стандартная неопределенность, обусловленная погрешностью измерения объема сливаемой жидкости, u(VΔ) вычисляется на основании информации о пределах погрешности измерения сливаемой жидкости используемой бюретки ±ΔV = ±0,1 см3, установленных ГОСТ 29251-91 (п. 6.1, таблица 1), в предположении треугольного распределения погрешности в установленных пределах по формуле и составляет

Стандартная неопределенность, обусловленная погрешностью при снятии показаний со шкалы бюретки, u(Vc) рассчитывается на основании информации о цене наименьшего деления используемой бюретки dV = 0,1 см3 в предположении прямоугольного распределения вероятностей погрешности снятия показаний в пределах ±dV/2 по формуле и составляет

Стандартная неопределенность, обусловленная отличием температуры, при которой проводятся измерения, от температуры, при которой нормируется погрешность измерения объема сливаемой жидкости бюретки, u(Vt) рассчитывается в предположении прямоугольного распределения исходя из возможных пределов изменения объема жидкости при условии колебания температуры в лаборатории при измерении в пределах (20 ± 2) °С (±Δt = ±2 °С) и коэффициента объемного расширения жидкости kV, ºС-1, (равного 0,00021 ºС-1 для воды) по формуле и составляет

Стандартная неопределенность измерений u(Vтр) вычисляется путем суммирования стандартных неопределенностей перечисленных выше вкладов по формуле и составляет

3.5 Объем пробы воды, взятой для анализа, Vпр

Значение оценки величины составляет Vпр = 100 см3 для двух параллельных проб питьевой воды и определяется по шкале мерной колбы номинальной вместимостью 100 см3 2-го класса точности по ГОСТ 1770-74«Посуда мерная лабораторная стеклянная. Цилиндры, мензурки, колбы, пробирки. Общие технические условия».

Стандартная неопределенность измерений u(Vпр), см3, состоит из следующих основных вкладов:

– стандартная неопределенность, обусловленная отклонением действительного объема мерной колбы от номинального, u(VΔ,4), см3;

– стандартная неопределенность, обусловленная влиянием температурных колебаний на изменение объема жидкости при проведении измерений, u(Vt,4), см3.

Стандартная неопределенность, обусловленная отклонением действительного объема мерной колбы от номинального, u(VΔ,4) рассчитывается на основании информации о пределах допускаемой погрешности номинальной вместимости используемой мерной колбы 2-го класса точности ±ΔVпр = ±0,20 см3, установленных ГОСТ 1770-74 (п. 2.4, таблица 7) при температуре 20 °С, в предположении треугольного распределения погрешности в установленных пределах по формуле и составляет

Стандартная неопределенность, обусловленная влиянием температурных колебаний на изменение объема жидкости при проведении измерений, u(Vt,4) рассчитывается в предположении прямоугольного распределения исходя из возможных пределов изменения объема жидкости при условии колебания температуры в лаборатории при измерении в пределах (20 ± 2) °С (±Δt = ±2 °С) и коэффициента объемного расширения жидкости kV, ºС-1, (равного 0,00021 ºС-1 для воды) по формуле и составляет

Стандартная неопределенность измерений u(Vпр) вычисляется путем суммирования стандартных неопределенностей перечисленных выше вкладов по формуле и составляет

3.7 Поправочный коэффициент, учитывающий расхождение между результатами параллельных определений, Fr

Значение оценки величины Fr принимается равным 1.

Стандартная неопределенность измерений поправочного коэффициента, учитывающего расхождение между результатами параллельных определений, u(Fr) определяется на основании установленного ГОСТ 31954-2012 (п. 4.7, таблица 1) предела повторяемости r = 0,1·Ж (см. ГОСТ Р ИСО 5725-6-2002, п. 4) для диапазона измерений жесткости св. 0,4 °Ж. С учетом того, что за измеренное значение жесткости пробы питьевой воды принимается среднее арифметическое двух параллельных определений, стандартная неопределенность u(Fr) рассчитывается по формуле и составляет

4 Корреляция

Все величины, входящие в функции измерений (1), (2), (3) рассматриваются как некоррелированные.

5 Оценка измеряемой величины

Оценка измеряемой величины – жесткость пробы питьевой воды Ж, рассчитывается по формуле (2) на основании результатов двух параллельных определений и составляет

Жесткость воды для каждого параллельного определения Жi рассчитывается по формуле (1) и составляет

Расхождение между результатами параллельных определений составляет 0,40 °Ж, что не превышает допускаемое значение (предел повторяемости)  °Ж, установленное ГОСТ 31954-2012 (п. 4.7, таблица 1).

6 Суммарная стандартная неопределенность

Стандартная неопределенность измерений жесткости пробы питьевой воды Ж вычисляется по закону распространения неопределенностей путем суммирования квадратов относительных стандартных неопределенностей влияющих величин, входящих в функции измерений (1) и (2), по формуле

Примечание – Для оценивания суммарной стандартной неопределенности измерений u(Ж) используются значения оценок влияющих величин (измеренные значения), полученные в рамках первого параллельного определения жесткости воды.

Значения относительных стандартных неопределенностей влияющих величин (измеряемой величины) рассчитываются как отношение стандартной неопределенности влияющей величины (измеряемой величины) к значению оценки влияющей величины (измеряемой величины).

Значения процентных вкладов неопределенностей влияющих величин в суммарную стандартную неопределенность рассчитываются как умноженный на 100 % квадрат отношения относительной стандартной неопределенности влияющей величины к квадрату относительной суммарной стандартной неопределенности измерений.

В таблице 2 представлен бюджет неопределенности для измеряемой величины.

Таблица 2 – Бюджет неопределенности измерений жесткости пробы питьевой воды

7 Расширенная неопределенность

Расширенная неопределенность измерений жесткости пробы питьевой воды U(Ж), °Ж, получается умножением суммарной стандартной неопределенности u(Ж), °Ж, на коэффициент охвата k, значение которого выбирается равным 2 при уровне доверия приблизительно 95 % в предположении нормального распределения вероятностей измеряемой величины (в бюджете неопределенности доминирует вклад от величины Fr, распределенной по нормальному распределению)

8 Результат измерения

Результат измерения жесткости пробы питьевой воды представляется в виде:

«Жесткость пробы питьевой воды составила (10,3 ± 0,5) °Ж, где число, следующее за знаком ±, является численным значением расширенной неопределенности, которая получена умножением суммарной стандартной неопределенности на коэффициент охвата k = 2, основанный на предполагаемом нормальном распределении, и определяет интервал, соответствующий уровню доверия приблизительно 95 %».

Библиография

ГОСТ 1770-74 Посуда мерная лабораторная стеклянная. Цилиндры, мензурки, колбы, пробирки. Общие технические условия

ГОСТ Р ИСО 5725-6-2002 Точность (правильность и прецизионность) методов и результатов измерений. Часть 6. Использование значений точности на практике

ГОСТ 25794.1-83 Реактивы. Методы приготовления титрованных растворов для кислотно-основного титрования

ГОСТ 29169-91 Посуда лабораторная стеклянная. Пипетки с одной отметкой

ГОСТ 29251-91 Посуда лабораторная стеклянная. Бюретки. Часть 1. Общие требования

ГОСТ 31954-2012 Вода питьевая. Методы определения жесткости

ТУ 2642-001-56278322-2008

Пример 3. Оценивание неопределенности измерений условной вязкости лакокрасочных материалов

Измеряемая величина – условная вязкость.

Методика измерений – ГОСТ 8420-74 «Материалы лакокрасочные. Методы определения условной вязкости» (п. 6.1).

Метод испытаний – испытания проводят с применением вискозиметра типа ВЗ-4 методом непрерывного истечения испытуемого материала. Пробу испытуемого материала доводят до температуры (20,0 ± 0,2) ºС. Время истечения измеряется секундомером типа СОПпр-2б-2-010.

Объект испытаний – эмаль.

1 Функция измерений

Условная вязкость эмали Х, с, согласно ГОСТ 8420-74 (п. 7) рассчитывается в соответствии с функцией измерений

где

t – среднее арифметическое двух измерений времени истечения, с;

K – поправочный коэффициент вискозиметра;

Fс – поправка на погрешность секундомера, с;

Fr – поправка, учитывающая допускаемый разброс значений времени истечения, F = 0,0 с.

Измеренное значение условной вязкости эмали округляется до целого числа.

2 Результаты испытаний

При определении условной вязкости эмали с применением вискозиметра ВЗ-4 и секундомера СОПпр-2б-2-010 было получено два значения времени истечения

t1 = 106,4 с, t2 = 104,2 с.

Относительная разность между двумя измерениями времени истечения составляет 2,1 % и не превышает допускаемого значения 5 %, которое установлено ГОСТ 8420-2022 (п. 7).

3 Анализ входных величин

3.1 Среднее арифметическое двух измерений времени истечения,

Значение оценки величины  рассчитывается по двум измерениями времени истечения по формуле и составляет

Возможное рассеяние значений единичных измерений времени истечения учитывается через влияющую величину Fr.

3.2 Поправка на погрешность секундомера, Fс

Значение оценки величины Fс принимается равным 0,0 с.

Стандартная неопределенность измерений поправки на погрешность секундомера u(Fс), с, рассчитывается на основании информации о пределах допускаемой относительной погрешности при температуре (20 ± 5) °С, δ = ±(1,7·A/Т + B), где А – значение скачка секундной стрелки, равное 0,2 с, В – составляющая относительной погрешности, определяемая отклонением частоты спускового регулятора от номинального значения, равная 4,3·10-4, Т – измеряемый интервал времени, с, и цене деления шкалы d = 0,2 c секундомера СОПпр-2б-2-010, используемого для измерения времени. В предположении прямоугольного распределения вероятностей значений погрешностей в границах ±δ и ±d/2 стандартная неопределенность измерений u(Fс) рассчитывается по формуле и составляет

3.3 Поправочный коэффициент вискозиметра, K

Значение оценки величины K берется из свидетельства о поверке на вискозиметр ВЗ-4 и составляет 0,96.

Стандартная неопределенность измерений поправочного коэффициента вискозиметра u(K) рассчитывается на основании информации о пределах значений основной относительной погрешности измерения времени истечения градуировочной жидкости ±δв = ±3 %, указанных в руководстве по эксплуатации на вискозиметр ВЗ-4, с учетом прямо пропорциональной зависимости поправочного коэффициента вискозиметра от времени истечения контрольной (градуировочной) жидкости (см. ГОСТ 8420-2022, Приложение В, п. В.3). В предположении прямоугольного распределения вероятностей значений погрешности в границах ±δв стандартная неопределенность измерений u(K) рассчитывается по формуле и составляет

3.4 Поправка, учитывающая допускаемый разброс значений времени истечения, Fr

Значение оценки величины Fr принимается равным 0,00 с.

Стандартная неопределенность измерений поправки, учитывающей допускаемый разброс значений времени истечения,u(Fr) рассчитывается на основании установленного ГОСТ 8420-2022 (п. 7) допускаемого относительного расхождения между двумя измерениями времени истечения = 5 %, которое рассматривается как предел повторяемости согласно положений ГОСТ Р ИСО 5725-6-2002, и с учетом того, что за измеренное значение времени истечения лакокрасочного материала принимается среднее арифметическое двух последовательных измерений, стандартная неопределенность u(Fr) рассчитывается по формуле и составляет

4 Корреляции

Все величины, входящие в функцию измерений (1), рассматриваются как некоррелированные.

5 Оценка измеряемой величины

Оценка измеряемой величины – условная вязкость эмали, X, рассчитывается по формуле (1) и составляет

6 Суммарная стандартная неопределенность

Стандартная неопределенность измерений условной вязкости эмали вычисляется по закону распространения неопределенностей путем суммирования квадратов произведений стандартных неопределенностей влияющих величин, входящих в функцию измерений (1), на соответствующие коэффициенты чувствительности по формуле и составляет

где коэффициенты чувствительности влияющих величин рассчитываются как частные производные функции измерений (1) по этим величинам и составляют

Вклад в неопределенность от влияющей величины рассчитывается как произведение стандартной неопределенности этой величины на соответствующий коэффициент чувствительности.

Значения процентных вкладов неопределенностей влияющих величин в суммарную стандартную неопределенность рассчитываются как умноженный на 100 % квадрат отношения вклада в неопределенность от влияющей величины к суммарной стандартной неопределенности измерений.

В таблице 1 представлен бюджет неопределенности для измеряемой величины.

Таблица 1 – Бюджет неопределенности измерений условной вязкости эмали

7 Расширенная неопределенность

Расширенная неопределенность измерений условной вязкости эмали U(Х), с, получается умножением суммарной стандартной неопределенности u(Х), с, на коэффициент охвата k, значение которого выбирается равным 1,65 при уровне доверия 95 % в предположении прямоугольного распределения вероятностей измеряемой величины (в бюджете неопределенности доминирует вклад от величины K, распределенной по прямоугольному распределению)

8 Результат измерения

Результат измерения условной вязкости эмали представляют в виде:

«Условная вязкость эмали составила (101 ± 3) с, где число, следующее за знаком ±, является численным значением расширенной неопределенности, которая получена умножением суммарной стандартной неопределенности на коэффициент охвата k = 1,65, основанный на предполагаемом прямоугольном распределении, и определяет интервал, соответствующий уровню доверия 95 %».

Библиография

ГОСТ Р ИСО 5725-6-2002 Точность (правильность и прецизионность) методов и результатов измерений. Часть 6. Использование значений точности на практике

ГОСТ 8420-2022 Материалы лакокрасочные. Методы определения условной вязкости

Пример 2. Оценивание неопределенности измерений плотности грунта методом режущего кольца

Измеряемая величина – плотность грунта.

Методика измерений – ГОСТ 5180-2015 «Грунты. Методы лабораторного определения физических характеристик» (п. 9).

Метод испытаний – методом режущего кольца, основан на взвешивании кольца-пробоотборника с грунтом и вычислении плотности грунта как отношения массы грунта к внутреннему объему кольца-пробоотборника.

Объект испытаний – проба глинистого грунта.

1 Функция измерений

Плотность пробы глинистого грунта ρ, г/см3, согласно ГОСТ 5180-2015 (п. 9.4) вычисляется в соответствии со следующей функцией измерений

где

m0 – масса кольца, г;

m1 – масса грунта с кольцом и пластинками, г;

m2 – масса пластинок, г;

V – внутренний объём кольца, см3, который рассчитывается по формуле

где

D – внутренний диаметр кольца, мм;

h – высота кольца, мм;

1000– коэффициент перевода мм3 в см3;

Fr – поправка, учитывающая допускаемый разброс результатов измерений плотности грунта, F = 0,00 г/cм3.

Измеренное значение округляют с точностью до 0,01 г/см3.

2 Результаты испытаний

При измерении плотности пробы глинистого грунта были получены результаты испытаний, представленные в таблице 1.

Таблица 1 – Результаты, полученные при испытаниях

3 Анализ входных величин

3.1 Масса грунта с кольцом и пластинками, m1

Значение оценки величины m1 определяется как показание весов электронных неавтоматического действия Pioneer РА213С, полученное при взвешивании грунта с кольцом и пластиной, и составило 120,335 г.

Стандартная неопределенность измерений массы грунта с кольцом и пластинками u(m1), г, оценивается на основании информации о пределах допускаемой погрешности в эксплуатации ±Δm, г, и дискретности отсчета d, г, установленных в описании типа на используемые весы электронные неавтоматического действия Pioneer РА213С (±Δ = ±0,02 г для диапазона измерений св. 50 до 200 г, d = 0,001 г). В предположении прямоугольного распределения вероятностей значений погрешностей в границах ±Δm и ±d/2 соответственно стандартная неопределенность u(m1) рассчитывается по формуле

3.2 Масса кольца, m0

Значение оценки величины m0 определяется как показание весов электронных неавтоматического действия Pioneer РА213С, полученное при взвешивании кольца-пробоотборника, и составило 60,058 г.

Стандартная неопределенность измерений массы кольца u(m0), г, оценивается на основании информации о пределах допускаемой погрешности в эксплуатации ±Δm, г, и дискретности отсчета d, г, установленных в описании типа на используемые весы электронные неавтоматического действия Pioneer РА213С (±Δ = ±0,02 г для диапазона измерений св. 50 до 200 г, d = 0,001 г). В предположении прямоугольного распределения вероятностей значений погрешностей в границах ±Δm и ±d/2 соответственно стандартная неопределенность u(m0) рассчитывается по формуле

3.3 Масса пластинок, m2

Значение оценки величины m2 определяется как показание весов электронных неавтоматического действия Pioneer РА213С, полученное при взвешивании пластинок, и составило 32,123 г.

Стандартная неопределенность измерений массы пластинок u(m2), г, оценивается на основании информации о пределах допускаемой погрешности в эксплуатации ±Δm, г, и дискретности отсчета d, г, установленных в описании типа на используемые весы электронные неавтоматического действия Pioneer РА213С (±Δ = ±0,01 г для диапазона измерений от 0,02 до 50 г, d = 0,001 г). В предположении прямоугольного распределения вероятностей значений погрешностей в границах ±Δm и ±d/2 соответственно стандартная неопределенность u(m2) рассчитывается по формуле

3.4 Внутренний диаметр кольца, D

Значение оценки величины D рассчитывается как среднее арифметическое четырех измеренных значений диаметра кольца в двух сечениях по высоте кольца (вблизи краев кольца) и в двух взаимно перпендикулярных направлениях (в каждом сечении) по формуле

где Di i-ое измеренное с применением штангенциркуля ШЦ-II-250-0,1 значение диаметра кольца, мм.

Измеренные значения Di, мм, представлены в таблице 1.

Стандартная неопределенность измерений внутреннего диаметра кольца u(D), мм, оценивается на основании информации о пределах допускаемой абсолютной погрешности ±Δш = ±0,10 мм и цене деления шкалы с = 0,1 мм, установленных в ГОСТ 166-89 на используемый штангенциркуль ШЦ-II-250-0,1 и разбросе измеренных значений Di. В предположении прямоугольного распределения вероятностей значений погрешностей в границах ±Δш и ±с/2 соответственно стандартная неопределенность u(D) рассчитывается по формуле

где sD– выборочное стандартное отклонение измеренных значений Di, которое рассчитывается по формуле

3.5 Высота кольца, h

Значение оценки величины h рассчитывается как среднее арифметическое четырех измеренных значений высоты кольца в четырех равномерно расположенных друг от друга точках по диаметру кольца по формуле

где hi i-ое измеренное с применением штангенциркуля ШЦ-II-250-0,1 значение высоты кольца, мм.

Измеренные значения hi, мм, представлены в таблице 1.

Стандартная неопределенность измерений высоты кольца u(h), мм, оценивается на основании информации о пределах допускаемой абсолютной погрешности ±Δш = ±0,10 мм и цене деления шкалы с = 0,1 мм, установленных в ГОСТ 166-89 на используемый штангенциркуль ШЦ-II-250-0,1 и разбросе измеренных значений hi. В предположении прямоугольного распределения вероятностей значений погрешностей в границах ±Δш и ±с/2 соответственно стандартная неопределенность u(h) рассчитывается по формуле

где sh– выборочное стандартное отклонение измеренных значений hi, которое рассчитывается по формуле

3.6 Внутренний объём кольца, V

Значение оценки величины V рассчитывается по формуле (2) на основании измеренных значений диаметра и высота кольца (см. пп. 3.4 и 3.5 соответственно) и составляет

Рассчитанное значение внутреннего объема кольца округляется до 0,1 см3 согласно положений ГОСТ 5180-2015 (п. 9.2.3).

Поскольку функция измерений (2) представляет собой произведения и частные влияющих величин, то стандартная неопределенность измерений внутреннего объема кольца u(V) будет вычисляться по закону распространения неопределенностей путем суммирования произведений относительных стандартных неопределенностей влияющих величин, входящих в функцию измерений (2), на соответствующие относительные коэффициенты чувствительности по формуле

3.7 Поправка, учитывающая допускаемый разброс результатов измерений плотности грунта, Fr

Значение оценки величины Fr принимается равным 0,00 г/cм3.

Стандартная неопределенность измерений поправки, учитывающей допускаемый разброс результатов измерений плотности грунта,u(Fr) рассчитывается на основании установленного ГОСТ 5180-2015 [1] (приложение A) допускаемого расхождения между результатами параллельных определений r = 0,03 г/cм3для глинистых грунтов, которое рассматривается как предел повторяемости согласно положений ГОСТ Р ИСО 5725-6-2002, и с учетом однократного определения плотности грунта рассчитывается по формуле

4 Корреляции

Все входные величины рассматриваются как некоррелированные.

5 Оценка измеряемой величины

Оценка измеряемой величины – плотность грунта ρ, рассчитывается по формуле (1) и составляет

6 Суммарная стандартная неопределенность

Стандартная неопределенность измерений плотности грунта ρ определяется по закону распространения неопределенностей путем суммирования квадратов произведений стандартных неопределенностей влияющих величин, входящих в функцию измерений (1), на соответствующие коэффициенты чувствительности по формуле

где коэффициенты чувствительности рассчитываются как частные производные функции измерений (1) по входным величинам

Значения вкладов неопределенностей влияющих величин в суммарную стандартную неопределенность рассчитываются как произведение коэффициента чувствительности влияющей величины на ее стандартную неопределенность.

Значения процентных вкладов неопределенностей влияющих величин в суммарную стандартную неопределенность рассчитываются как умноженное на 100 % отношение квадрата вклада стандартной неопределенности влияющей величины к квадрату суммарной стандартной неопределенности:

В таблице 2 представлен бюджет неопределенности для измеряемой величины.

Таблица 2 – Бюджет неопределенности измерений плотности грунта

7 Расширенная неопределенность

Расширенная неопределенность измерений плотности грунта U(ρ), г/см3, определяется путем умножения суммарной стандартной неопределенности u(ρ), г/см3, на коэффициент охвата k, значение которого выбирается равным 2 при уровне доверия приблизительно 95 % в предположении нормального распределения вероятностей измеряемой величины (в бюджете неопределенности доминирует вклад от величины Fr, распределенной по нормальному распределению)

8 Результат измерения

Результат измерения плотности грунта представляется в виде:

«Плотность пробы глинистого грунта, определенная методом режущего кольца, составила (0,22 ± 0,02) г/см3, где число, следующее за знаком ±, является численным значением расширенной неопределенности, которая получена умножением суммарной стандартной неопределенности на коэффициент охвата k = 2, основанный на предполагаемом нормальном распределении, и определяет интервал, соответствующий уровню доверия приблизительно 95 %».

Библиография

ГОСТ 5180-2015 Грунты. Методы лабораторного определения физических характеристик

ГОСТ 166-89 (ИСО 3599-76) Штангенциркули. Технические условия

ГОСТ Р ИСО 5725-6-2002 Точность (правильность и прецизионность) методов и результатов измерений. Часть 6. Использование значений точности на практике

Пример 1. Оценивание неопределенности измерений влажности пнеополистирола

Измеряемая величина – влажность пнеополистирола.

Методика измерений – ГОСТ 17177-94 «Материалы и изделия строительные теплоизоляционные. Методы испытаний».

Метод испытаний – основан на определении массы образца до и после высушивания до постоянной массы; высушивание образцов до постоянной массы проводится при температуре (105 ± 5) ℃; образцы материала считают высушенными до постоянной массы, если потеря их массы после повторного высушивания в течение 0,5 ч не превышает 0,1 % согласно ГОСТ 17177-94 (п. 3.5); перед каждым взвешиванием образцы материалов охлаждают в эксикаторе над хлористым кальцием.

Объект испытаний – образцы пнеополистирола в количестве 9 шт.

1. Функция измерений

Влажность пнеополистирола Wi, %, при испытании i-го образца вычисляется в соответствии с ГОСТ 17177-74 (п. 8.5) согласно функции измерений

где m4 – масса образца до сушки, г;

m5 – масса образца после сушки, г.

За измеренное значение влажности пнеополистирола W, %, принимается среднее арифметическое девяти параллельных определений

где Fr – поправка, учитывающая рассеяние результатов параллельных определений, %.

Результат измерения влажности пнеополистирола округляется с точностью до 0,1 %.

2 Результаты испытаний

При измерении влажности девяти образцов пнеополистирола были получены данные, представленные в таблице 1.

Таблица 1 – Исходные данные

3 Анализ входных величин

3.1 Масса образца до сушки, m4

Значения оценки величины m4 определяются как показание весов электронных Adventurer RV-64 в граммах, полученные при взвешивании i-го образца до сушки, и представлены в таблице 1.

Стандартная неопределенность измерений массы образца до сушки u(m4), г, оценивается на основании информации о пределах допускаемой погрешности в эксплуатации ±Δ = ±0,001 г (в диапазоне измерений от 0,01 до 50 г) и дискретности отсчета d = 0,0001 г используемых для взвешивания весов электронных Adventurer RV-64, установленных в описании типа на указанные весы. В предположении прямоугольного распределения вероятностей значений погрешностей в границах ±Δ и ±d/2 соответственно стандартная неопределенность u(m4) рассчитывается по формуле

3.2 Масса образца после сушки, m5

Значения оценки величины m5 определяются как показание весов электронных Adventurer RV-64 в граммах, полученные при взвешивании i-го образца после сушки, и представлены в таблице 1.

Стандартная неопределенность измерений массы образца после высушивания u(m5), г, состоит из следующих основных вкладов: стандартной неопределенностью, обусловленной погрешностью и дискретностью используемых весов, u(mw), г, и стандартной неопределенностью, обусловленной допускаемым возможным отклонением массы образца после высушивания, когда образец считается высушенным до постоянной массы, u(mdis), г.

Стандартная неопределенность u(mw) оценивается аналогично стандартной неопределенности измерений u(m4) по формуле (3) и составляет 0,000578 г.

Стандартная неопределенность u(mdis), г, оценивается на основании установленного ГОСТ 17177-94 (п. 3.5) требования к моменту, когда образец считается высушенным до постоянной массы: потеря массы после повторного высушивания в течение 0,5 ч не должна превышать 0,1 % от массы образца после высушивания m5. Это значение можно рассматривать как 95-% доверительный интервал для разности двух оценок величины, распределенной по нормальному закону распределения (предел повторяемости по ГОСТ Р ИСО 5725-6). В предположении нормального распределения вероятностей стандартная неопределенность u(mdis) будет равна стандартному отклонению, рассчитанному на основании указанного интервала по формуле

Примечание – Для расчета стандартной неопределенность u(mdis) используется значение оценки величины m5, полученное в рамках первого параллельного определения влажности (первого образца).

Стандартная неопределенность измерений u(m5) рассчитывается по формуле

3.3 Поправка, учитывающая рассеяние результатов параллельных определений, Fr

Значение оценки величины Fr принимается равным 0,0 %.

Стандартная неопределенность измерений поправки, учитывающей рассеяние результатов параллельных определений, u(Fr), определяется в предположении нормального распределения по формуле

где среднее арифметическое рассчитывается по формуле

4 Корреляции

Все входные величины рассматриваются как некоррелированные.

5 Оценка измеряемой величины

Оценка измеряемой величины «влажность пенополистирола (W,%)» определяется по формуле (2) на основании результатов девяти параллельных определений (см. результаты вычисления по формуле (1), представленные в таблице 1)

6 Суммарная стандартная неопределенность

Стандартная неопределенность измерений влажности пенополистирола получается по закону распространения неопределенностей путем суммирования квадратов произведений стандартных неопределенностей всех влияющих величин, входящих в функции измерений (1) и (2), на соответствующие коэффициенты чувствительности по формуле

где коэффициенты чувствительности определяются по формулам

Примечание – Для расчета коэффициентов чувствительности используются значения оценок влияющих величин (измеренные значения), полученные в рамках первого параллельного определения влажности (первого образца).

Значения вкладов неопределенностей влияющих величин в суммарную стандартную неопределенность рассчитываются как произведение коэффициента чувствительности влияющей величины на её стандартную неопределенность.

Значения процентных вкладов неопределенностей влияющих величин в суммарную стандартную неопределенность рассчитываются как умноженное на 100 % отношение квадрата вклада неопределенности влияющей величины к квадрату суммарной стандартной неопределенности.

Бюджет неопределенности измерений влажности пенополистирола представлен в
таблице 2.

Таблица 2 – Бюджет неопределенности измерений влажности пенополистирола

7 Расширенная неопределенность

Расширенная неопределенность измерений влажности пенополистирола U(W) рассчитывается по формуле

где коэффициент охвата k выбран равным 2 в предположении нормального распределения вероятностей, поскольку число эффективных степеней свободы более 10. Число эффективных степеней свободы рассчитывается по формуле Уэлча-Саттертуэйта (с округлением полученного значения до ближайшего наименьшего целого числа)

8 Результат измерения

Окончательный результат измерения влажности пенополитирола представляют в виде:

«Влажность пенополитирола составила (1,3 ± 0,3) %, где число, следующее за знаком ±, является численным значением расширенной неопределенности, которая получена умножением суммарной стандартной неопределенности на коэффициент охвата k = 2, основанный на предполагаемом нормальном распределении, и определяет интервал, соответствующий вероятности охвата 95 %».

Библиография

ГОСТ 17177-94 Материалы и изделия строительные теплоизоляционные. Методы испытаний

ГОСТ 34100.3-20017/ISO/IEC Guide 98-3:2008 Неопределенность измерения. Часть 3. Руководство по выражению неопределенности измерения

ГОСТ Р ИСО 5725-6-2002 Точность (правильность и прецизионность) методов и результатов измерений. Часть 6. Использование значений точности на практике

Ближайшее обучение

Успейте записаться на ближайшее мероприятие

Популярные услуги