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

Учебное пособие 800292

.pdf
Скачиваний:
7
Добавлен:
01.05.2022
Размер:
1.35 Mб
Скачать

 

 

 

 

2

1

6 ( 1)n 1

 

 

 

 

 

(4y 1)dsinn y

 

 

2

 

2

n sh2n

n

 

sh2n

0

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

(4y 1)sinn y

 

 

 

 

 

 

 

 

 

 

 

sinn y d(4y 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

n

 

 

sh2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

n

 

 

 

sh2 n

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6 ( 1)n 1

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sinn y dy

 

 

 

 

 

 

 

 

 

 

 

 

 

n

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n sh2n

 

 

sh2n 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

cosn y

 

 

 

 

 

 

 

 

 

 

 

 

 

6 ( 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n sh2n

 

 

 

 

 

n3 3

sh2n

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6 ( 1)n 1

 

 

 

8 ( 1)n

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

( 1)n(8 6n2 2) 8

.

 

 

n3 3 sh2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n sh2n

 

 

 

n3 3

sh2n

 

 

n3 3

sh2n

Получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 1)n(8 6n2 2) 8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(x,y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

shn x sinn y.

 

 

 

 

 

 

 

 

n

3

 

3

sh2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И тогда окончательно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,y) v(x,y) w(x,y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n y

 

 

 

 

 

 

n y

 

 

 

n x

 

 

 

 

 

 

 

 

 

 

anch

bnsh

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

Anchn x Bnshn x sinn y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 1)n(8 6n2 2) 8

 

1

 

 

n x

n y

 

 

 

 

 

 

 

 

 

 

sin

 

sh

 

 

n

3

 

3

 

 

 

 

2

2

n 1

 

 

 

n

 

 

 

 

 

 

 

sh

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

50

1 sh2n

shn x sinn y .

Замечание. Уравнение Лапласа хотя и является двумерным, но так как в него входят только пространственные переменные (отсутствует время), то также представляется возможным построить график его решения. Так график полученной при решении задачи 4.1 функции имеет следующий вид (рис. 4.1):

Рис. 4.1. График решения задачи Дирихле 4.1

ГЛАВА 5. ОДНОМЕРНОЕ УРАВНЕНИЕ ТЕПЛОПРОВОДНОСТИ

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

51

Рассмотрим тонкий однородный стержень длины l (рис. 5.1). Будем полагать, что боковая поверхность стержня теплоизолирована, а температура u во всех точках поперечного сечения стержня в каждый момент времени одинакова.

 

M

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

~

 

 

 

 

 

l

 

 

x

 

 

 

x =x + x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 5.1. Тонкий однородный стержень

Известно, что в неравномерно нагретом теле происходит движение тепла от более нагретых частей к менее нагретым. В теории теплопроводности принято, что количество тепла Q, прошедшего за время t через площадку S поперечного сечения стержня, вычисляется по формуле (закон Фурье):

Q k

u

S t,

 

(5.1)

 

 

 

x

 

 

где k ‒ коэффициент теплопроводности, u u(x,t) ‒ температура в

момент t в сечении с абсциссой x. Если Q1,

Q2

- количество

тепла, прошедшего за время t через

 

сечения M (x ),M (x ) стержня

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

соответственно, то приток тепла в элемент

 

М'М'' стержня за время

t будет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q1 Q2

kS t

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

x

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x x

 

 

 

 

 

x x

 

 

Применяя теорему Лагранжа, получим

2u

Q1 Q2 kS t x2 x ~x x.

Это тепло пойдет на повышение температуры u u t. Если

t

‒ плотность, c ‒ теплоемкость (т.е. количество тепла, которое надо потратить, чтобы нагреть единицу массы на один градус), то

52

 

Q Q с xS

u

t,

 

 

 

 

1

 

 

2

 

 

t

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

kS t

 

2u

 

 

x с xS

u

t.

 

 

 

 

x2

 

 

 

 

 

~

 

 

 

t

 

 

 

 

 

 

x x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как это соотношение должно быть выполнено для любых

x, t, то

 

 

 

 

 

 

 

 

 

 

 

u

a2

2u

,

где a2

k

.

(5.2)

 

 

 

 

 

 

t

x2

 

 

 

 

 

c

 

Уравнение (5.2) будем называть одномерным уравнением теплопроводности.

Согласно принятой нами классификации, уравнение (5.2) является уравнением параболического типа, так как для него

a11a22 a122 a2 0 02 0.

5.1.Решение задачи Коши для одномерного уравнения теплопроводности методом интегрального преобразования Фурье

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

Всегда можно считать, что стержень совпадет с осью Ox. Тогда математически задача будет сформулирована следующим образом:

найти решение однородного одномерного уравнения теплопроводности

53

ut(x, t) a2uxx(x, t)

( x

,

t 0),

(5.3)

удовлетворяющее начальному условию

 

 

 

u(x, 0) f(x).

 

 

(5.4)

Данная задача называется задачей Коши для однородного уравнения теплопроводности.

Применим для нахождения решения метод разделения

переменных, т.е. будем искать

частное решение

уравнения

ut (x, t) a2uxx(x, t) в виде произведения двух функций:

 

u(x,t) X(x) T (t).

(5.5)

Подставляя в уравнение, будем иметь

 

 

 

 

2

 

 

X(x) T(t) a X

 

(x) T(t),

 

или

 

 

 

 

 

 

 

 

 

 

 

T (t)

 

 

 

 

 

X (x)

.

(5.6)

a2T (t)

 

 

X(x)

 

 

При независимых переменных x и t равенство (5.6) возможно лишь тогда, когда левая и правая его части сохраняют постоянное значение. Обозначим эту постоянную через . Тогда исходное уравнение распадается на два уравнения:

 

 

 

 

 

 

T (t)

 

 

X

 

,

 

(x)

,

 

 

 

 

 

a2T (t)

 

X (x)

или

 

 

 

 

 

T (t) a2T (t) 0,

(5.7)

 

 

 

 

(5.8)

 

X (x) X(x) 0.

Решая уравнение (5.7), найдем его общее решение

 

T (t) C1e a2 t .

(5.9)

 

 

 

54

 

 

 

Поскольку ни в одном сечении стержня (т.е. ни при каком фиксированном x) температура u X (x) T (t) не может

неограниченно возрастать по абсолютной величине при t , то должно быть отрицательно. Положим поэтому 2. Тогда из (5.9)

имеем

 

 

T (t) C1e a2 2 t .

(5.10)

Уравнение (5.8) в этом случае примет вид

 

 

2

 

X (x) X(x) 0

 

и будет иметь общее решение

 

 

X(x) C2 cos x C3 sin x .

(5.11)

Тогда из (5.5), (5.10) и (5.11) получим общее решение уравнения

(5.1):

u(x,t) (A cos x B sin x)

e a2 2t ,

(5.12)

где A C1 C2, B C1 C3. Здесь C1, C2,

C3 , а следовательно,

A и B

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

Полученная функция является при любом фиксированном решением рассматриваемого уравнения, и мы можем, конечно, для каждого значения выбирать произвольные постоянные A и B. Это означает, что A и B можно считать произвольными функциями от :

A A( ), B B( ). Таким образом, в итоге

мы имеем

однопараметрическое семейство частных решений

 

u (x,t) A( )cos x B( )sin x e a2 2t ,

(5.13)

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

.

Полученная функция (5.13), вообще говоря, не удовлетворяет начальному условию. Ранее, при решении методом Фурье волнового

55

уравнения, чтобы удовлетворить заданным начальным условиям, мы суммировали все решения, придавая все допустимые значения параметру , суммируя по всем n , соответствующим допустимым . Естественно и теперь поступить аналогично. Но теперь допустимые значения – вся положительная полуось, а допустимые значения – вся числовая ось. Тогда логично вместо суммирования по дискретному параметру n использовать, как это делалось ранее, интегрирование по непрерывному параметру .

Рассмотрим функцию

u(x,t) u (x,t) d

 

A( )cos x B( )sin x e a2 2td .

 

 

(5.14)

 

 

 

 

Очевидно, что если выбором A( ) и

B( ) будет

обеспечена

законность дифференцирования под знаком интеграла (один раз по t и два раза по x ), то функция (5.14) будет удовлетворять дифференциальному уравнению (5.3).

Для того чтобы функция u(x,t), определяемая равенством (5.14), удовлетворяла и начальному условию (5.4), должно выполняться равенство

 

 

 

 

u(x, 0) A( )cos x B( )sin x d f (x).

(5.15)

 

 

 

 

Предположим, что заданная функция f (x) представима в виде

интеграла Фурье:

 

 

 

 

1

 

 

 

f (x)

d f ( )cos( x) d .

(5.16)

2

 

 

 

 

 

 

 

Так как cos ( x) cos cos x sin sin ,

то (5.16)

можно переписать в виде

 

 

56

 

 

 

1

 

 

f(x)

 

 

 

 

f( )cos d cos x

 

 

 

 

2

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

f( )sin d sin x d .

 

 

2

 

 

 

 

 

 

 

 

Сравнивая это разложение с (5.15), можно сделать вывод, что

 

1

 

 

 

A( )

 

f( )cos d ,

(5.17)

2

 

 

 

 

 

1

 

 

B( )

f( )sin d .

(5.18)

2

 

 

 

 

Подставляя (5.17) и (5.18) в (5.14), получим:

 

1

 

 

 

u(x,t)

f( ) d (cos cos x sin sin x) e a2 2td

2

 

 

 

 

 

 

 

1

 

 

 

 

 

 

f( ) d cos ( x) e a2 2td .

 

 

2

 

 

 

 

 

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

 

1

 

 

 

a2

2t

 

 

 

 

 

 

 

 

 

u(x,t)

 

 

f( ) d

e

 

 

cos ( x) d .

(5.17)

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

( x)2

 

 

 

 

 

 

 

 

 

 

 

e

 

e a2 2t cos ( x) d

 

 

 

 

(5.18)

 

 

 

4a2t

,

2a

 

0

 

 

 

 

 

 

 

 

t

 

 

 

и подставляя (5.20) в (5.19), окончательно получим, что решение задачи Коши для одномерного уравнения теплопроводности

57

ut(x, t) a2uxx(x, t)

( x , t 0),

u(x, 0) f(x)

определяется по формуле

u(x,t)

1

 

 

 

( x)2

d .

(5.21)

 

 

f( ) e

 

4a2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

t

 

 

 

 

 

 

Замечание. Интеграл, стоящий в формуле (5.21), часто бывает неберущимся, но при этом решение выражается через элементарные функции и специальную функцию Лапласа

 

1

 

x

 

2

 

(x)

 

e

 

 

 

 

2 d ,

(5.22)

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

или нормированную функцию Лапласа

 

1

 

x

 

2

 

0(x)

 

e

 

 

 

 

2 d ,

(5.23)

 

 

 

 

 

 

 

 

 

 

 

2 0

 

 

 

 

при этом (x) 0,5 0(x).

5.2. Решение краевой задачи для однородного одномерного уравнения теплопроводности методом Фурье

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

ut a2uxx

(0 x

l,

t 0),

(5.24)

u(x,0) f(x),

 

 

(5.25)

u(0,t) u(l,t) 0.

 

 

(5.26)

 

58

 

 

 

Здесь стержень считается совпадающим с отрезком [0,l ] оси Ox . Решим её методом Фурье разделения переменных. Положим,

как и ранее,

u(x,t) X(x) T(t).

 

(5.27)

Подставляя (5.27) в уравнение (5.24), получим

 

 

2

 

 

T

X

2

 

T X a T X

 

 

 

k .

a2 T

X

 

Откуда

 

 

 

 

 

 

X k 2X 0,

T k 2a2 T 0 .

 

Решим граничную задачу

 

 

 

 

 

 

X k 2X 0 ,

 

 

(5.28)

X(0) X(l) 0 ,

 

 

(5.29)

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

X(x) C coskx D sin kx .

Подставляя первое из граничных условий (5.26), получим

X(0) C cos 0 D sin 0 C 0

 

X(x) D sin kx .

Подставляя второе граничное условие из (5.26), получим

 

X(l) D sin kl 0

 

sin kl 0

 

k k

n

n

,n .

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

В результате получим совокупность решений:

 

 

 

 

X

(x) D sin n x,

n

,

 

 

 

n

 

n

l

 

 

 

 

 

 

 

 

 

 

 

 

 

удовлетворяющих требуемым граничным условиям.

59