CC BY
95
К вопросу об инженерной методике расчета лепестковых газодинамических
подшипников турбокомпрессоров
д.т.н. проф. Меркулов В.И., Плыкин М.Е., к.т.н. Тищенко И.В.
Университет машиностроения, ОАО НПО «Наука» 8 (495) 223-05-23, доб. 1297, 1573, tgtd@mami.ru, IV.Tishenko@npo-nauka.ru
Аннотация. В статье рассмотрена инженерная методика расчета лепестковых газодинамических подшипников. Выполнены расчеты по данной методике и выполнено сравнение с экспериментальными результатами, полученными другими авторами. Рассмотрены области применения таких подшипников.
Ключевые слова: лепестковый газодинамический подшипник, турбомашина, микрогазотурбинная установка
Постановка задачи
Возрастание скорости и энергонапряженности современных турбохолодильников для систем кондиционирования воздуха (СКВ) летательных аппаратов (ЛА) с одновременным ужесточением требований к ним по надежности, габаритным размерам и массе вызывает необходимость совершенствования подшипниковых узлов и в первую очередь за счет совершенствования форм поверхности качения.
В этой связи оправдан неуклонный рост интереса к разработкам и внедрению газодинамических подшипников в быстроходных лопаточных машинах различного назначения в РФ и за рубежом.
Это в первую очередь касается турбохолодильников для СКВ гражданских и военных самолетов, микрогазотурбинных установок, малошумных вентиляторов, компрессоров различного назначения, в т.ч. для тепловых насосов и т.д.
Этот интерес вызван принципиальными особенностями газовых подшипников:
• низкой вязкостью газа, обеспечивающей низкие потери на трение, увеличение срока службы и энергоэффективности;
• инертностью газа по отношению к перекачиваемой среде;
• температурной стабильностью газа (работоспособность подшипника при высокой температуре зависит только от материалов компонентов).
Газовые подшипники по принципу работы разделяются на:
• газостатические - с внешним наддувом;
• газодинамические - без внешнего наддува;
• гибридные.
В статье рассматриваются лепестковые газодинамические подшипники радиальные и упорные без внешнего наддува применительно к работе в СКВ гражданских самолетов и перспективы их использования в других областях энергетики.
В 2010-2011 г. сотрудниками МГТУ «МАМИ» при участии сотрудников КБ ОАО НПО «Наука» был выполнен ряд работ по разработке инженерной методики расчета несущей способности газодинамических подшипников в отсутствии внешних нагрузок и оценке виброустойчивости ротора при внешних нагрузках.
Инженерный метод расчета лепестковых газодинамических подшипников Основные допущения теории газовой смазки:
1) динамическая вязкость не зависит от давления;
2) течение в газовом слое ламинарное;
3) течение изотермическое;
4) инерционные и гравитационные силы пренебрежимо малы по сравнению с силами вязкости.
5) результатами расчета после выбора конструктивных параметров (длина, диаметр цапфы и количество лепестков в подшипнике) является определение зависимости несущей способности и жесткости от режимов работы подшипника.
На основе полученных данных выбирается жесткость пакета лепестков, исходя из
Серия 1. Наземные транспортные средства, энергетические установки и двигатели. предположения, что эта жесткость больше жесткости газового слоя в подшипнике.
Эпюры полного давления в газодинамическом слое , несущая способность Wrf и
жесткость Grf определяется как функция безразмерных параметров X, ^ и £ .
f ^ ст) = Pf 1 s, ст) + Prf 2 (x, S, (1)
Wf (х, в, а) = 4W2 + Wn2, (2)
= Wrf / ^, (3)
где: Рг^ - суммарное давление в газовом слое, [Па] Wf, W, Wn - несущие способности, [Н] Grf - жесткости, [Н/м]. Безразмерные параметры
6 u-со- r2 2r L Lf
X =-2—, a = —, s = _ , s f = , (4)
p0 • h L h h
где: ц - динамическая вязкость,
CO - угловая частота вращения, СО = 2%f, 1/с,
r - радиус вала, м,
p0 - статическое давление, Па,
h - зазор между лепестками и валом, м,
L - длина подшипника, м,
e - эксцентриситет подшипника, м,
ef - амплитуда изменения эксцентриситета, м.
Методика расчета была апробирована на примере турбохолодильника с ротором, состоящим из трех лопаточных колес и двух валов между ними. Параметры:
• общая масса вращающегося узла - 3,44 кг;
• масса, приходящаяся на левый подшипниковый узел 0,88 кг; на правый - 2,56 кг;
• диапазон частот вращения вала 1000 - 37000 об/мин;
• диаметр цапфы подшипника 50 мм;
• длина подшипника выбрана из условия равенства диаметру, как наиболее стабильное сочетание параметров;
„2 / 3
• - количество лепестков z ~ r ;
• - монтажный зазор в узле на радиус 50 - 100 мкм;
• - полуамплитуда кривизна лепестка 10 мкм;
• - эксцентриситет по условию всплытия принимается равным 50 - 10 - 3 = 37 мкм, где 3 мкм условно нулевой зазор при всплытии ротора.
При этих условиях несущая сила опоры в зависимости от частоты вращения при условии всплытия данной конструкции приведена на рисунке 1.
Находим на рисунке 1 пересечение зависимости и статической нагрузки на опору. Для данного варианта они равны 4000 об/мин для первой опоры и 13000 об/мин для второй опоры. Если нас устраивает частота всплытия, то переходим к определению жесткости опоры и рабочих зазоров. Если частота выше желаемой, то требуется изменить диаметры или длину опор. При изменении монтажного зазора стоит учитывать условия сборки и изготовления деталей опор. Минимальный монтажный зазор не рекомендуется задавать меньше 20 мкм. Увеличение зазора свыше 200 мкм может приводить к нежелательным вибрационным нагрузкам, пониженному демпфированию опоры и снижению ресурса опор.
О 2000 4000 6000 8000 10000 12000 14000 16000
Рисунок 1. Несущая сила опоры в зависимости от частоты вращения вала при условии
всплытия (Е, Н; п, об/мин)
Для определения жесткости газового слоя в опоре необходимо, задав частоту вращения на интересующем режиме, изменяя эксцентриситет, получить зависимость несущей способности от него. Жесткость газового слоя опоры определяется отношением разности несущих способностей к разности эксцентриситетов. При увеличении частоты вращения жесткость газового слоя увеличивается. Определим жесткость опор на режиме 35000 об/мин и получим зависимость несущей способности от эксцентриситета.
Таблица 1.
Несущая сила опора в зависимости от эксцентриситета
2. МКМ Жесткость. Шмкм
5 16 1.6
10 25 1,8
15 35 2 ш
20 46 12
25 59 2.6
30 74 3
35 96 4.4
Как видно из таблицы 1 и рисунка 2, жесткость увеличивается с увеличением эксцентриситета более чем в 2 раза. Также видно, что на максимальном режиме обеспечивается более чем 2-х кратный запас по несущей способности опор, же сте ость, Ы/ткт
0 5 10 16 20 25 30 35 40
есс, ткт
Рисунок 2. Жесткость опоры в зависимости от эксцентриситета
В дальнейшем, используя значения жесткости, можно провести уточненный расчет
критических частот роторной системы с учетом упругости опор.
В силу определенных ограничений данная расчетная методика является поверочной и имеет ряд допущений. Помимо допущений об однородности, несжимаемости и изотропности среды в опоре, есть и допущения конструктивного характера. Предполагается, что форма лепестка в расчете носит сформировавшийся характер под воздействием упругих сил и давления смазочного слоя. Это допущение вполне применимо для лепестков без отдельных подкладных пружин. С их внедрением необходимо более точно определять форму лепестков в процессе работы. Одним из расчетных методов, позволяющим учесть этот эффект, является численное моделирование газодинамики опоры с учетом упругости лепестков. Подобные расчеты имеет смысл проводить для отработки новых конструкций, при анализе особенностей конструкции, из-за гораздо большего машинного времени расчета, по сравнению с указанной выше методикой.
Разработанная методика инженерного расчета газодинамических опор позволяет, используя основные геометрические и режимные параметры установки, определить несущую способность опор, их жесткость и частоту вращения вала, при которой происходит всплытие вала, критические частоты системы в зависимости от жесткости опор и массы пакетов лепестков.
Данная методика апробирована на быстроходных лопаточных машинах с весом ротора от 100 г до 14 кг и рекомендована к применению в текущих проектах ОАО НПО «Наука».
Анализ экспериментальных данных лепестковых газодинамических подшипников
С целью верификации расчетной модели упорного газодинамического подшипника по результатам зарубежных и отечественных исследований был выполнен обзор современных работ. Сопоставление результатов расчета с экспериментальными данными приведено в таблице 2. Отклонение расчетных значений от экспериментальных составило от 5 % до 20 %.
Таблица 2.
Данные проверки аналитической методики и сравнения с экспериментальными
данными
Геометрия подшипника и режимные условия Аналитические результаты Численные и экспериментальные данные
Радиальный подшипник: r=2.4cm, L=S.4cm, h=16.9-10-'m, Ро=705Рз Щ,0, = 78.7Ы (£=0.335, у=10500грт( КУ1о1 = 66.7Ы (£=0.5, у=4700грт) W1M = 84.3N (experiment} Wl(tl = 74,5N (experiment)
Радиальный подшипник: г*2.4ст, L=5.4cm. h=e.25 1Cr6m. P„=1-10sPa \Ы1а = 125.1Ы (£=0.335, 7=4700грт| №,„, = 244.ЗМ (£=0.5, «=7400грт| W,o1= 115.6N (experiment) W,ct = 236.2N (experiment)
Радиальный статический подшипник: г=1.6ст, 1=4.8ст 1Л=1.5), Л,=0.75, d=0,1mm, с=0.2, h=5Wm, v=50000rpm. Po=5-10'JPa, па=12(2 rows х 6 outlets) = 157М - стат. компонент М« = 221N - динам, компонент Л1М = 37Ш W„ = 184N
Радиальный статический подшипник: г=1,6ст. !=2.4ст (А=0.8), А,=0, d=0.1mm с=0.2, h=5-10-'m, v=50000rpm, Pc=5105Pa. п„=6 (1 raw* 6 outlets) \М„ = 106N - стат. компонент = - динам, компонент \Л/Ш =202М W„ = 58N
Радиальный статическии подшипник: r=3.2cm. 1=9.6ст (А=1.5), А,=0.75, d=0.2rvm, t=0.1, h=10-10tm„ v=50000rpm, Po=510sPa. па=12 (2 rows x 6 outlets) 'Л., = 472Ы - стат. компонент \Л/И = 423И - динам, компонент =895М W,s = 324N
Радиальный лепестковый подшипник: r=1.6cm, l=4.8cm (A=1.5). h=10-1t>im, frO S, nr4 \Л/Н = 105М (Р0=4-Ю5Ра. £=0.1. у=50000грт! - 228N {Р0=4 10*Ра, £=0.2, *=50000грт1 = 403М 1Р0=4-105Ра, Е=0.3, \г=50000грт| \Л/„ = 343М {Р„=6-10'Ра, £=0.2, 1г=75000грт) 375И (Р„=6-Ю'Ра, £=0.2, *=90000грт) W,, = 96N W„ = 194N Wrf = 416N W,< = 350N W„ = 380N
Осевой статический подшипник: г,„=13ст, r„„,=J7cm IВ=2ст. 1*11.8cm), h*12W'm. л„=4. P0=5 105Pa Wr, = 4D58N (Г-0.4тт, д=0.4тт) \Л1„ = 9207Ы (1=1.4тт, Ь=1.4тт) W,,= 5816N Wrf = 9464N
Осевой лепестковый подшипник: г.п-1.8ст, г„„,=3.8ст (г„=2.дст. L0=2cm). h=22.5 10-'m. с^0.778. Р„-110sРа Wr, = 94И (п„=6, у=30000грт) 1Л/Н = 183№ <п^=6, *=5ООООгрт) ччп = 7эм (п^=а. v=зoooorpm) W,= 216N = 306N W„= 104N
При этом основные характеристики газодинамических подшипников (радиальных и осевых) менялись в следующих диапазонах:
• радиальные подшипники газостатические и лепестковые газодинамические с диаметром вала от 32 до 64 мм при частотах вращения от 50000 об/мин до 90000 об/мин обеспечивали несущую способность от 7 кг до 40 кг;
• осевые газостатические подшипники размером йвн / й = 130 мм/170 мм выдерживали
осевую нагрузку 405-920 кг; • осевой лепестковый подшипник размером / = 180 мм /380 мм соответственно
выдерживал осевую нагрузку 10,4 - 30,6 кг. В источнике [5] сравниваются различные методики оценки несущей способности осевых ЛГП в зависимости от частоты вращения (рисунок 3). 1
?о ео
п х (000 Т/мин
Рисунок 3. Зависимость экспериментальной относительной предельной несущей способности осевого ЛГП в зависимости от частоты вращения при различных
вариантах вычисления
Как следует из рисунка 3, в настоящее время имеются достаточно представительные методы расчета ЛГП. Однако последнее «слово» за экспериментом, определяющим настоящие ресурс газодинамических подшипников и технологию их изготовления.
Расчет влияния конструктивных факторов на виброустойчивость ротора на
газодинамических подшипниках Расчеты выполнены для турбохолодильника разработки ОАО НПО «Наука» для режимов, указанных в таблице 3.
Таблица 3
Создана аналитическая модель опоры в программном комплексе ANSYS Mechanical, позволяющая проводить расчет критических частот с учетом особенностей данного типа опор с учетом внешних синусоидальных колебаний при жесткости пакета лепестков 200 кН/м и расчет траекторий центра вала для различных частот колебаний и внешних нагрузок (рисунок 4).
-3.2 -1,6 О 1.6 Э.2
-.8 .8 z.i 4
TIME
Рисунок 4. Траектория центра вала для частоты 100 Гц и нагрузки 0,5g
Современные применения газодинамических опор в конструкциях высокооборотных
лопаточных машин
1. Авиационные СКВ
В США в течение 25 последних лет велись и ведутся работы по реализации газодинамических подшипников в агрегатах самолетов различного назначения.
В источнике [1] изложен обзор расчетных и экспериментальных работ по отработке лепестковых газодинамических подшипников турбомашин в СКВ на самолетах в течение последних 25 лет.
Из этого материала следует, что разработаны конструкции ЛГП для широкого спектра гражданских и военных самолетов.
Ведущая фирма США - разработчик и изготовитель подшипников и высокоскоростных турбомашин для СКВ - Hamilton Sundstrand реализовала свою продукцию для самолетов Boeing 777. Подшипники выдерживают 36000 пусков и остановок, что соответствует времени эксплуатации 30 лет.
2. Микрогазотурбинные установки
В микрогазотурбинных установках мощностью 100 кВт и менее применяются газодинамические подшипники.
Примерами являются разработки фирмы Capstone (США) и ЦИАМ им. Баранова (РФ). Если вопрос высокотемпературных антифрикционных покрытий в США решен, то в РФ находится в стадии разработок. Необходимая рабочая температура покрытия 400 - 500°С. В работе [3] приведены данные о покрытиях с коэффициентом трения 0,1, работоспособные при 500 пусках-остановах и максимальной температуре 810°С с нагрузкой 9,1 Н. Испытания покрытий были проведены от 370°С до 810°С.
3. Турбогенераторные установки и бытовые вентиляторы.
В работе [2] изложена практика внедрения ЛГП в промышленных турбовентиляторах.
Бытовые вентиляторы обладают повышенным шумом, увеличивающимся в процессе износа шарикоподшипников и вследствие срывающихся с лопастей вихревых потоков. Интенсивность шума зависит от частоты вращения ротора.
Так при числе лопастей z = 5 и частоте вращения ротора n = 3000 об/мин = 50 об/с, частота генерируемого звука составит f = z ■ n = 5 • 50 = 250 Гц, что находится в области высокой чувствительности человеческого уха.
Для радикального снижения шума вентиляторов предлагается существенно увеличить частоту вращения ротора электродвигателя и применить бесконтактные газодинамические подшипники. Так как человеческое ухо слабо слышит звук с частотой выше 2 кГц, то при z = S и n = 15GGG об/мин воздушный шум вентилятора с частотой f = z • n = 25G • S = 24GG Гц
> 2GGG Гц. Особенно эта задача актуальна для космонавтов в замкнутом пространстве.
Увеличение частоты вращения позволит уменьшить радиальные размеры колеса вентилятора, высокооборотный электродвигатель при одной и той же мощности имеет существенно меньшие габариты и массу, чем низкооборотный.
Недостатком вентиляторов с газодинамическими опорами по сравнению с вентиляторами с опорами качения является большая сложность изготовления высокоточных деталей опор и, соответственно, повышенная стоимость. Однако, учитывая уникальность межпланетных проектов, повышенная стоимость при достижении новых качеств вентиляторов является оправданной.
Выводы
1. Проведенные расчетные исследования газодинамических подшипников дают возможность реализации подобной инженерной методики проектирования радиальных и упорных газодинамических подшипников в рамках проектных работ.
2. Для разработки газодинамических опор необходимо создание испытательных стендов для проверки их работоспособности на различные температуры от 2GG0C до 45G0C с целью расширения областей применения подобных конструкций как в системах СКВ, так и микрогазотурбинных установках, вентиляторах, высокооборотных компрессорных и вентиляционных установках.
В этой связи МГТУ «МАМИ» подготовил обзор экспериментальных стендов для низких и высоких температур по материалам, опубликованным в Internet, которые могут быть использованы для проектных работ при создании стенда (стендов) по исследованию лепестковых газодинамических подшипников.
Несмотря на позитивную роль представленных численных экспериментов, как отмечено выше, значение прямых экспериментальных работ на испытательных стендах велико для успешного внедрения газодинамических подшипников для конкретных разработок ОАО НПО «Наука».
3. С учетом зарубежного опыта применения газодинамических подшипников не только на пассажирских самолетах, но и на военных самолетах (F-14, F-16) следует начать подобные разработки агрегатов СКВ для боевых самолетов и вертолетов.
4. Перспективной работой по внедрению газодинамических подшипников является микрогазотурбинная установка мощностью 5 кВт для энергетики как альтернатива паротурбинной установки израильской фирмы «Ормат».
5. С учетом того факта, что высокотемпературное антифрикционное покрытие лепестков и вала с малыми коэффициентами трения уже работают в ряде зарубежных фирм, создание отечественного высокотемпературного антифрикционного покрытия дело времени с помощью ведущих химических институтов РАН.
Литература
1. Agrawal. Foil air/gas bearing technology. An overview. ASME, 345E. New York, lGGl, publication 9l-GT-34l
2. GT-2GG3-3S6G6. Development of high-speed industrial turbo blowers with foil air bearing. ASME, Turbo Expo, 2GG3, USA
3. Low-friction wear-resistant coating for high-temperature foil bearings, report N.Y, 122G5, USA
4. Технические справки и отчеты МГТУ «МАМИ» за 2GlG-2Gllrr.
5. Ермилов Ю.И. Осевые лепестковые газодинамические подшипники, диссертация на соискание ученой степени кандидата технических наук, 2GGSr.
