Численное исследование температурных неоднородностей при нагреве металлов Текст научной статьи по специальности «Физика»

УДК 536.2

ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ТЕМПЕРАТУРНЫХ НЕОДНОРОДНОСТЕЙ ПРИ НАГРЕВЕ МЕТАЛЛОВ

И.А. КОНЦЕВОЙ

Учреждение образования «Гомельский государственный технический университет имени П. О. Сухого»,

Республика Беларусь

Введение

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

В работах [2-4] для рассмотрения динамического теплового гистерезиса (ДТГ) применялась классическая модель теплопроводности Фурье. Получены сведения о том, что в плоскости «тепловой поток - градиент температуры» необходимым образом существует петля гистерезиса, если: 1) интенсивность поверхностного источника энергии (^) немонотонно зависит от времени; 2) тепловой процесс происходит в

таком интервале температур, что эффективным образом проявляет себя температурная зависимость Я(Т) коэффициента теплопроводности. Данная работа является продолжением исследований [2-4]. Цель работы: 1) изучить нестационарные свойства теплового поля металлического образца при воздействии на него высокоинтенсивного поверхностного источника энергии; 2) исследовать влияние параметров поверхностного источника энергии и характера зависимости Я (Т) металлов на качественные и количественные особенности поведения петель ДТГ в плоскости «тепловой поток - градиент температуры».

Постановка задачи

Металлическая пластина толщиной И имеет в начальный момент времени I = 0 температуру Т = Т0; при ^ > 0 правая граница х = И изотермическая Т(ЪЛ) = Тк = Т0. На левую границу пластины х = 0 действует тепловой поток ^(0,1) = ц0 (I) , моделирующий воздействие концентрированного потока энергии. Нагрев осуществляется в широком интервале температур, на котором нелинейные теплофизические свойства металла проявляют себя в полной мере (фазовые превращения не происходят).

Теплофизические свойства металла с(Т), Я(Т) на интервале температур [Т1, Т2 ] описываются полиномами третьей степени с постоянными коэффициентами. Построение этих полиномов осуществляется на основе справочных данных [5].

Математическая модель поставленной задачи приведена, например, в [2-4]. Рассматриваются поверхностные импульсы энергии треугольной, колоколообразной и трапециевидной временных форм. Все расчеты проводятся в безразмерных величинах. Приведение размерных величин к безразмерным приведено в [2].

Метод решения

Решение поставленной задачи выполнено численным методом интегральных соотношений А.А. Дородницына [6]. Область интегрирования делится на пять полос: хг = ¡И/5, г = 0,1,..., 5. Алгоритм решения позволяет учитывать поглощательную

способность металла А(Т) = А0 + А1Т. Подробности построения схемы вычислений

даны в [2, 3].

Производство энтропии подсчитываем по формуле

ст = « Тх (Т ] •

Для характеристики нестационарных и нелинейных свойств теплопереноса применяем следующие критерии:

1 дТ)

\д Х у

к Ь ]=

1 + (р (1)[ Ь ])

3/2

Р |1) Ь ] = Ь(Х2) - Ь1 Х0) Р (2) Ь ] = Ь1 Х2) - 2 Ь1 Х1) + Ь1 Х0)

Х2 - Х0 (Х1 - Х0)2

Ь (Хі) = Ь (Т)/Ь (Т0), Тг = Т (Х, і);

О а( Хг,І) _0 Л(Т °)

0 = а0 , ^ = ~И^ •

а и

Здесь Q - безразмерный тепловой поток, g - безразмерный градиент температуры, О - безразмерное производство энтропии. В качестве /(Т) применяем теплофизические параметры Я(Т) и а(Т) = Я(Т) /с(Т) . Параметр кривизны К[] относится к левой границе области. Параметры Р(1) Ь ] , Р(2) Ь ] представляют собой конечно-разностные формулы для первой и второй производных изучаемых функций.

Результаты расчетов

Расчеты выполнены для нескольких металлов (железо, молибден, вольфрам, никель, ванадий) с различным характером температурных зависимостей теплофизических свойств Я(Т), с(Т). В качестве граничного теплового потока q0(і) рассмотрены импульсы различных временных форм (рис. 1).

Основными характеристиками этих импульсов являются:

- максимальный тепловой поток qmax;

- время і = т0 достижения qтах= q0(т0);

- длительность импульса т .

а) б) в)

Рис. 1. Импульсы треугольной (а), колоколообразной (б) и трапециевидной (в) временных форм

Представим результаты моделирования теплового воздействия на пластины из никеля Ni и ванадия V, теплофизические свойства которых на интервале [, T2] аппроксимируются следующими полиномами:

Ni: T, K е [300, 630],

c(T) = -6703290 + 75309,2T-173,51456 T2 + 0,135664 T3 Дж/(м3 • К);

A(T) = 223,65-0,8425T + 0,00169T2 -1,2217-10-6T3 Вт/(м• К);

V: T, K е [300, 2000],

c(T) = 2689498 + 996,3834 T - 0,33259 T2 + 0,0001506 T3 Дж/(м3 • К);

A(T) = 29,11 - 0,000264 T +1,39734-10-5 T2 - 4,23507-10-9T3 Вт/(м • К).

В основной серии расчетов (h = 5 • 10-3 м, Tw = T0 = 300 К ) для каждого вида q0 (t) рассмотрены три различные величины т :

I. Треугольный импульс (т0 = 0,125 т)

а) т = 0,5 т* = 1,6 с; б) т = 0,75 т* = 2,4 с; в) т = т* = 3,2 с;

II. Колоколообразный импульс (т0 = 0,5 т)

а) т = 0,25 т* = 0,25 с; б) т = 0,5 т* = 0,5 с; в) т = т* = 1,0 с;

III. Трапециевидный импульс (т0 = т- т1 = 0,125 т)

а) т = 0,5 т = 1,6 с; б) т = 0,75 т* = 2,4 с; в) т = т* = 3,2 с .

Для вариантов а-в всех видов импульса рассмотрены три величины qmax:

1 2

1) q max = 3 q* ; 2) qmax = 3 q ; 3) q max = q.

Величина q* для каждого вида импульса подобрана такой, что при т = т* температура T(0, t) левой границы пластины в процессе нагрева изменяется во всем интервале [, T2], т. е. температурные зависимости теплофизических свойств металла A(T), c(T) проявляют себя в полной мере. Для импульсов треугольной, колоколо-

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

Вт Вт Вт

N1: 6,411-Ю6 — ; 9 -106 —; 4,968 -106 —;

мм м

Вт Вт Вт

V: 20,71-106 — ; 31,72-106 —; 14,43-106 —.

м м м

Ниже приведены характерные гистерезисные зависимости при треугольной форме импульса (т = т*, цтах = ц ) для левой границы (а - материал пластины Мт, б -материал пластины V ): между тепловым потоком и градиентом температуры (рис. 2), между относительной температурой Т0 / Т0 и параметром кривизны К[Х\ (рис. 3), между производством энтропии и параметром кривизны К[Х\ (рис. 4).

а) б)

Рис. 2. Формирование петли гистерезиса в плоскости ^, £)

Рис. 3. Гистерезисные зависимости в плоскости (Т0 / Т0, К0[А\)

ВД

а)

Рис. 4. Гистерезисные зависимости в плоскости (К0[А\, 50)

б)

Как видно из рис. 2, формирование петли гистерезиса в плоскости (0, £■) для никеля и ванадия происходит различным образом. При подводе тепла (ц0 / & > 0) к пластине никеля формируется нижняя часть петли, а при отводе (&Ц0 /& < 0) - верхняя (начальное состояние соответствует началу координат). Для ванадия - наоборот: при &Ц0 /& > 0 формируется верхняя часть петли, а при &ц0 / & < 0 - нижняя. Такой

результат обусловлен тем, что на всем интервале температур [Т, Т2 ] зависимость Л(Т) никеля есть функция монотонно убывающая, а для ванадия функция Л(Т) монотонно возрастающая.

Мерой неоднозначной зависимости между градиентом температуры и тепловым потоком является относительное удлинение в (относительная площадь) петли ДТГ:

в=К. . V.,

Ъ0 Ъ0Ъ1

где Ъ - «длина петли», Ь0 - «ширина петли» [4]. Во всех рассмотренных вариантах форма лепестка такая, что К1 /Ъ0 > 1. Чем меньше в, тем сильнее выражен ДТГ. Результаты основной серии расчетов представлены в таблице 1.

Таблица 1

Влияние формы, длительности и максимальной интенсивности импульса на параметр в

Длительность импульса т , с Материал пластины

Никель Ванадий

Максимальная величина теплового потока д тах

д*/3 2д /3 * д д*/3 2д* /3 * д

Параметр в Параметр в

Треугольная форма импульса

1,6 39,27 24,52 15,46 21,64 11,46 9,21

2,4 37,33 22,06 11,33 17,90 10,44 9,02

3,2 36,55 21,16 8,77 16,09 9,79 9,20

Окончание табл. 1

Длительность импульса т , с Материал пластины

Никель Ванадий

Максимальная величина теплового потока д тах

д */3 2д */3 * д д */3 2д */3 * д

Параметр в Параметр в

Треугольная форма импульса

Колоколообразная форма импульса

0,25 75,89 46,07 35,77 45,29 19,70 14,02

0,5 62,41 37,68 21,39 29,80 13,71 10,95

1,0 47,66 28,24 11,01 20,75 10,84 10,33

Трапециевидная форма импульса

1,6 30,19 18,49 9,90 14,05 8,06 6,97

2,4 32,27 18,74 8,77 13,00 7,92 7,20

3,2 34,63 19,99 8,77 12,79 8,03 7,59

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

Неоднозначные зависимости между тепловым потоком и градиентом температуры, между относительной температурой Т0 /Т0 и параметром кривизны К[Л\, между производством энтропии и параметром кривизны К[Х] (см. рис. 2-4) обусловлены двумя факторами: 1) нелинейной зависимостью Л(Т) коэффициента теплопроводности; 2) немонотонным характером зависимости источника энергии ц0(Т) от времени.

Выводы

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

Импульсы энергии треугольной, колоколообразной и трапециевидной временных форм инициируют многовариантные эволюционные процессы в металлах. Свойства этих процессов зависят от функциональных зависимостей Л(Т), с(Т), а также от параметров импульса.

Работа выполнена под научным руководством профессора О.Н. Шабловского.

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

1. Дьюли У. Лазерная технология и анализ материалов. - М.: Мир, 1986. - 502 с.

2. Шабловский О.Н., Кроль Д.Г., Концевой И. А. Импульсный нагрев металла в широком интервале температур //Машиностроение. - Минск, 2002. - Вып. 18. -С. 516-520.

3. Шабловский О.Н., Кроль Д.Г., Концевой И. А. Нестационарные свойства поверхностного нагрева металлов //Вістник Запорізького держ. університету. Фіз.-мат. науки. - 2002. - № 1. - С. 148-152.

4. Шабловский О.Н., Концевой И.А. Гистерезисные свойства нестационарного температурного поля //Математические модели в образовании, науке и промышленности: Сб. науч. трудов. - С.-Пб.: Санкт-Петербургское отделение МАН ВШ, 2003. - С. 262-266.

5. Зиновьев В.Е. Теплофизические свойства металлов при высоких температурах: Справ. изд. - М.: Металлургия, 1989. - 384 с.

6. Белоцерковский О.М., Грудницкий В.Г. Исследование нестационарных течений газа со сложной внутренней структурой методами интегральных соотношений //Журнал вычисл. матем. и матем. физики. - 1980. - Т. 20. - № 6. - С. 1400-1415.

Получено 28.11.2003 г.