УДК 532.525.2 : 532.517:532 : 59
К.С. АХВЕРДИЕВ, М.А. МУКУТАДЗЕ, В.А. КОНСТАНТИНОВ
ГАЗОВЫЙ УПОРНЫЙ ПОДШИПНИК ПОВЫШЕННОЙ НЕСУЩЕЙ СПОСОБНОСТИ С ПОДАТЛИВОЙ ОПОРНОЙ ПОВЕРХНОСТЬЮ
Теоретически исследованы основные рабочие характеристики деформируемых газодинамических упорных подшипников скольжения, обладающих повышенной несущей способностью. Показано, что характеристики стационарной работы такого подшипника отличаются от характеристик жесткого подшипника.
Ключевые слова: упорный подшипник, несущая способность, адаптированный профиль, газовая смазка.
Введение. Рост скорости и рабочей температуры в машинах вращательного движения делает почти неизбежным использование подшипников с газовой смазкой. Существенным недостатком самогенерирующихся газовых подшипников является их низкая несущая способность. Каждый раз, когда нагрузка оказывается высокой или даже средней, требуется внешний наддув, усложняющий конструкцию. Так как размеры подшипника обычно ограничены, то возникает необходимость исследования конструкции подшипников, которые выдерживают максимальную нагрузку при заданных размерах подшипника.
В области гидродинамических подшипников появилось новое направление - подшипники с нежесткой опорной поверхностью. Подшипники с нежесткой поверхностью имеют явные преимущества, по сравнению с подшипниками с нежесткой опорной поверхностью. Податливость эластомера под действием давления смазки приводит к появлению своего рода губы, ограничивающей утечку сдавливаемой смазки из подшипника [1]. Таким образом, в условиях уменьшенного питания смазкой подшипник с податливой поверхностью может сохранять гидродинамическую или гидростатическую сплошную пленку смазки.
Цель наших исследований - расчет упорного газодинамического подшипника повышенной несущей способности с учетом деформации упругого слоя переменной толщины на рабочей поверхности ползуна.
Постановка задачи. Рассмотрим установившееся движение газовой смазки в зазоре упорного подшипника (между ползуном и направляющей). Ползун с нелинейным профилем опорной поверхности предполагается неподвижным, а направляющая движется со скоростью и*. Также предполагается наличие упругого слоя на рабочей поверхности ползуна (см. рисунок).
Схематическое изображение упорного подшипника с нелинейным податливым контуром опорной поверхности:
1 - упругий слой переменной толщины; 2 - контур упругого слоя, прилегающий к жесткой опорной поверхности; 3 - деформированный контур упругого слоя; 4 - недеформированный контур упругого слоя
1 / 2
3
4
Уравнения недеформированного внешнего контура упругого слоя, внутреннего деформированного контура упругого слоя, а также направляющей можно записать в виде:
y' = h+ x 'tga, y' = h0 + x 'tga- a sin ю' x' + Xf (x'), y' = 0. (1)
Здесь h0 - толщина пленки в начальном сечении; a - угол наклона линейного профиля ползуна
к оси Ox'; a и ю' - соответственно амплитуда и частота контурного возмущения на поверхности ползуна, которые характеризуют степень отклонения профиля ползуна от линейного; l - длина подшипника.
Предполагается, что /tga, a и X одного порядка малости; ю=ю' l в дальнейшем определяется из условия максимума несущей способности подшипника; Xf (x') - характеризует деформацию упругого слоя; f' (x') - безразмерная функция, подлежащая определению.
Основные уравнения и граничные условия. В качестве исходных уравнений берется безразмерная система уравнений движения вязкого газа для случая «тонкого слоя», уравнение неразрывности и уравнение состояния:
d"u =1Ф d(Pu) + d(pu) = 0 =р ду2 Л dx ’ dx dy ,
Здесь размерные величины x',y',u', и',p',p ' связаны с безразмерными x,y,u, и,p,p соотношениями:
f i ti f * f * f
x = /x, y = hoy, u = u u, U = SU и, p = PaP,
Г * * P A ^/u h0
P =PP, P =TÍ, Л = ^, S=^ (3)
RT Paho l
где u', и ' - компоненты вектора скорости; p' - гидродинамическое давление; p ' - плотность; ц - динамический коэффициент вязкости; ра - атмосферное давление; Л - параметр сжимаемости; Т - температура; R - газовая постоянная.
Система уравнений (2) решается при следующих граничных условиях:
u = 0, и = 0 при y = h(x) = 1 + rx - r|j sin rnx + r2ф(x),
u = 1, и = 0 при y = 0; p(0) = p(1), r = /tga, Г1 = a, Г2 =y~, ф(x) = f (ex). (4)
h0 h0 h0
Здесь ф( x) - ограниченная функция, подлежащая определению, max |ф x) = 1|.
xe[0,1]
К системе уравнений (2) необходимо добавить следующую безразмерную систему уравнений Ламэ для «тонкого слоя»:
d_u/ = 0 d_uxL = 0 (5)
^ *2 ’ ^ *2 ' '
dy dy
Здесь в упругом слое переход к безразмерным переменным осуществлен по формулам:
u'd, = ü*uy,, u'x, = u*ux,, y' = h0 +5*y*, 5* = h1-h0,
f f w o *
где ux,,uy, - компоненты вектора перемещений; 5 - толщина упругого слоя в начальном сече-
*
нии; u - характерная величина компоненты вектора перемещений. В переменных x,y и x, y* - уравнение недеформированного контура упругого слоя, прилегающего к смазочному слою, соответственно запишем в виде:
y = H1(x) = 1 + rx-r1sinrax, y* = H2(x) = r3x-r4sinrax, r3 = ltga/5*, r4 = X/5*. (6)
Уравнение внешнего контура упругого слоя, прилегающего к жесткой опорной поверхно-
* 1
сти подшипника, запишем в виде у = 1 + г3 х.
Систему уравнений (5) решаем при следующих граничных условиях:
= 5и
ду у*=я2(х)~ ду
ды,
У=Я2(х), М^
дУ
У* =Н2( х) =- Р Ых
* , = 0, ы,
У =1 У
.,= 0. (7)
У =1
1 + а Gы , т Gы 5 Здесь М = - —^—, N = -
(1 -а )5 рд цы 5
где G - модуль сдвига; а* - постоянная Мусхелишвили.
Граничные условия (4) означают прилипание смазки к поверхности ползуна и направляющей, и то, что давление равно атмосферному давлению в начальном и конечном сечениях. Граничные условия (7) означают равенство касательных и нормальных напряжений на недефор-мированной поверхности подшипника, прилегающей к смазочному слою и равенство нулю компонентов вектора перемещений на внешней поверхности упругого слоя, прилегающей к жесткой
опорной поверхности ползуна. Перейдем к решению задачи (2)-(4).
Интегрируя систему уравнений (2), с учетом (4), получаем:
1 dp у2 ( 1 dph 1 1
ы =---—— -\-------—— + — I у +1. (8)
Л dх 2 ^ Л dх 2 h)
Интегрируя уравнение неразрывности от 0 до h, с учетом (5), получаем:
ЦрИ' —р'\ = бЛ^. (9)
dх ^ dх) —х
Это уравнение однократным интегрированием приведем к виду:
—тр>-<')’ (10)
где с1 - постоянная интегрирования.
При очень высоких скоростях, как ы* ^да(Л^да), ёр/—х может оставаться ограниченным, если только ph ^ с1. Это дает асимптоту, соответствующую высокой скорости для подшипника бесконечной длины с газовой смазкой
р = С-. (11)
h
При очень низкой скорости уравнение (9) преобразуется к виду [2]:
—| У? — | = 6Л—, ^х) = 1 + гх -г^бш юх + Г2ф(х). (12)
ёх ^ ёх) ёх
В случае, когда ы* ;А<х>, ы* ^ 0, точное автомодельное решение системы (2), удовлетворяющее граничным условиям (4), будем искать в виде [3, 4]:
ры = дУ + и(х,у)р, 0 = -^ + V(х,у), у = у(5), и(х,у) = ы(5,х), ду дх
У(х,у) = -0(5,х)^, 5 = -^, Р—Р = Щ
^ 7л/-\/*\ Л ✓-/*/* 1л2 /-\/*\ 1л3 / -\/*\
, 5 =-------,--------= 9--------+ ^—, с2 = const. (13)
—х ^х) Л —х h (х) h (х)
Подставляя (13) в (2) и в граничные условия (4), получаем:
—3у ~ д2ы ~ ч до Цх)ды гды Л (,„ ч dh 1 ГЛЛ^
= с2, —^ = с(х),-у-— + 5— = 0, \ h'(х) = —I, (14)
д5 д5 h (х)дх дх ^ ёх
ж
d^(О) = О, ^У(1) = о, u (1, х) = О, и(1, x) = О, u (О, x) = -1, и(О, х) = О,
d £, d £,
i
р(О) = p(1) = 1, ju(^, x)d^ = О. (15)
О
Интегрируя систему уравнений (14), с учетом (15), получаем:
~ ^ 2 Í ^ \
dy-y(fч), u = £у-|í+р\ £ = 6р.
2
Для определения безразмерного гидродинамического давления в смазочном слое имеем уравнение:
р— = Л dx
f 6Р + С2 ^ vh2(x) h3(x),
(16)
с граничными условиями р(0) = р(1) = 1.
Из уравнений (11), (12) и (16) следует, что для определения гидродинамического давления предварительно необходимо найти функцию гф( х), обусловленную деформацией упругого слоя на поверхности ползуна.
Для определения этой функции воспользуемся системой уравнений (5). Интегрируя эту систему, будем иметь:
ыу' = £ (х) у* + £>2 (х), ых, = еъ (х) у* + еА( х).
Используя граничные условия (7) для ыу,, в принятом нами приближении, получаем сле-
дующее выражение:
uy' =“77 У* + 7т(1 + rx )• (17)
y M M
Воспользовавшись приближенной формулой, получим:
(18)
|Уд - Ун| = u
где уд = h(х) - безразмерное уравнение деформированного контура упругого слоя; ун = Н1(х) -
безразмерное уравнение недеформированного контура, прилегающего к смазочному слою.
С учетом (6) и (8) для ^2ф( х), в принятом нами приближении, получим выражение:
. . p p ._п.
л2ф( х) =---1---ЛлЭт юх. (19)
2 M M 4
В экстремальном случае (при Л ^ да) для определения р, приходим к уравнению
р2 1
—(1 + sin юх) + р (1 + ^х - ^1 sin юх) -1 -~ = О. (20)
Решение уравнения (20) будем искать в виде:
1 1
Р = РО +—Р1 +—^ Р2 +... (21)
О M 1 M2 2
Подставляя (21) в (20) для определения коэффициентов разложения (21) с точностью до
членов О | —1— |, будем иметь:
IM2 )
р0 (1 + rx -r^sin qx ) = 1, (22)
р1 (1 + rx - r1 sin Qx) + p2 (1 + r4 sin ШІ) = 1. (23)
Решая (22) и (23) с точностью до членов 0(ц2),0(ц2),0(ц2), получаем:
p0 = 1 - цх + r|j sin юх,
Pj = 2цх - 2ц1 sin ах - ц4 sin юх.
(24)
(25)
Во втором экстремальном случае (т.е. при малых скоростях скольжения направляющей), а также при промежуточных значениях скоростей скольжения направляющей, решение уравнений (12) и (16), с учетом (19), находится в виде ряда Тейлора в окрестности точки х=0:
р = р(0) + р' (0) X + р" (0)у +... (26)
При решении уравнения (12) будем исходить из следующих выражений:
р(0) = 1, Л(0) = 1 + -1, h"(0) = ^-^1®+ Р-(0) + ^4ю,
М ММ
"" =- 3h" (0) р" (0) + 6ЛИ' (0)
h(0) h3(0) ,
где р" (0) определяется из условия р(1)=1. При решении уравнения (16) исходными выражениями являются:
р(0) = 1, К0) = 1 + -1, h"(0) = л-Л1® + + р(0)^4®,
М ММ
Р ' (0) =
Л
Р'' (0) = -Л-РД Р (0)
6
h2(0) hi{0)
+
Р(0)
Л
6
Л
+
h (0) h (0)
Р(0)
6 p' (0)h2 (0) - 2h(0)h ' (0)6 p(0) - 3h ' (0)c2 h^Ó)
1
Здесь в принятом нами приближении с2 является корнем уравнения р'(0) + р"(0)— = 0.
Перейдем к определению основных рабочих характеристик подшипника в первом экстремальном случае, т.е. при Л^да.
При значениях xє
0,1
3
смазочный слой имеет форму сужающегося клина, а при
x є
i,i
3
форму расширяющегося клина. Нагрузка, вносимая сужающимся клином единичной
длины для подшипника с жесткой опорной поверхностью, определяется выражением:
1 1
3 . 3л ^ , л 2л1
J( p - ^ = Ц -^x+л sin^2~x Idx = -18+-3^.
18 3л
(27)
Воспользуемся основными результатами работы [5], где показано, что оптимальное значение ю (по несущей способности) составляет ю=3л/2.
Нагрузка, определяемая расширяющимся клином
Wp = {(Р (3)-P(x) dx= ^ -3л+л
: = |1 --Л + Л |2-2 + Л^ = 2Л + 2П, 1 3 13 3 2 18 9 3
(28)
Таким образом, нагрузка вносимая одновременно сужающимся и расширяющимся клином, определяется следующей суммой:
Л 2л 2 2 л 2 2л1
w = w + w =-------L + —+ — л + — Л1 = —+ —Л +—-
с Р 18 3л 9 3 6 3 3л
2
3
При ю = 0, т.е. в случае ползуна с линейным контуром опорной поверхности, будем
иметь:
Wo =1(1 - р( X) )ск = -. (30)
0 2
Из полученных выражений (29) и (30) следует, что w « 2w0, т.е. предлагаемая конструкция подшипника с жесткой опорной поверхностью обладает свойством подшипника двойного действия. С учетом (8) для безразмерной силы трения получим выражение:
^ <31>
0 1 + ЦХ - Б1П — ПХ
Перейдем к определению добавочной безразмерной составляющей поддерживающей си-
3
лы и силы трения при Ю = ^ .
С учетом (25), будем иметь:
і
М0 У 18 3л 3л )М
» к -
рдоб М -1
- -
М
Ріі з і-р(х)
, і 4 4 2л1 2 ) 1
Ух -'—л—л1-1 —л4 і—, (33)
^'93" л 3 4|М
і 5 4 10 2 . .
».до, + »„о, -|-3л-3Лі-3Л-3л< IМ■ (34)
Таким образом, с учетом (29) и (34):
л 2 2л, 1 і 1 4 10 2 .
- 6+^л1+ 1ЛГ- МІ3П+^+ 3Л+Ї^І, (35)
, -ч+П-М+^-З+М+Зц (36)
трп 2 3л М У 2 3л 3л у
Численный анализ полученных зависимостей (35) и (36) для основных рабочих характеристик позволяет сделать следующие выводы.
1. При больших значениях скорости скольжения направляющей, таких, как Л^-да при
3
ю-—л и при М ^да, подшипник с жесткой опорной поверхностью обладает свойством подшипника двойного действия.
2. С увеличением значения упругогидродинамического параметра М значения безразмерной несущей способности и безразмерной силы трения увеличиваются, оставаясь меньше, чем соответствующие значения для подшипника с жесткой опорной поверхностью.
3. Как и предполагалось, при М ^да найденные выражения для основных рабочих характеристик подшипника с деформируемой опорной поверхностью стремятся к соответствующим выражениям для подшипника с жесткой опорной поверхностью.
4. Найденные закономерности влияния параметра осцилляции ю, параметра сжимаемости Л, а также упругогидродинамического параметра М на основные рабочие характеристики остаются такими же и при малых и промежуточных значениях скоростей скольжения направляющей.
Заключение. По результатам теоретических исследований разработан метод расчета упорного газодинамического подшипника с повышенной несущей способностью с учетом деформации его опорной поверхности. Дана оценка влияния параметра осцилляции ю, параметра сжимаемости Л,
3
а также упругогидродинамического параметра М на основные рабочие характеристики подшипника. Установлено, что с увеличением значения параметра М несущая способность и сила трения возрастают, оставаясь по своему значению меньше, чем соответствующие значения для подшипников с жесткой опорной поверхностью.
Библиографический список
1. Heuer D.F. Dynamic and Environmental Evaluation of Compliat Foil Gas Lubricated Bearings / D.F. Heuer, R.A. Collins // Tech. Rep. AFAPL-TR-73-56. Air Rorce Aero Propulsion Lab. - June, 1973.
2. Bisson and Anderson. Advanced Bearing Technology // NASA SP-38, 1965. - Ch.5. - Gas Lubricated Bearings by J.S. Ausman. - Р.109-138.
3. Математическая модель гидродинамической смазки бесконечно широких опор, работающих в турбулентном режиме на микрополярной смазке / К.С. Ахвердиев [и др.] // Трение и смазка. - 2007. - №6. - С.278-284.
4. Математическая модель гидродинамической смазки бесконечно широких опор, работающих в нестационарном турбулентном режиме на микрополярной смазке / К.С. Ахвердиев [и др.] // Вестн. Рост. гос. ун-та путей сообщения. - 2007. - №4. - С.18-24.
5. Ахвердиев К.С. Газовый упорный подшипник с повышенной несущей способностью / К.С. Ахвердиев, В.А. Константинов, С.А. Солоп // Новые материалы и технологии в машиностроении: сб. науч. тр. по итогам междунар. науч.-техн. конф. - Брянск, 2009. - С.3-8.
Материал поступил в редакцию 17.11.10.
References
1. Heuer D.F. Dynamic and Environmental Evaluation of Compliat Foil Gas Lubricated Bearings / D.F. Heuer, R.A. Collins // Tech. Rep. AFAPL-TR-73-56. Air Rorce Aero Propulsion Lab. - June, 1973.
2. Bisson and Anderson. Advanced Bearing Technology // NASA SP-38, 1965. - Ch.5. - Gas Lubricated Bearings by J.S. Ausman. - R.109-138.
3. Matematicheskaya model' gidrodinamicheskoi smazki beskonechno shirokih opor, rabotayu-schih v turbulentnom rejime na mikropolyarnoi smazke / K.S. Ahverdiev [i dr.] // Trenie i smazka. -2007. - №6. - S.278-284. - In Russian.
4. Matematicheskaya model' gidrodinamicheskoi smazki beskonechno shirokih opor, rabotayu-
schih v nestacionarnom turbulentnom rejime na mikropolyarnoi smazke / K.S. Ahverdiev [i
dr.] // Vestn. Rost. gos. un-ta putei soobscheniya. - 2007. - №4. - S.18-24. - In Russian.
5. Ahverdiev K.S. Gazovyi upornyi podshipnik s povyshennoi nesuschei sposobnost'yu / K.S. Ahverdiev, V.A. Konstantinov, S.A. Solop // Novye materialy i tehnologii v mashinostroenii: sb. nauch. tr. po itogam mejdunar. nauch.-tehn. konf. - Bryansk, 2009. - S.3-8. - In Russian.
K.S. AKHVERDIYEV, M.A. MUKUTADZE, V.A. KONSTANTINOV
GAS THRUST BEARING OF INCREASED BEARING CAPACITY WITH YIELDING BEARING SURFACE
Basic operation factors of the deformed gas-lubricated thrust bearing with increased bearing capacity are theoretically analysed. It is shown that steady-state operation characteristics of this bearing differ from those of the rigid bearing.
Key words: thrust bearing, bearing capacity, adapted profile, gas lubrication.
АХВЕРДИЕВ Камил Самедович (р. 1938), заведующий кафедрой «Высшая математика-2» Ростовского государственного университета путей сообщения, доктор технических наук (1984), профессор (1984). Окончил Азербайджанский государственный университет (1962).
Область научных интересов: трение и износ в машинах.
Автор 460 научных публикаций.
МУКУТАДЗЕ Мурман Александрович (р. 1963), доцент кафедры «Высшая математика-2» Ростовского государственного университета путей сообщения, кандидат технических наук (1995), доцент (2001). Окончил Ростовский государственный университет путей сообщения (1986). Область научных интересов: трение и износ в машинах.
Автор 41 научной публикации.
КОНСТАНТИНОВ Владимир Артёмович (р. 1986), аспирант кафедры «Высшая математика-2» Ростовского государственного университета путей сообщения. Окончил Ростовский государственный университет путей сообщения (2008).
Область научных интересов: трение и износ в машинах.
Автор 6 научных публикаций.
Kamil S. AKHVERDIYEV (1938), Head of the Higher Mathematics-2 Department, Rostov State Transport University. PhD in Science (1984), Professor (1984). He graduated from Azerbaijan State University (1962).
Research interests: friction and machine wear.
Author of 460 publications.
Murman A. MUKUTADZE (1963), Associate Professor of the Higher Mathematics-2 Department, Rostov State Transport University. Candidate of Science in Engineering (1995), Associate Professor (2001). He graduated from Rostov State Transport University (1986).
Research interests: friction and machine wear.
Author of 41 publications.
Vladimir A. KONSTANTINOV (1986), Postgraduate student of the Higher Mathematics-2 Department, Rostov State Transport University. He graduated from Rostov State Transport University (2008). Research interests: friction and machine wear.
Author of 6 publications.