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

Lumf_p1&2(2010)

.pdf
Скачиваний:
38
Добавлен:
21.05.2015
Размер:
1.03 Mб
Скачать

Упражнение. Проверить, что (193) удовлетворяет уравнению (180) и соответствующему начальному условию.

Далее необходимо вернуться к исходной переменной t: = a2t и, подставляя в (193), получим

 

1

 

1

 

 

 

 

 

(x )2

 

 

 

 

 

 

 

Z

f( )e

 

 

 

 

 

 

u(x; t) =

2ap

 

4a2t d

(194)

 

t

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

Можно проверить, что функция

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

(x )2

 

 

 

 

 

 

 

 

 

 

 

 

4a2t

 

(195)

 

 

 

 

 

2ap

 

 

 

G(x; ; t) =

t

e

 

 

 

 

также является решением исходного уравнения и ее называют

фундаментальным решением уравнения теплопроводности.

Физическим тепловым импульсом называется начальное распре-

101

деление температуры

(0;0

 

xj x0

 

0j> ":

(196)

f"(x) =

;

 

 

u

x

x

< ";

 

 

 

j

 

 

j

 

В этом случае решение задачи будет иметь вид

 

u

x0+"

(x )2

 

 

2ap0

 

Z

e

 

 

 

u(x; t) =

 

4a2t d

(197)

t

x0 "

и по теореме о среднем оно может быть записано следующим образом

 

 

 

 

 

~ 2

 

 

2"u0

 

 

(x )

 

 

 

4a2t

(198)

 

2ap

 

u(x; t) =

t

e

 

Точечный тепловой импульс соответствует " ! 0. Количество теплоты, переданное стержню, пропорционально произведению

2"u0

102

и при " ! 0 должно оставаться конечным. Полагая

2"u0 = 1

получаем, что u0 ! 1 при " ! 0. Т.о., точечный тепловой импульс может быть записан в виде -функции Дирака:

f(x) = (x x0):

Подставляя записанное в таком виде начальное условие в (194), получаем решение

 

1

 

 

(x x0)2

 

 

 

2

 

 

(199)

u(x; t) =

2ap

 

e

4a

t ;

t

которое есть фундаментальное решение G(x; ; t) при = x0. Т.о., мы можем утверждать, что функция

 

 

1

 

 

 

(x )2

 

 

 

 

 

 

 

2

(t t0);

 

 

G(x; ; t

 

t0) =

 

 

 

 

e

4a

(200)

 

 

 

 

 

 

 

 

2ap (t t0)

 

 

 

 

 

 

103

дает температуру в точке x в момент времени t, если в начальный момент времени t = t0 в точке возникает точечный тепловой импульс. Функция G(x; ; t t0) носит название функции влияния точечного источника для неограниченной области или функции Грина, с ее помощью решение задачи записывается в виде:

 

1

 

 

u(x; t) =

Z

f( )G(x; ; t) d

(201)

 

1

 

 

104

4.3Решение задачи о теплопроводности для конечного отрезка

Рассмотрим задачу о теплопроводности на отрезке:

 

ut = a2uxx + g(x; t);

(0 < x < l; t > 0)

(202)

Начальное условие

 

(203)

u(x; 0) = f(x)

и однородные граничные условия

 

 

u(0; t) = 0;

u(l; t) = 0:

(204)

4.3.1 Однородная задача

 

Рассмотрим сначала однородную задачу

 

ut = a2uxx:

(205)

Будем искать решение в виде

 

u(x; t) = X(x)T (t)

105

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

1 T

0

 

X00

 

 

 

 

=

 

 

=

a2 T

 

X

В результате получаем два обыкновенных ДУ:

X00 + X = 0; T 0 + a2 T = 0:

(206)

(207)

(208)

Из граничных условий для u получаем граничные условия для

X:

X(0) = 0; X(l) = 0:

В результате для функции X(x) мы получили задачу о собственных значениях (задачу Штурма-Лиувилля):

X00 + X = 0; X(0) = 0; X(l) = 0:

(209)

Ранее было показано, что собственные значения этой задачи

n = n 2 (210) l

106

соответствующие собственным функциям

 

 

 

X (x) = sin

n

x = sin

n

x

(211)

 

n

 

 

l

 

Далее находим функцию T (t):

 

 

 

 

 

 

 

 

 

 

Tn(t) = Cne a2 nt

 

 

(212)

Таким образом, мы нашли частные решения однородной задачи:

un(x; t) = Cne a2 nt sin

n

x

(213)

 

 

l

 

Общее решение нашей задачи запишем как суперпозицию част-

ных

1 Cne

n

2a2t sin lnx

 

 

 

u(x; t) =

l

(214)

 

X

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

Из начального условия получаем

 

f(x) =

1 C sin nx

(215)

 

 

X

 

n

 

l

n=1

 

107

Последнее выражение есть разложение функции f(x) в ряд Фурье по синусам на интервале (0; l). Для нахождения Cn домножим уравнение (215) на sin lmx и проинтегрируем:

l

 

1

 

l

 

 

 

 

 

 

 

 

Z

f(x) sin

m

x dx = n=1 Cn Z

sin

n

x sin

m

x dx (216)

l

l

l

0

 

 

X

0

 

 

 

 

 

С учетом формулы

sin sin = 12(cos( ) cos( + ))

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

Z

l

lnx sin l x dx =

2 nml:

sin

 

 

 

 

m

1

0

 

 

 

 

 

 

 

 

108

В результате для коэффициента Cn имеем

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cn = l Z

f( ) sin ln d :

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставим в решение найденное значение Cn:

 

 

 

 

1

2

 

2

 

l f( ) sin

n

d 3e

n

2a2t sin

n

x

u(x; t) =

 

l

 

 

 

 

 

n=1

 

 

l

Z

l

 

 

 

 

 

 

 

 

 

l

 

X

 

 

 

 

0

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

Поменяем порядок суммирования и интегрирования

 

 

 

l

2

2

 

 

 

1 e

n

2a2t sin

n

sin

n

x3f( ) d

u(x; t) =

Z

 

 

 

l

 

 

 

 

 

 

 

l

n=1

 

 

 

l

 

l

 

 

 

 

0

6

 

 

X

 

 

 

 

 

 

 

 

 

7

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

(217)

(218)

(219)

109

Введем функцию

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G(x; ; t) = 2l

1 e

n

2a2t sin ln

sin lnx

 

 

l

(220)

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

функцию мгновенного точечного источника или функцию температурного влияния мгновенного точечного источника тепла. С ее использованием решение нашей задачи будет иметь вид

u(x; t) = Z

l

 

G(x; ; t)f( ) d

(221)

0

Покажем, что функция G(x; ; t) представляет собой распределение температуры в стержне в момент времени t, если в начальный момент температура равна нулю и в этот момент в точке x = мгновенно выделяется некоторое количество тепла, при том что на краях стержня поддерживается нулевая температура.

110

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