Алгоритм определения угла опережения зажигания при переключении с бензина на газ для контроллера электронной системы управления двигателем
Автор: Шишков В.А.
Журнал: Известия Самарского научного центра Российской академии наук @izvestiya-ssc
Рубрика: Механика и машиностроение
Статья в выпуске: 4-1 т.13, 2011 года.
Бесплатный доступ
На основании исследований, проведённых автором, предложен алгоритм определения угла опережения зажигания при переключении двигателя с бензина для работы на газовом топливе в соответствии со скоростью горения топливной смеси, зависящей от многих параметров, и приближенного к реальным физическим процессам, происходящим в камерах сгорания ДВС. Сделан ориентировочный расчет изменения угла опережения зажигания для компримированного природного газа (КПГ) и сжиженного углеводородного газа (СУГ). Показаны возможные погрешности и пути их минимизации, а также порядок реализации алгоритма при определении коэффициентов влияющих параметров.
Угол опережения зажигания, двигатель внутреннего сгорания, электронная система управления двигателем, алгоритм управления, газовое топливо, скорость горения
Короткий адрес: https://sciup.org/148200108
IDR: 148200108
Текст научной статьи Алгоритм определения угла опережения зажигания при переключении с бензина на газ для контроллера электронной системы управления двигателем
так и на переходных режимах при поставленных начальных условиях. Этими условиями могут быть максимальные значения крутящего момента и мощности двигателя или минимальный расход топлива с учетом отсутствия детонации на всех режимах работы ДВС. Обычно определяются обе эти характеристики и используются в зависимости от условий движения, нагрузки автомобиля и манеры езды конкретного водителя.
Как видно, все эти варианты требуют проведения дополнительных специальных испытаний. Но если были проведены испытания для одного вида топлива по первому пункту, то целесообразно перейти к четвёртому пункту, что позволит сократить на порядок объем проведения испы- таний при доводочных работах для другого вида топлива.
На процесс воспламенения топливовоздушной смеси влияет величина времени накопления энергии в катушке зажигания, которая может быть как одинаковой, так и различной для бензина и газового топлива. Увеличение угла опережения зажигания при работе на газовом топливе обычно приводит к уменьшению времени накопления энергии в катушке зажигания, если момент начала ее накопления неизменен по углу положения коленчатого вала двигателя. При работе на газовом топливе момент начала зажигания наступает раньше, чем при работе на бензине. В этом случае необходимо раньше начинать накопление энергии в катушке зажигания. Это выполняется для того, чтобы сохранить, а в лучшем случае и увеличить время накопления энергии для повышения энергии разряда при работе на газовом топливе.
Увеличение времени накопления энергии при работе на газовом топливе можно, например, определить в той же зависимости, что и изменение угла опережения зажигания. Повышенная энергия разряда на свече зажигания нужна потому, что для воспламенения газового топлива требуется большая энергия, чем для смеси паров бензина с воздухом.
Скорость горения зависит от химического состава топлива и, в частности, от содержания молекул водорода, который увеличивает число активных центров в начале процесса окисления. При одинаковом химическом составе топлива на скорость горения будет оказывать влияние и размер его молекул. Чем меньше молекула, тем меньше энергии необходимо на ее разложение с последующим окислением атомов водорода и углерода, а также тем проще доступ радикалов и молекул кислорода к атомам углерода. В работе [2] изложено недостаточно полное теоретическое описание горения бимолекулярных реакций, которое не может быть применено в настоящее время для контроллеров систем управления двигателем внутреннего сгорания.
Кроме этого, скорость горения зависит также и от качества приготовления смеси – от однородности топливовоздушной смеси, от степени турбулизации потока в камере сгорания ДВС и других параметров. В работах [3-5] приведены математические модели расчетов параметров при горении топлив. Они требуют значительного машинного времени для расчета параметров и неприемлемы для современных электронных блоков управления двигателем внутреннего сгорания, так как применение более мощных процессоров приведет к их удорожанию. Поэтому требуется более простой алгоритм, описываю- щий основные физические процессы при воспламенении и горении топливной смеси.
Для двигателей с одинаковыми степенью сжатия и степенью наполнения цилиндров топливовоздушной смесью давление в конце процесса сжатия как бензиновоздушной, так и газовоздушной смесей можно считать одинаковым, так как объем поступившей смеси примерно одинаков. Но он отличается по количеству воздуха и топлива, поступившего в цилиндр.
Скорость нарастания давления в процессе сгорания зависит от количества активных центров горения или количества водорода в смеси. Чем выше эффективное давление в процессе сгорания, тем ниже будет скорость горения. Но чрезмерное увеличение давления в камере сгорания может привести к детонации, когда скорость горения на два порядка выше обычной скорости горения.
Возникает вопрос: насколько отличается угол опережения зажигания при работе ДВС на бензине, СУГ или КПГ? Угол опережения зажигания – это функция многих параметров. Выделим основные параметры, которые оказывают на него максимальное воздействие, и пренебрежем остальными, влияние которых незначительно.
α з = f ( α ; T; P; W гор ; N Н2 ; V турб ;n; G в ; V кс и т.д.) , (1) где α з – угол опережения зажигания;
-
α – состав смеси;
T – температура топливной смеси перед началом воспламенения, зависит от начальной температуры всасываемой топливной смеси, а также от степени ее сжатия и температуры стенок цилиндров, то есть от температуры двигателя;
Р – давление топливной смеси перед началом воспламенения и в процессе горения, изменяется с течением времени, то есть Р = f(t) ;
Wгор – скорость горения топливовоздушной смеси;
NН2 – количество атомов свободного водорода перед началом горения;
Vтурб – турбулизационная составляющая скорости топливной смеси перед воспламенением и в процессе горения, изменяется с течением времени, то есть Vтурб = f(t) ;
n – частота вращения коленчатого вала (КВ) двигателя;
Gв – расход воздуха;
Vкс – объем камеры сгорания, м3 (для расчета примем 0,0004 м3).
Температуру двигателя можно не учитывать, так как ее влияние определяется через температуру топливной смеси, которую можно выразить зависимостью от температуры двигателя.
Преобразуем зависимость (1) к другому виду:
α з = А(n,Gв) f (Wfгор; n; Gв ) , (2)
где W fгор = f (W н ; α ; T; P; N Н2 ; V турб ; V кс ; К т ; К о2 ; t г ) (3)
– функция скорости горения топливовоздушной смеси в зависимости от влияющих параметров;
Wн – скорость нормального горения топливовоздушной смеси, определяется эксперимен- тальным путем для каждого вида топлива при стехиометрическом составе смеси по формуле Аррениуса, нормальном давлении, отсутствии в топливе свободного молекулярного водорода, а также при ламинарном движении топливовоздушной смеси и свободном объеме;
А(n,Gв) – матрица коэффициентов коррекции угла опережения зажигания для заполнения при проведении испытаний по калибровке параметров по конкретному двигателю (в начале проведения работ все ячейки матрицы заполнены единицами);
Кт – концентрация паров топлива в смеси в процессе горения, изменяется с течением времени, то есть Кт = f(t) ;
КО2 – концентрация кислорода в смеси в процессе горения, изменяется с течением времени, то есть КО2 = f(t) ;
tг – время горения топливовоздушной смеси.
Для упрощения зависимости (3) введем функцию концентрации топлива и кислорода в сме- си от времени fк = f(Кт; КО2; tг) .
Тогда выражение (3) можно представить в
следующем виде:
Wfгор = f (Wн; б; T; P; NН2; Vз; Vтурб;Vкс; fк) . (5) В более полном виде скорость горения топ- ливовоздушной смеси:
Wfго = (А e - E/(R T) +Vт б) fб P m (1+ + lnр (NН2)) ( Кva ln(Vукрс+Kvb)) fк ,
где А – постоянный коэффициент реакции горе- ния;
Vтурб – турбулентная составляющая скорости горения топливной смеси (для простоты расчета примем одинаковой для всех трех видов топлив, равной 20 м/с);
f α – функция влияния состава смеси, в первом приближении ее можно принять в виде функции синуса f α = sin β , где β = 90· α ;
Е – энергия активации при окислении кислородом воздуха для паров бензинов 78,5-80 кДж/моль, для КПГ 86-89 кДж/моль, для СУГ 82,5-84 кДж/моль , для водорода примерно 67 кДж/моль;
R – газовая постоянная, 8,31441 кДж/(кмоль К);
е – основание натурального логарифма, величина постоянная, 2,71828…;
ln – натуральный логарифм;
Кva – множитель влияния масштабности системы, для всех трех видов топлив примем равным 6;
Кvb – коэффициент масштабности системы, при- мем одинаковым для всех трех видов топлив 1,181;
Кva ln(Vкс+Kvb) – функция влияния масштабности системы на процесс воспламенения.
Функцию концентрации топлива и кислорода в смеси в первом приближении можно взять для нормальной скорости горения. Реально функция концентрации топлива и кислорода в смеси зависит также от скорости горения, поэтому необходимо сделать несколько последовательных приближений.
Как видно из уравнения (6), функция скорости горения топливовоздушной смеси зависит не только от параметров ее состояния, но и от времени протекания процесса, то есть изменяется с течением времени. Решение уравнения (6) в нестационарном процессе представляет собой сложную задачу, которая в настоящее время неприемлема для систем управления двигателем внутреннего сгорания, так как недостаточно быстродействие процессора контроллера для оперативного вычисления в процессе работы двигателя. Для упрощения решения можно задаться некоторыми условиями, которые позволяют решить уравнение (6) для стационарного процесса при некоторой средней скорости горения или ее максимальном значении.
Угол опережения зажигания для различных частот вращения КВ ДВС:
з газ з бензин w ,
где K w = (W fгор a / W fгор г ) .
Подставим выражение (6) для скорости горения бензина и газа в выражение Kw и сократим одинаковые параметры:
! [(А бенз e ER T +V tp6 ) f a P m (l+hW К hV+K b )).Я бн !
K w =
[(А газ e ER T) +У турб ) f a P m- (1+1п(Ы н= ))(К, а ln(V K +K b )) fj a
Нормальная скорость горения паров бензина и газа, энергия активации Е , начальная температура топливной смеси (из-за отсутствия энергии на испарение капель топлива при работе на газе), влияние давления (показатель степени), скорость движения заряда и турбулентная составляющая (из-за различного объемного содержания топлива), коэффициент масштабности – все эти параметры различны для бензовоз-душных и газовоздушных топливных смесей.
Примем в первом приближении следующие условия:
- функции fб влияния состава смеси на скорость горения одинаковы для паров бензина и газа;
-
- начальная концентрация молекул водорода в бензине и газовом топливах одинакова;
-
- функции концентрации fк топлива и кислорода в смеси в процессе горения от времени одинаковы;
-
- давление Р в конце процесса сжатия на всех видах топлива одинаково;
-
- турбулентную составляющую скорости V
турб для известной конструкции камеры сгорания и системы впуска топливовоздушной смеси при работе на бензине и газе считаем одинаковой.
Тогда уравнение (8) можно переписать в следующем виде:
( тбенз ~ тгаз ) Г [ (А бенз e +^ турб) ( Kva ln(V.с +Kvb, )) ] бен]
-
K.. Р
[ [ (А газ e - E/(RT)+V myp6 ) (К ln(V kC +K vb ))] газ ]
(9) где mбенз – показатель степени влияния давления при работе двигателя на бензине при турбулентном горении ( mбенз = -0,1 ... -0,4) примем - 0,2 ;
mгаз – показатель степени влияния давления при работе двигателя на газовом топливе примем из условия его пропорциональности детонационному числу топлива. Тогда для СУГ он будет равен mгаз с г = -0,2 · 95/105 = -0,18095, для КПГ mгаз кпг =у -0,2 · 95/120 = -0,15833. В данном расчете принято детонационное число бензина – 95, СУГ – 105, КПГ – 120. Например, давление в конце процесса сжатия 1200 кПа, тогда коэффициент влияния давления на угол опережения зажигания при работе на КПГ составит 0,744219, а при работе на СУГ 0,873658. То есть угол опережения зажигания при переходе с бензина на газовое топливо из-за различного влияния давления уменьшится на 25,57% для КПГ и на 12,63% при работе на СУГ;
Eбенз , Eгаз кпг , Eгаз с г – энергии активации для паров бензина, КПГ у , СУГ с воздухом приведены выше;
Tбенз и Tгаз – отличаются на значение, обусловленное энергией парообразования капелек бензина, то есть начальная температура смеси воздуха и газовых топлив незначительно выше, чем для смеси бензина и воздуха. Например, зная массу воздуха и топлива, поступившего в цилиндр, а также количество энергии (по изменению энтальпии смеси воздуха с топливом), затраченное на парообразование всего количества бензина для прогретого двигателя или его части для непрогретого двигателя, можно определить разницу (Tгаз – Tбенз) . Пример расчета для двигателя 1,5 л, мощностью 54 кВт при n = 5600 об/ мин и составе стехиометрической смеси, равном 1, приведен ниже в табл. 1.
В табл. 1: Токр – температура окружающего воздуха;
Gц возд – цикловая подача воздуха в двигатель;
Gц топл – цикловая подача топлива в двигатель; r – теплота парообразования бензина (среднее значение для различных составов бензина).
Примем разность температур топливной смеси при работе на бензине и газе, равной изменению температуры смеси бензина и воздуха до и после испарения бензина:
газ бенз нач кон , где примем в первом приближении Тнач = Токр ;
(Србн G топл+ С рвн Gц возд ) Т нач r Gц топл
Т кон = ------------------------------------------ , (10)
(С рбк С ц топл + С рвк Сц возд )
где Србн и Србк – начальная и конечная теплоемкость бензина соответственно, кДж/(кг·град), примем среднее значение;
Срвн и Срвк – начальная и конечная теплоемкость воздуха соответственно, кДж/(кг·град), примем среднее значение.
Как видно из таблицы, разность температур газовой и бензиновой смесей с воздухом для различных условий работы двигателя составляет примерно 20-21оС. По этим значениям можно определить предельные значения коэффициента влияния температуры топливной смеси и энергии активации на угол опережения зажигания. Как следует из формулы (9), при переходе с бензина на газ температура смеси для испаренного в газовом редукторе СУГ и для КПГ влияет одинаково. Отличия коэффициентов влияния температуры топливной смеси и энергии активации для СУГ и КПГ связаны с их различными энергиями активации, которая, в свою очередь, зависит от температуры топливной смеси в процессе горения. При повышении температуры энергия активации снижается.
Определим влияние температуры и энергии активации, например, для степени сжатия 11. Зададим температуру в конце сжатия Tбенз = 750 К , 800 К , 825 К и Тгаз = 770 К , 820 К , 845 К, коэффициенты Абенз и Агаз примем равными 7000000 . Вычислим коэффициент влияния для КПГ и СУГ:
K w =
[(А бенз e E/(R T +V myp^X K va ln(V кс +K vb))]бенз [(А газ e — E/(R T) +Vm yp6 )(K va ln(VK c +K vb ))] газ
.
Диапазон изменения коэффициента Kw’ влияния температуры, энергии активации и масштабной функции составляет при Тгаз = 770 К, 820 К,
Таблица 1. Определение разности температур газовой и бензиновой смесей, поступивших в цилиндр
Т окр , °С |
vц возд , г/цикл |
ц топл , г/цикл |
г, кДж/кг |
С рб , кДж/(кг°С) |
рв , кДж/(кг°С) |
Т газ - Т бенз , °С |
-30 |
0,3702 |
0,0254 |
348 |
1,73 |
1,0048 |
21,25 |
+45 |
0,2851 |
0,0195 |
328 |
1,72 |
1,0056 |
19,97 |
845 К соответственно для КПГ 1,467; 1,712; 1,671 , а для СУГ 1,1395; 1,223; 1,189. Функция данного влияния имеет экстремум из-за различного влияния температуры и энергии активации. Как видим, увеличение угла опережения зажигания при переходе с бензина на газовое топливо в зависимости от температуры, энергии активации топлива и масштабного коэффициента составит для КПГ 46,7-71,2%, а для СУГ 13,95-22,3%.
Тогда суммарное влияние всех параметров, а именно давления и коэффициента Kw’ по соотношению (11), на увеличение угла опережения зажигания по уравнению (9) составит для КПГ Kw’ = 1,092-1,274, а для СУГ Kw’ = 0,996-1,068. Например, угол опережения зажигания на бензине был 25 град, тогда на КПГ составит 27,331,85 град, то есть фактическое увеличение составит от 2,3 до 6,85 град по углу положения коленчатого вала двигателя. Аналогично для СУГ он составит 25-26,7 град, то есть фактическое увеличение составит от 0 до 1,7 град по углу положения КВ двигателя.
Величина погрешности вычисления изменения угла опережения зажигания при переходе с одного вида топлива на другой складывается из следующих погрешностей:
-
- экспериментального подбора коэффициентов в формуле Аррениуса и эмпирических зависимостях, используемых в данном алгоритме;
-
- выполнения калибровочных работ на конкретном двигателе при подборе коэффициентов влияния для получения максимальной мощности, максимального крутящего момента и минимальной токсичности отработавших газов;
-
- измерения или вычисления температуры и давления в камере сгорания перед процессом воспламенения;
-
- задания состава топливной смеси (погрешность измерения и вычисления расхода воздуха и топлива);
-
- определения турбулентной составляющей скорости в процессе горения;
-
- определения количества атомарного водорода в топливной смеси и т.д.
При осмысленной постановке экспериментов и выполнении калибровочных работ все эти погрешности можно минимизировать до необходимого уровня для получения физической картины протекающих процессов.
Предлагается такой порядок определения необходимых для расчета коэффициентов. Последовательно определяются:
-
. нормальная скорость горения паров бензина и газа с воздухом при одинаковых условиях;
-
. коэффициенты влияния состава смеси на нормальную скорость горения при одинаковых условиях;
-
. показатель влияния давления на нормальную скорость горения при одинаковых условиях для смесей паров бензина и газа с воздухом;
-
. по нормальной скорости горения - энергия активации и коэффициент А в формуле Аррениуса;
-
. количество атомарного водорода в топливе по химическому составу топлива или по расходу дополнительного водорода;
-
. турбулентная составляющая скорости (определяется оптическим методом);
-
. коэффициенты в функции влияния объема камеры сгорания на нормальную скорость горения;
. матрица коэффициентов А(n,Gв) (определяется в процессе проведения калибровочных работ электронного блока управления двигателем); при правильно выбранных коэффициентах в формуле (9) коэффициенты матрицы А(n,Gв) будут близки к единице.
Для ускорения процесса адаптации часть коэффициентов и функций, используемых в данной методике, можно найти в соответствующей специализированной литературе по теории горения паров бензина и газовых топлив в смеси с воздухом. Часть коэффициентов в первом приближении можно задать, при этом не нарушая физических процессов, а погрешности их задания компенсировать матрицей калибровочных коэффициентов А(n,Gв) . В случае использования более мощного процессора в электронном блоке управления двигателем алгоритм вычисления угла опережения зажигания при переходе с одного вида топлива на другое можно расширять введением алгоритмов вычисления соответствующих параметров влияния.
ЗАКЛЮЧЕНИЕ
-
1. Данная работа позволяет определить угол опережения зажигания при переходе с бензина на газовое топливо с учетом физического влияния основных факторов.
-
2. Предложенный алгоритм изменения угла опережения зажигания можно использовать в электронном блоке управления в системах управления ДВС, работающего на газовом топливе.
Список литературы Алгоритм определения угла опережения зажигания при переключении с бензина на газ для контроллера электронной системы управления двигателем
- Злотин Г.Н., Захаров Е.А., Кузьмин А.В. Регулировка бензинового ДВС при переводе его на сжиженный нефтяной газ//Двигателестроение. 2007. №2. С.29-31.
- Гардинер У., Диксон#Льюис Г., Цельнер Р. и др. Химия горения: Пер. с англ [под ред. У.Гардинера]. М.: Мир, 1988. С. 464.
- Zaslavsky M.Yu., Pergament A.Kh., Plushchenkow B.D. Dynamics and Stability of One Dimensional Combustion Problems. Preprint, Inst. Appl. Math. Naim M.V. Keldysh, the Russian Academy of Science. Moscow, 2002.
- Сабденов К.О. Теплофизические и гидрогазодинамические эффекты при горении газов и ракетных топлив. Автореф. дисс.... докт. физ. мат. наук. Томск, 2007. С. 41.
- Матиевский Д.Д., Сеначин П.К., Свердлов М.Ю., Ильина М.А. Новый подход к проблеме моделирования сгорания смеси в ДВС с искровым зажиганием. 20 c. URL: http://aomai.secna.ru:8080/Books/Files/1999-02/HTML/12/pap_12.html (дата обращения 18.01.2011)