CC BY
29
УДК 678-1:303.094.7
ИМИТАЦИОННОЕ ИССЛЕДОВАНИЕ ТЕПЛОПЕРЕНОСА В МЕТОДЕ НЕРАЗРУШАЮЩЕГО КОНТРОЛЯ СТРУКТУРНЫХ ПЕРЕХОДОВ В ПОЛИМЕРАХ
Н.Ф. Майникова1, Н.П. Жуков2, А.А. Балашов2, Д.Г. Бородавкин2
Кафедры: «Теория машин, механизмов и детали машин» (1), «Гидравлика и теплотехника» (2), ГОУВПО «ТГТУ»
Представлена членом редколлегии профессором В. И. Коноваловым
Ключевые слова и фразы: неразрушающий контроль; плоское полупространство; структурный переход; теплофизические свойства; численное моделирование.
Аннотация: Проведено имитационное исследование метода неразрушаю-щего контроля структурных переходов в полимерных материалах.
Обозначения
с - удельная теплоемкость, Дж/(кг-К); Т - избыточная температура, °С;
д - тепловой поток, Вт/м2; е - тепловая активность, Вт-с0,5/(м2К);
К - радиус м; т - время, с.
г - расстояние, м;
Разработанный авторами метод [1-3], входящий в группу методов термического анализа, позволяет определять неразрушающим способом теплофизические свойства (ТФС) твердых материалов и фиксировать температурные характеристики структурных переходов в полимерных материалах (ПМ) по изменениям тепловой активности с ростом температуры. Теоретическую основу данного метода неразрушающего контроля (НК) составляют аналитические закономерности распространения одномерного плоского температурного поля от поверхностного круглого источника тепла постоянной мощности при локальной регуляризации тепловых потоков в твердых телах, которые могут быть отнесены к классу плоского полупространства [4-6]. Одномерное температурное поле было использовано в связи с тем, что двух- и трехмерность полей обычно не обеспечивает заметного улучшения эксплуатационных показателей метода и при этом, как правило, ухудшает его метрологические возможности.
Согласно измерительной схеме метода (рис. 1, а) тепловое воздействие на исследуемое полимерное тело, имеющее равномерное начальное температурное распределение, осуществляется с помощью нагревателя постоянной мощности, выполненного в виде тонкого диска радиусом К, встроенного в подложку измерительного зонда (ИЗ). Начальное температурное распределение контролируется одновременно несколькими (не менее трех) термоэлектрическими преобразователями (ТП), расположенными в центре нагревателя и на расстояниях г, от центра в плоскости контакта подложки ИЗ и исследуемого тела. ТП в ходе эксперимента
фиксируют термограммы - зависимости избыточной температуры Т (или темпе*
ратуры исследуемого изделия Т ) от времени.
Рис. 1. Измерительная (а) и тепловая (б) схемы
На рис. 2 представлены термограммы, полученные экспериментально на изделии из политетрафторэтилена (ПТФЭ) при следующих условиях: начальная
* *
температура Тн = 23 °С; Т - температура изделия, Т = Тн + Т; радиус нагревателя Я = 4 мм; мощность нагревателя Ж = 0,6 Вт; временной интервал измерения температуры Дх = 0,25 с; расстояния от центра нагревателя т\ = 7 мм и г2 = 9 мм. В общем случае на каждой из термограмм можно выделить несколько участков, соответствующих различным состояниям температурного поля исследуемой системы. Так, для термограммы, зафиксированной центральным ТП на изделии из ПМ вне области структурных переходов, характерны восемь участков (рис. 3).
Первому участку соответствует одномерное температурное поле в локальной области исследуемого тела (вблизи нагревателя). Тепловые потоки, поступающие в изделие д\ и зонд q2, изменяются во времени, так как между нагревателем и исследуемым телом имеется термическое сопротивление, нагреватель обладает инерционностью. Второму участку отвечает одномерное температурное поле, но процесс проходит стадию регуляризации в локальной области исследуемого тела, расположенной вблизи нагревателя. Третьему участку соответствует двухмерное температурное поле в образце, поскольку нельзя пренебречь распространением тепла в радиальном направлении. Четвертому участку соответствует тепловой процесс, вышедший на стадию регуляризации. В локальной области исследуемого ПМ формируется температурное поле, близкое к одномерному полусферическому. На пятом участке нарушаются условия полуограниченности исследуемого тела.
После отключения нагревателя (т > тоткл), на стадии остывания, можно выделить шестой участок, когда тепловые потоки ql и q2 изменяются во времени, седьмой участок термограммы, тепловой процесс в котором проходит стадию регуляризации, и восьмой участок, где тепловой процесс изменяется.
Участки термограммы II, IV и VII - рабочие [1-3], так как вне области структурного перехода в ПМ возможно однозначно определить значения ТФС в зависимости от параметров аналитических моделей, описывающих термограмму на данных температурно-временных интервалах, используя регулярные тепловые режимы на моделях плоского и сферического полупространств [4]. Авторами разработаны методики определения границ рабочих участков с использованием статистических критериев [1].
Известные решения краевых задач нестационарной теплопроводности, описывающие процесс распространения тепла в твердых телах с различными ТФС [4-6] от плоского круглого нагревателя постоянной мощности радиуса Я, имеют
T
200 -
150 -
50-
\
200
600
1000
Рис. 2. Термограммы для ПТФЭ:
1 - в центре зонда; 2, 3 - на расстояниях 7 и 9 мм от центра
Рис. 3. Термограмма 1 (рис. 2) с выделенными участками
весьма сложный вид и малопригодны для их использования при реализации метода НК. Однако, учитывая, что распределение температуры от плоского круглого источника тепла постоянной мощности радиуса Я при малых значениях времени т близко к распределению температуры в плоском полупространстве [1-4], для получения математической модели, описывающей процесс распространения тепла в исследуемой системе (изделие-зонд) на втором участке термограммы, измерительная схема с плоским круглым нагревателем постоянной мощности радиусом Я заменена схемой с бесконечным плоским нагревателем [1-3]. Это позволило получить расчетное соотношение, описывающее термограмму на втором (рабочем) участке
Ti(0, т) =
2дл/т
дсн
(ei +62Ь/П (si +S2)2
(1)
где 7^1, 61, е2 - избыточная температура, тепловые активности материалов исследуемого изделия и подложки зонда, соответственно; q - удельная поверхностная мощность; т - время; сн - теплоемкость единицы площади нагревателя. Индексы 1, 2 отнесены к материалам изделия и подложки ИЗ соответственно.
Для четвертого и шестого (рабочих) участков термограммы расчетные соотношения приведены в работах [1-3].
Численным моделированием температурных полей методом конечных элементов с помощью пакета БЬСЦГ [7], определена адекватность полученной математической модели (1) реальному тепловому процессу на втором (рабочем) участке термограммы [9].
В данной работе авторами представлены результаты ряда вычислительных экспериментов, по которым получено уравнение регрессии для определения температуры середины второго (рабочего) участка термограммы Тр в зависимости от тепловой активности 61 полимерного материала и удельной поверхностной мощности q для измерительного зонда с подложкой из рипора.
Данная зависимость имеет следующий вид
Tp = 0,003q + 0,018s1 - 26,24 .
(2)
150 < б! < 1550 Вт-с0,5/(м2-К);
2 -з
1000 < q < 29000 Вт/м . Радиус нагревателя Я = 4-10 м. Результаты регрессион-
Диапазоны варьирования переменных: т/м2. Радиус нагре ного анализа представлены на рис. 4.
100
0
0
т. с
т, с
2
Для значения множественного коэффициента детерминации Яд = 0,999, свидетельствующего о том, что вариация исследуемой переменной Тр (температуры середины рабочего участка) на 99,9 % объясняется изменчивостью включенных в модель переменных, проверена значимость уравнения регрессии (2). Если Е (фактическое значение критерия Фишера) превышает некоторое критическое значение выбранных данных п, к (п - размер выборки; к - число коэффициентов модели) и принятой доверительной вероятности, то величина Яд считается существенной
Е П-к = °,999 15 - 2 = 12987
1 _ Я2 к -1 1 - 0,999 2 -1
1 Дд
Фактическое значение Е больше табличного Е0 05; 2; 12 = 3,89, определенного на уровне значимости а = 0,05 при к\ = к = 2 и к2 = п - к - 1 = 15 - 2 - 1 = 12 степенях свободы [8]. Таким образом, уравнение регрессии (2) значимо. Следовательно, исследуемая зависимая переменная Тр достаточно хорошо описывается включенными в регрессионную модель переменными q и б!.
На рис. 5 а, б представлены зависимости Тр = /(Я), Тр = /^), полученные
численным моделированием.
По уравнению регрессии (2) и данным имитационного исследования, представленным на рис. 5, определяют температуру Тр в области второго участка
термограммы, в котором возможна регистрация структурных переходов в полимерных материалах по отклонениям зависимостей тепловой активности с ростом температуры.
Пример реализации метода неразрушающего контроля температурных характеристик структурных переходов в полимерах представлен на рис. 6. Использованы термограммы, зафиксированные центральным ТП после калибровки ИЗ на материалах с известными ТФС [10]. Условия проведения опытов, выбранные с учетом выполненного имитационного исследования, следующие: Я = 4 мм; Ж = 0,7 Вт; Дх = 0,5 с; Тн = 21 °С. Исследовали твердофазные переходы в политетрафторэтилене (см. рис. 6, а) и капролоне (см. рис. 6, б).
60
40
20
60
40
20
3 4 5 6 Я-10 " 0 10000 20000 30000 q
а) б)
Рис. 5. Зависимости Гр = /(Я) и Гр = /(</). Условия моделирования:
а) е1 = 758,1; q = 30250 (1); q = 20500 (2); q = 10750 (3); q = 1000 (4); б) е1 = 758,1; Я = 6,5 (1); Я = 5 (2); Я = 4,5 (3); Я = 4 (4); Я = 3,5 (5); Я = 3 (6); Я = 2,5 (7); Я = 2 (8)
e*, Вт-с0,5/(м2 К)
e*, Вт-с0,5/(м2 К)
1 500
500
26 28 30 32 Ts, °С 24 25 26
а) б)
*
Рис. 6. Зависимость £ = / (Г ) для изделия из ПТФЭ (а) и капролона (б)
Ts, °C
T.
0
2000-
0
0
Известно, что в ПТФЭ при температуре ниже 19,6 °С элементарная ячейка ПТФЭ имеет триклиническую структуру. В интервале от 19,6 до 30 °С существует гексагональная элементарная ячейка, выше 30 °С стабильной становится псевдогексагональная решетка. Удельные теплоты переходов составляют соответственно: 4,0 ± 0,5 и 1,2 ± 0,3 кДж/кг. Твердофазный переход в ПТФЭ, происходящий при ^ = 30 °С и сопровождающийся эндотермическим эффектом, зарегистрирован при ^ = 29,8 °С (см. рис. 6, а). Твердофазный переход в капролоне, соответствующий переходу при 25 °С гексагональной решетки в моноклинную, проявился в интервале температур 24,7...25,3 °С (см. рис. 6, б).
Твердофазные переходы в политетрафторэтилене и капролоне проявились в температурных интервалах, соответствующих вторым участкам на термограммах.
Полученные результаты хорошо согласуются со справочными данными.
Таким образом, найденное уравнение регрессии и результаты имитационного исследования позволили обнаружить структурные переходы в ПТФЭ и капролоне по отклонениям на зависимостях тепловой активности от температуры.
Список литературы
1. Пат. 2167412 Российская Федерация, МПК в 01 N 25/18. Способ комплексного определения теплофизических свойств материалов / Н.П. Жуков, Н.Ф. Май-никова, Ю.Л. Муромцев, И.В. Рогов ; заявитель и патентообладатель Тамб. гос. техн. ун-т. - № 99103718/28 ; заявл. 22.02.99 ; опубл. 20.05.2001, Бюл. № 14. - 8 с.
2. Об одном методе термического анализа для неразрушающего контроля теплофизических свойств полимеров / Н.Ф. Майникова [и др.] // Пластические массы. - 2001. - № 2. - С.30-33.
3. Об одном методе исследования твердофазных переходов в полимерах / Н.Ф. Майникова [и др.] // Пластические массы. - 2002. - № 6. - С. 23-26.
4. Платунов, Е.С. Теплофизические измерения и приборы / Е.С. Платунов. -Л. : Машиностроение, 1986. - 256 с.
5. Лыков, А.В. Теория теплопроводности / А.В. Лыков. - М. : Высшая школа, 1967. - 599 с.
6. Козлов, В.П. Двумерные осесимметричные нестационарные задачи теплопроводности / В.П. Козлов. - Минск : Наука и техника, 1986. - 392 с.
7. БЬСиТ: Моделирование двумерных полей методом конечных элементов. Версия 5.1. Руководство пользователя. - СПб. : Производственный кооператив ТОР, 2003. - 249 с.
8. Химмельблау, Д. Анализ процессов статистическими методами / Д. Хим-мельблау - М. : Мир, 1973. - 960 с.
9. Определение условий адекватности модели распределения тепла в плоском полупространстве реальному процессу при теплофизическом контроле / Н.П. Жуков [и др.] // Вест. Тамб. гос. техн. ун-та. - 2006. - Т. 12, № 3А. -С. 610-616.
10. Жуков, Н.П. Многомодельные методы и средства неразрушающего контроля теплофизических свойств материалов и изделий / Н.П. Жуков, Н.Ф. Майникова. - М. : Машиностроение-1, 2004. - 288 с.
Simulated Research into Heat Transfer in the Method of Nondestructive Control over Structural Transformations of Polymers
N.F. Mainikova1, N.P. Zhukov2, A.A. Balashov2, D.G. Borodavkin2
Departments: "Theory of Machines, Mechanisms and Machine Parts" (1), "Hydraulics and Heat Engineering" (2), TSTU
Key words and phrases: flat semi-space; non-destructive control; numerical modeling; structural transformation; thermo-physical properties.
Abstract: Simulated research into the method of non-destructive control over structural transformations in polymer materials is carried out.
Imitationsforschung der Wärmeübertragung in der Methode der nicht zerstörenden Kontrolle der strukturellen Übergänge in den Polymeren
Zusammenfassung: Es ist die Imitationsforschung der Methode der nicht zerstörenden Kontrolle der strukturellen Übergänge in den Polymerstoffen durchgeführt.
Etude d'imitation du transfert de la chaleur dans la métodes du contrôle non-destructif des passages structuraux dans les polymères
Résumé: Est réalisée une étude d'imitation du transfert de la chaleur dans la métode du contrôle non-destructif des passages structuraux dans les polymères.
