Файл: Несенчук, А. П. Тепловые расчеты пламенных печей для нагрева и термообработки металла учеб. пособие.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) при граничных условиях третьего рода при мет вид
ПО