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

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

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

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

Добавлен: 22.10.2024

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

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

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

Подставляя это выражение в уравнение (1.92), находим

( - §

( C i c o s - | - e

+

C2 sin-^-e)+ 2 - ^ - ( c a c o s - | - е —

 

/->

со

\

F

со

 

Сг

sin -g- ej =

-g- cos -g- e.

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

г ( я2

со8 . Л , о гко r _ F .

разрешая которую, определяем

 

 

 

г

а — coa + 6a )8f

r

2nM6F

 

1

~ (ha — со2 + е а ) а + 4ла ш2

'

2 ~ (п2 —соа + 9а )а + 4 я а « в

Представим частное решение (1.93) в одночленном виде. Для

этого

введем замены

 

 

 

 

 

С%

— а sin р, Сх = а cos р.

(1.94)

Теперь решение (1.93)

принимает

вид

 

 

 

п

 

 

 

 

z — ae

cos ( - l h - р ) .

(1.95)

Разделив почленно выражения (1.94), определим сдвиг фаз между колебаниями и возмущением:

Возводя в квадрат и суммируя равенства (1.94), имеем С\ + С\ = a2 (sin2 р + cos2 р) = а2 .

Отсюда находим

а ~

C z ~ У (п2 — со2 + е а ) а + 4л2 ш2 '

( 1 - 9 7 ^

Подставляя это выражение в решение (1.95) и переходя к старым переменным в соответствии с формулами (1.15) и (1.86), находим ре­ шение для стационарных колебаний в виде

^

9 f c o s M - p )

=

v

/ ( п а — со2 + 92 )2 + 4ла ш2

Полагая здесь д; = A, cos (cof — р) = 1„ получаем уравнение амп­ литудно-частотной характеристики

/ (А) = ±

B f

.

(1.99)

' v '

V V — со2 + 92 )2 + 4naco2

v

32


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

В отличие от систем без трения, амплитудные кривые, построен­ ные по формуле (1.99), не имеют разрывов вдоль скелетных кривых (рис. 21). Как показывают исследования [7, 35, 57, 62, 72], место перескока зависит от направления изменения частоты возмущения.

Допустим, что в случае жесткой системы мы начали увеличивать частоту возмущения от некоторого ее значения ©А, которому соот­ ветствует точка / на кривой. В точке 2 амплитуда получает наиболь­ шее значение и начинает уменьшаться до точки 3 (с частотой со2),

А

Рис. 21. К вопросу о перескоке с резонансных на нерезонансные колеба­ ния:

а — ж е с т к а я система; б — мягкая система.

где происходит перескок (срыв) амплитуды до ее значения в точке 4, после чего опять следует ее плавное уменьшение. Наоборот, при уменьшении частоты возмущения от значения со2 амплитуда дости­ гает точки 5, где совершает скачок до ее значения в точке 6, а далее плавно уменьшается. Таким образом, в точках 3 и 5 нарушается не­ прерывное изменение амплитуды. Это явление свойственно только нелинейным системам. Аналогично обстоит дело и в случае мягкой упругой системы, где разрывы имеют место в точках 2 и 5.

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

Как видно из рис. 21, амплитудно-частотные характеристики всегда имеют зону частот со, при которых возможны колебания с боль­ шими амплитудами (участки 6—3 на рис. 21, а и 5—3 на рис. 21, б), именуемые резонансными, а также колебания с малыми амплитудами (участки 5—4 на рис. 21, а и 6—2 на рис. 21, б), именуемые нерезо­ нансными колебаниями. Как показывают исследования [17, 62], реализация тех или иных колебаний зависит от начальных условий. Так, при нулевых начальных условиях всегда реализуются нере­ зонансные колебания. Резонансные колебания.реализуются при не­ нулевых начальных условиях.

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

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

для уравнения Дуффинга. Заштрихованные зоны

представляют

3 4-5

33


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

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

оно возможно только при рассмотрении нестационарных

колебаний.

 

 

в процессе

установления

и

 

 

поэтому

здесь

не

рассматри­

 

 

вается.

 

 

 

 

 

 

 

 

 

 

 

Для

определения

 

макси­

 

 

мумов

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

 

 

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

 

необходимо

 

 

найти

геометрическое

место

 

 

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

 

 

кривые имеют горизонтальные

 

 

касательные. Для этого следу­

 

 

ет продифференцировать урав­

 

 

нение

(1.99)

по со и положить

 

 

dA/da

=

0. Имеем

 

 

 

 

 

 

-±- f(A)

[(«2 — со2

+

в2 )2

+

 

 

 

 

 

i _

 

 

 

 

 

 

 

Рис. 22. К вопросу о реализации

резо­

+ 4/i2 co2 j

2 [ 2 ( / г 2 - с о 2

+ 6 2

)

X

нансных и нерезонансных колебаний в

зависимости от начальных условий:

 

X (— 2ш) +

8/г2со] =

0.

 

а — ж е с т к а я

система

( а = 1 « к - 2 ,

р = 1

 

см~-X

X с е к - 2 ,

л =

0,125 « к -

1 ,

F=l

см

• сек~2);

6—

мягкая

система

( а =

1

сек

2 ,

0

=

= — 1,8

см 2

• с е / с - 2

,

л

=

0.03

сек-1,

 

F =

Это равенство удовлетворяет­ ся при условии

 

 

со2 =

0 2 — п 2 .

(1.100)

Итак, получено

уравнение кривой, которая оказывается

близкой

к скелетной кривой со = 0, так как /г ^

1.

 

 

Если подставить условие (1.100) в уравнение (1.99), то получим

простую формулу для определения максимальных значений

ампли­

тудно-частотных

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

 

 

 

 

/ ( Л т а х ) = ^ Г .

 

(1-101)

Подставляя условие (1.100) в (1.96),

находим tgp

= со,

р = у .

Итак, стационарные колебания с максимальной амплитудой (1.101) отстают от возмущения на 90°.

Для определения точек срыва колебаний необходимо найти гео­ метрическое место точек, в которых амплитудно-частотные харак­ теристики имеют вертикальные касательные. Для этого следует продифференцировать уравнение (1.99) по Л и положить dcaIdA = 0.

J Эти результаты получены на ЭЦВМ «Урал-3» М. И. Казакевичем, Э. Н.Квашей, С. Ф. Редько.

34


Тогда получаем

уравнение

кривой

 

 

 

 

 

 

/' (А) ] / ( л 2 — со2 +

 

0 2 ) 2

+

4 лV

=

±-^-F 1

262 (га2 ш2

+ б2 )

 

 

(га2 — со2 62 )2 +

2 ш2

 

 

 

 

 

 

 

 

 

dA

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.102)

 

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

 

построение

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

характе­

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

конкретными примерами.

Сначала рассмотрим стационарные колебания маятника

(см.

рис. 1). Добавляя

к

уравнению

(1.25)

 

член, учитывающий

вяз­

кое

трение,

имеем

ij; + 2/ггр +

^

град

 

 

 

 

+ -у- sin т|э =

М cos at.

Подстав­

 

 

 

 

 

" 1

ляя

выражения

(1.27)

и

(1.28)

 

 

 

 

 

120]

 

 

 

 

 

 

град

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

80

 

 

 

 

 

 

 

 

ВО

 

 

4

 

 

 

ft

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

\ д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

40

 

 

 

 

 

 

 

 

40

 

 

 

 

 

 

О

 

 

 

 

 

 

 

 

0

 

 

 

 

с

-;

0,4

0,8

 

1,1

 

1,6

 

о,сей'

 

 

 

 

 

 

 

 

 

Рис. 24. Зависимости максимальных ам­

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

харак­

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

плитуд колебаний маятника от коэффи­

для различных

М, сек.—2

 

 

 

 

циента затухания для

различных М,

 

 

 

 

 

 

 

 

 

сек.—2

 

 

 

 

 

 

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

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

ки

для маятника:

 

 

 

 

 

 

 

 

 

 

 

 

 

2 sin А

{ [ " 2

-

со2 + -

f ( l -

4 " sin2

 

+ 4 n 2 0 2 } ~

=

 

=± м ( 1 - 4 - 5 т 2 4 ) .

Амплитудные

кривые, построенные по этой

формуле 1

при gll

=

= 1 сект-2, п =

0,15 сект1,

приведены на рис. 23.

 

 

Подставляя

выражение

 

(1.27) в формулу

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

за­

висимость максимальной

амплитуды

от коэффициента

затухания:

 

 

I

in

A a x

м

 

 

 

 

 

w l n

2

~

2га •

 

 

 

 

V-

 

 

 

 

 

 

 

На рис. 24 изображены кривые, построенные по этой формуле при qll = 1 сект2.

Далее, рассмотрим колебания осциллятора с кубической харак­ теристикой, описываемые уравнением Дуффинга [17, 69]:

х + 2пх + ах + рх3 = F cos at, a > 0, f> ^ 0.

1 Для удобства вычислений это уравнение следует разрешить относительно частоты возмущения.

3*

35


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

 

 

 

 

Л ^ а + ^ б Л ^ п2

co2 + a ( l +

 

 

 

 

+

0,756-^-Л2

+

4/г2 со2 }~Г

 

 

 

 

1

'

а

 

 

 

 

 

 

 

= ± / ? ] / A a ( l

+ 0,756 i

Л 2 ) .

 

0,6

1.0

1,4 Цсе/С*

Построенные по этой формуле1

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

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

осциллятора с

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

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

кубической характеристикой при

ки изображены на рис. 25—27. Точка­

вязком трении

для a =

1 сек

ми на

рисунках отмечены

результа­

 

 

 

.—2

 

 

 

 

 

 

р =

1 см~2 -сек-2, F=0,25 смХ

ты *, полученные на ЭЦВМ «Урал-3».

X

сек~2.

 

 

В случае

решения

при нулевых на­

 

 

 

 

чальных условиях

реализованы нере­

зонансные колебания (см. рис. 25, 26), при различных начальных условиях реализованы как нерезонансные, так и резонансные коле­ бания (см. рис. 27). Совпадение аналитических и машинных решений

можно признать

хорошим, за исключением

участка со = 0,4

-f-

- г - 0,8 сек~1

для

F = 2 см • сект2 на рис. 27.

О причине этого

не­

соответствия

будет сказано ниже (см. § 5 данной главы).

 

см

/ /

1,2

1

 

 

1

 

 

 

1,2

 

2.4

а.сек'

°'1

 

 

 

 

 

Рис. 26.

Амплитудно-частотная кривая

Рис. 27. Амплитудно-частотные кривые

осциллятора

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

осциллятора

с

кубической

характери­

кой

при

вязком

трении

для

a =

стикой

при

вязком трении

для a

=

= 1

сек~2,

Р = 1

см~2

сек 2,

п <

1 сек 2 ,

р =

1 сек.—2

..—2

 

* 0,13 сек~\

F =

1 см •

секГ2.

 

= 0,16

сек 1

и различных F,

см-сек

-.

1

Рис. 25

—решения получены Л. В. Велик и А. С. Третьяковой.

 

 

4

Рис. 26 — результаты получены

В. И.

Глухих

и Т. И. Тимошенко,

рис. 27 — Н. Г. Новиковой и Л. И. Бобровой.

36