close

Вход

Забыли?

вход по аккаунту

?

Влияние зазора в подшипниках крепления на конические прецессии неуравновешенного ротора.

код для вставкиСкачать
УДК 534
Вестник СПбГУ. Сер. 1, 2005, вып. 1
И. А. Пасынкова
ВЛИЯНИЕ ЗАЗОРА В ПОДШИПНИКАХ КРЕПЛЕНИЯ
НА КОНИЧЕСКИЕ ПРЕЦЕССИИ
НЕУРАВНОВЕШЕННОГО РОТОРА
Было установлено [1–4], что динамически (моментно) неуравновешенный жесткий
ротор, укрепленный вертикально в нелинейных упругих опорах так, что его центр масс
находится на равном расстоянии от опор, может иметь режим вращения, который представляет симметричную коническую прецессию. При этом движении центр масс остается неподвижным, а ось вращения ротора зачерчивает в пространстве поверхность
прямого кругового конуса. Модель ротора представляла собой механическую систему с
четырьмя степенями свободы. Наличие зазора в подшипнике не учитывалось, т. е. предполагалось, что ось ротора во все время движения находится в контакте с поверхностью
подшипника. В [5] рассматривалась упрощенная модель гибкого статически неуравновешенного ротора с двумя степенями свободы. Исследовано прецессионное движение с
учетом зазора в подшипниках. Предложенные ранее (см. [1, 3]) идеи и методы в данной
работе развиваются на случай крепления ротора в подшипниках с зазором. Проводится
сравнение с моделью ротора без учета зазора.
Рассматривается динамически симметричный ротор массы M и длины L, имеющий
момент инерции относительно оси симметрии A и экваториальный момент инерции B.
Под динамической (моментной) неуравновешенностью понимается угол δ между осью
динамической симметрии и осью вращения ротора (осью геометрической симметрии
ротора). Предполагается, что податливость нелинейных упругих опор обладает центральной симметрией, в силу этого реакция опоры имеет только радиальную составляющую. Рассматривается существенно нелинейная зависимость реакции опоры от величины смятия в точке контакта, задаваемая формулой Герца. Кроме того, учитывается
наличие зазора в подшипнике, величина которого равна ρ0 . Предполагается, что ротор
вращается с постоянной угловой скоростью ω, а смещение вдоль оси ротора считается пренебрежимо малым. В качестве обобщенных координат выбираются координаты
центров опор в плоскостях, перпендикулярных оси вращения ротора в положении его
равновесия.
В подшипнике с зазором (радиальным люфтом) ротор может совершать движения
как находясь в контакте с подшипником, так и не имея такого контакта, поэтому сила
реакции опоры предполагается кусочно-непрерывной и задается в виде
⎧
Si
⎨
c0 (|Si | − ρ0 )3/2
, если |Si | > ρ0 ,
Pi =
(1)
|Si |
⎩ 0,
если |S | ρ ,
i
0
где Si — смещение центра соответствующей опоры от равновесного положения (коплексная переменная), а c0 — вещественная постоянная, характеризующая нелинейную упругость опоры. Кроме того, учитываются силы сопротивления, пропорциональные скорости движения центра опоры: Ri = −μ̃vi , i = 1, 2.
c
И. А. Пасынкова, 2005
103
Система дифференциальных уравнений движения ротора будет иметь вид (см. [6, 7])
M
(S̈1 + S̈2 ) + (R1 + R2 ) + (P1 + P2 ) = 0,
2
L2
L2
B(S̈2 − S̈1 ) − iωA(Ṡ2 − Ṡ1 ) +
(R2 − R1 ) +
(P2 − P1 ) =
2
2
= (B − A)ω 2 Lδ exp (iωt).
(2)
Последовательно
перейдем в системе (2) к безразмерному времени τ = ω0 t, где
√
ω02 = 2c0 Lδ/M, и безразмерным относительным переменным si (i = 1, 2). По существу
задачи величины |Si | и зазор ρ0 малы, поэтому в качестве характерного линейного
размера возьмем величину Lδ, введем безразмерный зазор ρ = ρ0 /(Lδ) и безразмерные
относительные смещения центров опор si = Si /ρ0 .
Как и раньше, введем безразмерные параметры
λ=
A
,
B
k=
M L2
M L2
=
,
4(B − A)
4B(1 − λ)
Ω=
ω
,
ω0
μ=
2μ̃
.
M ω0
(3)
Параметры k и λ всегда удовлетворяют соотношению k(1 − λ) > 0. Следовательно,
если ротор представляет собой динамически сжатое тело, то λ > 1 и k < 0. Если ротор
является динамически вытянутым телом, то λ < 1 и k > 0.
Система (2) в новых переменных примет вид
s̈1 + s̈2 + μ(ṡ1 + ṡ2 ) + (p1 + p2 ) = 0,
s̈2 − s̈1 + (μk(1 − λ) − iΩλ)(ṡ2 − ṡ1 ) + k(1 − λ)(p2 − p1 ) =
1
= (1 − λ)Ω2 exp (iΩτ ),
ρ
(4)
где точкой обозначено дифференцирование по безразмерному времени τ. Безразмерные
силы pi определяются по формулам:
%
si
, если |si | > 1,
ρ1/2 (|si | − 1)3/2
|si |
(5)
pi =
0,
если |si | 1.
Система (4) допускает решение вида
sj = rj exp (i φj ) exp (i Ω τ ),
j = 1, 2,
(6)
где rj , φj — вещественные постоянные и rj > 0. Это решение представляет собой прямую
синхронную прецессию ротора, которая может быть гиперболоидальной, конической
или цилиндрической. Характер прецессии можно установить по соотношению между
фазами φj . Для гиперболоидальной прецессии фазы могут быть произвольными, при
равных фазах гиперболоидальная прецессия вырождается либо в цилиндрическую (если r1 = r2 ), либо в коническую (если r1 = r2 ) с вершиной конуса вне опор. Наконец,
если фазы различаются на π, то будет иметь место коническая прецессия с вершиной
конуса между опор, при этом соотношение между радиусами может быть любым.
Рассмотрим сначала систему без сопротивления, соответствующую движению с контактом, чтобы определить ситуации, «опасные» с точки зрения возникновения резонансов. Система без сопротивления также допускает решения вида (6). Подставим решение
104
(6) в систему (4), положив μ = 0, и разрешим ее относительно «комплексных амплитуд»
rj exp (i φj ):
A2 Ω2
A1 Ω2
, r2 exp (i φ2 ) =
.
(7)
r1 exp (i φ1 ) = −
ρΔ
ρΔ
Здесь приняты следующие обозначения:
Δ = A1 B2 + A2 B1 ,
1/2
Aj =
ρ
(rj − 1)3/2 − Ω2 ,
rj
(8)
1/2
Bj = k
ρ
(rj − 1)3/2 − Ω2 .
rj
(9)
Множество Δ = A1 B2 +A2 B1 = 0 в пространстве переменных (Ω, r1 , r2 ) определяет
резонансное множество или множество нелинейных резонансов. Если рассматривать
симметричные прецессии, когда r1 = r2 = r, то это множество преобразуется к виду
Δ = A1 B1 = (ρ1/2 (r − 1)3/2 − r Ω2 ) (k ρ1/2 (r − 1)3/2 − r Ω2 ) = 0.
(10)
На плоскости (Ω2 , r) это две исходящие из одной точки Ω = 0, r = 1 непересекающиеся
кривые. Вдоль одной из них (A = 0) «резонирует» цилиндрическая прецессия, вдоль
другой (B = 0) — коническая.
Среди решений (6) системы (4) можно выделить решение вида s2 = −s1 = s, что соответствует симметричной конической прецессии ротора, когда его центр масс неподвижен, а центры опор движутся по окружностям одинакового радиуса. При этом первое
уравнение в (4) обращается в тождество, а второе принимает вид
s̈ + (μk(1 − λ) − iΩλ)ṡ + k(1 − λ)p =
1
(1 − λ)Ω2 exp (iΩτ ),
2ρ
(11)
решением которого будет функция s = r exp(iφ) exp(iΩτ ). Здесь r — безразмерная относительная амплитуда, а φ — безразмерная фаза смещения опоры в системе координат,
вращающейся с угловой скоростью ротора.
Значение r = 1 соответствует смещению центра опоры на величину, равную зазору
в подшипнике, и переходу от движения без контакта к движению с контактом, или
наоборот. При этом для движения ротора с контактом в подшипнике мы имеем уравнение, соответствующее ненулевому значению силы p, а для движения без контакта,
когда величина смещения центра опоры меньше величины зазора, сила p = 0. Переход
с одного уравнения на другое происходит при значении |s| = r = 1. После подстановки
решения в (11) и несложных преобразований мы получим зависимости амплитуды и фазы симметричной конической прецессии от частоты, соответствующие двум режимам
движения ротора, с контактом и без контакта. Удобно перейти к новым переменным
(x, y) по формулам x = Ω2 , y = r. Амплитудно- и фазово-частотные характеристики
конической прецессии в переменных (x, y) будут
3
4
x
, y>1
(12)
(xy − k ρ(y − 1)(y − 1))2 + k 2 μ2 xy 2 =
2ρ
√
4
x
y x + k 2 μ2 =
,
y1
(13)
2ρ
√
kμy x
4
tg φ =
.
(14)
xy − k ρ |y − 1||y − 1|
105
При увеличении частоты вращения ротора предельное значение амплитуды y, как
следует из (12) и (13), отлично от нуля и равно y∞ = 1/(2ρ), как для движения с
контактом, так и для движения без контакта. Это означает, что при большой угловой
скорости вращения происходит самоцентрирование ротора, при этом ось динамической
симметрии ротора стремится занять положение оси вращения ротора в равновесии. В
зависимости от величины зазора это движение будет либо с контактом, либо без него.
Если y∞ > 1, т. е. ρ < 1/2, то вращение ротора в состоянии, близком к центрированному, будет происходить с контактом в подшипнике. В противном случае самоцентрирование произойдет в состоянии без контакта. Если вернуться к размерным величинам,
то неравенство ρ < 1/2 означает, что зазор в подшипнике меньше отклонения оси динамической симметрии ротора от оси его вращения в плоскости опоры.
На рис. 1 приведены АЧХ конической прецессии ротора с зазором (толстая линия)
и без зазора (тонкая линия) при следующих значениях параметров: ρ = 0.4, k = 2.0,
μ = 0.09. На правом рисунке представлены те же характеристики в крупном масштабе. Как видно, влияние зазора самым существенным образом сказывается при малых
частотах. Так, если зазор мал (ρ < 1/2), то диапазон частот, в котором возможно движение без контакта, также мал. Координата точки перехода с одного типа движения
на другой определяется по формуле xtr = 4 ρ2 μ2 k 2 /(1 − 4 ρ2 ) и для данных значений параметров xtr = 0.0576. Кроме того, амплитуда прецессии ротора при движении с
контактом существенно больше (левая ветвь АЧХ), хотя максимальное значение амплитуды несколько меньше, чем при вращении ротора в подшипнике без зазора, и «срыв»
происходит при меньших значениях частот. Эта же особенность отмечена в работе [5].
Рис. 1.
В зависимости от параметров ротора возможны от одного до трех различных режима симметричной конической прецессии. Проведем исследование устойчивости по линейному приближению сначала для той части АЧХ, которая соответствует движению
с контактом, т. е. для y > 1. Пусть безразмерной частоте вращения ротора Ω отвечают
значения r, φ безразмерной амплитуды и фазы, причем r > 1. Введем для каждой опоры малые возмущения zi , αi и построим систему в вариациях, сохранив обозначения
x = Ω2 , y = r. Удобно ввести следующие обозначения:
106
L1i
=
L2i
=
L3i
=
+
L4i
=
+
√
√
34
z̈i + μ(żi − y xαi ) − 2y xα̇i + (
ρ(y − 1) − x)zi ,
2
4
√
√
y α̈i + μ(y α̇i + xzi ) + 2 xżi + ( ρ(y − 1)3 − x y)αi ,
√
√
z̈i + μk(1 − λ)(żi − y xαi ) + (λ − 2)y0 xα̇i +
3 4
(1 − λ)( k ρ(y − 1) − x)zi ,
2
√
√
y α̈i + μk(1 − λ)(y α̇i + xzi ) + (2 − λ) xżi +
4
(1 − λ)(k ρ(y − 1)3 − x y)αi .
Тогда уравнения в вариациях запишутся следующим образом:
L11 − L12 = 0,
L21 − L22 = 0,
L31 + L32 = 0,
L41 + L42 = 0.
(15)
Характеристический полином системы (15) имеет восьмую степень, но как и в случае подшипников без зазора, приводится к произведению двух полиномов, каждый из
которых имеет четвертую степень. Опуская промежуточные вычисления, выпишем эти
полиномы:
M 1 = a0 p 4 + a1 p 3 + a2 p 2 + a3 p + a4 ,
a0 = 1,
a1 = 2μ,
4
1 4x y + (5 y − 2) ρ(y − 1) + 2μ2 y ,
a2 =
2y
4
μ a3 =
4x y + (5 y − 2) ρ(y − 1) ,
2y
4
4
1 (2 x − 3 ρ(y − 1))(x y − ρ(y − 1)3 ) + 2μ2 x y .
a4 =
2y
(16)
Здесь l = 1 − λ и l > 0, так как рассматривается динамически вытянутый ротор:
M 2 = b 0 p4 + b 1 p3 + b 2 p2 + b 3 p + b 4 ,
b1 = 2 μ k l,
b0 = 1,
4
1 2(1 + l2 )x y + k l (5 y − 2) ρ(y − 1) + 2μ2 k 2 l2 y ,
b2 =
2y
4
μkl 4 x y + k l (5 y − 2) ρ(y − 1) ,
b3 =
2y
4
4
1 (2 x − 3 k ρ(y − 1))(x y − k ρ(y − 1)3 ) + 2 k 2 μ2 x y .
b4 =
2y
(17)
Как и для ротора без зазора, устойчивость различных режимов конической прецессии определяется знаком коэффициентов a4 и b4 . Все остальные коэффициенты полиномов M1 , M2 и определители Рауса—Гурвица первого порядка всегда положительны.
Определители Рауса—Гурвица второго порядка имеют вид
μ
(4y(x + μ2 ) + c) > 0,
2y
(18)
μkl
(4 l y(x + μ2 k 2 ) + k c) > 0,
2y
(19)
det2 (M1 ) =
det2 (M2 ) =
107
4
где использовано обозначение c = (5 y − 2) ρ (y − 1) и учтено, что y > 1 и l = 1 − λ > 0.
Определители Рауса—Гурвица третьего порядка приводятся соответственно к виду
μ2
(d + 4 y c (4x + μ2 )) > 0,
4y 2
(20)
μ2 k 3 l3
(k l d + 4 y c ((1 + l)2 x + k 2 l2 μ2 )) > 0,
4y 2
(21)
det3 (M1 ) =
det3 (M2 ) =
где использовано обозначение d = ρ(y 3 − 1) + 3ρ(y 2 − 1).
Определители Рауса—Гурвица четвертого порядка равны соответственно det4 (M1 ) =
a4 det3 (M1 ), det4 (M2 ) = b4 det3 (M2 ), т. е. условия Рауса—Гурвица нарушаются только
когда a4 < 0 или b4 < 0.
Рис. 2.
На рис. 2 показаны ассимптотически устойчивые (толстая линия) и неустойчивые
(тонкая линия) участки АЧХ, построенной при тех же значениях параметров. Пунктирными линиями проведены границы выполнения условий устойчивости Рауса—Гурвица
a4 = 0 и b4 = 0, которые представляют собой криволинейные гиперболы. Одна из
асимптот для a4 = 0 — это нелинейный резонанс A1 = 0 из (10), а для b4 = 0 асимптотой
является нелинейный резонанс B1 = 0, т. е. неустойчивость возникает в зоне влияния
нелинейных резонансов. Неустойчивость участков в диапазоне a4 < 0 левой ветви и в
диапазоне b4 < 0 правой ветви АЧХ имеет различный характер. Так, при b4 = 0 (вблизи резонанса для конической прецессии) наблюдается классическое явление «срыва» и
при этом характер прецессии (симметричная коническая) сохраняется. Неустойчивость
при a4 < 0 (вблизи резонанса для цилиндрической прецессии) приводит к качественному изменению характера движения. В этом диапазоне частот вовлекается в движение
центр масс ротора, и прецессия перестает быть конической. Такое явление качественного изменения характера движения носит название «пространственной неустойчивости»
(по терминологии [8]). Для разных значений параметра k взаимное расположение нелинейных резонансов может быть разным, и явление пространственной неустойчивости
может иметь место как в докритическом диапзоне частот (рис. 2), так и в закритическом.
Можно предположить, что коническая прецессия ротора при вращении без контакта
не будет асимптотически устойчивой. И действительно, характеристический полином
108
системы линейного приближения в этом случае имеет вид
P = P1 P2 P3 ,
(22)
где полиномы
P1 = (p2 + x)2 ,
P2 = p2 + 2μ p + x + μ2 ,
P3 = p2 + 2μ k(1 − λ)p + (1 − λ)2 (x + μ2 k 2 )
(23)
имеют чисто мнимые корни.
Рис. 3.
Что касается динамически сжатого ротора (λ > 1), то вращение с контактом возможно только при малых зазорах (ρ < 1/2), причем резонанса не будет ни прикаких значениях частоты, так как tg φ всегда ограничен, что следует из формулы (14).
На рис. 3 приведена АЧХ конической прецессии динамически сжатого ротора при
ρ = 0.4, k = −2.0, μ = 0.09 для значений y > 1. Кривая асимптотически стремится к
своему предельному значению y∞ = 1/(2 ρ). Исследование устойчивости показало, что
при всех частотах имеет место асимптотическая устойчивость, за исключением очень
узкого диапазона малых частот, когда при очень малых значениях ρ может возникнуть
пространственная неустойчивость. Если зазор ρ > 1/2, то ротор при своем вращении не
имеет контакта с подшипником, и такая прецессия не является асимптотически устойчивой.
Аналогичные результаты справедливы для цилиндрической прецессии статически
неуравновешенного ротора.
Summary
I. A. Pasynkova. An influence of bearings clearance on the rigid unbalanced rotor conical precessions.
The existence, characteristics and stability of the conic precessions of a rigid dynamically unbalanced rotor supported by bearings with radial clearance are investigated. It is assumed that the
restoring forces are non-linear and vary according to the Hertz formula and the resistance forces are
linear functions of velocity. The resonance curves are obtained and compared with the characteristics like these in the case with non-clearance bearings. The stability conditions are obtained with
the Hurwitz criteria.
109
Литература
1. Пасынкова И. А. Прецессии жесткого неуравновешенного ротора в нелинейных упругих
опорах // Проблемы и перспективы прецезионной механики и управления в машиностроении:
Материалы международной конференции. Саратов. 1997. С. 83–85.
2. Пасынкова И. А. Гиперболоидальная прецессия ротора в нелинейных упругих опорах //
Вестн. С.-Петерб. ун-та. Сер. 1. 1997. Вып. 4 (№ 22). С. 88–95.
3. Пасынкова И. А. Устойчивость конической прецессии жесткого неуравновешенного ротора // Вестн. С.-Петерб. ун-та. Сер. 1. 1998. Вып. 1 (№ 1). С. 82–86.
4. Архипова И. М., Пасынкова И. А. Исследование прецессионного движения неуравновешенного ротора // В сб.: Вторые Поляховские чтения: Избранные труды. СПб. 2000. С. 65–72.
5. Pascal M. Non linear vibrations of an unbalanced rotor with radial clearance // Третьи
Поляховские чтения: Избранные труды. СПб. 2003. С. 123–131.
6. Тимошенко С. П. Колебания в инженерном деле. М., 1959. 440 с.
7. Кельзон А. С., Циманский Ю. П., Яковлев В. И. Динамика роторов в упругих опорах.
М., 1982. 280 с.
8. Ганиев Р. Ф., Кононенко В. О. Колебания твердых тел. М., 1976. 432 с.
Статья поступила в редакцию 27 мая 2004 г.
Документ
Категория
Без категории
Просмотров
5
Размер файла
251 Кб
Теги
неуравновешенных, крепления, влияние, прецессия, подшипники, конические, ротора, зазоре
1/--страниц
Пожаловаться на содержимое документа