Материал: Теоретические основы теплотехники 2

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

При отсутствии внутренних источников тепла дифференциальное урав-

нение теплопроводности имеет вид

t

a

t

2

 

 

2

t

 

 

 

a

 

 

2

 

x

 

 

 

 

2

t

 

 

 

 

y

2

 

 

 

 

2tz 2

  

.

(54)

Уравнение (54) является линейным, однородным дифференциальным уравнением второго порядка в частных производных. Решения такого урав-

нения обладают свойством наложения аналогично решениям обыкновенного

однородного дифференциального уравнения: если t1 и t2 — частные решения

уравнения, то выражение C1t1 C2t2

является также его решением при про-

извольных значениях постоянных С1

и С2. Поскольку у постоянных С1, и С2

возможны различные значения, уравнение типа (54) может иметь бесконечно большое количество частных решений.

Для решения уравнения теплопроводности, удовлетворяющего задан-

ным условиям однозначности, берут сумму частных решений, в которых по-

стоянные Сi имеют определенные значения. Соответствующим подбором по-

стоянных Ci, удовлетворяют решение исходному дифференциальному уравнению и условиям однозначности.

К классическим методам решения уравнения теплопроводности отно-

сятся метод разделения переменных и метод источников.

Метод разделения переменных. По этому методу решается уравнение теплопроводности, а затем, исходя из начальных и граничных условий, опре-

деляются постоянные в общем решении. Частное решение t выражается про-

изведением двух функций, одна из которых U(τ) зависит только от времени τ,

а другая P(x,y,z) зависит только от координат

 

t CU P x, y,z ,

(55)

где С – произвольная постоянная.

Подставляя решение (55) в уравнение (54) получим

26

U P x, y,z aU

P x,

2

 

Уравнение (56) можно переписать так

y,z

.

(56)

U

 

 

2

P x, y,z

 

 

 

 

 

 

 

 

U

a

P x, y,z

.

(57)

 

 

 

Левая часть уравнения(57) может зависеть только от

или быть посто-

янным числом; она не зависит от координат. Правая часть может зависеть от координат или быть постоянным числом; она не зависит от времени. По-

скольку уравнение (57) справедливо при любых значениях времени и коор-

динат, то правая и левая части его равны постоянной величине, которую обо-

значим через D.

Таким образом, мы получим два дифференциальных уравнения для определения вида функций U(τ) и P(x,y,z):

U

 

U

 

D

;

 

 

2

P x, y,z

 

a

 

 

P x, y,z

 

 

D

.

(58)

Решением уравнения (58) является

U CeD ,

(59)

где С – постоянная интегрирования.

Постоянная величина D выбирается из физических соображений. В

большинстве случаев при нагревании или охлаждении тел по истечении дли-

тельного времени температура распределяется в теле определенным образом.

Для тепловых процессов, стремящихся к тепловому равновесию, величина D

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

что физически невозможно. Величина D не может равняться нулю, так как при D=0 функция U(τ) в уравнении (59) имела бы постоянное значение, а тем-

27

пература тела не зависела бы от времени, как это следует из уравнения (55),

что не реально.

Таким образом, из физических соображений следует, что величина D

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

гда экспонента (59) будет периодической функцией времени.

Рассматривая случай, когда D < 0, предположим, что

D аm

2

,

(60)

 

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

С учетом (60) имеем выражение для функции U

Уравнение (58) для

P x,

 

 

U Ce

 

2

 

 

 

 

am

 

.

 

 

 

 

 

 

y,z становится следующим:

 

2

P x, y,z m

2

P x, y,z 0.

 

 

(61)

(62)

Методы решения уравнения (62) излагаются в курсах высшей математи-

ки.

Исходя из того, что при заданных уравнение (62) найдено, и вид функции уравнения (54) примет вид

t Ce

am

 

условиях однозначности решение

P x, y,z известен, частное решение

2

 

 

 

 

 

P x, y,z .

(63)

 

 

Для общего решения уравнения (54) по принципу наложения берут сум-

му частных решений. Постоянная m определяется из граничных условий, а по-

стоянная C – из начальных условий.

28

Метод источников. Метод источников заключается в замене процесса распространения теплоты в теле теплопроводностью совокупностью процес-

сов выравнивания температуры от большого количества элементарных источ-

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

Сущность метода источников покажем на примере неограниченного те-

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

G x, ,

b

 

4a

exp

 

 

x

 

 

2

4a

 

 

.

(64)

Функция G представляет температуру в точке x, если

мент времени в точке

выделяется теплота в количестве

 

ство теплоты на бесконечной прямой равно

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

x

 

 

 

 

dx

 

bc

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

Q c

 

 

exp

 

 

 

 

 

 

 

 

 

 

 

e

 

 

4a

 

2

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

u

;

 

e

u

du

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в начальный мо-

Q bc . Количе-

u

2

du bc ,

(65)

 

 

 

(65а)

Таким образом, количество теплоты Q не зависит от времени. Оно рав-

но произведению площади, ограниченной кривой G и осью абсцисс x, на объ-

емную теплоемкость cp.

Функцию G называют фундаментальным решением уравнения тепло-

проводности, поскольку она удовлетворяет этому уравнению. В самом деле,

для неограниченного тела при одномерном потоке теплоты уравнение (54)

имеет вид

29

t

 

2

t

 

a

x

.

(66)

 

 

2

 

 

Если функция G является решением уравнения (66), его можно запи-

сать так

G

 

G

 

2

 

 

 

 

 

a

x

.

 

 

 

 

 

 

 

 

(67)

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

 

 

2

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пользуясь уравнением (64), найдем выражения для

 

и

x

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

G

 

b

 

x

 

 

1

 

x

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

4 a

 

4a

 

 

2

exp

 

4a

 

;

 

 

(68)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

G

 

b

 

 

 

 

 

x

2

 

4 a

 

 

 

 

 

 

x

 

 

 

2

 

4a

 

 

 

 

 

1

 

 

 

exp

2

 

 

 

x

 

 

2

4a

 

 

.

(69)

Сопоставление последних двух выражений показывает, что действи-

тельно справедливо уравнение (67).

Преобразование Лапласа. Преобразование Лапласа приводит к опера-

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

извольных постоянных.

 

 

 

 

 

Преобразование Лапласа функции f x ,

обозначаемое символом L u ,

называется операцией умножения

f x на

e

ux

с последующим интегрирова-

 

 

 

нием по в интервале от 0 до

 

 

 

 

 

 

 

 

 

 

 

 

L u e ux f

x dx .

(70)

 

0

 

 

 

 

30