Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Уравнения математической физики..pdf
Скачиваний:
29
Добавлен:
15.11.2022
Размер:
5.66 Mб
Скачать

Рассмотрим задачу Коши:

 

 

и„ =а2ихя +

t> 0 ,x e R ,

(3.1)

w|,=o =<Р(Х), «*|,=0

x e R

(3-2)

Классическим (регулярным) решением задачи (3.1)-(3.2) называется функция u(x,t) е С2 (/ > 0) П С2 (t > 0), удовлетворяющая всюду при t > 0 уравнению (3.1) и начальным условиям (3.2).

Теорема 3.1. Если <peC2( R ) ,i/s e C \R ) ,fe C l(t> 0), то классическое решение задачи (З.1.), (3.2) существует, единственно и выражается формулой Даламбера

1

 

1

х+ш

 

 

u{x,t) = -

[<р{х+ at) + <р(х - at)] + —

fy ffld g +

 

 

 

. .

° x~at

(3.3)

J

 

. t x+a{t-x)

 

 

v

+ Г- i J/«,r)d £ /r.

l a 0 x-a(l-r)

Для указанных функций <p,y/,f определена правая часть (3.3), поэтому можно формально найти функцию по формуле Даламбера.

Оказывается, найденная функция u{x,t) при каждом фиксированном t> 0 удовлетворяет уравнению (3.1), за исключением разве что конечного числа точек разрыва, а также удовлетворяет условиям (3.2). Более того, найденная функция оказывается обобщенным решением задачи.

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

3.1. Распространение волн в бесконечной струне. Задача Коши

Заметим, что задача Коши (3.1), (3.2) служит математической моделью различных физических явлений. В качестве одного из примеров рассмотрим задачу о малых поперечных колебаниях упругой однородной струны.

Будем считать, что колебания струны совершаются только в

направлении перпендикулярном оси струны х и искомая функция смещения точек струны u(x,t) зависит только от координаты х и времени t . Считается, что постоянная Т и плотность струны р постоянны. Тогда в указанном случае задача (3.1), (3.2) есть математическая модель задачи

о малых

непрерывных колебаниях

упругой однородной струны, где

Т

Fix t)

F(x,t)

- линейная плотность внешних сил,

а2 = — ,f(x ,t) = ------—,где

<р(х) - начальное смещение струны, а <//(х) - начальная скорость.

Хотя в природе бесконечной струны не существует, тем не менее, задача Коши адекватно описывает процесс распространения волны вдоль реальной струны до того момента времени, пока волна не достигла границ. Полезно графически представлять себе профиль струны в задачах о поперечных колебаниях струны при различных начальных условиях.

Задача 3.1. Бесконечная струна с линейной плотностью, равной 3, и натянутая с силой, равной 27, в начальный момент времени находится в покое и смещена относительно положения равновесия, причем функция

 

х е [0,11

смещения определяется равенством <р(х) = • 2 - х ,

х е (\,2\

0,

х R \ [0,2}

В момент времени I = 0 струну отпускают,

в результате чего она

начинает движение. Изобразить графически профиль струны в моменты времени / = — , / = 0,5. Схематично указать процесс дальнейшего

распространения волны, указать его скорость.

1г Решение. Т = 27, р = Ъ=> а - — = 3. Математическая модель колебаний

струны - задача Коши:

 

 

 

 

 

и„ =9и ^ ,

 

t > 0,

хе R,

 

I

,

ч

I

п

(ЗЗЛ )

% =0=Ф)>

*4=0 =0’

 

где <р(х) = Ж[о,1](*) • х + *(1>2](*) • (2 - х ).

 

 

Решение задачи (3.1.1)

есть

частный

случай задачи

(3.1), (3.2) при

f( x ,t ) = y/{t) = 0 и кусочно-гладкой функции . Решение (3.1.1) согласно

формулам Даламбера (3.3)

примет вид

 

и(х, 0

= 1[<Р(Х+ 3t) + <p(x- 3/)].

(3.1.2)

1) t - 0 и(х,0) = <р(х) - начальное положение струны.

 

2) / =

12

1 Удобно сначала построить вспомогательные графики —р(х) сжатием

графика <р(х) вдоль оси ординат, затем сдвигом влево (вправо) на -

построить графики функций — х + - (соответственно — <р X). 2 V

После этого сложить оба графика.

 

1 Г

Р

х

 

Г ’б Г г *

X + — + -(р

3 ) 4

■) 2

 

Построения проводим аналогично предыдущему пункту, с той лишь

разницей, что сдвиг графика функции —<р(х)

проводим на \ влево и

2

2

вправо.

 

О

 

1

 

2

4 ), = 1

5 ) 4

0 I

1

 

з)

1

3)

 

х , - \ =—(рf

+ -<рf

 

4,1

2

1

4J

2 г 1

4)

 

 

 

/К

 

. 1.

 

 

 

 

 

 

Ч(х,~)

 

 

 

 

 

 

4

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

--------------------------------X .

7

 

 

 

 

 

X

 

п

= ^ p ( * - l ) + ^ ( * + l)

Ь

4\

, 1,

6) / = —

5 | 1 [

5 ] 1 /

и

х + - .

12

I

 

7) / > | . Ясно, что <р{х - 3/) • (р{х + 3/) = 0, следовательно, графики

u+(x,t) и u_(x,t) не будут пересекаться.

1

/\

 

1

 

I

 

!

 

1

 

 

 

 

1

 

 

 

 

1

 

 

 

 

1

 

 

 

 

1

 

 

 

 

1

 

 

 

 

1

 

 

 

 

1

 

 

 

 

1

 

 

 

 

----- 1-----

J,

.

J

>

»

v

у

х \

 

1

at

at

|

 

1

 

 

i

 

1

 

 

 

1

 

 

i

 

Две волны: м+(х,/) и u_(x,t) распространяются вдоль оси ох в противоположные стороны с одинаковой скоростью. За время At

абсолютная

величина

смещения

AS1вдоль

оси

ох

функции

^(р{х + си)

 

 

ds

..

AS

..

aAt

 

TI

 

равна aAt, поэтому v = — =

lim —

= lim

---- = a . Итак, очевидно, что

 

 

st

At->o At

Al->о

At

 

 

 

скорость

распространения

волн

равна

а =

f f

прямо

— , т.е.

 

 

 

 

 

 

 

 

\ Р

 

пропорциональна

и обратно пропорциональна -fp В частности, чем

больше натянута струна, тем больше скорость распространения воли вдоль струны. В нашем примере скорость v = 3.

Задача 3.2. Бесконечная струна с линейной плотностью, равной 4, и натянутая с силой, равной 64, находится в прямолинейном положении. В момент времени / = 0 к участку [0,l] струны прикладывается импульс, в результате которого струна приобретает начальную скорость, равную 1 на

этом участке. Изобразить графически профиль струны в момент времени

f = -M i = 0,4).

Схематично

изобразить

процесс

дальнейшего

16

 

 

 

 

 

 

 

 

 

 

распространения волны, указать ее скорость.

 

 

 

 

Решение.

В

обозначениях

из

 

(3.1)-(3.2) а -

= 4, <р(х) = О,

у/(х) = Z[o,i](x) • Математическая модель задачи о колебаниях струны:

 

 

 

 

ии = 1 бм^,

t > 0, х е R ,

 

 

 

 

 

ML o =0’

Ч « о = ^*> = *[о,!](*)•

 

 

Согласно формуле Даламбера (3.3)

 

 

 

 

 

 

 

 

 

 

,

*+4/

 

 

 

 

 

 

 

u(x,t) = -

 

 

 

 

 

 

 

 

 

 

 

х-М

 

 

 

 

Для удобства графического изображения профиля струны в

различные

моменты

времени

введем

первообразную функцию

'F(JC) = JV(£)d£, где a eR . Тогда решение u(x,t) примет вид

 

 

, N 1

х+4t

 

х -4 t

=;Ч'(1 +4/)-1ч'(1-4/).(3.1.3)

 

u(x,t) = -

\w (S)d4 -

 

 

О

_ а

 

а

 

О

 

о

 

 

 

 

 

 

 

 

 

 

 

Разность ¥(х + At) - Т(х - 4/) не зависит от выбора

a e R . Рекомендуется

выбрать

а = inf € R | 'Р(х) * 0}.

В

нашем

примере

а = 0,

т.е.

X

 

 

 

 

 

 

 

 

 

 

4 4 * )= JiK £ )^

Для

нахождения

аналитического выражения

для

о

функции Ч'Сх) удобно воспользоваться непрерывностью функции (как известно из математического анализа, интеграл с переменным верхним пределом от непрерывной функции есть функция непрерывная) и равенством

у

Т ^ ) = 'Р ( х ) + |^ Ж

х

1)

x< 0

'F(JC)=

]У(£)</£ = 0.

 

 

 

 

О

 

 

 

 

x

x

2)

0<х< 1 'F(x) = 4'(0) + f ?(£№=* fld£ = x .

 

 

 

о

о

 

 

 

Jf

 

3)

х > 1

¥(х) = 'F(l) + jy(4)d£ = х .

 

 

 

 

i

 

 

 

0,

х < О,

 

Итак, ТДх) = • х,

0 < х < 1, Теперь с помощью формулы (3.1.3) можно по

 

 

1,

х>1.

 

аналогии с предыдущей задачей графически изобразить профиль струны

вразличные моменты времени.

0)/ = 0 м(х,0) = 0 .

 

 

 

 

 

 

/

 

 

 

 

 

 

 

 

 

1 -

 

 

 

 

 

 

 

 

 

 

 

 

 

___________^

 

 

 

 

 

 

 

 

 

 

X

1

и

(

1 "1

1 >т./

п

1 IT,(

п

"16

1

х,—

= -Ч /

X +

8

1

х —

у

 

16J

8

1

4>

4

Изображаем

графически

функции

-Ч'(х)

и - - Т ( х ) , симметричные

 

 

 

 

 

 

 

 

8

 

8

относительно оси ох. При этом на оси ох и ои удобно изображать разный масштаб. Далее из графика функции ^ЧДх) сдвигом влево на ^

получаем график функции

x + ^-j

и из графика -^Т Д х) сдвигом

1

1

(

О

Затем складываем

вправо на — график функции — Т

х —

4

8

\

 

получившиеся графики.

3 )f = 16

п

4) , Л и *>т

. 4)

I

 

8

 

 

----->

 

x

5) Ясно, что при

i//(x-4t) y/(x + 4t) = 0 , поэтому на рисунке

получим следующее:

 

Скорость распространения волн точно такая же, как и в предыдущей задаче V = 4.

X

Задача 3.3. Найти Ч^х) = ji//(^)d<^, если а = inf {х е /?| if/{x) * о}.

а

 

-2 ,

х е (—3;—l],

1,

х е ( - 1;1),

¥(х) = <-1,

хб[2;3),

4х,

х е [3;4j

О,

JCе Л \ {(- 3;1)U [2;4]}.

Решение. а = - 3. Воспользуемся формулой из предыдущей задачи и найдем функцию Ч'(х).

1)

х < -3

¥(х) = J o ^ = 0

 

 

 

-3

X

 

 

X

2)

- 3 < х < -1

Ч-'(х) = 'Р(-З) + Jyf(4)di; = J(-2)d£ = -2x + 6

 

 

-3

-3

 

 

Щ х) = ^(-1) + fys(£)d<* = -2 • (-1) + 6 +

3)

-1 < х < 1

-1

 

X

 

 

 

+ Jl-^ = 9 + x

 

 

 

-1

 

 

 

X

 

4)

1 < х < 2

vF(x) = vF(l)+ J 0 - ^

= 10

 

 

1

 

 

 

X

 

5)2 < х < 3 vF(x) = vi/(2)+ J(-l)-d£ = 12-

 

 

 

2

 

 

 

ДГ

6)

3 < х < 4

4/(x) = vF(3)+ j4 g •d£ = 2x2 -

 

 

 

3

 

 

 

-X

7)

х > 4

vF(x) = 4/(4)+ |0-tf£ = 14

 

 

 

3

 

 

0,

x e (- oo;-3]

 

 

- 2x + 6,

x e (- 3;-l]

 

 

x + 9,

x E (- l;l]

Ответ: Ч/(х) = - 10,

xe(l;2]

 

 

- x + 12,

xe(2;3]

 

 

2x2 -18,

xe(3;4]

 

 

14,

xe(4;+oo)