книги / Физика колебаний
..pdfε > 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 04−2δεA 02+ε2 ) 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(1−cos ϕ |
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 |
|
|
|
2∆C |
|
|
|
∆W = q0 |
|
|
|
− |
|
|
= q0 |
|
|
|
|
. |
|||||
C0 |
−∆C |
C0 + ∆C |
|
|
|
−(∆C ) |
2 |
||||||||||
2 |
|
|
|
|
2 C |
2 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|||
Полагая ∆C << C0 , |
можно записать |
|
|
|
|
|
|
|
|
||||||||
|
|
∆W ≈ |
q2 |
|
2∆C |
=W |
2∆C |
|
|
|
|
|
|||||
|
|
|
0 |
|
|
|
|
, |
|
|
|
|
|||||
|
|
2C |
|
C |
C |
|
|
|
|
||||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|||||
|
|
|
|
|
0 |
0 |
|
|
|
0 |
|
|
|
|
|
|
где W0 – энергия, запасенная в конденсаторе до скачка емкости, W0 = q02 / 2C0. Если ввести глубину модуляции параметра
µ = ∆C ,
C0
то общее приращение энергии в контуре за один период колебаний при параметрической «накачке» составит
|
q2 |
|
|
2∆W = |
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 на величину |
||||||
2∆L (∆L << L0 ) |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
1 |
|
1 |
|
|
|
∆W = |
Φ0 |
|
|
− |
|
|
≈W0 2µ, |
||
2 |
L0 |
−∆L |
|
|
|||||
|
|
|
L0 + ∆L |
|
где µ = ∆L / L0.
Отметим, что в линейной колебательной системе при выполнении условия параметрического возбуждения колебаний (условия параметрического резонанса) происходит неограниченное нарастание амплитуды возбужденных колебаний. Связано это с тем, что и поте-
130