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

3405

.pdf
Скачиваний:
1
Добавлен:
15.11.2022
Размер:
4.89 Mб
Скачать

3.6 Определение количества теплоты, отданного цилиндром в процессе охлаждения

Так же как и для пластины, количество теплоты Qп, Дж, которое отдается или воспринимается поверхностью цилиндра за время от = 0 до = , должно равняться изменению внутренней энергии цилиндра за период полного его охлаждения:

 

 

 

 

Q

п

 

r 2 l

c t

0

 

t

ж

.

(3.63)

 

 

 

 

 

 

0

 

 

 

 

 

За любой промежуток времени от

 

 

= 0 до

1 внутренняя

энергия цилиндра изменится на величину

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

Qп 1

Q .

 

 

 

(3.64)

где по-прежнему, как и для пластины,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t1

t ж

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 0

t ж

 

 

 

 

 

 

 

Средняя безразмерная температура цилиндра найдется из

уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

R

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

2 R dR

2

 

R dR

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R 2

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(R изменяется от 0 до 1).

Если в это уравнение подставить значение согласно уравнению (3.53) и проинтегрировать в указанных ранее пределах, то получим:

 

 

 

 

 

 

 

4J 21

 

n

 

exp

2 Fo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 n J02

 

 

J12

 

 

n

 

 

 

 

 

 

n 1

n

 

n

 

 

или, учитывая,

что J0(

)/J1(

) = /Bi,

 

 

 

 

 

 

 

 

4Bi

 

 

exp

n2 Fo .

 

(3.65)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

2 n

2 n

 

Bi 2

 

 

 

 

 

 

 

 

 

 

 

 

При расчете средней температуры цилиндра

в случае

Fo 0,25 также можно ограничиться одним первым членом ряда (3.65):

4Bi

2

exp

12 Fo .

(3.65')

 

 

21 21

Bi 2

 

 

 

198

Функцию

 

4Bi

2

M Bi можно заранее рас-

 

 

 

21

21

Bi 2

 

 

считать для соответствующих значений Bi и свести в таблицы.

3.7 Охлаждение шара

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

Если обозначить избыточную температуру для любой точки шара = t . tж, то дифференциальное уравнение теплопроводности шара в сферических координатах запишется:

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

.

(3.66)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r 2

 

 

r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а) Граничные условия:

 

 

 

 

 

 

 

 

 

 

 

на поверхности шара при r = r0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r .

 

 

 

 

 

 

 

r r r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

Из условий симметрии задачи в центре шара при r = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 .

 

(3.67)

 

 

 

 

 

 

 

 

 

 

r

r 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) Начальные условия: при

= 0

 

 

 

 

 

 

=

0 = t0 – tж для 0

r

r0.

 

Решая уравнение (3.66) методом разделения переменных

и подчиняя полученное решение условиям (3.67), получим:

 

2 sin

n

 

 

n cos

n sin

n R

exp

2 Fo ; (3.68)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin n cos

 

 

 

 

n R

 

n

n 1

n

 

n

 

 

 

198

здесь = / 0; R = r/r0.

Постоянная в уравнении (3.68) является корнем харак-

теристического уравнения, которое для шара имеет вид:

 

tg

Bi 1 .

(3.69)

Уравнение (3.69) является трансцендентным, имеет бесчисленное множество корней при заданном значении Bi и решается аналогично уравнению (3.14). Значения шести корней уравнения (3.69) для различных Bi приведены в [Л.111, табл.

6.5].

При Bi согласно характеристическому уравнению

(3.69)

n n ,

при этом «начальная тепловая амплитуда» уравнения (3.68)

D n

2 sin

n

n

cos

n

 

2

1

n 1

.

 

n

sin

n cos

n

 

 

 

 

 

 

 

 

С учетом

последнего

уравнения

 

формула

(3.68)

.принимает вид:

 

 

 

 

 

 

 

 

 

 

 

2

1 n 1

1

sin n R exp

 

n 2 Fo .

(3.70)

n R

 

n 1

 

 

 

 

 

 

 

 

 

 

При Bi = 1 согласно уравнению (3.69):

 

2n 1

 

 

 

,

 

 

Dn

2 1 n 1

1

,

 

n

2

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и уравнение (3.68) запишется:

 

 

 

 

 

 

 

1

n

1 2 sin

n R

exp

2

Fo .

(3.71)

 

 

 

 

 

 

 

 

n R

n

n

1

 

 

 

 

n

 

 

 

 

 

При малых значениях Bi (Bi < 0,l) начальные амплитуды (Dn) всех членов ряда (3.68), за исключением первого, стремятся к нулю. Начальная амплитуда первого члена ряда D1 = 1, a 21 = 3Bi При этих условиях соотношение (3.68) запишется

так:

 

 

 

 

 

 

 

 

 

 

sin

3BiR

exp 3Bi Fo .

(3.72)

 

 

 

 

3BiR

 

 

 

198

Из анализа уравнения (3.68) следует, что при значениях Fo 0,25 ряд становится настолько быстросходящимся, что для выражения температурного поля можно ограничиться первым членом ряда:

2 sin

1

1 cos

1 sin

1R

exp

2 Fo .

(3.73)

 

 

 

 

 

 

 

 

sin

1 cos

 

 

 

1R

1

 

1

1

 

 

 

 

Так как n в уравнении (3.68) зависят только от числа Bi, то уравнение температурного поля может быть записано в ви-

де

 

 

 

 

F R, Bi, Fo .

(3.74)

Для центра шара

 

 

 

R 0

F1

Bi, Fo .

(3.74 )

Для поверхности шара

 

 

R 1

F2

Bi, Fo .

(3.74")

Функции, определяемые выражениями (3.74 ) и (3.74) для различных значений чисел Bi и Fo, представлены на рис.

3.13 и 3.14.

198

Рис. 3.13 Зависимость = F1 (Fo, Bi) для центра шара

198

Рис. 3.14 Зависимость = F2 (Fo, Bi) для поверхности шара

198

Рис. 3.15. Зависимость Q/Qп = F (Fo, Bi).

198

Аналогично, как для пластины и цилиндра, количество теплоты, которое отдается или воспринимается шаром за про-

межуток времени от

 

= 0 до

, найдем по формуле

 

 

Q

 

6

 

sin

 

 

 

n cos

 

2

 

 

 

 

 

 

n

 

 

n

 

1

e n2 F0 .

(3.75)

 

Qп

3

 

 

 

sin n cos

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1 n

 

n

 

n

 

 

 

 

 

В уравнении (3.75)

Q

 

4

r 3

c t

 

t

 

. начальная из-

п

 

 

0

ж

 

 

 

 

 

 

 

3

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

быточная внутренняя энергия шара. Из рассмотрения (3.75) следует, что

Q

F Bi, Fo .

(3.76)

 

Qп

 

 

Значения функции (3.76) для различных значений чисел Fo и Вi представлены на рис. 3.15.

3.8 Охлаждение (нагревание) тел конечных размеров

а) Охлаждение параллелепипеда

Рассмотрим охлаждение параллелепипеда в среде с постоянной температурой и с постоянным коэффициентом теплоотдачи на всех его гранях. В начальный момент времени ( = 0) все точки параллелепипеда имеют одинаковую температуру t0. Параллелепипед с размерами 2 x 2 y 2 z является однородным и изотропным. Требуется найти распределение температуры в параллелепипеде для любого момента времени, а также среднюю температуру, необходимую для определения количества подведенной (отведенной) теплоты.

Поместим начало координат в центре параллелепипеда

(рис. 3.16).

198

z

y

-x

 

 

x

 

0

 

 

2 z

 

 

 

-y

-z

2

y

 

 

 

 

2 x

 

 

Рис. 3.16 К охлаждению параллелепипеда

При этом днфференциальное равнение запишется следующим образом:

 

t x, y, z,

a

2 t x, y, z, .

(3.77)

 

 

 

 

 

 

Начальные условия ( = 0)

 

 

t x, y,z

 

t0 const .

(3.78)

При заданных условиях задача симметрична относитель-

но центра параллелепипеда. Если ввести обозначение

= t –

tж, то граничные условия запишутся так:

 

а) для поверхности при

 

> 0

 

x

y

z

x x ;

x x

y y

;

(3.79)

y y

z z .

z z

б) в центре параллелепипеда при > 0

198

x

y

z

0 y y

x 0, 0 z

z

0 x x

y 0, 0 z z

0 x x

z 0, 0 x y

0;

0;

(3.80)

0.

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

Можно доказать, что решение таких задач представляется произведением безразмерных температур для тел неограниченных размеров, в результате пересечения которых образовалось рассматриваемое тело.

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

ператур для трех безграничных пластин:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

y

 

z ,

 

 

 

 

 

 

(3.81)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t x,

t ж ;

 

 

 

t y,

 

t ж ;

 

 

 

 

t z,

t ж .

x

 

 

t 0

t ж

 

 

 

y

 

t 0

 

t ж

 

 

 

z

 

t 0

t ж

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение (3.81) в развернутом виде запишется

следующим образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t x, y, z,

t ж

 

t x,

t ж

 

t y,

t ж

 

t z,

 

t ж

.

(3.81 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 0

t ж

 

 

t 0

t ж

 

t 0

t ж

 

t 0

t ж

 

 

Приведенное решение удовлетворяет как дифференциальному уравнению, так и граничным условиям, описывающим процесс теплопроводности в параллелепипеде.

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

198

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