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

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

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

А.И. Данилушкин, В.А. Данилушкин, Е.А. Никитина

Самарский государственный технический университет

443100, г. Самара, ул. Молодогвардейская, 244

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

Ключевые слова: индуктор, внутренние источники тепла, температурное распределение.

Рассматривается задача идентификации процесса индукционного нагрева цилиндрического ферромагнитного изделия во вращающемся магнитном поле, создаваемом обмоткой трехфазного индуктора. Индуктор выполнен в форме статора асинхронного двигателя. Реализация такой нестандартной конструкции индукционной установки требует решения ряда задач, связанных с исследованием взаимосвязанных электромагнитных и тепловых полей системы. Сложность моделирования электротепловых процессов обусловлена неоднородной составной структурой системы «индуктор - футеровка - изделие», различными условиями теплообмена между сопряженными телами различной физической природы, наличием внутренних источников тепла, неравномерно распределенных по объему нагреваемого изделия. Моделирование процесса сводится к решению двух взаимосвязанных задач - электромагнитной и тепловой.

В работе [1] приведены методика расчета и результаты исследования электромагнитных полей системы «трехфазный индуктор - цилиндрическая заготовка» и определен характер распределения внутренних источников тепла.

В настоящей работе исследуется температурное распределение в системе «индуктор - заготовка» при нагреве внутренними источниками тепла во вращающемся магнитном поле трехфазного индуктора.

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

Эскиз индукционного нагревателя представлен на рис.1.

Рис. 1. Геометрическая модель исследуемой установки: 1 - магнитопровод индуктора; 2 - воздушный зазор; 3 - футеровка; 4 - изоляция; 5 - нагреваемое изделие; 6 - медная трубка специального профиля

индуктор заготовка трехфазный теплообмен

Исходная математическая модель нестационарной теплопроводности представлена системой дифференциальных уравнений в частных производных вида [2]:

; (1)

, (2)

, (3)

соответственно для цилиндрической заготовки (1), футеровки (2) и магнитной системы индуктора (3).

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

, (4)

на торцевых поверхностях заготовки:

, (5)

на торцевых поверхностях футеровки:

, (6)

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

(7)

. (8)

Граничные условия на внешней поверхности магнитопровода:

; (9)

на торцевых поверхностях:

. (10)

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

При исследовании температурных полей принимается ряд допущений.

1. Контакт между футеровкой и магнитопроводом индуктора можно считать идеальным.

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

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

По результатам решения электромагнитной задачи [1] получены функции распределения внутренних источников тепла в нагреваемом цилиндре вида

, (11)

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

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

Следуя МКЭ, дифференциальному уравнению (1) ставится в соответствие вариационная формулировка о минимизации энергетического функционала, характеризующего тепловое состояние тела:

(12)

где - граница с конвективным теплообменом; .

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

При рассмотрении теплообмена с внешней средой по двум граничащим со средой сторонам принимается, что имеет место общий случай граничных условий вида (4).

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

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

Ниже приведены некоторые результаты исследования процесса индукционного нагрева стальной ферромагнитной цилиндрической заготовки в трехфазном индукторе с вращающимся магнитным полем. Параметры индукционной системы: диаметр заготовки - 140 мм, толщина футеровки (шамот группы В) - 20 мм, величина воздушного зазора между заготовкой и индуктором - 5 мм. Число пазов индуктора - 6. Обмотка индуктора выполнена из стандартной медной трубки специального профиля со смещенным отверстием. Индуктор охлаждается водой. Источником тепла служит объемная плотность тепловыделения, обусловленная вихревыми токами в заготовке. График распределения поверхностной плотности тока по внешнему радиусу заготовки, полученный по результатам решения электромагнитной задачи, представлен на рис. 2.

Рис. 2. График распределения мгновенного значения плотности полного тока по внешнему радиусу заготовки

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

Распределение температурного поля футеровки (рис. 3) показывает незначительный перепад температуры между областями под трубкой и под зубцом магнитопровода. Большее значение температуры под трубкой объясняется тепловым потоком с граничащей с футеровкой поверхности трубки, который обусловлен джоулевыми потерями в индуктирующем проводе. Максимум температуры магнитопровода наблюдается в точке, находящейся на стыке футеровки, магнитопровода и трубки, и составляет 66 єС.

Неравномерность температурного распределения по окружности заготовки наблюдается только на расстоянии от поверхности, соответствующем глубине проникновения тока в металл. С приближением к центру изотермы принимают вид концентрических окружностей с центром, совпадающим с осью симметрии заготовки. Максимальный перепад температур по окружности заготовки, представленный на рис. 4, составляет 25 К.

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

Рис. 3. Температурное поле футеровки на участках под пазом и под трубкой

Рис. 4. Распределение температуры по контуру поверхности заготовки

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

Библиографический список

1. Базаров А.А., Данилушкин А.И., Никитина Е.А. Моделирование и расчет внутренних источников тепла в трехфазном индукторе c вращающимся магнитным полем // Вестник СамГТУ. Сер. Технические науки, 2009. - №2 (24). - С. 120-127.

2. Лыков А.В. Тепломассообмен (Справочник). - М.: Энергия, 1978. - 480 с.

Размещено на Allbest.ru