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

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

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

ε > 0. При дальнейшем возрастании расстройки амплитуда падает и при ε = ε1 происходит скачкообразный срыв колебаний до существенно меньшей амплитуды. При обратном ходе расстройки скачок происходит при ε = ε2 < ε1 и амплитуда при этом резко возрастает.

Значения амплитуд A01 и A02 ,

при которых происходит срыв коле-

баний, определяется из уравнения

 

 

 

 

 

 

 

dA0 (ε)

= ∞

 

 

dε

 

 

 

 

 

 

 

(в точках A01 и A02 касательная к резонансной кривой вертикальна).

Для определения значений ε1, ε2 ,

A01, A02 и A0m запишем урав-

нение (3.31) для резонансной кривой несколько иначе:

 

(ε−δA02 )

2

 

f

0

 

2

 

 

=

 

 

 

−β2 ,

 

 

 

 

 

 

 

 

 

2A0ω

 

или

 

 

 

 

 

 

 

 

 

 

(β2 2 A 042δεA 022 ) A 02=

f0

2 .

(3.32)

 

 

 

 

 

 

 

 

2ω

 

Найдем теперь уравнение для касательной dA0 / dε. Для этого продифференцируем по ε уравнение (3.32), помня, что от ε зависит только A0:

dA

=

−εA A 3

 

 

0

0

0

 

.

dε

ε2 2 4δεA 02+

3δ2 A 04

 

 

Приравнивая выражение для производной dA0 / dε нулю, находим A 02m= ε/ δ. Тогда из (3.31) следует

A 02m= 2βωf0 .

121

Значения A01 и A02 (они лежат на вертикальных касательных

dA0 / dε = ∞ ) можно найти из совместного решения уравнения

 

ε2 2 4δεA 02+3δ2 A 04= 0

(3.33)

и уравнения (3.32) для резонансной кривой.

 

Величина внешней силы f0 , при которой на резонансной кри-

вой появляется гистерезис (неоднозначность), определяется из усло-

вия равенства корней уравнения (3.33)

ε1 = ε2 , т.е. обращения в нуль

его дискриминанта

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4δ2 A4 3δ2 A4

−β2 = 0.

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

Откуда A2

/ δ. При этом ε = ε

2

= 2β и из (3.31) следует

 

 

0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f0 =

 

8ω2β3

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ

 

 

 

 

где значение δ определяется выражением (3.30).

 

 

 

 

 

 

R(I)

 

Рассмотрим теперь влияние нелинейного за-

 

 

 

 

 

 

кона диссипации на резонанс в колебательной

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

системе. Для этого обратимся к электрическому

 

L

 

 

 

 

 

 

 

 

 

колебательному контуру с постоянными реактив-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

ными элементами L, C и нелинейным затуханием

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

(рис. 3.15). Пусть сопротивление зависит от тока

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U0 cosωt

 

по закону R(I ) = R0 (1+ γ0 I 2 ). Это

соотношение

 

 

Рис. 3.15

 

качественно

передает зависимость

омического

 

 

 

 

 

 

 

 

 

 

 

сопротивления проводников от протекающего че-

рез них тока за счет их нагрева. Тогда уравнение, описывающее поведение данного контура при внешнем гармоническом воздействии, можно записать в виде

x +2βx(1+ γx2 )02 x = f0 cos ωt,

(3.34)

122

где β – коэффициент затухания, β = R0 / 2L; γ – коэффициент нели-

нейности; f0 – амплитуда внешней силы, f0 =U0 / L;

ω – ее частота;

ω0 – частота собственных колебаний, ω0 =1/ LC.

Как и ранее вве-

дем в рассмотрение линейную расстройку частоты ε = ω−ω0 << ω0.

В этом случае для достаточно малых значений расстройки уравнение (3.34) преобразуется к виду

x 2 x = f0 cosωt 2βx 2βγx3 +2εωx.

(3.35)

Для применимости метода медленно меняющихся амплитуд потребуем, чтобы все слагаемые в правой части (3.35) были малы по сравнению с членами в его левой части. Решение уравнения (3.35) ищем в виде

x = A(t)cos[ωt (t)], x = −A(t)ωsin[ωt (t)].

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

 

 

A0 =

 

f0

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ω

ε2 2 1+

3

γωβA02

2

 

 

 

 

 

 

 

4

 

 

 

 

В

случае

линейного

контура

 

 

 

 

A0

 

 

(γ = 0)

это

уравнение

совпадает

 

 

 

 

 

 

 

 

 

 

 

 

 

с уже известным нам уравнением (3.9)

 

γ=0

 

 

 

f01 > f02

для амплитуды вынужденных колеба-

 

 

 

 

ний. На рис. 3.16 приведено несколько

 

 

 

 

 

 

 

резонансных кривых для разных зна-

 

 

 

 

 

 

 

чений

f0 (пунктиром отмечены резо-

 

 

 

 

 

 

ε

 

 

0

 

нансные кривые для линейного кон-

 

 

 

 

 

 

 

Рис. 3.16

 

 

тура). Видно, что резонансные кривые

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

123

для контура с нелинейным затуханием являются более плоскими в области расстроек, близких к нулю. Это изменение тем больше, чем больше резонансная амплитуда. Связано это с ростом эффективного затухания. Вдали от области малых расстроек резонансные кривые линейного и нелинейного контуров практически совпадают. В то же время эти кривые нигде не пересекаются и резонансная кривая нелинейного контура при γ > 0 даже вдали от резонанса всегда располо-

жена ниже резонансной кривой линейного контура.

124

Глава 4 ПАРАМЕТРИЧЕСКИЕ КОЛЕБАНИЯ

4.1. Параметрическое возбуждение

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

Простейшая механическая параметрическая система – математический маятник с изменяющейся со временем длиной нити l = l(t)

или с перемещающейся точкой подвеса. Электрический аналог такой системы – колебательный контур с изменяющейся со временем емкостью C = C(t). Математический анализ таких систем приводит

к обыкновенным дифференциальным уравнениям, коэффициенты которых зависят от времени.

Будем различать два класса задач, соответствующих двум классам параметрических систем.

1. Резонансные параметрические системы. К ним относятся системы, для которых характерное время изменения параметров того же порядка, что и характерное время изменения переменных в системе. Например, если частота гармонического осциллятора, описы-

ваемого уравнением x 02 x = 0, зависит от времени и время изменения параметра порядка tхар ~ 2π/ ω0 , то такой осциллятор относится

кклассу резонансных параметрических систем.

2.Нерезонансные параметрические системы. Им соответствуют системы, параметры которых меняются очень быстро или очень медленно по сравнению с характерным временем изменения переменных.

125

 

 

 

 

Классический пример параметрического

 

F

резонанса – раскачивание качелей ребенком,

 

 

регулярно

приседающим

и поднимающимся

 

ϕ0

и тем самым периодически изменяющим поло-

l

жение центра тяжести системы (правда, ребе-

 

 

нок об этом не подозревает!). Для выяснения

 

a

механизма этого способа возбуждения колеба-

 

ний обратимся к математическому маятнику,

mg

 

длину

которого

можно

менять,

подтягивая

 

 

и

опуская

нить,

переброшенную

через блок

 

 

(рис. 4.1). Пусть в момент каждого прохожде-

 

Рис. 4.1

ния через равновесное положение, когда нить

 

 

натянута наиболее сильно, маятник подтягива-

ется внешней силой

 

F

на некоторую небольшую высоту a << l.

А в каждом крайнем положении нить опускается на ту же длину a.

Таким образом, в течение каждого периода колебаний маятник будет

дважды удлинен и укорочен; другими словами, частота периодиче-

ского

изменения параметра

системы

(длины l ) будет в два раза

больше частоты собственных колебаний.

 

 

Поскольку удлинение нити происходит при наклонном положе-

нии маятника, то в эти моменты он опускается на высоту a cos ϕ0

( ϕ0

угловая амплитуда колебаний),

что меньше высоты подъема

в моменты подтягивания нити. Поэтому за каждое подтягивание

и опускание нити внешняя сила производит против силы тяжести ра-

боту

 

 

 

 

 

 

 

 

 

 

 

A

 

= mga(1cos ϕ

0

) 1 mgaϕ2

 

 

mg

 

 

 

 

2

0

 

 

 

 

 

 

 

 

 

 

 

(угол ϕ0 предполагаем малым,

т.е.

cos ϕ0 1−ϕ02 / 2 ). Кроме того,

внешняя сила производит работу против центробежной силы (растя-

гивающей нить) в нижнем положении маятника:

 

 

A 0 a ,

= mv2

цб l

126

где v0 – максимальная скорость маятника. Таким образом, полная работа внешней силы за период колебаний маятника составит

A = 2 1 mgaϕ02 + mv02

a .

2

2

 

 

С учетом того, что v0 = lϕ0ω0 ,

где ω0 =

g / l – собственная час-

тота колебаний маятника, имеем

 

 

 

A = 6 a mv02 .

 

 

l

2

 

 

Видно, что работа, производимая внешней силой над маятником, положительна и пропорциональна его энергии. Поэтому энергия маятника будет систематически возрастать, получая за каждый период небольшое приращение, пропорциональное самой этой энергии и величине µ = a / l. В этом и заключается механизм параметрическо-

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

dWdt W.

Это соотношение того же вида, как и при затухающих колебаниях, с тем, однако, отличием, что производная dW / dt положительна, а не отрицательна. Это значит, что энергия (а с нею и амплитуда) колебаний экспоненциально возрастает со временем.

Точно такая же ситуация может наблюдаться и в электрическом колебательном контуре. Рассмотрим линейный колебательный контур, состоящий из последовательно соединенных L, R и С. Пусть

емкость конденсатора меняется во времени с помощью какого-то

127

внешнего устройства по закону, изображенному на рис. 4.2, а. Предположим, кроме того, что заряд на конденсаторе меняется по закону, близкому к гармоническому.

а

б

Рис. 4.2

При каждом уменьшении емкости С0 на 2∆С энергия конденсатора увеличивается. Заряд же q при скачкообразном изменении емкости не меняется, так как является инерционной величиной. Пусть частотные и фазовые соотношения между q(t) и C(t) таковы, что ем-

кость конденсатора уменьшается в те моменты, когда заряд на емкости проходит через экстремум (рис. 4.2, б), а обратное увеличение емкости – при разряженном конденсаторе. Тогда вложение энергии в систему будет максимальным, так как при раздвижении обкладок конденсатора для уменьшения его емкости совершается максимальная работа против кулоновских сил. Частота изменения параметра (емкости) в этом случае в два раза выше частоты собственных колебаний в контуре.

128

Рассчитаем приращение энергии в контуре W в момент скачка емкости:

2

 

 

1

 

 

1

 

 

 

2

 

 

 

2C

 

 

W = q0

 

 

 

 

 

= q0

 

 

 

 

.

C0

−∆C

C0 + ∆C

 

 

 

(C )

2

2

 

 

 

 

2 C

2

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

Полагая C << C0 ,

можно записать

 

 

 

 

 

 

 

 

 

 

W

q2

 

2C

=W

2C

 

 

 

 

 

 

 

 

0

 

 

 

 

,

 

 

 

 

 

 

2C

 

C

C

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

0

 

 

 

 

 

 

где W0 – энергия, запасенная в конденсаторе до скачка емкости, W0 = q02 / 2C0. Если ввести глубину модуляции параметра

µ = C ,

C0

то общее приращение энергии в контуре за один период колебаний при параметрической «накачке» составит

 

q2

 

2W =

0

4µ =W 4µ.

(4.1)

 

 

0

 

 

2C0

 

Заметим, что, как и в случае раскачки качелей, величина вложения энергии пропорциональна самой энергии.

Определим теперь потери энергии в контуре. Если считать ко-

лебания заряда приближенно гармоническими, т.е.

q = q0 sin ωt, то

средняя мощность потерь составляет 1 Rq2 =

1 Rω2q . Тогда энергия,

 

 

2

2

0

 

 

 

теряемая системой за один период колебаний Т,

 

W

= 1 Rω2q2T = πRωq2 .

(4.2)

дис

2

0

0

 

Сравнивая (4.1) и (4.2), получаем условие, при выполнении которого вкладываемая энергия превосходит потери, и в системе происходит нарастание колебаний

129

 

q2

4µ > πRωq2

→µ > µ

 

=

1

πR

C

=

λ

 

 

 

0

порог

 

0

 

,

(4.3)

 

2C

2

L

2

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

где λ – логарифмический

декремент

затухания

контура,

λ =

= πR C0 / L. Если емкость меняется с той же периодичностью, но по

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

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

Ω = 2nω,

где – частота изменения параметра; ω – частота возбуждаемых колебаний, n =1, 2, 3,... При этом, конечно, скорость вложения энер-

гии в систему будет тем меньше, чем больше n.

Аналогичные соотношения мы получим, если скачком меняется индуктивность L. При этом WL = 12 LI 2. Однако поскольку при изме-

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

изменении индуктивности.

Тогда

при

изменении L на величину

2L (L << L0 )

 

 

 

 

 

 

 

 

 

 

2

 

 

1

 

1

 

 

 

W =

Φ0

 

 

 

 

W0 2µ,

2

L0

−∆L

 

 

 

 

 

L0 + ∆L

 

где µ = ∆L / L0.

Отметим, что в линейной колебательной системе при выполнении условия параметрического возбуждения колебаний (условия параметрического резонанса) происходит неограниченное нарастание амплитуды возбужденных колебаний. Связано это с тем, что и поте-

130