Файл: Бондарь, Н. Г. Нелинейные стационарные колебания.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 22.10.2024

Просмотров: 111

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Порядок вычисления по формулам (1.237) и (1.238) такой же, как для (1.233) — (1.235).

На рис. 56 изображены

амплитудные кривые ультра- и субгар­

монических

колебаний

для

системы с параметрами а = 1 сек-2,

Р = 1 см-2

• сект2, п =

0,1

сек-1.

Сопоставляя кривые на рис. 7, 29, 56, устанавливаем, что, в от­ личие от системы с кубической нелинейностью, амплитуды суб- и ультрагармонических колебаний в системах с перескоком могут иметь тот же порядок, что и гар­ монические колебания с часто­ той возмущения со. Заметим, что

1.6 о,сен'

о

0.6

1.6

а.сеи'

Рис. 56. Амплитудно-частотные

характе-

Рис. 57.

Зависимость

амплитудно-

ристики ультрагармонических (штриховые

частотной

характеристики

системы

линии) и субгармонических (сплошные ли-

с перескоком от субгармонических

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

колебаний,

 

 

вязком трении для различных F,

см.сек~2.

 

 

 

 

кривая F = 1/4 на рис. 56 имеет разрыв, так как формулами (1.238),

как показано ниже, нельзя

пользоваться для В >• 0,5.

 

 

Возвратимся к рис. 30 и попытаемся объяснить расхождение меж­ ду аналитическим и машинным решениями, влиянием субгармони­ ческих колебаний. Для этого рассмотрим рис. 57, где сплошной линией изображена характеристика суммы амплитуд .колебаний с частотами со и (со — 20). Как видим, эта характеристика лучше со­ ответствует машинному решению, представленному точками, чем амплитудно-частотная характеристика колебаний с частотой со, ко­ торая изображена штриховой линией.

Высшие тона. Рассмотрим

случай

<

1. Тогда

\2В cos 20*| <

< 1 и можно воспользоваться биноминальным

разложением

1 + 2BCOS2Q/ = 1 - 2 8 C 0 S Ш

+ < 2 5 C 0

S Ш ) * ~

{

2 В C 0 S Ш

) 3 + " '

(1.239) Отсюда видно, что полученные выше результаты соответствуют ис­ пользованию двух членов ряда (1.239). Если воспользоваться тремя членами разложения, то уравнение (1.228) примет вид

2 " (е ) -J. z (е ) = .|_Cos со* [1 — cos 20* + (2В cos 20/)2] =

73


=

{(1 +

 

2Ba ) cos at —В [cos (со — 26) t + cos (со +

26) t] +

 

+

В2

[cos (со — 49) t +

cos (со - f 48) t]} .

 

(1.240)

Как видно из этого выражения,

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

| о» ± 201,

возникают ультра-

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

колебания с

частотами | со ±

40 [. Легко видеть, что для этих колебаний можно

воспользоваться

 

полученными

выше

результатами,

если

заменить

В на В 2 и (со ±

20) на (со ± 40). В частности, для Ау

и Л с без нали­

чия трения имеем для систем с кубической

характеристикой

 

 

 

 

А

_

 

ee*F

.

 

 

 

 

 

 

У , °

у ^ б 2

(со ± 4G2 )2 ] '

 

 

для систем

с перескоком

 

 

 

 

 

 

А

1 /

Р

л

g _

 

 

WB*

 

2 .

 

у.с г

 

 

~ у . с

 

 

2 _ ( ш ± 4 е

 

 

 

 

 

2

 

е

б ) 2

 

 

 

Г

2

У-с

у

 

4

е£_(в>± М )2

^

 

Используя формулы (1.236) и (1.238), при наличии вязкого тре­ ния совершенно аналогично получаем следующие выражения амп­

литудно-частотных характеристик колебаний с частотами | со ±

401:

для систем с кубической характеристикой

 

 

 

 

 

д

_

 

 

OFB2

 

 

.

 

 

 

У'°~

Уа {[л2 — (со ±

46)2 +

б 2 ) 2 +

4п2 (со ±

4в)2 } '

 

для

систем

с

перескоком

 

 

 

 

 

 

 

 

 

 

 

 

 

еб ра2

 

 

 

 

 

 

 

 

V 2

— (со ±

49б )2

-f- е^]2

+

2 (со ± 4вб )2

 

 

 

 

 

а

 

т > 2 ;

 

 

 

 

 

_

£

Л 2

 

1

 

 

в ^ Д »

 

 

 

 

2

' у > с

УЩ,

4

_ ( ш ±

29м )2

+ б2 ]2

+

2 (со '+ 29м )2

'

 

 

 

 

 

1

< т < 2 .

 

 

 

 

Эти формулы по сравнению с (1.236) и (1.238) имеют значительно

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

с

частотами

| с о ± 4 0 |

значительно меньше

амплитуд колебаний

с

частотами

| со ± 20|.

Отметим, что,беря в

ряду (1.239) четвертый и

последующий члены, можно обнаружить колебания с частотами | со ±

± 601; | со ± 801; \ со ± 100 |; .... Амплитуды этих колебаний опре­ деляются аналогично изложенному. Очевидно, что по мере возра­ стания порядка колебаний значения амплитуд будут уменьшаться.

Как видно из выражения (1.240), учет суб- и ультрагармониче­ ских колебаний с частотами | со ± 401 приводит к корректировке колебаний с частотой со на коэффициент (1 + 2В2 ).

74


Поэтому в формулах для амплитудно-частотных характеристик, приведенных выше, следует вместо F подставлять F (1 + 2). Вычисления показывают, что эта корректировка для уравнения Дуффинга практически никакого значения не имеет. Однако для систем с перескоком корректировка оказывает заметное влияние на большие колебания и несущественна для малых колебаний. В ка­

честве примера на рис. 58 приведены амплитудно-частотные

харак­

теристики систем с перескоком при

а — 1 сек-2; р = 3 см.—2

сек-2;

п = 0,05 сек-1, заимствованные из

работы [32]. Устойчивые

ветви

см

 

 

 

 

 

 

 

/;

 

 

 

 

 

 

 

 

-/Л

V

 

 

 

 

 

 

 

/

 

 

 

 

F=Q,

 

•ЛЛг

/ /

 

 

 

 

 

 

у/

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

0J

 

 

 

-

 

 

 

 

0.2

 

 

 

 

 

 

 

 

 

 

 

 

0,5

 

 

1.5

1,0

2,5

о, сек'

Рис. 58.

Скорректированные

амплитудно-частотные

характеристики

системы

с

перескоком:

 

 

 

 

 

 

 

О — F —

.—2.

Ф

- F .

0.2

см

• сек~

 

 

0,5 см * сек'

 

 

изображены

сплошными

 

линиями,

неустойчивые

— штриховыми;

точками и кружками представлены результаты решения задачи на ABM МН-7. Для больших колебаний наблюдается заметное несовпа­ дение машинных и аналитических данных, полученных по формулам

(1.103) и (1.105). Это несовпадение

устраняется, если провести кор­

ректировку амплитудных кривых,

как показано штрих-пунктир­

ной линией. Несовпадение с аналитическим решением

точки F =

= 0,5 см - сект2 и со = 1 сект1

объясняется влиянием

субгармони­

ки с частотой J со — 201 подобно

тому, как это показано на рис. 57.

§ 6. Устойчивость симметричных колебаний

Как известно [2, 3, 21, 24, 35, 36], стационарные коле­ бания нелинейных систем не всегда устойчивы. Понятие устойчи­ вости колебаний поясним на элементарном примере маятника (см. рис. 1). Пока амплитуда колебаний маятника меньше 180°.

75


колебания будут устойчивыми. Однако, как только амплитуда'; до­ стигнет 180°, колебательное движение перейдет во вращение вокруг точки подвеса. Как видно из рис. 2, а, в этом случае амплитуда-до­ стигает значения ненулевого корня характеристики, превышение которого меняет знак восстанавливающей силы, т. е. превращает ее в толкающую силу, приводящую к апериодическому движению, именуемому иногда в литературе [22,48] вращательным режимом. Этот термин связан с устойчивостью колебаний маятника. Таким образом, используемое здесь понятие устойчивости колебаний сле­ дует понимать как сохранение колебаний (орбитальная устойчи­ вость), а понятие критического состояния — как границу между ко­ лебательным и вращательным режимами.

С позиций качественной теории колебательному режиму (устой­ чивость) соответствуют фазовые траектории, заключенные внутри петли сепаратрисы, вращательному (апериодическому) режиму (не­ устойчивость) соответствуют на фазовой плоскости траектории, лежащие вне сепаратрисы [22, 48]. Переход колебательного движе­ ния в апериодическое иногда называют неустойчивостью «в малом» [41 ]. В нелинейных системах может иметь место также неустой­ чивость «в большом». Это касается систем с перескоком, в которых, как показано выше, большие колебания относительно неустойчи­ вости положения равновесия могут смениться малыми колебаниями относительно одного из двух положений устойчивого равновесия. Условия перехода одних колебаний в другие в системах с переско­ ком легко определяются из анализа амплитудно-частотной характе­ ристики. Поэтому вопрос устойчивости в большом можно считать решенным. Определить'условия перехода колебательного движения в апериодическое для устойчивости в малом из анализа амплитудночастотных характеристик не представляется возможным. Необходи­ мы специальные исследования, которые и составляют содержание настоящего параграфа.

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

Найдем корни характеристики уравнения (1.192):

(1.241)

Отсюда видно, что ненулевые корни будут действительными только для мягких характеристик (0 > В = — |В|):

(1.242) Если амплитуда колебаний достигнет значения (1.242), т. е. станет

(1.243)

76


то, как видно

из формул

(Г.208) и (1.209),

Имеют место

равенства

k = 1; K{k) = оо; 9 = -jjL-

У 2-

= 0.

-(1.244)

Следовательно,

функция

sin 20/ в приближенном

равенстве (1.210)

вблизи от состояния неустойчивости будет медленно изменяющейся.

Это обстоятельство позволяет заменить ее первым членом

разложе­

ния в степенной ряд. Тогда

выражение (1.210) принимает

вид е =

= Ф (t) = (1 + 2В) Qt,

откуда находим

 

 

Ф (0 =

(1 +

25)0;

t=

( 1

+ e 2 f l ) 6 .

(1.245)

Подставляя формулы (1.245) в уравнение

(1.12), получаем

 

г » + г ( 6

) =

е ( 1 ^ 2

Д )

cos

е ( 1 + 2 Д ) g -

< L 2 4 6 >

Частное решение уравнения (1.246), определяющее стационарные колебания, может быть получено аналогично (1.18), т. е.

9(1

+2B)F

 

в3 (1 +

2В)2 — о*

6.(1 + 2В)

Переходя здесь к старым переменным в соответствии с формулами (1.3) и (1.245), находим приближенное решение для стационарных колебаний вблизи от состояния неустойчивости при отсутствии тре­ ния:

.

9(1 +2B)F

.

fM=

e * ( i + W - o *

c o s a L

Полагая x = A и cos со/ = ± 1 и принимая во внимание формулу

(1.32), получаем выражение для амплитудно-частотной

кривой мяг­

кой системы с кубической характеристикой:

 

Л]Д--±-|В|Л* = ± yf/V+fff^ .

(1-247)

Из приведенных соображений следует, что полученные ранее результаты можно использовать также для колебаний вблизи от неустойчивости, если заменить в полученных ранее формулах 0 на 9 (1 + 2В). В частности, используя формулы (1.32) и (L99), получаем выражение для амплитудно-частотной кривой мягкой си­ стемы с кубической характеристикой при наличии вязкого трения:

А / « _ - ' | Р | * = ±

В ( 1 + 2 5 ) / ?

. (1.248)

77