УДК 519.248+621.431.75
В01: 10.18287/2541-7533-2022-21-2-7-15
ОБОСНОВАНИЕ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ ПЫЛЕЗАЩИТНОГО УСТРОЙСТВА ГАЗОТУРБИННОГО ДВИГАТЕЛЯ
© 2022
доктор технических наук, профессор кафедры авиационных двигателей; Уфимский государственный авиационный технический университет; kafedra.ad@mail.ru
А. Х. Рахимов ст. преподаватель кафедры авиационных двигателей;
Уфимский государственный авиационный технический университет; mr.abdusattor@list.ru
Пылезащитное устройство считается важнейшим элементом для сохранения работоспособности вертолётных двигателей при эксплуатации в условиях запылённого воздуха. Интенсивность эрозии элементов газовоздушной проточной части газотурбинного двигателя зависит от многих факторов: от твёрдости и химического состава частиц пыли, от фракционного состава, концентрации песка, режима работы двигателя и др. Поэтому разработка методики выбора оптимальных значений параметров пылезащитного устройства является сложной многокритериальной и многофакторной задачей. В данной статье рассматривается методика выбора оптимальных значений геометрических размеров пылезащитного устройства авиационного газотурбинного двигателя. Методика включает шесть основных этапов. Выбор оптимальных значений геометрических параметров пылезащитного устройства проводится с учётом трёх частных критериев эффективности: степени очистки воздуха, потери давления в пылезащитном устройстве и массы на основе применения метода Парето. Рассматривается пример выбора оптимальных значений геометрических параметров циклонного пылезащитного устройства.
Пылезащитное устройство; показатели эффективности пылезащитного устройства; обобщённая функция цели; Парето-оптимальное множество
Цитирование: Гишваров А.С., Рахимов А.Х. Обоснование геометрических параметров пылезащитного устройства газотурбинного двигателя // Вестник Самарского университета. Аэрокосмическая техника, технологии и машиностроение. 2022. Т. 21, № 2. С. 7-15. DOI: 10.18287/2541-7533-2022-21-2-7-15
Пылевая эрозия элементов проточной части авиационных газотурбинных двигателей (ГТД) и отложения пыли вызывают снижение мощности, недопустимый рост температуры газа перед турбиной, поломки ослабленных эрозией лопаток, помпаж и другие отказы. Попавшие в двигатель частицы пыли не только изменяют форму профилей лопаток, но и оседают во входной части компрессора, что затрудняет надёжную работу двигателя вертолёта. В результате налипания пыли на лопатках компрессора их собственная частота колебаний может уменьшиться на 45%, что приводит к недопустимому увеличению действующих в лопатках напряжений и к их обрыву [1; 2].
Эрозионный износ лопатки компрессора приводит также к утонению передней и задней кромок, выработке поверхностей корыта, уменьшению значений хорды и максимальной толщины. Такие повреждения лопаток компрессора ведут к изменению режима их обтекания и, соответственно, к изменению параметров, определяющих условия работы профиля в потоке:
• угла между средней геометрической скоростью и хордой профиля;
• коэффициента подъёмной силы;
• угла набегания, определяемого как угол между вектором скорости воздуха перед решёткой и касательной к средней линии профиля у передней кромки;
• угла отставания потока;
• угла поворота потока в решётке.
Введение
Увеличение угла набегания за счёт эрозионного износа входной кромки может приводить к росту местных скоростей на поверхности профиля и к появлению сверхзвуковых скоростей потока.
Гидравлические потери при движении воздуха через канал решётки с небольшими скоростями до и после повреждения поверхностей профиля лопатки сводятся к трём видам:
1) профильные потери, связанные с образованием вторичных течений;
2) концевые потери, связанные с образованием пограничного слоя на стенках, ограничивающих канал решётки по высоте;
3) вторичные потери, связанные с возникновением протеканий в пограничном
слое.
Увеличение суммарных потерь полного давления, обусловленное ухудшением характеристик отдельных ступеней при эрозионном износе лопаток, приводит к падению напорности ступеней и уменьшению плотности воздуха по тракту компрессора. Происходит смещение характеристик в сторону меньших расходов, граница устойчивости смещается к рабочей линии. На границе помпажа при малых частотах вращения ротора работа первых ступеней компрессора неустойчива, а режим работы последних ступеней близок к режиму «запирания».
Увеличение удельного расхода топлива и температуры газов за турбиной происходит за счёт ухудшения аэродинамики проточной части вследствие увеличения радиальных зазоров между рабочими лопатками и корпусом компрессора, износа и шероховатости поверхности лопаток.
Снижение запаса газодинамической устойчивости из-за эрозионного износа на определённую величину может способствовать возникновению помпажа у двигателей. Одной из причин этого является разброс по характеристикам как новых, так и прошедших ремонт ГТД.
При эрозионном износе элементов конструкции двухконтурных двигателей с увеличением наработки также ухудшается их газодинамическая устойчивость. Наблюдается смещение линии рабочих режимов высокого давления, смещение напорных ветвей по расходу воздуха, изменение границы устойчивой работы. Например, изменение КПД компрессора или турбины каскада высокого давления двигателя типа Д-30КУ на 1% вызывает смещение линии рабочих режимов на 0,72%.
Поэтому на практике используются пылезащитные устройства (ПЗУ) авиационных двигателей, в которых для отделения пыли (песка) от воздуха используется их инерционность.
Целью данной работы является разработка методики выбора оптимальных значений геометрических параметров ПЗУ с учётом многокритериальности и многофакторности решаемой задачи.
Методика выбора параметров ПЗУ
Методика включает шесть основных этапов.
Этап 1. Формирование исходных данных. На данном этапе проводится:
- обоснование программного обеспечения, используемого для моделирования двухфазного потока «воздух - частицы песка». При этом рассматриваются Лагранже-вый, Эйлеровый и комбинированный методы [5];
- выбор режима работы двигателя;
- задание характеристик песка (количества, скорости, диаметра частиц и др.), области определения геометрических размеров ПЗУ и др.
Этап 2. Выбор частных параметров и критериев эффективности ПЗУ. На практике в качестве частных параметров эффективности рассматриваются два показателя [1; 2]:
- степень очистки воздуха - п, %;
- потеря давления в ПЗУ - Ap, Па.
В данной работе авторами вводится ещё один параметр эффективности - масса ПЗУ (М, кг). Таким образом, выбор оптимальных значений геометрических параметров ПЗУ Xopt = (, x2,..., xk ) проводится с учётом трёх критериев эффективности:
Ф (X ) = V ( X)= max п ( X);
Ф2 (X) = Apopt (X) = min Ар (X); (1)
Фз (X)= Mopt (X) = min Мг (X),
где X = (, x2,..., xk ) - вектор геометрических параметров ПЗУ.
Значения параметров п, Ар и М определяются моделированием двухфазного потока «воздух - частицы песка» с применением программного комплекса (ПК) Ansys Fluent для ПЗУ с различной геометрией. Область значений X = (xl,x2,...,xk) задаётся в
виде системы двухсторонних ограничений xmini £ xi < xmaxi; г = 1, k.
Этап 3. Нормирование частных параметров и критериев эффективности ПЗУ. При нормировании вместо «натурального» параметра вводится его отношение к некоторой нормирующей величине, измеряемой в тех же единицах, что и сам параметр [7 - 9]. В результате операции нормирования и приведения к единой шкале все частные критерии эффективности приобретают безразмерный вид с диапазоном измерения от 0 до 1:
Ф1 (X) = nOPt (X) = max п. (X)е 0...1;
Фн2 (X) = ApOHpt (X) = min DP (X)е 0.1; (2)
Фнъ (X)=Mp (X) = min Мг (X)е 0.1.
Критерии эффективности вида (2) можно преобразовать к виду, по которому эффективность ПЗУ будет тем выше, чем больше значение показателя (в %):
Ф1н (X) = max h (X) е (0... 100) %;
4ft (X)
Ф2н (X) = max
Фр (X) = max
Ap^ (X) Мг (X)
Мбаз (X)
100 е (0.100)%; (3)
100 е (0.100)%,
где Арбаз и Мбаз - базовые значения показателей Ар и М, соответствующие наиболее желательным значениям, к которым необходимо стремиться.
Представление критериев в виде (3) удобно для выбора оптимальных геометрических параметров ПЗУ в виде обобщённой целевой функции аддитивного вида [8; 9]:
ФЕ (X) = тах^ фН (X) =тах
/=1
V* (X) +
Ар. (X)
АРбаз. (X)
100 +
(X)
Мбаз. (X)
100
. (4)
Этап 4. Задаются двухсторонние ограничения на геометрические размеры ПЗУ, внутри интервала которых проводится
а.:
х^ ^ х Л ^ х^ ,
1тт 1 1тах'
ХК тш — ХК — ХК тах'
где х.тах, х.■ - максимальное и минимальное значения /-го фактора.
г тах ^ г тш а а
Внутри области Ох строится план эксперимента, в котором варьируемыми (независимыми) параметрами являются геометрические параметры ПЗУ (х1,х2,...,хк ), а выходными (зависимыми) параметрами являются значения п, Ар и М. По результатам эксперимента определяются уравнения регрессии, прошедшие проверку на значимость коэффициентов и адекватность уравнения:
П = /1 (( х2,..., хк ); АР = /2 (X1, х2,..., хк ) ; М = /з (X1, х2,..., хк ) .
Перед планированием эксперимента натуральные значения независимых факторов х1, х2,..., хк кодируются по формуле:
X = х/ х/ 0
Ах..
где X. - кодовое значение /-го фактора; х. - натуральное текущее значение /-го фактора; х.0 - начальный (нулевой) уровень фактора; Ах. - интервал варьирования /-го фактора:
Ах.. =
х — х
г тах г тт
Этап 5. Строится множество Парето-оптимальных значений функции Фх (X) = /\_Ф\ (X),Ф2Н (X),Ф3Н (X) с применением обобщённой целевой функции аддитивного вида (4) [10]:
Ф^ ) =
тах
^ПГ ^)+ а2
АР. (X)
Ар баз. (x)
100 + а3
Мг (X)
Мбаз. (X )
100
(6)
где аг - веса (коэффициенты важности) показателей эффективности:
2
0 £ а,. £ 1; Е а,. = 1.
,=1
Этап 6. Экспертно или используя дополнительный количественный критерий, назначаемый исследователем (разработчиком) ПЗУ, выбирается окончательный вариант геометрии ПЗУ из множества Парето-оптимальных значений.
Пример определения оптимальных значений геометрических параметров циклонного ПЗУ
Принципиальная схема циклонного ПЗУ приведена на рис. 1.
Рис. 1. Принципиальная схема циклонного пылезащитного устройства
Сепарация пыли в циклонах и их гидравлическое сопротивление зависят от основных конструктивных параметров, к которым относятся угол закрутки лопаточного аппарата а, диаметр корпуса циклона Б и ряд линейных размеров, характеризующихся относительными величинами (отнесёнными к диаметру Б): внутренним диаметром центральной трубки ё, внешним диаметром центрального обтекателя БА и др.
При увеличении угла закрутки а возрастают тангенциальные скорости воздуха и твёрдых частиц. Вследствие этого увеличиваются центробежные силы и радиальные скорости частиц. Ввиду того, что центробежная сила пропорциональна V2 (V - скорость твёрдых частиц), степень очистки п находится в обратной зависимости от tg2a . Однако увеличение угла закрутки а приводит к увеличению турбулентности потока внутри циклона, что, в свою очередь, приводит к потерям давления за циклоном.
С увеличением диаметра корпуса Б при постоянном расходе воздуха уменьшаются осевые скорости воздуха и твёрдых частиц, что ведёт к пропорциональному увеличению продолжительности пребывания частиц в поле центробежных сил и способствует сепарации. Однако при этом уменьшается тангенциальная скорость воздуха, влияющая на силы инерции.
Уменьшение диаметра центральной трубки ё приводит к уменьшению степени очистки воздуха п, так как при этом меньшая часть траекторий твёрдых частиц заканчивается в пределах входного сечения трубки. Это приводит к значительному уменьшению поступающего в ГТД воздуха и потерям давления на выходе из циклона.
В процессе исследования эффективность ПЗУ оценивалась степенью очистки воздуха (п, %), потерей давления в ПЗУ ( Ар, Па) и массой ПЗУ (М, кг) [2].
Область определения параметров ПЗУ соответствовала следующим ограничениям:
• угол закрутки лопаточного аппарата 60° < а < 90°;
• внешний диаметр центрального обтекателя 6 < БА < 12 мм;
• внутренний диаметр центральной трубки 18 < ё < 24 мм;
• внутренний диаметр внешней трубки 16 < Б < 30 мм.
Построение уравнений регрессий (п,Ар,М) = /(х1,х2,...,хк) проводилось по
данным центрального композиционного плана, включающего 25 экспериментов (опытов) [10].
В кодированных значениях переменных х1, х2,..., хк уравнения регрессии имеют
вид:
П = 94,1 — 1,37х4 -0,56х1х2 + 2,823х1х3 + 0,8х1х4 +
+ 0,8 х2 х3 х — 2,8 8 х2 х4 + 3,7 х1 — 3,4 х^ - 0,8 7 х^;
(погрешность оценки п по уравнению (7) равна 1,35%);
Ар = 629 — 737,88х1 — 772,68х2 — 372,7х 3 + +1275,5х 1 х 2+11441,28х 1 х 3 + 1252,68х 1 х 4 + +1103,2 х 2 х 4—1172,2 х 3 х 4+1750,6х 32;
(7)
(8)
(погрешность оценки Ар по уравнению (8) равна 1,43%);
М = 9856,68 + 157,848х2 + 4434,316х 4 —30,283х 3 х 4 +645,6х 2; (9)
(погрешность оценки М по уравнению (9) равна 0,5%).
В натуральных единицах измерения параметров а, ЮА, ё и Ю уравнения регрессии имеют вид:
П = 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,ЫЮ + 6,87Ю2, кг.
Результаты выбора параметров ПЗУ приведены в таблице, где под № 1 приведены результаты исходного ПЗУ, с которым проводилось сравнение ПЗУ с Парето-оптимальными значениями геометрических параметров (варианты № 2 ... № 5). Оценка эффективности выбора проводилась по формуле:
5 = 1р"сх ~ р 100,%,
рисх
где рисх - исходное значение показателя; р - текущее значение показателя.
Таблица. Значения Парето-оптимальных решений
№ Ар, Па М, кг S,
8л 8АР 8м
1 94,5 630 0,0115 0 0 0
2 98,4 501 0,0138 +4,13 +20,4 -20,4
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,13% имеет лучший показатель по параметру п, на 20% имеет лучший показатель по параметру Ар, но на 20% уступает исходному ПЗУ по параметру M. Оптимальные значения геометрических параметров ПЗУ следующие:
а opt =90^ DAopt =6 мм; dopt =15 мм; Dopt = 25 мм.
Библиографический список
1. Степанов Г.Ю., Зицер И.М. Инерционные воздухоочистители. М.: Машиностроение, 1986. 181 с.
2. Гишваров А.С., Аитов Р.Р., Айтумбетов А.М. Исследование эффективности пылезащитных устройств вертолётных газотурбинных двигателей // Вестник Уфимского государственного авиационного технического университета. 2015. Т. 19, № 2 (68). С. 100-110.
3. Рахимов А.Х., Салимзянова А.А. Исследование влияния геометрии вертолётного пылезащитного устройства на его эффективность // Материалы XI Всероссийской молодёжной научной конференции «Мавлютовские чтения». Т. 1. Уфа: УГАТУ, 2017. С. 62-67.
4. Железина Г.Ф., Соловьёва Н.А., Макрушин К.В., Рысин Л.С. Полимерные композиционные материалы для изготовления пылезащитного устройства перспективного вертолётного двигателя // Авиационные материалы и технологии. 2018. № 1 (50). С. 58-63. DOI: 10.18577/2071-9140-2018-0-1-58-63
5. Волков К.Н., Емельянов В.Н., Курова И.В. Лагранжев подход к моделированию турбулентных течений газовзвеси в приложении к проблемам внутренней газодинамики РДТТ // Сб. трудов шестой Всероссийской конференции «Внутрикамерные процессы и горение в установках на твёрдом топливе и ствольных системах» (8-10 сентября 2008 г., Санкт-Петербург, Россия). Ижевск: ИПМ УрО РАН, 2008. С. 320-332.
6. Рысин Л.С. Ещё раз о пыли // Двигатель. 2009. № 5 (65). http://engine.aviaport.ru/issues/65/page10.html
7. Ларичев О.И. Теория и метод принятия решений. М.: Логос, 2006. 392 с.
8. Мушкин Э., Мюллер П. Методы принятия технических решений. М.: Мир, 1990. 208 с.
9. Гишваров А.С. Совмещённые ресурсные испытания технических систем. Уфа: АН РБ, Гилем, 2002. 268 с.
10. Гишваров А.С., Рахимов А.Х. Планирование эксперимента при исследовании авиационных двигателей и энергетических установок: уч. пособие. Уфа: УГАТУ, 2019. 243 с.
JUSTIFICATION OF GEOMETRIC PARAMETERS OF THE ENGINE AIR PARTICLE SEPARATOR
© 2022
A S Gishvarov Doctor of Science (Engineering), Professor, Department of Aircraft Engines;
- Ufa State Aviation Technical University, Ufa, Russian Federation;
kafedra.ad@mail.ru
A. Kh. Rakhimov Senior Lecturer, Department of Aircraft Engines;
Ufa State Aviation Technical University, Ufa, Russian Federation; mr.abdusattor@list.ru
The Engine Air Particle Separator is considered to be the most important element for maintaining the performance of helicopter engines when operating in dusty air conditions. The intensity of erosion of the elements of the gas-air passage of the gas turbine engine depends on many factors: on the hardness and chemical composition of dust particles, on the fractional composition, sand concentration, engine operating condition etc.Therefore, the development of a method for selecting the optimal values of the Engine Air Particle Separator parameters is a complex multi-criteria and multi-factor task. The present document addresses the selection method of optimal values (geometric dimensions) of the Engine Air Particle Separator parameters in aviation gas turbine engines. The methodology includes six main steps. The choice of optimal values of the Engine Air Particle Separator geometrical parameters is accomplished taking into account three partial criteria for the Engine Air Particle Separator efficiency: the degree of air purification, pressure loss and mass. The Engine Air Particle Separator parameters are selected using the Pareto method. An example of choosing the optimal values of the geometric parameters of a cyclone Engine Air Particle Separator is considered.
Particle separator; optimal parameters of a dust protection device; generalized target function; Pare-to optimal set
Citation: Gishvarov A.S., Rakhimov A.Kh. Justification of geometric parameters of the engine air particle separator.
Vestnik of Samara University. Aerospace and Mechanical Engineering. 2022. V. 21, no. 2. P. 7-15.
DOI: 10.18287/2541-7533-2022-21-2-7-15
References
1. Stepanov G.Yu., Zitser I.M. Inertsionnye vozdukhoochistiteli [Inertial air cleaner]. Moscow: Mashinostroenie Publ., 1986. 181 p.
2. Gishvarov A.S., Aitov R.R., Aytumbetov A.M. Study the effectiveness of the dust devices. Vestnik UGATU. 2015. V. 19, no 2 (68). P. 100-110. (In Russ.)
3. Rakhimov A.Kh., Salimzyanova A.A. Issledovanie vliyaniya geometrii vertoletnogo pylezashchitnogo ustroystva na ego effektivnost'. Materialy XI Vserossiyskiy molodezhnoy nauchnoy konferentsii «Mavlyutovskie Chteniya». V. 1. Ufa: UGATU Publ., 2017. P. 62-67. (In Russ.)
4. Zhelezina G.F., Solovieva N.A., Makrushin K.V., Rysin L.S. Polymer composite materials for manufacturing engine air particle separation of advanced helicopter engine. Aviation Materials and Technologies. 2018. No. 1 (50). P. 58-63. (In Russ.). DOI: 10.18577/2071-9140-2018-0-1-58-63
5. Volkov K.N., Emel'yanov V.N., Kurova I.V. Lagranzhev podkhod k modelirovaniyu turbulentnykh techeniy gazovzvesi v prilozhenii k problemam vnutrenney gazodinamiki
RDTT. Sb. trudov shestoy Vserossiyskoy konferentsii «Vnutrikamernye Protsessy i Gorenie v Ustanovkakh na Tverdom Toplive i Stvol'nykh Sistemakh» (September, 8-10, 2018, Saint-Petersburg, Russia). Izhevsk: IPM UrO RAN Publ., 2008. P. 320-332. (In Russ.)
6. Rysin L.S. Once more about dust. Engine. 2009. No. 5 (65). (In Russ.). Available at: http://engine.aviaport.ru/issues/65/page10.html
7. Larichev O.I. Teoriya i metod prinyatiya resheniy [Theory and method of decision making]. Moscow: Logos Publ., 2006. 392 p.
8. Mushkin E., Myuller P. Metody prinyatiya tekhnicheskikh resheniy [Techniques for making technical decisions]. Moscow: Mir Publ., 1990. 208 p.
9. Gishvarov A.S. Sovmeshchennye resursnye ispytaniya tekhnicheskikh sistem [Combined life cycle tests of technical systems]. Ufa: AN RB, Gilem Publ., 2002. 268 p.
10. Gishvarov A.S., Rakhimov A.Kh. Planirovanie eksperimenta pri issledovanii aviatsionnykh dvigateley i energeticheskikh ustanovok: uch. posobie [Planning an experiment in the study of aircraft engines and power plants]. Ufa: UGATU Publ., 2019. 243 p.