Рекомендовано к изданию решением секции N 4 научно-технического совета НИИСФ Госстроя СССР.

Описаны методика и измерительная аппаратура для комплексного определения теплофизических характеристик материалов импульсным методом. Изложены физические концепции метода, причины появления и способы компенсации систематических ошибок. Для обработки полученных результатов приведена программа для ЭВМ.

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

Табл. 2, ил. 6.

В НИИСФ проведен комплекс теоретических и экспериментальных исследований, связанных с физическим обоснованием импульсного метода комплексного определения теплофизических характеристик строительных материалов, выявлением систематических ошибок измерений, их оценкой и уравновешиванием. Метод позволяет применять стандартную измерительную аппаратуру в совокупности с несложной дополнительной оснасткой. Для обработки получаемого массива экспериментальных данных на ЭВМ разработана программа NNK1, составленная на языке Фортран-IV, а также для работы в режиме диалога с ЭВМ на языке Бейсик-II.

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

Рекомендации разработаны в лаборатории долговечности ограждающих конструкций НИИСФ канд. техн. наук Ю.Д. Ясиным и инж. Н.Н. Кузнецовой и содержат конкретные указания по использованию методики и аппаратуры.

1.1. К основным теплофизическим характеристикам материалов, определяемым импульсным методом, относятся:

- объемная теплоемкость, Дж/(м³·К);

a - коэффициент температуропроводности, м²/с;

- коэффициент теплопроводности, Вт/(м·К);

b - коэффициент тепловой активности, Дж/(м²·с^1/2·К).

Для нахождения этих характеристик достаточно определить только две любые из них, остальные рассчитывают с помощью известных уравнений связи

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

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

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

1.5. Разработанные НИИСФ способ и устройство для определения теплофизических характеристик строительных материалов <1> обладают рядом неоспоримых преимуществ. Они позволяют экспрессно и комплексно определять теплофизические характеристики твердых, в том числе влажных капиллярно-пористых строительных материалов на образцах или изделиях без их нарушения и при одностороннем доступе к ним.

--------------------------------

<1> А.с. N 1122956 (СССР). Способ определения теплофизических характеристик материалов конструкций/Ю.Д. Ясин, Н.Н. Кузнецова. Опубл. в Б. И., 1984, N 41.

А.с. N 1107036 (СССР). Устройство для определения теплофизических характеристик материалов конструкций/Ю.Д. Ясин, Н.Н. Кузнецова. Опубл. в Б. И., 1984, N 29.

2.1. Метод основан на частном решении теплофизической задачи об изменении избыточной температуры во времени в центре плоского импульсного источника тепла центральносимметричной формы и конечных размеров, расположенного в плоскости соприкосновения двух (эталонного и исследуемого) полубесконечных тел. Это решение получено в виде

где - изменяющаяся во времени избыточная температура в центре плоского импульсного источника тепла, К;

Q - количество тепла, выделенное единицей площади импульсного источника, Дж/м²;

R - радиус плоского кругового импульсного источника, эквивалентного источнику центральносимметричной формы, м; индексы "э" и "м" относятся к эталонному и исследуемому телам соответственно.

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

2.3. Практически выражения (1) и (2) равны, если невязка в значениях избыточных температур, рассчитанных по ним, не превышает предельной абсолютной погрешности измерения температуры , величина которой определяется по паспортным данным измерительного прибора. В предельном случае при это может быть представлено следующим образом:

Таким образом, - момент времени, исчисляемый с момента возникновения импульса и разграничивающий время действия рассматриваемого теплового процесса на два интервала: и , т.е., когда этот процесс описывается выражениями (2) и (1) соответственно.

2.4. Значение определяют из уравнения (3), например, методом последовательных приближений

где i = 1...n - номер итерации; n определяется из условия , где - предельная абсолютная погрешность измерения , которая определяется по паспортным данным измерительного устройства.

2.5. Расчетные формулы, полученные из выражений (1), (2) с учетом (4), имеют вид:

для коэффициента тепловой активности

для коэффициента температуропроводности

В выражении (6) используют значение b_м, найденное предварительно по формуле (5).

2.6. Коэффициент теплопроводности и объемную теплоемкость с использованием уже найденных коэффициентов тепловой активности в температуропроводности определяют по формулам:

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

где x₁, ..., x_n - независимые аргументы (они же измеряемые параметры), предельная абсолютная ошибка определяемого параметра может быть выражена через предельные абсолютные ошибки измеряемых параметров по формуле

3.2. В общем случае предельная абсолютная ошибка может включать в себя систематическую и случайную составляющие:

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

3.3. После устранения систематических ошибок выражение (9) имеет вид

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

3.4. Кроме выражения (10) для априорных оценок могут быть использованы выражения:

или

Здесь - предельная относительная ошибка параметра; - его среднее значение в определяемом диапазоне изменения; , - среднеквадратические отклонения, .

3.5. Расчетные выражения (5), (6) получены на основе решения идеализированной математической модели. Это может быть причиной образования ряда систематических ошибок, наиболее важные из которых предопределены:

толщиной и теплофизическими свойствами плоского нагревателя;

отсутствием идеального теплового контакта между плоскостями нагревателя и исследуемого материала;

фактической длительностью теплового импульса;

конечными размерами соприкасающихся тел и т.д.

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

где H = (b_м - b_н)/(b_м + b_н), н - индекс принадлежности нагревателю.

Пренебрегая наличием нагревателя, т.е. полагая b_н = b_м = b, из выражения (13) получим

Относительная ошибка , обусловленная наличием нагревателя, определяется с учетом (13) и (14) по формуле

где .

3.7. Графическая зависимость представленная на рис. 1, позволяет выбрать время измерения или подбирать параметры нагревателя таким образом, чтобы ошибка, обусловленная его влиянием, не превышала наперед заданной величины.

Пример. Нагреватель выполнен из фольги толщиной , теплофизические характеристики которой в принятых единицах измерений находятся в диапазонах: ; 5·10^-6 < a < 1·10^-4; 7·10³ < b < 25·10³. Этим условиям удовлетворяют металлы с большим удельным сопротивлением (Ом), которые и целесообразно использовать при изготовлении нагревателей. Для такого нагревателя относительная ошибка , монотонно убывая, уже через 2 с после подачи импульса становится меньше 2·10^-2.

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

3.9. Влияние воздушных зазоров можно свести до минимума путем их заполнения консистентной смазкой с известными теплофизическими свойствами <1>, наличие которой в виде дополнительного слоя учитывают используя решение

где H = (b_м - b_см)/(b_м + b_см), .

Относительная ошибка , обусловленная наличием дополнительного слоя смазки толщиной , с учетом (2), (15) и (16), определяется по формуле

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

3.10. Поскольку значения тепловой активности консистентных смазок (солидол, графитная и контактная смазка, унитол, литол и др.) <1> находятся в диапазоне 500 - 800 Дж/(м·с^1/2·К) и близки к значениям тепловой активности большинства строительных материалов, относительная ошибка, обусловленная наличием слоя смазки, будет, как правило, незначительна.

--------------------------------

<1> Фройштетер Г.В., Трилиский К.К., Ищук Ю.А., Ступак П.П. Реологические и теплофизические свойства пластичных смазок. - М.: Химия, 1980.

Пример. При исследовании керамзитобетонов с коэффициентом тепловой активности 450 < b_м < 750 <1> с использованием в качестве заполнения воздушных прослоек солидола (b_см = 559 Дж/(м·с^1/2·К), согласно графической зависимости, представленной на рис. 1, для |H| < 0,15, (что соответствует рассматриваемому случаю) величина относительной ошибки будет меньше , начиная с , или при толщине слоя смазки 2·10^-4 (м) и температуропроводности a_см = 7·10^-8 (м²/с), начиная с момента времени .

--------------------------------

<1> Глава СНиП II-3-79*. Строительная теплотехника. - М.: Стройиздат, 1982.

3.11. Когда относительная ошибка превышает допустимое значение, а и H > 0, коэффициент тепловой активности материала с учетом влияния смазки рассчитывают по формуле

Здесь значение H определяют из выражения

которое получено из (16) при учете одного члена ряда.

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

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

3.14. Математически это можно выразить как включение в момент времени дополнительного отрицательного источника (стока) тепла равной (по абсолютной величине) мощности.

Избыточная температура нагревателя изменяется во времени под действием источника постоянной мощности по выражению

где q - удельная мощность нагревателя, Вт/м².

То же, но с учетом воздействия отрицательного источника той же мощности, включенного в момент времени ,

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

3.15. Относительная ошибка обусловленная фактической длительностью импульса, определяется по формуле

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

откуда следует, что для , а для ошибка определяется выражением

где .

Зависимость (19) графически представлена на рис. 2. С помощью последнего можно получить представление о том интервале времени, в котором предельная относительная ошибка не будет превышать заданной величины.

3.16. Используемое в методе решение для системы полубесконечных тел будет справедливо в течение некоторого времени - времени сохранения условия полубесконечности исследуемого образца. Эти условия определяют, используя известное решение, описывающее распределение избыточной температуры в толще материала по нормали к плоскости нагревателя

где Q_м - количество тепла, выделенное единицей площади импульсного источника в тело исследуемого материала, Дж/м². Путем решения уравнения при z = z_м находят значение , соответствующее появлению максимального значения температуры в плоскости z = z_м

где z_м - толщина исследуемого образца или максимальная глубина проникновения импульса.

Полагая в (20), что при в плоскости z = z_м уровень полезного сигнала не превышает порога чувствительности прибора , получим уравнение, определяющее условия полубесконечности исследуемого образца

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

Пример. Определить время , в течение которого удовлетворяется условие полубесконечности для облицовки из декоративного камня толщиной z_м = 0,015 м с коэффициентами тепловой активности b_м = 1000 Дж/(м·с^1/2·К) и температуропроводности a_м = 4·10^-7 м²/с при нагревателе размером R = 0,02 м, выделяющем единицей площади количество тепла Q_м = 800 Дж/м², при ошибке измерения температуры . Искомая величина с учетом условий задачи и расчетной формулы (22) равна 98 с (с округлением до целой секунды), что удовлетворяет условию полубесконечности, поскольку , вычисленное по формуле (21), равно 156 с.

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

3.18. Следуя (11), выражение для нахождения предельной относительной ошибки определения коэффициента тепловой активности исследуемого материала с учетом (5) будет иметь вид

где , - предельные относительные ошибки определения коэффициентов тепловой активности исследуемого и "эталонного" материала соответственно;

, - предельные относительные ошибки измеряемых параметров: температуры и времени соответственно.

3.19. В свою очередь предельная относительная ошибка определения коэффициента тепловой активности "эталонного" материала с учетом (5), в котором положено b_м = b_э, может быть найдена из выражения

С учетом (24) выражение (23) преобразуется к виду

3.20. Поскольку избыточная температура нагревателя изменяется во временном интервале , ее среднее значение рассчитывается по выражению

тогда предельная относительная ошибка определяется по формуле

Формула для определения предельной относительной ошибки , полученная аналогичным образом, имеет вид

3.21. Вводя безразмерное время в выражения (26) и (27), получим, соответственно:

для предельной относительной ошибки 

для предельной относительной ошибки 

Здесь

и

3.22. Исходя из характера изменения предельных относительных ошибок и (рис. 3) их суммарная ошибка будет иметь минимум, который можно определить исходя из условия . Учитывая последнее, будем иметь уравнение

решая которое относительно , например, итерационным методом при заданных значениях K₁ и K₂ можно определить оптимальное отношение .

3.23. Поскольку количество тепла, выделенное единицей площади импульсного источника Q, определяется из выражения

где W_н - мощность нагревателя, Вт; U - напряжение на клеммах нагревателя, В; I - сила тока, протекающего через нагреватель, А; S_н - площадь нагревателя, м²; r_н - сопротивление нагревателя, Ом, предельная относительная ошибка в соответствии с измеряемыми параметрами определяется по формулам

3.24. Предельная относительная ошибка определения коэффициента температуропроводности эталонного материала, согласно уравнению (1), в котором положено b_м = b_э, a_м = a_э, может быть определена по формуле

где

Предельная относительная ошибка определения коэффициента температуропроводности исследуемого материала с учетом (33) может определяться из выражения

где

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

Блок-схема программы приведена на рис. 4.

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

4.3. Программа работает следующим образом. После ввода исходных данных, в качестве первого приближения назначается величина , соответствующая , начиная с , и задаются начальные средние значения и . Затем для фиксированных избыточных температур нагревателя и соответствующих моментов времени по формулам (5) и (6) вычисляют значения b_м, a_м и их средние , . Полученные значения сравниваются с их первым или последующими приближениями и и, если условия

не выполняются, то происходит уточнение значений по итерационной формуле (4) в цикле, выход из которого осуществляется по условию , где присваивается значение предыдущего приближения . Таким же образом значениям и присваиваются значения и .

4.4. В случае выполнения условий (35), значения и выводятся на печать и вычисляются их среднеквадратические отклонения и относительные ошибки. Далее по формулам (7) и (8) определяются значения и . Все вычисленные значения также выводятся на печать.

4.5. С целью более компактной записи программы использованы две подпрограммы-функции FUN1 и FUN2, реализующие общие части выражений (4), (5) и (6).

Программа NNK1, написанная на языке Фортран-IV, приведена в прил. 1.

5.1. Для комплексного определения теплофизических характеристик необходимы: первичный преобразователь и вторичные устройства для создания теплового импульса и для измерения и регистрация изменения во времени избыточной температуры нагревателя первичного преобразователя.

5.2. Для измерения и регистрации температуры может быть использован ампервольтомметр Ф30К, в комплект которого входят измерительный прибор, транскриптор и цифропечатающая электроуправляемая машина. Аналогичное устройство можно собрать из отдельных приборов, например вольтметра универсального В7-21, транскриптора Ф5033 и цифропечатающей машины ЭУМ-23Д. Возможны и другие комбинации и замены аналогичных приборов.

5.3. Устройство для создания теплового импульса состоит из стабилизированного источника питания, например Б5-21, и таймера, задающего временной интервал нагрева. В качестве таймера может быть использовано устройство, разработанное Ройфе В.С., принципиальная электрическая схема которого и перечень основных элементов к этой схеме представлены в прил. 2.

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

где - относительная ошибка измерения температуры в момент времени .

Пример. Определить параметры первичного преобразователя для исследования материала с априорно заданными характеристиками a = 4·10^-7 (м²/с), b = 1000 Дж/(м²·с^1/2·К). Материалом для тела первичного преобразователя выбран полистирол со следующими теплофизическими характеристиками: a_э = 3·10^-7 (м²/с), b_э = 200 Дж/(м²·с^1/2·К). Измерение избыточной температуры в центре нагревателя производят измерительно-регистрирующим устройством, например Ф30К, с ее регистрацией через 5 с. При проведении эксперимента необходимо обеспечивать, с одной стороны - представительность экспериментального массива для статистической оценки искомых характеристик, а с другой - минимизировать время эксперимента, что обусловлено требованием экспрессности метода, а также условиями сохранения полубесконечности исследуемого образца и тела первичного преобразователя. С этой целью выбираем количество измерений N = 30, а . Поскольку значения тепловой активности определяются с большей надежностью, принимаем , что обеспечивает десять результатов измерений для вычисления коэффициента b_м и двадцать результатов измерений для вычисления коэффициента температуропроводности a_м. Наконец, для получения достоверных характеристик при абсолютной ошибке измерения температуры , потребуем, чтобы относительная ошибка не превышала 1·10^-2, при и 2·10^-2 при . Радиус нагревателя оценим по формуле (36), при этом получим R = 2·10^-2 м. Далее, с учетом заданных значений , , a_м, b_м, R по формуле (1) найдем значения Q, обеспечивающие соответствующие точности измерения и : и . Положим Q = 16000 Дж/м² и по итерационной формуле (4) уточним значение . Чтобы убедиться в сохранении условий полубесконечности исследуемого образца при заданных параметрах, предварительно оценим по формуле (20) при z = 0 и Q_м - количество тепла, выделенное единицей площади нагревателя в тело материала, которое составит 13400 Дж/м², а также определим по формуле (21) - время наступления максимального значения температуры в плоскости z = 0,05 м, равное 1122 с. Затем, пользуясь выражением (22), получим время, в течение которого сохраняются условия полубесконечности образца - , превышающее предполагаемое в эксперименте .

Для определения толщины эталонного тела первичного преобразователя воспользуемся уравнением (22), заменив в нем индексы материала (м) на индексы эталона (э) и приняв Q_э = Q - Q_м и определим глубину проникновения импульса в тело эталона за время , которая составит z_э = 0,026 м.

С целью расширения возможностей разрабатываемого первичного преобразователя выберем z_э = 0,05 м, тогда в этой плоскости время наступления составит 1496 с, а величина не будет превышать , что обеспечит сохранение условий полубесконечности при любом .

6.1. Для определения коэффициентов тепловой активности b_м и температуропроводности a_м эталонного тела первичного преобразователя последний устанавливают на поверхности образца из того же материала, что и тело первичного преобразователя так, как это схематически показано на рис. 5. При этом плоский нагреватель с датчиком температуры и поверхность исследуемого образца должны иметь тепловой контакт.

6.2. Далее включают устройство для задания теплового импульса, который через нагреватель передается в тело первичного преобразователя и исследуемый образец. Избыточную температуру в центре нагревателя регистрируют с помощью термопары и измерительного устройства с выводом результатов на печать.

6.3. Поскольку в данном случае исследуемый образец имеет те же теплофизические характеристики (b_м = b_э, a_м = a_э), расчетные формулы, полученные из выражений (1), (2) с учетом (4), имеют вид:

для коэффициента тепловой активности

для коэффициента температуропроводности

при этом в выражении (38) используют значения b_э, найденные предварительно по формуле (37).

Остальные теплофизические характеристики и рассчитывают по формулам (7) и (8).

6.4. Для определения теплофизических характеристик исследуемого материала (общий случай) производят испытание образца или изделия согласно пп. 6.1 и 6.2. Затем полученный массив экспериментальных данных обрабатывают на ЭВМ по программе NNK1.

6.5. Исходный массив экспериментальных данных набивают на стандартных перфокартах в следующей последовательности:

на первой карте - численное значение экспериментального массива N форматом I3, равное количеству измерений избыточной температуры нагревателя, проведенных с начала подачи импульса;

на второй карте - значения Q, AE, BE, DT, DTAU, R;

на последующих картах - массивы значений T и TAU по семь значений на каждой карте, причем начало массива совмещают с началом новой карты. Значения на второй и последующих картах набивают форматом E11.4.

6.6. Результаты расчета представляются на дисплее или, в виде распечатки ЭВМ, снабжены соответствующими заголовками, обозначениями и не требуют дополнительных разъяснений.

Пример. Результаты расчета на ЭВМ.

6.7. Необходимо регулярно проводить опытную проверку метода и устройства на образцовых мерах, согласно пп. 6.1 и 6.2, расчет их теплофизических характеристик, согласно п. 6.4.

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

где , - средние значения сравниваемых (известной и экспериментально определенной) теплофизических характеристик, образцовой меры соответственно; , - их выборочные дисперсии; n₁, n₂ - количество измерений; - число степеней свободы.

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

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

стекло ЛК-5 ;

оптическое стекло ;

оргстекло ;

коэффициент теплопроводности которых представлен метрологической службой. Из-за отсутствия других статистических характеристик образцовых мер положим S₁ = 0 и . Средние значения теплофизических характеристик этих материалов, экспериментально полученных предлагаемым методом на устройстве, представлены в таблице, где, кроме того, даны среднеквадратические отклонения и объем выборок.

Проверка согласия между имеющимися по данным метрологической службы значениями теплопроводности образцовых мер и соответствующими значениями, полученными по методу и на устройстве, показала, что сравниваемые значения статистически неразличимы, поскольку - критерий оказался меньше при общепринятом уровне значимости p = 0,05.