CC BY
83
Южно-Уральский государственный университет (национальный исследовательский университет) [email protected]
МОДЕЛИРОВАНИЕ СВЯЗАННЫХ ПРОЦЕССОВ В СИСТЕМЕ «УПРУГИЙ КОЛЕНЧАТЫЙ ВАЛ, ВРАЩАЮЩИЙСЯ В ПОДШИПНИКАХ ГИДРОДИНАМИЧЕСКОГО ТРЕНИЯ»
В работе предложена математическая и алгоритмическая модель, а также численные методы для моделирования системы коренных подшипников двигателя внутреннего сгорания. Новизна предлагаемого в работе подхода заключается в моделировании одновременного влияния нелинейных демпфирующих свойств смазочных слоев, обусловленных, в том числе, особенностями реологических свойств современных смазочных материалов; характеристик упругих связей шипов и подшипников; перекосов осей опорных и подвижных элементов трибосопряжений, вызванных комплексом внешних воздействий различной природы (динамических силовых, кинематических, тепловых).
Ключевые слова: коренной подшипник, связанные процессы, динамическая модель, смазка, неньютоновские свойства, коэффициент податливости, коленчатый вал, картер, двигатель внутреннего сгорания.
Анализ динамики многих механических систем осложняется наличием смазочных слоев, разделяющих движущиеся относительно друг друга поверхности упругих элементов конструкций и обеспечивающих передачу механического движения при минимальных потерях на трение. Это в первую очередь относится к осям и валам (прямолинейным, коленчатым) таких машин массового применения как двигатели внутреннего сгорания, поршневые компрессоры, насосы, турбоагрегаты. Учет действия смазочных слоев переводит проблему расчета таких объектов в класс нелинейных задач динамики упругих конструкций.
Кроме того, сама возможность существования смазочного слоя достаточной толщины при силовом взаимодействии упругих элементов конструкций обеспечивается соответствующим выбором формы поверхностей трения, положения и конфигурации источников для маслообес-печения, вязкостно-температурных и неньютоновских свойств современных смазочных материалов, обусловленных их структурной неоднородностью из-за наличия противоизносных и энергосберегающих присадок и другими факторами.
Исходя из этого, нелинейная задача динамики упругих конструкций должна быть дополнена гидродинамической, с учетом реологии современных смазочных материалов, и тепловой задачами смазки для каждого из сопряжений жидкостного трения. Такая постановка позволяет, при условии использования адекватных
методов моделирования, получать с приемлемой для инженерных приложений точностью результаты для отдельных (автономных) три-босопряжений.
Однако достаточно много трибосопряжений, например, подшипники поршневых и роторных машин, являются неавтономными и составляют единую трибосистему, характеризуемую устойчивыми связями трибологического характера. В такой трибосистеме совокупность трибосопряжений и рабочих сред находятся в механическом, гидравлическом, тепловом и других видах взаимодействия и характеризуются существенными обратными связями. Гидродинамические подшипники скольжения валов в поршневых и роторных машинах составляют именно такие трибосистемы, для которых наиболее характерны упругие связи. Причем связанными оказываются не только шипы трибо-сопряжений (подвижные элементы - цапфы валов), но и подшипники (неподвижные элементы), установленные в корпус турбокомпрессора или двигателя, жесткость которого конечна по величине и неравномерно распределена.
Таким образом, задача динамики и смазывания системы «вал - подшипники жидкостного трения» представляется связанной задачей теорий упругости, колебаний, термогидродинамической теории тонкого смазочного слоя, реологии жидких смазочных сред. Термин «связанный» предполагает два контекста использования: первый - как связанность процессов, происходящих в сопряжении жидкостного трения;
второй - как связанность сопряжений между собой.
Известные в настоящее время подходы к моделированию связанных процессов течения жидкости в тонких смазочных слоях и нелинейной динамики системы «упругий вал, вращающийся в подшипниках гидродинамического трения» сводятся к решению отдельных частных задач.
Первый подход заключается в решении задач динамики упругого вала на упругих опорах без промежуточного смазочного слоя. Такая постановка используется при оценке динамической нагруженности и прочности вала и элементов корпуса [1]-[5]. В ее рамках в настоящее время возможно достаточно точное трехмерное моделирование процессов, происходящих в упругих элементах конструкций. Однако, очевидно что нелинейные демпфирующие свойства подшипников влияют на динамическое поведение и прочность конструкции вала и корпуса.
Второй подход состоит в оценке работоспособности гидродинамических трибосопряже-ний как автономных подшипников. При таком подходе исключается взаимное влияние опор, а его целью является детальное исследование протекающих в смазочном слое процессов и разработка методов расчета подшипников различных конструкций [7], [8], [9]. Для тяжелонагру-женных трибосопряжений кривошипно-шатунного механизма ДВС такой подход является наиболее распространенным [10], хотя в его рамках невозможно оценить влияние многих существенных для оценки общей надежности ДВС конструктивных и технологических факторов: различной податливости опор картера, несоосностей их расположения относительно шеек коленчатого вала и др.
Третий подход характеризуется учетом упругих взаимосвязей между трибоопряжениями, однако при этом используются существенные допущения. Приближенные методы расчета траекторий движения коренных шеек коленчатого вала в подшипниках представлены в работах Мааса, Никитина А.М., Захарова С.М. и др. [11], [12] , для расчета траектории используется идея Холланда, в рамках которой принципиально невозможно учитывать особенности конструкции трибосопряжения, тем более корректно моделировать неньютоновские свойства современных смазочных материалов. В работах
Ветрова М.К. [6], Прокопьева В.Н. разработаны аппроксимирующие выражения на основе обобщенного уравнения Рейнольдса для конечной опоры и использован метод подвижностей. Использование современных средств вычислительной техники позволяют в рамках концепции решения связанной задачи динамики и смазки, предложенной в этих работах, использовать более точные модели и методы расчета гидродинамических процессов. Там же показано, что учет нелинейных демпфирующих свойств подшипников существенно влияет не только на динамическое поведение и прочность конструкции вала и корпуса, но и значительно изменяет гидромеханические характеристики подшипников. Однако данный метод не учитывает перекос осей вала и подшипников, а также разницу в величинах податливостей различных опор вала.
В работах Захарова С.М. [11], [13] рассматривается наиболее полная модель системы «коленчатый вал - подшипники - блок цилиндров», с учетом девиаций шеек и изгибной жесткости опор, хотя авторы для получения результатов вынуждены были использовать метод Холланда, недостатки которого приведены выше.
Несмотря на наличие большого количества исследований динамики и смазки систем «упругий вал, вращающийся в подшипниках гидродинамического трения», их несомненную значимость, в настоящее время отсутствуют общие подходы к решению таких задач с учетом перекосов осей подвижных и неподвижных элементов трибосопряжений, упругих свойств опор одновременно с учетом особенностей реологического поведения смазочных материалов, кинематического нагружения.
Обоснование системы уравнений
для моделирования связанных процессов
Математическая модель динамики связанной системы «упругий многоопорный вал, вращающийся в однослойных подшипниках гидродинамического трения», основана на механической модели, представленной на рисунке 1. Упругие взаимосвязи шипов коренных подшипников реализуются при помощи конструкции коленчатого вала, а упругие связи опор коренных подшипников - посредством конструкции картера двигателя.
Методика идентификации упругих свойств вала, а также линейной к^Х, к^У и угловой к,к®’У податливостей ] опоры представлена в работе [15]. При помощи нелинейно упругих элементов с податливостями к ^ (т) к (т) и к^ ’Х (т) к^’У (т) моделируется влияние смазочных слоев, разделяющих поверхности трения коренных шеек коленчатого вала и картера двигателя. Эти значения податливостей зависят от толщин смазочных слоев, которые изменяются во времени т , при этом центры шипов движутся по некоторым траекториям под действием системы сил г1’^, где I - номер колена вала, 8=1,..,7 - номер силы (см. рисунок 1). Индекс 8 соответствует номеру внешней нагрузки, действующей на коленчатый вал со стороны шатунов (8=3,5), а также инерционные нагрузки от противовесов (8=1,7) и массовых элементов коленчатого вала 8=2,4,6).
Мгновенные значения податливостей могут быть определены из соотношений
кЦ (т)= Х] (г)/ ГХ (т) ; к У (т) = (т)/ Г/ (т) ;
к
9,1 L ,х
(т)=а(т)/-
мх(т);
к
9,]
4 У
(г)=р](г)/
мі (т)
Здесь Т]Х (т) Г/ (т), (т) М/(т) мгновенные
значения сил и моментов, действующих на^'-й подшипник; Х7(т) /7 (т), а-'(т)в-'(т) - положение центра каждой/-й коренной шейки коленчатого вала.
Упругая подзадача
Так как механическая модель (рисунок 1) представляет собой статически неопределимую стержневую пространственную раму на упругом основании, то для определения сил и моментов, действующих на подшипники
ГХ (т) Г/ (т), Ы}х(т) М/(т) ,необходимо раскрыть
статическую неопределимость. Для полноопорных валов для раскрытия статической неопределимости наиболее эффективно использование уравнения пяти моментов
[5] {М}={*} (2)
где [5] - ленточная матрица коэффициентов влияния,
х х ( 1г^’]=1,п „Х,]=1,п „Х,]=1,п
5у=5а[ к*,Х ,к*,у ,кь,х ,
(1)
к9,] ,к9Х ,к?’] к9,]
4, у’ ’ 5,у ’ ’ 4, у
Рисунок 1. Динамическая модель коленчатого вала (а) и подшипника (б)
{M } - вектор-столбец искомых моментов,
{r} - вектор-столбец правых частей - перемещений в каждом направлении от внешнего силового и кинематического нагружения,
R = R (Ff1 n’(s), Fj=1 n’(s), K^j1’n, K*’j=1’n,
K k,j = hn KV,j д д 4
KL,x , KL,y’ case ^ necks ) •
Здесь Д necks, Д case - несоосности шеек коленчатого вала и отверстий под подшипники в картере двигателя, вызванные технологическими факторами и тепловым нагружением. Методика оценки теплового состояния картера двигателя представлена в работе [15], [16]. При помощи величин Дnecks, Дcase моделируется кинематическое нагружение коренных опор и коленчатого вала. Методика раскрытия статической неопределимости вала на опорах не только с линейной, но и с угловой податливостью, более подробно рассмотрено в работе [17]. Там же представлены выражения для определения сил и моментов, действующих на подшипники FJX(т) (т) MX(т) MY(т) исходя из значений
вектора-столбца моментов {M}.
Динамическая подзадача
Положение центра каждой j-й коренной шейки коленчатого вала XJ (т) YJ (т), aJ (т) в1 (т) определяется из системы уравнений движения (индекс j опущен)
го-диссипативных свойств опор с жидкостном смазкой; ис - гироскопический момент.
По аналогии с автономными тяжелонагру-женными подшипниками силами инерции шипов можно пренебречь в виду их малости по сравнению с действующими силами.
В отличие от системы уравнений движения для автономного подшипника, в системе (3) векторы внешних нагрузок г, ыр зависят не только от времени, но также и от характеристик демпфирования всех опор, а, следовательно, от координат и скоростей центров всех коренных шеек коленчатого вала (1).
Учитывая особенности рассматриваемой конструкции и малость величины инерционных членов, система (3) относится к классу жестких систем дифференциальных уравнений. Для ее решения использованы неявная численная схема и метод функций дифференцирования назад, хорошо зарекомендовавшие себя при решении системы уравнений движения для автономного подшипника.
Гидродинамическая подзадача
Для определения гидродинамических реакций смазочного слоя каждого подшипника необходимо знать поля гидродинамических давлений р(р,z,т) в тонком смазочном слое. Функция р(р,z,т) определяется решением дифференциального уравнения Рейнольдса для неньютоновской жидкости р(р,z,т) [18]-[20]
mU = F (U ,U ,у,у,т) + R (U ,U ,J,f,-) - ( ті I - - ( т2 I -
JY = Mp (U ,U ,yj-) + Mr (U ,U ,y,f,-) + д др hk+2 1 J, ^ 1 г 1-^ V _ др Рдр д + — дz hk+2 ф2 V _ф1_ ф0 _ др
+Mg (U ,U ,у,у,т). (3)
Здесь U = {X, Y}, у = {а, в} = {-у-cosе, у- sinе},
, и ,У - векторы линейных и угловых координат
и ускорений шипа; у - угол перекоса между осями шипа и подшипника; е - угол между плоскостью перекоса и координатными плоскостями; m, J - инерционные характеристики шипа; F = {FX, Fy }, MF = {Mx , My} - силы и моменты, приложенные шипу, определяемые решением системы (2); R, MR - результирующие векторы сил и моментов, полученные непосредственным интегрированием полей гидродинамических давлений в смазочном слое без использования каких либо упрощающих допущений о характере упру-
_Э_
др
о ph ґ \ |1 Ф 1
J 0 ~в~ llV
+-— (ph). дт
(4)
Здесь р - плотность смазки; р = (р - ра )/2 / ь -
функции безразмерного гидродинамического давления и толщины смазочного слоя соответственно, р, ра - размерные значения гидродинамического и атмосферного давлений, /, о - относительные значения зазора и угловой скорости движения шипа, і - безразмерная координата по ширине подшипника, ф - угловая координате подшипника; к - коэффициент, характеризующий степень неньютоновского поведения смазки;
— ^2 — к I —* і—
фк = 1 у /V “У, к = 0,1,2 ; у - безразмерная коор-у1
дината по нормали к поверхности шипа; Ц - безразмерная вязкость смазки как функция скорости сдвига, температуры и давления [21]
ц* = (!2)(к~1')!2■ с1 ■ е(С2/(те+сз')')+в(те)■р (5)
Здесь 12 = (дУх/ Эу)2 + (ЪУ21 Эу)2 - второй инвариант скоростей сдвига; Те -температура смазочного слоя Те; С, С2, С3 - константы смазки, определяемые экспериментально; в(Те) - пьезокоэффициент смазочного материала, который зависит от температуры и химического состава.
Методика определения характеристик смазочного материала представлена в работе [22]. Соотношение (5) позволяет учитывать структурную неоднородностью смазочного материала, вызванную наличием противоизносных и энергосберегающих присадок.
Размерная функция толщины смазочного слоя с учетом перекоса осей шипа и подшипника записана в виде
к(ф,Z\,т) = к (ф,Z\) -в(т)^(ф - 8(т)) +
^ 2S(т)
+ Z1— ------с°^ф - £(т))
В
(6)
Здесь е(т) 5(т) - эксцентриситет и угол положения линии центров в срединном сечении опоры; 8(т) - расстояние между геометрическими центрами шипа и подшипника на торцах опоры; В - ширина подшипника; 21 е [-в/2;я/2]; к*(р, 21) - зазор в опоре при совпадении осей шипа и подшипника, учитывающий геометрию поверхностей трения (некруглость, конусность и другие отклонения формы поверхностей трения от идеальных кругло-цилиндрических) и пространственную геометрию смазочного слоя.
Уравнение Рейнольдса интегрируется с учетом граничных условий Свифта-Штибе-ра, а также наличия источников смазки на поверхностях трения.
Решение обобщенного уравнения Рейнольдса выполняется численно с использованием оригинальной версии многосеточного метода определения давлений в тонких смазочных слоях неньютоновской жидкости с произвольной геометрией. Метод позволяет учитывать конструктивные особенности трибосопряжений, способы подачи смазки, вязкостно-температурные и простейшие неньютоновские свойства смазочного материала в широком диапазоне скоростных и нагрузочных параметров.
Тепловая подзадача
Теплопряженность каждого коренного подшипника оценивается решением уравнения теплового баланса за цикл работы двигателя тс [24]. Среднеинтегральное значение температуры определяется как
Те = То + N*/0*рст . (7)
* 1 Тс 1 Т с
Здесь N * = — | N (т)1т , О = — | 0(т)1т - сред-Тс о Тс о
ние за цикл значения потерь мощности на трение в слое смазки и расход смазки через торцы опоры; N(т) <2(т) - мгновенные значения потерь на трение и расхода смазки соответственно; для их расчета используется диссипативная функция [22]; р и ст - плотность и теплоемкость смазки соответственно.
Уравнение теплового баланса^'-го подшипника решается итерационным методом, причем на каждой итерации при постоянной температуре определяется вся траектория центра )-го шипа и многократно решается уравнение (4) с учетом изменения вязкости от температуры, вычисляемого по соотношению (5). Некоторое ускорение расчетов может быть достигнуто, если уравнение теплового баланса решать на каждом шаге определения траектории движения центра шипа. Тогда под Тс понимается длительность шага расчета траектории центра шипа. По окончанию цикла нагружения для проверки сходимости решения по температуре выполняется расчет цикла нагружения при постоянной температуре. Если изменение среднеинтегрального значения температуры за цикл не превышает 1%, то сходимость считается достигнутой.
Таким образом, решение задачи моделирования связанных процессов течения жидкости в тонких смазочных слоях и нелинейной динамики системы «упругий вал, вращающийся в подшипниках гидродинамического трения» основывается на совместном решении упругой (2), динамической (3), термогидродинамической подзадач (4)-(7). В рамках представленной математической модели возможно одновременно учитывать: а) нелинейные демпфирующие свойства смазочных слоев, обусловленные особенностями реологических свойств современных смазочных материалов; б) характеристики упругих связей шипов и подшипников; в) перекосы осей опорных и подвижных элементов трибосопряже-
ний, вызванных внешними воздействиями различной природы.
Вычислительные алгоритмы и программный комплекс для численного моделирования связанных процессов в системе «упругий многоопорный вал, вращающийся в подшипниках гидродинамического трения» основаны на итерационной схеме решения системы уравнений (2)-(4), (7). Итерационный процесс организован по нагрузкам, действующим на коренные опоры. Расчет начинается в предположении, что к£’Х (т)=к^у (т)= = КрХ (т ) = кру (т )= о, у = 1, п. Последовательным решением систем уравнений (3), (4), (7) для каждого подшипника определяются траектории движения центров шипов. Затем потоотношению (1) вычисляются величины К'Ь’Х (т),К^у (т),КрХ^ Кру (т), раскрывается статическая неопределимость (уравнение (2)) и определяются новые значения нагрузок, действующих на все опоры коленчатого вала в каждый момент
цикла нагружения. Если относительная интегральная ошибка между вычисленными «новыми» и «старыми» значениями нагрузок отличаются менее чем на 5%, то итерационный процесс прекращается. В противном случае динамическая, гидродинамическая и тепловая подзадачи решаются заново. Практика показывает, что количество итераций, обычно, не превышает 7-10 шт.
Область применения предложенной методики - моделирование тяжелонагруженных коренных подшипников коленчатых валов, оценка допусков на расположение отверстий, выполненных в картере двигателя для установки подшипников скольжения, оценка допусков на биения коренных шеек коленчатых валов в двигателях внутреннего сгорания. Ее использование позволит комплексно подойти к оценке надежности коренных подшипников, коленчатого вала и картера двигателя.
10.11.2013
Работа выполнена при поддержке Российского фонда фундаментальных исследований
(грант № 13-08-00875)
Список литературы:
1. Рыжов В.А. Разработка дизелей нового поколения на коломенском заводе // Двигателестроение. - 2009.- № 2. - С. 18-20.
2. Раенко М.И., Рыжов В.А. Применение иерархической системы моделей для оценки прочности деталей двигателей // Двигателестроение. - 2009. - № 2. - С. 21-26.
3. Антропов Б.С., Соцкая И.М., Ананьин М.Ю. Работоспособность вкладышей коленчатого вала // Тракторы и сельскохозяйственные машины. - 2008.- №8. - С. 44-45.
4. Абрамишвили М.М., Чистяков В.К. Динамическая напряженность коленчатого вала и его опор V -образного восьмицилиндрового двигателя, работающего по внешней скоростной и нагрузочной характеристикам // Двигателестроение. -1990. - № 1. - С. 10-12.
5. Краснокутский А.Н., Трифонов Ю.Ю. Расчет коленчатого вала на прочность по неразрезной схеме // Сборник научных трудов по проблемам двигателестроения, посвященный 175-летию МГТУ им. Н.Э. Баумана. - М.: МГТУ им. Баумана. -2005. - С. 96-102.
6. Ветров М.К. Разработка метода расчета параметров, характеризующих нагруженность подшипников многоопорных коленчатых валов поршневых машин: дис. ... канд. техн. наук. - Челябинск. - 1984. - 201 с.
7. Задорожная Е.А. Совершенствование и расширение области применения метода расчета динамики и гидромеханических характеристик опор скольжения с плавающими втулками: дис.. канд. техн. наук. - Челябинск. - 2002. - 171 с.
8. Ли, Роде. Статические и динамические характеристики радиальных подшипников с плавающей втулкой // Труды американского общества инженеров-механиков. Проблемы трения и смазки. - М.: Мир. 1982. - №3. - С. 64 - 70.
9. Прокопьев В.Н., Бояршинова А.К., Задорожная Е.А. Характеристики устойчивости подвижных элементов опор ротора турбокомпрессора систем наддува двигателей транспортных машин // Вестник Уральского межрегионального отделения Академии транспорта. - Курган. - 2001. - С. 180-186.
10. Гоц А.Н., Эфрос В.В. Расчет гидродинамических характеристик подшипников на неустановившихся режимах // Тракторы и сельскохозяйственные машины. - 2007. - №1. - С. 41-43.
11. Захаров С.М., Тарсис Ю.Л., Шорох Е.А. Совместный расчет многоопорного коленчатого вала и подшипников скольжения // Вест-ник машиностроения. -1985. - №1. - С. 5-7.
12. Семенов Г.И. Определение усилий в коленчатом валу с учетом его деформаций // Автомобильная промышленность. -1975. - № 4. - С. 20-30.
13. Захаров С.М., Сиротенко В.И., Жаров И.А. Моделирование работы трибосистемы «коленчатый вал-подшипники -опоры блока цилиндров» двигателей внутреннего сгорания // Трение и износ. - 1995. - Т16. - № 1. - С. 47-54.
14. Мыльников А.А. Экспериментальные исследования упругих характеристик коленчатого вала и блок-картера четырехцилиндрового рядного двигателя типа ЧН 13/15. - Вестник ЮУрГУ, серия «Машиностроение». - 2012. - № 12 (271). -С. 218-222.
15. Rozhdestvensky Y, Khozenjuk N, Mylnikov A and et al. Modeling of a multisupporting cranckshaft tribo-system of internal combustion engine. // Proceedings of the 15th Nordic Symposium on Tribology «NORDTRIB 2012», Trondheim, Norway, 1215 June 2012. CD - Р. 163.
16. Rozhdestvensky Y., Khozeniuk N., Mylnikov A., Levanov I., Romanov V. Modeling of the Main Bearings of a Multi-Supporting Crankshaft of the Internal Combustion Engine. // World Tribology Congress 2013, 8 - 13 September 2013, Torino, Italy. WTC2013 Proceedings, ISBN 9788890818509.
17. Мыльников А.А., Хозенюк Н.А. Методика оценки нагруженности опор коленчатого вала с учетом упругих свойств картера двигателя. // Альманах современной науки и образования. - 2013. - № 12. В печати
18. Прокопьев В.Н., Задорожная Е.А., Караваев В.Г. и др. Cовершенствование методики расчета тяжелонагруженных подшипников скольжения, смазываемых неньютоновскими маслами. // Проблемы машиностроения и надежности машин. - 2010. - № 1.- С. 63-67.
19. Paranjpe R, Tseregounis S and Viola M. Comparison between theoretical calculations and oil film thickness measurements using the total capacitance method for crankshaft bearings in a firing engine // J of Tribol Trans 2000. - № 43. Pp. 345-356.
20. Zhang R, Xueming H, Yang S and et al. Perturbation solution of non-newtonian lubrication with the convected Maxwell model // Trans of the ASME 2005. - № 127. - Pp. 302-305.
21. Задорожная Е.А., Леванов И.Г., Мухортов И.В. Применение неньютоновских моделей смазочных жидкостей при расчете сложнонагруженных узлов трения поршневых и роторных машин. // Трение и смазка в машинах и механизмах. - 2011. -№ 7. - С. 22-30.
22. Леванов И.Г. Методика расчета гидромеханических характеристик тяжелонагруженных подшипников скольжения поршневых и роторных машин, смазываемых неньютоновскими маслами // Вестник ЮУрГУ, Серия «Машиностроение». - 2011. - № 18. - С. 34-44.
23. Mukhortov I, Zadorojznaya E and Levanov I. Reological Model of a boundary layer of lubricant. In: STLE Annual Meeting & Exhibition, Hilton Atlanta, USA, 15-19 May 2011. - New York: ASME. - Pp. 235-241.
Сведения об авторе:
Хозенюк Надежда Александровна, доцент кафедры автомобильный транспорт и сервис автомобилей Южно-Уральский государственный университет (национальный исследовательский университет), кандидат технических наук,
e-mail [email protected]
UDC 62-144.3; 621.822.173
Khozeniuk N.A.
MODELING OF INTERCONNECTION PROCESS IN THE SYSTEM «ELASTIC CRANKSHAFT ROTATING ON THE HYDRODYNAMIC FRICTION BEARINGS»
In this paper the mathematical model and numerical methods for modeling of the main bearing system of internal combustion engines are suggested. Novelty of approach offered in work consists in 1) modeling of simultaneous influence of nonlinear damping properties of lubricant layers depending from rheological properties of the modern lubricant layers, 2) elastic links of the journals and the supports, 3) skewness of the axis of movable and unmovable parts of bearings. These skewness are caused by external forces of the various nature (dynamic forces or/and displasements, thermal factors).
Key words: main bearing, interconnected process, dynamic model, lubricant, non- Newtonian properties, elastic compliance coefficient, crankshaft, crankcase, internal combustion engine.
