№6
2008
ТЕХНОЛОГИЯ И ТЕХНОЛОГИЧЕСКИЕ МАШИНЫ
669.02/09
РАСЧЕТ ТЕПЛОВОГО ПРОФИЛЯ РАБОЧИХ СТЕНОК КРИСТАЛЛИЗАТОРА МАШИНЫ НЕПРЕРЫВНОГО
ЛИТЬЯ ЗАГОТОВОК
Канд.?пехн.наук,доу.С.В.Лукин} д-р.техн.наук,проф,В.Р.Аншелес, д-р,техн.наук,проф. П.Г Русаков,
аспТ.Н. Шесмаков, аспЛХ\ Поверим
Описана методика расчета оптимального с точки зрения теплопередачи профиля рабочей стенки кристаллизатора машины непрерывного литья заготовок при котором обеспечивается плотный контакт слитка с рабочей стенкой и высокая интенсивность теплообмена по всей высоте кристаллизатора. Опьчичие данной расчетной методики состоит в том, что процессы усадки, затвердевания и теплообмена в кристаллизаторе рассматриваются в тесной взаимосвязи друг с другом.
При разливке стали на машинах непрерывного литья заготовок (МНЛЗ) в кристаллизаторе формируется твердая оболочка, отделяющая жидкий металл от рабочей стенки кристаллизатора. Если рабочая стенка выполнена без учета усадки слитка, то отход оболочки от рабочей стенки происходит, когда скорость усадки превышает скорость пластической деформации оболочки под действием ферростатического давления жидкого металла [1].
Для обеспечения высокой интенсивности теплообмена слитка с кристаллизатором и исключения деформации его оболочки необходим плотный контакт слитка с рабочей стенкой, для чего профиль рабочих стенок должен быть выполнен с учетом усадки слитка.
При расчете усадки слитка и выборе профиля рабочей стенки процессы усадки, затвердевания и теплообмена слитка с кристаллизатором нужно рассматривать одновременно, поскольку они взаимно влияют друг на друга, однако в существующих методиках процессы затвердевания и усадки рассматривают в отрыве от процесса теплообмена.
58_Известия вузов. МАШИНОСТРОЕНИЕ__
№ 6 2008
В [2], [3] разработана математическая модель охлаждения и затвердевания слитка в кристаллизаторе, предполагающая, что по всей высоте кристаллизатора оболочка слитка в целом прижата к рабочей стенке ферростатическим давлением. При контакте горячей поверхности металла с относительно холодной рабочей поверхностью стенки поверхностный слой оболочки локально охлаждается и, в результате его термического сжатия, между металлом и стенкой образуется зазор, заполняемый смазкой. Зазор характеризуется эффективным коэффициентом теплопроводности Л39 который, при достаточной текучести смазки практически равен молекулярному коэффициенту теплопроводности смазки. При образовании зазора теплообмен между слитком и рабочей поверхностью резко ухудшается, оболочка разогревается и, в результате термического расширения металла, поверхность слитка прижимается к рабочей поверхности кристаллизатора. Рассмотренный процесс повторяется многократно, и протекает во всех точках рабочей поверхности.
В математической модели [2], [3] средняя стохастическая толщина зазора между поверхностью слитка и рабочей поверхностью кристаллизатора в зависимости от времени затвердевания т определяется выражением:
Чс(х) - ' ч ^
*/ ч V V п &(*>*) Г
V 0
• б/т
(1)
5т
где Р/ = Р/ (?) - температурный коэффициент линейного расширения твердой стали в зависимости от температуры; = (т) - толщина оболочки слитка по температуре солидус в зависимости от т; - температурное поле слитка в зависимости от расстоянии л* от его поверхности в момент т.
При расчете теплоотдачи от поверхности слитка к охлаждающей воде следует учитывать наличие гарнисажа из застывшей смазки на рабочей стенке вблизи уровня мениска жидкого металла. Пусть <50 > 0 - толщина гарнисажа. Термическое сопротивление от поверхности слитка к охлаждающей воде определяется выражениями:
Д(<М0Аом4АЛс + 1/Ч, 5(т)<50; (2,а)
Я(т) = б(т)Дз + 5С Дс + 1/<хв, 5(т) > 50, (2,6)
где Хш - коэффициент теплопроводности смазки; 8С - толщина рабочей стенки; Лс - коэффициент теплопроводности материала стенки; ав - эффективный коэффициент теплоотдачи от рабочей стенки к воде, учитывающий форму охлаждаемых каналов. Для щелевого кристаллизатора методика определения аъ приведена в [4]. Коэффициент теплоотдачи от поверхности слитка к охлаждающей воде равен:
а(т) = 1/Л(т). (3)
__Известия вузов. МАШИНОСТРОЕНИЕ_59
№ б 2008
Температурное поле оболочки слитка в кристаллизаторе с достаточной точностью описывается уравнением теплопроводности:
/ л \
, 0<х<5,0<т , (4)
эф " Р ' Л "" Л
ох дх
к* ах
Ч У
где СЭф~ эффективная теплоемкость стали; Хэф - эффективный коэффициент теплопроводности стали; р - плотность стали; В - половина толщины слитка. Выражения для расчета СЭф и Х^ приведены, например, в [5].
Начальное условие имеет вид:
*(х,0) = /0, 0 <х<В, (5)
где - температура жидкой стали, заливаемой в кристаллизатор. Граничные условия имеют вид:
&/дх\х=в = 0; дфх\х^ = - а(т)Д• (г(0, т) - ),
(б)
где Х- коэффициент теплопроводности стали; 4 - средняя температура охлаждающей воды.
Система уравнений (1) - (6) составляет математическую модель теплообмена слитка с кристаллизатором при условии, что обеспечивается плотный контакт слитка с рабочей стенкой. В [2] и [3] показана адекватность данной модели. На основе уравнений (1) - (6) можно численно рассчитать процесс затвердевания слитка в кристаллизаторе и определить толщину оболочки слитка по температуре ликвидус в зависимости от расстояния от мениска г. При стационарной скорости разливки V имеется простое соотношение между г иг: 2 = V- Т .
Рассмотрим усадку слитка в кристаллизаторе, при условии, что оболочка слитка не деформируется под действием ферростатического давления жидкого металла. Массовый расход металла через произвольное сечение кристаллизатора определяется так:
G = V■(ST■pr+Sж■pж), (7 ,а)
где V - скорость разливки; 1?т, — площади поперечного сечения слитка, занятые, соответственно, твердой и жидкой фазой; рт и рж - средние плотности твердой и жидкой фаз слитка.
Площадь сечения слитка равна: £ = + . В верхней части кристаллизатора, где имеется одна жидкая фаза, уравнение массового расхода запишем так:
где - площадь сечения кристаллизатора в верхней части.
Поскольку расход Э и скорость V являются одинаковыми для всех сечений кристаллизатора, то из выражений (7а) и (76) получим:
№6
1+^тЛ»
1 + (5т/5Ж)-(Р,/Р«)'
2005 (8)
Так как рж < рт, то из (8) следует, что 5 < , то есть сечение слитка в нижней части
кристаллизатора меньше, чем в верхней части.
Профиль рабочей стенки следует выбирать так, чтобы проходное сечение кристаллизатора на данном уровне равнялось сечению слитка.
Рассмотрим прямоугольный кристаллизатор, у которого площадь сечения на уровне мениска равна 2£0*2Й0, где Ь0 и Во - полуширина и полутолщина кристаллизатора на уровне мениска. Пусть 21 и 2В - ширина и толщина слитка на расстоянии г от мениска. Схема сечения кристаллизатора и слитка показана на рис. 1.
21
жидкая фаза
2 В
¡4--— /
л II ^......................
- —
| 1 — — —
у 11111111! Ш
V
<- \ , ►
2Во
21 о ^
твердая фаза
Рис. 1 Схема сечения слитка и кристаллизатора.
В силу симметрии достаточно рассмотреть четвертую часть сечения кристаллизатора и слитка. Площади сечения слитка и кристаллизатора на уровне мениска и на расстоянии г определяются выражениями:
^о = ¿о *5о; £ = (9)
Изменение полуширины А В и полутолщины ЛЬ слитка и кристаллизатора в результате усадки составит;
АВ = В0-В; = (10)
Можно считать, что А В намного меньше В0, а ЛЬ намного меньше Ь0. Если принять, что толщина твердой фазы слитка нарастает равномерно по периметру кристаллизатора, то Л В £ АЬ. С учетом этого, а также соотношений (9) и (10), изменение площади сечения слитка равно:
__Известия вузов, МАШИНОСТРОЕНИЕ_ 61
№ 6 2008
= = Ь0-В-1 = АВ-10 + М'В = АВ-(1 + В). (П)
Пусть £лО) - толщина твердой фазы по температуре ликвидус на отметке г. Площади сечения твердой и жидкой фаз слитка определяются так:
+ = + (12) На основе (8), (9), (11) и (12) получено выражение:
где рж - плотность жидкого металла; рТ(г) - средняя плотность твердой фазы на отметке г. Величина рТ(г) определяется так:
Ш
(14)
| V /
где р(г) - зависимость плотности твердой фазы от температуры.
Как следует из (13), величина усадки слитка по ширине и толщине пропорциональна толщине твердой фазы. Толщина твердой фазы по высоте кристаллизатора существенно зависит от скорости разливки, поэтому профиль рабочих стенок должен соответствовать максимальной скорости.
Для расчета усадки слитка по формулам (13) и (14), и выбора оптимального профиля рабочей стенки, необходимо рассчитать процесс теплообмена и затвердевания слитка в кристаллизаторе по математической модели, представленной уравнениями (1) - (6), в которой процессы теплообмена и затвердевания слитка в кристаллизаторе рассчитываются при условии плотного контакта между слитком и рабочей стенкой, что как раз и должно обеспечиваться при оптимальном профиле рабочих стенок.
В [6] представлены данные о проведенных в 1999 году эксплуатационных испытаниях на сортовой МНЛЗ. Для оптимизации теплопередачи здесь впервые был применен кристаллизатор с параболической конусностью. Испытания показали, что кристаллизаторы с параболическим профилем допускают чрезвычайно высокую скорость литья (порядка 5 м/мин и выше) при достаточно хорошем качестве сортовых заготовок. Рассчитаем процесс усадки квадратного слитка 100* 100 мм в сортовом кристаллизаторе для скоростей разливки 3,5 и 5 м/мин на основе уравнений (1) - (6), (13), (14). Теплофизические параметры зададим такими же, как для средне-углеродистой стали с содержанием углерода 0,45 %; начальная температура жидкой стали Г0 = 1550 °С; толщина медной рабочей стенки д0 = 11 мм; коэффициент теплопроводности меди Лс = 380 Вт/(м*°С); температура охлаждающей воды = 20 °С; эффективный коэффициент теплоотдачи аъ = 35000 Вт/(м2оС); в сортовом кристаллизаторе в качестве смазки применяется рапсовое масло с коэффициентом теплопроводности Хсм= 0,1 Вт/(м-°С); с учетом хорошей те-
62_Известия вузов. МАШИНОСТРОЕНИЕ__
№ 6 2008
кучести рапсового масла можно положить Хз = толщина гарнисажа принималась равной За - 0,015 мм. Расчет производился численно методом конечных разностей.
На рис. 2 показано уменьшение толщины слитка 2АВ в результате усадки в зависимости от координаты г для скоростей разливки 3,5 и 5 м/мин. Также на рис. 2 показана зависимость профиля рабочей стенки сортового кристаллизатора от координаты г, построенная по данным [б].
2ЛВА .....~................... ........."1 .....
мм ^■у
V =3,5 м/мщк ^^^^^^^
расчет
У профиль рабочей стенки
О ------1
О 0,2 0,4 0,6 0,8 ^»М
Рис.2 Усадка сортовой заготовки и профиль рабочей стенки
Сравнение расчетных данных, приведенных на рис. 2, и расчетных данных, приведенных в [6], показало, что они достаточно хорошо согласуются друг с другом. Чем больше скорость разливки, тем меньше усадка слитка в кристаллизаторе. Из рис. 2 следует, что при скорости 5 м/мин теоретическая усадка слитка больше, чем изменение профиля рабочей стенки. В действительности, под действием ферростатического давления оболочка слитка будет прижиматься к рабочей стенке. В случае, если теоретическая усадка слитка будет меньше, чем изменение профиля кристаллизатора, то будет происходить интенсивное истирание рабочей поверхности кристаллизатора и ухудшение теплопередачи в верхней части кристаллизатора.
ВЫВОДЫ
1. Профиль рабочих стенок кристаллизатора должен выбираться с учетом марки стали* теплофизических свойств смазки и скорости разливки.
2. Параболический профиль рабочих стенок кристаллизатора является достаточно близким к оптимальному.
-—
V =3,5 1/м V =5 М/^ЙЙГ"
расчет профиль абочей стенк I
1
№6
2008
СПИСОК ЛИТЕРАТУРЫ
1. Ефимов В.А., Эльдарханов A.C. Современные технологии разливки и кристаллизации сплавов. - М.: Машиностроение. 1998. - 360 с.
2. Лукин С.В., Шестаков H.H., Страшко Т.И. Математическая модель теплообмена сляба с рабочей стенкой кристаллизатора машины непрерывного литья заготовок // Изв. вузов. Черная металлургия. - 2007. - № 3. - С. 13-16.
3. Лукин С.В.; Шестаков Н.И., Страшко Т.И.. Зверев A.B. Охлаждение и затвердевание металла в кристаллизаторе машины непрерывного литья заготовок // Металлы, 2007. - № 3. - С. 20-26.
4. Шестаков Н.И., Калягин Ю.А., Манько О,В,, Аукин С.В. Теплообмен в рабочей стенке щелевого кристаллизатора машины непрерывного литья заготовок // Проблемы машиностроения и надежности машин, 2004. - №3. - С. 78 - 81.
5. Цаплин А.И. Теплофизика внешних воздействий при кристаллизации стальных слитков на машинах непрерывного литья. Екатеринбург: Изд-во УрО РАН. 1995. - 238 с.
6. Хорбах У., Коккентидт Й., Юнг В, Аитье сортовых заготовок с высокой скоростью через кристаллизатор параболического профиля // МРТ. - 1999. - С. 42 -51.
7.
621.825.52.001,2
МЕХАНИЗМ УПРОЧНЕНИЯ АУСТЕНИТНЫХ СТАЛЕЙ ПРИ ПЛАСТИЧЕСКОМ ДЕФОРМИРОВАНИИ
Докт. техн. наукпроф. В.А. Бубнов, канд. техн. паук, доц. С.Г. Костенко
Установлено влияние пластической деформации аустенитных сталей на изменение их механических характеристик прочности\ пластичности, ударной вязкости, а также магнитных свойств.На основании экспериментальных данных получена зависимость для определения предела выносливости аустенитных сталей, прошедших деформационное упрочнение.