Р.Ф. Калимуллин, Н.Н. Якунин
РАСЧЕТНАЯ ОЦЕНКА УСЛОВИЙ СМАЗКИ КОРЕННЫХ ПОДШИПНИКОВ АВТОМОБИЛЬНЫХ
ДВИГАТЕЛЕЙ
Разработана модель условий работы подшипников скольжения автомобильных двигателей в переходном смазочном процессе, для оценки которого использован параметр Рж, характеризующий продолжительность существования смазочного слоя. Рассмотрена связь данного параметра с нагрузочно-скоростными режимами, а также с конструктивным и технологическим состоянием кинематических пар.
Актуальность проблемы оценки условий смазки коренных подшипников. Коренные подшипники являются одними из трибосопряже-ний, ограничивающих долговечность автомобильного двигателя. Наиболее благоприятным в условиях эксплуатации режимом смазки для этих подшипников является жидкостный, когда изнашивание шеек и вкладышей незначительно. Поэтому оценка режима смазки до сих пор является актуальной, несмотря на имеющиеся решения в этой области / 1, 2, 3, 5 /.
В Оренбургском государственном университете первые работы по изучению условий смазки в подшипниках скольжения начали проводиться около пятнадцати лет назад. Экспериментально исследовались на машине трения СМЦ-2 подшипники скольжения в условиях переходного смазочного процесса в зависимости от радиальной нагрузки, частоты вращения вала и диаметрального зазора. Были определены оптимальные по условиям смазки области нагрузочно - скоростного режима подшипников /6, 7, 8, 9, 10/. В настоящее время на кафедре автомобильного транспорта ОГУ в лаборатории двигателей продолжается аналогичная работа, но уже объектом исследования являются коренные подшипники автомобильных двигателей.
Теоретические представления. Коренные подшипники являются нестационарно - нагруженными гидродинамическими подшипниками скольжения. Изнашивание вкладышей и шеек предотвращается разделением трущихся поверхностей тонким слоем масла. Под действием гидродинамического давления возникает реакция со стороны смазочного слоя -несущая способность смазочного слоя Ыж, которая противодействует внешней нагрузке на подшипник N.
В эксплуатации происходит разрушение смазочного слоя, которое ведет к локальному
контактированию трущихся поверхностей. Разрушение слоя происходит вследствие многих причин: недостаточной подачи, высокой температуры или разжижения масла; малой частоты вращения коленчатого вала; высоких нагрузок, передаваемых подшипником, и других.
Расчетные оценки условий смазки ведутся путем сравнения минимальной толщины смазочного слоя с критической толщиной (Ък = 2...7 мкм), средней температуры смазочного слоя с допустимой (1 = 1150С), оп-
ределением мест сближения траектории центра шейки за цикл. При этом внешняя нагрузка на смазочный слой от действия газовых и инерционных сил определяется по разным схемам расчета коленчатого вала /1/. Гидродинамический расчет подшипника /3, 5/ ведется на основе частных решений уравнения Рейнольдса для ламинарного течения жидкости при принятых упрощениях. Известны расчеты подшипников с учетом влияния кавитационных разрушений, смазочных канавок, девиации оси вала, макроотклонений подшипника.
На основе теоретического анализа авторами разработана методика расчетной оценки условий смазки коренных подшипников. Для этого построена математическая модель условий смазки нестационарно-нагруженного подшипника скольжения. Условия смазки оценены параметром продолжительности существования смазочного слоя Р.
Математическая модель базируется на теории рабочих процессов автомобильных двигателей /1/ и гидродинамической теории подшипников скольжения /3, 5/. Оценка параметра Рж производится путем сравнения внешней нагрузки N ( величина воздействия) на смазочный слой от рабочих процессов в двигателе и предельной несущей способности смазочного слоя Nж (способности сопро-
тивления).
Смазочный слой разрушается, если внешняя нагрузка N превышает его предельную несущую способность Nж. Предельная несущая способность характеризует способность смазочного слоя воспринимать внешнюю нагрузку не разрушаясь при достижения критической толщины Нк . При 50 % - ной вероятности неразрушения слоя “Nж - ^ = 0, где “Nж и - средние значения внешней нагрузки и предельной несущей способности смазочного слоя. При заданной вероятности неразруше-ния Р предельная несущая способность смазочного слоя “Nж должна быть больше ^ на величину, равную среднему квадратическому отклонению, помноженному на квантиль и :
^ж - ^ + и Б= 0, Р
гдГ£ = + &/'5 (Б^ж + SN - средние
квадратические отклонения предельной несущей способности и внешней нагрузки).
В расчетах Nж и N случайные величины входят в разных видах /4/:
а) У = Х1 + Х2 +... Х„ Б = (Б21 + Б22 + ... +
Б2/'5, У = “У + ир Б;
б) У = Х1 ЧХ2 Ч ... ЧХ, коэффициент вариации V = Б/“У = [(Б 1/Х1)2 + [(Б/Х)2 + ... +
[(Б„/Хп)2)]0’5 = (V2, + V22 + ... + V2/'5, У = “У (1+ ир ^;
в) У = Ха11 ЧХа22 Ч ... ЧХа” V = [(а^)2 + (а^2)2 + ... +(ау/)]°’5, У = “У ”(1+ ир V).
Случайные величины X можно принимать детерминированными при малости Б. по сравнению с Б, или V. по сравнению с V.
Нами выдвигается гипотеза, что вероятностной интерпретацией физического смысла параметра Рж является вероятность неразрушения смазочного слоя.
Так, если Рж =1, то смазочный слой разделяет поверхности и в этом случае формируется режим бесконтактного взаимодействия трущихся поверхностей (жидкостное трение); в случае, если Р = 0 - возникает режим контактного взаимодействия. Общим случаем в подшипниках является переходный смазочный процесс, характеризующийся поочередным возникновением контактного и бесконтактного взаимодействия, для которого 0 <Рж <1.
Таким образом, математическая модель условий смазки подшипников учитывает следующие конструктивные, технологические и эксплуатационные факторы: вязкостно-температурные свойства смазочного материала, геометрические параметры сопряжения, качество обработки поверхностей трения, условия работы двигателя, а также их вероятностный характер.
На первом этапе разработки методики принимаем, что внешняя нагрузка N и предельная несущая способность смазочного слоя Nж являются детерминированными величинами.
Этапы расчетной методики оценки условий смазки коренных подшипников, например, карбюраторного четырехцилиндрового, рядного двигателя сводятся к следующему.
1. Исходя из внешней скоростной характеристики рассчитываемого двигателя задаются максимальные значения среднего эффективного крутящего момента на маховике коленчатого вала M и частоты вращения колен-
max
чатого вала n . В координатном поле M - n
max
определяются узловые точки M i - n. для дальнейшего расчета.
2. В каждой узловой точке определяется среднее эффективное давление за цикл:
pe = I6 Ч M/(D2 S i), Па, где D - диаметр цилиндра, м;
S - ход поршня, м; i - число цилиндров, и среднее индикаторное давление: p. = pe + pm, Па,
гдеpm= 0,0034 + 0,00038ЧБ n - среднее давление механических потерь, Па.
Строится индикаторная диаграмма зависимости давления газов в цилиндре рг от угла поворота коленчатого вала a.
3. Определяются силы, действующие на коренные шейки (по разрезной схеме коленчатого вала).
На первую (пятую) шейку N = 0,5 [N2 1 + (К 1 - R 1)2]0,5, Н,
1(5) 7 L ш.ш.1 v рк1 пр1у J 7 7
на вторую (четвертую)
N2(4) = 0,5 [(Т1 - Т2)2 + (Крк1 - Kpк2)2]0,5, Н
Н,
2(4)
на третью N = 0,5 [(Т2
Тз)2 + (Крк2
рк2'
К- - М2]0^
ркЗ
0,5
где N . = (І2. + К21 + К2к
ш.ш-I I I „ R-шг. -
S ^ = (К2І + К2,ш + К\.)0'5’
I10=I~4) = [0,25 (p2 - p0 + p) p D Sin(a + Ъ)/ cos Ъ,
К (i=I...4) = [0,25 (p2 - p0 + p) p D2 cos(a +
b)/cos b,
p. = - m. R (p n/30)2(cos a + l cos 2a),
m. = m 1+ 0,275 m ,
1 п 7 ш7
l = R/Lui,
К„ = m R (p n/30)2,
R-ш шк vjr '
ККк = mK R (p n/30) ,
R = L m r (p n/30)2/l,
npI I пр УГ /
R З = 2 R
прЗ прІ
Здесь:
NmmI - сила, действующая на первую шатунную шейку, Н;
К . - сила, действующая на i- е колено вала по кривошипу, Н;
Т. - тангенциальная сила, направленная по касательной к окружности i - го кривошипа, Н;
К. - сила, направленная по радиусу i- го кривошипа, Н;
К - сила инерции вращающихся масс шатуна, Н;
К - сила инерции вращающихся масс кривошипа, Н;
КПР1, Кпр3 - силы инерции противовесов, действующие на первый (пятый) и третий коренные шейки соответственно, Н;
р0 - атмосферное давление, Па; р. - удельная сила инерции возвратно - поступательно движущихся масс, Па;
т. - возвратно - поступательно движущаяся масса, кг;
Ьш - длина шатуна, м; тп - масса поршневой группы, кг; тш - масса шатуна, кг; тшк - масса шатуна, сосредоточенная на оси кривошипа, кг;
тк - масса неуравновешенных частей одного колена вала без противовесов, кг; тп - масса противовеса, кг;
Ь - угол отклонения оси шатуна от оси цилиндра;
Я - радиус кривошипа, м; г - радиус центра масс противовеса, м;
I, I - расстояния между плоскостью противовеса и смежной коренной шейкой и смежными коренными шейками соответственно, м.
Таким образом, в результате расчета получается множество значений сил Ы(к, М, п, а), действующих на к -ю коренную шейку, в зависимости от крутящего момента М, частоты вращения коленчатого вала п, угла поворота коленчатого вала а .
4. Определяется предельная несущая способность смазочного слоя (подшипник с идеальной геометрией, абсолютно жесткий)
Ыж = р т п I й Ф/(30у2), Н, где т - динамическая вязкость масла, ПаЧс;
I - опорная длина шейки вала, м;
у - относительный зазор - разность между диаметром В подшипника и диаметром й шейки вала, отнесенная к диаметру шейки вала, у = (В - й) /й;
Ф - коэффициент нагруженности подшипника.
Коэффициент нагруженности подшипника определяется исходя из критической толщины смазочного слоя к .
кР
В результате расчета получается множество значений Ыж (п).
5. При текущем значении угла поворота коленчатого вала а за цикл от 0 до 7200 проводится сравнение Ы(к, М, п, а) и Ыж (п) (рисунок 1). В результате такого сравнения формируется функция
Р (к, М, п) = 8а /7200,
ж ' / жі 7
где аж. - протяженность і - го участка смазочного слоя в градусах.
Рисунок 1 Схема формирования продолжительности существования смазочного слоя Р в коренном подшипнике.
Результаты компьютерного моделирования. На основе математической модели разработан алгоритм расчетной оценки условий смазки в подшипниках автомобильного двигателя (рисунок 2).
Рисунок 2 Блок - схема алгоритмаоценки условий смазки в коренных подаипникахавгомобипшого двигателя
Данный алгоритм реализован на языке Turbo Pascal и применен для оценки условий смазки в коренных подшипниках рядного четырехцилиндрового карбюраторного автомобильного двигателя с рабочим объемом 1,5 дм3 в зависимости от нагрузочного и скоростного режимов работы двигателя.
После обработке данных программой STATISTICA 5.0 StatSoft, Inc на рисунке 3 представлены трехмерные графики зависимостей Px: = | (М, n) для каждого коренного подшипника (при монтажных зазорах 5О мкм) вида:
а)
Px1(5) = 0,957-10 -3 M+4Ч10 -5 n - 6,4Ч10 -9 п2+2,4Ч10 -7 M n - 1,6Ч10 -7 M2,
P = 0,789-2Ч10 -3 M-5m0
-7 M п- 5 Ч10 -7 M2,
п2+5,4Ч10
ж1(5)
= 0,915-2Ч10 -3 M+7,3i\10 -5 п-1,1Ч10
п2+4,8Ч10 -7 M n - 1,6Ч10 -6 M2,
для всех подшипников:
Px = 0,619-4Ч10 M n -1,7Ч10 -6 M2.
M-7H10 -8 п2+6,7Ч10
Анализ графиков показывает, что параметр Рж уменьшается практически по линейной зависимости при увеличении крутящего момента М, а зависимость от частоты вращения вала п имеет вид параболы с выраженным максимумом. Условия смазки третьего подшипника существенно отличаются от условий смазки первого (пятого) и второго (четвертого).
Дальнейший анализ позволяет определять области нагрузочно-скоростного режима работы двигателя, соответствующие определенному уровню параметра Рж в подшипнике (рисунок 4). Это дает возможность, например, выделить оптимальный нагрузочно - скоростной режим работы двигателя по критерию минимального изнашивания (Р ® Р тах) как конкретного подшипника, так и всей совокупности подшипников.
Выводы. Разработанная методика расчетной оценки условий смазки коренных подшипников открывает перспективы определения зависимости режима смазки от вероятностной природы параметров и нагрузочного, скоростного, теплового режима работы двигателя, вязкостно-температурных свойств моторного масла, геометрических параметров подшипников, уравновешенности коленчатого вала, условий подачи масла. Это позволит провести многопараметрическую оптимизацию по условию обеспечения жидкостной смазки в коренных подшипниках.
> -У
в)
г)
Рисунок 3. Зависимости параметра Рж в коренных подшипниках от нагрузочно-скоростного режима работы двигателя: а) первом (пятом), б) втором (четвертом), в) третьем, г) всей совокупности.
а)
б)
; : /
г/? / /" ' 1 < 1 И! Л 1 ! \ 1 : ! ■
! / /
я : / / I 1 г 1 1
•1: 1 / 1 ■ 1 ■
Шп М! ■ 1 П !11
П! 7/ И ; ■ |
и ■ [ 1 1 ь , ь
нч
Ь:/1. . 1 1 :\\
■3:И1
■ 1П
■ 1ййа
в) г)
Рисунок 4. Изопараметрические области для а) первого (пятого), б) второго (четвертого),
в) третьего и г) всей совокупности коренных подшипников.
Список использованной литературы
1.Двигатели внутреннего сгорания. В 3 кн. Кн. 2. Динамика и конструирование: Учеб. / В.Н. Луканин, И.В. Алексе-
ев, М.Г. Шатров и др.; Под ред. В.Н. Луканина. - М.: Высш. шк., 1995, 319 с.
2.Машино - ориентированные методы расчета комбинированных двигателей / Б.И. Иванченко, В.И. Каплан, К.Б. Цыреротов и др. - М.: Машиностроение, 1978, 168 с.
3.Расчет опорных подшипников скольжения: Справочник/ Е.И. Квитницкий, Н.Ф. Киркач, Ю.Д. Полтавский, А.Ф. Савин. - М.: Машиностроение, 1979, 70 с.
4. Решетов Д.Н. Детали машин: Учебник для студентов машиностроительных и механических специальностей вузов. - М.: Машиностроение, 1989, 496 с.
5.Справочник по триботехнике: В 3 т. Т. 2: Смазочные материалы, техника смазки, опоры скольжения и качения/
Под общ. ред. М. Хебды, А.В. Чичинадзе. - М.: Машиностроение, 1990. - 416 с.
6.Повышение долговечности транспортных машин: Учебное пособие для вузов/ В.А. Бондаренко, К.В. Щурин, Н.Н. Якунин и др.; Под ред. В.А. Бондаренко. - М.: Машиностроение, 1999.-144 с.
7. А.С. 1312444 СССР, МКИ3 О 01 N 3/56. Способ определения нарушения жидкостного режима трения подшипников скольжения / Абдрашитов Р.Т., Шевченко А.И., Якунин Н.Н.// Открытия. Изобретения. 1987. № 19.
8. А.С. 1523941 СССР, МКИ3 О 01 М 13/04. Способ динамической оценки состояния пары трения механизма / Абдрашитов Р.Т., Шевченко А.И., Якунин Н.Н.//Открытия. Изобретения. 1989. № 43.
9. Якунин Н.Н., Калимуллин Р.Ф., Баловнев С.В. Переходный смазочный процесс в подшипниках скольжения машин // Справочник. Инженерный журнал. 1999, № 12, С. 30-32.
10. Якунин Н.Н., Калимуллин Р.Ф., Баловнев С.В. Переходный смазочный процесс в подшипниках скольжения машин (продолжение) // Справочник. Инженерный журнал. 2000, № 1, С. 16-19.
Статья поступила в редакцию 20.12.99г. 5 8 ВЕСТНИК ОГУ 1'(4)2000