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

Lumf_p1&2(2010)

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

4.2Решение задачи о теплопроводности в бесконечном стержне методом Фурье

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

@u

= a2

@2u

(178)

@t

@x2

 

 

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

u(x; 0) = f(x)

(179)

что соответствует задаче Коши. Сделаем замену переменных

= a2t

91

тогда

@u@t = @u@ @@t = a2@u@

и наше уравнение принимает вид

@u @2u @ = @x2;

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

u(x; 0) = f(x):

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

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

подставляя его в (180), получаем

X(x)T 0( ) = X00(x)T ( )

или

T 0( ) = X00(x): T ( ) X(x)

(180)

(181)

92

Так как левая часть этого уравнение зависит только от , а правая

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

 

T 0( )

= ;

X00(x)

= :

(182)

 

T ( )

X(x)

 

 

 

 

В результате для Т( ) получаем

T ( ) = Ce :

Так как температура стержня должна оставаться конечной при t ! 1, то должно быть < 0, т.е. мы можем положить

= 2:

и

T ( ) = e 2 :

Уравнение для X(x) принимает вид

X00(x) + 2X(x) = 0

93

и его общее решение

 

X(x) = D cos x + E sin x:

 

Тогда частное решение уравнения (180) запишется в виде

 

u(x; ) = (A cos x + B sin x)e 2

(183)

В общем случае в (183) A = A( ), B = B( ) и семейство частных

решений уравнения (180) имеет вид

 

u (x; ) = (A( ) cos x + B( ) sin x)e 2 ; 1 < < 1

 

(184)

Общее решение уравнения (180) записывается как суперпозиция

частных

1

 

u(x; ) =

Z

u (x; ) d

 

1

 

94

или

 

1

 

 

u(x; ) =

Z

(A( ) cos x + B( ) sin x)e 2 d

(185)

1

Неизвестные функции A( ) и B( ) подбираются так, чтобы удовлетворить начальному условию:

u(x; 0) = f(x)

которое примет вид

1

Z

(A( ) cos x + B( ) sin x) d = f(x)

(186)

1

Равенство (186) представляет собой разложение функции f(x) в интеграл Фурье, которое в общем случае имеет вид:

11

f(x) = 2 Z

d Z

f( ) cos ( x) d

(187)

1

 

 

 

 

1

1

 

 

95

или с учетом

 

 

получаем

cos ( x) = cos cos x + sin sin x

(188)

 

 

 

 

 

 

 

 

f(x) =

1

2 Z 1 f( ) cos d cos x+

 

Z

 

 

 

1

 

 

 

 

 

 

 

1

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

Z 1 f( ) sin d sin x d (189)

 

 

 

 

2

 

 

 

 

 

 

 

1

 

Сравнивая (186) и (189), находим

 

 

 

1

f( ) cos d ;

 

 

 

 

A( ) = 2 Z

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

f( ) sin d :

(190)

 

 

 

B( ) = 2 Z

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

96

Подставляя (190) в (185) получаем

11

u(x; ) = 2

Z

d

Z

f( ) cos (x )e

d

(191)

1

 

 

 

2

 

 

 

1

 

1

 

 

 

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

u(x; ) =

1

 

1 f( )

8

1 e 2 cos (x

 

) d9d

(192)

 

 

2

Z

 

>

 

>

 

 

 

 

 

< Z

 

=

 

 

 

 

1

 

 

1

 

;

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

 

интеграл в (192). Для этого сделаем за-

Преобразуем внутренний>

 

 

>

 

 

 

 

 

=

 

 

 

x

= !

 

мену переменной

 

p

 

, введем обозначение

p

 

 

 

, и в

 

 

 

 

97

результате получим

1

e

cos (x ) d = p

1

e

 

cos ! d = p I(!)

Z

Z

 

 

2

1

 

 

2

1

 

1

 

 

1

 

 

 

 

Для вычисления

1

Z

I(!) = e 2 cos ! d

1

найдем его производную

1

Z

I0(!) = e 2 sin ! d

1

98

и выполним интегрирование по частям

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I0(!) = Z

e 2

sin ! d =

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

e sin ! 1

 

 

e cos ! d =

 

 

=

 

 

!

 

!I(!)

2

2

2

1

2

1

1

Z

2

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

Т.о., для функции I(!) получаем дифференциальное уравнение

 

 

 

 

1

 

 

 

 

I0(!)

!

 

 

 

 

 

I0(!) =

 

!I(!)

)

 

 

=

 

 

 

 

2

 

I(!)

2

Отсюда находим

ln I(!) = !2 + ln C 4

I(!) = Ce !2=4

99

Т.к.,

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I(0) = Z

e 2 d = p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(интеграл Пуассона), то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I(!) = p

 

 

e !2=4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и возвращаясь к старым переменным находим

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x )2

2

 

 

 

1

 

 

 

 

 

Z

e

cos (x ) d =

p

 

I(!) = r

 

 

e

4

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

(x )2

 

 

 

2p1

Z

 

f( )e

 

 

 

 

 

 

(193)

 

u(x; ) =

 

4

 

 

 

 

d

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

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