CC BY
72
ЭЛЕКТРОТЕХНОЛОГИИ
УДК 669.187.013.621
ИССЛЕДОВАНИЕ ЭЛЕКТРОДИНАМИЧЕСКИХ СИЛ, ДЕЙСТВУЮЩИХ НА ДУГИ В ТРЁХФАЗНОЙ ДУГОВОЙ СТАЛЕПЛАВИЛЬНОЙ ПЕЧИ
Н.В. Евсеева, Л.А. Лазуко, Ю.Б. Черкасова, С.У. Хасанов г. Челябинск, ОАО «Челябинский металлургический комбинат»
RESEARCH OF ELECTRODYNAMIC FORCES INFLUENCING THE ARCS OF THREE-PHASE EAF
N. Yevseyeva, L. Lazuko, Y. Cherkasova, S. Khasanov Chelyabinsk, Joint Stock Company “Chelyabinsk Metallurgical Plant”
Приведён принцип расчёта электродинамических сил, действующих на электрические дуги в трёхфазной ДСП, с построением графического изображения этих сил. Дано объяснение экстремальному характеру зависимости потребляемой из сети активной мощности от тока. Доказано, что образование мениска под дугой наблюдается только при токах, соответствующих нисходящей ветви зависимости мощности дуги от тока. Визуализация дуги позволяет выбрать оптимальную рабочую точку для любого периода плавки.
Ключевые слова: электродинамические силы, визуализация (графическое изображение), газовая струя ореола, столб дуги, восходящая и нисходящая ветвь зависимости мощности дуги от тока.
The article describes the calculation principle for electrodynamic forces influencing the electric arcs of a three-phase EAF and their graphic representation. Extreme type of the relationship between current and active power consumed from the network is explained. The article proves that metal meniscus is formed under the arc when the currents correspond with descending branch of graphic chart showing the relationship between current and arc power only. Graphic representation of the arc allows choosing the effective operating point for any period of melting process.
Keywords: electrodynamic forces, imaging (graphic representation), circle of flame efflux, arc column, ascending and descending branches of the graphic chart showing the relationship between current and arc power.
Исследование электродинамических сил, действующих на электрическую дугу в трёхфазной дуговой сталеплавильной печи (ДСП), выполнено по методике, применяемой в теории электрических аппаратов, и изложено в работах [1-4].
Все электродинамические силы, действующие на дуги в трёхфазных ДСП, имеют чёткую математическую интерпретацию, и при современном состоянии программирования и вычислительной техники, применяемой на ДСП, могут быть изображены графически для действующей печи в реальном или ином масштабе времени. Величина и направление электродинамических сил (э.д.сил) определяют поведение электрической дуги, её столба и ореола. При наличии в ДСП холодной шихты магнитная проницаемость среды ц возрастает в десятки раз (на жидкой ванне ц=1), э.д.силы,
пропорциональные величине ц, увеличиваются настолько, что выдувают не только газовую оболочку дуги (ореол), но и сам столб дуги. Поскольку в трёхфазной ДСП э.д.силы вращаются и направлены наружу в сторону стен печи, то в процессе расплавления шихты образуются колодцы, диаметр которых зависит от длины дуги и величины э. д.сил, определяемых магнитной проницаемостью ц, т. е. температурой шихты. В процессе работы печи электрический режим подбирается таким образом, чтобы проплавляемые колодцы имели форму перевёрнутого конуса.
В ходе расплавления шихты образуется жидкая ванна и магнитная проницаемость среды даже при наличии нерасплавленных, но нагретых выше точки Кюри (700 °С) кусков шихты, снижается до д=1. Э.д.силы также снижаются по величине и
даже при достаточно больших токах не могут изменить положение столба дуги, но ореол дуги сдувается этими силами в сторону стен печи и вверх, а при токах, больших максимального значения, - в сторону стен печи и вниз.
В характеристике Рд(1) зависимости мощности дуги Рд от тока дуги I будем выделять восходящую ветвь - до достижения критического значения тока 1КР, экстремальный (весьма «размытый») участок и нисходящую ветвь - до тока короткого замыкания.
Визуальное наблюдение за дугой на печах от 8 МВА до 80 МВА показало, что при работе на восходящей ветви зависимости Рд(1) не наблюдается мениска под дугой.
При токе 1кр происходит лавинообразное увеличение сечения дуги вследствие перехода дуги из режима работы с катодным пятном в режим работы без катодного пятна. При уменьшении тока происходит внезапный обратный переход к режиму с катодным пятном, однако при более низком значении тока, чем 1КР [6].
Как установлено в работе [4], плотность тока в катодном пятне для всех ДСП постоянна и равна около 82А/см2. При резком увеличении сечения дуги и расплывании её по всему сечению электрода плотность тока, естественно, падает.
Режим работы на нисходящей ветви характеристики Рд(1) широко применялся на 100-тонных ДСП, выпускающих готовые марки стали, и рекомендован институтом металлургии (НИИМ) для периода доплавления ^ф > 1).
При лавинообразном увеличении катодного пятна электрическое сопротивление столба дуги Гд длиной 1д и сечением 8д падает: Гд=Рд1д/8д, где рд - удельное электрическое сопротивление дуги, рд=1,1710-30мм [4, с. 49].
Определим величину и направление э.д.сил при токах, соответствующих восходящей и нисходящей ветви характеристики Рд(1).
На дугу первой фазы (фА) действуют э.д.силы:
- со стороны параллельных проводников двух других фаз р!/2||, Рш||;
- со стороны токов, протекающих в жидкой ванне по сторонам треугольника, Р1/2д, Р1/3д;
- аксиальная электродинамическая сила, возникающая при взаимодействии тока дуги с собственным магнитным полем, Ракс.
Запишем мгновенные значения этих сил:
Р1/2ц=2 11 12 ^10-7зт(ю1;) зт(ю1;-2/зя;) к1/2ц, Н,
Р1/3 ц =2 11 13 ^10-7зт(ю1;) зт(ю1+2/зл) к1/3 ц, Н, где 11, 12, 13 - действующие значения фазных токов, к1/2||, к1/3ц - коэффициенты контура сил Р1/2||, Р1/31 [6, с. 14].
Вычисление коэффициентов контура подробно рассмотрено в [1].
Суммарная сила Р1 ц, действующая на дугу фазы 1 со стороны параллельных проводников фаз 2 и 3, равна векторной сумме
1
г 1/2II
1/3
Электродинамические силы, действующие на дугу фазы 1 со стороны условных проводников 2 и 3 (рис. 1), образованных токами в ванне ДСП, равны:
Р1/2Л=М- 1А 1АВ 1 0 7k1/2A, H,
Р1/3Л=^ 1А 1СА 1 0 7k1/3A, H,
где 1А, 1АВ, 1СА - действующие значения токов, к1/2д, к1/3д - коэффициенты контуров для проводников 1-2 и 1-3 [1].
Рис. 1. Направление э.д.сил F1/2д, F2/1д, F1/з4, Fз/14 для положительных значений и, ІAB, ІCA
Мгновенные значения токов:
I— л/2"
и= л/2 1А зт(ю1;), иВ= 1А вт(ю1+я/6),
л/3
А
л/3
іСА= —^ ІА зіп(ю1+5/6 п)
Мгновенные значения сил :
2
Р1/2Д=_/=ІА2' 10 7к1/2дзіп(ю1:) зіп(ю1+п/6),
л/3
2 2 _7
Рі/3А= ^ * ІА2'10 кі/зд зіп(ю1:) зіп(ю1+5/6 п).
л/3
Суммарная сила Р1Д, действующая на дугу фазы 1 со стороны токов, протекающих в ванне печи по сторонам треугольника, равна векторной сумме сил Р1/2Д, Ршд:
Рід=Р,/2Д+Р,/ЗД .
На рис. 1 показано положительное направление векторов тока ІА, ІАв, ІСА. В соответствии с [6, с. 19], если ток проходит через угол от одного проводника к другому, то силы, действующие на эти проводники, направлены наружу (силы Р1/2Д, Р2/1Д,). Если токи направлены навстречу, к углу (как ІА и ІСА на рис. 1), то силы, действующие на проводники 1 и 3, стремятся уменьшить угол и направлены внутрь (силы Р1/3Д, Р3/1Д на рис. 1).
Аксиальная электродинамическая сила возникает в результате взаимодействия тока дуги с соб-
Исследование электродинамических сил, действующих на дуги...
ственным магнитным полем и определяется по формуле [2]:
F^m-8 I2 ц ln
Н,
где I - действующее значение тока дуги, А; д - магнитная проницаемость среды; Бэ - площадь сечения электрода, см2; Бр - площадь растекания тока в жидкой ванне, см2.
Площадь растекания тока Бр в жидкой ванне печи пропорциональна сечению дуги Бр=кр Бд, где кр - коэффициент растекания тока.
Статистическое сечение дуги (в см2) при токе I < 1КР может определяться по формуле Бд=12,2 I, где I - ток дуги, кА [7].
В общем случае сечение дуги определяется по формуле
БД=рД 1д /Гд, м2, где Рд - удельное сопротивление плазмы дуги, Рд=1,17- 10-3 Омм [7]; 1д - длина дуги, м.
Приняв, что градиент падения напряжения вдоль вертикальной оси дуги равен 1 В/мм [9, с. 49], получим эмпирическое соотношение: длина дуги 1д (мм) численно равна напряжению дуги ид (В). Тогда
Бд=11,7 Цд/гд (см2).
Для агрегата комплексной обработки стали (АКП) с печным трансформатором (ПТ) 32 МВА для 9-й ступени (Цст=306 В) принято кр=3,2 и формула для определения аксиальной силы примет вид
F акс=5^10-8 !д2 ц ln
8Э гд
3,2-11,7 U
д
Н.
Запишем эту формулу для мгновенного значения тока дуги 1д= >/2 Iд 8шю1:
(
Faкс=10 7 Iд2 sin2(rat) ц ln
8Э гд
л
3,2-11,7 Uц
v ’ ’ д /
, Н.
При токах Iпр, когда Бэ=Бр, выражение под знаком 1п равно 1, а Р^ =0. При токах I > ^р аксиальная сила меняет знак, становится отрицательной и направлена вниз. Результирующая э.д.сил, действующая на дугу первой фазы Рь равна векторной сумме сил Р1 | , Р1Д , Р !акс:
Р = Р1| + Р1Д + Р1акс .
Угол поворота результирующей силы Р^шг) относительно оси Ъ равен
(
11акс F1
Л
ч |^1акс | • |Fli
При cos ^1>0 результирующая сила F1 направлена вверх и в сторону стен печи, при cos ^1=0 сила F1 перпендикулярна вертикальной оси и выдувает ореол дуги в сторону стен печи, при cos ^1<0 сила F1 направлена вниз и в сторону стен печи.
На базе рассмотренного выше анализа электродинамических сил разработана программа рас-
чёта и визуализации этих сил в трёхмерном пространстве. Графическое изображение э.д.сил позволяет получить представление о положении ореола дуги и её столба по ходу плавки в реальном масштабе времени, так как непосредственное наблюдение за дугой на мощных ДСП с использованием газокислородных горелок просто невозможно. Ориентируясь на графическое изображение э.д.сил технологи могут выбрать наиболее рациональный режим работы ДСП. В программу можно ввести момент изменения ц магнитной проницаемости среды в зависимости от потребленной электроэнергии или связав уменьшение ц с исчезновением высших гармоник. Установка программы визуализации дуги особенно желательна на агрегатах ковш-печь, так как обработка металла на них ведётся, как правило, на нисходящей ветви характеристики Рд(1). Эта программа позволяет рассчитывать и отображать силы как динамически - в реальном масштабе времени по текущим значениям тока, напряжения и положения переключателя ступеней напряжения, так и статически - по заданным пользователем значениям. Переключение режимов работы и задание необходимых параметров производится в специальном графическом интерфейсе.
Программный продукт состоит из двух приложений: программы расчёта сил, действующих на дугу, «ForceWind» и программы визуализации векторов сил «ForceControl».
Программа «ForceWind» написана на языке Watcom C v10.6. Программа «ForceControl» реализована в пакете визуального программирования Photon Application Builder графической оболочки Photon v1.14. Обе программы работают в среде операционной системы реального времени QNX 4.25.
В качестве примера приводим расчёт э.д.сил в статическом режиме программы по данным пользователя, соответствующим экспериментальным характеристикам печи.
На рис. 2 приведены электрические характеристики дуговой печи с ПТ 32 МВА (агрегат ковш-печь, АКП). Зависимости активной мощности Р0 и реактивной мощности Q0 от тока сняты экспериментально, зависимости напряжения дуги ид и активного сопротивления гд от тока получены расчётным путём. Управляющий комплекс АКП содержит программу расчёта всех характеристик дуги, которые используются в программе визуализации дуги в динамическом режиме.
На рис. 3 показан вид сверху для режима, в котором токи по фазам выбраны одинаковыми, и диаграммы мгновенных значений э.д.сил образуют симметричную систему. Направление вращения поля ю0 - по часовой стрелке, а электродинамические силы вращаются навстречу полю, что предотвращает раскручивание ниппельных соединений электродов при работе на жидкой ванне.
э
р
тт
МВАр Go
мОм /4 % 6 yCPo
6' ■ !?- ■ loo-
4- ■ ю- \A
5- - 8 ■
2 ■ - 6- se-
і • 4-
/0 20 30 40 ^ к/7
Рис. 2. Электрические характеристики ДСП (АКП), ПТ 32 МВА, ступень 9, ист=306 В, короткая сеть: гкс=1,033 мОм, xкс=3,146 мОм
Рис. 3. Электродинамические силы (вид сверху) при !=20 кА, P0=8,2 МВт, Q0=5 МВАр
На рис. 4 изображены те же силы в изометрии. Показаны направление трёхфазного магнитного поля ю0, зеркало ванны и длина дуги, численно равная ид (117 мм). Векторы мгновенных значений э.д.сил для каждой фазы лежат в одной плоскости, силы при заданных токах направлены вверх и в сторону стен печи. Под действием этих сил происходит выдувание ореола дуги вверх, в результате чего может гореть футеровка свода.
При увеличении тока до 35 кА (рис. 5) э.д.сил меняют направление, но угол у лишь незначительно больше 90°, поэтому силы направлены почти горизонтально в сторону стен печи, что может
способствовать прогоранию футеровки в области шлакового пояса.
Выводы
Использование вычислительной техники позволяет не только выполнять расчёты электродинамических сил, действующих на дуги в трёхфазной ДСП, но и получать изображение этих сил как в статике, так и на действующей печи в реальном масштабе времени. Исследование э.д.сил в разные периоды плавления дает возможность анализировать механизм их воздействия на электрическую дугу и параметры технологического процесса. В результате исследований установлено.
Рис. 4. Электродинамические силы в изометрии при !=20 кА, P0=8,2 МВт, Q0=5 МВАр
^Шеепетд
0Я«* р-п | Г ^ Ьуіс | Г Заданные значения
" 1 Э5.00 АІкА ПС
-1 35.00 ий
131 Э5.00 0]кА Р0
ти 700 АІ ▼1 сю
Рис. 5. Электродинамические силы при !=35 кА, P0=13,8 МВт, Q0=12,4 МВАр
1. Вращающийся вектор э.д.сил для каждой фазы направлен вверх или вниз (в зависимости от величины тока) и в сторону стен печи.
2. Вектор э.д.сил вращается навстречу магнитному полю дуговой печи.
3. Э.д.силы резко возрастают при наличии в печи холодной шихты (ц>>1) и выдувают вверх и в сторону стен печи не только ореол, но и столб
дуги, что способствует интенсивному плавлению металла.
4. При достижении критического значения тока, соответствующего началу экстремального участка зависимости активной мощности дуги от тока, сечение дуги Бд лавинообразно растёт, растекаясь по всему торцу электрода, а результирующая электродинамическая сила меняет направление (сверху вниз).
5. Графическое изображение э.д.сил, выведенное на пульт сталевара, позволяет выбрать значение тока, оптимального для каждого периода плавки.
Литература
1. Определение электродинамических усилий, действующих на электрические дуги и электроды в ДСП / Н.В. Евсеева, И.В. Гончаренко, В.И. Смирнов и др. // Инжиниринг, инновации, инвестиции: сб. науч. тр. - Челябинск: Изд-во ЧНЦ РАЕН, 2006. - Вып. 8. - С. 83-98.
2. Аксиальные электродинамические силы, действующие на дугу в электродуговой печи / Н.В. Евсеева, И.В. Гончаренко, В.И. Смирнов и др. // Инжиниринг, инновации, инвестиции: сб. науч. тр. - Челябинск: Изд-во ЧНЦ РАЕН, 2006. -Вып. 8. - С. 76-82.
3. Мгновенные электродинамические силы, действующие на дугу в трёхфазной дуговой сталеплавильной печи / Н.В. Евсеева, И.В. Гончаренко, В. И. Смирнов и др. // Современные технологии и бизнес: сб. науч. тр. - Челябинск: Изд-во ЧНЦ РАЕН, 2006. - Вып. 1. - С. 81-95.
4. Анализ характера вращения электродинамических сил /Н.В. Евсеева, И.В. Гончаренко, В.И. Смирнов и др. // Наука и производство: сб. науч. тр. - Челябинск: Изд-во ЧНЦ РАЕН, 2007. - С. 172-180.
5. Сисоян, Г.А. Электрическая дуга в электрической печи /Г.А. Сисоян. - М.: Металлургия, 1974.
6. Арцимович, Л. Плазменные ускорители / Л. Арцимович. - М., 1973. - С. 158.
7. Холявский, Г.Б. Расчёт электродинамических усилий в электрических аппаратах / Г.Б. Холявский. - М.; Л.: Госэнергоиздат, 1962.
8. Анализ электрических параметров дуги переменного тока в дуговой сталеплавильной печи // Н.В. Евсеева, И.В. Гончаренко, В.И. Смирнов и др. // Интеллектика, логистика, системология: сб. науч. тр. - Челябинск: Изд-во ЧНЦ РАЕН, 2005. -Вып. 15. - С. 44-51.
9. Спелицин, Р.И. Влияние электрических режимов работы высокомощных ДСП на заглубление дуги в жидкую ванну /Р.И. Спелицин, Н.А. Пирогов, В.Д. Смоляренко // Производство электростали: темат. сб. - М.: Металлургия, 1977. - № 6.
- С. 46-50.
Поступила в редакцию 02.03.2011 г.
Евсеева Нина Васильевна - инженер, ОАО «Челябинский металлургический комбинат». Контактный телефон: (351) 725-41-71.
Yevseyeva Nina Vasilievna - an engineer of Join Stock Company “Chelyabinsk Metallurgical Plant”. Tel.: (351)725-41-71.
Лазуко Леонид Анатольевич - начальник участка, ОАО «Челябинский металлургический комбинат». Контактный телефон: (351) 725-41-71.
Lazuko Leonid Anatolievich is the head of the workshop of Joint Stock Company “Chelyabinsk Metallurgical Plant”. Tel.: (351) 725-41-71
Черкасова Юлия Борисовна - ведущий программист, ОАО «Челябинский металлургический комбинат». Контактный телефон: (351) 725-36-89.
Cherkasova Yulia Borisovna - a senior programmer of Joint Stock Company “Chelyabinsk Metallurgical Plant”. Tel.: (351) 725-36-89.
Хасанов Сабиржан Урунбасарович - ведущий математик, ОАО «Челябинский металлургический комбинат». Контактный телефон: (351) 725-36-89.
Khasanov Sabirzhan Urunbasarovich - a senior mathematician of Joint Stock Company “Chelyabinsk Metallurgical Plant”. Tel.: (351) 725-36-89.
