CC BY
61
Бутиков Евгений Иванович
МАЯТНИК С МОДУЛИРУЕМОЙ ДЛИНОЙ. II. Области параметрического возбуждения
Аннотация
В первой части этой статьи [1] были рассмотрены основные особенности параметрического возбуждения простого жесткого маятника путем периодической модуляции его длины. Был найден порог возбуждения и установлена роль нелинейных эффектов в ограничении размаха резонансных колебаний. Во второй части статьи исследуется зависимость частотных границ интервалов параметрической неустойчивости от глубины модуляции длины и от добротности маятника. Исследованы стационарные колебания, происходящие на границах интервалов неустойчивости.
Ключевые слова: параметрический резонанс, интервалы неустойчивости, фазовая синхронизация, компьютерное моделирование.
КПГ1П|1|ПТРр
13 унеБнам
ПрОЦЕССЕ
□ □□□□
УДК 534
1. ВВЕДЕНИЕ. ИССЛЕДУЕМАЯ ФИЗИЧЕСКАЯ СИСТЕМА
В данной статье рассматривается хорошо известный способ параметрического возбуждения маятника путем периодического изменения его длины. Во многих учебниках и статьях (см., например, [2-6]) такая физическая система используется в качестве простой модели детских качелей. Действительно, качели можно рассматривать как физический маятник, приведенная длина которого периодически изменяется из-за того, что стоящий на них ребенок сгибает ноги в коленях в положениях наибольшего отклонения и распрямляется каждый раз, когда качели проходят через положение равновесия.
В статье предполагается, что периодическая модуляция приведенной длины маятника происходит по прямоугольному (кусочно-постоянному) закону. Такой (несколько идеализированный) способ модуляции дает
© Бутиков Е.И., 2013
возможность сравнительно простого физического объяснения явления параметрического резонанса, вплоть до количественного анализа (см. [1]). Приближение прямоугольной модуляции применимо к реальному маятнику при условии, что изменения длины происходят в течение промежутков времени, составляющих малую долю периода собственных колебаний. Будем считать, что резкие, почти мгновенные уменьшения и увеличения длины маятника происходят последовательно через промежутки времени, равные Т/2 , так что Т есть полный период изменения длины маятника (период модуляции). Разработанная автором статьи специальная компьютерная программа [7] моделирует поведение этой системы. Компьютерное моделирование существенно помогает в исследовании изучаемого явления.
Легко понять, что изменения длины не могут вывести покоящийся маятник из положения равновесия: в противоположность прямому силовому возбуждению, параметрический резонанс возможен только если
маятник совершает хотя бы небольшие собственные колебания. Будем считать, что изменения длины I маятника происходят между значениями 11 = 1° (1 + т1) и 12 = 10 (1 - тг), где т г - безразмерная глубина модуляции. Параметрический резонанс происходит при следующих значениях частоты (О (или периода Т) модуляции параметра:
2(° г = г = пТ°, п = 2, ... (!)
Т = Т =
± ± п
п 2
Но параметрический резонанс возможен не только при частотах (п, удовлетворяющих уравнению (1), но также в некоторых интервалах частот ( , лежащих по обе стороны от значений оп (в интервалах параметрической неустойчивости). Эти интервалы тем шире, чем больше глубина модуляции параметра. Далее будут получены аналитические выражения для частотных границ интервалов параметрической неустойчивости на основе математически более строгого описания параметрического резонанса при прямоугольной кусочно-постоянной модуляции параметра.
2. ДИФФЕРЕНЦИАЛЬНОЕ УРАВНЕНИЕ
ПАРАМЕТРИЧЕСКИХ КОЛЕБАНИЙ И СРЕДНИЙ СОБСТВЕННЫЙ ПЕРИОД ПРИ БОЛЬШОЙ ГЛУБИНЕ МОДУЛЯЦИИ
В течение интервалов (0, Т/2) и (Т/2, Т длина маятника неизменна, и его движение можно рассматривать как собственное колебание, подчиняющееся соответствующему дифференциальному уравнению. Однако коэффициенты этого уравнения для соседних интервалов времени (0, Т/2) и (Т/2, Т) различаются:
ф + 2уф + (28т ф = 0,
( =
л/1 + т ф + 2уф + (>2 этф = 0,
для 0 < г < Т/ 2, (2)
(2 =
1
т
для - Т/2 < г < 0. (3)
Здесь (о° =^1 - собственная частота малых колебаний маятника со средней длиной 1°, а у - постоянная затухания,
обусловленного вязким трением. При медленном движении маятника в воздухе зависимость силы сопротивления от скорости с хорошей точностью можно считать линейной. Когда сила сопротивления приложена к грузу на конце стержня маятника, и эта сила пропорциональна линейной скорости груза ф I, момент силы сопротивления относительно оси вращения пропорционален ф 12. Поскольку момент инерции в нашей модели (легкий стержень с грузом на конце) также пропорционален 12, постоянная затухания не изменяется при изменении длины маятника (при смещении груза вдоль стержня). Это значит, что ее значения в уравнениях (2) и (3) одинаковы.
Каждый раз в моменты tn = пТ/2 (п = 1, 2, ...), когда длина маятника резко изменяется, нужно совершить переход от одного из уравнений (2)-(3) к другому. На протяжении каждого полупериода Т/2 график движения маятника представляет собой отрезок некоторого собственного колебания. Аналитическое исследование параметрического возбуждения можно выполнить должным «сшиванием» известных решений уравнений (2)-(3) для примыкающих интервалов времени.
Начальные условия для каждого интервала, в соответствии с принятой моделью, выбираются следующим образом. Начальное значение угла отклонения ф берется равным значению угла ф^), которого маятник достигает к концу предыдущего интервала. Начальное значение угловой скорости ф связано с угловой скоростью в конце предыдущего интервала законом сохранения момента импульса:
(1 + т1 )2ф1 = (1 - т1 )2ф2. (4)
В уравнении (4) ф1 - угловая скорость в конце предшествующего интервала, на протяжении которого момент инерции маятника был равен J1 = Т°(1 + т1 )2, а ф2 - начальное значение для следующего интервала, на протяжении которого момент инерции равен J2 = J0 (1 - т1 )2. Скачкообразное изменение угловой скорости при резком изменении момента инерции от значения /2 до /1 находится аналогично. Здесь допустимо
использовать закон сохранения момента импульса в виде (4), потому что за короткое время перемещения груза вдоль стержня влиянием силы тяжести можно пренебречь. Иными словами, можно считать, что маятник свободно вращается вокруг своей оси. Это предположение справедливо, если сдвиг груза вдоль стержня происходит за время, составляющее малую часть периода собственных колебаний.
Определяя условия, при которых уравнения (2)-(3) дают решения с нарастающей амплитудой, можно найти интервалы частоты со вблизи резонансных значений (0„ = 2ю0/ п, в пределах которых состояние покоя маятника неустойчиво при данной глубине модуляции Ш1. В таких интервалах параметрической неустойчивости сколь угодно малого отклонения маятника из положения равновесия достаточно для прогрессивного роста первоначально малых колебаний.
В первой части статьи [1] был найден порог параметрического возбуждения маятника в условиях резонанса, когда два цикла модуляции параметра происходили на протяжении одного периода или трех периодов собственных колебаний. Полученная оценка справедлива при малых значениях глубины модуляции т1 длины маятника. Для больших значений глубины модуляции т г само понятие собственного периода требует
уточнения. Пусть Т0 = 2л/щ = 2лЛ//0/g -период собственных колебаний маятника, когда массивный груз закреплен в среднем положении, то есть эффективная длина равна 10. Когда груз смещен вниз, период ста-
новится несколько больше: Т1 = Т0Л/1 + т1 » » Т0 (1 + тг/ 2). Напротив, период укорачивается, когда груз смещен вверх к оси:
Т2 = Т0л11 -т »Т0(1 -т//2).
Для дальнейшего удобно ввести средний собственный период не как арифметическое
среднее ^ (Т1 + Т2), а как период, соответствующий среднему арифметическому значению частоты (Уау = 1((+ю2), где ( = 2л/Т и со2 = 2л/Т2 . Таким образом,
средний период Тау определим следующим соотношением:
2л = 2Т1Т2
Т =-
ау
(5)
®ау (Т1 + Т2)
Действительно, период Т модуляции параметра, соответствующий точной настройке на один из резонансов, зависит не только от порядкового номера п резонанса, но также и от глубины модуляции т1. Для выполнения резонансных условий приращение фазы собственных колебаний за один цикл модуляции должно быть равно л, 2л, 3л, ..., пл, ... На протяжении первого полупериода фаза возрастает на со1 Т/2, а на протяжении второго полупериода - на со2 Т/ 2. Поэтому вместо приближенного условия (1) получаем:
( + (02
2
Т = пл,
л
или Т = Тп = п-= п
Т
од
(ау 2
(6)
Таким образом, для параметрического резонанса некоторого порядка п условие точной настройки можно выразить через значения двух собственных периодов Т1 и
Тау
Т2. Это условие имеет вид Т = п^-, где Тау
определяется соотношением (5). Для малых и умеренных значений глубины модуляции т I можно использовать следующие приближенные выражения для средних значений частоты и периода:
ю0
1
+
1
1+т л/1 - т
с
, з 2
1 + -т, 8 1
Тау » Т0
1 - 3 т} 8 1
(7)
Различие между Тау и Т0 проявляется в членах, квадратичных по глубине модуляции.
3. ЧАСТОТНЫЕ ИНТЕРВАЛЫ ПАРАМЕТРИЧЕСКИХ РЕЗОНАНСОВ НЕЧЕТНЫХ ПОРЯДКОВ
Для нахождения границ интервалов параметрической неустойчивости, содержащих резонансные значения Т = Тау/2,
(ау =
2
T = Tav, T = 3Tav/2, ..., нужно рассмотреть стационарные колебания бесконечно малой амплитуды, происходящие при совпадении периода модуляции с одной из этих границ. Такие стационарные периодические колебания можно представить как чередование собственных колебаний с периодами T и T2.
3.1. ОСНОВНОЙ ИНТЕРВАЛ ПАРАМЕТРИЧЕСКОЙ НЕУСТОЙЧИВОСТИ
Рассмотрим сначала окрестность основного резонанса, происходящего при T = Tav/2. Предположим, что период модуляции T соответствует левой границе интервала неустойчивости, то есть что он немного меньше резонансного значения. В таком случае между последовательными резкими увеличениями и уменьшениями длины маятника проходит чуть меньше четверти среднего собственного периода Tav. Стационарный режим колебаний постоянной амплитуды в отсутствие трения может реализоваться только при условии, что резкие увеличения и уменьшения угловой скорости равны по величине. Графики угла отклонения j(t) и угловой скорости j(t) для такого периодического процесса имеют характерный вид, показанный на рис. 1. Сегменты графиков собственных колебаний (происходящих на интервалах, соответствующих неизменной длине маятника) представляют собой чередующиеся отрезки синусоид и ко-
синусоид с периодами Т1 и Т2. Эти сегменты симметрично обрезаны с двух сторон.
Примем момент, когда груз смещается вниз, за начало отсчета времени t = 0. В этот момент величина угловой скорости резко уменьшается (см. рис. 1). Затем на протяжении промежутка времени ( 0, Т/ 2) график описывает отрезок собственного колебания
с частотой щ = щ / л/1 + m. Поскольку график симметричен относительно момента времени Т/ 4, соответствующие зависимости от времени для j(t) и j (t) можно записать в виде:
j1 (t) = - A1 cos щ (t - T/ 4), j1 (t) = Ащ sin w1 (t - Т/ 4),
0 < t < T12. (8)
Аналогично на интервале (- Т/ 2, 0) график представляет собой отрезок собственного колебания с частотой щ2 = (О0/л11 - m: j2(t) = - А2 sin W2 (t + Т/ 4), j2 (t) = - А2щ2 cos щ2 (t + Т/ 4),
- Т/ 2 < t < 0. (9)
Для определения постоянных А1 и А2 воспользуемся условиями, которые должны выполняться при «сшивании» отрезков графика, и учтем периодичность стационарного процесса колебаний. При t = 0 угол отклонения одинаков для j1 и j2, то есть j1(0) = j2(0). Угловая скорость при t = 0
Рис. 1. Фазовая траектория и графики стационарных параметрических колебаний на нижней границе основного интервала неустойчивости (вблизи Т = Тау /2)
претерпевает скачок, величина которого следует из закона сохранения момента импульса: (1 + т1 )2<^1(0) = (1 - т1 )2 ф2 (0), см. (4). Из этих условий получаем следующие уравнения для А1 и А2:
А1ео8(ю1 Т/ 4) = А281п(ю2 Т/ 4), (10)
А1 (1 + т1 )2 ( 81п( ( Т/ 4) = = А2(1 - т1)2ю2оо8(ю2 Т/ 4). (11)
Однородные уравнения (10)—(11) для А1 и А2 совместны только при равенстве нулю определителя, то есть при выполнении следующего условия:
(1 + т1)2 ю1^1п(ю1ТI4) 81п( ю2 Т/ 4) =
= (1 - т1)2ю2со^(ю1Т14) ео8(ю2 Т/ 4) .(12)
Это уравнение определяет период модуляции Т (при заданном значении глубины модуляции тг), который соответствует левой границе интервала параметрической неустойчивости. Преобразуем уравнение (12) к форме, удобной для получения численного решения относительно неизвестной величины Т:
(Ч +1) ео8(юау Т/2) = (ч -1) Ш5(Д®Т/4),
(13)
где юау = (ю1 + со2)/2и Дю = ю2 - ю1. В уравнении (13) введена безразмерная величина д, зависящая от глубины модуляции т г:
ч3/2
Ч =
(14)
1+т 1 - т
V У
Чтобы найти левую границу Т- интервала неустойчивости, содержащего основной параметрический резонанс, нужно искать решение Т уравнения (13) в окрестности Т0 = Т/2. Заменим Т в аргументе косинуса в левой части (13) на Тау/2 + ДТ. Поскольку юауТау = 2л, косинус можно представить в виде - 8т(юау ДТ/ 2). Тогда уравнение (13) принимает вид:
81п(®ау ДТ/2) = ■
Ч -1 соч Дю(Та^2 + ДТ)
соь .
Ч +1 4
(15)
Это уравнение для неизвестной величины ДТ можно решать методом итераций. В качестве приближения нулевого порядка можно взять значение ДТ = 0 и подставить его в правую часть уравнения (15). Тогда левая часть (15) дает значение ДТ первого порядка. Подставляя ДТ первого порядка в правую часть (15), в левой части получаем значение ДТ во втором порядке. Эту процедуру итераций можно повторять до тех пор, пока не будет получено самосогласованное значение ДТ. Выполняя такие вычисления для разных значений глубины модуляции тг, можно найти всю левую границу Т- (т1) первого интервала параметрической неустойчивости. Далее будет показано, как можно получить правую границу этого интервала, а также границы других интервалов неустойчивости.
Интервалы неустойчивости в плоскости Т - тц для первых шести параметрических резонансов, рассчитанные численно с помощью вышеописанной процедуры, показаны на рис. 2. Здесь можно увидеть аналогию с известной диаграммой Айнса-Стрэтта для интервалов параметрической неустойчивости системы, описываемой уравнением Ма-тье, в частности, для маятника с вертикальными осцилляциями точки подвеса.
Чтобы наблюдать стационарные колебания, соответствующие левой границе интервала неустойчивости (см. рис. 1) в моделирующем эксперименте, недостаточно выбрать в качестве периода модуляции Т самосогласованное решение уравнения (15)для заданного значения глубины модуляции тг: необходимо также должным образом задать начальные условия. Это можно сделать на основе выражения (8), согласно которому при произвольном (но малом) начальном отклонении ф1(0) начальная угловая скорость должна быть равна фф1(0) = (1ап(( Т/4)ф1(0).
Для правой границы основного интервала параметрической неустойчивости период модуляции Т немного больше резонансного значения Т = Тау /2. В этом случае между последовательными увеличениями и уменьшениями длины маятника происходит
Рис. 2. Интервалы параметрической неустойчивости при прямоугольной кусочно-постоянной модуляции длины маятника в отсутствие трения
чуть больше четверти среднего собственного периода Tav. Графики угла отклонения и угловой скорости для такого периодического процесса показаны на рис. 3.
Соответствующие зависимости от времени j(t) и j(t) на интервале (0, T/2) можно записать в следующем виде:
j1(t) = B1 sin col(t—T/4), j1 (t) = B1o1 cos w (t - T/4),
0 < t < T12 . (16)
На интервале (-T/2 , 0) график (см. рис. 3) представляет собой отрезок собствен-
ного колебания с частотой w2 =
V1 —
m
j2 (t) = — B2 cos 0)2 (t + T / 4), j2(t) = B2O2 sin 02(t + T/4),
-T12 < t < 0. (17)
Дальнейшие вычисления аналогичны тем, что описаны после уравнений (8)-(9) для левой границы. Оказывается, что DT для правой границы интервала находится решением уравнения, которое отличается от (15) противоположным знаком в правой части. Решая его численно методом итераций для разных значений m l, получаем правую границу основного интервала параметрической неустойчивости (см. рис. 2).
Для получения приближенных аналити-
Рис. 3. Стационарные параметрические колебания на верхней границе основного интервала
неустойчивости (вблизи Т = Тау/2 ).
ческих решений уравнения (15), справедливых при малой глубине модуляции тг, можно упростить его правую часть, учитывая, что д »1 + 3т,, д -1» 3 т,. Значение косинуса можно приближенно положить равным 1. В левой части (15) синус можно заменить его малым аргументом, в котором юау = 2к/ Тау. В результате получаем следующие приближенные выражения для обеих границ основного интервала, справедливые вплоть до квадратичных членов по т1:
Т+ = 1
+ 2
1 - 3т, 1 + —-
к
Т = 1 Ту 2
1 - 3т,
1 + —1-
3т,2 Л
к
Тн
(18)
3.2. ИНТЕРВАЛ ТРЕТЬЕГО ПОРЯДКА
Границы интервалов неустойчивости вокруг параметрических резонансов высших
порядков можно найти аналогично. При резонансе третьего порядка (п = 3) два цикла изменения длины маятника происходят на протяжении приблизительно трех периодов собственных колебаний (Т » 3Тау/2). Фазовые траектории и графики стационарных колебаний на левой и правой границах третьего интервала показаны на рис. 4.
Фазовая траектория этих периодических колебаний замыкается после двух циклов модуляции. Она образована двумя концентрическими эллипсами, соответствующими малым колебаниям маятника с частотами ю1 ий2 . Изображающая точка обходит эту траекторию по часовой стрелке, перескакивая с одного эллипса на другой каждый раз, когда груз смещается вдоль стержня маятника. Цифры на рис. 4 помогают проследить движение изображающей точки: эквивалентные точки фазовой траектории и графика
Рис. 4. Фазовые траектории и графики угловой скорости и угла отклонения для стационарных колебаний на левой и правой границах интервала неустойчивости вблизи Т = Тгу/2
8
угловой скорости отмечены одинаковыми цифрами.
Рассматривая условия сшивания графиков собственных колебаний с частотами щ и со2 для левой границы, получаем такие же уравнения (10)-(11) для А1 и А2 и такое же уравнение (15) для периода модуляции. Так будет для всех интервалов нечетных порядков. Аналогично, для правой границы получаем такие же уравнения для В1 и В2, как и в случае п = 1, а также уравнение (15) с противоположным знаком в правой части для нахождения соответствующего периода модуляции Т. Но для интервала третьего порядка нужно искать решение этих уравнений в окрестности Т = 3 Тау /2, как и для любого другого интервала нечетного порядка п - в окрестности Т = п Тау /2. Границы интервалов третьего и пятого порядков, полученные путем численного решения методом итераций, также приведены выше на рис. 2.
При малых значениях глубины модуляции Ш1 можно найти приближенные аналитические выражения для нижней и верхней границ третьего интервала, справедливые вплоть до квадратичных по тг членов:
г о л
Т+ = 3
+ 2
1 - т
1+—-р
Т = 3 Тау 2
тI 3т2
1 + ^ -р
8
Т
о >
тг <<1. (19)
В этом приближении третий интервал характеризуется такой же шириной (3тг /р)Т0 что и интервал неустойчивости в окрестности основного резонанса. Однако этот интервал отличается большей асимметрией: его середина смещена влево (в сторону меньших Т ) от резонансного значения
3 9 о
Т = -Т0 на —т2Ъ 2 0 16 1
0
4. ПАРАМЕТРИЧЕСКИЕ РЕЗОНАНСЫ ЧЕТНЫХ ПОРЯДКОВ
При малой и умеренной глубине прямоугольной модуляции длины маятника параметрический резонанс порядка п = 2 (один цикл модуляции на протяжении одного периода собственных колебаний) оказывается сравнительно слабым по сравнению с рас-
смотренными выше резонансами п = 1 и п = 3 . В случае п = 2 на протяжении каждого собственного колебания сдвиги груза приводят по одному разу и к увеличению, и к уменьшению энергии маятника. Рост колебаний происходит тогда, когда вложение энергии в момент сдвига груза вверх будет больше, чем уменьшение энергии в момент сдвига вниз. Это возможно лишь при условии, что груз сдвигается вверх тогда, когда угловая скорость маятника больше, чем когда груз смещается вниз. При Т » Тау такое условие может выполняться только благодаря существованию (малого) различия собственных периодов Т1 и Т2 , где
Т1 = Т^1 + тг - период при смещенном вниз
грузе, а Т2 = Т0лу1 - т1 - при смещенном вверх грузе маятника. Это различие собственных периодов пропорционально т г.
Рост колебаний при параметрическом резонансе второго порядка иллюстрируется рис. 5. Отметим асимметричный характер этих колебаний: угловое отклонение маятника в одну сторону больше, чем в противоположную. При этом результирующее вложение энергии за период пропорционально квадрату глубины модуляции т1 в отличие от резонансов п = 1 и п = 3, для которых вложение энергии пропорционально первой степени т1. Поэтому при одной и той же постоянной затухания у (одной и той же добротности О) для превышения порога параметрического возбуждения при п = 2 необходима значительно большая глубина модуляции. Рост амплитуды здесь также ограничивается нелинейными свойствами маятника.
Интервал неустойчивости в окрестности резонанса п = 2 (для малых значений т 1) значительно уже интервалов п = 1 и п = 3 . Его ширина также пропорциональна лишь квадрату глубины модуляции тг.
Для нахождения границ этого интервала неустойчивости можно представить стационарные колебания при Т » Т0 как образованные чередующимися отрезками собственных колебаний с периодами Т1 и Т2. Фазовая траектория и графики угловой ско-
Рис. 5. Фазовая траектория и графики угловой скорости ф(/) и угла отклонения ф(/) при колебаниях, соответствующих параметрическому резонансу второго порядка п = 2 (Т » Тау)
рости ф^) и угла отклонения ф^) стационарных колебаний на одной из границ показаны на рис. 6. При колебаниях на границе интервала неустойчивости резкие увеличения и уменьшения угловой скорости в отсутствие трения в точности компенсируют друг друга.
Для описания таких стационарных колебаний малой амплитуды можно воспользоваться следующими выражениями для ф^) и ф(I) на интервале (0, -Т/2) (см. рис. 6):
( (t) = - A1 cos w1 (t - T/ 4), (p1 (t) = A1w1 sin w1 (t - T/ 4),
0 < t < T/ 2, (20)
а на интервале (- T/2 , 0)
(2 (t) = A2 cos®2(t + T/4), ((2 (t) = - А2Ю2 sin W2 (t + T/ 4),
- T/ 2 < t < 0. (21)
Условия «сшивания» графиков при t = 0 такие же, как и для других резонансов,
Рис.6. Стационарные параметрические колебания на левой границе интервала неустойчивости второго порядка п = 2 (вблизи Т = Тау » Т0)
а именно, при t = 0 должно быть (р1 (0) = j2(0), а угловая скорость изменятся скачком в соответствии с законом сохранения момента импульса (см. уравнение (4)). Из этих условий следуют уравнения для A1 и A2:
A1 cos(( T/4) = -A2 cos(2 T/4). (22) A1 (1 + ml) ■2 ( sin(( T/4) = = A2(1 - ml )2®2sin(ffl2 T/4). (23)
Однородные уравнения (22)-(23) для A1 и A2 имеют нетривиальное решение только при выполнении следующего условия:
(1 + mt)2 ( sin(( T/4) cos(2 T14) = = -(1 - ml)2w2 sin(2 T14)cos(( T/4). (24)
Это уравнение определяет период модуляции T (при заданном значении ml), соответствующий левой границе 2-го интервала параметрической неустойчивости. Уравнение (24) можно преобразовать к следующему виду, удобному для численного решения итерациями:
(q + 1)sin(av T/2) = (q - 1)sin(A®T/4),
(25)
где q зависит от глубины модуляции mt в соответствии с (14). Далее можно заменить T в аргументе синуса в левой стороне (25) на Tav + AT. Поскольку (OavTav = 2p, синус можно представить как - sin(wav AT/ 2). Тогда уравнение (25) принимает вид:
• ( AT/-» q - 1 • A ((Tav +AT)
sin(wav AT/ 2) = —-—г si^ av
q +1
4
(26)
Это уравнение относительно ЛТ можно решать численно с помощью описанной выше процедуры итераций. Самосогласованные решения для разных значений глубины модуляции т1 дают левую границу интервала п = 2. Для рассчитанного таким образом периода модуляции Т начальные условия, обеспечивающие стационарные колебания, находятся из уравнения (20), в соответствии с которым при произвольном
начальном отклонении (р1 (0) начальная уг-
ловая скорость должна быть равна (1(0) = щ\ш{щТ/ 4) (1(0).
Правая граница 2-го интервала дается уравнением, отличающимся от (26) противоположным знаком его правой части. Обе границы показаны на диаграмме 2 вместе с интервалами более высоких четных порядков, полученными путем аналогичных численных расчетов.
Отметим, что интервалы четных порядков (п = 2, 4, 6) при малой глубине модуляции т1 намного уже резонансов нечетных порядков. С ростом тг четные интервалы расширяются и становятся сравнимыми с интервалами нечетных порядков.
При малых и умеренных значениях глубины модуляции т1 << 1 из решения уравнения (26) (и уравнения с противоположным знаком для другой границы) можно получить приближенные аналитические выражения для обеих границ 2-го интервала:
T+ =
2
1 m ml2 4l
Tav = T0 +
9
- 3 - 3
+ 4 8
^i2To.
(27)
3
то есть Т- = Т0(1 - 8 т/), Т+= Т(1 + 8 т/). Как уже отмечалось выше, ширина этого
3
интервала неустойчивости Т+ - Т- =—т2Т)
пропорциональна квадрату глубин^1 модуляции.
5. ПЕРЕСЕЧЕНИЯ ГРАНИЦ ИНТЕРВАЛОВ ПРИ БОЛЬШОЙ ГЛУБИНЕ МОДУЛЯЦИИ
Как видно из рис. 2, при определенных значениях тг обе границы интервалов с п > 2 совпадают (можно считать, что они пересекаются). Это значит, что при таких т1 данные интервалы параметрической неустойчивости исчезают. Можно показать, что такие значения т1 соответствуют собственным периодам колебаний Т1 и Т2, которые относятся как 2 : 1, 3 : 1, и 3 : 2.
При первом пересечении (отношение 2:1) маятник совершает в точности половину собственного колебания с периодом Т1 в течение первой половины цикла модуля-
Рис. 7. Фазовая траектория и графики угловой скорости ф(/) и угла отклонения ф(/) для стационарных колебаний при пересечении границ третьего интервала неустойчивости
ции (см. рис. 7). На фазовой диаграмме изображающая точка за это время описывает половину меньшего эллипса (1-2) и затем перескакивает вниз на больший эллипс (2-3). За вторую половину цикла модуляции маятник совершает полное колебание с периодом Т2 = 71/2, так что изображающая точка проходит вдоль всего большего эллипса (34) и затем перескакивает на меньший эллипс по тому же самому вертикальному отрезку (4-5).
На протяжении следующего цикла модуляции длины маятника изображающая точка проходит сначала вдоль другой половины меньшего эллипса (5-6), а затем снова вдоль всего большего эллипса (7-8). Таким образом, на протяжении любых двух последовательных циклов модуляции изображающая точка проходит один раз вдоль всего меньшего эллипса и дважды вдоль большего эллипса, возвращаясь в итоге в начальную точку. Легко видеть, что такое колебание будет периодическим при любых начальных условиях. Это значит, что для соответствующих значений глубины модуляции Ш1 и периода модуляции Т рост амплитуды невозможен даже в отсутствие трения (интервал неустойчивости исчезает).
Аналогичные объяснения возможны и для других случаев пересечения границ интервалов неустойчивости на рис. 2.
6. ИНТЕРВАЛЫ ПАРАМЕТРИЧЕСКОГО ВОЗБУЖДЕНИЯ ПРИ НАЛИЧИИ ТРЕНИЯ
Если в системе есть трение, интервалы параметрического возбуждения сужаются и при достаточно сильном трении (ниже порога) исчезают совсем. Выше порога возбуждения приближенное выражение для границ основного интервала можно получить из соотношения (18), заменив в нем тг на
д/т,2 - (т1 с пороговым значением
(т1 )ш;п = р/(60 (см. [1]). Доказательство приведено в Приложении. Для третьего интервала можно использовать выражение
(19), подставляя в него ^т2 - (т1 вместо т1 с (т1 )тп = р/(20. Когда т1 уменьшается до своего порогового значения (т1 , соответствующий интервал параметрической неустойчивости исчезает.
Границы второго интервала параметрического возбуждения при наличии трения приближенно даются выражением (27), если
заменить в нем
т
на
т,4 - т )тш с по-
роговым значением (т1 )шЬ 2/(30 , соответствующим параметрическому резонансу второго порядка (см. Приложение).
Диаграмма на рис. 8 показывает границы первых трех интервалов параметричес-
кого возбуждения для Q = 3, Q = 5 и Q = 7 (а также в отсутствие трения). Отметим «островок» параметрического резонанса третьего порядка (п = 3) при Q = 7. Этот резонанс исчезает, когда глубина модуляции становится больше 48 %, и появляется снова, когда тг превышает приблизительно 66 %.
При наличии трения для данного значения глубины модуляции т1 могут существовать лишь несколько первых интервалов параметрического возбуждения (когда тг превышает пороговое значение). Отметим также, что, когда положение равновесия маятника неустойчиво из-за модуляции параметра, параметрическая раскачка возможна лишь при условии, что маятник уже совершает колебания хотя бы малой амплитуды. Если же начальные значения ф и ф в точности равны нулю, маятник будет оставаться в покое несмотря на модуляцию длины. Этим параметрический резонанс существенно отличается от обычного резонанса при непосредственном силовом воздействии, когда амплитуда растет, даже если первоначально система покоится в положении равновесия.
7. ЗАКЛЮЧИТЕЛЬНЫЕ ЗАМЕЧАНИЯ
В этой статье мы рассмотрели маятник, длина которого подвержена модуляции по
прямоугольному кусочно-постоянному закону. Такая модель дает удобный пример, в котором все особенности явления параметрического резонанса в нелинейной системе допускают ясное физическое истолкование. Рост амплитуды при параметрическом резонансе ограничен нелинейными эффектами, благодаря которым собственный период зависит от амплитуды. Увеличение размаха колебаний маятника при резонансе приводит к удлинению собственного периода. Система выбивается из резонанса, размах колебаний уменьшается, и условия резонанса восстанавливаются. Возникают переходные биения, постепенно затухающие из-за трения, и в конце концов устанавливаются колебания конечного постоянного размаха. Компьютерное моделирование движения системы существенно облегчает понимание физических причин ограничения роста амплитуды из-за нелинейных свойств маятника.
В статье показано, что границы интервалов параметрической неустойчивости маятника при прямоугольной кусочно-постоянной модуляции его длины можно рассчитать с помощью весьма скромных математических средств. Для этого достаточно рассмотреть стационарные периодические колебания бесконечно малой амплитуды, происходящие при совпадении периода
Рис. 8. Интервалы параметрического возбуждения при прямоугольной модуляции длины маятника в отсутствие трения при Q = 7, Q = 5 и при Q = 3
модуляции с одной из этих границ. Такие периодические колебания можно представить как чередование собственных колебаний с периодами Т1 и Т2 , соответствующими длине маятника 11 = 10 (1 + ш1) и длине 12 = /0(1 - т1) соответственно. «Сшивание» известных аналитических выражений
для этих собственных колебаний на примыкающих интервалах времени в моменты скачкообразных изменений длины маятника позволяет определить значение периода модуляции Т( тг) для соответствующего стационарного режима как функцию глубины модуляции.
Приложение
ГРАНИЦЫ ИНТЕРВАЛОВ НЕУСТОЙЧИВОСТИ ПРИ НАЛИЧИИ ТРЕНИЯ
РЕЗОНАНСЫ НЕЧЕТНЫХ ПОРЯДКОВ
Стационарные колебания, происходящие на левой границе основного интервала параметрической неустойчивости при наличии трения показаны на рис. 9 (сравните с рис. 1). Дважды на протяжении одного цикла модуляции угловая скорость резко увеличивается и дважды уменьшается. Увеличения больше уменьшений, так что в целом за счет модуляции длины получаемая маятником энергия превосходит отдаваемую энергию. Этот избыток компенсирует рассеяние энергии из-за трения, которое происходит при собственных колебаниях маятника в промежутках между резкими сдвигами груза вдоль стержня маятника.
Чтобы найти условия, при которых происходят такие стационарные колебания, запишем выражения j(t) и j(t) для соседних интервалов, на протяжении которых маятник совершает затухающие собственные колебания, и затем произведем «сшивание» этих выражений на границах интервалов. В противоположность случаю отсутствия трения (см. рис. 1), теперь фазовая траектория не симметрична относительно оси ординат (рис. 9). Выберем начало отсчета времени t = 0 в момент, когда груз сдвигается вниз и угловая скорость уменьшается. На протяжении интервала ( 0, Т/ 2) происходит затухающее собственное колебание с частотой w =w 0/-д/1 + m¡ . Его можно представить как суперпозицию колебаний синусоидального и косинусоидального типа с некоторыми константами A1 и B1:
j(t) = (A1sinw1t + B1cosw1t)e, ji(t)»(A1wcosw1t -B1w1 sin wt)e-g. (28)
Выражение для ji(0) приближенно справедливо при сравнительно слабом трении (g << w0), так как оно соответствует дифференцированию j (t) по времени при постоянном экспоненциаль-
Рис. 9. Стационарные колебания на левой границе интервала неустойчивости при наличии трения
ном множителе e ~gt. Действительно, при слабом затухании основной вклад в производную обусловлен быстро осциллирующими множителями sin (Of и cos wt в j1(t). Аналогично, на протяжении интервала (- Т/ 2, 0) график на рис. 9 представляет собой сегмент затухающего собственного колебания с частотой a:
<р2(t) = (A2 sin a + B2 cosa2t)e, (29)
<р2(t)»(A2a2 cos a -B2a2 sina2t)e~g . (30)
Для определения постоянных A1, A2, и B1, B2 воспользуемся условиями сшивания сегментов графиков и примем во внимание периодичность всего процесса. При t = 0 угол отклонения совпадает для <р1 и <р2, то есть j1(0) = р2(0). Из этого условия следует, что B1 = B2 . В дальнейшем эти совпадающие постоянные будем обозначать как B. Угловая скорость при t = 0 изменяется скачком в соответствии с законом сохранения момента импульса: (1 + mi )2<р1 = (1 - mi )2<р2 (см. уравнение (4)). Это условие дает следующее соотношение между A1 и A2: A2 = qA1 = qA (далее обозначаем A1 как A), где множитель q зависит от глубины модуляции mi в соответствии с (14).
В условиях стационарных колебаний на границе основного резонанса, как и других резонансов нечетных порядков n = 1, 3, ... в уравнении (6), условия периодичности имеют вид:
<1(Т/2) = -<2(-Т/2), (1 + ml)2р>1(Т/2) = -(1 -ml)2<<2(-Т/2). (31)
Подставляя р1 и <р1, р2 и <р2 в (31), получаем систему однородных уравнений для A и B:
(pS1 - qS2) A + (p€l - C 2) B = 0,
q(pCx+C2)A - (pqSx - S2)B = 0. (32)
где p = exp(-gt). В (32) использованы следующие обозначения:
C1 = cos(a Т/2), C2 = cos(a2 Т/2),
S1 = sin(a Т/2), S2 = sin(a2 Т/2). (33)
Система однородных уравнений (32) для A и B имеет нетривиальное решение только при равенстве нулю ее определителя:
2qC1C2 - (1 + q2 )S1S2 + q(p +1/p) = 0. (34)
Это условие существования ненулевого решения системы (32) дает уравнение относительно неизвестного периода модуляции T, который входит в (34) как аргумент синусов и косинусов в S1, S 2 и C1, C2, а также как аргумент экспоненты в p = e ~Jt. Искомые границы интервалов неустойчивости Т- и Т+ находятся как корни уравнения (34). Для получения решений этого трансцендентного уравнения преобразуем его к более удобной форме. Представим сначала произведения C1C2 и S1S 2 в (34) следующим образом:
„„ 1. ДюТ 1 ДюТ
C1C2 = 2 (cos 2 + cosa.^), S1S2 = ^(co^-^--cosюаvT). (35)
Затем, используя тождество cos a = 2 cos2 (a/2) -1, приводим (34) к следующему виду:
(д + Т/2) = +Л/(д -1)2 соб2 (ДдоТ/4) - р + ]/р - 2) . (36)
Чтобы найти границы основного интервала параметрической неустойчивости, нужно искать решение уравнения (36) относительно неизвестного периода Т в окрестности Т = Т0/2 » Тау/2. Если для заданного значения добротности Q ^ входит в р = ) глубина модуляции т1 превышает пороговое значение, уравнение (36) имеет два решения, которые соответствуют искомым границам Т- и Т+ интервала неустойчивости. Такие решения существуют, если подкоренное выражение в (36) положительно. Равное нулю значение соответствует порогу:
^соб^ДюТ/4) = р + - - 2. (37)
д Р
Для оценки порогового значения Q при малой глубине модуляции ш1 << 1 можно положить здесь q »1 + Ъш1 (см. уравнение (14) и соя(ЛаТ/4)»1. В правой части (37), в р = в~7', можно считать уг »770/2 = л/(20 << 1, так что р + — - 2 »(уТ)2 = (к/20)2 . Таким образом, для порога основ-
Р
ного параметрического резонанса получаем:
г> л / л к
^ » 6ш , (Ш^ » 60 . (38)
В условиях порога подкоренное выражение в уравнении (36) равно нулю. Оба его корня (границы интервала неустойчивости) совпадают. Так происходит, если косинус в левой части (36) равен
нулю, то есть когда его аргумент равен л/2:
„ Т л л 1
®ау- = -, или Т =-= -Тау ,
22 ®ау 2
так что пороговые условия (38) соответствуют точной настройке на резонанс: Т = Тау/2.
Для нахождения границ Т- и Т+ интервала неустойчивости представим Т в аргументе косинуса в левой части (36) как Т^/2 + ЛТ . Поскольку (отТш = 2л , этот косинус можно записать в виде -зт(юау ЛТ/2). Тогда уравнение (36) принимает вид:
ЛТ/2) = ±Л-1Й- 1)^ Л^2>г" +ЛТ)-q<Р=»1. (39)
q +1 у 4 р - '
При нулевом трении р = 1, и (39) совпадает с (14). Диаграмма на рис. 8 получена путем численного решения уравнения (39) относительно неизвестного ЛТ методом итераций. Границы интервалов неустойчивости высших нечетных порядков п = 3, 5, ..., рассчитываются аналогично. Для этого Т в уравнении (36) представляется как иТау/2 + ЛТ . Эти интервалы также показаны на рис. 8 для нескольких значений добротности Q. При больших значениях глубины модуляции ш1 эти границы почти сливаются с соответствующими границами в отсутствие трения.
Для получения приближенного решения уравнения (39), справедливого при малых значениях
глубины модуляции ш1 << 1 вплоть до квадратичных по членов, можно упростить подкоренное выражение в правой части (32), полагая q »1 + 3ш,, ^ -1)2 » 9ш 2, а значение косинуса принять рав-
1 ( ^
ным 1. Последний член в подкоренном выражении можно представить как —
160,
части синус можно заменить его малым аргументом, где юау = 2л/Тау . В результате получаем:
(ш, )2 . В левой
' Ш1П
ЛТ _ 1 , .2 или т = Т^'
■■ +—. ш,2 - (ш1 )2. , или Тт = Тау 2л V 1 Ш1П ' + 2
1 + — Ш - (ш, )2
л V Ш1П
(40)
/
При отсутствии трения (ш, )Ш1П = 0, и эти приближенные выражения для границ интервала неустойчивости сводятся к выражению (18). В условиях порога ш1 = (ш, )ш1п , и обе границы интервала сливаются при Т = Тау/2, то есть интервал неустойчивости исчезает.
При подстановке любого из двух корней Т- или Т+ уравнения (18) в (32) уравнения для А и В становятся эквивалентными, и из них можно найти только отношение А В. Тем не менее, соответствующие колебания характеризуются вполне определенной формой, которая определяется отношением амплитуд А и В синусоидальных и косинусоидальных функций, сегменты которых образуют рисунок стационарных параметрических колебаний (см. рис. 1 и 3).
РЕЗОНАНСЫ ЧЕТНЫХ ПОРЯДКОВ
Для описания стационарных колебаний, происходящих на границах интервалов неустойчивости четных порядков, можно использовать такие же выражения для р(г) и ф(г), как в (28) и (30). Условия сшивания графиков при г = 0 также не отличаются. Однако условия периодичности для момен-
тов - T/2 и T/2 в случае резонансов четных порядков отличаются от (31) противоположным знаком. В результате вместо уравнения (36) получаются следующие уравнения для границ интервалов неустойчивости:
(д +1) япф^ Т/2) = ±д/(д -1)2 э1п2 (ДюТ/4) - д(р + 1/р - 2). (41)
Для интервала 2-го порядка нужно искать решение в окрестности Г0 » Тау. Если для заданного значения добротности Q (Q входит в р = е-у) глубина модуляции Ш1 превышает пороговое значение, уравнение (41) имеет два корня Т- и Т+ , которые соответствуют границам интервала неустойчивости. Эти корни существуют, если подкоренное выражение в (41) положительно. Нулевое значение соответствует порогу, то есть значению (ш; )шп для заданного Q или Qш1n для заданного ш1:
М2 э1п2(ДюТа^4) = ^^. (42)
др
Пороговые условия выполняются при точной настройке на второй резонанс, то есть при Т = Тау. Чтобы оценить пороговое значение Q при малой глубине модуляции ш 1, можно положить здесь
д »1 + 3ш1, э1п(ДюТ/4)» ДюТ0/4 и Дю » ш1ю0. В правой части уравнения (42) в р = е~у можно
1 ( Р — 1)2
положить //» уТ0 = л^ << 1, так что р +---2 = —-»(уТ)2 = (л/Q)2 . В результате для порога
рр
2-го резонанса получаем:
2 [Г
^п > (Ш )ш;п . (43)
3ш2 у 3Q
Чтобы найти границы Т— и Т+ 2-го интервала неустойчивости, можно представить Т в аргументе синуса в левой части уравнения (41) как Тау + ДТ . Поскольку юауТау = 2л , этот синус записывается как - э1п(юау ДТ/2). В результате уравнение (41) принимает вид:
. юДТ _ 1 I ~ . 2 ДюТ, + ДТ) (р -1)2
^— = +- (д -1)2 эш2-^-- - д ——— (44)
2 д+1 4 р
Такая форма уравнения удобна для численного решения итерациями. Она использовалась для получения соответствующих кривых диаграммы на рис. 8. В случае нулевого трения р = 1, и уравнение (44) совпадает с (26). Для получения приближенного решения уравнения (44) при малой глубине
модуляции Ш1 можно упростить подкоренное выражение в правой части (44), полагая д »1 + 3ш;, (д -1)2 »(3ш; )2 и э1п[Дю (Тау + ДТ у4]» Д (оТт/4 = лш1/ 2. Последний член подкоренного выражения можно представить как (2/3 Q)2 » (ш1 )4 . . В левой части синус можно заменить его малым аргументом, где юау = 2л/Тау . В результате для границ 2-го интервала неустойчивости получаем следующие приближенные выражения:
—» m з m - m )4. , или Tm J* m 3 l.4-tm л4
T 4 V l u mrn m
av
1 m-. m - m )4.
4 V l min
Tv. (45)
Литература
1. Е.И. Бутиков. Маятник с модулируемой длиной. I. Параметрический резонанс // Компьютерные инструменты в образовании, 2013. № 3. С. 33-44.
2. S. M. Curry. How children swing // Am. J. Phys, 1976. № 44. P. 924-926.
3. P.L. Tea Jr., H.Falk. Pumping on a swing // Am. J. Phys, 1968. № 36. P. 1165-1166 .
4. W. Case. The pumping of a swing from the standing position // Am. J. Phys, 1996. № 64. P. 215-220.
5. M. A. Pinsky, A. A. Zevin. Oscillations of a pendulum with a periodically varying length and a model of swing // Int. J. Non-Linear Mech, 1999. № 34. P. 105-109.
6. D. Stilling, W. Szyskowski. Controlling angular oscillations through mass reconfiguration: a variable length pendulum case // Int. J. Non-Linear Mech, 2002. № 37. P. 89-99.
7. E. I. Butikov. Pendulum with a square-wave modulated length (simulation program) / http:// facultv.ifmo.ru/butikov/Applets/PendParSquare.html (дата обращения 28.08.2013).
PENDULUM WITH A MODULATED LENGTH. PART II.
REGIONS OF PARAMETRIC EXCITATION
Abstract
The first part of this paper [1] dealt with the principal peculiarities of parametric excitation of a rigid planar pendulum caused by a square-wave modulation of its length. The threshold of excitation was found and the role of nonlinear properties of the pendulum in restricting the resonant swinging was considered. In the second part of the paper the frequency boundaries of parametric instability are determined as functions of the modulation depth and the quality factor. Stationary oscillations at these boundaries are investigated.
Keywords: parametric resonance, ranges of excitation, phase locking, computer simulations.
Бутиков Евгений Иванович, профессор физического факультета СПбГУ,
ещене. butikov@gm.ail. сот.
© Наши авторы, 2013. Our authors, 2013.
