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

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

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

и вычислим ее производную, учитывая уравнение (2.2.14):

 

 

&M

&H

1

D!.

 

2.2.18

 

 

!

ω

 

Из последнего выражения видно, что энергия системы не будет

возрастать, если

 

&lnω

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.2.19

 

 

 

ω º 0.

 

 

 

В свою очередь, если полная энергия

системы в процессе ее

устойчивым.

остается

 

|&H/&τ|

то

|H|

должны

оставаться

движения

ограниченной,

 

ограниченными

и величины

 

 

 

и

 

т.е.движение

будет

Вернувшись к исходной независимой переменной

, в конечном

итоге убеждаемся в справедливости следующего

утверждения

 

 

 

 

Теорема

2.2.1

Для

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

(2.2.14)

достаточно, чтобы для любого момента времени коэффициенты этого уравнения удовлетворяли неравенству

&lnω

D º 0.

2.2.20

&

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

Возвращаясь к примерам из предыдущего параграфа, легко

проверить, что коэффициенты уравнения (2.2.11), в котором

D 1 ;

131

ω 1 ,

2.2.22

удовлетворяют условию (2.2.20), а в уравнении (2.2.5) с неустойчивым решением оно нарушается.

Применим теорему 2.2.1 об устойчивости к уравнению

 

 

 

 

 

HL Q+H Q+H 0,

2.2.23

в котором

 

 

 

 

ω Q

 

 

 

 

;

 

2.2.24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.2.25

 

 

 

 

 

 

D Q+.

 

Для этого уравнения условие устойчивости принимает вид

 

 

 

 

2 Q

+

p 0

2.2.26

 

 

 

 

 

α

 

 

 

 

 

 

и, начиная с момента времени

1

 

 

 

α

 

 

 

 

 

 

 

W

 

 

 

 

2.2.27

 

 

 

¤

 

 

2 ln g2 h,

Второй

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нарушается, что означает потерю системой устойчивости, если она

работает время

 

W.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

решая уравнение, связан с уравнением вида

2.2.28

 

 

 

 

HL

1 φ H 0,

 

|φ| N 1 φ » 0

 

» ∞

 

 

в котором

 

 

и

H &

1

 

 

 

 

 

 

 

.

 

 

 

 

 

 

при

 

 

 

 

 

 

 

 

Умножим его на , перепишем в виде

 

 

 

 

 

& ,2 H H - φHU

2.2.29

 

 

 

 

 

 

 

132

 

 

 

 

 

и проинтегрируем по времен

 

H

 

H

 

 

1 φ

 

 

 

1

 

 

 

 

? ,

 

2

2

 

 

 

ž 2

;

φ H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

H

 

H

 

 

 

 

 

 

 

 

 

 

 

ž

@ 2

2 1 φ A ¤ 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из выражения (2.2.30) следует неравенство

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

1

|

 

|

 

? ,

 

 

2

 

1 φ

 

ž

2

;

φ H

которое можно усилить, переписав в виде

 

 

 

 

 

 

 

4

 

ž 2 ;

φ 4

 

? .

 

 

 

 

 

H

 

 

 

 

 

 

|

 

| U

 

 

 

 

Если ввести новые величины

¢ H4 ;

P 2|φ|,

то неравенство (2.2.33) примет вид

¢ ž ; q r? .

133

2.2.30

2.2.31

2.2.32

2.2.33

2.2.342.2.35

2.2.36

Воспользуемся

леммой,

которую

приведем

без

доказательства [12].

 

 

 

 

 

 

 

Лемма 2.2.1. Пусть

и

 

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

функции,

а С –

положительная постоянная,

причем для любого

 

выполняется

 

¢

P

 

 

¤ 0

 

 

неравенство (2.2.36), тогда

 

 

 

 

 

 

 

 

 

 

 

2.2.37

 

 

 

 

 

 

 

¢ žexp š; P& ›.

 

Из леммы и неравенства (2.2.33), которое сводится к неравенству

(2.2.36), следует, что

 

 

 

 

 

 

 

H2

 

 

X

 

2.2.38

 

 

 

 

4 e

 

 

0

 

žexp sM 2|z'|& t.

 

Таким образом, получен еще один критерий устойчивости системы в виде следующей теоремы

Теорема 2.2.2. Для устойчивости системы, описываемой

уравнением (2.2.28), достаточно, чтобы интеграл

2.2.39

¬ ;|φ|&

 

 

cходился при » ∞.

2.2.3. Общие соотношения для решений линейного дифференциального уравнения с переменным коэффициентом

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

 

&

 

H

 

H 0,

2.2.40

&

 

 

 

 

 

 

 

в котором произвольная функция.

134

Можно показать [13], что существует такая фундаментальная пара решений уравнения (2.2.40)

χ |χ | ;

2.2.41

χ |χ | ,

2.2.42

что любое частное решение уравнения (2.2.40) можно представить в виде

 

 

 

 

 

 

|

|

 

 

 

 

, 2.2.43

 

 

χ

 

χ

 

χ

 

cos

ψ

 

 

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

 

|

|

 

 

 

 

1

 

 

χ

 

 

|

|

 

 

 

 

|χ|3

0;

2.2.44

 

 

χ

 

 

 

 

ψ

 

1

 

 

 

 

 

 

 

|χ| .

 

2.2.45

 

 

 

 

 

 

Эти уравнения играют ключевую роль в одном из самых красивых и эффективных разделов теории пучков заряженных частиц, а именно: в теории огибающих пучков [13], элементы которой будут изложены ниже.

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

 

 

χ

 

 

 

 

;

2.2.46

 

 

 

χ

2.2.47

 

 

 

 

 

 

 

 

135

 

 

или в матричной форме

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&k

 

 

 

 

 

2.2.48

где

 

 

 

 

& <k,

 

 

 

 

 

 

 

 

 

uv

 

Ê

 

 

 

 

 

 

 

 

 

 

 

 

wPx ;

 

 

2.2.49

 

 

 

 

 

 

 

k

1

 

 

 

 

 

 

<

 

0

 

 

 

 

2.2.50

с начальным условием

 

 

0!,

 

uv

 

 

 

Ê

 

 

 

uuuv

 

 

 

 

 

 

 

 

 

 

0

 

 

2.2.51

 

 

 

 

 

y 0 w P

 

x y .

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

Решение

матричного уравнения

 

2.2.48

также

может быть

представлено в матричном виде:

 

 

 

 

 

 

 

 

 

 

 

 

uv

 

 

 

 

 

 

uuuv

 

2.2.52

 

 

 

(

 

 

y

( y ,

 

,

 

 

 

 

2.2.48

 

 

 

 

 

 

Ê

 

 

 

 

 

 

 

 

 

 

где

матрица

 

формируется из независимых частных решений

 

 

исходного уравнения

 

 

 

χ

 

χ !

 

2.2.53

 

 

 

 

 

(

 

 

 

 

 

 

 

 

 

 

 

 

χ

χ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с начальным условием

( 0 g1

 

 

0h ¬.

 

2.2.54

 

 

 

 

 

 

 

 

 

Несложно убедиться в том, что

0

 

 

1

 

 

 

 

 

 

 

 

 

&

det(

0

 

2.2.55

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

136

 

 

 

 

 

 

и, следовательно,

 

 

det( const,

2.2.56

 

 

 

 

а из начального

условия

 

находим значение

постоянной

величины

 

 

приходим к следующему замечательному

 

так что 2.2.54

 

 

свойству

исходного уравнения (2.2.40): у этих уравнений

 

const 1,

 

 

 

 

 

 

 

 

det( 1.

2.2.57

2.2.4. Уравнение Хилла

Вернемся к линеаризованному уравнению маятника с

колеблющейся точкой

подвеса

его

 

 

 

и

отметим

одну важную

особенность

этого уравнения,

 

коэффициент

 

периодическая

2.2.2

 

 

 

 

 

 

функция.

 

 

 

 

 

(

 

 

 

 

 

В общем

случае,

уравнение

 

:

,

в

котором

функция

 

 

 

является периодической с

периодом

 

 

 

 

 

 

 

 

 

2.2.40

 

,

 

 

 

 

 

 

 

(

 

 

 

 

 

2.2.58

 

 

 

 

 

 

 

 

 

 

 

 

называется уравнением Хилла.

Хотя первая публикация об этом уравнении появилась в 1868 году, до сих пор выходят посвященные ему работы, а количество опубликованных к настоящему времени монографий и статей об уравнении Хилла огромно. Это объясняется, с одной стороны,

обширными приложениями уравнения Хилла, а с другой стороны

с

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

 

многочисленными приложениями изложена в монографии [14].

 

 

также

χ

 

 

Непосредственно из определения уравнения Хилла, в частности,

χ (

 

 

есть решение этого уравнения, то функция

следует, что если

 

должна быть его решением.

В основе теории уравнения Хилла лежит теорема Флоке [14]. Теорема Флоке. Существуетφ ( zединственная( , пара независимых

решений уравнения Хилла и удовлетворяющая условиям

137

 

 

 

 

 

 

 

 

φ ( λ φ );

 

 

 

(2.2.59)

 

λ

£

 

 

 

 

φ ( λ φ ),

 

 

 

(2.2.60)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

постоянные,

которые

 

можно было

бы назвать

собственными значениями функций

 

 

 

 

и

 

 

 

.

функциями

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Хилла называется

 

Эта пара решений уравненияφ

 

φ

 

 

 

 

Флоке, через которые можно выразить любое решение уравнения

 

Константы

 

 

 

χ С φ С φ .

 

2.2.61

 

 

и

 

связаны

 

друг

 

с

другом,

в чем несложно

убедиться.

Действительно, если уравнения

 

 

 

 

 

 

 

 

λ

£

φL

 

 

 

 

 

φ

 

0;

 

 

 

2.2.62

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

0

 

 

 

 

2.2.63

 

 

 

 

 

 

 

 

φL

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

умножить соответственно на

φ

и

 

φ

и вычесть одно из другого, то

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

φ φ

?

 

 

 

φ'φ

0,

2.2.64

 

 

 

 

 

φ φ

 

 

φ'φ

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

из которого следует, что

 

 

?

 

φ φ' const.

2.2.65

 

 

 

 

 

 

3 φ'φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С другой стороны, согласно теореме Флоке,

 

 

 

 

 

 

 

 

3 ` λ λ 3 λ λ const,

2.2.66

так что

 

 

 

 

 

 

 

 

λ

138λ 1.

 

 

 

 

 

2.2.67

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из теоремы Флоке следует, что

 

 

2.2.68

 

φ, %(

λ, φ, ,

и поэтому решения уравнения Хилла будут устойчивыми, если

 

 

 

 

 

 

Ìλ,Ì 1.

 

2.2.69

сопряженными:

2.2.67

 

являются комплексными и, как это

В общем случае, числа

£

и

 

видно

из выражения

 

 

£, они

должны быть

комплексно-

 

 

 

 

 

£ £W ,

 

2.2.69

2.2.70

будут

располагаться

на

единичной

 

£

£

 

так что при выполнении условия устойчивости

 

числа 1 и

 

2

 

 

 

 

 

 

 

окружности комплексной

плоскости, рис. 2.2.2,а. В частности, если

 

2.2.71

 

 

 

 

λ

λ

1,

 

то решения уравнения Хилла будут описывать устойчивые периодические движения с периодом Т при λ=1 и 2Т при λ= 1.

jy

jy

 

 

 

 

 

 

 

j

j

 

 

 

°λ

 

 

 

 

°

 

°1 x

 

λ

1

x

λ

 

 

λ

 

 

°

 

а

 

б

 

 

Рис. 2.2.2. Случаи с устойчивыми (а) и неустойчивыми (б) решениями уравнения Хилла

139

2.2.73

Если же числа λ и λ «соскакивают» с единичной окружности, как на рис. 2.2.2,б, то система теряет устойчивость. Потеря устойчивости системы, описываемой уравнением Хилла, называется параметрическим резонансом.

Название резонанса объясняется тем, что он может инициироваться периодическим изменением во времени параметра системы, как в случае с уравнением (2.2.2), описывающим в линейном приближении маятник, точка подвеса которого периодически колеблется.

2.2.5.Диаграмма устойчивости

Вэтом параграфе решается задача построения диаграммы устойчивости уравнения Хилла, точнее, одной из его разновидностей,

имеющей вид:

&

 

H

 

 

 

 

δ εcosτ H 0

2.2.72

 

 

 

 

 

 

 

и известной под названием уравнения Матье.

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

Заметим также, что в общем случае период коэффициента в уравнении Хилла или, как в данном частном случае, в уравнении Матье, может быть произвольным, однако с помощью простого

преобразования независимого переменного

τ (

легко перейти к уравнению вида 2.2.72 с периодом ` 2π. Напомним, что диаграммой устойчивости принято называть

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

140