Термодинамический расчет цикла ДВС

1. Расчет цикла двигателя внутреннего сгорания

Краткое описание процессов, составляющих цикл карбюраторного двигателя

Идеализированный цикл карбюраторного двигателя представлен циклом Карно. В этом цикле подвод и отвод теплоты реализуется в процессах V=const, а сжатие свежего заряда и расширение продуктов сгорания – в политропических процессах с отводом теплоты (с постоянными значениями показателей политроп).

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

Состав топлива

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

Основные реакции при горении топлива имеют вид:

Под реакциями подписаны молярные массы веществ, участвующих в реакциях, а в правых частях в общем виде записано количество теплоты, выделяющейся в этих реакциях. На основании этих записей можно составить формулу для расчета теоретически необходимого количества воздуха для сгорания 1 кг топлива. Следует учесть количество кислорода, содержащегося в топливе, и массовую долю кислорода в воздухе (0,23):

где M₀ – масса воздуха, необходимая для сгорания 1 кг топлива, кг; C, H, O – массовые доли углерода, водорода и кислорода в топливе.

Последнюю формулу можно записать в виде: (1) подставив значения получим кг

Действительное количество воздуха, подаваемое для сгорания 1 кг топлива

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

, (2)

где a – коэффициент избытка воздуха; в карбюраторных двигателях обычно a=0,8…1,15. Учитывая, что у нас a=1,14, получим кг.

Количество теплоты, выделяющееся при сгорании топлива

Если известны основные химические реакции, протекающие при сгорании топлива, и тепловые эффекты этих реакций, то легко записать формулу для вычисления суммарного количества теплоты, МДж/кг, выделяющейся при сгорании 1 кг топлива (формула Менделеева):

. (3)

При сгорании топлива часть теплоты уносится с водяными парами и не дает вклада в суммарное количество теплоты (низшая теплота сгорания топлива). Подставим значения: МДж/кг.

Расчет процесса сжатия

Параметры начальной точки

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

T₁=(350…430) K;

p₁=(0,9…0,95)*10⁵ Па (в тихоходных двигателях);

p₁=(0,75…0,85)*10⁵ Па (в быстроходных двигателях);

Сравнительно высокие значения температуры в начальной точке связаны с нагревом воздуха во входных каналах двигателя.

Расчет процесса сжатия свежего заряда

4.2.1. Молекулярная масса свежего заряда определяется по формуле

, (4)

здесь m_б, m_в – массовые доли паров бензина и воздуха; m_б, m_в-молярные массы паров бензина и воздуха.

Масса свежего заряда – M_с.з.= 1 кг паров бензина + 16,9 кг воздуха = 17,9 кг. Массовая доля паров бензина m_б==0,06, массовая доля воздуха m_в==0,94. Подставляем эти значения в (4): кг/кг*моль.

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

Среднее значение молярной теплоемкости для изохорического процесса в интервале температур 0-T рассчитывается по формуле (5), где .

Задаемся значением Т₂=625 К. ДЖ/кмоль*К, теперь можно определить величину удельной массовой теплоемкости (6) Дж/(кг*К).

Показатель адиабаты для процесса сжатия. Газовая постоянная для свежего заряда вычисляется по формуле (7) Дж/(кг*К)

Среднее значение теплоемкости при постоянном давлении (8) . Дж/(кг*К)

Показатель адиабаты для процесса сжатия (9) =1,378.

Показатель политропы для процесса сжатия. В задании приводится значение (n₁-k₁)=-D₁, поэтому n₁= k₁-D₁=1,378 – 0,009=1,37.

P₁*V₁₌RT; =>

Теперь можно определить параметры в конце процесса сжатия: м³ /кг, Па, К. Полученное значение температуры отличается от изначально принятого на 207К.

Зададимся другим значением Т₂.

Среднее значение молярной теплоемкости для изохорического процесса в интервале температур 0-T рассчитывается по формуле (5), где .

Задаемся значением Т₂=832 К. ДЖ/кмоль*К, теперь можно определить величину удельной массовой теплоемкости (6) Дж/(кг*К).

Показатель адиабаты для процесса сжатия. Газовая постоянная для свежего заряда вычисляется по формуле (7) Дж/(кг*К)

Среднее значение теплоемкости при постоянном давлении (8) . Дж/(кг*К)

Показатель адиабаты для процесса сжатия (9) =1,373.

Показатель политропы для процесса сжатия. В задании приводится значение (n₁-k₁)=-D₁, поэтому n₁= k₁-D₁=1,373 – 0,009=1,364.

P₁*V₁₌RT; =>

Теперь можно определить параметры в конце процесса сжатия: м³ /кг, Па, К. Полученное значение температуры отличается от изначально принятого на 8К.

Итерация: Взяли Т₂=832, получили 824 после второй подгонки.

Расчет процесса сгорания

Состав продуктов сгорания

Из основных реакций и следует, что в результате реакций на 1 кг С приходится 44/12=3,67 кг CO₂, а на 1 кг Н приходится 36/4=9 кг Н₂О.

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

Общая масса продуктов сгорания, кг:

М_п.с.=3,67С + 9Н + 0,77М₀ + (a – 1) М₀=3,14+1,305+11,51+1,94=17,89 кг

Массовые доли веществ, составляющих продукты реакции горения:

Молярная масса продуктов сгорания

Вычисляется по формуле (12):

кг/моль.

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

В интервале температур (Т₁, Т₂) для a1 определяется по формуле (14)

, где .

Задаемся значением Т₃ = 2850К кг/моль. Удельная массовая теплоемкость вычисляется по формуле

ДЖ/(кг*К).

Параметры в конце процесса сгорания

Температура в конце сгорания вычисляется по формуле (15) , где q_2,3 – количество теплоты выделившейся при сгорании 1 кг свежего заряда. Её можно вычислить по формуле (16) , где x_Z – коэффициент подвода теплоты, его значение – для карбюраторных двигателей находится в пределах 0,85–0,95, выбираем 0,9, x_a – учитывает меньшее выделение теплоты – x_a=1,4a-0,4, при α1

К. Полученная температура отличается от первоначально принятой на 5 К, что находится в пределах допустимого.

, V₃=V₂, .=> Па.

Итерация: Взяли Т₃=2850, получили 2845,3

Расчет процесса расширения продуктов сгорания

Показатель адиабаты

Задаемся значением температуры в конце процесса расширения Т₄=1610 К: К. Вычисление средних значений молярных теплоемкостей (в интервале температур) производится по формулам (13) и (14).

; ДЖ/кмоль

ДЖ/(кг*К),

, ДЖ/(кг*К)

.

Показатель политропы

,

Расчет процесса выхлопа газа

Па, V₄=V₁,

.

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

Итерация: Взяли Т₄=1610, получили 1620

Энергетические характеристики цикла

Уравнение теплового баланса

Для рассмотренного цикла можно записать баланс в виде: q_2,3 + q_1,2 + q_3,4 + q_4,1 =l_3,4 + l_1,2, или q_2,3 = q_и, где q_и – энергия, полученная в цикле q_и = l_3,4 + l_1,2 – q_1,2 – q_3,4 – q_4,1; (17)

l_1,2 – работа сжатия, (18)

Дж/кг,

l_3,4 – работа расширения, (19)

Дж/кг,

q_1,2 – теплота, отведенная в процессе сжатия, (20)

Дж,

q_3,4 – теплота, отведенная в процессе расширения, (21)

Дж,

q_4,1 – теплота, отведенная с выхлопными газами, (22) ;

Дж/(кмоль *К),

, ДЖ/(кг*К)

Дж.

Полезная работа

Дж. Проводим сопоставление количества теплоты, выделившейся в процессе сгорания 1 кг свежего заряда q_2,3 и суммы полезной работы и отведенной теплоты q_и в цикле. Эти величины не совпадают, поэтому подсчитываем относительную величину дисбаланса по формуле (25).

Дисбаланс не превышает пяти процентов, поэтому делаем вывод, что подсчет был правильным.

Среднее теоретическое индикаторное давление вычисляется по формуле (26)

Индикаторная диаграмма

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

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

Действительная индикаторная диаграмма меньше теоретической за счет отличия действительных процессов от теоретических. Уменьшение площади индикаторной диаграммы можно учесть с помощью коэффициента полноты диаграммы V=0,95, а механические потери – относительным механическим КПД h_м =0,95. Среднее эффективное давление цикла (27) МПа

Термический КПД цикла

(28) ®

Геометрические характеристики двигателя

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

(30) ® л

Определение диаметра цилиндра и рабочего хода поршня

При заданном значении .

=>

Расчет теплообменной поверхности радиатора

1. Исходные данные

Мощность двигателя Р_e =60Вт

Температура воды на входе t=90С

Температура воздуха на входе t=30

Скорость обдува =25 м/с

Высота радиатора Н=300 мм

Ширина В=50 мм

Размер трубки ba 245

Размещение трубок двухрядное

Шаг трубок s=15 м

Ребра стальные

Толщина 0,2 мм

Теплопроводность 53,6 Вт/(м К)

2. Расчет радиатора

Определение количества элементов n:

n = = =30

принимаем 42 шт.

Уточняем тепловой поток, отводимы одним элементом Q:

Вт

Расчет коэффициента теплоотдачи от воды к стенке трубки. Теплофизические свойства воды принимаем при температуре входа 95С:

кг/м; м/c; ; Р =1,95.

Определяем эквивалентный диаметр трубки:

а) площадь сечения трубки f

б) Смачиваемый периметр

в) Эквивалентный диаметр

мм

Вычисляем критерий Рейнольдса для течения воды в трубке, задавшись скоростью м/c:

Вычисляем критерий Нуссельта:

Определяем коэффициент теплоотдачи:

Расчет коэффициента теплоотдачи от стенки трубки к воздуху. Теплофизические свойства воздуха принимаем при температуре 28С: кг/м; м/c; Вт/(м*К); .

Вычисляем критерий Рейнольдса для течения воздуха в межтрубном пространстве, за характерный размер принимаем ширину радиатора В:

Вычисляем критерий Нуссельта:

Определяем коэффициент теплоотдачи:

Определение средней температуры теплоносителей:

Определяем массовый расход воды :

кг/с

Определяем массовый расход воздуха:

кг/с

Определяем среднюю температуру теплоносителей если теплоемкость воды и воздуха соответственно Дж/кг*К; Дж/кг*К:

Определение коэффициента эффективности оребрения.

Вычисляем длину ребра:

мм

Определяем безразмерный параметр х:

Находим коэффициент эффективности оребрения :

=th x / x=th 0,828 / 0,828=0,82

Предварительное определение площади оребрения.

Площадь боковой поверхности трубки :

м

Определим среднюю температуру стенки трубки:

Площадь поверхности оребрения:

Количество ребер :

Расстояние между ребрами:

Уточненный расчет.

Определяем критерий Рейнольдса, за эквивалентный диаметр принимаем 2h:

Вычисляем критерий Нуссельта:

Уточняем коэффициент теплоотдачи α_в от оребренной стенки к воздуху:

Уточняем температуру, для чего определяем живое сечение радиатора S и пересчитываем расход воздуха G_в:

Уточняем коэффициент эффективности оребрения:

Определяем свободную поверхность трубки между ребрами:

Уточняем площадь ребер:

Оцениваем погрешность:

и увеличиваем высоту трубки пропорционально недостающим процентам:

Определяем длину радиатора L, полагая двухрядное расположения трубок:

Определяем окончательные габариты радиатора, мм: