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

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

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

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

Добавлен: 22.10.2024

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

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

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

Используя формулы (1.32) и (1.139), совершенно аналогично полу­ чаем амплитудно-частотную характеристику при турбулентном сопротивлении

А

уа-

21

| Р И 2

=

_L ,

e ( 1 +

2

g ) F

(1.249)

В (1.247) — (1.249) 9 и В определяются по формулам (1.209) и (1.212). Следует заметить, что полученные выражения легко обобщают­ ся на случай любой мягкой системы, если заменить левые части фор­

мул (1.247) — (1.249) соответствующим значением амплитудной функции / (А) и пользоваться точным значением 9.

Критерии неустойчивости систем без трения. В критическом со­ стоянии, когда амплитуда стационарных колебаний . определяется выражением (1.243), по формулам (1.212) и (1.244) находим

(1.250)

Подставляя выражения (1.243) и (1.250) в равенство (1.247), прихо­ дим к уравнению критического состояния

Попытаемся оценить точность этой приближенной формулы. Для этого рассмотрим стационарные колебания, описываемые урав­

нением

"x + ax — \fi\x* = Fsn(at, k) = Fsnu,

(1.252)

где snu — эллиптический синус. Будем искать частное

решение

уравнения

(1.252) в виде

 

 

x = Asn(at,k).

(1.253)

Подставляя решение (1.253) в уравнение (1.252) и группируя члены, имеем А [—со2 (k2 2k2 sn2 и + 1) + а — ] f\\ А2 sn2 и — F/A} sn и = 0. Приравнивая нулю квадратную скобку и группируя члены, по­ лучаем (2k2a2 — ВЛ2) sn2 и + (a— a2k2 со2 F/A) = 0. Это уравнение удовлетворяется при условиях 2&2со2 — ВЛ2 = 0; а —

— со2 (1 + k2) — F/A = 0, откуда находим 1

А = « - ( i + » ) a , * = ( ] - 2 5 4 > Подставляя формулы (1.254) в равенство (1.253), приходим к

точному решению:

Р

 

х =

sn (at, k).

j

a - ш 2 - — | | 3 | Л 2

1 Этот результат приведен

в работе

[1].

78


Полагая здесь х = А и sn (cor', k) = 1, получаем точное выражение

для амплитудно-частотной

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

Л =

F

 

a - c o * - i - | P | Л2

Подставляя сюда формулу (1.243), приходим к тому же уравнению (1.251) критического состояния.

Итак, уравнение (1.251), будучи приближенным для тригономет­ рического возбуждения, оказалось точным для эллиптического воз­ буждения. Поскольку эллиптические функции разлагаются в быстросходящиеся тригонометрические ряды [23, 28, 55], точность формулы (1.251)

оказывается удовлетворительной. Этот вывод подтверждается решениями на ABM МН-7.

 

j/2 yS-- /

Z a,air<

 

и.сек

Рис.

59.

График

критических

Рис. 60. К построению кривой / /

графика

состояний

симметричных

коле­

критических состояний.

 

баний

без трения.

 

 

 

 

На рис. 59 кривая

/ построена по формуле (1.251) для системы

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

| р | = 0,2 см~2 сек—2. Заметим, что

для левой ветви кривой / в формуле (1.251) перед радикалом

берет­

ся знак

«плюс», для правой

ветви — знак «минус». Кривая

/ / на

этом же рисунке построена на основании следующих соображений. Как известно [7, 35, 36, 47, 57, 62], в нелинейных системах возможно

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

колебаний.

На рис. 60 изображены два типа амплитудно-частотных

характери­

стик мягких систем *. В первом типе (кривые III) перескок с нере­

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

колебани­

ям, во втором типе (кривые II) — к апериодическому движению (не­

устойчивость колебаний).

Как видно из рис. 60, перескок колеба­

ний происходит в точке,

где касательная к

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

1 Для примера принята система с параметрами a =

1 сек~2, | р ) = 0,2 см~2Х

X сек~2 и F = 0,5; 0,05 см •

сек~2.

 

79



характеристике становится вертикальной. Следовательно, Для1 'оп­ ределения условий перескока колебаний необходимо продифферен­

цировать уравнение амплитудной кривой по Л и положить

= 0.

Поскольку точки срыва колебаний находятся вдали от критической амплитуды колебаний, то можно воспользоваться уравнением ам­ плитудно-частотной характеристики в виде (1.23), т. е.

 

f(A)[Q2(A) — со2] =FB(A).

(1.255)

Дифференцируя

выражение (1.255) по Л и полагая

= 0, полу­

чаем 1

/'(82 — со2 )+ 2/0 • 9' = F9\

(1.256)

 

Выражения

(1.255) и (1.256) представляют собой уравнения кри­

вой / / критического состояния в параметрической форме. Вычисле­ ния по ним производятся в следующем порядке. Для мягкой систе­ мы с заданными параметрами задаются амплитудой Л и, вычисляют / (Л), 8 (Л), f (А) и 9' (Л). Затем решают систему уравнений (1.255) и (1.256) и находят со и F. Если эти величины вещественны и поло­ жительны, то на графике наносится точка с координатами со и F. Следует заметить, что эти вычисления весьма громоздки. Оказывает­ ся, что результаты могут быть аппроксимированы простым выраже­ нием |

Это выражение получено из следующих соображений. Принимаем линейное приближение для частоты свободных колебаний 92 = а. Тогда£0' = 0 и уравнение (1.256) принимает вид f (a — со'2) = 0. Поскольку a —1 со2 Ф 0, то f = 0. Отсюда f = с = const, т. е. ли­ нейное приближение для частоты влечет за собой осреднение ампли­ тудной функции. Теперь уравнение (1.255) принимает вид

c(a — a2) = VaF.

(1.258)

Ниже будет показано, что кривая

/ /

пересекается

с осью орди­

нат в точке

 

 

 

со = 0; F - F u

-

^ y ^ .

d-259)

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

—§-«Vw- . (L260)

Подставляя выражение (1.260) в (1.258), получаем уравнение (1.257). На рис. 59 кривая / / построена по формуле (1.258) для тех же параметров, что и кривая /. Точками и кружками показаны

1 Штрихами обозначены производные по А.

80


результаты решения 1 на ABM МН-7. Как видим,

совпадение

ана­

литического и машинного решений удовлетворительно 2 .

 

 

 

Заметим, что кривые / и / / н а рис. 59 пересекаются. Это проис­

ходит потому, что кривая

/, как видно из выражения

(1.251), при

F

= 0 и со =

0 пересекается

с горизонтальной осью в

точке

со =

=

"j/' - y - и с

вертикальной

осью

в точке Ft =

|/"-гуг- •

Кри­

вая / / пересекается с

горизонтальной

осью в точке со = ] / а .

Этот

результат вытекает из следующих

соображений.

Полагая

в

урав­

нениях

(1.255) и (1.256)

F = 0,

получаем

 

 

 

 

 

9 = со;

2/09'=0;

299' = 0; _dQ2/dA = 0;

02 = const.

 

Поскольку для малых F при со =

У а (см. рис. 59) будут иметь место

малые линейные колебания, то 03

= а.

 

 

 

 

 

 

Определим точку пересечения

кривой

/ / с вертикальной

осью.

Полагая в формулах (1.255) и (1.256) со =

0, находим

 

 

 

 

 

 

 

/9 =

/?//;

/'92

+

2/00' =

FnQ'.

 

 

(1.261)

Исключая отсюда Fu,

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/'0 +

/0' =

0.

 

 

 

(1.262)

 

При статическом воздействии

на

систему (со = 0,

х =

0,

х =

=

Л„)

уравнение (1.192)

принимает

вид

 

 

 

 

 

 

 

 

 

 

« A 0 - | B | A o W / 7 .

 

 

(1.263)

Дифференцируя это равенство по Л0 , имеем

 

 

 

 

 

 

 

 

 

а 3 |Р \А\ = F'u.

 

 

(1.264)

В то же время, дифференцируя

по А0 первое выражение

(1.261),

находим f'0

-f- /8' = F/j.

Сопоставляя это равенство с формулой

(1.262), определяем F'u =

0. Далее, из выражения (1.263) получаем

А0

=

" j / " 3 "p j

. Подставляя это значение в равенство (1.263), по­

лучаем формулу

(1.259).

 

 

 

 

 

 

 

 

 

 

 

На рис. 59 точки, лежащие выше сплошных линий,

характери­

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

которых реализуются

нерезонансные

колебания

(см. § 2 данной

главы). Штриховая часть кривой / / характеризует

параметры

воз­

буждения, при которых происходит перескок с нерезонансных

коле­

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

резонансные

(см. кривые

/ / / на рис. 60).

1 Решения на АВМ, приведенные в этом параграфе, получены В. С. Горбато­

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

резо-

нансов субгармонических колебаний, т. е. при со = 6/3, 6/5...

 

2 Следует иметь в виду, что погрешность вычислений на АВМ МН-7 может достигать 10—15%. Причем особенность решения задач устойчивости на АВМ МН-7 состоит в том, что результат для кривой / всегда занижен.

6 4-5

81