3405
.pdf3.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