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

АВИАЦИОННАЯ И РАКЕТНО-КОСМИЧЕСКАЯ ТЕХНИКА

УДК 621.52

П. В. Кащеева

О ВЫБОРЕ ОПРЕДЕЛЯЮЩИХ ПАРАМЕТРОВ МНОГОСТУПЕНЧАТОЙ ТУРБИНЫ НИЗКОГО ДАВЛЕНИЯ ГТД

Рассмотрен аналитический метод выбора определяющих параметров, назначаемых на начальной стадии проектирования многоступенчатой неохлаждаемой газовой турбины. Метод основан на использовании количественных зависимостей оптимальных параметров, дающих необходимые в процессе проектирования представления о закономерностях изменения КПД от определяющих параметров в конкретных условиях проектируемого объекта. Газотурбинные двигатели; турбины низкого давления; оптимизация; выбор параметров

Выбор определяющих параметров турбины ГТД осуществляется на этапе одномерного газодинамического расчета и представляет собой начальную стадию многодисциплинарного мно-гоуровнего итерационного процесса проектирования. Для газодинамического расчета многоступенчатой турбины необходимо задаться рядом определяющих параметров, от правильного выбора которых зависят газодинамические потери в турбине и ее КПД. Необходимо определить оптимальное количество ступеней и распределить теплоперепады между ступенями, назначить степени реактивности, выбрать форму проточной части при соблюдении всех заданных ограничений. От эффективности решения этих задач зависит степень совершенства турбины, поэтому особенно важно не сделать на этом этапе ошибки.

На практике задача решается методом вариантных расчетов для конкретных условий проектируемого объекта, базирующемся на рекомендациях по выбору параметров и личном опыте проектировщика. В этом случае результат выбора становится зависимым от квалификации проектировщика и носит субъективный характер. К тому же при переборе ограниченного числа вариантов остается сомнение в том, что выбран действительно наилучший вариант решения. Автоматизация процесса проектирования и развитие методов оптимизации способствовали попытке внедрения процедуры оптимизации на всех этапах проектирования турбины. Использование метода численной оптимизации, работающей непосредственно с математической моделью, позволяет определить точные значения оптимальных параметров, обеспечивающие экстремум функции КПД турбины при задан-

Контактная информация: p.kascheeva@rambler.ru

ных ограничениях. Такой подход мог бы существенно повысить эффективность и сократить сроки проектирования, однако не получил широкого применения на практике вследствие трудностей, связанных с корректной постановкой задачи, определением ограничений и нереализуемыми полученными решениями.

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

Для достижения этой цели использовалась математическая модель одномерного газодинамического расчета многоступенчатой дозвуковой неохлаждаемой турбины, в которой исходные параметры представляются в безразмерном виде. Такое представление исходных данных позволило получить универсальные зависимости, применимые к турбинам с различными исходными данными. Достоверность такой методики во многом зависит от совпадения истинных и принятых в расчете потерь в проточной части. В этой связи проанализированы результаты выполненных разными авторами исследований профильных и вторичных потерь и потерь в радиальном зазоре. На основе результатов анализа и полученных статистических данных по геометрическим параметрам проточной части сделаны эмпирические и полуэмпирические обобщения, которые составляют основу универсальной методики расчета потерь. Математическая модель расчета также учитывает особенности преобразования энергии в диагональных ступенях турбины, обусловленные инерционными силами, порождаемыми кориолисовым

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

Алгоритм расчета турбины состоит из основных элементов, описанных в работе [1]. Задаются следующие исходные параметры: число ступеней г, степень понижения полного давления в турбине п , параметр газодинамической

нагруженности У = Л/^и,2 /сад (и7 - окружная

скорость на среднем диаметре 7-й ступени, сад -условная адиабатическая скорость многоступенчатой турбины), приведенная скорость и углы потока в абсолютном движении на входе и выходе из турбины - соответственно А0, а0 и Ат, а,., доли адиабатного теплоперепада, срабатываемого в ступенях Над 7 = Над 7 / Над т, степени реактивности ступеней р7, отношения среднего диаметра к высоте лопаток последней ступени (Б / И)т, отношение средних радиусов на выходе и входе в турбину гт / г0, относительные величины радиального зазора 5рз 7 = (5рз / Н)7, число Яе на входе в турбину, коэффициент расхода радиального зазора Црз и показатель адиабаты к.

Параметр относительного расхода газа для у-го венца:

— 1.П(1.) Iк + 1 1

О =-^. 8ш(а.) = ^.л1к—1Я-

.)

обобщает параметры состояния газа.

Раскрытие проточной части определяется изменением аксиальной площади:

F_

F

---- nt +1

Gn * 2 nt

GT

где п - показатель политропы процесса расширения газа в ступени находится методом последовательных приближений по формуле [2]:

ln

— = 1 + ■

( \

1 - 1 кт * к Рт V “т у * Лт

nt ln pT

Форма проточной части определяется из условия равномерного приращения аксиальной площади по венцам:

F

F„

= 1 +

F

Л

F

\ ax 0

1

j /2z

и изменением радиальности г7 + 1 / г7 в соответствии со среднестатистическим данными и всегда удовлетворяет условию ее плавности.

Используемые зависимости для расчета газодинамических потерь в решетках, описаны в работе [3]. Для определения мощностного КПД турбины по параметрам торможения с учетом возврата тепла использовалась формула, полученная в работе [2]:

( V1

Лт =

1 -

1

Лт

к-1 * к

1-

к-1 * к lT J

П

1-

1-

к-1 * к

Л

Мощностной коэффициент полезного действия турбины вычислялся по формуле [2]:

Лт

Л«т г

1 -«)к

1 - (Р )

Поиск оптимальных параметров Над 7, р7 при условии а = const (рис. 2, 3) осуществляется методом скользящего допуска [5], который является методом прямого поиска в задачах с ограничениями. При решении задач безусловной оптимизации этим методом используется метод деформируемого многогранника Нелдера-Мида. При этом одномерная модель дает возможность охватить широкие диапазоны изменения большого числа варьируемых параметров, что особенно важно на первых этапах оптимизации турбины.

Турбина низкого давления с тремя ступенями, согласно статистическим данным, - наиболее часто встречающийся случай многоступенчатой турбины в современной практике проектирования. Причем количество исследуемых параметров в трехступенчатой турбине невелико по сравнению с турбиной с большим количеством ступеней, что позволяет без особых сложностей анализировать их влияние на эффективность. Поэтому объектом расчетного исследования была выбрана трехступенчатая не-охлаждаемая турбина с умеренными числами Маха в проточной части, основные параметры которой типичны для ТНД ТРДД.

Выбор параметров многоступенчатой турбины лимитируется рядом ограничений. Весьма серьезным с точки зрения работы затурбинного устройства двигателя является ограничение по углу выхода потока из турбины ат. Ограничение по степени реактивности в корневых сечениях лопаток рвт 7 продиктовано желанием повысить

X

т

1

1

7=1

n=0

т

1-к

устойчивость потока против отрывных явлений в прикорневых сечениях. Учет этих взаимосвязанных требований налагает строгие ограничения на выбор параметров последней ступени. Поэтому целесообразно первоначально назначать Над ш и рш. Для выбора этих параметров предлагается использовать зависимости, показывающие их влияние на КПД при оптимальном распределении теплоперепадов и степеней реактивности в предшествующих ступенях. В качестве примера на рис. 1 приведены такие зависимости, полученные для турбины со следующими исходными данными:

П* = 4, У = 0,5, а0 I = 90°, А0 I = 0,3, ^

От = 0,448, (Б / Л)т = 5, гт/ г0 = 1,18. ( )

Исследуемый диапазон Над 111 ограничивался значениями скоростей на выходе из решеток СА и РК 0,5 < Аи (А^О < 0,9 (отклонение от ограничения обозначалось штрих-пунктирной линий), чтобы избежать некорректности в определении профильных потерь, а нижние значения степеней реактивности р7 достижением отрицательных значений у корня лопаток рвт 7 > 0 (отрицательные значений рвт 7 обозначено штриховой линией). Однако, как показывает практика, использование ступени с пониженной на среднем диаметре степенью реактивности оказывается возможным, если применить какой-либо из известных способов снижения радиального градиента реактивности [4].

Условие осевого выхода из турбины при умеренных значениях У не позволяет назначить оптимальные параметры III ступени, чтобы получить наибольший Пи т .

Как видно из рис. 1, а, б, чтобы обеспечить заданный угол с увеличением нагрузки Над ш требуется снижение рш. Очевидно, что меньший угол выхода потока из турбины аг позволяет назначить большую нагрузку на последнюю ступень Над ш, а также получить больший КПД Пи т . Однако это приводит к необходимости спрямления потока в спрямляющем аппарате (часто эту функцию выполняют профилированные силовые стойки за ТНД). Очевидно, что для определенного ат существует оптимальное сочетания параметров Над ш и pIII. Зависимости, показывающие влияние Над III и pIII на Дци т (а), Дп т (б) позволяют выбрать рациональные параметры нагрузки и реактивности последней ступени, обеспечивающие максимальный из возможных КПД турбины и удовлетворяющие заданным ограничениям по ат и рвт III. Кроме того, по ним можно оценить степень влияния этих параметров на КПД и ат.

К основным газодинамическим параметрам турбины относятся параметр нагруженности У и приведенная скорость на выходе из турбины вследствие их существенного влияния на КПД, массу и прочностные характеристики турбины. Оценить влияние этих параметров на КПД и возможность его увеличения позволяют зависимости Дци т = ^(У, Ат). КПД получен при оптимальных нагрузках и степенях реактивностей ступеней и ограничении по углу выхода потока из турбины. На рис. 2 для примера представле-

* л

ны такие зависимости, полученные для спт = 4 и ат = 90°. Выбор У и Ат согласуется с конструктивными и прочностными требованиями, предъявляемыми к турбине. По зависимостям, аналогичным показанным на рис. 2, определяются более точные значения оптимальных нагрузок и степеней реактивности ступеней.

Достаточно часто при проектировании и доводке турбины невозможно назначить оптимальные Над 7 и р7. Причиной тому может являться необходимость снижения температуры газа перед рабочими лопатками с целью улучшения условий работы ротора и сокращения числа охлаждаемых решеток, либо уменьшена глубина охлаждения и, как результат, снижены соответствующие потери. К тому же иногда требуется воспрепятствовать появлению сверхзвуковых скоростей на выходе из соплового аппарата, изменить осевую силу, воздействующую на ротор турбины или перепад давлений, под которым работает лабиринтное уплотнение между сопловым аппаратом и рабочим колесом. В этом случае необходимо знать влияние этих параметров на КПД турбины. В качестве примера на рис. 3 представлены зависимости, полученные для пт* = 4, У = 0,5, а0 I = 90°, А0 I = 0,3, ат = 90°, Ат = 0,5.

Рассмотрим применение представленного метода и возможные варианты увеличения КПД трехступенчатой турбины вентилятора ТРДД с исходными данными на проектном режиме, указанными выше (1), при следующих ограничениях: ат = 90°, Ат = 0,5. Исходное распределение нагрузок и реактивностей по ступеням следующее: Над_! = 0,43, Р! = 0,32, Над II = 0,3, рп = 0,35, Над III = 0,27, рш = 0,4. Расчетное значение КПД турбины ци т = 0,9329. Сопловой аппарат первой ступени выполнен охлаждаемым и его изменение требует значительных материальных затрат. Поэтому увеличение КПД достигается перераспределением нагрузок и реактивностей 2-й и 3-й ступеней с сохранением формы проточной части.

Р|ор1

0,44

0,4

0,36

0,32

0,6

0,2^ 0,5 0,4

0,3 * Л

*' * « ш . * *

0,2 0,25 0,3 0,35 0,4 0,45

д е

Рис. 1. Влияние величин НадШ и рш на КПД турбины (а, б) и оптимальные параметры I (в, д) и II (г, е) ступеней: ◊ - ат = 90°, ♦ - ат = 100°, ♦ - ат = 80°, ( - • - •-0,5 > Аи (А^) > 0,9,

---------- рШвт < 0)

а

в

а

в

д

б

г

е

Рис. 2. Влияние У и Ат на КПД турбины (а) и оптимальные параметры I (б, в), II (г) и III (д, е) ступеней

АПы*т

%

0

-0,5

-1

-1,5

0,5 0,4 0,3

0,6 / \

Pl К4 = 0.2 '

hu*max= 0,9483

0,25 0,3 0,35 0,4 Над!

б

в

г

Рис. 3. Влияние величин Над! и р! на КПД турбины (а) и оптимальные параметры II (б) и III (в, г) ступеней ( - ■ - ■-----------------0,5 > Хи О'^ц) > 0,9,------рПьт < 0)

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

Как следует из рис. 3, оптимальные параметры рассматриваемой турбины: Над п = 0,265, рп = 0,445, Над ш = 0,305, рш = 0,34. Результаты газодинамического расчета показывают, что можно получить более высокий (на Дци т = = 0,114 %) КПД, чем в исходном варианте, обеспечив на выходе осевое направление потока.

Из рис. 1 видно, что, уменьшая выходной угол турбины до 80°, можно выиграть Дци т » » 0,2 %. Согласно аналогичным, представленным на рис. 3, зависимостям для варианта с возможной закруткой потока на выходе из турбины (а = 80°) оптимальными параметрами будут Над п = 0,24, Pii = 0,45, Над ш = 0,33, рш = 0,42 и Дци т = 0,24 %. Рассмотрение зависимостей рис. 2 показывает возможное увеличение КПД за счет увеличения параметра на-груженности турбины Y и уменьшением выходной скорости.

Таким образом, последовательное применение основополагающих универсальных зависимостей, аналогичных описанным выше, образует аналитический метод выбора параметров назначаемых на начальной стадии проектирования многоступенчатой турбины. Используя данный подход, проектировщик имеет возможность не только определить оптимальные параметры, но и получает наглядное представление о степени влияния отклонения от них на КПД турбины. Применение метода позволяет повысить эффективность и сократить сроки проектирования газовой турбины ГТД.

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

1. Абианц В. Х. Теория авиационных газовых турбин. М.: Машиностроение, 1979. 246 с.

2. Богомолов Е. Н. Основы теории и выбор параметров авиационных газовых турбин. Ярославль, 1986. 88 с.

3. Богомолов Е. Н., Кащеева П. В. Оптимизация основных параметров, назначаемых при проектировании ступени осевой газовой турбины ГТД // Авиационно-космическая техника и технология. 2006. № 8(34). С. 129-132.

4. Мамаев Б. И., Сандимирова Т. А. Распределение теплоперепада по ступеням турбины вентилятора ТРДД. Лопаточные машины и струйные аппараты. Сборник статей. Вып. 8. Тр/ ЦИАМ № 1179, 1987. 140 с.

5. Химмельблау Д. Прикладное нелинейное программирование. М.: Мир, 1975. 536 с.

ОБ АВТОРЕ

Кащеева Полина Витальевна, старший преподаватель Рыбинск. гос. авиац. технологическ. акад. (РГАТА). Дипл. инж.-мех. (РГАТА, 2003). Канд. техн. наук по тепл. двигателям летательн. аппаратов (2009). Иссл. в обл. проектирования газовых турбин ГТД.