CC BY
134
ПЕРСПЕКТИВНЫЕ ДВИГАТЕЛИ ВНУТРЕННЕГО СГОРАНИЯ
УДК 621.43.052
КРИТЕРИЙ ДЕТОНАЦИИ В ГАЗОВОМ ДВИГАТЕЛЕ С ВЫСОКОЭНЕРГЕТИЧЕСКОЙ СИСТЕМОЙ ЗАЖИГАНИЯ
Ф.И. Абрамчук, профессор, д.т.н., А.Н. Кабанов, доцент, к.т.н.,
А.П. Кузьменко, ассистент, М.С. Липинский, аспирант, ХНАДУ
Аннотация. Предложен подход к выбору предельной степени сжатия в цилиндре газового двигателя с высокоэнергетической системой зажигания. Выполнено сравнение результатов расчета с экспериментальными данными.
Ключевые слова: методика расчета, процесс сгорания, детонация, газовый двигатель, степень сжатия.
КРИТЕРІЙ ДЕТОНАЦІЇ В ГАЗОВОМУ ДВИГУНІ З ВИСОКОЕНЕРГЕТИЧНОЮ СИСТЕМОЮ ЗАПАЛЮВАННЯ
Ф.І. Абрамчук, професор, д.т.н, О.М. Кабанов, доцент, к.т.н.,
А.П. Кузьменко, асистент, М.С. Липинський, аспірант, ХНАДУ
Анотація. Запропоновано підхід до вибору граничного ступеня стиску в циліндрі газового двигуна з високоенергетичною системою запалювання. Виконано порівняння результатів розрахунку з експериментальними даними.
Ключові слова: методика розрахунку, процес згоряння, детонація, газовий двигун, ступінь стиску.
KNOCK CRITERION IN GAS ENGINE WITH HIGH-ENERGY
IGNITION SYSTEM
F. Abramchuk, Professor, Doctor of Technical Science,
O. Kabanov, Assosiate Professor, Candidate of Technical Science,
A. Kuzmenko, assistant, M. Lipinskiy, postgraduate, KhNAHU
Abstract. The approach to choosing limiting compression value in gas engine cylinder with high-energy ignition system has been proposed. The comparison of calculated results with experimental data has been performed.
Key words: calculation method, combustion process, knock, gas engine, compression ratio.
Введение
Известно, что повышение степени сжатия (е) в двигателе с искровым зажиганием приводит к улучшению мощностных и экономических характеристик. Однако при этом сдерживающим фактором является возникновение детонации при достижении определённого порогового значения е.
Определить этот порог можно двумя путями - экспериментальным и расчетным. Первый путь более надежный, однако требует очень больших затрат времени, человеческих и материальных ресурсов в связи с изготовлением большого количества комплектов поршней (в некоторых случаях - и головок цилиндров) с разной геометрией поверхности камеры сгорания и необходимостью много-
кратной разборки-сборки двигателя с целью замены деталей.
двигателе с высокоэнергетической системой зажигания.
Второй путь был бы намного менее затратным. Оптимальным подходом было бы использование математической модели, позволяющей предсказать появление детонации при работе двигателя на данном режиме, и проверка результатов расчета при помощи эксперимента при одновременном уточнении расчетных данных. Поэтому выбор и применение на практике адекватной методики расчета детонации в газовом двигателе с высокоэнергетической системой зажигания в настоящее время является актуальной задачей.
Анализ публикаций
Обзор литературы по теме исследования показал, что основная масса расчетных соотношений для степени сжатия представлена эмпирическими уравнениями вида [1]
Описание критерия
Предельное значение степени сжатия в двигателе с искровым зажиганием лимитируется в основном возникновением детонации при превышении определённого значения в.
В настоящее время экспериментально доказано, что детонация возникает вследствие самовоспламенения несгоревшей части смеси в процессе сгорания. Исходя из этого, для расчета момента появления детонации предлагается использовать критерий детонации, предложенный в [2]. Данный критерий выглядит следующим образом
К, =
^ зт
е?
(2)
в = 17,0-0,0420• + 5,6• (а -1) < 16, (1) где - энергия 13 единице объёма, нако-
где в - предельное значение степени сжатия газового двигателя, соответствующее пределу детонации; ^ - температура топливовоздушной смеси на входе в двигатель, °С; а -коэффициент избытка воздуха (а > 1).
Однако такие соотношения для других двигателей, как правило, оказываются не применимыми.
В [2] впервые предложено использовать так называемый критерий детонации на основе разницы энтальпий несгоревшей смеси в цилиндре в начале и в процессе сгорания. Данный критерий может применяться ко всем типам двигателей с принудительным воспламенением смеси, хотя и создавался изначально для расчета детонации в газовых ДВС. Он нашел применение в ряде работ, где авторы отметили хорошую сходимость расчетных данных с экспериментом [3, 4]. Данный критерий предлагается для расчета детонации в газовом двигателе с высокоэнергетической системой зажигания.
Цель и постановка задачи
Целью исследования является выбор и обоснование методики, позволяющей выполнить выбор предельной степени сжатия в газовом
пившаяся в несгоревшей смеси, кДж/м ; ег -энергия в единице объёма, выделившаяся при сгорании, кДж/м3.
Энергия в единице объёма, выделившаяся в процессе предпламенных реакций в несгоревшей смеси, кДж/м3
Ез
Ус
(3)
где Езт - полная энергия, выделившаяся в процессе предпламенных реакций в несгоревшей смеси, кДж; Ус - объём камеры сго-
3
рания, м .
Энергия в единице объёма, выделившаяся при сгорании, кДж/м3
Ее
Ук
(4)
где Е? - полная энергия, выделившаяся при сгорании, кДж; УК - рабочий объём, м3.
Полная энергия несгоревшей смеси в любой момент времени, кДж
Езт =( К'- К )■
т„
(5)
где кс' - удельная энтальпия несгоревшей смеси в момент начала сгорания, кДж/кг; к -удельная энтальпия несгоревшей смеси в любой момент времени, кДж/кг; ти - масса несгоревшей смеси, кг.
Энергия, выделившаяся в процессе сгорания, кДж
ЕЕ Н и • тР ■
(6)
Для выполнения данной задачи расчетным методом использовалась двухзонная модель, в которой принималось допущение, что температура несгоревшей смеси в процессе одинакова по всему объему несгоревшей зоны. Температура несгоревшей смеси, К
Т = Т •
(11)
где Ни - низшая теплота сгорания топлива, кДж/кг; тР - масса топлива, поступившего в цилиндр за цикл, кг.
где г - номер шага расчета; р - давление в цилиндре, МПа; пи - показатель политропы сжатия несгоревшей смеси
Таким образом, уравнение критерия детонации (2) можно представить в виде
пи = Ь
А.
с,„.
(12)
т„
Ка =-
Ни • т-.
где А - характеристическая газовая посто-(7) янная смеси, Дж/(кг-К); Суи - массовая изо-хорная теплоемкость смеси, Дж/(кг-К).
^ Ук 1 ти
С учетом того, что — = в -1, а-------------= 1 - х
Ус
тр
К, = ( к)- ти-(в-1), (8)
Ни • т¥
Для расчета коэффициента теплоотдачи из зоны несгоревшей смеси в стенки камеры сгорания использовалась формула В. Аннан-да [7]
р0,7 с 0,7
а,ер1 = 0,26-Л,т-^“0,3
и°п • £>и
г^агпзт с
или
Ка =%ГК)( -х)-(в-1).
(9)
где кс' - удельная энтальпия несгоревшей смеси в момент начала сгорания, кДж/кг; к -удельная энтальпия несгоревшей смеси в любой момент времени, кДж/кг; Ни - низшая теплота сгорания, кДж/кг, х - доля выгоревшего топлива; е - степень сжатия.
Для расчета энтальпии несгоревшей смеси использовались зависимости из [6]. Данные полиномы имеют следующий вид
к = / (Т ),
(10)
-2,1-10-
• (т4 + Т4)
\ зт 4 }
(13)
где Хзт - теплопроводность смеси в цилиндре, Вт/(м-К); рзт - плотность смеси, кг/м3; Ст - средняя скорость поршня, м/с; ^.агпзт -динамическая вязкость смеси, Па-с; Вс - диаметр цилиндра, м.
Величины, входящие в уравнение (13), рассчитываются с использованием следующих зависимостей [7]
Xт = 0,000361-Т,0:75.
Рзт = 3,49 -10
-3 Рз
Т„
(14)
(15)
где / (Т) -
полином пятой степени.
^агпзш = 0,56-10-6 • Т^2, (16)
То есть для расчетов необходимо знать изменение температуры несгоревшей смеси в процессе сгорания.
где рзт - давление смеси, Па; Тзт - температура смеси, К.
13
Экспериментальная проверка модели
Экспериментальный стенд, при помощи которого осуществлялась экспериментальная проверка модели, создан на кафедре ДВС ХНАДУ на базе газового двигателя 4ГЧ7,5/7,35.
Для нахождения экспериментальных значений температур несгоревшей смеси использовалась методика обработки индикаторных диаграмм, предложенная в [5]. Для снятия и обработки индикаторных диаграмм использовался измерительно-вычислительный комплекс с использованием программного комплекса PowerGraph.
В процессе эксперимента определить появление детонации в двигателе можно двумя способами: восприятием человеком специфических акустических колебаний и появлением волн давления на индикаторной диаграмме. Первый метод является слишком субъективным, поэтому более рациональным является применение второго метода.
Критерием появления детонации является появление волн давления на индикаторной диаграмме в процессе сгорания. При статистической обработке индикаторных диаграмм детонация считается присутствующей на режиме, если из 200 диаграмм более чем на 20 % и больше присутствует ясно различимая детонация [3]. Если в выборке 10-20 % диаграмм с детонацией, такой режим считается пограничным.
На основании полученных экспериментальных данных определяется максимальное значение критерия детонации Катах, соответствующее появлению детонации. Изменение критерия детонации К^тах в зависимости от относительного угла сгорания в процессе сгорания представлено на рис. 1.
Рис. 1. Изменение критерия детонации в процессе сгорания: 1 - в = 14, 0 = 35 град. пкв до ВМТ, п = 3600 мин-1; 2 -в = 14, 0 = 30 град. пкв до ВМТ, п = =3600 мин-1; 3 - в = 12, 0 = 39 град. пкв до ВМТ, п = 4500 мин1
Во время снятия этого режима отчётливо был слышен характерный «стук», доля индикаторных диаграмм с детонацией на данном режиме составила 73 %.
На рис. 2 приведена индикаторная диаграмма режима со следующими параметрами:
а = 1,17, в = 14, 0 = 35 град. пкв до ВМТ, п = 3600 мин-1, Ктах = 2,1. На рисунке отчётливо виден признак детонации - «пила» на индикаторной диаграмме, появляющаяся вследствие многократного столкновения волн давления и разрежения с чувствительным элементом датчика давления.
Рис. 2. Копия экрана процедуры снятия индикаторной диаграммы (а = 1,17, в = 14, 0 = 35 град. пкв до ВМТ, п = 3600 мин1)
На рис. 3 приведена индикаторная диаграмма режима а = 1,02, в = 14, 0 = 30 град. пкв до ВМТ, п = 3600 мин1 (Ктах = 1,48). Данный режим представляет собой пограничный случай. Детонация появляется поздно и проявляется очень слабо, при снятии режима характерный «стук» не был слышен. Доля
индикаторных диаграмм с детонацией на данном режиме составила 31 %.
Рис. 3. Копия экрана процедуры снятия индикаторной диаграммы (а = 1,02, в = 14, 0 = 30 град. пкв до ВМТ, п = 3600 мин-1)
На рис. 4 приведена индикаторная диаграмма режима а = 1,1, в = 12, 0 = 39 град. пкв до ВМТ, п = 4500 мин-1 (Ктх = 1,38). Признаки детонации на индикаторных диаграммах на данном режиме отсутствуют.
На основании полученных экспериментальных данных установлено, что появлению на режиме признаков детонации, указанных выше, соответствует значение критерия детонации К^тах = 1,42-1,57, в зависимости от режима. С учетом погрешности измерений, рекомендуется считать, что на режиме присутствует детонация, если значение критерия детонации на данном режиме равно Ктх = 1,5, что соответствует предложенному в [3].
Рис. 4. Копия экрана процедуры снятия индикаторной диаграммы (а = 1,1, в = 12, 0 = 39 град. пкв до ВМТ, п = 4500 мин1)
Таблица 1 Сравнение результатов расчета максимального значения критерия детонации с экспериментальными данными
в а 0 п Кйшакт К^шах
град. пкв до ВМТ мин-1
12 1,1 23 1600 1,18 1,25
12 1,1 27 2000 1,21 1,29
12 1,1 32 3000 1,22 1,31
12 1,1 35 3500 1,27 1,35
14 1,02 30 3600 1,37 1,48
14 1,12 35 3600 2,37 2,54
14 1,01 37 3600 1,94 2,15
14 1,17 35 3600 1,87 2,1
12 1,1 39 4500 1,30 1,38
12 1,1 42 5000 1,25 1,41
Для сравнения результатов расчёта критерия детонации с экспериментом использовались данные режимов, приведенных в табл. 1. Погрешность определения расчетного максимального значения критерия детонации К^тахт в сравнении с экспериментальным значением Клтях обозначена как А. Как видно из таблицы 1, погрешность расчета критерия детонации Ктах по предлагаемой методике составила А = 8,7-11,4 %.
Наиболее склонный к детонации режим -режим внешней скоростной характеристики двигателя (ВСХД) с частотой вращения, соответствующей максимальному крутящему моменту (в данном случае п = 3600 мин1). Соответственно выбор оптимальной степени сжатия необходимо выполнять путём анализа этого режима.
Регулировочные параметры (а и 0) на режиме выбирались исходя из условий обеспечения максимальной мощности на конкретном режиме. Значение степени сжатия перебиралось в пределах в = 9,8-14 с шагом Ав = 0,1. Расчёты с использованием критерия детонации Ка показали, что максимальное значение степени сжатия, обеспечивающее бездетона-ционное сгорание на режиме, наиболее склонном к детонации, является в = 13,5.
Выводы
1. Для газового двигателя с искровым зажиганием предложен подход к расчету максимального значения степени сжатия, обеспечивающего бездетонационное сгорание топлива.
2. Сравнение результатов расчета максимального значения критерия детонации Kdmax с экспериментальными данными показало, что погрешность расчетов по приведенной модели составляет А = 8,7.. .11,4 %.
3. Максимальное значение степени сжатия для газового двигателя 4ГЧ7,5/7,35 с искровым зажиганием, обеспечивающее бездето-национное сгорание на всех режимах, рассчитанное с использованием предложенного подхода, составило в = 13,5.
Литература
1. Быков Г.А. Детонационные ограничения
при использовании альтернативных топлив в двигателях с искровым зажиганием / Г.А. Быков // Экотехнологии и ресурсосбережение. - 1995. - № 3. - С. 3-9.
2. Gao J. Knock Modelling in S.I. Engines: PhD
thesis: Mechanical Engineering / J. Gao. -Calgary, Canada, 1993. - 230 p.
3. Attar A. A. Optimization and Knock Model-
ling of a Gas Fueled Spark Ignition Engine: PhD thesis: Mechanical Engineering /
A.A. Attar. - Calgary, Canada, 1997. - 248 p.
4. Bade Shrestha S.O. A Predictive Model for
Gas Fueled Spark Ignition Engine Applications / S.O. Bade Shrestha, G.A. Karim. -Calgary, Canada, University of Calgary, 1999. - 18 p. (Preprint / SAE, University of Calgary, № 1999-01-3482).
5. Иноземцев Н.В. Процессы сгорания в дви-
гателях / Н.В. Иноземцев, В.К. Кошкин. -М. : Машгиз, 1949. - 342 с.
6. Гурвич Л.В. Термодинамические свойства
индивидуальных веществ : справочное
издание: в 4 т. / Л.В. Гурвич, И.В. Вейц,
В.А. Медведев и др. - 3-е изд., перераб. и расшир. - М. : Наука, 1978. - Т. 1 : Элементы O, H (D, T), F, Cl и их соединения, 1978. - 496 с.
7. Кавтарадзе Р.З. Локальный теплообмен в
поршневых двигателях : учеб. пособие для вузов / Р.З. Кавтарадзе. - М. : Изд-во МГТУ им. Баумана, 2001. - 502 с.
Рецензент: В.П. Волков, профессор, д.т.н., ХНАДУ.
Статья поступила в редакцию 27 мая 2011 г.
