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

Вынужденные колебания неуравновешенного ротора в неизотропных опорах Текст научной статьи по специальности «Физика»

CC BY
462
95
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ДИНАМИКА РОТОРОВ / ROTOR DYNAMICS / ПРЯМАЯ ПРЕЦЕССИЯ / DIRECT WHIRLING / ОБРАТНАЯ ПРЕЦЕССИЯ / REVERSE WHIRLING / АНИЗОТРОПИЯ ОПОР / ANISOTROPY OF BEARINGS

Аннотация научной статьи по физике, автор научной работы — Пасынкова Инна Анатольевна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Вынужденные колебания неуравновешенного ротора в неизотропных опорах»

УДК 517.933:534.1:531.36

Вестник СПбГУ. Сер. 1. Т. 1 (59). 2014. Вып. 2

ВЫНУЖДЕННЫЕ КОЛЕБАНИЯ НЕУРАВНОВЕШЕННОГО РОТОРА В НЕИЗОТРОПНЫХ ОПОРАХ

И. А. Пасынкова

Санкт-Петербургский государственный университет,

Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7/9

В работе исследуется влияние анизотропии упругих опор на вынужденные колебания статически и моментно неуравновешенного ротора, укрепленного как на жёстком, так и на гибком валу. Рассматривается ротор с четырьмя степенями свободы. Предполагается, что вал укреплен в линейно упругих неизотропных опорах. Выписаны дифференциальные уравнения вращения ротора в комплексных переменных и найдено точное решение системы, соответствующее эллиптической синхронной прецессии. Точное решение представляет сумму двух векторов, один из которых параметризует прямую прецессию, а другой — обратную. Построены амплитудно-частотные характеристики прямой и обратной прецессий и эллиптические траектории концов оси ротора. Показано, что в случае неизотропных опор могут быть реализованы как прямая, так и обратная прецессии, а также движение оси сложного типа, когда один конец её движется в прямом направлении, а другой — в обратном. Влияние анизотропии упругих опор проявляется также в сдвиге критических частот в сторону уменьшения и появлении дополнительных критических частот в нижней части спектра, что существенно усложняет динамику высокоскоростного ротора в период выхода на рабочую угловую скорость. Библиогр. 7 назв. Ил. 6.

Ключевые слова: динамика роторов, прямая прецессия, обратная прецессия, анизотропия

опор.

Введение. При исследовании вынужденных колебаний ротора, вызванных его неуравновешенностью, часто возникает необходимость учитывать влияние анизотропии упругих опор. В работе [1] были изучены особенности колебательного движения гибкого вала с диском при различии упругих свойств опор в двух взаимно-перпендикулярных направлениях. Были рассмотрены простейшие системы с двумя степенями свободы при симметричном расположении диска относительно одинаковых опор. Была показана возможность существования обратной прецессии вала, что невозможно в случае изотропных опор. В монографии [2] также исследуется влияние ортотроп-ных упругих опор на вращение ротора с двумя степенями свободы, получено точное решение уравнений, представляющее цилиндрическую эллиптическую прецессию.

В настоящей работе изучаются колебания полностью неуравновешенного ротора с четырьмя степенями свободы в двух случаях, именно укрепленного на жёстком и на линейно упругом валу. В обоих случаях с учётом анизотропии опор составлены уравнения движения в комплексных переменных и найдено точное решение этих уравнений. Отметим, что для жёсткого ротора, укрепленного в линейно-упругих неизотропных опорах и вращающегося с постоянной угловой скоростью, уравнения движения в декартовых координатах и без учёта сил трения получены в [3].

Описание модели жесткого ротора и уравнения движения. Рассматривается жёсткий ротор, укрепленный в двух неподвижных упругих опорах. Ротор установлен вертикально, что позволяет не учитывать влияние силы тяжести.

Введены следующие обозначения: масса и длина ротора соответственно М и Ьг, расстояние между опорами — Ь. Точка С — центр масс ротора. Пусть ротор обладает динамической симметрией и его моменты инерции суть 1р (полярный или осевой) и ^ (трансверсальный или экваториальный). Ротор статически и моментно неуравновешен. На рис. 1 показана модель ротора. Статический эксцентриситет, т. е. смещение

центра масс ротора от его оси вращения, равен е. Моментный эксцентриситет характеризуется двумя углами, 6 и е. Угол 6 — угол между осью динамической симметрии ротора и прямой, параллельной оси вращения, проходящей через центр масс, а угол е — это угол между плоскостью I, проходящей через ось вращения ротора и его центр масс, и плоскостью II, в которой лежит угол 6. Угол е определяет фазовый сдвиг моментного эксцентриситета по отношению к статическому.

Считаем, что ротор приводится во вращение двигателем, способным поддерживать заданный закон изменения угловой скорости ^ = у>(£). Тогда, пренебрегая перемещением ротора вдоль оси вращения, его можно рассматривать как механическую систему с четырьмя степенями свободы.

Обозначим через Q точку пересечения оси вращения ротора и плоскости, перпендикулярной к оси вращения и проходящей через центр масс (рис. 1). Пусть ротор расположен не центрально между опорами, так что расстояние от первой (нижней) опоры до точки Q равно е\ Ь, а расстояние от второй опоры QQ2 равно е2 Ь.

Если ротор укреплен консольно, то е2 < 0, так что всегда выполняется е\ + е2 = 1.

Введем следующие системы координат: Ох у г — инерциальная система координат с осью Ог, совпадающей с направлением оси вращения ротора в его неподвижном состоянии; С С — связанная с ротором система главных центральных осей инерции (рис. 2).

В качестве обобщенных координат можно выбрать либо декартовы координаты точки Q: (xQ, yQ) и углы (а, в), характеризующие направление оси вращения Ql Q2 в пространстве, либо декартовы координаты (х^, у^) (_?' = 1, 2) ]-го конца оси ротора (точки Qj) в плоскости, перпендикулярной оси подшипников. Именно эти координаты используются в дальнейшем.

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

считать, что главные направления обеих опор совпадают, и коэффициенты упругости равны соответственно cxj, cyj ,j = 1, 2.

Абсолютные координаты центра масс C равны

xc = xq + e cos <p, ye = yq + e sin (1)

Координаты (xq, yq, а, в) можно выразить через координаты (xj, yj) (j = 1, 2):

xq = e2 xi + ei x2, yq = e2 yi + ei y2,

. x2 - xi . y2 - yi (2)

i3=—

Решение задачи ведется в геометрически линейной постановке, т. е. справедливы оценки xj/L ^ 1, yj/L ^ 1, и углы (а, в) являются малыми.

Удобно перейти к комплексным переменным Sj и ввести полусуммы и полуразности коэффициентов упругости:

, , __. _ C-xj I Cyj _ Cxj Cyj

= Xj + cmj = ^ I сф' = ^ •

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

Rje) = -/j Sj (j = 1,2). (4)

Применяя теоремы о движении центра масс и теорему моментов, составим уравнения движения ротора в комплексных переменных:

-jSj + /tj Sj + cmjSj + CrfjSj j = Mew2 exp(wi),

j=1,2

y^ (-1 y (jt Sj -ijpùj Sj + L2 ej /TtjS'j + L2 ej (cmjSj + cdjSj= (5) j=1,2

= ( Jt — Jp) LSw2 exp(i (wt — e)),

где S'j — комплексно-сопряжённые переменные.

В системе уравнений (5) перейдем к безразмерному времени т: т = wo t, где wo —характерная угловая скорость, и к безразмерным переменным Sj = Sj/h, где

масштабный коэффициент Н может быть равным либо статическому эксцентриситету е, либо величине Ь 6. Получим

+ Ч + + = ехр(гГ2г),

5 = 1,2

Е + - + - + = (6)

5 = 1,2

= (1 — Х)(Ш2 ехр(г(Пт - £)).

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

■7Р МЬ2 Л Л -Ы ..е_ Щ ,,2 _ . ,2 _ СФ

Тг' к-\у

Jt ¿г(1 — Л) Н Н МШ0 МШ0 МШ0

(7)

Заметим, что всегда справедливо к (1 — Л) > 0.

Прямые и обратные синхронные прецессии. При исследовании вращения ротора в изотропных опорах было показано, что при выполнении определённых условий существует прямая синхронная круговая прецессия ротора как в случае линейных упругих опор, так и в случае нелинейных опор. Покажем, что в рассматриваемом случае существуют как прямые, так и обратные эллиптические синхронные прецессии. Уравнения (6) допускают точные решения вида

= аЛз ехр(Ют) + аГ5 ехр(—Шт), ] = 1, 2. (8)

Здесь аф- и аг^ — комплексные амплитуды точки ( ротора в прямой и обратной прецессии.

После подстановки решения в виде (8) в систему (6) получим два равенства вида

где

А+ ехр(Шт) + А- ехр(—Шт) = 0, В+ ехр(Шт) + В- ехр(—Шт) = 0,

А+ = ^ ((ш^ - ез^-П2 + + - с^П2,

5 = 1,2

А~ = ~ - щЩагз + ^фвф) ,

5=1,2

В+ = ^ (-1-У + г^О.) - П2)а,ф + - ¿2£12 ехр(-ге),

5=1,2

В- = (-1)5 + - ЩЩ - ■

(9)

(10)

(11)

Равенства (9) должны тождественно удовлетворяться при любом значении т. Функции ехр(Шт) и ехр(—¿Пт) линейно независимы, следовательно, коэффициенты А+, В+, А-,В- должны равняться нулю. Заменим выражения А-, В- на комплексно-сопряжённые и приравняем их, а также величины А+,В+ нулю, получим неоднородную систему линейных алгебраических уравнений относительно величин аф-, агу:

Рг = (, (12)

где г = (а^ъ «гЪ вгг) , Я = (<¿1^2 ,0, ¿2^1 ехр(—ге), 0) , а элементы матрицы Р равны (а = 1, 2)

Р1,з = »Шз - ез-з+ Р1,2+з = »|,

Р2,з = »ф, Р2,2+з = »Шз - ез-зО2 + г^еп,

рз,з = (-1)3(кез(»Ш + - П2), рз,2+ = (-1)3^»2, (13)

= ("1 р^+о = (-1У (ке^ш^ + г^П) - ) •

Для вычисления критических частот, на которых резонируют прямая или обратная прецессии ротора, приравняем нулю определитель матрицы Р при отсутствии сопротивления. Пусть = 0. Для произвольных значений остальных параметров определитель системы (12) имеет громоздкое выражение, и найти критические (резонансные) угловые скорости возможно только численно. Сделаем некоторые упрощающие предположения. Пусть опоры одинаковые, т. е. »Ш1 = »Ш2 = »ш, »41 = »42 = »4, и ротор укреплен центрально, т. е. его центр масс находится точно посередине между опорами и е1 = е2 = 1/2. Тогда уравнение Д = 0, где Д — определитель матрицы Р в отсутствии сопротивления, распадается на три уравнения:

X - 2(»ш + »2) = о, X - 2(»ш - »2) = о,

(14)

4(1 + А)Х2 - 4к»ШX + (1 - А)к2 (»Ш - »4) =0.

Если в этих уравнениях от величин »ш, »4 вернуться к частотам, связанным с главными направлениями »х, »у, то получим систему

X - 2»2 = 0, X - 2»,2 = о,

х у (15)

4(1 + А)X2 - 2к(»2 + »„Д + (1 - А)к2»2»2 =0.

Как видно из этих уравнений, для динамически вытянутого («длинного») ротора, когда А < 1, существуют четыре критические угловые скорости, а для динамически сжатого («короткого») ротора при А > 1 — три. Этот результат согласуется со случаем вращения ротора в одинаковых изотропных опорах [4].

Решение системы (12) позволит определить комплексные амплитуды прямой и обратной прецессии точки (. Абсолютные величины этих амплитуд представляют амплитудно-частотные характеристики (АЧХ) соответствующих движений. На рис. 3 показаны АЧХ прямой (жирная сплошная линия) и обратной (штриховая линия) прецессий точек (1 и (2 для следующих значений параметров: е1 = е2 = 1/2, е = п/3, А = 0.7, к = 3.8, ^ = 2, й2 = 0.5, »2 = 1.2, »^ = 0.4, ^ = ^ = 0.05. Тонкими сплошными линиями показаны линейные резонансы, соответствующие критическим частотам, которые при данных параметрах равны X! = 0.19, X2 = 0.80, Xз = 1.60, X4 = 2.40. Отметим, что при выбранном малом внешнем сопротивлении доминируют резонансы, соответствующие собственным частотам в двух главных направлениях. На тех частотах, где амплитуда прямой прецессии превосходит амплитуду обратной прецессии, соответствующая точка ( движется по эллиптической траектории в направлении вращения ротора, если же больше амплитуда обратной прецессии, то результирующее движение происходит противоположно вращению ротора.

На рис. 4 показаны траектории точки (1 (жирная сплошная линия) и точки (2 (штриховая линия) для разных значений частот. На частоте X =1.5 имеет место

обратная прецессия, а на частоте X = 3.5 — прямая. Это хорошо видно, если использовать функцию анимации в математическом пакете Maple.

цессия.

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

|Si| = + Щ - cose, |S'2| = + Щ + 4Ьлк2 cose, (16)

а амплитуды обратной прецессии равны нулю.

Уравнения движения ротора, укрепленного на гибком валу. Рассмотрим ротор с четырьмя степенями свободы, который представляет собой абсолютно твёрдое динамически симметричное тело, прикрепленное к упругому валу, массой которого

можно пренебречь по сравнению с массой тела. Гибкий вал укреплен в линейно-упругих анизотропных опорах. Модель ротора приведена на рис. 5.

Рис. 5. Модель гибкого ротора с 4 степенями свободы, укрепленного в упругих опорах.

Введем следующие системы координат: Ох у г — инерциальная система координат с осью Ог, совпадающей с направлением оси вращения ротора в его неподвижном состоянии; система координат Q£nC, жёстко связанная с ротором и с осью QZ, направленной вдоль касательной к изогнутой оси вала.

Положение ротора можно охарактеризовать с помощью 8 параметров: (х, у) — декартовы координаты точки Q; (а, в) —углы, определяющие направление оси QZ; (х3-, уз) —декартовы координаты точки Qj (] = 1, 2).

Кинетическая энергия ротора с точностью до линейных членов относительно параметров дисбаланса е, 6 и квадратичных членов относительно малых величин а, в и их производных а, в может быть записана в виде

Т = ^М (ж2 + у2 + еш(усов(шг) - хвт(шг))) + - 2ш$а)+

+ ^ (а2 + /З2) + (1р - 6ш (а эт^ - е) - ¡3 сов(шг - е)^ . (17)

Потенциальная энергия изогнутого вала и упругих сил со стороны опор может быть записана в виде

Пв = ^ сп ((ж - ж0)2 + (у - г/о)2) + ^ с22 ((а - а0)2 + (/3 - /Зо)2) +

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

+ С12 ((ж - х0)(а - а0) + (у - у0)(в - /Зо)) + \ ^ (сХэж2 + сю-г/2). (18)

2 3=1,2

Здесь введены обозначения C = { c;m} (l,m = 1, 2) — матрица жесткости упругого вала, закреплённого в жестких подшипниках; (хо, yо) —декартовы координаты точки Qo; углы (ао, во) определяют направление прямой QiQ2. Величины (xo, yo, «о, во) характеризуют перемещение ротора как твердого тела (рис. 5) и вычисляются по формулам (2).

Уравнения Лагранжа II рода относительно комплексных переменных (3) могут быть записаны в форме

МБ + 11еБ + сп(Б - во) + 012(7 - 70) = Меи2 ехр(шЬ),

■¿а - ирШ'у + ^еЬ2л/ + С12(Б - Бо) + 022(7 - 7о) = (19)

= (■ - ■р)6и2 exp(г(wt - е)), (сце2 - ^Ч (Б - Б0) + (с 12е2 - (7 - 70) = с^Бх + с,цБъ

(20)

(сце 1 + (Б - Б0) + (с12ел + (7 - 70) = ст2Б2 + сЛ2Б2.

Алгебраические уравнения (20) отражают баланс упругих сил со стороны изогнутого вала и восстанавливающих сил со стороны деформированных опор. Эти уравнения можно рассматривать также как уравнения связей между координатами Б, 7 и Б1, Б2. Любой набор переменных (Били (Б1, Б2) может быть выбран в качестве обобщенных координат. Пусть (Б1, Б2) будут обобщенными координатами.

Уравнения (20) являются линейными относительно (Б, 7), точное решение системы будет

5= 53 ((ез-5 + К1уСту)Бу +К1уСфБу) ,

5 = 1,2

7= Е(~1у((1 + к^т+К2^)' (21)

5=1,2 ^ '

к15 = сИ + (-1)5 с12 е5 Ь к25 = с22 е5 Ь + (-1У с12'

где С2 = {с*1т} — матрица податливости. Компоненты этой матрицы сЦт зависят от способа крепления вала. Для различных видов опор эти функции можно найти во многих монографиях, например [5].

Подставим решение (21) в уравнения (19), затем перейдем к безразмерному времени т = ио Ь и безразмерным переменным 35 = Б5/к. Получим следующие безразмерные уравнения:

53 (Мт ({е3-з + + егу^в,-) + + ш^) = дл П2 ехр(гПт),

5=1,2

53 (-1)^ (<?т ((1 + (¡23^)83 + (723^3) + к I + ш2а^з)) = (22)

5=1,2

= Ы2 П2 ехр(г (Пт - е)).

Здесь

¿2 ё ¿2 ё МЛ*) = + М*) = ТТ(*) + (-ША + ^/Ме) — (*). (23)

ат2 ат ат2 ат

Основные обозначения соответствуют формулам (7) и

<715 = к 15 Ми>2, 025 = К25ЬМ и>2. (24)

Прямые и обратные синхронные прецессии. Как и в случае жесткого ротора, уравнения (22) допускают точное решение вида (8). Проведем вычисления, как

для жесткого ротора, получим линейную систему алгебраических уравнений вида (12) для определения критических частот и амплитуд прямой и обратной прецессий. Получим следующие элементы матрицы Р будут (_?' = 1, 2):

= + (ез-з + »2тз)(-П +

Р1,2+з=¿3(1+(-^2+ Р2,з = ¿2 (1 + (-П2 + 3)),

Р2,2+з = »2тз + (е3-3 + 0-11^3-)(-П +

Рз,3 = (-1)3 (кез »2т3 + (1 + 023 ¿т з)(-02 + гкез 3)), Рз,2+з = (-1)3 ¿ф (ке3 +023 (-П2 + гке3 ^П)),

р4,з = ("1)^3 + ^23 ("+ гке^р^ ) , (26)

РА,2+3 = (-1)' (ке^ + (1 + а2з"шз) + ■

Рассмотрим в качестве примера модель турбомолекулярного насоса, которая представляет собой твёрдое тело, консольно укрепленное на гибком валу. В работе [6] предполагалось, что вал установлен в жёстких опорах, и были найдены две критические частоты: ¿1 = 10496 об/мин, ¿2 = 67080 об/мин. В работе [7] было показано, что учёт упругости опор в случае их изотропности приводит к значительному сдвигу критических частот в сторону уменьшения. Так, для той же модели при коэффициенте упругости для стального подшипника соответствующие частоты равнялись ¿1 = 7832 об/мин, ¿й2 = 17926 об/мин. В этой работе предположим, что вал укреплен в анизотропных одинаковых опорах, в которых жёсткости в двух главных направлениях различаются в два раза. Возьмём в качестве наибольшей жёсткости её значение из [7]: стах = 1.08 • 107 н/м. Наименьшая жёсткость будет ст¡п = 0.54 • 107 н/м. Безразмерные параметры этой модели примут значения е1 = 1.2, е2 = -0.2, е = п/3, Л = 0.64, к = 12.0, Н1 = 1, Н2 = 0.25, ¿т = 1, ¿2 = 0.33, ^ = ^ = 0.05. Для этих параметров найдены четыре критических значения: = 0.173, Х2 = 0.299, Х3 = 1.126, Х4 = 2.360. Этим значениям соответствуют угловые скорости ротора ¿1 = 4616 об/мин, ¿2 = 7010 об/мин, ¿3 = 11779 об/мин, ¿4 = 17052 об/мин. Для выбранных значений параметров на рис. 6, а и Ь показаны АЧХ прямой (жирная сплошная линия) и обратной (штриховая линия) прецессий для точек Ql и Q2. Тонкими сплошными линиями показаны линейные резонансы, соответствующие критическим частотам. Отметим, что прямая прецессия резонирует на частотах, близких к частотам в случае изотропных опор. На тех частотах, где амплитуда прямой прецессии превосходит амплитуду обратной прецессии, соответствующий шип движется по эллиптической траектории в направлении вращения ротора, если же больше амплитуда обратной прецессии, то результирующее движение шипа происходит противоположно вращению ротора.

На рис. 6, с показаны траектории точки Ql (жирная сплошная линия) и точки Q2 (штриховая линия) для значения частоты X = 1.2. При этом наблюдается интересная ситуация, когда точка Ql движется по эллипсу в направлении вращения ротора, а точка Q2 —в противоположном.

Заключение. На основании полученных результатов можно сделать вывод о значительном влиянии анизотропии опор на динамику неуравновешенных роторов,

(25)

Рис. 6. АЧХ прямой и обратной прецессии: а — точки Ql; Ь — точки с— траектории точек Ql и Q2.

KI

1,5 1,0 0,5

\ /

\ / \ /

\ / \j

Д / \

/ \

j \

/

1 2

У

1,5-

0,5

-i,o\

\

4-0,5

N

\ -1,0-1,5-

3 х

С

ч

\

\

Ж

укрепленных как на жёстком, так и на гибком валу. Происходит сдвиг критических частот в сторону уменьшения и появляются дополнительные критические частоты, причём в нижней части спектра. Это будет существенно влиять на выход роторной машины на рабочую угловую скорость. Кроме того, могут быть реализованы прямая и обратная прецессии, а также движение оси сложного типа, когда один конец её движется в прямом направлении, а другой — в обратном. Все это подтверждает высокую практическую значимость учёта даже незначительной анизотропии опор.

Литература

1. Диментберг Ф. М. Изгибные колебания вращающихся валов. М.: Изд-во АН СССР, 1959. 248 с. (Dimentberg F. Flexural Vibrations of Rotating Shafts. Butterworth, London, 1961.)

2. Genta G. Dynamics of Rotating Systems. Springer-Verlag, New York, 2005. 658 p.

3. Тимошенко С. П. Колебания в инженерном деле. М.: Физматгиз, 1959. 440 с. (Timo-shenko S. P. Vibration Problems in Engineering. Van Nostrand, Toronto, 1928.)

4. Пасынкова И. А. Бифуркации прецессионного движения неуравновешенного ротора // Прикладная математика и механика. Т. 70. Вып. 4. 2006. С. 605—616. (Pasynkova I.A. Bifurcations of the precessional motion of an unbalanced rotor // Journal of Applied Mathematics and Mechanics, 2006. Vol. 70. Issue 4. P. 549-559.)

5. Гробов В. А. Асимптотические методы расчета изгибных колебаний валов турбомашин. М.: Изд-во АН СССР. 1961. 166 с.

6. Genta G. Vibration of structure and machines: practical aspects. Springer-Verlag, Berlin, 1999. 599 p.

7. Pasynkova I. A. Whirling motion of an Unbalanced Rotor in Linear and Nonlinear Elastic Bearings // Den 7. Magdeburger Maschinenbau-Tage. 11-12. Oktober 2005. Tagungsband, 2005, 143-148.

Статья поступила в редакцию 26 декабря 2013 г.

Сведения об авторах

Пасынкова Инна Анатольевна — доктор физико-математических наук, профессор; [email protected]

FORCED OSCILLATION OF AN UNBALANCED ROTOR IN NONISOTROPIC BEARINGS

Inna A. Pasynkova

St.Petersburg State University, Universitetskaya nab., 7/9, St.Petersburg, 199034, Russian Federation; [email protected]

An influence of bearings anisotropy on forced oscillation of an unbalanced rotor with four degrees of freedom was considered. Two cases of the rotor mounting on the rigid shaft and on the flexible one were thoroughly examinated. The differential equations of the rotor rotation were obtained in complex variables. The exact solution of the system was found. Appearance of the additional critical frequencies especially in the bottom part of the spectrum was shown. The dynamic responses and the elliptical trajectories of the shaft's ends were obtained It was also shown that the direct and reverse whirling motion could be realized. Refs 7. Figs 6.

Keywords: rotor dynamics, direct whirling, reverse whirling, anisotropy of bearings.

ХРОНИКА

24 апреля 2013 г. на заседании секции теоретической механики им. проф. Н. Н. Поляхова в Санкт-Петербургском Доме Ученых РАН выступил кандидат физ.-мат. наук, доцент А. Б. Бячков (Пермский государственный национальный исследовательский университет) с докладом на тему «Дифференциально-алгебраические модели динамики несвободного движения».

Краткое содержание доклада:

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

i Надоели баннеры? Вы всегда можете отключить рекламу.