CC BY
42
УДК 658.512.011; 519.711.3; 669.18
ОПТИМИЗАЦИЯ КОНСТРУКЦИИ ЭЛЕКТРОДУГОВОЙ ПЕЧИ ЭНЕРГОМЕТАЛЛУРГИЧЕСКОГО КОМПЛЕКСА МЕТОДОМ КОМПЬЮТЕРНОГО ИНЖЕНЕРНОГО АНАЛИЗА
А.А. Арсеньева
Представлены результаты математического моделирования процесса получения жидкого железа прямого восстановления на базе уравнения термодинамического состояния веществ. Описаны результаты исследования влияния конструктивных элементов печи на движение расплава и распределение концентраций и температуры по объёму расплава.
Ключевые слова: физико-математическое моделирование, электродуговая сталеплавильная печь, железо прямого восстановления, газификация угля.
Физическая и математическая модель исследования. Объектом оптимизации является электродуговая печь для плавки стали непрерывного действия (ДСП). Исходными материалами процесса являются кислород и порошки руды и угля. В печи протекает процесс восстановление оксидов железа растворенным в железе углеродом. Так же осуществляется процесс получения газообразного топлива из находящегося в растворе с железом угля, реагирующего с расплавленными оксидами железа, поступающими в ДСП из железной руды, и жидкими оксидами железа, образующимися при подаче газообразного кислорода внутрь объема железоуглеродного расплава рис.1.
Дуговая сталеплавильная печь имеет корпус, в котором наводится металлическая и шлаковая ванны, три графитовых электрода, три донные фурмы для подачи шихты (железорудного концентрата) и три донные фурмы для подачи порошкообразного угля. Порошкообразный уголь подается в большом количестве для осуществления реакции восстановления железа из поступающего в ДСП железорудного концентрата. Три кислородные фурмы ориентированы для подачи кислорода в объем железоуглеродистого расплава. Электрические дуги используются для начального создания объёма расплава и для управления процессом. Железная руда и уголь имеют меньшую плотность, чем расплав и перемещаются вверх к поверхности расплава, создавая конвективное течение расплава, которое усиливается кислородными струями. Это течение обеспечивает распределение оксидов железа и углерода по объёму и возможность реакции между ними, при которой выделяется значительное количество монооксида углерода. По мере подачи железной руды и угольного порошка уровни металлической и шлаковой ванн поднимаются. В ходе гибридного процесса положение графитовых электродов регулируют, поддерживая постоянную длину (напряжение) электрических дуг.
142
Рис. 1. Устройство дуговой плавильной печи энергометаллургического комплекса: 1 - корпус печи; 2 -электроды; 3 - фурмы для подачи железной руды; 4 - фурмы для подачи угольного порошка, 5 - летка слива железа прямого восстановления, 6 - летка слива шлака,
7 - кислородная фурма
Задачей является определение оптимального расположения конструктивных элементов печи, при котором обеспечивается равномерное распределение температуры и концентрации реагентов по объёму печи. Решение этой задачи возможно только методами компьютерного инженерного
анализа, основанными на физико-математическом моделировании процессов в изучаемом объекте. Критериями оценки являются термодинамическое состояние материалов в печи и химический состав расплава. Так как процесс протекает во времени и пространстве, то термодинамическое состояние непрерывно изменяется вследствие выделения теплоты электрической дугой и химическими реакциями. Химический состав также непостоянен вследствие протекающих реакций и подачи элементов во время плавки.
Физико-математическое моделирование процесса заключается в решении системы дифференциальных уравнений энергии и массопереноса.
Так как дуговая печь близка по форме к телу вращения, то при моделировании будет использоваться цилиндрическая система координат: радиального расстояния от оси симметрии г , угла поворота ф относительно плоскости оси электрода и расстояния от пода печи г.
Пространство моделирование условно разделено на области, рис.1: Е - область угольных электродов; В - область электрических дуг; Я - область расплавленного шлака; М - область расплава металла; Р - область футеровки печи; О - область газовой среды. Каждая из областей в зависимости от свойств веществ и процессов, происходящих в ней, описывается специальной системой уравнений.
Во всех указанных областях ДСП протекает нестационарный термодинамический процесс, который описывается изменением энтальпии Н (?) множества точек пространства во времени ?. Нестационарное линейное уравнение теплопроводности в цилиндрической системе координат г, ф, г имеет вид [1]
Г ЗтЛ
1г— дг
ЭН _ 1 _д_
д? г дг дН
+уг ~дг+% - У*,
+
1 д
1——
д
+ —
дН дН
г2 дф ^ дф) дг \ дг ) дг у дф (1)
где Т - температура точек пространства; 1 - коэффициент теплопроводности среды, зависящий от координат расположения точки в пространстве, типа вещества и температуры в этой точке; уг , Уф, Уг - скорости движения вещества в направлении соответствующих координат, д^, д* - удельные значения мощности выделения и поглощения теплоты в данной точке пространства.
В жидком металле действует давление, возникающее вследствие подачи через донные фурмы Р^Оз, С и выделения СО, а также меньшей плотности данных веществ по сравнению с расплавом. Кроме того, струи кислорода создают давление на поверхности металлической ванны. Суммарное давление от перечисленных факторов вызывает перемещение расплава.
Течение жидкости описывается уравнением Навье - Стокса, которое в цилиндрической системе координат имеет вид [1]
Р
Э*г _ Эр Эх Эг
'э v
Э 2*
Э V
Э*
Р
ф _
Р
Эх
Э*
Эр
Эф 1
г Эр Эх Эг 1
V
Эг 2 г 2Эф2 Эг 2 ?
^ Э 2*ф , Э2*Ф Э 2* _1_ Л ф
Эг 2 г 2Эф2 Эг 2 У
э ч Э Ч э ЧЛ
1 Эг 2 1 г2Эф2 ' Эг2 ,
(2)
где р - плотность расплава; р - давление в данной точке пространства; ^ - динамическая вязкость.
Распределение давления в расплаве определяется при решении уравнения неразрывности
х
С í
р _ Е
Э*ф Э*г Э*„
+ ■
гЭф Эг Эг
(X +
(3)
где Е - модуль упругости; р^ - гравитационное давление.
Концентрация элементов в расплаве изменяется при плавлении компонент шихты и перехода элементов металлошихты в расплав, химического взаимодействия элементов в расплаве и перемещения продуктов взаимодействия потоком жидкости. Изменение концентрации С, 1 -го эле-
мента расплава описывается уравнением переноса
^ _ ^ Эх 1
Э
гЭг
эс
Эг
Э
г 2Эф
эс
эф
э (эс +— 1
Эг
V
Эг
у
+
С
Эг
+ V,,
ЭС, ^ ЭС
(4)
+ V
ф гЭф г Эг
где Б, - коэффициент диффузии 1-го элемента в жидком железе; *г, , *ф
- скорости движения расплава, определяемые из решения уравнения Навье
- Стокса.
Результаты моделирования представлены текущими распределениями по объёму расплава: интенсивности объёмного тепловыделения; энтальпии и температуры; плотности расплава; гравитационного давления; скоростей движения расплава в направлении каждой из координат; концентраций каждого из веществ.
Различие давлений создаёт потоки в расплаве. Кроме того, в следствии подачи струй кислорода осуществляется вращательное движение
0
г
расплава в направлении ортогональных координат в плоскости поперечного сечения ух + Уу . Подача порошка углерода и оксида железа вызывает
движение расплава в виде трёх вертикальных вихрей Уг . Движение расплава вызывает перемещение углерода и оксида железа в объёме расплава и определяет распределение их концентрации СсСрео.
Электрические дуги и окисление углерода нагревают приповерхностную область расплава, соответственно температура поверхностных слоёв выше. В придонные слои поступают холодные углерод и оксид железа, реакция взаимодействия между которыми поглощает тепло, что снижает температуру расплава.
Теплоперенос осуществляется исключительно потоками расплава, которые выравнивают температуру по объёму.
Влияние конструктивных характеристик печи на процесс плавки.
Было проведено исследование влияния расположения донных фурм для подачи угля и железной руды и расположения электродов на движение расплава и распределение концентраций и температуры по объёму расплава.
Так как реакция восстановления железа протекает в центральной зоне массы расплава с большим поглощением тепла, а источники теплоты - реакция окисления углерода потоками кислорода и электрические дуги действуют на поверхности расплава, температура в печи распределена очень неравномерно. Расположение фурм оказывает значительное влияние на неравномерность распределения температуры, рис. 2.
При малом диаметре центров фурм, рис. 2, а, в нижней центральной зоне реакция поглощения теплоты, рис. 2, в, понижает температуру расплава до недопустимого уровня (до температуры кристаллизации) несмотря на высокие значения температуры на поверхности расплава. По мере удаления центров фурм от центра, распределение температуры становится более равномерным. Но уменьшение неравномерности температур существенно при увеличении диаметров центров фурм до значения, равного 0,71 внутреннего диаметра печи, рис. 2, в. Более равномерного распределения температуры удалось добиться при уменьшении диаметра центров электрических дуг, рис. 3, б, т. е. за счёт более равномерного распределения тепловыделения на поверхности расплава.
Но наиболее равномерное распределение температуры в печи обеспечивается при сдвиге осей электрических дуг и зон окисления углерода струями кислорода относительно центров донных фурм на 30о, рис. 3, в. Это объясняется тем, что теплоперенос в печи осуществляется преимущественно движением расплава.
На рис. 4 показано распределение температур и давлений на поверхности расплава, и концентрации моноокиси углерода, которая выделя-
146
ется при восстановлении железа и которая соответствует интенсивности поглощения теплоты.
а) 6) в)
Рис. 2. Распределение температур в разных сечениях ДСП при диаметре центров фурм: а - 200 см; б - 400 см; в - 420 см
а) 6) в)
Рис. 3. Распределение температур в разных сечениях ДСП при диаметре центров электрических дуг: а - 420 см; б - 200 см; в - 200 см и повороте электродов и кислородных фурм на 30°
Рис. 4. Распределение температур (а), давления (б) на поверхности расплава и концентрации монооксида углерода (в) на глубине 1/3 толщины расплава при диаметре центров электродов (1) 200 см, диаметре центров кислородных фурм (1) и фурм для подачи угольного порошка (3) и руды (4) 420 см
Под действием подачи железной руды, углерода через донные фурмы 3, 4 (рис. 4, в) и струй кислорода, касательно направленных на поверхность расплава 2 (рис. 4, а, б), возникает давление, которое на поверхности имеет ярко выраженные зоны растяжения А (рис. 4, б) и сжатия Б (рис. 4, б). Это давление создаёт сложное вихревое движение расплава, направленное вверх от донных фурм и вниз в зоне воздействия струй кислорода. При сдвиге кислородных фурм на 30о зона тепловыделения реакции окисления углерода кислородом оказывается в зоне нисходящего потока расплава. Этот поток уносит расплав, нагретый этой реакцией и электрическими дугами, в зону эндотермической реакции восстановления железа.
Разработанная компьютерная модель гибридного процесса получения жидкого железа прямого восстановления и газификации углерода в дуговой электропечи энергометаллургического комплекса позволяет решить задачу оптимизации конструкции элементов ДСП и технологии ведения данного процесса. Полученные результаты численного эксперимента показывают, что наиболее рациональной является конструкция дуговой плавильной печи, при которой центры нижних фурм для подачи измельчённых руды и угля расположены от оси печи на расстоянии 0,71 её внутреннего диаметра. Центры зон воздействия кислородных струй и электрических дуг сдвинуты относительно центров донных фурм на угол 30о. Центры электрических дуг желательно расположить ближе к оси печи, на расстоянии, равном 0,35 - 0,5 внутреннего радиуса печи. Выводы
1. Разработана физико-математическая модель функционирования энергометаллургического процесса получения железа прямого восстановления, в котором измельчённые железная руда и уголь подаются в расплав железа через шесть равнорасположенных донных фурм, а поверхность расплав нагревается реакцией окисления углерода струями кислорода и
электрическими дугами.
2. Разработаны метод и алгоритм численного решения для выполнения численного исследоваля влияния расположения фурм для подачи руды, угля, кислорода и электродов на распределения по объёму расплава скоростей его движения, концентраций окиси железа, углерода и моноокиси углерода, интенсивности выделения и поглощения теплоты и температуры расплава.
3. Показано, что равномерное распределение температуры в расплаве железа, необходимое для нормального протекания процесса, обеспечивается при расположении зон выделения теплоты на поверхности расплава над зонами нисходящих вихревых потоков, возникающих в расплаве при подаче руды и угля и струй кислорода.
4. Установлено, что наиболее рациональной является конструкция дуговой плавильной печи энергометаллургического комплекса, при которой центры нижних фурмы для подачи измельчённых руды и угля расположены от оси печи на расстоянии 0,71 её внутреннего диаметра. Центры зон воздействия кислородных струй и электрических дуг сдвинуты относительно центров донных фурм на угол 30о. Центры электрических дуг желательно расположить ближе к оси печи, на расстоянии, равном 0,35 - 0,5 внутреннего радиуса печи.
Работа представлена на 3-й Международной Интернет - конференции по металлургии и металлообработке, проведенной в ТулГУ 1 мая - 30 июня 2014 г.
Список литературы
1. Тихонов А. Н., Самарский А. А. Уравнения математической физики. М.: Наука, 1972. 735 с.
2. Тихонов, А. Н., Калько В. Д., Гласко В. Б. Математическое моделирование технологических процессов и метод обратных задач в машиностроении. М.: Машиностроение, 1990. 264 с.
3. Рябов А. В., Чуманов И. В., Шишимиров М. В. Современные способы выплавки стали в дуговых печах. М.: Теплотехник, 2007. 192 с.
4. Richard MacRosty Modelling, Optimization and Control of an Electric Arc Furnace. Hamilton: McMaster University, 2005. 160 p.
Арсеньева Алина Алексеевна, асп., [email protected], Россия, Тула, Тульский государственный университет
DESIGN OPTIMIZATION OF ELECTRIC ARC FURNACE OF A POWER-USING COMPUTER ENGINEERING ANALYSIS
A.A. Arseneva 149
The results of mathematical modeling of the process of producing liquid DRI based on the thermodynamic state equation substances. The results of studies of the effect of structural elements on the motion of the melt furnace and the concentration distribution and temperature of the melt volume.
Key words: physical-mathematical modeling, electric arc furnace, direct reduced iron, coal gasification.
Arseneva Alina Alekseevna, postgraduate, Silabykv@yandex. ru, Russia, Tula, Tula State University
УДК 621.833
ИНСТРУМЕНТАЛЬНЫЕ НАЛАДКИ ПРИ ИЗГОТОВЛЕНИИ
ЧЕРВЯЧНЫХ КОЛЕС
В.Н Скрябин, И.И. Феофилова, С. Л. Рахметов
Рассмотрены сборные червячные фрезы, инструментальные и станочные наладки, необходимые при изготовлении червячных колес ортогональных передач. Показано, что локализация пятна контакта в передаче способствует компенсации погрешностей изготовления и сборки, а также повышению технологичности изготовления червячной фрезы.
Ключевые слова: фреза, передача, профиль, колесо.
При производстве червячных ортогональных передач строгая сопряженность активных поверхностей зубьев является скорее недостатком, чем достоинством. С целью компенсации неизбежных погрешностей изготовления и сборки, а также силовых и температурных деформаций зацепление следует проектировать как приближенное, локализуя контакт в окрестности некоторой расчетной точки. Приближенность вводится искусственно путем модификации одной сопрягаемой поверхности относительно другой так, чтобы при любых допустимых погрешностях пятно контакта не выходило ни на одну из кромок зуба, а точность передачи по нормам плавности оставалась в пределах допуска.
Вопросам приближенного червячного зацепления посвящены работы Зак П.С., Лагутина С.А., Сандлера А.И., Верховского А.В. и других. Особое внимание они уделяли восстановлению крупномодульных червячных передач путем изготовления обоих звеньев передачи, т. е. червяка и колеса с учетом возможностей инструментального производства и наличия червячных фрез с параметрами близкими к параметрам червяка вышедшей из строя червячной передачи.
В производственных условиях эксплуатации сложных технологиче-
