АМПЛИТУДНЫЙ ПОЛИГАРМОНИЧЕСКИЙ МЕТОД ТЕМПЕРАТУРНЫХ ВОЛН ДЛЯ КОНТРОЛЯ ТЕМПЕРАТУРОПРОВОДНОСТИ ТВЕРДЫХ ИЗОТРОПНЫХ МАТЕРИАЛОВ
Е.Л. Артюхина
Кафедра «Управление качеством и сертификация», ФГБОУВПО «ТГТУ»;
artyukhina@yahoo.com
Ключевые слова и фразы: амплитудный полигармонический метод; коэффициент температуропроводности; краевые задачи теплопереноса; математическое моделирование; температурные волны.
Аннотация: Разработан амплитудный полигармонический метод температурных волн. Получено решение инверсной коэффициентной задачи теплопереноса для цилиндрического однородного изотропного образца, на плоской поверхности которого задана периодически изменяющая температура, интеграл от которой по периоду равен нулю, а длина волны меньше высоты цилиндра.
Введение
Широкое применение методов температурных волн в практике измерения тепловых свойств материалов объясняется рядом неоспоримых преимуществ, среди которых высокая помехозащищенность, возможность проводить эксперименты при низком уровне температурных возмущений, независимость результатов измерений от начального распределения температур, высокая информативность [1].
Существующие технические средства контроля температуропроводности, реализующие метод температурных волн, позволяют проводить измерения тепловых свойств материалов в широком диапазоне их изменения (от теплоизоляторов до металлов).
Однако в абсолютном большинстве применяемые для контроля температуропроводности методы температурных волн являются моногармоническими. При этом используется два основных подхода к решению поставленной задачи. Первый заключается в зондировании исследуемого образца тепловым потоком одной частоты. В этом случае о тепловых свойствах материалов судят по реакции на тепловое возмущение. Во втором случае из зондирующего теплового потока и реакции на тепловое возмущение выделяют разложением в ряд Фурье одинаковые гармоники [2, 3], по которым восстанавливают искомый коэффициент температуропроводности. Очевидными недостатками описанных методов являются: использование только части информации о тепловом процессе и значительные погрешности, связанные с гармоническим анализом температур и тепловых потоков. Альтернативой изложенным способам определения температуропроводности являются полигармонические методы температурных волн [4, 5], лишенные описанных недостатков.
В значительном числе случаев при реализации метода температурных волн исследуемые образцы представляют собой цилиндрические диски, радиус которых существенно превосходит высоту. Для таких образцов рассмотрим задачу теплопереноса. При описании процесса теплопереноса в исследуемом образце предположим однородность и изотропность материала, отсутствие в нем физикохимических превращений, приводящих к зависимости свойств от времени. Будем
Т0
считать периодические процессы установившимися. Кроме того, | Tdт = 0, где
0
То - период изменения температуры.
Температура на плоской поверхности образца представляет собой периодическую функцию.
Математическая модель описанного процесса имеет вид:
дT(х, т) д T(х, т)
—^J- = a-----------(т > О, О < х <то),
дт дх2
T(х,О) = TG = О ,
T (О, т) = TAf (т),
где
f (т) =
1, кто < т < I k +1 іто ;
-1,
к = О, 1, 2,
(1)
(2)
(3)
(4)
к + — |то < т < (к + 1)то ,
Т(х,т) - температура тела в точке с координатой х в момент времени т, К; а - коэффициент температуропроводности, м2/с; Т0 - начальная температура тела, К; Тл - амплитуда колебания температуры, К; f (т) - некоторая периодическая функция времени, определенная ниже.
Для решения сформулированной задачи (1) - (3) воспользуемся результатом работы [6], полученным для случая f(т) = cos(ют - ф) в установившемся состоянии,
T(х, т) = Ta exp
cos
(5)
Этот результат позволяет получить решение задачи (1) - (3) для произволь-
2п
ной периодической функцииf(г) с периодом —, раскладывая эту функцию в ряд
ю
Фурье,
f (т) = ao + Iak cos(krox -фк). к=1
(6)
Общее решение, полученное на основании принципа суперпозиции, является суммой частных решений (5) для каждого члена ряда (6)
T (х, т) = ao + I aк exp
к=1
- хл
, 1кю ^
cos кют-фк - хJ—
V 2a
(7)
В случае, когда температура поверхности исследуемого образца изменяется по закону (4), решение краевой задачи теплопереноса (1) - (4) в соответствии с (7) имеет вид
— х
ют- х
то
Рис. 1. Плоская температурная волна в образце
T(х,т) = 4Гл у_±_ п k=02k +1
(
ехр
(2k + 1)п
- х-„
v . ато
sin
(4k + 2)пт (2k + 1)п
" ато
(8)
ато ~ т
Перейдя в уравнении (8) к безразмерным переменным Fo = ——, т= —
х2 то
^ ~ Т (х, т)
&( х, т) =----------, получим
&(Fo, ~) = — X 2гГ-
п t^n 2k +1
^ 1 f \lk + 1)п^ (
4
' exp
Fo
Sin
(4k + 2)п~ -
(2k + 1)п
Fo
(9)
Распределение относительных температур $(х,~) =---------,— по глубине исследуе-
ТА
мого образца иллюстрируется графиком (рис. 1).
Постановка и решение коэффициентной задачи теплопроводности
Уравнение температурной волны в безразмерном виде (8) является функцией двух переменных - числа Фурье и относительного времени. Для каждого заданного числаГо,, находим т/тахтакие что $(го/ ,~тах )=&тах/, где &тах/ - амплитуда /-й температурной волны. Таким образом, получаем двухмерный массив данных {Гоу,&тах/}. Аппроксимируя сформированный массив, получим функцию
Го = Го(^тах).
Поскольку период колебания т0 и расстояние от источника задания теплового возмущения до места расположения датчика температуры х в эксперименте заданы, искомый коэффициент температуропроводности а определяется по формуле
Го($ ) х 2
а =------тах---, где Го(^тах) - аппроксимация таблично заданной функции.
т0
Ниже рассмотрен пример реализации описанного метода. Основываясь на решении (8), рассчитываем $(Го/,~) для фиксированных значений Го в диапазоне Го [1...70] и ~ [0...1]. Строим графики зависимостей $(Го/ = со^^т/) (рис. 2).
s
Рис. 2. Температурные кривые й^о;, т;):
1 - 10; 2 - 25; 3 - 35; 4 - 55; 5 - 70
С использованием пакета МаШетайса определяем координаты максимумов функций температурных кривых &таХ;- (таблица).
Аппроксимируем массив данных {Fo;, smax;} выражением
Fo = 0,925 +
47,409
s
4
max; J
43,0423 15,3267 3,3964
—^---------+ -
sm
sm
sr
(10)
На рисунке 3 представлена аппроксимация таблично заданной зависимости амплитуды от числа Фурье йтах. ^ог-).
Так для стандартных образцов из полиме-
Координаты максимумов функций smax t
i Fo; s max;
1 1 0,2384
2 3 0,4866
3 5 0,5890
4 7 0,6477
5 10 0,7021
6 15 0,755
7 20 0,7874
8 25 0,8068
9 30 0,8271
10 35 0,8407
11 40 0,852
12 45 0,8506
13 50 0,8654
14 55 0,8689
15 60 0,8779
16 70 0,8843
тилметакрилата при x = 2,840
80 с, зна-
чение йтах составило 1,675, соответствующее ему Fo - 1,224, а найденное значение коэффициента температуропроводности а = 1,2-10-7 м2/с.
и 1ЛЯУ
Рис. 3. Аппроксимация таблично заданной зависимости йтах^о)
тах
м, Т
0
Разработанный амплитудный полигармонический метод температурных волн, обладая всеми преимуществами методов температурных волн, отличается от моногармонических методов на порядок более высокими значениями регистрируемых температур и существенно большим объемом информации, что позволяет повысить точность контроля температуропроводности и в то же время существенно сократить время измерений.
Работа выполнена при финансовой поддержке Министерства образования и науки Российской Федерации (государственный контракт № 16.526.11.6010 от 28.10.2011 г.).
Список литературы
1. Кравчун, С.Н. Метод периодического нагрева в экспериментальной теплофизике / С.Н. Кравчун, А.А. Липаев. - Казань : Изд-во Казан. ун-та, 2006. -208 с.
2. Angstrom, A.J. Neue Methode das Warmeleiltungs-vermogen der Korper zu bestimmen / A.J. Angstrom // Ann. d. Physik. - 1881. - Bd. 14. - S. 513-530.
3. Филиппов, Л.П. Измерение теплофизических свойств веществ методом периодического нагрева / Л.П. Филиппов. - М. : Энергоатомиздат, 1984. - 105 с.
4. Артюхина, Е.Л. Фазовый полигармонический метод температурных волн для контроля температуропроводности твердых изотропных материалов / Е.Л. Артюхина, С.В. Мищенко // Вопр. соврем. науки и практики. Ун-т им. В.И. Вернадского. - 2013. - № (45). - С. 48-52.
5. Артюхина, Е.Л. Теоретическое обоснование полигармонических методов температурных волн для контроля температуропроводности твердых изотропных материалов / Е.Л. Артюхина, С.В. Мищенко // Вестн. Тамб. гос. техн. ун-та. -2013. - Т. 19, № 1. - С. 30-37.
6. Лыков, А. В. Конвекция и тепловые волны / А.В. Лыков, Б.М. Берковский. - М. : Энергия, 1974. - 286 с.
Amplitude Polyharmonic Method of Temperature Control for Thermal Waves of Solid Isotropic Materials
E.L. Artyukhina
Department “Quality Management and Certification ", TSTU; artyukhina@yahoo. com
Key words and phrases: coefficient of thermal conductivity; heat transfer boundary value problems; mathematical modeling; polyharmonic amplitude method; temperature waves.
Abstract: A method of temperature polyharmonic amplitude waves has been developed. A solution of the inverse problem of heat transfer coefficient for a homogeneous isotropic cylindrical specimen on a flat surface has been obtained. Periodically changing temperature is given, its integral by period is zero and the wavelength is less than the height of the cylinder.
Polyharmonische Amplitudenmethode der Temperaturwellen fur die Kontrolle der Temperaturleitfahigkeit der festen Isotropmaterialien
Zusammenfassung: Es ist die polyharmonische Amplitudenmethode der Temperaturwellen entwickelt. Es ist die Losung der inversen Koeffizientaufgabe der Warmeubertragung fur das zylindrischen gleichartigen isotropen Muster erhalten, auf dessen flachen Oberflache die periodisch andernde Temperatur aufgegeben ist, das Integral von deren nach der Periode der Null gleich ist, und die Lange der Welle ist kleiner als Hohe des Zylinders.
Methode d’amplitude polyharmonieuse des ondes de temperature pour le controle de la conductibilie de temperature des materiaux solides isotropes
Resume: Est elaboree la methode d’amplitude polyharmonieuse pour les ondes de temperature. Est obtenue la solution du probleme inverse de coefficient du transfert de temperature pour un echantillon cylindrique isotrope sur la superficie plate duquel est donnee la temperature qui change periodiquement et dont l’integrale par periode est egal au zero et la longueur de l’onde est moindre que la hauteur du cylindre.
Автор: Артюхина Екатерина Леонидовна - аспирант кафедры «Управление качеством и сертификация», ФГБОУ ВПО «ТГТУ».
Рецензент: Пономарев Сергей Васильевич - доктор технических наук, профессор, исполняющий обязанности заведующего кафедрой «Управление качеством и сертификация», ФГБОУ ВПО «ТГТУ».