Влияние технологических факторов на окисленность конвертерного шлака

ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2013р. Серія: Технічні науки Вип. 27

ISSN 2225-6733

УДК 621.385.6:669

© Сущенко А.В.1, Попов Е.С.2, Сидорчук Р.С.3, Гриценко А.С.4

ВЛИЯНИЕ ТЕХНОЛОГИЧЕСКИХ ФАКТОРОВ НА ОКИСЛЕННОСТЬ

КОНВЕРТЕРНОГО ШЛАКА

Выполнена проверка адекватности известных выражений для определения содержания оксидов железа в конвертерном шлаке после продувки плавки. Проведен статистический анализ и получены зависимости для прогноза окисленности шлака применительно к условиям работы 350-т конвертеров ПАО «МК «Азовсталъ».

Ключевые слова: конвертерная плавка, параметры шлака, содержание оксидов железа, технологические параметры, статистический анализ, регрессия.

Сущенко А.В., Попов Є.С., Сідорчук Р.С., Гриценко О.С. Вплив технологічних факторів на окисненість конвертерного шлаку. Виконана перевірка адекватності відомих виразів для визначення вмісту оксидів заліза у конвертерному шлаку після продування плавки. Проведений статичний аналіз та отримані залежності для прогнозування кисненості шлаку стосовно до умов работи 350-т конвертерів ПАТ «МК «Азовсталъ».

Ключові слова: конвертерна плавка, параметри шлаку, вміст оксидів заліза, технологічні параметри, статистичний аналіз, регресія.

A. V. Sushchenko, E.S. Popov, R.S. Sidorchuk, A.S. Gritsenko. Influence technological factors on converter slag oxidization. Examination of adequacy of known expressions for determination of iron oxide consistence in BOF slag after purging was completed. Statistic analysis was performed and dependences to forecast slag oxidization for the 350-t converters of PAS “MK “Azovsteel” operating conditions were obtained.

Keywords: converter melting, slag parameters, iron oxides consistence, technological parameters, statistical analysis, regression.

Постановка проблемы. Содержание оксидов железа в шлаке ^LFexOy) определяет его

окислительную способность, консистенцию и процесс шлакообразования в целом. От динамики изменения параметра ?,FexOy по ходу, а также его величины в конце продувки плавки, существенным образом зависят такие технико-экономические показатели производства стали, как: удельный расход металлошихты, степени усвоения шлакообразующих материалов, стойкость футеровки конвертера, качество полученного металла и др. Регулирование содержания оксидов железа в шлаке по ходу продувки, поддержание его на оптимальном уровне, а также прогнозирование конечной окисленности шлака, являются важнейшими задачами совершенствования технологии конвертирования чугуна и развития автоматизированных систем управления технологическими процессами в конвертере.

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

1 канд. техн. наук, доцент, ГВУЗ «Приазовский государственный технический университет», г. Мариуполь; sushchenko. andrei@gmail. com

2 директор по технологии и качеству ПАО «МК «Азовсталъ», входящий в ООО «Метинвест холдинг», г. Мариуполь; evseniy. popov@azovstal. com. ua

3 начальник технологического управления ПАО «МК «Азовсталъ», входящий в ООО «Метинвест холдинг», г. Мариуполь sidorchuk-rs@azovstal. com. ua

4 специалист первой категории ГВУЗ «Приазовский государственный технический университет», г. Мариуполь; Gritcenko. a. s@mail. ru

40

ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2013р. Серія: Технічні науки Вип. 27

ISSN 2225-6733

- из [1, 2] при [C]ме < 0.10 %:

(FeO )шл =

(

12 +

0,9

\

2 %

3, %

[C] ме У

где [C ] ме - содержание углерода в металле после продувки, %; - из [3]:

ИВ

(FeO L = 5,15

[C ] м

%

где В - основность шлака после продувки; - из [4] - в таблице 1:

Таблица 1

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

[C ] ме , % 0,03 0,04 0,05 0,06 0,07 0,08 0,09 0,1

(FeO )шл .% 19,00 16,95 15,50 14,50 13,85 13,75 13,65 13,58

- из [5, 6]:

(FeO )шл =

(

1,25 + 4 • В +

0,3

\

[C ] м

+10 ~6 • t 2

2 %

3, %

где Тме - температура металла после продувки, 0С.

Для проверки адекватности указанных выше зависимостей были выполнены сравнения расчетных (по [1-6]) значений (FeO)mn с фактическими данными по массиву плавок, проведенных на 350-т конвертерах ПАО «МК «Азовсталь» с ноября 2012 г. по август 2013 г. При этом в обработке использовали показатели только тех плавок, на которых выполнялись отборы проб для определения химических составов жидкого чугуна, металла и шлака после продувки, а также были выполнены замеры температур чугуна и металла. Были исключены плавки с возвратом металла в конвертер, промежуточным скачиванием шлака и проведенные на рельсовую сталь (с науглероживанием металла жидким чугуном в ковше). В итоге количество плавок, участвующее в обработке составило 2513. Результаты сравнения расчетных значений (FeO)mn по приведенным выше выражениям с промышленными данными приведены на рис. 1. Как видно, все указанные выше известные расчетные зависимости являются не пригодными для прогноза окисленности шлака после продувки плавок в реальных условиях работы 350-т конвертеров МК «Азовсталь».

Цель данной работы - изучить влияние производственно-технологических факторов на величину окисленности шлака после продувки плавки, и получить адекватные зависимости для ее определения в условиях конвертерного цеха МК «Азовсталь».

Изложение основного материала. Окисленность конвертерного шлака выражают различными способами: через суммарное содержание в нем железа (T.(Fe) или Fe^), концентрацию закиси железа (FeO ), суммарную концентрацию оксидов железа (FeO + Fe2O3), или приведенную суммарную их концентрацию, пересчитанную через эквивалентное количество FeO . В ККЦ МК «Азовсталь» в качестве показателя окисленности шлака используется содержание в нем FeO. Величина Fe^ в конвертерных шлаках определяется в экспресс-лаборатории на рентгеновском спектрометре ARL 9900 в соответствии с требованиями ДСТУ 3564-97. Для пересчета на содержание (FeO)mn используется эмпирическая зависимость, полученная путем статистической обработки результатов анализа проб шлака, выполненных химическими методами. Учитывая, что указанная зависимость периодически корректируется (в связи с изменением соотношения долей FeO и Fe2O3 в шлаках), для минимизации погрешности определения показателя окисленности шлака в качестве последнего была выбрана величина Fe^.

41

ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2013р. Серія: Технічні науки Вип. 27

ISSN 2225-6733

б)

Фактическое содержание FeO в шлаке после продувы:. 0 о

в)

30

<а

Факті рте ское содержание FeO в шлаке после продувки, %

г)

Фактическое содержание FeO в шпаке после продувки, %

Рис. 1 - Сравнение результатов расчета FeO в шлаке после продувки по методикам: а) - [1, 2]; б) - [3]; в) - [4]; г) - [5, 6] с данными промышленных плавок

42

ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2013р. Серія: Технічні науки Вип. 27

ISSN 2225-6733

Для определения влияния технологических параметров на величину Fe0^ был проведен статистический анализ (парная и множественная регрессии) показателей плавок по указанному выше массиву. Уравнения парной регрессии, имеющие коэффициент корреляции R > 0.32, приведены в таблице 2. При этом не значимыми оказались такие аргументы: Тме, [S]ме, [Р]ме,

Тчуг, [Мп]ч

[S ]ч

[ Р]ч

Мчуг • [Мп]чуг .10-

Мчуг • [S]чуг -10-

Мчуг • [Р]чуг -10-

Мчуг [Мп]чуг -10 /G0, Мчуг

[S]чуг • 103 /G0, Му

[Р]чуг' 103 /G0, В, ВВ, ^gO)^, Мизв -103 /G0

М_ • 103/ G, МБР • 103/ G, V

Е

. V Т

0 ’ О2 ’ прод ’

103 • М,,

ічуг

/ М

2

2

2

Таблица 2

Парные корреляционные зависимости величины Fe^ от некоторых технологических факторов

№ п/п Фактор Ед. изм. R Уравнение регрессии

1 [Si]4yz % 0.33 y = 13.499x2 - 30.163x + 34.706

2 Му • ту -10 -2 т 0.39 y = 1.4626x2 - 10.406x + 34.994

3 Мчуг ' Тчуг т*0С 0.40 у = - 4 • 10-10 x2 + 0.0001x + 21.777

4 Му • my -103 / G0 кг Si/ т МШ 0.40 у = 0.1973x2 - 3.8659x + 35.057

5 Мчуг  Тчуг / G0 т*0С/ т МШ 0.46 у = 7 -10-5x2 - 0.1969x + 145.14

6 [С ] ме % 0.36 у = 1637.3x2 - 300.03x + 30.686

7 [Мп ] ме 0/ % 0.49 у = 1135.5x2 - 280.64x + 34.14

8 О И о кг/т МШ 0.40 у = 0.0002x2 - 0.4397x + 212.79

9 Мприв у V чуг ж ) кг/т МШ 0.46 у = 0.0001x2 - 0.2715x + 151.14

10 Мприв У1 чуг ж кг/т МШ 0.53 у = 0.0003x2 - 0.4631x + 220.89

11 Мприв У11 чуг ж кг/т МШ 0.54 у = 0.0003x2 - 0.4933x + 232.52

12 Мприв Y чуг ж кг/т МШ 0.53 у = 0.0004x2 - 0.6229x + 284.87

13 Мприв У чуг ж кг/т МШ 0.50 у = 0.0001x2 - 0.2756x + 149.19

14 Мприв У1 чуг ж кг/т МШ 0.49 у = 0.0002x2 - 0.3321x + 173.13

При этом были приняты следующие обозначения: [Si] , [Мп]чуг, [S]чуг, [Р]ме - содержание в чугуне кремния, марганца, серы и фосфора соответственно, %; Тчуг - температура чугуна, 0С; Му - расход чугуна, т; [С] ме [Мп] ме [S] ме [Р] ме - содержание в металле после про? ? ?

дувки углерода, марганца, серы и фосфора соответственно, %; Тме - температура металла после 0~ тэ CaO „„ CaO + MgO

продувки, С; В =------ и ВВ = ---------5— - основности шлака после продувки соответствен-

SiO 2 SiO 2 + P2O5

но; (MgO)mji - содержание MgO в шлаке, %; Мизв, Мугля, М

- суммарные расходы извести

угля и брикетов раскисляющих на плавку, т; G0 - масса металлошихты (садка), т; Vq - суммарный расход технического кислорода на плавку, тыс. м ; Т прод - продолжительность продувки,

мин; Мме - масса жидкого металла после продувки, т; М^Ж - приведенный удельный расход жидкого чугуна на плавку, соответственно равный (варианты), кг/т МШ:

БР

43

ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2013р. Серія: Технічні науки Вип. 27

ISSN 2225-6733

(мприв Y

V чуг ж )

Мчуг • (61.9 + 0.88• Тчуг) чуг v чуг’ • 1000

чуг

G0 • (61.9 + 0.88• 1300) Мчуг •( (б19 + 0.88• Тчуг)+ 26970•

МпривVі

V чуг ж /

G0 •( (61.9 + 0.88• 1300)+ 26970•

[SQ 100 0.7 100

•1000

(мприв YIJ

V чуг ж )

—

(61.9 + 0.88 • Тчуг)+ 26970 • + 7000 •[—П]чуг

V чуг/ 100 100

—п]ч

чуг

0.7 0.2

G0 •! (61.9 + 0.88 1300) + 26970• — + 7000• — 0 1 v 7 100 100

•1000

Мприв Y'

X чуг ж /

Мприв)

X чуг ж /

—чуг •( (61.9 + 0.88• Тчуг)+ 26970• + 7000•

чуг \ чуг 100

[Мйічу

100

+147.753^ 4.25

f 0 7 0 2

G0 •! (61.9 + 0.88• 1300) + 26970• — + 7000• — +147.753^4.25 0 ^ 7 100 100

Мчуг + —21)• f 26970• + 7000• \М,і]чуг І + (61.9 + 0.88• Тчуг)• М,

х чуг чуг> у 100 100 J v чу ’

f 0 7 0 2

G0 •! (61.9 + 0.88• 1300)+ 26970• — + 7000• —

0 Г 7 100 100

—чш + MZ)• I 26970• + 7000•[ ^ +147.753^ 4.25 1 + (61.9 + 0.88• Г )• Муг

t VI X чуг чуг ’ у 100 100 J v чуг’ чуг

—Твж) =---------------Y-----------------------------------J-----------.---------1000

х чуг ж} / 07 02

G0-I (61.9 + 0.88• 1300) + 26970• —+ 7000• —+147.753^4.25 0 0 100 100

где (61.9 + 0.88 • Тчуг) - физическая теплота (энтальпия) жидкого чугуна, кДж/кг;

14775,3, 26970 и 7000 - приведенные тепловые эффекты реакций окисления углерода, кремния и марганца соответственно, кДж/кг элемента; 4,25 - среднее содержание углерода в чугуне, %; —Тув - суммарная масса твердого чугуна, т; 1300, 0.7 и 0.2 - средние (за 2012 г.) значения тем-

—]ч

1000

•1000

пературы, содержания кремния и марганца в жидком чугуне.

Наибольшее влияние на величину Fe^ оказывает комплексный параметр - приведенный удельный расход жидкого чугуна, определяемый по выражению (6), учитывающий физический и химический тепловые потенциалы жидкого чугуна. Это указывает на то, что в подавляющем большинстве случаев содержание Fe^ определяется тепловым балансом процесса, особенно в условиях «передува» плавок (окисление железа «в шлак» и «в дым» компенсирует недостаток приходной части теплового баланса).

В отличие от известных литературных данных, влияние [С ] ме и Тме на величину Fe0^ выражено в гораздо меньшей степени, причем с увеличением Тме имеет место тенденция снижения Fe^. Это, очевидно, объясняется следующим. Момент окончания продувки плавки машинистом дистрибутора зачастую определяется «по факелу» горящих отходящих газов над горловиной конвертера и показаниям системы газового анализа (т. е. фактически по содержанию углерода в металле); при этом, более высокая температура на повалке имеет место при большем тепловом потенциале шихтовки плавки. Следует также отметить, что величина [—n ] ме сильнее влияет на Fe^, чем [С ] ме . Наблюдается устойчивая тенденция увеличения окисленности шлака с ростом [S]чуг и [S]ме . В качестве примера, на рисунках 2, 3 приведены

графики зависимостей Fe^ от [С]ме и —прГж , определенному по (6).

При проведении множественной регрессии, незначимые параметры (аргументы) автоматически исключались из рассмотрения и не отражались в итоговом уравнении. Максимальное значение коэффициента корреляции (R=0.74) получено при совокупности следующих исходных параметров с соответствующими коэффициентами их влияния (рис. 4):

44

ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2013р. Серія: Технічні науки Вип. 27

ISSN 2225-6733

Feo4 = 56,5618 -2,87086 •[Si]4y! + 42.86821 • [S]чуг +118.0113 • [Р]чуг -40.0498 • [С]-34,9923 • [S]ж -- 75,2349 • [Р]ме - 2,01533 • (MgO)шл + 28,62067 • ВВ - 29,879 • В - 0,03505 • M

прив чуг. ж

Рис. 2 - Влияние содержания углерода в металле после продувки на содержание Fe0^ в шлаке

Пріведенный удельный расход жидкого чугуна с учетом его фіпігческой II XIВПРІЄ скоп {Si, Мп) теплоты, кг/т МШ

Рис. 3 - Зависимость Fe^ от приведенного удельного расхода жидкого чугуна, с учетом его физической теплоты и теплоты окисления, содержащихся в нем Si и Mn

Рис. 4 - Сравнение результатов расчета Fe^ в шлаке после продувки по зависимости (10) с данными промышленных плавок

45

ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2013р. Серія: Технічні науки Вип. 27

ISSN 2225-6733

Уравнение (10) является не удобными для прогнозирования величины Fe^ по известным (до начала плавки) исходным данным. Обычно задаются: параметры жидкого чугуна, одно значение основности шлака (как правило ВВ) и (MgO)mn, а также величина [С]ме. Параметры [Мп]ме, [S]*, [Р]ме и В являются не известными. При этом, если они будут рассчитаны не точно, то, не смотря на высокие коэффициенты корреляции, это приведет к значительной погрешности при прогнозировании величины Fe^. После исключения из рассмотрения в (10) указанных неизвестных параметров получена следующая статистическая зависимость (R=0.60):

Feotбщ = 65,4973 + 45,7757 • [5^ + 46.8525 • [Р]чуг - 97.4827 • [С]ме - 0,205 •

С целью увеличения точности прогнозирования величины Fe^ рассматриваемый массив плавок были разбит на части в соответствии с определенными диапазонами по [С]ме (0,02; 0,03; 0,04; 0,05; более 0,05 %). Однако это не привело к увеличению коэффициента корреляции.

Диапазоны изменения параметров, входящих в уравнение (11), а также показатели их влияния на изменение величины Fe^ приведены в таблице 3. Как видно, наибольшее влияние оказывает изменение теплового потенциала шихтовки плавки и степень ее «передува». Зависимость от величины [5]чуг связана с характерным «холодным» ходом продувки плавки при повышенном содержании серы в металле. Влияние [Р]чуг объясняется, очевидно, намерениями

машинистов дистрибутора дополнительно «умягчить» продувку или «передуть» металл с целью уменьшения величины [Р]ме при высоком содержание фосфора в исходном чугуне. Влияние основности ВВ связано с «разбавлением» концентрации оксидов железа в шлаке основными оксидами CaO и MgO (прямое влияние), а также с тем, что при более высоком тепловом потенциале шихтовки расход извести на плавку, как правило, увеличивают, при этом содержание Fe^ уменьшается, а содержание CaO возрастает (косвенное влияние). Зависимость от величины (MgO)mn связана с повышенной вязкостью шлака (особенно при больших значениях

(MgO)mjl), что приводит к усилению угара металла «в дым» и «в шлак», а также к увеличению количества корольков металла в нем.

Влияние аргументов зависимости (11) на изменение величины Fe^ в шлаке после продувки

Параметр Диапазон изменения аргумента Знак влияния Максимальная амплитуда колебания Fe^

МпРи 1 чуг. ж 685 - 854 - 9,73

[С ] ме 0,03 - 0,09 - 5,85

[5 ]чуг 0,006 - 0,08 + 3,39

Основность ВВ 1,7 - 4,9 - 3,07

(MgO)шл 2,35 - 14,25 + 2,44

[ Р]чуг 0,046 - 0,08 + 1,59

Таблица 3

Результаты сравнения расчетных значений Fe^ в шлаке после продувки по зависимости (11) с промышленными данными приведены на рис. 5. Отклонение расчетных значений от фактических связано как с точностью определения исходных параметров (масса, хим. состав, температура чугуна, хим. составы шлака и металла после продувки и др.), так и со следующими не учтенными объективными и субъективными факторами: качество металлоохладителей (различные виды лома, скрапа зашлакованного и т.д.), колебания хим. состава (в т.ч. величины ППП) в шлакообразующих и других сыпучих материалах, ход продувки плавки (режимы дутья, присадки сыпучих материалов и т.п.).

46

ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2013р. Серія: Технічні науки Вип. 27

ISSN 2225-6733

у = 0.3б24х+ 13.223

еС Фактическое содержание Feo6iit в шпаке после

продувки,0 о

Рис. 5 - Сравнение результатов расчета Feo64 в шлаке после продувки по зависимости (11) с фактическими данными плавок

Выводы

1. Показано, что известные и широко используемые в настоящее время, как на практике, так и в учебном процессе, зависимости для определения содержания оксидов железа в шлаке после продувки являются не пригодными для условий работы 350-т конвертеров МК «Азов-сталь».

2. На основании проведенного статистического двухфакторного анализа получено, что наибольшее влияние на величину Fe^ оказывает приведенный удельный расход жидкого чугуна (6), учитывающий его физический и химический тепловые потенциалы.

3. При проведении многофакторного регрессионного анализа получены уравнение (10) для прогнозирования окисленности конечного конвертерного шлака с коэффициентом корреляции R = 0,74, а также более удобная для практического использования зависимость (11), включающая только известные (до начала плавки) аргументы (R = 0.60).

4. Рассмотрены механизмы влияния значимых параметров (МПриеж, [C] ме, [«S']чуг, основность шлака ВВ, {MgO)mj, [Р] ) на содержание оксидов железа в шлаке после продувки.

5. Задачей дальнейших исследований является получение адекватных зависимостей, учитывающих дополнительное влияние на величину Fe^ видов и химических составов металлоохладителей, шлакообразующих и других сыпучих материалов, а также режимов продувки плавки и присадки сыпучих материалов.

Список использованных источников:

1. Бигеев А.М. Математическое описание и расчеты сталеплавильных процессов /

A. М. Бигеев. - М.: Металлургия, 1982. - 160 с.

2. Бигеев А.М., Металлургия стали. Теория и технология плавки стали / А.М. Бигеев,

B. А. Бигеев. - Магнитогорск: МГТУ, 2000. - 544 с.

3. Бойченко Б.М. Методичні вказівки до виконання розрахунків киснево-конвертерної плавки / Б.М. Бойченко, С.В. Бочка. - Дніпропетровськ: НМетАУ, 2001. - 25 с.

4. Димніч А.Х. Методичні вказівки до курсової роботи «Матеріальний і тепловий баланс конверторної плавки» / А.Х. Димніч [та ін.]. - Донецьк: ДонНТУ, 2005. - 44 с.

5. Харлашин П.С. Методические указания к выполнению курсовой работы «Расчет материального и теплового балансов кислородно-конвертерной плавки при переделе обыкновенных чугунов» / П.С. Харлашин, В.Я. Бакст. - Мариуполь: ПГТУ, 2005. - 20 с.

6. Бигееев A.M. Основы математического описания и расчеты кислородно-конвертерных процессов производства стали / А.М. Бигееев, Ю.А. Колесников. - М.: Металлургия, 1970. -232 с.

47

ВІСНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХНІЧНОГО УНІВЕРСИТЕТУ 2013р. Серія: Технічні науки Вип. 27

ISSN 2225-6733

Bibliography:

1. Bigeev A.M. Mathematic description and calculations of steelmaking processes / A.M. Bigeev. -

M.: Metallurgy, 1982. - 160 p.

2. Bigeev A.M., Metallurgy of steel. Theory and technology of steelmaking / A.M. Bigeev, V.A. Bigeev. - Magnitogorsk: MGTU, 2000. - 544 p.

3. Boichenko B.M. Methodic pointing for students of any specialties to making calculations of converter melting / B.M. Boichenko, S.V. Bochka. - Dnipropetrovsk: NMetAU, 2001. - 25 p.

4. Dimnich A.H. Methodic pointing for the term paper “Material and thermal balance of converter melting”

5. Harlashin P.S. Methodic pointing for the term paper “Material and thermal balance of oxygen-converter melting for redistribution of usual cast-irons” / P.S. Harlashin, V.Y. Bakst. - Mariupol: PSTU, 2005. - 20 p.

6. Bigeev A.M. The base of mathematic descriptions and calculations of BOF process / A.M. Bigeev,

Y.A. Kolesnikov. - M.: Metallurgy, 1970. - 232 p.

Рецензент: В.А. Маслов

д-р техн. наук, проф., ГВУЗ «ПГТУ»

Статья поступила 30.10.2013

48