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

СЕКЦИЯ 4

ЭНЕРГЕТИКА И ЭНЕРГЕТИЧЕСКИЕ ТЕХНИКА И ТЕХНОЛОГИИ

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

Кибардин Антон Николаевич

аспирант кафедры «Промышленная теплоэнергетика», ЧГУ,

г. Череповец Е-mail: anton [email protected]

Лукин Сергей Владимирович

канд.техн. наук, доцент, профессор ЧГУ, г. Череповец

Е-mail: [email protected]

THE CALCULATION OF HEAT TRANSFER IN THE SLIT MOULD WALL OF THE CONTINUOUS CASTING MACHINE

Anton Kibardin

post-graduate student of “Industrial heat-power engineering " department,

Cherepovets State University, Cherepovets

Sergey Lookin

candidate of technical sciences, professor of Cherepovets State University,

Cherepovets

АННОТАЦИЯ

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

ABSTRACT

In the paper is presented the engineering method of calculation of heat transfer in the slit mould wall of the continuous casting machine, allowing to take rather exactly account of the wall ribbing effect on the side of cooling water.

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

Key words: Heat transfer; mould; continuous casting machine.

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

Схема рабочей стенки щелевого кристаллизатора и основные геометрические размеры показаны на рис. 1, где приведены следующие обозначения: 1 — медная плита; 2 — стальная пластина; 3 — прямоугольный паз, заполненный охлаждающей водой; 4 — рабочая поверхность стенки кристаллизатора; 5 — толщина стенки от рабочей поверхности до основания паза; 2s — расстояние между пазами; h — глубина паза (длина ребра); 21 — ширина паза (канала); q — плотность теплового потока, поступающего к рабочей стенке от затвердевающего слитка.

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

Расчетная схема элемента рабочей стенки щелевого кристаллизатора показана на рис. 2. На рис. 2 показано, что поверхностью 1 рабочая стенка контактирует с отливаемым слитком, а

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

Рисунок 2. Расчетная схема элемента рабочей стенки.

Для расчета эффективного коэффициента теплоотдачи от рабочей стенки щелевого кристаллизатора к охлаждающей воде,

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

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

где а — обычный коэффициент теплоотдачи от поверхности к охлаждающей воде;

X — коэффициент теплопроводности материала стенки;

5, 5, I, И — геометрические размеры рабочей стенки, показанные на рис. 1.

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

г

3

4

5

\

/

5

*

6

ь

а.

(1а)

'эф

5 +1

(1б)

Введем избыточную температуру 3 = t -1 , где t — температура

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

д3+£|=а (2)

дх ду

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

д3/дх|х=0 = 0, 0 < у < s; Э3/Эу|у=0 = 0, 0 < х < h; д3ду\у=s = - (аА) •3Ц, 0 < х < h; 3L=h = 30, 0 < у <s,

где х, у — текущие координаты;

30 — избыточная температура у основания ребра.

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

3(х'у)=23 •£ •cos(у)•ch(• Iх) (3)

sin ц к

= - .

- к + вт - к • соэ - к

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

(-) = -, (4)

где Bi = а • я/Х — критерий Био.

Известно, что корни цк уравнения (4) представляют бесконечную

последовательность возрастающих чисел: ц1<ц2<...<ци<........ Поэтому

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

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

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

q. (К у) = 20. ■$„ ( ivf)sh ( Н.- ■ к) <«

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

1 s да D Is

qx(hI=-jqx(уIdy = 2'9o'Z h( Isin(h.I'sh(.-VsI

s 0 k=1 ch (k'h/sI

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

да D

q,P = qx(hI's = 2'3o'Zch( k.h/sIsin(h.I'sh(.-VsI (б)

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

qt0 = а'30 'l. (7)

Введем линейное термическое сопротивления ребра рабочей

стенки кристаллизатора:

Rlp = -^ = -f2'j h ( Dk Isin(h.I'sh. (В)

q,P І t=-ch (-k•h/s I f

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

q, = ( о + qp I2m = 3о '(V R,p + а •11'2m’ (9)

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

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

аэф =, (10)

эф P' 30

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

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

а эф =(1Яр + “•1 )/( 5 +1), (11)

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

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

ми = І------------*£*-------------, (12)

8 1 + 900 +12,7 ^ •Рг23 -1)

Яе \8 У ’

где Яе = w■dэ| Ув — критерий Рейнольдса;

Ми = а ■ dэ /Хв — критерий Нуссельта;

Рг — число Прандтля для воды;

Е — коэффициент гидравлического сопротивления;

w — средняя скорость движения воды в охлаждаемых каналах;

dэ — эквивалентный диаметр канала;

ув — коэффициент кинематической вязкости для воды;

Хв — коэффициент теплопроводности для воды.

Величины уъ, Хв и Рг выбираются по средней температуре охлаждающей воды 4. Коэффициент Е можно рассчитать по формуле [3]:

І = ( 0,79 •Іп (е/8))-2.

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

, 4 / 21 ■к

< = — = -т-г,

р I + к

где f=2l • 2к — площадь сечения канала; р=2 ^(21+2к) — периметр канала.

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

каналах — м>=6 м/с. При таких параметрах воды и размерах канала коэффициент теплоотдачи от стенки канала к воде, рассчитанный по формуле (12), равен а=20000 Вт/(м- К).

Определяем число Био: й = а-5/1 = 20000-0,013/385 = 0,675. Первые четыре корня характеристического уравнения (4) при этом Ы [2]: ц1=0,74; ц2=3,34; ц3=6,39; ц4=9,49. Расчет линейного

термического сопротивления ребра по формуле (8) дает значение Яр=0,00382 м К/Вт. Бесконечный ряд в (8) является быстросходящимся, и большая точность расчета достигается при использовании первых трех членов ряда. Расчет эффективного коэффициента теплоотдачи от стенки к воде по формуле (11) дает значение аэф=20110 Вт/(м2 - К), что ненамного больше, чем обычный коэффициент теплоотдачи а, несмотря на то, что на 1 м2 плоской поверхности рабочей стенки приходится 1,69 м2 оребренной поверхности, взаимодействующей с водой. Это объясняется достаточно низким коэффициентом эффективности оребрения Е, который, рассчитанный по формуле (1б), в данном примере равняется

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

Е = и = 2(“Т1'^ Ь(ВкЫ ^^). (13)

а-30-и ^ а-и 1 сЬ (цк-и/5)

Формула (13) дает для данного примера Е=0,52, т. е. существенно меньше, чем формула (1б). Для сравнения известной инженерной методики (формулы (1а) и (1б)) и разработанной инженерной методики (формулы (8) и (11)), был проведен численный расчет теплопередачи через рабочую стенку щелевого кристаллизатора при указанных условиях. Численный расчет позволяет учесть двухмерность температурного поля не только в ребре, но также в основании рабочей стенки, т. е. является более точным. Было получено, что численный расчет дает аэф = 19800 Вт/(м- К), известная инженерная методика - аэф=21840 Вт/(м К), разработанная инженерная методика - “эф=20110 Вт/(м К). Таким образом, “эф, полученный по известной инженерной методике, на 10,2 % превышает аэф, полученный численным методом. Значение “эф, полученное по разработанной инженерной методике, лишь на 1,5 % превышает аэф, полученное численно.

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

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

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

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

'раб. пов = 'осн. кан + Я - V1 = 'в + Я - (V1 + V“эф ) . (14б)

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

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

Выводы: Представлена инженерная методика расчета

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

охлаждающей воды.

Список литературы:

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

2. Лыков, А.В. Тепломассообмен: Справочник. — М.: Энергия, 1978. — 480 с.

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