Матфизика Мурга Е.В
.pdfПідрахуємо інтеграл, що стоїть справа:
|
x |
|
1 |
|
|
2 |
x |
1 |
|
1 |
|
|
|
1 2tg |
|
||||||
sin2 |
k |
dx |
|
1 cos |
k |
|
dx |
|
|
|
|
|
sin2 k |
|
|
1 |
|
|
k |
|
|
|
2 |
|
|
2 |
4 k |
|
|
2 |
|||||||||||||
0 |
|
0 |
|
|
|
|
|
2 |
|
2 k 1 tg k |
і пригадаємо, що tg k k . Тоді p
|
|
k |
|
|
|
p2 2 |
p |
||
sin2 |
|
dx |
|
|
|
|
k |
|
|
|
|
2 |
|
p |
2 |
|
2 |
||
0 |
|
|
|
k |
Знаходимо
|
|
2 |
|
p |
2 |
2 |
|
|
n x |
|
ak |
|
|
|
k |
|
(x)sin |
dx . |
|||
|
|
|
|
2 |
|
|||||
|
|
|
p(p 1) k |
0 |
|
Залишилося знайдені значення аk підставити в ряд і записати відповідь.
Відповідь:
|
|
|
|
2 |
2 |
|
|
a2 k2 |
|
|
|
|
|
u(x,t) |
2 |
|
p |
|
k |
|
e |
|
sin |
k x |
(x)sin |
n x |
dx. |
|
|
2 |
|||||||||||
|
|
|
|
2 |
|
|
|||||||
|
k 1 |
p(p 1) k |
|
|
0 |
|
Приклад 3. Поставити і вирішити крайову задачу про визначення температури стрижня 0 ≤ х ≤ з теплоізольованою бічною поверхнею, якщо його початкова температура є довільною функцією х,
а на кінці стрижня подається ззовні заданий тепловий потік.
Постановка задачі. Хай на бічній поверхні стрижня відбувається конвективний обмін з навколишнім середовищем,
111
температура якого є заданою функцією часу. Нехтуючи деформацією ізотермічних поверхонь, поставимо крайову задачу про визначення температури в стрижні коли на кінці стрижня подаються ззовні постійні теплові потоки.
Рівняння теплопровідності в даному випадку має вигляд:
u |
|
|
1 |
u |
|
|
p |
(u u |
|
) , |
(1) |
t |
|
xx |
|
0 |
|||||||
|
|
c |
|
c |
|
|
При 0 t + , 0 x , де р – периметр поперечного перетину стрижня; - коефіцієнт теплообміну між поверхнею стрижня і навколишнім середовищем; температура якої рівна u0 (у нас - u1 = u0)
с – питома теплоємність; - коефіцієнт температуропровідності матеріалу, з якого виготовлений стрижень; - густина маси; - площа поперечного перетину.
U(x, ) = f (x) = - 4, 0 |
x , |
- u x(0,t)= q1 (t) = 0 |
||||
|
|
|
ux( ,t) |
= q2(t), |
0 t + |
|
Хай, |
x |
a2 , |
p |
b2 |
тоді |
|
c |
|
|
||||
|
|
c |
|
|
ut = a2uхx-b2u + b2u0; u0= 1; = 1 u(x, 0)= f(x)= -4; ux(0,t)= ux(1,t)= 0.
Рішення даного рівняння шукатимемо у вигляді u = v + w, де v = v(x) рішення рівняння
112
|
|
|
a2 |
2v |
b2v b2 0 |
(2) |
||||
|
x2 |
|||||||||
|
|
|
|
|
|
|
|
|
||
при |
v x(0) = v x(1) = 0, |
aw – рішення рівняння |
||||||||
|
|
|
|
dw |
a2 |
d2w |
b2w , |
(3) |
||
|
|
|
|
|
dx2 |
|||||
|
|
|
|
dt |
|
|
|
|||
за умов wx(0,t) = wx (0,1) = 0 |
|
|
|
|
||||||
|
w(x,0)= f(x) – v (x) = - 4 – v (x). |
|||||||||
|
Дійсно, хай u = v + w, тоді |
|
|
|||||||
|
ut = v t + wt = 0 + wt; uxx = v xx + wxx. |
|||||||||
|
Складаючи (2) і (3), отримаємо |
|
||||||||
|
0 |
w |
a2 (wxx vxx ) b2 (w v) b2 |
|||||||
|
|
|||||||||
|
|
t |
|
|
|
|
||||
|
де wxx vxx uxx , w v= u . |
|
|
|||||||
|
Перевіримо ще початкові і крайові умови |
|||||||||
|
u(x,0)= v (x)+ w (x,0) = v (x)+f(x) – v (x)= f(x) |
|||||||||
|
ux(x,t)|x=0 = v (0)+wx(0,t)=0 ; |
ux( ,t)= |
v ( )+wx( ,t) |
|||||||
|
Вирішуємо рівняння |
|
|
|
|
113
a2 |
d2v |
|
b2v b2 |
0, |
|
|
v (0) v ( ) 0, |
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
dx2 |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
d2v |
|
b2 |
v |
b2 |
, |
|
|
v A |
|
, |
|
|
2 |
b2 |
0, |
|
b |
. |
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
0 |
a2 |
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
dx2 |
a2 |
|
|
|
|
|
a2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
||||||||||||
Підставляючи в рівняння, маємо: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
0 – |
b2 |
A |
|
|
|
|
b2 |
|
A |
|
|
|
|
1 |
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
a2 |
0 |
|
a2 |
0 |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Вирішуємо відповідне однорідне рівняння другого порядку |
|
||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
d |
2 |
v |
|
|
b |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
b |
x |
|
|
|||
|
|
|
|
|
|
|
|
|
|
v 0, |
v ea , v |
|
e a |
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
dx2 |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
a2 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|||||||
Загальне рішення однорідного рівняння: |
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
x |
|
|
|
|
|
|
|
|
b |
x |
|
, |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
v 00 = C1e a C2e a |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
v |
|
|
|
|
b |
|
b |
x |
|
|
b |
|
|
|
|
|
|
|
b |
x |
|
, |
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
C |
|
|
|
|
ea |
|
|
|
|
|
C |
ea |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 a |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
(C1 C2 ) 0 , |
|
|
|
|
|
|
|
(4) |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
v |
(0) a |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
(C1e |
b |
|
C2e |
|
b |
|
) 0. |
|
|
|
|
(5) |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
a |
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
v |
(1) a |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
Отримаємо з (4) C = C |
|
та з (5) – C |
|
|
|
b |
|
|
|
|
|
|
|
b |
|
0, |
, звідки C |
C |
0. |
||||||||||||||||||||||||||||||||||
2 |
|
|
ea e |
a |
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 = |
2 = |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
114
Отже, рішенням рівняння (2) є функція
v (x) 1.
Займемося тепер рішенням рівняння (3)
2
dw a2 d w b2w . dt dx2
wx(0,t) = wx ( ,t)= 0; w(x,0) = -4 -1= -5.
Рішення шукатимемо у вигляді
w(x,t)=X(x)T(t)
w |
|
2w |
|
2X |
T(t). |
t |
X(x)T (t); |
x2 |
x2 |
Підставляючи в (3), одержуємо
X(x)T (t) = a2 |
d2X(x) |
b |
2 |
X(x)T(t) . |
dx2 |
|
Розділимо обидві частини рівності на аХТ, отримаємо
b2 T X a2 a2T X
Ліва частина не залежить від x, а права – від t. Рівність
115
можливо лише у тому випадку, коли ці частини постійні
T |
|
|
b |
2 |
|
|
|
X |
|
||
|
|
|
|
|
2 |
|
2 , |
||||
|
|
|
|
|
|
|
|||||
a2T a2 |
|
|
|
X |
|||||||
T a2T( |
b2 |
2 ) 0, |
|
T (b2 a2 2 )T 0, |
|||||||
|
|
||||||||||
|
|
|
|
|
|
|
a2 |
|
|
|
|
x 2x 0, |
|
|
|
||||||||
X A cos x B sin x, |
|||||||||||
dT |
(b2 a2 2 ), |
|
ln |T | (b2 a2 2 )t ln |C|, |
||||||||
|
|
||||||||||
T |
|
|
|
|
|
|
|
|
|
|
|
ln | T | (b2 a2 2 )t, C
2 2 2
T C e (b a )t.
Підставляючи X(x) і T(t) в w(t,x), отримаємо
w(x,t) Ce (b2 a2 2 )t (Acos x Bsin x).
Оскільки T(t) ≠ 0, бо було б w(x,t) 0, то знаходячи похідну х,
одержуємо
wx Сe (b2 a2 2 )t ( A sin x B cos x).
|
Поставляючи сюди х = 0 і х = , маємо Сe |
b2 |
a2 2 |
t |
B 0 |
|
|
|
|
|
|||
звідки В = 0 |
|
|
|
|
|
|
|
|
Ce (b2 a2 2 )tA sin 0, , |
|
|
|
|
звідки |
sin =0, |
n= n; n = 0, 1, 2.., |
|
|
|
|
116
х 0 и w 0
Таким чином, Хn(х)= А cos nx. Враховуючи n, запишемо
T (t) С |
n |
e (b2 na2 )t |
n 0,1,2... |
n |
|
|
Таким чином, власні функції задачі Штурма-Ліувілля знайдені
wn (x,t) Сne (b2 na2 )t Acos nx
і рішення нашого рівняння запишеться у вигляді
2 2
w(x,t) Dne (b na )t cos nx
n 1
Тут Dn=CnА. Залишилося задовольнити початковим умовам:
w(x,0)= f(x) – v(x)= - 5,
w(x,0)= -5 = Dn cos nx .
n 0
Розкладаючи Р(х) = -5 на відрізку 0,1 в ряд по косинусах
(тобто продовжуючи її парним чином), маємо
|
|
|
a0 |
|
|
|
|
|
5 |
an cos nx, |
|
|
|||
|
|
|
|
||||
|
|
2 |
n 1 |
|
|
||
|
1 |
|
|
1 |
sin nx |
|
|
де |
a0 2 5dx 10; |
an 2 5cos nxdx 10 |
|10 0 |
||||
|
|||||||
|
0 |
|
|
0 |
n |
117
таким чином, Dn=an=0 для |
n 1, і рішення має вид |
w(x, |
||
t)=D0e b2t 5e b2t оскільки D0 |
|
a |
5 . |
|
|
|
|||
|
2 |
|
|
Таким чином, рішенням початкової задачі є функція
u(x,t)=v(x)+w(x,t)=1 5е b2t
Перевірка:
2 |
b2t |
; |
uxx = 0; |
2 |
|
2 2 |
b2t |
||
ut=5b е |
|
|
|
-b u= -b +5b е |
|
||||
|
|
|
|
|
ut=a2uxx - b2u + b2; |
|
|
||
2 |
е |
b2t |
2 |
|
b2t |
2 |
|
||
5b |
|
|
|
= 0 – 5 b (1-5е |
|
) + b , |
|
тобто рішення знайдено вірно, ця функція задовольняє також початковим і крайовим умовам.
Задачі для самостійного розв’язування |
|
|||||||||
1. Знайти рішення рівняння |
u |
a2 |
2u |
|
(0 x , |
t 0) , |
||||
|
|
|||||||||
|
t |
|
|
|
x2 |
|
|
|
||
яке задовольняє умовам u(0,t)= u( , t) = 0, t 0 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
||
x при 0 |
х |
|
|
|
|
|||||
2 |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|||
u(x,0) |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|||||
х при |
2 |
х |
|
|
|
|||||
|
|
|
|
|
|
|
118
Відповідь:
|
4 |
|
( 1) |
n |
|
|
e |
(2n 1)2 2a2 |
t sin |
(2n 1) x |
|
|||
u(х,t) |
|
|
|
|
|
2 |
. |
|||||||
2 |
(2n 1) |
2 |
|
|
||||||||||
|
|
n 0 |
|
|
|
|
|
|
|
|||||
2. Даний тонкий однорідний ізольований стрижень завдовжки |
||||||||||||||
, початкова температура |
|
якого |
(x) |
cx( x) |
. Кінці стрижня |
|||||||||
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
підтримуються при температурі, рівній нулю. Визначити температуру стрижня у момент часу t 0.
Відповідь:
|
8c |
1 |
|
|
2n 1)2 |
2a2 |
2n 1 x |
||||
|
|
( |
|
|
t |
||||||
u(x,t) |
|
|
|
|
|
e |
|
|
sin |
|
. |
|
3 |
(2n 1) |
3 |
|
|
|
|||||
|
|
n 0 |
|
|
|
|
|
l |
3. Даний тонкий однорідний стрижень завдовжки , бічна поверхня якого теплоізольована. Початкова температура стрижня відома. На обох кінцях стрижня (х = 0 і х = ) відбувається теплообмін з навколишнім середовищем, температура якого вважається рівною нулю. Визначити температуру стрижня у момент часу t 0.
Відповідь:
|
2 |
|
|
|
n2a2 |
|
|
|
n |
cos |
n x |
sin |
n x |
|
|
||||||||
|
|
|
|
|
|
|
|||||||||||||||||
u(x,t) |
Ane |
|
2 |
t |
|
|
|
|
|
|
|
|
|
|
; |
||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|||||||||||
|
|
n 1 |
|
|
|
|
|
|
|
|
|
|
( 2) n |
|
|
|
|||||||
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
z |
|
|
|
|
|
||||
|
|
n |
|
|
|
|
|
n |
|
|
|
|
|
||||||||||
де An (z) |
n cos |
|
|
|
|
sin |
|
|
dz |
; |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1, 2, 3…. позитивні корені рівняння
119
2ctg |
|
|
|
, |
де |
h |
. |
|
|
|
|||||
|
|
|
|
k |
Вказівка. Крайові умови:
u |
hu| |
|
0, |
u |
hu| |
|
0 . |
|
x 0 |
|
x |
||||
x |
|
x |
|
||||
|
|
|
|
4. Вирішити задачу про охолодження однорідного стрижня з теплоізолйованою бічною поверхнею, якщо його початкова
температура |
u(x,0)= (x). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
Один кінець теплоізолйований, а інший підтримується при |
||||||||||||||||||
постійній температурі u0. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
Відповідь: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u(x,t) u0 ane a2 n2t cos n x; |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
n 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
(2n 1) x |
|
( 1) |
2 |
4u0 |
|
|
|
|
(2n 1) |
|
||
|
|
|
|
|
an |
(x)cos |
dx |
|
; |
|
n |
|
. |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
0 |
|
2 |
(2n 1) |
|
|
|
2 |
||||||
|
|
|
Вказівка. |
Задача |
приводиться |
до інтегрування рівняння |
|||||||||||||||
|
|
2 |
u |
|
|
|
|
|
u 0,t |
|
|
|
|
|
|
|
|
|
|||
|
u |
a2 |
|
за |
умов |
|
|
0; |
u l,t u |
|
; |
u x,0 x . |
|||||||||
|
t |
x2 |
|
|
|
|
x |
|
|
|
|
0 |
|
|
|
|
Рішення слід шукати у вигляді u(x,t)=u0+v(x,t), де v(x,t) - невідома функція.
5. Знайти стаціонарне розподілення температури всередині нескінченного циліндру радіусу , якщо на лівій половині поверхні циліндру (0≤φ<π) підтримується температура –Т, а на правій половині
120