Научная статья на тему 'Влияние качества процесса топливоподачи на экономические и экологические показатели транспортного дизеля'

Влияние качества процесса топливоподачи на экономические и экологические показатели транспортного дизеля Текст научной статьи по специальности «Механика и машиностроение»

CC BY
519
66
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по механике и машиностроению, автор научной работы — В. А. Марков, С. Н. Девянин, В. И. Мальчук

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по механике и машиностроению , автор научной работы — В. А. Марков, С. Н. Девянин, В. И. Мальчук

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Влияние качества процесса топливоподачи на экономические и экологические показатели транспортного дизеля»

УДК 621.436

В. А. Марков, С. Н. Д е в я н и н, В. И. Мальчук

ВЛИЯНИЕ КАЧЕСТВА ПРОЦЕССА ТОПЛИВОПОДАЧИ НА ЭКОНОМИЧЕСКИЕ И ЭКОЛОГИЧЕСКИЕ ПОКАЗАТЕЛИ ТРАНСПОРТНОГО ДИЗЕЛЯ

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

Система топливоподачи относится к наиболее важным системам дизельного двигателя. Процессы впрыскивания топлива, его распыли-вания и смесеобразования предопределяют эффективность рабочего процесса дизеля, его мощно стные, экономические и экологические показатели [1, 2]. Причем, в связи с ухудшающейся экологической обстановкой в крупных городах показатели токсичности отработавших газов (ОГ) дизелей становятся приоритетными.

собенностями дизельных двигателей транспортного назначения являются широкий диапазон скоростных и нагрузочных режимов (в условиях реальной эксплуатации транспортных средств) и сложная зависимость между показателями топливной экономичности и токсичности ОГ. В частности, при повышении эффективности процесса сгорания удельный эффективный расход топлива де и выбросы с ОГ продуктов неполного сгорания (монооксида углерода С О, несгорев-ших углеводородов СНх и сажи С) снижаются, а выбросы одного из важнейших токсичных компонентов ОГ — оксидов азота N0^ — напротив, увеличиваются. Поэтому необходима оптимизация процессов топливоподачи и смесеобразования в целях достижения компромисса меду указанными показателями работы дизеля [2, 3].

а современном этапе развития двигателестроения в качестве непременного условия получения высоких экономических и экологических показателей дизелей с объемным смесеобразованием рассматривается интенсификация процесса топливоподачи. Максимальные давления впрыскивания топлива уже перешли рубеж 100 МПа и планируется выпуск топливной аппаратур с давлением впрыскивания порядка 200 МПа [4, 5]. Основной целью повышения давления впрыскивания является улучшение мелкости распыливания. Однако при достигнутом уровне давлений впрыскивания средний диаметр капель слабо зависит от интенсивности процесса топливоподачи, поэтому эффективность

дальнейшего повышения давления впрыскивания невелика. Кроме того, рост давлений впрыскивания ограничивается допустимым уровнем динамики сгорания топлива (величиной максимального давления сгорания pz и интенсивностью нарастания давления при сгорании dp/dp), возможным увеличением эмиссии оксидов азота NOx, а таюке возмож-ностью попадания топлива на относительно холодные стенки камеры сгорания. Следует отметить, что длина струи топлива увеличивается при повышении давления рфmax до 80МПа, а при дальнейшем его повышении остается постоянной [6]. При этом существенно ухудшается гидравлический КПД процесса топливоподачи и увеличивается мощ-ность, необходимая для привода насоса [7]. Все эти факторы вынуждают двигателестроителей искать пути улучшения качества процесса топливоподачи, в том числе и без существенного повышения давления впрыскивания.

Задача совершенствования процессов топливоподачи и смесеобразования наиболее эффективно решается путем отработки конструкции проточной части распылителя форсунки. Заметное влияние на топливоподачу оказывает расположение распыливающих каналов. В большинстве отечественных и зарубежных распылителей распыли-вающие отверстия выходят в центральный канал под конусом иглы распылителя [8-10]. Таковыми являются, в частности, распылители типа 145.1112110 производства Ногинского завода топливной аппаратуры (НЗТА). В этих распылителях устья распыливающих отверстий 1 (рис. 1, а), располагаются в колодце 2 ниже запорного конуса 3 иглы 4 форсунки. Такая конструкция более технологична, однако ее недо-

Рис. 1. Конструктивные схемы распылителей типа 145.1112110 производства НЗТА (а) и распылителей типа DOP 119 S 534 фирмы Motorpal (Чехия) (б):

1 — распыливающее отверстие; 2 — колодец; 3 — седло; 4 — игла

статком является подача топлива, остающегося после закрытия иглы в колодце, в камеру сгорания дизеля под низким давлением, что способствует повышению выбросов сажи и других продуктов неполного сгорания с ОГ. Этот недостаток устраняется при уменьшении объема колодца или при выходе распыливающих отверстий на поверхность запирающего конуса распылителя [8-10]. По такой схеме выполнены, в частности, распылители типа DOP 119 S 534 фирмы Motorpal (Чехия), в которых устья распыливающих отверстий 1 (рис. 1, б) размещены на конической поверхности 3 выше колодца 2, выполненного с уменьшенным объемом. В этом случае при посадке иглы 4 на конус 3 в заключительной фазе топливоподачи объем топлива, сосредоточенный в полости под иглой (в колодце 3), оказывается изолирован от распыливающих отверстий 1.

Исследования распылителей с выходом распыливающих отверстий в колодец (распылители, выполненные по варианту I) и на поверхность запирающего конуса иглы (распылители, выполненные по варианту II), проведенные фирмами R. Bosch (Германия) и AVL (Австрия), а в отечественном двигателестроении — в ГТУ "МАДИ" [9,10], позволили отметить больший объем ("пучистость") струи топлива у распылителей, изготовленных по варианту II. Основной причиной этого является относительно низкая эффективная площадь проходного сечения распылителей II (распылителей Motorpal, кривая 2 на рис. 2) по сравнению с распылителями I (распылителей НЗТА, кривая 1, см. рис. 2), особенно при малых подъемах иглы форсунки. Следует отметить увеличение максимальной скорости истечения топлива в сопловых каналах

распылителей, выполненных по варианту II. Имеет место также более быстрое открытие иглы у этого варианта форсунки из-за более низкой пропускной способности распылителя, которая приводит к повышению давления впрыскивания, действующего на поверхность иглы в направлении ее откр -тия.

Указанное улучшение качества процессов топливопо-дачи и смесеобразования при использовании распылителей (вариант II) позволило повысить эффективность сгорания

ibfp,

О ОД 0,2 0,3 h„, мм

Рис. 2. Зависимости эквивалентного проходного сечения ^р/р от хода Нп иглы форсунки распылителей:

1 — конструкция НЗТА; 2 — фирма Мйогра!

и уменьшить эмиссию продуктов неполного сгорания топлива. В частности, при переходе от работы с распылителями, выполненными по варианту I, к работе с распылителями, выполненными по варианту II, отмечено снижение содержания в ОГ несгоревших углеводородов в 34 раза. Причиной высокой эмиссии СНх при использовании форсунки (вариант I) является то, что топливо, находящееся в ее подыгольной полости, "выпаривается" или выдувается в камеру сгорания в конечной фазе впрыскивания и не успевает полностью сгореть, а затем попадает в ОГ. Меньшие выбросы с ОГ продуктов неполного сгорания топлива в двигателе с форсункой, изготовленной по варианту II, объясняются также тем, что при посадке иглы на седло в камеру сгорания впрыскивается меньше топлива, и при больших давлениях (топливо получает дополнительный импульс и при топливоподаче лучше распыливает-ся). Форсунка (вариант II) отличается меньшим количеством топлива, впрыскиваемого в период задержки воспламенения, что способствует снижению выброса с ОГ оксидов азота NOx. Кроме того, как отмечено в работе [9], при рациональном выборе отношения длины распылива-ющего отверстия /р к его диаметру dp, распылитель, выполненный по варианту II, обеспечивает более мелкое распыливание топлива.

помянутые результат исследований получен либо при стационарном истечении топлива через распыливающие отверстия форсунок (ГТУ "МАДИ"), либо на ограниченном числе режимов работы дизеля (фирм R. Bosch и AVL), что не позволяет в полной мере оценить влияние конструкции распылителя на характеристики процессов впрыскивания и распыливания топлива, а таюке на показатели работа дизеля. Для проведения такой оценки авторами выполнены дополнительные экспериментальные исследования как на безмоторных установках, так и в условиях работы топливной системы на дизеле.

лияние конструкции проточной части распылителя форсунки на параметры и характеристики процессов впрыскивания и распыливания топлива при безмоторных исследованиях оценено с использованием дизельной топливоподающей системы (производство НЗТА), состоящей из топливного насоса высокого давления (ТНВД) типа 4 УТНМ c диаметром плунжеров ¿пл = 9 мм и их полным ходом Л,пл = 10 мм, топливопроводов высокого давления с внутренним диаметром 2,0 мм и длиной Ьт = 690 мм и форсунок ФДМ-22, которые отрегулированы на давление начала впрыскивания рф0 = 22, 0 МПа. Такая система топливоподачи используется в дизеле Д-245 (411/12,5), устанавливаемом на тракторы "Беларусь" (Минский моторный завод), малотоннажные автомобили ЗиЛ-5301 "Бычок" (завод им. Лихачева) и автобусы "ПАЗ" (Павловский автобусный завод). При испытаниях форсунки поочередно оснащали упомянутыми распылителями типа 145.1112110

(производство НЗТА) и распылителями типа БОР 119 Б 534 (фирма МоШгра1). Конструктивные схемы исследуемых распылителей приведены на рис. 1, зависимости эквивалентного проходного сечения ^р/р распылителей от хода Ни иглы форсунки — на рис. 2, а некоторые характеристики этих распылителей приведет: в табл. 1 и 2. При этом распылители конструкции НЗТА выполнены с диаметром иглы йи = 6, 0 мм (по направляющей), а распылители фирмы Мо1;огра1 — с йи = 5, 0 мм.

Таблица 1

Параметры распылителей

№ распыли-теля Изготовитель, маркировка Диаметр распыливающих отверстий d>, мм Число распыли-вающих отверстий гр Максимальный ход иглы Ки, мм Суммарная эффективная площадь распылителя в сборе ^р/р, мм2

1 НЗТА, 145.1112110 0,32 5 0,26 0,278

2 Motorpal, DOP 119 S 534 0,34 5 0,30 0,250

Примечание. Величины ^р/р приведены при максимальном подъеме иглы форсунки; указанные значения Ни и ^р/р — средние для комплекта распылителей.

Таблица 2

Расположение распыливающих отверстий

№ отвер- Распылитель № 1 Распылитель № 2

стия (НЗТА, 145.1112110) (Motorpal, DOP 119 S 534)

Угловое Угол наклона Угловое Угол наклона

расположение отверстия отно- расположение отверстия от-

отверстия сительно оси отверстия носительно оси

относительно распылителя, относительно распылителя,

штифта, град град штифта, град град

1 8 62 8 63,5

2 90 71,5 90 72

3 172 62 172 63

4 237 52 237 52,5

5 303 52 303 53,5

Безмоторные исследования указанных распылителей проведены на топливном стенде венгерской фирмы Мег^. Исследуемый ТНВД, вращающийся от электродвигателя стенда, подавал топливо через нагнетательный топливопровод к исследуемой форсунке, которая впрыскивала топливо в среду без противодавления, причем одна из струй топлива была направлена на датчик, измеряющий силовое воздействие струи топлива.

Стенд оборудован комплексом необходимой измерительной аппаратуры. В процессе испытаний измеряли частоту вращения птн кулачкового вала ТНВД с помощью стендового тахометра аналогового типа. Объемная цикловая подача дц топлива определялась путем замера его расхода через форсунки при помощи мерных колб (мензурок) стенда за фиксированное число впрыскиваний. При испытаниях регистрировали давление топлива рф в нагнетательном топливопроводе у штуцера форсунки при помои пьезометрического датчика давления типа Т6000 и опорный сигнал (угловую метку) положения вала привода ТНВД с помощью индукционного датчика, фиксирующего момент прохождения "флажка" мимо датчика.

ак отмечено ранее, стенд был дооборудован датчиком, регистрирующим силу гидродинамического давления струи распыливаемого топлива на воспринимающую поверхность датчика. Датчик создан на основе пьезодатчика давления ЛХ-601, имеющего высокую чувствительность и хорошую частотную характеристику (диапазон измеряемых частот до 90 кГц). На чувствительную мембрану датчика была наклеена воспринимающая пластинка массой 0,2 г и диаметром 18 мм. При испытаниях на стенде можно изменять расстояние Ьд от одного из распыливающих отверстий форсунки до воспринимающей поверхности датчика. Датчик устанавливался в двух положениях — при Ьд = 5 мм и = 50 мм. Схема установки датчика силового воздействие струи топлива приведена на рис. 3.

Сигналы от датчика углового положения вала привода ТНВД, датчика давления рф топлива у штуцера форсунки и датчика гидродина-

Рис. 3. Схема установки датчика гидродинамического давления струи топлива на воспринимающую поверхность датчика:

Ьд — расстояние от распыливающего отверстия форсунки до датчика

Рис. 4. Осциллограммы процессов впрыскивания и распиливания топлива исследуемым распылителем конструкции НЗТА на режимах с частотой вращения птп = 1200 мин-1 и полной нагрузкой, полученные при установке датчика силового воздействия струи топлива на расстоянии Ьд = 5 мм от распиливающего отверстия форсунки:

Рф — давление топлива на входе в форсунку; — сила гидродинамического давления струи топлива на поверхность датчика

мического давления струи распиливаемого топлива подавались на измерительно-вычислительный комплекс (ИВК-9), имеющий девять входных независимых каналов с синхронной частотой опроса датчиков, проводимого через 6,1 мкс. Информация в ИВК-9 обрабатывалась и сохранялась в ПЭВМ на жестком диске в виде файла. При необходимости эту информацию можно вывести на экран монитора в графическом или цифровом виде, на принтер и т.д.

римеры получаемых при исследованиях осциллограмм давления Рф (кривая 1), силы (кривая 3) и сигнала метки (кривая 2) приведены на рис. 4. Сигнал угловой метки на полученных осциллограммах позволяет связать процесс изменения давления рф и реакции струи с положением вала ТНВД при изменениях частоты вращения, цикловой подачи топлива, расстояния Ьд от датчика давления струи до распы-ливающего отверстия. Перед началом экспериментальных исследований датчики тарировали. Статическая тарировка датчика силы — это определение величины электрического сигнала, регистрируемого ИВК-9, при нагружении воспринимающей поверхности датчика гидродинамического давления струи гирьками различной масс . Тарировку проводили гирьками, создающими силу тяжести в диапазоне от 0 до 8Н. Датчик имел линейную тарировочную зависимость с величиной доверительной вероятности, равной 0,9998.

Исследуемую систему топливоподачи испытывали на двух скоростных режима с частотой вращения кулачкового вала ТНВД птн = 1200 и 800 мин-1. На каждом из этих скоростных режимов устанавливали пять нагрузочных режимов (пять цикловых подач топлива дц), соответствующих нагрузочным режимам 13-ступенчатого цикла, используемого при испытаниях транспортных дизелей на токсичность их ОГ. При этом полной нагрузке двигателя (крутящему моменту двигателя Ме, равному 100 %) соответствовала цикловая подача топлива дц = 75 ... 83 мм3 (разброс значений дц при неизменном положении дозирующей рейки Т обусловлен различными значениями суммарной эффективной площади ^р/р распыливающих отверстий исследуемых распылителей, см. табл.1). Частичным нагрузкам двигателя при Ме =75, 50, 25 и 10% соответствовали цикловые подачи топлива дц=62... 74, 40... 50, 30... 35 и 21.. .29 мм3. Испытания проводили на дизельном топливе с плотностью рт=830 кг/м3 и вязкостью ^т=3,0 мм2/с.

Качественный анализ полученных при экспериментальных исследованиях результатов позволил отметить следующие особенности осциллограмм силы гидродинамического давления струи топлива (см. рис. 4). Характерной особенностью этих осциллограмм является наличие выраженных пульсаций силы свидетельствующих о значительных возмущениях струи распыливаемого топлива. Таюке заметны деформация переднего и заднего фронтов характеристики силы и значительное уменьшение величины этой силы при увеличении Ьд из-за рассеяния энергии струи топлива, взаимодействующей с воздухом.

Геометрия проточной части исследуемых распылителей оказывает заметное влияние на основные параметр процессов впрыскивания и распыливания топлива. Обработка полученных при экспериментальных исследованиях осциллограмм (см. рис. 4) позволила оценить величины двух важных параметров топливоподачи — среднего давления впрыскивания рфср и его продолжительности тв при использовании исследуемых распылителей. Эти параметры заметно влияют на эффективность процессов впрыскивания и распыливания топлива и, в конечном итоге, — на экономические и экологические показатели транспортного дизеля [1,2,4].

Среднее давление впрыскивания рфср характеризует суммарную энергию, реализуемую за период впрыскивания. При повышении рфср изменяется характер протекания процесса смесеобразования. Распад струи топлива начинается непосредственно на выходе из распыливаю-щего отверстия, возрастает объем струи, и достигается более полный охват камеры сгорания струями топлива. Кроме того, при интенсификации впрыскивания увеличиваются скорости протекания топлива по распыливающим отверстиям (приводящие к турбулизации потока)

и скорости распространения струй топлива в камере сгорания. Эти факторы приводят к уменьшению диаметра капель топлива и выравниванию значений коэффициента избытка воздуха а по объему камеры сгорания.

При увеличении давлений впрыскивания сокращается продолжительность впрыскивания твпр (или выраженная в градусах поворота кулачкового вала (п.кул.в.) — ^Впр). В современных быстроходных дизелях транспортного назначения величина ^впр обычно не превышает 10... 15° п.кул.в. ТНВД (20.. .30о п.к.в. — поворота коленчатого вала двигателя). В большинстве случаев, если нет нарушений согласования ^впр и параметров воздушного заряда, уменьшение продолжительности впрыскивания улучает качество процесса смесеобразования, увеличивает эффективность процесса сгорания и снижает выбросы сажи (углерода С) и других продуктов неполного сгорания топлива (монооксида углерода СО и легких углеводородов СНх). Но при этом может увеличиваться эмиссия оксидов азота N0х. Однако при сокращении продолжительности впрыскивания появляется возможность работа с поздними углами опережения впрыскивания топлива, позволяющими снизить эмиссию N0х.

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

арактеристики среднего за период впрыскивания давления топлива рф ср на входе в форсунку, полученные при экспериментальных исследованиях распылителей (конструкции НЗТА и фирмы Мо1;огра1 (см. табл. 1,2), представлены на рис. 5. Анализ характеристик Рфср, полученных на режимах с различными частотой вращения птн кулачкового вала ТНВД и нагрузкой (при подачах топлива дц, соответствующих различным значениям крутящего момента Ме коленчатого вала двигателя), показал, что на всех исследованных режимах наибольшие давления рфср наблюдаются при установке распылителей фирмы Мо1;огра1. В частности, на режиме номинальной мощности при птн=1200 мин-1 и Ме=100% распылитель фирмы Мо1;огра1 обеспечивает среднее за время впрыскивания давление на входе в форсунку рфср = 23, 2 МПа, а распылитель конструкции НЗТА — рфср = 21, 2 МПа. Преимущество

22

20

18

16

14

У К* У

п тн =1200 ми! ! 1 У У ** ' 2

\ \ У ✓ \ * \ ^ ' ^ ^^^ \

> ■и/

✓ > / / ^ / ✓

/ ?тн=800 мин \\ ✓ У / У 1

\ \ * / Г

< У" / //

Рф ср, МПа

■22

20

18

16

20

40

60

80

14

N1 , %

Рис. 5. Зависимость среднего за время впрыскивания давления рфср топлива на входе в форсунку от скоростного и нагрузочного режимов при использовании распылителей:

1 — конструкция НЗТА; 2 — фирма Мо1огра1

распылителей фирмы Мо1огра1 по величине рфср объясняется меньшими значениями суммарной эффективной площади ^р/р, имеющими место при всех подъемах иглы форсунки Ни (см. рис. 2).

Тип используемого распылителя влияет и на продолжительность процесса впрыскивания, о чем свидетельствуют характеристики продолжительности воздействия тв струи топлива на датчик силы (рис. 6). Наибольшая продолжительность впрыскивания отмечена при установке распылителей фирмы Мо1огра1. На номинальном режиме работы при частоте вращения птн = 1200 мин-1 и нагрузке Ме = 100 % использование этих распылителей дает продолжительность впрыскивания тв = 1,13 мс, а распылителя конструкции НЗТА — тв = 0, 93 мс. Большая продолжительность впрыскивания распылителями фирмы Мо1огра1 отмечена и на других исследуемых режимах. Как и в случае со средним давлением рфср, такая зависимость тв от типа распылителя объясняется меньшими значениями эквивалентного проходного сечения распыливающих отверстий ^р/р у распылителей фирмы Мо1огра1.

Геометрия проточной части исследуемых распылителей оказывает заметное влияние на ряд других параметров процессов впрыскивания и

Тв,

MC 1,1

1,0

0,9

0,8

0,7

0,6

2 Ч--

Лтн= 1200 Л МИН"1 ■i» ж»

/> х' f 1

\

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

г*

2

У У é

п Т н= 800 мин л * 1 /

/ \ ^' / \s / / у \

У у ✓ /»" \

У к»

* Г

20

40

60

80 Me , %

Тв, MC

1,3 1,2 1,1 1,0 0,9 0,8 0,7 0,6 0,5

Рис. 6. Зависимость времени тв впрыскивания от скоростного и нагрузочного режимов при использовании распылителей: 1 и 2 — см. рис. 5

распиливания топлива. Одним из таких параметров является силовое воздействие F струи топлива, характеризующее интенсивность процесса впрыскивания топлива. Средняя за время впрыскивания сила ^сср гидродинамического давления струи топлива на воспринимающую поверхность датчика может быть определена по известному аналитическому соотношению [11,12]

^сср = Рт^с^с sin Y, (1)

где рт — плотность используемого топлива; ис — средняя скорость истечения топлива из распыливающего отверстия; — объемный расход топлива через распиливающее отверстие; y — угол между осью струи и воспринимающей поверхностью датчика. Из выражения (1) следует, что величина ^сср характеризует энергию струи топлива, которая, в свою очередь, определяет другие параметры процессов впрыскивания и распыливания.

Как отмечено ранее, осциллограммы силы F гидродинамического давления струи топлива на поверхность датчика (см. рис. 4) имеют ярко выраженный колебательный характер. Поэтому величина си-

лы Fc не в полной мере характеризует возмущенность струи топлива (ее турбулизацию) и склонность ее к распаду в камере сгорания дизеля. Степень возмущенности струи топлива (наличие в ней пульсаций давления) мо но оценить при частотном анализе характеристики силы F;, т.е. оценивая величины амплитуд низко-, средне- и высокочастотных составляющих указанной характеристики. Для этого были проведены расчетные исследования на персональном компьютере в программной среде операционной системы Windows 2000 с использованием встроенного пакета математических функций электронных таблиц Excel.

При расчетном исследовании было выделено 256 точек, по которым анализировали осциллограммы силы F (см. рис. 4), равномерно расположенные во временном диапазоне Д т = 1, 56 мс (интервал между смежными точками составил 8т = 0, 0061 мс = 6,1 мкс). На каждом исследуемом режиме работы характеристика силы описываемая этой совокупностью точек, разлагалась в ряд Фурье, позволявший представить исходну негармоническу функци в виде суммы бесконечного множества гармонических функций с различными частотой f (периодом T) и амплитудой A. По результатам проведенных расчетных исследований построены амплитудно-частотные спектры (АЧС), вид которых представлен на рис. 7.

Сравнительный анализ полученных АЧС показал, что при установке датчика силы F на расстоянии Ьд = 5 мм от распыливающего отверстия форсунки диапазон частот f колебаний силы в котором амплитуды имеют значимые величины, составил 0,64... 80 кГц (см. рис. 7, а). Для экспериментов с установкой датчика силы F на расстоянии Ьд = 50 мм от распыливающего отверстия форсунки характерно уменьшение максимальных частот колебаний со значимыми амплитудами до f = 50 кГц (см. рис. 7, б). Такая закономерность изменения характеристик силы F при увеличении Ьд наблюдается на всех режимах работы обоих типов исследуемых распылителей.

Для оценки степени турбулизации струи распыливаемого топлива с использованием полученных АЧС были рассчитаны значения площади Sa под АЧС, характеризующие интегральную для всего диапазона частот f (0,64... 80 кГц) возмущенность струи топлива. Анализ полученных характеристик площади Sa под АЧС (рис. 8) показал, что при удалении фронта струи топлива от распыливающего отверстия форсунки степень турбулизации струи топлива значительно уменьшается. Например, при исследовании распылителя НЗТА на номинальном режиме при птн = 1200 мин-1 и Ме = 100 % площадь под АЧС уменьшается с SA = 3,60 Н-кГц при Ьд = 5 мм до SA = 2, 35 Н-кГц при Ьд = 50 мм (см. рис. 8).

0,3

0,2

1

1

1 1 1

к м .л/У V U)

0,1

10 20 30 40 50 60 70 f, кГц

0,2

1 1 1 1 1 1

1 1 1 к2 1

1 L1 'II

% ц /I I 1 \ 1V- /*»Г|

0,1

50

70 f, кГц

Рис. 7. Амплитудно-частотные спектры колебаний силы гидродинамического давления струи топлива на поверхность датчика на расстоянии Ьд = 5 мм (я) и Ьд = 50 мм (б) на режимах с частотой вращения итл = 1200 мин-1 и нагрузкой Ме = 10 % при использовании распылителей:

1 и 2 — см. рис. 5

Значительно изменяется величина площади Б а и при смене скоростного и нагрузочного режимов работы топливоподающей аппаратуры. При использовании распылителя НЗТА, полной нагрузке Ме = 100 % (полной подаче топлива) и Ьд = 5 мм снижение частоты вращения птн

приводит к уменьшению

кулачкового вала ТНВД от 1200 до 800 мин 1

Sa, H кГц

3,5 3,0 2,5 2,0 1,5 1,0 0,5

3

"тн = 1200 МИ1 - л---- 2

У 1

/ л \ // 4

/ ' t /

У ' / / 3^_

/ / У 2

* Л - ' '

^ ^ V \ 1

Птн1 800 УЛИН"1

20

40

60

80

Ме

Н кГц

2,5

2,0 1,5 1,0

0,5

%

Рис. 8. Зависимость площади За под АЧС колебаний силы Ре гидродинамического давления струи топлива на поверхность датчика на расстоянии Ьд = 5 мм (кривая 3) и Ьд = 50 мм (кривые 1 и 2) от скоростного и нагрузочного режимов при использовании распылителей:

1 и 3 — конструкция НЗТА; 2 — фирма Мо1огра1

площади Ба от 3,60 до 2,65 Н-кГц. Падение нагрузки Ме со 100 до 10 % (при птн = 1200 мин-1) приводит в этом случае к уменьшению Б а от 3,60 до 1,15 Н-кГц (см. рис. 8). Таким образом, нагрузочный режим работа топливоподающей аппаратуры (цикловая подача топлива дц) оказывает большее влияние на возмущенность струи распыливаемого топлива, чем ее скоростной режим (частота вращения птн). Этим и обусловлено значительное ухудшение качества процессов распыливания топлива и смесеобразования на режимах с малой нагрузкой.

Улучшение качества распыливания топлива в широком диапазоне режимов работы топливоподающей аппаратуры может быть достигнуто путем оптимизации геометрии проточной части распылителя форсунки, оказывающей заметное влияние на величину площади Б а под АЧС. При этом следует заметить, что полученный при экспериментальных исследованиях распылителя фирмы Мо1огра1 на режимах п = 800 мин-1 (см. рис. 8) график зависимости Ба = /(Ме) оказался более пологим, чем для распылителя НЗТА. Это объясняется следую -щим. Известно, что с уменьшением нагрузки увеличивается доля топ-

лива, подаваемого в камеру сгорания дизеля при частичных подъемах иглы Ни. Анализ показывает [13], что при частичных подъемах иглы (при прочих равных условиях) располагаемый напор перед распиливающими отверстиями у распылителей типа II значительно выше, чем у распылителя типа I. Это объясняется наличием дополнительных гидравлических потерь при внезапном расширении потока топлива, втекающего в колодец из щели между запирающими конусами, что имеет место у распылителя, выполненного по варианту I.

При Ьд = 50 мм на всех исследуемых режимах наилучшие показатели по величине Sa дает использование распылителя фирмы Motorpal (см. рис. 8). Это свидетельствует о том, что при продвижении фронта струи топлива вглубь камеры сгорания наибольшую возмущенность струи обеспечивает именно этот распылитель, что создает предпосылки к быстрому распаду струи, повышению качества процесса смесеобразования, улучшению экономических и экологических показателей транспортного дизеля.

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

ценка влияния конструкции проточной части распылителя форсунки на экономические и экологические показатели проведена с использованием результатов экспериментальных исследований дизеля Д-245 на моторном стенде АМО "ЗиЛ". Дизель был оснащен топливной системой с ТНВД фирмы Motorpal (Чехия) типа PP4M10U1f c диаметром плунжеров = 10 мм и их полным ходом Л,пл = 10 мм, топливопроводами высокого давления с внутренним диаметром 2,0 мм и длиной Ьт = 540 мм и форсунками ФДМ-22 производства АО "Куро-аппаратура" (г. Вильнюс), которые были отрегулированы на давление начала впрыскивания рфо = 22, 0 МПа.

Программа исследований дизеля Д-245 с распылителями различ-нх типов предусматривала работу двигателя на различных установившихся режимах, с частотой вращения коленчатого вала n от 750 до 2400 мин-1 и нагрузкой — эффективным крутящим моментом дизеля Ме — от 0 до 360 H•м. Дизель исследовали при неизменном положении упора дозирующей рейки ТНВД (упора максимальной подачи топлива) с постоянным штатным для исследуемого дизеля установочным углом опережения впрыскивания топлива, равным в = 13° п.к.в. до верхней мертвой точки (ВМТ).

GT, кг/ч

20

15

10

п = 240I ) мин"1 V / X /

я < 2

* / V/' / / / о / / /

/ / / / / / / / / / / / X

/> / / / / г/ / / / / / / / / / '/ /

/ / X / / / / / п / = 1500 м ин'1 / / / / / / X / / s'' >о

/V XV о' / / / / У

1 X о-' ' уУ^ У У о У у

л = 760 / мин"1

100

200

300 Ме, Н м

Рис. 9. Зависимость часового расхода топлива Gт от скоростного и нагрузочного режимов (частоты вращения п и эффективного крутящего момента Ме) дизеля Д-245 при использовании распылителей:

1 и 2 — см. рис. 5

Оценка показателей двигателя при его работе на режимах с различными частотой вращения и нагрузкой проведена по экспериментальным данным, полученным при работе дизеля на режимах 13-ступенчатого испытательного цикла. Такие данные по часовому расходу топлива Ст, приведенные на рис. 9, свидетельствуют о том, что на всех исследуемых режимах 13-ступенчатого цикла наибольший расход топлива наблюдается при использовании распылителей конструкции НЗТА. Сравнительная оценка интегральной топливной экономичности дизеля с распылителями различных типов проведена по условному среднему на режимах 13-ступенчатого цикла эффективному расходу топлива, определяемому из работы [3] по выражению:

geусл = ^(GiK)/ ),

(2)

i=1

i=1

где Стг — часовой расход топлива на г-м режиме; — мощность дви-

гателя на этом режиме; K — коэффициент, отражающий долю времени каждого режима; k — число режимов. Полученные при расчетах значения де усл для дизеля с исследуемыми распылителями приведены в табл. 3. Из табл. 3 видно, что при замене распылителей конструкции НЗТА на распылители фирмы Motorpal условный средний на режимах 13-ступенчатого испытательного цикла расход топлива де усл снижается на 5,9%. Такое улучшение топливной экономичности при использовании распылителей фирмы Motorpal обусловлено улучшением качества процесса смесеобразования по сравнени с применением распылителей конструкции НЗТА.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Таблица 3

Показатели дизеля Д-245 с распылителями различных типов

Распылитель Условный средний расход топлива на режимах 13-ступенчатого цикла, де усл, г/(кВт-ч) Удельные массовые выбросы токсичных компонентов при работе по 13-ступенчатому циклу, г/(кВт-ч)

eNOx ecO eCHx

Конструкции НЗТА 268,3 6,439 4,464 2,958

Фирмы Motorpal 252,4 6,741 3,067 0,945

Примечание. Требование EURO-2 — eNOx = 7,0; eCO = 4,0; eCHx =1,1 г/кВт^ч соответственно

От типа применяемых распылителей в значительной степени зависит и токсичность ОГ двигателя, т.е. выбросы нормируемых токсичных компонентов ОГ — оксидов азота NOx, монооксида углерода СО, легких углеводородов СНх, а также эмиссия твердых частиц или выбросы сажи (углерода С).

Приведенные на рис. 10, а графики концентрации в ОГ оксидов азота NOx от типа распылителя свидетельствуют о достаточно сложном характере этой зависимости. В частности, на режиме номинальной мощности при n = 2400 мин-1 наименьшее содержание оксидов азота в ОГ (CNOx = 0, 0545 %) имеет место при установке на исследуемый дизель распылителей конструкции НЗТА (кривая 1, см. рис. 10, а). Использование распылителей фирмы Motorpal приводит к увеличению CNOx на этом режиме соответственно до 0,064% (кривая 2, см. рис. 10,а). На режимах с n = 2400 мин-1 и малой нагрузкой концентрации NOx в ОГ примерно одинаковы. На режиме максимального крутящего момента при n = 1500 мин-1 наименьшая концентрация NOx в ОГ отмечена при использовании распылителей конструкции НЗТА (CNOx = 0, 06 %), а наибольшая — при установке распылителей

250

я

2 / V / '' У

/ У ■ V'

п = 2400 V ЛИН"1 к и / f * "ч

У / , / / ,у\ ' / \ / / \ 4

/ 1 J 9/ л \ = 1500 X МИН-1

/ / / / л f

\ X Г / / / / f

2ч У ✓ У

/7 = 7 60 мин г1

*10l %

750

500

250

100

200

300 Ме , Н м

Рис. 10 (начало). Зависимость объемной концентрации в ОГ оксидов азота (а), монооксида углерода Ссо (б) и легких углеводородов СсНх (в) от скоростного и нагрузочного режимов (частоты вращения п и эффективного крутящего момента Ме) дизеля Д-245 при использовании распылителей:

1 и 2 — см. рис. 5

фирмы Motorpal (Скох = 0, 07 %). На этом же скоростном режиме, но при минимальной нагрузке с Ме = 36 Н•м установка распылителей конструкции НЗТА, напротив, приводит к максимальной концентрации N0^ (Скох = 0, 025 %) по сравнению с распылителями фирмы Motorpal (^оГ = 0,015 %).

Расчет интегральных удельных массовых выбросов оксидов азота на режимах 13-ступенчатого испытательного цикла проведен по формуле из работы [3]

13

13

eNOx = Ki) Ki),

i=1 i=1

(3)

° со

»104, %

500

О \ * 4 V

\ч N п- = 2400 м \н ИН"1 1

■ ^--о ~ . 0 4 2

N

V X \ п = 1! 500 N ЛИН"1 / °

»4. •V л / / / ✓ / / / г

Л /7=76 0 мин""1

^со

<104, %

1000

500

100

200

300 Ме, Н м

Рис. 10 (продолжение)

где Е]0хг — массовый часовой выброс оксидов азота на г-м режиме. Полученные при расчете данные свидетельствуют о том, что оба типа исследуемых распылителей обеспечивают примерно одинаковые значения е]0х, меньше предельно допустимого выброса е]0хпред = = 7, 0 г/(йЗт-ч), установленного европейским стандартом ЕиК.О-2 на выброс этого токсичного компонента ОГ (см. табл. 3).

Тип применяемого распылителя оказывает значительное влияние и на эмисси с Г продуктов неполного сгорания топлива — монооксида углерода СО и углеводородов СНх. Влияние типа распылителя на выбросы СО неоднозначно. Из приведенных на рис. 10, б данных по содержанию в ОГ монооксида углерода СО следует, что на режимах максимальной (номинальной) мощности (при п = 2400 мин-1) и максимального крутящего момента (при п = 1500 мин-1) минимальные концентрации этого токсичного компонента в ОГ (соответственно ССО = 0, 02 и 0,054%) обеспечивают распылители фирмы Мо1югра1 (кривые 2, см. рис. 10, б). Максимальная концентрация СО в ОГ (ССО = 0, 09 %) отмечена на режиме с п = 1500 мин-1 и Ме=85 Н-м при установке распылителей НЗТА (кривые 1, см. рис. 10, б).

На режимах 13-ступенчатого испытательного цикла интегральные удельные массовые выбросы монооксида углерода еСО рассчитаны по следующему выражению из работы [3]:

13 13

еСо = К)/ (4)

г=1 г=1

где ЕСОг — массовый часовой выброс монооксида углерода на г-м режиме. Сравнение расчетных данных по удельным массовым выбросам СО показывает, что лишь использование распылителей фирмы Мо1;огра1 позволяет выполнить требования стандарта ЕЦКО-2 по выбросам этого компонента ОГ. При установке распылителей фирмы Мо1»гра1 удельный массовый выброс монооксида углерода на режимах 13-ступенчатого цикла составил еСО = 3,067г/(кВт-ч), что существенно ниже предельно допустимого выброса (еСОпред = 4, 0 г/(кВт-ч)) по стандарту ЕЦКО-2.

Наибольший эффект от использования распылителей фирмы Мо1;огра1 получен по выбросам с ОГ легких углеводородов. Как следует из анализа данных рис. 10, в, на котором приведены зависимости концентрации в Г легких углеводородов С х от типа используемых распылителей, замена распылителей конструкции НЗТА на распыли-

тели фирмы Motorpal приводит к сокращению содержания СНх в ОГ в среднем в 2-3 раза на всех исследуемых режимах. Максимальное снижение концентрации СНх в ОГ отмечено на режиме с n = 1500 мин-1 и Ме = 85 Н-м, на котором переход от распылителей конструкции НЗТА к распылителям фирмы Motorpal приводит к снижению ССНх с 0,152 до 0,031 %, т.е. в 4,9 раза. Такое уменьшение содержания СНх в ОГ объясняется улучшением качества (мелкости) распыливания топлива, увеличением объема (ширины) топливных струй, сокращением их длины и уменьшением количества топлива, попадающего на стенки камеры сгорания, при использовании распылителей с устьями распыливающих отверстий, расположенных на конусе седла иглы распылителя.

роведенные расчеты интегральных удельных массовых вбросов легких углеводородов еСНх по выражению [3]

13 13

еснх = ^(Еснх, K)/^(Ne i Ki), (5)

i=1 i=1

где ЕСНхг — массовый часовой выброс углеводородов на i-м режиме, также подтвердили значительное уменьшение еСНх при переходе от распылителей конструкции НЗТА к распылителям фирмы Motorpal. При этом требования стандарта EURO-2 по выбрoсам СНх (ССНхпред = = 1,1 г/(кВт-ч)) обеспечивались лишь при использовании распылителей фирмы Motorpal (ССНх = 0, 945 г/(кВт-ч), см. табл. 3).

Таким образом, при анализе интегральных показателей токсичности ОГ, соответствующих работе двигателя на режимах 13-ступенчатого испытательного цикла, отмечено, что исследуемый дизель имеет показатели, удовлетворяющие требованиям стандарта EURO-2 к токсичности ОГ лишь при оснащении его распылителями фирмы Motorpal (см. табл. 3). Установка распылителей конструкции НЗТА не обеспечивает получение требуемых показателей токсичности Г по выбросам монооксида углерода СО и легких углеводородов СНx.

Результаты исследований, проведенных на моторной установке, подтвердили преимущества распылителей, входные кромки распыливающих отверстий которых расположены на запирающем конусе седла игл форсунки, а не в объеме подыгольной полости форсунки (в колодце распылителя). Следует отметить, что моторные испытания двигателя, оснащенного распылителями различных типов, с определением показателей топливной экономичности и токсичности ОГ, являются достаточно трудоемкими, дорогостоящими и потому не всегда возможными, особенно при исследовании большого числа различных комплектаций распылителей. Поэтому было бы очень желательно прогнозировать экономические и экологические показатели транспортного дизеля на

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

Такое прогнозирование экономических и экологических показателей дизеля может быть проведено путем сопоставления названных показателей с величиной площади БА под АЧС силы гидродинамического давления струи топлива на поверхность датчика этой силы. Следует учитывать, что в быстроходных дизелях с объемным смесеобразованием начальные очаги воспламенения распыливаемого топлива возникают в камере сгорания на расстоянии 12... 15 мм и более от распиливающего отверстия [14]. Поэтому при сопоставлении экономических и экологических показателей дизеля с величиной площади БА под АЧС силы ^си целесообразно использовать данные, полученные при установке датчика силы ^сна расстоянии = 50 мм от распыливающего отверстия форсунки.

Приведенные на рис. 11 данные, построенные с использованием результатов проведенных безмоторных (см. рис. 7, 8) и моторных (см. рис.9, 10) испытаний топливной системы дизеля Д-245, свидетельству т о наличии зависимости показателей топливной экономичности и токсичности ОГ от площади Б а под АЧС силы Так, с увеличением значения Б а (с ростом возмущенности струи распыливаемого топлива, вызванным увеличением цикловой подачи топлива) удельный эффективный расход топлива де, концентрации в ОГ монооксида углерода ССО и легких углеводородов ССНх монотонно уменьшаются, а содержание в ОГ оксидов азота С]Ох, напротив, увеличивается (см. рис. 11). Такая зависимость показателей дизеля от Б а отмечена как при частоте вращения кулачкового вала ТНВД птн = 1200 мин-1 (рис. 11, а), так и при птн = 800 мин-1 (рис. 11, б). Повышение концентрации ССО в ОГ при больших значениях Б а на режиме максимального крутящего момента при птн = 800 мин-1 и Ме = 100 % (две точки при Ба = 1, 75 и 1,90Н-кГц на рис. 11, б). обусловлено увеличением температур сгорания топлива на режиме максимального крутящего момента и диссоциацией образовавшегося диоксида углерода СО2 с образованием монооксида СО и кислорода О2 в высокотемпературных зонах топливного факела [2].

Проведенный по методике работ [15, 16] корреляционный анализ приведенных на рис. 11 зависимостей показал, что при птн = 1200 мин-1 (см. рис. 11, а) выборочный коэффициент корреляции для зависимости де = f (Ба) оказался равным г = 0, 776, для зависимости С]Ох = f (БА) — г = 0, 915, для зависимости ССО = f (БА) — г = 0, 888, для зависимости СсНх = f (Ба) — г = 0,863. Для режимов с птн = 800 мин-1

де,г I (кВт-ч) ge, г/(кВтч)

Рис. 11 (начало). Зависимость удельного эффективного расхода топлива де, объемной концентрации в ОГ оксидов азота CNox , монооксида углерода Ссо и легких углеводородов Сснх от площади Sa под амплитудно-частотным спектром колебаний силы Fe гидродинамического давления струи топлива на поверхность датчика, установленного на расстоянии Ьл = 50 мм от распыливающего отверстия на скоростных режимах птн = 1200 мин-1 (а) и птн = 800 мин-1 (б), при использовании распылителей конструкций НЗТА (1) и фирмы Motorpal (2): х — 1; о — 2

(см. рис. 11, б) коэффициенты корреляции r исследуемых зависимостей параметров дизеля от SA равны соответственно 0,850; 0,877; 0,825 и 0,908. Сравнительно большие значения коэффициентов корреляции, полученные при корреляционном анализе, подтверждают наличие тесной связи между экономическими и экологическими показателями работы дизеля и качеством процесса топливоподачи, характеризуемым величиной Sa.

Рис. 11 (окончание)

Характерно, что при равных значениях Ба для зависимостей С]Ох = f (Ба) и ССО = f (Ба) комплектация II имеет существенные преимущества. В то же время по концентрации ССО это не наблюдается, т.е. можно предположить, что у данного дизеля имеются значительные резервы, позволяющие проводить дальнейшее совершенствование процесса топливоподачи.

Предложенная методика оценки качества процессов впрыскивания и распыливания топлива в транспортном дизеле по результатам безмоторных испытаний позволяет определить конструктивные соотношения топливоподающей аппаратуры без проведения моторных испытаний. Это создает предпосылки для совершенствования качества рабочего процесса дизеля в широком диапазоне его рабочих режимов. Результаты исследований подтвердили значительную зависимость расхода топлива дизеля и токсичности его ОГ от типа используемых распылителей.

СПИСОК ЛИТЕРАТУРЫ

1. Топливные системы и экономичность дизелей / И.В. А стахов, Л Л. Голубков, В.И. Трусов и др. - М.: Машиностроение, 1990. - 288 с.

2. Марков В. А., Баширов Р.М., Га бито в И.И. Токсичность отработавших газов дизелей. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2002. - 376 с.

3. Марков В. А., Кислов В. Г., Хватов В. А. Характеристикитопливоподачи транспортных дизелей. - М.: Изд-во МГТУ им. Н.Э. Баумана, 1997. - 160 с.

4. Исследование возможности получения повышенных давлений впрыскивания топливоподающей аппаратурой разделенного типа в автотракторных дизелях / Л.В. Грехов, В.А. Марков, В.А. Павлов и др. // Вестник МГТУ. Серия "Машиностроение". - 1997. - № 1. - С. 92-103.

5. Грехов Л. В. Топливная аппаратура с электронным управлением дизелей и двигателей с непосредственным впрыском бензина. - М.: Изд-во "Легион-Автодата", 2001. - 176 с.

6. Kamimoto T., Yokota H., Kobayashi H. Effect of High Pressure Injection Soot Formation in a Rapid Compression Machine to Simulate Diesel Flames // SAE Technical Paper Series. - 1987. - № 871610. - P. 9.

7. Некоторые результаты исследований энергетического баланса системы топливоподачи быстроходного дизеля / В.И. Мальчук, Л.П. Музыка, В.А. Павлов и др. // Автотракторные двигатели внутреннего сгорания: Сб. науч. трудов МАДИ.

- М.: Изд-во МАДИ, 1978. - С. 60-66.

8. Ф а й н л е й б Б. Н. Топливная аппаратура автотракторных дизелей: Справочник.

- Л.: Машиностроение, 1990. - 352 с.

9. Трусов В. И., Мальчук В. И., Зрячкин М. В. К выбору конструктивных параметров распылителя по заданным характеристикам впрыска и распыливания топлива // Труды МАДИ. - 1979. - Вып. 178. - С. 58-62.

10. Повышение эффективности подачи и распыливания топлива в дизелях /

B.А. Марков, В.И. Мальчук, С.Н. Девянин и др. // Грузовик &. - 2003. - № 6. -

C. 30-32; - № 7. - С. 23-27; - № 8. - С. 50-51.

11. Башта Т. М. Машиностроительная гидравлика. - М.: Машиностроение, 1971.

- 676 с.

12. Г и д р а в л и к а, гидромашины и гидроприводы / Т.М. Башта, С.С. Руднев, Б.Б. Некрасов и др. - М.: Машиностроение, 1982. - 123 с.

13. Луканин В. Н., Мальчук В. И. Метод расчета гидравлических параметров корректирующих распылителей // В сб.: "Двигатели внутреннего сгорания: проблемы, перспективы развития". - Труды МАДИ (ТУ). - 2000. - С. 104-113.

14. Frolov S. M., Skripnik A. A., Kavtaradze R. Z. Modeling of Diesel Spray Ignition // Combustion and Atmospheric Pollution. - М.: Torus Press Ltd., 2003.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

- P. 220-227.

15. Крамер Г. Математические методы статистики. - М.: Мир, 1975. - 648 с.

16. Митропольский А. К. Техника статистических вычислений. - М.: Наука, 1971.-576 с.

Статья поступила в редакцию 24.06.2004

i Надоели баннеры? Вы всегда можете отключить рекламу.