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

Матфизика Мурга Е.В

.pdf
Скачиваний:
29
Добавлен:
06.02.2016
Размер:
996.67 Кб
Скачать

Підрахуємо інтеграл, що стоїть справа:

 

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