Научная статья на тему 'Метод учета потерь в элементах проточной части газотурбинных двигателей'

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

CC BY
190
43
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ГАЗОТУРБИННЫЙ ДВИГАТЕЛЬ / ПРИВОД / ЭЛЕМЕНТ ПРОТОЧНОЙ ЧАСТИ / КОМПРЕССОР / ТУРБИНА / ГАЗОДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ / ПОТЕРИ / КОЭФФИЦИЕНТ ПОЛЕЗНОГО ДЕЙСТВИЯ / ПАРАМЕТРЫ / РЕЖИМЫ / ТЕХНИЧЕСКОЕ СОСТОЯНИЕ / ДИАГНОСТИКА / GAS TURBINE ENGINE / DRIVE / ELEMENT OF THE FLOW PART / NODES-COMPRESSOR / TURBINE / GAS DYNAMIC CHARACTERISTICS / LOSSES / EFFICIENCY / PARAMETERS / MODES / TECHNICAL CONDITION / DIAGNOSTICS

Аннотация научной статьи по механике и машиностроению, автор научной работы — Кривошеев Игорь Александрович, Рожков Кирилл Евгеньевич, Симонов Николай Борисович, Багаутдинов Марат Закиевич

Рассматриваются методы учета потерь в элементах проточной части (ПЧ) силовых установок (СУ) самолетов, газоперекачивающих агрегатов (ГПА), газотурбинных двигателей (ГТД) и приводов (ГПА) в их составе. Предложен новый метод аппроксимации газодинамических характеристик, расчета коэффициентов полезного действия (КПД) узлов, в т.ч. компрессоров и турбин оптимизации их геометрии и режимов. Для верификации метода использованы результаты заводских испытаний, 2D и 3D CAD/CAE-моделирования течения и газодинамических характеристик элементов ПЧ, включая продувки решеток профилей. Показано, что данный метод может эффективно использоваться для оптимизации геометрии ПЧ и выбора режимов работы входных и выходных устройств, лопаточных венцов в составе компрессоров и турбин, переходных каналов между каскадами. Предложено использовать метод при контроле изменения характеристик узлов ГТД в эксплуатации. Это позволяет существенно ускорить проектирование и доводку, повысить достоверность параметрической диагностики ГТД и ГТУ.

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

Похожие темы научных работ по механике и машиностроению , автор научной работы — Кривошеев Игорь Александрович, Рожков Кирилл Евгеньевич, Симонов Николай Борисович, Багаутдинов Марат Закиевич

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

Method for accounting losses in the elements of flowing part of gas turbine engines

The methods of accounting for losses in the elements of the flow part of the power plants (SU) of aircraft, gas pumping units (GPU), gas turbine engines (GTE) and drives (HPA) in their composition are considered. A new method of approximation of the gas-dynamic characteristics, the calculation of the efficiencies (efficiency) of the nodes, including compressors and turbines optimize their geometry and modes. For verification of the method, the results of factory tests, 2D and 3D CAD / CAE modeling of the flow and gas-dynamic characteristics of the elements of the inverter, including the purge of the gratings of profiles, were used. It is shown that this method can be effectively used to optimize the geometry of the inverter and select the operating modes of the input and output devices, blade crowns in the composition of compressors and turbines, transition channels between stages. It is proposed to use the method when monitoring changes in the characteristics of the CCD units in operation. This allows to significantly increase the reliability of parametric diagnostics of gas turbine engines and gas turbines. The proposed method can be used to monitor the state of the units and systems in the energy sector, in transport and in other areas.

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

ISSN 1992-6502 (Print)_

2019. Т. 23, № 1 (83). С. 40-48

Вестник УГАТУ

ISSN 2225-2789 (Online) http://journal.ugatu.ac.ru

УДК 621.452.3:621.48.01:681.518.54

Метод учета потерь в элементах проточной части

газотурбинных двигателей

и. а. Кривошеев1 , к. е. Рожков2 , н. б. Симонов3 , м. з. Багаутдинов4

1 krivosh777@mail.ru, 2 rke@mail.ru, 3 sventigo@yandex.ru,4 baganet@mail.ru ФГБОУ ВО «Уфимский государственный авиационный технический университет» (УГАТУ)

Поступила в редакцию 26.02.2019

Аннотация. Рассматриваются методы учета потерь в элементах проточной части (ПЧ) силовых установок (СУ) самолетов, газоперекачивающих агрегатов (ГПА), газотурбинных двигателей (ГТД) и приводов (ГПА) в их составе. Предложен новый метод аппроксимации газодинамических характеристик, расчета коэффициентов полезного действия (КПД) узлов, в т.ч. компрессоров и турбин оптимизации их геометрии и режимов. Для верификации метода использованы результаты заводских испытаний, 2D и 3D CAD/CAE-моделирования течения и газодинамических характеристик элементов ПЧ, включая продувки решеток профилей. Показано, что данный метод может эффективно использоваться для оптимизации геометрии ПЧ и выбора режимов работы входных и выходных устройств, лопаточных венцов в составе компрессоров и турбин, переходных каналов между каскадами. Предложено использовать метод при контроле изменения характеристик узлов ГТД в эксплуатации. Это позволяет существенно ускорить проектирование и доводку, повысить достоверность параметрической диагностики ГТД и ГТУ.

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

ВВЕДЕНИЕ

В настоящее время при выборе геометрии [1-5] и режимов работы элементов ПЧ СУ самолетов и ГПА, ГТД и ГТП в их составе, а также их узлов-компрессоров, турбин, переходных каналов между каскадами [6-11] производится условная оптимизация. При этом допустимые уровни КПД и потерь, запасов ГДУ (газодинамической устойчивости) и прочности (Кз) в компрессоре и турбине, требуемых пропускной способности (Свпр и Ат), степени повышения давления пк (в компрессорах) и мощности Ыт (в турбинах) рассматриваются как ограничения, а габариты (диаметры и длина) и массу - рассматривают как функции цели (при их минимизации). В ряде случаев - при проектировании лопаточных венцов (ЛВ) в составе турбин и

компрессоров (при выборе степени реактивности р, углов атаки /, распределения работы по ступеням), геометрии переходных каналов между каскадами компрессоров и турбин наоборот, минимизируются потери, максимизируется коэффициент восстановления полного давления и КПД.

При этом потери Ьг в МЛК (межлопаточных каналах) компрессоров учитывают с использованием коэффициента гидравлических потерь с.

2

В РК (рабочем колесе) с w

рк 2

-=V -

удельная работа (изменение статической энтальпии) в МЛК РК. Потери (на трение) в РК

Ар**

Pi

■. Поэтому

Ар

(* * \

Р„1 - Р*2)

'"рк

р1 *1 2

Р1 *1

Аналогично в НА (направляющем аппара-

2

те) $ * — = Ьна - удельная работа (изменена 2

ние статической энтальпии) в МЛК НА.

арн

Потери (на трение) в НА Ьг = Поэтому

$ на

* ^на

Ар

Г Н

Р 2

(* * \

Р2 - Рз /

2

Р с2

2

2

2

Р с2

Потери в компрессоре (рис. 1) используют при расчете КПД. Так КПД РК:

Прк = 1 -— = 1 -

= 1 —

¿рк

А/

рк 1 - £ *

1-

8т в

Ар

где Ь г =-

Грк Р1 т.д.) в РК.

рк

- потери энергии (на трение и

Аналогично КПД НА рассчитывается по формулам:

п = 1-

• ня

ь

ь.

ь

1__Гна

=1

£ на С22 2 - 2 С2 С3

А/..

н

1-

Бта

\2

2

ч Бта 3 у

Ар,

где Ьг = —— - потери энергии (на трение

нЯ Р2

и т.д.) в НА.

Соответственно КПД ступени определяется

как:

Пст =

(ЬркПрк + ЬнаПна ) = 1 (£*рк + £наС2 )

2 Н,,

= 1 - ПркР1-Пна(1-Р2).

Рис. 1. Диаграмма Хауэлла по продувке решеток профилей РК с отмеченными точками нулевых углов атаки, оптимального и номинального режимов, линиями коэффициента потерь и КПД

рк

$

2

2

на

В МЛК турбины используют коэффициенты скорости (в СА - сопловом аппарате у^/с^, в РК - рабочем колесе фРк=н2/н28), энергии £ (в СА £ СА=с 2 /с ^, в РК - рабочем колесе £РК=н 2/н \5). Кроме этого, для всех

элементов ПЧ используют коэффициент

* *

восстановления полного давления о=р,+1/р,. Здесь р - давление, с - скорость в абсолютном движении, н - скорость в относительном движении, * - заторможенное значение параметра, индексы: в компрессоре 1 - вход в РК, 2 - выход из РК и вход в НА, 3 - выход из НА; в турбине 1 - выход из СА, 2 -выход из РК. Индекс г - вход в МЛК, индекс г +1 - выход из МЛК. С использованием этих коэффициентов при оптимизации производится расчет КПД (п) узлов и элементов ПЧ. Однако начальные приближения при такой оптимизации до сих пор выбираются на основе недостаточно четких рекомендаций, с достаточно широким диапазоном параметров. Кроме того, для полноценной оптимизации требуются не значения этих коэффициентов «в точке» только для конкретного режима или комбинации параметров, а в виде газодинамической характеристики (ГДХ) элемента ПЧ.

ПРЕДЛАГАЕМЫЙ МЕТОД УЧЕТА

ПОТЕРЬ В ЭЛЕМЕНТАХ ПРОТОЧНОЙ ЧАСТИ ГТД

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

В данной работе для анализа потерь в элементах ПЧ и расчета КПД узлов предлагается использовать предложенный авторами метод отслеживания струй (МОС) [12]. Он основан на одновременном рассмотрении траекторий в относительном и абсолютном движении при движении частиц вдоль поверхностей тока в ПЧ СУ, ГТД и ГТП, узлов в их составе. Все элементы ПЧ в составе СУ самолета и ГПА, ГТД и ГТП, а также узлов в их составе работают в режиме как диффузоров (с торможением потока), так и конфузоров (с ускорением потока). При этом входное устройство (ВУ) самолета при работе ГТД на неподвижном самолете и в начале разгона при взлете, в конце торможения при посадке и при рулении на аэродроме работает в режиме конфузора. Во время полета самолета ВУ выполняет функцию диффу-

зора. Все остальные элементы ПЧ СУ самолетов и ГПА (входное и выходное устройства), ГТД и ГТП в их составе, а также их узлов - компрессоров, турбин, переходные каналы между каскадами на всех режимах выполняют функции либо диффузоров, либо конфузоров. Для диффу-зорного элемента ПЧ, в котором входное сечение обозначается индексом «г», коэффициент гидравлических потерь с записывается (в абсолютном движении)

_ 2 (р,;, - р)

(ре2),

Для конфузорного элемента ПЧ, в котором входное сечение обозначается индексом «г» коэффициент гидравлических потерь с записывается (в абсолютном движении)

2 (рм - Р*)

5= (рс )м .

Например, в абсолютном движении в ВНА компрессора или СА турбины

= 2 (р - р*)

5са(вна) / 2 \ ,

(рс А

в НА компрессора

=2 (р*- р.*) 5на / 2 \ .

(рС )2

Соответственно в относительном движении в РК компрессора в относительном движении

(* * \

р*1 - Р^

рк

(РН 2 )1 '

в РК турбины

(* * 1

рн1 - рн2 /

(РН 2 )2 '

Для таких элементов ПЧ переход к коэффициенту восстановления полного давления о производится по формулам: в абсолютном движении для диффузора

о = 1 - ^ Я2 к+1 * 1 - 5 [1 - п ( ^ )];

для конфузора

Ь *1 - 5 - п(^+■)] •

о = 1 - +1 )^2+1

Например, в абсолютном движении для ВНА компрессора и СА турбины

к

0са(вна) =1 ?са(вна)£(^1)^1 * 1 - Сса(вна) [1 - П (X )] ,

для НА компрессоров

°на = 1 - Сня^Х^ —

к +1

к +1

а - [1 - п (^2)].

В относительном движении для РК компрессора

ь

о =1 - С *я

* к + 1 * 1 - С[1 - П (^* )] ,

для РК турбины

ь

о = 1 - С б(Х

к +1

1 - С^ [1 - п (X*2)].

Для построения ГДХ таких элементов ПЧ предлагается использовать аппроксимацию вида: в абсолютном движении для диффузора

о = 1 - 8(КЛ /) [кл / ]2 для конфузора

к

к +1

1 - п (кл /),

о = 1 - в(КЛ/+,) [КЛг+, ]2

к +1

: п

( КЛ г+,).

Например, в режиме конфузора работает входное устройство (ВУ) ГПА и ГТУ в энергетике, а также входное устройство самолета перед взлетом. При расчете потерь в ПЧ ГТД эти зависимости также удобно использовать. Например, в абсолютном движении для ВНА компрессора и СА турбины

Оса(вна) = 1 - 8^) [^ ]2

для НА компрессора Она = 1 - 8(КЛ) []2

к +1

п ( ^Х),

к +1

: п

() .

В относительном движении для РК компрессора

*рк V х )[к хх *1]

рк

* мк хХ * ),

к+1

для РК турбины

V X *2 )КXх*2 ]2

к

к + 1

"П(КXх*2 ).

Здесь К^Х - «кажущаяся» приведенная скорость в элементе ПЧ. Очевидно, что чем меньше значение Кх (0< Кх<1), тем лучше ГДХ элемента ПЧ (меньше потери). Значение коэффициента Кх определяется по значению с и приведенной скорости X (на входе в диффузор или на выходе из диффузора, соответственно, в абсолютном движении для элементов статора и в относительном движении для элементов ротора). Для этого «в точке» решается уравнение

1 - Св ( X /) X = 1 - 8 (КхХ/ )[КхХ/ ]2 '

к +1

к + 1

или к

1 - Св (X/) х2 — = п (КхХ/) .

Если по аналогии с обратными тригонометрическими функциями (типа Arcsin...) ввести понятие обратной газодинамической функции агсп(...), то

К, =— агсп

х X,.

1 - Св (X/) X!

к +1

Очевидно, что для МЛК РК и НА показатели К зависят от степени диффузорно-сти В (например, от фактора диффузорно-сти ¥В по Либляйну) и от угла атаки /. Обобщенные зависимости ^(В,/) получаются путем обработки результатов продувок решеток профилей [1-3].

АНАЛИЗ ПОТЕРЬ В ЭЛЕМЕНТАХ ПРОТОЧНОЙ ЧАСТИ ГТД С ИСПОЛЬЗОВАНИЕМ МЕТОДА ОТСЛЕЖИВАНИЯ СТРУЙ

Предложенный авторами метод отслеживания струй (МОС) предполагает одновременное рассмотрение движения частиц в ПЧ компрессоров и турбин вдоль траекторий в относительном и абсолютном движении. При этом изменение вдоль ПЧ заторможенной температуры в относительном

гт*

движении Т * рассматривается не только в рабочих колесах, но и для лопаточных вен-

к

2

цов статора, оно определяется коэффициентом теоретического напора в виртуальном вращении, которое, в свою очередь, определяется геометрией лопаточных венцов (ЛВ) НА и ВНА в компрессоре, СА в турбине. Полагается, что в виртуальном относительном движении статор вращается с той же частотой, что и ротор, но в обратном направлении, в каждом сечении виртуальная окружная скорость статора ин=-и. Рассматриваются траектории движения вдоль поверхностей тока вдоль втулки ротора (сечения «вт»), вдоль внутренней поверхности корпуса (сечения «к») и вдоль промежуточных поверхностей тока, например, эффективной поверхности «эф», для которой на входе в 1 ступень

Бф = Б

эф к

У1 + й7

2

Вдоль каждой поверхности тока для компрессора о=По,. Здесь о, - показатели для отдельных ступеней. Для каждой ,-ступени

о =

0вна \/0оз10оз1» 0ркн

4

0оз20оз2 н0на

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

Здесь 0вна , 0оз1 , 0оз2 , 0на - показатели в

абсолютном движении для элементов ПЧ статора (МЛК ВНА, осевых зазоров между

ВНА и РК, РК и НА, НА), о

оз1н ■

о

оз2н '

оРКн - показатели в относительном движении для элементов ПЧ ротора (МЛК РК, осевых зазоров между ВНА и РК, РК и НА). При этом ооз1 и ооз2 определяются в

абсолютном движении по и о

оз1н

и

ооз2н определяются в относительном движении по и Аналогично в турбине для каждой ступени

о =

.4

о о о

оз оз н ркн

Здесь оса, ооз - показатели в абсолютном движении для элементов ПЧ статора (МЛК СА, осевого зазоров между СА и РК), оозн, оркн - показатели в относительном

движении для элементов ПЧ ротора (МЛК РК, осевого зазора между СА и РК). При

этом ооз определяется в абсолютном движении по ^1. оозн определяется в относительном движении по

Традиционно при определении КПД в компрессоре и турбине рассматривается идеализированное (виртуальное) значение изоэнтропической температуры за РК. За РК компрессора

к-1 к

Т* = Т (прк)

за НА ступени

т* = т; (<т)"

для многоступенчатого компрессора (МК)

к-1 к

77* гт1*1 * \

к; = Т (пк )

к-1 к

В турбине за СА

Т * = Т0

■45 к г -1

(П™ т )

кг

за РК ступени

Т * =■

Т *

Т

п

к г -1

(П * ) кг (П * )

\ рк т / \ ст т /

кг -1

для многоступенчатой турбины (МТ)

Т * = То

1TS кг-1 •

( п* У

В отличие от этого, авторами предложено по фактическому изменению затормо-

*

женной температуры Т в ЛВ, ступени и многоступенчатого компрессора или турбины определять идеализированное изменение

*

заторможенного давления Р ид. При этом в РК компрессора

(т*Л

к-1

ПРКид

ТМ к

2 Т *

К11 У

В НА

п„

С Т * V

Т

К1 2

к-1

= 1.

В ступени компрессора

к-1

п,

СТид

т

2 Т *

V у

= п

РКид'

В МК

ГТ

к-1

П^ =

Кид

т

V11 у

Тогда фактическое значение степени повышения полного давления в РК определяется как

прк пРКид°.

Если перед РК находится ВНА, то следует определять степень повышения давления в системе ВНА-РК, для чего умножать

п

'РКид

на произведение Овнал10оз1Ооз1* 0рк* .

В общем случае степень повышения давле-[ в стуг Здесь

ния в ступени компрессора пст = пстидост.

0ст 0вн^Т°оз10оз1* 0рк^0оз20оз2* 0на .

В МК (многоступенчатом компрессоре)

* *

П К =П Кид а К •

Здесь

ат

=п

а

Аналогично в турбине: в СА

п,

САид

В РК турбины

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

гт *т 10

т

V у

кг -1

к- -1

= 1.

п

РКид

С Т * V кг Т *

V12 * у

В ступени турбины

п,

СТид

* V" к

10

Т '*

V12 у

кг-1

С Т * V1Г

10 Т *

V12 у

кг -1

С учетом потерь полного давления фактические отношения давления: в СА

пса пСАид0са 0са ,

в РК

п

прк =■

'РКид

РК*

в ступени турбины

п,

пСТт

СТид

в многоступенчатой турбине

п* = ^

Здесь

0СТт 0с^0ОЗ0ОЗ* 0РК* .

В многоступенчатой турбине

о =П о ..

т X X ст г

В компрессоре

гр* Гр* кГ -1

= ТК о кг

Т* ~ т К ' Т1 Т1

Аналогично в ступени

В РК

гр* гр* кГ -1

= ^^ о кг гр* гр* ст *

Т1 Т1

кг-1

= ^^ о кг

гр* ~ гр* ° РК* • Т1 Т1

С учетом этого КПД вдоль выбранной поверхности тока: для компрессора

Пк =1-

к-1

1-а

Г гр* N

с Т

+ 1

Н„

V тк у

для ступени

* 1

Пст = 1 -

С к-1V С Т* V

^+1 Н

V Тст у

1 - о к

ст

V у

= 1 -

С к-1 V 1 - о к

ст

V у

гр*

Срт1и

2 Л

у—+1 Н

V Тст у

для РК

к

*

т

ст г

*

к

К

* 1

Прк =1 -

С к-1V С т* V

^+1 н

V Тст у

1 - о к

рк *

у

^ к-^ С С т%,2 V

= 1 -

1 - о к

рк *

7-г* 2

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

р_ 1й +1

Н.

Тст

у

для НА

С к-1 V

* 1

Пна = 1 -

1 - о к

на

V у

к-1 С

= 1 -При этом

1 - о к

на

V у

т*

С рТ2* + 1 Н

V т* у 2 V

с Т„ и

р 2* + 1 Н

V т*

у

Н т = 1-

С 2а

где

1

С2а

С3а

( ^ Р 2 + 2 )' 1

(^2 + С^3 )

Нт = (с^а2 - С^а3 )Саи . Аналогично в турбине

гр* гр* кг-1

= о кг Т0 Т0

В ступени турбины

гр* гр* кг -1 Т25 = ^^ о кг т* ~ т* СТт '

В СА турбины

гр* гр* кг-1

= о кг

тг* Т7* Са *

Т0 Т0

В РК турбины

гр* кг-1 Т25 = о кг

гр* ~ гр* РК» • Т0 Т0

С учетом этого КПД вдоль выбранной поверхности тока: для турбины

пТ =1/

1+

С кг- V 1- окг

гр* СрТ0

ь,,

-1

для ступени

С = 1/

для РК

1+

С кг-1 VС т* V

1-о кг

ст

V у

р 0 -1 ь

V и

у

Прк = 1/ для СА

Пса = 1/

1 +

( кг-Л^С Т* ^

1-о

РКw

"У 0

ь

V и

-1

1+

С кг-1 V 1-о кг

са

V у

гр* СрТ0 *

ь

V и*

-1

Как для компрессора, так и для турбины внутренний изоэнтропический КПД определяется с учетом концевых потерь (потерь в радиальном зазоре и потерь на трение дисков, трения потока рабочего тела о поверхности корпуса и втулки ротора): в * * * — —

компрессоре П К =Л Ки ■ Пконц = П Ки ■ Врад •П/, в турбине П Т = Пт-и ' Пконц = •П *Ти ■ Прад ' П/ . Методы учета радиальных потерь рассмотрены в ряде работ.

Приведенные зависимости позволяют строить линии КПД на характеристиках компрессоров и турбин, ступеней и ЛВ в их составе, обоснованно выбирать значения степени реактивности р, углов атаки /, густоты решеток профилей ¿Л и т.д. Анализ по предложенной методике результатов 3Б СЛО/СЛЕ-моделирования, испытаний компрессора позволяет организовать и ускорить доводку компрессоров и турбин. Например, требуется обеспечить постоянство значений ог для всех ЛВ РК и НА в составе ступеней МК (многоступенчатого компрессора). Тогда параметр о=Пог для всего МК, т.е. произведение о/ всех ЛВ и ОЗ (осевых зазоров) максимально, и КПД компрессора максимален. Это позволяет оптимально распределять работу по ступеням и по высоте ПЧ. Из этого же следует, что для промежуточных ступеней компрессора следует выбирать степень реактивности р не р-0,5, а даже меньше чтобы обеспечить условие X*!- XС2.

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

Предложенный метод позволяет получить ГДХ ВУ самолета при работе двигателя на стоянке (при нулевой скорости движения). Аналогично может быть получена характеристика ВУ ГПА. Обычно при нулевой наработке изначально по результатам ПСИ (приемосдаточных испытаний) для ГТД со стендовым (с лемнискатным) воздухозаборником (у которого предполагается о=1) определяется ДрХ (дроссельные характеристики), среди которых Свпр(п1пр) и q(X1в)=/(n1пр). После установки ГТД в СУ самолета или на ГПА может быть определено значение оВУ и построена ГДХ ВУ самолета. Для этого ГТД выводится на режим с некоторым значением п1пр и измеряется значение рв перед входом в двигатель (перед ВНА или РК 1 ступени). Одновременно по ДрХ ГТД (т.к. наработка на новом ГТД все еще близка к нулю т~0) определяется значение ^ на входе. Значение ^(п1пр) должно быть такое же, как при ПСИ. С учетом того, что

(рв/РН )= п(^1в )аВУ * п(^1в МКАв У

то существует возможность определить значение К%. Для упрощения решения задачи получения ГДХ ВУ предложена аппроксимация

(рвР 0>5 МКр«Ч,/2у+1|

В этом случае подбирается значение Кр. Тогда на последующих этапах эксплуатации появляется возможность определять значения Х1в = агссо8 [2(рв/ри) -1] • 2 / Кр путем замера статического давления рв перед входом в двигатель. При этом определяется исходное значение показателя Х1в/п1пр и сравнивается с исходным значением, полученным при ПСИ. Кроме того, определяется теку-

щее значение комплекса

i *

К ) k -1

и

сравнивается также с исходным значением, полученным при ПСИ. В данном случае предполагается, что по мере увеличения наработки изменяются (деградируют) характеристики ГТД, компрессоров и турбин

в его составе, но характеристика ВУ изменяется мало. Это позволяет каждый раз, выводя двигатель на какой-либо режим, измерять п1пр и рв. Решая уравнение

(рв/ри МЧ, )я(К^ У или более простое

(рвР 0>5 МКрП^1в/2)+l], имеется возможность определить значение Х1в и значение Х1в/п1пр. Уменьшение значения этого показателя является признаком загрязнения ПЧ. Уменьшение значения показателя

v. \м " ) k -1

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

а:

является признаком эрозии лопаток компрессора. Аналогичным образом в анализ может быть включен параметр - скольжение роторов п2/п1.

ЗАКЛЮЧЕНИЕ

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

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

1. Бунимович А. И., Святогоров A. A. Аэродинамические характеристики плоских компрессорных решеток при большой дозвуковой скорости // Сборник статей «Лопаточные машины и струйные аппараты». М.: Машиностроение, 1967. Вып. 2. С. 5-35. [ A. I. Bunimovich, A. A. Svyatogo-rov, Aerodynamic characteristics of flat compressor grids at high subsonic speed, (in Russian), in "Lopatochnie mashini I stryinie apparati". Moscow: Mashinostroenie, 1967. ]

2. Emery J. C., Herring J. I. and Erwin J. R., "Systematic twodimensional cascade test of NACA-65 serias compressor blades at low speeds," NACA RML51G, 31, 1951. [ J. C. Emery, J. I. Herring and J. R. Erwin, "Systematic twodimensional cascade test of NACA-65 serias compressor blades at low speeds," in NACA RML51G, 1951. ]

3. Дейч М. Е., Самойлович Г. С. Основы аэродинамики осевых турбомашин / М. Е. Дейч, Г. С. Самойлович. М.: Машгиз, 1959. 428 с. [ M. E. Deutsch, G. S. Samoylovich, Fundamentals of aerodynamics of axial turbomachines, (in Russian). M.: Mashgiz, 1959. ]

4. Федоров Р. М. Характеристики осевых компрессоров. Воронеж, Изд-во Научная книга, 2015. 220 с. [ R. M. Fedorov, Characteristics of axial compressors, (in Russian). Voronezh: Publishing House Scientific book, 2015. ]

5. Белоусов А. Н., Мусаткин Н. Ф., Радько В. М. Теория и расчет авиационных лопаточных машин / А. Н. Белоусов, Н. Ф. Мусаткин, В. М. Радько. Самара: Сам. Дом печати, 2003. 344 с. [ A. N. Belousov, N. F. Musatkin, V. M. Radko, Theory and Calculation of Aircraft Blade Machines, (in Russian). Samara: Sam. House of Printing, 2003. ]

6. Чичков Б. А. Модели и параметрическая диагностика авиационных двигателей: учебное пособие. М.: МГТУ ГА, Часть 1. 2004. 96 c. Часть 2. 2005. 60 c. [ B. A. Chichkov, Models and parametric diagnostics of aircraft engines, (in Russian). Moscow: MGTU GA. ]

7. Технология эксплуатации, диагностики и ремонта газотурбинных двигателей: учебное пособие / Ю. С. Елисеев и др. М.: Высшая школа. 2002. 355 с. [ Y. S. Eliseev, et. al., Technology of operation, diagnostics and repair of gas turbine engines, (in Russian). Moscow: Higher School, 2002. ]

8. Оценка технического состояния авиационных ГТД / Лозицкий Л. П. и др. М.: Транспорт, 1982. 167 с. [ L. P. Lozitsky, et. al., Assessment of technical condition of aviation GTE, (in Russian). Moscow: Transport, 1982. ]

9. Практическая диагностика авиационных газотурбинных двигателей / под редакцией В. Л. Степаненко. М.: Транспорт, 1985. 116 c. [ V. L. Stepanenko, Practical diagnostics of aviation gas turbine engines, (in Russian). Moscow: Transport, 1985. ]

10. Шереметьев А. В. Анализ опыта эксплуатации зарубежных газотурбинных двигателей по техническому состоянию / А. В. Шереметьев // Авиационно-космическая техника и технология. Х.: ХАИ. 2003. Вып. 40/5. С. 5-9. [ A. V. Sheremetev, "Analysis of operational experience of foreign gas turbine engines according to the technical state", (in Russian), in Aviachionno-kosmicheskai technika I technology, HAI, issue 40/5, pp. 5-9, 2003. ]

11. Тарасенко А. В., Кондратьев Ю. В. Оценка технического состояния устройств компенсации реактивной мощности в тяговой сети // Успехи современной науки и образования. 2016. №10. Т. 5. C. 155-157. [ A. V. Tarasenko, Yu. V. Kondratiev, "Evaluation of the technical condition of reactive power compensation devices in the traction network", (in Russian), in Progress in Modern Science and Education, vol. 5, no. 10, pp. 155-157, 2016. ]

12. Кривошеев И. А., Рожков К. Е., Симонов Н. Б. Основы метода отслеживания струй // Вестник УГАТУ. 2015. Т. 19, № 3 (69). С. 103-110. [ I. A. Krivosheev, K. E. Rozhkov, N. B. Simonov, "Streamline tracking method basics", (in Russian), in Vestnik UGATU, vol. 19, no. 3 (69), pp. 103-110, 2015. ]

ОБ АВТОРАХ

КРИВОШЕЕВ Игорь Александрович, д-р техн. наук, проф.

каф. АД УГАТУ.

РОЖКОВ Кирилл Евгеньевич, канд. техн. наук, доц. каф.

АТиТ.

СИМОНОВ Николай Борисович, асс. каф. АД УГАТУ.

БАГАУТДИНОВ Марат Закиевич, асп. каф. АД УГАТУ.

METADATA

Title: Method for accounting losses in the elements of flowing part of gas turbine engines

Authors: I. A. Krivosheev1, K. E. Rojkov2, N. B. Simonov3, M.Z. Bagautdinov.

Affiliation: Ufa State Aviation Technical University (UGATU), Russia.

Email: 1 krivosh777@mail.ru, 2 rke@mail.ru, 3 sventigo@yandex.ru, 4 baganet@mail.ru/

Language: Russian.

Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), vol. 23, no. 1 (83), pp. 40-48, 2019. ISSN 2225-2789 (Online), ISSN 1992-6502 (Print).

Abstract: The methods of accounting for losses in the elements of the flow part of the power plants (SU) of aircraft, gas pumping units (GPU), gas turbine engines (GTE) and drives (HPA) in their composition are considered. A new method of approximation of the gas-dynamic characteristics, the calculation of the efficiencies (efficiency) of the nodes, including compressors and turbines optimize their geometry and modes. For verification of the method, the results of factory tests, 2D and 3D CAD / CAE modeling of the flow and gas-dynamic characteristics of the elements of the inverter, including the purge of the gratings of profiles, were used. It is shown that this method can be effectively used to optimize the geometry of the inverter and select the operating modes of the input and output devices, blade crowns in the composition of compressors and turbines, transition channels between stages. It is proposed to use the method when monitoring changes in the characteristics of the CCD units in operation. This allows to significantly increase the reliability of parametric diagnostics of gas turbine engines and gas turbines. The proposed method can be used to monitor the state of the units and systems in the energy sector, in transport and in other areas.

Key words: Gas turbine engine; drive; element of the flow part; nodes-compressor; turbine; gas dynamic characteristics; losses; efficiency; parameters; modes; technical condition; diagnostics.

About authors:

KRIVOSHEEV, Igor Aleksandrovich, prof. Aviation Engines Dept. Dipl. engineer (USATU 1976). Dr. Tech. Sciences (USATU, 2000).

ROJKOV, Kirill Evgenevich, cand. of Tech.Sci. Associate Professor of the Dept. of Aviation Heat and Power Engineering.

SIMONOV, Nikolai Borisovich, ass. of the Dept. of Aircraft Engine.

BAGAUTDINOV, Marat Zakievich, postgraduate student of the Dept. of Aircraft Engine.

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