Engineering method of calculation of gasdynamic foil bearings of turbocompressors

Abstract

In this article the engineering method of calculation of gasdynamic foil bearings is considered. Calculations according to this method were made and comparison with the experimental results, obtained by other authors, was done. Ranges of application are also considered.

Full Text

Постановка задачи Возрастание скорости и энергонапряженности современных турбохолодильников для систем кондиционирования воздуха (СКВ) летательных аппаратов (ЛА) с одновременным ужесточением требований к ним по надежности, габаритным размерам и массе вызывает необходимость совершенствования подшипниковых узлов и, в первую очередь, за счет совершенствования форм поверхности качения. В этой связи оправдан неуклонный рост интереса к разработкам и внедрению газодинамических подшипников в быстроходных лопаточных машинах различного назначения в РФ и за рубежом. Это, в первую очередь, касается турбохолодильников для СКВ гражданских и военных самолетов, микрогазотурбинных установок, малошумных вентиляторов, компрессоров различного назначения, в т.ч. для тепловых насосов и т.д. Этот интерес вызван следующими принципиальными особенностями газовых подшипников: · низкая вязкость газа обеспечивает низкие потери на трение, увеличение срока службы и повышение энергоэффективности; · инертность газа по отношению к перекачиваемой среде; · температурная стабильность свойств газа обеспечивает работоспособность подшипника при высокой температуре, которая зависит только от используемых материалов компонентов. Газовые подшипники по принципу работы разделяются на: o газостатические – с внешним наддувом; o газодинамические – без внешнего наддува; o гибридные. В статье рассматриваются лепестковые газодинамические подшипники радиальные и упорные без внешнего наддува применительно к работе в СКВ гражданских самолетов и перспективы их использования в других областях энергетики. В 2010-2011г. сотрудниками МГТУ «МАМИ» при участии сотрудников КБ ОАО НПО «Наука» был выполнен ряд работ по разработке инженерной методики расчета несущей способности газодинамических подшипников в отсутствии внешних нагрузок и оценке виброустойчивости ротора при внешних нагрузках. Инженерный метод расчета лепестковых газодинамических подшипников Основные допущения теории газовой смазки: 1) динамическая вязкость не зависит от давления; 2) течение в газовом слое ламинарное; 3) течение изотермическое; 4) инерционные и гравитационные силы пренебрежимо малы по сравнению с силами вязкости. 5) результатами расчета после выбора конструктивных параметров подшипника (длины, диаметра цапфы и количества лепестков в подшипнике) является определение зависимости его несущей способности и жесткости от режимов работы. На основе полученных данных выбирается жесткость пакета лепестков, исходя из предположения, что эта жесткость больше жесткости газового слоя в подшипнике. Эпюры полного давления в газодинамическом слое , несущая способность и жесткость определяются как функции безразмерных параметров , и . , [Па] (1) , [Н] (2) , [Н/м] (3) где: – суммарное давление в газовом слое, , , – несущие способности, – жесткость. Используемые здесь безразмерные параметры, равны: , , , (4) где: – динамическая вязкость, – угловая частота вращения , 1/с, – радиус вала, м, – статическое давление, Па, – зазор между лепестками и валом, м, – длина подшипника, м. Методика расчета была апробирована на примере турбохолодильника с ротором, состоящим из трех лопаточных колес и двух валов между ними. Параметры: · общая масса вращающегося узла – 3,44 кг; · масса, приходящаяся на левый подшипниковый узел 0,88 кг, на правый – 2,56 кг; · диапазон частот вращения вала 1000 – 37000 об/мин; · диаметр цапфы подшипника 50 мм; · длина подшипника выбрана из условия равенства диаметру как наиболее стабильное сочетание параметров; · количество лепестков ; · монтажный зазор в узле на радиус 50 – 100 мкм; · полуамплитуда кривизны лепестка 10 мкм; · эксцентриситет по условию всплытия принимается равным 50 – 10 – 3 = 37 мкм, где 3мкм условно нулевой зазор при всплытии ротора. При этих условиях несущая сила опоры в зависимости от частоты вращения при условии всплытия данной конструкции приведена на рисунке 1. Рисунок 1 – Несущая сила опоры F [H] в зависимости от частоты вращения вала n [об/мин] при условии всплытия Находим на рисунке 1 пересечение зависимости и статической нагрузки на опору. Для данного варианта они равны 4000 об/мин для первой опоры и 13000 об/мин для второй опоры. Если нас устраивает частота всплытия, то переходим к определению жесткости опоры и рабочих зазоров. Если частота всплытия выше желаемой, то следует изменить диаметр или длину опор. При изменении монтажного зазора стоит учитывать условия сборки и изготовления деталей опор. Минимальный монтажный зазор не рекомендуется задавать меньше 20 мкм. Увеличение зазора свыше 200 мкм может приводить к нежелательным вибрационным нагрузкам, пониженному демпфированию опоры и снижению ресурса опор. Для определения жесткости газового слоя в опоре необходимо, задав частоту вращения на интересующем режиме, изменяя эксцентриситет, получить зависимость несущей способности от него. Жесткость газового слоя опоры определяется отношением разности несущих способностей к разности эксцентриситетов. При увеличении частоты вращения жесткость газового слоя увеличивается. Определим жесткость опор на режиме 35000 об/мин и получим зависимость несущей способности от эксцентриситета. Таблица 1 Несущая сила опора в зависимости от эксцентриситета Как видно из таблицы 1 и рисунка 2, жесткость увеличивается с увеличением эксцентриситета более чем в 2 раза. Также видно, что на максимальном режиме обеспечивается более чем 2-х кратный запас по несущей способности опор. Рисунок 2 – Жесткость опоры в зависимости от эксцентриситета В дальнейшем, используя значения жесткости, можно провести уточненный расчет критических частот роторной системы с учетом упругости опор. В силу определенных ограничений, данная расчетная методика является поверочной и имеет ряд допущений. Помимо допущений об однородности, несжимаемости и изотропности среды в опоре, есть и допущения конструктивного характера. Предполагается, что форма лепестка в расчете носит сформировавшийся характер под воздействием упругих сил и давления смазочного слоя. Это допущение вполне применимо для лепестков без отдельных подкладных пружин. С их внедрением необходимо более точно определять форму лепестков в процессе работы. Одним из расчетных методов, позволяющих учесть этот эффект, является численное моделирование газодинамики опоры с учетом упругости лепестков. Подобные расчеты имеет смысл проводить для отработки новых конструкций, анализе особенностей конструкции из-за гораздо большего времени машинного расчета по сравнению с указанной выше методикой. Разработанная методика инженерного расчета газодинамических опор позволяет, используя основные геометрические и режимные параметры установки, определить несущую способность опор, их жесткость и частоту вращения вала, при которой происходит всплытие вала, критические частоты системы в зависимости от жесткости опор и массы пакетов лепестков. Данная методика апробирована на быстроходных лопаточных машинах с весом ротора от 100 г до 14 кг и рекомендована к применению в текущих проектах ОАО НПО «Наука». Анализ экспериментальных данных лепестковых газодинамических подшипников С целью верификации расчетной модели упорного газодинамического подшипника по результатам зарубежных и отечественных исследований был выполнен обзор современных работ. Сопоставление результатов расчета с экспериментальными данными приведено в таблице 2. Отклонение расчетных значений от экспериментальных составило от 5 % до 20 %. При этом основные характеристики газодинамических подшипников (радиальных и осевых) менялись в следующих диапазонах: · радиальные подшипники газостатические и лепестковые газодинамическиес диаметром вала от 32 до 64 мм при частотах вращения от 50000 об/мин до 90000 об/мин обеспечивали несущую способность от 7 кг до 40 кг; · осевые подшипники – газостатический размером = 130 мм/ 170 мм выдерживал осевую нагрузку 405-920 кг; осевой лепестковый подшипник – = 180 мм / 380 мм – соответственно 10,4 – 30,6 кг. Таблица 2 Данные проверки аналитической методики и сравнения с экспериментальными данными В источнике [5] сравниваются различные методики оценки несущей способности осевых ЛГП в зависимости от частоты вращения (рисунок 3). Рисунок 3 – Зависимость экспериментальной относительной предельной несущей способности осевого ЛГП в зависимости от частоты вращения при различных вариантах вычисления Как следует из рисунка 3, в настоящее время имеются достаточно представительные методы расчета ЛГП. Однако последнее слово за экспериментом, определяющим настоящие ресурс газодинамических подшипников и технологию их изготовления. Расчет влияния конструктивных факторов на виброустойчивость ротора на газодинамических подшипниках Расчеты выполнены для турбохолодильника разработки ОАО НПО «Наука» для режимов, указанных в таблице 3. Таблица 3 Создана аналитическая модель опоры в программном комплексе ANSYS Mechanical, позволяющая проводить расчет критических частот с учетом особенностей данного типа опор с учетом внешних синусоидальных колебаний при жесткости пакета лепестков 200 кН/м и расчет траекторий центра вала для различных частот колебаний и внешних нагрузок (рисунок 4). Рисунок 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 об/с частота генерируемого звука составит = 250 Гц, что находится в области высокой чувствительности человеческого уха. Для радикального снижения шума вентиляторов предлагается существенно увеличить частоту вращения ротора электродвигателя и применить бесконтактные газодинамические подшипники. Так как человеческое ухо слабо слышит звук с частотой выше 2 кГц, то при z = 8 и n = 15000 об/мин = 250 об/с воздушный шум вентилятора с частотой = 2400 Гц превышает 2000 Гц. Особенно эта задача актуальна для космонавтов в замкнутом пространстве. Увеличение частоты вращения позволит также уменьшить радиальные размеры колеса вентилятора, и высокооборотный электродвигатель одной и той же мощности будет иметь существенно меньшие габариты и массу, чем низкооборотный. Недостатком вентиляторов с газодинамическими опорами по сравнению с вентиляторами с опорами качения является большая сложность изготовления высокоточных деталей опор и, соответственно, повышенная стоимость. Однако, учитывая уникальность межпланетных проектов, повышенная стоимость при достижении новых качеств вентиляторов является оправданной. Выводы 1. Проведенные расчетные исследования газодинамических подшипников дают возможность реализации подобной инженерной методики проектирования радиальных и упорных газодинамических подшипников в рамках проектных работ. 2. Для разработки газодинамических опор необходимо создание испытательных стендов для проверки их работоспособности при различных температурах (от 200˚С до 450˚С) с целью расширения областей применения подобных конструкций как в системах СКВ, так и микрогазотурбинных установках, вентиляторах, высокооборотных компрессорных и вентиляционных установках. В этой связи МГТУ «МАМИ» подготовил обзор экспериментальных стендов для низких и высоких температур по материалам, опубликованным в Internet, которые могут быть использованы для проектных работ при создании стенда (стендов) по исследованию лепестковых газодинамических подшипников. Несмотря на позитивную роль представленных численных экспериментов, как отмечено выше, значение прямых экспериментальных работ на испытательных стендах велико для успешного внедрения газодинамических подшипников для конкретных разработок ОАО НПО «Наука». 3. С учетом зарубежного опыта применения газодинамических подшипников не только на пассажирских самолетах, но и на военных самолетах (F-14, F-16) следует начать подобные разработки агрегатов СКВ для боевых самолетов и вертолетов. 4. Перспективной работой по внедрению газодинамических подшипников является микрогазотурбинная установка мощностью 5 кВт для энергетики как альтернатива паротурбинной установке израильской фирмы «Ормат». 5. С учетом того факта, что высокотемпературное антифрикционное покрытие лепестков и вала с малыми коэффициентами трения уже работают в ряде зарубежных фирм, создание отечественного высокотемпературного антифрикционного покрытия с помощью ведущих химических институтов РАН – дело времени.