Магнитогидродинамические перемешиватели жидких металлов с несинусоидальными токами Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

Journal of Siberian Federal University. Engineering & Technologies 2 (2013 6) 166-177

УДК 621.643; 537.8; 532.5

Магнитогидродинамические перемешиватели жидких металлов с несинусоидальными токами

В.Н. Тимофеев*, Г.ф. Лыбзиков, М.Ю. Хацаюк, М.А. Еремин, С.П. Тимофеев

Сибирский федеральный университет, Россия 660041, Красноярск, пр. Свободный, 79

Received 15.03.2013, received in revised form 22.03.2013, accepted 31.03.2013

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

Ключевые слова: магнитогидродинамический перемешиватель, жидкий металл, электромагнитные силы, индуктор, периодическое напряжение, периодический ток.

Введение

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

В настоящее время широкое распространение получили индукционные магнитогидродинамические (МГД) перемешиватели, позволяющие бесконтактно воздействовать на расплав и полностью автоматизировать этот процесс.

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

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

© Siberian Federal University. All rights reserved

* Corresponding author E-mail address: [email protected]

работы источника питания МГД-перемешиватель, питаемый несинусоидальными периодическими токами, имеет явное преимущество перед МГД-перемешивателями с синусоидальным питанием, а наличие пульсирующих электромагнитных сил способствует устранению микро-неоднородностей в многокомпонентных расплавах.

По принципу действия индукционные МГД-перемешиватели аналогичны линейным индукционным машинам [2]. На рисунке 1 представлен эскиз печи-миксера с МГД-перемешивателем, установленным под подиной ванны [3]: 1 - ванна печи-миксера; 2 - жидкий металл (расплав); 3 - многофазная обмотка индуктора; 4 - сердечник индуктора; 5 - источник питания.

Индуктор МГД-перемешивателя состоит из многофазной обмотки 3 и сердечника 4. Длина сердечника Ь = 2рт, где р - число пар полюсов; т - полюсное деление.

При подключении многофазной (т = 2, 3 - число фаз) обмотки к источнику периодического синусоидального или несинусоидального напряжения в зазоре между индуктором и расплавом образуется бегущее магнитное поле. Скор) ость движения амплитуды магнитной индукции в направлении осиХ равна и1 = 2и£ где f— частота тока. Пре взаимодействии бе гущего магнитного поля с наведенными в расплаве электрическими токами расплав начинает двигаться со скоростью и также в направлении оси х, осуществляя перемешивание расплава в ванне.

Эффективность перемешивания зависит от глубины проникновения электромагнитного пооя в расплав. Синусоидальное электромагнитное поле по мере распространения в расплав (ось X) затухает по экспо неециальному зокону [2]. Так, зависимость комплексной величины магнитной индукции от координаты (X) описывается выражением

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

(1)

где Б0 - величин;! магнитной иидукции на поверхности расплава (Т). Комплексное значение ^ равно

(2)

у,

ИП

где т - полюсное деление индуктора (m); со = 2nf— угловая частота (rad/s); f- частота тока (1/s); fio = 4п10'7 - магнитная постоянная (H/m); у - удельная электропроводность расплава (1/Q m);

Расстояние, на котором электромагнитное поле уменьшается в е~2,72 раза (е - основание натурального логарифма), называется глубиной проникновения электромагнитного поля. В соответствии с (1-2) глубина проникновения электромагнитного поля Д обратно пропорциональна ф и при не изменнык свойетвах расплава оависит от частотыы токаf и полюсного деления индуктора а. При т-—о глубина празикновения эзектромагнитного поля равна глубине про-никновеная плоской элеатромакнитной 1золны>1 и при неизменныых свойствзх расплава полностью зевааит от частоты! тока /[4]:

В реальнызх индукторая величина полюсного деления т имеет ограниченное значение и обыачно выбирается из услолиа т3=Л(33—где <5 -а веоичина расстояния ет магнитопровода ин-дукторв до расплааленноео металлн "Чем меньше т, тема быстрее зату хает электромагнитное по ле в расславе. При фиксированном о глубино проникновения электромагнитного по ля также зависит от частот ы тока / Еслиf=0 (постоянный ток в обмотке индуктора), то в соответствии с (2) глубину проникнооения магазтаого поля (также постояннтго во времени) будет максимальной и рсвной

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

Известно, что источниками электромагнитного поля индукаора являются электрические токи в облотке индукалрл. В связи с этим целесообразно сравнить электромагнитные характеристика двух МГЛ-перемеши вателей . Первый питается синусоидальными периодическими токами, второй - несинусоидальными, которые в течение полупериода Т/2 имеют один интервал времени, когда д1дф0, а в течение другого интервала времени д1д = 0.

При рассмотрении второго варианта предположим, что на каждую фазу обмотки индуктора илдается периодическои напряжлние прямоугольной формы1 (рис. 2). Поскольку каждая фаза имеет активное сопротивления Iе и индуктивность Ь, т ней возникает переходный процесс при каждом изменении полярности напряжения. Длительность переходных процессов tп зависит от параметров индуктора II и Ь, а именно определяется постоянной времени тв и приблизительно рввна

(3)

д =

т

(4)

tп * 3т,

(5)

где [5].

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

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

и

I =+_»_

у Я

(6)

Электромагнитное поле будет повторять форму токов. Условно процесс взаимодействия такой формы поля с жидким металлом ванны можно раздевить на два. Первый - когда поле не меняется (вр емя tй и проникает в толщу металлй максимально глубоко, и второй - когда поле меняется во времени с большой! скоростью (время () и в толгме металла 5 наводятся кольцевые (вихревыб) токи, которые, взаимодействуя с полем, создают силу, приводящую металл 5 в движение. Для эффективного перемешивания жидкого металла в ванне целесообразно на расплав воздействовать бегущим электромагнитным полем. Для этого в индукторе устанавливают две и более обмотки, которые запитываются от источника одинаковым по форме и частоте напряжением, но при этом напряжения сдвинуты по (разе . Например, для двухфазного индуктора (т = 2) фазивый сдвиг между напряжениями может быть Т/С (рис. 2).

В получе нных выражениях для ущкощения взаимная инддктивность между обмотками не учитывалась.

В трехфазном индукторе (т = 3) фазовый сдвиг между питающими напряжениями может составлять Т/3.

Так как длительность переходного процесса tп~3L/R, а длительность установившегося состояния

/ = Т/2 - 3 — у Я

(7)

можно получить различные значения ty путем изменения параметров индуктора L и R при постоянном значении периода Т.

Далее определим электромагнитные характеристика МГД-перемешивателей с синусоидальным и несинусоидальным периодическим питанием при различных значениях относительного параметра ty* = ty/(0,5T).

Математическое и экспериментальное моделирование

На рис. 3 представлена расчетная 2D-модель системы «индуктор - расплав», построенная в расчетном комплексе ANSYS. Для построения модели были выделены следующие подобласти: 1, 2 - индуктор; 3 - расплав; 4 - окружающее пространство. Подобластям были присвоены следующие материалы с соответствующими электрофизическими свойствами: 1 - медь (р = 1,7e-8 Qm, /r = 1); 2 - шихтованная трансформаторная сталь (рх = 14e-8 Qm, ру = 14e-8 Qm, pz = 1e32 Qm, /r = 300); 3 - алюминий при t = 730 °C (р = 24,76e-8 Qm, /r = 1).

При постановке задачи приняты следующие допущения:

1. Магнитная проницаемость магнитопровода имеет постоянное значение.

2. Расчет электромагнитного поля проводится без учета поля скоростей.

3. В качестве источника задается источник тока, имеющий линейные переходы.

Система уравнений, описывающих электромагнитное пол«;, состоит из уравнения на векторный магнитный потенциал и уравнения непрерывности. С учетом принятых допущений эта система имеет вид:

V2 A = -Ма Ô, (8)

V Ô = 0 , (9)

где A - магнитный потенциал, (Fs/m); /а - абсолютная магнитная проницаемость; ô - вектор плотности электрического тока, (Л/т2) :

Г . .S ■ 4 J

Г Y'

' ' В:-'

ИДЯ

б = -у

дА

д(

(10)

Чтобы электромагнитное по ле в расчетной области определялось однозначно, дополним систему уравнений (8 - 9) граничными условиями. Считая, что расчетную область окружает среда с идеальными магнитными свойствами (м = да), на границах расчетной области для касательной и нормальной составляющих векторного магнитного потенциала будут выполняться условия:

дА,

— = 0.

(11)

дп

А=0, (12)

где Ал и Аи - тангенсальная и нормаленая составляющие ве кторного мегнитного потенциала; п -а нормаль к поверхности.

В качестве источникового члена задавался источник тока. функция двух фаз синусоидального периодическооо тока имеет вид:

п (') = 1т 8ш(2П?), П (') = + п / 2),

(13)

несиносоидаюьнаго периодического тока:

—) =

п 1

т' -

у

- п

- 1т

, 0< 1<Гп , + << < '„+'+

" У,, < 2 < + 2,,

' т(=) =

(

2

V-

+ - 1т

2'п + 'у <'< 2'п + 2'у

Т/4 < 2 <774 + 'п Т/4 + 'п<7 <Т 14 +'п+'у 77/4 + Пп + '< <1тТП4 + 2+ + 'у ,Т/4 + 2'п +'у < 7 <Т/4 + 2/п + 2'у.

(14)

Верификации результатов численного моделирования

Доя проверки построенной мятематической модуло была разработана и с обранк экспериментальная установка (рис. 44) 2-фазного индуктора 1 с импульсным источником питания 2. Частота питания индуктора 3±0,3 И током 20 А.

Измерение электромагнитного поля проводили с помощью двух датчиков Холла, расположенных ортогонально друг к другу. Такие датчики дают две пространственные компоненты вектора напряженности магнитного поля. Результаты измерений выводили на дисплей осциллографа 3 и в текстовый файл. Воздействие периодическим несинусоидальным электромагнитным полем осуществлялось на алюминиевую заготовку 4.

Рис. 4

Датчик для выполнения измерений помещали между центральным зубцом магнитопрово-да (рис. 4 и 5, точка А) и алюминиевой заготовкой и в той же плоскости в толщине заготовки (рис. 4и 5, точка Б).

С аналогичными параметрами было выполнено численное моделирование с использованием ргзработанной математической модели [6], в результате которого была получена динамика электромагнитного поля. Так же из точки в расчетной области между центральным зубцом магнитопровода и заготовкой и из точки в толщине заготовки были выведены в текстовый файл две пространственные компоненты вектора нгпряженности магнитного поля.

График-сравнение по лученных данных представлен на рис. 5. Некоторые расхождения результатов математического моделирования и физического эксперимента объясняются тем, что, во-первых, при расчете использовалась 2Б-модель, во-вторых, при проведении эксперимента частотг прямоуаольно го периодического напряжения отличалась от расчетной из-за специфики источника питания. Вместе с тем из приведенных графиков видно, что математическая модель правильно воспроизводит характер формы магнитного поля и может быть использована для анализа электромагнитных процессов в МГД-перемешивателях с несинусоидальными периодическими токами.

Анализ результатов численного моделирования

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

На рис. 6 представлены картины распределения мгновенных значений напряженности магнитного поля Н, магнитной индукции В и мгновенных значений объемных электромагнитных сил f в расплаве.

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

- 172 -

0.7 точка А

0,6 0.7

точка В

Рис. 5

картин распределения Ни В позволяет правильно проектировать магнитопровод и выбрать линейную токовую нагрузку.

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

На рис. 7 изображены зависимости относительных составляющих электромагнитной силы F* = ехЕ*х+е1Е*2 от относительного параметра 2t/T для двухфазной и трехфазной обмоток индуктора соответственно. Здесь Е*х и Е*1 - тангенциальная (ось х) и нормальная (ось ¿) составляющие электромагнитной силы соответственно; е и е2 - единичный вектор, направленный по оси х и г соответственно.

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

А

Рис. 6

| ? [

А

НТШШТааЭ нттмтнк »^с.^ТлТ.вм от*т.1н> ■ ■

В

Г\ «IV. 3

П

\

ч

* л/ — 11

г. 1 - ™

N

1,1 о.* ф.с я, т

от о.4

ол л,

Рис. 7

р* =.

р * =.

(15)

Пунктирные линии на ¡рис;. 7 соответствуют тангенциальным, а сплошные - нормальным составляющим электромагнитных сил.

Как следует из приведкнных графиков, при

Т„ /0 ДТ >0,4

(16)

Рис. 8

значения как тангенциальной, так и нормальной составляющих электромагнитной силы при несинусоидальном питании обмоток индуктора превосходят аналогичные значения при питании обмоток индуктора синусоидальными токами. Следовательно, если форма питающих обмотку индуктора несинусоидальных токов имеет такое значение длительности установившегося состояния 0, при заданном периоди Ту при котором выполняется неравенство (16), перемешивание будет более эффективным по сровнению с вариантам, I! котором обмотка индуктора будет питаться синусоидальными токами.

Величину периода Т питающего напряжения можно выбрать из следующих соображений. Известно, что зависимость относительной электромагнитной0 силы от безразмерной уеличины е - коэффициента добротномти пру различныт значениях дКт имеет вид, представленный на рис. 8. Здесь д - расстояние от индуктора до расплава металла. Коэффициент добротности определяется как [2]

(17)

При известном значении S/т можно определить значение коэффициента добротности е = ет, при котором относитеиьная электромагнитная сила будет иметь максимальное значение F*3M = F*MM max. Здесь F*3M = FJF0, где за базисную силу принята Fo = p<>2pvaA, (Ж), где p - число пар) полюсов индуктора ; <а - ширина индуктора, m; A - линейная тоиовая нагрузка, (A/тг).

Зная ет, в соответствии с (17) можно определить период (T = 1/f), а именно

T =

2ц0ут2

(18)

Следовательно, индукционная! установка, обмнлкл индуктора которой выполнена таким образом, чтобы соблюдалось условие ( 16), паоволяет подо браль чаятоту многофазного (т = 2, 3) источника несинусоидаиьного периодического напряжения (выражение 18), при которой будет возможно максимально эффективное перемешивание металла.

Следует отметить, что при наличии в электрической це пи с последовательным соединением R и L элементов периодического тока i напряжение на ее зажимах определяется выражением

di

u = Ri + L —. (19)

dt

Во время установившегосн состояния ty изменение электричелкого тока отсутствует (di/dt = 0)и тогда в соответствии с (19) u = Ri.

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

Заключение

Таким образом, в результате сравнения двух вариантов МГД-перемешивателей вариант МГД-перемешивателя с несинусоидальными периодическими токами при определенных параметрах индуктора и режимах работы осуществляет более эффективное перемешивание жидких металлов, при этом отмечено снижение потребления реактивной мощности (повышение коэффициента мощности) и сокращение времени перемешивания.

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

[1] Тимофеев В.Н., Корчагин А.И., Павлов Е.А., Тимофеев Н.В. // Journal of Siberian Federal University 1 (2011 5) 28-38.

[2] Вольдек А.И. Индукционные магнитогидродинамические машины с жидкометалличе-ским рабочим телом. Л.: Энергия, 1970. С. 272.

[3] Тимофеев В.Н., Христинич Р.М., Бояков С.А., Темеров А.А. // Международный журнал для алюминиевой индустрии. 2009. № 1. С. 16-20.

[4] Демирчян К.С., Нейман Л.Р., Коровкин Н.В., Чечурин В.Л. Теоретические основы электротехники: учебник для вузов. Т. 2. 4-е изд. СПб.: Питер, 2003. 576 с.

[5] Бессонов Л.А. Теоретические основы электротехники: учебник для студентов энергетических и электротехнических вузов. М.: Высшая школа, 1973. С. 752.

[6] Pervukhin M.V., Minakov A.V., Sergeev N.V., Khatsayuk M.U. // Mathematic simulation of electromagnetic and thermal-hydrodynamic processes of the system "inductor-ingot" of an electromagnetic mould. 2011. No. 1. P. 79-88.

Magnetohydrodynamic Stirrer Liquid Metal with a Non-Sinusoidal Currents

Viktor N. Timofeev, Gennady F. Lybzikov, Maxim Yu. Khatsayuk, Mikhail A. Eremin and Sergey P. Timofeev

Siberian Federal University 79 Svobodny, Krasnoyarsk, 660041 Russia

The article compared the induction MHD stirrer liquid metals with sinusoidal and non-sinusoidal periodic currents. On the basis of numerical simulation software product ANSYS, define differential and integral electromagnetic characteristics magnetohydrodynamic stirrer with two-phase and three-phase windings. It is shown that under certain conditions the inductor and the frequency of the supply current, magnetohydrodynamic stirrers with non-sinusoidal currents are more efficient and less energy-intensive than the magnetohydrodynamic stirrers with sinusoidal currents.

Keywords: magnetohydrodynamic stirrer, liquid metal, electromagnetic force, inductor, periodic voltage, periodic current.