Тепловой и динамический расчет двигателя внутреннего сгорания
Костанайский социально
- технический университет
имени академика З.
Алдамжар
ТЕХНИЧЕСКИЙ факультет
Кафедра ТРАНСПОРТ И
ТЕХНОЛОГИИ
Курсовая работа по
дисциплине
Энергетические установки транспортной техники
ТЕМА:
Тепловой и динамический расчет двигателя внутреннего сгорания
Выполнил:
Смагул
Жанболат
Дулатулы
курса
специальности
ТТТиТ,
заочной формы обучения
Научный
руководитель
Сатыбалдин
Т. Т.
ст.
преподаватель кафедры «ТиТ»
Костанай
2013
ЗАДАНИЕ
на курсовую работу студенту
Смагулу
Жанболату Дулатулы
Тема курсовой работы: Тепловой и динамический
расчет двигателя внутреннего сгорания
Целевая установка: Произвести расчеты
четырехтактного карбюраторного двигателя предназначенного для легковых
автомобилей (1-вариант). Эффективная мощность карбюраторного двигателя Ne=47,1
кВт при частоте вращения коленчатого вала nN=5600 мин-1.
Двигатель четырехцилиндровый, i=4
с рядным расположением. Система охлаждения жидкостная закрытого типа. Степень
сжатия ε=8,5
для карбюраторного двигателя.
С учетом приведенных рекомендаций и
заданий nN=5600 мин-1
тепловые расчеты последовательно проводятся для карбюраторного двигателя при
п=1000, 3200, 5600 и 6000 мин-1.
Объем курсовой работы - не
менее 25-30 стр.
Срок доклада руководителю о
ходе разработки курсовой работы:
а) доклад о собранном
материале и ходе разработки курсовой работы
до « »_____ 20___г.
б) доклад о ходе написания
курсовой работы до «____»___________ 20____г.
Срок сдачи курсовой работы -
«____» ________20_ г.
Руководитель курсовой работы:
______________ /Сатыбалдин Т.Т./
«___»_____________20___ г.
СОДЕРЖАНИЕ
Введение
.
Тепловой расчет
1.1 Топливо
1.2
Параметры рабочего тела
.3
Параметры окружающей среды и остаточные газы
.4
Процесс впуска
.5
Процесс сжатия
.6
Процесс сгорания
.7
Процессы расширения и выпуска
.8
Индикаторные параметры рабочего цикла
.9
Эффективные показатели двигателя
.10
Основные параметры цилиндра и двигателя
.11
Построение индикаторных диаграмм
.
Тепловой баланс
.
Расчет внешних скоростных характеристик
.
Кинематика
.
Динамика
.1
Силы давления газов
.2
Приведение масс частей кривошипно-шатунного механизма
.3
Удельные и полные силы инерции
.4
Удельные суммарные силы
.5
Крутящие моменты
.6
Силы, действующие на шатунную шейку коленчатого вала
.7
Силы, действующие на колено вала
.8
Равномерность крутящего момента и равномерность хода двигателя
Заключение
Список
использованной литературы
ВВЕДЕНИЕ
двигатель тепловой динамический газ
Карбюраторные двигатели находят широкое
применение в современной жизни. Их используют в основном на транспортных средствах
(из-за высокой стоимости топлива которые данные виды двигателей используют), к
таким транспортным средствам относятся: мотоциклы, автомобили, а также катера,
моторные лодки и т. п.
Используются карбюраторные двигатели и в
современной автомобильной промышленности.
Самый распространенный тип современного
теплового двигателя - двигатель внутреннего сгорания. Двигатели внутреннего
сгорания могут работать на жидком топливе (бензин, керосин и т. п.) или на
горючем газе, сохраняемом в сжатом виде в стальных баллонах или добываемом
сухой перегонкой из дерева (газогенераторные двигатели).
Особое место в конструкции автомобиля занимает
силовая установка, т. е. двигатель внутреннего сгорания (ДВС).
На автомобильном транспорте применяются
карбюраторные и дизельные двигатели, а также бензиновые двигатели с впрыском
топлива и принудительным воспламенением рабочей смеси.
Бензиновые двигатели с впрыском топлива и
принудительным воспламенением рабочей смеси в зависимости от организации
процесса смесеобразования и их конструктивных особенностей могут сочетать в
себе положительные свойства и карбюраторных двигателей и дизелей.
Актуальность исследования:
грамотная эксплуатация автомобиля предполагает не только знание его устройства
и технического обслуживания, но и понимание физических процессов, посредством
которых работает данное транспортное средство. Современный этап развития теории
автомобиля требует углубленное изучение отдельных его составляющих и
эксплуатационных свойств автомобиля, оптимизацией их показателей и технических
параметров, что позволяет еще на стадии проектирования создавать наиболее
рациональные конструкции автомобилей и обеспечить максимальную эффективность их
применения.
Объект исследования:
тепловой и динамический расчет двигателя внутреннего сгорания.
Предмет исследования:
бензиновый двигатель внутреннего сгорания.
Цель исследования:
произвести тепловой и динамический расчет двигателя внутреннего сгорания.
Задачи исследования:
изучение научно-методической и справочной литературы, применение теории на практике
расчета.
Методы исследования:
анализ научно-методической и справочной литературы, тепловой и динамический
расчет двигателя внутреннего сгорания.
1. ТЕПЛОВОЙ РАСЧЕТ
1.1 Топливо
В соответствии с заданной степенью
сжатия ε=8,5 можно
использовать бензины марок АИ-91, Регуляр-91 [1].
Средний элементарный состав и
молекулярная масса бензина: С=0,855; Н=0,145 и тт=114
Низшая теплота сгорания:
Ни=33,91С+125,60Н-10,89(О-S)-2,51(9H+W)=
=33,91·0,855+125,60·0,145-2,15·9·0,145=43930
кДж (1.1)
1.2 Параметры рабочего тела
Теоретически необходимое количество
воздуха для сгорания 1 кг топлива
L0==0,517 кмоль
возд/кг топ.
l0==14,957 кг
возд/кг топл.
Коэффициент избытка воздуха
устанавливается на основании следующих соображений. На современных двигателях
устанавливают многокамерные карбюраторы, обеспечивающие получение почти
идеального состава смеси по скоростной характеристике. Возможность применения
для рассчитываемого двигателя двухкамерного карбюратора с обогатительной
системой и системой холостого хода позволяет получить при соответствующей
регулировке как мощностной, так и экономичный состав смеси. Стремление получить
двигатель достаточно экономичный и с меньшей токсичностью продуктов сгорания,
которая достигается при α=0,95-0,98, позволяет принять α=0,96 на
основных режимах, а на режимах минимальной частоты вращения α=0,86. Далее
непосредственный числовой расчет будет проводиться только для режимов
максимальной мощности, а для остальных режимов окончательные значения
рассчитываемых параметров приводятся в табличной форме [2].
Количество горючей смеси для
карбюраторного двигателя:
М1=α·L0+1/m=0,96·0,517+1/114=0,5051
кмоль гор.см/кг топл.;
Количество отдельных компонентов
продуктов сгорания
при К=0,5 и принятых скоростных режимах [3]:
=0,0770 кмоль СО2/кг
топл.
МСО=2··0,208·L0=2··0,208·0,517=0,0057
кмоль СО/кг топл.
МНО=-2К··0,208·L0=-2·0,5··0,208·0,517=0,0696
кмольН2О/кг топл.
МН=2К··0,208·L0=2·0,5··0,208·0,517=0,0029
кмольН2 /кг топл.
МN=0,792·α·L0=0,792·0,96·0,517=0,3931
кмоль N2/кг топл.
Общее количество продуктов сгорания
для карбюраторного двигателя:
М2=МСО+Мсо+МНО+МН+МN=
0,077+0,0057+0,0696+0,0029+0,3931=0,5483
кмоль пр.сг/кг топл.
Результаты расчетов при п=1000,
3200, 5600 и 6000 мин-1 сведены в таблицу 1.
Таблица 1
Параметры рабочего тела
Параметры
|
Рабочее
тело
|
n
|
1000320056006000
|
|
|
|
α
|
0,860,960,960,96
|
|
|
|
M1
|
0,4530,50510,50510,5051
|
|
|
|
МСО2
|
0,0910,0770,0770,077
|
|
|
|
MCO
|
0,02010,00570,00570,0057
|
|
|
|
МН2О
|
0,06250,06960,06960,0696
|
|
|
|
МН2
|
0,01000,00290,00290,0029
|
|
|
|
МN2
|
0,35210,39310,39310,3931
|
|
|
|
M2
|
0,53570,54830,54830,5483
|
|
|
|
1.3 Параметры окружающей среды и
остаточные газы
Давление и температура окружающей
среды при работе двигателей без наддува: рk=р0=0,1
МПа и Tk=T0=293 К
Температура остаточных газов.
При постоянных значениях степени сжатия ε=8,5
температура остаточных газов практически линейно возрастет с увеличением
скоростного режима при α-const,
но уменьшается при обогащении смеси. Учитывая уже определенные значения п
и α,
можно принять значения Тr
для расчетных режимов карбюраторного двигателя [4]. При номинальных режимах для
карбюраторного двигателя Тr=1060
К.
Давление остаточных газов
рr
за счет расширения фаз газораспределения и снижения сопротивлений при
конструктивном оформлении выпускных трактов рассчитываемых двигателей можно
принять на номинальном скоростном режиме для карбюраторного двигателя:
рrN=1,18·p0=1,18·0,1=0,118
МПа;
Тогда величины давлений на остальных
режимах работы двигателей можно подсчитать по формулам:
рr=р0(1,035+Ар·10-8·n2)=0,1(1,035+0,4624·10-8·5600)=0,1180,
где Ар=(рrN -p0·1,035)·108/(nN2 p0)
При nN=5600 мин-1
Ар=(0,118-0,1·1,035)·108/(56002·0,1)=0,4624
1.4 Процесс впуска
Температура подогрева свежего заряда.
С целью получения хорошего
наполнения двигателей на номинальных скоростных режимах принимается для
карбюраторного двигателя ΔTN=8°С.
Тогда на остальных режимах значения ΔT
рассчитываются по формуле:
ΔT=AT(110-0,0125n),
где AT=ΔTN
/(110-0,0125·nN)
ΔT=ΔTN
/(110-0,0125·nN) (110-0,0125n)=8 (1.4.1)
Плотность заряда на впуске:
ρо=ро·106/(RвТо)=0,1·106/(287·293)=1,189
кг/м3
где Rв=287
Дж/(кг·град)-удельная газовая постоянная для воздуха
Потери давления на впуске. В
соответствии со скоростным режимом n=5600 мин-1 и при учете
качественной обработки внутренних поверхностей впускных систем можно принять
для карбюраторного двигателя β2+ξвп=2,8 и ωвп=95 м/с.
Тогда ΔPa на всех
скоростных режимах двигателей рассчитывается по формуле: ΔPa=(β2+ξвп)An2·n2·pk·10-6/2,
где An=ωвп/nN
Потери давления на впуске
карбюраторного двигателя при n=5600 мин-1,
Аn=95/5600=0,01696,
тогда ΔPa=2,8 0,016962·56002·1,189·10-6/2=0,0150
МПа.
Давление в конце впуска в карбюраторном
двигателе при nN=5600 мин-1
а=P0-ΔPa=0,1-0,015=0,0850
МПа;
Коэффициент остаточных газов. При
определении γr для
карбюраторного двигателя без наддува принимается коэффициент очистки φоч=1, а
коэффициент дозарядки на номинальном скоростном режиме, φдoз=1,10 что
вполне возможно получить при подборе угла опаздывания закрытия впускного
клапана в пределах 30-60°. При этом на минимальном скоростном режиме (n=1000
мин-1) возможен обратный выброс в пределах 5%, т. е. φдоз=0,95. На
остальных режимах значения φдоз можно
получить, приняв линейную зависимость φдоз от
скоростного режима [5]
Тогда при nN=5600
мин-1:
γr=·==0,0495(1.4.2)
Температура в конце впуска:
Та=(Т0+ΔТ+γrТr)/(1+γr)=(293+8+0,0495·1060)/(1+0,0495)=337
К (1.4.3)
Коэффициент наполнения:
ηv=(φдоз·ε·ра-φоч·рr)=
=(1,1·8,5·0,085-1·0,118)=0,8784
Результаты расчетов при п=1000,
3200, 5600 и 6000 мин-1 сведены в таблицу 2.
Таблица 2
Параметры процесса впуска газообмена
Параметры
|
Процесс
впуска газообмена
|
n
|
1000320056006000
|
|
|
|
α
|
0,860,960,960,96
|
|
|
|
Tr
|
900
|
1000
|
1060
|
1070
|
pr
|
0,10400,10820,11800,1201
|
|
|
|
ΔT
|
19,514,08,07,0
|
|
|
|
ΔРa
|
0,00050,00490,01500,0172
|
|
|
|
Рa
|
0,09950,09510,08500,0828
|
|
|
|
φдоз
|
0,950
|
1,025
|
1,10
|
1,12
|
γr
|
0,05160,04610,04950,0504
|
|
|
|
Ta
|
341338337337
|
|
|
|
ηυ
|
0,87440,91670,87840,8701
|
|
|
|
1.5 Процесс сжатия
Средний показатель адиабаты сжатия k1
при ε=8,5,а также
рассчитанных значениях Та определяется по номограмме, а
средний показатель политропы сжатия n1 принимается несколько
меньше k1.
При выборе n1
учитывается, что с уменьшением частоты вращения теплоотдача от газов в стенки
цилиндра увеличивается, а n1 уменьшается по сравнению с k1
более значительно для карбюраторного двигателя при nN=5600
мин-1, Ta=337 К и ε=8,5
показатель адиабаты сжатия определен по номограмме k1=1,3775.
Давление в конце сжатия для
карбюраторного двигателя при nN=5600 мин-1: Pс=Pаεn1=0,085·8,5
1,377=1,6189 MПа, где n1=1,377 принят несколько
меньше k1=1,3775.
Температура в конце сжатия: Tc=Taεn-1=337·8,51,377-1=755
К;
Средняя мольная теплоемкость:
а) свежей смеси (воздуха):
(mcv)=20,6+2,638
10-3 tc, где tc=Tc-273°С=755-273=482;
(mcv)=20,6+2,638
10-3 482=21,872 кДж/(кмоль·град) (1.4.5)
б) остаточных газов (mc"v)
определяется методом интерполяции по табл.3.8
для карбюраторного двигателя при nN=5600
мин-1, α=0,96 и tc=482°С
(mc"v)=23,586+(23,712-23,586)·=23,611 (1.4.6)
где 23,586=ср.тепл.при α-0,95(t-400) и
23,712=ср.тепл.при α-1,00(t-400)-значения
теплоемкости продуктов сгорания при 400°С соответственно при α=0,95 и α=1,00,
взятые по табл.3.8;
(mc"v)=24,014+(24,15-24,014)·=24,041
где 24,014=ср.тепл.при α-0,95(t-500) и
24,150=ср.тепл.при α-1,00(t-500)-значения
теплоемкости продуктов сгорания при 500°С соответственно при α=0,95 и α=1,0, взятые
по табл. 3.8.
Теплоемкость продуктов сгорания при tc=482°С
и α=0,96
(mc"v)=23,611+(24,041-23,611)=23,964
кДж/(кмоль·град);
в) рабочей смеси (mc'v)=[(mcv)+γr(mc"v)]:
(mc'v)=[21,872+0,0495·23,964]=21,971
кДж/(кмоль·град);
Результаты расчетов при п=1000,
3200, 5600 и 6000 мин-1 сведены в таблицу 3.
Таблица 3
Параметры процесса сжатия
Параметры
|
Процесс
сжатия
|
n
|
1000320056006000
|
|
|
|
k1
|
1,3770
|
1,3778
|
1,3775
|
1,3775
|
n1
|
1,370
|
1,377
|
1,377
|
1,377
|
Рc
|
1,86691,81131,61891,5770
|
|
|
|
Tc
|
753757755755
|
|
|
|
tc
|
480484482482
|
|
|
|
(mcv)21,86621,87721,87221,872
|
|
|
|
|
(mc"v)23,95523,97223,96423,964
|
|
|
|
|
(mc'v)21,96921,96921,97121,972
|
|
|
|
|
1.6 Процесс сгорания
Коэффициент молекулярного изменения
горючей μ0=М2/М1
и рабочей смеси μ=(μ0+γr)/(1+γr) для
карбюраторного двигателя при nN=5600 мин-1: μ0=0,5483/0,5051=1,0855
и μ=(1,0855+0,0495)/(1+0,0495)=1,0815
количество теплоты, потерянное
вследствие химической неполноты сгорания, и теплота сгорания рабочей смеси
ΔHи=119950(1-α)L0 и Нраб.см= (1.6.1)
Для карбюраторного двигателя при nN=5600
мин-1
ΔHи=119950(1-0,96)·0,517=2481
кДж/кг,
Нраб.см==78191
кДж/кмоль раб.см.;
Средняя мольная теплоемкость
продуктов сгорания
(mC"v)=(1.6.2)
определяется по эмпирическим
формулам, приведенным в табл.3.6 для интервала температур от 1501 до 2800°С
=22,490+0,001430·tZ
=39,123+0,003343·tZ
=19,678+0,001758·tZ
=26,670+0,004438·tZ
=23,723+0,001550·tZ
=21,951+0,001457·tZ
(mC"v)=(1/0,5483)·[0,077·(39,123+0,003343·tz)+0,0057·(22,490+0,001430·tZ)+0,0696·(26,670+0,004438·tZ)+0,0029·(19,678+0,001758·tZ)+0,3931·(21,951+0,001457·tZ)]=1,8238[3,012+0,00026tZ+0,128+0,00008tZ+1,85623+0,00031tZ+0,057066+0,000005tZ+8,62894+0,00057tZ]=5,4933+0,0004742tZ+0,23345+0,0001459tZ+3,3853923+0,0005654tZ+0,10407697+0,0000091tZ+15,7374608+0,0001040tZ=24,954+0,0012986tZ
кДж/(кмоль·град).
Коэффициент использования теплоты ξz зависит от
совершенства организации процессов смесеобразования и сгорания топлива. Он
повышается за счет снижения потерь теплоты газов в стенки цилиндра и
неплотности между поршнем и цилиндром. При увеличении скоростного режима ξz снижается.
При проведении расчетов двигателя ξz выбирается
по опытным данным в зависимости от конструктивных особенностей двигателя. В
зависимости от реальной зависимости ξz от
скоростного режима карбюраторного двигателя принята величина коэффициента
использования теплоты для карбюраторного двигателя ξz=0,91 при nN=5600
мин-1. Так же определены значения ξz для всех
расчетных режимов.
Температура в конце видимого
процесса сгорания:
ξzНраб.см+(mc'v)·tc=μ(mC"v)·tZ (1.6.3)
Для карбюраторного двигателя при nN=5600
мин-1
0,91·78191+21,971·482=1,0815·(24,954+0,0012986tZ)tZ или 0,001405tZ2+26,993tZ-81757=0
откуда tZ=(-26,993+)/(2·0,001405)=2660°С;
Тz=tZ+273=2660+273=2933
К
Максимальное давление сгорания
теоретическое:
z=pc·μ·Tz/Tc=1,6189·1,0815·2933/755=6,8016
МПа;
Максимальное давление сгорания
действительное:
pzд=0,85·pz=0,85·6,8016=5,7814
МПа;
Степень повышения давления: λ=pz/pc=6,8016/1,6189=4,201;
Результаты расчетов при п=1000,
3200, 5600 и 6000 мин-1 сведены в таблицу 4.
Таблица 4
Параметры процесса сгорания
Параметры
|
Процесс
сгорания
|
n
|
1000320056006000
|
|
|
|
μ0
|
1,18261,08551,08551,0855
|
|
|
|
μ
|
1,17361,08171,08151,0814
|
|
|
|
ΔHи
|
8682248124812481
|
|
|
|
Нраб.см
|
73992784457819178124
|
|
|
|
(mC"v)25,291+0,0012675tZ24,925+0,0021065tZ24,954+0,0012986tZ24,954+0,0021065tZ
|
|
|
|
|
ξz
|
0,81
|
0,91
|
0,90
|
0,88
|
tZºC
|
2129250226602129
|
|
|
|
Тz
К
|
2402277529332402
|
|
|
|
pz
|
7,04277,20756,80165,4461
|
|
|
|
pzд
|
5,98636,12645,78144,6292
|
|
|
|
λ
|
3,7723,9794,2013,453
|
|
|
|
1000(mC"v)=(1/0,5357)·[0,091·(39,123+0,003343·tz)+0,0201·(22,490+0,001430·tZ)+0,0625·(26,670+0,004438·tZ)+0,0100·(19,678+0,001758·tZ)+0,3521·(21,951+0,001457·tZ)]=1,8667[3,560+0,00030tZ+0,396+0,00003tZ+1,66688+0,00028tZ+0,196780+0,000018tZ+7,72895+0,00051tZ]=6,6455+0,0005600tZ+0,73921+0,0000560tZ+3,1115649+0,0005227tZ+0,36732923+0,0000336tZ+14,4276310+0,0000952tZ=25,291+0,0012675tZ
кДж/(кмоль·град).
(mC"v)=(1/0,5483)·[0,077·(39,123+0,003343·tz)+0,0057·(22,490+0,001430·tZ)+0,0696·(26,670+0,004438·tZ)+0,0029·(19,678+0,001758·tZ)+0,3931·(21,951+0,001457·tZ)]=1,8238[3,012+0,00026tZ+0,112+0,00001tZ+1,85623+0,00031tZ+0,057066+0,000005tZ+8,62894+0,00057tZ]=5,4933+0,0004742tZ+0,20427+0,0000182tZ+3,3853923+0,0005654tZ+0,10407697+0,0000091tZ+15,7374608+0,0010396tZ=24,925+0,0021065tZ
кДж/(кмоль·град).
(mC"v)=(1/0,5483)·[0,077·(39,123+0,003343·tz)+0,0057·(22,490+0,001430·tZ)+0,0696·(26,670+0,004438·tZ)+0,0029·(19,678+0,001758·tZ)+0,3931·(21,951+0,001457·tZ)]=1,8238[3,012+0,00026tZ+0,128+0,00001tZ+1,85623+0,00031tZ+0,057066+0,000005tZ+8,62894+0,00057tZ]=5,4933+0,0004742tZ+0,23345+0,0000182tZ+3,3853923+0,0005654tZ+0,10407697+0,0000091tZ+15,7374608+0,0010396tZ=24,954+0,0021065tZ
кДж/(кмоль·град).
1000 ξzНраб.см+(mc'v)·tc=μ(mC"v)·tZ
0,81·73992+21,969·480=1,1826·(25,291+0,0012675tZ)tZ
или 0,001499tZ2+29,909tZ-70479=0
откуда tZ=(-29,909+√29,9092+4·0,001499·70479)/(2·0,001499)=2129°С
3200 ξzНраб.см+(mc'v)·tc=μ(mC"v)·tZ
0,91·78445+21,969·484=1,0855·(24,925+0,0021065tZ)tZ
или 0,002287tZ2+27,056tZ-82018=0
откуда tZ=(-27,056+√27,0562+4·0,002287·82018)/(2·0,002287)=2502°С
6000 ξzНраб.см+(mc'v)·tc=μ(mC"v)·tZ
0,88·78124+21,972·482=1,0855·(24,954+0,0021065tZ)tZ
или 0,001499tZ2+29,909tZ-70479=0
откуда tZ=(-29,909+√29,9092+4·0,001499·70479)/(2·0,001499)=2129°С
1.7 Процессы расширения и выпуска
Средний показатель адиабаты
расширения k2
определяется по номограмме при заданном ε для соответствующих значений
α и Tz,
а средний показатель политропы расширения п2 оценивается по
величине среднего показателя адиабаты для карбюраторного двигателя при ε=8,5, α=0,96 и Tz=2933 К k2=1,25,
что позволяет принять п2=1,251.
Давление и температура в конце
процесса расширения: pb=pz/εп2 и Тb=Tz/εn2-1
Для карбюраторного двигателя при nN=5600
мин-1: pb=6,8016/8,51,251=0,4676
МПа и Тb=2933/8,51,251-1=1714
К;
Проверка ранее принятой температуры
остаточных газов: Тr=;
При nN=5600 мин-1:
Тr==1088 К,
Тr=1714/(0,4676/0,118)^0,33=1088
ΔTr=100(1088-1060)/1060=2,64
%;
где ΔTr-погрешность
расчета.
На всех скоростных режимах
температура остаточных газов принята в начале расчета достаточно удачно, так
как ошибка не превышает 6 %.
Результаты расчетов при п=1000,
3200, 5600 и 6000 мин-1 сведены в таблицу 5.
Таблица 5
Параметры процесса расширения и
выпуска
Параметры
|
Процесс
расширения и выпуска
|
n
|
1000320056006000
|
|
|
|
k2
|
1,2625
|
1,2525
|
1,251,2582
|
|
n2
|
1,262
|
1,252
|
1,251
|
1,258
|
0,47290,48320,46760,3689
|
|
|
|
Tb
|
1371161817141383
|
|
|
|
Tr
|
86710161088949
|
|
|
|
ΔTr, %
|
-3,671,602,64-10,47
|
|
|
|
1.8 Индикаторные параметры рабочего
цикла
Теоретическое среднее индикаторное давление:
р= (1.8.1)
для карбюраторного двигателя при nN=5600
мин-1
р==1,4835 МПа
р=(1,6189/8,5-1)[4,201/1,251-1(1-8,51-1,251)-1/1,377-1(1-8,51-1,377)]=1,4835
Среднее индикаторное давление: рi=φдоз р
карбюраторного двигателя pi=0,95·1,4835=1,4093
МПа.
Индикаторный КПД: ηi=;
Индикаторный удельный расход
топлива: gi =;
карбюраторного двигателя при nN
=5600 мин-1
ηi=1,4093·14,957·0,96/43930·10-3·1,189·0,8784=0,4410;
gi =3600/0,441·43930·10-3=186
г/(кВт·ч);
1.9 Эффективные показатели двигателя
Среднее давление механических потерь
для бензиновых двигателей с числом цилиндров до шести и отношением S/D≤1:
рм=0,034+0,0113vп.ср.
Для карбюраторного двигателя,
предварительно приняв ход поршня S равным 78 мм, получим значение
средней скорости поршня при nN=5600 мин-1
vп.ср=S·n /(104·3)=78·5600/(104·3)=14,56
м/с
Тогда рм =0,034+0,0113·14,56=0,1985
МПа
Среднее эффективное давление: ре
=рi-рм =1,4093-0,1985=1,2108 МПа
Механический КПД: ηм=ре /рi
карбюраторного двигателя: ηм=1,2108/1,4093=0,8591
Эффективный КПД: ηе=ηм·ηi
Эффективный удельный расход топлива:
ge=
карбюраторного двигателя: ηe=0,8591·0,441=0,3789
и ge=3600/43,93·0,379=311 г/(кВтч);
Результаты расчетов при п=1000,
3200, 5600 и 6000 мин-1 сведены в таблицу 6.
Таблица 6
Индикаторные и эффективные параметры
двигателей
Параметры
|
Индикаторные
и эффективные параметры двигателей
|
n
|
1000320056006000
|
|
|
|
|
1,48831,48351,48351,4835
|
|
|
|
рi
|
1,41391,40931,40931,4093
|
|
|
|
ηi
|
0,44250,44100,44100,4410
|
|
|
|
gi
|
185,2185,8185,8185,8
|
|
|
|
vп.ср
|
2,608,3214,5615,60
|
|
|
|
рм
|
0,06340,12800,19850,2103
|
|
|
|
ре
|
1,35051,28131,21081,199
|
|
|
|
ηм
|
0,95520,90920,85910,8508
|
|
|
|
ηe
|
0,42270,4010,37890,3752
|
|
|
|
ge
|
346329311307
|
|
|
|
1.10 Основные параметры цилиндра и
двигателя
Литраж
карбюраторного двигателя: Vл=, (1.10.1)
где τ-тактность (τ=4)
Vл=30·47,1·4/1,2108·5600=0,83357
л;
Рабочий объем одного
цилиндра карбюраторного двигателя: Vh=
Vh=0,83357/4=0,20839
л;
Диаметр цилиндра. Так как
ход поршня предварительно был принят S=78 мм, то: D=2·103·=2·103·0,20839/(3,14·78)=58,34
мм;
Окончательно принимается для
карбюраторного двигателя D=65 мм и S=65 мм.
Основные параметры и показатели
двигателей определяются по окончательно принятым значениям D и S:
-площадь поршня: Fп=π·D2/(4·100)=3,14·652/(4/100)=33,17
см2;
-литраж двигателя: Vл==3,14·652·65·4/(4·106)=0,86
л;
-мощность двигателя: Ne==1,2108·0,86·5600/30·4=48,59
кВт;
-литровая мощность двигателя: Nл=Ne /Vл=54,767 л;
-крутящий момент: Ме==·47,1/5600=82,9
Н·м;
-часовой расход топлива: GT=Ne·ge·10-3=48,59·311·10-3=15,11
кг/ч
Результаты расчетов при п=1000,
3200, 5600 и 6000 мин-1 сведены в таблицу 7.
Таблица 7
Основные параметры и показатели
двигателя
Параметры
|
Основные
параметры и показатели двигателя
|
n, мин
|
1000320056006000
|
|
|
|
Fп
|
33,17
|
Vл
|
0,86
|
Nл
|
54,767
|
Ne
|
8,6827,7748,5952,06
|
|
|
|
Ме
|
82,9382,9182,982,9
|
|
|
|
GT
|
3,009,1415,1115,98
|
|
|
|
1.11 Построение индикаторных
диаграмм
Индикаторная диаграмма
карбюраторного двигателя построена для номинального режима работы двигателя, т.
е. при Ne=47,1 кВт и n=5600 мин-1,
аналитическим методом.
Масштабы диаграммы: масштаб хода
поршня Ms=l мм в мм;
Масштаб давлений Мр=0,05
МПа в мм.
Величины в приведенном масштабе,
соответствующие рабочему объему цилиндра и объему камеры сгорания:
АВ=S/Ms=65/1,0=65
мм; ОА=АВ/(ε-1)=65/(8,5-1)=8,67 мм.
Максимальная высота диаграммы (точка
z): Pz/Mp=6,8016/0,05=136,0 мм.
Ординаты характерных точек: рa/Mp=0,085/0,05=1,7
мм;
pc/Mp=1,6189/0,05=32,4
мм
рb/Mp=0,4676/0,05=9,4
мм;
pr/Mp=0,118/0,05=2,4
мм
р0/Mp=0,1/0,05=2
мм
Построение политроп сжатия и
расширения аналитическим методом:
а) политропа сжатия px=pa(Va/Vx)n
Отсюда рx/Mp=(pа/Mp)(OB/OX)n1=1,7(73,7/OX)1,377
мм;
где OB=OA+AB=8,67+65=73,7
мм;
б) политропа расширения рх=рb(Vb/Vx)n2
Отсюда рx/Mp=(pb/Mp)(OB/OX)n2=9,4(73,7/OX)1,251
мм;
Результаты расчета точек политроп
приведены в таблице 8.
Таблица 8
Координаты точек политроп
№
точек
|
OX,
мм
|
OB/OX
|
Политропа
сжатия
|
Политропа
расширения
|
|
|
|
(OB/OX)1,377
|
рx/Mp, мм
|
рх, МПа
|
(OB/OX)1,251
|
рx/Mp, мм
|
рх, МПа
|
1
|
10,4
|
7,0914,8425,21,26 (точка с)11,59108,95,45
(точка z)
|
|
|
|
|
|
|
2
|
11,0
|
6,713,7223,31,1710,80101,55,08
|
|
|
|
|
|
|
3
|
12,6
|
5,811,2519,10,969,0284,84,24
|
|
|
|
|
|
|
4
|
17,7
|
4,27,2112,30,626,0256,62,83
|
|
|
|
|
|
|
5
|
22,1
|
3,35,188,80,444,4541,82,09
|
|
|
|
|
|
|
6
|
29,5
|
2,53,536,00,33,1529,61,48
|
|
|
|
|
|
|
7
|
44,2
|
1,72,083,50,181,9418,20,91
|
|
|
|
|
|
|
8
|
58,9
|
1,31,442,40,121,3913,10,66
|
|
|
|
|
|
|
9
|
88,4
|
0,80,741,30,07 (точка а)0,767,1
(точка b)0,36
|
|
|
|
|
|
|
Теоретическое среднее индикаторное давление:
pi'=F1Mp/AB=1929·0,05/65=1,484
МПа,
где F1=1929 мм2-площадь
диаграммы aczba
Величина pi'=1,484
МПа, полученная планиметрированием индикаторной диаграммы, очень близка к
величине pi'=1,4835 МПа, полученной в тепловом расчете.
Скругление индикаторной диаграммы
осуществляется на основании следующих соображений и расчетов. Так как
рассчитываемый двигатель достаточно быстроходный (n=5600 мин-1),
то фазы газораспределения необходимо устанавливать с учетом получения хорошей
очистки цилиндра от отработавших газов и обеспечения дозарядки в пределах,
принятых в расчете. В связи с этим начало открытия впускного клапана (точка r')
устанавливается за 18° до прихода поршня в в.м.т., а закрытие (точка а')
- через 60° после прохода поршнем н.м.т.; начало открытия выпускного клапана
(точка b') принимается за 55° до прихода поршня в н.м.т., а закрытие
(точка а')-через 25° после прохода поршнем в.м.т. Учитывая
быстроходность двигателя, угол опережения зажигания в принимается равным 35°, а
продолжительность периода задержки воспламенения Δφ1=5°.
В соответствии с принятыми фазами
газораспределения и углом опережения зажигания определяют положение точек r', а',
а", с', f и b' по формуле для перемещения
поршня: АХ=[(1-cos φ)+(1-cos 2φ)], где λ-отношение
радиуса кривошипа к длине шатуна.
Выбор величины λ
производится при проведении динамического расчета, а при построении
индикаторной диаграммы предварительно принимается λ=0,285.
Расчеты ординат точек r', a',
а", с', и b' сведены в таблице 9.
Таблица 9
Координаты точек политроп
Обозначение
точек Положение точек φº (1-cos
φ)+(1-cos 2φ)Расстояние
точек
от
в.м.т. (АХ), мм
|
|
|
r'
|
18º
до
в.м.т.
|
18
|
0,0655
|
2,1
|
a'
|
25º
после
в.м.т.
|
25
|
0,1223
|
4,0
|
а"
|
60º
после
в.м.т.
|
120
|
1,6069
|
19,7
|
с'
|
35º
до
в.м.т.
|
35
|
0,2313
|
7,5
|
f
|
30º
до
в.м.т.
|
30
|
0,1697
|
5,5
|
b'
|
55º
до
в.м.т.
|
125
|
1,6667
|
21,7
|
Положение точки с" определяется из
выражения
рс" =(1,15÷1,25) рс=1,25·1,6189=2,024
МПа;
р"с=2,024/0,05=40,5
мм.
Действительное давление сгорания: pzд=0,85pz=0,85·6,8016=5,7814
МПа;
zд/Мр=5,7814/0,05=116
мм.
Нарастание давления от точки с"
до zд составляет 5,7814-2,024=3,757 МПа или 3,757/12=0,34
МПа/град п.к.в., где 12° положение точки zд по горизонтали
(для упрощения дальнейших расчетов можно принять, что действительное
максимальное давление сгорания pzд достигается
через-10° после в.м.т., т. е. при повороте коленчатого вала на 370°).
Соединяя плавными кривыми точки r, с, а',
с', с, с" и далее с zд и кривой
расширения b' с b" (точка b" располагается
обычно между точками b и а) и линией выпуска b"r'r,
получим скругленную действительную индикаторную диаграмму ra'ac'fc"zдb'b"r.
Координаты точек действительной
индикаторной диаграммы сведены в Таблицу 10. Построена индикаторная диаграмма
по данным точкам (в соответствии с рисунком рис. 1).
Таблица 10
Координаты точек действительной
индикаторной диаграммы
ex
|
pr
|
pa
|
pсж
|
pz
|
pрас
|
10,1180,0851,6191,619
|
|
|
|
|
|
10,1180,0851,6196,801
|
|
|
|
|
|
1,0630,1180,0851,4885,7815,781
|
|
|
|
|
|
1,1890,1180,0851,2765,781
|
|
|
|
|
|
1,7000,1180,0850,7803,696
|
|
|
|
|
|
2,8330,1180,0850,3861,951
|
|
|
|
|
|
4,250,1180,0850,2211,175
|
|
|
|
|
|
8,5000,1180,0850,0850,494
|
|
|
|
|
|
Рисунок 1. Индикаторная диаграмма
карбюраторного двигателя.
2. ТЕПЛОВОЙ БАЛАНС
В реальном двигателе потери тепла
возрастают из-за трения, теплообмена, неполноты сгорания и других причин. В
связи с этим эффективный КПД ηe цикла имеет
меньшее значение по сравнению с величиной ηt
Распределение тепловой энергии
топлива, сгорающего в двигателе, наглядно иллюстрируется составляющими внешнего
теплового баланса, которые определяются при установившемся тепловом состоянии
двигателя в процессе его испытаний. Приближенно составляющие теплового баланса
можно найти аналитически по данным теплового расчета двигателя.
Тепловой баланс позволяет определить
тепло, превращенное в полезную эффективную работу, т. е. установить степень
достигнутого совершенства теплоиспользования и наметить пути уменьшения
имевшихся потерь. Знание отдельных составляющих теплового баланса позволяет
судить о теплонапряженности деталей двигателя, рассчитать схему охлаждения,
выяснить возможность использования теплоты отработавших газов и т. д.
В общем виде внешний тепловой баланс
двигателя может быть представлен в виде следующих составляющих:
Qо=Qe+Qr+Qв+Qн.с+Qост=НиGт/3,6 (2.1.1)
где Q0-общее
количество теплоты, введенной в двигатель с топливом;
теплота, эквивалентная эффективной
работе двигателя за 1 с: Qe=1000Nе;
теплота, потерянная с отработавшими
газами:
Qr=(GТ /3,6)·{M2 [(mc"V)+8,315]·tr-M1[(mcV)+8,315]·to} (2.1.2)
Теплота, передаваемая охлаждающей
среде:
Qв=c·i··nm·(Hu -ΔHu)/(α·Hu) (2.1.3)
где с=0,45-0,53-коэффициент
пропорциональности для четырехтактных двигателей. В расчетах принято с=0,5;
i-число цилиндров; D-диаметр цилиндра, см; n-частота
вращения коленчатого вала двигателя, мин-1; m=0,5-0,7-показатель
степени для четырехтактных двигателей. В расчетах принято для карбюраторного
двигателя при n=1000 мин-1 m=1,6, а на всех остальных
скоростных режимах m=0,65.
Теплота, потерянная из-за химической
неполноты сгорания топлива:
Qн.с.=ΔНи·GТ / 3,6; (2.1.4)
неучтенные потери теплоты: Qост.=Qо-(Qе+Qr+Qв+Qн.с.)
Знание абсолютных значений
составляющих теплового баланса позволяет осуществить количественную оценку
распределения теплоты в двигателе. Если же необходимо сравнить распределение
теплоты в различных двигателях или оценить степень теплоиспользования
конкретного двигателя, то составляющие теплового баланса удобнее представлять в
относительных величинах, например, в процентах по отношению ко всей теплоте,
подведенной с топливом: qо=qe+qr+qв+qн.с+qост=100%
Величины отдельных составляющих
теплового баланса двигателя не являются постоянными, а изменяются в процессе
его работы в зависимости от нагрузки, быстроходности и других факторов.
Характер распределения теплоты,
подводимой в цилиндр с топливом, в процессе превращения в полезную эффективную
работу наглядно может быть представлен в виде кривых теплового баланса.
Графические зависимости строятся на основании определения каждой составляющей в
зависимости от частоты вращения, нагрузки, качества смеси и т. д. Необходимые
для построения указанных кривых теплового баланса данные получают при
проведении специальных испытаний двигателя либо путем использования результатов
ранее выполненных экспериментов. Тепловой баланс может быть также построен по
данным теплового расчета двигателя с использованием формул.
Общее количество теплоты, введенной
в двигатели при номинальном скоростном режиме (все данные взяты из теплового
расчета) карбюраторного двигателя: Qo=43930·15,11/3,6=184384
Дж/с;
Теплота, эквивалентная эффективной
работе за 1 с карбюраторного двигателя: Qe=1000·48,59=48590
Дж/с;
Теплота, передаваемая охлаждающей
среде карбюраторным двигателем:
B=0,5·4·6,51+2·0,65·56000,65·(43930-2083,9)/(0,96·43930)=79348
Дж/с;
Теплота, унесенная с отработавшими
газами карбюраторного двигателя:
r=(15,11/3,6){0,5483[25,25+8,315]·815-0,5051[20,775+8,315]20}=50979
Дж/с;
где (mc"V)=25,250
кДж/(кмоль·град)-теплоемкость отработавших газов при α=0,96 и tr=Tr-273=1088-273=815°С);
(mcV)=20,775
кДж/(кмоль·град)-теплоемкость свежего заряда;
Теплота, потерянная из-за химической
неполноты сгорания топлива карбюраторного двигателя: Qн.с.=2481·15,11/3,6=10413
Дж/с;
Неучтенные потери теплоты
карбюраторного двигателя:
ост=184384-(48590+50979+79348+10413)=-4946
Дж/с; (2.1.5)
Составляющие тепловых балансов
карбюраторного двигателя представлены в таблице 11.
Таблица 11
Составляющие теплового баланса
Составляющие
теплового баланса
|
Частота
вращения двигателя, мин-1
|
|
1000320056006000
|
|
|
|
|
Q, Дж/с
|
q, %
|
Q, Дж/с
|
q, %
|
Q, Дж/с
|
q, %
|
Q, Дж/с
|
q, %
|
Теплота,
эквивалентная эффективной работе
|
868023,72777024,94859026,45206026,7
|
|
|
|
|
|
|
|
Теплота,
передаваемая охлаждающей среде
|
2589570,75515249,47934843,08298742,6
|
|
|
|
|
|
|
|
Теплота,
унесенная с отработавшими газами
|
48101131,44812843,25097927,64811524,7
|
|
|
|
|
|
|
|
Теплота,
потерянная из-за химической неполноты сгорания топлива
|
20685,662995,6104135,6110135,6
|
|
|
|
|
|
|
|
Неучтенные
потери теплоты
|
-48136-131,5-25816-23,1-4946-2,78250,4
|
|
|
|
|
|
|
|
Общее
количество теплоты, введенной в двигатель с топливом
|
36608100111533100184384100195000100
|
|
|
|
|
|
|
|
Из приведенных таблиц и рисунков видно, что
основная часть теплоты топлива расходуется на эффективную работу, нагрев
охлаждающей среды и потери с отработавшими газами. На величину теплоты,
отводимой охлаждающей средой, оказывают влияние многие эксплуатационные и конструктивные
факторы. С увеличением частоты вращения двигателя и температуры охладителя, а
также коэффициента избытка воздуха величина Qв уменьшается, а
с увеличением размеров охлаждающей поверхности и отношения хода поршня к
диаметру цилиндра двигателя возрастает. Из всего количества теплоты, отдаваемой
охладителю, наибольшая часть (до 60-65%) воспринимается стенками цилиндра и
камеры сгорания во время процессов сгорания и расширения, остальная часть - в
течение выпуска.
Уменьшение доли теплоты, отводимой с охладителем,
повышает долю теплоты, отводимой с выпускными газами. В двигателях с
газотурбинным наддувом такое перераспределение теплоты приводит к повышению
располагаемой энергии выпускных газов и, следовательно, к увеличению работы в
газовой турбине. Этот принцип реализуется в комбинированных двигателях с
уменьшенным теплоотводом от рабочего тела («адиабатный дизель») за счет
тепловой изоляции деталей камеры сгорания и выпускного трубопровода
керамическими материалами.
3. РАСЧЕТ ВНЕШНИХ СКОРОСТНЫХ
ХАРАКТЕРИСТИК
На основании тепловых расчетов, проведенных для
четырех скоростных режимов работы бензиновых двигателей, получены и сведены в
таблицу 12 необходимые величины параметров.
Таблица 12
Параметры внешней скоростной характеристики
Частота
вращения коленчатого вала, мин-1
|
Параметры
внешней скоростной характеристики
|
|
Ne,
кВт
|
ge г/(кВтч)
|
Ме, Н м
|
GT, кг/ч
|
ηυ
|
α
|
10008,6834682,933,000,87440,86
|
|
|
|
|
|
|
320027,7732982,919,140,91670,96
|
|
|
|
|
|
|
560048,5931182,915,110,87840,96
|
|
|
|
|
|
|
600052,0630782,915,980,87010,96
|
|
|
|
|
|
|
Коэффициент приспособляемости по скоростным
характеристикам карбюраторного двигателя: К=Меmax/МeN=82,93/82,91=1;
Для сравнения различных методов построения
скоростных характеристик и проверки правильности выполнения теплового расчета
для нескольких скоростных режимов двигателя дополнительно приведен расчет
изменения мощности и удельного расхода топлива на основе процентных соотношений
между параметрами относительной скоростной характеристики карбюраторного
двигателя. Результаты расчета сведены в таблице 13. Внешние скоростные
характеристики (в соответствии с рисунком. 2, 3, 4).
Таблица 13
Параметры внешней скоростной характеристики
Частота
вращения коленчатого вала, пх
|
Мощность,
Ne
|
Удельный
расход топлива,
ge
|
%
|
мин-1
|
%
|
кВт
|
%
|
г/(кВт-ч)
|
20
|
1120209,72115358
|
|
|
|
|
40
|
22405019,44100311
|
|
|
|
|
60
|
33607329,1597302
|
|
|
|
|
80
|
44809238,8795295
|
|
|
|
|
100
|
560010048,59100311
|
|
|
|
|
120
|
67209258,31115358
|
|
|
|
|
Рисунок 2. Внешние скоростные характеристики ge
и Ме.
Рисунок 3. Внешние скоростные характеристики Ne
и GT.
Рисунок 4. Внешние скоростные характеристики
ηυ и α.
На основе сравнения полученных данных с кривыми Ne
и gе,
построенными по результатам теплового расчета, можно сделать следующие выводы:
Точки относительной характеристики практически
полностью совпадают с внешней скоростной характеристикой мощности
рассчитываемого двигателя.
Точки относительной характеристики
удельного расхода топлива несколько отличаются от кривой gе,
построенной по данным теплового расчета, в сторону увеличения gе и особенно
при малых значениях частоты вращения коленчатого вала. Максимальное расхождение
составляет при п=1000 мин-1 около 23% [350 и 284 г/(кВт·ч)].
4. КИНЕМАТИКА
Выбор λ и длины Lш
шатуна.
В целях уменьшения высоты двигателя
без значительного увеличения инерционных и нормальных сил отношение радиуса
кривошипа к длине шатуна предварительно было принято в тепловом расчете λ=0,285 [6]
При этих условиях Lш=R/λ=32,5/0,285=114,0
мм.
Принятые значения Lш
и λ
обеспечивают движение шатуна без задевания за нижнюю кромку цилиндра.
Следовательно, перерасчета величин Lш и λ не
требуется.
Перемещение поршня
x=R·[(1-соs
φ)+ (1-соs 2φ)]=32,5·[(1-соs
φ)+ (1-соs 2φ)] мм.
(4.1.1)
Расчет sx
производится аналитически через каждые 10° угла поворота коленчатого вала.
Угловая скорость вращения
коленчатого вала:
ω=πn/30=3,14·5600/30=586
рад/с.
υп=ωR(sin φ+sin 2φ)=586·0,033×(sin φ+0,285/2sin 2φ) м/с.
(4.1.2)
Ускорение поршня:
j=ω2R(cos φ+λcos 2φ)=5862·0,033
(соs φ+0,285 соs 2φ) м/с2
(4.1.3)
Параметры перемещения поршня даны в
таблице 14.
Построены графики (в соответствии с
рисунком 5, 6, 7).
Таблица 14
Параметры перемещения поршня
φ°
|
[(1-cosφ)+(1-cos2φ]sx
мм(sinφ+sin2φ)υп м/с(соsφ+0,285 соs2φ)j, м/с2
|
|
|
|
|
|
1
|
2
|
3
|
4
|
5
|
6
|
7
|
0
|
0,0000
|
0,0
|
0,0000
|
0,0
|
+1,2860
|
+14352
|
30
|
+0,1697
|
5,5+0,6234+11,9+1,0085+11255
|
|
|
|
|
60
|
+0,6069
|
19,7+0,9894+18,8+0,3575+3990
|
|
|
|
|
90
|
+1,1425
|
37,1+1,0000+19,0-0,2850-3181
|
|
|
|
|
120
|
+1,6069
|
52,2+0,7426+14,1-0,6425-7171
|
|
|
|
|
150
|
+1,9017
|
61,8+0,3766+7,2-0,7235-8075
|
|
|
|
|
180
|
+2,0000
|
65,00,00000,0-0,7150-7980
|
|
|
|
|
210
|
+1,9017
|
61,8-0,3766-7,2-0,7235-8075
|
|
|
|
|
240
|
+1,6069
|
52,2-0,7426-14,1-0,6425-7171
|
|
|
|
|
270
|
+1,1425
|
37,1-1,0000-19,0-0,2850-3181
|
|
|
|
|
300
|
+0,6069
|
19,7-0,9894-18,8+0,3575+3990
|
|
|
|
|
330
|
+0,1697
|
5,5-0,6234-11,9+1,0085+11255
|
|
|
|
|
360
|
+0,0000
|
0,0
|
-0,0000
|
0,0
|
+1,2850
|
+14341
|
По табличным данным строем графики sx
в масштабе Мz=2
мм в мм, υп-в
масштабе Му=1 м/с в мм, j-в
масштабе Mj=500
м/с2 в мм. Масштаб угла поворота коленчатого вала Мφ=3°
в мм.
При j=0
υп=±υmax,
а на кривой sx-это точка перегиба.
Рисунок 5. График перемещения поршня.
Рисунок 6. График скорости поршня.
Рисунок 7. График ускорения поршня.
5. ДИНАМИКА
5.1 Силы давления газов
Индикаторную диаграмму, полученную в тепловом
расчете, развертывают по углу поворота кривошипа по методу Брикса (в
соответствии с рисунком 8).
Рисунок 8. Развернутая индикаторная диаграмма.
Основные параметры двигателя даны в таблице 15.
Таблица 15
Основные параметры двигателя
Частот
вращения n, об/мин
|
5600D65n21,251
|
|
|
|
|
Радиус
кривошипа R, мм
|
32,5ε8,5Pb0,4676
|
|
|
|
|
Длина
шатуна L, мм
|
114,0Pa0,085Pr0,1180
|
|
|
|
|
Постоянная
механизма λ
|
0,285n11,377Fп33170
|
|
|
|
|
Угловая
скорость коленчатого вала
|
586Pz6,8016Vh0,20839
|
|
|
|
|
Поправка Брикса:
Rλ/(2Ms)=32,5·0,285/(2·1)=4,63
мм,
где Мs-масштаб
хода поршня на индикаторной диаграмме.
Масштабы развернутой диаграммы:
давлений и удельных сил Мр=0,05
МПа в мм;
полных сил МF=МрFп=0,05·0,003317=0,000166
МН в мм, или Мр=1659 Н в мм, угла поворота кривошипа Мφ=3° в мм,
или М'φ=4π/ОВ=4·3,14/240=0,0523
рад в мм, где ОВ - длина развернутой индикаторной диаграммы, мм.
По развернутой диаграмме через
каждые 10° угла поворота кривошипа определяют значения Δрr и заносят в
таблицу 16.
Таблица 16
Значения давлений и удельных сил
φо Δрr.
МПа j,
м/с2 pj,
МПа p, МПа tg
β pN,
МПа ps,
МПаpк,
МПарТ,
МПаТ,
кНМкр.ц,
Н·м
|
|
|
|
|
|
|
|
|
1
|
2
|
3
|
4
|
5
|
6
|
7
|
8
|
9
|
10
|
11
|
12
|
13
|
14
|
15
|
0
|
0,18
|
+14352-2,029-2,2090,0000,0001,000-2,209+1-2,2090000
|
|
|
|
|
|
|
|
|
|
|
|
|
30
|
0,18
|
+11255-1,591-1,771+0,144-0,3181,010-1,789+0,794-1,406+0,625-1,107-3,672119,3
|
|
|
|
|
|
|
|
|
|
|
|
|
60
|
0,18
|
+3990-0,564-0,630+0,253-0,1591,031-0,650+0,281-0,177+0,993-0,626-2,07667,5
|
|
|
|
|
|
|
|
|
|
|
|
|
90
|
0,18
|
-3181+0,450+0,630+0,295-0,1861,043+0,657-0,295-0,186+1-0,630+2,09067,9
|
|
|
|
|
|
|
|
|
|
|
|
|
120
|
0,18
|
-7171+1,014+1,194+0,253+0,3021,031+1,231-0,719-0,858+0,740+0,884+2,93295,3
|
|
|
|
|
|
|
|
|
|
|
|
|
150
|
0,18
|
-8075+1,142+1,322+0,144+0,191,010+1,335-0,938-1,240+0,376+0,497+1,64953,6
|
|
|
|
|
|
|
|
|
|
|
|
|
180
|
0,18
|
-7980+1,128+1,3080,0000,0001,000+1,308-1-1,7110000
|
|
|
|
|
|
|
|
|
|
|
|
|
210
|
0,03
|
-8075+1,142+1,172-0,144-0,1691,010+1,184-0,938-1,099-0,376-0,441-1,46347,5
|
|
|
|
|
|
|
|
|
|
|
|
|
240
|
0,03
|
-7171+1,014+1,044-0,253-0,2641,031+1,076-0,719-0,751-0,740-0,773-2,56483,3
|
|
|
|
|
|
|
|
|
|
|
|
|
270
|
0,05
|
-3181+0,450+0,500-0,295-0,1481,043+0,522-0,295-0,148-1-0,500-1,65953,9
|
|
|
|
|
|
|
|
|
|
|
|
|
300
|
0,12
|
+3990-0,564-0,684-0,253+0,1731,031-0,705+0,281-0,192-0,993+0,679+2,25273,2
|
|
|
|
|
|
|
|
|
|
|
|
|
330
|
0,49
|
+11255-1,591-2,081-0,144+0,31,010-2,102+0,794-1,652-0,625+1,301+4,315140,2
|
|
|
|
|
|
|
|
|
|
|
|
|
360
|
4,15
|
+14341-2,028-6,1780,0000,0001,000-6,178+1-6,1780000
|
|
|
|
|
|
|
|
|
|
|
|
|
370
|
6,9
|
+14352-2,029+8,929+0,050+0,4461,001+8,938+0,976+8,715+0,222+1,982+6,574213,7
|
|
|
|
|
|
|
|
|
|
|
|
|
390
|
2,28
|
+11255-1,591+3,871+0,144+0,5571,010+3,910+0,794+3,074+0,625+2,419+8,024260,8
|
|
|
|
|
|
|
|
|
|
|
|
|
420
|
0,63
|
+3990-0,564+1,194+0,253+0,3021,031+1,231+0,281+0,336+0,993+1,186+3,934127,9
|
|
|
|
|
|
|
|
|
|
|
|
|
450
|
0,3
|
-3181+0,450+0,750+0,295+0,2211,043+0,782-0,295+0,221+1+0,750+2,48880,9
|
|
|
|
|
|
|
|
|
|
|
|
|
480
|
0,2
|
-7171+1,014+1,214+0,253+0,3071,031+1,252-0,719-0,873+0,740+0,898+2,97996,8
|
|
|
|
|
|
|
|
|
|
|
|
|
510
|
0,16
|
-8075+1,142+1,302+0,144-0,1871,010+1,315-0,938-1,221+0,376+0,490+1,62552,8
|
|
|
|
|
|
|
|
|
|
|
|
|
540
|
0,27
|
-7980+1,128+1,3980,0000,0001,000+1,398-1-1,3980000
|
|
|
|
|
|
|
|
|
|
|
|
|
570
|
0,17
|
-8075+1,142+1,312-0,144-0,1891,010+1,325-0,938-1,231-0,376-0,493-1,63553,1
|
|
|
|
|
|
|
|
|
|
|
|
|
600
|
0,17
|
-7171+1,014+1,184-0,253-0,3001,031+1,221-0,719-0,851-0,740-0,876-2,90694,4
|
|
|
|
|
|
|
|
|
|
|
|
|
630
|
0,17
|
-3181+0,450+0,620-0,295-0,1831,043+0,647-0,295-0,183-1-0,620-2,05766,9
|
|
|
|
|
|
|
|
|
|
|
|
|
660
|
0,17
|
+3990-0,564-0,734-0,253+0,1861,031-0,757+0,281-0,206-0,993+0,729+2,41878,6
|
|
|
|
|
|
|
|
|
|
|
|
|
690
|
0,17
|
+11255-1,591-1,761-0,144+0,2541,010-1,779+0,794-1,398-0,625+1,101+3,652118,7
|
|
|
|
|
|
|
|
|
|
|
|
|
720
|
0,17
|
+14341-2,028-2,1980,0000,0001,000-2,198+1-2,1980000
|
|
|
|
|
|
|
|
|
|
|
|
|
5.2 Приведение масс частей
кривошипно-шатунного механизма
С учетом диаметра цилиндра, отношения S/D
устанавливаются:
-масса поршневой группы (для поршня
из алюминиевого сплава принято т'п=100 кг/м2): тп=т'п·Fп=100·0,003317=0,332
кг;
масса шатуна (для стального кованого
шатуна принято т'ш=150 кг/м2):
тш=т'ш·Fп=150 0,003317=0,498
кг;
масса неуравновешенных частей одного
колена вала без противовесов (для литого чугунного вала принято m'к=140
кг/м2):
тк=m'к
Fп=140 0,003317=0,464 кг.
Масса шатуна, сосредоточенная на оси
поршневого пальца:
тш.п.=0,275тш=0,275·0,498=0,137
кг.
Масса шатуна, сосредоточенная на оси
кривошипа:
mш.к.=0,725mш=0,725·0,498=0,361
кг.
Массы, совершающие
возвратно-поступательное движение:
mj=тп+тш.п=0,332+0,137=0,469
кг.
Массы, совершающие вращательное
движение:
mR=тк+mш.к.=0,464+0,361=0,825
кг.
5.3 Удельные и полные силы инерции
Удельная сила инерции
возвратно-поступательно движущихся масс:
pj=-jmj/Fn=-j·469·10-6/0,003317=-j141·10-6
МПа.
Центробежная сила инерции
вращающихся масс:
R=-mRRω2=-0,825 0,033-5862·10-3=-9,207
кН.
Центробежная сила инерции
вращающихся масс шатуна:
КRш=-mш.к·Rω2=-0,361·0,033
5862·10-3=-4,029 кН.
Центробежная сила инерции
вращающихся масс кривошипа:
КRк=-mк·Rω2=-0,464-0,033·5862·10-3=-5,178
кН.
5.4 Удельные суммарные силы
Удельная сила (МПа), сосредоточенная
на оси поршневого пальца:
=Δpr+pj
Удельная нормальная сила (МПа): pN=p·tg
β
Удельная сила (МПа), действующая
вдоль шатуна: ps=p (1/cos β)
Удельная сила (МПа), действующая по
радиусу кривошипа:
рх=р
cos (φ+β)/cos β
Удельная и полная тангенциальные
силы (МПа и кН):
T=p·sin(φ+β)/cos β и T=pT ·Fп=pT·0,003317·103
По табличным данным строят графики
изменения удельных сил pj, р, рS, pN, рК
и рТ в зависимости от изменения угла поворота коленчатого
вала φ
(в
соответствии с рисунком 9, 10, 11).
Рисунок 9. Графики изменения pj
и p.
Рисунок10. Графики изменения рN и ps.
Рисунок 11. Графики изменения рК
и pТ.
Среднее значение тангенциальной силы
за цикл:
по данным теплового расчета:
Тср=рiFп=1,4093·0,003317=744
Н; (5.4.1)
по площади, заключенной между кривой
рт и осью абсцисс:
рТср=Мр=0,05=0,215 МПа, (5.4.2)
а Тср=рТср·Fп=0,215 0,003317·106=713
Н
ошибка Δ=(744-713)100/744=4,2%.
5.5 Крутящие моменты
Крутящий момент одного цилиндра: Мкр.ц=ТR=T·0,033·103
Н·м
Период изменения крутящего момента
четырехтактного двигателя с равными интервалами между вспышками: θ=720/i=720/4=180°
Суммирование значений крутящих
моментов всех четырех цилиндров двигателя осуществляется табличным методом
(Таблица 17) через каждые 10° угла поворота коленчатого вала и по полученным
данным строится кривая Мкр в масштабе Мм=10
Н·м в мм. (в соответствии с рисунком 12).
Рисунок 12. График изменения
крутящего момента.
Таблица 17
Значения крутящих моментов
φ°
|
Цилиндры
|
Мкр, Н·м
|
|
1-й
|
2-й
|
3-й
|
4-й
|
|
|
φ° кривошипа
|
Мкр, Н·м
|
φ° кривошипа
|
Мкр, Н·м
|
φ° кривошипа
|
Мкр, Н·м
|
φ° кривошипа
|
Мкр, Н·м
|
|
0
|
0
|
0
|
180
|
0
|
360
|
0
|
540
|
0
|
0
|
10
|
10
|
-132,4
|
190
|
-28,0
|
370
|
+125,5
|
550
|
-29,8
|
-64,7
|
20
|
20
|
-203,2
|
200
|
-65,2
|
380
|
+161,3
|
560
|
-67,1
|
-174,2
|
30
|
30
|
-223,3
|
210
|
-94,6
|
390
|
+176,6
|
570
|
-97,4
|
-238,7
|
40
|
40
|
-209,7
|
220
|
-123,0
|
400
|
+156,6
|
580
|
-126,8
|
-302,9
|
50
|
50
|
-174,3
|
230
|
-154,7
|
410
|
+128,6
|
590
|
-155,7
|
-356,1
|
60
|
60
|
-127,6
|
240
|
-165,2
|
420
|
+124,8
|
600
|
-337,7
|
70
|
70
|
-65,2
|
250
|
-159,4
|
430
|
+142,6
|
610
|
-165,9
|
-247,9
|
80
|
80
|
+19,6
|
260
|
-139,8
|
440
|
+185,5
|
620
|
-146,3
|
-81,0
|
90
|
90
|
+97,4
|
270
|
-103,9
|
450
|
+234,3
|
630
|
-103,5
|
+124,3
|
100
|
100
|
+140,7
|
280
|
-36,3
|
460
|
+248,9
|
640
|
-39,1
|
+314,2
|
110
|
110
|
+160,3
|
290
|
+32,6
|
470
|
+244,2
|
650
|
+61,5
|
+498,6
|
120
|
120
|
+165,2
|
300
|
+97,0
|
480
|
+229,3
|
660
|
+121,4
|
+612,9
|
130
|
130
|
+152,9
|
310
|
+123,0
|
490
|
+190,1
|
670
|
+183,6
|
+649,6
|
140
|
140
|
+127,7
|
320
|
+138,0
|
500
|
+147,3
|
680
|
+208,8
|
+621,8
|
150
|
150
|
+94,6
|
330
|
+137,8
|
510
|
+115,3
|
690
|
+219,6
|
+567,3
|
160
|
160
|
+63,4
|
340
|
+120,2
|
520
|
+76,4
|
700
|
+201,3
|
+461,3
|
170
|
170
|
+29,8
|
350
|
+71,8
|
530
|
+30,8
|
710
|
+139,8
|
+272,2
|
180
|
180
|
0
|
360
|
0
|
540
|
0
|
720
|
0
|
0
|
Средний крутящий момент двигателя:
по данным теплового расчета
Mкp.cp=Mi=Mе/ηм=82,9/0,8591=96,5
Н·м;
по площади, заключенной под кривой Мкр:
Мкр.ср=Мм=10=92,3 Н м; (5.5.1)
ошибка Δ=100=4,4 %.
Максимальный и минимальный крутящие
моменты:
Мкр.max=5 Н·м; Мкр.min=-1 Н·м. (5.5.2)
5.6 Силы, действующие на шатунную
шейку коленчатого вала
Для проведения расчета
результирующей силы, действующей на шатунную шейку рядного двигателя,
составляют таблицу 18.
Таблица 18
Полные силы
φº
|
Полные
силы, кН
|
|
T
|
K
|
Pк
|
Rш.ш.
|
КРк
|
Rк
|
0
|
0
|
-7,327-11,35611,356-2,9698,387
|
|
|
|
|
30
|
-3,672-4,664-8,69312,365-3,7728,593
|
|
|
|
|
|
60
|
-2,076-0,587-4,6166,692-5,0011,691
|
|
|
|
|
|
90
|
+2,090-0,617-4,6466,736-4,9921,744
|
|
|
|
|
|
120
|
+2,932-2,846-6,8759,807-4,325,487
|
|
|
|
|
|
150
|
+1,649-4,113-8,1429,791-3,9385,853
|
|
|
|
|
|
180
|
0
|
-5,675-9,7049,704-3,4676,237
|
|
|
|
|
210
|
-1,463-3,645-7,6749,137-4,0795,058
|
|
|
|
|
|
240
|
-2,564-2,491-6,529,084-4,4274,657
|
|
|
|
|
|
270
|
-1,659-0,491-4,526,179-5,031,149
|
|
|
|
|
|
300
|
+2,252-0,637-4,6666,918-4,9861,932
|
|
|
|
|
|
330
|
+4,315-5,480-9,50912,365-3,5268,839
|
|
|
|
|
|
360
|
0
|
-20,492-24,52124,521-125,521
|
|
|
|
|
370
|
+6,574+28,908+24,87931,453-3,53734,990
|
|
|
|
|
|
390
|
+8,024+10,196-6,16714,191-2,10412,087
|
|
|
|
|
|
420
|
+3,934+1,115-2,9141,020-4,842-3,822
|
|
|
|
|
|
450
|
+2,488-0,733-4,7627,250-4,9572,293
|
|
|
|
|
|
480
|
+2,979-2,896-6,9259,904-4,3055,599
|
|
|
|
|
|
510
|
+1,625-4,050-8,0799,704-3,9575,747
|
|
|
|
|
|
540
|
0
|
-4,637-8,6668,666-3,784,886
|
|
|
|
|
570
|
-1,635-4,083-8,1129,747-3,9475,800
|
|
|
|
|
|
600
|
-2,906-2,823-6,8529,758-4,3275,431
|
|
|
|
|
|
630
|
-2,057-0,607-4,6366,693-4,9951,698
|
|
|
|
|
|
660
|
+2,418-0,683-4,7127,130-4,9722,158
|
|
|
|
|
|
690
|
+3,652-4,637-8,66612,318-3,788,538
|
|
|
|
|
|
720
|
0
|
-7,291-11,3211,320-2,988,340
|
|
|
|
|
Суммарная сила, действующая на
шатунную шейку по радиусу кривошипа: Рк=К+КRш=(К-4,029)
кН, где К=рк Fп=рк
0,003317·103 кН.
Результирующая сила Rш.ш,
действующая на шатунную шейку, подсчитывается графическим сложением векторов
сил Т и Рх при построении полярной диаграммы (в
соответствии с рисунком 13).
Рисунок 13. Полярная диаграмма.
Масштаб сил на полярной диаграмме
для суммарных сил МР=0,1 кН в мм. Значения Rш.ш для
различных φ> заносят
в таблицу и по ним строят диаграмму Rш.ш в
прямоугольных координатах (в соответствии с рисунком 14).
Рисунок 14. Диаграмма нагрузки на
шатунную шейку.
По развернутой диаграмме Rш.ш определяют:
Rш.ш.cp=FMP/OB=24030·0,05/74=16,303
кН;
Rш.ш.max=31,453 кН; Rш.ш.min=1,02 кН,
где ОВ - длина диаграммы, мм;
F-площадь под кривой Rш.ш, мм2.
По полярной диаграмме строят
диаграмму износа шатунной шейки. Сумму сил Rш.ш.j,
действующих по каждому лучу диаграммы износа, определяют с помощью таблицы 19.
Таблица 19
Значение сил, действующих по каждому
лучу
α
|
Значения
Q, кН для точек
|
|
1
|
2
|
3
|
4
|
5
|
6
|
7
|
8
|
9
|
10
|
11
|
12
|
10
|
44
|
44
|
|
|
|
|
|
|
|
|
|
44
|
40
|
39
|
39
|
|
|
|
|
|
|
|
|
|
39
|
70
|
36
|
36
|
|
|
|
|
|
|
|
|
|
36
|
100
|
37
|
37
|
37
|
|
|
|
|
|
|
|
|
37
|
130
|
32
|
32
|
32
|
|
|
|
|
|
|
|
|
32
|
160
|
22
|
22
|
22
|
|
|
|
|
|
|
|
|
22
|
190
|
19
|
19
|
|
|
|
|
|
|
|
|
19
|
19
|
220
|
29
|
29
|
|
|
|
|
|
|
|
|
29
|
29
|
250
|
39
|
39
|
|
|
|
|
|
|
|
|
39
|
39
|
280
|
40
|
40
|
|
|
|
|
|
|
|
|
40
|
40
|
310
|
36
|
36
|
|
|
|
|
|
|
|
|
36
|
36
|
340
|
102
|
|
|
|
|
|
|
|
102
|
102
|
102
|
102
|
370
|
|
|
|
|
427
|
427
|
427
|
427
|
|
|
|
|
400
|
|
|
144
|
144
|
144
|
144
|
|
|
|
|
|
|
430
|
73
|
73
|
73
|
73
|
|
|
|
|
|
|
|
|
460
|
62
|
62
|
62
|
62
|
|
|
|
|
|
|
|
|
490
|
52
|
52
|
52
|
|
|
|
|
|
|
|
|
52
|
520
|
37
|
37
|
37
|
|
|
|
|
|
|
|
|
37
|
550
|
22
|
22
|
|
|
|
|
|
|
|
|
22
|
22
|
580
|
28
|
28
|
|
|
|
|
|
|
|
|
28
|
28
|
610
|
36
|
36
|
|
|
|
|
|
|
|
|
36
|
36
|
640
|
36
|
36
|
|
|
|
|
|
|
|
|
36
|
36
|
670
|
38
|
38
|
38
|
|
|
|
|
|
|
|
|
38
|
700
|
43
|
43
|
43
|
|
|
|
|
|
|
|
|
43
|
ΣQi
|
902
|
800
|
278
|
571
|
571
|
427
|
427
|
102
|
102
|
388
|
768
|
qi
|
1
|
0,54
|
0,36
|
0,19
|
0,38
|
0,38
|
0,29
|
0,29
|
0,07
|
0,07
|
0,26
|
0,52
|
По данным таблицы в масштабе МР=50
кН в мм по каждому лучу откладывают величины суммарных сил ΣRш.ш.j,
от окружности к центру (в соответствии с рисунком 15). По диаграмме износа
определяют расположение оси масляного отверстия (φм=68°).
Рисунок 15. Диаграмма износа шатунной шейки
карбюраторного двигателя.
5.7 Силы, действующие на колено вала
Суммарная сила, действующая на колено вала по
радиусу кривошипа:
Pк=Pк+KRк=Рк-8,960
кН (5.7.1)
Результирующая сила, действующая на колено вала Rк=Rш.ш+KРк,
определяется по диаграмме Rш.ш. Векторы из полюса Ок
до соответствующих точек на полярной диаграмме в масштабе МР=0,1
кН в мм выражают силы Rк, значения которых для различных φ
заносят в таблицу.
5.8 Равномерность крутящего момента
и равномерность хода двигателя
Равномерность крутящего момента:
μ=(Mкр.max-Mкр.min)/Mкp.cp=[5-(-1)]/92,3=0,04 (5.8.1)
Избыточная работа крутящего момента:
Lизб=Fabc·ММ·М′φ=841·10·0,3622=3046,1
Дж, (5.8.2)
где Fabc-площадь
над прямой среднего крутящего момента мм2; М′φ=4π/(iОА)=4·3,14/(4·9)=0,3622
рад в мм-масштаб угла поворота вала на диаграмме Мкр.
Равномерность хода двигателя
принимаем δ=0,01.
Момент инерции движущихся масс
двигателя, приведенных к оси коленчатого вала: J0=Lизб/(δω2)=3046,1/(0,01
5862)=0,887 кг·м2
ЗАКЛЮЧЕНИЕ
Курсовая работа “Тепловой и динамический расчет
двигателя” выполнена в соответствии с заданием на основе методической и учебной
технической литературы.
Рассчитанные показатели рабочего цикла, работы,
размеров, кинематики и динамики проектируемого двигателя отличаются от
прототипа топливной экономичностью и габаритными размерами. Снижение удельного
расхода топлива на 17 г/кВт ч достигнуто уменьшением хода и диаметра поршня,
т.е. снижением габаритов, скорости и потерь на трение.
В целом из выполненной работы следуют выводы:
1. Обоснованы исходные данные для проектирования
эффективного двигателя по заданию с учетом прототипа и методических
рекомендаций.
2. Рассчитаны с применением ЭВМ рабочий
цикл, работа и размеры двигателя, его удельные мощности и топливные показатели,
кинематика и динамика, регуляторная (нагрузочная) характеристика. Проектируемый
двигатель отличается повышенной топливной экономичностью и меньшими габаритами.
. Получены навыки расчета и опыт
оформления материалов по проектированию автомобильного двигателя, отвечающего
современным техническим требованиям.
СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ
1. Краткий
автомобильный справочник НИИАТ, М.: Транспорт, 1984.
2. Колчин
А. И., Демидов В. П. «Расчет автомобильных и тракторных двигателей», М.: Высшая
школа, 1980.
. Ховах
М. С., Вихерт М. М., Воинов А. Н. и др. «Автомобильные двигатели», М.:
Машиностроение, 1977.
. Ленин
И. М., Попык К. Г, Малашкин О. М. и др. «Автомобильные и тракторные двигатели»,
М.: Машиностроение, 1969.
. Орлин
А. С., Вырубов Д. Н., Круглов М. Г. и др. «Двигатели внутреннего сгорания.
Конструкция и расчет.», М.: Машиностроение, 1972.
6.Вахламов
В. К. Автомобили, основы конструкции М: Асадема, 2004 г.