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

Journal of Siberian Federal University. Engineering & Technologies, 2017, 10(1), 74-84

УДК 517.9; 537.8; 004.94

Analytical Calculation of Differential and Integral Electromagnetic Characteristics of MHD-Stirrer

Viktor N. Timofeev and Maksim Yu. Khatsyuk*

Siberian Federal University 79 Svobodny, Krasnoyarsk, 660041, Russia

Received 14.11.2016, received in revised form 26.12.2016, accepted 19.01.2017

The analysis of electromagnetic field in the system «inductor - bath of molten (liquid metal)" by different winding characteristics of the inductor due to analytical method is shown in this paper. Integral and differential characteristics of the system are defined and compared for the most common coil inductors MHD-stirrers. It provides a framework for the analysis offield harmonics on the surface of melt for various types of inductors. Also made the comparative assessment of influence of longitudinal edge effect on system performance.

Keywords: MHD-stirrers, electromagnetic characteristics, analytical calculation.

Citation: Timofeev V.N., Khatsyuk M.Yu. Analytical calculation of differential and integral electromagnetic characteristics of MHD-stirrer, J. Sib. Fed. Univ. Eng. technol., 2017, 10(1), 74-84. DOI: 10.17516/1999-494X-2017-10-1-74-84.

Аналитический расчет дифференциальных и интегральных электромагнитных характеристик МГД-перемешивателя

В.Н. Тимофеев, М.Ю. Хацаюк

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

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

© Siberian Federal University. All rights reserved

Corresponding author E-mail address: maxhac@ya.ru, viktortim0807@mail.ru

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

Введение

В процессе приготовления многокомпонентных сплавов важной технологической операцией является перемешивание расплава с целью выравнивания химического состава и температуры во всем объеме ванны миксера. Перемешивание расплава в плавильных печах позволяет повысить интенсивность плавления и сократить на эту операцию потребление энергии. Основным достоинством электромагнитного перемешивания признано отсутствие контакта с высокотемпературным и агрессивным расплавом и возможность автоматизировать операцию перемешивания. В металлургической практике наибольшее распространение получили индукционные магнитогидродинамические (МГД) перемешиватели, воздействующие на расплав бегущим магнитным полем. Источником бегущего магнитного поля служит индуктор, представляющий собой катушки, расположенные на магнитопроводе, сдвинутые в пространстве и питаемые да-фазной системой переменных токов. Как правило, электромагнитные, тепловые и гидродинамические процессы в миксерах с МГД-перемешивателями моделируются численными методами с использованием таких программных пакетов, как А№У8 и т.д. [1] На ранней стадии проектирования МГД-перемешивателей требуется выбрать тип обмотки индуктора, количество фаз, частоту питающего напряжения, токовую нагрузку и т.д. Для этого требуется простая методика проектирования этих устройств. В [2] решена одномерная задача анализа электромагнитных процессов в МГД-перемешивателе. Интегральные и дифференциальные характеристики получены при наличии только первой гармоники электромагнитного поля в зазоре. В данной статье аналитически решена двухмерная задача по анализу электромагнитного поля в МГД-перемешивателе с применением рядов Фурье в комплексной форме. Полученные выражения позволяют учитывать тип обмотки, количество и размеры пазов, распределение токовой нагрузки и влияние продольного краевого эффекта.

Расчетная модель системы «МГД-перемешиватель - ванна с расплавом»

Эскиз системы, состоящей из ванны с жидким металлом (расплавом) и индуктора бегущего магнитного поля, представлен на рис. 1.

Индуктор расположен на расстоянии 3 от расплава и имеет длину Ь=2рт, где р - число пар полюсов; т - полюсное деление.

Для анализа электромагнитного поля примем следующие допущения:

1. Магнитопровод индуктора имеет магнитную проницаемость ¡л=ж и электропроводность у=0. В М пазов магнитопровода уложены проводники с током, каждый да-й паз заменен поверхностным током с комплексной линейной плотностью

w I

/ , (1)

т 2Д К>

т

где Vда - число проводников; 1т - комплексный ток; 2Дда - ширина да-го паза.

2. Воздушный зазор равномерный и равен

? = (2) где кд и кт - коэффициенты, учитывающие влияния пазов и магнитного сопротивления ярма и зубцов соответственно.

3. Система имеет бесконечную длину в направлении оси у, т.е. электромагнитное поле является плоскопараллельным. Вектор комплексной напряженности электрического поля Е=ёу-Еу, вектор комплексной напряженности магнитного поля Н=ёхНх+ё2Н2, где ёх, ёу, е2 -единичные векторы прямоугольной системы координат.

4. Систему координат считаем жестко связанной с индуктором, а все элементы металла в ванне двигаются относительно индуктора с одинаковой скоростью и.

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

Дифференциальные уравнения и граничные условия

Решение будем искать для комплексной напряженности электрического поля в области

1 (0 < г < 5') и области 2 (5' < г < да) Е1=ёуЕ1у; Е2=ёуЕ2у..

В расчетных областях справедливы дифференциальные уравнения и граничные условия [3]:

д2 Е. д2 Ёу д2 Е 2 у д2 Е 2 у . дЕ 2 у -Г- +-г- = 0; -2^ +-2^-' -Мо-7-ю-Е2у-н>-Г-и-^ = 0; (3, 4)

&2 5х2 522 5х2 0 2у 0 5х

Рис. 1. Индуктор бегущего магнитного поля вблизи ванны с расплавом: 1 - прямоугольная ванна; 2 - расплав; 3 - многофазная обмотка; 4 - магнитопровод индуктора

Рис. 2. Расчетная модель МГД-перемешивателя расплава в ванне, xm - расстояние до центра m-го паза от начала системы координат

dK

dz

■Mo

E1y E 2 y

s \z=5 J \z=5

Jm ; Xm ~A m * X * Xm +A m

0

dE.

dE, y

dz

2 y

dz

IE,

2 У

= 0.

(5)

(6, 7)

Здесь 1=^-1; ц0=4 ж-Ш Гн/м; т=2 ж/с"1; /- частота тока, Гц; у - удельная электропроводность расплава, 1/(Ом-м).

Комплексная линейная плотность тока определяется выражением (2), где

J = J в""".

m m

(8)

Величина фазы да-го тока зависит от количества фаз, числа пар полюсов р и количества пазов М.

Решение дифференциальных уравнений

Полагая, что электромагнитное поле периодично по координате х с периодом Т, будем искать решения в виде рядов Фурье в комплексной форме [4]:

El,2 (x, z) E№ (z)eiknX,

1 T

Eia„ (z) = - J El,2 (x, z)e-ik'xdx;

kn =

2 пж

(9) (10) (11)

Умножив уравнения (3) и (4) на функцию е п и проинтегрировав получившееся выражение по х от 0 до Т, получим:

d 2 E1n к 2E

dz

■2- - К2En = 0;

d2 E.

dz

2r -Vn2 E 2n = 0;

(12)

(13)

Ф„ = К +1 -Мо +1 ■ к„ -Мо 'Г'и- (14)

Общие решения обыкновенных дифференциальных уравнений (12) и (13) определяются выражениями

EEln = Aeknz + Л2в-кп*;

Егп = Л3в^"2 + Л4в-

(15, 16)

Постоянные интегрирования Л] (]=1..4) определяются из граничных условий (3-5), которые также необходимо преобразовать в соответствии с (10). При этом условие (3) заменяется выражением

dE,,

dz

= Wn

(17)

z=0

z=S

z=S

z

z=0

где

¥n =

imu M Xm »

£ Jm j e-"'xdx = i ^ £ J sm(k„ Am )e

T m=l

m

m=1 x...-Д.

KT m=1

(18)

а условия (4) и (5) сохраняются для Е1п и Е2п.

Определив постоянные интегрирования из граничных условий и подставив их в (15) и (16), получим:

E _ Vn Gn(z).

1Я~ К Qn '

n(S-z)

E2n = ~Wn -

Qn

где

Gn(z) = kn cosh^n(z - 8)) - Vn sinh(K(z - £)).

Qn = kn sinh(kn8) + Vn coSh(K8)-

(19, 20)

(21) (22)

В соответствии с (6) запишем искомые решения:

¿1, (x, z) = -i ^ £ Л ^ C,_eil!"x'

T Qn

где

ee2..(x, z) = -, -L^p Cmek

1 n=-a kn Qn

M

Cin =E Jm sin(k„ Am )e-iKxm.

m=1

Из уравнения Максвелла rotE=-ifflB, имеем: txx -

ia dz

Подставив (23) и (24) в (26), получим:

t =- 2BL у _1G'n(z) C e'k,x.

x1 rp ¿„I 7 1n ■

1 n=—& n Sin

B z = — ;

irn dx

^ - (^-z)

B 2 =zl-l у 2

9n e

n=-OT kn Qn

C e^'

iknx.

B 1 = i IbL V1 GnSzlc e'k"x;

z1 t h kn Qn ln

XX ='у1 e_C eik«x.

nz2 1 T k Q ''

1 n=-W kn Qn

где

(23)

(24)

(25)

(26)

(27, 28) (29, 30)

(31)

(2) = кп втЬ^ (2 - 8)) - фп С0&(кп (2 -

Плотность тока в области 2 определяется законом Ома с учетом движения расплава со скоростью и относительно системы координат:

. . . 2п -A ev"(s~z)

¿2, = Г(Е 2 + uB2 z ) = i— X п C 1nP(S)e

Тт nt^ nQn

iknx.

P(s) = e| n(1 - s)-

T

(32)

(33)

-ikx

e

где s = оау1лйтг /л2; T* = T/r.

Интегральные и дифференциальные характеристики МГД-перемешивателя

Комплексная электромагнитная мощность, отдаваемая обмоткой индуктора в зазор, равна

4aMols ^ 1 Оп (0)

^ =-l,Z J Ei(x,0) • J'mdx = i£

C C

m=1 x —Л

Г/ ^ , 3 w1n 2n'

n=-«kn Qn

(34)

где С2п = ^ 3 т 8т(кпДт )е 'к"х-; 15 - ширина индуктора по оси у.

т=1

Электромагнитная сила, действующая на обмотку индуктора в направлении движения, равна силе, действующей на расплав, и определяется выражением

Fem = Rejl.S JДг(x,0)J>URe-ji^ £

4Mols ^ 1 Gn (0)

2 ^ w1n 2n

m=1 xm —Д m | [ J n=-« kn

CinC2n k (35)

Т п=-Лп Qn

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

f = Re{¡>2 х В2 )= ef + ef,

где

fx = ReR В J

fz =- Re{4 X)

(36)

(37, 38)

2 ^ 2 ^ Jz (.2^

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

р =ojj4s2_ j2.

л

Pb 2

Fb = Mols2r- J2,

2т-f

(39, 40)

M W I где J = Y^Lm.

m=1 T

В соответствии с (34), (36), (39) и (40) будем иметь:

о = Sm = ^ у 1 (Gn(0)

em Ps 4л hn3 Qn

C1nC2n; Fm = FFm = Refc £4-^cc К (4l, 42)

n Qn

где

Gn(0) = kncosh (K (z - - Ф„smh (K (z - Qn = knsmh(hs) + Фп cosh(hs)

2n

<Pn I + ie

(1 - Tо -.));

(43)

x„ +Д

x... +A

M

Cn = I - sink A m )e-kx-; С2и = £ К sin(k„ Am )e"

m Упт^ > 2n / > m V n m г

m=1 m=1

k =Lk = 2П; T = -; J = J-mm-; - = -

n * n - t m - m - "

Относительные активная и реактивная мощности соответственно таковы:

Pmm = Re{^em J Qem = (44)

Коэффициент мощности МГД-перемешивателя и электромагнитный коэффициент полезного действия устройства в двигательном режиме:

cosp = n=Pem ; r/em = (1 - s). (45, 46)

JP 2 + Q 2 Pem

V em ^em

Анализ расчета электромагнитных характеристик

Как следует из полученных выражений, безразмерные характеристики индукционного МГД-перемешивателя определяют следующие параметры:

1. Количество пазов M, размеры пазов 2Дт и их расположение xm, амплитудно-фазовые характеристики синусоидальных токов в пазах Ут=Уте'¥т, соотношение размеров индуктора (L=2pт=2т) и ванны (Т=Т/т)..

2. Безразмерный коэффициент добротности е=ую^дт2 /л2.

3. Относительная величина зазора между индуктором и расплавом - 5/т.

4. Скорость движения расплава (как твердого тела) относительно индуктора - u или величина скольжения 5=1-u / т/ Для упрощения сравнения характеристик будем полагать, что ширина всех пазов одинакова (Дт=соп5/) и они занимают половину длинны индуктора, т.е. Т!Мп=12кт =т.. Условие Т=2т позволяет анализировать электромагнитное поле и электромагнитные характеристики без учета влияния продольного краевого эффекта. Рассмотрим три типа обмоток (рис. 3): трехфазная обмотка с фазной зоной 120° (рис. 3а); двухфазная обмотка с фазной зоной 90° (рис. 36); трехфазная обмотка с фазной зоной 60° (рис. 3е). Очевидно, что линейная токовая нагрузка во всех трех вариантах обмоток одинакова.

Рис. 3. Исследуемые типы обмоток индуктора

a

0 0.3 113 2 0 0.5 1 5 2 0 05 15 2

a б в

Рис. 4. Относительная напряженность электрического поля и нормальная составляющая относительной магнитной индукции для трех вариантов обмоток

На рис. 4 представлены огибающие кривые относительной напряженности электрического поля (E1y / m^Jt) и нормальной составляющей относительной магнитной индукции (B1z /¡л J) на поверхности индуктора (z=0) - сплошные линии - и поверхности расплава (z=d') - пунктирные линии. По оси абсцисс показаны относительные значения. Кривые, представленные на рис. 4a, 46, 4в, соответствуют обмоткам на рис. 3a, 36, 3в. При расчете принималось Т=2т; s=1; 8/т=0,2; s=2, т.е. скорость расплава u=0. Если на поверхности индуктора зубцово-пазовые гармоники явно проявляются, то на поверхности расплава их проявление ослабевает. При этом с увеличением числа пазов проявление гармоник также уменьшается.

Картина распределения электромагнитного поля изменяется при наличии движения расплава относительно индуктора (пф0). На рис. 5 показаны те же огибающие при наличии движения. При расчете принималось Т= 6т; s=0,1; е=2; д/т=0,2. От начала координат до входа в индуктор предусмотрено расстояние т, а от выхода индуктора до конца периода Т - расстояние 3т. Вследствие проявления продольного краевого эффекта на входе индуктора электромагнитное поле сжимается, а на выходе вытягивается. Неравномерность распределения электромагнитного поля по длине индуктора приводит к несимметрии токов в фазах при их подключении к источнику фиксированного напряжения. С уменьшением зазора S неравномерность возрастает.

На рис. 6 и 7 представлены интегральные характеристики, соответствующие рассматриваемым обмоткам. Зависимости относительных активной Рет и реактивной Qem мощностей от е при M=3, 4, 6, показаны на рис. 6a и 66 соответственно. На рис. 6в показана зависимость

cos <р = /(е).

0 15 3 45 6 0 1j 3 4.5 6 0 15 J 4.5 6

0 1.5 3 45 6 о 1 5 3 4.5 6 0 :.5 3 4.5 6

а б в

Рис. 5. Относительная напряженность электрического поля и нормальная составляющая относительной магнитной индукции для трех вариантов обмоток при Тф2т и 5=0,1

о i.:s з.н j и ¡.ы ;.з з.н j "о S :.j i.is 'j

а б в

Рис. 6. Относительные активная (а), реактивная мощности (б) и коэффициент мощности от е (в)

На рис. 7а изображены зависимости относительной электромагнитной силы Еет от добротности е, а на рис. 7б представлены зависимости коэффициента полезного действия п от скольжения 5.

Как следует из представленных графиков, индуктор с трехфазной обмоткой с фазной зоной 120° (М=3) потребляет из сети большие активную и реактивную мощности, чем индуктор с обмотками М=4 и 6. Вместе с тем индуктор с этой обмоткой развивает электромагнитную силу Еет , действующую на расплав, меньшую, чем индукторы с другими обмотками (рис. 7а). Этому соответствуют графики \\=/5), представленные на рис. 7б. С увеличением величины

- 82 -

B)

Рис. 7. Относительная электромагнитная сила от е (сверху), КПД от 5 (снизу): а - д/т=0,2; б - д/т=0,3; в - ¿/т=0,4

зазора (д/т от 0,2 до 0,4) разность электромагнитных характеристик представленных обмоток уменьшается.

Выводы

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

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

Исследование выполнено при финансовой поддержке Российского фонда фундаментальных исследований, Правительства Красноярского края, Красноярского краевого фонда поддержки научной и научно-технической деятельности в рамках научных проектов № 16-48-242018 и № 16-43-242013р_офи_м.

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

[1] Khatsayuk M, Timofeev V, Michaylov D, Research of MHD influence during preparation and casting aluminum melt in electromagnetic field, Applied Mechanics and Materials, 2015, 698, 181187.

[2] Baake E, Jakovics A, Khatsayuk M, Lupi S, Timofeyev V, MHD technologies in metallurgy. Intensive Course Specific IV. Study guide, Publishing house of ETU, St. Petersburg, 2013. 232.

[3] Вольдек А.И. Индукционные магнитогидродинамические машины с жидкометалличе-ским рабочим телом, Ленинград, Энергия, 1970. 272 с. [Vol'dek A.I. The induction MHD machines with liquid metal working fluid, Leningrad, Energy, 1970, 272 p. (in Russian)]

[4] Гринберг Г. А. Избранные вопросы математической теории электрических и магнитных явлений, Л.: Издательство АН СССР, Москва, 1948. 727 с. [Grinberg G.A. Selected problems of mathematical theory of electrical and magnetic phenomena. Moscow, USSR Academy of Sciences, 1948, 727 p. (in Russian)]