Повышение точности метода измерения теплофизических свойств за счет выбора рациональных параметров проведения эксперимента и обработки опытных данных Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

Автоматика. Информатика. Управление. Приборы

УДК 519.6:536.21.022/023

ПОВЫШЕНИЕ ТОЧНОСТИ МЕТОДА ИЗМЕРЕНИЯ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ ЗА СЧЕТ ВЫБОРА РАЦИОНАЛЬНЫХ ПАРАМЕТРОВ ПРОВЕДЕНИЯ

ЭКСПЕРИМЕНТА И ОБРАБОТКИ ОПЫТНЫХ ДАННЫХ

С.В. Пономарев1, П.В. Балабанов1, В.Ф. Сорочинский2,

А.С. Щекочихин1

Кафедра «Автоматизированные системы и приборы», ГОУВПО «ТГТУ» (1); ГНУ «Всероссийский научно-исследовательский институт зерна и продуктов его переработки», Россельхозакадемия, г. Москва (2); pav-balabanov@yandex.ru

Ключевые слова и фразы: задача Штурма-Лиувилля; повышение точности измерения; регулярный режим первого рода.

Аннотация: Развита теория метода двух альфа для измерения комплекса теплофизических характеристик. Предложены способы обработки экспериментальных данных в зависимости от вариантов проведения экспериментов. На основе оценки погрешности определения теплофизических характеристик даны предложения по повышению точности измерения.

Введение

Целью настоящей работы являляется разработка метода и измерительной установки для исследования теплофизических свойств (ТФС) зерна нута (ЗН), представляющего собой ценную пищевую культуру семейства бобовых. Информация о свойствах требуется, в частности, для расчета режимов сушки зерна.

Плод нута имеет близкую к шаровой форму диаметром 6,5...7,5 мм, поэтому для измерения ТФС отдельных зерен применимы известные методы и устройства [1, 2], разработанные для измерения свойств тел шаровой формы. Наиболее простыми с точки зрения приборной реализации и дающими информацию о комплексе ТФС (теплопроводность, температуропроводность, объемная теплоемкость) являются методы, основанные на теории регулярного режима первого рода [3]. Для реализации таких методов наиболее часто используют установки, в состав которых входят калориметры различной конструкции (а-калориметры, ламбда-калориметры, бикалориметры) и термостаты (жидкостные или воздушные). При измерениях калориметр термостатируют, а затем помещают в среду с температурой, отличной от начальной температуры, и с постоянным значением коэффициента теплоотдачи. В процессе остывания (нагрева) определяют темп охлаждения (нагрева), а затем вычисляют искомую температуропроводность, теплопроводность или комплекс ТФС.

Калориметр, как правило, представляет собой металлическую оболочку, в которую помещают исследуемый материал с измерителями температуры. Однако при изготовлении калориметра для измерения ТФС ЗН практически невозможно обеспечить хороший тепловой контакт между металлической оболочкой и поверхностью зерна (так как зерно не является идеальным шаром). Поэтому в качестве прототипа для вновь разрабатываемого метода измерения ТФС ЗН был взят метод, известный в литературе как метод двух альфа [1, 4], характерной особенностью которого является необходимость проведения двух экспериментов при различных значениях коэффициентов теплоотдачи на внешних поверхностях исследуемых образцов.

Теоретические основы метода измерения

Для измерения ТФС ЗН в центр испытуемого зерна помещалась термопара медь - константан, сваренная встык (рис. 1). Диаметр проводов 0,1 мм. Для испытаний использовались зёрна диаметром (7 ± 0,3) мм. Каждое зерно разрезалось на две доли, в центр зерна 1 помещался спай термопары (см. рис. 1), после чего доли склеивались по периметру, а для уменьшения утечек тепла через термопарные провода их помещали внутрь зёрен нута 2 того же размера, что и испытуемое.

Дополнительно к испытуемой системе изготавливалась эталонная система, отличающаяся от испытуемой тем, что вместо образцов из нута в ней использовались образцы из эталонного материала (полиметилметакрилата) диаметром 7 мм с известными ТФС.

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

Сформулируем задачу теплопереноса в испытуемом и эталонном образцах.

Пусть два образца шаровой формы радиуса R с начальной температурой 70 помещаются в начальный момент времени в среду с постоянной температурой 7с Ф То. Тогда температурные поля в образцах будут описываться дифференциальными уравнениями:

дТ} (г, т) _ ^ д_

дт ] г 2 дг

с начальными условиями

Т} (г,0) _ То,

Рис. 1. Конструкция испытуемой системы:

1 - испытуемый образец; 2 - вспомогательные образцы; 3 - провод из меди; 4 - провод из константана; 5 - рабочий спай термопары

dTj (r, т) dr

j = 1,2, 0 < r < R, т> 0,

2

r

с граничными условиями:

dT, (0, т)

, dTJ т) (т T (R ))

х J dr =а(Тс - TJ (R,т))

дг

дТ} (Я, т)

дг

где а у, Xj - температуропроводность и теплопроводность у-го образца; а - коэффициент теплоотдачи.

Замена переменных © у _ (Ту (г, т)-Тс )/(Т0 - Тс), Бо у _ ау- т/я2, г _ г / Я, Ыу _ аЯ/А,у позволит перейти к безразмерным краевым задачам:

d© j (r ,Fo j ) = ^д dFo r2 dr

r2 J

g© j (r ,Fo j )' dr

© J (r,0) = 1 d©, (0,Fo,)

j = 1,2, 0 < r < 1, Fo j > 0;

dr

d© j (1,Fo j )

J/ + Bi J © J (1,Fo J )= 0,

dr J J

общее решение которых

© у (° у )_Е Спу Хпу е"Е * Р°;, (1)

п _1

где Спу - постоянные, определяемые из начальных условий; е^-, хпу - собственные значения и собственные функции краевой задачи Штурма-Лиувилля:

2

х'Пу (г)+■= х'п,-)+еПу Хпу(г)_0; (2)

хПу (0)_ 0; (3)

Х'пу (1) + В1 уХпу (1)_ 0. (4)

Известно [4, 5], что ряды (1) являются знакопеременными, а их собственные 2

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

*

р° у > Б° у всеми членами рядов (1), кроме первых, можно пренебречь и при наступлении в образцах регулярного теплового режима первого рода считать, что температурные поля в них вычисляются по формулам:

© у (,р° у)_ С1 уХ1у(г)е ^ Р°1, (5)

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

Лиувилля (2) - (4).

Прологарифмировав выражения (5), получим

1П © у (Г , Бо) _ 1П С1 у Х1 у (Г ) - Бо у .

Пусть индекс у' _ 1 соответствует эталонному образцу, а у' _ 2 - опытному образцу. Если по данным эксперимента с эталонным образцом построить график зависимости 1п©1(0,Б°1)_ /(Б°1) (рис. 2), то можно вычислить первое собственное значение задачи Штурма-Лиувилля как тангенс угла наклона его прямолинейного участка

о / ** ( ** \ * ( * \\ / / ** * \

— 8^1 _(1п©* (0,Б°* )- 1п©1(0,Р°1 ))/(Б°1 -Бо*). (6)

После чего из решения задачи Штурма-Лиувилля (2) - (4), можно определить Б^, а затем по формуле а_ Я вычислить коэффициент теплоотдачи.

По данным эксперимента с опытным образцом можно построить график зависимости 1п ©2 (0, т)_ /(т), из которого определить величину - 8^2/Я2 как тангенс угла наклона его прямолинейного участка

-81^2/Я2 _(1п©2* (0, т**)- 1п©2 (0,т*))/(т** -т*).

Обозначим

2

^.2 L _**/ **\ _ */ *W// ** *\

A = R (in©2 (0,т )-In©2(0,т ))(т -т ).

(7)

Если эксперимент проводить при Ы ^ да, то граничное условие (4) для опытного образца примет вид

Х12 (1)_ 0,

а первое собственное значение будет равно п2.

Из эксперимента с опытным образцом (при Ы ^ да) определяется А, после чего вычисляется искомая температуропроводность опытного образца

$2 — —A /822 — —A /п .

(8)

Для определения теплопроводности X 2 опытного образца необходимо провести эксперимент при Ы << 101 (например, в условиях естественной конвекции). По данным эксперимента с эталоном определяют коэффициент теплоотдачи а, по

данным эксперимента с опытным образцом определяют А, вычисляют 8ц _ -А / а2

и подбирают значение ^2, удовлетворяющее задаче (2) - (4). Затем вычисляют теплопроводность

^2 _ аЯ/В12 . (9)

Одним из недостатков изложенного выше метода измерения (применительно к зёрнам нута) является необходимость проведения эксперимента при В1 ^ да, что на практике достигается путем интенсивного обдува образцов воздухом, за счет чего изменяется их влажность (на 3-5 % в ходе опыта). Поэтому возможен

подход, при котором проводят два эксперимента при двух различных значениях

коэффициентов теплоотдачи а 1 и а 2 (при условии, что В1 < 101).

Рассмотрим краевую задачу (2) - (4) для первого и второго экспериментов с

опытным образцом:

2 А, —

Х'к(г) + - Х'хк(г)--L %1к(г)_ 0, к _1,2 (10)

Г а2

х1к (0)_ 0; (11)

х'1 к (1)+^^кЯ Х1к (1)“ 0, (12)

Я 2

где к - номер эксперимента.

Систему граничных условий (12) преобразуем к виду

а1 Х±(1)_а2 _ В. (13)

Х11(1) Х12 (1)

Решая численным методом уравнения (10) можно подобрать значение а2, удовлетворяющее условию (13). Далее вычисляют теплопроводность по формуле

Я 2 _-Я/В. (14)

Проведение экспериментов и обработка данных

Для вычисления комплекса ТФС (температуропроводности и теплопроводности) необходимо провести два эксперимента. Возможно два варианта проведения экспериментов, но, не зависимо от вариантов, перед началом каждого эксперимента эталонную и опытную системы выдерживают в воздушном термостате до достижения постоянной температуры 70 .

1-й вариант проведения экспериментов. В первом эксперименте эталонную и опытную системы помещают в среду с 7с Ф 7) при В1 ^ да. На практике разность температур 7) -7с = 7...10 °С, а В1 и 102. Вычисление температуропроводности осуществляют по формуле (8). Во втором эксперименте эталонную и опытную системы помещают в среду с температурой 7с Ф 70, но при В1 и 0,2...10. Вычисление теплопроводности осуществляют по формуле (9).

2-й вариант проведения экспериментов. В первом эксперименте эталонную и опытную системы помещают в среду с 70 - 7с = 7.10 °С и при В1 и 50, а во втором эксперименте при В1 и 0,2...10. По полученным в экспериментах данным вычисляют А1 и А2 , решают задачи (10) - (12), подбирая значение а2, удовлетворяющее условию (13), и вычисляют теплопроводность по формуле (14).

Выбор рациональных параметров проведения эксперимента и обработки опытных данных

При выборе рациональных параметров проведения эксперимента и обработки опытных данных по известным методикам оценивались теоретические погрешности определения ТФС, а затем отыскивались критерии рационального окончания эксперимента [6] и диапазон экспериментальных данных, использование которых при вычислениях ТФС дает наименьшие погрешности вычисления.

Для оценки теоретической погрешности определения температуропроводности прологарифмируем выражение (8) и после дифференцирования получим

^ (г?2 )

da2

$2

dA

A

Заменим в последнем выражении дифференциалы на конечные разности, а знак «-» на «+»

Да2 _ АА . )

a2

A

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

А$2

а2 /ск

1

АЛ

A

^ ) W.

(15)

Проведя аналогичные преобразования с выражением (7), получим среднеквадратичную оценку погрешности измерения величины А

AA

~A

1

( ©** Л2

in ©*

V © /

^А©** Л

©**

*2 ' А© Л

©*

( AR

+1 "R~

2 (

А

** * (т*-т )

** * т -т

Из полученных выражений видно, что на погрешность измерения температуропроводности оказывают влияние погрешности измерения температур, геометрических размеров, времени, а также погрешность А(е22). Погрешности измерения геометрических размеров и времени являются небольшими по величине, поэтому пренебрегаем ими. Оценим вклад члена д(е22) в общую погрешность измерения величины А. При В1 _ 51 первое собственное значение 8ц = 3,0801, а при В1 = 101 - 8ц = 3,1105. Поэтому, даже если эксперимент будет проведен при 51 < В1 < 101, то при вычислении температуропроводности по формуле (8) погрешность А^ц )/(822) составит менее 2 %, и с увеличением числа Био будет

уменьшаться. Поэтому, пренебрегая погрешностью ), формулу (15) шем в виде

запи-

Аа

___2

а2 /ск

AA

~

А©

©**

2

(А© Л

©*

in

©_

©*

(16)

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

Из формулы (9) видно, что теплопроводность Я2 является функцией коэффициента теплоотдачи, радиуса и числа Био. Значение В12 подбирается путем

2

+

2

+

+

ск

2

2

+

численного решения задачи (2) - (4), которое решается с высокой степенью точности современными прикладными программами такими, например, как МаШСай Поэтому погрешностью определения Б12 можно пренебречь. Погрешность измерения радиуса Я является небольшой по величине, поэтому пренебрегаем и ей. Следовательно, среднеквадратичная погрешность теплопроводности определяется, главным образом, погрешностью измерения коэффициента теплоотдачи

АХ; X 2

Аа

а

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

является функцией измеряемых в эксперименте температур 0 грешностей измерения этих температур

0

а также по-

АХ

X 2

і

А0

0**

Л

2

2

0*

1п

0_

0*

(17)

* ** , ,

Найдем пару значений 0 , 0 , при которых (ДХ2/ X2 )ск ^ шт.

Для этого оценим погрешности вычисления безразмерных температур

А0

0

1

А(Т(г, т)-Тс ) /А(7Ъ - Тс)

+ 5;.

(Т(г, т)-Тс) 1(Т0 - Тс)

В последнем выражении учтена погрешность 5о, обусловленная отбрасыванием всех членов ряда Фурье (1) кроме первого.

Выполненные расчеты (Д0 / 0)ск = /(0) в диапазоне 0,25 < Ы < 101 и

5 < То - Тс < 10 показали (рис. 3), что при 0,2 < 0 < 0,6 погрешность вычисления безразмерной температуры минимальна, а следовательно можно утверждать, что выполняется условие (ДХ2/X2 )ск ^ шт.

Таким образом, для автоматизации процесса обработки экспериментальных данных на ЭВМ можно программно отслеживать значение безразмерной темпера-

*

туры 0. Значение 0 = 0,2 служит критерием окончания эксперимента, а для расчетов ТФС следует использовать данные из диапазона 0,2 < 0 < 0,6.

(А0/0)с

Рис. 3. График зависимости (А0/0)ск = / (0) при ВІ = 101 и Г0- Тс = 7 °С

ск

ск

+

ск

Опытная установка и результаты измерения ТФС зерен нута

Для реализации предложенного метода измерения была собрана опытная установка (рис. 4), состоящая из воздушного термостата ТС-80М, двух систем образцов с термопарами, подключенными к плате сбора данных РСІ 6010. Результаты эксперимента выводятся на монитор компьютера.

Рассмотрим пример обработки данных одного из экспериментов.

Температурные кривые, полученные по результатам экспериментов с эталонной системой при различных значениях числа Ві, представлены на рис. 5, а.

1

У

РСІ

6010

ПК

5

Рис. 4. Схема опытной установки:

1 - термостат ТС-80М; 2 - опытная система; 3 - эталонная система; 4 - термопары (медь - константан); 5 - плата сбора данных РСІ-6010

©і

Рис. 5. Графики зависимостей 01(0, Fo1) =/(Ро^ (а) и 1п 01(О, Ро1) = /(Ро1) (б) для эталонной системы

По данным экспериментов с эталоном по формуле (6) вычисляют первое собственное значение задачи Штурма-Лиувилля е21 как тангенса угла наклона прямолинейного участка графика функции 1п ©1(0,Ро1 ) = / (Ро1), представленного на рис. 5, б. Полученные в ходе экспериментов данные показывают, что в диапазоне 0,2 <©< 0,6 находится прямолинейный участок зависимости 1п ©1(0, Бо1 ) = = / (Бо^), что подтверждает выполненные выше теоретические расчеты.

Из решения задачи (2) - (4) определяют Ы1, а далее коэффициент теплоотдачи а.

Результаты экспериментов с опытной системой представлены на рис. 6, а.

По полученным данным были построены графики зависимости (рис. 6, б) 1п ©2 (0, т) = / (т), из которых определены значения параметра А.

Сводные результаты обработки экспериментальных данных представлены в табл. 1.

В результате расчетов получили следующие теплофизические свойства исследуемых образцов нута X = 0,24 Вт/(м-К), а = 1,8-10 7 м2/с.

02

1п @2(0, т) = f (т) (б) для опытной системы

Таблица 1

Сводные результаты экспериментов

№ опыта 2 S11 Bi1 а AIR2 A-107

1 0,64 0,25 14 - 0,009 1,10

2 2,83 1,2 67 - 0,036 4,41

3 9,30 50,0 2785 - 0,14 17,15

Результаты измерения ТФС образцов нута при параметрах эксперимента 0,25 < Ы < 51 и 5 < Т00 - Тс < 10 показали, что при использовании для расчетов ТФС экспериментальных данных из диапазона 0,2 < © < 0,6 погрешности вычисления теплопроводности и температуропроводности будут минимальными.

Таким образом, для автоматизации процесса обработки экспериментальных данных на ЭВМ можно программно отслеживать значение безразмерной температуры ©. Значение ©** = 0,2 служит критерием окончания эксперимента, а для расчетов ТФС следует использовать данные из диапазона 0,2 < © < 0,6.

Список литературы

1. Осипова, В. А. Экспериментальное исследование процессов теплообмена /

B. А. Осипова. - 2-е изд., перераб. и доп. - М. : Энергия, 1969. - 392 с.

2. Чиркин, В.С. Теплопроводность промышленных материалов / В.С. Чир-кин. - М. : Машгиз, 1962. - 245 с.

3. Лыков, А. В. Теория теплопроводности / А.В. Лыков - М. : Высшая школа, 1967. - 599 с.

4. Теоретические и практические основы теплофизических измерений /

C.В. Пономарев [и др.]. - М. : Физматлит, 2008. - 408 с.

5. Карташов, Э.М. Аналитические методы в теории теплопроводности твердых тел : учеб. пособие / Э.М. Карташов. - 2-е изд., доп. - М. : Высшая школа, 1985. - 480 с.

6. Балабанов, П.В. Повышение точности измерения теплопроводности путем введения критерия управления ходом измерения / П.В. Балабанов, А.В. Трофимов // Метрология, стандартизация, сертификация и управление качеством продукции : программа, материалы шк.-семинара молодых ученых / Тамб. гос. техн. ун-т. -Тамбов, 2003. - С. 137.

Improvement of Accuracy of Measuring Technique of Thermo-Physical

Properties through Selection of Rational Parameters of Operation and Experimental Data Processing

S.V. Ponomarev1, P.V. Balabanov1, V.F. Sorochinsky2, A.S. Shchekochikhin1

Department “Automated Systems and Devices ", TSTU (1);

All-Russian Research Institute of Grain and Products of its Processing, Russian Agricultural Academy (2); pav-balabanov@yandex.ru

Key words and phrases: measurement accuracy improvement; regular model of 1st type; Sturm-Liouville problem.

Abstract: The paper develops the theory of two-alpha method for measuring a set of thermo-physical characteristics. The ways of processing experimental data depending on the procedure of experiment running are proposed. On the basis of error estimation for determining thermo-physical characteristics the ways of improving measurement accuracy are put forward.

Erhohung der Exaktheit der Methode der Messung der warmephysikalischen Eigenschaften infolge der Auswahl der rationelen Parameter der Experimentdurchfurung und der Bearbeitung der Experimentangaben

Zusammenfassung: Es ist die Theorie der Methode der zwei Alfa fur die Messung des Komplexes der warmephysikalischen Charakteristiken entwickelt. Es sind die Arten der Bearbeitung der Experimentalangaben je nach der Experimentdurchfuhrung vorgeschlagen. Auf Grund der Einschatzung des Fehlers der Bestimmung der warmephysikalischen Charakteristiken sind die Vorschlage fur die Erhohung der Messungsexaktheit angegeben.

Augmentation de la methode de la mesure des proprietes thermophysiques au depends du choix des parametres rationnels de la mise en oeuvre de l’experiment et du traitement des donnees

Resume: Est evaluee la theorie de deux alphas pour la mesure du complexe des caracteristiques thermophysiques. Sont proposes les moyens du traitement des donnees experimentales compte tenue des variants la mise en oeuvre des experiments. A la base de revaluation de l’erreur de la definition des caracteristiques thermophysiques sont donnees les propositions sur l’augmentation de l’exactitude de la mesure.

Авторы: Пономарев Сергей Васильевич - доктор технических наук, профессор, заместитель заведующего кафедрой «Автоматизированные системы и приборы»; Балабанов Павел Владимирович - кандидат технических наук, доцент кафедры «Автоматизированные системы и приборы», ГОУ ВПО «ТГТУ»; Соро-чинский Владимир Федорович - доктор технических наук, заместитель директора ГНУ ВНИИЗ Россельхозакадемии; Щекочихин Андрей Сергеевич - магистрант группы НГ-61, ГОУ ВПО «ТГТУ».

Рецензент: Коновалов Виктор Иванович - доктор технических наук, профессор кафедры «Химическая инженерия», ГОУ ВПО «ТГТУ».