Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
2018.pdf
Скачиваний:
25
Добавлен:
07.01.2021
Размер:
2.61 Mб
Скачать

Глава 4. ПРИМЕНЕНИЕ РЯДОВ ФУРЬЕ К ИССЛЕДОВАНИЮ УРАВНЕНИЙ КОЛЕБАНИЙ СТРУНЫ

В гл. 3 описывалось применение рядов Фурье к некоторым ис- Сследованиям упругого изгиба балок [13]. В данной главе мы рассмот-

рим метод решения уравнения колебаний струны методом разделения переменных, который также называется методом Фурье (хотя для случая колебан й струны был предложен еще Даниилом Бернулли).

Фурье.

иущественным для этого метода является также использование рядов

Рассмотр м задачу о сво одных колебаниях струны.

Вывод волнового уравнения.

Струной называется ги кая упругая натянутая нить, не оказы-

вающая сопрот влен я зги у.

Рассмотрим натянутую струну, которая в начальный момент времени совмещена с отрезком [0, l] оси 0х. Пусть концы струны закреплены неподвижно. Если струну тем или иным способом отклонить от первоначального положения и затем предоставить самой себе,

бА§ 1. Уравнения гиперболического типа

то она начнёт совершать колебанияДотносительно положения равновесия, называемые свободными колебаниями (см. рис. 4.1). Требуется найти закон этих колебаний, то есть зависимость перемещения точек

струны от времени x. И Рис. 4.1

149

Будем предполагать, что струна однородна, то есть имеет постоянную линейную плотность ρ [кг·с24] (ρ = γ/g – удельная плотность материала [кг/м3], g = 9,81 м/с2 – ускорение силы тяжести).

Пусть длина струны равна l, а в состоянии равновесия струна прямолинейна и расположена вдоль оси между точками х = 0 и х = l. Если вывести струну из состояния равновесия, подвергнув ее действию какой-нибудь силы, то струна начнет колебаться. Будем считать, что дв жен е всей струны происходит в одной плоскости и что каждая ее точка движется перпендикулярно оси 0х. Смещение точки струны с коорд натой х в момент времени t будем обозначать

через

(x,t) ли просто через и. Предположим далее, что все деформа-

С

удем понимать, что малы как сме-

ции струны малы. Под этим мы

щения

каждого з элементов струны, так и их повороты и'х.

Рассмотр м элемент струны (см. рис. 4.1), который в положении

равновес

меет концами точки х и хх. Пусть в результате откло-

 

струны в некоторый момент времени этот элемент переходит в

нения

 

 

положен е MM′.

 

 

Очев дно, дл на элемента MM′ равна

 

 

x x

1 uх

2dx,

 

 

 

 

 

x

 

 

 

 

бА

что в предположении малости угла поворота элемента (и тем самым тангенса этого угла) приближенно равно х.

и φ + φ. Тогда вертикальная составляющаяДравнодействующей этих двух сил натяжения будет равна

Рассмотрим воздействие на элемент ММ' равнодействующей вертикальных составляющих сил натяжения Т, действующих на его концы. Эти силы действуют в направлении касательных к струне.

Обозначим углы, образуемые этими касательными с осью 0х, через φ И

Т sin ( + ) sin .

Ввиду малости углов φ и φ + φ мы можем синусы заменить тангенсами:

Тtg ( + )-Т tg .

Но тангенсы углов наклона касательных равны производным:

 

 

u

 

 

 

u

 

 

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

x

 

 

 

 

 

 

M

 

 

 

M .

 

 

 

 

150

Сила инерции, действующая на элемент ММ', очевидно равна

2u

t2 x,

где ρ – масса единицы длины струны. Составляя на основании принципа Даламбера (Даламбер Жан Лоран (16 ноября 1717 – 29 октября 1783) – французский механик, физик и математик) уравнение равновесия с л, действующ х на элемент струны, в проекции на ось 0y,

приход м к равенству

 

 

 

 

 

 

 

 

 

 

 

 

2и

 

 

 

 

или 2u

 

 

 

 

 

 

 

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

u

 

 

 

 

 

u

 

 

и

 

Т

 

 

 

 

 

 

 

 

х,

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

x

 

M

 

 

x

 

 

M

х

 

 

 

 

 

 

 

 

, деля обе части этого равенства на

х и переходя к пределу при

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

 

 

 

 

.

 

 

 

 

 

 

 

 

 

x2

t2

 

 

 

Здесь мы воспользовались теоремой Лагранжа о конечных при-

ращениях (Лагранж Жозеф Луи (25 января 1736 – 10 апреля 1813) –

французский механик и математик). Частное приращение производ-

ной

u

при переходе от аргументов (х,t) к аргументам (хх,t) заме-

 

 

x

бАx2

 

 

ной её частным дифференциалом, то есть величиной

2u

.

 

 

 

 

Обозначив, наконец, отношение Т/ρ через а2, мы получим урав-

нение

 

 

 

 

 

 

 

 

 

 

 

a

2 2u

 

2и

,

 

 

(4.1)

 

 

 

 

 

 

 

 

 

 

 

Д2 2

 

 

 

 

 

x

 

t

 

 

 

 

которое и называется уравнением свободных колебанийИструны.

Это уравнение и описывает процесс малых свободных поперечных колебаний струны и называется одномерным волновым уравнением. Это уравнение и его решение впервые были получены Ж.Даламбером в 1743 году.

151

§ 2. Начальные и граничные условия

Уравнению (4.1) удовлетворяет всякое свободное колебание струны, независимо от его физического происхождения, а также от

способов закрепления концов струны в точках х = 0 и х = l.

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

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

.

Из сказанного следует, что для определения движения струны, кроме уравнен я (4.1), нео ходимо еще задать начальные условия, описывающ е поведен струны в начальный момент времени t = 0,

т. е. ту форму, которую струна

при выводе ее из положе-

ниям

 

 

ния равновес я,

 

 

приобретает

 

 

А

(4.2)

и(х, 0) = f(х),

 

итескорости,которыесоо щаютсяточкамструныпри «отпускании»ее:

и х,0 (х).

 

(4.3)

t

 

 

Д

 

Кроме того, необходимо задать граничные условия задачи, т. е.

описать характер поведения концов струны в процессе ее колебаний.

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

концы струны закреплены неподвижно:

И

u(0,t) 0,

 

(4.4)

u(l,t) 0.

 

(4.5)

Разумеется, в частности, может оказаться, что в начальный момент времени струна не имеет отклонения от равновесного состояния (f(х) = 0) или же неподвижна ( (х) 0).

Граничные условия задачи вместе с начальными ее условиями иногда называются краевыми условиями.

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

152

Требуется решить однородное волновое уравнение

 

 

 

a2

2u

 

2и

 

 

 

x2

t2

 

 

 

 

 

при однородных граничных условиях

С

 

u(0,t) 0,

 

 

 

 

 

 

u(l,t) 0.

начальными услов ями

 

 

 

и

 

 

 

 

 

 

(х, 0) = f(х),

 

 

 

и х,0

(х).

 

 

 

t

 

 

 

 

бА

 

§ 3. Решен е однородного волнового уравнения методом

 

разделения переменных (методом Фурье)

В

соответств

с методом

разделения переменных сначала

ищется

частное ненулевое решение однородного уравнения (4.1),

удовлетворяющее лишь однородным граничным условиям (4.4)–(4.5), в виде произведения двух функций, каждая из которых зависит только от одной переменной:

 

 

u(x,t) X (x)T(t).

 

(4.6)

Подставив (4.6) в (4.1), получим

 

 

 

 

X(x)T

 

 

а

2

X

 

(x)T(t).

 

 

(t)

 

 

 

 

Разделив обе части уравнения на

а2X(x)T(t) 0

, приходим к ра-

венству

 

Д

 

 

 

X (x)

 

T2

(t) .

 

(4.7)

 

 

X(x)

 

 

 

 

 

а T(t)

 

 

В этом равенстве при изменении t леваяИчасть, не зависящая от t, остаётся постоянной, поэтому будет постоянной и равная ей правая часть, то есть обе части равенства (4.7) не зависят от t. С другой стороны, при изменении x правая часть равенства, не зависящая от x, будет оставаться постоянной, значит, будет постоянной и не зависеть от x и равная ей левая часть. Таким образом, обе части равенства (4.7) не зависят ни от x, ни от t. Следовательно, они являются постоянными.

153

Обозначая эту постоянную через μ (её называют постоянной разделения), то есть принимая

X (x)

 

T (t)

,

 

а2T(t)

X(x)

 

С

 

 

 

 

получим два независимых обыкновенных линейных однородных

дифференциальных уравнения 2-го порядка

 

 

T

 

2

(4.8)

 

(t) а T(t) 0 ,

 

Х

 

(4.9)

 

(х) Х(х) 0 .

Подстав в далее (4.6) в граничные условия (4.4) – (4.5), получим

 

 

X(0) Х(l) 0.

(4.10)

бА

 

иВ результате для определения функции X(x) приходим к задаче на собственные значен я (4.9), (4.10) (задаче Штурма Лиувилля).

Эта задача меет нулевое (тривиальное) решение 0≡X, не представляющее ф з ческого нтереса, так как тогда

u(x,t) 0.

(Штурм Жак Шарль Франсуа (29 сентября 1803 – 18 декабря 1855) – французский математик, Лиувилль Жозеф (24 марта 1809 – 08 сентября 1882) – французский математик). Однако при некоторых значениях параметра μ, называемых собственными значениями, задача (4.9), (4.10) имеет решения, не равные тождественно нулю. Эти решения называются собственными функциями.

Покажем сначала, что ненулевые решения (4.9), (4.10) существуют при μ˂0. Доказательство от противного. Пусть μ=0. Тогда урав-

нение (4.9) примет вид

Д

 

Х

 

(х) 0.

Его общее решение будет иметь вид И

Х(х) С1х С2.

Определяя произвольные постоянные из граничного условия (4.10), получаем

С1 С2 0.

154

Тогда

Х(х) 0.

Пусть μ > 0 ( 2, 0). Характеристическое уравнение, со-

ответствующее уравнению (4.9), имеет вид

 

 

С

 

 

 

k2 2

0.

 

 

 

 

 

 

 

 

 

 

 

Его корни действ тельны и различны: k1,2

. Поэтому общее

решен е уравнен я (4.9) запишется в виде

 

 

 

 

Х(х) С e х С

e х .

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

Подставляя его в граничные условия (4.10), получаем систему

однородных л нейных алге раических уравнений

 

 

бА

 

 

 

С1

С2

0;

 

 

 

 

 

С е С

е 0,

 

 

и1

 

 

 

 

2

 

 

 

 

которая меет ед нственное нулевое решение С1 = С2 = 0. Соответст-

венно задача (4.9), (4.10) удет иметь только нулевое решение

 

 

 

 

 

 

Х(х) 0.

 

 

 

Пусть μ < 0 ( 2,

0). Тогда корни характеристического

уравнения

 

 

 

 

Д

 

 

 

 

 

 

 

 

Х

 

 

 

2

Х(х) 0

 

(4.11)

 

 

(х)

 

 

будут мнимыми:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i2 2 i.

 

 

 

k

2

 

 

 

 

1,2

 

 

 

 

 

 

 

 

И

 

 

 

 

 

 

 

 

 

 

 

Поэтому частные линейно независимые решения будут

 

 

 

X1(x) cos x,X2(x) sin x,

 

и общее решение уравнения (4.11) запишется в виде

 

 

 

X(x) C1 cos x C2 sin x.

 

(4.12)

Из первого граничного условия (4.10) получаем С1 = 0, а из вто-

рого граничного условия следует

 

 

 

 

 

 

 

 

 

 

 

 

C2 sinlx 0.

 

(4.13)

155

Если С2 = 0, то опять получим нулевое решение

X (х) = 0.

Поэтому, чтобы существовало нетривиальное решение задачи, необходимо принять

 

 

 

 

 

sinlx 0.

 

 

 

 

 

 

 

 

(4.14)

 

Из уравнения (4.14) следует

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l n ;

n 1,2,... (n 0, так как 0).

 

 

Поэтому собственные значения параметра для задачи (4.9),

 

 

(4.10) будут

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С

 

 

 

 

 

 

n

 

 

n

.

 

 

 

 

 

 

 

 

(4.15)

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

Так как со ственные значения будут различными для разных n,

то

пр п сывается соответствующий индекс.

 

 

 

 

им

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соответствующ е им со ственные функции с точностью до по-

стоянного множ теля С2

определяются по формуле

 

 

 

 

 

 

 

 

Х

n

х sin

nx

.

 

 

 

 

 

 

(4.16)

 

 

 

 

 

 

 

 

 

 

 

 

 

бl

 

 

 

 

 

 

 

 

 

 

Здесь уместно о ъяснить причину,

по которой выше было при-

нято рассматривать только положительные значения параметра .

Это было сделано потому,

что при отрицательных значениях или

номера n будут получатьсяАсобственные функции, отличающиеся

лишь постоянным множителем.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим, что собственные функции (4.16) ортогональны на от-

резке [0, l].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С учетом (4.15) уравнение (4.8) запишется в виде

 

 

 

 

 

 

 

 

 

 

Д

 

 

 

 

 

 

an

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

(t)

 

 

 

 

 

 

T (t) 0.

 

 

 

(4.17)

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

l

 

 

 

n

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T (t) A

cos

ant

 

B sin

ant

;

n 1,2,...

(4.18)

 

 

 

 

 

 

n

n

 

 

 

 

l

 

 

 

 

 

n

 

 

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

Подставляя (4.16) и (4.18) в (4.6) и суммируя частные решения

линейного однородного уравнения (4.3), получим

 

 

 

 

 

 

 

 

 

 

 

 

ant

 

 

 

ant

 

nx

(4.19)

 

u(x,t) A

cos

 

 

 

 

 

 

 

B

sin

 

sin

 

.

 

 

 

 

 

 

 

 

 

 

 

n 1

n

 

 

 

 

 

l

 

n

 

 

l

 

l

 

156

Произвольные постоянные An и Bn находим далее из начальных условий. Подстановка (4.19) в (4.5) приводит к равенствам

 

u(x,0) f x A

sin nx .

(4.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

n 1

 

 

 

 

(4.21)

 

u(x,0) x B an sin nx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

t

 

 

 

n 1

 

l

 

 

l

 

 

Если функц

f x и x

удовлетворяют условиям Дирихле,

то про звольные постоянные An и Bn

могут быть определены как ко-

эффиц l

0

 

 

 

l

 

 

 

 

 

енты Фурье для соответствующих функций при разложении

Сих в ряды Фурье по с нусам на промежутке [0, l], равном длине стру-

ны. Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

.

бАu (x,t) A cos B sin sin

 

 

 

2 l

 

 

 

n

 

 

 

 

 

 

An

 

f (x)sin

 

 

xdx;

(4.22)

 

 

B

2

l

 

 

 

 

n

xdx.

(4.23)

 

 

 

(x)sin

l

 

 

n

 

ak 0

 

 

 

 

 

 

 

Выражение (4.19) с учетом (4.22) и (4.23) и даёт окончательное решениезадачи омалыхсо ственныхпоперечныхколебаниях струны.

В полученном решении разные значения n соответствуют различным собственным формам колебаний струны. Функцию

 

n

 

 

n

ant

 

n

 

ant

 

nx

 

 

 

 

 

l

 

 

 

 

 

 

 

l

 

 

 

l

 

 

введением вспомогательного угла

 

 

 

 

 

 

 

легко преобразовать

n

arctg

An

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

к виду

 

 

 

 

 

 

 

Д

 

 

 

u

 

(x,t) F

 

ant

 

 

 

 

 

nx

,

 

(4.24)

 

 

 

sin

 

 

 

 

 

sin

 

 

 

 

 

 

 

l

 

l

 

 

 

 

n

 

 

n

 

 

 

 

 

n

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И

гдеFn

An2 Bn2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из (4.24) следует, что все точки струны совершают гармонические колебания относительно положения равновесия с одинаковыми

частотами n

 

ant

и амплитудой

 

 

 

l

F sin nx , зависящей от про-

n

l

 

дольной координаты x точки струны. Такие колебания называются

157

стоячими волнами. При этом максимальное отклонение от положения

равновесия будет достигаться при sin nx 1, то есть в точках с абс- l

циссами x = (2k +1)l/(2n), (k = 0,1,2,...,n-1) на отрезке [0, l]. Точки, в

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

СвремениПервые три формы колеблющейся струны в разные моменты показаны на р с. 4.2.

sin nx 0. Так х точек на отрезке [0, l] будет (n +1) с абсциссами

l

x = kl/n, (k = 0,1,2,..., n).

Рис. 4.2

Результирующее отклонение u(x, t) произвольной точки струны, как следует из (4.19), равно сумме отклонений, соответствующих раз-

Она соответствует основному тону колебаний струны. Как видно, частота основного тона колебаний тем выше, чем сильнее натянута струна и чем она короче и легче. Высшие тона колебаний называются обертонами.

ным формам коле аний.

 

 

ant

 

 

 

 

 

 

Частоты колебаний

n

 

 

называются собственными час-

 

 

 

бАl

тотами. Наименьшая собственная частота колебаний соответствуетn = 1

иравна

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

Т

.

 

 

 

 

 

 

1

 

 

l

 

l

 

 

 

 

 

 

 

Д

 

 

 

 

 

 

 

 

 

 

 

И

Пример. Найти закон колебания струны длиной l, если в начальный момент струне придана форма кривой

u x l x ,

8l

а затем струна отпущена без начальной скорости. Струна закреплена на концах. Внешние силы отсутствуют.

158

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