THE SOLUTION OF FIELD PROBLEMS AT DESIGNING INDUCTOR WITH A MAGNETIC CIRCUIT

Full Text

Системы индукционного нагрева применяются довольно длительное время, с большим успехом составляя конкуренцию другим видам нагрева. Вместе с тем такие особенности технологического процесса, как высокая температура и повышенная частота питающего напряжения, существенно ограничивают применение магнитопроводов при изготовлении индукторов. Это сказывается отрицательно на энергетических характеристиках устройств индукционного нагрева, главным образом за счет больших полей рассеяния магнитного поля. В последнее время несколько расширилось применение магнитопроводов из новых материалов, таких как магнитодиэлектрики, что положительно отражается на управлении распределением магнитного поля и внутренних источников тепла в загрузке. Предлагаемая работа направлена на расширение возможностей индукционных устройств за счет использования замкнутых магнитопроводов. Конечной целью такого решения является повышение управляемости процесса формирования внутренних источников. На современном этапе ставится задача разработки математических моделей и алгоритмов проектирования процессов и устройств. Особенность рассматриваемых задач заключается в сложности математического аппарата и значительной ресурсоемкости вычислений. Ориентация на использование персонального компьютера предполагает оценку возможностей моделирования процессов, погрешностей вычислений и получения приемлемых в проектной практике результатов. Конструкция индуктора с магнитопроводом (рис. 1) содержит магнитопровод, аналогичный статору двигателя переменного тока, с вложенными в пазы изолированными витками обмотки, футеровку и загрузку. Работы в данном направлении были посвящены исследованию отдельных вопросов, связанных с разработкой конструкции магнитопровода, выбором систем охлаждения проводников и железа [1]. В предлагаемой работе ставится задача более детального исследования характеристик. Для простоты рассматривается обмотка с комбинацией прямого и обратного чередования фаз, чтобы исключить проявление вращающего момента на загрузку. Исследование характеристик индуктора с немагнитной загрузкой также позволяет не опасаться притяжения загрузки к футеровке за счет электромагнитных сил с возможностью ее быстрого разрушения ввиду ее хрупкости и низкой стойкости к истиранию. Рис. 1. Трехфазный индуктор с магнитопроводом: 1 - магнитопровод; 2 - водохлаждаемый виток обмотки; 3 - изоляция витка; 4 - футеровка; 5 - загрузка; 6 - воздушный зазор Проектирование индуктора включает определение не только конструктивных параметров, таких как внутренний и внешний диаметр магнитопровода, длина индуктора, толщина футеровки, но и мощности индуктора, тока, напряжения, коэффициента мощности. Кроме того, необходимо согласовать функционирование оборудования не только с другими технологическими узлами, например прессом, но и с питающей сетью. Такие режимы, как пусковой, могут серьезно усложнить работу. Анализ всех вопросов одновременно невозможен, требуется многократное уточнение параметров при более детальном рассмотрении различных аспектов. Сходная по сложности картина наблюдается при исследовании влияния конструктивных параметров на энергетические характеристики индуктора. В качестве исходных данных взяты физические свойства для алюминиевого сплава марки Д16. Размеры слитка: диаметр 0,42 м, длина 0,8 м. Температура максимально нагретой точки - 480 , температура центра - 430 . Допустимое отклонение температуры - ±12 . Физические свойства: плотность 2700 ; температура плавления 660 °С; коэффициент теплопроводности ; коэффициент теплоемкости ; удельное сопротивление . В отличие от электрических машин в рассматриваемой задаче нет жестких ограничений на величину допустимой индукции магнитопровода. В то же время отсутствие магнитных свойств у загрузки делает непростой процедуру выбора ширины паза и зубца магнитопровода. Решение этой задачи достигается построением зависимостей от величины соотношения размеров паза и зубца. Для выполнения этой процедуры вполне возможно обойтись двумерной электромагнитной задачей. Перейдя от системы уравнений Максвелла [2] к формулировке с использованием векторного магнитного потенциала, можно записать дифференциальное уравнение с соответствующими граничными условиями: ; . (1) Здесь: A - векторный магнитный потенциал; - относительная магнитная проницаемость материала по осям x, y; x, y - координаты; j - мнимая единица; - круговая частота тока; g - удельная электрическая проводимость; L - граница расчетной области; - плотность стороннего тока. На базе уравнений (1) строится конечно-элементная формулировка [3] для плоской двумерной области. В качестве программной реализации используется пакет Elcut. В версиях начиная с 5.6 реализована поддержка нелинейной электромагнитной задачи, что позволяет корректно решать подобные задачи. Конечная система уравнений имеет вид . (2) Здесь - матрица жесткости и векторы источников и неизвестных величин. При решении задачи моделирования электромагнитных процессов в системе «индуктор - загрузка» были приняты: число пазов, равное 12; суммарная ширина паза и зубца 123 мм; диаметр индуктора внутренний 470 мм; диаметр индуктора внешний 640 мм; диаметр загрузки 420 мм. В результате расчетов были получены характеристики индукторов для пяти вариантов соотношения размеров паза и зубца (табл. 1). Макcимальное значение индукции в магнитопроводе находится в пределах 1,19 - 1,33 Тл. Полученные результаты отражают не только зависимость различных величин от соотношения размеров паза и зубца, но и в некоторых пределах явное влияние соотношения размеров паза и зазора между индуктором и загрузкой. На втором этапе рассмотрим влияние длины индуктора на его характеристики при неизменной длине загрузки. Таблица 1 Параметры индуктора при различной ширине паза № 1 2 3 4 5 Ширина паза, мм 60 65 70 75 80 Ширина зубца, мм 63 58 53 48 43 Фазный ток, А 34 700 35 300 35 800 36 900 36 450 Фазное напряжение, В 30,6 28,6 27,6 27 27,2 Мощность в загрузке, кВт 152 150 149,8 150,8 149,8 Потери мощности в индукторе, кВт 178 168 160 159 236 Мощность индуктора, кВт 330 318 309,8 309,8 385,8 Полная мощность, кВА 3 185 3 028 2 964 2 989 2 974 КПД 0,46 0,47 0,483 0,487 0,388 Cos ϕ 0,103 0,105 0,105 0,104 0,13 Для классической конструкции индуктора, представляющей собой соленоидальную катушку, имеют место краевые эффекты. По краям загрузки в зависимости от заглубления, или, другими словами, от соотношения длин индуктора и загрузки может наблюдаться как значительное снижение удельной мощности, так и повышение [4]. Регулируя величину заглубления, можно добиться желаемого уровня повышения мощности на краях, чтобы скомпенсировать тепловые потери на торцах, где не предусмотрена тепловая изоляция. Решение задачи моделирования электромагнитных, а также еще и тепловых процессов требует использования трехмерной постановки. Это существенно повышает требования к ресурсам компьютера. Если в двумерной постановке электромагнитная задача может содержать 2 млн элементов, то это не будет чрезмерной нагрузкой для компьютера. В случае трехмерных моделей соотношение числа кончных элементов и оперативной памяти совершенно другое. Так, для компьютера с оперативной памятью в 32 Гб для решения задачи с 100 000 элементами потребуются все ресурсы. Кроме того, подобная задача требует времени около 2 суток. Для описания электромагнитного поля в трехмерной постановке требуется использовать не только векторный магнитный потенциал , но и электрический скалярный потенциал , чтобы обеспечить замыкание магнитных потоков и токов при изменении направления движения. (3) Здесь: - оператор набла, где - единичные векторы по осям x, y, z соответственно. Оператор набла при применении к вектору записывается в виде векторного произведения , а при применении к скалярной величине - в виде скалярного произведения . Расчетная система уравнений составляется аналогично (2), но с отличием в векторе неизвестных величин: ; (4) . (5) Как видно, размерность вектора существенно больше, чем у аналогичного в выражении (2). Конечно, изменения в размерах сказываются не только на векторах, но и на матрице жесткости. Как отмечалось выше, значительно возрастающие требования к вычислительным ресурсам заставляют корректировать подход к построению модели. Вместо рассмотрения всего объекта приходится ограничиться его фрагментом. В данном случае с некоторыми допущениями можно произвести расчет для сектора (рис. 2). Граничные условия в этом случае приходится задавать не только на боковой поверхности вокруг индуктора и на торцах, но и на плоскостях среза. Из имеющихся вариантов наиболее подходящими являются условия антипериодичности: (6) Здесь индексы dst и src означают «назначение» и «источник», фактически отражая значения функций слева и справа от поверхности раздела. Условия магнитной изоляции задаются на боковой поверхности и на торцах (7) Более правильным вместо (6) было бы задание магнитного потока с фазовым сдвигом, но такое в программе не предусмотрено. При анализе расчетной модели (см. рис. 2) напрашивается аналогия с линейным двигателем. При кажущейся равноценности этот вариант намного предпочтительнее. Дело в том, что формируемые матрицы жесткости имеют отличие: для конструкции в виде сектора они являются несимметричными. Это приводит к тому, что метод решения является более трудоемким и ресурсоемким. На примере рассматриваемой задачи прямоугольная модель может содержать вместо 100 000 элементов 300 000. Кроме того, скорость решения несравнимо выше. Таким образом, при достаточном анализе результатов решения можно существенно сэкономить ресурсы при повышении точности расчетов за счет изменения модели. Моделирование тепловых процессов само по себе не столь ресурсоемко. Однако при увеличении векторов источников и матриц жесткости аналогично выражениям (5, 6) при объединении электромагнитных и тепловых задач такая процедура существенно усложняет расчет. Гораздо экономичнее и быстрее выявить траекторию протекания тока и определить мощности тепловыделения, чтобы затем перенести в тепловую задачу. Конечно, это невыгодно делать при одноразовом решении или при двумерной постановке. Рис. 2. Сектор индуктора с двумя витками: 1 - магнитопровод; 2 - паз; 3 - виток; 4 - загрузка; 5 - воздух Уравнение теплопроводности для трехмерной области с внутренними источниками тепла имеет вид . (8) Здесь: - плотность материала; - коэффициенты теплопроводности по осям x, y, z; Т - температура; w - удельная мощность тепловыделения. В качестве граничных условий на боковой поверхности и торцах используется конвективный теплообмен в форме . (9) Здесь индекс x модет принимать значения других координатных осей или нормали n по отношению к рассматриваемой поверхности. На поверхностях среза задается условие тепловой изоляции . (10) Таким образом, сформирован комплекс математических моделей для трехмерной постановки, позволяющий рассчитывать электромагнитные и тепловые процессы при исследовании свойств индукционного нагревателя с магнитопроводом. Для определения влияния соотношения длин загрузки и индуктора был произведен ряд расчетов. Результаты представлены в табл. 2 и на рис. 3. Нужно иметь в виду, что края витков индуктора находятся вне пределов магнитопровода (по аналогии с конструкцией электрических машин - лобовые части витков). Распределение мощности для загрузки длиной 0,9 и 1 м совпадают. Отличие в том, что для загрузки длиной 0,9 м кривая начинается на расстоянии 0,05 м от края. Таблица 2 Параметры индуктора при разной длине загрузки № Длина загрузки, м Ток индуктора, кА Мощность в загрузке, кВт Удельная мощность, Вт/м3, 1 0,6 43,2 9,9 1,35 2 0,7 43,1 16,8 2 3 0,8 43,2 24,8 2,5 4 0,85 43,2 28,7 2,8 5 0,9 43,16 39,4 4 6 1 43,2 38,9 4 Рис. 3 Диаграммы распределения удельной мощности в загрузке под витком индуктора: 1 - при длине загрузки 0,6 м; 2 - при длине загрузки 0,7 м; 3 - при длине загрузки 0,8 м; 4 - при длине загрузки 0,85 м; 5 - при длине загрузки 1 м Как и следовало ожидать, мощность в загрузке за пределами магнитопровода быстро снижается до нуля. Поэтому нагрев такой загрузки будет неэффективным, хотя полезная мощность возрастает до значения длины 0,9 м. Более детальный ответ может быть получен после моделирования тепловых процессов. Для определения распределения мощности возле торцевой поверхности загрузки построен ряд диаграмм на линиях, параллельных торцу (рис. 4). В торцевой зоне загрузки распределение плотности тока и удельной мощности является более размытым, чем в регулярной зоне. Это в конечном счете приводит к уменьшению мощности. Диаграммы температуры вдоль осевой линии под витком и посередине между витками приведены на рис. 5. При моделировании принято время нагрева, равное 1 000 с. Затем моделируется процесс охлаждения в естественных условиях в течение 200 с. Благодаря высокой теплопроводности сплава перепад температуры по длине даже для длинной загрузки (1 м) составляет 70 в момент отключения нагрева. Выдержка в течение 50 с приводит к снижению перепада до 45, а через 100 секунд - до 30 . При уменьшении длины загрузки до 0,8 м перепад температуры снижается до 25 , что укладывается в заданные отклонения. Рис. 4. Диаграммы распределения удельной мощности на поверхности загрузки на линии, параллельной торцу: 1 - на торце; 2 - на расстоянии 10 мм от торца; 3 - на расстоянии 20 мм от торца; 4 - на расстоянии 100 мм от торца Результаты моделирования электромагнитных и тепловых процессов позволяют сделать несколько выводов, касающихся алгоритма поиска варианта, обеспечивающего максимум КПД и заданную точность нагрева. 1. Основная масса расчетов, связанная с поиском числа пазов, размеров паза и зубца, осуществляется на базе двумерных моделей. На этом этапе предварительно определяются КПД и интегральные характеристики индуктора. 2. С помощью трехмерных моделей рассчитываются краевые эффекты и определяется температурное распределение для различных длительностей интервалов управления. В случае необходимости производится корректировка размеров паза. Рис. 5. Диаграммы температур на поверхности по длине загрузки: 1 - под витком в момент отключения нагрева; 2 - между витками в момент отключения нагрева; 3 - под витком через 50 с после отключения нагрева; 4 - под витком через 100 с после отключения нагрева; 5 - под витком через 150 с после отключения нагрева; 6 - под витком через 200 с после отключения нагрева 3. На основании уточненных данных по размерам паза и уровню мощности производится расчет параметров индуктора и поиск оптимального режима нагрева. Представленный алгоритм помимо поставленных задач позволяет решить проблему сокращения затрат времени на поиск параметров индуктора.

About the authors

Alexander A Bazarov

Samara State Technical University

Email: aleksbazarov@yandex.ru
(Dr. Sci. (Techn.)), Professor. 244, Molodogvardeyskaya st., Samara, 443100, Russia

Alexander I Danilushkin

Samara State Technical University

Email: idanilushkin@yandex.ru
(Dr. Sci. (Techn.)), Professor. 244, Molodogvardeyskaya st., Samara, 443100, Russia

Danil N Pimenov

Samara State Technical University

Graduate Student. 244, Molodogvardeyskaya st., Samara, 443100, Russia

References

Supplementary files