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

книги из ГПНТБ / Бондарь, Н. Г. Нелинейные стационарные колебания

.pdf
Скачиваний:
30
Добавлен:
20.10.2023
Размер:
7.67 Mб
Скачать

колебаний, описываемых уравнением

(1.54):

 

 

 

 

cos at

— В

cos (ш — 20) /

+

,

cos (со + 29) /

•со2

• (со — 28)2

} + A - ( l + - | - S c o s 2 0 / ) .

(1.319)

е а (со+29)2

Сопоставляя решения (1.231) и (1.319), устанавливаем, что в не­ симметричных системах при несимметричном возбуждении возника­ ют суб- и ультрагармонические колебания такие же, как и в сим­ метричных системах. Следовательно, для определения амплитуд этих колебаний можно полностью воспользоваться результатами § 5 на­ стоящей главы, только 0 и В следует определять по формулам данного параграфа. Наличие несимметрии приводит к появлению ко­ лебаний с частотой 20, амплитуды которых, как видно из формулы (1.319), определяются уравнением

Л е Y a * + ~Т Hi e =

.

(1.320)

Вычисления показывают, что имеет вместо неравенство сх^ ^> Это позволяет упростить (1.320) к виду

Л 2 0

д а - | | £ ^ .

 

30 у а .

1 ал?

-^-РЛ

(1.321)

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

У | / « * + 4 " = -|р + a cos (at — р) + ас cos [(со — 20) t — рс ] —

— а у cos [(со -+- 20) г? —Р у ] + а2о cos (20с" — р2е),

где амплитуды колебаний и сдвиги фаз определяются по формулам, приведенным в § 5 данной главы. Кроме того,

<z28 =

г

206,5

\

 

Р29 =

,

4я0

 

'

arctg

л * - 3 0 * •

 

W - 3 e » ) * + 1 6 i i * e »

 

V M

S

Сопоставляя полученные данные с результатами § 5 данной главы, видим, что в несимметричных системах при вязком тре­ нии возникают суб- и ультрагармонические колебания, амплитуды которых можно определять по формулам § 5. Кроме того, возникают колебания с частотой 20, амплитуды которых определяются из приближенного уравнения, полученного аналогично (1.321):

Л 2 9

да *

= -

2 6 6 * 5

.

(1.322)

 

Уа*

 

/ а Л ( л 2 - 3 0 2 ) а +

 

V

Полагая в формуле (1.322) п = 0, как и следовало ожидать, при­ ходим к выражению (1.321).

100

В случае несимметричного возбуждения симметричных систем в приведенных выше формулах следует использовать равенства (1.291).

Более подробно1 рассмотрим стационарные колебания при турбу­ лентном сопротивлении. Используя замены (1.3) и формулы (1.133)

и (1.134),

приведем уравнение (1.303) к виду, аналогичному (1.12):

 

 

г" (Б) + 2 (е) =

- Л -

ек «-*) (F0 + F cos

at).

 

 

 

 

ф(<)

 

 

Подставляя сюда формулу (1.227), приближенно получаем

 

 

2

(г) + це)

- е

в (I + 2В cos 2Ш)

 

~

 

(FQ +

Fcosat)(\

 

„к (1-х)

{F0 + Fcosat —

е

— 25 cos26/) = -=—g

— 2BF0 cos 26i! — BF [cos (ю — 26) t + cos (со + 26) *]}.

Используя приближенное равенство (1.15), запишем это урав­ нение в виде

г"(е) + 2(е) = -^-е (Т ~*) \pQ + Fcos~е

2BF0cos2е —

c o s ( 4 — 2 ) е - b c o s ( - r + 2 ) e }•

Частное решение этого уравнения, полученное аналогично из­ ложенному выше (см. § 3, 5, 7), имеет вид

2(e) = Л °

Ti{4

+ a c o s ( - f e _ p ) +

a c cos[( - f — 2 )

+

 

+

aycos

-g- +

2) е — р у

+

а2 е cos (2s — р2 е)},

(1.323)

где

 

 

 

 

 

 

 

 

 

а —

 

8F

 

 

 

 

 

 

 

 

 

 

 

Л 2 я 2

 

 

 

 

 

 

 

+

16

 

 

р =

arctg

(4Л2 л2 — л2 ) со2

+ 82 it2

'

 

/ [ ( c o - 2 e ) 2 ( - i ^ n 2 _ i )

+

e2

+

16- АС*П?(со—26)4

 

 

 

 

 

 

 

 

 

я а

 

 

 

 

ААспк

(со

28)2

 

 

 

 

рс --- arctg • (4Лс 2 л2 — я2 ) (со — 28)2

+

8 2 я 2 '

 

1 В § 5 влияние сопротивления не рассматривалось.

101

 

 

 

 

 

 

QFB

 

 

 

 

 

 

 

 

V

 

 

 

У . я 2

— 1

+

0 2 1

+16

У

 

(со +

20)1

 

 

я 2

/

 

 

J

 

 

я 2

 

 

 

 

 

 

 

 

 

у пл (ш-f 29)2

 

 

 

 

 

 

 

р у = arctg • (4Л2 я2 — я2 ) (со +

2G)2 +

0 2 я 3

'

 

 

 

а 2 е

=

 

 

 

 

Рге

 

 

 

 

16Л2дЛЛ

 

 

0 ( 1 6 Л 2 2 0 л 2

—Зя 2 ) *

a r

c

t g

16AW— Зл2

 

 

Возвращаясь в решении (1.323) к старым переменным в соответ­

ствии с заменами (1.3) и (1.15) и учитывая приближенные

равенства

(1.134) и (1.143), для стационарных

колебаний

получаем выражение

/ (х) =

- А . _|_а

cos (at — р) +

ас

cos [(со 2В) t — рс ]

+

 

+

ау cos [(со + 28) t — ру ] +

а2 9

cos (20t р2о).

 

 

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

суб- и

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

 

колебаний

определяются

 

по

формулам

/ (Лс ) = ас ;

/ (Лу ) = ау;

/ (Л2в) = а2в.

 

Подставляя сюда

приве­

денные выше выражения для а^, ау

и а2 в

и используя приближенное

равенство / (Л) =а Л У~а, для амплитудных

кривых

получаем урав­

нения

 

 

 

 

ЛС | У

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

±

 

 

 

 

 

9FS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/

 

со Т 20). 1

4 <

У

1

+ 0 2

 

+ 16

с, у

л 2 (со ¥

20)4

 

 

 

 

 

^ 2 0

= =F

^ " 5

л 3

 

 

 

.

 

 

 

(1.324)

 

 

 

 

 

QVa.(\&A2n>

Зл2 )

 

 

 

 

v

Знак «минус» в выражении

(со =р 29) соответствует

амплитудам Л с

субгармоник,

знак «плюс»—амплитудам Л у

ультрагармоник. Пола­

гая в формуле (1.324) п =

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

выражению

(1.321).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следует

отметить,

что

если

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

разложение

(1.240),

то аналогично изложенному выше (см. § 5) можно обнаружить коле­

бания с частотами | со ±

49|; | со ± 691; | со ± 801;

а также коле­

бания с частотами 2iB,

I =

1, 2, 3, 4,

амплитуды которых зна­

чительно меньше амплитуд

колебаний

с частотами

| со ± 20 \ и 26.

Устойчивость колебаний. Сначала

рассмотрим

простейший слу­

чай отсутствия трения, т. е. когда колебания описываются уравне­ нием (1.54), где параметры и \ определяются по формуле (1.282).

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

102

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

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

9 ( 1 + 2 5 ) = - ^ .

= 2ц.

Из третьей формулы (1.296) видно, что величина ц. выражается через корни уравнения энергии (1.295). Выразить ее аналитически через параметры системы невозможно. Однако, как показывают вы­ числения, хорошую сходимость дает весьма простое выражение

е(1 + 2В) = 2 | х » | / Г - ^ - ,

(1.326)

аналогичное формуле (1.250) для симметричных систем. Причем не­ симметрия системы учитывается параметром а), определяемым из уравнения (1.287) и первого равенства (1.285). Точность приближенно­ го выражения (1.326) тем выше, чем меньше несимметрия, т. е. чем меньше \у^\.

Подставляя формулу (1.326) в равенство (1.325), после простых преобразований получаем уравнение критических состояний

F =

(1.327)

 

m[n

Входящая в эту формулу критическая амплитуда Л к р равна ненулевым корням характеристики уравнения (1.281), которые свя­ заны с ненулевыми корнями характеристики уравнения (1.284) пер­ вой зависимостью (1.292), т. е.

Ар = У = У* — - gjjj-r •

Ненулевые корни характеристики уравнения (1.284), как извест­ но [13], имеют вид

У * = ± ТГРГ ( ^ v 2 + 4 a : | P | ± У*

Следовательно, для критических амплитуд получаем два выра­ жения:

А«р = "QTRT f"Г" ± l / V ; + 4a:|P|).

(1.328)

2IPI

 

Из этих-двух значений в формулу (1.327) подставляется то, которое дает минимальную по абсолютному значению амплитуду возбужде­ ния F.

Заметим, что принимая в формуле (1.327) Л к р согласно выраже­ нию (1.243) и полагая а* = а и б* = 0, получаем, как следовало

103

ожидать, уравнение (1.251) критических состояний симметричной системы. Полагая в формуле (1.327) попеременно F = 0 и со = О, получаем координаты точек пересечения кривой / критических со­ стояний с осями координат (рис. 71):

< и Fi =

 

V а

 

 

а .

min

Для определения

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

места точек срыва

колебаний

продифференцируем

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

по

Л и положим

 

= 0.

 

Имеем

 

 

 

 

(/'0 +

/ 0 ' ) ( 0 2 - с о 2 ) +

2 ( / 0 -

 

 

— б*) 60' = ± 200'F. (1.329)

о.сек

Рис. 71. Кривые критических состоя­ ний несимметричной системы без тре­ ния.

Это уравнение должно разре­ шаться совместно с уравнением амплитудно-частотной характе­ ристики (1.54), которое запишем так:

(/0 — SJ (92 — со2) = ± 0 2 F . (1.330)

Формулы (1.329) и (1.330) представляют собой параметри­ ческие уравнения кривой / / кри­ тических состояний (см. рис. 71). Координаты точек пересечения кривой // критических состоя­ ний с координатными осями, най­ денные так же, как в § 6 данной главы, следующие:

(1.331)

min

Использование выражений (1.329) и (1.330) затруднительно из-за громоздкости вычислений. Поэтому построим приближенную форму­ лу, приняв линейное приближение для частоты свободных колеба­ ний 02 = а;. Тогда 0' = 0 и выражение (1.329) принимает вид f У о ; (а; — со2) = 0.

Поскольку а ; ^ 0 и а ; — со2 Ф 0, то /' = 0, откуда / = с = = const, т. е. линейное приближение для частоты 0 (так же, как и для симметричных систем) влечет за собой осреднение амплитудной функции. Теперь уравнение (1.330) принимает вид

Val - 6J (о! - со2) = ± <F.

(1.332)

104

Используя граничные условия со = О, F = Fц, получаем

Здесь использована формула (1.331). Теперь уравнение критических состояний (1.332) принимает вид

F =

( ± - г а * / з Т р Т + б *

1 — со

(1.333)

 

а*

 

На рис. 71 кривые критических состояний построены по форму­ лам (1.327) и (1.133) для системы с параметрами а = 1 се/с- 2 ; у =

= 0,15 см~х • сект"1; В = —0,2 см~2 • сект2; 6„ = —0,25 см • сект2; F0 = 0,25 см • сект2.

Кружками (кривая /) и точками (кривая //) представлены ре­ зультаты решения 1 на ABM МН-7. Как видим, совпадение аналити­ ческих и машинных решений можно признать удовлетворительным, если учесть, что машинные решения для кривой / занижены на 10—

15%. Физический смысл кривых и порядок пользования

рис. 71

такой же, как на рис. 59 для симметричных систем (см. § 6

данной

главы).

 

Далее рассмотрим устойчивость колебаний с учетом вязкого трения.

Пользуясь формулой (1.110), записываем уравнение амплитудно-

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

в виде

 

F = ± - gr (/в -

6J V(n* - со2 + 82 )2 + 4га2

(1.334)

Используя выражения (1.32) и (1.334) и заменяя частоту 9 на

9 (1 +

+ 25), аналогично (1.248) получаем следующее уравнение для амп­ литудно-частотных характеристик колебаний, близких к неустой­ чивости:

± Q2 (1 + 2 S ) 2

=FB(I

+ 2 5 ) Л ] / " а # - ^ - | В | Л 2 - (1.335)

- б , У [п2

— со2 +

92 (1 +- 2В)2 ]2 -+ 4/г2со2.

Подставляя сюда приближенное выражение (1.326), приходим к уравнению критических состояний (кривая /)

F=

А,

(±-^rV2a*-\V\A2Kp-

 

 

 

. (1.336)

1 Результаты решения

на АВМ, приведенные в этом параграфе, получены

В. С. Горбатовым.

105

В формулу (1.336) подставляется то значение (1.328) критической амплитуды /4к р , которое дает минимальное абсолютное значение F.

Исследуя выражение (1.336) на экстремум, определяем координа­ ты минимума кривой / критических состояний:

со - / 4 - « : - / . «

 

F =nV2 ±AK?V2am-\^\Alp-y^-

(1.337)

 

min

Заметим, что приближенное равенство (1.326) справедливо в случае отсутствия трения. Поэтому формулой (1.336) можно поль­ зоваться для малых сопротивлений; при больших сопротивлениях естественнее воспользоваться приближенным равенством

2 + 6 2 ( 1 + 2 В ) 2 ^ - ^ - .

Подставляя это выражение в формулу

С1.335),

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

критических состояний для больших значений

/г и со <С Уа; г:

± - ^ 1 / 2 а , - | 6 | 4 р

 

со2 +4п2 со2

V V а*

 

 

 

 

(1.338)

Для приближенного определения геометрического места точек срыва колебаний в соответствии с изложенным выше примем в

уравнении (1.334) линейное приближение для частоты 02 +

пг = ai

и осредним

амплитудную

функцию /

с = const.

Имеем

F

=

• (с У а,

— бJ У (а: со2)2 +

4n2co2.

(1.339)

 

 

а;

 

 

 

 

 

 

 

Постоянную

с найдем

из

граничных

условий

со =

О, F — F/t.

Используя формулу (1.231), получаем такое же значение с,

как для

случая отсутствия трения. Подставляя его в формулу

(1.339), прихо­

дим к следующему уравнению критических состояний

(кривая / / ) :

F =

(1.340)

1 Ограничение со < J^a* вытекает из сопоставления аналитических и машин­ ных решений.

106

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

ний п. При

больших значениях п лучшую

сходимость с машинными

решениями

дает

уравнение.

 

 

 

 

F = - 4 ( ± - ! - « , / V r p T +

 

 

+

l / ( a : со2)2 + 12шо

- со)2

(1.341)

На рис. 72 приведены построенные по формулам (1.336) — (1.341) графики критических состояний для системы с теми же параметрами,

 

 

 

 

 

 

 

 

а.сен

 

 

1

/

о,сек'

Рис. 72.

Кривые

критических

состоя­

Рис. 73.

Кривые

критических

состоя­

ний несимметричной системы при вяз­

ний несимметричной системы при тур­

ком трении:

 

 

 

 

булентном

сопротивлении:

 

 

 

 

 

 

 

 

 

 

 

/

л

=

0,05 сек-1;

2

л =

0,20

сек—

t — п —

0,05

cm~1;2

— л = 0,20

см~1.

3

п

=

0,30

4

л =

0,50 сек~1

 

 

 

 

 

что и на рис. 71. Здесь же условными обозначениями представлены результаты решения на ABM МН-7. Как видно, совпадение резуль­ татов можно признать удовлетворительным. Порядок пользова­ ния графиками такой же, как и для симметричных систем (см. § 6 динной главы).

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

на

устой­

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

Корректируя путем

замены 0 на

0 (1 + 2В),

а на а.^ и F0 на

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

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

(1.308),

для несимметричных колебаний мягких систем, близких

к неустой­

чивости, имеем

107

0(1 + 2B)F

 

 

6 (1 + 2В) ^

'

Л2 Я2 д

] / [ с о 2 ( ~ - л - 1 ) + 8* (1 + 2Bf

+ 16

— ш4

Подставляя сюда приближенное выражение (1.326), получаем урав­ нение критических состояний (кривая /)

+ 1 6

(1.342)

В эту формулу подставляется то значение (1.328) критической ам­ плитуды Л к р , которое дает наименьшее абсолютное значение F.

Параметрические уравнения геометрического места точек срыва колебаний (кривая // ) весьма громоздки и неудобны для использо­ вания, так как они содержат амплитуду колебаний. Вычисления по­ казывают, что хорошее совпадение результатов дает простое вы­ ражение 1

 

F

=

 

 

+ 8^ У

( о ; - со2)2 + 6Д, (У< - ш)3>

(1.343)

где Ал

абсолютная величина минимума характеристики

системы,

3|PI Imin"

На рис. 73 изображены построенные по формулам (1.342) и (1.343) кривые критических состояний для системы с теми же параметрами, что и на рис. 71, при различных значениях п. Дискретными ус­ ловными обозначениями представлены результаты решения на ABM МН-7. Как видно, совпадение результатов удовлетворитель­ но. Порядок пользования графиками такой же, как и для сим­ метричных систем.

1 Эмпирические формулы, приведенные в этом параграфе, подобраны В. С. Гор­ батовым.

Г Л А В А I I

БИГАРМОНИЧЕСКОЕ

ВОЗБУЖДЕНИЕ

§ 1. Симметричное

возбуждение

 

систем без трения

 

 

 

Будем рассматривать

(§ 1,2)

возмущение в виде двух

гармоник: F (t) = Fx cos coji + ^2 c o s

F (0 =

^1 s ' n W J +

-f- F2 sin co2*. Как известно, эти

функции

будут

периодическими

только в случае кратности частот, т. е. со2 = р.©!, ц = 0, 1,2,3, ...

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

F (t) = Fx cos со* + F2 cos LICO*;

F (t) = ^ s i n at +

F2 sin p,co*. (II. 1)

Частота со называется основной (низшей) частотой возмущающей

силы F (t). При р, = 0 первая формула (11-1) дает

несимметричное

гармоническое возбуждение (см. § 7 гл. I), вторая формула

( I I . 1) —

гармоническое возбуждение

(см. § 1—6 гл. I). При

L I =

 

1 формулы

( I I . 1) описывают

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

возбуждение

с частотой

со и ам­

плитудой Fx -f- F2.

Несмотря на аналогичность

выражений

(II.1),

описание

стационарных колебаний

 

для них

несколько

различно.

Поэтому

будем

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

и

именовать

соответ­

ственно

четным

(косинусоидальное)

и нечетным (синусоидальное)

возбуждениями.

 

 

 

 

 

 

 

 

 

 

 

 

Четное возбуждение. Рассмотрим

[34] стационарные

вынужден­

ные колебания

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

 

 

x(t)

+ R(x)

= F x

cos со* + F2 cos LICO*.

 

 

 

(II.2)

Используя замены

(1.3), преобразуем нелинейное уравнение

(П.2)

к линейному, аналогичному (1.12), т. е.

 

 

 

 

 

 

 

г" (е) +

г (е) = —^— (F±

cos со* + F2 cos ixco*).

 

 

 

 

 

 

 

 

ф(9

 

 

 

 

 

 

 

 

 

Линеаризуя

фазовую функцию в соответствии с выражением (1.13)

и используя

зависимости (1.14) и (1.15), получаем

 

 

 

 

 

г" (е) + г(е) = -L

(F± cos

 

е + F2 cos ц -f- е).

 

 

(П.З)

109

Соседние файлы в папке книги из ГПНТБ