Файл: Несенчук, А. П. Тепловые расчеты пламенных печей для нагрева и термообработки металла учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 115

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Частные производные в выбранной точке заменим через разно­ стные отношения (6.29) — (6.32):

д іп,

 

tn, т+ 1 tn, ■

(6.34)

дх

 

Ат

 

 

d2t„ , т

Іп—І,

т2tn

(6.35)

дх2

 

(Аху

 

 

Рис. 6.16. Определение производной

Рис. 6.17. Сетка прямоугольного типа,

функции f(x).

 

Тогда

дифференциальное

уравнение (6.33)

для узла А

(рис. 6.17) заменится соотношением

 

 

 

 

tn, лі+1 tn,

 

 

tn—1,

т 2 / „ ,

 

(6.36)

 

 

Ат

 

 

 

(Ах)2

 

или

 

 

 

 

 

 

 

2аДт

 

 

 

аАх

 

 

і

Г1

1

*

I

 

 

2аАт

(tn—l,m~ytn+i,m)-

(6.37)

7п, 7П+1

 

( Д х ) 2

J

Г" - т +

 

 

 

 

 

 

хУ

 

 

Выбирая различным

образом

соотношения между Ах и Ат,

формулу (6.37)

можно значительно упростить. Так,

приняв

Д т=

= (Дх)2/2а, получим

 

 

 

 

 

 

 

 

 

 

tn—i, ni + ^n+1. т

 

(6.38)

 

 

tn, т - Ц = =

 

Д

 

или

 

 

 

tn+l, m—l~\~tn—l, m—1

 

 

 

 

tn, тп

 

 

(6.39)

 

 

 

 

 

 

 

Последнее выражение, называемое уравнением Э. Шмидта, имеет большое практическое применение как при численном, так и при графическом решении задач нестационарной теплопровод­ ности.

106


Формула (6.39) позволяет найти температуру для всех узлов горизонтального ряда (например, ряда т) по известной температуре в узлах предшествующего ряда (т— 1). Так как начальными усло­ виями (при т = 0) задается распределение температуры по сечению тела (известна температура в узлах, находящихся на оси абсцисс), то можно последовательно найти температуру в узлах первого, вто­ рого и других рядов.

Рассмотренная сетка (рис. 6.17) удобна для численного интег­ рирования дифференциального уравнения (6.33) при граничных условиях первого рода (в любой момент времени известна темпера­ тура на поверхности тела), так как в этом случае граничные прямые х — 0 и х = Х принадлежат самой сетке.

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

 

^2экп—Хо Лі

Т^ЗЭКВ—Хз

7-1

 

(6.40)

 

 

Т Г '

 

 

 

 

Аз

 

 

где Х2 и Лз— действительная

толщина

соответственно

второго

Аі,

и третьего слоев, м;

 

 

 

 

 

 

и А3 — соответственно коэффициенты теплопроводности пер­

 

вого,

второго

и третьего

слоев,

ккал/м-ч-° С

 

(вт/м-° К).

 

 

 

 

 

 

 

 

В этом случае

расчетные

участки

(по толщине)

для

второго

и третьего слоев:

 

 

 

 

 

 

 

 

 

 

/ а \

(Ах)іАі ~і/ а2

 

 

 

 

 

(А-Ѵ) 2экв—

^

 

У

ßl

 

 

 

 

 

 

Л2

 

'

 

 

 

 

 

( Д х ) з э к в = (Д.ѵ)і

Л3

V

 

 

( 6- 41)

 

 

 

 

г

а1

 

 

При этом справедлива зависимость (6.39) для расчета температуры в каждом из слоев эквивалентной стенки.

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

tP, п= ^0(Д.т)2энв, т — ^(Дк)2экв. т-1+|7р-1, т - і ~ ^(Дк)2экв, т-і] *1, (6.42)

где

(Ах) 2ЭКВ

ki =

( Д х ) і + ( Д х ) 2экв

р— номер расчетного участка (A-c)i первого слоя, примыкаю­ щего ко второму слою многослойной стенки.

107


Температура на границе раздела второго и третьего слоев

ts, т — І0(Ах)3эка, т — і ц & х ) 3,л .й, m -l+[^-1, го-1—^(Дх)3э[(п, m -l] &2, (6 .4 3 )

где

(Дя)зэкв

W ) зэкп-)- (А'Ѵ) зэки

s — порядковый номер расчетного участка (ДлОгэіт второго слоя, примыкающего к третьему слою составной стенки.

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

iq, т ( t q —l, т - 1 t в) ^3, (6.44)

где

Ач/аз

ka

( А л ) з э и в - Ь ^ ч / а з

— коэффициент теплопроводности первого (основного) слоя составной стенки, ккал/м-ч-° С (0г/лг-°К);

аз — коэффициент теплоотдачи от наружной стенки к воздуху,

ккал/ м2■ч-° С (вт/м2-° К ); (Д.ѵ)зэкв — толщина расчетного слоя материала, прилегающего к на­

ружной поверхности, л і ;

q — порядковый номер расчетного участка (Дя)зэкв наружно­ го слоя составной стенки.

Если необходимо решить дифференциальное уравнение (6.33) при граничных условиях третьего рода (известны температура окру­ жающей среды и условия теплообмена между телом и окружающей

средой) и повысить точность определения температуры

на поверх­

ности, надо вводить дополнительные

узловые точки, не

лежащие

в области твердого тела.

 

 

Уравнение (6.33) решаем при граничных условиях третьего рода

. _ 0 — [Ф, % UI —К

Х=х = а [^л', т)—^c], (6.45)

где /с — температура окружающей среды.

Сетка строится так, как показано на рис. 6.18.

Таким образом, в рассмотрение вводятся точки, не лежащие в области твердого тела. Определяем температуру в дополнительных точках, отстоящих от рассматриваемой области на расстоянии Д.ѵ/2,

(to, m ГІ t h + l , m ) -

108


Производную ( - ^ г ) г.= ѵ>входящую в условие (6.45), заменяем

разностным соотношением

 

dt

\

=

tk+1, т tk, :

(6.46)

(

)

Ах

дх

1Х=х

 

 

Температуру на поверхности тела (в точке В) берем как среднее арифметическое температур в точках Л и С:

t (х, тАх) = tfi, т+ 4 + 1 , :

(6.47)

Рис. 6.18. Расчетная схема для гранич­ ных условии третьего рода.

Граничное условие (6.45) запишем так:

4 + 1 , т 4 , 1

= - —

( 4 , m + 4 + i , :

(6.48)

Ах

1

V

 

После некоторых преобразований (6.48) получим выражение для расчета температуры в дополнительном узле (точка С)

4+1, Vi— 4, ш+ (4 4, т) ki,

(6.49)

где

2Ах

4— 2К/а+Ах

Температура на поверхности плоской стенки (в узле В)

t (х, тАх) = 4 , тп+ (4 4, m )h,

(6.50)

где

Дх

5= 2%/а+Ах '

109

Температура в дополнительном узле D (рис. 6.18)

to, m = tit тЧ- (^c—ti' т)&4,

(6.51)

где titm — температура в узле F. Температура левой поверхности (узел Е)

t(0, тДт) =*!,,»+(/«.-*,, m)As.

(6.52)

Если рассматривается задача по расчету процесса нагрева (охлаждения) многослойной стенки, то температура полуслоев эле­ ментарных участков, прилегающих к границе раздела, выразит­ ся так:

для области первого слоя стенки

^(Р—0,5) (Дх)і, m 1о,5(Дд)2, (m-l)+[^(p-l,5) (Д.г-)і, (ni-1)—^0,5 (Дж)2, (m-1)] k', (6.53)

где

k, _

(Ах) 2

 

3 (Дх) і+ (Дх)г

для половины толщины элементарного участка второго слоя сложной стенки

to,5(Дх)2, т—^1,5 (Дэс)2, (т-і)+[^(р-0,5) (Д.ф, (т-1) А),5 (Д.ф, (т-1)] k", (6.54)

где

2(Дх)2

3 ( Д х ) 2+ ( Д х ) і '

Температура на границе раздела двух слоев составной стенки

tp, т = to,5(Дгс)2, m=to,5 (Дф, т+[^(р-0,5) (Д.ѵ)і, m~to,5 (Дж)2, т] &*, (6.55)

где

(Ах) 2

k*

(Дх) 1+ (Дх)2

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

^1“Мп+1 + ^2+- • •-Мп-1

tcpm

С (°К).

(6.56)

Выражение (6.56) при граничных условиях третьего рода при­ мет вид

ПО