Расчетная модель радиального подшипника повышенной несущей способности, работающего на микрополярной смазке с учетом ее вязкостных характеристик от температуры Текст научной статьи по специальности «Строительство и архитектура»

МАШИНОСТРОЕНИЕ И МАШИНОВЕДЕНИЕ MACHINE BUILDING AND MACHINE SCIENCE

УДК 51:621.891 DOI 10.12737/23784

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

К. С. Ахвердиев1, Е. О. Лагунова2, К. С. Солоп3

1 2- 3 Ростовский государственный университет путей сообщения, г. Ростов-на-Дону, Российская Федерация

Design model of heavy-duty radial bearing running on micropolar grease with account of its viscosity-temperature characteristics *

K. S. Akhverdiev1, E. O. Lagunova2, K. S. Solop3**

1 2- 3 Rostov State Transport University, Rostov-on-Don, Russian Federation

Целью данного исследования является описание метода расчета радиального подшипника скольжения с адаптированным профилем опорной поверхности, работающего на микрополярном смазочном материале. Кроме того, предполагается, что вязкостные характеристики микрополярной жидкости зависят от температуры. Рассматривается частичное заполнение смазочным материалом рабочего зазора (т. е. наличие свободной поверхности). Кроме обычных безразмерных параметров теории ньютоновской жидкости вводятся и другие параметры для микрополярной жидкости: взаимодействие N и N (характеризует зависимость от размера молекул смазки). В результате численного анализа показано, что с увеличением теплового параметра К несущая способность подшипника уменьшается. С увеличением параметра N и параметра ю, характеризующего адаптированный профиль опорной поверхности, несущая способность увеличивается. При значении параметра ю = 0,5 достигается максимум несущей способности. При N ^ ж значение несущей способности стремится к соответствующему значению для случая ньютоновского смазочного материала.

The study objective is to develop a method for calculating the plain journal bearing with an adapted profile of the bearing surface operating on the micropolar lubrication. Furthermore, it is assumed that the micropolar fluid viscosity characteristics are temperature dependent. Partial filling of the front gap with lubricant (i.e. availability of a free surface) is considered. Except regular dimensionless parameters of Newtonian fluid theory, other options for micropolar liquid are introduced, namely the interaction of N and N1 parameters characterizing the dependence on the size of the grease molecules. The numerical analysis shows that with the heat parameter (K) increment, the carrying capacity of the bearing decreases. With the increment of the parameters N1 and ra characterizing the adapted profile of the bearing surface, the bearing capacity increases. At the parameter point of ra =0.5, maximum load capacity is achieved. At N1 the bearing capacity value

tends to the corresponding value for the case of Newtonian lubricant

Ключевые слова: радиальный подшипник, несущая способность, микрополярный смазочный материал, зависимость вязкостных характеристик от температуры, адаптированный профиль

Keywords: radial bearing, load capacity, micropolar lubricant, temperature dependent viscous characteristics, adapted profile.

Введение. Как известно, в настоящее время в качестве модели гидродинамической смазки подшипников скольжения широко используется микрополярный жидкий смазочный материал. Поэтому разработка методов расчета 5 подшипников скольжения, работающих на микрополярном жидком смазочном материале, требует учета зависимости -2 вязкости не только от давления, а также зависимость вязкостных характеристик от температуры. Значимый недоста-§ ток существующих методов расчета подшипников скольжения, работающих на микрополярном жидком смазочном 2 материале, состоит в том, что в большинстве случаев или вовсе не учитывается зависимость вязкости от давления и 'й температуры [1-7], или учитывается эта зависимость только от давления [8-11].

¡и Постановка задачи. Рассмотрим установившееся движение микрополярного смазочного материала в зазоре

радиального подшипника. Предполагается, что подшипник с адаптированным профилем опорной поверхности непо-^ движен, а вал вращается с угловой скоростью О (рис. 1).

*Работа выполнена по гранту ОАО «РЖД» на развитие научно-педагогических школ в области железнодорожного транспорта.

**E-mail: vm@rgups.ru, lagunova@rambler.ru

The research is done on OJSC "RR" grant on the development of academic schools in the sphere of railway transport.

Рис. 1. Расчетная схема радиального подшипника

Предполагается, что имеет место частичное заполнение смазочным материалом рабочего зазора (т. е. наличие свободной поверхности). Кроме того, предполагается, что вязкостные характеристики микрополярной жидкости зависят от температуры.

-р 'Г

-р 'Т'

-р 'Т'

ц' = ц0е - , к' = к0е Т , Y = Yo* С1)

Здесь ц0 — характерная вязкость ньютоновского смазочного материала; к 0 и Yo — характерные вязкости микрополярного смазочного материала;

T' — температура, в' — экспериментальная постоянная величина.

Уравнение контуров вала и подшипниковой втулки в полярной системе координат r' , 9 в центре вала запишутся в виде:

r = r0, r = r1 + e cos9- а sinю9 , (2)

где го — радиус шипа; г1 — радиус подшипника; е — эксцентриситет; —, — — малая величина одного порядка (

5 5

5 = г1 - г0); ю — параметр контактного профиля подлежит определению.

Основные уравнения и граничные условия. В качестве исходных уравнений берется система безразмерных уравнений движения микрополярной жидкости для случая «тонкого слоя» с учетом (1), а также уравнение неразрывности:

д и т до 1 dp --+ N2 — =--—

dr2 dr e-вт d9

д2 и и 1 ди ди dv o dr2 = N1 N1 дr ' dr д9 = '

(3)

Здесь размерные величины и',V, и', р', ц', к', у',т',Т', г' связаны с соответствующими безразмерными и, V, и, р, к, у, Т, г следующими соотношениями:

и' = иОг0; V' = Q5v; у' = у*у; р' = р* р;

д'=до д; к'=ко к; у'=7о у; Т' = ТоТ; г' = г +5г;

* '0 и = 0

r Q * r02Q( 2^0 + К0)

25

Р =■

252

5 = r1 -r0;

(4)

N2 =-

N1 = /_2И0; 12 = J0_

1 52 кп

2До + к о 5 2 Ко 4Цо

Здесь и' , V — компоненты вектора скорости; и' — скорость вращения микрочастиц, То — характерная температура.

Как видно из системы (3), помимо обычных безразмерных параметров, встречающихся в теории ньютоновской жидкости, здесь вводятся и другие параметры для микрополярной жидкости. Параметр взаимодействия N ^ о при ко ^ 0 .

(U

К X <и

(U

и о X

к

3

й

(U

К X <и О Л

н о о X

к

3

й

Безразмерный параметр Ы1 зависит от параметра I, имеющего размерность длины, и его можно трактовать как характеристику, зависящую от размера молекул смазки. Система уравнений (3) с точностью до членов О

решается при следующих граничных условиях:

и = 0, V = 0, и = 0 при г = И (е) = 1 + sinю6;

и = 1, V = 0, и = 0 при г = 0;

р = (о) при е = е1 и е = е2; (5)

е а

П = 8; П = 8;

91 и 02 — соответственно координаты начала и конца свободной поверхности. Осредним второе уравнение системы (3) по толщине смазочного слоя. Получим:

1 h d2v , 1 h , 1 h du , — J ——dr = —- J vdy Л--- J —dr

h о dr

N h о

N1 h 0 dr

Найдем решение уравнения (6) в виде:

и=а1 (е) г2+а2 (е) г+а (е).

Из граничных условий (5) следует, что

A3 = 0; A2 = -A1 h.

С учетом (8) для и получим выражение:

и = A1 (е)-(r2 -rh).

(6)

(7)

(8) (9)

Подставляя (9) в (6) с точностью до членов O

N

по[„.±

\ i л \

O

1 У

V N1 У

получим:

и = —

— (r2 - 2h), 2 N v '

do

1

dy 2n

(2r - h).

(10)

С учетом (10) система уравнений (3) в принятом нами приближении преобразуется к следующему виду:

- И ) = евт^Р,

' dе

d2u N 2

dr 2 2Nj

1 i

и = - --(

2Nj у

du dv ,

— + — = (

де dr

(11)

Точное автомодельное решение системы (11) с граничными условиями (5) будем искать в виде:

u = ё1 + u (r, е), V = Л v (r, е),

dr де v '

с о тз

"¡3

-м

ХЛ

(U &

С Л

F(r.e) = -v(^)/U(r,Q) = u(l-)), 4 = dp _ Сх С2

h(ey

v(r,e) = m ^ =

d е

и = —

—(52 - 5).

2 N \ '

h2 (е) h3 (е)

(12)

С учетом (12) из системы (11), получим:

г

Г = с2

ы1

йр = е_Рт й 0

2ЫХ

С| С*2 к2 (0)+к3^)

Граничные условия (5) примут вид:

у'(0) = 0, у'(1) = 0; й (0) = 1,

/7(1) = 0, У(1) = 0,

?(0) = 0,

(13)

(14)

Решая задачу (13) - (14) непосредственным интегрированием, получим:

у ' 2 2ЫХ I 3 2

^ ( N2

12 N1 2

и = —

- — (§2 " §). 2 N У '

(15)

~ N ~

Здесь С] =6--, а С2 в дальнейшем определим из условия ) = /'("г ) = " •

N1

Как и ожидалось, при N ^ ж, V ^ ж полученные результаты полностью совпадают с результатом для случая ньютоновской смазки.

Для определения гидродинамического давления имеем:

С,

СРГ Ф _ С\

й0 к2(0) к3(0)

(16)

1

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

скорости диссипации энергии под действием сил сдвига. Получим

2ц0цЛ2г02к 1

И00~ 8 0

V

к2(0) к(0)

Повышенная температура будет определяться выражением.

(17)

й0 й0 ' с,0 с,0 ' 8 0[к2 (0) к(0)

Здесь б — расход смазки в единицу времени; с — теплоемкость при постоянном давлении.

1 -80г С

0 12

Продифференцируем по 0 зависимость ц от температуры Т:

йц _ йТ — = -цр— •

й0 й0

Комбинируя (18) - (20), получим:

1 ф_24РцоОго/7(0)1/ ф'(4) |

ц2(0) й0 Введем следующие обозначения:

С2ср82Т0 1[/?2(0) ' А(0)

К =

24РЦ0

к1СрТ0

(18)

(19)

(20)

(21)

щ

5 I

и «

и т

0

1 5

3

ей

Щ

5 I

(и О

а н о

0

1 5

3

ей

е

Ai = = % д2 = Ь^К^И = jcA,

12

N2 . N2 N4

2 '

Д, = \(й'У (11 =-+ 4--+

о 720 12Л^

Здесь к — тепловой параметр. С учетом (22) уравнение (21) преобразуется к виду:

1 ф_ К ( ¿1 + А2 | Аз ^

n2(e)rf0 сЦл3(е) л2(е) л(е)

Интегрируя уравнение (23), получим:

ц (0) = -

1 - К [ J3 (0)Д1 + J2 (0)Д2 + J (0)Д3 ]

Jк (0)=i

0 d0

С

1 hK (0)

Гидродинамическое давление p с учетом (24) определяем из уравнения (16).

p = q) = qj2 (0)+c2J3 (0) =

01 h2 (0) 01 h3 (0)

6 —

N

2

(22)

(23)

(24)

N j{(0-01 )-л(0-01 )2-^cos®(0-01 )-

-01 )2 +^(02 -01 )3 + ^ (02 -01 ) cosю(02 -01 ) 2 ю

)-01)-3^(0-01 )2 -— cosю(02 -01) 2 ю

Перейдем к определению основных рабочих характеристик радиального подшипника. Для безразмерной несущей способности и безразмерной силы трения с учетом (25) получим:

rn3 Q( 2цп + к0 ) е2 гЩ Q( 2цп + к0 ) е2 Ry = - 0 V -0J- J p sin 0d0, Rx = - 0 -^ J p cos

(25)

2 5

L = ^0

^ 5

2 5

!? C9 ¿TO

i-

h2 (0)

r7N 2 12 N

- 4

i hM

(26)

Входные параметры для расчета несущей способности и силы трения, определяемые выражением (22) - (26), приведены в таблицах №1 и №2.

Таблица 1

Входные параметры для расчета несущей способности и силы трения

J1 J 2 J3 П=П1 N N1 ю в k 02 -01

7,36 7,12 2,83 0,01 0,95 100 0,5 0,7 0,8 15

7,34 7,10 2,82 0,01 0,95 200 0,5 0,8 0,9 15

7,38 7,14 2,84 0,01 0,95 300 0,5 0,6 0,7 15

14,02 13,22 5,377 0,01 0,95 100 0,5 0,7 0,8 30°

14,006 13,21 5,358 0,01 0,95 200 0,5 0,8 0,9 30°

14,04 13,24 5,396 0,01 0,95 300 0,5 0,6 0,7 30°

с о тз

"¡3

-м

ХЛ

<и &

£ Л

Входные параметры для расчета несущей способности и силы трения

Таблица 2

Д1 А 2 A3 N N1 ю П=П1 в k

0,003 0,1997 3,99 0,95 100 0,5 0,01 0,7 0,8

0,0033 0,1998 3,995 0,95 200 0,5 0,01 0,8 0,9

0,0033 0,1999 3,997 0,95 300 0,5 0,01 0,6 0,7

П

г0 Р

0.7-

0.6-

Рис. 2. Зависимость безразмерной несущей способности от параметра связи N и параметра ю

Я

г0 Р

0.65-

0.60-

0.35

0.50-

0.45-

100 то 1 ' ' "ТТт- -

Пи 200 250 0'Я 0.6 0.4 0.2 и

N1 о

Рис. 3. Зависимость безразмерной несущей способности от параметра N1 и параметра ю

Я

г0 Р 0,5

гит 0 4 0>2 0,8 { 0,8 °'6 •

К со

Рис. 4. Зависимость безразмерной несущей способности от теплового параметра К и параметра ю

(и 5 I

и «

и т

0

1 5

3

ей

Щ

5 I

(и О

а н о

0

1 5

3

ей

=<

=<

Выводы. Результаты численного анализа, приведенные в табл. 1, 2 и на рис. 2-4 позволяют сделать следующие выводы.

1) Основные рабочие характеристики подшипника существенно зависят от:

— параметра ю , характеризующего адаптированный профиль опорной поверхности;

— теплового параметра К;

— экспериментального параметра р , обусловленного зависимостью вязкости от температуры Т.

2) С увеличением теплового параметра К несущая способность подшипника уменьшается. Наиболее резкое уменьшение отмечается при к<0,5. При к>0,7 несущая способность стабилизируется. С увеличением параметра N и параметра ю , характеризующего адаптированный профиль опорной поверхности, несущая способность увеличивается. Параметр связи N не оказывает существенного влияния на несущую способность. При значении параметра ю =0,5 наблюдается ярко выраженный максимум несущей способности. При N ^ ж значение несущей способности стремится к соответствующему значению несущей способности для случая ньютоновского смазочного материала.

Библиографический список

1. Мукутадзе, М. А. Математическая модель сжимаемой микрополярной гидродинамической смазки радиального подшипника с адаптированным профилем его опорной поверхности / М. А. Мукутадзе // Вестник Дон. гос. техн. ун-та. — 2011. — Т. 11, № 8 (59). — С. 1400-1404.

2. Математическая модель гидродинамической смазки радиального подшипника, работающего в нестационарном режиме на микрополярной смазке / К. С. Ахвердиев [и др.] // Вестник Ростов. гос. ун-та путей сообщения. — 2008. — № 1 (29). — С. 147-151.

3. Гидродинамический расчет радиального подшипника, работающего в нестационарном режиме на вязкопла-стичной смазке, обладающей микрополярными свойствами / К. С. Ахвердиев [и др.] // Вестник Ростов. гос. ун-та путей сообщения. — 2008. — № 4 (32). — С. 131-138.

4. Аналитический метод прогнозирования значений критериев микрополярной смазки, обеспечивающих устойчивый режим работы радиального подшипника скольжения / К. С. Ахвердиев [и др.] // Трение и износ. — 2008.

— Т. 29, № 2. — С. 184-191.

5. Математическая модель микрополярной смазки подшипников скольжения с податливой опорной поверхностью / К. С. Ахвердиев [и др.] // Трение и смазка в машинах и механизмах. — 2012. — № 6. — С. 22-25.

6. Вовк, А. Ю. Математическая модель прогнозирования значений безразмерных критериев микрополярной смазки, обеспечивающих рациональный режим работы радиального подшипника скольжения / А. Ю. Вовк, М. А. Мукутадзе, М. А. Савенкова // Вестник Ростов. гос. ун-та путей сообщения. — 2007. — № 1 (25). — С. 5-8.

7. Расчетная модель радиального подшипника скольжения с повышенной несущей способностью, работающего на микрополярной смазке с учетом ее вязкостных характеристик от давления [Электронный ресурс] / К. С. Ахвердиев [и др.] // Инженерный вестник Дона. — 2013. — № 4. — Режим доступа: http://ivdon.ru/magazine/archive/n4y2013/2200 (дата обращения: 20.11.14).

8. Расчетная модель упорного подшипника скольжения с повышенной несущей способностью, работающего на неньютоновских смазочных материалах с адаптированной опорной поверхностью [Электронный ресурс] / К. С. Ахвердиев [и др.] // Инженерный вестник Дона. — 2013. — № 4. — Режим доступа: http://ivdon.ru/magazine/archive/n4y2013/2201 (дата обращения: 20.11.14).

9. Ахвердиев, К. С. Разработка расчетных моделей подшипников скольжения на основе усовершенствования упругодинамической теории смазки / К. С. Ахвердиев, М. А. Мукутадзе, А. Ч. Эркенов. — Ростов-на-Дону : Изд-во РГУПС, 2012. — 371 с.

10. Мукутадзе, М. А. Разработка системы расчетных моделей подшипников скольжения на основе развития гидродинамической и реодинамической теории смазки : дис. ... д-ра техн. наук / М. А. Мукутадзе. — Ростов-на-Дону : РГУПС, 2015. — 476 с.

§ 11. Лагунова, Е. О. Расчетная модель с учетом зависимости вязкости от давления радиального подшипника,

обладающего повышенной несущей способностью / Е. О. Лагунова, Т. С. Черкасова, М. А. Мукутадзе // Тр. Между-'й нар. науч.-практ. конф. «Транспорт-2014». Ч. III. Технические и естественные науки. — Санкт-Перербург, 2014. —

£ С. 218-220.

>

-тз

а £ .й

References

1. Mukutadze, M.A. Matematicheskaya model' szhimaemoy mikropolyarnoy gidrodinamicheskoy smazki radial'nogo podshipnika s adaptirovannym profilem ego opornoy poverkhnosti. [Mathematical model of compressible micropolar hydro-dynamic lubrication for radial bearing with adapted profile of its surface.] Vestnik of DSTU, 2011, vol. 11, no. 8 (59), pp. 1400-1404 (in Russian).

2. Akhverdiev, K.S., et al. Matematicheskaya model' gidrodinamicheskoy smazki radial'nogo podshipnika, rabotay-ushchego v nestatsionarnom rezhime na mikropolyarnoy smazke. [Mathematical model of hydrodynamic lubrication for radial bearing operating in an unsteady mode on micropolar lubrication.] Vestnik RGUPS, 2008, no. 1 (29), pp. 147-151 (in Russian).

3. Akhverdiev, K.S., et al. Gidrodinamicheskiy raschet radial'nogo podshipnika, rabotayushchego v nestatsionarnom rezhime na vyazkoplastichnoy smazke, obladayushchey mikropolyarnymi svoystvami. [The hydrodynamic calculation of radial bearing operating in nonstationary conditions with viscous-plastic lubricant having micropolar properties.] Vestnik RGUPS, 2008, no. 4 (32), pp. 131-138 (in Russian).

4. Akhverdiev, K.S., et al. Analiticheskiy metod prognozirovaniya znacheniy kriteriev mikropolyarnoy smazki, o-bespechivayushchikh ustoychivyy rezhim raboty radial'nogo podshipnika skol'zheniya. [Analytical method for prognosis of values of micropolar lubrication criteria providing stable operation of radial sliding bearing.] Friction and Wear, 2008, vol. 29, no. 2, pp. 184-191 (in Russian).

5. Akhverdiev, K.S., et al. Matematicheskaya model' mikropolyarnoy smazki podshipnikov skol'zheniya s podatlivoy opornoy poverkhnost'yu. [Mathematical model of micropolar greasing of bearings of sliding with a pliable basic surface.] Friction & Lubrication in Machines and Mechanisms, 2012, no. 6, pp. 22-25 (in Russian).

6. Vovk, А.Y., Mukutadze, М.А., Savenkova, M.A. Matematicheskaya model' prognozirovaniya znacheniy bez-razmernykh kriteriev mikropolyarnoy smazki, obespechivayushchikh ratsional'nyy rezhim raboty radial'nogo podshipnika skol'zheniya. [Mathematical model for predicting values of dimensionless criteria of micropolar lubrication ensuring rational operation of the radial plain bearing.] Vestnik RGUPS, 2007, no. 1 (25), pp. 5-8 (in Russian).

7. Akhverdiev, K.S., et al. Raschetnaya model' radial'nogo podshipnika skol'zheniya s povyshennoy nesushchey spo-sobnost'yu, rabotayushchego na mikropolyarnoy smazke s uchetom ee vyazkostnykh kharakteristik ot davleniya. ^omputa-tional model of the radial plain bearing with high load carrying capacity operating at micropolar grease with account of its viscous characteristics from pressure.] Engineering Journal of Don, 2013, no. 4. Available at: http://ivdon.ru/magazine/archive/n4y2013/2200 (accessed: 20.11.14) (in Russian).

8. Akhverdiev, K.S., et al. Raschetnaya model' upornogo podshipnika skol'zheniya s povyshennoy nesushchey spo-sobnost'yu, rabotayushchego na nen'yutonovskikh smazochnykh materialakh s adaptirovannoy opornoy poverkhnost'yu. [Computational model of thrust sliding bearing with high load-bearing capacity operating on non-Newtonian lubricants with adapted supporting surface.] Engineering Journal of Don, 2013, no. 4. Available at: http://ivdon.ru/magazine/archive/n4y2013/2201 (accessed: 20.11.14) (in Russian).

9. Akhverdiev, K.S., Mukutadze, M.A., Erkenov, A.C. Razrabotka raschetnykh modeley podshipnikov skol'zheniya na osnove usovershenstvovaniya uprugodinamicheskoy teorii smazki. [Development of computational models of plain bearings on the basis of improving elastic dynamical theory of lubrication.] Rostov-on-Don: RSTU Press, 2012, 371 p. (in Russian).

10. Mukutadze, M.A. Razrabotka sistemy raschetnykh modeley podshipnikov skol'zheniya na osnove razvitiya gidrodinamicheskoy i reodinamicheskoy teorii smazki: dis. ... d-ra tekhn. nauk. [Development of calculation models system of sliding bearings on the basis of the hydrodynamic lubrication theory and rheodynamic theory of lubrication: Dr.Sci. (Engs.) diss.] Rostov-on-Don: RSTU, 2015, 476 p. (in Russian).

11. Lagunova, Е.О.,Cherkasova, T.S., Mukutadze, M.A. Raschetnaya model' s uchetom zavisimosti vyazkosti ot

davleniya radial'nogo podshipnika, obladayushchego povyshennoy nesushchey sposobnost'yu. [Computational model with ^

account of dependence of viscosity on the radial bearing pressure which has higher load-bearing capacity.] Tr. Mezhdunar. ¡^

nauch.-prakt. konf. «Transport-2014». Ch. III. Tekhnicheskie i estestvennye nauki. [Proc. Int. Sci.-Pract. Conf. "Transport- m

2014". Part III. Engineering and Science] St. Petersburg, 2014, pp. 218-220 (in Russian). g

s

a

<u

<u

Поступила в редакцию 24.06.2016 Received 24.06.2016

Сдана в редакцию 24.06.2016 Submitted 24.06.2016 g

Ж

Запланирована в номер 30.09.2016 Scheduled in the issue 30.09.2016 g

о, н о о к к

3

й