ISSN 1992-6502 (Print)
2021. Т. 25, № 4 (94). С. 38-45
Вестник УГАТУ
http://journal.ugatu.ac.ru
ISSN 2225-2789 (Online)
УДК 519.248
doi 10.54708/19926502_2021_2549438
Методика выбора оптимальной геометрии циклонного
пылезащитного устройства авиационного ГТД
а. с. Гишваров 1, а. х. Рахимов 2
2
1 kafedra.ad@mail.ru, 2 mr.abdusattor@list.ru ФГБОУ ВО «Уфимский государственный авиационный технический университет» (УГАТУ)
Поступила в редакцию 22.04.2021
Аннотация. Рассматривается методика выбора оптимальной геометрии циклонного пылезащитного устройства (ПЗУ) авиационного двигателя. Методика включает следующие основные этапы: формирование исходных данных, выбор частных показателей эффективности ПЗУ и их нормирование, построение регрессионных уравнений связи вида «параметры эффективности ПЗУ - геометрические размеры ПЗУ», формирование обобщенной функции цели, определение множества Парето-оптимальных решений и выбор окончательного варианта оптимальной геометрии ПЗУ. Приводится численный пример выбора оптимальной геометрии циклонного ПЗУ.
Ключевые слова: пылезащитное устройство (ПЗУ); показатели эффективности; обобщенная функция цели; Парето-оптимальное множество; методика.
Пылевая эрозия элементов проточной части авиационных ГТД и отложения пыли вызывают снижение мощности, недопустимый рост температуры газа перед турбиной, поломки ослабленных эрозией лопаток, помпаж и другие отказы. Попавшие в двигатель частицы пыли не только изменяют форму профилей лопаток, но и оседают во входной части компрессора, что затрудняет надежную работу двигателя вертолета. В результате налипания пыли на лопатках компрессора их собственная частота колебаний может уменьшиться на 45 %, что приводит к недопустимому увеличению действующих в лопатках напряжений и к их обрыву [1, 2].
Эрозионный износ лопатки компрессора приводит к утонению передней и задней кромок, выработке поверхностей корыта, уменьшению значений хорды и максимальной толщины [3].
Такое повреждение лопаток компрессора ведет к изменению режима их обтекания и, соответственно, к изменению параметров, определяющих условия работы профиля в потоке:
- угла между средней геометрической скоростью и хордой профиля;
- коэффициента подъемной силы;
- угла набегания, определяемого как угол между вектором скорости воздуха перед решеткой и касательной к средней линии профиля у передней кромки;
- угла отставания потока;
- угла поворота потока в решетке.
Увеличение угла набегания за счет эрозионного износа входной кромки может приводить к росту местных скоростей на поверхности профиля и к появлению сверхзвуковых скоростей потока.
Гидравлические потери при движении воздуха через канал решетки с небольшими скоростями до и после повреждения поверхностей профиля лопатки сводятся к трем видам:
а) профильные потери, связанные с образованием вторичных течений;
ВВЕДЕНИЕ
б) концевые потери, связанные с образованием пограничного слоя на стенках, ограничивающих канал решетки по высоте;
в) вторичные потери, связанные с возникновением протеканий в пограничном слое.
Увеличения суммарных потерь полного давления, обусловленное ухудшением характеристик отдельных ступеней при эрозионном износе лопаток, приводит к падению напорно-сти ступеней и уменьшению плотности воздуха по тракту компрессора. Происходит смещение характеристик в сторону меньших расходов, граница устойчивости смещается к рабочей линии. На границе помпажа при малых частотах вращения ротора работа первых ступеней компрессора неустойчива, а режим работы последних ступеней близок к режиму «запирания» [3].
Увеличения удельного расхода топлива и температуры газов за турбиной происходит за счет ухудшения аэродинамики проточной части вследствие увеличения радиальных зазоров между рабочими лопатками и корпусом компрессора, износа и шероховатости поверхности лопаток.
Снижение запаса газодинамической устойчивости из-за эрозионного износа на определенную величину может способствовать возникновению помпажа у одних двигателей, в то время как у других этого явления при остальных одинаковых условиях не возникает. Одной из причин этого является разброс по характеристикам как новых, так и прошедших ремонт ГТД.
При эрозионном износе элементов конструкции двухконтурных двигателей с увеличением наработки также ухудшается их газодинамическая устойчивость. Наблюдается смещение линии рабочих режимов высокого давления, смещение напорных ветвей по расходу воздуха, изменение границы устойчивой работы. Например, изменение КПД компрессора или турбины каскада высокого давления двигателя типа Д-30КУ на 1 % вызывает смещение линии рабочих режимов на 0,72 %.
Поэтому на практике используются пылезащитные устройства авиационных двигателей, в которых для отделения пыли (песка) от воздуха используется их инерционность.
АКТУАЛЬНОСТЬ
ПЗУ считается важнейшим элементом для сохранения работоспособности вертолетных двигателей при эксплуатации в условиях запыленного воздуха.
Оценка эрозии элементов газо-воздушного тракта ГТД является сложной задачей, поскольку зависит от многих факторов: от твердости и химического состава частиц пыли, от фракционного состава, концентрации песка, режима работы двигателя и др. Эффективность ПЗУ оценивается, как правило, двумя параметрами: степенью очистки воздуха и потерей давления в воздушном тракте ПЗУ. Поэтому разработка методики выбора оптимальной геометрии ПЗУ является сложной многокритериальной и многофакторной задачей [4].
ПОСТАНОВКА ЗАДАЧИ
Целью данной работы является разработка метода выбора оптимальной геометрии ПЗУ с учетом многокритериальности и многофакторности решаемой задачи.
МЕТОДИКА ВЫБОРА ОПТИМАЛЬНОЙ ГЕОМЕТРИИ ПЗУ
Методика включает шесть основных этапов.
Этап 1. Формирование исходных данных. На данном этапе проводится:
- обоснование программного обеспечения, используемого для моделирования двухфазного потока «воздух - частицы песка». При этом рассматриваются Лагранжевый, Эйлеровый и комбинированный методы [5];
- выбор режима работы двигателя;
- задание характеристик песка (количества, скорости, диаметра частиц и др.);
- области определения геометрических размеров ПЗУ и др.
Этап 2. Выбор частных параметров и критериев эффективности ПЗУ.
На практике в качестве частных параметров эффективности рассматриваются только два показателя [6-8]:
- степень очистки воздуха - п, %;
- потеря давления в ПЗУ - АР, Па.
В данной работе авторами введено в рассмотрение еще одного параметра эффективности - массы ПЗУ (М, кг), то есть выбор оптимальной геометрии ПЗУ Хсрг = (х1, хг, ...,хк) проводится с учетом трех частных критериев эффективности:
Ф, (X) = nopt (X)= max n (X); Ф2 (X) = APopt (X) = min Ар ( X); Фз (X)= Mopt ( X) = min Ыг (X),
(1)
гдеX = (xi, Х2, ..., xk) - вектор геометрических параметров ПЗУ.
Численные значения параметров n, ЛР и M определяются моделированием двухфаз- ного потока «воздух - частицы песка» с применением программного комплекса (ПК) Ansys Fluent. Область значений X = (xi, Х2, ..., xk) задается в виде системы двухсторонних ограничений вида:
jC ^ JC • ^ JC тгш — г — тахг
г = I к,
где Хгшах, хгшш - максимальное и минимальное значения /-го фактора.
Этап 3. Нормирование частных параметров и критериев эффективности ПЗУ. При нормировании вместо «натурального» параметра вводится его отношение к некоторой нормирующей величине, измеряемой в тех же единицах, что и сам параметр [9-11]. В результате операции нормирования и приведения к единой шкале измерения все частные критерии эффективности приобретают безразмерный вид с диапазоном изменения от 0 до 1:
Ф f
Ф
Ф
X =т£, X = тахг|г X G0...1; 2 X =A/J X = ттД/> X G0...1; н х =M"pt X = minМг X G0...1.
(2)
Критерии эффективности вида (2) можно преобразовать к виду, при котором эффективность параметров п, АР и М многокритериальное^ по принципу чем больше значения показателя, тем выше его эффективность в виде:
Of X = maxг|г X е(0...100)%;
•100 £ (0...100) %;
Of X =тах
Of X =тах
1
X
^^баз X
(3)
М, X
Мблз. X
юое(о..лоо)%,
где АРбаз. и Мбаз. - базовые значения показателей АР и М, соответствующие наиболее наилучшим значениям, к которым необходимо стремиться.
Представление критериев в виде (3) удобно для выбора оптимальной геометрии ПЗУ с применением обобщенной целевой функции аддитивного вида [8, 9]:
ФЕ X = maxOf X = max r|fW +
Z=1
1-
ЛР X
A/f, X
M X
•100 + 1 ! •100
M6». X
(4)
3
Этап 4. Задаются двухсторонние ограничения на геометрические размеры ПЗУ, внутри интервала которых проводится поиск оптимальных значений:
О
х
х . < х < х
1шт 1 1тах
ктш к ктах
Внутри области Ох строится план эксперимента, в котором варьируемыми (независимыми) параметрами являются геометрические параметры ПЗУ (х1, Х2,..., хк), а выходными (зависимыми) параметрами являются значения п, АР и М. По результатам эксперимента определяются уравнения регрессии, которые проверяются на значимость коэффициентов регрессии и адекватность самих уравнений:
п = /1( х1, Х2, ----, ХК )' = /2( х1, х2, ----, хк )' М = /э( х1, Х2, ----, хк )• (5)
Перед планированием эксперимента натуральные значения независимых факторов XI, Х2,хк кодируются по формуле:
X =
х хю
где Х - кодовое значение /-го фактора; х/ - натуральное текущее значение /-го фактора; х/о - начальный (нулевой) уровень фактора; Ах/ - интервал варьирования /-го фактора:
Л„ _ хгшах хгшт ^ " 2
Этап 5. Строится множество Парето-оптимальных значений функции Ф (X) = / [(ФГ (X) ,Ф^ (X) ,Ф^ (X)] с применением обобщенной целевой функции аддитивного вида (4) и коэффициентов важности а/ [9]:
ФЕ X = шах
1-
АР, X
^^баз X
100 + аг.
1 —
М, X
Л^баз. X
100
(6)
где а/ - веса (коэффициенты важности) показателей эффективности:
3
0 <а <1; = 1-
;=1
Этап 6. Экспертно или используя дополнительный количественный критерий, назначаемый исследователем (разработчиком) ПЗУ, выбирается окончательный вариант геометрии ПЗУ из множества Парето-оптимальных значений.
ПРИМЕР ВЫБОРА ОПТИМАЛЬНОЙ ГЕОМЕТРИИ ЦИКЛОННОГО ПЗУ
Принципиальная схема циклонного ПЗУ приведена на рис. 1.
Рис. 1. Принципиальная схема циклонного ПЗУ
Сепарация пыли в циклонах и их гидравлическое сопротивление зависят от основных конструктивных параметров, к которым относятся угол закрутки лопаточного аппарата а, диаметр корпуса циклона В и ряд линейных размеров, характеризующихся относительными величинами (отнесенными к диаметру В): внутренним диаметром центральной трубки ё, внешним диаметром центрального обтекателя В а и др.
При увеличении угла закрутки а возрастают тангенциальные скорости воздуха и твердых частиц. Вследствие этого увеличиваются центробежные силы и радиальные скорости частиц. Ввиду того, что центробежная сила пропорциональна V2 (V - скорость твердых частиц), степень очистки п находится в обратной зависимости от 1§2а. Однако, увеличение угла закрутки а приводит к увеличению турбулентности потока внутри циклона, что, в свою очередь, приводит к потерям давления за циклоном.
С увеличением диаметра корпуса В при постоянном расходе воздуха уменьшаются осевые скорости воздуха и твердых частиц, что ведет к пропорциональному увеличению продолжительности пребывания частиц в поле центробежных сил и способствует сепарации. Однако при этом уменьшается тангенциальная скорость воздуха, влияющая на силы инерции.
Уменьшение диаметра центральной трубки ё приводит к уменьшению степени очистки воздуха п, так как при этом меньшая часть траекторий твердых частиц заканчивается в пределах входного сечения трубки. Это приводит к значительному уменьшению поступающего в ГТД воздуха и потерям давления на выходе из циклона.
В процессе исследования эффективность ПЗУ оценивалась степенью очистки воздуха (п, %), потерей давления в ПЗУ (АР, Па) и массой ПЗУ (М, кг) [1-3].
Область определения параметров ПЗУ соответствовала следующим ограничениям:
- угол закрутки лопаточного аппарата 60° < а < 90°;
- внешний диаметр центрального обтекателя 6 < В а < 12 мм;
- внутренний диаметр центральной трубки 18 < ё < 24 мм;
- внутренний диаметр внешней трубки 16 < В < 30 мм.
Построение уравнений регрессии (п, АР, М) = / (х1, Х2, ..., хк) проводилось по данным центрального композиционного плана, включающего 25 экспериментов (опытов) (табл. 1)
[4].
Таблица 1
Матрица плана и результаты эксперимента
№ опыта хг х2 хз х4 Уг У2 Уз
(а, град.) (ВА, мм) (ё, мм) (В, мм) (П, %) (АР, Па) (М, кг)
1 -1 -1 1 1 91,5 581 0,0117
2 1 -1 -1 1 95,7 621 0,0148
3 -1 1 -1 1 89,5 535 0,0139
23 1 -1 -1 -1 93,0 609 0,0162
24 -1 1 -1 -1 97,9 653 0,0181
25 0 0 0 0 98,3 709 0,0201
В кодированных значениях переменных (х1, Х2, ..., хк)уравнения регрессии имеют вид:
П = 94,1-1,37х4 -0,56хх + 2,823х1х3 + 0,8хх +
2 2 2 + 0,8х2х3х - 2,8 8х2х4 + 3,7х1 - 3,4х2 - 0,8 7х3 ;
(погрешность оценки п по уравнению (7) равна 1,35 %);
АР = 629,84 -737,88х -772,68х2 -372,7х 3+1275,5^х 2+1441,28х Xx3 + +1252,68х: х4+ +1103,2 х2 х4 -1172,2х3 х4+1750,6x2;
(погрешность оценки АР по уравнению (8) равна 1,43 %);
M = 9856,68 +157,848х2 + 4434,316х 4 -30,283х 3 х 4 +645,6х 2
(8)
(9)
(погрешность оценки М по уравнению (9) равна 0,5 %).
В натуральных единицах измерения параметров а, Па, d и В уравнения регрессии имеют вид:
П = 104,7 -1,02 • В-1,035 • а • В А+1,89 • а • ё + 0,01-а • В + 56,98 • В А-й + + 0,91 • ВА•В-0,006• а2 + 76,27• В2Л-150,051 й2 %
ДР = 12230,3 - 20,32 • а -1681,15 • В А-26355,87 • ё + 21,18 • а • ВА +
+12,56• а• ё + 0,45• а• В + 36,44• ВА•В-41,1 ё• В + +15424,07• й2 Па;
М = 1861,5 + 3124,5ВА + 167,8В + 44,1 ё-В + 6,87 •В2 кг.
Результаты выбора параметров ПЗУ приведены в табл. 2, где под № 1 приведены данные исходного ПЗУ, а под (№ 2 ... № 5) 4 - Парето-оптимальных значений.
Оценка эффективности выбора оптимальной геометрии проводилась по формуле:
5 = Рисх - Р • 100 % Р
исх.
где Рисх- исходное значение показателя; Р - текущее значение показателя.
Парето-оптимальные значения
Таблица 2
№ п, % AP, Па М, кг д,
п/п 8 V ÖAP ¿м
1 94,5 630 0,0115 0 0 0
2 98,4 501 0,0138 +4,13 +20,4 -20,0
3 98,6 546 0,0132 +4,36 +13,2 -14,5
4 99,0 543 0,0159 +4,76 +13,7 -39,0
5 99,9 643 0,0173 +5,70 +2,2 -50,4
Примечание: «-» - показатель эффективности ухудшается; «+» - показатель эффективности улучшается.
В качестве окончательного было выбрано ПЗУ № 2, которое по сравнению с исходным ПЗУ:
- на 4,1 % имеет лучший показатель по параметру п;
- на 20 % имеет лучший показатель по параметру АР;
- но на 20 % уступает исходному ПЗУ по параметру M.
При этом оптимальные значения геометрии ПЗУ следующие: aopt = 90°; DAopt = 6 мм; dopt = 15 мм; Dopt = 25 мм.
Рис. 2. Влияние на степень очистки и потери давления ПЗУ:
а - угла закрутки а; б - внутреннего диаметра центральной трубки ё; в - внутреннего диаметра корпуса циклона D; г - внешнего диаметра центрального обтекателя DA
8 10
а б
Рис. 3. Влияние на массу ПЗУ:
а - внешнего диаметра центрального обтекателя DA; б - внутреннего диаметра корпуса циклона D
б
а
в
г
СПИСОК ЛИТЕРАТУРЫ
1. Полимерные композиционные материалы для изготовления пылезащитного устройства перспективного вертолетного двигателя / Г. Ф. Железина [и др.] // Авиационные материалы и технологии. 2018. № 1 (50). С. 58-63. [ G. F. Zhelezina, et al., "Polymer composite materials for the manufacture of a dustproof device for a promising helicopter engine", (in Russian), in Aviacionnye materialy i tehnologii, no. 1 (50), pp. 58-63, 2018. ]
2. Рысин Л. С. Еще раз о пыли // Журнал «Двигатель». 2009. № 2. [Электронный ресурс]. URL: http://engine.aviaport.ru/issues/65/page10.html (дата обращения 15.04.2021) [ L. S. Rysin (2021, Apr. 15), "One more time about dust" [Online], (in Russian), in zhurnal "Dvigatel", no. 2, 2009. Available: http://engine.aviaport.ru/issues/65/page10.html ]
3. Основы конструирования, производства и эксплуатации авиационных газотурбинных двигателей и энергетических установок в системе CALS технологий. Кн. 3: Эксплуатация и надежность ГТД и ЭУ / Н. Н. Сиротин [и др.]. М.: Наука, 2012. 616 с. [ N. N. Sirotin, et al., Fundamentals of design, production and operation of aircraft gas turbine engines and power plants in the CALS technology system. Book. 3: Operation and reliability of the gas turbine engine and power plant, (in Russian). Moscow: Nauka, 2012. ]
4. Гишваров А. С. Повышение эффективности многокритериального планирования многофакторного эксперимента. М: Машиностроение, 2014. 215 с. [ A. S. Gishvarov, Improving the efficiency of multi-criteria planning of a multi-factor experiment, (in Russian). Moscow: Mashinostroenie, 2014. ]
5. Волков К. Н., Емельянов В. Н., Курова И. В. Лагранжев подход к моделированию турбулентных течений газовзвеси в приложении к проблемам внутренней газодинамики РДТТ. Санкт-Петербург: «Военмех», 2008. С. 320-332. [ K. N. Volkov, V. N. Emel'ianov, I. V. Kurova, Lagrangian approach to modeling turbulent flows of a gas suspension as applied to the problems of internal gas dynamics of solid propellants, (in Russian). St. Petersburg: "Voenmeh", 2008. ]
6. Степанов Г. Ю., Зицер И. М. Инерционные воздухоочистители. М.: Машиностроение, 1986. 184 с. [ G. Y. Stepanov, I. M. Zicer, The inertial air cleaner, (in Russian). Moscow: Mashinostroenie, 1986. ]
7. Гишваров А. С., Аитов Р. Р., Айтумбетов А. М. Исследование эффективности пылезащитных устройств вертолетных газотурбинных двигателей // Вестник УГАТУ. 2015. Т. 19, № 2 (68). С. 100-110. [ A. S. Gishvarov, R. R. Aitov, A. M. Aytumbetov, "Modeling and optimization features dustproof device helicopter turboprop", (in Russian), in Vestnik UGATU, vol. 19, no. 2 (68), pp. 100-110, 2014. ]
8. Гишваров А. С., Салимзянова А. А., Рахимов А. Х. Исследование влияния геометрии вертолетного пылезащитного устройства на его эффективность // Мавлютовские чтения: материалы XI Всероссийский молодежной научной конференции УГАТУ. Уфа: УГАТУ, 2017. Т. 7. С. 62-67. [ A. S. Gishvarov, A. A. Salimzyanova, A. A. Rakhimov, "Investigation of the influence of the geometry of a helicopter dustproof device on its effectiveness", (in Russian), in Mavlyutov Readings: Materials of the XI All-Russian Youth Scientific Conference UGATU, 2017, vol. 7, pp. 62-67. ]
9. Ларичев О. И. Теория и метод принятия решений. М.: Логос, 2006. 392 с. [ O. I. Larichev, Theory and method of decision making, (in Russian). Moscow: Logos, 2006. ]
10. Мушкин Э., Мюллер П. Методы принятия технических решений. М.: Мир, 1990. 208 с. [ E. Mushkin, P. Miuller, Techniques for making technical decisions, (in Russian). Moscow: Mir, 1990. ]
11. Гишваров А. С. Совмещенные ресурсные испытания технических систем. Уфа: Гилем, 2002. 268 с. [ A. S. Gishvarov, Combined life tests of technical systems, (in Russian). Ufa: Gilem, 2002. ]
ОБ АВТОРАХ
ГИШВАРОВ Анас Саидович, проф., зав. каф. авиац. двиг. Дипл. инж.-мех. по авиац. двиг. (УАИ, 1973). Д-р техн. наук по тепл. двиг. летательных аппаратов (УГАТУ, 1993). Иссл. в обл. надежности, ресурса, испытаний и прогнозирования состояния техн. систем.
РАХИМОВ Абдусаттор Хасанович, асп. каф. авиационных двигателей. Дипл. маг. по авиастр. (УГАТУ, 2017). Иссл. в обл. надежности и ресурса авиац. двиг.
METADATA
Title: Study the effectiveness of the dust devices. Authors: A. S. Gishvarov1 A. H. Rahimov2
Affiliation: Ufa State Aviation Technical University (UGATU), Russia. Email: 1 kafedra.ad@mail, 2 mr.abdusattor@list.ru.ru Language: Russian.
Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), vol. 25, no. 4 (94), pp. 38-45, 2021.
ISSN 2225-2789 (Online), ISSN 1992-6502 (Print). Abstract: The present document addresses the selection method of optimal values (geometric dimensions) of the Engine Air Particle Separator (EAPS) parameters in aviation gas turbine engines (GTE). The methodology includes six main steps. The choice of these optimal values is carried out taking into account three particular criteria for the EAPS efficiency: the degree of air purification, pressure loss and mass. The EAPS parameters are selected using the Pareto method. Key words: particle separator; aircraft engine; method; optimal parameters. About authors:
GISHVAROV, Anas Saidovich, Prof., Dept. of Aircraft Engines. Dipl. engineer (USATU, 1973). Dr. of Tech. Sci. (USATU, 1993). RAHIMOV, Abdusattor Hasanovich, PhD Stud., Dept. of Aircraft Engines. Accident (USATU, 2017).