CC BY
72
ISSN 1992-6502 (P ri nt)_
2017. Т. 21, № 1 (75).С. 72-79
Ъьомшь QjrAQnQj
ISSN 2225-2789 (Online) http://journal.ugatu.ac.ru
УДК 621.45.037
Анализ влияния малоразмерности
на параметры эффективности осевого компрессора
а. а. Бобрик 1, а. ю. Ткаченко2, в. с. Кузьмичев3
1ЬоЬпк000а1@таМ.ш, 2 [email protected], [email protected] ФГАОУ ВО «Самарский национальный исследовательский университет имени академика С. П. Королева»
Поступила в редакцию 01.12.2016
Аннотация. В настоящее время малоразмерные газотурбинные двигатели находят все большее применение. Для двигателей средних и больших размеров накоплен опыт использования в процессе проектирования математических моделей разного уровня сложности, в частности, одномерной модели осевого компрессора. На основе численных расчетов оценивается возможность применения этой модели для малоразмерных газотурбинных двигателей. Переход к малоразмерности осуществляется за счет масштабирования исходной модели, исходя из сохранения положения рабочей точки на напорной характеристике компрессора. При последовательном масштабировании оцениваются границы применимости используемой модели и ее адекватность в сравнении со статистическими данными по уровню эффективности работы малоразмерных газотурбинных двигателей.
Ключевые слова: одномерная модель, малоразмерный газотурбинный двигатель, осевой компрессор, масштабирование.
ВВЕДЕНИЕ
В настоящий момент малоразмерные газотурбинные двигатели (МГТД) находят все большее применение в различных отраслях. В авиации МГТД используются в качестве силовых установок (СУ) беспилотных летательных аппаратов (ЛА), а также в качестве вспомогательных силовых установок (ВСУ) самолетов и вертолетов. В перспективе планируется использовать МГТД в составе распределенных силовых установок (РСУ) средне- и дальнема-гистральных самолетов [1].
В современной практике на этапе концептуального проектирования двигателя используются одномерные (1-0) термодинамические модели узлов и двигателя в целом. Данные модели позволяют быстро оценить эффективность двигателя и его составных узлов до начала их детального проектирования.
Математические модели МГТД, используемые на этапе концептуального проектирования, по существу не отличаются от аналогичных моделей двигателей средней и большой размерности. Однако в них необходимо учесть особенности, которые присущи МГТД.
В данном исследовании рассматриваются следующие актуальные вопросы, связанные с
разработкой эффективных методов проектирования МГТД:
• оценка адекватности 1-0 модели осевого компрессора (ОК), созданной применительно к полноразмерным двигателям, в диапазоне размерностей, характерных для МГТД;
• сравнение расчетных результатов по 1-0 модели ОК с обобщенными статистическими данными;
• разработка способов учета особенностей МГТД, необходимых для использования существующей 1 -0 модели ОК.
ОПИСАНИЕ РАСЧЕТНОЙ МОДЕЛИ
В работе исследуется влияние размерности на параметры эффективности осевого компрессора при уменьшении габаритов его проточной части за счет пропорционального масштабирования осевых и диаметральных размеров. В компрессоре при уменьшении габаритов увеличиваются относительные величины радиальных зазоров, толщина пограничного слоя в лопаточных венцах. Критерий Рейнольдса выходит за пределы автомодельности. Все эти факторы приводят к снижению КПД компрессора. Чаще всего на этапе концептуального проектирования двигателя это учитывается статистиче-
скими поправками на КПД в зависимости от его размерности.
Исследования проводились на одномерной расчетной модели трехступенчатого ОК с входным направляющим аппаратом, созданной в САЕ - системе АСТРА [2].
Компрессор работает при стандартных атмосферных условиях (САУ) на земле:
рн = р* = 101325Па; Тн = Т* = 288К .
Исходная модель компрессора характеризуется параметрами рабочего процесса, представленными в табл. 1.
Таблица 1
Параметры рабочего процесса компрессора исходной размерности
Физический расход воздуха через компрессор, Св, кг/с 100
Расход воздуха, приведенный к параметрам рабочего тела за последней ступенью компрессора, Св пр к, кг/с 53,15
Частота вращения ротора, об/мин 5000
Суммарная степень повышения давления в компрессоре, пк 2,115
Изоэнтропический КПД компрессора, Пк 0,902
Геометрические параметры компрессора, характеризующиеся радиусами (втулочным, средним и периферийным) на входе в соответствующий лопаточный венец (ЛВ), а также, параметры выходного сечения представлены в табл. 2. Радиусы на входе в каждый последующий ЛВ равны выходным радиусам предыдущего:
я.вх (г+1) = явых (0 .
Таблица 2 Геометрические параметры компрессора исходной размерности
Обозначе- Втулоч- Сред- Периферий-
ние ный ний ный радиус,
лопаточ- радиус, радиус, ЯК, мм
ного вен- Явт, мм ЯСР, мм
ца
ВНА 196 350,5 505
РК-1 246,5 366,9 487,3
НА-1 310 403,1 496,2
РК-2 310 403,1 496,2
НА-2 348 421,1 494,3
РК-3 355 424,4 493,8
НА-3 356 422,1 488,2
Выходное
сечение 340 407,8 475,6
Геометрия компрессора масштабировалась в несколько этапов - от исходной до уменьшенной в 20 раз. В качестве параметра масштабиро-
вания был принят коэффициент масштабирования, Км. По результатам масштабирования получено несколько моделей компрессора со следующими значениями КМ: Км ={1; 0,8; 0,6; 0,4; 0,2; 0,15; 0,1; 0,05} .
Для примера на рис. 1 представлены проточные части исходного варианта компрессора и компрессора, уменьшенного в 10 раз.
\ии
I I I иг
Рис. 1. Проточные части исходного компрессора и уменьшенного компрессора (КМ = 0,1)
Абсолютные значения радиальных зазоров у рабочего колеса (РК) для всех моделей составляют:
= 1 мм; дрз2 = 1,2 мм; Зрз3 = 0,8 мм.
ПРИНЦИП МАСШТАБИРОВАНИЯ
Масштабирование компрессора проводилось из принципа сохранения положения рабочей точки на его напорной характеристике, то есть из постоянства треугольников скоростей для каждой из рассчитываемых моделей.
Далее приведены основные зависимости для параметров, характеризующих компрессор:
F = я(Я2К -ЯВТ ) = я(Як -ЯВТ )(Як + ЯВТ ) =
= ж- И-л • 2ЯСР , (1)
с = ш• р* • д(X)-^
и =
4г
я Я • п
(2)
30 (3)
где ^ - осевая площадь рассматриваемого сече* гг*
ния; р , Т - полные давление и температура рабочего тела; - относительная плотность тока; т - численный коэффициент, для воздуха т = 0,0405 (кг х К/Дж)0,5; и - окружная скорость на >том радиусе РК.
Пусть исходная модель компрессора масштабируется в к раз (Км = к). Радиусы всех ступеней, а также входного устройства изменяются пропорционально Км. Исходя из зависимости (1), площадь изменяется пропорционально изменению высоты лопатки и среднего радиуса, то есть пропорционально к2 раз.
Расход воздуха (2) на входе в компрессор изменяется пропорционально площади, так как полные давление и температура зависят только от внешних условий, которые не изменяются. Таким образом, расход изменяется пропорционально квадрату Км в k2 раз.
Чтобы обеспечить постоянство треугольников скоростей при масштабировании моделей, необходимо сохранить неизменной величину окружной скорости, которая изменяется пропорционально радиусу и частоте вращения ротора. Тогда, исходя из (3), частота вращения должна меняться обратно пропорционально радиусу, то есть в 1/k раз.
На следующей схеме проиллюстрировано описанное изменение параметров: KM i (k) ^ R i (k) ^ Ft i( k ) ^ GBX i (k2);
R i (k) ^ Ui i (k); U = const при n i ^ 1 j ot (k).
Важно отметить, что постоянство треугольников скоростей для компрессоров разной размерности не выдерживается из-за влияния дополнительных потерь, за счет которых изменяется картина течения потока, в том числе и углы атаки, в результате чего происходит рассогласование ступеней, что влияет на параметры ОК.
ОПТИМИЗАЦИЯ ПАРАМЕТРОВ КОМПРЕССОРА
Расчет каждого из масштабируемых вариантов компрессора в САЕ - системе АСТРА включал первичный расчет, в котором учитывалось только пропорциональное изменение геометрии, и оптимизированный, в котором дополнительно подбирались оптимальные величины характерных углов ЛВ с целью обеспечения максимального КПД. Расчетные параметры ЛВ, используемые в одномерной модели, приведены на рис. 2.
Рис. 2. Геометрические параметры лопаточного венца РК На рис. 2 представлены следующие параметры:
• с - хорда лопатки;
• 5 - шаг лопаточного венца;
• вш - конструктивный лопаточный угол на входе в венец;
• в2л - конструктивный лопаточный угол на выходе из венца;
• У - угол между продольной осью компрессора и касательной к средней линии профиля лопатки на входе в венец;
• 0 - угол между касательными к средней линии профиля на входе и выходе из венца.
Оптимизация проводилась путем подбора значений углов у и 0:
7 = 90о,
© = г-(90о-Ал) = Г+А "90.
(4)
(5)
Подставляя (4) в (5), получаем:
© = 90°-%л +Ал - 90° = Ал -Ьл. (6)
Таким образом, подбор значения угла у позволил определить оптимальный угол атаки, а подбор угла 0 - оптимальную работу сжатия, подводимую к ступени. Данные обозначения соответствуют РК. Для направляющих аппаратов (НА) вместо углов в используются углы а.
РЕЗУЛЬТАТЫ РАСЧЕТОВ
Серия расчетов при различных коэффициентах масштабирования проводилась в следующей последовательности:
М 1:1 ^М 0,8:1
(исх / первич)
М 1:1
(опт)
0,05: 1
(первич) ^.....^ \ji\j-j . -i (опт)•
Значения углов у и 0 для первичного расчета каждой следующей масштабной модели брались из оптимизационного расчета предшествующей модели большего масштаба. Затем проводился пересчет углов отмасштабированной модели. Так как при оптимизации изменяются параметры ЛВ, то изменяется и кинематика потока и, соответственно, треугольники скоростей.
Численные результаты расчетов представлены в табл. 3, 4, а сравнение параметров, полученных двумя видами расчетов, в зависимости от Км проиллюстрировано на рис. 3-7.
Рис. 3. Зависимость приведенного расхода воздуха Gb пр к от Км
Таблица 3
Результаты первичных расчетов для
моделей ОК различной размерности
Коэффициент масштабирования, Км 1 0,8 0,6 0,4 0,2 0,15 0,1 0,05
Расход воздуха через ОК, Gв, кг/с 100 64 36 16 4 2,25 1 0,25
Частота вращения ротора, п, об/мин 5000 6250 8333 12500 25000 33333 50000 100000
Приведенный расход воздуха за ОК, Gк пр, кг/с 53,154 34,226 19,469 8,584 2,171 1,246 0,571 0,152
Степень повышения давления в К, пК 2,115 2,101 2,074 2,095 2,071 2,024 1,96 1,842
Изоэнтропический КПД компрессора, п к 0,902 0,899 0,897 0,891 0,874 0,867 0,844 0,756
Закрутка потока в РК-1, ДСи, м/с 135,45 135,35 138,95 141,51 133,69 133,97 135,9 134,82
Число Рейнольдса в РК-3, Яе 1817162 1450731 1089135 690428 349972 244939 157818 72641
Результаты расчетов с оптимизацией для моделей ОК различной р Таблица 4 азмерности
Коэффициент масштабирования, Км 1 0,8 0,6 0,4 0,2 0,15 0,1 0,05
Расход воздуха через К, Gв, кг/с 100 64 36 16 4 2,25 1 0,25
Частота вращения ротора, п, об/мин 5000 6250 8333 12500 25000 33333 50000 100000
Приведенный расход воздуха за ОКДк пр, кг/с 53,221 34,37 19,368 8,396 2,19 1,252 0,555 0,147
Степень повышения давления в К, пК опт 2,112 2,089 2,086 2,15 2,05 2,012 2,024 1,904
Изоэнтропический КПД компрессора, п к опт 0,904 0,903 0,900 0,893 0,878 0,868 0,847 0,772
Закрутка потока в РК-1, ДСи, м/с 135,47 139,13 141,74 134,24 134,29 136,48 136,83 137,92
Число Рейнольдса в РК-3, Яе 1816791 1456061 1040395 706453 328236 238873 150303 74081
Рис. 4. Зависимость степени повышения давления от коэффициента масштабирования
Рис. 6. Зависимость критерия Рейнольдса от коэффициента масштабирования
Рис. 5. Зависимость закрутки потока от коэффициента масштабирования
Рис. 7. Зависимость изоэнтропического КПД от коэффициента масштабирования
Как видно из таблиц и рисунков, приведенный расход воздуха для обоих вариантов расчета фактически не различается, степень повышения давления зависит от углов ЛВ, определяемых в результате оптимизации. Рис. 5 иллюстрирует изменение подобия треугольников скоростей при масштабировании.
КПД (см. рис. 7) существенно уменьшается при Км< 0,2. Данную величину (Км= 0,2) при исходном значении расхода воздуха через компрессор Gв = 100 кг/с можно принять за границу малоразмерности.
СРАВНЕНИЕ РАСЧЕТНЫХ РЕЗУЛЬТАТОВ СО СТАТИСТИЧЕСКИМИ ЗАВИСИМОСТЯМИ
Размерность двигателя в целом принято оценивать по различным параметрам (тяге, расходу воздуха на входе в двигатель, приведенному расходу воздуха по параметрам за компрессором). В данном исследовании в качестве критерия размерности принят расход воздуха, приведенный по параметрам рабочего тела за компрессором [3]. Для выбранного критерия существует статистическая зависимость (7), полученная на основе обобщения экспериментальных данных [4].
Данная зависимость определяет поправку на политропический КПД компрессора Ппол.к исх в зависимости от его размерности.
Таблица 5
Сравнение величин расчетного и статистического КПД для серии расчетов без оптимизации
Политропический КПД, Ппол. к 0,912 0,909 0,907 0,901 0,886 0,879 0,858 0,776
Статистическая поправка на политропический КПД, ДПпол. к х 102 0,565 0,589 0,641 0,791 1,1585 2,2374 4,4564 15,706
Статистический КПД, Ппол. к стат 0,906 0,906 0,905 0,904 0,896 0,888 0,866 0,755
§Ппод. к, % 0,62 0,34 0,17 0,28 1,15 1,03 0,98 2,71 Таблица 6
Сравнение величин расчетного и статистического КПД для серии расчетов с оптимизацией
Политропический КПД, Ппол. к 0,913 0,913 0,910 0,904 0,889 0,881 0,861 0,791
Статистическая поправка на политропический КПД, ДПпол. к х 102 0,565 0,589 0,641 0,797 1,576 1,2365 4,6810 16,223
Статистический КПД, Ппол. к стат 0,908 0,907 0,907 0,905 0,897 0,890 0,866 0,751
§Ппол. к, % 0,62 0,60 0,31 0,11 0,91 1,02 0,63 5,08
0,02308
Л'л.к =—-+ °,0522.
°в пр К (7)
Взаимосвязь изоэнтропического и политропического КПД определяется зависимостью (8).
( — 1п л к
Лп
v
пол.К / к-1 Л
' л -1 1п 1 + л-1
Лк
(8)
Рассчитанные по этой зависимости изоэн-тропические КПД были пересчитаны на политропические.
Полученные значения политропического КПД для расчетных моделей и по обобщенной статистической зависимости приведены в табл. 5-6.
Относительная разность между КПД вычисляется по следующей формуле:
пол.к '
'пол.к расч (9)
По результатам расчетов построены графики (рис. 8-9), где сплошными линиями изображаются расчетные значения КПД, а пунктирными - определенные по статистическим зависимостям.
Рис. 8. Зависимости расчетного и статистического КПД от приведенного расхода воздуха без оптимизации ЛВ (логарифмическая шкала)
Рис. 9. Зависимости расчетного и статистического КПД от приведенного расхода воздуха с оптимизацией ЛВ (логарифмическая шкала)
Из анализа полученных результатов видно, что для всех масштабируемых моделей (кроме модели наименьшего масштаба) статистические данные не отличаются от расчетных более чем на 1 %.
Таким образом, созданная математическая модель при введении поправочных коэффициентов может использоваться для расчета компрессоров малоразмерных ГТД.
ЗАКЛЮЧЕНИЕ
С помощью одномерной термодинамической модели проведен численный анализ влияния малоразмерности на параметры эффективности осевого компрессора при пропорциональном масштабировании его осевых и диаметральных размеров. В качестве параметра масштабирования был принят коэффициент масштабирования КМ, изменяющийся от 1 до 0,05.
Исходный вариант компрессора масштабировался из условия сохранения рабочей точки на напорной характеристике, однако, вследствие оптимизации углов у и 9, корректирующих углы атаки и закрутку потока в лопаточных венцах с целью получения максимального КПД, положение рабочей точки изменялось.
По результатам расчетов определено влияние размерности компрессора на величину его КПД. Показано, что изоэнтропический КПД компрессора при уменьшении его размеров в 20 раз и оптимизации уменьшается на 13% (от 0,9 до 0,77).
Показано, что расчетные значения КПД осевого компрессора, определенные по его одномерной модели, согласуются с обобщенными статистическими зависимостями, расхождение не превышает 1%.
Одномерная математическая модель осевого компрессора, применяемая при проектировании полноразмерных двигателей, может использоваться и для расчета компрессоров малоразмерных газотурбинных двигателей.
СПИСОК ЛИТЕРАТУРЫ
1. Скибин В. А., Соломин В. И., Палкин В. А. Работы ведущих авиадвигателестроительных компаний в обеспечении создания перспективных авиационных двигателей. М.: ЦИАМ, 2010. 677 с. [ V. A. Skibin, V. I. Solonin, V. A. Pal-kin, The works of the leading aircraft engine companies in ensuring the development of advanced aircraft engines, (in Russian). Moscow: Central Institute of Aviation Motors "CIAM", 2010. ]
2. Формирование виртуальной модели рабочего процесса газотурбинного двигателя в CAE системе «АСТРА» / В. С. Кузьмичев [и др.] // Электронный журнал «Труды
МАИ». 2013. №67. С. 1-15. [ V. S. Kuz'michev, et al., "Development of the virtual model of the working process of gas turbine engine using the computer-aided system ASTRA", (in Russian), in Elektronnyj zhurnal «Trudy MAI», 2013, no. 67, pp. 1-15. ]
3. Скибин В. А. Машиностроение. Энциклопедия. Самолеты и вертолеты. Том 4-21. Книга 3. М.: Машиностроение, 2010. 720 с. [V. A. Skibin, Engineering. Encyclopedia. Airplanes and helicopters. Vol. 4-21. Book 3, (in Russian). Moscow: Mashinostroenie, 2010. ]
4. Григорьев В. А., Зрелов В. А., Игнаткин Ю. М. и др. Вертолетные газотурбинные двигатели. М.: Машиностроение, 2007. 491 с. [ V. A. Grigor'ev, V. A. Zrelov, Y. M. Ignatkin, et al., Helicopter gas turbine engines, (in Russian). Moscow: Mashinostroenie, 2007. ]
ОБ АВТОРАХ
БОБРИК Алексей Александрович, асп. каф. теории двигателей летательных аппаратов. Дипл. инж. (Самарский Университет им. академика С. П. Королева, 2016). Готовит дисс. о методах и средствах концептуального проектирования распределенных силовых установок летательных аппаратов на базе малоразмерных газотурбинных двигателей.
ТКАЧЕНКО Андрей Юрьевич, доц. каф. теории двигателей летательных аппаратов. Дипл. инж.-констр. (СГАУ, 2005). Канд. техн. наук по тепловым, электроракетным двигателям и энергоустановкам ЛА (СГАУ, 2009). Иссл. в обл. разработки автоматизированных средств проектирования ГТД.
КУЗЬМИЧЕВ Венедикт Степанович, проф. каф. теории двигателей летательных аппаратов. Д-р техн. наук по тепловым, электроракетным двигателям и энергоустановкам ЛА (СГАУ, 2000). Иссл. в обл. разработки автоматизированных средств проектирования ГТД.
METADATA
Title: Analysis of influence of reducing axial compressor's size.
Authors: A. A. Bobrik1, A. U. Tkachenko2, V. S. Kuz'michev3
Affiliation: Samara University, Russia.
Email: [email protected], [email protected], [email protected].
Language: Russian.
Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), vol. 21, no. 1 (75), pp. 72-79, 2017. ISSN 2225-2789 (Online), ISSN 1992-6502 (Print).
Abstract: Small-size gas turbine engines are finding larger application in present. Experience of using mathematical models of engines and parts of engines has accumulated preference for medium and big-size gas turbine engines. Appraise of possibility using mathematical models axial compressor of medium and big-size engines for small-size engines is a main point of this article.
Key words: one-dimensional model; small-size gas turbine engine; axial compressor; scale.
About authors:
BOBRIK, Alexey Aleksandrovich, postgraduate student of Department theory of aircraft engine. Diploma design engineer (Samara University, 2016). Prepares diss. methods and means of the conceptual design of distributed power plants which based on small-size gas turbine engines.
TKACHENKO, Andrey Yurievich, associate professor of Department theory of aircraft engines. Diploma design engineer (SSAU, 2005). PhD for thermal, electric propulsion and power installations of aircrafts (SSAU, 2009). Research in the development of automated design tools GTE.
KUZ'MICHEV, Venedikt Stepanovich, professor of Department theory of aircraft engines. Doctor of technical science for thermal, electric propulsion and power installations of air-crafts (SSAU, 2000). Research in the development of automated design tools GTE.
