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

Баев Теория колебаниы 2015

.pdf
Скачиваний:
25
Добавлен:
12.11.2022
Размер:
2.85 Mб
Скачать

периодическому процессу. Кривая 2 на рис. 2.3.1 – пример неустойчивого предельного цикла.

Помимо рассмотренных выше особых точек фазовой плоскости типа «центр» и типа «седло», изображенных на рис. 1.6, часто встречаются и другие их разновидности.

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

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

x

y

 

y

 

 

 

t

x

x

 

а

б

в

 

 

x

y

 

y

 

t

x

x

 

г

д

е

 

 

Рис. 2.3.2. Особые точки фазовой плоскости типа «фокус»

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

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

161

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

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

Допустим, что для некоторой системы по мгновенным значениям построены фазовые траектории, изображенные на рис. 2.3.1.

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

x

y

y

 

 

t

x

x

а

б

в

 

 

y

x

y

t

x

x

г

д

е

Рис. 2.3.3. Особые точки фазовой плоскости типа «узел»

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

162

163
ω ´ ( ,
где и

констатируем, что на фазовой плоскости есть один устойчивый предельный цикл (кривая 1) и один неустойчивый (кривая 2).

Если начальные условия в системе будут таковы, что изображающая

(кривая

I H

 

ее точка, определяемая начальными значениями координаты xи

скорости

,

попадает внутрь неустойчивого предельного цикла

 

2), то

в системе возникает затухающий колебательный

устойчивого равновесия. Если же начальные значения x и I будут таковы, что изображающая точка попадет в область, находящуюся

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

снаружи кривой 2, то в результате переходного процесса возникнут автоколебания – изображающая точка будет двигаться по устойчивому предельному циклу (кривая 1).

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

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

2.3.2. Метод усреднения по быстрым осцилляциям

Рассмотрим движение системы, находящейся в силовом поле, определяемом потенциальной функцией U и, кроме того,

испытывающей воздействие силы f, меняющейся с частотой ω:

2.3.1

cosω sinω ,

1 функции, зависящие только от координат.

Пусть T период колебаний, которые совершает система, двигаясь в силовом поле U. Сделаем два допущения, в рамках которых рассмотрим движение системы. Во-первых, будем считать частоту изменения силы достаточно большой, чтобы выполнялось условие

2.3.2

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

относительно этих возмущений. Заметим, что при этом величина не

предполагается малой по сравнению с силой, действующей в поле U.

Сначала ограничимся одномерным движением описываемым

уравнением

 

 

&6

 

 

 

 

2.3.3

 

 

#HL &H .

Представим решение этого уравнения в виде суммы двух

функций

 

H j ξ

2.3.4

и потребуем, чтобы

Õ

2.3.5

 

 

 

 

ξ 0;

 

 

HÕ j ;

2.3.6

здесь символы

и означают величины, усредненные по интервалу

времени

Õ

 

 

 

 

ξ HÕ

 

2.3.7

 

 

a

 

 

(

ω .

Выражение (2.3.4) и условия (2.3.5), (2.3.6) означают, что движение системы представляется состоящим из двух компонент.

Одна компонента,( , описываемая функцией j , мало меняется за время > так что с достаточной степенью точности считается справедливым выражение (2.3.6). Можно сказать, что функция

j

описывает усредненное по быстрым осцилляциям, плавное движение системы. Вторая компонента решения (2.3.4), представляемая

функцией ξ , как раз и является теми быстрыми осцилляциями, о которых шла речь выше и которые накладываются на усредненное

плавное движение системы. Эта |компонентаξ | q 1 вызвана силой , причем согласно принятому допущению .

164

Подставив функцию (2.3.4) в уравнение (2.3.3), разложив его

правую часть в степенной ряд по величине

 

и ограничившись

линейными членами

разложения,

преобразуем

его к

 

 

 

 

 

¨

 

 

следующему виду:

L

6

 

 

 

 

 

 

L

 

 

 

2.3.8

#j #ξ

j

ξ ƒ

 

j, ξ ƒ.

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

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

Выделим из уравнения (2.3.8) часть, описывающую быстро

осциллирующее движение:

L

j, .

2.3.9

 

 

В правую часть этого уравнения вошла только одна из трех быстро осциллирующих функций правой части исходного уравнения (2.3.8),

поскольку по сравнению с функцией

j,

две

другие,

пропорциональные

являются величинами

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

более высокого порядка

малости и поэтому ξ,оказывают значительно меньшее влияние на

быстро осциллирующую компоненту, чем функция

 

 

.

которой

Интегрируя уравнение (2.3.9)

с функцией

 

(2.3.1), в

 

j,

 

 

 

(

 

 

 

 

 

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

величину

 

 

можно считать постоянной, так как она мало меняется за

время

 

>, j

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

найдем функцию

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

.

 

 

 

 

 

 

 

2.3.10

Õ Õ

 

 

 

 

:

 

 

ω

 

 

 

 

 

 

 

ξ 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω . Поскольку

Усредним теперь уравнение (2.3.8) по времени

 

 

 

 

то после

операции усреднения получим(уравнение для

функции

j

 

6

 

„„„„„

6

 

1

 

„„„„„„

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ƒ j

ξ ƒ j

ω

ƒ

 

 

2.3.11

165

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

2.3.12

где

 

 

ƒ ƒ

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

2

 

$

1

 

2

%.

2.3.13

 

6 2

ω

6

4 ω

2

2

 

 

 

 

 

 

 

 

 

 

 

Функцию

можно

интерпретировать

 

как

 

эффективную

потенциальную6cэнергию

системы, включающую в

 

себя

быстро

осциллирующие силы, действующие на систему. Если выражение (2.3.13) для эффективной энергии объединить с выражением (2.3.10),

то получим следующее выражение:

# 2

 

 

 

 

 

 

 

 

6 2

'

2.3.14

 

 

ξ .

Это

выражение означает,

что усредненное

по быстрым

 

6

 

 

 

 

осцилляциям движение системы происходит так, как если бы помимо

поля

 

на нее действовало

бы еще

и дополнительное поле,

определяемое кинетической энергией быстро осциллирующего движения.

 

С помощью выражения (2.3.14) полученные результаты легко

обобщаются на системы со многими степенями свободы: если

2.3.15

то

( 2

,

 

1

 

 

 

,

 

 

 

 

,

 

„„„„„„„

 

6 2

 

 

' '

2.3.16

 

$ξ ξ

%.

 

1

 

 

 

 

В качестве примера рассмотрим систему, представляющую маятник с колеблющейся вертикально точкой подвеса. Вертикальные

166

167
12 N 2DBγ N 1,

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

подвеса приводят к появлению в уравнении ее движения переменной силы

 

 

#B γ cosγ sinφ;

φ

2.3.17

здесь и

 

масса и длина маятника соответственно;

– его угол

отклонения от положения равновесия.

 

#

B

 

 

В соответствии с выражением (2.3.13) эффективная потенциальная функция рассматриваемой системы принимает вид

 

#DB cosφ

2γ2

2

2.3.18

 

4DB

sin φ .

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

определяемом функцией . Из выражения

 

 

& $

#DB ,1

γ

cosφsinφ 0

2.3.19

2DB

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

– это точка φ 0

и, если выполнено условие

 

 

 

точка

 

2DB

N 1,

 

 

 

 

 

2.3.20

 

 

 

γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

φ

arccos

2DB

!.

2.3.21

 

 

 

γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

при условии

2.3.22

а ее координата лежит в интервале

 

 

 

 

 

 

^ π.

2.3.23

 

 

 

4 N φ

 

φ

Таким

образом, при вертикальных колебаниях

точки подвеса

0

 

сохраняется очевидная

 

точка устойчивого равновесия

маятника

 

 

 

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

условии (2.3.22) появляется еще одна точка устойчивого равновесия, которая в зависимости от значения параметров маятника может иметь значения из интервала (2.3.23). Задача решена.

2.3.3. Укороченные уравнения

Пусть движение системы описывается уравнением вида

2.3.24

 

 

|ε| q 1

 

 

HL • H Ð H, H ,

Если

 

.

 

 

 

в котором

 

 

 

 

 

положить

 

то уравнение (2.3.24) перейдет в уравнение

гармонического

осциллятора с известным общим решением. Для

 

ε 0,

 

 

аналогий, рассуждая следующим образом. Из-за малости параметра ε можно допустить, что общее решение данного уравнения внешне

анализа этой математической модели снова воспользуемся методом

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

величина

 

. То же самое скажем и о фазе колебаний, представив ее как

вторую

неизвестную функцию искомого решения.

 

 

ε

 

 

Итак, решение уравнения (2.3.24) будем искать в виде

2.3.25

 

 

 

H cosφ .

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

168

осциллятора, а именно: потребуем, чтобы искомое решение удовлетворяло требованию

H ωsinφ.

2.3.26

Дифференцируя функцию (2.3.25) и учитывая условие (2.3.26),

придем к выражению cosφ φsinφ ωsinφ 0.

2.3.27

Дифференцируя далее функцию (2.3.26) и подставив получившуюся вторую производную в исходное уравнение (2.3.24),

получим еще одно уравнение:

ωsinφ ωφcosφ ω cosφ ε cosφ, ωsinφ . 2.3.28

Рассматривая

выражения (2.3.27)

и

(2.3.28)

как

 

уравнения

 

 

 

 

 

 

 

 

и

 

и

решая

 

 

H

 

 

 

относительно

производных

 

 

 

 

 

их,

 

найдем

дифференциальные

уравнения для искомых функций

 

 

и

 

 

:

 

 

 

φ'

 

 

'

 

ε

1

 

 

 

 

 

 

 

 

 

'

ε

 

 

 

 

 

 

 

 

 

 

, φ ;

2.3.29

ω

cosφ, ωsinφ sinφ ω

 

ε

 

 

 

 

 

 

 

 

 

 

 

ε

2

, φ ; (2.3.30)

φ' ω ω

cosφ, ωsinφ cosφ ω ω

 

Первое из этих уравнений, как и следовало ожидать, показывает,

что амплитуда

колебаний

 

 

 

 

меняется

медленно, так

как

ее

производная имеет

порядок

 

. Следовательно, можно ожидать,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что за одно колебание, т.е.

за время, в течение которого фаза

 

 

 

ε

 

 

 

 

 

 

 

 

 

 

 

движения

изменится на величину, равную 2

 

, амплитуда а и характер

 

 

 

φ

изменятся мало. Поэтому

не

 

должен привести к большой ошибке

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

переход в правых частях выражений (2.3.29), (2.3.30) к их

усредненным за период значениямε

1

 

 

' ω

„„„„„„„

2.3.31

 

;

169

 

 

 

 

 

 

 

ε

2

 

 

 

 

 

φ' ω ω

„„„„„„„

где

 

 

 

 

,

„„„„„„„

 

0

cosφ, ωsinφ sinφ?φ;

1

2e M

 

1

2d

 

 

 

„„„„„„„

 

0

cosφ, ωsinφ cosφ?φ.

2

2e M

 

 

1

2d

 

 

 

2.3.32

2.3.332.3.34

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

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

uuuuuuu

2.3.35

 

0.

 

 

Если оно не имеет действительных решений, то это, очевидно, означает, что в системе стационарные колебания невозможны.

Проиллюстрируем сказанное на примере следующего уравнения

HL ω H εαH2.

2.3.36

Его укороченные уравнения имеют вид

2.3.37

' 0;

2

 

2.3.38

φ' ω .

170