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

книги / Прочность, устойчивость, колебания. Т. 3

.pdf
Скачиваний:
6
Добавлен:
13.11.2023
Размер:
24.79 Mб
Скачать

Колебания нелинейных систем с одной степенью свободы

261

в которой уг (/), уг (/), — должным образом выбранные функции

вре­

мени; а1г а2,

— параметры, определяемые

из системы уравнений

 

 

 

I

[ту + Р (У) — Р 0 51П 0)/] у1(/) = 0

 

о

(1 = 1, 2, .

.,я).

(96)

 

При сохранении одного члена суммы (95)

получается соотноше­

ние (93).

 

 

 

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

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

Р = Рх51П < 0 -\~Р 2 5Й1 00а/,

на систему с нелинейной характеристикой

Р (у) = су + V!/3

приближенное решение имеет вид

У = а1 з т по>]/ + а2 з т соа/.

Амплитуды а1 и а2 определяют из нелинейной системы уравнений:

=____________ Рг_____________ .

(р“ — ©?) + ~4"

(а? + 2аг)

________________ Р з________________

(р2- ш 1 ) + - | - | - ( а2 + 2а?)

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

ту + ку + Р {у) — Р0 з!п о /

(97)

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

у = а з т (ш/ — у).

(98)

Амплитуда колебаний может быть определена из уравнения

 

Р (а) = VРр — (Ласо)2 + таю1

(99)

или из уравнения

 

а = -

( 100)

йаша

 

/ [ 1- -й)Т+т2р* (а)

262

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

в котором р- (а) определяется формулой (78). На рис. 20, а схематически показана амплитудно-частотная характеристика, а на рис. 20, б — из­ менение амплитуды а при изменении частоты ю возмущения от нуля до значения со* и затем от значения со* до нуля. В системах с демпфирова­ нием срыв амплитуды неизбежен даже при монотонном увеличении

частоты.

Для построения высших приближений используют метод БубноваГалеркина. Закон движения принимают в виде суммы (95) и вместо уравнений (96) получают систему уравнений

(I)

 

|

[ту + Ьу 4- Р (у) — Р„ 5’щ со/] У1

 

о

 

 

(/ = 1, 2, . . . . л)

(101)

с неизвестными

параметрами

 

Нелинейные диссипативные системы

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

Свободные колебания. С и с т е м а с к в а д р а т и ч н ы м з а ­

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

Дифференциаль­

ное уравнение движения может быть приведено к виду

У±-% -У*+РгУ - 0 ,

(102)

где Р — коэффициент, зависящий от вязких свойств системы. Если о0 = = 0 и с0 — начальное смещение системы, то следующее наибольшее по величине отклонение ах (достигаемое через полупериод колебаний)

определяют из

трансцендентного

уравнения

 

 

;п (1 + рД1) -

ра1 =

1п (1 + ра0) - &а0.

(103)

Следующее

наибольшее

по величине отклонение а2 (достигаемое

по истечении еще одного полупериода) определяют из трансцендентного уравнения

1п (1 — ра2) раа = 1п (1 — РаО + Ра,

(104)

и т. д.

Для графического последовательного решения уравнений типа (102) и х.103) удобно воспользоваться стандартной кривой (рис. 21, а)

т] =1п(1 +Е)-6-

(Ю5)

Вычислив по данным задачи

= ра0, следует по стандартной

кривой найти соответствующее

значение Г|0 (точка

А 0 справа от оси

ординат). Проведя через точку

А 0 горизонтальную

прямую, находим

точку А 1ъ определяющую значение | |р и вычисляем

(106)

Колебания нелинейных систем с одной степенью свободы

263

Для определения: а2 нужно повторить то жепостроение, вновь от­ ложив модуль ^ справа от оси ординат и т. д. (см. рис. 21, б).

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

С и с т е м а с о с т е п е н н ы м з а к о н о м , н е у п р у г о г о

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

К = к у | У \п~

 

Дифференциальное уравнение движения приводится к виду

У + - ^ 1) \ У |Л-1 + Р*У = о,

(107)

где п — показатель

нелинейности; к — постоянный

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

Вместо точного решения этого сложного уравнения приближенно при­ нимают

у = а $ т (р* + ф)

(108)

и находят переменную амплитуду колебании а = а (/) из дифференци-

щие значения:

п . .

0

0,5

1

1.5

2,0

2,5

ф (л)

2,000

1,750

1,571

1,437

1,333

1,249

п . .

3,0

4,0

5,0

6,0

7,0

 

ф (л)

1,178

1,087

0,982

0,914

0,857

 

Уравнение (109) следует из энергетического соотношения: работа силы трения К за один цикл равна уменьшению энергии системы за этот цикл. При п Ф 1 решение дифференциального уравнения имеет вид

а

а0

(1 1 1 )

(п— \)к{ра0)п 1Ф{п) ( пт

264Основы теории колебаний механических систем

Вчастности, при п = 2 затухание колебаний следует гиперболиче­ скому закону

а

До

 

( 112)

■ ,

4/грдр

г

 

“*■

З т л

 

 

С и с т е м а с к у л о н о в ы м

т р е н и е м .

Дифференциальное

уравнение движения приводится к виду

 

 

0 + Л ± - § - =

О,

<ПЗ)

где # — сила кулонова трения.

График движения в этом случае состоит из отрезков синусоид, имеющих одинаковый период, но различные амплитуды (рис. 23, о, фазовая траектория показана на рис. 23, б).

Связь между двумя последовательными максимальными отклоне­ ниями щ и а 1, разделенными интервалами времени, равными полу-

Х* + Т

Т п

периоду = — , имеет вид

*Р

а

! = —щ +

2Р,

(114)

где величина

/? _

*

 

 

(115)

 

тр2

с

 

 

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

силой трения 7?. Этот же результат следует из формулы (111) при

п = 0.

Последовательность амплитуд образует арифметическую прогрес­

сию и огибающие графика движения имеют вид прямых

 

Д = ± (о 0 — - ^ г - ) .

(116)

Как только отклонение а станет меньше, чем {5, движение прекра­ щается.

Колебания нелинейных систем с одной степенью свободы

265

С и с т е м а с н е л и н е й н ы м в н у т р е н н и м т р е н и е м . Внутреннее трение характеризуется рассеянием энергии Ф за один цикл, для многих материалов и целых конструкций можно принять, что ве­ личина У не зависит от частоты процесса и определяется формулой

У =

ааЛ+1,

(117)

где а , п — постоянные системы; а — амплитуда цикла.

В частности, при п = 0 и а =

4/? имеем V =

4/?а; это соответствует

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

добно дифференциальному уравнению

движения

(109)

а +

0.

(118)

где с — коэффициент жесткости системы; Т — период колебаний (при­ ближенно равный периоду свободных колебаний недемпфированной

системы). Решением дифференциального уравнения (118)

при п =^= 1

будет

 

 

°о

 

 

 

а

 

 

(119)

 

 

(п— \)а(

лП— 1

 

 

 

+

------“о

 

 

 

 

 

В

частности, при п ■- 0

и а =

(кулоново трение)

огибающая

имеет

вид прямой

Ш

 

40/

 

 

а — а0 -

= а0

( 120)

 

 

сТ

 

 

 

При п = 2 (что приближенно соответствует закономерностям вну­ треннего трения для многих сталей) огибающая имеет вид кривой ан-

перболического

типа

До

 

 

 

а =

 

( 121)

 

! +

•*-’

 

 

^

сТ

 

В случае п =

1 дифференциальное уравнение (118)

имеет решение

 

 

_

(122)

 

а — а^г

и колебания убывают по экспоненциальному закону (как и в случае

вязкого трения).

Обработка опытных виброграмм на основе выражения (119) позволяет определить значения п и а, необходимые для расчета вынужденных ко­

лебаний.

В системах с умеренным нелинейным трением логарифмический де­

кремент определяется

соотношениями

 

6 =

1п Л ( _ *

=

(123)

 

Д/+1

ДI

ЫЦ шах

где

и я/+1 — две последовательные амплитуды процесса свободных

266

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

затухающих колебаний; Да* и Д/7, — уменьшение амплитуды и потен­ циальной энергии за один /-и цикл; П[ тах — максимальная потенциаль­ ная энергия системы в *-м цикле. В отличие от линейных систем лога­ рифмический декремент колебаний не является постоянным для данной системы, а зависит от амплитуды колебаний.

Вынужденные колебания (случай гармонического возмущения). При умеренном нелинейном демпфировании пользуются линеаризацией сил трения и приходят к дифференциальному уравнению (20). Коэф­ фициент к (или п) эквивалентного линейного трения определяют из условия равенства энергии, рассеиваемой за один цикл в нелинейном (заменяемом) и линейном (заменяющем) элементах трения, при этом коэффициент оказывается зависящим от частоты и амплитуды колеба­ ний (табл. 17).

17. Коэффициенты к и п

при замене нелинейного трения

 

экви валентнымлинейным

 

 

 

Закон нелинейного трения

Эквивалентные коэффициенты

 

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

Квадратичное трение

.

ВВао)

 

40аш

Я = ±

 

*‘ = Ч е г - : п' = - т г

 

 

 

 

 

 

Нелнпсйно-вязкое трение

,

2*Ф (л) (аш)л—1

 

 

 

я

 

Л = ку 1у |Л-1

 

 

АФ(л> (ай))л—1

 

 

 

л. = -------------------

 

 

 

*

 

пт

 

Кулоново трепне

к

■ л

 

2*

± я

 

*

 

*

яа(0

'

птаа*

Внутреннее нелинейное трение

и

аал—1 _ _

аал—1

ЧГ = аал-И

 

 

лсо

П*

2пат

Выражения для п* следует подставить в формулу (67):

/т _______________Уст_____________

(124)

/

о ^2 + 4ш2л2------ ----------(а,

ш)

 

 

Это соотношение следует рассматривать как уравнение относительно неизвестной амплитуды вынужденных колебаний. В табл. 18 приведены найденные таким способом резонансные амплитуды при ш = р.

Колебания нелинейных систем с одной степенью свободы

267

18. Резонансные амплитуды при нелинейном трении

 

Закон нелинейного трения *

Резон

некая амплитуда

 

Квадратичное трение

1 - | /

ЗПР0

 

Р

У

 

 

 

Нелинейно-вязкое трение

1

ЛЯ0

 

р

у

2кФ (л)

 

 

 

Кулоново трен и

При малом трении

 

 

со

 

Внутреннее нелинейное трение

 

л

?

 

 

 

 

* См. табл. 17.

Фрикционные автоколебательные системы

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

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

чения о0.

19. Фрикционные автоколебания

Схема

1

1ЛЛЛММЧ111

^ — г 1

2

3

! Ц м м м /ууу|

|

Возможные автоколебания

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

Вертикальные автоколебания упруго закрепленного тела, прижатого к вращающемуся диску

Горизонтальные автоколебания те­ ла около режима равномерного дви­

жения

268 Основы теории колебаний механических систем

Существенным условием появления фрикционных автоколебаний

сИ?

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

<3 0 (см. табл. 4, сухое трение, нижний график). Отличительной чертой автоколебании является независимость их амплитуды и частоты от на­ чальных условий. На рис. 24 показаны огибающие кривые у = у {() при различных начальных условиях.

Дифференциальное уравнение возмущенного движения имеет вид

а{(П

 

ту +

ДЯ к » У) +

СУ = 0;

(125)

 

 

здесь

перемещение

у

отсчиты­

 

 

вают от невозмущенного уровня,

 

 

т. е. от

положения

равновесия

 

 

схем

/

и 2 табл.

19 или от дви­

 

 

жущегося со скоростью

о0 на­

П р о ц е сс и сгп а н о б л сн и я У т а н о б и б ш и й с я

чала

координат

для

схемы 3

 

проц есс

Рис.

24

табл.

19.

(125)

разность

Д# (ь>о« У)—К

 

В

уравнении

(уо—У) представляет

собой приращение силы

трения, возникающее из-за изменения скорости скольжения; на падаю­ щем участке характеристики эта разность отрицательна ^отрицатель­ ное» демпфирование, способствующее раскачке колебаний).

Точное решение уравнения (119), как правило, затруднительно. Для упрощения решения в зависимости от параметров системы

пользуются одним из двух приближенных способов.

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

Я(о0) — я (Но —

-у ) + Ч/ = 0. (126)

Закон автоколебания такой системы существенно отличается

от гармонического, и возможны интервалы полного сцепления (отсут­ ствия скольжения, рис. 25); такие автоколебания называют разрыв­ ными (из-за разрывов скорости), или релаксационными.

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

< /= асо зр /.

(127)

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

р

 

Г ДК у (II = 0,

(128)

о

Колебания нелинейных систем с одной степенью свободы

269

выражающее равенство нулю работы силы трения за один цикл коле­ баний.

В уравнение (128) вместо скорости у следует подставить ее выраже­

ние,

вытекающее из

соотношения (127):

 

 

 

 

у = ар 51П р1.

 

(129)

Пример Б. Определить амплитуду автоколебаний

в случае,

когда ' при

V > 0

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

трения имеет вид, показанный

на рис.

26;

 

 

 

 

(130)

причем Я, и у# — сила трения и скорость, соответствующие точке минимума характеристики. В этом случае

Я (у„) =ЗЯ, X

„3

Я(у0 у) = 3 Я' X

Уо— У , (Ур — У)8

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

аЯ = 3я .

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

Р

1Н - 4Ь *-ь У 41 = 0

или после подстановки выражения (120)

ар 51П* ф — агр251П4 ф <*ф = 0.

0 *

Отсюда амплитуда автоколебаний

( 131)

При этом наибольшая скорость автоколебаний

270

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

не должна превосходить номинальную скорость скольжения о0 [при Утак —

= у0 произойдет остановка и процесс уже недопустимо описывать зависи­ мостью (130)]; поэтому условия существования квазилинейных автоколебаний

имеют вид двух неравенств

Г

 

 

V -

<V0 < О*,

 

 

 

0,895о, < ов <

(132)

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

Л1

у У

У

0

Рис. 27

Свойства автоколебательных систем общего вида могут быть исследованы по фазовому портрету системы путем изучения характера фазовых траекто­

рий на плоскости у, у.

Для автономной нелинейной систе­ мы, движение которой описывается

Удифференциальным уравнением

У+ Р2У= I {у> У),

(133)

фазовые траектории являются интегральными кривыми дифференци­ ального уравнения первого порядка

Ау

Н у. у) — р*у

(134)

 

 

у

Для графо-аналитического построения фазовых траекторий удобен делыпа-метод, согласно которому траектории строят на фазовой пло­ скости у, V, причем

V = — •

(135)

Р

Далее вводят функцию

(136)

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

 

V2 +

+ б)а = сопз1.

(137)

Это уравнение

описывает

окружность, центр которой расположен

в точке у = —б;

V = 0. Построение фазовой

траектории начинают

с точки с координатами у 0=

у (0), \ 0=

определяемыми началь­

ными условиями. В соответствии с выражением (136) по этим значениям вычисляют б0 = б о, V0); найденное значение определяет центр ок­ ружности, описываемой уравнением (137). Далее по ходу часовой стрелки проводят малую дугу окружности, начинающуюся в точке Уо» (рис. 27). С чертежа снимают новые значения ух, вновь под-