Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

tot_book

.pdf
Скачиваний:
40
Добавлен:
12.03.2015
Размер:
4.69 Mб
Скачать

141

Рис. 77.

(во всех одномерных стационарных задачах следовало бы использовать обозначения dT, dх и т. д.; сохраним знак частной производной лишь для однообразия и в силу традиции). Задача (2.20) допускает прямое интегрирование:

 

T

= C

(2.21)

 

 

 

x

1

 

 

 

 

T = C1 x + C2

(2.22)

 

где С1, С2 — постоянные интегрирования.

 

При х = 0 из равенства (2.22) следует, что Тw1 = С2; при х = δ

 

Тw2

= С1

+ Тw1, откуда C

=

Tw2 Tw1

. Получаем

 

 

 

 

 

 

1

 

 

δ

 

 

 

 

 

 

 

 

 

T =

Tw2 Tw1

x +T

(2.23)

 

 

 

 

 

 

 

δ

 

 

w1

 

 

 

 

 

 

 

 

 

Ранее были известны лишь температуры Тw1 и Тw2 уравнение (2.23) позволяет определить температуру в любом сечении пластины и показывает, что по координате х температура меняется линейно.

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

подставляя в уравнение (2.4) значение T x и уравнения (2.21):

 

 

 

 

142

 

 

 

 

 

q = λ

Tw1 Tw2

 

(2.24)

 

 

 

δ

 

 

 

 

 

 

 

 

Величину

δ

= R

 

называют

внутренним

термическим

λ

 

 

λ

 

 

 

 

сопротивлением4 ; ее размерность — м2 К/ Вт; очевидно, что

q = Tw1 Tw2 Rλ

Решения (2.23) и (2.24) можно обобщить для многослойной пластины (рис. 78).

 

 

 

Рис. 78.

4 Заметим, что

R

=

Tw1 Tw2

— отношение разности

 

 

λ

 

q

 

 

 

температур на двух изотермических поверхностях к плотности теплового потока. Для пластины q = idem, в более общем случае q может задаваться для одной из изотермических поверхностей.

143

Пусть в пластине n слоев с толщинами δ1, δ2,...,δn и теплопроводностями λ1, λ2,...,λn, соответственно; температуры границ слоев равны Тw1, Тw2,...,Тwn+1. Физически ясно, что в стационарном режиме плотность теплового потока q во всех сечениях пластины остается одинаковой:

q = Tw1 Tw2 ;

δ1 λ1

q= Tw2 Tw3 ;

δ2 λ2

.......................

q= Twn Twn+1

δn λn

Разрешим эти уравнения относительно разностей температур, а затем сложим почленно:

Tw1 Tw2 = qδ1 λ1 ;

Tw2 Tw3 = qδ2 λ2 ;

.......................

 

 

 

 

 

 

 

 

 

 

T

 

T

 

= qδn

;

 

 

 

 

 

 

 

 

wn

 

wn+1

 

 

 

λn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

T

 

 

 

δ

+

δ

2

+... +

δ

n

 

=

 

= q

1

 

 

 

 

w1

 

wn+1

 

 

 

λ1

 

 

 

λ2

 

 

λn

 

= q(Rλ1 + Rλ2 +... + Rλn ), поэтому

 

 

 

 

q =

 

 

Tw1 Twn+1

 

 

 

 

 

 

 

 

 

R

+ R

 

+...

+ R

 

 

 

 

 

 

 

(2.25)

 

 

λ1

λ2

 

 

λn

 

 

 

 

 

Величину

Rλ1 + Rλ2 +... + Rλn = RλΣ

 

 

назовем суммарным

внутренним термическим сопротивлением многослойной пластины (здесь видна аналогия с сопротивлением электрических цепей постоянного тока: суммарное сопротивление

144

последовательно соединенных проводников равно сумме составляющих). В осях Тх распределение температур в сечении пластины представляет ломаную линию; наклон любого ее звена тем больше, чем меньше теплопроводность соответствующего слоя. Усложним задачу (рис. 79): пусть обе поверхности пластины омывают жидкости с различными температурами (Тf1 > Тf2), а коэффициенты теплоотдачи на этих поверхностях равны соответственно α1 и α2. Рассмотрим, таким образом, задачу теплопроводности при граничных условиях III рода.

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

q

 

 

= α

(T

 

T

), внутри пластины

q

 

 

=

Tw1 Tw2

и на

 

 

 

 

 

 

 

 

 

 

 

 

 

x=0

1

 

f 1

w1

 

 

 

0<x<δ

 

δλ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

правой границе q x=δ = α2 (Tw1 Tf 2 ).

Рис. 79.

Выполнив те же преобразования, что и в предыдущем случае, получим

q =

 

Tf 1 Tf 2

 

=

Tf 1

Tf 2

,

(2.26)

1

 

+ δ + 1

 

R

+ R + R

 

 

 

 

 

 

 

 

α1

λ

α2

 

α1

 

λ α2

 

 

 

 

 

 

 

 

 

 

145

где

R

= 1

α1

, R

= 1

α2

внешние

термические

 

α1

 

 

α2

 

 

 

 

 

 

сопротивления, (м2 К)/Вт — величины,

обратные коэффициентам

теплоотдачи α1 и α2.

 

 

 

 

 

 

 

 

 

 

Величину

k =

 

1

 

 

,

которая

характеризует

 

1 α +δ λ +1 α

2

 

 

 

 

1

 

 

 

 

 

 

интенсивность передачи теплоты, называют. Из формулы (2.26) следует, что значение k численно равно отношению плотности теплового потока, передаваемого через пластину от одной жидкости к другой, к разности между температурами жидкостей

(теплоносителей)5.

Если величина α характеризует только конвективный теплообмен на поверхности пластины, то коэффициент

теплопередачи k является характеристикой всей теплопередающей системы, поскольку кроме коэффициентов теплоотдачи α1 и α2 она включает параметры δ и λ!

Величина 1k = Rk — ясно, что Rk = Rα1 + Rλ + Rα2 . Распределение температур в многослойной пластине при

граничных условиях III рода в пределах каждого слоя остается линейным, а в жидкости или газе вблизи границ пластины приобретает нелинейный характер (рис. 80); плотность теплового потока

q =

 

 

Tf 1 Tf 2

 

 

=

Tf 1 Tf 2

= k(T

f 1

T

f 2

).

(2.27)

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

Rk

 

 

1

 

+ δ

 

+

1

 

 

 

 

 

 

α1

λ

α2

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 Подробнее о теплоносителях см. разд. 2.12.

146

Рис. 80.

Нередко требуется определить температуру в произвольном сечении пластины 0 < х < δ. Запишем значения перепадов температур, используя прежнее условие q = idem:

Tf 1 Tw1 = q α1 ;

Tw1 Tw2 = qδ1 λ1 ;

.......................

Twn Twn+1 = qδn λn ;

Twn+1 Tf 2 = q αn+1 .

Температуры на границах слоев

Tw1 =Tf 1 q 1α1

Tw2 =Tf 1 q 1α1 +δ1 λ1

.....................................................

Twn+1 =Tf 1 q 1α1 +δ1 λ1 +...+δn λn ;

Tfn+1 =Tf 1 q 1α1 +δ1 λ1 +... +δn λn + 1αn+1 .

Пусть сечение х находится в j-м слое стенки. Тогда расстояние от этого сечения до границы с температурой Тwj будет

равно xj = х – (δ1 + ... + δj–1), а температура

147

T (x)=T

 

 

1

 

+

δ

+...+

δj1

+

x j

 

(2.28)

f 1

q

α1

1

 

 

.

 

 

 

 

 

λ1

 

λj1

 

λj

 

Рассуждения остаются в силе, если расчет начинать не с Тf1, а с Тf2 перед q поставить знак “+”, а все члены суммировать в обратном порядке.

Теплопроводность цилиндрической стенки при граничных условиях I рода рассчитать чуть сложнее, поскольку в этом случае плотность потока q зависит от текущего радиуса r (рис. 81).21.10

Рис. 81.

Обычно цилиндрические стенки достаточно “длинные” (трубы, каналы и т. п.); их поле температуры вполне точно описывают одномерным уравнением Фурье, т. е. считают, что температура вдоль оси цилиндра постоянна.

Краевая задача Фурье для осесимметричной цилиндрической стенки без внутренних источников имеет вид

2T

+

1 T

= 0;

r2

 

 

r r

 

(2.29)

Tr=r1 =Tw1;

 

Tr=r2 =Tw2.

Введем новую переменную U = Tr , чтобы понизить порядок производных, тогда

 

 

 

 

 

148

 

 

 

2T

=

U

;

1 T

= U

;

U

+U

= 0;

r2

 

r

 

r r

r

 

r

r

 

разделим переменные и проинтегрируем; получим последовательно

lnU +ln r = C1,

 

∂ =

r

,

(2.30)

T

C1 r

 

T = C1 ln r +C2.

Постоянные С1 и С2 найдем из граничных условий задачи

(2.29):

при

r = r1

 

T =Tw1;

Tw1 = C1 ln r1 +C2;

 

 

 

при

r = r2

 

 

T =Tw2;

Tw2 = C1 ln r2 +C2;

 

 

 

C =

Tw1 Tw2

 

; C

 

=T

(T

T

)

ln r1

 

 

.

 

 

 

ln(r r

)

1

ln(r r

)

 

 

2

 

w1

w1

w2

 

 

 

1

2

 

 

 

 

 

 

 

 

 

1

2

 

 

После подстановки получим значение температуры Т на произвольном радиусе r:

T (r)=T

(T

T

)

ln r r1

 

.

(2.31)

ln(r r

)

w1

w1

w2

 

 

 

 

 

 

 

1

2

 

 

 

Уравнение (2.31) описывает убывающую логарифмическую функцию, что физически вполне объяснимо: с увеличением значения r плотность теплового потока снижается, хотя тепловой поток, подводимый на поверхности r =r1, и отводимый с

поверхности r = r2, остается неизменно равным Q = −λ Tr 2πrl

(здесь l — длина расчетного участка).

Поскольку в соответствии с равенством (2.30)

T

=

C1

=

(Tw1

Tw2 )

,

Q =

2πλl(Tw1 Tw2 )

 

, Вт.

r

r

 

 

r ln(r r )

 

 

r ln(r r

)

 

 

 

 

 

 

 

 

1

2

 

 

 

1

2

 

 

Для цилиндрической стенки плотность теплового потока удобнее рассматривать либо на поверхности 2πr1l, либо на поверхности 2πr2l:

 

 

 

 

 

149

 

 

 

 

q1 =

 

Q

=

(Tw1 Tw2 )

,

Вт м2 ,

(2.32)

 

2πr1l

 

 

 

 

(r1 λ)ln(r2 r1 )

 

 

 

 

q2 =

 

Q

 

=

(Tw1 Tw2 )

 

,

Вт м2 ,

(2.33)

2πr2l

 

 

 

 

 

(r2 λ)ln(r2 r1 )

 

 

 

Можно также рассчитать тепловой поток, протекающий через единицу длины цилиндра,

ql =

Q

=

π(Tw1 Tw2 )

, Вт м,

(2.34)

l

 

 

(1 2λ)ln(r2 r1 )

 

 

эту величину назовем линейной плотностью теплового потока.

Легко заметить, что

ql = 2πr1q1 = 2πr2q2;

(2.35)

иначе говоря, линейная плотность теплового потока ql , не зависит от текущего радиуса r. Для n-слойной цилиндрической стенки все рассуждения проведем так же, как для многослойной пластины; при этом линейная плотность теплового потока окажется равной

ql =

π(Tw1 Tw2 )

, Вт м.

 

n

 

1

ln(ri+1 ri )

(2.36)

 

 

 

 

 

i=1

2λi

 

 

 

 

Величину

R

 

=

1 ln(r

r )

назовем

линейным

 

lλi

 

i+1

i

 

 

 

 

 

 

2λi

 

 

 

внутренним термическим сопротивлением i-го слоя, а сумму

n

1 ln(ri+1 ri ) суммарным линейным внутренним

RlλΣ =

i=1

2λi

термическим сопротивлением многослойной цилиндрической стенки, (м·К)/Вт.

Для граничных условий III рода (рис. 82)

ql =

 

π(Tf 1 Tf (n+1))

 

 

, Вт м,

(2.37)

n

 

 

 

 

2α1 r

21λ

ln(ri+1 ri )+

 

1

 

+

 

 

 

2α

r

 

1 1

i=1

 

i

 

n+1 n+1

 

150

Рис. 82/

где

1

= R

 

,

1

= R

 

линейные

внешние

 

 

 

 

 

 

lα

 

 

lα

 

 

 

 

2α1r1

1

 

2α2rn+1

2

 

 

 

 

 

 

 

 

 

 

термические сопротивления, (м·К)/Вт. Равенство (2.37) можно представить в виде

 

 

 

 

ql

= kl π(Tf 1 Tf (n+1)),

 

Вт м,

 

 

 

(2.38)

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

kl =

 

 

 

1

 

 

 

 

=

1

 

,

Вт

 

 

n

ln(ri+1 ri )+ 2α

 

1 r

 

Rlα +RlλΣ+Rlα

 

м К

 

 

2α1 r + 21λ

 

1

n+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 1 i=1 i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n+1 n+1

 

 

 

 

 

 

 

линейный

коэффициент

 

теплопередачи,

а

величина

1

= R ,

м К

полное линейное термическое сопротивление.

 

Вт

 

 

lk

 

 

 

 

 

 

 

 

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если

расчет ведут для

плотностей теплового

потока

q1 = k1(Tf 1 Tf 2 )

и q2 = k2 (Tf 1 Tf 2 ),

 

где k1, k2 — коэффициенты

теплопередачи, Вт/(м2·К), рассчитанные для соответствующих поверхностей однородного цилиндра, то справедливо соотношение kl = 2r1k1 = 2r2k2 .

Температура на радиусе r, отдаленном на расстояние rj от границы (j – 1)-го и j -го слоев (см. рис. 82), равна

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]