УДК 531.36
ОСНОВНОЕ ДВИЖЕНИЕ ДВУХМАЯТНИКОВОГО АВТОБАЛАНСИРА НА ГИБКОМ ВАЛУ С УПРУГИМИ ОПОРАМИ
В.А. Дубовик, Г.Р. Зиякаев
Томский политехнический университет Е^У: [email protected]
Получены условия балансировки маятниками неуравновешенного ротора на гибком валу, вращающегося в упруго-податливых опорах.
Ключевые слова:
Ротор, маятниковый автобалансир, дисбаланс, устойчивость, критическая угловая скорость, уравнения движения системы. Key words:
Rotor, balancing device, unbalance, dynamic stability, critical speed, equations of motion.
Изучению динамики маятниковых и шаровых автобалансиров посвящены работы [1-4]. В [5] исследуется в первом приближении устойчивость основного движения маятникового автобалансира. В этих работах считается, что опоры вала, на котором крепится ротор с балансировочным устройством, являются абсолютно жесткими. Практически опоры вала, как правило, оказываются упруго деформируемыми, что приводит к изменению частоты колебаний вала [6] и, как следствие этого, условий балансировки. Ниже изучается влияние упругой податливости и массы опор вала с неуравновешенным ротором на основное движение маятникового автобалансира. Считается, что на ротор в месте крепления его к валу помимо упругой силы действуют внутреннее трение, пропорциональное скорости деформации [6], и сила внешнего сопротивления, пропорциональная скорости точки крепления. Целью работы является при сформулированных выше положениях установление условий существования основного движения маятникового автобалансира, т. е. движения, при котором отсутствуют прогиб вала и смещение опор (рис. 1, а).
На (рис. 1, б) приведена схема возмущенного движения системы. Ротор, центр масс которого составляет с его геометрическим центром С эксцентриситет СС'=е, закреплен посредине вала в точке С так, что он может двигаться только в плоскости, перпендикулярной оси опор АВ - хОу. При вращении вала вокруг неподвижной оси Ог с постоянной угловой скоростью со он прогибается на величину О'С, а ось опор АВ смещается на - ОО Предполагается, что ось вращения Ог вертикальна, либо силы инерции намного больше сил тяжести. Это позволяет пренебречь влиянием последних на движение системы.
За обобщенные координаты принимаем х,у и х0,у0 перемещения центра ротора и опор вала относительно неподвижной системы координат хОу, а также уьу2 - углы поворота маятников относительно вектора дисбаланса ротора СС' (рис. 1, б). Массу ротора обозначим через тр; массу опоры - т0; массу маятника - т; моменты инерции относительно центральных осей, параллельных оси Ог ротора и маятника соответственно - 1р и /; изгибную жесткость вала - с; жесткость опоры - с0; длину маятника - I.
Упругие силы, действующие на опору и на ро-т—р со _стороны вала, равны ¥С0=-с0-ОО' и
- =-с0-О'С; сила вязкого сопротивления среды ¥т=-к1-УС; моменты сил вязкого сопротивления, препятствующие .повороту м.аятников относительно вала Ыс1=-к1у1, Мс2=-к2уг- Скорость деформации вала в точке С равна разности скоростей точек Си О'в системе координат £,О-, вращающейся вместе с валом, -гС0=УгС-Ую=УС-Уо-тхО,С. Тогда сила внутрен-его трения, приложенная к ротору в точке С равна ¥т^в=-кУгС0< [6], а её проекции на неподвижные оси будут
РТЕх =~к[Х ~К + Ш(У - Уо )]>
РТЕу = -[У - Уо - Ш(Х - Хо )]-
Здесь кък2 - коэффициенты вязкого трения ротора, маятника; к - коэффициент внутреннего трения; точка сверху означает производную по времени.
В принятых обозначениях кинетическая, потенциальная энергии и диссипативная функция рассматриваемой механической системы имеют вид
Т = 1 тр [(х - гю$т ю?)2 + (у + ею сов ю?)2] +
+ Т ^рШ + 2Т то(х0 + у1) +
где
М = mp + 2m, h = h21 ( J + ml2),
+1 ¿i 2
Гp 2 ['x -1 (o + y) sin(ot + y¡)]2 + [ m I ^ +
|+[y +l (a + y ) cos(ot + y ¡)]2 J +J (o + y¡ )2
П = 2 C[( -V - x. )2 + (y - y. )■] + 2Í c. (x.2 + y.2);
ф = i k |[ V - x +"(У - y)]’ +1 +
2 [+[У - y. -®(V - )]2 J
+т hi( x 2 + y2) + 2 h2 X y <22 2 /=1
Подставляя эти выражения в уравнения Лагранжа 2-го рода, приходим к дифференциальным уравнениям движения системы
Mx + c(x - xo) + \Х + k(X - Xo) + km(y - yo) -
2
-mlX [y sin(mt + y/) + (m + y )2 cos(mt + y)] =
/=i
= mpem2 cos mt; (1)
My + c(У - У.) + hiy + k(y- yo) - km(x - x.) +
2
+mlXy cos(mt + y/) - (m + y )2 sin(mt + y )] =
/=i
= mpem2 sin mt; (2)
2m.H. + (2co + c)x. + kx. - cx - kx - km(y - y.) = 0; (3)
2m.yo+(2c. +c)y. + ky.- cy - ky + km(x - x)=0; (4) y i + hy / = ¡лф (x,y, y), i =12, (5)
(7)
ц= ml / (J + ml2) < 1,
4>i ( x,У, У ) = x sin(ot + У ) - У cos(ot + y). (6)
Для достижения поставленной цели воспользуемся теорией синхронизации объектов с равномерными вращательными движениями [3]. Принимаем л за малый параметр. Соответствующая (1-6) порождающая система уравнений (л=0) допускает семейство стационарных периодических решений вида у = const1 =а1, у2 = const2 =а2; x(t) = U cos cot - V sin cot, y(t) = U sin ot + V cos ot; xo (t) = Uo cos ot - Vo sin ot, yo (t) = Uo sin ot + Vo cos ot.
Здесь введены обозначения
U = aR(al,a1) + bl (apa2),
V = al (a1,a2) - bR(a1,a2); a = a1a2 / (a^ + b22), b = a1b2 / (a^ + b22); a1 =-2moa>2 +2co +c, a2 =a1(c-Mo2)-c2; b2 = a1h1o;
U0 = aRaa) + boI (a1,a2),
Vo = aoI(a1,a2) - boR(a1,a2);
ao = ca2 / (a^ + b22), bo = cb2 / (a^ + b22);
R(a1,a2) = mpeo2 + mlo2(cosa1 + cosa2),
1 (a1,a2) = mlo2(sina1 + sina2). (8)
Эти решения зависят от двух постоянных a1 и a2, которые фиксируют маятники относительно вектора дисбаланса ротора. Для определения этих постоянных составим основные уравнения. Решения (7) подставляем в (6), усредняем по времени за период 2я/ю, полученные порождающие функции приравниваем нулю
p (a1,a2) =
2 [[aR(a1,a2) + bI(a1,a2)]sina--]
= -a i ^ = 0
[-'aI(a1, a2 ) - bR(a1, a2 )] cos a ¡
¡ = 1,2. (9)
Система уравнений (9) имеет решение
a1 = a1 =n-a0, a2 =a2 =-(л-a0),
mpe
cosa0 =^— = X < 1, 0 2ml
(10)
при котором R(a1*,a2*)=I(a1*,a2*)=0 и, как следует из (7), x(t)=y(t)=0, xo(t)=yo(t)=0, т. е. движение автобалансира, соответствующее этому решению, является основным. Решение (10) можно получить из равновесия сил инерции при основном движении.
Достаточные условия асимптотической устойчивости решения (10) имеют вид [3]
(
Ь1(а1,а2) = -
8Р дР2
—L + —-8а1 да
>0,
2 J
Ь2(а1,а2) =
8Р2
> 0.
(11)
дР1 дР2 дР1 ^ 2
да1 да2 да2 да1
Подставляя в (11) производные от порождающих функций (9), для решения (10) получаем
¿1(а1*,а:*) = -2т1а>4 а,
¿2(а;,а2‘) = (т/ю4)2(а2 + 62)4Я2(1 -Я2). (12)
В соответствии с (10) величина Х2(а1*,а2*) всегда положительна. Условие (12) для Х1(а1*,а2*) приводит к требованию а<0 или с учетом обозначений (8) к неравенству, записанному в безразмерных переменных
2 (Р) = -у2р3 + [2у(1 + 7)+у2]Р2 -
-[у(1 + 27) + (1 + 7)2]Р+7(1 + 7) < 0, (13)
со > (01,
при наличии одной критической скорости и
в случае трех критических скоростей.
Если две последние скорости совпадают о2=о3=о*, то движение будет устойчивым при о>о1, кроме случая о=о*.
Аналогичные результаты для жестких опор получены в [5].
На рис. 2 в плоскости (о/о0, 2co/c) изображены графики критических скоростей в зависимости от жесткости опор y=2cjc при v=1. Здесь область А между кривыми 1 и 2 и область В над кривой 3 являются областями устойчивого движения независимо от дисбаланса ротора.
ю/га0
¿т„ ю м ю 2с I
где у =—0, Р =-= —-, 7 = —°, ю0 = V с / М
М с ю0 с
- собственная частота системы с неподвижными маятниками.
Решение (7) и условие устойчивости (13) не содержат коэффициента к, отсюда рассматриваемое движение не зависит от внутренних сил, что очевидно, т. к. при стационарном вращении вала скорость деформации его равна нулю.
При отсутствии массы опор (у=0) условие устойчивости основного движения будет
о>4Ш. (14)
Здесь с*=2с„с/(2с„+с) - приведенная жесткость системы. В случае с„^да (опора становится абсолютно жесткой) с*^с, и неравенство (14) переходит в условие устойчивости маятникового автобалансира на гибком валу с жесткими опорами [2].
Исследование функции г(Р) показывает, что г(0)=7(7+1)>0, г(Р)^-да при Р^+да. В интервале [0,+да) в зависимости от значений параметров у, 7 она имеет либо один корень Р1, либо три действительных корня Д, Р2, Р3, обозначенные в порядке возрастания, при этом два последних могут совпадать. Неравенство (13) выполняется, что соответствует устойчивости движения, в первом случае при Р>Д, а во втором - при Р__3<Р2 и _>_3. Угловые скорости ш=Р0, ю=Р2ю0, ш=\РШо разделяют прямую 0<ю<+да на устойчивые и неустойчивые области, т. е. являются критическими. Отсюда условия существования и устойчивости основного движения маятникового автобалансира имеют вид
2
JL ©
0
'XF
О
2 4 2^/с
Рис. 2. Графики критических скоростей в зависимости от жесткости опор и области устойчивости движения
Из проведенных расчетов следует, что при у<1 область А увеличивается, а область В уменьшается; при у>1 - наоборот. С ростом параметра у область неустойчивости между кривыми 2 и 3 увеличивается. При 7>5920 для у=0,5; 7>8090 для у=1; 7>10927 для у=2 существует только одна критическая скорость Ш~ю0. В этом случае движение устойчиво при всех ю, удовлетворяющих неравенству ю>ю0.
Выводы
Уравновешивание ротора на гибком валу с податливыми опорами малой жесткости маятниковым автобалансиром возможно при вращении вала с угловой скоростью в диапазоне критических скоростей от первой до второй, а также свыше третьей. При большой жесткости опор уравновешивание происходит при угловой скорости, превышающей собственную частоту устройства. За счет выбора жесткости опор можно добиться включения автобалансира на угловой скорости, значительно меньшей, чем собственная частота.
СПИСОК ЛИТЕРАТУРЫ
1. Гусаров А.А. Автобалансирующие устройства прямого действия. - М.: Наука, 2002. - 119 с.
2. Пановко Я.Г Основы прикладной теории колебаний и удара. -Л.: Машиностроение, 1976. - 320 с.
3. Блехман И.И. Синхронизация в природе и технике. - М.: Наука, 1981. - 352 с.
4. Нестеренко В.П. Автоматическая балансировка роторов приборов и машин со многими степенями свободы. - Томск:
Изд-во Том. ун-та, 1985. - 84 с.
5. Филимонихин Г.Б. К устойчивости основного движения двухмаятникового автобалансира // Доповид НАН Укра'ши. -1996. - № 8. - С. 74-78.
6. Диментберг Ф.М. Изгибные колебания вращающихся валов. -М.: Изд-во АН СССР, 1959. - 246 с.
Поступила 18.03.2010 г.