CC BY
219
Автоколебания в одномерных упругих системах с трением
д.т.н. проф. Пожалостин А.А., к.ф.-м.н. доц. Паншина А.В.
МГТУ им. Н.Э. Баумана 8-(499)-263-63-75 panalv@mail.rH Аннотация. Исследуется возможность появления автоколебаний в случае поперечных колебаний упругих однородных балок при наличии в системе сухого трения. Поставленная задача решается на примере балки, подвешенной одним концом на вращающемся валу, и на примере балки, закрепленной шарнирно на своих концах. Применяется метод замены исходной системы эквивалентной механической колебательной системой.
Ключевые слова: балка, сухое трение, упругие колебания, форма колебаний, автоколебания.
Известно [1], что в маятнике Фроуда-Жуковского при определенных условиях возможны нарастающие угловые колебания (автоколебания) и установление периодического колебательного движения с амплитудой Лак. В данной статье показано, что режим автоколебаний возможен в случае поперечных колебаний упругих балок при наличии в системе сухого трения.
Приняты основные допущения:
• балка прямолинейная, однородная, изгиб балки прямой;
• колебания предполагаются малыми;
• материал балки идеально упругий, подчиняется закону Гука;
• сечения балки при ее колебаниях остаются плоскими;
• инерция вращения сечений не учитывается.
Рассмотрим случай 1: изгибные колебания консольной балки. Пусть однородная балка массой m и длиной l подвешена верхним концом на валу О, который вращается с постоянной угловой скоростью Q (рисунок 1).
/
Рисунок 1. Балка на вращающемся валу Рисунок 2. Зависимость момента сухого
трения от угловой скорости
При движении балки на нее действует момент от сил сухого трения. Зависимость момента сухого трения Мтр (ю), действующего на консоль при вращении вала в центре О, от
угловой скорости скольжения СО принята такой же, как и в [1] (рисунок 2). Угловая скорость скольжения СО = С1 — ф, где ф - угол поворота балки как абсолютно твердого тела. Зависимость Мтр (ю) имеет падающий участок АВ . Величина угловой скорости О выбирается на
участке АВ , чтобы продемонстрировать эффект упругих автоколебаний балки.
Дифференциальное уравнение свободных изгибных колебаний однородной балки имеет вид [2]:
ех4 +т° елА
Примем граничные условия закрепления балки в следующем виде:
у(о,г) = о, у ' (о,г) = о,
у "(¡,1) = о, у'"(\,г) = о.
+ = 0 . (1)
(2)
„ , еу ,, е2у ,,, е3у
Здесь производные у = , у , у = —3
ех ех ех3
Решение, согласно [2], ищем в виде:
у = /(х) ■ э(г). (3)
Из (1) получим дифференциальное уравнение для формы колебаний функции /(х) :
Е1о/" + По®2/ = 0. (4)
Представим решение уравнения (4) в следующем виде:
/ (х) = (Хх) + С2Т (Хх) + Сзи(1х) + С^(Хх), (5)
где: £, Т, и, V - функции Крылова [1].
Граничные условия (2) для формы колебаний /(х) [2]:
/(0) = 0, /"(0) = 0, Г (I) = 0, /" (I) = 0. (6)
Удовлетворяя первым двум соотношениям граничных условий (6), получим, что С1 = С2 = 0. Поэтому искомое решение (5) примет вид:
/ (х) = си+C4V. (7)
Из третьего соотношения граничных условий (6) следует, что С4 = -С3 ^^. Поэтому форма колебаний I -го тона:
Дх) = СзХи&х) - ^Т^М)], (8)
Т (Хг1)
где: С31 - произвольная константа.
Из третьего и четвертого соотношений граничных условий (6) получаем трансцендентное уравнение для нахождения собственных значений X1 (I = 1,2,...) краевой задачи и частот
^2
свободных колебаний балки юг- = Х^^
V ^ о
ч2
& о
82(Х/) - Т (XI) ■V (XI) = 0. (9)
Известно [2], что функции / (х) удовлетворяют условию ортогональности:
I
\mofr (х)1} (х)йх = 0 (г ф 7). (10)
0
С учетом условий ортогональности (10) построим эквивалентную маятниковую систему [3], [4] (рисунок 3) (I = 1,2,...).
При построении эквивалентной маятниковой системы (механического аналога) постулируется равенство частот собственных колебаний исходной и эквивалентной систем.
Выпишем кинетическую энергию исходной системы с учетом условий ортогонально-
сти (10) форм колебаний и приравняем ее кинетической энергии эквивалентной системы [2]:
/
1 ' 00 1 00 т=- =%' (п)
2 о 1=1 2 1=1
I
где приведенная масса = |цо/г2 (х)ёх.
Рисунок 3. Эквивалентная маятниковая система
Пусть временная функция из (3) имеет вид = , тогда приведенный момент
инерции эквивалентной системы (маятника) равен Joj = ш°!2, где I - длина балки.
Потенциальная энергия системы:
1 ¥
п = 2 £
2 1=1
с ь2.
I
где: сг =|Ы0(/")2ёх, а крутильная жесткость эквивалентной системы (маятника) равна
с = ■ — о О
со = 0,12 .
Здесь принято во внимание второе условие ортогональности форм колебаний:
I
| г'; ах = о (1 * ;). (12)
о
Проверим равенство частот ю, собственных колебаний упругой балки и эквивалентной системы (механического аналога):
ю2 =
с° с.
¿о* т°
В дальнейшем будем учитывать только первый тон колебаний (, = 1). Дифференциальное уравнение колебаний эквивалентной системы (механического аналога) (/ = 1) имеет вид [1]:
= (13)
Характеристика момента сухого трения Мтр (ю) представлена на рисунке 2. Угловую
скорость вращения вала О выбираем на падающем участке.
Рассмотрим случай 2: однородная балка массой ш и длиной I шарнирно закреплена на концах О и А (рисунок 4).
Балка расположена горизонтально и опирается в середине на шероховатую горизон-
тальную плоскость, вращающуюся с угловой скоростью О. Тогда в середине (х = 2 ) балка имеет участок, где сила трения Етр убывает при возрастании скорости скольжения V.
Рисунок 4. Шарнирно закрепленная балка Рисунок 5. Эквивалентная
система с грузом
Граничные условия закрепления концов балки для формы колебаний f (x) имеют вид
[1]:
f (0) = f (l) = 0, f'(l) = f'(l) = 0. (14)
Собственная форма fi(x) представляется в виде:
f (x) = A sin^, i = 1,2,... (15)
Система с грузом (i = 1,2,...), эквивалентная рассмотренной горизонтальной упругой балке, представлена на рисунке 5.
Для номера i = 1, как и в первом случае, будем иметь приведенные массу и жесткость для эквивалентной системы:
i i nx , l
0 С • 2 nx j l m7 = luo sin — dx = m —,
1 { l 2
I (16)
СО = ] Е.10 (Р )4вт2 ^х = Е,1о (Р )4.
0 II I
Дифференциальное уравнение для координаты ^ груза в эквивалентной системе имеет
(17)
У 8Л(5*1) ртр Ь &1
Здесь обобщенная сила Q1 = —-=-2-.
2). И так как А(1-
Ql = Ртр (V).
Уравнение (17), аналогичное уравнению (13), принимает вид:
(18)
Введем обозначение:
) в случае 1.
д(*) = 1 /А о
^(л) в случае 2.
Так как возможное перемещение ф 0, то Q1 = Fmp (v)f1 (—). И так как f (—) = 1, то
Тогда дифференциальное уравнение для координаты q(t) имеет вид:
2 №)
2 + = (19)
здесь т2 - - Т в случае 1. с Г-0 в случае 1. _{мтр(в случае 1.
3деСЬ: ^ _ т' Т К в слУчае 2. - _{<% в случае 2.' Р(у) _\РщРМ в случае 2.
Разложим функцию ¥(V) в ряд по степеням V в первом и во втором случаях, соответственно. В результате уравнение (19) запишем в виде [1]:
'4 + со^ = -а 4 + уд3. (20)
Используя метод Ван-дер Поля [5] и полагая q — А(t)cos(Ш + в) , получим из (20) для A(t) «уравнение установления» в следующем виде [1]:
_ йА .3 2 л3 ^
2--аА + -ую2 А3 — 0. (21)
йг 4
Проинтегрируем уравнение (21) при начальном условии q(0) — qo, затем возьмем предел при t ® ¥. В результате получим амплитуду автоколебаний установившегося периодического колебательного режима [1]:
А _ 2 Га
ак \ /-> ,
сох у 3у
С
где _у!I •
В [1] доказана устойчивость периодического режима, то есть при qo < Аак, — > 0, а
dt
при qo > Аак, — < 0. Таким образом, это утверждение доказано. dt
Заключение
В работе показана возможность существования автоколебаний при упругих изгибных
колебаниях балок при наличии сухого трения.
Литература
1. Стрелков С.П. Введение в теорию колебаний. СПб.: Лань, 2005. 438 с.
2. Тимошенко С.П. Колебания в инженерном деле. Москва: КомКнига, 2006. 439 с.
3. Шиманский Ю.А. Динамический расчет судовых конструкций. Л.: Судпромгиз, 1969. 408 с.
4. Пожалостин А.А., Кулешов Б.Г., Паншина А.В. Колебания упругих одномерных систем с трением//Инженерный журнал: наука и инновации, 2013, вып. 12. URL: http://engiournal.ru/catalog/eng/teormech/1136.html (дата обращения 15.01.2014).
5. Лампер Р.Е. Введение в теорию нелинейных колебаний авиаконструкций. М.: Машиностроение, 1985. 88 с.
