CC BY
25
УДК 004.422.8 Дата подачи статьи: 14.12.16
DOI: 10.15827/0236-235X^30.3.537-545 2017. Т. 30. № 3. С. 537-545
КОМПЛЕКС ПРОГРАММ ДЛЯ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ ЭЛЕКТРИЧЕСКИХ ДУГ ТРЕХФАЗНОГО ПЕРЕМЕННОГО ТОКА, ГОРЯЩИХ НА ГОРИЗОНТАЛЬНУЮ ПОВЕРХНОСТЬ
И.М. Ячиков, д.т.н., профессор, [email protected]; Е.М. Костылева, аспирант, [email protected] (Магнитогорский государственный технический университет им. Г.И. Носова, просп Ленина, 38, г. Магнитогорск, 455000, Россия)
В металлургической промышленности эксплуатируется множество агрегатов, в которых в качестве источников тепловой энергии используются электрические дуги, горящие на горизонтальную тепловоспринимающую поверхность. Знание электромагнитных сил, действующих на дуги, их положение и поведение важны при конструировании новых печей и автоматизации их работы.
Обзор литературы показывает, что многие вопросы, связанные с электрическими дугами трехфазного переменного тока, освещены недостаточно полно, поскольку на промышленных или опытных установках изучение большинства проблем технически сложно, дорого и часто невыполнимо по самым разным причинам.
Возможным направлением исследования является создание адекватных математических моделей, их реализация в виде прикладного ПО и дальнейшее компьютерное моделирование.
В данной работе показаны математическое модели и разработанный комплекс программ, включающий в себя модули для определения напряженности магнитного поля, электромагнитных сил, формы дуг и теплового потока вблизи трех дуг трехфазного переменного тока, горящих на горизонтальную поверхность.
Результатом работы созданного программного комплекса являются расчет и визуализация векторных полей -напряженности магнитного поля и среднеинтегральных электромагнитных сил, а также визуализация средней за период формы осей столбов электрических дуг и распределение тепловых потоков по горизонтальной тепловосприни-мающей поверхности.
Ключевые слова: трехфазный переменный ток, напряженность магнитного поля, электромагнитные силы, форма оси столба дуги, тепловой поток, математическая модель, программное обеспечение.
В настоящее время существует множество агрегатов, в которых в качестве источников тепловой энергии используются трехфазные электрические дуги, горящие на горизонтальную поверхность. Это дуговые печи для выплавки металлов и сплавов, рудотермические, рудовосстановительные печи и другие. Знание ряда параметров электрических дуг, таких как магнитные поля, электромагнитные силы, действующие на дуги, положение и динамика поведения дуг, а также тепловые потоки излучения вблизи горящих дуг, является важным при эксплуатации существующих, конструировании новых печей и автоматизации их работы. На основе литературного анализа установлено, что многие вопросы, связанные с электрическими дугами трехфазного переменного тока, горящими на горизонтальную поверхность, освещены недостаточно полно. На промышленных или опытных установках исследование этих и других подобных проблем технически сложно, дорого и часто невозможно по самым разным причинам. Одним из возможных методов их изучения являются создание адекватных математических моделей, их реализация в виде ПО и дальнейшее исследование с использованием компьютерных моделей.
Для электрических дуг трехфазного переменного тока наиболее практически значимыми являются такие параметры, как напряженность магнитного поля в районе горения дуг, электромагнитные силы, действующие на дуги, форма дуг и радиаци-
онные тепловые потоки, падающие на горизонтальную поверхность.
В работах [1, 2] рассмотрена математическая модель для определения результирующей мгновенной напряженности магнитного поля в произвольной точке пространства A от токов, проходящих через три параллельных электрода. Находятся векторы напряженности магнитных полей, вызываемые токами в j-м электроде, а затем определяется результирующая сумма векторов с использованием принципа суперпозиции магнитных полей.
В выбранной полярной системе координат относительно центра распада электродов центры первого, второго и третьего электродов имеют координаты C\(rp, 0), C2(rp, 2л/3) и C3(rp, 4л/3) соответственно. Выбранная полярная система координат и совмещенная с ней декартовая система координат Oxy, а также схематическое расположение токов и точки А показаны на рисунке 1.
Напряженность магнитного поля НЛ в произвольной точке A(r, ф) определяется суммой трех векторов, координаты которых определяются по формуле
Haj = Haj ((-Xa + yc¡ ) • d, (Xa - xCJ ) • d) , (1) í0 • sin (2%vt + ф .)
где d =---^—--^-— ; fp - радиус
2 ^•(^уа + ycj ) +(xa - xcj ) )
окружности, на которой расположены токи; т -
время; г0 - амплитуда тока; ф/ = 2л(/-1)/3, ' = 1, 2,
3; V = 50 Гц - частота тока; Т = 1/у - период; И л] -напряженность магнитного поля в точке А от тока '-го электрода; А(ха, у0) - координаты точки А; С'(хд', уд) - координаты точек Су.
Результирующая напряженность магнитного поля в точке А от действия трех токов, проходящих через графитированные электроды
ил - ил ( Хи ! + хи 2 + Хнъ, ут + у и 2 + Уи 3 ) , (2)
где хщ, уну - координаты х и у вектора напряженности магнитного поля от действия '-го тока.
На основе математической модели в среде разработки и моделирования MathCad создан программный модуль, позволяющий определять магнитные поля вблизи трех параллельно горящих дуг трехфазного переменного тока. В качестве входных параметров задаются диаметр распада электродов, амплитудное значение тока дуг и полярные координаты точки, в которой ищем магнитное поле. Программа позволяет строить графики изменения напряженности магнитного поля в зависимости от безразмерного расстояния до центра распада электродов в различные моменты времени или при разных значениях угловой координаты, нормированные векторные поля напряженности магнитного поля вблизи трех дуг трехфазного переменного тока, годографы векторов напряженности магнитного поля. В качестве примера были построены нормированные векторные поля напряженности магнитного поля в различные моменты времени. Так, для моментов времени т=пТ/2, где п=0, 1, 2...,
магнитное поле вблизи первого электрода отсутствует. В последующие моменты времени оно постепенно увеличивается, принимает максимальное значение при т = Т/4 + пТ/2, после чего постепенно уменьшается и к моменту времени т = Т/2 + пТ/2 пропадает. Аналогично ведут себя напряженности магнитных полей вблизи второго и третьего электродов, возле которых поле отсутствует в моменты времени т = Т/6 + пТ/2 и т = Т/3 + пТ/2 и максимально т = 5Т/12 + пТ/2 и т = 7Т/12 + пТ/2 соответственно. Установлено, что вблизи трех электрических дуг, горящих на горизонтальную поверхность, имеет место бегущее магнитное поле. На рисунке 2 приведен пример работы программы, показаны нормированные векторные поля напряженности магнитного поля вблизи трех параллельно горящих дуг трехфазного переменного тока в разные моменты времени, а также годографы вектора напряженности магнитного поля внутри и вне диаметра распада электродов.
Для электрометаллургической теории и практики важными являются методы расчета локальной электромагнитной силы, действующей на элемент дуги, и среднеинтегральной силы, действующей на всю дугу. При постоянных токах для локальной силы, действующей на элемент дуги от тока, протекающего через жидкий металл, в [3] получено выражение
_ г / \ _ УрН-р__1_
Л - Л (г)- 4я(д + 21 2 + (д + 2)2 ,
а для среднеинтегральной электромагнитной силы [4, 5] зависимость
F =
^gVG 4я
ln
+ Д2 + I)(Д + I )
Афд +А)2
+I2 +I)
(3)
где Д = 3,2-3,6 см - глубина электрической оси то-копроводящего слоя от поверхности жидкого металла [6]; /0 - ток дуги; /1 - ток, протекающий через жидкий металл со стороны другой дуги; 1д - длина дуги; I - длина горизонтального отрезка, по которому протекает ток Л; г - расстояние от точки привязки дуги на горизонтальной поверхности до рассматриваемого элемента дуги.
В работе [4] получены выражения для электромагнитной силы, действующей на ток дуги /1 = /А со стороны трех фазных токов /ав, /вс, /са, протекающих через металл. Установлено, что модуль равнодействующей силы пропорционален 4312 ¡4 . С учетом действующего значения линейного тока /да - л/3/ / л/2 получено локальное среднее значение электромагнитной силы, действующей на единицу длины дуги:
/м ( Z ) =
^12
1
/
(Д + z 12 + ( д + z)2
(4)
CK
4*3
y
СУ! Ф>)
2r 2-2rp' я
C21 rp
о
x
1.6-
. ¡12>?J
C1 (Гр,o) К,o)
в)
г)
Рис. 2. Результаты работы программного модуля: нормированные векторные поля напряженности магнитного поля для разных моментов времени: a) t=T/7; б) т=Т/5; годографы векторов напряженности магнитного поля в точках, расположенных: в) внутри диаметра распада электродов, г) за пределами диаметра распада электродов
Fig. 2. The results of the software module: normalized vector fields of a magnetic field for different moments of time: a) t=T/7; б) t=T/5; intensity vector hodographs of a magnetic field at points located: в) inside a pitch circle diameter,
г) outside a pitch circle diameter
и среднеинтегральное значение силы, действующей на дугу:
г 2
F ln
м 8^3-
(V l2 +А2 +1 )(Д + /д )
А фд +А )2 +12 +1 )
(5)
В работе [7] рассмотрено взаимодействие двух вертикальных дуг постоянного тока, получена локальная электромагнитная сила, действующая на единицу длины дуги, на ток столба дуги /0 со стороны тока /1, протекающего через дугу (длиной 1д)
и графитированный электрод (длиной 4д):
( \
di ~ f д -
МоН- о
4-l
L + К„ - г
Vl2 + ( 1д + 1эд - г)2
и среднеинтегральная сила, действующая на дугу:
F -
H-о¥о 4-l
L + L. - г
dz .
Ф2 + ^2 V2 + (1д + 'эд - ^ ,
Среднюю результирующую силу взаимодействия одной из дуг в трехфазной дуговой печи переменного тока с двумя другими дугами с учетом действующего значения линейного тока I = I / -\/2 и локальную электромагнитную силу,
действующую на единицу длины дуги [4, 7], можно определить следующим образом:
II
Fd ---^Г--1Wl2 +(1д + 1эд)2
/д ( z ) -
4-l
2
о1 дл
4-1
L + lm - z
Ф2 + z2 Vi2 + ( 1д + 1эд - z)2
,(6)
. (7)
4
3
На основе разработанных математических моделей в среде разработки и моделирования Ма1ЬаЬ создана компьютерная программа «Расчет электромагнитных сил, действующих на дуги в трехфазной дуговой печи», позволяющая рассчитывать основные электромагнитные силы (Свид. о гос. регистр. № 2016618499), ее главное окно показано на рисунке 3. Исходными данными для программы являются расстояние между дугами, длина дуги, длина электрода, действующее значение тока. Все они проходят проверку на корректность ввода [8]. Далее рассчитываются значения среднеинтегральных
сил от токов, протекающих по жидкому металлу, и от токов, протекающих через параллельные дуги и электроды. Значения этих сил одинаковы для каждой из дуг и различаются только направлениями. Далее при помощи методов аналитической геометрии визуализируются расположение электродов, среднеинтегральных электромагнитных сил и их результирующее направление. Также в программе имеется возможность построения графиков локальных и среднеинтегральных электромагнитных сил. При нажатии на кнопку «Графики» на экране появляется рабочее окно для вывода графиков (рис. 4).
а)
б)
Рис. 3. Программа «Расчет электромагнитных сил, действующих на дуги в трехфазной дуговой печи»: а) главное окно, б) результаты моделирования основных векторов среднеинтегральных электромагнитных сил, действующих на дугу
Fig. 3. A program "Calculation of electromagnetic forces affecting arcs in a three-phase arc furnace ": a) the main window, б) simulation results of the main vectors of average integral electromagnetic forces affecting the arc
Рис. 4. Рабочее окно программы «Расчет электромагнитных сил, действующих на дуги в трехфазной дуговой печи» для вывода графиков локальных и среднеинтегральных электромагнитных сил
Fig. 4. The working window of the program "Calculation of electromagnetic forces affecting arcs in a three-phase arc furnace "for the charts of local and average integral electromagnetic forces
Строится график локальных сил, представляющий собой зависимость ¥ /12 - значения локальной
дл
силы, деленной на квадрат тока от координаты точки г, и график среднеинтегральных сил ¥ / /2л
от безразмерной координаты точки ///д.
Рассмотрим математическую модель формы параллельно горящих на горизонтальную поверхность электрических дуг постоянного или переменного тока при их электромагнитном взаимодействии под действием силы ¥д . Для определения
формы двух дуг постоянного тока при их электромагнитном взаимодействии использовалась математическая модель, предложенная в работе [9]. Форма оси первой дугиу(Г) (рис. 5), через которую протекает ток /1, при электромагнитном взаимодействии с дугой, по которой протекает ток /2, описывается следующим интегро-дифференциальным уравнением:
. 1 '■2 J
Г J
dy (Ç) d Ç
(Ç-z )-(-y (Ç)-y ( z ))
(2 \ 3/ 2
(|y (Ç)| + У (z )) +(Ç-z )2 )
( "fVz )-( У (0-У ( z ))' (( У (0 + y ( z ))2 +(Ç- z )2 + y 2 )3
d Ç -
-J
-d ç-
Ц oß( z )
4л (l + (dy / dz)2 )
d y
dz
= 0.
(8)
где $(г) - функция, определяющая параметры дуги; уо - среднее геометрическое расстояние площади сечения дуги в ее цилиндрической части от самой себя.
Уравнение (8) относится к классу интегро-диф-ференциальных уравнений. Свойства таких уравнений и методы их решения рассмотрены в работах [1о-16]. В отдельных частных случаях найдены аналитические решения, но в общем случае они мо-
Рис. 5. Основные размеры, используемые в математической модели взаимодействия дуг
Fig. 5. The basic dimensions in an arc interaction mathematical model
гут быть решены только приближенно. Численное решение получается в результате итерационного процесса, точные оценки погрешности получены не для всех случаев. При численных расчетах в качестве критерия для завершения итерационного процесса применяется малое различие результатов, полученных на последовательных шагах процесса. Задача сводится к нахождению функции, которая с приемлемой точностью удовлетворяет уравнению (8). Решение затрудняется тем, что оно нелинейное, так как искомая функция входит в знаменатель правой части и в дифференциальном уравнении имеется модуль.
Для решения подобных интегро-дифференци-альных уравнений обычно используются приближенные итерационные алгоритмы, например метод Пикара [12]. Однако его численная реализация показала, что для описываемой задачи полученный итерационный процесс является расходящимся.
Для нахождения приближенного решения этого уравнения предложены и реализованы три алгоритма [17]. По первому алгоритму предполагалось, что по дугам текут равные постоянные токи i1 = i2 (дуги зеркально симметричны относительно оси Oz), а неизвестная функция положения оси первой дуги имеет форму параболы y(z) = aiz2 + biz + c.
dy
Применяя начальные условия >>(0) = a, —
dz
где а - тангенс угла наклона оси дуги в ее начальном сечении, получено c = a, b1 = а. Решая нелинейное уравнение (8) при разных z, находился неизвестный коэффициент a1, который затем усреднялся.
По второму алгоритму также рассматривался случай одинаковых токов i1 = i2, при этом искомая функция имела вид полинома n-й степени y(z) = = a0 + a1z + a2z2 + a3z3 + ...+ a„z". Неизвестные коэффициенты a2, a3, ... an находились как минимум
функции нескольких переменных: 91(a2, a3, ...) = i i = JФ2dz = min или 92(a2, a3, ...) =J|ф|dz = min, 0 0 где Ф - левая часть уравнения (8).
По третьему алгоритму рассматривалось нахождение формы осей двух дуг при разных токах, проходящих через них, i1 Ф i2. В этом случае для нахождения двух неизвестных функций решалась система двух дифференциальных уравнений (8), записанных для каждой дуги, с граничными условиями:
dy,
yi(0) = a, -f1 dz
dy
= ai, У2(0) = - a, — dz
= a,
Для получения приближенного решения предполагалось, что искомые функции формы осей столбов дуг имеют вид уДТ) = а' + а' + ауТ + + а3,]2ъ +... + ап/, где ' = 1, 2 - номер дуги. Из начальных условий, принимая а,1 = а2 = 0, имеем а11 = а12 = о ао1 = ао2 = о. Искомые коэффициенты
z =0
z 2 = 0
многочленов a2j = a3g = ... = arg находились как минимум функции нескольких переменных:
i
Ф3 = ДФf + Ф^) dz = min, где Ф1, Ф2 - левые части
0
уравнения (8), записанные для первой и второй дуг соответственно.
Уравнение равновесия элемента столба дуги расширено и применено для нахождения формы осей двух или трех дуг постоянного или трехфазного переменного тока и может быть записано как
к. dF. , + dF. . + dF. = О,
1 1-2 1-1 1-м ?
(9)
где к\ - коэффициент, учитывающий вид электромагнитного взаимодействия; с1Р1 2 - сила взаимодействия дуги с токами, протекающими по другим дугам; - сила взаимодействия дуги с собственным током; с!Р^ - центростремительная сила, выпрямляющая дугу. Значение коэффициента к1 и характер электромагнитного взаимодействия зависят от рода тока и количества дуг [17].
На основе разработанных математических моделей в среде разработки и моделирования Ма1ЬаЪ создан программный модуль «Взаимодействие дуг», позволяющий в зависимости от количества дуг, рода и значения протекающего по ним тока, расстояния между дугами, физических свойств плазмообразующего газа рассчитывать и визуализировать форму дуг (Свид. о гос. регистр. № 2016618499); главное окно программы показано на рисунке 6. Исходными данными являются расстояние между дугами, длина дуги, токи первой и второй дуг, степень полинома, при помощи которого аппроксимируется форма дуги, род тока, участвующего во взаимодействии дуг. Далее с учетом физических свойств плазмообразующего газа и
Рис. 6. Рабочее окно программы «Взамодействие дуг» для расчета и построения формы осей столбов дуг
Fig. 6. The working window of the program "Arc interaction " for calculating and constructing shapes of arc axes
свойств окружающей среды находятся коэффициенты для построения формы оси столба дуги. Программа строит график формы осей столбов дуг у(2) в сантиметрах.
В дуговой сталеплавильной печи (ДСП) трансформация электрической энергии в тепловую происходит в разрядном промежутке. За счет электромагнитного взаимодействия электрических токов, протекающих по графитированным электродам и в поверхностных слоях расплава с током дуги, происходит ее выдувание от центра печи. Средний угол наклона столба дуги 9 по отношению к нормали зеркала ванны колеблется для ДСП различной мощности и емкости в пределах 45°-65° (см., например, [18, 19]).
В работе [20] была получена модель теплового потока на произвольно расположенную площадку ёБ, находящуюся на поверхности ванны от электрической дуги, расположенной под углом к ней:
Чизл (r> A) = ■
К • P •
2 л • l • r
- B + A • arcth
B
BJ
(10)
B • r(lA + r) - A(r2 + 2A • r • l +12 )arcth
1д + rA rB
B (r2 + 2A• r• 1д + 1Д2)
где А = соБу-Бшб; В = V А2 -1 ; г - расстояние от начала координат до площадки ёБ; у - угол, под которым расположена площадка ёБ к началу координат; 9 - угол между направлением излучения и нормалью к поверхности ванны; Кэ - поправочный экспериментальный коэффициент, учитывающий долю мощности дуги, идущую на излучение; Р -полная мощность излучения дуги; 1д - длина дуги (рис. 7).
Рис. 7. К расчету плотности теплового излучения дуги на площадку dS
Fig. 7. To the calculation of an arc thermal radiation density on dS area
+
Тепловой поток от трех наклонных горящих электрических дуг длиной /д на произвольно расположенную горизонтальную площадку на -ходится по принципу суперпозиции как сумма тепловых потоков, падающих в данную точку А (рис. 8). Тепловой поток, падающий на элементарную площадку, определяется как
4Мл К а. 0) = Z qj
j=1
■Ч изл (Tj. V j. 0), (11)
где д шл(г , V/, 0) - тепловой поток, падающий на
горизонтальную поверхность от у-й дуги; гу - расстояние оту-й дуги до площадки ёБ; - угол, под которым расположена площадка ёБ ку-й дуге; 0 -угол между направлением излучения и нормалью к поверхности ванны [21].
Vv. , У
С2\\ A(a,a)
Vi
) 0 С x
Сз А
Рис. 8. К расчету плотности теплового излучения дуги на площадку dS от трех дуг
Fig. 8. To the calculation of an arc thermal radiation density on dS area from three arcs
С использованием созданной математической модели (10)-(11) разработан компьютерный модуль, позволяющий рассчитывать и визуализировать тепловые потоки, падающие на горизонтальную поверхность, вблизи горения трех наклонных электрических дуг. В качестве входных параметров задаются диаметр распада электродов, действующее значение тока, координаты точек, в которых ищется и визуализируется тепловой поток. На рисунке 9 изображено распределение теплового потока по азимутальной координате при разном расстоянии от центра распада электродов в трехфазной дуговой печи ДСП-180 (ток дуги - 50 кА, /д = 0,5 м, диаметр распада электродов - 1,2 м ), построенное при помощи разработанного модуля.
Таким образом, на основе оригинальных математических моделей разработан комплекс компьютерных программ для определения основных параметров электрических дуг трехфазного переменного тока. Комплекс программ позволяет рассчитать напряженность магнитных полей, опре-
делить и визуализировать среднеинтегральные электромагнитные силы, действующие на дуги, и локальные силы, действующие на элементы дуги: от токов, протекающих через жидкий металл, и от токов, проходящих через параллельно горящие дуги и графитированные электроды. Также специальные модули позволяют рассчитать и визуализировать формы осей столбов дуг в зависимости от их количества и рода тока и для определения и визуализации теплового потока от трех горящих наклонных электрических дуг по азимутальной координате. Структурная схема комплекса программ представлена на рисунке 10. Разработанный комплекс компьютерных программ может использоваться при конструировании многоэлектродных
3
Рис. 10. Структурная схема комплекса программ Fig. 10. The software complex structural diagram
дуговых или плазменных печей для определения положения электрических дуг в системах автоматизированного регулирования их тепловой работы, для оценки стойкости футеровки вблизи горящих дуг при разных технологических режимах, для определения воздействия электромагнитных полей на приборы и человека.
Литература
1. Леушин А.И. Дуга горения. М.: Металлургия, 1973. 240 с.
2. Ячиков И.М., Зарецкая Е.М. Анализ поведения магнитного поля вблизи электродов дуговых печей посредством математического моделирования // Изв. вузов: Черная металлургия. 2011. № 1. C. 18-20.
3. Миронов Ю.М. Электрическая дуга в электротехнологических установках: монография. Чебоксары: Изд-во Чуваш. ун-та, 2013. 290 с.
4. Yachikov I.M., Kostyleva E.M. Electromagnetic forces on the arc in a three-phase arc furnace. Steel in Translation. 2015. vol. 45, no. 7, pp. 467-472. DOI: 10.3103/S0967091215070141.
5. Егоров А.В. Расчет мощности и параметров электроплавильных печей. М.: Изд-во МИСиС, 2000. 272 с.
6. Егоров А.В. Электроплавильные печи черной металлургии. М.: Металлургия, 1985. 280 с.
7. Никольский Л.Е., Смоляренко В.Д., Кузнецов Л.Н. Тепловая работа дуговых сталеплавильных печей. М.: Металлургия, 1981. 320 с.
8. Лазарев Ю. Моделирование процессов и систем в MATLAB. СПб: Питер, 2005. 512 с.
9. Бортничук Н.И., Крутянский М.М. Плазменно-дуговые плавильные печи. М.: Энергоиздат, 1981. 120 с.
10. Азбелев Н.В., Максимов В.П., Рахматуллина Л.Ф. Элементы современной теории функционально-дифференциальных уравнений. Методы и приложения. Ижевск: Изд-во ИКИ, 2002. 304 с.
11. Дмитриев С.С., Кузнецов Е.Б. Численное решение систем интегро-дифференциально-алгебраических уравнений с запаздывающим аргументом // Журн. вычислит. матем. и матем. физики. 2008. № 3. С. 430-444.
12. Вержбицкий В.М. Численные методы. Математический анализ и обыкновенные дифференциальные уравнения. М.: Высш.шк., 2001. 383 с.
13. Федоренко Р.П. Введение в вычислительную физику. Долгопрудный: Интеллект, 2008. 504 с.
14. Зализняк В.Е. Основы научных вычислений. Введение в численные методы для физиков. М.: Editorial URSS, 2002. 296 с.
15. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы М.: Бином. Лаборатория знаний, 2003. 632 с.
16. Лёзин И.А. Решение систем полиномиальных уравнений на ЭВМ // Программные продукты и системы. 2012. № 3. С. 22-25.
17. Ячиков И.М., Костылева Е.М. Математическое моделирование формы дуг при их электромагнитном взаимодействии. Сообщ. 2: Форма электрических дуг постоянного или трехфазного переменного тока, горящих между электродами и токопод-водящей поверхностью // Изв. вузов: Черная металлургия. 2014. № 5. С. 56-61.
18. Макаров А.Н. Теплообмен в дуговых сталеплавильных печах. Тверь: Изд-во ТГТУ, 1998. 184 с.
19. Миронов Ю.М., Михадаров Д.Г. Влияние межэлектродного промежутка на отклонение столба дуги // Вестн. ЧГУ. 2009. № 2. С. 139-142.
20. Yachikov I.M., Zaretskaya E.M. Radiant heat flux from an electric arc inclined to the plane of the receiving surface. Steel in Translation, 2011, vol. 41, no. 7, pp. 550-554. DOI: 10.3103/S0967091211070151.
21. Макаров А.Н., Свенчанский А.Д. Оптимальные тепловые режимы дуговых сталеплавильных печей. М.: Энергоатом-издат, 1992. 96 с.
Software & Systems Received 14.12.16
DOI: 10.15827/0236-235X.030.3.537-545 2017, vol. 30, no. 3, pp. 537-545
A SOFTWARE SUITE TO DETERMINE THE PARAMETERS OF THREE-PHASE ALTERNATING CURRENT ELECTRIC ARCS BURNING IN A HORIZONTAL SURFACE
I.M. Yachikov 1, Dr.Sc. (Engineering), Professor, [email protected] E.M. Kostyleva l, Postgraduate Student, [email protected]
1 G.I. Nosov Magnitogorsk State Technical University), Lenin Ave. 38, Magnitogorsk, 455000, Russian Federation
Abstract. The metallurgical industry uses a lot of equipment with electrical arcs as the sources of thermal energy. They burn in the direction of the horizontal heat-absorbing surface. Knowing electromagnetic forces affecting the arcs, their location and behavior is significant when designing new furnaces and automating their operation.
The literature review shows that many issues connected with electrical arcs of three-phase alternating current are under-served due to technical complexity, high cost and often impossible on industrial or pilot plants.
A possible method of research is the development of appropriate mathematical models, their implementation as applied software and further computer modeling.
The paper shows mathematical models and developed software, which includes modules for determining the magnetic field intensity, electromagnetic forces, arc and heat flow forms, near three arcs of three-phase alternating current burning in the direction of a horizontal surface.
The result of the created software complex is calculation and representation of vector fields including magnetic field strength and average-integral electromagnetic forces, as well as visualization of the average form of electric arc column axis by the period and the heat flux distribution on a horizontal heat-absorbing surface.
Keywords: three-phase alternating current, magnetic field, electromagnetic forces, form of arc column axis, heat flux, mathematical model, software.
References
1. Leushin A.I. Duga goreniya [The Burning Arc]. Moscow, Metallurgiya Publ., 1973, 240 p.
2. Yachikov I.M., Zaretskaya E.M. The analysis of a magnetic field behavior near electric arc furnace electrodes using mathematical modeling. Izvestiya vysshikh uchebnykh zavedeny. Chernaya metallurgiya [Izvestiya. Ferrous Metallurgy]. 2011, no. 1, pp. 18-20 (in Russ.).
3. Mironov Yu.M. Elektricheskaya duga v elektrotekhnologicheskikh ustanovkakh [An Electric Arc in Electrotechnolog-ical Installations]. Monograph. Cheboksary, Chuvash. un-t Publ., 2013, 290 p.
4. Yachikov I.M., Kostyleva E.M. Electromagnetic forces on the arc in a three-phase arc furnace. Steel in Translation. 2015, vol. 45, no. 7, pp. 467-472.
5. Egorov A.V. Raschet moshchnosti i parametrov elektroplavilnykh pechey [The Calculation of Power and Parameters of Electric Smelting Furnaces]. Moscow, MISiS Publ., 2000, 272 p.
6. Egorov A.V. Elektroplavilnyepechi chernoy metallurgii [Electric Melting Furnace for Ferrous Metallurgy]. Moscow, Metallurgiya Publ., 1985, 280 p.
7. Nikolsky L.E., Smolyarenko V.D., Kuznetsov L.N. Teplovaya rabota dugovykh staleplavilnykhpechey [Thermal Performance of Electric Arc Furnaces]. Moscow, Metallurgiya Publ., 1981, 320 p.
8. Lazarev Yu. Modelirovanie protsessov i sistem v MATLAB [Modeling Processes and Systems in MATLAB]. St. Petersburg, Piter Publ., 2005, 512 p.
9. Bortnichuk N.I., Krutyansky M.M. Plazmenno-dugovye plavilnye pechi [Plasma-arcs Melting Furnace]. Moscow, Energoizdat Publ., 1981, 120 p.
10. Azbelev N.V., Maksimov V.P., Rakhmatullina L.F. Elementy sovremennoy teorii funktsionalno-differentsialnykh uravneny. Metody iprilozheniya [Elements of the Modern Theory of Functional Differential Equations. Methods and Applications]. Moscow, IKI Publ., 2002, 304 p.
11. Dmitriev S.S., Kuznetsov E.B. Numerical solution to systems of delay integrodifferential algebraic equations. Zhurn. vychislitelnoy matematiki i matematicheskoy fiziki [Computational Mathematics and Mathematical Physics]. 2008, vol. 48, no. 3, pp. 406-419.
12. Verzhbitsky V.M. Chislennye metody. Matematichesky analiz i obyknovennye differentsialnye uravneniya [Numerical Methods. Mathematical Analysis and Ordinary Differential Equations]. Moscow, Vyssh. shk. Publ., 2001, 383 p.
13. Fedorenko R.P. Vvedenie v vychislitelnuyufiziku [An Introduction to Computational Physics]. 2nd ed., Intellekt Publ., 2008, 504 p.
14. Zaliznyak V.E. Osnovy nauchnykh vychisleny. Vvedenie v chislennye metody dlya fizikov [Fundamentals of Scientific Computing. Introduction to Numerical Methods for Physicists]. URSS Publ., 2002, 296 p.
15. Bakhvalov N.S., Zhidkov N.P., Kobelkov G.M. Chislennye metody [Numerical Methods]. Moscow, Binom. Labora-toriya znany Publ., 2003, 632 p.
16. Lezin I.A. Solving simultaneous polynomial equations using a computer. Programmnye produkty i sistemy [Software & Systems]. Tver, 2012, no. 3, pp. 22-25 (in Russ.).
17. Yachikov I.M., Kostyleva E.M. The mathematical modeling of arcs' form during their electromagnetic coupling. Message 2. Izvestiya vysshikh uchebnykh zavedeny. Chernaya metallurgiya [Izvestiya. Ferrous Metallurgy]. 2014, vol. 57, no. 5, pp. 56-61.
18. Makarov A.N. Teploobmen v dugovykh staleplavilnykhpechakh [Heat Exchange in Arc Steel-Melting Furnaces]. Tver, TSTU Publ., 1998, 184 p.
19. Mironov Yu.M., Mikhadarov D.G. The influence of the interelectrode gap on the deviation of an arc column. Vestn. ChGU [Cherepovets State Univ. Bulletin]. 2009, no. 2, pp. 139-142 (in Russ.).
20. Yachikov I.M., Zaretskaya E.M. Radiant heat flux from an electric arc inclined to the plane of the receiving surface. Steel in Translation. 2011, vol. 41, no. 7, pp. 550-554.
21. Makarov A.N., Svenchansky A.D. Optimalnye teplovye rezhimy dugovykh staleplavilnykh pechey [Optimum Thermal Conditions of Electric Arc Furnaces]. Moscow, Energoatomizdat Publ., 1992, 96 p.
