CC BY
40
Для каждого утвержденного типа транспортного средства в системе качества предприятия-изготовителя должны быть предусмотрены документированные процедуры периодических проверок или испытаний по оценке соответствия выпускаемой продукции требованиям нормативных документов, используемых при подтверждении соответствия. Для выполнения этого требования изготовитель должен:
• обеспечить наличие процедур эффективного контроля соответствия транспортных средств официально утвержденному типу;
• иметь доступ к необходимому контрольному оборудованию для проверки соответствия каждого официально утвержденного типа;
• обеспечить регистрацию и хранение данных о результатах испытаний;
• анализировать результаты испытаний и принимать меры по поддержанию стабильности характеристик продукции;
• обеспечить проверку всех официально утвержденных типов транспортных средств;
• обеспечить в случае обнаружения несоответствия производства на любой выборке продукции, проведение новых испытаний и принятие необходимых мер для восстановления соответствия производства.
Заключение
Сформулированные рекомендации по оценке безопасности и модели рисков могут найти применение в автотранспортном комплексе в следующих практических случаях:
• для обнаружения опасных сочетаний факторов риска при оценивании безопасности эксплуатации высоконадежных систем;
• при обосновании эксплуатационных требований по безопасности к конструкциям автомобилей и системам, при оценивании приспособленности АТС к условиям эксплуатации.
Список литературы
1. Комаров В.В. Методические основы оценки безопасности автотранспортных средств с помощью моделей рисков «катастоф» [Текст]/ В.В. Комаров, Е.А. Куклев // Автомобильная промышленность. - 2008. -№ 5. - С.26-29.
2. Комаров В.В. Методология оценки безопасности автотранспортных средств по техническому уровню и сроку эксплуатации [Текст]/В.В. Комаров//Известия МГТУ «МАМИ». -2007. - № 2(4).- С.114-122.
3. Комаров В.В. Интегральные критерии безопасности автотранспортных средств [Текст]/ В.В. Комаров // Автотранспортное предприятие. -2007.- № 12. С.46-48.
4. Рябинин И.А. Надежность и безопасность структурно-сложных систем [Текст]/ И.А. Ря-бинин -СПб.:Политехника, 2000.-248 с.
5. Куклев Е.А. Оценивание рисков на основе цепей случайных событий [Текст]/ Е.А. Куклев - М.:РАТИ, 2003.- Серия «Наука и техника».
6. ISO/IEC Guide 2: 1996 Standardization and related activities. General Vocabulary, p.9.
7. Малинецкий Г.Г. Управление риском. Риск. Устойчивое развитие. Синергетика [Текст]/ Г Г. Малинецкий, В.В. Кульба, С.А. Косяченко, М.Г. Шнирман - М.: Наука, 2000,- 431с.-Серия «Кибернетика», РАН.
8. Хэнли Э. Надежность технических систем и оценка риска [Текст]: пер. с англ./ Э. Дж. Хэнли, Х. Кумамото. - Машиностроение, М.: 1984. -526 с.
Исследование кинетического механизма окисления азота в автомобильном двигателе с послойным смесеобразованием
д.т.н. проф. Фомин В.М., Платунов А.С.
МГТУ «МАМИ» 8 (495) 369-90-48, mixalichDM@mail.ru
Аннотация. Предлагается методика расчета, которая позволяет определить локальные образования оксидов азота в камере сгорания автомобильного бензино-
вого двигателя с внутренним смесеобразованием для режимов с послойной организацией заряда. С учетом резко неоднородного состава расслоенной смеси рабочее пространство камеры сгорания разделяется по пути распространения фронта пламени на зоны с различными (локальными) коэффициентами избытка воздуха. Для двигателя типа AUX фирмы Volkswagen проведены тестовые расчеты, по результатам которых установлено, что выход оксидов азота при сгорании расслоенной смеси снижается в 2,2 раза по сравнению с традиционным вариантом однородной гомогенной смеси. Расхождение результатов расчета с экспериментом не превышает 8%, подтверждая тем самым удовлетворительную адекватность расчетной методики.
Ключевые слова: автомобильный двигатель, непосредственное впрыскивание бензина, послойное смесеобразование, гетерогенная смесь, гомогенная смесь, расчетная методика.
Основным стимулом, определяющим постоянно растущее применение на легковом автотранспорте двигателя с непосредственным впрыскиванием бензина (НВБ), является его высокий уровень топливно-экономических показателей, которые достигаются благодаря возможности создавать в камере сгорания (КС) систему горения с послойным распределением смеси.
При работе автомобильного двигателя с НВБ на режиме послойного смесеобразования обеспечивается высокая полнота сгорания заряда, и суммарная токсичность продуктов сгорания определяется, главным образом, содержанием в них оксидов азота (в основном -
NO- ) [1].
Общепринятое математическое описание процесса сгорания в бензиновом ДВС предполагает наличие однородных полей температуры, коэффициента избытка воздуха а, а следовательно, и концентрации компонентов смеси в объеме КС. Такая идеализация процесса суживает диапазон возможностей по поиску мер снижения эмиссии токсических веществ, в частности, оксидов азота. В условиях системы горения расслоенных смесей целевой задачей
моделирования процесса образования NO- является выявление пространственного поля температур и коэффициента избытка воздуха а в объеме КС и его влияния на этот процесс. Отсюда вытекает необходимость разработки нового подхода к математическому представлению полей температур, а и горения в этих полях.
Для кинетики окисления азота характерна сильно выраженная нелинейная зависимость от соотношения топлива и воздуха в исходной смеси. В связи с этим важно рассмотреть возможность решения пространственных задач горения для изучения особенностей локальных кинетических процессов окисления азота, происходящих в отдельных областях (зонах) КС с различным коэффициентом избытка воздуха в смеси.
С учетом резко неоднородного состава (расслоенных) смесей предлагается разделить весь объем КС по пути распространения сферического фронта пламени на зоны с различными (локальными) коэффициентами избытка воздуха, при этом количество выделенных расчетных зон может быть принято практически неограниченно. Область распространения фронта пламени охватывает последовательно одну за другой зоны выделенных объемов КС.
Подобный подход позволяет решать задачу по определению локальных температур и кинетики окисления азота для каждой из зон (по мере подхода фронта горения) и получить представление о распределении концентрации оксидов азота по объему КС.
Такая задача может быть традиционно описана уравнением энергии для определения температуры зоны по тепловым эффектам процесса сгорания топлива с учетом тепловых потерь через локальные поверхности теплообмена стенок КС. Далее по полученным данным
NOx
можно спрогнозировать состав продуктов сгорания и уровень эмиссии x .
Из зарубежных публикаций известно [2, 3], что на режиме работы ДВС с послойным
Раздел 1. Наземные транспортные средства, энергетические установки и двигатели. распределением заряда процессы впрыска, испарения, перемешивания и горения в небольшом объеме камеры сгорания протекают интенсивно и, практически, одновременно. Поэтому при разработке модели было принято, что продукты сгорания оказываются полностью перемешанными уже в начале расширения [2]. В результате во время расширения в цилиндре находится однородно перемешанная (гомогенная) смесь.
Продолжительность интервала активного горения и время его подхода к каждой выделенной зоне КС могут быть определены по данным предварительного расчета закона тепловыделения (в варианте - на основе обработки индикаторной диаграммы), после чего внутри этого интервала можно выделить макроточки начала горения в центрах каждой из выделенных зон.
Учитывая сильную нелинейность зависимости скорости химических реакций окисления азота от температуры, необходимо использовать не среднюю («индикаторную») по объему КС температуру, а ее локальное значение в выделенной зоне. Очевидно, что именно локальная температура в макроточке отдельной зоны определяет местную кинетику процесса
образования . Расширяя, в первом приближении, макроточку до размеров зоны считаем, что температура внутри каждой зоны одинакова. Это позволяет решать задачу кинетики окисления азота для центральной точки каждой из выделенных зон (по мере подхода фронта
N0
горения) и получить представление о распределении концентрации х по объему КС и её изменении вплоть до открытия выпускного клапана.
Т
Значение локальной текущей (по времени тепловыделения) температуры л для каж-
дой выделенной зоны КС с различным локальным коэффициентом избытком воздуха смеси можно рассчитать согласно зависимости, предложенной в работе [4]:
а
ёТ
л
ёт
1
ш„
ди л
дТ„
и
1--л-
Н
ёу
ёт ёт
■ р--т
ёт
ди л
л
да
ёт
( \ ётвп ётвьш ивп - ил + КвпТвп )—:--
ёт
ёт
Зависимость выведена на основе рассмотрения уравнения баланса энергии с учетом характера процесса тепловыделения (сгорания) и тепловых потерь через поверхности теплообмена стенок КС. Следует отметить, что приведенная выше зависимость, в отличие от исходной [4], упрощена: изменение внутренней энергии рабочего тела локальной зоны ил в зависимости от давления не учитывалось (с допускаемой ошибкой не более 1,5%, [4]), и она рассматривалась только как функция локальных значений температуры и коэффициента избыт-
ил ?(Тл алЛ
ка воздуха: л = I ( л, л).
Процесс сгорания происходит при закрытых органах газораспределения, поэтому:
ёт„„ ёт
' вп
вып
ёт
ёт
0
; тогда окончательно имеем:
ёТл 1
ёт ди л
тл
л дТ
и
Л
Н
ёЯх ёО^
и /
ёт ёт
Р
ёу_
ёт
т.
ди л Кдал У
ёа
ёт
(1)
тп г „ тл = т + тлв
где: л - масса рабочего тела в локальной зоне; ^ ;
т
л.ц.
т
- масса топлива, впрыснутого в локальную зону за цикл;
- масса воздуха локальной зоны.
в
1
- количество теплоты, отведенной от рабочего тела локальной зоны;
Ни
и - низшая теплота сгорания топлива. Для определения текущей температуры рабочего тела отдельной локальной зоны
Т (Т) а О (Т)
л у ' при известном значении л необходимо задать закон тепловыделения х 4 . Изменение объема цилиндра ^^ / определяется с учетом конструкции двигателя и кинематической схемы его КШМ:
V, пП2
V = + (т)
1 4
п V,
где: п, " - диаметр и рабочий объем цилиндра;
5х (т) - кинематическая функция перемещения поршня.
„ т = тц + тк
Масса рабочего тела для всего объема КС ц В для конкретного ДВС может
быть определена по опытным данным (без учета массы остаточных газов) следующим образом. Располагая опытным значением массового секундного расхода воздуха через двигатель
О т т
в , с учётом частоты рабочих циклов и числа тактов за цикл ц , масса воздуха В , поступившего в цилиндр, определяется из формулы:
30 • О В Тц
т =-
1Пф
пр
тц О
а цикловая подача топлива ц при известном часовом расходе топлива Т :
Ст • ТД
тц =-—
ц 120•п•\ .
Здесь I — число цилиндров; п — частота вращения вала двигателя; ^пр. — коэффициент продувки (для двигателя с НВБ может быть принят, равным 1).
С учетом того, что по определению объемы зон одинаковы, массы содержащегося в них
т„ / 2
воздуха соответственно будут равны В , где: ъ - число выделенных зон.
а
В соответствии с предварительно заданным составом смеси в каждой зоне ( л ) распределение массы впрыснутого за цикл топлива по выделенным зонам КС к началу сгорания т
лц. определится по формуле:
1
т тв'-
л.в 2
т =-=-—
/ / ал 10 ал • 10
Таким образом, масса рабочего тела в отдельно взятой локальной зоне будет равна:
т = т + т „ = т +1/г • т„
л л.ц. л.В. л.ц. В
После поэтапного (позонного) расчета изменения температуры по уравнению (1) можно определить изменение давления, используя уравнение состояния Клапейрона:
р(т) = Т. Кл (т)тл (Т) Р Т'= V (т) .
В соответствии с зависимостью (1) температурное состояние рабочего тела в каждой
Раздел 1. Наземные транспортные средства, энергетические установки и двигатели. выделенной зоне в значительной степени зависит от интенсивности двух конкурирующих процессов: тепловыделения и теплоотвода. В этих условиях особую значимость приобретает точность определения теплоты, воспринимаемой теплообменными поверхностями стенок КС, соприкасающихся с рабочим телом отдельной зоны. Все известные формулы для расчета нестационарного, но усредненного по поверхности КС коэффициента теплоотдачи содержат усредненную по объему КС нестационарную температуру. При принятом разделении объема КС на зоны с различным составом смеси требуется рассчитать эту температуру и теплообмен обособленно для каждой зоны, что несколько усложняет задачу.
Изменение суммарного теплового потока для всех зон в общем виде запишем в виде уравнения теплового баланса:
йт 1 йт
Количество теплоты, переданной через поверхность стенок КС рабочим телом одной из локальных зон, может быть определено традиционно численным интегрированием уравнения Ньютона-Рихмана:
А0^ = а(тл -КАг,
р Т
где: л№ , № - локальная площадь и температура тепловоспринимающей поверхности для отдельной выделенной зоны КС;
Т
Л - текущая температура рабочего тела в объеме локальной зоны;
- коэффициент теплоотдачи. Точное определение текущих площадей теплообмена для каждой из выделенных зон представляет собой трудную, еще не решенную в теории теплообмена ДВС задачу. Поэтому в первом приближении можно принять, что величины локальных площадей теплообмена
р / р = V / V
пропорциональны соответствующему объему локальной зоны: л № № л ; или с учетом того, что объем КС был разделен изначально на ъ равных частей, теплообменные площади для каждой локальной зоны равны:
р =1/г■р л .
р
Суммарная площадь тепловоспринимающей поверхности КС находится как сумма площадей поверхностей головки поршня, головки цилиндра и цилиндра:
= рп + рГ + рЦ
Площадь теплообменной поверхности цилиндра является функцией времени:
рц М = ^ + пОБх (г)
где - объем пространства сжатия цилиндра.
Одним из важнейших параметров, который определяет условия теплообмена в КС дви-
а
гателя, является текущий коэффициент теплоотдачи . Для его определения целесообразно использовать зависимость проф. Р.З.Кавтарадзе [5], с помощью которой, на наш взгляд, можно наиболее корректно описать условия нестационарного теплообмена в отдельно взятой локальной зоне КС.
На основе решения уравнения Фурье-Кирхгофа, описывающего изменение температуры по толщине пограничного слоя, с учётом ряда допущений автором было получено выра-
а
жение для коэффициента теплоотдачи № [5]:
а =
ЬЛ.
4П
1 + 24т Ах Ср (( - Т№) * АТ
(2)
Ьл -\1^рт.л ' Ср.л ' рл
л рт.л р.л г л 11
где: у - коэффициент проникновения теплоты для пограничного слоя
теплообменной поверхности локальной зоны;
к к = 0,000361 • Т0'75
рт.л - теплопроводность рабочего тела в локальной зоне: рт.л ср .
С
р.л - изобарная теплоемкость рабочего тела в локальной зоне:
Вт/(м-К);
л —
С = 0,0762 • Т + 837,4
(для периодов сг
Рл =
Кл • Тср
' р.л ср '' " (для периодов сгорания-расширения);
Р
Рл = —
Рл - плотность рабочего тела локальной зоны:
К
л - газовая постоянная для рабочего тела локальной зоны;
Т
л и р - значения локальной температуры и давления рабочего тела. Теплофизические свойства рабочего тела локальной зоны в каждом интервале времени
Т = 0,5(Т + Т )
рассчитывались по средней температуре пограничного слоя ср . Относи-
Ах
тельная доля тепловыделения в локальной зоне для каждого расчетного шага определяется с учетом принятого при моделировании закона тепловыделения.
Согласно выражению (2), коэффициент теплоотдачи рассматривается как функция тепловыделения. Это позволяет исследовать взаимосвязь процессов выделения тепла и теплоотдачи в каждой отдельно взятой локальной зоне КС.
При выводе (2) учтено, что в течение малого (расчетного) промежутка времени Ат температура рабочего тела изменяется незначительно. Это позволяет считать теплофизиче-
ские свойства рабочего тела в данном интервале Ат постоянными величинами. Таким образом, выражение (2) позволяет определить значение коэффициента теплоотдачи для каждой
Т, Т +1 =Т + Ат]
выделенной зоны в интервале времени .
Численное интегрирование зависимости (1) позволяет определить значение локальной температуры рабочего тела для каждой выделенной зоны. С учетом принятого значения а равновесный состав продуктов сгорания для каждой зоны может быть определен по известным методикам, например, [6].
По полученным данным производится расчет локальной концентрации N° в составе продуктов сгорания каждой зоны с учетом следующих соображений.
В вы
азота N°. В зависимости от процессов, приводящих к образованию оксидов азота, различают следующие виды : термические , быстрые и топливные .
Термические образуются в результате окисления атмосферного азота со свободным кислородом в зоне продуктов сгорания при высокой температуре. Быстрые образуются непосредственно в зоне горения углеводородных топлив в результате связывания молекул азота радикалами и . Их содержание в продуктах сгорания бензина пренебрежительно мало. Топливные образуются в процессе горения в результате частичного окисления до NO азотсодержащих соединений топлива. При сгорании бензинов, получае-
N°
В выпускных газах двигателей более 90% всего количества х составляет монооксид
Раздел 1. Наземные транспортные средства, энергетические установки и двигатели. мых из нефти, в силу малости концентраций этих соединений в топливе можно также пренебречь актом образования N° из топливного азота.
К настоящему времени механизм образования окиси азота в пламёнах изучен достаточно подробно. Общепринятой теорией образования N° из атмосферного азота и кислорода является термическая теория, разработанная Я.Б.Зельдовичем, П.Я. Садовниковым и Д.А. Франк-Каменецким [7].
Применительно к условиям бензинового двигателя обычно используют модель окисления азота, включающую в себя в качестве основных следующие реакции цепного механизма [8]:
(3)
N + NO ^ N2 + O
N + °2 ^ N° + О (4)
N + ОН ^ N° + Н. (5)
Известно [8], что в богатых смесях, в условиях нехватки азота и кислорода, выход ^° определяется кинетикой разложения и температура «замораживания» снижается. Поэтому для богатых смесей определяющими реакциями являются (3) и (5), т.е. допустимо использование упрощенной модели с использованием этих реакций. В случае бедных смесей «замораживание» содержания окиси азота происходит на уровне максимальной концентрации
NO
. И лишь для больших значений а мож-
смеси, что обусловлено кинетикой образования но ожидать преобладания реакции (4).
Учитывая широкий диапазон изменения состава смеси по отдельным выделенным зонам КС двигателя с НВБ, необходима обобщенная модель. С этой целью система реакций (3) - (5) дополняется бимолекулярной реакцией [8]:
N2 + О2 ^ N° + N° (6)
Выражения для констант прямых реакций (3) - (6) приведены в табл. 1 [1].
Константы скорости прямых реакций:
Таблица 1.
Реакции^ к^.л /(моль- е)о
(3)о kiH 1,32 1О10п
(4)о к2н 1,81 Ю^Т1-5 ехр (-3000 / Т)о
(5)о кзн 4,2 1О10й
(6)о кфз 9 1 1021 Т-2^ ехр (-64800-/-Т)о
Следуя [2], уравнение для расчета концентрации окиси азота запишется в виде:
dx
dz
где:
ß
XNO ! rNO
- = -ß(5S, + S2 + S3) + S + 5(S2 + S3) + (1 -ß2)S4 5 = (S + ß)/(S ß +1) S = Si/(S2 + S3) Si = кi • rm • rn.
(7)
х и г - неравновесные и равновесные мольные доли №; индексы т, п обозначают компоненты смеси в уравнениях (3) - (6). Используя принцип квазистационарности промежуточных компонентов, уравнение [1]:
получено
dx
NO
dr
■К,
2S1
1 + ßS
+ S,
(8)
Равновесные доли r для каждой зоны могут быть определены по результатам традиционного термодинамического расчета состава продуктов сгорания при заданном локальном значении коэффициента избытке воздуха для текущих значений давления и температуры [6]. Применительно к двигателю AUX фирмы Volkswagen проведена апробация данной ме-
Раздел 1. Наземные транспортные средства, энергетические установки и двигатели. тодики расчета. Как известно, в двигателях с НВБ в зависимости от условий эксплуатации предусмотрены режимы работы с различными способами организации рабочего процесса. На высоких нагрузках ДВС работает на гомогенных смесях стехиометрического состава, а на частичных - на обедненных смесях с расслоением заряда. Для указанного ДВС в соответствии со штатной программой управления (Motronic MED 7.1.1. фирмы Bosch) установлен предельно максимальный режим работы с послойной организацией заряда, который реализуется
Р
при значении среднего эффективного давления е = 0,4 МПа при частоте вращения вала ttd/^ n = 3000мин-1 ^
ДВС . Тестовые расчеты проводились для данного режима, который характе-
NOx
ризуется максимальным уровнем эмиссии x .
Предполагалось, что для более наглядного понимания влияния характера расслоения
NOx КС
заряда на эмиссию x пространство КС двигателя достаточно разделить на три характерные локальные зоны (рисунок 1). Оксиды азота образуются в трех выделенных зонах (Н, С и К) КС с различными значениями коэффициентов избытка воздуха (рисунок 1). Объемы выделенных зон одинаковы.
Для выбранного расчетного режима среднее по объему КС значение коэффициента из-
а
бытка воздуха известно и равно cp = 2,3. Поэтому, учитывая, что горение во фронте пламени начальной (Н) выделенной зоны протекает при оптимальном избытке воздуха а = 0,9 (соответствующем максимальной скорости горения), можно принять следующее распределе-
а
ние состава смеси по зонам: а = 0,9 / 2,0 / 4,0; обеспечивающее cp « 2,3. Принято, что сгорание начинается в обогащенной начальной (Н) зоне (рисунок 1), что не противоречит известным положениям теории сгорания в ДВС, и далее последовательно распространяется на
другие зоны.
Рисунок 1 - Схема разделения на зоны камеры сгорания двигателя AUX VW: Н - начальная; С - средняя; К - конечная; S - свеча зажигания.
Кривые расчетных температур для трех выделенных зон, полученные по уравнению (1), изображены на рисунке 2. Процесс горения смеси продолжается всего ~ 20 мкс и изображается вертикальным участком температурной кривой. Пунктиром изображены расчетные кривые температур горения - максимальные температуры соответствуют окончанию горения в каждой из зон. Начало этих кривых соответствует времени прихода фронта горения в центры выделенных зон Н, С, К.
Конечная температура горения для богатых смесей (а = 0,9) достигается быстрее, чем для бедных, что связано с быстрым сгоранием горючих компонентов топлива. В начале процесса ( т < 10 с) скорость изменения температуры горения непрерывно растет и при достижении максимальной температуры горения стабилизируется на уровне, определяемом ис-
T
ходным коэффициентом избытка воздуха. Максимальное значение max наблюдаются вбли-
зи стехиометрического соотношения горючего и окислителя, при а ~ 0,9, что соответствует известным литературным данным.
Рисунок 2 - Изменение локальных температур рабочего тела в процессе сгорания и расширения для двигателя AUX VW (режим работы с послойным смесеобразованием
Pe = 0,4 МПа, П = 3000 мин-1 )
После окончания горения (т « 1,7 мс) считаем, что продукты сгорания в камере сгорания быстро перемешиваются, при этом усредняются температура и состав смеси. Эти данные могут быть использованы в качестве исходных при решении задачи для однородно перемешанной смеси в процессе ее расширения в цилиндре при известном изменении давления. Результаты этого решения изображены сплошной линией на рисунке 2.
Вследствие недостатка кислорода и несмотря на высокую температуру в начальной зоне, скорость реакций окисления азота незначительна. В зонах же с бедным составом смеси
температура сгорания существенно ниже, что способствует пониженному содержанию NOx в составе продуктов сгорания, даже при избытке окислителя. Совокупное влияние температурного и концентрационного факторов в данном случае проявилось следующим образом.
Содержание
NO„
в начальной зоне составляет
1100 млн-
, что в 1,64 раза превышает кон-
центрацию
NO
в средней зоне и в 2,9 раза для конечной зоны (рисунок 3).
ЫОх,
900
6 00
300
Зона H
г
Зона С
1 Зона К
\
/—
\ 1 г / Начал о сгорани! в зонах
1 [L
I 1
342 378 414 ф, град, пкв
Рисунок 3 - Текущая концентрация оксидов азота в локальных зонах КС двигателя
AUX VW (режим работы с послойным смесеобразованием
Pe = 0,4 МПа, n = 3000 мин-1 )
Раздел 1. Наземные транспортные средства, энергетические установки и двигатели. На рисунке 4 представлены результаты расчета совокупного, осредненного по объему
N0
КС, образования х . Так как объем трех локальных зон принят одинаковым, осредненные температура и доли оксидов азота определялись как средние арифметические.
Температура смеси при расширении резко падает. После перемешивания (т > 1,7 мс)
NОх й
концентрация оксидов азота х достигает максимума и в дальнейшем при температурах ниже 2000...2200К практически не изменяется (рисунок 4). Причиной является уменьшение скорости ведущих реакций окисления азота практически на два порядка, поэтому для процесса окисления азота имеет смысл говорить о температуре «замораживания» реакций обра-
NО_
зования
на последующем участке такта расширения.
Осредненное по объему КС расчетное содержание
в продуктах сгорания после их
перемешивания соответствует
710 млн"
(рисунок 4). Здесь же приведены данные по экспе-
й N0x ДВС
риментально замеренной концентрации х в выпускных газах на выходе из ДВС, откуда видно, что расхождение результатов расчета с экспериментом не превышает 8%, подтверждая тем самым удовлетворительную адекватность расчетной методики.
Эксперимент
NOx
450
300
150
1
350 360 370 380 390 <р, ПКВ
Рисунок 4 - Текущая концентрация NOX в продуктах сгорания двигателя AUX VW
P
0,4 МПа, П = 3000 мин-1 )
(режим с послойным смесеобразованием
Заметим, что полученные данные для расслоенной смеси не соответствуют традиционному решению задачи при горении однородно перемешанной горючей смеси. В варианте системы горения однородной смеси концентрация оксидов азота оказалась в 2,2 раза выше, чем в варианте горения с послойным распределением заряда.
Выводы
Поскольку состав смеси в КС двигателя с послойным зарядом характеризуется существенной неоднородностью, предложено разделить весь объем КС на зоны с различным распределением состава смеси. На этой основе разработаны алгоритм модели и методика расчета текущей температуры рабочего тела и динамики окисления азота в КС двигателя с НВБ.
Предложенная методика расчета позволяет определить локальные образования оксидов азота в зависимости от таких реальных факторов, как распределение топлива по объему КС, локальных скоростей сгорания и теплообмена в локальных зонах, а также локальной температуры, существенно отличавшейся от «индикаторной» температуры рабочего тела.
Для двигателя типа AUX фирмы Volkswagen для режима с послойной организацией за-
а
ряда со средним по объему КС значением коэффициента избытка воздуха cp = 2,3 проведены тестовые расчеты процесса сгорания в отдельных зонах. Доля оксида азота для начальной
зоны сгорания (а =0,9) достигает 1100 млн-1, что в 1,64 раза превышает концентрацию NOx в средней зоне ( а =2) и в 2,9 раза для конечной зоны ( а =4). При этом совокупный по всему объему КС выход оксидов азота снижается в 2,2 раза по сравнению с традиционным вариантом сгорания однородной гомогенной смеси. Установленное расхождение результатов расчета с экспериментом не превышает 8%, подтверждая тем самым удовлетворительную адекватность расчетной методики.
Литература
1. Чесноков С.А. Mоделирование высокотемпературных реакций горения. - Тула: Изд-во Тульского гос. ун-та, 2002. - 163 с.
2. Teodorczyk A., Rychter T.J. Mathematical model of nitric oxide formation in an SI piston engine // J. Techn. Phys. - 28. - 1987. - Р. 47-65.
3. Eichlseder Н., Baumann E., Mtiller P., Neugebauer S. Chancen und Risiken von Ottomotoren mit Direkteinspritzung // MTZ. N 3. 2000. S. 144-152.
4. Кавтарадзе Р.З. Теория поршневых двигателей. Учебник для вузов.- M.: Изд-во MГТУ им. Н.Э. Баумана, 2008.-720 с.
5. Кавтарадзе Р.З. Локальный теплообмен в поршневых двигателях. 2-е изд. испр. и доп. M:. Изд-во MF^ им. Н.Э. Баумана, 2007.-592 с.
6. Хзмалян ДМ., Коган Я.А. Теория горения и топочные устройства. - M.: Энергия, 1976. -488с.
7. Зельдович Я.Б., Садовников П.Я., Франк-Каменецкий Д.А. Окисление азота при горении. - M.: Изд-во Академии наук СССР, 1947. - 147 с.
8. Baulch D.L. et all. High Temperature Reaction Rate Data / D.L. Baulch, D.D. Drysdale, D.D. Horne, A.C. Lloyd // Report University of Leeds.- 1969.- №4.- 58 р.
Коэффициент усиления вакуумногоусилителя тормозов
к.т.н. доц. Шуклинов С.Н.
ХНАДУ
Аннотация. Рассмотрено влияние на коэффициент усиления вакуумного усилителя вида его статической характеристики. Получены зависимости для определения параметров статической характеристики вакуумного усилителя тормозов, коэффициента усиления с учетом нелинейностей статической характеристики.
Ключевые слова: вакуумный усилитель тормозов, коэффициент усиления, параметры статической характеристики
Введение
Коэффициент усиления вакуумного усилителя тормозов является важным критерием для оценки его энергопреобразующих свойств. По существу значение коэффициента усиления определяется как отношение выходного сигнала к входному, т. е. это отношение усилия на штоке усилителя к управляющему усилию, приложенному к поршню управления (толкателю) [4, 5]. Однако в работе [2} в качестве коэффициента усиления используется соотношение давления в тормозном приводе к управляющему усилию на педали тормоза. По сути это отношение является не коэффициентом усиления вакуумного усилителя, а коэффициентом передачи педального узла с главным тормозным цилиндром и вакуумным усилителем. В работах [3, 6] для определения коэффициента усиления вакуумного усилителя тормозов используется отношение давления в тормозном приводе при работающем усилителе к давлению в приводе без подвода энергии к усилителю. В указанных работах очевиден аффинный перенос определения коэффициента усиления гидровакуумного усилителя тормозов для определения данного коэффициента вакуумного усилителя.
