Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
2155.pdf
Скачиваний:
62
Добавлен:
07.01.2021
Размер:
3.48 Mб
Скачать

3.4. Измерение коэффициента наполнения цилиндра

По массе всасываемого в цилиндр ДВС воздуха дозируется топливо. Для косвенного измерения массового расхода воздуха (пошедшего через впускной клапан цилиндра ДВС) используют два способа: регистрируют скорость прохождения воздуха через входной коллектор, используя анемометрический датчик; измеряют давление в нем, используя датчик давления. Причем датчик давления устанавливают в области разряжения между дроссельной заслонкой и впускным клапаном.

На рис. 31, 32 показаны дроссели с различным диаметром отверстия, которые используются в учебном процессе для экспериментов.

РисСибАДИ. 31. Дроссели с различным диаметром отверстий

а

б

Рис. 32. Дроссели, установленные на макет: а дроссель с постоянным диаметром; б дроссель с изменяемым сечением отверстия

35

На рис. 33 показана кривая изменения давления во впускном коллекторе после дроссельной заслонки от времени за один цикл движения поршня в цилиндре макета.

СибАДИ

 

Рис. 33. Разбиение площади над кривой графика

 

функции давления на n участков

 

Если эту кривую раз ить на временные участки

t, то,

воспользовавшись формулой (3), можно вычислить расход воздуха за цикл, взяв интеграл по времени вычислив площадь замененного участка над кривой (см. р с. 33), корень квадратный и умножив на соответствующий коэфф ц ент.

цикл = K

,

(3)

или

 

 

цикл = K

,

(4)

где K = 0,004μ√ρв– постоянный коэффициент; Р – разница между давлением во впускном коллекторе и атмосферным по абсолютной величине; di – диаметр дросселя, который может изменяться; если он не изменяется, то его можно включить в постоянный коэффициент K;

tизмеряется в часах.

Вслучае, если t постоянно (не изменяется со временем), то формулу (4) можно записать в виде

36

цикл = K’ • , (5)

где K’ = К t.

Вычисление расхода воздуха по формуле (5) можно сделать программно. При этом контроллеры Адруино позволяют получить число отсчетов за цикл n около 1000, т. е. обеспечить очень высокую точность.

Теоретические зависимости между массой засосанного воздуха (или коэффициентом наполнения цилиндра) и этими давлениями и скоростями довольно сложны, но нам требуются относительные диапазоны изменения параметра, пропорционального массе засосанного воздуха, который мы перенесем на диапазон изменения параметра, регулирующего массу впрыскиваемого топлива (например, функцией map()). В случае учета нелинейности нам потребуются корректирующие коэффициенты, а также коэффициенты, учитывающие влияние внешних факторов (температуры, атмосферного давления и т.п.). Эти коэффициенты могут быть заданы аналитически илиСибАДИв виде таблиц, сохраняемых в EEPROM контроллера.

На рис. 34 показаны упрощенные кривые зависимости величины сигнала с датчика давления от времени (убраны помехи и колебания давления из-за турбулентности, амплитуда которых около 0,5 В) при изменении скорости вращения и положения дроссельной заслонки.

Рис. 34. Сигналы на входе (с датчика давления) и выходе компаратора

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

37

Схема электрическая принципиальная для превращения сигналов с датчика давления в прямоугольные импульсы представлена на рис. 35. По этим импульсам можно точно определить моменты времени открывания и закрывания клапанов.

СибАДИРис. 35. Схема компаратора для получения прямоугольных импульсов

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

Рис. 36. Датчик давления с платой компаратора, установленный на макете ДВС

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

38

задаются коэффициентом наполнения цилиндра. По определению, этот коэффициент определяется следующим образом:

Совершенство процесса наполнения оценивается коэффициентом наполнения ηV, равным отношению количества свежего заряда, находящегося в цилиндре, к началу действительного сжатия, то есть к моменту закрытия органов газообмена, к тому количеству свежего заряда, которое могло бы заполнить рабочий объем цилиндра при условиях на впуске [7].

ηV =

Gk.cж = ,

(6)

 

Vh ρk

 

где Gk.сж – массовое количество свежего заряда, находящееся в цилиндре к

началу сжатия, кг; Vk – объем, занимаемый свежим зарядом и приведенный

к условиям на впуске (ρk и Tk),

м3.

Объемы Vk и Vh отмечены на рис. 37.

Рис. 37. Полный Vh всасываемый Vk объемы цилиндра

Если ηV =1, то в цилиндр

при работе ДВС попадает ровно столько

воздуха, сколько бы уместилось в объеме, равном этому цилиндру при

данном атмосферном давлении и температуре.

Предположим, что клапаны цилиндра безынерционны и не

СибАДИ

испытывают трения. Тогда они

будут переключаться из закрытого

состояния в открытое и наоборот при атмосферном давлении с двух сторон клапана. Моменты переключения клапанов при вращении маховика и возвратно-поступательном движении поршня совпадают с моментами появления фронтов импульсов на выходе «1» компаратора (см. рис. 35 и рис. 39). Зная время в моменты переключения клапанов и время прохождения нижней (или верхней) мертвой точки, можно точно вычислить коэффициент наполнения цилиндра для нашей модели. В верхней мертвой точке начинается всасывание воздуха в цилиндр и эти моменты совпадают со спадающими фронтом импульсов Uвых1 на рис. 39.

39

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

Это соответствует нарастающему фронту импульсов на выходе Uвых2 (см. рис. 38).

Зная перечисленные моменты времени, коэффициент наполнения цилиндра можно определить как

ηV = Vk/Vh = 0,5 (1 сos a);

(7)

если угол а в градусах, то

 

a = 180(1 – ΔT1/T) = 180ΔT2/T,

(8)

если в радианах –

 

a = π(1 – ΔT1/T) = πΔT2/T,

(9)

 

 

 

 

И

 

где a – угол, показанный на рис. 37, соответствующий переключению

впускного клапана; T половина периода вращения (поворот на 180о

или

на π рад).

 

 

Д

 

T = T1 + ΔT2; ηV = 1

при ΔT1 = 0

или ΔT2 = Т.

 

 

 

А

 

 

 

б

 

 

 

и

 

 

 

 

С

 

 

 

 

 

Рис. 38. Кнопка и схема для регистрации нижней мертвой точки

 

Окончательно формула для вычисления коэффициента наполнения

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

 

 

 

 

 

ηV = Vk/Vh = 0,5 (1 сos (πΔT2/T)).

(10)

40

А

И

Рис. 39. Временные диаграммы напряжений

Ниже приведен код программы, в которой определяются скорость

б

 

вращения маховика и коэффициент наполненияДцилиндра по длительности

импульсов в последовательностях импульсов с выходов 1 и 2, показанных

 

 

и

на рис. 37 при постоянном диаметре дросселя.

 

Код:

С

//** Измерение коэфф ц ента наполнения цилиндра

/* Используемые штырьки:

А0

– вход с переменного резистора,

А1

– вход от датчика давления,

D2

– с выхода компаратора (выход 1, см. рис. 35)

D3

– с выхода RS-триггера (выход 2, см. рис. 38).

D9 – выход ШИМ

D10 – с выхода делителя на 2

*/

 

int rPin = A0;

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

int deltaT1;

// длительность импульса

int deltaT2; ;

//………..

int hpT;

 

double knc;

//коэффициент наполнения цилиндра

float pi;

//число пи

int cycleT; // время одного цикла хода цилиндра float speedR;

41

unsigned int rValue = 0; // цифровое значение напряжения с переменного резистора

void setup() { pinMode(2, INPUT); pinMode(3, INPUT);

Serial.begin(9600); // старт последовательного вывода данных pi= 3,14;

}

 

void loop() {

 

rValue = analogRead(rPin);

// считываем значение с переменного

резистора

 

rValue = map(rValue, 0, 1023, 0, 255); //масштабируем

analogWrite(9, rValue); // Pin9 - выход ШИМ

cycleT = pulseIn(10, HIGH); // Pin10 – с выхода делителя на 2 (см. рис. 25)

deltaT1 =pulseIn(2, HIGH);

 

 

 

И

 

 

 

 

deltaT2 =pulseIn(3, HIGH);

 

 

 

 

hpT = deltaT1+deltaT2;

 

 

 

Д

if(cycleT < 3*hpT)

 

 

 

 

 

 

 

 

 

 

{

 

 

 

 

 

А

 

Serial.print("Error!!!");

 

 

 

 

 

 

 

 

}

 

 

 

 

б

 

 

else

 

 

 

 

 

 

 

 

 

 

 

 

{

 

 

и

 

 

 

knc = (1-cos(pi* deltaT2 / hpT))/2;

 

 

 

}

speedR = 60000 / cycleT;

 

 

 

 

 

 

 

 

 

 

 

С

 

 

 

 

 

//Вывод на монитор

 

 

 

 

 

 

 

Serial.print("Speed = ");

 

// скорость вращения

 

Serial.print(speedR);

 

 

 

 

 

 

 

Serial.print(", Koef.nap.cil = ");

 

 

 

 

Serial.println(knc);

 

//коэффициент наполнения цилиндра

delay(100);

 

 

 

 

 

 

 

}

 

 

 

 

 

 

 

 

Данная программа выводит на экран дисплея две колонки: в первой скорость вращения, во второй коэффициент наполнения цилиндра. Это позволяет зарегистрировать зависимость коэффициента наполнения цилиндра от скорости вращения коленвала ДВС при различных диаметрах дросселя (см. рис. 31), т. е. получить семейство механических характеристик для данного макета без нагрузки (на холостом ходу).

42

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]