Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Кочо, В. С. Физико-химические и теплофизические особенности современного мартеновского процесса

.pdf
Скачиваний:
6
Добавлен:
21.10.2023
Размер:
10.39 Mб
Скачать

Г Л А В А 7

НАГРЕВ КИПЯЩЕЙ ВАННЫ МАРТЕНОВСКОЙ ПЕЧИ

Несмотря на существенные успехи, достигнутые в изучении про­ цесса нагрева ванны мартеновской печи1 [44, 45, 119, 157, 205— 208], качественное и количественное описание особенностей тепло­ передачи в жидкой ванне нельзя считать завершенным.

Основная трудность заключается в том, что в кипящей ванне идет сложный массотеплообмен, зависящий от многих, пока не кон­ тролируемых параметров. Исследованиями [157, 207, 208] установ­ лено, что решающая роль в процессе теплопередачи в ванне при­ надлежит механическому перемешиванию ванны, а виртуальный коэффициент теплопроводности ванны [206] зависит непосредст­ венно от мощности процесса кипения ванны. Мощность процесса кипения определяется скоростью обезуглероживания, глубиной за­ легания основного фронта реакций окисления и температурой ме­ талла. Как показали лабораторные исследования на холодной мо­ дели [209], при значительном увеличении скорости окисления пропорциональное увеличение мощности перемешивания не дости­ гается, поэтому рост скорости нагрева металла в результате сопри­ косновения его с более нагретыми шлаком и факелом мартеновской печи при увеличении скорости обезуглероживания носит ограничен­ ный характер [206].

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

СКОРОСТЬ ДВИЖЕНИЯ ЭЛЕМЕНТАРНЫХ ОБЪЕМОВ КИПЯЩЕЙ ВАННЫ

Проведенная киносъемка процесса кипения мартеновской ванны показала возможность определения глубины основного (наиболее вероятного) фронта реакции окисления углерода. Зная глубину фронта реакции окисления углерода, можно достаточно точно рас­ считать мощность процесса кипения ванны, например, по уточнен­ ной формуле Кочо [210]. Зная мощность кипения всего объема ванны, можно рассчитать скорость движения элементарных объемов кипящей ванны.

Запишем в несколько измененной форме первое начало Даламбера для объема ванны.

F Tpd h - ^ - m ^ - d h ~ W d z = 0 ,

■ К а п у с т и н Е. А. Исследование тепло- и массообмена в мартеновских печах. Автореф. докт. дис. М., 1966.

158

где ■FTp dh — работа по определению сил трения на пути dh, эрг; FTp — сила трения, дин;

d2h

т 2 dh — работа по определению сил инерции на пути dh;

т — масса ванны, г;

Wdx — энергия, подведенная за время dx, эрг;

W — мощность процесса кипения ванны, эрг/с.

Принимая во внимание, что

 

 

дин,

 

(127)

где &тр = 'Птр/см2 — коэффициент трения,

дин•с

1

 

см

см^

г)тр— коэффициент внутреннего трения,

дин • с см

FTp — поверхность трения, см2; ю — скорость движения, см/с.

Запишем уравнение (126) после подстановки в него уравнения

(127)

 

 

 

 

 

 

 

 

 

 

£тр5тр«о d h + m - ^ - d h - W d x = 0.

(128)

Поскольку необходимо определить скорость

движения

элементар­

ного объема ванны, то разделим уравнение

(128) на массу ванны

 

 

k7t>—

u d fi Jr - ^ - d h - —

dx=0.

(129)

Разделив

переменные

и сделав

ряд

преобразований,

уравнение

(129) запишем

 

 

 

 

 

 

 

 

 

Л } р5 уд+ ^ -

У

- Г

уд= 0 ,

(130)

где k'тр

kj-p

приведенный коэффициент внутреннего трения, с;

g?M

 

 

 

 

 

 

 

 

 

рм — плотность металла, г/м3;

 

 

 

Fуд-

5

удельная поверхность трения элементарного объ­

“У

 

ема, см-1;

 

 

 

 

 

 

 

 

 

 

 

 

 

V — объем ванны, см3;

 

см/с2;

 

 

g — ускорение сил тяжести,

 

 

Гуд — удельная мощность процесса кипения, отнесенная

 

 

к единице элементарного объема ванны, см/с.

Скорость движения элементарного объема ванны (x>= dh/dx опре­ делим из условия устойчивого процесса кипения, т. е. когда в урав­

нении (130) можно допустить: d&/dx = 0

 

 

W,уд

м/с.

(131)

со =

^тр^уд

159

Из уравнения (131) следует, что скорость движения элементар­ ного объема ванны зависит, кроме мощности процесса кипения, от температуры ванны (изменяется k' ) и от геометрии элементарного

объема (шар, цилиндр, куб и т. д.), которая учитывается через £уд. Мощность всплывающих пузырей окиси углерода затрачивается на перемещения элементарных объемов ванны в вертикальном и го­ ризонтальном направлениях. Например, по результатам производ­ ственных исследований процесса массообмена в ванне большегруз­ ной мартеновской печи с применением изотопа Со60 средние значе­ ния соответственно для вертикальной составляющей скорости ыу равны 9 см/с, а горизонтальной сож 5,7 см/с, т. е. результирующее

значение скорости движения элементарных объемов со = УаЯ +

=

= 10,5 см/с (напомним, что, по данным тех же исследований,

ско­

рость всплывания пузырей окиси углерода в сталеплавильной ванне равна 20 см/с).

Поскольку формула (131) дает результирующее значение ско­ рости движения элементарных объемов со, а в процессах передачи тепла от греющего факела в глубь металлической ванны решающая роль принадлежит вертикальной составляющей скорости соу, то в дальнейших расчетах и выводах будем оперировать понятием вер­ тикальной составляющей со у, величину которой можно определить из выражения (131)

(132)

где k y — вероятностный коэффициент, учитывающий соотношение вертикальной и результирующей скоростей движения эле­ ментарных объемов (ky~0,9).

Расчеты, основанные на сопоставительном анализе процесса ки­ пения на холодной модели сталеплавильной ванны и в реальных ус­ ловиях выплавки стали в большегрузной мартеновской печи, пока­

зывают, что под элементарным

объемом

сталеплавильной

ванны

в первом приближении можно

понимать

шарик диаметром

1 мм

[ 211].

В дальнейших выводах значение W1Vi [см. формулу (132)] будем рассматривать как отношение общей мощности кипения к тому объ­ ему ванны (массы), в котором происходит основной процесс окис­ ления углерода (например, при донном характере кипения берется во внимание весь объем ванны, а при поверхностном кипении — лишь ее верхние слои).

ФИЗИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА НАГРЕВА ВАННЫ

СУЧЕТОМ ПЕРЕХОДНОГО СЛОЯ ШЛАК—МЕТАЛЛ

Врезультате движения элементарных объемов ванны обра­ зуется слой шлак—металл, роль которого для протекания физико­ химических процессов сталеплавильной ванны очень велика (см. гл. 2). Здесь рассматривается влияние переходного слоя шлак—

160

металл на процесс теплопередачи в жидкой ванне, когда основное количество тепла в ванну поступает от греющего факела, а не от экзотермических реакций окисления примесей расплава, как это имеет место при интенсивной продувке металла кислородом. Таким образом, будут рассмотрены закономерности процесса нагрева ванны в условиях интенсификации плавки подачей кислорода в фа­ кел или в ванну, но в умеренных количествах.

В процессе передачи тепла из рабочего пространства в объем жидкой ванны переходный слой шлак—металл играет большую роль. Во-первых, переходный слой шлак—металл как бы уменьшает шлаковый покров, экранирующий металл от греющего факела и, во-вторых, циркулирующие объемы металла пронизывают высоко­ температурный переходный слой и, нагреваясь в нем, вносят тепло

во внутренние

объемы ванны, ускоряя

процессы

теплопередачи

в жидкой ванне. Сознавая

 

 

 

всю сложность условий теп­

 

 

const

лопередачи в кипящей ван­

 

 

 

 

 

не и практическую невоз­

 

 

 

можность учета всего много­

1

 

 

образия процессов, проте­

“■■> tuj

 

кающих в ней, была сделана

V '

 

 

попытка

создания

физичес­

 

 

 

кой модели процесса, более

 

 

 

близко

отражающей реаль­

 

 

 

ный процесс.

математичес­

 

 

 

При

этом

А .

 

 

кая модель позволит про­

 

 

анализировать сложный про­

 

 

 

цесс теплопередачи

не для

Рис. 70. Физическая модель процесса нагрева

обнаружения

количествен­

кипящей ванны

мартеновской

печи

ных закономерностей и свя­ зей, на основе которых можно было бы дать практические рекомен­

дации,

выраженные в количественных

координатах

процесса,

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

влияния многих важ­

ных, но

не контролируемых параметров,

определяющих

скорость

и направление сложного процесса тепло- и массообмена в кипя­ щей ванне, чтобы принять во внимание некоторые качественные закономерности при необходимости управления тепловым и техно­ логическим режимами мартеновской плавки [212].

На рис. 70 показана’ физическая модель процесса и градиент температур в ванне. Температура металла /м. Индексом 0 обозна­ чены начальные условия. Теплопроводность жидкой ванны принята такой, что температура практически сразу выравнивается во всем объеме ванны, на подине всегда в большей или меньшей степени идет зарождение пузырей окиси углерода и ванна подвижна на­ столько, что вследствие тепло- и массообмена градиент температур по глубине практически отсутствует или им можно пренебречь при постановке данной задачи. Закон распределения температур в пе­ реходном слое шлак—металл принимаем линейным и температурное

И Зак. № 603

161

Рис. 71. Кривая изменения темпе­ ратуры жидкой стали и шлака при медленном перемещении спая тер­ мопары в вертикальной плоскости
(240-т печь)

состояние слоя может быть охарактеризовано tn. с (средним зна­ чением температур граничных слоев /ш_м и tM, где (ш- м температу­ ра переходного слоя на границе со спокойным слоем шлака).

Температура поверхности шлака, контактирующая с греющим факелом, принята постоянной (£ra= const). В процессе нагрева ванны изменяется лишь градиент в спокойном слое шлака. Таким образом, тепловой поток, пронизывающий спокойный слой шлака (уменьшением теплосодержания последнего пренебрегаем), обеспе­ чивает выполнение условий: tu. c= const, нагрев жидкого металла и

покрытие постоянных тепловых потерь, в статье которых учтен (для удобства математических преобразований) ре­ зультирующий эффект экзо- и эндотермичных реакций, протекающих в объ­ емах шлака и металла.

Характер распределения темпера­ тур по глубине ванны, положенный в основу физической модели, отвечает действительному распределению тем­ ператур в сталеплавильной ванне, ко­ торое удалось обнаружить в период доводки плавки при помощи установ­ ки для непрерывного измерения темпе­ ратуры стали (рис. 71).

Поясним физическую модель про­ цесса нагрева кипящей ванны (см. рис. 70). Переходный слой шлак—ме­ талл является одновременно теплооб­ менником, в котором подогреваются циркулирующие через него элементар­ ные объемы металла, и непосредствен­ ным проводником теплового потока, направленного от шлака к металлу, ве­ личина которого определяется гради­ ентом температур на границах слоя и его теплопроводностью. По указанной причине тепло интенсивно подводится

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

РАСЧЕТ НАГРЕВА КИПЯЩЕЙ ВАННЫ

Принимаем, что в переходном слое содержится металл в коли­

честве

 

Оп. c=-j- KFpu,

(133)

где F — площадь пода, м2.

 

162

Под переходным слоем находится металл в количестве

(Jм G п. (

( 134)

где G*— садка печи, т.

Продолжительность пребывания элементарного объема метал­ лической ванны в переходном слое тц (время цикла) принимаем

V 2/?ш_м/(by с.

(135)

Коэффициент 2 означает, что элементарному объему необходимо преодолеть двойной путь в переходном слое, чтобы снова вернуться в объем металлической ванны. Не исключена возможность отклоне­ ния величины этого времени для некоторых объемов. В этом смы­ сле величина цикла носит вероятностный характер.

Расход металла, циркулирующего через переходный слой шлак—металл в единицу времени

 

Ga=jFpMuy т/с.

(136)

Температура, с которой элементарные объемы металла покидают

переходный слой,составляет

 

 

<„-ш= *п.с —(*п.с— 4 )ех р (—Атц) °С,

(137)

где

tM— температура, с которой поступают

элементарные

 

объемы в переходный слой из объема металличе­

 

ской ванны,°С;

 

А

1

 

постоянный коэффициент, с-1;

 

7?РмСм

R — тепловое сопротивление элементарного объема ме­ таллической ванны, зависящее от теплопроводно­ сти металла, геометрии элементарного объема и условий теплопередачи в слое, м2°С/Вт;

см — теплоемкость металла, ккал/ (кг • °С). Особенность математического описания процесса теплообмена

согласно предложенной модели заключается в том, что после каж­ дой циркуляции объемов металла через переходный слой (продол­ жительность цикла) система характеризуется новыми начальными условиями, вследствие этого вывод уравнений, описывающих про­ цесс нагрева, будет основываться на методе математической де­ дукции.

Запишем изменения температуры металлической ванны для ряда циклов.

Цикл I. Начальные условия: ^м = ^м„; ^н-м, =Ка-м0. Необходимо

найти температуру металла

в конце цикла.

11*

163

Уравнение баланса тепла для металлической ванны (под пере­ ходным слоем шлак—металл) за период одного цикла:

 

 

 

 

-

<?„) F \ = G McJ»-

 

(138)

 

 

 

1

пш-м

 

) Вт/м2>

 

 

( 1 3 9 )

 

 

 

 

 

 

 

 

 

 

 

Ч = Ч , - Ч = Р м %

см ( ^ . с - у [ 1 - ( ~

е х Р Ч ) ] в т/м 2,

 

(140)

где q

— тепло, подведенное, в металлическую

ванну вследствие

 

теплопроводности переходного слоя;

 

 

 

Уш-м — коэффициент

молярной

теплопроводности переходного

Aq

слоя;

 

 

 

 

 

 

 

 

 

 

— тепло, подведенное в металлическую ванну в результате

 

циркуляции элементарных объемов через переходный

 

слой;

 

 

 

циркулирующими

элемен­

А<7^ — тепло, унесенное из ванны

 

тарными

объемами, которые

имеют

температуру

;

А<7" — тепло, подведенное в ванну теми же объемами металла,

 

подогретыми в переходном слое

 

за

время цикла

т до

 

температуры 1м-ш [см. формулу (137)];

 

 

qn— постоянные потери с учетом

результирующего эффекта

 

экзо- и эндотермической реакций,

протекающих в объ­

 

емах шлака и металла.

 

 

 

температуры

металла

Из уравнения (138)

находим значение

в конце первого цикла:

 

 

 

 

 

 

 

 

 

 

 

 

 

^М,— ^М0—|—Д^м,> °С,

 

 

 

 

(141)

ДА,

 

- ( * ш - м „ -

У

+ " у РмСМ

с -

У

[ !

-

( -

е х Р А \ ) ] -

Яп

Ъи-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(142)

Цикл II. Начальные условия: ^М= ^М|;

 

 

=

— А7

. Напи­

шем аналогичные (138) — (140) балансовые уравнения для цикла II:

 

 

G j t

+

( g l2+ А<7ц2-

<7П)

 

=

GMc J u -

 

(143)

 

 

 

 

Аш_» •(^ш_м2— 4,,)

Вт/м2;

 

 

(144)

 

Д^ц,=“уРмСи(^п.с—^м.) [1 — (—ехр Лхц)]

Вт/м2.

 

(145)

Значение температуры металла в конце этого цикла t

 

 

 

 

 

 

^м2 14 .,“4 Д^Мг>

 

 

 

 

 

(146)

 

ДШ-М

/jt

4,,) + “ уР м См (* п .с

 

 

 

 

 

 

 

 

Д4

"Uш -м

« Ш...........- M j -

- * м .)

П

-

( ~

е Х Р А Н 1 -

<7П

 

 

 

 

 

 

 

 

 

 

 

^Ц.

 

 

 

 

 

 

 

 

 

 

 

 

(147)

164

Цикл III. Начальные условия:

tM= t

;

t

=t

— At

— At . Балансовые уравнения для цикла III:

 

 

 

0»CJ* ,+

К +

Ч , , -

?„) F\

= G*cJ*j

(148)

^з =

Т ~

(*ш-м,-*м,)

Вт/м2;

 

(149)

 

 

— (— expАхц)]

Вт/м2.

(150)

Температура металла в конце третьего цикла будет:

 

 

^Мз

j

^^М3'

 

 

 

(151)

/Т Т Л (*ш-м, - ^м2) + %Рмем (V

с - 'м,)

Р -

( -

ехР ^ ц)1

- qn

Д*« = -

 

 

 

 

 

 

F^.

 

 

 

 

 

 

 

(152)

Пользуясь методом математической дедукции, для любого цикла запишем значения всех интересующих нас параметров.

Температура металла после n-го цикла tMn:

Al„-- ^М,Д-А^М,-|‘^^М2Ч_^^Мз“Ь • • •

°С,

(153)

 

 

i = п

 

 

 

 

L = t

+ У

ДА,. °С,

 

(154)

 

Мп

Мп 1

 

 

 

t

-- t

 

 

 

 

Hin

 

 

 

/и£щ—M(?ш-мп

T- wypMcM(^n c

[1 (■

exp Л-Сц)]

gn

 

 

 

 

 

/% . (155)

Приращение температуры металла после n-го цикла:

д Ш М (Аи-м„

^ п - \ ) + “ уРмсм (Ai.c — Al„_!)

[1

,4ш—м

ц ( - е*рЛтц>1 - g "

£

/гТц. (156)

ДА, ------------------

Температура переходного слоя шлак—металл на границе с псевдоспокойным слоем шлака после га-го цикла Ад-м :

о

"^ м 3

■ • *

_ 1 °С; (157)

 

i =

п

 

 

t

t - У

А/

°г

(158)

i ==1

165

В уравнении (155) для t№ji раскроем значение

(/ш- м„ —

при помощи уравнений (154)

и (158):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i= п —1

 

 

 

 

 

 

 

 

 

 

 

 

 

—^м„ + 2

^^мг;

 

 

 

 

 

(159)

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1=п —1

 

 

/= я —1

 

 

 

 

 

ш-мп ^m,;_i)

 

м„

 

1

 

А^м;

^м0"— ^ а/м.;

 

 

 

 

 

 

 

 

 

 

 

 

 

i—ii—1

1

 

 

 

 

 

 

 

 

 

 

— f

 

 

 

 

 

 

 

 

 

 

 

 

/

— /

 

_ /

 

—9

^

 

°С.

 

 

(160)

 

‘'ш—глп

 

^м/2_ 1

—Мд *м0

^

 

 

 

 

Решая одновременно уравнения (155), (159), (160), запишем

общее выражение для температуры металла после n-го

цикла tMn-,

 

 

 

 

 

 

 

 

i=n —1

 

 

 

 

 

 

 

 

 

 

 

tUn—4i0+

^

 

 

^м/ +

 

 

 

 

 

 

и

MI i

 

_

1м0~

2

|

 

 

)

+

°\РмСм (

с -

Ч

-

 

 

 

I 1ш-м0

 

 

 

 

 

 

 

-

 

2

ДЧ )

[! -

(-ехр^Тц)] — да

 

 

 

 

 

 

 

 

 

 

 

 

 

GuCM

 

 

 

 

-------

Лц,

(161)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1Мп--

(^ш-м0 - ^mq)

+

%РмСм (*„■ с -

*„,) t 1 ~

( ~ ехр Л тд)1

"

? -

 

 

 

 

 

 

 

с м<?м

 

 

 

 

 

 

■X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/=Я—1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2 —

 

 

«уРмСм[1 —(—ехрЛхц)]

\

 

 

 

 

 

X I

 

 

Лш"м

 

 

а мсм

 

 

 

-/ “Чц I.

 

 

(162)

Принимая во внимание,

что

^•0= 4^ ш -М о + ^М„)

(СМ‘

РИС' 70’

 

 

 

 

 

 

 

 

71), запишем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(*ш-м0 -

г‘м„) {

 

 

+

0,5мурмсм [1 - ( —ехрЛхц) ] |

-

дп

 

CL

 

 

 

\

Лш-м

 

 

 

 

 

 

 

)

 

 

I 1ц,

---------------------------------------------

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Gy,c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

мсм

 

 

 

 

 

 

(163)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ь = 1

 

 

 

 

+

шуРм^м [1 “

( —ехр Лтц)]

Ц*

 

 

(164)

 

 

 

 

 

 

 

'п-‘мсмг

 

 

 

1

 

 

 

 

 

 

п ш - м ____________________________________________________ Р -

 

 

 

 

 

166

После этого уравнение

(162) запишется в виде

 

 

('= /!—1

 

 

-\- Ь ^

Atu;.

(165)

 

1

 

 

Принимая во внимание, что

 

 

/= « —1

 

 

 

2

' ^мо> ^мл_1==^мл

 

1

 

 

 

уравнение (165) можно записать

 

 

^ п= а

(^мл

^ м п) °С.

(166)

Учитывая, что при п-> оо Д(Мп->0

(см. рис. 72,

где представ­

лена геометрическая интерпретация закона изменения температуры металла во времени), уравнение (166) запишется

 

 

*мп= а + ЫМп- Ы щ °С,

 

(167)

откуда и определится значение tMn для п -у оо

 

 

 

 

а —btМ.0

о

с.

 

(168)

 

 

1

 

 

После подстановки значений а и b из уравнений (163) и

(164)

получаем

 

 

 

 

 

 

I

t'1"-*- ~ '“•> (~ fa ~ +

П - (~ ехр АчУ1} -

т

(169)

 

2 ^ш~м- + (йурмсм [1 -

(—ехрЛтц)]

 

 

 

 

ЛШ-М

 

 

 

 

Для 1

оо уравнение (169) запишется в виде:

 

 

 

 

4«„= *M„+(l---- —) х

 

 

(*ш-м0- ^м0) (

+ °.5“уРмС [1 - (—ехрЛтц)] j-

qn

(170)

х -

Am—i

 

 

 

 

 

м- + « jfPmCm[1 -- (—ехр^Тц)]

 

 

где jV— постоянный коэффициент — определяется из граничных ус­

ловий при л =

1, /

= t

+ At .

г

1 Ml

Mo

Ml

Принимая во внимание уравнение (142) для определения Д(м ,

находим значение N

 

Омсм —Т'хц | 2

+ юуРмсм [1 —(-ехр .4тц)] |

N=

(171)

167

Соседние файлы в папке книги из ГПНТБ