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

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

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

 

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

2.3.39

где

Žγ ω φ ,

 

 

2.3.40

 

γ 3εα.

Если по аналогии с решениемφ уравнения гармонического осциллятора под производной понимать частоту колебательного движения системы, то выражение (2.3.39) представляет ее амплитудно-частотную характеристику. Эта характеристика изображена на рис. 2.3.4.

a

T 0

^ 0

 

 

j jj

ω φ

Рис. 2.3.4. Амплитудно-частотная характеристика

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

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

гак

называемой

мягкой

εαH

 

,

 

 

. У такой системы

 

ω H

системы, у которой

 

внешняя сила, равная

 

2

 

вычитаясь из

собственной силы

системы, равной

,

 

 

 

α ¤ 0

 

 

 

 

ослабляет жесткость системы, так как уменьшает силу,

возвращающую систему в положение ее устойчивого равновесия. В

области II располагается фрагмент, представляющий амплитудно-

силой

α N 0

 

частотную характеристику так называемой жесткой системы, у

которой

 

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

системы, увеличивая в сумме силу, возвращающую систему в положение равновесия.

171

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

(2.3.24) к новой независимой переменной

τ ω

. Тогда уравнение

примет следующий вид

HL H µ

 

 

2.3.41

 

H, H ,

где

ε

,

 

2.3.42

 

µ ω

 

 

 

τ

 

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

переменной .

 

 

Представим решение уравнения (2.3.41) в виде

2.3.43

 

H cosτ sinτ,

 

 

Ì Ì q Ì Ì

 

где

τ .и τ - медленно меняющиеся функции, так что | | q | | и

 

Потребуем, чтобы выполнялось условие

2.3.44

 

H sinτ cosτ,

тогда, если продифференцировать функцию (2.3.43) и объединить получившееся выражение с выражением (2.3.44), то получим выражение вида

 

2.3.45

coa‡ sin‡ 0.

Подставив (2.3.43) в уравнение (2.3.41), после несложных преобразований, аналогичных проделанным выше, и с учетом выражений (2.3.44) и (2.3.45) придем к системе из двух

дифференциальных уравнений первого порядка

2.3.46

µ cosτ sinτ, sinτ cosτ si%τ;

' µ cosτ sinτ, sinτ cosτ cosτ,

2.3.47

172

 

из которых, после операции усреднения, получаем укороченные

уравнения следующего вида: ' φ , ;

2.3.48

'

 

 

2.3.49

ψ , ,

где

 

 

 

1

e

, , τ sinτ&τ;

2.3.50

φ ,

; µ

1

 

 

 

e

, , τ cosτ&τ.

2.3.51

ψ ,

; µ

Система укороченных τуравнений (2.3.48), (2.3.49) не содержит в явном виде переменную , что во многихτ случаяхτ позволяет ее проинтегрировать и получить функции и , являющиеся медленно меняющимися амплитудами искомого решения.

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

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

173

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

2.3.4. Метод последовательных приближений

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

аппроксимируя функцию следующими суммами:

 

 

H H H

.

9

1

#

H H ~

 

H H D H H H b

#

 

 

f

 

 

# f #

 

#

 

 

2

,

 

 

 

 

, ,#

, ,#

 

 

 

 

 

 

 

 

 

 

 

 

(2.3.52)

Эта функция состоит

из

двух

частей. Часть,

представленную

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

Перейдем

от

 

произвольных координат

 

 

 

 

(i=

1,2,…,n)

к

нормальным координатам

 

 

( = 1,2,…,n),

введенным в линейном

 

 

H

 

 

 

новых

приближении

(см.

 

параграф

2.1.1).

Функция

Лагранжа в

 

 

 

 

5

α

 

 

 

 

 

 

 

 

 

координатах будет иметь вид

 

 

 

 

 

 

 

 

 

 

 

 

9 g 5

 

ω 5 h λ 5 5 5 µ 5 5 5 .

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

g

 

 

 

g

 

g

2

 

 

 

 

 

 

 

 

 

, ,g

 

 

 

 

, ,g

 

(2.3.53)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Переходя далее к уравнениям Лагранжа, получим следующую

систему из n нелинейных дифференциальных уравнений движения:

 

 

 

 

L

ω

 

5

 

 

 

 

 

 

 

 

2.3.54

 

 

5

 

 

 

5, 5 α 1,2, … , % .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

174

Как это уже делалось ранее, в упомянутом выше параграфе,

представим

искомые

 

решения уравнений

2.3.54

в

виде сумм

функций:

 

 

,

2.3.55

 

5

5

 

5

5

l

 

α 1,2, … , %

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( ) ( )

(2)

 

 

 

 

 

 

 

считая, во-первых, что каждое следующее слагаемое суммы (2.3.55), являясь уточнением предыдущего приближения, удовлетворяет условию

³5( )³ q ³5( )³ α 1,2, … , % ,

2.3.56

и, во-вторых, что функции 5(1)α (α= 1,2,…,n) являются решениями

однородных линейных уравнений

 

 

 

ω

5

α

0

α 1,2, … , %

2.3.57

@

& '

(1)

 

2

 

 

 

 

%

α

 

 

 

 

 

т.е. имеют вид

 

5α(1) \%cos ω% θ% .

2.3.58

 

 

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

В соответствии с условием (2.2.56), оставляя в правых частях уравнений (2.2.54) лишь члены первого порядка малости, получим неоднородные линейные дифференциальные уравнения для функций

второго приближения:

 

α

 

α

 

 

'

 

 

ω

5

 

w5

 

x α 1,2, … , % , 2.3.59

@

& '

(2)

 

(1)

, 5

(1)

 

2

 

 

 

 

 

 

%

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

175

 

 

 

правые части которых будут содержать тригонометрические функции, например, такого вида:

5

 

5

 

cos ω θ cos ω

θ

2

€ C ,

 

( )

 

( )

 

α

α

β

β

1

 

 

 

 

 

 

 

 

3 cos g ωα ωβ

θα θβh ;

 

(2.3.60)

где

 

 

 

 

2.3.61

 

 

 

 

 

 

 

 

 

 

 

 

 

 

O cos g ωα ωβ

θα θβh.

 

2.3.62

Таким

 

образом, в

правые

части уравнений

(2.3.59) войдут

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

ее линейного приближения.

 

 

 

 

 

 

 

ω

 

 

 

α

 

 

Решив

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

убеждаемся

в

том,

 

 

 

 

 

 

 

 

что

во

×

 

втором

 

 

 

 

 

Öω

 

s ω

 

 

 

α, β

 

приближении на нормальные колебания

с

частотами

 

 

 

 

 

(

=

 

1,2,…,n) наложатся колебания с частотами

 

ω

 

 

 

 

 

 

(

 

=

 

 

 

 

 

 

 

 

 

 

 

1,2,…,n),

в частности, с

удвоенной частотой

2

 

 

 

и

 

 

 

нулевой

Частоты колебаний ω s ω называются комбинационными.

частотой, соответствующей постоянному смещению системы.

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

Суть ее состоит в том, что уже в первых приближениях среди комбинационных частот могут появлятьсяω частоты, совпадающие с собственными частотами системы . Из теории линейных систем известно, что тригонометрические функции в правой части уравнений движения с такими частотами вызывают резонансы. Если рассматриваемая система консервативная, то в ней не может происходить самопроизвольное нарастание амплитуды колебаний. Такую неустойчивость можно назвать математической, так как она вызвана не физическими процессами, протекающими в системе, а погрешностями метода решения уравнений движения. Члены приближений, вызывающие математические неустойчивости, принято называть сингулярными. Следовательно, чтобы получить верное

176

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

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

сингулярных членов.

 

 

 

Проиллюстрируем

сказанное

конкретным

примером,

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

сведено

к рассматриваемому

путем разложения

нелинейной

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

2.3.63

 

HL ω H αH βH2.

Ограничимся случаем малых значений величин и

Представим

искомое решение уравнения в виде суммы трех последовательных

приближений:

α

β.

 

 

считая, что

H H( ) H( ) H(2) ,

2.3.64

ÌH( ) Ì q ÌH(1) Ì

2.3.65

 

и

ÌH(2) Ì q ÌH( ) H( ) Ì.

2.3.66

 

Тогда в качестве первого приближения можно взять решение уравнения гармонического осциллятора (еще раз хочется обратить внимание на очередной случай использования этого решения в очередной методике исследования математической модели)

177

 

 

 

H&1' cosω ,

2.3.67

в котором

 

собственная частота колебаний

рассматриваемой

 

 

системы.

 

нелинейной ω

 

 

ω

В соответствии с вышесказанным, собственную частоту колебаний

 

 

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

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

 

 

 

ω ω ω( ) ω( ),

 

2.3.68

и также считая, что

( )

 

и

 

 

( )

 

 

ω

( ).

первыми двумя

Сначала

ограничимся в суммах

 

 

,

 

 

Ìω

Ì q ω

 

Ìω

 

Ì q ω

 

 

 

слагаемыми. Введя их в исходное

уравнение (2.3.63) и собрав члены

 

2.3.64

 

2.3.68

 

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

 

 

HL

 

ω H

 

 

2

 

 

2

 

cos2ω 2 ω ω

 

cosω . 2.3.69

 

 

( )

 

 

( )

 

 

α

 

 

α

 

 

 

 

 

 

 

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В правой части этого уравнения, т.е. уже во втором приближении,

появился сингулярный элемент, пропорциональный функции

 

вызывающий

 

 

расходимость

 

 

 

решения

 

(математическую

 

 

 

 

 

 

 

 

cosω ,

ω

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

неустойчивость). Чтобы ее устранить, очевидно, достаточно положить

 

( )

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решив при этом условии уравнение (2.3.69), найдем функцию

второго приближения:

 

 

 

 

α

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

 

cos2ω .

 

2.3.70

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Перейдя к полным суммам из трех слагаемых (2.3.64), (2.3.68), введя их в исходное уравнение (2.3.63) и собрав элементы третьего порядка малости, получим уравнения для функции третьего приближения:

178

HL ω H

 

 

 

 

 

 

 

 

 

β

 

 

α

 

- cos3ω

 

 

 

 

 

 

,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

(2)

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

,2ω ω

 

 

 

 

 

 

 

3

 

 

 

 

 

 

2.3.71

 

 

 

 

 

4 β- cosω .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

пропорциональный функции

 

 

 

 

 

 

 

 

 

. Для его устранения приравняем

нулю коэффициент

 

 

при

 

этой функции, что дает поправку к

 

 

 

 

 

cosω

 

 

 

 

 

 

 

 

 

 

 

 

собственной частоте колебаний нелинейной системы

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

- .

 

2.3.72

 

 

 

 

 

 

,

 

 

12ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

Решив уравнение (2.3.71) при условии (2.3.72), найдем функцию

третьего приближения

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

β

 

 

 

 

 

H

&3'

 

 

 

2

 

 

 

2

 

 

 

 

2.3.73

 

 

16ω0

3•0 2 cos3ω .

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

При необходимости этот процесс уточнения вида функции

 

и

значений собственной

 

 

частоты

 

колебаний системы

 

можно

 

 

 

 

H

 

продолжить и дальше. Если же ограничиться третьим

приближением,

 

ω

 

 

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

уравнения (2.3.63):

α

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos2ω

 

H r

 

cosω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.3.74

 

16•

 

,

 

 

2- cos3ω ,

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

ω r ω

 

 

 

 

 

- .

2.3.75

 

 

,

12ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

179

 

 

 

 

 

 

 

 

 

2.3.5. Нелинейные резонансы

Резонансные явления в нелинейных системах оказываются более сложными по сравнению с аналогичными явлениями в линейных системах. Проиллюстрируем их проявления на примере одномерной

математической модели следующего вида:

2.3.76

 

 

 

 

 

 

HL 2λH ω H # cosγ αH βH2.

 

Проанализируем поведение системы при различных значениях

частоты

 

внешнего на нее воздействия и начнем с хорошо знакомого

случая,

когда

 

 

 

S

 

 

γ ω ε,

2.3.77

где

|ε| q 1

.

 

 

 

 

 

 

 

 

 

В этом случае в линейных системах с постоянными параметрами

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

 

 

 

 

ω ω ,

 

 

 

2.3.78

где

 

постоянный коэффициент,

определяемый параметрами

системы, а b – амплитуда вынужденных колебаний.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

 

ε γ ω

,

 

 

 

2.3.79

 

 

 

 

 

 

и если сохранить для символа

тот же смысл, что вкладывался в него

применительно

к линейным

системам, а именно:

 

 

,

то

 

 

 

 

ε

 

 

 

 

 

 

 

 

расстройка по частоте будет иметь вид

 

 

, и

выражение (2.1.76)

для

линейной

 

 

 

 

 

ε γ ω

 

 

для

системы трансформируется в

 

выражение

 

 

 

 

 

 

ε

 

 

 

 

 

 

 

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

180