CC BY
43
УДК 681.32
РАСЧЕТ ТЕМПЕРАТУРНЫХ ПОЛЕЙ В ПОЛЫХ СЛИТКАХ ПРИ ЭЛЕКТРОШЛАКОВОМ ПЕРЕПЛАВЕ
А.Н. Суров, В.И. Потапов, М.С. Бугаев
Одним из перспективных направлений производства труб большого диаметра для особо ответственных трубопроводов является изготовление их из полых слитков путем раскатки. Полые слитки из нержавеющих и плоходеформируемых сталей производятся электрошлаковой технологией в специальных печах с двумя кристаллизаторами (внешним и внутренним). При раскатке таких слитков очень важно, чтобы стенка полого слитка была однородной структуры, не содержала дефектов в виде расслоений, наплывов и т.д. Поэтому при наплавлении слитка методом электрошлаково-го переплава (ЭШП) требуется строго выдерживать комплекс технологических параметров - ток переплава, глубину ванны, температуру охлаждающего теплоносителя и др. Поэтому для оперативного управления процессом наплавления слитка, оптимизации в данной работе рассматривается задача по разработке математической модели теплофизических процессов печи ЭШП с двумя кристаллизаторами. В работах [1] и [2] рассмотрена подобная задача, но с одним кристаллизатором.
Печь ЭШП с внутренним и внешним водоохлаждаемыми кристаллизаторами схематично изображены на рис. 1. В коаксиальном зазоре между кристаллизаторами наплавляется слиток. Расходуемые электроды (восемь штук) равномерно распределены в коаксиальном пространстве печи. Диаметры кристаллизаторов внешнего и внутреннего, соответственно равны 580 мм и 400 мм. Качество металла в слитках вторичного переплава во многом определяется направленной кристаллизацией, управление которой представляется возможным путем изменения теплофизических условий.
Рис. 1. Схема печи ЭШП для получения полого слитка в разрезе
Математическое моделирование теплофизических процессов и рассматривается в данной работе.
Печь со сплошным слитком была симметрична относительно оси А, поэтому было три дифференциальных уравнения в системе, описывающих связь между соответственно: вода, стенка и /-й слой, где i= 1, 2, 3,4; 1 - слой слитка твердого, 2 -слитка жидкого, 3 - шлака, 4 - электрода. Печь с полым слитком симметрична относительно оси В, поэтому система состоит из пяти дифференциальных уравнений энергии, описывающих связь между средами соответственно: вода внешнего слоя (в0, стенка^), /'-й слой, стенка(с2) и вода внутреннего слоя (в2).
При выводе уравнений приняты допущения: теплофизические параметры воды, стенки кристаллизатора, электрода, слитка не зависят от температуры, отсутствует диссипация энергии. Теплообмен аппарата ЭШП с внешней средой отсутствует. Ввиду цилиндрической формы печи ЭШП рассматривается осесимметричная задача.
Итак, с учетом принятых допущений, уравнения энергии для области Д = {х, V. 0 < х < г, 0 < ? < /*} запишем в виде:
<4,
д1
Эй.
ЗА
+ ув, ^ ^в,.с, (^с, ^В|) + ав
8%
дх2
= ^с1.в1К-^)+^с1,стК-^1)+«с дх
д2в(
с1 .
дх~
ьв„+и<л^=*ст_с1 (*ч - )+*ст,С2 (вС2 -в„)+
+а„
8% ' дх2
- + Я_
1^г
г дг
_8%
дг2
(1)
д(
дх*
061.
^в2,с2(#с2 в,г)+‘
3/ ”2 дх
/)2 = {*, ^ 2 Р,^> 0}:
дв* 8вл , х Г + °В1 Яг = (^с1 ~ ^)+ а1
д\
дх2
д2в.
дI
дА
д1
дв„
дх “1'"1' 4 дх
= К^(ввГвч) + К^{в„-вч) + а,
^сж,с,К-^) +
32й.
дх1
- + о.
дв„
д1 ™ дх
"*"-^-сж,с2 (^сж ~ вс2 ) ^ ас
двг
8% : дх2
- + я„
дг дг
(2)
___£2_ _ гг
с2’в2
К2-^2) + ^с2,«(0сж-^с2) + йс
80», дд. , , д\
± % Ох ~ (^2 ~ ^>2 ) + ЙВ;
8%
дх2
*2
5/
дх*
£>з = {х, /: р < х < сі, ґ > 0}:
за
ЗА
32А
■ + и„
а а*‘ *В,'С1^С| &2 ’
= *С],, К -*ч) + Ксиш(вш-вСі) + % ^;
+ иш ^ = ХС1>Ш (вш - ^) + КС2>Ш (#ш - <?С2) +
З/
+а„
д%
дх2
1 ЗА,, 32А Л
г дг дг
(?)
/
" К<2 ,в2 (*.2 - ^ ) + *с2,ш (*ш - »с2 ) + % ■
ЗА
±ив2^Г = ^в2,с2К-^в2) + «в
ЗЧ;
&2
д(
Зг 82 Зх
/)4 ={х, ґ. сі < х< І, х >0}:
два, ЗА , , д%
—- + Ч, ~^ = КВ с (е. - А ) + а ----------------------------Ь-
Я/ І с>Х ВЬС1 \ С1 В1/ В1 0^.2
дх
З/
ЗА
32А
^-иэ^=Кэ,С1(вС1-еэ) + Кэ,С2(вС2-еэ) +
СІ .
сЬ^
з2а
Зх2
дх
■ + а.
Л1 ЗА 32А А ---------2. +-----2
г дг дг
V
(4)
ч
ді
ЗА
326>
:^с2,в2(^в2 ^с2) + ^с2,з(^ 0с2)+ас2
ЗА
& --«2 &2 К*2.‘2\"С2 “*2І ' “В2 &2 ■
Начальные условия:
в'(х,г,0) = ір'(х,г), ф =[%]>%2,ірСі,рС2,^]7', / = 1,2,3,4; Г| <г < г2', 2(0 ,г) = г0, 20 = [0,0,0,2(г)]Г;
Д(0,г) = Д0, До=[0Д0,Д(г)]Г,
на подвижной границе фазового перехода слиток твердый - слиток жидкий
Рст<?2/ = Кг [<?ст (2 - 0.0], - Ясж |Чж (г + М], >
на подвижной границе шлак - электрод Рш9дг ^ЛЛ^р + Д-О’О^-^йДс + л + о.О]* ■ Условия на границах областей Ц : тг,0) = Тш;
в'(0,г,1) = (/(г. О, у = ,у/,2 ,Ц/Ч ,Ц/Ч ,^ст] ;
^"(/,г,0 = 0; ^[0ст(О,г,О]я = «ст,с, [^ст(0,г,0 - вч ],« = 1,2; Єсг(г-0,г,0 = Йсж(г + 0,г,/) = 7ш;
^еж(.р^,0 = тш\
Є'г(х,Гі,ґ) = [0], (Гг(х,г2,0 = [0], А’ = [0СТ ДЖ,0ШД]г;
Л[А"(х,г,0]л = «у [(^"(лг./г,/) - <9С], і = 1, 2;
К-*.2)
ЗА
+ <з».
в2
;/=12;
X 0 0 0' аст,с, 0 0 0 *
0 0 0 0 ^сжх. 0 0
Л= ; ц =
0 0 0 0 0 °тс, 0
0 0 0 л. 0 0 0
6|"(^,г,0 = 7ш;
^Ш(Р + А - О.г.О = бэ(Р + д + °^.0 = Гщ,; Є = /2^,;
х стр4 (Й3 + 273)4 -(ЙС + 273)4, г
2,4 /ул ^э~ес ’ 4 лсА
В уравнениях приняты обозначения: индексы вь в2, сь с2, ст, сж, ш, э у соответствующих параметров относятся к воде, стенкам, слитку твёрдой фазы, слитку жидкой фазы, шлаковой ванне, электроду; г = 1, _/ = ст; г = 2, у = сж; г = 3, ] = ш; і = 4, ] = э ;
<?„, > <?в2 . 0с, , 0с2 » 0ст . ^сж і вш, вэ ~ соответствующие температуры сред, °С;
р,, с„ Л,-, а, - соответственно плотность, кг/м3; удельная теплоёмкость, Дж/(кгград); теплопроводность, Вт/(мград); температуропроводность /-Й среды, м2/с;
аи, р„ 8, - соответственно коэффициент теплоотдачи между і-й и 7-й средой, Вт/(м2-град); периметр раздела, м; площадь поперечного сечения г-й среды, м2;
V,- - скорость движения г-й среды, м/с; х, г, і, г, А, р - соответственно текущие координаты по длине, м; радиусу аппарата, м; время, с; координата подвижной границы в слитке, м; координата подвижной границы в электроде и граница между электродом и жидкой ванной, м;
Тш, Тш - соответственно температуры плавления и шлаковой ванны, °С;
<7 - теплота кристаллизации, Дж/кг; щ, у/,, - известные распределения температур по соответствующим координатам;
тп - номер коаксиального слоя для слитка; а - коэффициент лучеиспускания. Термодинамические соотношения, замыкающие эту систему, такие, как зависимость коэффициента теплоотдачи от режима течения воды, коэффициента излучения и др., взяты из экспериментальных данных, полученных разными авторами.
Системы уравнений (1)—(4) решены численными методами с использованием метода конечных разностей.
На рис. 2 приведены некоторые результаты расчетов температурного поля в шлаковой ванне по полученной модели.
о
н-
гі
О.
*5*1
Н
«
О.
4>
С
2
ф
1800
1600
1400
1200
1000
800
600
400
200
0
у у N N
V N
I
\
\
А N
V N
0.4 0,43 0,46 0,49 0,509 0,535 0,565 0,58
радиус, м
Рис. 2. Распределение температуры в коаксиальном зазоре для шлаковой ванны при *1=0,4 и г2=0,58 для времени # = 2,5 ч по радиусу
74
Вестник ЮУрГУ, № 10, 2006
Суров А.Н., Потапов В.И.,
Бугаев М.С.____________________________________
Решение систем (1)—(4) позволяет определять температуру в процессе ЭШП в любой точке аппарата и в каждый момент времени.
Полученная в работе математическая модель тепловой работы печи ЭШП может быть основой разработки системы управления вторичным переплавом с целью экономии электроэнергии и повышения качества металла.
Литература
1. Потапов В.И., Торопов Е.В., Суров А.Н., Андрианова Р.А. Теппофизические процессы при электрошлаковом переплаве // Материалы Всероссийской научно-технической конференции «Те-
Расчет температурных полей в полых слитках ___________ при электрошлаковом переплаве
плофизика технологических процессов». - Рыбинск: Изд-во РГАТА, 2005.
2. Потапов В.И Математические модели теплофизических процессов в объектах многослойной структуры. - Челябинск: Изд-во ЮУрГУ, 2004. -270 с.
3. Потапов В.И., Рощин В.Е., Торопов Е.В. и др. Математическое моделирование теплофизических процессов при электрошлаковом переплаве трубчатых электродов // Труды второй всероссийской научной конференции «Математическое моделирование и краевые задачи». - Самара: Изд-во САМГТУ, 2005. -42. -208-211 с.
4. Йоши С., Митчелл А. Тепловые процессы при ЭШП// Электрошлаковый переплав. 1973. -138-180 с.
