CC BY
22УДК 621.822.575: 621.892.9 - 032.1 DOI: 10.12737/20238
Л.А. Савин, А.В. Сытин, В.О. Тюрин
СМАЗКА ЛЕПЕСТКОВЫХ ГАЗОДИНАМИЧЕСКИХ ПОДШИПНИКОВ МАЛОВЯЗКИМИ КРИОГЕННЫМИ ЖИДКОСТЯМИ
Рассмотрена возможность применения маловязких криогенных жидкостей в качестве смазочных материалов для лепестковых газодинамических подшипников. Приведена информация о теп-лофизических свойствах криогенных жидкостей. Рассмотрено формирование математической модели в расчете характеристик однолепестковых газо-
динамических подшипников с гофрированными элементами. Рассмотрены вопросы аппроксимации теплофизических свойств криогенных жидкостей.
Ключевые слова: лепестковые газодинамические подшипники, маловязкие криогенные жидкости, гофрированный элемент, теплофизические свойства, математическая модель.
L.A. Savin, A.V. Sytin, V.O. Tyurin LUBRICATION OF TULIP GAS BEARINGS BY LOW-VISCOUS CRYOGENIC LIQUIDS
The efficiency increase at minimum dimensions of movable turbo-sets is achieved at the expense of maximum admissible values of rotation frequencies (up to hundreds of thousands of rpm). At present as promising bearings with an elastic surface are tulip gas bearings. The advantages of these tulip gas bearings are the following: the capacity to self-adjustment, capacity for work under conditions of thermal deformations and at a certain gas pollution fed for lubrication. The computation of tulip gas bearings supposes a joint solution of the problems of gas dynamics, thermal physics, and theory of elasticity. The function of bending included in the equation of gap is calculated through common equations of elasticity theory for corresponding solids
and in the totality with boundary conditions gives a complete description of elastic surfaces in a friction bearing. As lubrication material in gas bearings with elastic elements air is mostly used. In addition, it is possible to use low-viscous cryogenic liquids, for instance, neon. The comparison of thermal physical properties of air and low-viscous cryogenic liquid (neon) witnesses of neon application possibilities, but only at temperature of below 150°K at which its viscosity is equal approximately to the viscosity of air under normal conditions.
Key words: tulip gas bearings, low-viscous cryogenic liquids, corrugated element, thermal and physical properties, simulator.
Введение
Высокоскоростные турбомашины в процессе проектирования предполагают расчет не только проточной части, состоящей из входных и выходных патрубков, улитки и рабочих колес, но также роторно-опорного узла. При этом рост производительности при минимальных размерах нестационарных турбоагрегатов достигается за счет максимально допустимых значений частот вращения (до сотен тысяч оборотов в минуту) [1]. С минимальными потерями на трение и нагревом данным требованиям удовлетворяют только подшипники с газовой смазкой, которые, как нестандартные узлы, в каждом конкретном случае требуют всестороннего анализа, математического моделирования и проверочных расчетов по основным показателям, таким как не-
сущая способность и потери мощности на трение. Главным условием расчета подшипников скольжения является определение полей давления в смазочном слое. Высокие частоты вращения служат источником значительных вибраций в опорах роторов, что может привести к снижению ресурса и разрушению не только роторно-опорного узла, но и рабочих колес, ремонт которых - сложный и дорогой процесс. Решением этой проблемы может быть применение подшипников с упругой поверхностью, среди которых перспективными на данный момент являются лепестковые газодинамические подшипники (ЛГДП). Разнообразие конструкций, материалов и покрытий позволяет подобрать такие подшипники в качестве опор рото-
ров высокоскоростных турбомашин в соответствии с требованиями.
Достоинства лепестковых газодинамических подшипников: способность к самоустановке, работоспособность в условиях тепловых деформаций и при некотором загрязнении газа, поступающего для смазки. Силы трения, возникающие при работе в зонах контакта между лепестками, между лепестками и корпусом, а также силы упругой деформации лепестков рассеивают энергию колебаний ротора и улучшают динамические характеристики машины, позволяют достигать высоких частот вращения. Податливость лепестка позволяет обеспечить чрезвычайно малые зазоры между цапфой и лепестком на большой площади и получить несущую способность, достаточную для многих высокооборотных турбомашин различного назначения [2].
На сегодняшний день запатентован и испытан ряд лепестковых подшипников учеными из США и других стран (пат. Ш № 4415280, 5634723), российскими учеными (пат. РФ № 2350794) и авторами данной статьи (пат. РФ № 2525497, 2558161, 2568005; рис. 1 [2]).
Среди разнообразия твердых смазочных покрытий, разработанных применительно к узлам трения механизмов, в лепестковых подшипниках нашли применение два вида покрытий: с органическим и неорганическим полимерным связующим. Эти покрытия наносятся методами напыления (из пульверизатора), электроосаждения, электростатического напыления. После нанесения и термообработки покрытие в некоторых случаях полируют. В массу связующего вводят наполнители двух типов: антифрикционные (графит, сульфиды, селениды и теллуриды тяжелых металлов) и износостойкие (порошки окиси хрома, стекла, окиси алюминия) [3].
Рис. 1. Лепестковый газодинамический подшипник
Расчет лепестковых газодинамических подшипников
Расчет лепестковых
газодинамических подшипников
предполагает совместное решение задач газодинамики, теплофизики и теории упругости (рис. 2). В зависимости от геометрии входящих в подшипник упругих элементов, которые можно представить в виде пластин, оболочек, пологих оболочек и т.п., функция прогиба, входящая в уравнение зазора, рассчитывается по классическим уравнениям теории упругости для соответствующих тел и в совокупности с граничными условиями дает полное описание упругих поверхностей подшипника скольжения [4].
Расчетная схема лепестка представляет собой тонкую незамкнутую цилиндрическую оболочку, защемленную по одной образующей и свободную по другой, нагруженную внутренним давлением (рис. 3). Торцы оболочки свободны, нагрузка не симметрична относительно оси цилиндра,
следовательно, деформация сводится к прогибу. Величину прогиба в подобных случаях можно получить с достаточной точностью, совершенно пренебрегая растяжением срединной поверхности оболочки.
Рис. 2. Схема решения комплексной задачи: м - прогиб упругого элемента, м; рт - касательные напряжения, Па; к0 - начальный зазор, м; е - эксцентриситет, м; 0 - угловая координата, рад; ф - расчетный угол, рад; для смазочного материала: р - давление, ат; Т- температура, К; р - плотность, кг/м3; Ср - теплоемкость, Дж/(кг-К); I - энтальпия, Дж/кг; р - вязкость, Па-с
Рис. 3. Расчетная схема лепестка
Для описания деформированного оболочечного элемента моментную цилиндрических
теорию оболочек
напряженно-состояния используем расчета [5], что
позволяет рассмотреть неосесимметричное нагружение.
Э2и 1 -уЭ2и 1 + у Э2Ф V Эw _
^т-т +-+-^т^ + —^ = 0;
Э72 2Я2 Э92 2Я Э7Э0 Я Э7
1 + у Э2и 1 Э2Ф 1 -V Э2Ф 1 Эw 1 -V
2Я Э7Э0 + Я2 Э62 + 2 Э72 + Я2 Э6 Е5 Рт(8,7):
V Эи 1 ЭФ w
--+ ^-+ —т +
Я Э7 Я2 Э6 Я2
+к.
Я
Я3 ^ + 2Я^ Э74 Э72Э92
.1ЭЧ
+ Я Э64
1 -V2
Е 5
Р (8,7
где р(0,^) - нормальное контактное усилие (давление); Рт(0,^) - касательное контактное усилие; z - осевая координата, м; 9 - перемещение в окружном направлении, м; и - перемещение в осевом направлении, м; 5 - толщина оболочки, м; V - коэффициент Пуассона; Я -внутренний радиус подшипника, м; Е -модуль упругости, Па; к] = 52 / 12Я2.
Рассмотрение кругового
гофрированного элемента (рис. 4) целесообразно начать с отдельного гофра. Каждый гофр рассматривается в виде незамкнутой цилиндрической оболочки со свободными торцами и свободным
Рис. 4. Расчетная схема кругового гофрированного элемента
опиранием образующих на корпус подшипника.
Деформация упругих элементов находится с учетом совместного прогиба
лепестка и гофра под действием приложенной к ним нагрузки (рис. 5).
Рис. 5. Схема совместного прогиба лепестка и гофра
В данной работе считается, что лепесток и гофр, в недеформированном состоянии представляющие две цилиндрические оболочки, контактируют по общей образующей. Главными условиями совместного прогиба являются условия контакта лепестка и гофра после
деформации и взаимного непроникновения контактирующих поверхностей [2].
В общем виде процесс совместной деформации лепестка и кругового гофрированного элемента под действием газодинамического давления представлен на рис. 6.
Рис. 6. Совместная деформация лепестка и кругового гофрированного элемента
Расчет системы сводится к последовательному расчету элементов с последующим перерасчетом с учетом их взаимного влияния. Этим достигается
максимальное приближение к реальным процессам, происходящим во время работы данной конструкции.
Применение маловязких криогенных жидкостей в качестве смазочного материала
Критериями выбора смазочных материалов для подшипников жидкостного трения могут быть различные факторы в зависимости от назначения машин и агрегатов. Так, для смазки опорных узлов агрегатов двигательных установок летательных аппаратов, например турбонасосных агрегатов ЖРД, используется перекачивание компонентов топлива. Смазка и охлаждение аэростатодинамических подшипников высокоскоростных электрошпинделей и турбокомпрессоров осуществляется
воздухом, подающимся под давлением или из окружающей среды. В качестве смазочного материала в лепестковых газодинамических подшипниках наиболее часто используется воздух [6; 7].
В последние годы находят широкое применение нетрадиционные смазочные материалы. В первую очередь это низкотемпературные и криогенные жидкости (жидкий кислород, водород, неон, аммиак и т. д.), сфера применения которых связана с использованием в энергетических силовых и
вспомогательных установках
транспортных и технологических систем. При смазке и охлаждении подшипников этими маловязкими жидкостями возникают специфические явления в
несущих слоях, обусловленные
повышенной турбулентностью,
критическими течениями и фазовыми превращениями [4].
Применение криогенных рабочих тел в транспортных системах связано с решением двух глобальных проблем: разработкой надежных экономичных двигателей длительного ресурса и созданием промышленных установок для получения дешевых низкотемпературных жидкостей. Одним из важных компонентов их реализации является использование высокоскоростных агрегатов с турбинным приводом (насосы, компрессоры, детандеры), работающих на лепестковых подшипниках при смазке маловязкими криогенными жидкостями, например неоном. Для большинства параметров один из факторов - давление или температура - является определяющим, поэтому большинство параметров зависят от одного главного фактора при определенном значении другого. Результаты аппроксимации свойств неона в диапазоне Т = 28...280 К, Р = 0,1...20 МПа и свойств воздуха в диапазоне Т = 270.380 К, Р = 0,1.30 МПа представлены в виде аналитических зависимостей /=/(Р,Т) и приведены в таблице [8; 9].
Таблица
Параметр Аналитическая зависимость
Неон
Вязкость ц, х10-6 Па-с 6,638 + 17,172Т +10,576Р - 0,035Т2 - 0,108Р2 - 0,061ТР + +4,68 ■ 10-5 Т3 - 0,014Р3 + 0,0029ТР2 - 2,6 ■ 10-5 Т2Р
Плотность р, х103 кг/м3 20,833 - 0,328Т + 0,0018Т2 -3,3210-6Т3 + 4,9761п(Р) + 1,144(1п(Р))2
1 + 0,008Т - 7,12 ■ 10-6Т2 - 4,5 ■ 10-9Т3 - 0,3811п(Р) -1,64 ■ 10-5 (1п(Р))2
Энтальпия I, х105 Дж/кг -1044,48 + 870,651п(Т) -14,319Р - 224,67(1п(Т))2 - 0,002Р2 + 5,626Р 1п(Т) + +20,42(1п(Т))3 -7,58■Ю-5Р3 + 0,8■ 10-3Р21п(Т)-0,552Р(1п(Т))2
Теплоемкость СР, х103 Дж/(кг-К) 2,54-0,6461п(Т) + 0,044Р + 0,068(1п(Т))2 -7,76■Ю-5Р2 -0,008Р 1п(Т)
Воздух
Вязкость ц, х10-6 Па-с 10,85 + 0,221Р - 0,0358Р2 + 0,37 ■ 10-3 Р3 - 0,02321Т - 0,24 10-3 Т2
1 + 0,0184Р - 0,00156Р2 + 2,2937 ■ 10-5 Р3 - 0,00657
Плотность р, х103 кг/м3 0,00718 - 0,0164Р - 2,251 ■Ю-5 Т
1 + 0,00191Р - 0,264 ■ 10-3Р2 - 0,009777 + 6,128 10-6 Т2
Энтальпия I, х105 Дж/кг 14,282 + 0,0276Р 6267,21 0,482 ■Ю"3 Р2 +13351?92,7 Т Т2 - 12,9Р + 3,53110-6РЗ _1,06 3108 _ 615,06Р + °,196Р2 Т у т Т Т
Теплоемкость СР, х103 Дж/(кг-К) 1,113-7,110-3Р+0,374■Ю-3Р2 +4,079Ю-6Р3 -7,4440-3Г+1,70540-5Г2
1+0,0138Р+0,528 10-3Р2 -6,6910-37+1,5810-572
Графические результаты аппроксимации представлены на рис. 7.
Рис. 7. Графические результаты аппроксимации теплофизических свойств воздуха и неона
Заключение
Сравнение теплофизических свойств воздуха и маловязкой криогенной жидкости (неона) свидетельствует о возможности применения неона, но только при температуре ниже 150 градусов Кельвина, при которой его вязкость приблизительно равна вязкости воздуха при нормальных условиях.
Полученные аппроксимационные зависимости неона позволяют рассчитывать лепестковые подшипники в неизотермиче-
ской постановке. Аппроксимационные зависимости с данной степенью точности (погрешности расчета не превышают 1%) при подстановке в уравнение Рейнольдса для определения полей давления и уравнение баланса энергии позволяют рассчитывать подшипники скольжения в неизотермической постановке, что дает более точное описание процессов в смазочном слое, приближая математическую модель к реальному объекту.
Статья подготовлена в рамках выполнения проекта № 9.101.2014/К проектной части государственного задания в сфере научной деятельности «Гидродинамические эффекты в напорно-сдвиговых течениях сред сложной реологии в каналах переменной геометрии» (2014 - 2016 гг.).
СПИСОК ЛИТЕРАТУРЫ
1. Гаевик, Д.Т. Подшипниковые опоры современных машин / Д.Т. Гаевик. - М.: Машиностроение, 1985. - 248 с.
2. Сытин, А.В. Решение комплексной задачи расчета характеристик радиальных лепестковых газодинамических подшипников: дис....канд. техн. наук / А.В. Сытин. - Орел, 2008. - 200 с.
3. Леонов, В.П. Лепестковые газовые подшипники турбомашин / В.П. Леонов, Т.И. Максимович. -М.: Изд-во МГТУ им. Н.Э. Баумана, 2002. - 60 с.
4. Савин, Л.А. Моделирование роторных систем с опорами жидкостного трения: монография / Л. А. Савин, О.В. Соломин. - М.: Машинострое-ние-1, 2006. - 444 с.
5. Хешмет, Х. Анализ газового ленточного радиального подшипника / Х. Хешмет, Дж. А. Уоло-
вит, О. Пинкус // Проблемы трения и смазки. -1983. -Т.105. -№4. - С. 124-132.
6. Пешти, Ю.В. Газовая смазка / Ю.В. Пешти. -М.: МГТУ, 1993. - 381с.
7. Самсонов, А. И. Подшипники с газовой смазкой для турбомашин / А.И. Самсонов. - Владивосток: Изд-во ДВГТУ, 1996. - 112 с.
8. Варгафтик, Н.Б. Справочник по теплофизиче-ским свойствам газов и жидкостей / Н.Б. Варгафтик. - М., 1972. - 720 с.
9. Рабинович, В.А. Теплофизические свойства неона, аргона, криптона и ксенона / В.А.Рабинович, А.А. Вассерман, В.И. Недоступ, Л.С. Векслер. - М.: Изд-во стандартов, 1976. -636 с.
1. Gaevik, D.T. Bearing Support of Modern Machinery / D.T. Gaevik. - M.: Mechanical Engineering, 1985. - pp. 248.
2. Sytin, A.. Solution of complex problem in computation of journal tulip gas bearing properties: Thesis for Can.Eng.Degree / A.V. Sytin. - Orel, 2008. -pp. 200.
3. Leonov, V.P. Tulip Gas Bearings for Turbomachines / V.P. Leonov, T.I. Maximovich. -M.: Publishing House of Bauman STU of Moscow, 2002. -pp. 60.
4. Savin, L.A. Modeling of Rotor Systems with Bearings of Liquid Friction: Monograph / L.A. Savin, O.V. Solomin. - M.: Mechanical Engineering-1, 2006. - pp. 444.
5. Heshmet, H. Analysis of gas foil journal bearing / H. Heshmet, J. A. Wollowit, O. Pinkus // Problems
of Friction and Lubrication. -1983. -T.105. -№4. -pp. 124-132.
6. Peshty, Yu.V. Gas Lubrication / Yu.V. Peshty. -M.: MSTU, 1993. - pp. 381.
7. Samsonov, A.I. Bearings with Gas Lubrication for Turbomachines / A.I. Samsonov. - Vladivostok: Publishing House of FESTU, 1996. - pp. 112.
8. Vargaftik, N.B. Reference Book on Thermal and Physical Properties of Gas and Liquids / N.B. Vargaftik. - M., 1972. - pp. 720.
9. Rabinovich, V.A. Thermal and Physical Properties of Neon, Argon, Krypton and Xenon / V.A. Rabinovich, A.A. Wasserman, V.I. Nedostup, L.S. Weksler. - M., Publishing House of Standards, 1976. - pp. 636.
Сведения об авторах:
Савин Леонид Алексеевич, зав. кафедрой «Ме-хатроника и международный инжиниринг» При-окского государственного университета, e-mail: [email protected].
Сытин Антон Валерьевич, доцент кафедры «Ме-хатроника и международный инжиниринг» При-окского государственного университета, e-mail: [email protected].
Тюрин Валентин Олегович, аспирант Приокско-го государственного университета, e-mail: [email protected].
Статья поступила в редколлегию 5.04.2016. Рецензент: д.т.н., профессор Брянского государственного университета Киричек А.В.
Savin Leonid Alexeyevich, Head of the Dep. "Me-chatronics and International Engineering" Prioksky State University, e-mail: [email protected]. Sytin Anton Valerievich, Assistant Prof. of the Dep. "Mechatronics and International Engineering" Prioksky State University, e-mail: [email protected]. Tyurin Valentin Olegovich, Post graduate student Prioksky State University, e-mail: [email protected].
