Расчёт скольжения в межвальных роликовых подшипниках ГТД

Автор: Макарчук Владимир Владимирович, Жильников Евгений Петрович

Журнал: Известия Самарского научного центра Российской академии наук @izvestiya-ssc

Рубрика: Научная жизнь

Статья в выпуске: 4-1 т.16, 2014 года.

Бесплатный доступ

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

Подшипники роликовые, скольжение, изнашивание

Короткий адрес: https://sciup.org/148203171

IDR: 148203171

Текст научной статьи Расчёт скольжения в межвальных роликовых подшипниках ГТД

центробежного расширения колец подшипника и деталей подшипникового узла.

Увеличение наружного диаметра центрирующей поверхности сепаратора определяется температурным расширением и центробежным расширения вращающегося сепаратора. Тогда зазор “плавания” сепаратора g пл в рабочих условиях будет равен разности диаметров центрирующих поверхностей.

Экспериментальные исследования показывают существенное влияние формы контакта торцов роликов и направляющих бортиков на проскальзывание в роликовых подшипниках. В настоящей работе рассмотрим стандартный вариант – контакт плоских торцов ролика и направляющего бортика.

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

Геометрические характеристики контактов бомбинированных роликов с беговыми дорожками колец рассмотрены в [4]. В настоящей работе рассмотрим характеристики контактов колец с роликами, образующая которых описывается радиусом.

По рекомендациям [4] расчёт толщин элас-тогидродинамических смазочных слоев в контактах выполняется по формуле

Безразмерные параметры в приведенных формулах определяются соотношениями

P s = M o U v( E np R np ) ; P u = « E np ;

P n = E np R np A lp m ( i ) ; Р т = M o ^^ c .

Здесь Д о = Д о t ) – динамическая вязкость смазки; « = Ор ) – пьезокоэффициент вязкости смазки; Л . = A . t ) – коэффициент теплопроводности смазки; Р - коэффициент зависимости вязкости от температуры; E пр – приведенный модуль упругости в контакте; F m ( i )/А l -нагрузка на единицу длины в среднем сечении линии контакта; t c – температура смазки в зоне контакта; U ^ скорость качения в зоне трения.

Степень жидкостного трения, характеризующая наличие металлического контакта микронеровностей шероховатых поверхностей, определится критерием, предложенным Т. Тальяном:

Здесь R aw – среднее арифметическое отклонение микронеровностей поверхности ролика. В качестве R а 1 принимается среднее арифметическое отклонение микронеровностей поверхности беговой дорожки внутреннего R или наружного R колец соответственно.

С использованием критерия Л относительная площадь металлического контакт микронеровностей вычисляется по формуле [3] П = 0,25 ( А )/ 2

Здесь Ф ( А ) – интеграл вероятности.

Коэффициент трения в контакте с учетом металлического контакта определим суммой f = f мет П + f см ( 1 - П ) .

Коэффициент трения f мет в контактах микронеровностей принимается по рекомендациям [5].

Коэффициент трения, определяемый сдвигом смазочного слоя, рассчитывается по отношению f = W А l/Fm ( i ) .

Здесь W – сила трения, отнесенная к единице длины линии контакта, определяемая по формуле

Выражение для интеграла J имеет вид

J = b m Ф ( х а ) exP ( O P о )V Р о , где Ф ( х а ) - интеграл вероятности, в котором предел интегрирования определяется по формуле х а = 7 2 а Р о

При расчётах для контакта с внутренним кольцом принимаем р о = СТ в ( i ) , для контакта с наружным кольцом - р о = СТ н ( i ) .

В формуле для расчёта W знак “+” принимается для поверхности, имеющей большую величину окружной скорости.

В приведенных расчётах вязкость смазки зависят от температуры. В этой связи расчёт коэффициента трения выполняется итерационным методом.

Температура смазки принимается равной температуре соответствующего кольца. Определяются вязкость смазки, толщина смазочного слоя и коэффициент трения. С использованием формулы Блока определяется приращение температуры в зоне трения:

А t = о,83, f ( F m ( i )/ A l )US ({^m() x

x(4 AnYncnVn + 4 Aw^A-Vw- У

Здесь Л п и A w - коэффициенты теплопроводности, У п и y w - плотность, cn и cw - удельная теплоемкость материалов колец подшипника и роликов соответственно; V n и V w – скорости поверхностей беговых дорожек колец и роликов; U S скорость скольжения.

Уточнённое значение температуры смазки определяется суммированием tcv = tn + 2 A t/3.

Расчёт повторяется до сходимости с заданной точностью принятого и уточнённого значений температуры смазки.

После уточнений температуры смазки и коэффициента трения определяется сила трения скольжения в контакте

,х Ff = fF(i), где F (i) - нагрузка в контакте ролика с внутренним или наружным кольцом.

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

Принимаем, что между каждым из торцов ролика и гнезда сепаратора устанавливается смазочный слой толщиной h« =А 6р

При вращении ролика касательные напряжения в слое смазки определяются по формуле Ньютона ^ = Д о V s6 б , где V – скорость скольжения.

При расчётах вязкость смазки принимаем в зависимости от температуры наружного кольца подшипника.

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

Усилие взаимодействия роликов с перемычкой сепаратора определяется суммой

Fs(i)=+ F;, ±Ff ±2Ff,.

Здесь Ffв , Ffн и Ffб силы трения в контактах ролика с внутренним и наружным кольцами и бортиком. Знак “+” выбирается, если скорость поверхности ролика меньше скорости соответствующей поверхности колец в контакте.

С полученным значением усилия взаимодействия определяется напряжение в контакте и полуширина площадки контакта.

Расчёт толщины эластогидродинамического смазочного слоя и силы трения в контакте выполняется как для контакта ролика с кольцом. Приращение температуры в зоне трения определяем при одной не движущейся относительно источника тепловыделения поверхности.

Момент трения между торцами ролика и гнезда сепаратора определяем, принимая, что между каждым из торцов ролика и гнезда сепаратора в рабочих условиях устанавливается смазочный слой толщиной, равный половине зазора между торцами.

При расчётах вязкость смазки принимаем в зависимости от средней температуры подшипника.

Расчёт контакта сепаратора с центрирующей поверхностью бортика наружного кольца выполним как узкого подшипника скольжения [6].

Усилие прижатия сепаратора к центрирующей поверхности определим по формуле

Здесь F cs – центробежная сила при смещении сепаратора.

Составляющие усилия прижатия роликами по координатным осям определяются суммированием проекций усилий взаимодействия роликов с перемычкой сепаратора:

FSB = L (Fs (i)sin(^)- Ffs (i)cos(^));

i = 1

i = 1

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

Приращение температуры в зоне трения определяем при одной не движущейся относительно источника тепловыделения поверхности.

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

При расчёте кинематики ролика скорости поверхностей колец (рис.1) в обращенном движении (при остановленном сепараторе) определяются по формулам:

V e = d ep ( H - Ю в )l2 – для внутреннего кольца;

Рис. 1. Схема скоростей поверхностей межвального подшипника

V h = D ,p ( ш. - to o )l 2 – для наружного кольца.

Ведущим в межвальном подшипнике является наружное кольцо. В этой связи поверхность ролика будет иметь проскальзывание некоторой величины относительно кольца. Тогда принимаем V w = V . ( 1 £ w ) , где £ w - относительная величина проскальзывания ролика в контакте с наружным кольцом. При этом угловая скорость вращения ролика относительно собственной оси будет равна H w = 2 V w/ D wp .

Скорости скольжения ролика в контактах с кольцами будут равны:

Vse = V w Ve - для внутреннего кольца;

  • Vs . = V . V w - для наружного кольца.

На рис. 2 приведена схема сил, действующих на ролик в нагруженной зоне.

Здесь F S и F fS – нормальное усилие и сила трения в контакте ролика с перемычкой сепаратора; F с – центробежная сила ролика.

Условие равновесия сил в радиальном направлении определяется выражением

  • F . = F e + F fs + F c .

Здесь F в , F н – усилия в контактах ролика с внутренним и наружным кольцами.

В связи с малостью величины силы трения F fS уравнение равновесия в радиальном направлении можно привести к виду

F = F A-F н в c.

Это упрощение позволяет решать задачу о распределении нагрузки по телам качения без учёта сил трения в контактах.

Силы трения F и F fS определяются, как показано ранее. Центробежная сила ролика определяется по формуле:

Fc = mw (D.p — Dwp 'ЩЦ2..

Как показано ранее, все силы трения опреде-

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

Третьим уравнением равновесия является уравнение моментов

S Т E MFfH - Ffe - Ffs)Dwp /2 + + Ff6 (Dwp - (DHp - D61 V2)- Тsc - Ткн -

  • - Т кв - 0

Здесь моменты сопротивления качению ролика определим по приближенным зависимостям Т кн - f кFH и Т кв - fк FB - где fк коэффициент трения качения.

Разделив на радиус ролика, получим

S Т E w ) - F fH - F f6 - F fs +

2 F f ( 1 - ( D h P - D 6 1 )/ 2 D wp )

  • 2 ( Т sc + Т кн + Т кв )/ D wp - 0. '

Полученное уравнение является функцией одной переменной E . Его решение выполняется итерационным методом Ньютона. При этом производную функции заменяем отношением приращений функции и аргумента.

В результате получаем величину угловой скорости вращения ролика to w и усилие Fs - нагрузки в контакте ролика с перемычкой сепаратора.

Рассматривается также кинематика ролика в не нагруженной зоне.

Таким образом, в межвальном подшипнике все ролики являются “ведущими” по отношению к сепаратору. В то же время в подшипнике с вращающимся внутренним и не вращающимся наружным кольцами ролики в разгруженной зоне являются “тормозящими” и “толкаются” сепаратором. Это объясняет установленное эксперимен- тально небольшое по величине скольжение сепаратора в межвальном подшипнике. Однако при этом наблюдается повышенное скольжение в контактах ролика с кольцами, что может привести к повышенному изнашиванию.

Расчёт сопротивления движению сепаратора в воздушно-масляной среде подшипника выполним по рекомендациям [7].

Температуру воздушно-масляной среды примем равной средней температуре подшипника

t6M = (te + tn )/ 2 .

вм в н

Масло в подшипнике не полностью заполняет весь свободный объем. В этой связи расчётное значение плотности воздушно–масляной среды рекомендуется в зависимости от объемной доли смазки О дсм в полости подшипника по формуле

Г ср = У„о Дсм Кол + 0,6 О дс,. ) .

Характер движения сепаратора в воздушно-масляной среде определяется числом Рейнольдса, которое в данном случае вычисляется по формуле

Re - Rc tos (v .

Момент сопротивления движению сепаратора в воздушно – масляной среде рекомендуется определять по формуле:

Здесь: R c и r c – радиусы боковой поверхности сепаратора, характеризующего сопротивление вращению его в масляной среде; B s – ширина наружной поверхности сепаратора; C n – коэффициент, зависящий от числа Рейнольдса.

Как показано выше, усилия взаимодействия роликов с перемычками сепаратора являются “ведущими”, обеспечивающими вращение сепаратора. Величина момента сил взаимодействия роликов с перемычками сепаратора определится по формуле i -1

Кроме того “ведущим” для сепаратора межвального подшипник является момент от силы трения сепаратора в контакте с центрирующими поясками вращающегося наружного кольца. Его величина определяется по формуле:

T s6 - fF Rs D s 1 .

Здесь F RS – усилия прижатия сепаратора и f – коэффициент трения в контакте, определяемые, как показано ранее.

Сопротивление вращению сепаратора межвального подшипника будет только от момента трения T ts от перемешивания воздушно-масля-ной среды в полости подшипника.

Величины всех моментов T S , T и T RS зависят от угловой скорости вращения сепаратора.

Величину угловой скорости сепаратора в свою очередь можно определить по формуле

  • У0 = У00 (1 - £0 ).

Здесь 8 о - коэффициент скольжения а У 00

- теоретическое (при отсутствии скольжения) значение угловой скорости вращения сепарато- ра, величина которого определяется по известной формуле

У00 = (dврЮв + dhp^h V(2(Dp - Dwp )).

Принимаем условие 0 < ( У н . Тогда с использованием выражений для у и 00 получим

£ 0 ^

(ув - Ун )(1 - Dwp /Dpp ) 2У00

Условие равновесия при равномерном вращении сепаратора будет иметь вид:

  • ■*0 ) = T s + T s6 - T rs = 0 .

Получено уравнение являющееся функцией одной переменной £ 0 . Его решение выполняется итерационным методом Ньютона. При этом производную функции заменяем отношением приращений функции и аргумента.

В результате получаем величину угловой скорости вращения сепаратора 0, величину скольжения £ 0 и усилие и температуры в контактах роликов с кольцами и сепаратором.

Исследование проскальзывания выполним для межвального подшипника 55-2672919Р5 изделия “99В”.

При расчетах скольжения сепаратора принята минимальная радиальная нагрузка F r = 50 Н . Перекос колец принят равным нулю. Скорости вращения колец приняты соответствующими режимам работы подшипника в двигателе.

Температура наружного кольца принята равной 43°С, внутреннего – 47,3°С. Монтажный радиальный зазор в подшипнике принят равным g = 0,043 мм. На рис. 3 и 4 приведены зависимости скольжения сепаратора от разности частот вращения наружного и внутреннего колец (скольжения роторов) при постоянной частоте вращения наружного или внутреннего колец.

При этом постоянная частота вращения наружного кольца (рис. 3) принята равной n н =13067 об/мин , а частота вращения внутреннего кольца (рис. 4) соответственно – n в =8152 об/мин .

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

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

Следовательно, уменьшается скольжение (увеличивается отрицательное значение проскальзывания).

Таким образом, зависимость скольжения сепаратора от частоты вращения наружного кольца имеет сложный характер.

На рис. 5 приведена зависимость скольжения сепаратора от одно-временного увеличения частот вращения внутреннего и наружного колец при постоянной разности скоростей.

Рис. 3. Зависимость проскальзывания от скольжения роторов при п = 13067 об/мин н

Рис. 4. Зависимость проскальзыванияот скольжения роторов при n =8152 об/мин в

Результаты расчётов показывают, что скольжение сепаратора зависит как от величины скольжения роторов, но и от величин скоростей вращения колец.

На рис. 6 приведена зависимость скольжения сепаратора от радиальной нагрузки при температуре наружного кольца 33 °С и частотах вращения колец: n в =5180 об/мин и п н = 11313 об/ мин . Результаты расчётов показывают известное уменьшение скольжения сепаратора с увеличением радиальной нагрузки на подшипник.

На рис. 7 приведены результаты испытаний подшипника 55-2672919Р5 при различных значениях радиального зазора при смазке маслом ИПМ – 10 [8]. Величины проскальзывания сепаратора приведены в зависимости от частоты вращения наружного кольца. При этом частоты вращения внутреннего кольца соответствовали режимам работы подшипника в изделии.

Рис. 5. Зависимость проскальзывания от частоты вращения наружного кольца при П н П в = 4300 об/мин

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

Экспериментальные зависимости проскальзывания в подшипнике от разности скоростей вращения колец при нагрузке Fr=50 Н и различных значениях радиального зазора приведены на рис. 8. Видно увеличение проскальзывания с увеличением “скольжения роторов”.

Приведенные результаты экспериментов показали, что на всех режимах испытаний частота вращения сепаратора была выше эпициклической, рассчитанной без учета сил трения в контак-

Рис. 6. Зависимость скольжения сепаратора от радиальной нагрузки

Рис. 7. Зависимость проскальзывания комплекта роликов в подшипнике при F =50 Н от частоты вращения и величины радиального зазора

Рис. 8. Зависимость проскальзывания комплекта роликов от разности частот вращения колец при различных значениях радиального зазора

тах и сопротивления вращению сепаратора.

На рис. 9 приведены результаты расчёта проскальзывания сепаратора на режимах испытаний при радиальном зазоре в сборе 0,043 мм и радиальной нагрузке 50 Н .

При выполнении расчётов температуры смазки принимались по экспериментальным данным. Для сравнения там же приведены экспериментальные значения проскальзывания комплекта роликов. Результаты расчётов имеют удовлетворительное соответствие с экспериментальными в связи с многофакторностью зависимости проскальзывания: от радиальной нагрузки, радиального зазора, как от скоростей вращения колец, так и от “сколь жения” роторов, температуры не только наружного, но и внутреннего колец, температуры смазки, ее свойств, объема и способа подачи в подшипник и другое.

К сожалению, многие факторы сложно оценить и учесть в расчётах.

Вместе с тем результаты экспериментов и теоретических расчётов позволяют сделать некоторые выводы о работоспособности межвальных подшипников.

Прежде всего, установлено, что термин “проскальзывание” для межвальных подшипников является условным, так как скорость вращения сепаратора оказывается не ниже, а выше теоретической (эпициклической). Это объясняется тем, что в межвальном подшипнике большинство роликов (даже в разгруженной зоне) являются “ведущими” так как прижимаются значительными по величине центробежными силами к вращающемуся наружному кольцу. Теоретический анализ показывает, что вследствие этого скольжение роликов относительно наружного кольца мини-

Разность частот вращ.роторов. (об./мин.)

Рис. 9. Зависимость проскальзывания комплекта роликов от разности частот вращения колец при g = 43 мкм и Fr = 50 Н мально и близко к нулю. В то же время наблюдается значительное по величине скольжение роликов относительно внутреннего кольца.

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

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

Список литературы Расчёт скольжения в межвальных роликовых подшипниках ГТД

  • Кузнецов Н.Д. Обеспечение надежности современных авиадвигателей//Проблемы надежности и ресурса в машиностроении. М.: Наука, 1986. С.51-68.
  • Акифьев В.И. Разработка методики расчета роликовых подшипников опор ГТД с учетом проскальзывания: Дисс. … канд. техн. наук. Самара, Самарск. гос. аэрокосм. ун-т, 1998. 165 с., ил.
  • Балякин В.Б., Жильников Е.П., Самсонов В.Н., Макарчук В.В. Теория и проектирование опор роторов авиационных ГТД. Самара: Изд-во Самар. гос. аэрокосм. ун-та, 2007. -254 с.: ил.
  • Макарчук В.В. Разработка методов расчета и проектирования высокоскоростных межвальных роликовых подшипников: Диссертация на соискание ученой степени канд. техн. наук. Самара, Самарск. гос. аэрокосм. ун-т, 2009. 165 с., ил.
  • Крагельский И.В., Добычин М.Н., Комбалов В.С. Основы расчетов на трение и износ. М: Машиностроение, 1997. 526 с.: ил.
  • Жильников Е.П., Самсонов В.Н. Трение и изнашивание в узлах авиационной техники: учеб. пособие. Самара: Изд-во Самар. гос. аэрокосм. ун-та, 2007. 144 с.: ил.
  • Силаев Б.М. Трибология деталей машин в маловязких смазочных средах: монография.-Самара: Изд-во Самар. гос. аэрокосм. ун-та, 2008. -264 с.
  • Исследование проскальзывания межвального подшипника 55 -2672919Р5 изделия 99В: Технический отчет/Руковод.: Г.М. Косинов, Н.И. Петров/ЦИАМ -Москва, 2002г. -18 с.
Еще
Статья научная