Научная статья на тему 'Инженерная методика расчета теплопередачи в рабочей стенке щелевого кристаллизатора машины непрерывного литья заготовок'

Инженерная методика расчета теплопередачи в рабочей стенке щелевого кристаллизатора машины непрерывного литья заготовок Текст научной статьи по специальности «Механика и машиностроение»

CC BY
196
49
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ТЕПЛОПЕРЕДАЧА / КРИСТАЛЛИЗАТОР / МАШИНА НЕПРЕРЫВНОГО ЛИТЬЯ ЗАГОТОВОК / HEAT TRANSFER / MOLD / MACHINE OF THE CONTINUOUS CASTING BILLETS

Аннотация научной статьи по механике и машиностроению, автор научной работы — Лукин Сергей Владимирович, Кибардин Антон Николаевич, Калягин Юрий Александрович, Славов Владимир Ионович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по механике и машиностроению , автор научной работы — Лукин Сергей Владимирович, Кибардин Антон Николаевич, Калягин Юрий Александрович, Славов Владимир Ионович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Инженерная методика расчета теплопередачи в рабочей стенке щелевого кристаллизатора машины непрерывного литья заготовок»

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

Литература

1. Калягин, Ю.А. Исследование теплообмена в кристаллизаторе МНЛЗ с круглыми и щелевыми каналами / Ю.А. Калягин, Н.И. Шестаков, О.В. Манько, С.В. Лукин // Заготовительные производства в машиностроении. - 2004.

- № 12. - С. 29 - 31.

2. Лукин, С.В. Тепловые процессы при разливке стали на машинах непрерывного литья заготовок / С.В. Лукин. -Череповец, 2008.

3. Самарский, А.А. Вычислительная теплопередача / А.А. Самарский, П.Н. Вабищев. - М., 2003.

УДК 621.746. 27

С.В. Лукин, А.Н. Кибардин, Ю.А. Калягин, В.И. Славов

ИНЖЕНЕРНАЯ МЕТОДИКА РАСЧЕТА ТЕПЛОПЕРЕДАЧИ В РАБОЧЕЙ СТЕНКЕ ЩЕЛЕВОГО КРИСТАЛЛИЗАТОРА МАШИНЫ НЕПРЕРЫВНОГО ЛИТЬЯ ЗАГОТОВОК

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

Теплопередача, кристаллизатор, машина непрерывного литья заготовок.

The paper presents the engineering method of calculation of heat transfer in the working wall of slot mold of the machine of the continuous casting billets, allowing to take into account the ribbing effect of the wall from the side of the cooling water.

Heat transfer, mold, machine of the continuous casting billets.

Современные машины непрерывного литья заготовок часто оборудуют щелевыми кристаллизаторами, в которых охлаждающая вода циркулирует в пазах, расположенных между медной стенкой и чугунным корпусом. В настоящее время отсутствуют достаточно точные инженерные методики расчета теплопередачи в рабочей стенке щелевого кристаллизатора. Схема рабочей стенки щелевого кристаллизатора и основные геометрические размеры показаны на рис. 1.

Рис. 1. Схема рабочей стенки щелевого кристаллизатора

1 - медная плита; 2 - стальная пластина; 3 - прямоугольный паз, заполненный охлаждающей водой; 4 - рабочая поверхность стенки кристаллизатора; 5 - толщина стенки от рабочей поверхности до основания паза; 2$ - расстояние между пазами; Н - глубина паза (длина ребра); 21 - ширина паза (канала); д - плотность теплового потока, поступающего к рабочей стенке от затвердевающего слитка

Расчетная схема элемента рабочей стенки щелевого кристаллизатора показана на рис. 2. Можно увидеть, что поверхностью 1 рабочая стенка контактирует с отливаемым слитком, а поверхностью 3 и 4 - с охлаждающей водой. Поверхность 5 соприкасается с чугунным корпусом. На поверхностях 2 и 6 выполняется условие геометрической и тепловой симметрии, поэтому здесь теплообмен отсутствует (адиабатические поверхности). Размеры ребра: Нхз. Размеры основания стенки: 5х(^ + I).

Для расчета эффективного коэффициента теплоотдачи от рабочей стенки щелевого кристаллизатора к охлаждающей воде, учитывающего эффект оребрения, можно использовать выражение, полученное в [1]:

а (Е • Н +1)

$ +1

(1а)

где Е - коэффициент эффективности оребрения, определяемый выражением:

Е = Ж ((к/5 )^/ а • 5/ X )Д( ^5 )^/ а • 5/ X), (1б)

где а - обычный коэффициент теплоотдачи от поверхности к охлаждающей воде; X - коэффициент теплопроводности материала стенки; 5, 5, /, к - геометрические размеры рабочей стенки, показанные на рис. 1.

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

Введем избыточную температуру 3 = / - /в, где / - температура стенки; /в - температура охлаждающей воды, выбранная за начало отсчета. Температурное поле стенки опишется уравнением Лапласа:

ЭЧ ЭЧ = 0

дх2 + ду2 ~ .

(2)

Рабочую стенку можно представить состоящей из двух элементов: основания стенки толщиной 5 и ребра. Ввиду малой толщины задачу расчета теплообмена в основании стенки считаем одномерной. Расчет температурного поля и термического сопротивления здесь не вызывает затруднений. Остановимся на расчете теплообмена в ребре длиной к и толщиной 25. Граничные условия будут следующими:

Э3/Эх|х0 = 0, 0 < у < 5;

Э3/Эу| =0 = 0, 0 < х < к;

Э 3/Эу|у=5 = -(а/Х) • 3|у=5, 0 < х < к;

3х=к =J0, 0 < У < 5

где х, у - текущие координаты; 30 - избыточная температура у основания ребра.

С учетом граничных условий методом разделения переменных получено решение уравнения (2), описывающее температурное поле в ребре рабочей стенки щелевого кристаллизатора:

Д„

і сЬ(ц* • V5)

хсов| ц* • • ЛI ц* • X

Дк =■

віп ц *

Коэффициенты т* находятся из характеристического уравнения:

(ц ) = В1/ ц :

(4)

где В1 = а • $/X - критерий Био.

Известно, что корни т* уравнения (4) представляют бесконечную последовательность возрастающих чисел: т < т2 < ■■■< Дя < ...Поэтому ряд в (3)

является быстро сходящимся для точек с координатой х < Н.

Плотность теплового потока определяется законом Фурье:

д = X • |Э^/дЯ ,

где я - текущая координата. Из выражения (3) получим, что плотность теплового потока в ребре рабочей стенки кристаллизатора в направлении х в основании ребра равна:

дх (Н, у ) = 2 • X Ч • X Д' ц* А, х

'Л(ц* ■ Н/э)

х Соя I ц* ■ у I •I ц* ■ ~ I •

(5)

Для определения средней плотности теплового потока qх в основании ребра функцию (5) необходимо проинтегрировать на отрезке от 0 до 5 и разделить на длину этого отрезка 5:

1 $

дх (Н) = - [ дх (Н у) Лу =

е ^

Д*/$

*=1 сЬ (ц* • )

•я1п (ц* )• яЬ (ц* • ).

Тепловой поток, проходящий через половину основания ребра длиной 1 м (Вт/м), определится выражением:

др = дх (Н> $ =2 •х Ч • X

д.

'сЬ(ц* ■ )

х я1п(ц* )•(ц* • V*).

(6)

Тепловой поток, проходящий через половину основание паза, определяется выражением теплоотдачи к охлаждающей воде:

д,о = а • Чо • I.

(7)

ц* +я1п ц* •Соя ц*

Введем линейное термическое сопротивления ре-(3) бра рабочей стенки кристаллизатора:

я =3 =

я/р

%

аэф =

о

*

д„

*=1 сЬ (ц* -Н^)

’я1п (ц* )•(ц*^/$ Н . (8)

Линейная плотность теплового потока рабочей стенки с учетом (6), (7) и (8) определяется выражением:

щ =( Чю + % )•2т = 30 \11я1р + «•1 )• 2т (9)

где т - количество ребер в стенке.

Эффективный коэффициент теплоотдачи от рабочей стенки к охлаждающей воде определится выражением:

Ч,

а л. = ■ эф Р • Ч

(Ю)

где Р = 2т • (5 + /) - периметр рабочей стенки кристаллизатора.

С учетом (9) и (10) получим выражение для расчета аэф:

аэф =(!/ К,р + а 1)/($ +1),

(11)

где величина Я/р определяется выражением (8).

Обычный коэффициент теплоотдачи а определяется по известным эмпирическим зависимостям. Достаточно точной является формула Петухова [3], полученная при турбулентном течении капельных жидкостей и газов при 0,5 < Рг < 200 и 4000 < Яе <

<Ю'

Ми = — •

Яе • Рг

8 1 + ?оо +12,7./— •( Рг“ -1)

Яе \8 [ >

(12)

где Яе = w•dэ/ Ув - критерий Рейнольдса;

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Ми = а •dз/Xв - критерий Нуссельта; Рг - число

Прандтля для воды; X - коэффициент гидравлического сопротивления; w - средняя скорость движения воды в охлаждаемых каналах; - эквивалент-

ный диаметр канала; ув - коэффициент кинематической вязкости для воды; 1 - коэффициент теплопроводности для воды. Величины Ув, 1 и Рг выбираются по средней температуре охлаждающей воды гв. Коэффициент X можно рассчитать по формуле [3]:

X = ( 0,79 1п (Яе/8))-2.

Эквивалентный диаметр прямоугольного канала определяется выражением:

где /= 2/ • 2Н - площадь сечения канала; р = 2 • (2/ + + 2Н) - периметр канала.

Рассмотрим пример. На вертикальной МНЛЗ сталеплавильного производства ЧерМК ОАО «Северсталь» применяется щелевой кристаллизатор со следующими геометрическими размерами рабочей стенки: к = 0,025 м; / = 0,0045 м; 5 = 0,013 м; 5 = = 0,02 м. Стенка выполнена из меди с коэффициентом теплопроводности X = 385 Вт/(м • К). Средняя температура охлаждающей воды - /в = 30 °С. Скорость воды в каналах - w = 6 м/с. При таких параметрах воды и размерах канала коэффициент теплоотдачи от стенки канала к воде, рассчитанный по формуле (12), равен а = 20000 Вт/(м • К).

Определяем число Био: Ы = 5 X =

= 20000 0,013/385 = 0,675. Первые четыре корня характеристического уравнения (4) при этом Ы [3]: т = 0,74; |12 = 3,34; Дз = 6,39; Д4 = 9,49. Расчет линейного термического сопротивления ребра по формуле (8) дает значение Я/р = 0,00382 (м • К)/Вт. Бесконечный ряд в (8) является быстросходящимся, и большая точность расчета достигается при использовании первых трех членов ряда. Расчет эффективного коэффициента теплоотдачи от стенки к воде по формуле (11) дает значение аэф = 20110 Вт/(м2 • К), что ненамного больше, чем обычный коэффициент теплоотдачи а, несмотря на то, что на 1 м2 плоской поверхности рабочей стенки приходится 1,69 м2 ореб-ренной поверхности, взаимодействующей с водой. Это объясняется достаточно низким коэффициентом эффективности оребрения Е, который, рассчитанный по формуле (1б), в данном примере равняется 0,59. Более точно величину Е можно рассчитать в рамках разработанной инженерной методики по выражению:

Е = -

др

= 2^ (V Н )• X

д.

*=1 сЬ (ц *'Н1$)

х я1п (ц * )•яЬ (ц *-НІї).

(13)

Формула (13) дает для данного примера Е = о,52, т.е. существенно меньше, чем формула (1б). Для

сравнения известной инженерной методики (формулы (1а) и (1б)) и разработанной инженерной методики (формулы (8) и (11)), был проведен численный расчет теплопередачи через рабочую стенку щелевого кристаллизатора при указанных условиях. Численный расчет позволяет учесть двухмерность температурного поля не только в ребре, но также в основании рабочей стенки, т. е. является более точным. Было получено, что численный расчет дает аэф = 198оо Вт/(м • К), известная инженерная методика -аэф = 2184о Вт/(м • К), разработанная инженерная методика - аэф = 2о11о Вт/(м • К). Таким образом, аэф, полученный по известной инженерной методике, на 1 о,2 % превышает аэф, полученный численным методом. Значение аэф, полученное по разработанной инженерной методике, лишь на 1,5 % превышает аэф, полученное численно.

Таким образом, разработанная инженерная методика позволяет с хорошей точностью рассчитывать аэф на основе выражения (11). Инженерный расчет

теплопередачи в рабочей стенке щелевого кристаллизатора проводится также, как для плоской стенки толщиной 8. Если известна плотность теплового потока д, проходящего через рабочую стенку кристаллизатора, то температура поверхности основания канала, заполненного водой, определяется выражением:

'осн. кан = *в + Ч/аэф . (14а)

Температуру рабочей поверхности стенки, контактирующей со слябом, можно оценить по выражению:

'раб. пов = 'осн. кан + Ч'= 'в + Ч' ( 5/Х + 1/ «эф )• (14б)

Выражения (14) позволяют достаточно точно рассчитывать ?осн.кан и /раб.пов, если плотность теплового потока ч незначительно изменяется по высоте рабочей стенки кристаллизатора. В действительности величина ч быстро изменяется в верхней части кристаллизатора, что приводит к перетеканию теплоты вдоль рабочей стенки, и действительные температуры 'осн.кан и /раб.пов получаются несколько другими, чем дают формулы (14 а, б). В верхней части кристаллизатора действительные температуры /осн.кан и 'раб.пов будут меньше, чем рассчитанные по формулам

(14 а, б), а в нижней части кристаллизатора - несколько больше. Так как надежность работы кристаллизатора во многом зависит от температур /осн.кан и 'раб.пов на уровне мениска жидкого метала, где плотность теплового потока д является максимальной, то определение величин /осн.кан и /раб.пов по формулам (14 а, б) дает некоторый запас при расчетах надежной работы кристаллизатора.

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

Литература

1. Калягин, Ю.А. Исследование теплообмена в кристаллизаторе МНЛЗ с круглыми щелевыми каналами / Ю.А. Калягин, Н.И. Шестаков, О.В. Манько, С.В. Лукин // Заготовительные производства в машиностроении. - 2004.

- № 12. - С. 29 - 31.

2. Лыков, А.В. Тепломассообмен: Справочник / А.В. Лыков. - М., 1978.

3. Петухов, Б.С. Теплообмен при смешанной турбулентной конвекции / Б.С. Петухов, А.Ф. Поляков. - М., 1986.

УДК 621.746. 27

С.В. Лукин, А.Н. Кибардин, В.В. Мухин

УТИЛИЗАЦИЯ НИЗКОПОТЕНЦИАЛЬНОЙ ТЕПЛОТЫ В СИСТЕМЕ ОХЛАЖДЕНИЯ МАШИНЫ НЕПРЕРЫВНОГО ЛИТЬЯ ЗАГОТОВОК

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

Утилизация, низкопотенциальная теплота, машина непрерывного литья заготовок.

The article gives proof of the possibility of utilization of low potential heat in the cooling system of the machine of the continuous casting billets. The way of using the heat in the heating supply system of the metallurgical plant and the nearby town is suggested in the paper.

Utilization, low potential heat, machine of the continuous casting billets.

Машина непрерывного литья заготовок (МНЛЗ) является современным высокотемпературным металлургическим агрегатом, в котором температура жидкой стали, подаваемой в кристаллизатор МНЛЗ, достигает значений 1520^1550 °С. Доля жидкой стали, разливаемой на МНЛЗ, непрерывно растет во всем мире. В ведущих странах (США, Япония) эта доля приблизилась к 100 %. Наибольшее количество стали разливается на слябовых криволинейных МНЛЗ, получивших распространение благодаря относительно небольшим капитальным затратам и большой производительности.

Основная технологическая задача МНЛЗ - превратить жидкую сталь в твердую заготовку правиль-

ной формы (сляб), для чего нужно отвести от металла значительное количество теплоты. Приблизительно 50 % теплоты жидкой стали, подаваемой в МНЛЗ, отводится элементами оборудования МНЛЗ - кристаллизатором и роликами, охлаждаемыми водой, водой, подаваемой непосредственно на поверхность сляба, и паровоздушной смесью, образовавшейся при испарении части воды и удаляемой из бункера МНЛЗ [1], [2].

В настоящее время теплота, отведенная от металла в пределах МНЛЗ, практически никак не используется и теряется в окружающей среде, несмотря на значительное количество. Так, на пяти слябовых криволинейных МНЛЗ ЧерМК ОАО «Северсталь»

i Надоели баннеры? Вы всегда можете отключить рекламу.