CC BY
62
УДК 669.1.004.86
ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ЭШП*
Д.А. Пятыгин, М.А. Порсев, А.Н. Суров, И.В. Чуманов
PARAMETRICAL IDENTIFICATION OF MATHEMATICAL MODEL OF ESR
D.A. Pyateegin, M.A. Porsev, A.N. Surov, I.V. Chumanov
Повышение энергетической эффективности электрошлакового переплава (ЭШП) является актуальной проблемой. Предложена математическая модель теплофизических процессов, позволяющая определить температурные поля в любой точке шлаковой ванны. Результаты данной работы могут служить базой для разработки систем управления вторичным переплавом с целью экономии электроэнергии и повышения качества металла.
Ключевые слова: электрошлаковый переплав, математическая модель, теплофизические процессы, распределение температур, шлак.
Increasing of the energetic effectiveness of the electroslag remelting (ESR) is an urgent problem. The authors of the article suggest a mathematical model of the thermophysical processes which enables to define the temperature fields in any point of the slag bath. Results of this work can be the basis for the development of the control system of second remelting with a view to electrical energy economy as well as metal refinement.
Keywords: electroslag remelting, mathematical model, thermophysical processes, temperature distribution, slag.
Экспериментальное определение теплофизических параметров в плавильном пространстве электрошлаковой печи в условиях высоких температур, химически агрессивной среды и оптически непрозрачных сред - задача достаточно сложная. В таких условиях наиболее эффективным решением является математическое моделирование.
При выводе уравнений приняты допущения: теплофизические параметры воды, стенки кристаллизатора, электрода, слитка не зависят от температуры, отсутствует диссипация энергии. Теплообмен плавильного пространства с внешней средой отсутствует. Ввиду цилиндрической формы печи ЭШП рассматривается осесимметричная задача.
Математическая модель теплофизических процессов при электрошлаковом переплаве [1] представлена уравнениями для областей
= {х, х : 0 < х < z, 0 < т < хк };
В2 = {х, х : г<х<с, х > 0};
В3 = {х, х : с <х < d, х > 0};
В4 = {х, х : d<x< 1, х > 0}:
5ЄВ 50в /Q Q ч д\
^Г + ив^—_^вс(ес _0в) + ав'Т77';
5т 0х дх
50,
д20.
. =Хсв(0в-ес) + 5Сс,,(е/-ес) + «с _ 2 .
ся дх
50, 50, (а а ч 520,
7Г7 + и; - Х/с (6° ~®,) + а1 ~ГТ дх дх дх
(1)
+
+<х
150, 5 0
-----+------
г дг дг
где %, =
aiipi
1 і =1,2, 3,4-
'J PM
Систему уравнений (1) дополним начальными
условиями:
0' (дг, г, 0) = <р' (дг, г), ф' = [срвФсФ, ]Г, /■ = 1,2,3,4;
г(о,г) = г0,г0=[оАг(г)]Г
на подвижной границе фазового перехода слиток твердый - жидкий металл
РспЯ^Х ~ ^СТ [®СТ ~ ^5 ~
-^сж[бсж(^+о,т)]д.; (3)
* Работа проведена в рамках научной программы Федерального агентства по образованию - «Развитие научного потенциала высшей школы (2009-2010)» № 2.1.2/657 и при финансовой поддержке гранта Президента РФ для молодых докторов наук № МД-2610.2008.8.
Пятыгин ДА., Порсее М.А., Суров А.Н., Чуманов И.В.
Параметрическая идентификация математической модели ЭШП
на подвижной границе шлак - оплавляющийся электрод
рш^Ат = А,ш [0Ш (с + А - 0, т)]л -
~К [еэ(с+А + 0,т)]л (4)
условиями на границах областей Д :
0"0СТ (х, г, 0) = 0СЖ (х, г, 0) = Тш;
0'(О,/-,т) = у'(>,г);
= Ус Усг]Г; = 0;
^■ст [®ст ^ст,с [®ст 00 ] ?
0СТ (г - 0, г, х) = 0СЖ (г + 0, г, т) = Тпл;
^сж{с,г,\) = Тш-, 0'(х,0,т) = [О];
Л [0К (х, Я, = а [0К (х, Я, т) - 0С ];
ЯСт 0 0 ■ гу ст,с 0 0
А = 0 ^сж 0 ; а = 0 ^сж,с 0
0 0 I а 0 0 аш,с_
0”(х,/-,о) = х"'(л:^); х"(^'*) = [хв хс хш]Г;
вт(с!,г,т) = Тш; 0эл(с + Л-О,/-,т) =
= 0Эл(с + А + О,г,т) = 1’пл;
0щл =^ШЛ-
В уравнениях (1)-(5) приняты обозначения: индексы в, с, ст, сж, ш, э у соответствующих параметров относятся к воде, стенке, слитку твердой фазы, слитку жидкой фазы, шлаковой ванне, элек-троду; где 0В, 0С, 0СТ, 0СЖ, 0Ш, 0Э - соответст-вующие температуры сред; р7, с,, Х^, а1 - соответственно плотность, теплоемкость, теплопроводность, температуропроводность /-й среды; ау-,
Р1, ^ - соответственно коэффициент теплопередачи между /-й и 7-й средой, периметр раздела и площадь поперечного сечения /-й среды; Р, - скорость движения /-й среды; х, г, т, г, А, с - соответственно текущие координаты по длине, радиусу кристаллизатора, время, координата подвижной границы в слитке, координата подвижной границы в электроде и граница между электродом и жидкой ванной; ТП}1, Тш - соответственно температуры плавления и шлаковой ванны; д - теплота кристаллизации; ср,, ц/у, - известные распределе-
ния температур по соответствующим координатам; т - номер коаксиального слоя для слитка; а - коэффициент лучеиспускания; А г - величина шага по переменной г.
Численная реализация модели, представленная системой дифференциальных уравнений (1)—(5), даёт возможность получить внутренние переменные модели в виде численных значений температур 0В , 0С, 0СТ, 0СЖ, 0Ш • Дня получения, например, 0Ш ис-пользовались значения Хш=100 ккал/м-ч °С коэффициента теплопроводности шлака и аш = Хш/сшрш ,
где ст и рш соответственно теплоемкость и плотность шлака.
В данной работе разработана методика параметрической идентификации математической модели теплового режима шлака.
Задача параметрической идентификации математической модели определяется как задача отыскания таких числовых значений параметров 0Ш , аш, Хш, при которых расчётные значения выхода модели наилучшим образом согласовались бы с экспериментально полученными данными. Необходимость решения такой задачи обусловлена тем, что существуют трудности по измерению Хш [2]. Поэтому такой параметр лучше всего определять по экспериментальным данным. Как видно из расчётных и экспериментальных данных на участке идентификации их близость может оказаться неприемлемой.
С целью повышения точности настройки модели на реальный процесс целесообразно решить задачу параметрической идентификации, которая сформулирована следующим образом:
Исходной информацией в данном случае является модель и экспериментальные данные. Находим зависимость 0(гДш) как функцию, зависящую от радиуса г и коэффициента Хш :
1. При Хш = 25 ккал/м-ч °С, радиусе г = (0,05;
0,10; 0,15; 0,20; 0,25) и температуре шлака 0^ = = (1444; 1553; 1561; 1519; 1397) получим зависимость, которая изображена на рис. 1. Получим значение функции
г.^фг-еї)1-
= ^(33489 + 5041 +1521 + 3249 + 34596) * 125 .
е, °С
- - Расчетная
—■— Экспериментальная
Рис. 1. Распределение температуры в зоне шлаковой ванны
Серия «Металлургия», выпуск 13
29
На рис. 1 видно, что расхождение между расчетными и экспериментальными данными увеличилось.
2. При Хш = 40 ккал/м-ч-°С, радиусе г= (0,05;
0,1; 0,15; 0,2; 0,25) и температуре шлака 0^ = = (1539; 1592; 1584; 1542; 1431) значение функции
3. При Хш = 50 ккал/м-ч °С, радиусе г = (0,05; 0,1; 0,15; 0,2; 0,25) и температуре шлака 0щ = (1575; 1609; 1596; 1556; 1451) значение функции Р3 « 65 .
4. При Хш = 75 ккал/м ч*°С, радиусе г= (0,05; 0,1; 0,15; 0,2; 0,25) и температуре шлака 0^ = (1624; 1635; 1617; 1581; 1489) значение функции Р4 « 43.
5. При Хш = 100 ккал/м ч*°С, радиусе г = (0,05; 0,1; 0,15; 0,2; 0,25) и температуре шлака 0щ = (1647; 1648; 1629; 1597; 1515) значение функции ^ « 36.
6. При Хш =110 ккал/м-ч-°С, радиусе г = (0,05; 0,1; 0,15; 0,2; 0,25) и температуре шлака 0щ = (1653; 1652; 1632; 1601; 1523) значение функции Рь « 36.
7. При Хш =115 ккал/м*ч*°С, радиусе г = (0,05; 0,1; 0,15; 0,2; 0,25) и температуре шлака 0^ = (1656; 1653; 1634; 1603; 1526) значение функции « 36.
8. При Хш = 125 ккал/м*ч*°С, радиусе г = (0,05; 0,1; 0,15; 0,2; 0,25) и температуре шлака 0щ = (1660; 1656; 1636; 1606; 1532) значение функции « 37 .
9. При Хш = 150 ккал/мч°С, радиусе г = (0,05; 0,1; 0,15; 0,2; 0,25) и температуре шлака 0щ = (1668; 1660; 1641; 1613; 1545) значение функции « 38 .
10. При = 200 ккал/м-ч-°С, радиусе г = (0,05;
0,1; 0,15; 0,2; 0,25) и температуре шлака 0щ = (1678; 1666; 1646; 1620; 1562) значение функции Р10 « 39.
График зависимости ^(^Ш,0Ш) от Хш показан на рис. 2.
По графику (рис. 2) видно, что функция /Г(А,Ш,0Ш) имеет достаточно выраженный минимум, то есть функция является унимодальной. Есть основания утверждать, что выбранный метод решения параметрической идентификации работоспособен.
На рис. 3 изображены точками экспериментальные данные, а сплошной линией построена кривая по данным значениям параметров математической модели.
Среднеквадратическая погрешность настройки модели равна 36 °С, что составляет 2 % расхождения с экспериментальными данными.
Рис. 2. Зависимость функции Р от теплопроводности Лш
—►— Расчетная - ■ Экспериментальная
Рис. 3. Распределение температуры в зоне шлаковой ванны
Построенная модель, представленная уравнениями (1)-(5), позволяет разработать метод определения поля температур в шлаковой ванне, а также установить с помощью параметрической идентификации, коэффициент теплопроводности X шлака, в данном случае флюс АНФ-6. Для данной модели наиболее точное значение Хш = = 110 ккал/м-ч-°С.
Литература
1. Суров, А.Н. Расчет температурных полей в полых слитках при ЭШП / А.Н. Суров, В.И. Потапов, М. С. Бугаев // Вестник ЮУрГУ. Серия «Металлургия». - Вып. 7. - № 10. - 2006. - С. 73-75.
2. Тепловые процессы при электрошлаков ом переплаве / Б. И. Медовар, В.Л. Шевцов, Г.С. Марийский и др. - Киев: Наукова думка, 1978. - 304 с.
Поступила в редакцию 31 августа 2009 г.
