Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

MODEL OF A TURBO-SUPERCHARGING IN A CYCLE OF THE ENGINE WITH VARIABLE PRESSURE OF AIR UPON AN ADMISSION

Gavrilov A.A. 1 Gots A.N. 1
1 Vladimir State University named after Alexander and Nikolay Stoletovs
The analysis of efficiency of teamwork of the piston engine and turbocharger is lead at any kind of pressurization (pulse, isobaric or intermediate) and is shown, that it depends on size of amplitudes of fluctuation of pressure of a fresh charge in the inlet pipeline. Because of cyclicity of processes of an admission pressure of a working body and factor of filling that affects not only on economic, but also ecological parameters of the engine changes. For reduction of influence of wave processes usually recommend to choose volumes of the inlet pipeline not less than three volumes of working process of the engine. For the coordination of account characteristics of the engine and the compressor on a design stage the mathematical model of the turbo-supercharging, taking into account features of processes in the cylinder and systems of gas exchange have been developed at the big amplitudes of fluctuation of the pressure, the admissions caused by non-uniform order (or by their absence) a fresh charge in the cylinders connected to the pipeline. The offered model is the compound block of model of a cycle of the piston engine.
the piston engine
turbocharger
fluctuation of pressure
working body
factor of filling
1. Gavrilov A.A., Efros V.V. Pulse system of pressurization four-cycle малоцилиндровых diesel engines // Tractors and agricultural machines. 1997. no. 10. pp. 16 18; no. 11. pp. 24 27.
2. Gavrilov A.A. Definition of the current charge of air in engines with a turbo-supercharging with use of the characteristic of the compressor // Perfection power, economic and ecological parameters ICE: Materials VIII Internal scientific-practical conference, Vladimir State University-Vladimir, 2001. pp. 41–45.
3. Gavrilov A.A., M.S.Ignatov M.S. Designing turbochargers for pressurization of piston engines of internal combustion: Teaching materialstudies. Vladimir state university. Vladimir: Publishing house Vladim. State university, 2009. 88 p. ISBN 978-5-9984-0003-2.
4. Gavrilov A.A., Lazarev V.M., Sharapov N.A., Efros V.V. Approximation of the universal characteristic of the compressor in mathematical model of cycle ICE. Fundamental and applied problems of perfection of piston engines. Materials of XII International scientific-practical conference.- Vladimir State University. Vladimir, 2010. pp. 100–104.
5. Dehovich D.A., Ivanov G.I., Kruglov M.G., Morgulis P.C., Perfilov V.G. The aggregates air the combined engines. Under edit. M.G.Kruglov. M.: Mechanical engineering, 1973. 296 p.
6. Kavtaradze R.Z. Theory of piston engines. Special the chapter. The textbook for high schools. M.: Publishing house of MGTU it N.E.Bauman, 2008. 720 p. (ISBN 978-5-7038-3086-4)
7. Khak G., Landkabel. Turbo and compressors: справ. The handbook. M.: Open Company «Publishing house Astrel»: Open Company «Publishing house of nuclear heating plant», 2003. 351 p.
8. Hanin N.S., Aboltin E.V., Ljamtsev B.F., Zajchenko E.N., Arshins H.p.: Automobile engines with a turbo-supercharging. The handbook. M.: Mechanical engineering, 1991. 336 р.

Эффективность совместной работы поршневого двигателя и турбокомпрессора (ТК) при любом виде наддува (импульсном, изобарном или промежуточном) в значительной мере зависит от величины амплитуд колебания давления свежего заряда ps во впускном трубопроводе. Именно в нём проявляется влияние цикличности процессов впуска на колебания давления рабочего тела. Чтобы уменьшить их негативное влияние на показатели совместной работы двигателя и компрессора, в справочной литературе рекомендуется принимать соотношение объёма трубопровода Vs и рабочего объёма Vh цилиндра Vs ≥ Vh. Основными параметрами, определяющими интенсивность колебаний давления, являются объём впускного трубопровода Vs, количество подсоединённых к нему цилиндров и очерёдность их работы [4]. Это подтверждено исследованиями дизеля 2ЧН 10,5/12 с турбонаддувом [1]. Предварительные расчёты показали, что для увеличения мощности от 22,1 до 30 кВт при α = 1,7 требуется расход воздуха Gk ≈ 0,049 кг/с и степень повышения давления в компрессоре πk = 1,6. На характеристике компрессора TD03-06G этому режиму соответствует КПД ηk = 0,73. При испытаниях с Vs = 0,91Vh ≈ 1 дм3 дизель показал неудовлетворительные результаты (Ne = 24 кВт, ge = 290 г/(кВт∙ч)). С увеличением объёма Vs до 6,5 дм3 показатели дизеля повысились до Ne = 30 кВт при ge = 230 г/(кВт∙ч). При этом было установлено, что чем больше колебания давления pk(πk) и расхода воздуха Gk на впуске, тем менее достоверными являются результаты согласования совместной работы двигателя и ТК по средним за цикл показателям.

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

Методы исследования. Аналитическое исследование.

Результаты исследования и их обсуждение

При моделировании цикла в двигателе используются только входные и выходные показатели: расходы, давления, КПД, частота вращения ротора, т.е. агрегат наддува представляется в виде «черного ящика». Передача энергии от выпускных газов к свежему заряду через турбокомпрессор происходит по схеме: EТ → EТК → EК. Здесь EТ – кинетическая энергия турбины; EТК – кинетическая энергия турбокомпрессора; EК – кинетическая энергия компрессора [5].

При работе двигателя на заданном режиме ротор турбокомпрессора имеет значительную величину кинетической энергии EТК. В течение цикла он получает от выпускных газов через турбину добавочную энергию EТ и через компрессор в количестве EК передает ее свежему заряду. Для промежутка времени dτ этот процесс описывается уравнением

Eqn40.wmf (1)

Кинетическая энергия ротора турбокомпрессора (ТК) EТК при установившемся вращении

Eqn41.wmf (2)

Момент инерции ротора JТК(кг∙м2) для консольного колеса радиально-осевой турбины определяется по эмпирической формуле

Eqn42.wmf

где dТ – диаметр рабочего колеса турбины, см; ωТК = πnТК/30 – угловая скорость вращения ротора ТК, рад/с.

Для конечного промежутка времени Δτ энергия Eт, передаваемая газами турбине, равна [7]

Eqn43.wmf (3)

где Hт – удельная адиабатная работа газа в турбине

Eqn44.wmf

pТО – давление газа за турбиной; pт – давление перед турбиной (принято равным давлению газов в выпускном трубопроводе pт = pp); kт – показатель политропы расширения газа в турбине; ηт – КПД турбины, работающей в условиях переменного давления pт var, для сложившихся к рассматриваемому моменту времени расходу газа Gт и работы Нт, определяется по экспериментальным характеристикам турбины [8].

Анализ работы ТК показывает, что энергия, получаемая турбиной от выпускных газов Eт за цикл, не превышает 2,0 % от кинетической энергии Етк. Вследствие этого переменное давление газа перед турбиной приводит к незначительным, менее 1 %, колебаниям частоты вращения ротора турбокомпрессора nтк. Незначительные изменения энергии Eтк за цикл обусловливают и небольшие колебания значений энергии EК, передаваемой компрессором свежему заряду. Это позволяет вести расчёт наддува в периоды отсутствия подачи газов на турбину и впуска в цилиндры в одно и двух цилиндровых двигателях [3, 4].

Для определения текущего значения энергии EК введем коэффициент пропорциональности (запаса кинетической энергии), вычисляемый по средним за цикл показателям наддува на расчётном режиме работы двигателя

Eqn45.wmf

Тогда количество энергии, полученной компрессором от ротора ТК, обладающего к данному моменту времени энергией Етк + Eт, можно вычислить по уравнению

Eqn46.wmf (4)

Уравнение для вычисления энергии, переданной в компрессоре свежему заряду, имеет вид

Eqn47.wmf (5)

где GК, ηК – расход воздуха и КПД компрессора; HК удельная работа сжатия воздуха

Eqn48.wmf

k – показатель политропы сжатия в компрессоре; pОК, pК – давления воздуха перед компрессором и после компрессора (принято равным давлению воздуха во впускном трубопроводе pК = ps).

Для определения текущих значений GК и ηК при работе двигателя с турбонаддувом уравнение (5) решается относительно отношения (GК/ηК)d

Eqn49.wmf (6)

pic_9.tif

Рис. 1. Характеристика центробежного компрессора

При расчёте цикла по углу поворота коленчатого вала (°пкв) вместо Δτ принимается интервал Δφ = ωΔτ, где угловая скорость ω,°п.к.в./с (ω = ωr180/π = 6nd, ωr, рад/с; nd – частота вращения коленчатого вала, мин–1.

По величине отношения (GК/ηК)d значения GК и ηК определяются с использованием характеристики компрессора.

Методика её применения базируется на уникальной особенности отношения (Gk/ηk)x, вычисленного при πk = const на характеристике компрессора (рис. 1) – возрастать с увеличением расхода воздуха Gk (рис. 2). Теоретически это отношение представляет предельный расход идеального газа на данном режиме работы центробежного компрессора. Для практического использования характеристика компрессора (например, на рис. 1) перестраивается в графики зависимости (Gk/ηk)x = f(Gk, πk).

Значения Gk и ηk берутся в точках пересечения линии πk = const с кривыми ηk = const, начиная с границы помпажа. Значения (Gk/ηk)x аппроксимируются полиномами не выше четвёртой степени (Gk/ηk)x = f(Gk, πk). Графики их для трёх πk (1,5; 1,8; 2,1) приведены на рис. 2.

pic_10.tif

Рис. 2. График зависимости Gк/ηк от Gк при πk равном 1,5; 1,8; 2,1 (см. рис. 1)

В программе расчёта цикла двигателя эти зависимости оформляются в виде матрицы коэффициентов полиномов. Для лучшего восприятия используемой методики приведём фрагмент ввода параметров с характеристикой при πk = 1,8. Для удобства оформления исходных данных значения расхода воздуха, снятые с характеристики, увеличены в соответствии с масштабом шкалы mG, т.е. GK/mG (для формата 00.00).

KG = 12 PIK = 1.8000

4.3 4.9 6.1 7.0 8.1 8.5 11.2 11.8 12.8 13.7 14.2 14.9

.64 .65 .68 .70 .72 .73 .73 .72 .70 .68 .65 .60

Здесь принято: KG – число парных значений расходов Gk (1-я строка) и ηk (2-я строка);

PIK – степень повышения давления в компрессоре πk. Результаты расчёта отношений Gk/ηk///, их аппроксимации функцией (Gk/ηk)x = f(Gk), коэффициенты полинома и оценка погрешностей (отклонения) для приведённого фрагмента имеют вид:

KP = 16 KG = 12 PIK = 1.8000

K Коэф.полин.

0 0.84661E-01

1 -0.26380E + 01

2 0.77620E + 02

3 -0.67986E + 03

4 0.21924E + 04

N Gk Gk/кпд (Gk/кпд) = f(Gk) Oтклонение

1 0.43000E-01 0.67188E-01 0.68190E-01 -0.10030E-02

2 0.49000E-01 0.75385E-01 0.74421E-01 0.96332E-03

3 0.61000E-01 0.89706E-01 0.88611E-01 0.10948E-02

4 0.70000E-01 0.10000E + 00 0.99791E-01 0.20889E-03

5 0.81000E-01 0.11250E + 00 0.11332E + 00 -0.82359E-03

6 0.85000E-01 0.11644E + 00 0.11817E + 00 -0.17288E-02

7 0.11200E + 00 0.15342E + 00 0.15270E + 00 0.72142E-03

8 0.11800E + 00 0.16389E + 00 0.16219E + 00 0.16954E-02

9 0.12800E + 00 0.18286E + 00 0.18148E + 00 0.13752E-02

10 0.13700E + 00 0.20147E + 00 0.20429E + 00 -0.28146E-02

11 0.14200E + 00 0.21846E + 00 0.21997E + 00 -0.15119E-02

12 0.14900E + 00 0.24833E + 00 0.24651E + 00 0.18218E-02

Здесь KP = 16 ‒ порядковый номер значения πk, при котором вводились Gk и ηk, а столбец «Отклонение» позволяет оценить корректность их ввода. Для границы помпажа дополнительно формируются полиномы

(Gk)pom = f(πk) и (Gk/ηk)pom = f(πk). (7)

В программе расчёта цикла реализован следующий алгоритм определения текущих значений Gkи ηk:

1. Вычисляется расход газов через турбину GT.

2. Рассчитываются последовательно параметры HT, ET, EТК, HК, EК и отношение (Gk/ηk)d при степени повышения давления πk на данный момент времени.

3. Определяется положение точки совместной работы двигателя и компрессора (совмещённая точка) относительно границы помпажа при данном πk. Если (Gk/ηk)2 ≤ (Gk/ηk)pom, то вычисляются значения расхода (Gk)pom и КПД (ηk)pom на границе помпажа и выдаётся соответствующее сообщение.

4. Проверяется положение совмещённой точки относительно нижнего и верхнего значений πk, принятых для аппроксимации характеристики компрессора. Если πk выходит за принятые значения, то Gkи ηk вычисляются на соответствующей границе.

5. Если совмещённая точка не вышла за пределы, установленные в пунктах 3 и 4, то для данной степени πk находятся ближайшие на характеристике компрессора πk1 ≤ πk ≤ πk2.

6. По значениям коэффициентов в полиномах функций (Gk/ηk)x = f(Gk) для πk1 и πk2 интерполируются коэффициенты полинома функции (Gk/ηk)x = f(Gk) для πk.

7. По сформированным коэффициентам для значений Gkx, начиная от границы помпажа, с принятым шагом вычисляются значения полинома (Gk/ηk)x = f(Gkx) и сравниваются со значением (Gk/ηk)d, определённым по уравнению (6). Как только будет достигнуто соотношение (Gk/ηk)x ≥ (Gk/ηk)d, вычисляется КПД компрессора ηk = Gkx/(Gk/ηk)x и текущий расход воздуха Gk = ηk(Gk/ηk)d.

Выводы

Практика использования предлагаемой методики показала, что в районе средних значений πk на характеристике компрессора погрешность определения текущих значений GKи ηK не превышает 2 %, а при πk, близких к πk1 или πk2 может повышаться до 5 % и более. Методика обеспечивает получение более достоверных результатов согласования совместной работы двигателя и турбокомпрессора при подсоединении к впускному трубопроводу менее трёх цилиндров, а при неравномерном порядке работы цилиндров – не более четырёх цилиндров.

Рецензенты:

Кобзев А.А., д.т.н., профессор, заведующий кафедрой, ФГБОУ «Владимирский государственный университет им. Александра Григорьевича и Николая Григорьевича Столетовых», г. Владимир;

Кульчицкий А.Р., д.т.н., профессор, главный специалист, ООО «Завод инновационных продуктов», г. Владимир.

Работа поступила в редакцию 24.05.2013.