Нестационарная расчетная модель неоднородного двухслойного пористого подшипника конечной длины при комбинированной подаче смазки

Нестационарная расчетная модель неоднородного двухслойного пористого подшипника конечной длины при комбинированной подаче смазки

К.С. Ахвердиев, И.М. Елманов, Н.С. Задорожная, А.М. Мукутадзе, Е.В.

Поляков

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

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

Ключевые слова: двухслойная пористая втулка; проницаемость пористых слоев; несущая способность; уравнение движения шипа; устойчивость.

Введение.

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

В большинстве опубликованных работ, посвященных расчету пористых подшипников, рассматриваются подшипники либо бесконечной длины [1 - 3], либо пористые подшипники конечной длины [4 - 7], работающие в стационарном режиме трения. В существующих нестационарных расчетных моделях пористых подшипников конечной длины не учитывается не только анизотропия проницаемости пористого слоя, но и влияние подачи источника смазки [8 - 10]. В расчетной модели, приведенной в работах [11, 12], хотя и учитываются анизотропия проницаемости пористого слоя и влияние подачи смазки либо в осевом, либо перпендикулярном оси подшипника, направлениях, однако здесь не рассматривается комбинированная подача смазки, и пористая втулка считается однослойной. Для повышения несущей способности подшипника, возникает необходимость комбинированной подачи смазки и конструкция подшипника с двухслойной пористой втулкой, каждая из которых имеет свое функциональное назначение.

Постановка задачи.

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

кг' = кгв V У, I = 1,2.

Здесь к У - проницаемости пористых слоев; к - характерные проницаемости; ф

- известная безразмерная функция г ; I - длина подшипника.

Уравнения, определяющие течения смазки в пористых слоях, согласно [8], представляются в виде

+ + «£^1 = 0, ! = 1,2,

ду дг дг дг

где у, г - прямоугольные координаты (рис. 1); р* - гидродинамическое давление в пористых слоях.

Давление в смазочном слое (между шипом и подшипником) определяется из модифицированного уравнения Рейнольдса [14] в рамках модели короткого подшипника

д_

ді

Н3 8-РЛ

V & у

= ЄЦ

СІф

&

12Ц)| у=0,

где Н = С(1 + єеоБв) - толщина пленки смазки; С - радиальный зазор; є -относительный эксцентриситет; в - угловая координата; р - давление в пленке смазки; ц - динамический коэффициент вязкости; оъ, (, ( - угловые скорости

подшипника, шипа и нагрузки; р - угол положения, ґ -время; у0 - координата

скорости в направлении Оу на внутренней границе пористого слоя, прилегающего к зазору,

ц

/л *\ др,

ду

V у

р !

і = 1,2.

у=о

Рис. 1 Радиальный подшипник конечной длины с пористой обоймой

Для решения системы уравнений

* * ^ * *

Р2 = Рі = р = Рн при 1 =—- ; р2 = Рі = р = Рк при 1

* * ~

Р1 = Р при у = 0; Р2 = Р§ при у = — Н ;

Ф* Л *

1 7 дР2 тг

Р1 = Р2,Л1^ = к2—^ при у = —Н1. ду ду

Здесь Н - толщина двухслойного пористого слоя; Н1

прилегающего к смазочному слою.

Перейдем к безразмерным переменным по формулам

Ь_

2

толщина слоя,

P

*

p,c

2

MR0 ] ^О^-'і

= Pjf2

g

P =

pc

2

кВ

C

2 z y

Z = —, Y = ^~, T = rn.t, L В ] '

2

juR2j * J JURo} ' * Л JuRbWj

Установим закономерность подачи смазки и изменения проницаемости пористых слоев по координате z следующим образом

к= к.к, к = ев , в< О, Pg = aZ + b + Pg

‘‘ g <5'/<5

С учетом

д2 P *

dY2

+ 4^ —1 -^4- + 4eZ — = О,

V L ) dZ2 dZ

д2 P

дZ2 (1

+ є cos

e)3

ґ 1 л

(р

V

2

slne + є cose

+

3Ф

ґ * \ dP Л

(l + є cos в)3 ((L)2

oy

v у

,

Y=О

где ^ = 2Р0; точкой обозначено дифференцирование по Т . Граничные условия (5) примут вид

* *

Р1 = Р при 7 = 0; Р2 = Р^ при 7 = -1;

*

Р = Р = Ря при X = -1; р = Р = РК при X = 1,

эр:

д/у

oy

Pl* = P*, k^ = k2 ^ при Y = -^ .

oy

Ві

В

Решение систем уравнений (S) и (9) будем искать в виде P1* = B1Y4 + B2Y3 + B3Y2 + B4Y + aZ + b + P1 (в, Z),

P2* = A1Y3 + A/2 + A3Y + aZ + b + р(в, Z),

P = aZ + b + P1(0,Z), Pg = aZ + b + Pg(z2 -1) Pg Подставляя (11) в (8) и (9), будем иметь Л = A1 + A3 - P1 + Pg (z2 - 1);

= const

4

В

\L J

22

д2 A

dZ 2

1 +

В1

3\

V В у

Ґ и \

1 +

В1

2

V

V В

+

+ 4pZ

dPg 1

dZ 3

-1+

В1

V В у

3Л

+ 4PZ

-1 +

Ві

В

a = О:

A

В

1+

+2

1+—1

В j

+ 2 A

1+—1 В

2P

1+—

В j

+

+2 p (z :

i) - і+

+-

д2 p

дZ2

+ 4PZ

г-1+в

В

—

Hi и j

f

-1+

\ Ґ u\

+ 4

H

L

д2 A.

1 +

4

1

+ —

з

г Hi

и

з

д2 P„ г

+ -

дЛ.

дZ

Ґ и \

1 +

H±

v и J

4

1

+ -

з

J_

/

v

дZ2

-1 +

г—

и

з

;> +

Ґ и \

1 +

v и j

з

+

1 +

v

Вл

и j

з

+

+

д4

дZ

1 г г—і 1 31 1 г Hi 1 21

-1 + + — 1 + +

3 v v H J J 2 v H J J

+ •

дЛ

дZ

1 +

- І +

г — і3 Л

v — j

+

дР.

дZ

- І +

Hi

— j

- І +

г — і3 Л

v — j

+

+

І дР

3 дZ

І +

и.

V—ly

з

0;

B - — A •

в4 7 A3;

k1

A1

Hi

v — j

+ Л

Hi v И J

И. _ „ г и.

- A3 — — B. з В 1

v — j

-B

Hi v — j

+ BQ

И. v И J

-B

И.. и ’

4 B.k.

ҐИ±'

и j

+3B2k

ґ Hl

и J

2B3k.

ґ Hl и J

з A1 k2

f Hi 1

v И J

2k

Ґ u\

з

- 4 B1

д2 B2 І

И.

V — y

Ґ TJ \

+ 3B0

v — j

- 2 B.

v — j

+ 4

I

v L J

В

1 fHi і5

5 дZ2 I И

v — j

+

V11 /

дZ2 4

Hi

v — j

д2 B3I

дZ2 3

f—i! v—J

д2 B4 1 г—i 1 2

дZ2 2 v И J

+

+ 4ytfZ

1 дв. г И.л

дв4 1

дZ 2

5 дZ

v — j

v — j

+

І дв.

4 дZ

дР.

В

дZ

v — j

-a

f и.!

v—J

и. I H J

4

І дв,

3 дZ

+

— 0

з

2

2

2

2

з

5

Из уравнения (13) определяется А3, функции А1 и А2 найдутся из уравнений (12) и (14), затем определяется гидродинамическое давление Р1. Функции В1, В2, В3 определяются из уравнений (16). При определении гидродинамического давления Р1 явный вид функций А1, А2, В1, В2, В3 нам не понадобится. Решая уравнение (13) учитывая (15) имеем

В 4 - -

Pse к2 g f!-f H1T1 v — у V v y у (z4 -1)+ 2Paf-1 + —11 V H у

Y 1 1V L у 3 f -1+ v f H1 Y ' tfjj 1 + 2 f f H \21 [1+f—1 1 V vH))

\z3 - Z)

Решая уравнение (9) с учетом граничных условий (10) и учитывая (17), для Р1 получим следующее выражение

Pi

12

Z + scosO)3

1

р — 2

sinO + s cosQ

z 3 -11

Фкп

6k1 (l + scosO)3 [ —

•X

L

вРг

-1 +

f H1 13 Л

v Ну

X-

f Z 6 - Z2 1

30- T

7

+-----

15

+ 2fia

Z5 Z Z3 Z

----------------------------1-----

20 20 6 6

у

-1 +

f H1 13 Л v Ну

1

+---

2

1 +

f H1 12 Л v Ну

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

определяется из условий

s cos 01 + sp sin O1 - 0,

s sin 62 - sp) cos 62 - 0, 62 - O1 + n.

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

F(e)-- 41 / f fPcosQdQdZ. 2C2 J1 f

F

-1 O

В случае полного заполнении зазора смазкой имеем

( ) ЦК(Л Ю Г Г

р(е) =_^ 0 у Г ГрсоМвсИ 2С } }

-1 0

2С -1 0

р - р

где Р = а2 + Ь + Р1 определяется уравнением (18), а = к ^ н .

Задача на устойчивость шипа в подшипнике.

Для определения движения шипа применим безразмерное уравнение записанное в следующем виде:

а є ре

- +

ат2 ю умс

(ю я 1 (арЛ 2

СОБ Р+ЄІ —

і ат

2

а2 р рР і(ю я | . 2 (ає у ар

ат2 ю умс є

2 (ає У ар^

бій р-

ю. в V ат у\ ат

V 1 У

где М - масса ротора, F(е) и F(р) - усилия масляной пленки в случае неполного заполнения смазкой зазора, определяемые формулами (18)-(20), соответственно, в случаях подачи смазки в направлении, перпендикулярном оси подшипника, и в осевом направлении.

Уравнения (21), определяющие движение шипа, решаются численно с учетом

полученных данных (19)-(20). Компоненты ускорения аТт, аТ* представляют собой явные функции параметров в и р, —, —, Р , Рн, Рк,вх,в2, Ф, в, —, а .

ат ат х к1

Уравнения (21) записывается в стандартной форме первого порядка и решается с помощью метода, разработанного Гиром [13].

Как и в работе [9], после получения решения уравнений движения, устойчивость рассматриваемого движения определяется визуально по графику. При заданных значениях выше указанных параметров, области устойчивости приведены на рис. 2. Здесь все точки, которые лежат ниже кривых устойчивости, соответствуют устойчивому движению шипа, а все точки, которые лежат выше кривых, соответствуют неустойчивому движению

(ю = 4ё/с), где g - ускорение силы тяжести.

Из приведенных на рис. 2 зависимостей следует, что:

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

направлениях подачи смазки.

2. Подшипник при комбинированной подаче смазки работает более устойчиво, чем при отдельно осевой или перпендикулярно оси направлениях подачи смазки.

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

4. В случае, когда отношение проницаемостей пористых слоев больше единицы

ґк "] к

— >1 , подшипник работает более устойчиво, чем при соотношении — <1.

\кі ) к

і. ф = одаь ^ = 1,5? 2, ф = ^ = і,:? ф = о,аз? ^ = ідг , ф = еда ^ = і*

Ь,

3.

Й.

А,

К

Рис. 2. Схематическое изображение границ устойчивости (двухслойный пористый подшипник при полном заполнении смазкой зазора и комбинированной подаче

смазки)

Литература:

1. Ахвердиев К.С., Шевченко А.И., Слоистый пористый подшипник бесконечной длины [Текст] // Вестник РГУПС. - 2000. - № 1. - С. 19-26.

2. Ахвердиев К.С., Мукутадзе М.А., Новгородова В.С., Черкасова Т.С.

Гидродинамический расчет двухслойного пористого подшипника бесконечной длины с учетом анизотропии проницаемости пористого слоя и сил инерции [Текст] // Вестник ДГТУ. - С. 211-215.

3. Мукутадзе М.А., Александрова Е.Е., Константинов А.А., Шевченко А.И. Гидродинамический расчет радиального пористого подшипника бесконечной длины с повышенной несущей способностью с учетом сил инерции [Текст] // Вестник РГУПС. - 2012. - № 2(46). - С. 194-197.

4. Дерлугян Ф.П., Щербаков И.Н. Обоснование процесса получения

композиционных антифрикционных самосмазывающихся материалов с

заданными техническими характеристиками методом химического наноконструирования. [Электронный ресурс] // «Инженерный вестник Дона», 2010 г., №4 - Режим доступа: http://ivdon.ru/magazine/archive/n4y2010/287

(доступ свободный) - Загл. с экрана. - Яз. рус.

5. Ахвердиев К.С., Прянишникова Л.И., Пустовойх Ю.И., Гидродинамический расчет пористых подшипников с переменной проницаемостью вдоль оси с учетом нелинейных факторов [Текст] // Трение и износ. - 1993. - Т. 14, № 5, -С. 813-821.

6. Ахвердиев К.С., Прянишникова Л.И. Об одном точном решении задачи о радиальном пористом подшипнике конечной длины [Текст] // Трение и износ.

- 1993. - Т. 12, № 1, - С. 24-32.

7. Ахвердиев К.С., Приходько В.М., Шевченко А.И., Казанчян О.Р. Математическая модель течения смазки в зазоре радиального подшипника конечной длины со слоистым пористым вкладышем переменной толщины [Текст] // Проблемы машиностроения. РАН М.: Наука. - 2000. - № 6, - С. 85 -91.

8. Кузано К., Фанк П.Е.. Исследование коэффициента передачи упругой опоры качения в демпфере со сдавливаемой пленкой и пористой обоймой [Текст] //Проблемы трения и смазки, изд-во «Мир». - 1974. - № 1, - С. 54.

9. Конри, Кузано К. Об устойчивости пористых радиальных подшипников [Текст] // Конструирование и технология машиностроения. - 1974. - № 2. - С. 206-216.

10.Ахвердиев К.С., Муленко О.В. Об устойчивости двухслойных пористых радиальных подшипников [Текст] / // Вестник РГУПС. - 2002. - № 3.- С. 5-7.

11.Ахвердиев К.С. Мукутадзе М.А., Задорожная Н.С., Флек Б.М., Поляков Е.В., Расчетная модель гидродинамической смазки неоднородного пористого подшипника конечной длины, работающего в устойчивом нестационарном режиме трения при наличии принудительной подачи смазки [Электронный ресурс] // «Инженерный вестник Дона», 2013 г., № 3. - Режим доступа: http://ivdon.ru/magazine/archive/n3y2013/1765 (доступ свободный) - Загл. С экрана. - Яз. Рус.

12.Ахвердиев К.С., Мукутадзе М.А., Задорожная Н.С., Флек Б.М. Аналитическое прогнозирование передаточных характеристик центрально нагруженного демпфера со сдавливаемой пленкой и пористой обоймой с учетом влияния анизотропии проницаемости пористого слоя и источника смазки [Текст] // Вестник РГУПС. - 2013. - № 4 (52). - С.131-142.

13.Gear C.W., Numarical Initial Value Problems in Ordinary Differential Equations, Prentice-Hall, Inc., Englewood Cliffs. - N.J., 1972.

14.Reynolds, O. On the theory of lubrication and its application to Mr. Beauchamp Towers experiments / O. Reynolds. - Phil. Trans. Roy. Soc. - London, 1886, vol. 177, pt. 1.