УДК 532.135; 621.822.114
МЕТОДИКА РАСЧЁТА ГИДРОМЕХАНИЧЕСКИХ ХАРАКТЕРИСТИК СЛОЖНОНАГРУЖЕННЫХ ПОДШИПНИКОВ СКОЛЬЖЕНИЯ ПОРШНЕВЫХ И РОТОРНЫХ МАШИН, СМАЗЫВАЕМЫХ НЕНЬЮТОНОВСКИМИ МАСЛАМИ
И.Г. Леванов
THE DESIGN PROCEDURE OF CHARACTERISTIC HYDROMECHANICAL OF COMPLEX-LOADED JOURNAL BEARINGS OF RECIPROCATOR AND ROTOR-TYPE MACHINES LUBRICATED BY NON-NEWTONIAN OILS
I.G. Levanov
Представлены основные положения методики расчёта, позволяющей учитывать влияние неньютоновских свойств моторных масел на работу сложнонагруженных подшипников скольжения; а также область её возможного применения и результаты оценки влияния указанных свойств на работу подшипников скольжения коленчатого вала двигателя 4ЧН 13/15.
Ключевые слова: гидромеханические характеристики, неньютоновские масла, подшипники скольжения.
The main positions of calculation methods which to take account of the influence of non-Newtonian properties of engine oils to work complex-loaded journal bearings, as well as the range of its possible application and an assessment of the influence of these properties to work journal bearings crankshaft engine 4ЧН 13/15 are presented in article.
Keywords: characteristic hydromechanical, non- Newtonian oils, journal bearings.
Введение. Решение проблемы повышения надежности таких машин массового применения, как двигатели внутреннего сгорания (ДВС), турбокомпрессоры, неразрывно связано с совершенствованием конструкции основных гидродинамических трибосопряжений (ТС): динамически нагруженных шатунных и коренных подшипников скольжения (ПС) коленчатого вала, ротора турбокомпрессора.
Разработка более совершенных методик расчета динамически нагруженных (сложнонагруженных) ПС остается актуальной задачей, особенно на современном этапе развития инженерной деятельности, для которого характерно широкое применение электронно-вычислительной техники и специального программного обеспечения, позволяющих повышать качество проектирования и надёжность подшипников скольжения и, в конечном итоге, увеличивать ресурс машины в целом.
Значительный вклад в разработку методик расчёта подшипников скольжения внесли многие отечественные и зарубежные исследователи: А.К. Дьячков, С.М. Захаров, В.Г. Караваев, М.В. Ко-ровчинский, В.Н. Прокопьев, О.И. Рабецкая, Ю.В. Рождественский, Л.А. Савин, В.И. Суркин, И.А. Тодер, Д.И. Фёдоров, Н.Н. Типей, T.W. Bates, J. F. Booker, P.K. Goenka, B.A. Gecim, S.D. Gul-wadi, D.R. Chen, R.S. Paranjpe, H.K. Hirani и др.
Тем не менее известные методики не полностью отражают физические процессы, происходящие в системе «шип - смазочный слой - подшипник», поскольку основаны на том допущении классической гидродинамики, что смазочный материал является ньютоновской жидкостью, вязкость которой зависит только от температуры и давления.
Известно, что все товарные марки смазочных масел получают путём загущения маловязкой основы различными добавками (присадками). Количество этих добавок определяется уровнем эксплуатационных свойств, предъявляемых к смазочному маслу. Например, главной целью
введения вязкостных присадок в современные моторные масла является получение улучшенной (пологой) вязкостно-температурной характеристики. То есть при низких температурах вязкость должна быть не слишком высокой, чтобы обеспечить прокачиваемость по системе смазки, доступ к узлам трения и минимальное сопротивление при проворачивании, а при высоких - достаточной, чтобы обеспечить смазочный слой способный нести нагрузку.
Сегодня широкое распространение при эксплуатации ДВС получили именно всесезонные моторные масла, загущенные полимерными присадками. Реологическое поведение таких масел имеет особенности, называемые в литературе неньютоновскими свойствами. К наиболее известным неньютоновским свойствам загущенных масел относятся: зависимость их вязкости от скорости сдвига (псевдопластичность, временное снижение или аномалия вязкости), вязкоупругие эффекты (релаксация касательных напряжений, появление нормальных напряжений при сдвиге).
Благодаря зависимости вязкости от скорости сдвига, такие масла также называют «энергосберегающими», поскольку они позволяют снизить потери мощности на трение в двигателе, а следовательно, и расход топлива.
А.Ю. Вовк, В.Н. Прокопьев, Л.А. Савин, B.A. Gecim, H.G. Elrod, A.V. Harnoy, H.K. Hirani, R.S. Paranjpe и другие авторы предпринимали попытки учесть неньютоновское поведение смазки в методиках расчёта подшипников скольжения. Теоретические и экспериментальные исследования указывают на то, что неньютоновские свойства масел оказывают влияние на характеристики ПС. В частности, установлено, что псевдопластичность смазки приводит к снижению минимальной толщины смазочного слоя, потерь мощности на трение; вязкоупругие эффекты способствуют увеличению толщины смазочного слоя [І].
Сегодня, очевидно, что особенности реологического поведения смазочных материалов необходимо учитывать на этапе проектирования ТС. В частности, эффективное применение всесе-зонных моторных масел в ДВС возможно только в том случае, если это допускает конструкция двигателя и его ТС (подшипников скольжения коленчатого вала).
Однако использование теоретических положений и методик, учитывающих неньютоновские свойства смазки, в инженерной практике проектирования сложнонагруженных гидродинамических ПС ограничено. Прежде всего, это связано с недостаточным объёмом информации о реологическом поведении современных смазочных масел при высоких температурах и скоростях сдвига в ПС; остаётся открытым вопрос о совместном влиянии неньютоновских свойств масел на характеристики сложнонагруженных гидродинамических ПС.
Таким образом, актуальность работы обусловлена, с одной стороны, широким применением для смазки поршневых и роторных машин масел, обладающих неньютоновскими свойствами; с другой стороны, отсутствием методики расчёта сложнонагруженных гидродинамических ПС поршневых и роторных машин, позволяющей оценивать влияние неньютоновских свойств смазочных жидкостей на характеристики подшипников.
В настоящей статье представлены основные положения разработанной методики и результаты её применения к расчёту гидромеханических характеристик подшипников скольжения коленчатого вала дизельного двигателя 4ЧН І3/І5.
Основные положения методики. При разработке методики приняты следующие допущения:
1. Трение в ПС обусловлено внутренним трением смазочной жидкости.
2. Смазочная жидкость «прилипает» к поверхностям шипа и подшипника.
3. Пространство между шипом и подшипником полностью заполнено смазочной жидкостью.
4. Течение смазочной жидкости ламинарное, силы инерции не учитываются.
5. Поверхности шипа и подшипника считаются идеально круглыми, гладкими и жёсткими.
6. Перекосы осей шипа и подшипника не учитываются.
7. Теплоотдача от смазочной жидкости в шип и подшипник не учитывается.
Задача расчёта динамики сложнонагруженного подшипника скольжения сводится к построению траектории движения центра масс шипа (например, шатунной шейки) под действием внешней периодической нагрузки в системе координат XOY (рис. І), закрепленной на втулке (вкладыше), и определению его гидромеханических характеристик (ГМХ), к которым относятся: мгновенные значения минимальной толщины смазочного слоя hmin (t) и максимального гидро-
динамического давления pmax (t) в момент времени t, а также их экстремальные inf hmin, sup pmax и средние h*mn, pmax за цикл нагружения t4 величины, эффективная температура смазочного слоя Тэ, мгновенные и средние потери мощности на трение N (t) , N ; расходы смазки
QB (t) , QB .
Траектория движения центра шипа строится по координатам, получаемым в результате решения уравнения равновесия:
1
h0 «0
L Ru (U ,U) + Fu (t) = 0.
(1)
Здесь и = и / к0 = (X, У) - безразмерный вектор перемещений центра шипа; и, и - производные по безразмерному времени ґ ; Р = (Рх, Ру ) - проекции на соответствующие оси внешней нагрузки; ^ (и ,и) - реакции смазочного слоя; И0 - радиальный зазор; характерная угловая скорость вращения шипа ю0 = &а / &.
Безразмерные реакции смазочного слоя на шип определяются формулой
RU =
R
R
Y1
=_ 4a іp
А
cos ф sin ф
(2)
Здесь р - безразмерные гидродинамические давления в смазочном слое; ф,- угловая и осевая координаты подшипника; 0.А - активная область смазочного слоя, где р > 0.
Эффект релаксации касательных напряжений в смазочном слое учитывается на основе подхода, описанного в работах Рагагуре [2]. Этот подход заключается в замене внешней нагрузки, действующей на подшипник, модифицированной. То есть в уравнение (1) подставляется выражение
ри = ¥и + ^®0 д1ги/5а, (3)
где X - время релаксации смазки.
Интегрирование уравнения движения (1) осуществляется методом, основанным на формулах дифференцирования назад для уравнений второго порядка и описанным в работах В.Н. Прокопьева, Ю.В. Рождественского, К.В. Гаврилова и др. [3].
Поле гидродинамических давлений, необходимое для вычисления реакции смазочного слоя Яц , определяется интегрированием уравнения Рейнольдса:
_д_
дф
h3 ( XL ^ ІФ2-Ф- -д p д h ( ф2 Т2 ^ _4l -д p д «21 p h І1 -І'
1 ф0 у дф д z V ф0 у дф ф0 у
(4)
Здесь р , p = (p - pa)у2/Ц0юо , V = h)lr , Z = z/r , -a < z < a, х = ю0? - безразмерные: плотность смазки, относительный зазор, координата по ширине подшипника, время; a = B/D - относительная ширина подшипника; ц0 - характерная вязкость смазки; pa - атмосферное давление; B,D = 2r, r - ширина, диаметр, радиус подшипника; Ю21 = (ю2 - Ю1 )/юо - безразмерная угловая
скорость вращения шипа; безразмерная толщина смазочного слоя h и её производная д h/ Эх определяются формулами:
h =1 -%cos (ф-8), dh / дх = — cos (ф-8)-х§ cos (ф-8), (5)
где х - относительный эксцентриситет; 5 - угол относительного положения линии центров;
У -к
фк = |Ц-'
0 Ц
(6)
где к = 0; 1; 2 ; Ц* - вязкость смазки, являющаяся функцией скорости сдвига, температуры и давления; у - безразмерная координата поперёк смазочного слоя.
Уравнение (4) интегрируется многосеточным методом, описанным в работах В.Н. Прокопьева, Е.А. Задорожной, А.К. Бояршиновой и др. [4], при граничных условиях Свифта-Штибера с учётом наличия на поверхности шипа источников подачи смазки (отверстий, канавок):
р(ф, 1 = ±а) = 0; р (ф, 7) = р(ф + 2л, 7); р(ф, 7) > 0; р(ф, 1) = р,
на (ф, 1) еО,, 5 = 1,2,..., 5 *. (7)
Здесь О, - область источника смазки, в которой давление постоянно и равно давлению подачи
*
р,; 5 - количество источников смазки.
Зависимость вязкости смазки от скорости сдвига аппроксимируется степенным законом Оствальда-Вейла [1]; от температуры - трёхконстантной формулой Фогеля [5]; от давления -формулой Баруса [6].
Таким образом, математическая модель вязкости смазочного материала имеет вид
Ц* = I (п~1)/2 . С • е^С2/(Тэ +С3 ))+Р(Тэ )'р (8)
Здесь Ц - параметр консистенции (вязкость при низкой скорости сдвига до 102 с-1);
12 =
8у
8^
8У
второй инвариант скоростей сдвига; п - параметр, характеризующий сте-
пень неньютоновского поведения; Р(Гэ) - пьезокоэффициент вязкости смазки, являющийся функцией температуры; р - гидродинамическое давление; Тэ - эффективная температура сма-
зочного слоя; С^, С2, Съ - константы.
Методика определения параметров Ц и п, а также результаты экспериментальных исследований реологических свойств моторных масел представлены в работе [7].
Градиенты скоростей сдвига определяются выражениями:
8^
8У
1
21 + кп+1 'о
У
др
8ф
8У„
1
• к
П+1
Л
8р 8Ї'
(9)
8У Ц . .
Для описания реологического поведения загущенных моторных масел предложено использовать следующую модель (рис. 2, где у = ^/12 ). На участке 1 масло ведёт себя как ньютоновская жидкость с вязкостью ц (Тэ, р). Для участка 2 характерно снижение вязкости по степенному закону (8). На участке 3 масло считается ньютоновской жидкостью с вязкостью ц 2 (Тэ, р).
В такой постановке реологическая модель не противоречит представлениям о поведении загущенного моторного масла, не вызывает трудностей при программной реализации.
Для оценки теплового состояния под-
шипника скольжения используется изотермический подход, основанный на уравнении теплового баланса:
(10)
отражающего равенство средних за цикл
*
значений теплоты ЛК, рассеянной в смазочном слое, и теплоты Л**, отведенной
смазкой, вытекающей в торцы подшипника скольжения.
Ц , Па-с Ц1 ц 2
1 2 /3
' ■ л— і і і 1 ^
1
102
106
у, с
-1
Рис. 2. Зависимость вязкости от скорости сдвига
•!<
Приращение температуры в смазочном слое на каждом временном шаге расчёта траектории центра шипа (за цикл нагружения) определяется по формуле N
ЛТЭ =---------- (11)
в -р-со
где N, вВ, р, Со - потери мощности на трение в смазочном слое, расход смазки через подшипник, плотность и теплоёмкость смазки, соответственно.
Использование такого подхода при решении задач динамики и смазки сложнонагруженных подшипников скольжения обосновано в работах В.Н. Прокопьева, В.Г. Караваева и др. [8, 9].
Потери мощности на трение определяются интегрированием диссипативной функции рассеивания Д по объему смазочного слоя:
N = .Ц д^ф ^ (12)
8а
оа
где О- = ОА - О5;
Д = 1
8У,
Су. (13)
ду ) [ду
Расход смазки в оба торца подшипника определяется выражением ББ 2ую0 2а
—I- Сф. (14)
8,|,=а
2
[
Совокупность выражений (1) - (14) представляет математическую модель смазочного слоя, учитывающую зависимость вязкости от скорости сдвига и эффект релаксации касательных напряжений в смазочном слое.
Алгоритм расчёта. Алгоритм расчёта гидромеханических характеристик подшипников скольжения реализован в комплексе программ:
1. Перед началом вычислений исследователь определяет тип объекта проектирования, метод расчёта и задаёт исходные данные:
- геометрические характеристики подшипника (диаметр, ширину, радиальный зазор, радиус кривошипа);
- параметры смазочного материала (коэффициенты динамической вязкости при характерных температурах, показатель степенного закона и параметр консистенции как функции температуры п (Т) , ц (Т), время релаксации);
- типы и расположение источников смазки (канавки, отверстия), а также давление подачи.
Начальные (стартовые) значения:
- температуры смазки на входе в подшипник (предполагаемое значение эффективной температуры смазочного слоя Тэ ); распределения давления в смазочном слое р0 (ф, 1); координат центра масс шипа; начальных скоростей смазки и0(ф, 1);
- нагрузочный режим: угловой скорости вращения вала, внешней нагрузки и её распределения во времени.
Для решения прямой гидродинамической задачи - расчета траектории движения центра шипа, т. е. определение координат под действием внешней периодической нагрузки, весь процесс дискретизируется по времени.
2. Цикл итерационного расчета траектории равен периоду нагружения То (для четырёхтактных ДВС составляет 720° поворота коленчатого вала). Интервал периода То разбивается на шаги: Атг- = тг-+1 — тг-; затем ищется приращение угла поворота шипа Афг- = ю • Атг-, где ю - угловая скорость шипа.
3. На каждом временном шаге цикла нагружения подшипника определяется функция зазора, интегрируется уравнение Рейнольдса, вычисляется поле гидродинамических давлений (ГДД) в тонком смазочном слое. При этом происходит обновление массивов приближения гидродинамических давлений и градиентов скоростей. По найденным значениям давлений определяем реак-
ции, действующие на шип со стороны смазочного слоя. В следующий момент времени (хг- + Ах) значения параметров, полученные на предыдущем временном шаге хг-, используются как начальные значения для периода времени хг-+!.
4. Рассчитываются мгновенные значения гидромеханических характеристик подшипника: расход смазки через торцы подшипника, диссипативная функция рассеяния энергии, потери мощности на трение и др.
5. Интегрированием уравнения движения определяются координаты и скорости движения центров подвижных элементов трибосопряжения.
6. На каждом временном шаге расчета в результате решения уравнения теплового баланса уточняется значение эффективной температуры Тэ , корректируется коэффициент динамической вязкости.
7. После завершения счета всего цикла нагружения х0 вычисляются среднеинтегральные характеристики гидромеханических характеристик за период Хо , а также их экстремальные значения.
8. На каждом шаге расчета и по завершении всего цикла нагружения Хо результаты записываются в систему файлов. По результатам расчета строятся графические зависимости в удобной для восприятия форме.
При этом необходимо отметить, что расчёт ведется в два этапа. Первый этап заключается в расчете полного цикла нагружения (пп. 2-8) без учёта особенностей реологического поведения смазочного материала (без учёта неньютоновских свойств). На втором этапе расчет повторяется с учетом реологических свойств смазочного материала, которые задаются в помощью набора ключей в исходных данных.
Главные различия приведённого алгоритма от традиционной схемы решения заключается в следующем
1. Внешняя нагрузка, действующая на подшипник, до начала расчёта заменяется на модифицированную с учётом времени релаксации смазки. Тем самым учитываются эффекты релаксации касательных напряжений в смазочном слое.
Вход в подпрограмму расчёта ГДД
Вычисление зазоров.
Первое приближение по давлениям и вязкости
А 'і
^1
Нет
Расчёт градиентов скоростей сдвига, неньютоновской вязкости
Решение уравнения Рейнольдса
I
Корректир овка вязкости от давления
■и
Уравнение Рейнольдса после его разностной аппроксимации сводится к системе алгебраических уравнений, которая с учетом граничных условий Свифта-Штибера решается итерационным методом Зейделя
Расчёт потерь мощности на трение и расхода смазки
Выход в подпрогр амму расчёта реакций смазочного слоя
Рис. 3. Схема вычисления гидродинамических давлений смазочного слоя с учётом зависимости вязкости от скорости сдвига
2. Процедура вычисления градиентов скоростей сдвига предшествует решению уравнения Рейнольдса (рис. 3). Это обусловлено тем, что вязкость смазочного материала является функцией скоростей сдвига.
Решение тестовой задачи. В качестве тестовой была выбрана задача расчёта гидромеханических характеристик третьего коренного подшипника коленчатого вала двигателя V6 производства компании General Motors. Исходные данные для расчёта представлены в работе [10]. Авторы этой работы рассчитали гидромеханические характеристики указанного подшипника с помощью программного пакета FLARE, а также экспериментально ёмкостным методом определили толщину смазочного слоя подшипника. Сравнение результатов тестовых расчётов с экспериментальными и теоретическими данными американских исследователей Paranjpe, Tseregounis и Viola представлено в табл. 1 и на рис. 4 и показывает хорошее качественное и количественное соответствие.
Таблица 1
Сравнение результатов тестовых расчётов
Источник Характеристика
inf hmin, мкм аh ■ , град "min Средняя нагрузка, Н Т, °С Скорость сдвига1 х106, с-1
Эксперимент [10] 1,1б 29б,0 - 104-10б -
Расчёт FLARE [10]2 1,25 285,0 4281 102,3 1,75
Расчёт автора3 1,43 294,0 4331 102,9 1,52
Расчёт автора4 1,81 29б,1 4327 10б,1 1,26
Расчёт автора5 2,12 291,4 4340 105,9 1,17
Расчёт автора6 1,82 290,1 4350 105,8 1,15
Расчёт автора7 1,44 288,2 4347 102,8 1,34
Расчёт автора8 1,11 289,1 4354 102,1 1,52
Примечания: 1средние значения за цикл нагружения; 2Paranjpe, Tseregounis, Viola, программный комплекс (США); 3с учётом зависимости вязкости от скорости сдвига; 4без учёта неньютоновских свойств масла; 5с учётом зависимости вязкости от давления; 6с учётом вязкоупругих свойств масла при X = 5 10~4 с ; 7с учётом зависимости вязкости от скорости сдвига и вязкоупругих свойств масла; 8с учётом зависимости вязкости от давления, скорости сдвига, и вязкоупругих свойств масла. Во всех расчётах учитывалась зависимость вязкости от температуры.
Из рис. 4 видно, что расчётная зависимость 3 минимальной толщины смазочного слоя от угла поворота коленчатого вала, полученная автором, хорошо согласуется с зависимостью 2 Paranjpe и др.
Расчётиая зависимость 4, полученная без учёта неньютоновских свойств масла, располагается выше 3, при этом значение inf hmin завышено (см. табл. 1) на 21 % относительно 3.
Учёт зависимости вязкости масла от давления приводит к повышению inf hmin на 14,6 % относительно 4; учёт вязкоупругих свойств масла при X = 5 -10_4 с-1 приводит к незначительному, около 0,5 % увеличению inf hmin по отношению к 4; совместный учёт зависимости вязкости от давления, скорости сдвига и влияния вязкоупругих свойств смазки приводит к снижению inf hmin на 16,5 % относительно 5. Отличие по величине температуры смазочного слоя не превышает 3 % относительно 3.
h ■
min’
мкм
10
, | Р езультаты эк :сперимента [10] ■ 1
4 /
3 2
0 180 360 540 а, град
Рис. 4. Минимальная толщина смазочного слоя в коренном подшипнике в зависимости от угла поворота коленчатого вала двигателя V6 GM
5
0
Описание программного комплекса. Разработанная методика реализована в комплексе программ «Неньютон-II» [11], который написан в среде Compaq Visual Fortran Professional Edition 6.6.G, имеет модульную структуру, ориентирован, прежде всего, на расчёт шатунных и коренных подшипников коленчатого вала двигателей внутреннего сгорания с учётом неньютоновских свойств смазочного материала; конструктивных особенностей (источников подачи смазки), скоростного и нагрузочного режимов работы поршневой или роторной машины. Однако, наряду с этим, может использоваться для расчёта других динамически нагруженных подшипников скольжения.
Результаты расчёта гидромеханических характеристик шатунного подшипника коленчатого вала двигателя 4ЧН 13/15. На первом этапе расчётных исследований была выполнена оценка влияния зависимости вязкости масла Mobil XHP Delvac 1GW-4 от давления и времени релаксации на ГМХ шатунного подшипника. Учёт зависимости вязкости масла от давления приводит к заметному увеличению hj^ и inf hmin на 18,5 % и 2G,G %, соответственно. Остальные
*
гидромеханические характеристики меняются незначительно - значение N увеличивается на 3,3 %, что сопровождается незначительным ростом температуры смазочного слоя Тэ на G,3 %. Значения
*
расхода смазки QB и максимального гидродинамического давления sup pmax увеличиваются на
1,8 % и 2,7 %, соответственно. Увеличение времени релаксации смазки приводит к увеличению
** среднеинтегральных значений потерь мощности на трение N и толщины смазочного слоя hmin .
Затем была выполнена оценка совместного влияния неньютоновских свойств масел на гидромеханические характеристики шатунного (табл. 2, рис. 5, б) и коренных подшипников коленчатого вала двигателя 4ЧН 13/15.
Таблица 2
Гидромеханические характеристики шатунного подшипника коленчатого вала двигателя 4ЧН 13/15
Масло Хаі рактеристика
inf hmin * мкм sup Pmax МПа * N , Вт h* ■ "min 5 мкм * pmax , МПа /с * * R С °
М8Г2к (SAE 2G) 1,7241 3G3,6 5б5,2 4,897 64,7G G,2124 Юб,3
Castrol SLX GW-3G 2,G461 1,88б2 284,3 383,2 б12,4 49G,G 5,3G5 4,784 6G,96 бб,44 G,2G49 G,19G5 1 G8,3 1 G5,7
GM 5W-3G 2,3391 1,9382 288,б 387,9 б32,б 477,3 5,337 4,714 6G,59 66,G8 G,2162 G,2G22 1G7,9 1G4,4
Shell Helix Ultra 5W-4G 2,б151 2,G992 284.7 367.8 б73,8 5G3,1 5,542 4,94G 59,34 б3,88 G,213G G,2GG8 1 G9,3 1 G5,3
Mobil XHP Delvac 1GW-4G 2,5891 2,1222 282,9 363,G б75,5 524,8 5,597 5,G88 59,G7 б2,99 G,2G96 G,2G17 1G9,7 1 G5,9
Shell Rimula R2 15W-4G 2,1G91 1,9192 282,7 3G2,6 б22,5 528,G 5,519 5,G19 59,95 б2,73 G,2G38 G, 1912 1 G8,7 Юб,9
Liqui Moly Molygen 5W-5G 2,66G1 2,G8G2 2б8,б 277,3 725.2 556.2 6,G38 5,534 5б,7б 59,27 G,1872 G,1767 113,7 1G9,2
Примечания: 'без учёта неньютоновских свойств (вязкость-функция температуры и давления); 2с учётом неньютоновских свойств (вязкость-функция температуры, давления, градиента скорости сдвига;
время релаксации смазки X = 5 • 10~4с ).
Анализ результатов свидетельствует о заметном влиянии неньютоновского поведения моторных масел на гидромеханические характеристики шатунного подшипника двигателя 4ЧН 13/15. В частности, в зависимости от марки и класса вязкости масла среднеинтегральные значения минимальной толщины смазочного слоя h*m[n снижаются на 8-12 %, потери мощности
*
на трение N - на 15-25 %. Эффективная температура смазочного слоя Тэ уменьшается на 1,5-4 %. Снижение значения inf hmin составляет 7-22 %. Среднеинтегральные значения максимальных
1
Л=
1 —j
V S vs и 1
1 \ V л \У Чч. 1
/ Л J
1) г
2) Mob il D slva< IP 1 0W- 40
0 60 120 180 240 300 360 420 480 540 00, град Рис. 5. Зависимость минимальной толщины смазочного слоя от угла поворота коленчатого вала
15,55
10 и I л aI I— V-I гидродинамических давлений pmax по-
»mm , х II___ >| I I Л 1_ вышаются на 4-9 %, при этом значения
мкмГ^^/ \ \ f sup pmax повышаются на 3-35 %. На
3-7 % уменьшается расход смазки Q*B , вытекающей в торцы подшипника.
Снижение потерь мощности на трение в шатунном подшипнике, работающем на загущенных маслах, относительно потерь, рассчитанных при работе подшипника на масле М8Г2к, составляет от 1,5 до 15,5 % (рис. 6).
С позиции гидродинамического режима смазки такое снижение потерь мощности на трение можно назвать «энергосберегающим эффектом».
С учётом того, что для шатунного подшипника коленчатого вала двигателя 4ЧН 13/15 допустимая минимальная толщина смазочного слоя составляет 1,9 мкм, то для надёжной работы необходимо использовать масла с вязкостью более 3 мПас (рис. 7).
Аналогичные расчёты были выполнены для коренных подшипников коленчатого вала двигателя 4ЧН 13/15.
Результаты расчётных исследований позволили сформулировать практические рекомендации по выбору класса вязкости масла для подшипников коленчатого вала двигателя 4ЧН 13/15. В частности:
- масла классов вязкости SAE 5W-30, 5W-40, 10W-40, 15W-40, 5W-50 обеспечат надёжную работу подшипников;
- наиболее предпочтительными с позиции энергосбережения двигателя являются классы SAE 5W-30, 5W-40, 10W-40, поскольку обеспечивают достаточную толщину смазочного слоя при наибольшем снижении потерь мощности на трение;
- масла классов SAE 20, 0W-30 не обеспечивают допустимого значения минимальной толщины смазочного слоя, то есть не могут быть рекомендованы для смазки двигателя 4ЧН 13/15.
мкм Заключение. Обобщая результаты
расчётов, необходимо отметить, что наибольшее влияние на ГМХ подшипников оказывает эффект временного
снижения вязкости. Влияние эффекта релаксации касательных напряжений в смазочном слое на работу подшипника требует дальнейших исследований.
1,5 Н---------------------------------------------- Реализация методики в виде ком-
плекса программ позволяет на ранней стадии проектирования сложнонагру-женных гидродинамических подшип-
10 2 0 3 0 4 0 Ц, мПа-с)
ников скольжения оценивать влияние Рис. 7. Зависимость р^етнм мчнчма^^ толщины конструктивных, режимных и экс-
смазочного слоя шатунного подшипника от вязкости, ’ г
измеренной при Т = 150 °С, у = 6,58103 с-1 плуатационных факт°р°в на их гвдр°-
15
10
0,99
7,15
6,58
1,59
0’№-30 5’№-30 5’№-40 10Ш-40 15Ш-40 5’№-50
Рис. 6. Снижение потерь мощности на трение в шатунном подшипнике, работающем на загущенных маслах, относительно потерь, рассчитанных при работе подшипника на масле М8Г2к
inf hm
2,5
2
R 2 Г = 0,64 0
h '\цоп = 1,9 мкм ♦ | 4
♦
механические характеристики, обосновывать рекомендации по совершенствованию их конструктивных параметров и подбору смазочного материала.
Представленная работа выполнена при поддержке Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России на 2009-2013 годы» и Российского фонда фундаментальных исследований (проект 10-08-00424).
Литература
1. Леванов, И.Г. Обзор реологических моделей моторных масел, используемых при расчётах динамики подшипников скольжения коленчатого вала / И.Г. Леванов // Вестник ЮУрГУ. Серия «Машиностроение». - 2010. - Вып. 15. - № 10 (186). - С. 54-62.
2. Paranjpe, R.S. Analysis of Non-Newtonian Effects in Dynamically Loaded Finite Journal Bearings Including Mass Conserving Cavitation / R.S. Paranjpe // Trans ASME Jour of Trib. - 1992. -V. 114. - P. 736-746.
3. Прокопьев, В.Н. Методы решения уравнений движения в задачах динамики опор коленчатого вала двигателей внутреннего сгорания / В.Н. Прокопьев, Ю.В. Рождественский, Н.В. Широбоков // Вестник Уральского межрегионального отделения Российской академии транспорта. - Курган: Изд-во КГУ, 1999. - № 2. - С. 79-82.
4. Прокопьев, В.Н. Многосеточные алгоритмы интегрирования уравнения Рейнольдса в задачах динамики сложнонагруженных подшипников скольжения / В.Н. Прокопьев, А.К. Бояршинова, Е.А. Задорожная // Проблемы машиностроения и надежности машин. - 2005. - № 5. -С. 16-21.
5. Влияние микрополярных свойств масел на динамику сложнонагруженных подшипников скольжения / Е.А. Задорожная, В.Г. Караваев, И.Г. Леванов, А.В. Чеснов // Вестник Самарского гос. аэрокосм. ун-та. - 2009. - № 3 (19). - С. 329-337.
6. Виссуссек, Д. Реологические свойства смазочных масел и их значение для практики / Д. Виссуссек; пер. Ф.Я. Гульбиса. - Свердловск, 1976. - 94 с.
7. Леванов, И.Г. Экспериментальные исследования реологических свойств всесезонных моторных масел /И. Г. Леванов, Е.А. Задорожная //Вестник ЮУрГУ. Серия «Машиностроение». -2011. - Вып. 17. - № 11 (228). - С. 70-76.
8. Караваев, В.Г. Совершенствование метода теплового расчета сложно-нагруженных опор скольжения / В.Г. Караваев, Ю.В. Рождественский, К.В. Гаврилов // Надежность и качество в промышленности, энергетике и на транспорте: тр. междунар. конф. - Самара, 1999. -
Ч. 2. - С. 76-77.
9. Прокопьев, В.Н. Динамика радиальных опор скольжения, смазываемых микрополярными жидкостями / В.Н. Прокопьев, Н.В. Анисимова // Техническая эксплуатация, надежность и совершенствование автомобилей: науч. тр. / Челябинск: ЧПИ, 1982. - № 276. - С. 48-65.
10. Paranjpe, R.S. Comparison between Theoretical Calculations and Oil Film Thickness Measurements Using the Total Capacitance Method for Crankshaft Bearings in a Firing Engine / R.S. Paranjpe, S.I. Tseregounis, M.B. Viola // Tribology Transaction. - 2000. - Vol. 43. - P. 345-356.
11. Совершенствование методики расчёта сложнонагруженных подшипников скольжения, смазываемых неньютоновскими маслами / В.Н. Прокопьев, Е.А. Задорожная, В.Г. Караваев, И.Г. Леванов // Проблемы машиностроения и надёжности машин. - 2010. - № 1. - С. 63-67.
Поступила в редакцию 10 июня 2011 г.
Леванов Игорь Геннадьевич. Аспирант, инженер кафедры «Автомобильный транспорт и сервис автомобилей», Южно-Уральский государственный университет. Область научных интересов -гидродинамическая теория смазки, неньютоновские жидкости, расчёт трибосопряжений поршневых и роторных машин. E-mail: IL74-LIG@rambler.ru, kafat@susu.ac.ru
Levanov Igor Gennadevich. The post-graduate student, the engineer of «Motor transport and service of cars» department of the South Ural state university. The area of scientific interests - the hydrodynamical theory of greasing, non-Newtonian fluids, calculation of tribounits piston and rotor-type machines. E-mail: IL74-LIG@rambler.ru, kafat@susu.ac.ru