Колебания нелинейных систем с одной степенью свободы |
261 |
в которой уг (/), уг (/), — должным образом выбранные функции |
вре |
мени; а1г а2, |
— параметры, определяемые |
из системы уравнений |
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ъ определяющую значение | |р и вычисляем
Колебания нелинейных систем с одной степенью свободы |
263 |
Для определения: а2 нужно повторить то жепостроение, вновь от ложив модуль ^ справа от оси ординат и т. д. (см. рис. 21, б).
По этим данным можно построить огибающие кривых затухающих колебаний, если известен период колебаний; можно принять, что он ие отличается от периода свободных колебаний недемпфированной системы (рис. 22).
С и с т е м а с о с т е п е н н ы м з а к о н о м , н е у п р у г о г о
с о п р о т и в л е н и я |
К = к у | У \п~ |
|
Дифференциальное уравнение движения приводится к виду |
У + - ^ 1) \ У |Л-1 + Р*У = о, |
(107) |
где п — показатель |
нелинейности; к — постоянный |
коэффициент. |
Вместо точного решения этого сложного уравнения приближенно при нимают
и находят переменную амплитуду колебании а = а (/) из дифференци-
щие значения:
п . . |
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 решение дифференциального уравнения имеет вид
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 |
и а = |
4К (кулоново трение) |
огибающая |
имеет |
вид прямой |
Ш |
|
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 |
|
|
|
л. = ------------------- |
|
|
|
* |
|
пт |
|
Кулоново трепне |
к |
4Я |
■ л |
|
2* |
± я |
|
* |
|
* |
яа(0 |
' |
птаа* |
Внутреннее нелинейное трение |
и |
аал—1 _ _ |
аал—1 |
ЧГ = аал-И |
|
|
лсо |
’ |
П* |
2пат |
Выражения для п* следует подставить в формулу (67): |
/т _______________Уст_____________ |
(124) |
/ |
о -Я^2 + 4ш2л2------ ----------(а, |
ш) |
|
|
Это соотношение следует рассматривать как уравнение относительно неизвестной амплитуды вынужденных колебаний. В табл. 18 приведены найденные таким способом резонансные амплитуды при ш = р.
Колебания нелинейных систем с одной степенью свободы |
267 |
18. Резонансные амплитуды при нелинейном трении |
|
Закон нелинейного трения * |
Резон |
некая амплитуда |
|
Квадратичное трение |
1 - | / |
ЗПР0 |
|
Р |
У |
8Р |
|
|
|
Нелинейно-вязкое трение |
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. При умеренно крутом падении характеристики и упругой связи небольшой жесткости принимают, что автоколебания носят гармониче ский характер и происходят с частотой р свободных колебаний той же системы (но без трения)
В этом случае систему называют квазилинейной. Для определения амплитуды установившихся автоколебаний используют уравнение энергетического баланса
2л
Колебания нелинейных систем с одной степенью свободы |
269 |
выражающее равенство нулю работы силы трения за один цикл коле баний.
В уравнение (128) вместо скорости у следует подставить ее выраже
ние, |
вытекающее из |
соотношения (127): |
|
|
|
|
у = — ар 51П р1. |
|
(129) |
Пример Б. Определить амплитуду автоколебаний |
в случае, |
когда ' при |
V > 0 |
характеристика |
трения имеет вид, показанный |
на рис. |
26; |
|
|
|
|
(130) |
причем Я, и у# — сила трения и скорость, соответствующие точке минимума характеристики. В этом случае
Я (у„) =ЗЯ, X
„3
Я(у0 — у) = 3 Я' X
Уо— У , (Ур — У)8
следовательно,
аЯ = 3я .
Подставляя это выражение в уравнение (128). имеем
2я
Р
1Н - 4Ь *-ь У 41 = 0
или после подстановки выражения (120)
2л
ар 51П* ф — агр251П4 ф <*ф = 0.
0 *
Отсюда амплитуда автоколебаний
( 131)
При этом наибольшая скорость автоколебаний
270 |
Основы теории колебаний механических систем |
не должна превосходить номинальную скорость скольжения о0 [при Утак —
= у0 произойдет остановка и процесс уже недопустимо описывать зависи мостью (130)]; поэтому условия существования квазилинейных автоколебаний
имеют вид двух неравенств
Г |
|
|
V - |
<V0 < О*, |
|
|
|
0,895о, < ов < |
(132) |
При нарушении правого неравенства автоколебания невозможны, при нарушении левого неравенства они будут носить релаксационный характер.
Свойства автоколебательных систем общего вида могут быть исследованы по фазовому портрету системы путем изучения характера фазовых траекто
рий на плоскости у, у.
Для автономной нелинейной систе мы, движение которой описывается
Удифференциальным уравнением
фазовые траектории являются интегральными кривыми дифференци ального уравнения первого порядка
у
Для графо-аналитического построения фазовых траекторий удобен делыпа-метод, согласно которому траектории строят на фазовой пло скости у, V, причем
Р
Далее вводят функцию
(136)
которую в малых интервалах времени можно считать постоянной, при этом дифференциальное уравнение (134) может быть проинтегрировано
|
V2 + |
(У+ б)а = сопз1. |
(137) |
Это уравнение |
описывает |
окружность, центр которой расположен |
в точке у = —б; |
V = 0. Построение фазовой |
траектории начинают |
с точки с координатами у 0= |
у (0), \ 0= |
определяемыми началь |
ными условиями. В соответствии с выражением (136) по этим значениям вычисляют б0 = б {уо, V0); найденное значение определяет центр ок ружности, описываемой уравнением (137). Далее по ходу часовой стрелки проводят малую дугу окружности, начинающуюся в точке Уо» (рис. 27). С чертежа снимают новые значения ух, вновь под-