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

книги / Физика колебаний

..pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
2.56 Mб
Скачать

Первое выражение можно рассматривать как меняющуюся со временем амплитуду колебаний с периодом 2π/ (ω2 −ω1 ). Максимальное значение амплитуды равно ϕ0 , а минимальное – ϕ0 (æ1 æ2 )/ (æ1 + æ2 ). Колебания второго маятника представляют собой биения

ϕ2 (t) = ϕ0

1

 

æ2

 

 

sin ω2

−ω1 t sin ω2

1 t .

 

 

æ1 +

 

æ2

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

При сильном различии парциальных частот, один из коэффици-

ентов распределения много

больше другого. Пусть,

например,

ν2

>> ν2. Тогда из (6.7)

и (6.8)

следует

æ ≈ ν2 / α ,

 

 

 

æ

2

 

≈ α

2

/ ν2.

 

 

1

2

 

 

 

 

 

1

1

1

 

 

 

 

 

 

1

И поскольку, как правило,

α α

2

<< ν2ν2 , то æ >>1,

 

æ

2

 

<<1. В этом

 

 

 

 

1

 

1

2

1

 

 

 

 

 

 

 

 

 

 

 

 

случае минимальное значение амплитуды первого маятника приближенно равно ϕ0 (12 æ2 / æ1 ), т.е. амплитуда колебаний первого маятника изменяется мало. Максимальная амплитуда второго маятника приближенно равна 2ϕ0 æ2 , т.е. много меньше амплитуды первого

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

ϕ1 ϕ0

t

ϕ2

2ϕ0 æ 2

t

Рис. 6.5

211

лены соответствующие этому случаю колебания обоих маятников. С приближением æ1 к æ2 минимальное значение амплитуды перво-

го маятника уменьшается, т.е. растет перекачка энергии от первого маятника ко второму.

При равенстве парциальных частот

ν1 = ν2 = ν коэффициенты

распределения амплитуд также

становятся

 

равными по модулю

æ1 =

 

æ2

 

= æ, а колебания маятников приобретают вид

 

 

 

 

 

 

ω −ω

 

 

ω +ω

 

,

 

 

 

 

ϕ1(t) = ϕ0 cos

2

1 t

cos

2

 

1 t

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

ω

−ω

 

 

ω

 

 

 

 

 

 

ϕ2 (t) = ϕ0æ sin

2

1 t sin

 

2

2

1 t

 

 

 

 

 

 

 

2

 

 

 

 

 

 

(для одинаковых маятников æ =1 ).

На рис. 6.6 представлены колебания маятников при равенстве парциальных частот ν1 = ν2 = ν. Оба маятника совершают биения,

огибающая которых сдвинута по времени на T / 4 (T = 4π/ (ω2 −ω1 )),

т.е. периодически происходит полная перекачка энергии от одного маятника к другому. В этом и заключается физический смысл поня-

тия сильной

связанности. Время

передачи энергии τ =T / 4 =

= π/ (ω −ω ),

и так как обычно α << ν2 , это время составляет

2

1

 

 

 

 

τ = πν

,

 

 

α

 

т.е. обратно пропорционально коэффициенту связи α.

В любой реальной системе имеет место затухание с постоянной времени τ0. Наличие даже малого затухания может кардинально из-

менить картину колебаний в системе с малой связью, потому что колебания в первом маятнике успеют затухнуть быстрее, чем раскачается второй маятник. Таким образом, сильная связанность в системе проявляется только в том случае, если τ << τ0. При известном затуха-

212

ϕ1 T

ϕ0

t

æϕ0 ϕ2

t

Рис. 6.6

хании в системе это условие накладывает ограничение на величину минимальной связи, при которой происходит сильное взаимодействие двух систем.

6.3. Вынужденные колебания в связанных системах

Рассмотрим теперь вынужденные колебания в системе из двух индуктивно связанных контуров, в один из которых включена внешняя ЭДС õ(t) = õ0 cos t (рис. 6.7). Уравнения колебаний токов I1

и I2 в этих контурах при отсутствии затухания следуют из закона Ома

1

I1dt

+ L1

dI1

=

õ0 cos t + M

dI2

,

C

dt

dt

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

(6.12)

 

1

 

 

 

dI2

 

dI1

 

 

 

I2dt + L2

= M

,

 

 

 

 

C

2

dt

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

213

 

 

 

 

M

 

 

 

 

 

где

 

L1, L2 , C1, C2

индуктивности

 

 

 

 

 

I2

 

 

 

 

и емкости контуров,

M – взаимная

 

 

 

 

 

 

 

 

 

 

 

I1

L1

 

 

 

 

 

индуктивность контуров (для нее

 

 

C1

L2

C2

 

 

всегда

 

выполняется

соотношение

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

M

 

 

L1L2 ). После дифференциро-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вания по времени уравнения (6.12)

 

 

õ0 cost

 

 

 

 

 

 

 

 

 

 

 

 

 

 

принимают вид

 

 

 

 

Рис. 6.7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

2 I

−α I

2

= − õ0 sin t,

 

 

 

 

 

1

1

 

1

1

 

 

L1

(6.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

2

2 I

1

−α I

1

= 0,

 

 

 

 

 

 

 

 

2

 

1

 

 

где ν1 =1/

L1C1 , ν2 =1/ L2C2 , α1 = M / L1,α2 = M / L2 .

 

Колебания в системе будут состоять из собственных колебаний с частотами ω1 и ω2 и вынужденных колебаний с частотой . Для

определения спектра собственных частот положим õ0 = 0 в (6.13) и будем, как и ранее, искать решение в виде I1 = Asin (ωt ), I2 = Bsin (ωt ) . Тогда приходим к уравнениям

(ν12 −ω2 ) A 1ω2 B = 0,

(ν22 −ω2 )B 2ω2 A = 0.

Из условия нетривиальности решения данной системы относительно величин A и B получаем уравнение для определения собственных частот

(ν12 −ω2 )(ν22 −ω2 )−α1α2ω4 = 0.

(6.14)

Собственные колебания нами уже рассматривались, поэтому ограничимся только вынужденными колебаниями. Будем искать решение уравнения (6.13) в виде

214

I1 = I01 sin t, I2 = I02 sin t.

Подставляя данные выражения в (6.13), получаем уравнения для определения амплитуд токов I01, I02:

 

 

 

(ν12 −Ω2 )I01 12 I02 = −

õ0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ν22 −Ω2 )I02 22 I01 = 0.

 

 

 

 

 

 

 

 

Из данной системы находим

 

 

 

 

 

 

 

 

 

 

 

 

 

I01 = −

 

 

 

õ0(ν22

−Ω2 )

 

 

 

,

 

 

 

 

 

 

 

L1

(ν12

−Ω2 )(ν22 −Ω2 )−α1α24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.15)

 

 

 

 

 

 

 

 

 

α2õ03

 

 

 

 

 

 

 

 

I02 =

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

L1 (ν12

−Ω2 )(ν22 −Ω2 )−α1α24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На рис. 6.8 приведены зависи-

 

I01

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

мости амплитуд токов от частоты

 

 

 

 

 

 

 

 

 

 

 

внешней

силы.

В

 

 

точках

Ω = ω1

 

 

 

 

 

 

 

 

 

 

 

и Ω = ω2

амплитуды I01 и I02 обра-

 

 

 

 

 

 

 

 

 

 

 

щаются в бесконечность (это видно

0

 

 

 

 

 

 

 

 

 

 

 

ω1

1

 

ν2

ω2

и без графиков, если сравнить зна-

 

 

 

 

 

 

 

 

 

 

 

менатели

 

в

выражениях

(6.15)

 

 

 

 

 

 

 

 

 

 

 

с уравнением (6.14) для собственных

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

частот). Таким образом, в системе

 

I02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с двумя степенями свободы резонанс

 

 

 

 

 

 

 

 

 

 

 

происходит

на

обеих

собственных

0

 

 

 

 

 

 

 

 

 

 

частотах системы. При

Ω < ν2 токи

 

ω1

1

 

ν2

ω2

 

I1 и I2 совершают противофазные

 

 

 

 

 

 

 

 

 

 

 

колебания,

а при Ω > ν2

синфаз-

 

 

 

 

 

 

 

 

 

 

 

ные колебания.

Точка

Ω = ν2 инте-

 

 

 

Рис. 6.8

 

 

 

 

ресна тем,

что

в

ней

ток

первого

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

215

контура, в который включена внешняя ЭДС, обращается в нуль. Данная ситуация уже обсуждалась нами в подразд. 3.3. Во втором контуре колебания не обращаются в нуль ни при каком конечном значении . В общем случае частота успокоения колебаний всегда равна парциальной частоте того контура, который получается при разрыве цепи в точке включения внешней ЭДС.

Физическая природа отсутствия колебаний в первом контуре связана с тем, что взаимная ЭДС õ1 = MdI2 / dt, наводимая в первом контуре колебаниями второго контура на частоте внешней силы Ω = ν2 в точности компенсирует внешнюю ЭДС. Нетрудно проверить, что

õ1 = M dIdt2 = M I2 cos t = −õ0 cos ν2t.

Поэтому вынужденные колебания в первом контуре на частоте ν2 и не совершаются.

Учет затухания в контурах при вынужденных колебаниях приводит к тому, что амплитуды колебаний при резонансе становятся конечными и при Ω = ν2 амплитуда колебаний в первом контуре не обращается в нуль (кроме тривиального значения Ω = 0 ). Таким об-

разом, демпфирование в системе с потерями никогда не бывает пол-

 

 

 

 

 

 

 

 

ным. Однако уменьшение ам-

 

I1

 

 

 

 

 

плитуды колебаний на частоте

 

 

 

 

 

 

 

 

 

 

 

 

ν2 при малом затухании весьма

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

значительно. На рис. 6.9 приве-

 

 

 

 

 

 

 

 

дена резонансная кривая в пер-

 

 

 

 

 

 

 

 

вом контуре при наличии по-

 

 

 

 

 

 

 

 

терь. Наличие потерь приводит

 

 

 

 

 

 

ω

к уменьшению амплитуды коле-

 

 

 

 

 

 

 

ω ν1

ν

 

1

 

 

 

2

2

 

баний вблизи резонанса и к уши-

 

 

 

 

 

 

 

 

Рис. 6.9

рению резонансной кривой. При

достаточной близости частот

ω1

 

216

и ω2 потери могут превратить резонансную кривую в одногорбую, полностью сгладив провал между частотами.

6.4. Колебания в системе с произвольным числом степеней свободы

В качестве последнего примера системы, состоящей из большого числа связанных осцилляторов, рассмотрим электрическую линию передачи, состоящую из N одинаковых индуктивностей L и конденсаторов C, соединенных как показано на рис. 6.10. Будем полагать, что концы этой цепи закорочены, т.е. напряжения на них равны нулю U0 =UN +1 = 0.

L

qn

L

qn+1

L

L

In1

 

In

 

In+1

 

Un

C

Un+1

C

 

C

Рис. 6.10

Запишем закон Ома для произвольно выделенного участка цепи с номером n, содержащего одну индуктивность и емкость и по кото-

рому протекает ток In:

L dIdtn =Un Un+1 .

Здесь Un = qn / C – напряжение на конденсаторе с зарядом qn. Из закона сохранения заряда следует

dqdtn = In1 In.

217

Из этих уравнений нетрудно получить одно уравнение, связывающее вторую производную qn с величинами заряда в данной точке и двух соседних

q =

1

(q

 

2q + q

 

),

n =1,2,3,..., N.

(6.16)

 

n1

n+1

n

LC

 

n

 

 

 

 

 

 

 

 

 

 

 

Точно так же в силу пропорциональности заряда и напряжения будет выглядеть и уравнение для напряжения Un:

 

Un =

1

 

(Un1 2Un +Un+1 ),

n =1,2,3,..., N.

 

 

LC

 

 

 

 

 

 

 

 

 

 

 

 

 

Предположим, что в случае нормальной моды колебаний с час-

тотой ω временнáя зависимость заряда qn

является гармонической,

т.е. в комплéксном виде

 

 

 

 

 

 

 

 

 

 

 

 

 

qn (t) = Qn exp(iωt),

 

 

(6.17)

где Qn – максимальный заряд. Аналогично

 

 

 

 

qn1 = Qn1 exp(iωt), qn+1 = Qn+1 exp(iωt)

 

(напомним,

что i

мнимая единица). Подставляя эти значения

в уравнение (6.16), получаем

 

 

 

 

 

 

 

−ω2Q eiωt

=

1

(Q

2Q +Q

)eiωt

 

 

 

 

 

 

 

 

n

 

LC

n1

 

n

n+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или после сокращения на exp(iωt )

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

Qn1 + 2

ω2

Qn Qn+1

= 0,

n

=1, 2, 3, ..., N,

(6.18)

 

 

 

ω0

 

 

 

 

 

 

 

где ω0 =1/

LC.

 

 

 

 

 

 

 

 

 

 

Уравнение (6.18) является рекуррентным соотношением, позво-

ляющим найти значение Qn по значениям Qn1

и Qn+1 в двух сосед-

218

них точках. Систему уравнений (6.18) необходимо дополнить соотношениями Q0 = QN +1 = 0. Таким образом, мы пришли к системе N линейных уравнений, решив которую можно найти N разных значений ω2 , причем каждое такое значение ωS будет частотой S-й нор-

мальной моды. Обычно решение подобных систем основано на теории матриц, но мы поступим несколько иначе. Перепишем (6.18) в виде

Qn1 +Qn+1

= 2ω02 −ω2 .

(6.19)

Q

ω2

 

n

0

 

Видно, что при любом заданном значении частоты нормальной

моды ωS правая часть этого уравнения не зависит от номера n,

а это

может быть выполнено только для вполне определенного закона изменения Qn. Предположим, что амплитуда заряда Qn выглядит следующим образом:

 

 

 

 

 

Qn = Asin (nθS ),

(6.20)

где A – постоянная, определяемая из начальных условий, а θS

– не-

которая величина, зависящая только от номера моды S. Тогда

 

 

Q

+Q

 

 

A sin (n 1)θS +sin (n +1)θS

 

 

n1

n+1

=

 

 

=

 

 

Qn

 

 

 

Asin nθS

 

 

 

 

=

2sin nθS cos θS

= 2cosθS ,

 

 

 

 

 

sin θS

 

 

 

 

 

 

 

 

 

т.е. наше предположение о выборе вида Qn полностью оправдалось (левая часть уравнения (6.19) при таком выборе Qn на самом деле не зависит от номера n ). Величину θS можно найти из граничных условий (заметим, что мы еще нигде ими не воспользовались). Условие Q0 = 0 выполняется автоматически в соответствии с (6.20). Из второго граничного условия QN +1 = 0 находим

219

 

 

 

Q

+1

= Asin (N +1)θ

 

=

0

→ θ

 

=

 

πs

 

 

 

 

S

S

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и тогда (6.17) для S -й моды колебаний частоты ωS

 

приобретает вид

 

 

 

 

 

 

 

q

 

=

Asin

 

nsπ

exp(iω

 

t ).

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

+1

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для нахождения разрешенных значений ωS

обратимся к урав-

нению (6.19), полагая в нем ω= ωS :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Qn1 +Qn+1

 

= 2ω02 2

ω2S

= 2cos θS

= 2cos

 

πs

 

.

 

 

 

 

 

 

 

 

N +1

 

 

 

Q

 

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

πs

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ωS = 2ω0

1cos

 

 

 

 

,

 

 

 

 

(6.21)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N +1

 

 

 

 

 

 

 

 

 

 

где s может принимать значения 1, 2, 3, ..., N,

 

а ω0 =1/ LC. Харак-

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

 

N представлен на рис. 6.11

(значения частот отображены стрелками).

 

 

 

 

 

 

 

 

 

 

 

 

 

0,2

 

 

 

0,6

1,0

 

 

 

1,4

 

 

 

1,8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N = 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N = 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N = 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N = 5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ωS / ω0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0,4

 

 

 

 

0,8

 

 

1,2

 

 

1,6

 

 

 

2,0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 6.11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Кроме того, из (6.21) видно, что существует максимальная частота собственных колебаний ωmax = 2ω0 , которая называется крити-

220