CC BY
136
УДК 621.436
МЕТОДИКА ОЦЕНКИ ГРАНИЧНЫХ УСЛОВИЙ ПРИ РЕШЕНИИ ЗАДАЧИ ОПРЕДЕЛЕНИЯ ТЕПЛОВОГО СОСТОЯНИЯ ПОРШНЯ БЕНЗИНОВОГО ДВИГАТЕЛЯ
Г.В. Ломакин, М.А. Мацулевич
Предложена универсальная методика оценки граничных условий при определении теплового состояния поршня бензинового двигателя 4Ч 8,2/7,56 методом конечных элементов, позволяющая учесть все составляющие теплового баланса и особенности газодинамического нагружения в ходе рабочего цикла.
Ключевые слова: поршень, граничные условия, газодинамическое нагружение, квазистационарный теплообмен.
Введение
Повышение мощности и снижение массы деталей поршневых двигателей в значительной степени ограничивается тепловой и механической нагруженностью. При решении проблем расчета и проектирования основных деталей, образующих внутрицилиндровое пространство двигателя, широкое применение нашел метод конечных элементов, позволяющий использовать закономерности термо- и газодинамики, теории упругости, усталостной прочности на начальной стадии проектирования при создании поршневых двигателей повышенного ресурса.
Граничные условия теплообмена в цилиндре двигателя
Оценка теплового состояния исследуемого поршня на установившемся режиме работы сводится к определению его температурного поля в выбранном критическом сечении. Нахождение температурного поля поршня предполагает решение дифференциального уравнения теплопроводности, что можно осуществить, имея в распоряжении условия однозначности, которые включают: распределение температуры в сечении поршня в начальный момент времени - начальное условие, теплофизические свойства материала и геометрическую форму исследуемого сечения поршня, взаимодействие поверхности поршня и окружающей среды - граничные условия теплообмена. Под окружающей средой понимаются рабочие газы, омывающие поверхность поршня в цилиндре двигателя, масло, движущееся в каналах и полостях поршня, и детали цилиндропоршневой группы, находящиеся в контакте с поршнем [1].
При оценке теплообмена со стороны рабочих газов использованы граничные условия третьего рода - температура окружающей среды (рабочих газов в цилиндре) и закон теплообмена между рабочими газами и поверхностью поршня (коэффициент сложной теплоотдачи). Последний учитывает процесс конвективной и радиационной теплоотдачи.
Рассматривая процесс теплообмена как квазистационарный, теплоотдачу можно характеризовать средним по времени коэффициентом теплоотдачи от рабочих газов к поверхности поршня и некоторой результирующей по теплоотдаче температурой рабочих газов.
Наибольшие трудности представляет определение количества теплоты, отданной поршню рабочими газами, поскольку каждый участок тепловоспринимающей поверхности поршня оценивается своим коэффициентом теплоотдачи и определяющей температурой рабочих газов.
Для оценки граничных условий теплового нагружения поршня рассмотрены тепловой баланс и особенности газодинамического нагружения основных деталей цилиндропоршневой группы (ЦПГ) в ходе рабочего цикла (рис. 1).
Количество теплоты Q1, полученное поршнем от рабочих газов в цилиндре дизеля, определится зависимостью
Ql =aF АТ х, (1)
где ^ - площадь тепловоспринимающей поверхности поршня, м2;
аэ - коэффициент теплоотдачи от рабочих газов к поршню, Вт/м2 • °С;
АТ - разность температур поршня и рабочих газов в цилиндре, °С;
х - время подвода теплоты, с.
йх
Рис. 1. Основные составляющие теплового баланса поршня: 01 - количество тепла, полученное поршнем от рабочих газов; 02 - количество тепла, отведенное от поршня через поршневые кольца, юбку и гильзу цилиндра в систему охлаждения двигателя; 03 - количество тепла, отведенное от поршня смазочным маслом
При определении количества теплоты, полученной поршнем, необходимо учесть неравномерность распределения коэффициента теплоотдачи по поверхности камеры сгорания.
Одним из основных затруднений при решении практических задач теплопроводности является достоверное определение текущего коэффициента теплоотдачи а от рабочих газов. При анализе теплообмена в быстроходных двигателях с уменьшенным отводом теплоты для определения коэффициента теплоотдачи широко применяется зависимость Г. Вошни, позволяющая учитывать различия в интенсивности движения рабочего тела в отдельные периоды рабочего цикла, размеры цилиндра и дополнительную турбулизацию, вызываемую процессом сгорания, а также учитывать конвективную и радиационную составляющие [1-3]:
а = СдР 0,8 Т “°,53 £-°,2 Ж 0,8, (2)
для процессов газообмена:
Ж = 2,28Ст + 0,417С, (3)
для процессов сжатия, сгорания, расширения:
Ж = 2,28Ст + 0,308С„ + СУ, Та (Р - РПр) / (Ра Уа), (4)
где Сч = 2,3 • 10-0,5(Тт - 600) + 0,005;
Ст - средняя скорость поршня, м/с;
С„ - тангенциальная скорость вихря на радиусе, м/с;
С1 - коэффициент, для двигателя с внешнем смесеобразованием равен 3,24 • 10-3 м/(с • К);
Та - температура рабочего тела в начале сжатия, К;
Т - температура рабочих газов в цилиндре, К;
У, - рабочий объем цилиндра, м3;
Уа - объем рабочего тела в начале сжатия, м3;
Ра - давление рабочего тела в начале сжатия, МПа;
Рпр - давление рабочего тела без подачи топлива, МПа;
Р - давление рабочих газов в цилиндре, МПа;
Б - диаметр цилиндра, м.
Зависимость Г. Вошни используется для определения коэффициента теплоотдачи а при наличии результатов теплового расчета или индикаторных диаграмм давления в цилиндре, обра-
ботка которых предоставляет необходимые давления и температуры рабочего тела для требуемого угла поворота коленчатого вала, что позволяет учесть различия в интенсивности движения рабочего тела в различные периоды рабочего цикла, размеры цилиндра и дополнительную тур-булизацию, вызываемую процессом сгорания.
Для расчета температурных полей действительные условия нестационарного теплового нагружения заменяют некоторыми эквивалентными стационарными, характеризуемыми параметрами, полученными из условия равенства количества теплоты, воспринимаемой участком поверхности поршня в действительном и условно эквивалентном процессах. Эквивалентные параметры теплообмена вычисляются по следующим зависимостям [4]: хс
01 • | аdх, (5)
аэ = аср =х
э ср 0
ас
Tэ = --------------------------------------------------------------------------------------------------------------------------------------------^, (6)
'ср
откуда следует, что эквивалентный коэффициент теплоотдачи равен его среднему значению за цикл, а эквивалентная температура значительно выше интегральной температуры за цикл.
Для четырехтактных двигателей она составляет
Tэ *Тф + (0,6...0,8)( -273). (7)
Для определения температурного состояния поршня и распределения температур по поверхности и в объеме составляющих его элементов необходимо учесть неравномерность распределения коэффициента теплоотдачи по поверхности камеры сгорания, обусловленную ее формой и характером движения рабочих газов в процессе сгорания. Характер распределения коэффициента теплоотдачи по поверхности камеры сгорания и внешнего днища поршня задается в зависимости от текущего радиуса и описывается известным уравнением Славинского [4]:
\2п
= A + (В - А)
R
1 - A-С - А)
R
x | [ x
Я) "IЯ
(8)
ац апер
где А =---- и В =----— - относительные значения коэффициентов теплоотдачи в центре (ин-
аа
^шах ^шах
декс «ц») и на периферии (индекс «пер») огневой поверхности поршня; х - расстояние от оси цилиндра, м;
Я - радиус цилиндра, м; п = 0,5 / lg р;
р - относительный радиус, соответствующий ашах; к - показатель степени, определяющий крутизну функции.
Для поршня рассматриваемой конструкции значения параметров в уравнении (8) составляют А = 0,7; В = 0,75; п = 1,0 и к = 2,2. Максимальное значение коэффициента теплоотдачи ашах определяется зависимостью
ашах = а /
шах
(9)
где а - среднее значение коэффициента теплоотдачи для рассматриваемой поверхности;
Фх = ах / ашах - отношение коэффициентов теплоотдачи.
Полученные результаты эквивалентных значений коэффициента теплоотдачи, учитывающие неравномерность распределения коэффициента теплоотдачи по поверхности днища поршня и температуры газа в цилиндре нефорсированного двигателя 4Ч 8,2/7,56 на режиме номинальной мощности, определялись с использованием программного комплекса «NewGRAN», разработанного на кафедре «Двигатели внутреннего сгорания» Южно-Уральского государственного университета (рис. 2).
а х
а
шах
Рис. 2. Индикаторная диаграмма давления, температуры и коэффициенты теплоотдачи от газов в стенки внутрицилиндрового пространства двигателя 4Ч 8,2/7,56 на режиме номинальной мощности (Ре = 0,895 МПа, п = 5200 мин-1)
Количество теплоты 02, переданное от поршня через поршневые кольца и гильзу цилиндров в охлаждающую жидкость, определится зависимостью
02 = ^2 Г (Гп - Тж )Т , (10) где к2 - коэффициент теплопередачи от поршня через поршневые кольца и гильзу цилиндра в систему охлаждения дизеля, Вт/(м2 • °С);
2
г 2 - площадь поверхности контакта, м ;
Тп - температура поршня, °С;
Тж - температура охлаждающей жидкости, °С.
Процесс передачи теплоты от поршня через поршневые кольца в гильзу цилиндров и, далее, в охлаждающую жидкость осуществляется через две контактирующие поверхности, обусловленные наличием промежуточных деталей: компрессионных и маслосъемного колец [4].
Количество теплоты 03, отведенное смазочным маслом, определится зависимостью
03 =«з Г (тп -Тм)т, (11)
где а3 - коэффициент теплопередачи от поршня в охлаждающее смазочное масло, Вт/(м2 • °С);
Г3 - площадь поверхности поршня, омываемая маслом, м2;
Тм - температура смазочного масла, °С.
Количество тепла 04 можно установить как алгебраическую сумму
04 = 01 + 02 + 03. (12)
С другой стороны, согласно закону Фурье, можно записать
04 = X Г т (& / ёп), (13)
где X - коэффициент теплопроводности материала поршня, Вт/м • °С,
Г - площадь сечения поршня, м2;
т - время, с;
к /ёп - градиент температуры в сечении поршня, °С / м.
Тогда температурное поле в сечении поршня можно определить, приравняв выражения (12) и (13)
к / ёп = (01 + 02 + 03) / (X Г т). (14)
Известно, что наибольшее значение имеет составляющая теплового баланса, учитывающая
теплоту, полученную поршнем от рабочих газов в цилиндре. При этом коэффициент теплоотдачи
от днища поршня в масло зависит от способов охлаждения [4]. Так, в неохлаждаемых поршнях теплоотдача в масляный туман от внутренних поверхностей поршня характеризуется значением коэффициента теплоотдачи ам = 60...500 Вт/(м2 • °С). Относительно высокая эффективность принудительного охлаждения поршня достигается использованием преимуществ локального охлаждения.
Среднее значение коэффициента теплоотдачи аср при орошении поверхности внутреннего днища поршня смазочным маслом в качестве первого приближения определено с учетом рекомендаций [4]. Для плоской поверхности орошения внутреннего днища поршня максимальное значение а
ашах = 524,42л/ХРср п , (15)
где X - коэффициент теплопроводности смазочного масла, Вт/(м • °С); п - частота вращения коленчатого вала, мин-1.
Для конусной и кольцевой поверхности орошения внутреннего днища поршня минимальное значение а составит
ашт ашах к (16)
где к - коэффициент пропорциональности, определяется по зависимости
к = 32,74 [(лЯ + tg2 Р) / Р 0 33] (т / п р tg Р)0 68, (17)
Р - угол конусности поверхности внутреннего днища поршня;
2
р - площадь поверхности, орошаемой маслом, м ; т - плотность орошения.
Среднее значение коэффициента теплоотдачи аср определится зависимостью
аср 0,5 (ашт + ашах). (18)
При установке в цилиндре двигателя поршень контактирует с гильзой цилиндра прежде всего через поршневые кольца и смазочный слой, а также имеет место контакт направляющей части поршня через слой смазки непосредственно с гильзой цилиндра. Теплообмен между этими двумя контактирующими телами характеризуется граничными условиями четвертого рода, которые могут быть записаны в следующей форме для случая равенства температуры контактирующих тел [1]:
X (д t / д п) = Х2 (д t / д п)п; (19)
Ьп (х) = t2n (х), (20)
где п - координаты пространства.
На гранях элементов, совпадающих с сопрягаемыми поверхностями, теплообмен рассматривается как контактный и, следовательно, описываемый параметрами: коэффициентом контактного теплообмена ак, численно равным проводимости контакта, и температурой Т поверхностей контактирующих элементов. Коэффициент контактного теплообмена в этом случае, как и проводимость контакта, являются величинами обратными термическому сопротивлению контакта Як, которое зависит от идентичности формы, чистоты исполнения поверхности и усилия, обеспечивающего сохранение контакта сопрягаемых поверхностей.
Коэффициент контактного теплообмена между горизонтальными поверхностями контакта поршневых канавок и поршневых колец может быть определен с использованием рекомендаций Г.Б. Розенблита [5] по выражению:
ак = [2ХС] / (М + к2) + 2,1[(рХм) / (Зав)] 104, (21)
где Хс, Хм - теплопроводность среды и приведенная теплопроводность материала, (Вт/м • °С); к1, к2 - высота неровностей профилей контактирующих поверхностей, м; р - контактное давление, Па;
ав - предел прочности менее пластичного материала, Па.
Применительно к контактирующим поверхностям поршневых колец с выточками в поршне следует отметить в соответствии с условиями их работы достаточно малое общее термическое сопротивление Япк, которое представляет собой часть последовательно-параллельных сопротивлений трущихся граней. При этом для приближенных оценок можно использовать зависимость
Япк = {[Япк2 ЯпкЗ / (Япк2 + Япкз)] + Япк1 + (Ьк / Хк)}у, (22)
где Ьк, Хк - соответственно, ширина и коэффициент теплопроводности материала колец;
у - относительное термическое сопротивление поршневых колец, равное у быстроходных двигателей для первого кольца 1,2.. .1,7, для второго кольца - 1,2.. .1,4;
Япк1, Япк2, Япк3 - сопротивления трущихся граней поршневого кольца.
В этой связи для контактирующих поверхностей за условное значение коэффициента теплоотдачи целесообразно принять коэффициент контактного теплообмена, а в качестве температуры окружающей среды принять температуру поверхности сопрягаемых элементов, т. е. описать процесс теплообмена эквивалентными граничными условиями третьего рода, причем при моделировании теплового состояния поршня необходимо учесть действительные теплофизические свойства используемых материалов, смазочного масла и газовой среды.
Заключение
Представленная универсальная методика оценки граничных условий при решении задачи определения теплового состояния поршня исследуемого двигателя методом конечных элементов предусматривает следующую последовательность действий:
• определение характера изменения давления и температуры газов в цилиндре дизеля посредством синтеза рабочего цикла или индикаторных диаграмм давления в цилиндре;
• определение параметров теплообмена со стороны рабочих газов в цилиндре двигателя, смазочного материала и сопряженных элементов цилиндро-поршневой группы;
• определение эквивалентных значений коэффициента теплоотдачи, учитывающих неравномерность распределения коэффициента теплоотдачи по поверхности днища поршня и температуры газа в цилиндре исследуемого двигателя с использованием программного комплекса «NewGRAN»;
• описание геометрических характеристик элемента, создание конечно-элементной модели и определение условий теплообмена на границах расчетной области;
• указание теплофизических свойств материалов и выполнение расчетных процедур по решению уравнения теплопроводности для нахождения распределения температур в анализируемом элементе.
Литература
1. Лазарев, Е.А. Методические указания к практическим занятиям по дисциплине тепловая и механическая напряженность поршневых двигателей: учеб. пособие / Е.А. Лазарев. - Челябинск: Издат. центр ЮУрГУ, 2011. - 40 с.
2. Двигатели внутреннего сгорания. Теория поршневых и комбинированных двигателей /
Д.Н. Вырубов, Н.А. Иващенко, В.И. Ивин и др.; под ред. А.С. Орлина и М.Г. Круглова. - 4-е изд.,
перераб. и доп. - М. : Машиностроение, 1983. - 383 с.
3. Woschni, О. Е1пв Methode 2ыг VorausЬerechnung der Anderung des Brennverlaufs тШв1-schnellaufender Dieselmotoren Ьei geanderten BetrieЬs-Ьedingungen / О. Woschni, р. Anisitis //МТ2. -1973. - № 4. - Б. 160-165.
4. Костин, А.К. Теплонапряженность двигателей внутреннего сгорания / А.К. Костин, В.В. Ларионов, Л.И. Михайлов. - Л.: Машиностроение, 1979. - 222 с.
5. Розенблит, Г.Б. Теплопередача в дизелях /Г.Б. Розенблит. - М. : Машиностроение, 1977. -216 с.
Ломакин Георгий Викторович. Кандидат технических наук, доцент кафедры «Двигатели внутреннего сгорания», Южно-Уральский государственный университет (Челябинск), Lgeorge@ уап<!ех.ги.
Мацулевич Михаил Андреевич. Ассистент кафедры «Двигатели внутреннего сгорания», Южно-Уральский государственный университет (Челябинск), шшacule-vich@шail.ru
Поступила в редакцию 6 марта 2014 г.
Bulletin of the South Ural State University Series “Mechanical Engineering Industry” ____________2014, vol. 14, no. 1, pp. 26-32
METHODS OF BOUNDARY CONDITIONS IN SOLUTION OF DETERMINING THERMAL CONDITION OF THE PISTON GASOLINE ENGINE
G.V. Lomakin, South Ural State University, Chelyabinsk, Russian Federation, Lgeorge@yandex.ru, M.A. Maculevich, South Ural State University, Chelyabinsk, Russian Federation, mmacule-vich@mail. ru
A universal method of evaluation the boundary conditions for determining the thermal condition of the piston gasoline engine 4Ch 8.2 / 7.56 finite element method, which allows to take into account all the components of the heat balance and features a gas-dynamic loading during the working cycle.
Keywords: piston, the boundary conditions, the gas-dynamic loading, quasi-stationary heat.
References
1. Lazarev, E.A. Metodicheskie ukazanija kprakticheskim zanjatijam po discipline teplovaja i me-hanicheskaja naprjazhennost' porshnevyh dvigatelej [Methodical Instructions for Practical Lessons on Discipline Thermal and Mechanical Tension Reciprocating Engines]. Chelyabinsk: South Ural State University Publ., 2011, 40 p.
2. Vyrubov D.N., Ivashchenko N.A., Ivin V.I., e.a. Dvigateli vnutrennego sgoranija. Teorija porshnevyh i kombinirovannyh dvigatelej [Theory of Piston Engines and Combined]. Moscow, Mashi-nostroenie Publ., 1983. 383 p.
3. Woschni, G. Anisitis F. Eine Methode zur Vorausberechnung der Anderung des Brennverlaufs mit-telschnellaufender Dieselmotoren bei geanderten Betriebs-bedingungen. MTZ, 1973, no. 4, pp. 160-165.
4. Kostin A.K., Larionov V.V., Mikhailov L.I. Teplonaprjazhennost' dvigatelej vnutrennego sgoranija [Thermal Stress Combustion Engines]. Leningrad, Mashinostroenie Publ., 1979. 222 p.
5. Rosenblit, G.B. Teploperedacha v dizeljah [The Heat Transfer in Diesel Engines]. Moscow, Mashinostroenie Publ., 1977. 216 p.
Received 6 March 2014
