CC BY
56
DOI: 10.18721/^.230406. УДК 621.165.001.2: 621.438.001.2
А. К. Нгуен, К.Л. Лапшин
Санкт-Петербургский политехнический университет Петра Великого,
Санкт-Петербург, Россия
ВЛИЯНИЕ ТАНГЕНЦИАЛЬНОГО НАКЛОНА РАБОЧИХ ЛОПАТОК НА ПОТЕРИ КИНЕТИЧЕСКОЙ ЭНЕРГИИ В ОСЕВОЙ ТУРБИННОЙ СТУПЕНИ
Исследовано влияние тангенциального наклона рабочих лопаток на характеристики, структуру потока и потери кинетической энергии в осевой турбинной ступени. Выбраны расчетные модели ступеней с тангенциальным наклоном рабочих лопаток и для них выполнены 3D газодинамические расчеты в программном комплексе ANSYS CFX. Представлены графики изменения КПД, степени реактивности и коэффициента суммарных потерь кинетической энергии на переменных по частоте вращения ротора режимах, распределения параметров потока (выходная скорость из ступени и коэффициенты потерь кинетической энергии в направляющем аппарате и рабочем колесе) по высоте лопатки. Выполнено сравнение результатов расчета для ступеней с различным тангенциальным наклоном рабочих лопаток. Дано заключение о влиянии тангенциального наклона рабочих лопаток на характеристики, структуру потока и потери кинетической энергии и о его рациональном применении в осевых турбинных ступенях.
ТУРБИННАЯ СТУПЕНЬ; ТАНГЕНЦИАЛЬНЫЙ НАКЛОН РАБОЧИХ ЛОПАТОК; ПОТЕРИ КИНЕТИЧЕСКОЙ ЭНЕРГИИ; СТЕПЕНЬ РЕАКТИВНОСТИ; АШУ8 СБХ.
Ссылка при цитировании:
А. К. Нгуен, К. Л. Лапшин. Влияние тангенциального наклона рабочих лопаток на потери кинетической энергии в осевой турбинной ступени // Научно-технические ведомости СПбПУ. Естественные и инженерные науки. 2017. Т. 23. № 4. С. 66-73. DOI: 10.10871ДЕБТ.230406.
A. Q. Nguyen, K. L. Lapshin
Peter the Great St. Petersburg polytechnic university, St.-Petersburg, Russia
THE INFLUENCE OF TANGENTIAL INCLINATION OF ROTOR BLADES ON KINETIC ENERGY LOSSES IN AN AXIAL TURBINE STAGE
The influence of tangential inclination of rotor blades on the characteristics, flow structure and kinetic energy losses in an axial turbine stage were explored. Calculation models of the stages with tangential inclination of the rotor blades were constructed and 3D-testing gas-dynamic calculations were performed in ANSYS CFX. The graphs of the efficiency index, the degree of reactivity and the coefficient of the total kinetic energy losses at variable rotor frequency were presented. The distributions of the flow parameters (output speed from the turbine stage and the coefficient of kinetic energy losses for the guide vanes and the impellers) were given by the height of the blade. The calculation results between turbine stages with different tangential inclination of the rotor blades were compared. The conclusion about the influences of tangential inclination of the rotor blades and the rational application on the characteristics, flow structure and kinetic energy losses of the axial turbine stages was given.
TURBINE STAGE; TANGENTIAL INCLINATION OF THE ROTOR BLADES; KINETIC ENERGY LOSSES; DEGREE OF REACTIVITY; ANSYS CFX.
Citation:
A. Q. Nguyen, K. L. Lapshin. The influences of tangential inclination of the rotor blades on the losses of kinetic energy in the axial turbine stage, Peter the Great St. Petersburg polytechnic university journal of engineerigsciences and technology, 23(04)(2017)66-73, DOI: 10.10871/JEST.230406.
Введение
Постановка задачи исследований
Как известно, тангенциальный наклон направляющих лопаток (ТННЛ) дает возможность эффективно снизить градиент степени реактивности в осевой турбинной ступени [1—9, 15]. Так, например, в [15] представлено сравнение результатов 3D проверочного расчета с помощью аэродинамического пакета ANSYS CFX с экспериментами, полученными на воздушной турбине в лаборатории аэродинамики кафедры турбиностроения Ленинградского политехнического института [1]. В расчетах и опытах исследовалось влияние ТННЛ на характеристики, структуру потока и потери кинетической энергии осевой турбинной ступени с отрицательным градиентом степени реактивности. Положительный опыт сравнения 3D расчетов с экспериментами [15] показал, что программу ANSYS СБХ можно рассматривать в качестве виртуального стенда и применять её для исследования характеристик и структуры потока дозвуковых осевых турбинных ступеней.
Тангенциальный наклон рабочих лопаток (ТНРЛ) по соображениям прочности применяется гораздо реже, в основном с целью компенсации изгибных напряжений, действующих на РЛ со стороны потока пара (или газа). Влияние ТНРЛ на характеристики, структуру потока и потери кинетической энергии осевой турбинной ступени практически не изучалось [16], и поэтому его исследование представляется актуальным.
Цели работы
Цель работы состоит в том, чтобы с помощью аэродинамического пакета ANSYS СБХ в качестве виртуального стенда исследовать влияние ТНРЛ на характеристики, структуру потока и потери кинетической энергии в осевой турбинной ступени. По результатам выполненных исследований планируется сформулировать выводы по рациональному применению ТНРЛ в осевой турбинной ступени.
Геометрические параметры исходной ступени взяты из книги [1]. Четыре варианта сту-
(
пени (рис. 1) средней веерности
^ = 6,34 1
\
имеют одинаковый направляющий аппарат. Угол а1эф постоянен вдоль радиуса и равен 18°10'. Линии центров масс направляющих лопаток ориентированы по радиусу. Рабочие лопатки не закручены, имеют постоянные профили и применены в разных вариантах (рис. 1) Вариант I (исходный): линия центров масс РЛ ориентирована по радиусу.
Вариант II: угол наклона РЛ на среднем радиусе (угол между осью г и линией центра масс РЛ в плоскости г-и) равен 10 градусам.
Вариант III: угол наклона РЛ на среднем радиусе (угол между осью г и линией центра масс РЛ в плоскости г-и) равен —10 градусам.
Вариант IV: угол наклона РЛ на среднем радиусе (угол между осью г и линией центра масс РЛ в плоскости г-и) равен —10 градусам; и линия центра масс на обводах проточной части в корневой и периферийной областях касается направлений радиусов.
Остальные геометрические характеристики ступеней следующие:
Радиус корневого сечения Щ . . . 159,1 мм
Длина НЛ 11.............58,2 мм
Хорда НЛ Ь1.............16,5 мм
Шаг НЛ г^..............11,86 мм
Число НЛ т1.............100 шт.
Длина РЛ 12 ............. 59,55 мм
Хорда РЛ Ь2с............. 27,85 мм
Шаг РЛ ^..............17,45 мм
Число РЛ т2 ............. 68 шт.
Осевой зазор между НЛ и РЛ Д^ . 8 мм Выполнялись 3D газодинамические расчеты с учетом рекомендаций, представленных в статьях [11—13]. Для упрощения решения задачи исследовались осевые турбинные ступени без радиального зазора над РК. Построение расчетных моделей выполнено с помощью программ, входящих в пакет ANSYS: Bladgen, TurboGrid и СБХ (ТигЬотасШпе). Расчетные модели являются
По "А-А"
О
Рис. 1. Схема осевых турбинных ступеней с различными вариантами установки рабочих лопаток Fig. 1. Scheme of axial turbine stages with variants of installation of rotor blades
секторами проточной части, состоящими из трех каналов НЛ и двух каналов РЛ. Вблизи поверхностей, обтекаемых потоком, обеспечено сгущение сетки, что необходимо для корректного расчёта параметров потока в пограничном слое. Диапазон частоты вращения ротора соот-
u
ветствует диапазону изменения параметра —
C0
от 0,46 до 0,62. Граничные условия следующие: Давление торможения потока
на входе в ступени P0*......0,137 МПа
Температура торможения
потока на входе в ступени T0*.....368 К
Давление в потоке на выходе
из ступеней P2 .........0,103298 МПа
Рабочее тело..........Air Ideal Gas
Модель турбулентности.........к-г
Постановка задачи...........Stage
При исследовании вариантов турбинных ступеней с помощью аэродинамического пакета ANSYS CFX использован способ решения Stage (стационарная задача), при котором на поверхности сопряжения между направляющим аппаратом и рабочим колесом выполняется осреднение параметров потока в окружном направлении. Применение Stage оправдано тем, что поток газа во всей области течения — дозвуковой, количество НЛ велико, осевой за-
зор между НЛ и РЛ также достаточно велик, поэтому степень неравномерности потока перед рабочими лопатками, вычисленная по известной формуле Г. Ю. Степанова, составляет примерно всего 3 %.
Результаты расчетов
После решения стационарной задачи с помощью аэродинамического пакета ANSYS CFX все полученные в 3D расчётах параметры газа усреднены в контрольных сечениях 0—0, 1—1 и 2—2 (см. рис. 1) по расходу. Затем определялись зависимости окружного КПД, степеней реактивности у корня р'т и у периферии р"т и потерь кинетической энергии ^ в направляющем аппарате и в рабочем колесе исследуемых вариантов осевых турбинных ступеней с ТНРЛ от
и
различных значений параметра —. Для вы-
С0
числения указанных величин использовались стандартные формулы, представленные, например, в [1].
Вычислялись также:
коэффициент потерь кинетической энергии для направляющего аппарата относительно располагаемого изоэнтропийного перепада энтальпий ступени Н0
ZHA =
С 2 - С 2 C1f C1 .
2H '
(1)
коэффициент потерь кинетической энергии для рабочего колеса относительно располагаемого изоэнтропийного перепада энтальпий ступени Н0
ZPK =
W2 - W22
2H
(2)
0
коэффициент потерь кинетической энергии с выходной скоростью относительно располагаемого изоэнтропийного перепада энтальпий ступени Н0
ZBC =
С2
2H
(3)
значении
С0
'opt
меньший КПД при значении
( \ u
С0
= 0,54
opt
0,84
коэффициент суммарных потерь кинетической энергии относительно располагаемого изоэнтропийного перепада энтальпий ступени Н0
<^=Сна +СРК +Свс . (4)
Полученные зависимости КПД, степеней реактивности у корня и у периферии, а также коэффициента суммарных потерь кинетиче-
и
скои энергии от параметра — представлены
на рис. 2—5.
На рис. 2 показаны зависимости КПД исследуемых ступеней от параметра —. Видно,
что при наличии как отрицательного, так и положительного угла наклона рабочих лопаток КПД исследуемых ступеней заметно уменьшается. Вариант I имеет наибольший КПД при
Рис. 2. Зависимости КПД от параметра —
C0
Fig. 2. Dependences of the efficiency
on the parameter — Cn
= 0,58, а вариант III — наи- -i
(см. рис. 5). Что важно, КПД варианта IV на 0,35 % больше КПД варианта III на оптимальном режиме [10]. Это обусловлено уменьшением коэффициента суммарных потерь кинетической энергии ступени варианта IV по сравнению с вариантом III (рис. 5), так как у корня и на периферии рабочая лопатка варианта IV приближается по свойствам к рабочей лопатке варианта I.
и/С0
Рис. 3. Зависимости степеней реактивности p'T
u
у корня от параметра — С0
Fig. 3. Dependences of the degree of reactivity p'T
u
on the parameter — Cn
Рис. 4. Зависимости степеней реактивности p"T
у периферии от параметра —
C0
Fig. 4. Dependences of the degree of reactivity p"T
—
on the parameter —
Рис. 5. Зависимости коэффициентов суммарных
с. и
потерь кинетической энергии от параметра —
Fig. 5. Dependences of coefficients of total losses of kinetic energy on the parameter —
На рис. 3 и 4 показаны зависимости степеней реактивности исследуемых ступеней рТ у
корня и рТ у периферии от параметра —. При
С0
положительном угле наклона рабочих лопаток (вариант II) степень реактивности у корня ступени увеличивается, а у периферии уменьшается, то есть уменьшается градиент степени реактивности по сравнению с исходным вариантом I. Наоборот, при отрицательном угле наклона рабочих лопаток (варианты III и IV) степень реактивности у корня уменьшается и увеличивается у периферии, то есть градиент степени реактивности увеличивается. Таким образом, тангенциальный наклон рабочих лопаток, как и тангенциальный наклон направляющих лопаток, является средством управления градиентом степени реактивности в осевой турбинной ступени.
Далее для более детального исследования влияния ТНРЛ на потери кинетической энергии в осевой турбинной ступени были рассмотрены графики изменения коэффициентов потерь кинетической энергии в направляющем аппарате и в рабочем колесе, а также выходной скорости по высоте проточной части на оптимальных режимах, на которых ступени достигают максимальных КПД. Соответствующие кривые для выходной скорости и коэффициентов потерь кинетической энергии в направляющем аппарате и в рабочем колесе представлены на рис. 6-8.
Следует подчеркнуть, что рабочие лопатки во всех вариантах ступеней не закручены. Поэтому даже в исходном варианте I на оптимальном режиме перед рабочими лопатками, в корневой зоне, неизбежно появятся положительные, а в периферийной зоне отрицательные углы атаки [14]. Эффективный угол р2эф в рабочих решетках нарастает вдоль радиуса, что приводит к увеличению выходной скорости с2 вдоль радиуса для всех исследованных вариантов (рис. 6).
Коэффициенты потерь кинетической энергии для этих вариантов несколько возрастают по сравнению с исходным вариантом I (рис. 7).
Понижение степени реактивности рТ у корня в вариантах III и IV увеличивает положительные углы атаки, вследствие чего коэффициенты потерь кинетической энергии увеличиваются в этой зоне по сравнению с вариантами I и II (рис. 8).
Рис. 6. Распределение выходной скорости по высоте лопатки Fig. 6. Distribution of the speed output by the height of the blades
Рис. 7. Распределение коэффициента потерь кинетической энергии для направляющего
аппарата по высоте лопатки Fig. 7. Distribution of the coefficient of losses of kinetic energy for the guide vanes by the height of the blades
Рис. 8. Распределение коэффициента потерь кинетической энергии для рабочего колеса
по высоте лопатки Fig. 8. Distribution of the coefficient of losses of kinetic energy for the impellers by the height of the blades
Коэффициенты потерь кинетической энергии Z2 вариантов III и IV почти совпадают по высоте лопатки, но у корня и на периферии для варианта IV они немного меньше, чем для варианта III. Этот различие объясняется тем, что в варианте IV на цилиндрических поверхностях тока сохраняются профили рабочих лопаток варианта III в срединной части, но меняются профили рабочих лопаток в корневой и периферийной областях, которые по конфигурации приближаются к варианту I. В итоге вариант IV развивает более высокий КПД, чем вариант III (см. рис. 2). По аналогии с «интегральными» сопловыми лопатками [9] назовём IV вариантом с «интегральными» рабочими лопатками.
Выводы
Влияние даже небольшого ТНРЛ на характеристики, структуру потока и потери кинетической энергии осевой турбинной ступени довольно существенно как при положительных углах наклона рабочих лопаток, так и, особенно, при отрицательных углах наклона рабочих лопаток.
С целью компенсации изгибных напряжений, действующих на рабочие лопатки со стороны газа, по соображениям экономичности рационально применять «интегральные» рабочие лопатки.
СПИСОК ЛИТЕРАТУРЫ
1. Афаснасьева Н. Н., Бусурин В. Н., Гоголев И. Г., Дьяконов Р. И. Аэродинамические характеристики ступеней тепловых турбин // Под общ. ред. В. А. Черникова. Л.: Машиностроение, Ленингр. отд-ние, 1980. 263 с., ил.
2. Кириллов И.И, Кириллов А. И. Теория турбо-машин. Примеры и задачи. Л.: Машиностроение (Ленинград. отд-ние), 1974. 320 с.
3. Кириллов И. И. Теория турбомашин. Л.: Машиностроение, 1972. 533 с.
4. Алексеев О. Н., Кириллов А. И., Лапшин К. Л., Черников В. А. Исследования турбинных ступеней с навалом сопловых лопаток // Труды ЛПИ. 1969. № 310. С. 35-40.
5. Кириллов И. И., Агафонов Б. Н., Афанасьева Н. Н., Кириллов А. И., Проскуряков Г. В., Черников В. А., Шварцман О. А. Исследование и отработка ступеней для высокотемпературных газовых турбин // Теплоэнергетика, 1983. № 3. С. 50-55.
6. Лапшин, К. Л. Оптимизация проточных частей осевых тепловых турбин: энциклопедия. М.: Машиностроение, 2015. Т. 4-19, глава 1.4. С. 49-75.
7. Топунов А. М., Тихомиров В. А. Управление потоком в тепловых турбинах. Л., 1979. 151 с.
8. Лапшин К. Л. Оптимизация проточных частей паровых и газовых турбин. СПб., СПбГПУ, 2013. [Эл. ресурс] < URL: http://dl.unilib.neva.ru/dl/2/3476.pdf>
9. Лапшин К. Л. Оптимизация проточных частей паровых турбин с применением «интегральных» сопловых лопаток // Научно-технические ведомости СПбГПУ. 2013. № 1(166). С. 61-66.
10. Осипов А. В., Анисин А. К. Теория и расчет одноступенчатой газовой турбины / Брянск, БГТУ, 2013. 73 с.
11. Смирнов Е. М., Кириллов А. И., Рис В. В. Опыт численного анализа пространственных турбулентных течений в турбомашинах // Научно-технические ведомости СПбГПУ. 2004. № 1(35). С. 56-70.
12. Кириллов А. И., Рис В. В., Смирнов Е. М. Опыт решения методических и практических задач численного моделирования течений и теплообмена в турбомашинах // Справочник. Инженерный журнал. Приложение: «Методы совершенствования энергетических установок». М.: Машиностроение. 2004. № 10. С. 13-19.
13. Епифанов А. А., Кириллов А. И., Рассохин В. А. Расчет трехмерного течения в ступенях малорасходных турбин // Научно-технические ведомости СПбГПУ. Наука и образование. 2012. № 1(142). С. 65-70.
14. Лапшин К. Л. Основы теории осевых тепловых турбомашин. СПб.: Изд-во СПбГПУ, 2016. 103 с.
15. Нгуен А. К., Лапшин К. Л. Характеристики и структура потока турбинной ступени с отрицательным градиентом степени реактивности // Научно-технические ведомости СПбГПУ. 2016. № 2(243). С. 163-173. БС1: 10.5962/ЛБ8Т.243.7.
16. Нгуен Ан Куанг, Лапшин К. Л. Влияние тангенциального наклона рабочих лопаток на характеристики осевой турбинной ступени// Неделя науки СПбПУ. Материалы научной конференции с международным участием.14-19 ноября 2016 года. Лучшие доклады. С. 12-15.
СВЕДЕНИЯ ОБ АВТОРАХ
НГУЕН Ан Куанг—аспирант Санкт-Петербургского политехнического университета Петра Великого. E-mail: [email protected]
ЛАПШИН Кирилл Леонидович — доктор технических наук профессор Санкт-Петербургского политехнического университета Петра Великого. E-mail: [email protected]
REFERENCES
[1] Afasnasyeva N. N., Busurin V. N., Gogolev I. G., Dy-akonov R. I. Aerodinamicheskiye kharakteristiki stupeney teplovykh turbin [Aerodynamic characteristics of thermal turbines stages] / Pod obshch. red. V. A. Chernikova. Leningrad: Mashinostroyeniye, Leningr. otd-niye, 1980. 263 p., il. (rus)
[2] Kirillov I.I, Kirillov A. I. Teoriya turbomashin. Primery i zadachi. [Theory of turbomachinery. Examples and exercises. Textbook for universities]. Leningrad: Mashinostroyeniye (Leningrad. otd-niye), 1974. 320 s. (rus)
[3] Kirillov I. I. Teoriya turbomashin. [Theory of turbomachinery]. Leningrad: Mashinostroyeniye, 1972. 533 s. (rus)
[4] Alekseyev O. N., Kirillov A. I., Lapshin K. L., Chernikov V. A. Issledovaniya turbinnykh stupeney s na-valom soplovykh lopatok [The research of turbine stages with the bulk of the nozzle vanes]. Trudy LPI. Leningrad. 1969. № 310. P. 35-40. (rus)
[5] Kirillov I.I, Agafonov B.N, Afanasyeva N.N, Kirillov A.I, Proskuryakov G.V, Chernikov V.A, Shvarts-man O. A. Issledovaniye i otrabotka stupeney dlya vysokotemperaturnykh gazovykh turbin. [The research and treatment stages for high temperature gas turbines]. Teploenergetika. 1983. № 3. S. 50-55. (rus)
[6] Lapshin, K. L. Optimizatsiya protochnykh chastey osevykh teplovykh turbin [Optimization of flow parts of axial thermal turbines]: entsiklopediya. M.: Mashinostroyeniye, 2015. T. 4-19, glava 1.4. S. 49-75.(rus)
[7] Topunov A. M., Tikhomirov V. A. Upravleniye po-tokom v teplovykh turbinakh. [Flow control in thermal turbines]. Leningrad, 1979. 151 s. (rus)
[8] Lapshin K. L. Optimizatsiya protochnykh chastey parovykh i gazovykh turbin [Optimization of flow parts of steam and gas turbines] / SPb., SPbGPU, 2013. [El. resurs] (rus) URL: http://dl.unilib.neva.ru/dl/2/3476.pdf
[9] Lapshin K. L. Optimizatsiya protochnykh chastey parovykh turbin s primeneniyem «integralnykh» soplovykh lopatok [Optimization of flow parts of steam turbines using "integrated" nozzle vanes]. NTV SPbGPU. 2013. № 1(166). S. 61-66. (rus)
[10] Osipov, A.V., Anisin, A. K. Teoriya i raschet odnos-tupenchatoy gazovoy turbiny [Theory and calculation of one-stage gas turbines] / Bryansk, BGTU, 2013. 73 s.(rus)
[11] Smirnov Ye.M., Kirillov A. I., Ris V. V. Opyt chislennogo analiza prostranstvennykh turbulentnykh techeniy v turbomashinakh [Experience of numerical analysis of spatial turbulent flows in turbomachines]. St. Petersburg State Politechnical Universiti Journal. 2004. № 1(35). S. 56-70. (rus)
[12] Kirillov A. I., Ris V. V., Smirnov Ye. M. Opyt resheniya metodicheskikh i prakticheskikh zadach chislennogo modelirovaniya techeniy i teploobmena v tur-bomashinakh[Experience of solving methodological and practical problems of numerical simulation of flows and heat exchange in turbomachines]. Spravochnik. Inzhen-ernyy zhurnal. Prilozheniye «Metody sovershenstvovaniya energeticheskikh ustanovok» [Application "Methods of improving power plants"]. Moscow: Machine building Publ, 2004. № 10. S. 13-19. (rus)
[13] Yepifanov A. A., Kirillov A. I., Rassokhin V. A. Raschet trekhmernogo techeniya v stupenyakh malor-askhodnykh turbin. Nauchno-tekhnicheskiye vedomosti SPbGPU. Nauka i obrazovaniye. 2012. № (142)1. S. 6570. (rus)
[14] Lapshin K. L. Osnovy teorii osevykh teplovykh turbomashin [Basics of theory of axial thermal turbines]. SPb.: Izd-vo SPbGPU, 2016. 103 s.(rus)
[15] Nguyen A. K., Lapshin K. L. Kharakteristiki i struktura potoka turbinnoy stupeni s otritsatelnym gradi-yentom stepeni reaktivnosti [Characteristics and structure of the flow in a turbine stage with a negative gradient of the degree of reactivity]. St. Petersburg State Politechnical Universiti Journal. 2016. № 2(246). S. 163-173. (rus), DOI: 10.5962/JEST.243.7
[16] Nguyen An Kuang, Lapshin K. L. Vliyaniye tan-gentsialnogo naklona rabochikh lopatok na kharakteris-tiki osevoy turbinnoy stupeni [The influences tangential inclination of the rotor blades on the characteristic of the axial turbine stage]. Nedelya nauki SPbPU. Materi-aly nauchnoy konferentsii s mezhdunarodnym uchastiyem. 14—19 noyabrya 2016goda. Luchshiye doklady. [Week of sciences SPbPU. Materials of the scientific conference with international participation. November 14-19, 2016. The best reports] S.12-15. (rus)
AUTHORS
NGUYEN An Quang — Peter the Great St. Petersburg polytechnic university. E-mail: [email protected]
LAPSHIN Kirill L. — Peter the Great St. Petersburg polytechnic university. E-mail: kirilUapschm@gmailxom
Дата поступления статьи в редакцию: 27 июня 2017 г.
© Санкт-Петербургский политехнический университет Петра Великого, 2017
