Файл: Гусаров А.А. Балансировка гибких роторов с распределенной массой.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 07.07.2024

Просмотров: 149

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

реакции опор, снизить изгибающие моменты в роторе и обеспе­ чить сбалансированность последнего в некотором диапазоне ■скоростей.

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

Решение дифференциальных уравнений колебаний неуравновешенного гибкого ротора

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

Важным обстоятельством, существенно облегчающим реше­ ние задачи балапсировки гибких роторов, является свойство ортогональности собственных форм колебаний, отмеченное еще А. Мельдалем [6]. В математической форме это условие запи­

сывается

так:

 

ъ

п

если г = п

о

 

если г =j= п

 

 

где Z,. (s) — упругая линия ротора при критической скорости сог; Zn (s) — то же, при критической скорости со„.

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

Отсюда следует, что любую распределенную по длине неуравно­ вешенность можно разложить на составляющие, соответствую-' щие формам упругой линии при собственных колебаниях. Каждая из этих составляющих вследствие условия ортогональности вызы­ вает вынужденные колебания только по своей форме. Баланси­ ровку гибкого ротора поэтому можно проводить раздельно по каждой из составляющих последовательно на каждой критиче­ ской скорости, где данная составляющая имеет преобладающее. значение п изгиб ротора происходит в основном по соответствую­ щей форме упругой линии. Балансировка, выполненная таким образом с помощью распределенных по длине масс, должна при­ вести к уравновешенности ротора на всех скоростях. Теорети­

9



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

Рассмотрим общую методику решения уравнений колебаний двухопорного ротора постоянного сечения (рис. 1.3) с равномерно распределенной массой и с непрерывно распределенной по длиненеуравновешенностью. При этом будем предполагать, что трение отсутствует, так как в работе [7] было показано, что в реальных машинах при отсутствии специальных демпферов на скоростях, не близких к критическим, практически можно не учитывать влияние трения.

Р п с. 1.3. Гибкий ротор с распределенным дисбалансом

Дифференциальное уравнение колебаний ротора в комплекс­ ной форме может быть получено на основании известного соотно­ шения ElzІѴ = q.

При колебательном

движении нагрузка q представляет силы

инерции элементарных

масс ротора

 

q — Tfiz -j- mco2p (s) exp i [coi + VF (s)],

 

и уравнение колебапий будет иметь вид

 

EIzlv/m + z = (о2р (s) exp i [coi + XF (s)].

(1.4)

Функцию распределения неуравновешенных масс можно раз­ ложить в ряды по формам колебаний ротора, которые для случая жестких шарнирных опор выражаются функциями синуса раз­ личной кратности

СО

 

р (s) exp i*F (s) = 2 (An + iBn) sm{nnsß).

(1.5)

П=1

 

Фаза каждой гармоники определяется фазой соответствую­ щего комплексного множителя.

С учетом разложения (1.5) дифференциальное уравнение коле­

баний ротора

будет

 

°с

Elzivjm +

z = со2ехр mt 2 {Ап + іВп) sin {mts/l).

 

п=1

10


Это уравнение решается. с помощью подстановки

оо

z = ехр Ш 2

Znsin (nns/l),

 

(1.6)

 

 

«=i

 

 

 

приводящей

его

к виду

 

 

ОО

 

 

оо

 

 

2 к

-

со2) Zn = со2 2 (4п +

в п),

(1.7)

п = 1

 

 

п= 1

 

 

где

 

 

 

 

 

соп =

(пя/1)2 'YEI/т.

 

 

Из выражения (1.7) находим значения Zn

Zn = г« (Ап +

іВпЖ1 — T n )

(я =

1, 2, 3 ,...) ,

подставляя которые в выражение (1.6) получаем

ОО^2

z = ехр Ш 2

TZZT {Ап +

iBn)sin (nns/l).

(1.8)

 

 

П=1

‘п

 

 

 

 

Решение (1.8) определяет вынужденные колебания рассма­

 

триваемого ротора. Зависящие от начальных условий собствен­

 

ные колебания в задаче балансировки нас не интересуют, так

 

как в реальных машинах они быстро затухают.

 

 

Уравнение (1.4) можно решать и не разлагая функцию распре­

 

деления

неуравновешенности

в ряды

по формам колебаний.

 

В этом случае для освобождения от временной функции приме­

няем подстановку

 

 

 

 

 

z = Z (s) ехр

ісо/,

 

 

(1.9)

'

которая приводит уравнение (1.4) к виду

 

 

ZIV — k:iZ =

/с*р (s) ехр PF (s),

 

(1.10)

где к4 ;= ігмй2ІЕІ.

 

 

 

(1.11)

Уравнение (1.10) является обыкновенным уравнением чет­

вертого порядка с правой частью. Общее решение его может

быть записано в форме [8]

 

 

 

 

Z =

/15 (ks) +

ВТ (ks) +

CU (Äs) +

DV (ks) + Ф (s).

(1.12)

Здесь А, В, С, D — произвольные постоянные, определяемые из ■соответствующих граничных условий; 5 (/cs), Т (ks), U {ks), V {ks) — функции А. Н. Крылова, представляющие линейные ■комбинации тригонометрических и гиперболических функций:

5 (/cs) = х/г (ch /cs + cos ks),

T (/cs)

=

1/a (sh ks -j- sin ks),

U (/cs) — V2 (ch ks — cos ks),

V {ks)

=

V2 (sh ks — sin ks).

 

 

 

(1.13)

11


Функции А. Н. Крылова обладают следующими свойствами:

dS/ds =

kV, dT/ds = kS, dU/ds = kT, dVIds = kU, S (0) = 1,

T (0) =

U (0) = V (0) = 0. В случае, когда аргумент к является

комплексной величиной к = р А іа и а мало, функции А. Іі. Кры­ лова с точностью до величин первого порядка малости можно вы­ разить через функции от действительного параметра в виде

S (ks) Ä ; S (ps) -f msF (pS),

T (ks) «

T (ps)

-f- iasS (ps),

U (ks) Ä ; U (ps) А icisT (ps),

V (ks) Ä

V (ps)

iasU (ps).

Функция Ф (s) в выражении (1.12) представляет частное реше­ ние при нулевых начальных условиях, равное

S

Ф(а)

о

где q (s) — внешняя возмущающая сила, равная в рассматривае­ мом случае

q (s) = k4p(s) exp i ¥ (s).

Подставив в выражение (1.12) значения найденных с учетом гра­ нитных условий произвольных постоянных и умножив получен­ ное выражение на временную функцию ехр (іш£), получим иско­ мое решение (1.9) дифференциального уравнения колебаний ро­ тора.

Рассмотренное уравнение (1.4) и его решение применимы для однопролетного ротора постоянного сечения с непрерывно распре­ деленной нагрузкой при отсутствии в пролете сосредоточенных сил.

Для ротора с изменяющимися по длине величинами попереч­ ных сечений или нагрузки, а также для многоопорных роторов не­ обходимо ротор разбивать на несколько участков. Границами' участков служат сечения, в которых либо меняется размер попе­ речного сечения гп (момента инерции І пи погонной массы тп), либо' расположена опора, либо приложена сосредоточенная сила или изменяется нагрузка. В пределах каждого участка величины поперечного сечения, погонной массы, момента инерции и нагруз­ ки полагаем неизменными. Длину »-го участка обозначаем Іп.

Для каждого участка составляем уравнения типа (1.4), причем начало координат каждый раз переносится на границу участка. Записывая в форме (1.12) решения уравнений (1.4) для каждого участка, получим п уравнений упругой линии ротора по участкам.

Zn = AjiiS (knsn) -f- BnT (knsn) -p CnU (knsn)

D nV (knsn) -j-

+ Фп (sn)>

(1.14)'

где К = mn(a2/EIn.

Произвольные постоянные A n, Bn, Cn, D n определяются из: условий сопряжения на границах участков и на опорах.

12