Lumf_p1&2(2010)
.pdf4.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