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

t, мин

Состав металла, %

T, К

nFe, %/мин

nSi, %/мин

nС, %/мин

nСr, %/мин

n T, К/мин

C Cr Si

Период I

0 1,160 17,29 0,400 1748,00 0,039 0,044 0,088 1,190 0,0095 0 11,84
1 1,150 17,29 0,356 1759,84 0,039 0,044 0,088 1,300 0,0095 0 11,82
2 1,141 17,29 0,312 1771,66 0,039 0,044 0,088 1,410 0,0095 0 11,80
7 1,085 17,29 0,090 1830,77 0,039 0,044 0,088 2,0 0,095 0 11,72
8 10,75 17,29 0,046 1842,49 0,039 0 0,442 2,210 0,0474 0 3,34
22 0,420 17,29 0,046 1889,06 0,039 0 0,442 1,020 0,0474 0 3,28
23 0,373 17,29 0,046 1892,34 0,038 0 0,425 0,910 0,1380 0,0456 0 2,91

Период II

24 0,327 17,29 0,046 1895,25 0,038 0 0,425 0,807 0,0560 0,0456 0 2,91
25 0,282 17,29 0,046 1898,16 0,038 0 0,425 0,700 0,0313 0,0313 0,046 6,38
26 0,251 17,24 0,046 1904,54 0,038 0 0,425 0,651 0,0250 0,0250 0,068 7,40
32 0,132 16,73 0,046 1959,25 0,038 0 0,425 0,525 0,0148 0,0148 0,1016 10,36
33 0,118 16,63 0,046 1969,61 0,020 0 0,221 0,507 0,0280 0,0237 0 –0,57

Период III

34 0,094 16,63 0,046 1969,04 0,020 0 0,221 0,400 0,00182 0,0182 0,018 0,07
35 0,076 16,61 0,046 1969,70 0,020 0 0,221 0,323 0,0130 0,0130 0,035 1,93
39 0,0375 16,42 0,046 1980,84 0,020 0 0,221 0,174 0,0057 0,0057 0,060 3,72
39,5 0,0346 16,39 0,046 1982,70                

Примечание. В таблицу внесены округленные значения; промежуточные значения в каждом из периодов не приводятся.


Как видно из приведенного ниже сопоставления, расчетные значения (знаменатель) [С], [Сr] и T в конце каждого периода достаточно близки к экспериментальным (числитель):

Н.п.

С, % Cr, %

Т, К

1,16 17,29 1748
I 0,41/0,37 16,9/17,29

1923/1892

II 0,12/,0,118 16,7/16,3

1998/1969

III 0,047/0,0346 16,45/16,39

2003/1982

Примечание. Н.п. - начало продувки, I–III - моменты окончания периодов продувки.

Примечание. Образование оксида с формулой Cr3O4 экспериментально не зарегистрировано. Формула принята только потому, что она обеспечила наилучшее согласование между собой равновесных опытных значений [Cr], [C], и Т при рСО = 1 атм. Значение DGº образования Cr3O4 получено обратным расчетом из 9.2.

83. Распределение серы
между металлом и шлаком

Задача. Определить количество шлака, необходимое для десульфурации 1 т стали в ковше при 1923 К.

Исходные данные.

1. Состав шлака: 50 % СаО; 50 % Al2O3. 2. Содержание алюминия и серы в конечном металле: 0,010 % [Аl], 0,005 % [S]. 3. Содержание серы в исходном металле 0,015 %.

Теория. В условиях сталеварения обычно валентность серы nS = -2 (см. задачу 24), поэтому растворимость серы в шлаке определяется соотношением

(S) = ,

где S - концентрация серы в шлаке.

Константу K(S) называют сульфидной емкостью шлака, и если массовое содержание серы выражено в процентах, то ее обозначают символом cS (см. задачу 29):

(% S) = . (9.27)

Окислительный потенциал системы ( pO2) и давление паров серы ( pS2), необходимые для расчета равновесного содержания серы в шлаке (% S), могут быть найдены по заданным значениям [Аl] и [S] из реакций

2[Аl] + О2 = (Аl2О3);

[S] = S2 .

Сульфидная емкость cS шлака заданного состава оценивается при помощи одного из теоретических методов или принимается на основании опубликованных экспериментальных данных.

Решение. На основании справочных данных принимаем значения DG реакций:

О2(г) = [% О]Fe ; DGº = -117200 - 7,89Т, Дж;

lg K[O] = + 0,151;

K[O] (1923 К) = 2154; (9.28)

S2 (г) = [% S]Fe ; DGº = -7200 – 10,26Т, Дж;

lgK[S] = + 0,535;

K[S] (1923 К) = 309; (9.29)

2[Al] + 3[O] = Al2O3 (т) ; DGº = -1199700 – 372,2Т, Дж;

lgK4 = - 19,431;

K4 (1924 К) = 1,36×1013. (9.30)

Чтобы не усложнять примера, коэффициенты активности fC , fO , fAl и gAl2O3 принимаем равными 1. Решая совместно уравнения
(9.27) – (9.30), определяем содержание серы в шлаке после обработки металла:

,

Пользуясь данными [4, с. 118], получим

lg cS = -2,37, cS = 4,3×10-3.

Вычисляем содержание серы в конечном шлаке:

(% S) = 43×10-3 = 0,24.

Из 1 т металла должно перейти в шлак [(0,015 - 0,005)×1000 кг]/100 =
= 0,1 кг серы. Для этого потребуется шлака (0,1×100)/(% S) =
= (0,1×100)/0,24 = 41,7 кг.

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

84. Дегазация стали продувкой
жидкого металла в ковше аргоном

Задача. Определить количество аргона, необходимого для продувки жидкой стали ШХ15 в ковше при 1873 К через фурму или пористые элементы футеровки для снижения концентрации водорода.

Исходные данные. 1. Масса металла в ковше g = 40 т. 2. Исходная и конечная концентрации водорода равны [H]исх = 6×10-4 %,
[Н]кон = 3×10-4 %. 3. Константа (KH) равновесия реакции Н2 = 2[Н] при 1873 К равна 0,73×10-5.

Теория. В соответствии с формулой Геллера [56] уравнение, описывающее зависимость концентрации водорода от количества продуваемого аргона VAr3), имеет вид

. (9.31)

После интегрирования (при 1873 К) получим

VAr = 0,82 . (9.32)

Решение. Коэффициент активности водорода в жидкой стали ШХ15 (1 % С, 1,5 % Сr) вычисляем по параметрам взаимодействия (см. задачу 62): fH = 1,07. Тогда VAr = [(0,82×10-3×40)/(1,072)][( ) -
– ( )] = 47 м3.

Примечания. 1. Расчет сделан с учетом того, что Н2 в пузыре аргона находится в равновесии с водородом, растворенным в жидкой стали. В связи с этим размер пузырей аргона, так же как и продолжительность продувки, роли не играют. Для водорода обычно наблюдают удовлетворительное соответствие между опытными и расчетными данными, что указывает на достижение равновесия ( ) и отсутствие диффузионных сопротивлений в объеме металла. Для азота в большинстве случаев равновесие не достигается (a = 0,2–0,4). Это связано с тем, что при продувке азотом процесс лимитируется массопереносом азота из объема расплава к поверхности пузыря аргона. Другими словами, равновесие между азотом в пузыре аргона с азотом, растворенным во всем объеме металла, не достигается, хотя между поверхностным слоем и пузырем равновесие может иметь место.

2. Чтобы ввести в металл требуемое количество аргона (47 м3) за реальное время продувки (10 мин), необходим расход аргона 4,7 м3/мин. Это соответствует интенсивности продувки
0,12 м3/(т×мин). Такие интенсивности при продувке в ковше недостижимы из-за слишком бурного кипения металла. Поэтому удаление в ковше водорода из металла без вакуума практически невозможно. Процесс можно реализовать только в аргоно-кислородном конвертере при продувке с интенсивностью 1 м3/(т×мин).

85. Условие пузырькового
выделения газов из стали

Задача, Вычислить концентрацию водорода в стали 30ХГСА, при которой в ходе вакуумирования возможен пузырьковый механизм выделения водорода с заданной глубины металлической ванны.

Исходные данные. 1. Состав стали: 0,3 % С; 1 % Сr; 1 % Мn; 1 % Si. 2. Плотность расплава r = 7,0×103 кг/м3. 3. Поверхностное натяжение металла s = 1,68 Дж/м2. 4. Начальный радиус пузыря r = 5×10-4 м. 5. Температура 1873 К. 6. Остаточное давление в вакуумной камере рост = 4×102 Па. 7. Глубина ванны h = 0,5 м.

Теория. Процесс растворения водорода в стали (или выделения его из металлической ванны) может быть описан следующим уравнением [3]:

H2 (г) = [H]Fe) ; lgK = (-1900/T ) - 1,577. (9.33)

Запишем выражение константы равновесия процесса (9.33) применительно к условиям данной задачи:

K = [H] , (9.34)

где рH2 - давление водорода в пузыре. Тогда искомая концентрация [H] равна

[H] = . (9.35)

Величина рH2 в пузыре

рH2 = pocт + pgh + 2s/r, (9.36)

где pgh - ферростатическое давление (g - ускорение свободного падения, равное 9,8 м/с2).

Решение. Вычислим значение K: lgK = (-1900/1873) - 1,577 = = -2,59; K = 0,0026. Определим fH по параметрам взаимодействия (см. задачу 16); lg fH = 0,04; fH = 1,1. По формуле (9.36) вычислим рH2 в пузыре: рH2 = 4×102 + 7,0×103×9,8×0,5 + 2×1,68/(5×10-4)2 = 41420 Па » » 0,41 атм. По формуле (9.35) определим [H]:

[H] = (0,0026/1,1) = 0,0015 %.

Таким образом, при вакуумировании ванны жидкой стали 30ХГСА глубиной 0,5 м водород может выделяться с пода в виде пузырей, если концентрация его в ванне [H] ³ 0,0015 %.

86. Удаление серы при обработке металла сульфидообразующими реагентами

Задача. Вычислить максимально возможную (равновесную) концентрацию серы в стали при обработке расплава карбидом кальция. Рассчитать науглероживание стали и минимально необходимый (теоретический) расход карбида кальция.

Исходные данные и рекомендуемые допущения. 1. Состав стали: 0,3 % С; 0,3 % Si; 0,4 % Mn. 2. Температура 1873 К. 3. Активности карбида кальция и образующегося сульфида принять равными 1.
4. Давление пара кальция принять равным 1 атм, что приблизительно соответствует условиям обработки стали порошком карбида кальция в струе инертного газа при атмосферном давлении.

Теория. Задача сводится к расчету равновесия реакции

СаС2 (т) + [S] = CaS(т) + 2[С]. (9.37)

Значение стандартной энергии Гиббса для этой реакции может быть найдено путем алгебраического суммирования DGº соответствующих реакций образования веществ из элементов (см. табл. П1) и реакций растворения элементов в жидком железе (см. табл. П2). Зная , определяем выражение для константы реакции (9.37):

lg K = , (9.38)

из которого вычислим равновесную концентрацию серы. Коэффициенты активности серы и углерода можно определить при помощи параметров взаимодействия (см. табл. П3). Для определения равновесной концентрации углерода в металле рассчитаем равновесие реакции

СаС2 (т) = Са(г) + 2[С]; DGº = -RT ln , (9.39)

приняв заданные значения aCaC2 и рCa . Для определения минимально необходимого (теоретического) расхода карбида кальция, необходимо воспользоваться стехиометрическими соотношениями реакции (9.37), из которых следует, что для удаления 1 моля серы (32 г) требуется 1 моль СаС2 (64 г). В металл при этом переходят 2 моля углерода (24 г).

Решение. Определим равновесную концентрацию углерода в металле. Для определения DGo реакции (9.39) воспользуемся данными табл. П1 и П2:

СаС2 (т) = Са(г) + 2С(гр); DGº = 213000 - 61,0T Дж

(гр) = 2[С]; DGº = 45200 - 84,6Т Дж

СаС2 (т) = Са(г) + 2 [С]; DGº = 258200 - 145,6T Дж (9.40)

Из уравнения (9.40), приняв Т = 1873 К, aCaC2 = 1 и рCa = 1 атм, получим = (-258200 + 145,6×1873)/(19,14×1873) = 0,4045; аC = = 1,6 %. С помощью табл. П3 получим lg fC = 0,033; fC = 1,08, следовательно, [С] = аC /fC = 1,48 %. Полученный результат показывает, что десульфурация заданной стали карбидом кальция должна сопровождаться науглероживанием. Соотношение этих процессов определим ниже.

Для нахождения представим реакцию (9.37) как сумму следующих процессов:

СаС2 (т) = Са(г) + 2С(гр); DGº = 213000 - 61,0T Дж (9.41)

[S] = S2 (г); DGº = 72000 + 10,25Т Дж (9.42)

Са(г) + S2 (г) = CaS(т); DGº = -692000 + 186,0Т Дж (9.43)

(гр) = 2 [С]; DGº = 45200 - 84,6Т Дж (9.44)

СаС2 (т) + [S] = CaS(т) + 2[С]; = -361800 + 50,65T Дж (9.45)

При 1873 К = -267000 Дж. Согласно выражению (9.38) с учетом равенства aCaC2 = aCaS = 1 получим [С]2/fS [S] = -267000/(19,14×1873) =
= 7,448. В соответствии с табл. П3 lg fS = 0,02; fS = 1,05.

Подставив найденные значения fC и fS в (9.39), получим

lg[(1,082×1,482)/(1,05×[S])] = 7,448; lg[S] = 7,448 + 2×0,033 + 2×0,17 -
- 0,02 = -7,067 [S] » 1×10-7 %.

Этот результат показывает, что десульфурация в рассматриваемых условиях может обеспечить очень низкий уровень концентрации серы в металле. Из приведенных выше стехиометрических соотношений видно, что удаление из металла каждой 0,01 % серы требует не менее 0,01×(64/32) = 0,02 % (0,2 кг/т) карбида кальция. При этом в металл переходит не менее 0,01×(24/32) = 0,0075 % углерода. Суммируя полученные результаты, можно сделать вывод, что карбид кальция является хорошим десульфуратором. Теоретически при расходе 0,2 кг карбида кальция на 1 т стали можно удалить 0,01 % серы при сравнительно небольшом науглероживании стали (на 0,0075 %). На практике в связи с угаром кальция и действием кинетических факторов, препятствующих достижению равновесия, расход карбида кальция и науглероживание должны быть выше.

87. Скорость удаления серы
из расплава в газовую фазу

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

Исходные данные и рекомендуемые допущения. 1. Масса расплава т = 1000 кг; состав: 0,1 % [С]; 0,02 % [S]; T = 1873 К. 2. Параметры вакуумной обработки: остаточное давление робщ = 1 Па (10-5 атм); производительность откачки 15×10–5 м3/с; [Н] = 0,0002 %. 3. Параметры обработки водородом: давление р = 101,325 кПа (1 атм), расход водорода 0,015 н×м3/с. 4. Принять, что максимальная скорость десульфурации соответствует насыщению газовой фазы серосодержащими газами (до их равновесных давлений). Концентрацию кислорода при вакуумной обработке определить из условия равновесия с углеродом при рCO = 10,1325 кПа (0,1 атм). При обработке водородом принять [О] = 0,0022 %.

Теория. Для решения задачи необходимо определить максимально возможное (равновесное) содержание серы в газовой фазе, а затем, по известной скорости газообмена, рассчитать максимально возможную скорость удаления серы из металла. Равновесное содержание серы в газовой фазе определяется термодинамическими характеристиками испарения серы и реакций образования серосодержащих газов: H2S, CS, CS2 , COS , SO2 и др. Эти характеристики могут быть определены по табл. П1 и П2. Результаты расчетов с использованием данных работ [9, 62] представлены в табл. 1 и для некоторых конкретных условий – на рисунке.

Таблица 1

Термодинамические характеристики реакций образования серосодержащих газов (по данным [9, 62], расчеты А.В. Богданова)

Реакция

Коэффициенты в уравнении
lgK = (А/Т ) + В

А В
1. [S] = 1/2S2(г) -3760 -0,5355
2. [S] + 1/2Н2 (г) = HS(г) -7718 0,2828
3. [S] + [H] = HS(г) -5812 1,874
4. [S] + H2 = H2S(г) 973 -3,12
5. [S] + 2[H] = H2S(г) 4786 0,066
6. [S] + [C] = CS(г) -13672 1,845
7. 2[S] + [C] = CS2(г) -5714 -3,002
8. [S] + [C] + [0] = COS(г) 2003 -2,397

Окончание табл. 1

Реакция

Коэффициенты в уравнении
lgK = (А/Т ) + В

А В
9. 2[S] + [O] = S2O(г) -4731 -4,510
10. [S] + [O] = SO(г) -6803 -0,426
11. [S] + 2[O] = SO2(г) 2662 -4,563
12. [S] + 3[O] = SO3(г) 1528 -9,450
13. [S] + [Si] = SiS(г) -12908 2,735
14. 2[S] + [Si] = SiS2(г) -8113 -1,873

Равновесное давление серосодержащих газов pi (цифры на кривых соответствуют номерам в табл. 1, для сплавов на основе железа
при 1873 К) при разных расчетных условиях:
2, 4 - рH2 = 1 атм; 3, 5 - [Н] = 0,001 %; 6, 7 - [С] = 0,1 %;8 - [С] = 0,1 %;
[О] = 0,23 % (условно на поверхности металла); 8 - [С] = 0,1 %; [О] = 0,022 %;
9, 10, 11, 12 - [O] = 0,23 %; 9', 10', 11', 12' - [O] = 0,022 %; 13, 14 - [Si] = 0,1 %

В качестве примера расчетов, проведенных при построении графика, определим pH2S при T = 1873 К и [H] = 0,001 %.

По данным табл. 1

lg K = lg[pH2S /(aS[H]2)] = (4786/1873) + 0,066 = 2,621.

Отсюда

lg pH2S (атм) = 2,621 + 2lg(0,001) + lg aS = -3,379 + lg aS . (9.46)

Полученное уравнение соответствует прямой 5 на рисунке. Определив таким способом парциальные давления серосодержащих газов (р i) в заданных условиях, можно вычислить суммарную мольную долю серы в газовой фазе {S}S по формуле

,

где nS (i) - число атомов серы в молекуле газа i.

По известным значениям расхода газов, удаляемых из системы (q, н×м3/с), массы металла (m, г) и pобщ (Па) можно определить скорость десульфурации металла, %×с-1:

.

Решение. Коэффициент активности серы в заданном расплаве равен lg fS = -0,028×0,020 + 0,11×0,1 = 0,01044; fS = 1,024. Активность серы aS = 0,020×1,024 = 0,0205. Содержание кислорода при вакуумной обработке определим, приняв величину т = [% С]×[% О] =
= 0,0022×0,1 = 0,00022 (в 10 раз ниже, чем при атмосферном давлении, поскольку эффективное давление СО при вакуумной обработке обычно равно 0,1 атм). Получим [% О] = 0,00022/[% С] = 0,0022 %.

Равновесные давления серосодержащих газов (р i) находим с помощью уравнений табл. 1 или по рисунку (с корректировкой на более низкие содержания кислорода и водорода). Для обработки в атмосфере водорода получим

lg pH2S = lg K + lgaS + lg рH2 = - 3,12 - 1,6882 + 0 = -4,2887;
pH2S = 5,142×10-5 атм.

Скорости десульфурации в заданных условиях:

при обработке водородом

=
= 1,181×10-7 % ×с-1(0,00042 %×ч-1);

при вакуумной обработке

=
= 2,4×10-6 %×с-1(0,00864 %×ч-1).

Результаты расчетов значений р i , {S}i и {S}S приведены в табл. 2.

Таблица 2

Результаты расчетов значений р i , { S } i и { S } S

Газ

Для атмосферы H2

Для вакуума

р i ×105, атм {S}i ×105 р i ×105 aтм {S}
H2S 5,142 5,142 0,0343 0,0343
COS 0,0514 0,0514 0,0514 0,0514
HS 0,2978 0,2978 0,0242 0,0242
SO 0,0004 0,0004 0,0004 0,0004
CS 0,0007 0,0007 0,0007 0,0007
S2 0,0004 0,0008 0,0004 0,0008
 

{S}S = 5,4931×10-5

{S}S = 0,1118

Примечание. Размерность {S}å – моль S/моль смеси газов.

Примечание. Скорости десульфурации в реальных условиях внепечной вакуумной обработки и даже вакуумной плавки при указанных параметрах существенно ниже, так как равновесие между расплавом и газовой фазой не достигается вследствие кинетических условий, поэтому на практике рассмотренные процессы не применяют. Для удаления серы разработаны другие, более эффективные методы: обработка шлаками, продувка аргоном в смеси с порошками (извести, силикокальция, РЗМ и др.).

88. Раскислительная способность шлака

Задача. Рассчитать количество синтетического шлака, необходимое для рафинирования 1 т стали от кислорода при 1873 К.

Исходные данные. 1. Состав шлака: 50 % СаО; 50 % Аl2О3 2. [О]нач =
= 0,05 %; [О]кон = 0,005 %.

Теория. Шлак, находящийся в равновесии с металлом, содержит определенное количество FeO в соответствии с реакцией

[Fe] + [О] = (FeO), DGo = -120998 + 52,37T Дж/моль,

K = a(FeO) /a[Fe] a[O] . (9.47)

Если металл, отделенный от равновесного шлака, привести в контакт с новым шлаком, не содержащим FeO, то в результате реакции (9.47), понизится содержание кислорода в металле, а в шлаке появится FеО, т.е. произойдет раскисление металла шлаком.

Решение. Чтобы не усложнять примера, принимаем gFеO , f[O] и a[Fe] равными 1: K = х(FeO)/[% О], lnK = -DGo/RT; lnK = (14553/T ) – 6,3; K1873 К = 4,3.

К концу обработки равновесное содержание FeO в шлаке окажется равным xFeO = K [% О] = 4,3×0,005 = 0,0215 моль-1. Для этого из металла в шлак должно перейти кислорода: тO = [(0,05 - 0,005)×1000]/100 = 0,45 кг или FeO: mFeO = тOМFeO /МO = (0,45 – 72)/16 = 2,025 кг, где МFeO и МO - молекулярные массы FeO и кислорода. В 100 г исходного шлака (50 % СаО и 50 % Al2O3) содержится молей СаО и Al2O3 : nCaO = 50/56, nAl2O3 = 50/102. Соотношение nCaO /nAl2O3 = (50×102)/(56×50) = l,82 остается постоянным в процессе обработки металла. После обработки шлак содержит молей FeO: nFeO = 2025/72 = 28,125 молей.

Решая совместно систему уравнений

x = = 0,0215;

nCaO /nAl2O3 = 1,82;

nFeO = 28,125.

получим nCaO = 826,507 молей; nAl2O3 = 454,125 молей.

Масса шлака

mш = nFeO mFeO + nCaO mCaO + nAl2O3 mAl2O3 ,
mш = 2,026 + 826,607×0,056 + 454,125×0,102 = 94,6 кг

или 9,5 % от массы металла, что заметно превышает количество синтетического шлака в практике обработки стали в ковше (3,5–5 %).

89. Расчет равновесного состава металла
и шлака с применением ЭВМ

Задача. В ковш выпустили из ДСП-100 низкоуглеродистую сталь. Во время выпуска ее раскислили ферросилицием и алюминием. Для дополнительной десульфурации металла добавили известь. Определить количество и состав конечного металла и шлака, считая, что в системе устанавливается равновесие.

Исходные данные. 1. Жидкая сталь перед выпуском в ковш содержит 0,18 % [О], масса 100 т. 2. Ферросилиций ФС65 содержит 65 % Si, 35 % Fe, масса 0,1 т. 3. Алюминий вторичный АВ92, содержит 92 % Al, 8 % Si, масса 0,1 т. 4. Известь 100 % СаО, масса 0,3 т. 5. Температура 1873 К. 6. Валентности (степени окисления) элементов в шлаке считать постоянными, их значения: nFe = +2, nCa = +2, nSi = +4, nAl = +3, nO = -2. 7. Изменения энергий Гиббса образования соединений по реакциям

[Ме] + (9.48)

и выражения для определения констант равновесия химических реакций:

= -120850 + 23,3Т; (9.49) lgKFe = (6317/T ) - 1,22; (9.50)

= -668230 + 165,1Т; (9.51) lgKCa = (34951/T ) - 8,63; (9.52)

= -701550 + 140,6Т; (9.53) lgKSi = (36672/T ) - 7,35; (9.54)

= -663810 + 121,9Т; (9,55) lgKAl = (34669/T ) - 6,37. (9.56)

8. За стандартные состояния в пункте 7 приняты чистые жидкие металл и его оксид. Если их температуры плавления больше 1873 К, то стандартное состояние для условий задачи гипотетическое. Для кислорода за стандарт принят его гипотетический бесконечно разбавленный раствор в железе с х[O] = 1, т.е. при х[O] ® 0 коэффициент активности g[O] = 1, a[O] = х[O] (см. [3, с. 23–24]; g[O] в тексте – аналог j[O] в [3]).

Теория. Решение основано на сочетании материального баланса и термодинамического расчета. Термодинамический расчет позволяет установить составы фаз в равновесии.

Распределение i-го элемента по кислороду может быть представлено константой (см. задачу 24).

(9.57)

где ni - валентность (степень окисления) i-го элемента в фазе; х[i] , x(i) - атомные доли i-го элемента в металле и шлаке; g[i] - коэффициент активности в металле; y(i) - коэффициент активности в шлаке, рассчитываемый по модели коллективизированных электронов (МКЭ).

Однако в литературе не приводятся значения константы Ki/O . Чтобы использовать справочные данные по константам распределения, необходимо учесть переход от одного стандартного состояния (когда компонентом шлаковой фазы является элемент, как принято в МКЭ) к другому (когда компонент шлака - оксид).

Тогда получим выражение


, (9.58)

которое эквивалентно константе реакции (9.48):

(9.59)

Значения константы Ki можно вычислить по температурной зависимости констант, приводимых в литературе (они собраны в пункте 7 исходных данных задачи). Атомные доли выразим через массы элементов

; (9.60)

где – массы i-го элемента в металле и шлаке; – атомная масса, – суммы числа молей элементов, которые определяются следующим образом:

(9.61)

Учитывая выражение (9.60) и вводя обозначение

, (9.62)

получим окончательное выражение

(9.63)

Алгоритм решения

1. Определение массы каждого элемента, поступившего в систему с исходными материалами по формуле

, (9.64)

где – масса S-гo исходного материала; – процентное содержание i-го элемента в этом материале; – общее число исходных материалов.

Например, кислород находится в стали и в извести , тогда его масса равна

т.

2. Составление и решение системы уравнений. В систему входят: балансовое уравнение (1), уравнение константы (2), уравнение электронейтральности (3):

.

(9.65)

Система уравнений (9.65) решается методом итераций. Для этого необходимо задать начальные значения переменных. Вначале принимаем, что , , задаем значение отношения , эти величины уточняются в процессе итераций, поэтому их значения определяются лишь соображениями удобства. Задаемся также начальным значением окисленности металла, оно по условию задачи равно .

Вначале итерации ведутся на основе уточнения содержания кислорода в металле методом половинного деления. Указанное в условии задачи содержание кислорода принимается за верхнее значение интервала, за нижнее принимается . Если решение находится в заданном интервале, то последующее значение определяется по формуле

. (9.66)

Атомная доля кислорода в металле определяется первоначально по формуле

, (9.67)

затем, после определения состава металла, она рассчитывается по формуле (9.60).

Критерием завершения итераций служит баланс по кислороду. Если не соблюдается неравенство

,

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

3. Вычисление коэффициентов активности компонентов. С помощью МКЭ вычисляем g(i) по формуле (9.62), учитывая стандартное состояние
i-го компонента. Рассчитываем g[i] с помощью параметров взаимодействия Вагнера. При этом приняли, что раствор является разбавленным, и тогда коэффициенты активности имеют постоянные значения: g[Fe] = 1; g[Ca] = 1; ; ; .

Затем вновь обращаемся к пункту 2, подставив в систему (9.65) новые значения и приняв за нижнее значение , а за верхнее – то, которое получилось в предыдущей итерации. Критерием окончания вычислений являются установившиеся значения составов металла и шлака.

К решению прилагаются блок-схема алгоритма (см. рисунок), и результаты расчета (табл. 1 и 2). Описанный алгоритм можно использовать для моделирования любого сталеплавильного процесса или отдельных его периодов, когда известны составы и количества исходных материалов в начале периода.

Блок-схема алгоритма расчета задачи 89

Таблица 1

Результаты расчета задачи 89

Номер итерации по составу

Номер итерации по [O]

Массы элементов, т

Фаза

Fe Са Si Аl O
    99,8550 0,2143 0,0730 0,0920 0,2657 100,5000 3
1 1 98,3007 0 0 0 0,1073 98,4780 М
    1,5543 0,2143 0,0730 0,0920 0,6950 2,6285 Ш
  2 99,0718 0 0 0 0,0892 99,1610 М
    0,7832 0,2143 0,0730 0,0920 0,4747 1,6372 Ш
  3 99,4618 0 0 0 0,0448 99,5066 М
    0,3932 0,2143 0,0730 0,0920 0,3632 1,1357 Ш
  4 99,5580 0 0 0 0,0224 99,6805 М
    0,1970 0,2143 0,0730 0,0920 0,3072 0,8834 Ш
  5 99,7564 0 0 0 0,0112 99,7677 М
    0,0986 0,2143 0,0730 0,0920 0,2791 0,7569 Ш
  6 99,8057 0 0,0001 0 0,0056 99,8114 М
    0,0493 0,2143 0,0729 0,0920 0,2649 0,6935 Ш
  7 99,8303 0 0,0003 0 0,0028 99,8334 М
    0,0247 0,2143 0,0727 0,0920 0,2576 0,6613 Ш
  8 99,8180 0 0,0001 0 0,0042 99,8224 М
    0,0370 0,2143 0,0729 0,0920 0,2613 0,6775 Ш
  9 99,8118 0 0,0001 0 0,0049 99,8169 М
    0,0432 0,2143 0,0729 0,0920 0,2631 0,6855 Ш
  10 99,8149 0 0,0001 0 0,0046 99,8196 М
    0,0401 0,2143 0,0729 0,0920 0,2622 0,6815 Ш
  11 99,8165 0 0,0001 0 0,0044 99,8210 М
    0,0385 0,2143 0,0729 0,0920 0,2618 0,6895 Ш
  12 99,8172 0 0,0001 0 0,1043 99,8217 М
    0,0378 0,2143 0,0729 0,0920 0,2616 0,6785 Ш
  13 99,8176 0 0,0001 0 0,0043 99,8220 М
    0,0374 0,2143 0,0729 0,0920 0,2615 0,6780 Ш
2 1 99,7875 0 0,0040 0,0001 0,0043 99,7959 М
    0,0675 0,2143 0,0690 0,0919 0,2655 0,7081 Ш
  2 99,8213 0 0,0137 0,0004 0,0021 99,8375 М
    0,0337 0,2143 0,0593 0,0916 0,2446 0,6435 Ш
  3 99,8044 0 0,0068 0,0002 0,0032 99,8146 М
    0,0506 0,2142 0,0662 0,0918 0,2574 0,6803 Ш
  4 99,7960 0 0,0051 0,0002 0,0037 99,8050 М
    0,0590 0,2143 0,0679 0,0918 0,2638 0,6948 Ш
  5 99,9717 0 0,0045 0,0002 0,0040 99,8004 М
    0,0633 0,2143 0,0685 0,0918 0,2637 0,7016 Ш
  6 99,7939 0 0,0048 0,0002 0,0039 99,8027 М
    0,0611 0,2143 0,0682 0,0918 0,2628 0,6982 Ш
  7 99,7949 0 0,0050 0,0002 0,0038 99,8038 М
    0,0601 0,2143 0,0680 0,0918 0,2623 0,6965 Ш
  8 99,7954 0 0,0050 0,0002 0,0038 99,8844 М
    0,0596 0,2143 0,0680 0,0918 0,2620 0,6957 Ш

Окончание табл. 1

3 1 99,7926 0 0,0054 0,0002 0,0038 99,8019 М
    0,0624 0,2143 0,0676 0,0918 0,2624 0,6985 Ш
  2 99,8238 0 0,0177 0,0002 0,0019 99,8438 М
    0,0312 0,2143 0,0553 0,0918 0,2392 0,63l6 Ш
  3 99,8082 0 0,0091 0,0002 0,0028 99,8204 М
    0,0468 0,2143 0,0639 0,0918 0,2537 0,6704 Ш
  4 99,8004 0 0,0069 0,0002 0,0033 99,8108 М
    0,0546 0,2143 0,0661 0,0918 0,2584 0,6852 Ш
  5 99,7965 0 0,6061 0,0002 0,0035 99,8063 М
    0,0585 0,2143 0,0669 0,0918 0,2605 0,6920 Ш
  6 99,7945 0 0,0057 0,0002 0,0036 99,8041 М
    0,0605 0,2143 0,0673 0,0918 0,2615 0,6953 Ш
  7 99,7936 0 0,0056 0,0002 0,0037 99,8030 М
    0,0614 0,2143 0,0674 0,0918 0,2620 0,6969 Ш
4 1 99,7932 0 0,0056 0,0002 0,0037 99,8027 М
    0,0618 0,2143 0,0674 0,0918 0,2620 0,6973 Ш

Таблица 2

Состав металла (числитель) и шлака (знаменатель), % масс.

Номер итерации по составу
1 0,0043
2 0,0038
3 0,0037
4 0,0037

90. Расчет распределения элементов
переменной валентности между металлом
и шлаком с применением ЭВМ

Задача. Подшипниковую сталь обрабатывают в вакууме под синтетическим шлаком состава 40 % CaO, 30 % SiO2 , 30 % Аl2О3 . Рассчитать равновесное содержание железа в шлаке и его валентность при различном общем давлении в системе: 1; 0,1; 0,01; 0,001; 0,0001; 0,00001 атм, пренебрегая восстановлением компонентов шлака углеродом, содержащимся в стали, пользуясь классической моделью (законом действующих масс) и МКЭ.

Исходные данные. 1. Состав стали: 1,0 % С; 0,27 % Si; 0,30 % Mn; 0,027 % Р; 0,020 % S; 1,52 % Сr; 0,30 % Ni. 2. Параметры взаимодействия ; ; ; ; ; , . 3. Давление СО в системе равно общему давлению. 4. Термодинамические характеристики реакции окисления углерода:

; ;

5. Термодинамические характеристики реакции образования оксида железа (II):

; ;

6. Температурные зависимости параметров для определения валентности железа :

;

где – коэффициенты распределения кислорода между фазами1.

7. Температурная зависимость константы интегрирования:

.

8. Стандартные состояния: для углерода – 1 %-ный раствор в железе, для железа и оксида железа – чистые жидкие вещества.

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

(9.68)

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

Равновесие элементов описывается дифференциальным уравнением распределения их между фазами

, (9.69)

где – коэффициенты распределения i-го и j-го элементов между фазами, – их валентности в шлаке.

Коэффициент распределения i-го элемента между фазами определяется по формуле

, (9.70)

где , если шлак контактирует с газом и , если шлак контактирует с металлом.

Из выражения (9.70) получается дифференциальное уравнение:

(9.71)

При постоянных уравнение (9.71) переходит в выражение закона действующих масс. Действительно, после интегрирования (9.71) получим, что отношение и является постоянной величиной.

Зависимость от , полученная на основе МКЭ с помощью квантовой статистики Ферми – Дирака, имеет вид

. (9.72)

Выражение (9.72) определяет валентность как среднюю статистическую степень заселенности электронных уровней (3d и 4s) железа в шлаке. Константы соответствуют уровням 3d и 4s соответственно. На основании выражений (9.71) и (9.72), учитывая, что , имеем

. (9.73)

Алгоритм решения

1. Определение активности углерода в металле с помощью параметров взаимодействия по формулам

; (9.74)

где – коэффициент активности углерода в металле. Численное значение активности углерода для приведенных исходных условий (для контроля).

2. Определение давления кислорода

. (9.75)