CC BY
49
УДК 621.365.5 Е.С. Кинев
ОПТИМИЗАЦИЯ ЭНЕРГЕТИЧЕСКИХ ПОКАЗАТЕЛЕЙ ТРЕХФАЗНОГО ИНДУКЦИОННОГО НАГРЕВАТЕЛЯ
В статье рассмотрено решение задачи оптимизации индукционной системы промышленного трехфазного нагревателя, предназначенного для нагрева цилиндрической алюминиевой загрузки с позиций повышения электрического КПД, а также приведены результаты численного и экспериментального исследования индукционной установки, полученные с применением комбинированной математической модели.
Ключевые слова: индукционный нагреватель, математическая модель, мощность, динамика, ток.
Ye.S. Kinev OPTIMIZATION OF THE POWER INDICATORS OF A THREE-PHASE INDUCTION HEATER
The article considers the solution of the task of induction system optimisation for industrial three-phase heater which is designed to heat the cylindrical aluminium loading from the position of electric CUA increase, and it also gives the results of the numerical and experimental research of the induction plant, which were received with the use of the combined mathematical model.
Key words: induction heater, mathematical model, capacity, dynamics, current.
Нагрев металлов с высокой удельной электропроводностью выполняют при помощи индукционных нагревателей промышленной частоты с трехфазными индукторами. При расположении секций разных фаз индуктора встык вследствие эффекта переноса мощности между обмотками разных фаз нарушается симметрия токов трехфазной системы [1]. Уменьшению взаимной индуктивности секций и ослаблению переноса мощности способствует увеличение осевого зазора между обмотками индуктора, однако это приводит к увеличению провалов кривой удельной мощности в загрузке [2]. Несимметричную трехфазную индукционную систему в динамике нагрева следует рассматривать не только с учетом неравномерности температурного поля, но и с учетом электромагнитных процессов в силовой схеме. Для построения энергетически эффективной установки оказывается недостаточно решения задачи анализа физических полей в сложных индукционных системах. Поэтому на завершающем этапе выполняют комплекс исследований по многокритериальной оптимизации индукционной системы.
Комбинированная математическая модель, использованная в работе, построена на базе метода дискретизации свойств сред (МДСС), предложенного и подробно рассмотренного в [3]. Эскиз расчетной области трехфазного индукционного нагревателя, реализованной в математической модели, представлен на рис. 1.
Рис. 1. Эскиз расчетной области модели для анализа ЭМП в индукционной системе
В модели в двухмерной постановке учтено взаимосвязанное изменение электромагнитного (ЭМП) и теплового полей в электромагнитной системе индукционного нагревателя совместно с электромагнитными режимами в схеме электропитания с учетом температурной зависимости физических свойств загрузки, нелинейности магнитных свойств магнитопровода, неодинакового распределения потерь мощности по слоям обмоток индуктора, эффекта переноса мощности между обмотками разных фаз и отдельными секциями обмотки каждой фазы [4].
Комплексные векторы напряженности магнитного и электрического полей имеют следующие проекции:
н = ерНр+ёгНг.,
(1)
Эскиз расчетной области в подсистеме модели, предназначенной для анализа электромагнитного поля обмотки индуктора, представлен на рис. 2.
Рис. 2. Эскиз расчетной области в модели для анализа ЭМП в трехслойной обмотке
Касательная составляющая вектора напряженности магнитного поля на краях расчетной области у боковых стенок индуктирующего провода (і = 0,2 = I) равна нулю Нр = 0. Касательные составляющие вектора напряженности магнитного поля на краях расчетной области определяют по закону полного тока при заданном токе (/„):
р—<4-1!2
1 -к
(2)
Н2
<УС.
-к
'Ї
(3)
где Л/сл - число слоев обмотки; к - номер слоя; с(„ -диаметр к-го слоя; / - ширина витка с изоляцией.
После решения краевой задачи и расчета распределения векторов напряженности электрического и магнитного полей определяют резистивное и реактивное сопротивление индуктора п и Х1 через полную электромагнитную мощность (4) по теореме Умова-Пойтинга либо путем интегрирования джоулевых потерь и потокосцеплений (5):
м
*і=Х
к=1
К.І я1-./ ■2прА- - Нм,] ■ 2прмА.
■ = л + А’
(4)
N с..
Г{= I
к=1
КМ 2
И\Е,,; У І, АД2 2% Рі
і=1 ]= 1_____________
- w
N М | . |2
м I 1^7Г®йо^.,ЛрЛ--2лР,
_______:__________и , (5)
Т2 сл
к-1
Ї
*
*
2
И
2
И
И
где Ар, Az - шаг дискретизации расчетной области по осям р и z соответственно; wcn - число витков в слое; ю - циклическая частота.
Тепловой расчет в комбинированной модели реализован для двухмерного температурного поля в цилиндрической загрузке на каждом шагу итерационного алгоритма. Эскиз расчетной области для подсистемы анализа температурного поля в математической модели представлен на рис. З.
Рис. 3. Эскиз расчетной области для модели анализа теплового поля
На рис. 3 приняты следующие обозначения: 1 - немагнитная загрузка; 2 - футеровка; 3 - обмотка индуктора; 0 - плотности тепловых потоков на границах; X - коэффициент теплопроводности; с - теплоемкость; та -тепловые источники; Т - температура. Граничные условия заданы исходя из свободного теплообмена излучением и конвекцией с окружающей средой для торцевых частей загрузки. Потери с поверхности цилиндрической загрузки к водоохлаждаемой обмотке с температурой 20оС через футеровку определяются ее теплопроводностью.
Решение дифференциального уравнения Фурье выполнено численным методом конечных разностей (МКР) с применением неявного четырехточечного алгоритма конечно-разностной аппроксимации. В подсистеме анализа теплового режима рассчитывается температурное поле в загрузке при заданных на каждом временном интервале источниках тепла та и корректируются значения удельной электропроводности в каждой точке загрузки. Далее расчет продолжается в цикле, когда после определения интегральных параметров обмоток уточняется распределение токов и электромагнитного поля.
Схема силовых цепей трехфазной УИН с учетом несимметричных магнитных связей между фазами для этапа расчета электромагнитного режима силовой установки представлена на рис. 4. Для упрощения рисунка резистивные элементы батарей конденсаторов (БК) и секций индуктора для каждой фазы на схеме не показаны. Использование силовых вольтодобавочных трансформаторов ВДТ, автотрансформаторная схема включения индукторов И1, И2, Из и инверсное включение второй катушечной группы (рис. 4) наряду с другими особенностями схемных решений создают определенные трудности в построении единой математической модели, связывающей динамику электромагнитных и тепловых процессов в поле с электромагнитными процессами в реальной силовой схеме. Для их преодоления применена максимальная детализация элементов с построением схемных моделей всех звеньев силовой установки с применением расширенного элементного базиса теории цепей.
Рис. 4. Схема силовых цепей трехфазной индукционной установки
Взаимные (Мт) и собственные индуктивности (Ц и) различных частей индуктора определены в модуле анализа ЭМП после расчета мощностей (§ Р, О) путем численного интегрирования потока вектора Пойтинга (П) по площади каждой из секций (К11, К22, К33) в соответствии с рис. 1, а также в смежных областях (К12, К23).
$кк ~ |Пб/л - Ркк + jOkk ; Нкт = = /Окш ■
Ккк К,
(6)
кт
При заданных комплексах токов поочередно в каждой паре обмоток т и численно рассчитаны интегральные значения Ц 1т и Мт
(2кк^^^к' ^1)1 (~тт ^ ^ ^ т '
кт ~ ((2кт (2 к к О-тт ) ^ ® ^к ^т 1
(7)
(8)
где 1кк - действующее значение тока в соответствующей катушке.
В комбинированной математической модели для силовой установки применено описание схемных моделей с управляемыми источниками по компонентным уравнениям в базисе узловых напряжений и схемотехническое моделирование [2]. Это исключает необходимость решения топологических уравнений цепей, автоматизирует исследование электромагнитных режимов силовой схемы при большом разнообразии схемотехнических решений для силовых установок, а также позволяет учитывать оборудование в реальных режимах работы. Такой подход делает возможным переход непосредственно от распределенных параметров поля к сосредоточенным интегральным параметрам элементов электрической схемы.
Пример схемной модели трехфазного индуктора представлен на рис. 5.
Рис. 5. Схемная модель несимметричного трехфазного индуктора
Это решение для совокупности цепно-полевых электромагнитных и тепловых режимов сложной многофазной индукционной системы позволило оценить распределение мощности в загрузке с учетом многих факторов. Однако на этом этапе остался нерешенным ряд вопросов конструктивного исполнения собственно
индукционной системы и оценки ее влияния на эффективность преобразования электроэнергии в тепло и возможности повышения коэффициента полезного действия устройства в целом.
Для определения путей повышения энергетической эффективности был использован метод полного перебора по нескольким параметрам оптимизации. Перебор был осуществлен для числа слоев (No) от 2 до 12, высоты проводника (Ли), одинаковой для всех слоев от 0,002 до 0,01 м, частоты (f) питающего напряжения от 50 до 1000 Гц. В расчетах принято, что диаметр заготовки составляет d = 0,05 м, диаметр индуктора с/и = 0,09 м, длина индуктора / = 1 м. На рис. 6, а, б изображены полученные в расчете семейства линий уровней, включающих точки с максимальным значением КПД для промышленной частоты (г| = 0,453 при 12 слоях с высотой витков 3 мм и напряжении 700 В, то есть при использовании индукционного регулятора с удвоением выходного напряжения) и частоты выше промышленной (г| = 0,55 при 5-ти слоях с высотой витков 2 мм, частоте в диапазоне от 200 до 350 Гц и напряжении 600 В).
По результатам исследования сделан вывод, что при сохранении существующей схемы включения индуктора и использовании обмотки из плотно намотанных плоских проводников возможно повышение КПД с 0,441 до 0,55.
Другим способом увеличения энергетической эффективности является применение неравномерного сечения плоской шины по слоям многослойного индуктора. Однако существенное увеличение количества варьируемых параметров ограничивает возможность применения метода полного перебора и приводит к необходимости применения методов численной оптимизации. При этом одной из основных задач становится выбор метода оптимизации [5].
Рис. 6. Семейства линий уровней максимальных значений КПД
Анализ методов оптимизации применительно к УИН позволил выделить два наиболее эффективных алгоритма: комплексный метод Бокса и генетический алгоритм. Очевидно, что для определения оптимума гладких унимодальных функций, например, представленных на рис. 6, необходимо использовать метод из детерминистской группы, например, комплексный метод Бокса. В случае, если функция имеет несколько локальных экстремумов и один глобальный, то следует обратиться к методам стохастической группы, например, генетическому алгоритму.
Численные эксперименты показали, что если в состав варьируемых параметров не входит ширина витка, то зависимость КПД от параметров оптимизации гладкая и на этом этапе применен комплексный метод Бокса. В исследовании с учетом ширины витка использован генетический алгоритм.
Важной задачей является выбор критерия оптимальности. Эффективность процесса нагрева алюминиевой заготовки в индукционном нагревателе можно косвенно оценить по коэффициенту полезного действия установки. Задача поиска оптимальной обмотки была сведена к нахождению минимума критерия оптимальности, поэтому в соответствии с выражением (9) минимизировался максимум относительного расстояния значения КПД от заданной цели. Однако высокое значение КПД - это не единственное техническое требование, которое предъявляют к разрабатываемой индукционной системе. Например, зачастую необходимо выполнить условие согласования параметров нагрузки и с параметрами источника питания по напряжению или по техническим требованиям, следует ограничить плотность тока в обмотке. Ограничение допустимой
плотности тока в индуктирующих проводниках обмотки нагревателя на уровне 8 А/мм2 позволяет обеспечить эффективный съем тепла с медной обмотки при интенсивном обдуве с помощью вентилятора специального назначения. Это, в свою очередь, предоставляет возможность отказаться от водяного охлаждения и, как следствие, существенно повысить эксплуатационную надежность индукционного нагревателя.
Для оптимизации индукционной системы индукционного нагревателя по техническим условиям, рассмотренным выше, использована многокомпонентная функция цели:
^ = ^(1 - л)/1 + к2\(ип -ии)/и„\ + кз (5Д -8)/8,
^■тт,
(9)
где к], кг, кз - весовые коэффициенты; г|и - альтернативные значения КПД при изменении оптимизируемых параметров; ин - напряжение на индукторе; Ц, - расчетное напряжение; 5 - расчетная плотность тока в индуктирующем проводе; 5Д - допустимая плотность тока по условию охлаждения обмотки.
Выражение (9) можно использовать для максимизации только одного фактора - КПД, приняв равными нулю весовые коэффициенты к2 и кз, или осуществлять поиск оптимума по двум слагаемым выражения (9), приравняв нулю соответствующий весовой коэффициент.
Разработанная методика параметрической оптимизации опробована на задаче оптимизации, решение которой получено выше при помощи метода полного перебора. Оптимизационный поиск для индукционной системы нагревателя заготовки алюминиевой катанки проведен по трем параметрам (частота, высота индуктирующего проводника, равная для всех слоев, количество слоев обмотки) с помощью комплексного метода Бокса и генетического алгоритма с использованием целевой функции (9) при кз = 0. Полученные данные подтвердили результаты в методе полного перебора.
История трехмерного оптимизационного поиска для обоих методов представлена на рис. 7, а, б.
Рис. 7. История оптимизационного поиска с помощью генетического алгоритма (а) и метода Бокса (б)
В рассмотренном случае поверхность отклика целевой функции гладкая унимодальная, что обеспечивает преимущество комплексного метода в скорости при сохранении надежности нахождения глобального экстремума.
Как показал численный эксперимент, достоинством первого метода является высокая скорость нахождения оптимального решения, второго - надежность нахождения глобального оптимума при наличии локальных минимумов. Тем не менее комплексный метод имеет некоторые существенные ограничения, свойственные для всех детерминистских методов.
В случае мультимодальной целевой функции независимо от количества оптимизируемых параметров, как показало исследование, он часто не в состоянии найти глобальный оптимум. Поэтому, чтобы увеличить вероятность нахождения искомого решения - глобального оптимума - процедуру оптимизации необходимо повторять несколько раз из различных стартовых точек.
Для оптимизации конструкции обмоток УИН критичным параметром является ширина индуктирующего провода. Так, например, для двухслойного индукционного нагревателя была получена поверхность функции КПД в зависимости от длины проводника в каждом слое (рис. 8).
Рис. 8. Поверхность функции КПД
Проекция поверхности искомой функции КПД приведена на рис. 9, а. Как видно из рис. 9, а, оптимизационная функция имеет несколько локальных экстремумов и один глобальный экстремум, укрупнено представленный на рис. 9, б.
После подтверждения работоспособности предложенной методики оптимизации конструкции УИН выполнены исследования путей повышения энергетической эффективности установки индукционного нагрева заготовки алюминиевой катанки перед волочением.
а б
Рис. 9. Проекция поверхности функции КПД (а) и глобальный максимум (б)
С помощью разработанного программного обеспечения было исследовано влияние высоты и ширины прямоугольного индуктирующего провода на потери в обмотке, как по отдельности, так и совместно. Исследованы диапазоны изменения: длины витка - от 10 до 70 мм, высоты витка - от 2 до 10 мм, частоты напряжения источника - от 50 до 1000 Гц.
Исследование функции КПД комплексным методом Бокса дало положительный результат только при проведении процедуры оптимизации несколько раз из различных стартовых точек. В свою очередь, исследование той же функции с помощью генетического алгоритма гарантирует высокую точность нахождения глобального экстремума, хотя и требует значительно больше итераций. Независимо от начальных условий, оптимизационная процедура, реализованная с помощью генетического алгоритма, во всех исследованных случаях находила глобальный экстремум.
По результатам оптимизационного поиска, представленным в таблице, сделано заключение, что использование существующей схемы электропитания установки (вариант № 1) ограничивает увеличение КПД на уровне п = 0,503 (вариант № 2).
Сравнение результатов параметрической оптимизации ИН
Наименование Вариант исполнения
№ 1 № 2 № 3 № 4
Весовой коэффициент - 0,4 0,3 0,3
0,6 0,7 0,4
0,0 0,0 0,3
Число фаз 1 1 3 3
Частота, Г ц 960 941 149 152
Число слоев обмотки 1 3 8 7
Мощность ИН, кВт 188 124 101 112
Электрический КПД 0,44 0,50 0,541 0,556
Напряжение ИН,В 760 750 380 380
Плотность тока, А/мм2 17,3 15,7 12,9 7,8
Тиристорный преобразователь частоты позволяет регулировать частоту питающего напряжения в пределах от 800 до 1000 Гц. Поэтому предусмотрен альтернативный вариант модернизации с применением транзисторного преобразователя частоты на трехфазную нагрузку, что позволило получить значение КПД П = 0,541 при водяном охлаждении (вариант № 3). Вариант исполнения № 4 предполагает использование принудительного воздушного охлаждения обмотки индуктора. Он имеет сравнимый КПД п = 0,556 при меньших текущих расходах вследствие отсутствия затрат на оборотную воду. Кроме того, вариант № 4 имеет значительно более высокую эксплуатационную надежность, чем вариант с водяным охлаждением обмоток.
По результатам оптимизации для трех лучших вариантов выполнены технические проекты модернизации установки индукционного нагрева заготовки алюминиевой катанки. Результаты работы переданы на Иркутский алюминиевый завод для промышленного внедрения.
Литература
1. Кинев, Е.С. Математическое моделирование физических процессов при индукционном нагреве / Е.С. Кинев, Е.А. Головенко // Оптимизация режимов работы электротехнических систем: межвуз. сб. науч. тр. - Красноярск: ИПЦ КГТУ, 2006. - С. 19-30.
2. Analysis of power redistribution between windings of multi-phase induction heaters / Е.А. Golovenko, E.S. Kinev, E.V. Kuznetsov [at el.] // International Symposium on Heating by Electromagnetic Sources HES - 07 (Section P ). - Italy, Padua, 2007. - P. 328-336.
3. Тимофеев, В.Н. Анализ электромагнитного поля методом дискретизации свойств сред: Математическое моделирование физических полей в алюминиевых электролизерах: моногр. / В.Н. Тимофеев. -Красноярск: ИПЦ КГТУ, 2002. - 264 с.
4. Исследование распределения мощности и температуры в трехфазных индукционных нагревателях / Е.А. Головенко, Е.В. Кузнецов, Е.С. Кинев [идр.] // Вестн. ВГТУ. - 2007 - Т. 3. - № 6. - С. 170-177.
5. Кинев, Е.С. К вопросу оптимизации несимметричных режимов индукционных установок в трехфазной сети / Е.С. Кинев, Е.А. Головенко, Е.В. Кузнецов // Оптимизация режимов работы электротехнических систем: межвуз. сб. науч. тр. - Красноярск: ИПЦ КГТУ, 2006. - С. 4-14.
'--------♦------------
