Раздел 1
УДК 669.2
Мартусевич Е.А., Рыбенко И. А.
ДЕТЕРМИНИРОВАННАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДИНАМИКИ ИЗМЕНЕНИЯ ХИМИЧЕСКОГО СОСТАВА И ТЕМПЕРАТУРЫ РАСПЛАВА ПРИ ПОРЦИОННОМ СМЕШИВАНИИ ПЕРВИЧНОГО АЛЮМИНИЯ-СЫРЦА В ПРОЦЕССЕ ПОЛУЧЕНИЯ АЛЮМИНИЕВЫХ СПЛАВОВ В ЭЛЕКТРИЧЕСКИХ МИКСЕРАХ СОПРОТИВЛЕНИЯ ДЛЯ ОПРЕДЕЛЕНИЯ ОПТИМАЛЬНЫХ УСЛОВИЙ РЕАЛИЗАЦИИ ПРОЦЕССА
Аннотация. Формирование алюминиевого расплава в электрических миксерах сопротивления является сложным физико-химическим процессом, который характеризуется последовательным смешиванием алюминия-сырца, полученного методом электролиза и его взаимодействием с различными лигатурами и флюсами. На этой стадии происходит формирование заданного химического состава и качества алюминиевых сплавов посредством рафинирования и легирования расплава. Этот этап производства связан с наличием избыточного количества корректирующих воздействий в связи с многозадачностью и многофакторностью металлургического процесса, что приводит к увеличению времени приготовления расплава, снижению производительности электрического миксера сопротивления и повышению затрат на единицу готовой продукции. Следовательно, совершенствование технологии формирования алюминиевого расплава в электрических миксерах сопротивления и разработка оптимальных ресурсосберегающих режимов на основе математических моделей является актуальной научно-практической задачей современной металлургии. В связи с этим важное значение приобретает математическое моделирование металлургических процессов с применением программно-инструментальных систем, обеспечивающее высокое качество алюминиевых сплавов в результате прогнозирования конечных параметров химического состава с учетом свойств исходного сырья. Поэтому в рамках интенсивной цифровой трансформации металлургической отрасли актуальна разработка и промышленное применение математических моделей, позволяющих оптимизировать технологический процесс получения алюминиевых сплавов, снизив энергетические и материальные затраты, а также повысить производительность миксеров литейных отделений.
Ключевые слова: алюминиевый расплав, электрический миксер сопротивления, литейный ковш, математическая модель, задача оптимизации
Выделение объекта исследования
Объектом исследования является технологический процесс формирования алюминиевого расплава в миксере литейного отделения, который относится к сложным многопараметрическим объектам [1-4]. При выделении объекта исследования можно выделить следующие параметры: входных потоков, выходных потоков, технологические параметры и технико-экономические показатели (рис. 1).
При формировании алюминиевого расплава в миксере параметрами входного потока являются: К - количество литейных ковшей с алюминием-сырцом из электролизного цеха; Ь - количество легирующих добавок; Е - количество флюсовых добавок; тк - масса расплава алюминия-сырца из к-го литейного ковша (к =1 ^ К), кг; т1 - масса 1-й легирующей добавки (I =1 ^ Ь), кг; т/ - масса /-й флюсовой добавки для коррекции текущего химического состава расплава (/=1 ^ Е), кг; [Е]к - исходный химический состав алюминия-сырца в к-м литейном ковше, %; [Е],-г - химический состав 1-й лигатуры, %; [Е]/ - химический состав /-го флюса, %; т^ - исходная масса расплава в миксере перед началом выполнения технологической операции формирования расплава, кг; - исходная масса расплава алюминия-сырца в к-м литейном ковше, кг;
- температура исходного расплава алюминия-сырца в миксере, °С; ^ - температура алюминия-сырца в к-м литейном ковше, °С.
© Мартусевич Е.А., Рыбенко И.А., 2022
К параметрам выходного потока относятся: тМ -масса сформированного расплава в миксере с учетом легирующих и флюсовых добавок, кг; [Е],- - химический состав расплава в миксере согласно заданным требованиям установленного заказа, %; /М - температура расплава в миксере, °С. К технологическим показателям и технико-экономическим показателям относятся: тобщ - общее время формирования алюминиевого расплава в миксере, с; СМ - затраты на формирование алюминиевого расплава в миксере, руб; Цк - стоимость исходного алюминия-сырца в к-м литейном ковше, руб/кг; Ц - стоимость 1-й лигатуры, руб/кг; Ц - стоимость /-го флюса, руб/кг; Цэл - стоимость тарифа электроэнергии, руб./кВтч.
Разработка математической модели
Процесс получения алюминиевого расплава в миксере можно разделить на два последовательных этапа: формирование основы расплава путем смешивания алюминия-сырца из литейных ковшей и последующая корректировка расплава легирующими добавками и флюсами. Первый этап процесса формирования алюминиевого расплава заключается в смешивании исходного алюминия-сырца из к-х литейных ковшей и его дальнейшем перемешивании в миксере [5, 6]. Полученная масса расплава тМ в миксере определяется количеством перелитого алюминия-сырца из литейных ковшей:
МЕТАЛЛУРГИЯ ЧЕРНЫХ, ЦВЕТНЫХ И РЕДКИХ МЕТАЛЛОВ
ПАРАМЕТРЫ ПРОЦЕССА ФОРМИРОВАНИЯ АЛЮМИНИЕВОГО РАСПЛАВА В МИКСЕРЕ
ПАРАМЕТРЫ ВХОДНЫХ ПОТОКОВ
К - количество литейных ковшей с
алюминием-сырцом;
Ь - количество легирующих добавок;
Р - количество флюсовых добавок;
тк - масса алюминия-сырца, переливаемого в
миксер из к-го литейного ковша (к = 1 ^ К), кг;
т° - исходная масса алюминия-сырца в к-ом литейном
ковше, кг;
т1 - масса 1-й легирующей добавки (I =1 ^ Ь), кг; т/ - масса/-й флюсовой добавки (/"=1 ^ Р), кг; [Д]а - химический состав алюминия-сырца в к-м литейном ковше, %;
[Д]^ - химический состав 1-й лигатуры, %;
[Л],/ - химический состав /-го флюса, %;
Хк - температура расплава в к-м литейном ковше, °С;
тЦ - исходная масса расплава в миксере, кг;
¿и - температура исходного расплава в миксере, °С.
ПАРАМЕТРЫ ВЫХОДНЫХ ПОТОКОВ
тМ - масса сформированного расплава в миксере, кг;
[Д] - сформированный химический
состав расплава в миксере, %;
Хм - температура расплава в миксере, °С.
ТЕХНИКО-ЭКОНОМИЧЕСКИЕ ПОКАЗАТЕЛИ
Тобщ - общее время формирования алюминиевого расплава в миксере, с; См - затраты на формирование алюминиевого расплава в миксере, руб./т; Цк - стоимость 1кг алюминия-сырца в к-ом литейном ковше, руб.; Ц1 - стоимость 1кг 1-ой лигатуры, руб.; Ц/ - стоимость 1кг/-ого флюса, руб.; Цэл - стоимость электроэнергии, руб./кВт-ч.
Рис. 1. Параметры технологического объекта
¡\
тм = тИ + ^
тк,
(1)
к=1
Ш, =
'тМШ! 100
+
г
к=1
100
(1-Чг). (2)
расплавов алюминия-сырца из к-х литейных ковшей и их начальными значениями:
где К - количество литейных ковшей;
тМ - исходная масса расплава в миксере перед началом технологической операции формирования расплава, кг;
тк - масса алюминия-сырца, поступающего из к-го литейного ковша, кг.
Масса ,-го вещества в миксере рассчитывается с учетом массы и химического состава алюминия-сырца из литейных ковшей, а также угара элементов:
*'т/с[Д к4
с йть
= — = ~Утк
Ум =
йту[ с1т
= К К + т £ - тм) ,т 0 ) = т £ ;
171,
,(0) = ш°м,
(3)
где К - количество литейных ковшей;
- исходная масса расплава в миксере, кг;
тк - масса алюминия-сырца, поступающего из к-го литейного ковша, кг;
[Д], - химический состав алюминиевого расплава в миксере, оставшегося с предыдущего технологического этапа плавки, %;
[Д],к - химический состав алюминия-сырца в к-м литейном ковше, %;
- коэффициент угара ,-го элемента.
Скорости изменения масс химических веществ при формировании алюминиевого расплава в миксере можно представить в виде задачи Коши, определяемой системой обыкновенных дифференциальных уравнений относительно скоростей изменения масс
где г^ - скорость расхода массы (истечения) расплава из литейного ковша, кг/с;
- скорость прироста массы расплава в миксере,
кг/с;
- текущая масса расплава в к-м литейном ковше, кг;
- начальная масса расплава в к-м литейном ковше, кг;
- текущая масса расплава в миксере, кг;
- начальная масса расплава в миксере, кг;
- текущее время, с;
- постоянная времени запаздывания инерционного процесса, 1/с.
При этом то есть прирост массы форми-
руемого расплава в миксере равен расходу массы алюминия-сырца в к-м литейном ковше, что соответствует соотношению тм + тМ = т £-т Прирост масс ,-го вещества в миксере при переливании исходного алюминия-сырца из литейных ковшей с учетом начальных условий можно представить в следующем виде:
С1тщ йт
= кО
'м
-т£); т£ ( 0 )
х0
(4)
Раздел 1
Уменьшение массы ,-го вещества из к-х литейных ковшей при переливании расплава алюминия-сырца в миксер с учетом начальных условий можно представить уравнением:
^тк у у X
= -ут&; т£{0) = тк°.
{ ^ = У«0 + - т£); т-(0) = <-
^ = уК0 + т£'0 - т£); т*(0) = т*'0
йт.
= у(т£ 0 + т['° - тЩ'); т£(0) = тЩ'0
» и т и — т™ I ■ — --- и
йт
^ = УК10 + т^0 - г<); г<(0) = г<0 <"п£и = уК"0 + т£"0 - т£"); тМ"(0) = т,
йт йтпи
Си0 М
йт
гп0 м
<1т
Мп
м
йт
= у(г<0 + т*"0 - т£п); т^"(0) = т,
= у«"0 + тГ0 - тГ); тГ(0) = г<"0
йт\
щ- = у(т7 + т^0 - т7); т^(0) = <<
йт
йт
= у(т£ь + т£& - тЩ&); г<(0) = т,
РЬ
м
йт
С1тре
= -Г<е; <е(0) = тРке°
йтк йт йтк1 йт
йтк1 йт
гСкь
= -утЦ; 0) = т°к10
5(П
= -утпЦ; тп11(0) = тк1°
Т(п
= -ут$; т?(0) =
йтЦи
= -утсми; шски{0) = т£и°
йт (1т2кп йт
йт йтмк9 йт
йтркь йт
ат1п
= —ут^1; ш|п(0) = т2п°
= -ут"п; шГ(0) = ш™По = <«(0) =
= -утрмь; трк\0) = т<&°
(¿Т
{ ^ = ^ = + - Хм) Хк (0) = Х°к
0) — Х°
(8)
(5)
Математическая модель динамики изменения концентрации химических элементов в миксере и литейных ковшах соответственно определяется системой обыкновенных дифференциальных уравнений следующего вида:
„Ре
где - градиент изменения концентраций химических элементов в миксере;
- градиент изменения концентраций химических элементов в к-м ковше.
Последующий этап корректировки расплава с использованием флюсов и лигатур можно представить в виде уравнения
и г
тм = тЩ + А т ,, + ^ т 1 + А
к=1 1=1 /=1
Ш
(9)
где К - количество литейных ковшей; Ь - количество применяемых лигатур; Е - количество применяемых флюсов;
- исходная масса расплава в миксере, кг; тк - масса расплава, поступающего из к-го литейного ковша, кг;
т1 - масса 1-й лигатуры, добавляемой в алюминиевый расплав, кг;
т/ - масса /-го флюса, добавляемого в алюминиевый расплав, кг.
Тогда масса ,-го вещества в миксере будет определяться количеством алюминия-сырца из литейных ковшей и введенных в расплав лигатур и флюсов:
^ = у(т+ гпГ - т^); т°?(0) = (6)
, 100
100 А 100 А
/=1
1т/И;/ 1 00
' (1-4,). (10)
= т1п(0) = т?0 . (7)
В результате исходная задача Коши в векторной форме имеет следующий вид:
где К - количество литейных ковшей;
Ь - количество применяемых лигатур;
Е - количество применяемых флюсов;
- исходная масса алюминиевого расплава в миксере, кг;
тк - масса жидкого алюминия-сырца, поступающего из к-го ковша в миксер, кг;
т1 - масса 1-ой лигатуры, добавляемой в алюминиевый расплав в миксере, кг;
т/ - масса /-ого флюса, добавляемого в алюминиевый расплав в миксере, кг;
[Е],- - химический состав алюминиевого расплава в миксере с предыдущей плавки, %;
[Е],к - химический состав алюминия-сырца в к-м литейном ковше, %;
[Е] л - химический состав 1-й лигатуры, %;
[Е]/ - химический состав /-го флюса, %; - коэффициент угара ,-го элемента.
Технологический процесс формирования алюминиевого расплава реализуется в заданном температурном режиме 900-1000°С. Для расчета текущей температуры формируемого алюминиевого расплава в результате смешивания исходного алюминия-сырца с разной температурой используется уравнение
МЕТАЛЛУРГИЯ ЧЕРНЫХ, ЦВЕТНЫХ И РЕДКИХ МЕТАЛЛОВ
Ыы + п М ~ '
(И)
материального и теплового балансов [7, 8]. Целевая функция представляет собой суммарные затраты на получение единицы продукции:
где СМ - теплоемкость расплава в миксере, кДж/(кгград);
Хм - текущая температура алюминиевого расплава в миксере, °С;
^ - исходная температура алюминиевого расплава в миксере перед началом смешивания алюминия-сырца, °С;
4 - температура жидкого алюминия-сырца в к-м литейном ковше, °С;
At - приращение температуры в миксере за счет нагревательных элементов, °С.
Блок-схема математической модели процесса формирования алюминиевого расплава представлена на рис. 2.
Решение задачи оптимизации
Для определения оптимальных расходов алюминия-сырца, поступающего из литейных ковшей, и корректирующих добавок была поставлена задача оптимизации, которая заключается в минимизации затрат на формирование алюминиевого расплава в миксере с заданными свойствами при соблюдении всех технологических ограничений и выполнении условий
_ (££=1 т?сЦк+2[=1 £/=1 ш/Ц/+Цзл)
тм
-> тт, (12)
где Цк, Цг, Ц/ - стоимость алюминия-сырца, лигатур и флюсов соответственно, руб./т; Цэл - стоимость электроэнергии, руб./кВтч.
Решение задачи оптимизации заключается в определении минимума целевой функции при соблюдении следующих ограничений:
а) масса формируемого алюминиевого расплава в миксере должна не превышать полезную емкость миксера:
т„
< тм < гп„
(13)
б) масса алюминия-сырца, переливаемая в миксер из к-го литейного ковша, должна не превышать исходную массу алюминия-сырца в к-м литейном ковше:
О <тк< т°к,
(14)
К, Ь, Р, тк, т
о о
mk, т ■
I
т/ И» И, Д/ к '<
1 этап. Смешивание алюминия-сырца из литейных ковшей в миксере
1. Масса расплава тм в миксере:
ти = т
М + Ё
тк
2. Масса ,-го вещества в миксере:
т> = 1Ч00" + Ё 1С00 )'(1 ^
V к = 1 /
3. Скорости изменения масс веществ при смешивании алюминия-сырца из литейных ковшей в миксере:
йтк
"к = -ТТ = -Утк;
йт йт.]
V™ =■
•м
йт
= у(т£+т°к-тш);
тк( 0 ) = тк; тм ( 0 ) = .
4. Прирост массы ,-го вещества в миксере:
. х V х х
-¿¡г = 7{тш° + тк° - т&); т^( 0) = тш°.
5. Уменьшение массы ,-го вещества в литейных ковшах:
йт\
йт
ут^; тк ( 0 ) = тхк°.
2 этап. Корректировка алюминиевого расплава легирующими добавками и флюсами
6. Масса сформированного расплава тМ в миксере: тм = т£, + .
7. Масса ,-го вещества в миксере с учетом введенных в расплав легирующих добавок и флюсов:
Щ = I
тм И,
10 0 Ё 100 Ё
т1 [й]н
-т,[й]
100 100 100 100
\ к=1 I=1 f=1
8. Температура алюминиевого расплава в миксере:
( 1-4 ,) .
^м =
+ 1тк^к)См
шмСы
+ М .
Т
тм, Хм, [Д],
к=1
к
ь
ь
Рис. 2. Блок-схема математической модели процесса формирования расплава в миксере
Раздел 1
в) массы лигатур и флюсов должны не превышать массы имеющихся в наличии корректирующих материалов:
ml min ^ ml ^ ml тах>
mfmin < mf < rnfmax; (15)
г) формируемый химический состав алюминиевого расплава должен соответствовать заданному химическому составу установленной марки сплава:
[Я] i min < [я l < [я ] l тах>
(16)
д) температура формируемого расплава должна соответствовать заданному диапазону значений согласно требованиям технологической инструкции:
t min — ^М — ^ тах- 0-7)
При такой постановке задачи применение формальных методов оптимизации не всегда позволяет получить решение при соблюдении всех ограничений. При этом возможно несколько вариантов. В первом случае оптимальное решение по формированию расплава со всеми необходимыми характеристиками может быть получено при заданных начальных условиях, то есть при исходном наборе литейных ковшей с алюминием-сырцом с учетом лигатур и флюсов. Ко второму случаю относится ситуация, когда можно получить лишь неполное решение, например сформировать заданный химический состав алюминиевого расплава в миксере в соответствии с установленной маркой, но меньшей массы. К третьему случаю относится ситуация, когда оптимальное решение, исходя из заданных начальных условий, получить невозможно.
Однако в связи с тем, что в литейный цех непрерывно поступают новые литейные ковши с алюминием-сырцом, имеется возможность скорректировать начальные условия и добиться полного решения поставленной задачи, то есть получить алюминиевый расплав заданной массы, химического состава и температуры. Такую задачу можно решить методом динамического программирования, который применяется для определения оптимальных параметров объекта при постоянно изменяющихся начальных условиях. Для поиска оптимального решения в рамках метода динамического программирования на каждом этапе предлагается использовать метод линейного программирования, который обеспечивает решение задач, где целевая функция и ограничения имеют линейный вид.
Выводы
Для решения задачи определения оптимальных технологических режимов процесса формирования алюминиевого сплава в элекрическом миксере выделен объект исследования, определены параметры входных-выходных потоков и выбраны технико-экономические показатели. Разработана математическая модель технологического процесса формирования алюминиевого расплава в миксере, позволяющая рассчитывать пара-
метры смешивания алюминия-сырца из литейных ковшей в миксере с учетом присадок лигатур и флюсов. Математическая модель включает в себя уравнения материального и теплового балансов, динамические уравнения скоростей изменения масс и химического состава расплава при смешивании алюминия-сырца из литейных ковшей в миксере, что позволяет определить необходимое количество алюминия-сырца из литейных ковшей с учетом расходов лигатур и флюсов для формирования расплава с заданной массой, химическим составом и температурой. Поставлена и решена задача оптимизации, заключающаяся в минимизации затрат на получение необходимого количества алюминиевого расплава в миксере с заданным химическим составом и температурой при ограничениях на параметры входных-выходных потоков и технологические параметры. В связи с тем, что решение оптимизационной задачи невозможно получить формальными методами условной оптимизации из-за постоянно меняющихся начальных условий, для поиска решения предложен метод динамического программирования с использованием симплекс-метода линейного программирования.
Разработанная математическая модель предназначена для определения оптимального количества используемого алюминия-сырца из литейных ковшей с учетом введения лигатур и флюсов, обеспечивающего формирование алюминиевого расплава заданной массы, химического состава и температуры.
Работа выполнена в соответствии с грантом РФФИ («Конкурс на лучшие проекты фундаментальных научных исследований, выполняемых молодыми учеными, обучающимися в аспирантуре», договор № 19-37-90087\19), а также в рамках государственного задания (шифр темы 0809-2021-0013).
Список литературы
1. Агеев Н.Г. Моделирование процессов и объектов в металлургии. Екатеринбург: Изд-во УрФУ, 2016. 108 с.
2. Федоткин И.М. Математическое моделирование технологических процессов. М.: Ленанд, 2015. 416 с.
3. Лаврентьев Г.В., Лаврентьева Н.Б., Неудахина Н.А. Инновационные обучающие технологии в профессиональной подготовке специалистов. Барнаул: Изд-во АГТУ, 2009. 203 с.
4. Митришкин Ю.В. Линейные математические модели динамических систем с управлением. М.: Ленанд, 2019. 464 с.
5. Металловедение алюминия и его сплавов / А.И. Беляев, О.С. Бочвар, Н.Н. Буйнов [и др.]. М.: Металлургия, 1983. 280 с.
6. Уткин Н.И. Производство цветных металлов. М.: Изд-во «Интермет инжиниринг». 2004. 442 с.
7. Математическое моделирование действующего технологического процесса формирования алюминиевого расплава в миксере литейного отделения / С.Н. Калашников, Е.А. Мартусевич, Е.В.
МЕТАЛЛУРГИЯ ЧЕРНЫХ, ЦВЕТНЫХ И РЕДКИХ МЕТАЛЛОВ
Мартусевич, И.А. Рыбенко // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2021. № 56. С. 20-28.
8. Проектирование и реализация комплекса программ для моделирования технологического процесса формирования алюминиевых сплавов
Сведения об авторах
на основе объектно-ориентированного подхода / С.Н. Калашников, Е.А. Мартусевич, Е.В. Мартусевич, В.Н. Буинцев, И.А. Рыбенко // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2021. № 57. С. 110-117.
Мартусевич Ефим Александрович - старший преподаватель кафедры прикладных информационных технологий и программирования, Сибирский государственный индустриальный университет «СибГИУ», г. Новокузнецк, Россия. E-mail: [email protected]
Рыбенко Инна Анатольевна - доктор технических наук, доцент кафедры прикладных информационных технологий и программирования, Сибирский государственный индустриальный университет «СибГИУ», г. Новокузнецк, Россия. E-mail: [email protected]
INFORMATION ABOUT THE PAPER IN ENGLISH
A DETERMINISTIC MATHEMATICAL MODEL OF THE DYNAMICS OF CHANGES IN THE CHEMICAL COMPOSITION AND TEMPERATURE OF A MELTING DURING BATCH MIXING OF PRIMARY RAW ALUMINUM IN THE PROCESS OF PRODUCING ALUMINUM ALLOYS IN ELECTRIC RESISTANCE MIXERS TO DETERMINE THE OPTIMUM PRODUCTION CONDITIONS
Martusevich Efim A. - Senior Lecturer, Department of Applied Information Technologies and Programming, Siberian State Industrial University, Novokuzneck, Russia. E-mail: [email protected].
Rybenko Inna A. - Dr. Sci. (Eng.), Associate Professor, Head of the Department of Applied Information Technologies and Programming; Siberian State Industrial University, Novokuzneck, Russia. E-mail: [email protected].
Annotation. The formation of aluminum melt in electric resistance mixers is a complex physical and chemical process, which is characterized by sequential mixing of raw aluminum obtained by electrolysis and its interaction with various master alloys and fluxes. At this stage, the formation of a given chemical composition and quality of aluminum alloys occurs through refining and alloying of the melt. This stage of production is associated with the presence of an excessive number of corrective actions due to the multitasking and multifactorial nature of the metallurgical process, which leads to an increase in the melt preparation time, a decrease in the productivity of the electric resistance mixer and an increase in costs per unit of finished product. Therefore, the improvement of the technology of aluminum melt formation in electrical resistance mixers and the development of optimal resource-saving modes based on mathematical models is an urgent scientific and practical task of modern metallurgy. In this regard, mathematical modeling of metallurgical processes using software and instrumental systems is of great importance, which ensures high quality of aluminum alloys as a result of predicting the final parameters of the chemical composition, taking into account the properties of the feedstock. Therefore, within the framework of the intensive digital transformation of the metallurgical industry, it is relevant to develop and industrially apply mathematical models that allow optimizing the technological process for producing aluminum alloys, reducing energy and material costs, as well as increasing the productivity of foundry mixers.
Key words: aluminum melt, electric resistance mixer, casting ladle, mathematical model, optimization problem.
Ссылка на статью:
Мартусевич Е.А., Рыбенко И.А. Детерминированная математическая модель динамики изменения химического состава и температуры расплава при порционном смешивании первичного алюминия-сырца в процессе получения алюминиевых сплавов в электрических миксерах сопротивления для определения оптимальных условий реализации процесса // Теория и технология металлургического производства. 2022. №4(43). С. 9-14.
Martusevich E.A. Rybenko I.A. A deterministic mathematical model of the dynamics of changes in the chemical composition and temperature of a melting during batch mixing of primary raw aluminum in the process of producing luminum alloys in electric resistance mixers to determine the optimum production conditions. Teoria i tecnología metallurgiceskogoproizvodstva. [The theory and process engineering of metallurgical production]. 2022, vol. 43, no. 4, pp. 9-14.