НЕСТАЦИОНАРНЫЕ КОЛЕБАНИЯ РОЛИКА, КОНТАКТИРУЮЩЕГО С ТВЕРДОЙ ПОВЕРХНОСТЬЮ, ПРИ НАЛИЧИИ СМАЗОЧНОГО СЛОЯ

Полный текст

Введение. Стационарная гидродинамическая задача контакта ролика с пластиной при наличии смазочного слоя, разделяющего поверхности, рассматривалась во многих публикациях [1-5] и достаточно хорошо изучена. В то же время нестационарные аспекты гидродинамического контакта ролика с твердой поверхностью остаются в значительной мере неисследованными. Эта тема важна и актуальна, так как именно нестационарный контакт свойственен роликоподшипникам, работающим при переменных нагрузках, которые возникают при стартах космических аппаратов. При такой работе подшипников происходят быстрые изменения зазоров между контактирующими телами, сопровождающиеся значительным ростом вертикальных скоростей, которые, в свою очередь, приводят к резкому росту пиковых значений давления в смазочном слое. В данной статье рассматривается асимптотический аналитический метод [6-8] решения нестационарной гидродинамической задачи с применением коэффициента демпфирования смазочного слоя. Схема контактного взаимодействия цилиндрического ролика с твердой поверхностью, покрытой слоем смазочного материала, представлена на рис. 1. G:\АСПИРАНТУРА\ОСНОВНОЕ\Статьи\вестник сибгау 3\рис 1.tif Рис. 1. Схема расположения ролика и слоя жидкого смазочного материала Расчет давления в смазочном слое. Рассмотрим идеализированную модель контакта ролика с поверхностью при постоянном коэффициенте вязкости в смазочном слое [9; 10]. В этом случае распределение давления определяется из уравнения Рейнольдса [11]: (1) где P - давление в смазочном слое; V - тангенциальная скорость ролика относительно поверхности; m - динамическая вязкость масла при нормальном давлении; h - толщина смазочного слоя, зависящая от деформации поверхностей. Ось x ориентирована вдоль поверхности контакта, как показано на рис. 1. Здесь используется система отсчета, в которой ролик имеет нулевую тангенциальную скорость, а поверхность, соответственно, движется в направлении X. В предположении, что площадка контакта цилиндра и плоскости мала по сравнению с радиусом кривизны R, имеем следующее выражение для толщины слоя смазочного материала [2]: (2) где - минимальная толщина смазочного слоя; - координата точки минимального зазора. Граничные условия в рассматриваемом случае имеют следующий вид [2]: (3) где и - входная и выходная границы смазочного слоя. Для удобства решения задачи вводим безразмерные переменные: (4) Используя (4), преобразуем исходное уравнение Рейнольдса (1) к более простому виду: (5) Положения входной и выходной границ будем характеризовать безразмерными параметрами a и c. Значение параметра a зависит от количества смазки. В случае обильной смазки полагаем [2; 3; 12]. Интегрируя уравнение (5) и используя нулевое граничное условие (3) для производной функции давления при , получаем дифференциальное уравнение первого порядка: (6) Решение уравнения (6) зависит от безразмерного параметра υ, который связан с нормальной скоростью перемещения ролика. Интегрируя уравнение (6) для различных значений υ, находим распределение давления в смазочном слое с учетом нормальной скорости ролика (рис. 2): (7) Определение коэффициента демпфирования смазочного слоя. По найденным распределениям давления определяем значения несущей способности, являющейся функцией параметра υ: (8) Согласно формуле (7), давление линейно зависит от параметра υ, поэтому несущую способность также представим в виде линейной функции следующего вида: (9) где постоянные коэффициенты и A равны 0,401 и 1,125 соответственно. Для перехода к размерному виду используем зависимость (10) С учетом соотношений (4) и (10), преобразуем выражение (9) к размерному виду: (11) Здесь первое слагаемое выражает зависимость стационарной несущей способности от зазора при нулевой нормальной скорости, а второе слагаемое учитывает влияние нормальной компоненты скорости. Коэффициент перед скоростью будем называть коэффициентом демпфирования l: (12) Исследование нестационарного процесса. Для исследования нестационарного процесса запишем уравнение движения ролика по нормали к поверхности: (13) где m - масса ролика; F - внешняя нагрузка. С учетом зависимостей (11) и (12) уравнение (13) принимает следующий вид: (14) Полагая равными нулю производные по времени, определяем равновесное значение зазора: (15) G:\АСПИРАНТУРА\ОСНОВНОЕ\Статьи\вестник сибгау 3\Рис 2.tif Рис. 2. Безразмерное распределение давления в смазочном слое при разных υ: 1 - υ = +0,1; 2 - υ = +0,05; 3 - υ = 0; 4 - υ = -0,05; 5 - υ = -0,1 Принимая h0 в качестве базового зазора, вводим безразмерные переменные: , , (16) Используя нормировки (16), приводим уравнение динамики к безразмерному виду: (17) Уравнение (17) определяет зависимость зазора от времени в процессе установления стационарного режима. Характерное время переходного процесса определяется параметром t0. Зададим начальные условия при t = 0: (18) Так как параметр e при старшей производной является малым, то уравнение называется жестким [13], и его решение может быть приближенно представлено в виде асимптотического разложения по сингулярному малому параметру [14; 15]: где первое слагаемое является регулярной частью асимптотики, а второе слагаемое представляет собой сингулярную часть асимптотики, называемую также в литературе погранслойной асимптотикой. На пограничные функции накладываются дополнительные граничные условия затухания (19) Для первых членов асимптотики имеем выражения (20) (21) Уравнения на функции регулярной и погранслойной частей асимптотики можно получить путем подстановки выражения (14) в уравнение колебаний (17), группировки членов с одинаковыми степенями малого параметра и приравнивая к нулю. В результате получаем следующую систему уравнений: (22) (23) (24) Интегрируя уравнение (22), получаем: (25) Функция (25) удовлетворяет начальному условию h(0), но дает ненулевую начальную скорость: (26) Для компенсации этой скорости далее вводится погранслойное решение. Далее переходим к определению второго члена регулярной части асимптотики. Для этого преобразуем уравнение (23) к виду (27) где (28) Подставляя (28) в (27) и выполняя алгебраические преобразования, получаем: (29) Введем замену: (30) Подставляя (30) в (29), получим следующее выражение: (31) Следовательно, (32) Выражение (31) с учетом (32) примет вид (33) Разделив уравнение (33) на коэффициент при старшей производной, получим дифференциальное уравнение следующего вида: (34) Общим решением уравнения (34) является функция (35) Преобразуем подынтегральное выражение: (36) Тогда уравнение (35) с учетом (36) примет вид (37) где G - неизвестная функция, определяющаяся интегрированием следующего уравнения: (38) Вычисляя интеграл, находим (39) В подынтегральном выражении второго слагаемого уравнения (39) применяем замену переменных и выполняем алгебраические преобразования: (40) Интегрируя (40) и производя обратную замену переменных, получим выражение: (41) С учетом уравнения (41) уравнение (37) примет вид (42) Далее, интегрируя уравнение (24) с учетом условия затухания при t → ¥, получаем выражение (43) Постоянная интегрирования С1 определяется из условия нулевой начальной скорости для полного решения, представляющего собой сумму регулярной и погранслойной асимптотик. Так как первый член регулярной асимптотики дает ненулевое начальное значение скорости (44) то подбираем постоянную С1 таким образом, чтобы погранслойное решение компенсировало значение (26) при t = 0: (45) Подставляя найденное значение в выражение (44) и повторно интегрируя, находим выражение погранслойной части асимптотики: (46) Константа C2 подбирается таким образом, чтобы решение удовлетворяло нулевому граничному условию при t → ¥ (С2 = 0). Найденная погранслойная функция дает ненулевое возмущение в начальный момент времени, которое можно скомпенсировать подбором постоянной интегрирования C0 в выражении второго члена регулярной асимптотики: (47) В результате итоговое выражение принимает вид: (48) После определения всех функций регулярной и погранслойной частей асимптотики уравнения (20) построим график зависимости зазора от времени для различных начальных значений (рис. 3). Рис. 3 показывает, что после любого внезапного скачка нагрузки, который характеризуется безразмерным параметром h′(0), величина зазора стремится к новому равновесному значению. При этом скорость изменения зазора тем больше, чем больше начальное отклонение зазора от равновесного значения. На рис. 4 представлен график изменения относительной скорости сближения поверхностей между контактирующими телами со временем для начального условия h′(0) = 2. Рис. 4 показывает, что увеличение параметра ε приводит к более плавному изменению вертикальной скорости (dh′/dt′). При этом для очень малых значений ε скорость изменения зазора возрастает практически мгновенно, что в свою очередь вызывает резкий и большой скачок давления. Зная зависимость зазора от времени, построим график изменения максимума давления (рис. 5). Из рис. 5 видно, что процесс установления можно охарактеризовать двумя временными интервалами. На первом интервале отмечается очень быстрый рост давления. Амплитуда и длительность роста давления определяется параметром ε. Чем меньше малый параметр ε, тем быстрее увеличивается со временем давление в смазочном слое. На втором временном интервале происходит процесс плавного спада давления к стационарному значению, соответствующему установившемуся состоянию смазочного слоя при постоянной нагрузке. На рис. 6 представлен график нелинейной зависимости пикового значения давления от величины начального зазора, отнесенного к равновесному значению. Данный график показывает, что чем больше малый параметр ε и чем больше отклонение от равновесного зазора, тем сильнее проявляется нелинейность увеличения пикового давления. Величина hm/h0 обратно пропорциональна скачку внешней нагрузки (отношению действующих сил F1/F0). Рис. 5 и 6 показывают, насколько важно учитывать нестационарные переходные процессы в узлах трения летательных аппаратов. Например, при медленном (квазистационарном) увеличении нагрузки в 2 раза максимальное по слою давление увеличивается в 2,2 раза. Однако после внезапного скачка нагрузки в 2 раза во время переходного процесса максимальное по слою давление кратковременно возрастает в 13 раз. Такой резкий скачок давления в смазочном слое между контактирующими поверхностями критически сказывается на ресурсе всего узла трения. G:\АСПИРАНТУРА\ОСНОВНОЕ\Статьи\вестник сибгау 3\рис 3.tif Рис. 3. Зависимость зазора от времени при e = 0,01 для различных начальных условий: 1 - h′(0) = 2; 2 - h′(0) = 1,7; 3 - h′(0) = 1,4 G:\АСПИРАНТУРА\ОСНОВНОЕ\Статьи\вестник сибгау 3\Рис 4.tif Рис. 4. Относительная скорость изменения зазора при начальном условии h′(0) = 2 для различных значений малого параметра: 1 - ε = 0,04; 2 - ε = 0,02; 3 - ε = 0,01; 4 - ε = 0,005 G:\АСПИРАНТУРА\ОСНОВНОЕ\Статьи\вестник сибгау 3\Рис 5.tif Рис. 5. Зависимость максимума давления от времени: 1 - ε = 0,04; 2 - ε = 0,02; 3 - ε = 0,01; 4 - ε = 0,005 G:\АСПИРАНТУРА\ОСНОВНОЕ\Статьи\вестник сибгау 3\Рис 6.tif Рис. 6. Влияние начального условия на изменение пика давления для различных e: 1 - ε = 0,04; 2 - ε = 0,02; 3 - ε = 0,01; 4 - ε = 0,005 Заключение. Построено асимптотическое аналитическое решение задачи нестационарного контактного взаимодействия ролика с твердой поверхностью при наличии смазочного слоя в зоне контакта. Показано, что процесс установления характеризуется двумя временными масштабами. Первый - определяет резкий рост максимума давления сразу после скачка нагрузки. Второй - отражает процесс плавной релаксации давления к стационарному значению, соответствующему возросшему значению нагрузки. Полученные результаты обосновывают важность учета нестационарных переходных процессов в узлах трения летательных аппаратов. Если при медленном (квазистационарном) возрастании нагрузки в 2 раза максимальное по слою давление испытывает примерно двукратное увеличение, то в результате аналогичного по величине, но внезапного скачка нагрузки максимальное по слою давление во время переходного процесса кратковременно возрастает более чем на порядок. Такой значительный и резкий скачок давления в смазочном слое может критически отразиться на ресурсе всего узла трения.

Об авторах

В. А. Иванов

Сибирский федеральный университет, Политехнический институт

Email: VinTextrim@yandex.ru
Российская Федерация, 660074, г. Красноярск, ул. Киренского, 26

Н. В. Еркаев

Институт вычислительного моделирования СО РАН

Российская Федерация, 660036, г. Красноярск, Академгородок, 50/44

Список литературы

Дополнительные файлы