Тепловой и динамический расчет двигателя

Размещено на http://www.allbest.ru/

МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ

УО «Полоцкий государственный университет»

Кафедра автомобильного транспорта

КУРСОВАЯ РАБОТА

по дисциплине «Силовые агрегаты»

на тему: «Тепловой и динамический расчет двигателя»

Разработал: студент гр. 7ЭТМ (с) 41

Ильясов И.Х

Руководитель: Павлюк А.С

Барнаул 2017

Содержание

• Введение
• 1. Тепловой расчет двигателя
• 1.1 Процесс наполнения
• 1.2 Процесс сжатия
• 1.3 Процесс сгорания
• 1.4 Процесс расширения
• 1.5 Процесс выпуска

1.6 Индикаторные показатели

1.7 Эффективные показатели

1.8 Определение основных размеров двигателя

• 2. Динамический расчет
• 2.1 Построение индикаторной диаграммы

2.2 Силы, действующие в КШМ

2.3 Построение полярной диаграммы нагрузки

• 2.4 Построение графика суммарного крутящего момента
• Заключение
• Список литературы

Введение

Тепловой расчет двигателя служит для определения параметров рабочего тела в цилиндре двигателя, а также оценочных показателей рабочего процесса, для оценки мощностных и экономических показателей, позволяющих оценить мощность и расход топлива.

В основе методики расчета лежит метод В.И. Гриневецкого, в дальнейшем усовершенствованный Е.К. Мазингом, Н.Р. Брилингом, Б.С. Стечкиным и др.

Задачей динамического расчета является определение сил, действующих в механизмах преобразования энергии рабочего тела в механическую работу двигателя.

В настоящей работе тепловой и динамический расчеты выполняются для режима номинальной мощности.

1. Тепловой расчет двигателя

Исходные данные для теплового расчета поршневого двигателя внутреннего сгорания.

марка автомобильного двигателя ЗИЛ-130

номинальная мощность двигателя N_e = 110,3 кВт.

частота вращения коленчатого вала n = 3200 об/мин.

число тактов ф = 4

число цилиндров и расположение цилиндров i = 8V

степень сжатия е = 7.4

диаметр цилиндра м.

ход поршня м.

рабочий объем цилиндров (литраж двигателя ) л.

максимальный крутящий момент М_max= 402 Нм.

минимальный удельный расход топлива g_emin = 326 г/кВт?ч

тип двигателя карбюраторный

Для расчета двигателя в качестве топлива принимаем бензин А-92 с элементарным составом по массе: g_c= 0,85; g_n= 0,15; g₀ = 0.

Низшая теплота сгорания данного топлива Н_и = 44000 кДж/кг.

Давление и температуру окружающей среды принимаем равными р₀ = 0,1 МПа, Т₀= 298 К.

В начале сжатия температура отработавших газов для бензиновых

ДВС изменяется в пределах от 900 до 1100 К, для расчета принимаем

Т_r = 1000 К. Давление остаточных газов:

р_r = (1,05...1,25) · p₀ =1,15·0,1 = 0,115 МПа

Температура подогрева свежего заряда для бензиновых ДВС изменяется в пределах: ДТ= 5...30 К, принимаем ДТ = 15 К.

Величина потери давления на впуске для бензиновых ДВС

Д р_а = 0,12·0,01 = 0,012 МПа

1.1 Процесс наполнения

Давление газов в цилиндре в конце впуска:

МПа.

Коэффициент остаточных газов вычисляется по формуле:

Температура газов в цилиндре в конце впуска:

К

Коэффициент наполнения вычисляется по формуле:

1.2 Процесс сжатия

Давление и температура газов в конце сжатия вычисляются по формулам:

МПа

К

где n₁ = 1.38 - показатель политропы сжатия.

1.3 Процесс сгорания

Теоретически необходимое количество воздуха для сгорания топлива:

кмоль/кг

где g_c= 0,85; g_n= 0,15; g₀ = 0 - средний элементарный состав топлива для бензина А-92

Количество свежего заряда (кмоль) для бензиновых двигателей определяется по формуле:

где m_I - молекулярная масса топлива. Для бензина m_I = 110…120 кмоль/кг

б - коэффициент избытка воздуха б =0,8…1,3 для карбюраторных двигателей;

принимаем б = 0,85 на основных режимах.

кмоль/кг

Количество продуктов сгорания при работе двигателя на бензине при б < 1:

кмоль/кг

Теоретический коэффициент молекулярного изменения:

Действительный коэффициент молекулярного изменения:

Низшая теплота сгорания топлива Н_и = 44000 кДж/кг.

Потеря тепла вследствие неполноты сгорания топлива:

кДж/кмоль·К

Средняя мольная теплоемкость свежего заряда:

кДж/кмоль·град

Средняя мольная теплоемкость продуктов сгорания:

Коэффициент использования теплоты для бензиновых ДВС изменяется в пределах Принимаем

Максимальная температура сгорания подсчитывается по уравнению:

Подставив в уравнение найденное значение , имеем:

Решаем квадратное уравнение:

,

где А = 2,766 ·10^-3; В = 20,96; С = - 78263,05

К

Теоретическое максимальное давление цикла:

МПа

Степень повышения давления

Действительное давление конца сгорания МПа

1.4 Процесс расширения

Степень предварительного расширения для бензиновых двигателей = 1.

Степень последующего расширения

Температура в конце расширения

,

где средний показатель политропы расширения для бензиновых двигателей. Принимаем .

К.

Давление в конце расширения:

МПа.

1.5 Процесс выпуска

Параметры процесса выпуска выбираем в начале расчета (см. стр. 4):

температура отработавших газов для бензиновых ДВС изменяется в пределах от 900 до 1100 К, для расчета принято Т_r = 1000 К.;

давление остаточных газов: р_r = (1,05...1,25) · p₀ принято р_r = 15·0,1 = 0,115 МПа

Правильность выбора температуры Т_r проверяем по формуле профессора Е.К. Мазинга

К

Погрешность расчета составляет , что меньше предельной - 10 %, следовательно, нет надобности уточнять произведенные расчеты. бензиновый двигатель кривошипный шатунный

1.6 Индикаторные показатели

Среднее индикаторное давление т е о р е т и ч е с к о г о цикла:

МПа

Среднее индикаторное давление д е й с т в и т е л ь н о г о цикла:

МПа

где коэффициент полноты индикаторной диаграммы,

принимаем

Индикаторный КПД:

Удельный индикаторный расход топлива:

г/кВт·ч

1.7 Эффективные показатели

Величина механического КПД двигателя выбирается исходя из того, что для бензиновых ДВС изменяется от 0,7 до 0,85. Принимаем .

Среднее эффективное давление:

Эффективный КПД:

Удельный эффективный расход топлива:

г/кВт·ч

1.8 Определение основных размеров двигателя

По эффективной мощности, частоте вращения коленчатого вала и среднему давлению определяем литраж двигателя по формуле:

л

Рабочий объем одного цилиндра:

л,

где i - число цилиндров.

Принимаем соотношение S/D в соответствии с прототипом:

Диаметр цилиндра

м.

Ход поршня

м.

Эффективная мощность

кВт

Эффективный крутящий момент

Н·м

Часовой расход топлива

кг/ч

Средняя скорость поршня

м/с

Литровая мощность

кВт/л

Таким образом, найдены основные параметры рабочего цикла двигателя, индикаторные и эффективные показатели его работы, также определены основные параметры двигателя.

Выводы: точность расчетов температуры Т_r составляет 1,215 %;

отношение максимального крутящего момента прототипа к расчетному эффективному соответствует характеристике карбюраторного двигателя;

максимальный крутящий момент при

об/мин,

что соответствует технической характеристике прототипа.

2. Динамический расчет

Динамический расчет кривошипно-шатунного механизма выполняется с целью определения суммарных сил и моментов, возникающих от давления газов и от сил инерции. Результаты динамического расчета используются при расчете деталей двигателя на прочность и износ.

В течение каждого рабочего цикла силы, действующие в кривошипно-шатунном механизме, непрерывно изменяются по величине и направлению. Поэтому для характера изменения сил по углу поворота коленчатого вала их величины определяют для ряда различных положений вала через каждые 30 град. поворота коленчатого вала.

2.1 Построение индикаторной диаграммы

Индикаторная диаграмма строится в координатах p-V. Построение индикаторной диаграммы двигателя внутреннего сгорания производится на основании теплового расчета.

В начале построения на оси абсцисс откладывают отрезок АВ, соответствующий рабочему объему цилиндра, а по величине равный ходу поршня в масштабе ms, который в зависимости от величины хода поршня проектируемого двигателя может быть принят 1:1,1,5:1. Принимаем 1:1.

Отрезок ОА, соответствует объему камеры сгорания, определяется из соотношения

мм.

При построении диаграммы выбираем масштаб давления m_р=0,04 МПа/мм.

Затем по данным теплового расчета на диаграмме откладывают в выбранном масштабе величины давлений в характерных точках а, с, z', z, b, r.

Таблица 1

точки	a	c	z'	z	b	r
p, МПа	0,088	1,393	5,092	4,328	0,377	0,115

По наиболее распространенному графическому методу Бауэра политропы сжатия и расширения строим следующим способом.

Из начала координат проводим луч ОК под углом а₀ = 20° к оси координат. Далее из начала координат проводим лучи ОД и ОЕ под углами ₁ и ₂ к оси ординат. Эти углы определяют из соотношений:

Политропу сжатия строим с помощью лучей ОК и ОД. Из точки С проводим горизонталь до пересечения с осью ординат; из точки пересечения - линию под углом 45° к вертикали до пересечения с лучом ОД, а из этой точки - вторую горизонтальную линию, параллельную оси абсцисс. Затем из точки С проводим вертикальную линию до пересечения с лучом ОК. Из этой точки пересечения под углом 45° к вертикали проводим линию до пересечения с осью абсцисс, а из этой точки - вторую вертикальную линию, параллельную оси ординат, до пересечения со второй горизонтальной линией. Точка пересечения этих линий будет промежуточной точкой 1 политропы сжатия. Точку 2 находим аналогично, принимая точку 1 за начало построения.

Политропу расширения строим с помощью лучей ОК и ОЕ, начиная от точки z, аналогично построению политропы сжатия. Критерием правильности построения политропы расширения является приход ее в ранее нанесенную точку b.

После построения политропы сжатия и расширения производим скругление индикаторной диаграммы с учетом предварения открытия выпускного клапана, опережения зажигания и скорости нарастания давления, а также наносим линии впуска и выпуска. Для этой цели под осью абсцисс проводим на длине хода поршня S как на диаметре полуокружность радиусом R=S/2. Из геометрического центра О' в сторону н.м. т. откладываем отрезок

мм.

Величина O'O'₁ представляет собой поправку Брикса. Из точки O'₁ под углом г₀= =45° (угол предварения открытия выпускного клапана, который выбирается из таблицы или по прототипу) проводим луч O₁B₁. По учебнику [4], стр. 44 угол предварения открытия выпускного клапана равен 42 - 47°. Полученную точку В₁ соответствующую началу открытия выпускного клапана, сносим на политропу расширения (точка b₁)/

Луч O'C₁ проводим под углом и₀, соответствующим углу опережения зажигания (и₀ = 27° поворота кривошипного вала до в. м. т.), а точку С₁ сносим на политропу сжатия, получая точку с'₁. По справочникам угол опережения зажигания и₀ =25 - 27°. Затем проводят плавные кривые c₁'c” изменяя линии сжатия в связи с опережением зажигания и b₁'b” изменения линии расширения в связи с предварением открытия выпускного клапана. При этом можно считать, что точка b” находиться на середине расстояния bа , а ордината точки с” находится из соотношения р_с”=1,2р_с и откладывается на линии AZ'.

Развертку индикаторной диаграммы в координаты р-ц выполняем справа от индикаторной диаграммы. Ось абсцисс развернутой диаграммы располагаем по горизонтали на уровне линии р₀ индикаторной диаграммы. Длина графика (720° поворота кривошипного вала) делим на 24 равных участка, которые соответствуют определенному углу поворота коленчатого вала. Каждую точку на линии абсцисс нумеруем (0є, 30є, 60° ПКВ). По наиболее распространенному способу Ф. А.

Брикса дальнейшее перестроение индикаторной диаграммы ведем в следующей последовательности.

Полученную полуокружность делим вспомогательными лучами из центра О' на 6 равных частей, а затем из центра Брикса (точка Oi') проводим линии, параллельные вспомогательным лучам, до пересечения с полуокружностью.

Вновь полученные точки на полуокружности соответствуют определенным углам ц ПКВ. Из этих точек проводим вертикали до пересечения с соответствующими линиями индикаторной диаграммы. Развертку индикаторной диаграммы начинаем, принимая за начало координат положение поршня в в.м.т. в начале такта впуска. Далее для каждого значения угла ц на индикаторной диаграмме определяем величину давления в надпоршневой полости и заносим в табл.1.

Для проверки правильности построений и дальнейших расчетов сил производим аналитический расчет политроп сжатия и растяжения через каждые 30є поворота кривошипа по следующим формулам: политропа сжатия

,

где АВ и АО - определены в разделе 3.1;

S_i - перемещение поршня в i-ое положение из положения начала сжатия.

политропа расширения

,

где S_i - перемещение поршня в i-ое положение из положения начала расширения.

Результаты расчетов сведены в таблицу 2.

АВ = 80 мм.; АО = 10,67 мм.; n₁ = 1.38; n₂ = 1.3

Таблица 2

сжатие	расширение
№№ пол.	цi град.	pi МПа	№№ пол.	цi град.	pi МПа
7	180	0,088	14*	376,59	4,328
8	210	0,096	15	390	3,202
9	240	0,123	16	420	1,470
10	270	0,192	17	450	0,788
11	300	0,373	18	480	0,518
12	330	0,852	19	510	0,408
13	360	1,393	20	540	0,377

2.2 Силы, действующие в КШМ

Определяем силу давления газов на днище поршня для положений коленчатого вала, отстоящих друг от друга на 30 град. поворота коленчатого вала в пределах (0...720)⁰ поворота коленчатого вала. За начало отсчета принимаем такое положение кривошипа, когда поршень находится в начале такта впуска.

Силу давления газов на днище поршня определяем по формуле:

(1)

где р_г - индикаторное давление в каждом рассчитываемом положении, определяемое по построенной индикаторной диаграмме.

• Результаты расчета заносятся в табл. 3.
• Определяем силу инерции от возвратно-поступательно движущихся масс:

(2)

где л = 0,25 - отношение длины кривошипа к длине шатуна

• Масса поступательно движущихся частей кривошипно-шатунного механизма определяется из выражения:
• , (3)
• где доля массы шатуна, отнесенная к возвратно-поступательно движущимся массам.
• Значения m_n и m_ш вычисляются:
• кг.
• кг., где
• м^2, Теперь
• кг.
• Угловая скорость кривошипа
• мин^-1 (4)
• Радиус кривошипа
• м. (5)
• Результаты расчета заносятся в табл. 3.
• Находим суммарную силу, действующую в кривошипно-шатунном механизме. Определение этой силы ведем путем алгебраического сложения сил давления газов и сил инерции возвратно-поступательно движущихся масс: (6)
• Результаты расчетов заносим в таблицу 3.
• Определяем нормальную силу К, направленную по радиусу кривошипа
• (7)
• где
• l - длина шатуна
• Определяем тангенциальную силу Т, направленную по касательной к окружности радиуса кривошипа
• (8)

Результаты расчетов К и Т заносим в таблицу 4.

Графики сил Р_г и Р_j строим в одинаковом масштабе и на той же координатной сетке. На основании полученных графиков Р_г и Р_j на той же координатной сетке и в том же масштабе строим график суммарной силы Р_У.

Определение модуля силы Р_У для различных значений угла ц выполняем путем суммирования в каждой точке ординат графиков и с учетом их знаков и модулей сил Р_г и P_j из табл. 3.

Таким же образом строим графики сил Т и К.

Масштабные коэффициенты т_Рг = т_Рj = т_РУ = т_К = т_Т =200 Н/мм.

Результаты расчета сил Р_г, Р_j и Р_У

Таблица 3

ц, град	Рг , МПа	Рг , Н	cosц+лcos2ц	РJ , Н	PУ , Н
0	0,115	118,0	1,250	-11155,9	-11037,9
30	0,111	82,6	0,991	-8844,6	-8762,0
60	0,106	47,2	0,375	-3346,8	-3299,6
90	0,102	11,8	-0,250	2231,2	2243,0
120	0,097	-23,6	-0,625	5578,0	5554,4
150	0,093	-59,0	-0,741	6613,5	6554,5
180	0,088	-94,4	-0,750	6693,6	6599,1
210	0,096	-35,0	-0,741	6613,5	6578,5
240	0,123	182,3	-0,625	5578,0	5760,2
270	0,192	726,0	-0,250	2231,2	2957,2
300	0,373	2145,4	0,375	-3346,8	-1201,4
330	0,852	5912,4	0,991	-8844,6	-2932,3
360	1,393	10174,1	1,250	-11155,9	-981,8
376,59	4,328	33263,4	1,168	-10420,6	22842,8
390	3,202	24406,9	0,991	-8844,6	15562,3
420	1,470	10781,2	0,375	-3346,8	7434,5
450	0,788	5414,8	-0,250	2231,2	7646,0
480	0,518	3289,6	-0,625	5578,0	8867,5
510	0,408	2422,6	-0,741	6613,5	9036,1
540	0,233	1042,4	-0,750	6693,6	7736,0
570	0,126	204,5	-0,741	6613,5	6818,0
600	0,124	187,2	-0,625	5578,0	5765,2
630	0,122	169,9	-0,250	2231,2	2401,1
660	0,119	152,6	0,375	-3346,8	-3194,2
690	0,117	135,3	0,991	-8844,6	-8709,3
720	0,115	118,0	1,250	-11155,9	-11037,9

Результаты расчета сил Т и К

Таблица 4

ц, град	cos(б+)/cos	К , Н	sin(б+)/cos	Т, Н
0	1,000	-11037,9	0,000	0,0
30	0,803	-7036,2	0,609	-5337,0
60	0,308	-1016,1	0,977	-3223,4
90	-0,258	-579,1	1,000	2243,0
120	-0,692	-3843,9	0,755	4194,3
150	-0,929	-6089,2	0,391	2562,1
180	-1,000	-6599,1	0,000	0,0
210	-0,929	-6111,5	-0,391	-2571,5
240	-0,692	-3986,4	-0,755	-4349,8
270	-0,258	-763,5	-1,000	-2957,2
300	0,308	-370,0	-0,977	1173,6
330	0,803	-2354,7	-0,609	1786,1
360	1,000	-981,8	0,000	0,0
376,59	0,938	21425,2	0,354	8088,8
390	0,803	12497,0	0,609	9479,1
420	0,308	2289,4	0,977	7262,8
450	-0,258	-1974,2	1,000	7646,0
480	-0,692	-6136,8	0,755	6696,2
510	-0,929	-8394,7	0,391	3532,1
540	-1,000	-7736,0	0,000	0,0
570	-0,929	-6334,1	-0,391	-2665,1
600	-0,692	-3989,8	-0,755	-4353,5
630	-0,258	-620,0	-1,000	-2401,1
660	0,308	-983,6	-0,977	3120,4
690	0,803	-6993,9	-0,609	5304,9
720	1,000	-11037,9	0,000	0,0

2.3 Построение полярной диаграммы нагрузки

Полярная диаграмма нагрузок на шатунную шейку строится для определения величин, направления и точек приложения сил, действующих на шатунную шейку при различных положениях коленчатого вала.

В разделе 3 определены нормальная сила К, направленная по радиусу кривошипа

и тангенциальная сила Т, направленная по касательной к окружности радиуса кривошипа. Нагрузка вызывается силой, равной геометрической сумме сил К и Т.

Отрицательные значения силы К откладываем вверх, а положительные значения силы Т откладываем вправо. Сила S образует с горизонтальной осью угол :

Расчеты сводим в таблицу 6.

Результаты расчета силы N и угла г

Таблица 6

ц, град	К , Н	Т, Н	S, Н	г, град
0	-11037,9	0,0	11037,9	90,0
30	-7036,2	-5337,0	8831,3	127,2
60	-1016,1	-3223,4	3379,7	162,5
90	-579,1	2243,0	2316,5	14,5
120	-3843,9	4194,3	5689,3	42,5
150	-6089,2	2562,1	6606,3	67,2
180	-6599,1	0,0	6599,1	90,0
210	-6111,5	-2571,5	6630,5	112,8
240	-3986,4	-4349,8	5900,2	137,5
270	-763,5	-2957,2	3054,2	165,5
300	-370,0	1173,6	1230,6	17,5
330	-2354,7	1786,1	2955,4	52,8
360	-981,8	0,0	981,8	90,0
377,13	21425,2	8088,8	22901,2	290,7
390	12497,0	9479,1	15685,3	307,2
420	2289,4	7262,8	7615,1	342,5
450	-1974,2	7646,0	7896,7	14,5
480	-6136,8	6696,2	9083,0	42,5
510	-8394,7	3532,1	9107,5	67,2
540	-7736,0	0,0	7736,0	90,0
570	-6334,1	-2665,1	6871,9	112,8
600	-3989,8	-4353,5	5905,3	137,5
630	-620,0	-2401,1	2479,9	165,5
660	-983,6	3120,4	3271,8	17,5
690	-6993,9	5304,9	8778,2	52,8
720	-11037,9	0,0	11037,9	90,0

Диаграмма изменения силы S и есть полярная диаграмма сил, действующих на шатунную шейку, но без учета центробежной силы массы шатуна К_rш, отнесенной к его нижней головке:

, где

кг.

При установившемся движении сила Н.

имеет постоянную величину. Ее действие изменяет величину силы К, что учитывается переносом начала координат, построенной полярной диаграммы сил вниз по оси К на величину К_rш/m_K мм. Полученная точка О_ш является новым полюсом, а ранее построенная кривая относительно него будет полярной диаграммой сил, действующих на шатунную шейку.

Вектор, направленный из полюса О_ш к любой точке кривой на диаграмме, определяет в выбранном масштабе величину и направление R_шш нагрузки на шатунную шейку для соответствующего угла поворота коленчатого вала.

2.4 Построение графика суммарного крутящего момента

Для построения кривой суммарного крутящего момента М_i = М_крср многоцилиндрового двигателя графически суммируем кривые крутящих моментов от каждого цилиндра, сдвигая влево одну кривую относительно другой на угол и поворота кривошипа между вспышками.

Для двигателя с равными интервалами между вспышками суммарный крутящий момент для четырехтактного двигателя будет периодически повторятся через:

Поскольку , а R = const, то кривая будет отличаться от кривой лишь масштабом.

Масштаб крутящего момента ,

Таблица 5

Угол поворота коленчатого вала, град.	Крутящий момент для цилиндра, Н·м	Суммарный крутящий момент, Н·м
	№№ цилиндров
	1	5	4	2	6	3	7	8
0	0,0	106,6	0,0	-140,6	0,0	363,5	0,0	-114,1	215,4
30	-253,7	199,4	-122,2	55,8	450,6	318,3	-126,7	148,3	669,8
60	-153,2	121,8	-206,8	84,9	345,3	167,9	-207,0	252,2	405,1
90	106,6	0,0	-140,6	0,0	363,5	0,0	-114,1	0,0	215,4

Средний крутящий момент М_кр.ср. определяется по площади, лежащей под кривой графика суммарного М_кр.:

Н·м,

где F₁ и F₂ - соответственно положительная и отрицательная площади под кривой суммарного М_кр, мм₂;

L - длина интервала между вспышками по диаграмме крутящего момента, мм.

Найденный момент М_кр.ср. представляет собой средний индикаторный момент двигателя.

Эффективный крутящий момент двигателя:

Значение з_М определено в разделе 2.7. Значение эффективного крутящего момента, полученное по данной формуле, должно совпадать с величиной М_е, вычисленной ранее.

Нм

Погрешность вычислений составляет:

Относительная погрешность вычислений М_кр.ср не превышает

Заключение

В результате выполнения курсовой работы был произведен тепловой и динамический расчет двигателя М - 2106 - 70.

Проведя тепловой расчет, определили параметры рабочего тела в цилиндре двигателя, а также произвели оценочные показатели процесса, позволяющие определить параметры двигателя и оценить его мощностные и экономические показатели.

При выполнении динамического расчета определили действующие на кривошипно-шатунный механизм силы, а также крутящий момент, развиваемый двигателем.

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

1. Методические указания к курсовой работе по дисциплине «Автомобильные двигатели» для студентов специальности 1-37 01 06 - Новополоцк, 2001.-86 с.

2. Термодинамика, теплопередача и двигатели внутреннего сгорания/ Железко Б.Е., Адамов В.М., Есьман Р.И.- Мн.:«Выш.шк.»,1985.- 272с.

3. Техническая термодинамика/ Кириллин В.А., Сычев В.В., Шейндлин А.Е. - М.: Энергоатомиздат, 1983. - 416с.

4. Автомобили. Основы конструкций/ Вахламов В.К.- М.:Издательский центр «Академия», 2008.- 528с.

Размещено на Allbest.ru

...