CC BY
28
В1СНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХН1ЧНОГО УН1ВЕРСИТЕТУ 2003 р. Вип. №13
УДК 669.184
Сущенко А.В.*
ОПРЕДЕЛЕНИЕ КРИТИЧЕСКОЙ КОНЦЕНТРАЦИИ УГЛЕРОДА В КИСЛОРОДНОМ КОНВЕРТЕРЕ С ВЕРХНЕЙ ПРОДУВКОЙ
С использованием масспектрометрической системы «ФТИАН-3» выполнены исследования динамики процесса обезуглероживания расплава в 350 т кислородных конвертерах верхнего дутья, получены полуэмпирические и регрессионные зависимости для определения критических концентраций углерода и относительного расхода расплава, циркулирующего через реакционную зону.
Прогнозирование наступления момента времени продувки конвертерной плавки т» (или достижения концентрации углерода в расплаве ванны [С]*), соответствующих началу снижения скорости обезуглероживания ис в конце периода преимущественного окисления углерода имеет большое значение для теории и практики кислородно-конвертерного производства стали [1-4 и др.].
Снижение величины ис при [С] < [С]* связывают с тем, что потоки металла в конвертерной ванне начинают подводить к реакционной зоне (РЗ) массу углерода меньшую, чем та, которая может быть окислена в соответствии с интенсивностью подвода окислителя 102 ,
м3/(т-мин), [1-10]. Неиспользованный на окисление углерода в РЗ кислород (в виде растворенного в металле и оксидов Бе) переносится в ванну и процесс обезуглероживания начинает дополнительно развиваться в ее объеме. При этом вследствие «улучшения перемешивания» [1] может в начале этого этапа несколько увеличиться, а затем снижается пропорционально [С].
В [5, 7] эту характерную концентрацию углерода ([С]») называют переходной - [С]пер, в [6, 8-10] - критической - [("],у), а в [1, 3, 11] - и критической и переходной одновременно. В ряде работ [4, 5, 7, 12] понятия переходной и критической концентраций углерода разделяют, принимая за [С]пер концентрацию углерода, соответствующую началу интенсификации процесса в ванне за пределами РЗ, а за \С\кр - концентрацяию углерода, соответствующую началу последовательного снижения ис в агрегате.
По вопросу о величине \%С\* также нет единого мнения. По различным данным для конвертеров верхнего дутья она равна: 0,5-1,0 (1,5) [1,3,7,13,14], 0,2-0,7 [15], не более 0,3 [16], а для конвертеров с донной и комбинированной продувкой: 0,3-0,6 [1,13,14,17] и 0,3-0,8 [4,18,19 и др.] соответственно.
Значительный "разбег" данных, особенно для ЬБ-конвертеров, связан как с многофакторностью и быстротечностью процесса объёмного (за пределами РЗ) обезуглероживания расплава, так и с различием методик определения ис. Наиболее достоверная информация о процессе может быть получена непосредственно на промышленных конвертерах с использованием методов непрерывного контроля ис (например, по газовому анализу отходящих конвертерных газов [20 и др.]) неподверженных влиянию дискретности отбора проб, плавления лома, промежуточных повалок конвертера и т.п.
В [21,22] получены статистические уравнения, связывающие величину [С]* с 102 в основном для малотоннажных конвертеров верхнего дутья, причём в [21] - зависимость квадратичная, а в [22] - линейная.
Целью настоящей работы являлось определение величины [С]* и функциональных связей её с параметрами продувки на 350 т конвертерах верхнего дутья (МК «Азовсталь») с использованием непрерывного контроля ис по расходу и составу отходящих газов. Величину \С\кр, отвечающую характерному изменению скорости окисления углерода в кон-
* ПГТУ, канд. техн. наук., зав. отделом ПНПЛ
вертерной ванне, определяли по показаниям масспектрометрической системы "ФТИАН-3" [20]. При этом пределы допустимой абсолютной погрешности измерения содержания углерода в расплаве для 90 % плавок составляли не более 0,04 , 0,06 и 0,1 % в области содержаний углерода 0,2-0,4 , 0,4-0,6 и 0,6-0,8 % соответственно. Погрешность измерения объёмных концентраций СО, С02 , 02, Л'2. Щ в конвертерных газах (дискретность отбора проб - 2,4 с) составляла не более 0,2 % [23]. Для анализа отбирали плавки, на которых результаты расчёта текущей концентрации углерода в конвертерной ванне (в точке, соответствующей [С]к;,) по статическому (на основе уравнения баланса углерода) и динамическому (на основе регрессионного уравнения кинетической кривой заключительного этапа процесса обезуглероживания) методам расчёта совпадали с абсолютным расхождением не более 0,1 %.
В качестве критерия, характеризующего влияние конструктивных и режимных параметров дутьевых устройств на процесс взаимодействия струй с ванной была принята относительная (по отношению к глубине ванны в спокойном состоянии Нв) глубина проникновения дутьевых струй в расплав Ь = Нсгпр I Нв( 1-4 и др.). Величину Нстр определяли с использованием зависимостей [24]. Исследования были проведены в 1988 г. с использованием кислородных фурм с пятью продувочными соплами (диаметр минимального сечения £/„„„ = 42 мм; число Маха Мс= 2,293; угол наклона к оси фурмы а = 15°). Учитывая относительное постоянство величины 1а^ для конкретного агрегата рабочий расход кислорода на продувку У(Х был ограничен диапазоном 1000-И200 м7мин.
Результаты обработки данных промышленных плавок с использованием стандартных методов математической статистики представлены на рис. 1. При этом весь массив контролируемых плавок был разбит на две части: при = 1000-^1100 м /мин и 1ги = 1100н-1200 м7мин.
Полученные функциональные зависимости
[%С]кр = 0,220 / Ь при 1<ь = 2,9+3,1 м7(т-мин); [%С]кр = 0,184 / Ъ при Зщ = 3,1*3,4 м7(т-мин); [%С\Щ} = 0,200 / Ъ при /01 = 2.9-3.4 м7(т-мин);
(1)
(2) (3)
и
описывают экспериментальные точки с корреляционными отношениями травными 0,781 (1), 0,671 (2) и 0,658 (3) соответственно.
Изменение величины г в выражениях (1)-(3) свидетельствует о дополнительной существенной связи [%С\кр с интенсивностью продувки.
При дальнейших исследованиях процесса обезуглероживания расплава в конвертере, связанных с оптимизацией параметров кислородных фурм в ККЦ МК "Азовсталь" (1994/95 гг.) было установлено, что кинетическая кривая процесса имеет более сложный
характер. Фактически величина [%С].кр , определённая с использованием системы "ФТИАН-3" в (1)-(3), отражает ход суммарной скорости окисления углерода во всей ванне - цх, а не только в районе РЗ. В кислородных конвертерах верхнего дутья, в отличие от агрегатов с донной продувкой, на
характер изменения величины и^, после достижения критической концентрации углерода для условий его
А ' / 0,220 1 Ь
. /. 1 1ч, = 0,200 ' / * Ь
"■--. 1 А гч, » ... / [%с / ® 7 .' ]кр- 0,184 / Ъ •
"а- А А • •/ — .'"*.-.- / А
сЗ КС
0
а. (о К 1-
т
й
Н
1
о
а
д
О 0.4 0.45 0.5 0.55
Я,
Относительная глубина проникновения дутьевых струй в расплав, Ь
Рис. 1 - Зависимость критической концентрации углерода в конвертерной ванне от относительной глубины проникновения дутьевых струй в расплав: ▲ - г = 0,781; Га е [1000; 1100] м3/мин;
• - г = 0,671; е (1100;1200] м3/мин;
Ж, • - г = 0,658; У0 е [1000; 1200] м3/мин
окисления в РЗ [25] - [С]^, оказывает существенное влияние состояние шлака (шлако-газо-металлической эмульсии - ШГМЭ), а также температура и окисленность ванны. В зависимости от указанных факторов в крайних случаях и^ либо резко и непрерывно уменьшается - при "свёрнутом" шлаке, либо резко увеличивается - при развитой ШГМЭ и высокой температуре ванны. В промежуточных случаях имеет место непродолжительное снижение и^ с последую-
щим "всплеском" и непрерывным уменьшением. На различный характер изменения и^ в области \CYkI указывалось в [15,25,26 и др.].
На основе данных газового анализа и наблюдений за ходом плавок было установлено, что величина [С]^ в рассматриваемых условиях (при содержании углерода в металле после продувки [%С\К = 0,04-0,06) достигается при относительном времени продувки / т„р= 0,70-0,80 (в среднем ~ 0,75). В отдельных случаях величина ? 3 уменьшалась
Т,.„ = т..
до 0,64, что было связано с аномальным ходом плавок - малой величиной Нстр при больших значениях и т.п. Такие плавки из анализа исключались. На рис.2, в качестве примера, представлены типичные кривые изменения и^, по ходу продувки плавок с различными 1а , где
стрелками отмечены моменты достижения [С]^3 для соответствующих плавок. На рис.3 приве-
J кр '
дены данные по динамике содержания углерода в расплаве ванны на промышленных плавках с промежуточными повалками, проведенными сотрудниками конвертерной лаборатории ЦЛМК МК "Азовсталь" (№№ 1-6) и с участием автора (№№ 7 -9). В соответствии с этими данными
величина | %С' | Ркр составляла в среднем 0.7-:-1.2 (при т''^ рошо согласуются с данными авторов [26-29 и др.].
= 80
О
и
5 I 60
5 и
га у
а в к ^
« Я
О Рч ,(1
ш и 4и § а
3 « ¡■ч ра
« х
х § 20 «
% г? I
¡к и и к
V
1г\ 1 *** Т / // -ч 1 1 V
1' - Уг-Х \ V, к В Ал л*
Тя
о
20
0,75). Полученные результаты хо-
Для более точного определения величины | С |:использо-
40
60
80
100
Относительное время продувки плавки, т = т / т %
Рис.2 - Характерные кривые динамики скорости обезуглероживания расплава в 350-т конвертере МК "Азовсталь" по ходу продувки плавок (непрерывный отбор газовых проб - на 56-ой отметке за ОКГ):
№
1
2
3
4
№ плавки
Дата
1528 1576 592 679
19.04.94г. 22.04.94г. 13.02.94г. 18.02.94г.
Расход 02 в области
/%с/:г ~[%С]кр
м/мин
1300 1080-1130 1300-1350 940-1000
[%С]К
0,04 0,04 0,04 0,04
Щ> 5
мин-с
14-40
15-30 13-43 15-43
вали систему "ФТИ-АН-3" с соответствующей корректировкой настроечных коэффициентов её расчётных моделей при проведении экспериментальных исследований. Плавки проводились (1994 г.) в основном с использованием кислородной фурмы с пятью основными (с/,,;;;; = 47 ММ; Мс= 2,215; а =15°) и одним центральным = 30 мм; Мс = 1,5) соплами. Для анализа отбирали только плавки с ярко выраженным
изменением и^, в об-
1 рз
ласти [С]кр .
Попытка установить функциональную
сЗ
ч
о
Сы I-
а и
(н 0\ «
и
н я
о
я
X
о
5,0 4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0,5 0,0
*
0 10 20 30 40 50 60 70 80 90 100
Относительное время продувки, г Рис.3 - Изменение концентрации углерода в расплаве конвертерной ванны по ходу продувки плавок: • - 1; о - 2; А - 3; А - 4:, ■ - 5; □ - 6: ♦ - 7; 0 - 8; * -9;
мч, и, |%Л7|Ч \%Мп\ч Ул,
т °С % м / (т-мин)
1 295 1320 0,74 0,66 23,4 3,2
2 290 1320 0,72 0,32 23,7 3,3
3 270 1315 0,97 0,51 25,0 3,4
4 290 1290 0,78 0,47 25,6 3,3
5 290 - 0,69 0,32 23,7 3,1
6 300 1350 0,41 0,47 23,0 3,7
7 290 1300 0,96 0,64 25,6 3,6
8 295 1270 0,80 0,42 20,3 3,4
9 295 1300 0,85 0,40 21,3 3,4
Л/.,. |.%Л7 [,., [%М?]Ч - параметры чугуна; у/, - доля лома в металлошихте; 1а^ - средняя (за плавку)
удельная интенсивность продувки.
связь между параметрами [С]^ и Ъ по
типу выражений (1)-(3) оказались неудачными вследствие низкого коэффициента парной корреляции (0,35). Приняв за основу зависимость
для определения через относи-
тельный расход расплава, циркулирующего через РЗ [9,25], и представив её в виде
[%С]Ц Щ -т2« .(УЪ),
(4)
где т - массовый расход дутья (кислорода), кг/с; к1 - эмпирический коэффициент,
с учетом, что Ъ ~ т1 3 [30 и
о2
др.], получили
[%С]1
(5)
В соответствии с результатами статистической обработки экспериментальных данных (рис.4) зависимость (5) была представлена в виде
[%С£ = 0,090-0,634-з/7
(приб = 350т), г = 0,721.
(6)
1.2.
св
КС
0 , „
а, 1.1 а
^ 1Д
1 &
0,9
га г—I
5 0,8, Я 3
О 0,7 0,6
•
Ц С]^3= 0,42 + 0,00046г = 0,725
■ \
V -
* 1
•
9 • 1
Щ г ч
X \
• \ ,090*2 г - 0 ,721
\
>10 ,028• К г о ,725
900
1000
1100
'1200
1300 3/
1400
Расход кислорода, , мУмин Рис.4 - Зависимость величины критической концен-
трации углерода от расхода кислорода на
1 кр
продувку Уа (в области /%("/,'' ')
кр
Сравнение результатов расчёта величины | %(" | '.у по выражению (6) с
известными зависимостями (рис.5) показало их хорошее качественное соответствие данным [31]; при этом значения величины по (63) оказались
на ~15 % меньше, что может быть связано с различием методик её определения.
Статистическая обработка полученных экспериментальных данных в виде степенной функции
~ к3 ■ \ позволила установить
показатель степени к4 = 0,52 и получить выражение, несколько лучше описывающее искомую связь
[%С]РК3р = 0,028-й0,524-^, г = 0,725.
(7)
и
cS
ч о р,
о «
ÎH
Ы s
и &
я
(а Я К о
M
1,4
и
1,0 0,8 0,6 0,4 ОД
t°/»CLf = 1,319 0,6 »Э/10а [Черн Ьгевич А.Г.]
1,13 + 0,179*1
fi V°cg » 0,634 зД^
—' ; H [Кричевцов S.A. и др.]
; £ H0,lJ5-IO3 |I огушевский и др.]
2 2,5 3 3,3 4 4,5 5
Интенсивность кислородной продувки, /0 , м7(т-мин) Рис. 5 - Зависимости критических концентраций углерода /%( '/'',; и [%С]кр от интенсивности подачи кисло-
рода в конвертер
вертикальными пунктирными
линиями ограничена область проведения исследований на 350 т конвертере верхнего дутья
лить относительный расход расплава, циркулирующего верхней продувкой {О = 350 т)
Для уточнения значений к4 требуются дополнительные исследования в более широком диапазоне изменения величины |i .
ui
Отметим, что отличие показателя степени в (5) от 1/3 может быть связано с тем, что в конвертерах верхнего дутья при изменении параметра Fq "автоматически" осуществляется коррекция высоты фурмы 11 ф (заглубления струй - Ъ). Если допустить, что при Fa = var b ~ const, то в соответствии с (4) показатель степени при
К0 составит 2/3. По-видимому
величина этого показателя на практике находится в диапазоне 1/3... 2/3 и близка к 1/2.
Используя полученные зависимости (6) и (7), а также положения работы [25], можно опреде-через РЗ кислородного конвертера с
qP3 = 833,3'3 « 118,2-10]13, кг/кг 02, qP3 = 2679- V0]'2 « 143,1 • /0У2, кг/кг 02 .
(8) (9)
щ -
Ш g
о
а о,е -о
4
Ь чо
сч '
5 а,
3 "
^
s
о
а £
о «
0,4 ■
°С]кр = -0,13 +0,00051'Vq2 ; г = 0,782
N
X
• Ч •
■s
• ( •
•
•
и**** m
•
•
•
Расход кислорода, I, м/мин
Рис.6 - Зависимость величины критической концентрации углерода \С\кр от расхода кислорода на продувку (в области
При определении величины критической концентрации углерода для условий окисления его во всей ванне и соответствующей последовательному снижению и^ - \%С\,, были исключены случаи, когда « |%( ' . т.к. та-
кие плавки являются, как правило, технологически неоптимальными (имеют место значительные теплопотери вследствие раннего "свёртывания" шлака, более высокая окисленность ванны в конце продувки и т.п.). Величина \%С]кр= 0,28-0,64 (* 0,30-0,60) достигалась при г = гкр/тпр» 0,80-0,90 (в
среднем ~ 0,84). Полученная статистическая зависимость для определения [%С\кр (рис.6) имеет вид
[%С\кР = -0,13 + 0,0005 1-Fq
-0.13
0,179- г = 0,782,
(10)
она оказалась ближе к данным работы [21] (см. рис. 5).
В области \%С\Р процесс одновременно протекает в различных зонах (РЗ, ШГМЭ,
основном объёме ванны) h может лимитироваться разными звеньями. Причем моменты времени, соответствующие достижению критических концентраций углерода \С\кр и началу последовательного снижения ис - К'|,7/ в этих зонах могут существенно различаться. Так, при окислении углерода на пузырях СО в объёме ванны [%С]кр = 0,05-0,08 , а [%С]Ш = 0,20-0,40 [25 и др.]. Поэтому получение аналитической (детерминированной) зависимости [%С]кр трудно
осуществимо.
Выводы
¡.Установлено, что для 350 т кислородных конвертеров верхнего дутья (при [%С]К = 0,04-0,06 , Iq2 = 2,4-4,0 м3/(т-мин) ) величина критической концентрации углерода для
условий окисления его в РЗ - [%С]^ , равная 0,7*1,2 , соответствующая смене основного лимитирующего звена процесса (с «кислородного» на «углеродное»), достигается при относительном времени продувки т^р = тI тпр = 0,70-0,80 (в среднем ~ 0,75). При этом критическая концентрация углерода для условий окисления его во всей ванне, соответствующая последовательному снижению Vç - |%С|,;/) « 0,30-0,60 достигается при т = ткр/ тпр « 0,80-0,90 (в среднем ~ 0,84).
2. Получены полуэмпирические и регрессионные зависимости параметров [%С]^, [%С]Кр и относительного расхода расплава, циркулирующего через РЗ - qP3 (кг/кг 02) от интенсивности верхней продувки. Показано, что величина [%С]^ пропорциональна расходу кислорода (интенсивности продувки) в степени 1/3...2/3 (на практике - ближе к 1/2).
3. Задачами дальнейших исследований по рассматриваемому вопросу являются: 1) проверка и уточнение полученных зависимостей на LD-конвертерах с различной садкой и в более широком диапазоне изменения I; 2) моделирование динамики процесса обезуглероживания в
конвертере между моментами продувки плавки, соответствующими достижению [%С]^ и
[%С]кр .
Перечень ссылок
1. Баптизманский В.И. Конвертерные процессы производства стали / В.И. Баптизманский, М.Я. Меджибожский, В.Б. Охотский. - Киев-Донецк: Вища школа, 1984. - 343с.
2. Явойский В.И. Научные основы современных процессов производства стали / В.И. Явой-ский, A.B. Явойский. - М.: Металлургия, 1987. - 184с.
3. Баптизманский В.И. Проблемы технологии и теории современной конвертерной плавки / В.И. Баптизманский II Сталь. -1989,-№9,- С. 18-23.
4. Арсентьев П.П. Конвертерный процесс с комбинированным дутьем / П.П. Арсентьев, В.В. Яковлев, C.B. Комаров. -М.: Металлургия, 1991. -176с.
5. Охотский В.Б. Модель процесса рафинирования металла при продувке его кислородом сверху / В.Б. Охотский, А.Г. Чернятевич II Известия вузов. Чёрная металлургия. - 1972. -№ 10. -С.61-64.
6. Лузгин В.П. Газы в стали и качество металла / В.П. Лузгин, В.И. Явойский. - М.: Металлургия, 1983. -232с.
7. Охотский В.Б. Динамика окисления углерода при продувке металла / В.Б. Охотский II Сталеплавильное производство. Тематич. отрасл. сб. - М.: Металлургия, 1976. - №4. - С.20-26.
8. Охотский В.Б. Перемешивание сталеплавильной ванны в конвертере / В.Б. Охотский II Изв. АН СССР. Металлы. - 1986. - № 6. - С.3-8.
9. Капустин Е.А. Влияние гидродинамики на процесс обезуглероживания при продувке расплава заглубленной струёй / Е.А. Капустин, A.B. Сущенко II Известия вузов. Чёрная металлургия. - 1986. -№ 8,-С. 149-150.
10. Баптизманский В.И. О критической концентрации углерода в металле сталеплавильных
ванн / В.И. Баптизманский // Металлургия и коксохимия. - К.: Техшка, 1988. - С.3-8.
11. Наливайко А. П. Исследование процессов аэрогидродинамики и тепломассообмена в кислородном конвертере с целью совершенствования энергосберегающих методов продувки низкомарганцевых чугунов /А.П. Наливайко : Дис...канд. техн. наук: 05.16.02. - Днепродзержинск, 1990.
12. Явойский В.И. Металлургия стали / В.И. Явойский, Ю.В. Кряковский, В.И. Григорьев и др. -М.: Металлургия, 1983.- 384 с.
13. Saigusa М. Metallurgical control in top-and bottom blown oxygen steelmaking process / M. Saigusa, J. Naga, F. Sudo е.a // Ironmaking and Steelmaking. - 1980. - V.7. - № 5. - P.242-248.
14. Kamo E. Движение потока стали и металлургические реакции в кислородном конвертере донного дутья / Е. Kamo II Тэцу то хаганэ. - 1979. -Т.65, №8. - С. А137 -А140.
15. Богушевский B.C. Основы математического описания технологических процессов конвертерного производства стали / B.C. Богушевский, H.A. Рюмшин, Н.А.Сорокин. - К.: НПО "КИА", 1992. - 168с.
16. Меджибожский М.Я. Основы термодинамики и кинетики сталеплавильных процессов / М.Я. Меджибожский. - К.: Вища школа. - 1986. - 280с.
17. Арсентьев П.П. Конвертерный процесс с донным дутьём / П.П. Арсентьев, М.П. Квитко. -М.: Металлургия, 1983. - 128с.
18. Китамура С. Влияние интенсивности перемешивания ванны и скорости подачи в ванну кислорода на характер обезуглероживания высокохромистой жидкой стали / С. Китамура, К. Окохира, А. Танака // Тэцу то хаганэ. - 1986. - Т.72, №1. - С. А47 - А54.
19. Перемешивание металла при конвертировании с верхним и донным дутьём / Т. Осэки, К. Хасэагава, Е. Тамура и др. // Тэцу то хаганэ. - 1981. - Т.67. - № 4. - С.269.
20. Тучкевич В.М. Контроль и управление конвертерным производством стали / В.М. Тучкевич II Вестник АН СССР. - 1985. - № 5. С 59-61.
21. Кричевцов Е.А. О критической концентрации углерода в ванне кислородных конвертеров / Е.А. Кричевцов, Л.А. Смирнов, В.ГЛалетин //Известия вузов. Чёрная металлургия. - 1974. -№ 1. -С.14-17.
22. Богушевский B.C. АСУ ТП производства стали в конвертерах / B.C. Богушевский, H.A. Рюмшин, H.A. Сорокин. - К.: Техшка, 1991. - 180с.
23. Система "ФТИАН-3". Техническое описание и инструкция по эксплуатации. - ФТИ им.
A.Ф. Иоффе АН СССР, 1985. - 55с.
24. Чернятевич А.Г. К вопросу образования газо-шлако-металлической эмульсии в кислородном конвертере / А.Г. Чернятевич, Г.М. Соломон II Производство стали в кислородно-конвертерных и мартеновских цехах: Науч. тр. МЧМ СССР. - М.: Металлургия, 1981,- № 9. - С.12-18.
25. Капустин Е.А. Критическая концентрация углерода и анализ процесса обезуглероживания в сталеплавильных агрегатах / Е.А. Капустин, A.B. Сущенко // Известия вузов. Чёрная металлургия. 1988. - № 9. - С.40-44.
26. Чернятевич А.Г. Наблюдение через прозрачную стенку за поведением конвертерной ванны при продувке / А.Г. Чернятевич, Е.Я. Зарвин, М.И. Волович II Известия вузов. Чёрная металлургия. - 1975. - № 2. - С.37-42.
27. Weniger J. Integrated Automation Systems for Converter Steel Prants / J. Weniger, J. Müller II International Steelmaking Conference. - Desing Center Linz / Austria, Linz. - 1998. - P. 16/5-19/5.
28. Охотский В.Б. Шлакообразование и окисление углерода в сталеплавильных процессах /
B.Б. Охотский II Изв. РАН. Металлы. - 1992. - №1. - С.20-23.
29. Харлашин П.С. Мышьяк и его роль в металлургических процессах / П.С. Харлашин. - К.: Вища школа, 1993. - 303 с.
30. Марков Б.Л. Физическое моделирование в металлургии / Б.Л. Марков, A.A. Кирсанов. - М.: Металлургия. 1984. - 119с.
31. Чернятевич А.Г. Разработка энергосберегающих методов продувки конвертерной ванны на основе высокотемпературного моделирования / А.Г. Чернятевич : Автореф. дис... д-ра техн.наук:05.16.02 / ДПИ,- Донецк, 1990,- 34с.
Статья поступила 15.01.2003
