Имитационное аналитическое и статистическое моделирование балансировки гибких роторов турбокомпрессоров Текст научной статьи по специальности «Механика и машиностроение»

УДК (621:62-755:62-251)(075.8)

ИМИТАЦИОННОЕ АНАЛИТИЧЕСКОЕ И СТАТИСТИЧЕСКОЕ МОДЕЛИРОВАНИЕ БАЛАНСИРОВКИ ГИБКИХ РОТОРОВ ТУРБОКОМПРЕССОРОВ

© 2010 Н.В. Корнеев

Поволжский государственный университет сервиса, г. Тольятти

Поступила в редакцию 10.03.2010

В статье рассматриваются основные принципы уравновешивания гибких роторов. Приводится аналитическая модель точной низкочастотной балансировки гибких роторов. При использовании метода Монте-Карло строится статистическая модель указанной балансировки на базе программной оболочки МаЛ^эгкБ Ма1ЬаЬ. Проводится анализ и указывается возможность применения полученных результатов в производстве.

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

В машиностроении находят широкое применение так называемые гибкие роторы, для которых выполнение условий уравновешенности уже не обеспечивает спокойной работы при минимальном уровне вибраций [1, 2, 4, 6-7]:

J U (z )dZ +Х Uk = 0;

0 k=1 l m

J U (z Xz - zc )dz +^Uk (zk - zc ) = 0

0 k=1

, (1)

Таким образом, значения и' и и'' отвечают условиям уравновешенности жесткого ротора. Этим, однако, не исключается возможность его прогиба, значение которого в точке крепления диска будет равно[1, 6-7]:

h =

mso

-(a11 -ai2 X

1 -

f

o

o j

где и (z) = m(z )• e (z) - погонный дисбаланс; е (z) - смещения центров масс элементов ротора от оси вращения (эксцентриситеты); m(z) -погонная масса ротора, U¿ - сосредоточенные дисбалансы.

Гибким, согласно ГОСТ 19534-74, называют ротор, который, будучи даже идеально сбалансированным на частотах вращения, значительно меньших первой критической скорости, на более высоких частотах (вплоть до наибольшей эксплуатационной) имеет динамические реакции и прогибы, недопустимые в условиях нормальной эксплуатации.

Сущность проблем, относящихся к этому случаю, покажем на примере одномассового ротора (рис. 1). Диск имеет радиальное смещение (эксцентриситет) s относительно осевой линии вала. Соответствующий этому эксцентриситету дисбаланс U = те скомпенсирован двумя симметрично расположенными корректирующими массами, создающими дисбалансы:

и '=и "=4 2 W.

Корнеев Николай Владимирович, доктор технических наук, доцент кафедры социально-культурного сервиса. E-mail: niccyper@mail.ru

где 0О =

ma

значение критическои ско-

11

рости; ац и а^ - статические коэффициенты податливости.

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

Итак, установка уравновешивающих грузов и' и 0" отнюдь не исключает прогиба вала, который неограниченно возрастает по мере приближения угловой скорости вращения к критической. Значит, выполнения условий (1) недостаточно для уравновешивания гибкого вала. Его уравновешенность будет достигнута только при том дополнительном условии, что ап = а^ . Но

Рис. 1. Схема одномассового ротора

2

1

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

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

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

Основные положения, в соответствии с которыми балансируются гибкие роторы, сводятся к следующему [1, 2, 4, 6-7].

Критерием гибкости ротора служит отношение максимальной рабочей частоты вращения к « « *

первой критической частоте п* ротора, установленного на абсолютно жестких (в радиальных направлениях) опорах. Если это отношение не превышает 0,5, ротор относится к категории жестких, в противном случае - к категории гибких.

Наиболее эффективными являются такие методы балансировки, в которых комбинируются условия статического равновесия и некоторые динамические условия. Таковыми могут быть: устранение определенного числа собственных форм изгибных колебаний из кривой динамического прогиба; устранение прогибов на заданных частотах и в заданных точках; устранение динамических реакций подшипников на заданных частотах и др. С точки зрения эффективности балансировки все они примерно равноценны, а чис-

ло потребных плоскостей коррекции на 2 превышает количество выбранных динамических условий. Для ротора, работающего в диапазоне 0 < п < 4,5п*, чаще всего достаточно четырех плоскостей; в диапазоне 0 < п < 3п* число плоскостей коррекции можно уменьшить до трех.

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

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

В работах [1, 4, 9] предложена технология низкочастотной балансировки гибких роторов компрессорных агрегатов "Синтезгаз", включающая следующие операции. _ _

1. Измерение дисбалансов А0 и В0, центрального вала в плоскостях опорных подшипников;

2. Сборка первого промежуточного узла, включающего вал и один из дисков, и изм ерение в тех же плоскостях дисбалансов А1 и В1;

3. Операция 2 последовательно повторяется для второго узла, включающего два диска, затем для третьего, четвертого, и, наконец, последнего п-го узла, представляющего уже окончательно собранный ротор; _

_4. По результатам измерения дисбалансов Ак и Вк (к=0...п) последовательно определяются:

а) проекции дисбалансов Аку, Аь и Вку, Вь на координатные плоскости УХ и 7Х, жестко связанные с ротором (ось Х направлена по оси вращения);

б) проекции статических и моментных дисбалансов самих дисков:

иу = ААку +АВку; иь = Мь +АВЬ; Уь =ААуХк -В«-Х); К =ААЛ -В(I-Х),

где AAky = \ - Ak-l y;

ABky = Bky- Bk-i,y; aAkz = Akz - Ak-u ; = Bkz - Bk_u; ]

; к - номер

балансируемого узла; хк - расстояние от левой опоры до центра масс к-го диска. Таким образом, функция распределения дисбалансов становится полностью определенной.

в) назначаются плоскости коррекции дисбалансов. Примем, что это плоскости х1, х4 , х7, х8, проходящие через центры масс соответствующих дисков.

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

д) определяются проекции дисбалансов корректирующих масс D1 у, Б4у, Б7у , В8у из системы уравнений:

N+2 r r r

I Dk = -(A + B);

k=1

N+2 r r

I DkXk = -B • l; k=1

N+2 r r

I Dkaik = -rio, k=1

(2)

где N - число динамических условий выполняемых при уравновешивании; A и B - дисбалансы ротора, измеряемые на его опорах; l - расстояние между опорами; a ik - динамические коэффициенты податливости; i=1, 2...N; r^ -прогибы в некоторых точках X = Xi , вызываемые начальным дисбалансом и подлежащие устранению посредством балансировки.

е) определяются проекции дисбалансов корректирующих масс Dlz,D4z,D7z,D8z из системы уравнений, аналогичной системе (2), но записанной для плоскости ZX.

ж) определяются модули и углы дисбалансов корректирующих масс:

i--Diz

Dt =VD2y + D¡ ; щ= arceos D,

где i = 1, 4, 7, 8 - номера плоскостей коррекции;

5. Выполняется коррекция дисбаланса и контроль качества балансировки: при правильно выполненном расчете и точной коррекции величины остаточных дисбалансов, замеренные на опорах ротора, не должны превышать установленные допуски.

Проведенное аналитическое моделирование по системе уравнений (2) с использованием специально разработанного программного обеспечения [2, 3, 8] позволяет сделать вывод о высокой эффективности подобной балансировки, так как с одной стороны, обеспечивается уравнове-

шенность ротора как твердого тела, а с другой -выполняются определенные динамические условия, гарантирующие работу ротора с достаточно низким уровнем вибрации по сравнению с неот-балансированным ротором (рис. 2).

На рис. 2 обозначено: ук = гк/гк0 , где гк , гк0 - динамические прогибы в точках х = хк отбалансированного и неотбалансированного роторов; т/=1...12 обозначают различные сочетания условий, по которым определялись дисбалансы корректирующих масс; кривые Б(3) и Б(4) показывают уровень значений ук для ротора, отбалансированного по предлагаемой технологии соответственно в трех и четырех плоскостях коррекции, а кривая А показывает средне производственное значение ук .

Полученные результаты (рис. 2) позволяют выбрать такой вариант балансировки, при котором обеспечивается минимальное значение у, то есть вариант оптимальный с точки зрения вибронадежности роторной системы.

Однако на выбор метода балансировки оказывают влияние и другие факторы, например, жесткие ограничения на значения корректирующих масс приводят к тому, что достаточно сложно обеспечить сочетаемость минимальных значений коэффициента ук с приемлемыми значениями дисбалансов корректирующих масс Di. Это особенно важно, когда коррекция дисбаланса достигается за счет съема металла, поскольку резерв съёма ограничен.

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

Для реализации статистического поиска нами разработан программный модуль на базе программной оболочки МаЙШо^ Ма^аЬ, реализующий статистическое моделирование указанной технологии балансировки. Алгоритм действия программного модуля представлен на рис. 3.

0.6

ОА

0,2

О

Б[3)

/А

лЛл БА)

1234567 89 10 11 12 v Рис. 2. Результаты расчета у ( v)

И

Ввод сттшюяъного анолитическ?г&условия г=/{аХТ,с)

±

РАСЧЕТ-

Генерация четырех случайных чисел а., а,, а,. а4 е{0,1] I

Оо' Ф-О > = а; 2е' <?7

Ш В71 07сояр7,

. ят ц>4 ; =£)7з гп р7

□

т

Определяем ук, О,, щ с помощью фшямтичесхой модели

111:7

номер колеса

Рис. 4. Значение у при моделировании [1. Т3, '<г-ч

Рис. 3. Алгоритм программного модуля статистического моделирования

Здесь согласно результатам аналитического моделирования, за оптимальный аналитический вариант взят вариант со следующими условиями балансировки: а =4 - количество плоскостей коррекции; г =1, 4, 7, 8; Т =2 - количество плоскостей, в которых устраняются динамические прогибы; с =4, 7 - номера плоскостей, в которых устраняются динамические прогибы. На рис. 3 обозначено: 0О = 2(|Лп| + \Вп|) ; ¿=1...1000 -число исследований, принятое нами для статистического моделирования.

Некоторые результаты расчета, проведенного с использованием разработанного программного модуля, по указанным оптимальным параметрам, аналитической модели, приведены на рис. 4 и 5, где обозначено: 1 - при аналитическом моделировании, 2 - при статистическом моделировании.

ВЫВОДЫ

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

7 8 плоскости коррекции

Рис. 5. Значение Ог при моделировании

го решения, возможно определить не только более оптимальные параметры балансировки, но и добиться более равномерного распределения уровня вибрации по поверхности ротора (рис. 4, кривая 2);

- значения корректирующих дисбалансов при статистическом моделировании не намного превосходят соответствующие значения при аналитическом и отвечают необходимому резерву съёма материала при коррекции;

- отмеченные особенности, свидетельствуют о высокой сходимости аналитической и статистической моделей;

- применение подобного анализа в производстве играет существенную роль, так как аналитическая модель позволяет добиться нулевых значений динамических прогибов в определенных точках, при указанных динамических условиях. Это вызывает необоснованное повышение корректирующих дисбалансов, что недопустимо ввиду ограниченности резерва съема металла при коррекции. Используя статистическое моделирование, возможно наиболее равномерно распределить корректирующие дисбалансы по поверхности ротора;

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

СПИСОК ЛИТЕРАТУРЫ

1. Корнеев Н.В. Методы прогнозирования и снижения вибрации гибких систем турбоагрегатов. Монография. М.: Спутник+, 2007, С. 157.

2. Корнеев Н.В. Многокритериальная оптимизация дисбаланса гибких роторных систем // Известия Самарского научного центра РАН. 2008 Т. 10. №3. С. 830-834.

3. Корнеев Н.В. Аналитическая и статистическая оптимизация уровня дисбаланса гибких систем турбоагрегатов // Машиностроитель. 2007. №12. С. 25-28.

4. Корнеев Н.В. Технология динамической балансировки роторов турбоагрегатов // Известия Самарского научного центра РАН. 2009. Т. 11. №5. Самара, 2009. С. 216-221.

5. Корнеев Н.В. Методология прогнозирования дисбаланса деталей и узлов турбоагрегатов //Машино-

строитель. 2006. №7. С. 19-21.

6. Корнеев Н.В. Алгоритмы прогнозирования эксплуатационного дисбаланса роторов турбоагрегатов // Экономика и производство. 2006. №3. С. 71-75.

7. Корнеев Н.В. Методология прогнозирования начального дисбаланса турбоагрегатов в условиях сборки // Техника машиностроения. 2006. №3. С. 72-75.

8. Корнеев Н.В. Многокритериальная параметрическая оптимизация динамических характеристик роторных систем турбоагрегатов // Наука-производству. 2006. №6. С. 44...46.

9. Korneev N.V. Forecasting of a vibration level of nonrigid rotary tables of compressor units from an out-of-balance with allowance for of series of the random mechanical, gaseous dynamic and operation factors // Programm Abstracts Internationaler Kongress&Fachmesse „Ökologische, Technologische und Rechtliche Aspekte der Lebensversorgung", EURO-ECO Hannover 2009, Europäische wissenschaftliche Gesellschaft e.V. in Hannover, Europäische Akademie der Naturwissenschaften e.V. in Hannover Leibniz Universität Hannover, Universität Bremen, P. 45-47.

IMITATIVE ANALYTICAL AND STATISTICAL SIMULATION OF TRIMMING OF NONRIGID ROTARY TABLES OF TURBOCOMPRESSORS

© 2010 N.V. Korneev

Volga State University of Service, Togliatti

In the article the philosophys of a counterbalance of nonrigid rotary tables are esteemed. The analytical model of precise low frequency trimming of nonrigid rotary tables is resulted. At usage of Monte Carlo method the statistical model of the indicated trimming on the basis of programmatic shell MathWorks MatLab is created. The analysis is carried conducted and the feasibility of the obtained result ins effecting is underlined.

Keywords: simulation, out-of-balance, trimming, turbocompressor, rotary table

Nikolay Korneev, Doctor of Technics, Associate Professor at the Socio-Cultural Service Department. E-mail: niccyper@mail.ru