В базу знаний

Колебания механических систем

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

12 мин чтения От разработчика

Материалы к статье

К разделу

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

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

5.gif,(1)

где 5.1.gif – постоянные коэффициенты;
5.2.gif – обобщенная координата перемещения системы;
5.3.gif – обобщенная масса, это может быть масса, момент инерции или иная характеристика;
5.4.gif – коэффициент упругой силы или коэффициент жесткости;
5.5.gif – коэффициент диссипативной или демпфирующей силы.

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

5.6.gif,(2)

где 5.7.gif – собственная угловая частота колебаний;
5.8.gif – коэффициент затухания колебаний.

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

5.9.gif.(3)

Периодическая частота колебаний определяется следующим образом:

5.10.gif.

Безразмерный коэффициент затухания колебаний определяется следующим образом:

5.11.gif.(4)

При 5.12.gif движение системы представляет собой затухающие гармонические колебания. При 5.13.gif движение системы перестает быть колебательным. Значение 5.14.gif называется критическим затуханием, поэтому значение коэффициента 5.8.gif иногда называют долей критического затухания. На рисунке 1 представлены графики колебаний при различных значениях коэффициента затухания 5.8.gif.

В некоторых случаях уравнение (1) принято записывать в виде:

 5.16.gif,

где 5.17.gif ; 5.18.gif,

или в виде:

5.19.gif.(5)

где 5.20.gif – коэффициент демпфирования.

При 5.21.gif общее решение дифференциального уравнения (5) имеет вид:

5.22.gif,

где 5.23.gif – постоянные, определяемые из начальных условий.

5.1.png Рис. 1. Характер колебаний при различных значениях коэффициента затухания (доли критического затухания)

Логарифмический декремент и другие характеристики колебаний

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

 5.24.gif,

где 5.25.gif – амплитуды двух последовательных колебаний.

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

5.26.gif,

где 5.27.gif – число периодов колебаний.

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

При слабом затухании разность 5.28.gif мала по сравнению с 5.29.gif, и тогда для оценки значения логарифмического декремента можно использовать следующую формулу:

5.30.gif.(6)

Для систем с линейными характеристиками типа (1) логарифмический декремент является постоянной величиной, его значение равно:

5.31.gif.(7)

При слабом затухании колебаний:

5.32.gif.(8)

В программном комплексе EULER для задания диссипативных свойств некоторых силовых элементов используется относительный коэффициент демпфирования:

5.33.gif,(9)

где 5.34.gif – коэффициент диссипативной силы линеаризованной модели; 5.4.gif– коэффициент упругой силы линеаризованной модели или коэффициент жесткости.

Используя (3), (4), (5), (8) выражение (9) можно записать в виде:

5.35.gif,(10)

где 5.36.gif – периодическая частота первой формы собственных колебаний.

Затухание колебаний при различных типах диссипативных сил

В реальных механических системах логарифмический декремент колебаний не всегда имеет постоянное значение. Это означает, что даже при линейных упругих силах диссипативные силы не являются линейными. Нелинейные диссипативные силы могут описываться различными зависимостями. Для многих практических случаев можно принять описание диссипативной силы в виде следующей степенной зависимости:

5.37.gif,(11)

где 5.38.gif – постоянный коэффициент; 5.39.gif – показатель степени зависимости силы от скорости.

Рассмотрим характер колебаний систем с линейными упругими силами и некоторыми типичными диссипативными силами.

При 5.40.gif зависимость (11) описывает линейное или вязкое трение, которое является наиболее распространенным. Как уже отмечалось, в процессе колебаний логарифмический декремент в этом случае имеет постоянное значение, а огибающая линия кривой процесса колебаний описывается экспонентой:

5.41.gif,(12)

где 5.42.gif – начальная ордината огибающей линии (начальная амплитуда);
5.43.gif
5.44.gif – период колебаний;
5.45.gif – время.

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

При 5.46.gif зависимость (11) описывает кулоновское или сухое трение, которое не зависит от величины скорости. Логарифмический декремент в этом случае увеличивается в процессе уменьшения амплитуды колебаний, а огибающая линия кривой процесса колебаний описывается линейной зависимостью:

5.47.gif(13)

где 5.51.gif – постоянный коэффициент.

При 5.49.gif зависимость (11) описывает квадратичное трение. Примером квадратичного трения является аэродинамическое сопротивление. Логарифмический декремент в этом случае уменьшается в процессе уменьшения амплитуды колебаний, а огибающая линия кривой процесса колебаний описывается гиперболой:

5.50.gif,(14)

где 5.51.gif – постоянный коэффициент.

На рисунках 2 и 3 представлен вид огибающих линий кривой процесса колебаний и зависимости логарифмического декремента от амплитуды колебаний при различных типах диссипативных сил.

5.2.jpgРис. 2. Огибающие линии кривой процесса колебаний при различных типах диссипативных сил (5.46.gif – сухое трение; 5.40.gif – вязкое трение; 5.49.gif – квадратичное трение)

5.3.jpgРис. 3. Зависимость логарифмического декремента от амплитуды колебаний при различных типах диссипативных сил (5.46.gif – сухое трение; 5.40.gif – вязкое трение; 5.49.gif – квадратичное трение)

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

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

Рассеяние энергии в материалах

В процессе колебаний механических конструкций некоторая часть энергии рассеивается в конструкционных материалах. Эти свойства связаны с несовершенной упругостью материалов и зависят от уровня напряжений, вида напряженного состояния, температуры, структуры материала, а также от размеров образцов, схемы их закрепления, типа колебаний (например, продольные или поперечные колебания стержня) и других факторов. В таблице 1 по данным [3, стр. 128] приведены значения логарифмических декрементов некоторых конструкционных материалов при средних уровнях нагружения. Данные для алюминиевых сплавов и стеклопластика получены при изгибных колебаниях свободно подвешенных образцов-стержней кругового сечения длиной до нескольких метров, для остальных материалов – при колебаниях небольших призматических образцов в условиях чистого изгиба.

Материал  3.gif
Алюминиевые сплавы 0.0002
Титановые сплавы 0.0002 – 0.0003
Легированные стали 0.001 – 0.01
Стеклопластик 0.01

Таблица 1. Логарифмические декременты конструкционных материалов (5.52.gif)

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

Конструкционное рассеяние энергии

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

По данным [3, стр. 122] значения логарифмических декрементов для конструкций ракет и космических аппаратов, как правило, не выходят за пределы диапазона 0.015-0.35. При этом декременты, величины которых меньше 0.025, а также больше 0.2, встречаются редко. Зависимость декрементов от амплитуды колебаний может иметь различный характер. В одних случаях декременты возрастают при увеличении амплитуд колебаний, в других – остаются практически постоянными. Иногда наблюдается также уменьшение декрементов с увеличением амплитуд колебаний.

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

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

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

5.4.pngРис. 4. Процесс колебаний с различными значениями логарифмического декремента

Система подрессоривания автомобиля

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

  • частота собственных колебаний подрессоренной массы автомобиля должна находиться в интервале 1.0-2.0 [Hz];
  • коэффициент затухания колебаний подрессоренной массы автомобиля должна находиться в интервале 0.2-0.3 (доля критического затухания).

Вынужденные колебания линейных и нелинейных систем

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

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

Введем обозначения:

5.53.gif,

где 5.54.gif – частота установившихся вынужденных колебаний системы равная частоте внешнего возмущения;
5.55.gif – собственная частота системы (для нелинейных систем при малых амплитудах).

5.56.gif

где 5.57.gif – амплитуда установившихся вынужденных колебаний системы;
5.58.gif – статическое значение амплитуды (амплитуда установившихся вынужденных колебаний при бесконечно малой частоте внешнего возмущения).

На рисунке 5 представлены амплитудно-частотные характеристики механической системы (1) с линейными характеристиками упругости и демпфирования в зависимости от коэффициента затухания колебаний (доли критического затухания). Для систем с малым демпфированием характерно наличие резонанса – резкого увеличения амплитуды колебаний при частоте внешнего возмущения близкой к собственной частоте системы.

5.4.jpgРис. 5. АЧХ линейной системы

В нелинейных системах возникают различные сложные эффекты: зависимость частоты свободных колебаний от амплитуды; неоднозначная зависимость амплитуды вынужденных колебаний от частоты внешнего возмущения; наряду с основными колебаниями (частотой 5.54.gif) возникновение супергармонических (частотой 5.59.gif) и субгармонических (частотой 5.60.gif) колебаний; хаотизация колебаний и другие эффекты.

В качестве простого примера сложностей динамического поведения нелинейной системы можно продемонстрировать вынужденные колебания системы с одной степенью свободы, линейным демпфированием, но нелинейной упругой силой. На рисунке 6 представлена амплитудно-частотная характеристика системы с жесткой характеристикой упругой силы. Штрихпунктирной линией обозначена зависимость амплитуды свободных колебаний от их частоты. При малом демпфировании АЧХ нелинейной системы демонстрирует наличие как бы «деформированного» или «кривого» резонанса. В диапазоне частот возбуждения 5.61.gif амплитуда вынужденных колебаний определена неоднозначно, причем режим между точками АЧХ 5.62.gif и 5.63.gif является неустойчивым. При постепенном увеличении частоты от достаточно малых значений при 5.64.gif происходит резкий срыв амплитуды колебаний, а при постепенном уменьшении частоты от достаточно больших значений при 5.64.gif происходит резкое увеличение амплитуды колебаний. На рисунке 7 представлена аналогичная амплитудно-частотная характеристика системы с мягкой характеристикой упругой силы.

5.5.jpgРис. 6. АЧХ с жесткой характеристикой упругой силы

5.6.jpgРис. 7. АЧХ с мягкой характеристикой упругой силы

В файле «resonans.elr» содержится проект для программного ком-плекса EULER, который демонстрирует наличие «кривого» резонанса в простейшей системе «груз на вибростенде». Груз соединен с колеблющей-ся платформой вибростенда упругим и демпфирующим силовыми элементами. Эти элементы воспроизводят типовые характеристики автомобильной подвески. В качестве исходных данных для расчета в проекте задаются следующие параметры:

//
scalar type_K "тип характеристики жесткости: 1 - линейная; 2 - нелинейная"=2;
scalar type_W "тип изменения вынуждающей частоты: 1 - нарастание; 2 - убывание"=2;
scalar T "общее время процесса"=100[s];
scalar A "амплитуда вынуждающих колебаний"=50[mm];
scalar f_max "максимальная вынуждающая частота"=5[Hz]; scalar mG "масса груза"=500[kg];
scalar fG "собственная частота колебаний груза"=1[Hz]; scalar ksiG "коэффициент затухания колебаний груза"=0.2;

//

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

a) увеличение частоты возбуждения       b) уменьшение частоты возбуждения 2026-10-16_16-29-04.pngРис. 8. Колебания груза с нелинейной характеристикой упругой силы

Приведенный пример показывает, что в диапазоне частот 1.5-3.5 Hz установившиеся колебания груза, имитирующего корпус автомобиля на штатной подвеске, могут иметь сильно отличающиеся амплитуды, которые зависят от истории прихода к текущей частоте колебаний.

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

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

5.65.gif,(15)

где 5.66.gif – номер собственной формы колебаний стержня;
5.67.gif – длина стержня;
5.68.gif – изгибная жесткость стержня;
5.69.gif – погонная масса стержня;
5.70.gif – 5.71.gif-е собственное значение.

Собственные значения являются корнями частотного уравнения стержня. Собственные значения и уравнения форм свободных поперечных колебаний однородного стержня для различных схем закрепления пред-ставлены в таблице 2 по материалам [4] стр. 334, [5] стр. 200-206.

Схема закрепления стержня 5.71.gif Собственное значение 5.70.gif Вид/уравнение формы свободных колебаний
5.9.jpg
1
2
3
…
5.71.gif
3.142
6.283
9.425
 
5.72.gif
2026-10-19_14-37-23.png
5.73.gif
5.9.1.jpg 
1
2
…
5.71.gif
3.927
7.069
 
5.74.gif
2026-10-19_14-40-49.png
5.75.gif
5.13.jpg
1
2
…
5.71.gif
4.73
7.853
 
5.76.gif
2026-10-19_14-41-27.png
5.77.gif
5.16.jpg
1
2
…
5.71.gif
1.875
4.694
 
5.78.gif
2026-10-19_14-44-53.png
5.79.gif
5.19.jpg
1
2
3
…
5.71.gif
4.73
7.853
11.00
 
5.76.gif
2026-10-19_14-45-11.png
5.80.gif

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

5.81.gif,(16)

где 5.82.gif – номера собственных форм колебаний для сторон пластины (длины, ширины);
5.83.gif – размеры сторон пластины (длина, ширина);
5.84.gif – толщина пластины;
5.85.gif – модуль упругости материала пластины;
5.86.gif – коэффициент Пуассона материала пластины;
5.87.gif – плотность материала пластины.

В соответствии с (16) первая собственная угловая частота квадратной пластины равна:

5.88.gif.

Литература

  1. Тимошенко С. П. Колебания в инженерном деле. – М.: Наука, 1967. 444 с.
  2. Пановко Я. Г. Введение в теорию механических колебаний. – М.: Наука. Гл. ред. физ.-мат. лит., 1991. 256 с.
  3. Микишев Г. Н. Экспериментальные методы в динамике космиче-ских аппаратов. – М.: Машиностроение, 1978. 248 с.
  4. Машиностроение. Энциклопедия / Ред. совет: К.В. Фролов (пред.) и др. – М.: Машиностроение, Динамика и прочность машин. Тео-рия механизмов и машин. Том I-3. В 2-х кн. Кн. 1 / К.С. Колесников, Д.А. Александров, В.К. Асташев и др.; Под общ. ред. К.С. Колесникова. 1994. 534 с.
  5. Ильин М.М., Колесников К.С., Саратов Ю.С. Теория колебаний. – М.: Издательство МГТУ имени Н.Э. Баумана, 2003. 271 с.

Материалы к статье