
- •1. Система общих уравнений моделирования физических явлений в проточных частях двс. Газодинамические и термодинамические методы моделирования.
- •2. Термодинамическая математическая модель расчета установившихся и переходных процессов дизельного двигателя
- •2.1. Моделирование процессов тепловыделения и теплообмена в цилиндрах
- •2.3. Определение состава и объема рабочего тела в расчетной зоне
- •2.3. Моделирование процессов во впускных и выпускных
- •2.4. Моделирование совместной работы дизеля и турбокомпрессора
2.3. Определение состава и объема рабочего тела в расчетной зоне
Текущий объем цилиндра и его увеличение вычисляется по формулам (1.43) [56] или:
; (46)
. (47)
Количество молей смеси в цилиндре можно определить как:
(48)
а массовое количество смеси как:
(49)
где
– коэффициент молекулярного изменения
при aц
= 1;
g – коэффициент остаточных газов;
aц – коэффициент избытка воздуха в цилиндре;
x – доля топлива, которое сгорело к данному моменту времени.
Молекулярная масса смеси равняется
(50)
где
- молекулярная доля в смеси продуктов
сгорания,
-
молекулярная доля воздуха в смеси.
Удельные мольные теплоемкости воздуха, продуктов сгорания и их смеси определяются с учетом их зависимости от температуры по формулам:
(51)
(52)
(53)
а массовая теплоемкость по формуле:
(54)
Численные значения коэффициентов “а” и “b” в формулах (51) и (52) приведенные в литературе, например, в [1, 3]. Наиболее точными эмпирическими выражениями теплоемкостей есть формулы (1.34)…(1...37) [56]. Аналогично вычисляются mСр пов и mСр г.с.
2.3. Моделирование процессов во впускных и выпускных
системах
Интенсивность потока массы через контрольные поверхности газораспределительных органов в зависимости от вида истечения определяется по формулам:
- для подкритического истечения:
(55)
- для надкритического истечения:
(56)
где р1 и r1 - давление и плотность рабочего тела со стороны контрольной поверхности, где они больше;
р2 - давление рабочего тела из той стороны контрольной поверхности, где оно меньше;
µiFi - эффективные проходные сечения газораспределительных органов и элементов проточных частей ДВЗ (расчетные, как в [1], или полученные экспериментально путем продувки);
k - показатель адиабаты (берется с учетом состава и температуры рабочего тела по формуле (1.33)).
Энтальпия рабочего тела, которое пересекает контрольную поверхность, определяется по заторможенным параметрам:
(57)
где Cpmi и Ti - средняя изобарная массовая теплоемкость и температура с той стороны контрольной поверхности, где давление больше.
В результате численного интегрирования системы уравнений (21)...(23) можно получить текущие значения давления, температуры и состава рабочего тела в цилиндре, а также показатели качества процесса газообмена: масса воздушного заряда, коэффициенты наполнения, утечки продувочного воздуха и остаточных газов, среднее давление насосных потерь (см. (1.1)...(1.4) [56]).
Параметры рабочего тела в выпускном коллекторе определяются путем численного интегрирования уравнений,
(58)
(59)
(60)
которые получены из основной системы (18)…(20).
Здесь qv=0 и dV=0, так как в выпускной системе отсутствующее тепловыделение и объем коллектора постоянный V=const, значение количества цилиндров в данном выпускном коллекторе z1 зависит от типа системы выпуска: одноколлекторная ли многоколлекторная. Для одноколлекторной системы, как источника массы выступают все цилиндры (z1=z), а как сток – входной патрубок турбины или сопло преобразователя импульсов, поэтому суммирование по количеству контрольных поверхностей для каждого выпускного коллектора идет к (z1 + 1). Величиной потерь теплоты q через поверхности, которые ограждают объем коллектора Fq, в теплоизолированных, как правило, выпускных коллекторах дизелей можно пренебречь, поэтому:
(61)
Если между выпускными коллекторами и турбиной имеются преобразователи импульсов, как это показан на схеме рис.1, их расчет можно вести по методу, изложенному в работе [41].