Научная статья на тему 'Анализ аэроупругого поведения лопаточного венца в полуторной ступени осевого компрессора'

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

CC BY
89
41
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ТРЕХМЕРНЫЙ ПОТОК ИДЕАЛЬНОГО ГАЗА / ЛОПАТОЧНЫЙ ВЕНЕЦ / 3D IDEAL FLOW / BLADE ROW

Аннотация научной статьи по механике и машиностроению, автор научной работы — Гнесин В. И., Колодяжная Л. В., Огурцов К. В.

Частично-интегральный метод использован для численного моделирования нестационарных нагрузок и аэроупругих колебаний компрессорного лопаточного венца в трехмерном потоке идеального газа с учетом аэродинамического взаимодействия направляющего аппарата), рабочего колеса и статора)

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

Похожие темы научных работ по механике и машиностроению , автор научной работы — Гнесин В. И., Колодяжная Л. В., Огурцов К. В.

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

THE NUMERICAL ANALYSIS OF AEROELASTIC BEHAVIOUR OF BLADE ROW IN STAGE AND A HALF COMPRESSOR STAGE

The partial integration method has been used to numerical modelling of the unsteady loads and aeroelastic oscillations of the compressor blade row in 3D ideal gas flow with taking into account the stator – rotor – stator 2 aerodynamic interaction

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

■а о

Частково-інтегральний метод використано для чисельного моделювання нестаціонарних навантажень та аеропружних коливань компресорного лопатевого вінця у тривимірному потоці ідеального газу з урахуванням аеродинамічної взаємодії направляючого апарату, робочого колеса і статора

Ключові слова: тривимірний потік ідеального газу, лопатковий вінець

□-----------------------------------□

Частично-интегральный метод использован для численного моделирования нестационарных нагрузок и аэроупругих колебаний компрессорного лопаточного венца в трехмерном потоке идеального газа с учетом аэродинамического взаимодействия направляющего аппарата), рабочего колеса и статора)

Ключевые слова: трехмерный поток идеального газа, лопаточный венец

□-----------------------------------□

The partial integration method has been used to numerical modelling of the unsteady loads and aeroelastic oscillations of the compressor blade row in 3D ideal gas flow with taking into account the stator - rotor - stator 2 aerodynamic interaction

Key words: 3D ideal flow, blade row ------------------□ □----------------------

УДК 621.165

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

В.И. Гнесин

Доктор технических наук, профессор, заведующий

отделом

Отдел нестационарной газодинамики и аэроупругости* Контактный тел.: (057) 94-18-67 E-mail: [email protected] Л.В. Колодяжная Кандидат технических наук, старший научный сотрудник*

Контактный тел.: (057) 94-18-67 E-mail:[email protected] К.В. Огурцов Аспирант*

*Институт проблем машиностроения им. А.Н. Подгорного Национальной академии наук

Украины

ул. Дм. Пожарского, 2/10, г.Харьков, Украина, 61046

Введение

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

Традиционный подход к расчету флаттера об-лопаченных дисков основан на частотном анализе [1,2], в котором колебания лопатки представляются гармонической функцией во времени с постоянным сдвигом по фазе угла колебаний между соседними лопатками. Однако этот подход игнорирует эффект обратной связи воздействия колеблющихся лопаток на поток газа.

В настоящее время развиты новые подходы для исследования взаимодействия основного потока с колеблющимися лопатками, основанные на одновременном интегрировании во времени уравнений движения газа и колебаний лопаток с обменом информацией на каждой итерации [3-5].

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

В работе представлены амплитудно-частотные характеристики нестационарных аэродинамических нагрузок на рабочих лопатках и упругих колебаний лопаток.

1. Постановка связанной аэроупругой задачи

Трехмерный трансзвуковой поток идеального нетеплопроводного газа через полуторную ступень осевого компрессора рассматривается в физической области (рис.1), включающей статор0 (С0), рабочее колесо (РК), вращающееся с постоянной угловой скоростью, статор1 (С1).

Учитывая непериодичность потока в окружном направлении, расчетная область включает все лопатки С0, РК и С1. Расчетная область в осевом направлении состоит из трех подобластей, содержащих по одному лопаточному венцу и имеющих пересечения в осевых зазорах. В каждой из подобластей геометрические и

э

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

Рис. 1. Расчетная область

В каждый момент времени структура потока характеризуется периодичностью на угловом шаге 2п (на одном обороте ротора).

Расчетная область содержит zC0 + zPK + zC1 - сегментов ^С0^РК^С1 - числа лопаток лопаточных венцов), каждый из которых включает одну лопатку и имеет протяженность в окружном направлении, равную шагу статора или ротора. Каждый из сегментов расчетной области дискретизируется с использованием гибридной Н-Н разностной сетки. Разностные сетки в периферийных сечениях межлопаточных каналов ротора и статоров показаны на рис. 2.

Трехмерный трансзвуковой поток идеального газа описывается в относительной вращающейся с постоянной угловой скоростью декартовой системе координат полной системой нестационарных уравнений Эйлера, представленной в интегральной форме законов сохранения массы, импульса и энергии [6].

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

Принимая, что нестационарные возмущения вызваны вращением ротора в неравномерном потоке, а поток на входе в компрессор и на выходе содержит лишь малые отклонения от основного потока, постановка граничных условий основана на одномерной теории характеристик [6].

Динамическая модель колеблющейся лопатки с использованием модального подхода приводится к системе независимых обыкновенных дифференциальных уравнений относительно модальных коэффициентов собственных форм [7]:

Ч (і) + 2Мі (і) +ю2 (і) = X і (і),

(1)

где ч; - модальный коэффициент і-й формы; Ь;

- коэффициент механического демпфирования і-й формы; Ю; - собственная частота і-й формы; X; - модальная сила, соответствующая перемещению по і-й форме, которая рассчитывается на каждой итерации по мгновенному распределению давления на поверхности лопатки

Цр и; • п do

х,- = -

ДОр

а)

б)

в)

Рис. 2. Разностная сетка ротора и статоров: а) - С0; б) - РК; в) - С1

где р, р - давление и плотность на поверхности лопатки; и; (х,у^) - вектор перемещения лопатки по 1

- й форме.

Определив модальные коэффициенты qi из системы дифференциальных уравнений (1), найдем перемещение и скорость лопатки на каждой итерации

и (х,у,^)=Х и! (х,у^) qi (^ ,

1

и (х,у,^)=Х и! (х,у^) q 1 (^,

1

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

Для подтверждения достоверности предложенного метода проведено сопоставление численных результатов с экспериментальными данными для 11-ой стандартной конфигурации [2], которое

и

показалоудовлетворительное соответствие как по осреднен-ным характеристикам, так и по нестационарным [7].

2. Численные результаты

Численное исследование выполнено для полуторной ступени осевого компрессора с учетом аэродинамического взаимодействия лопаточных венцов и колебаний лопаток.

Соотношение чисел лопаток ротора и двух статоров zco:

Zpк: Zcl = 44:28:33.

Разностная сетка в каждом межлопаточном канале содержит 11x43x85, 11x67x79 и 11x57x69 сеточных узлов в радиальном, окружном и осевом направлениях для статора С0, ротора РК и статора С1 соответственно (рис. 2).

Расчет проведен для номинального режима работы компрессорной ступени:

- число оборотов ротора

п=12000 об/мин* Рис* 3* Распределение изолиний чисел Маха: а) — корневое сечение; б) — среднее

- полное давление в абсо- сечение; в) - периферийное сечение

лютной системе координат на

нагрузок. Ниже представлены результаты численного расчета аэродинамических нагрузок, действующих на рабочую лопатку, и колебаний лопатки в течение 12 оборотов ротора.

Анализ нестационарных характеристик проведен с использованием разложения в ряды Фурье исследуемых величин:

F(t) = Fo + ^Г Fli ■ соз(2п^ії) + F2І бш(2^й),

І=1

где F(t) - физическая нестационарная характеристика* Fo- среднее значение* FliИ F2i - коэффициенты Фурье* і-номер гармоники* V - 1-я гармоническая частота.

Так как угловая частота вращения ротора V рот=200 Гц, а время расчета соответствует времени 12 оборотам ротора ^=0.06 сек), первая гармоническая частота 1

равна V =— урет =16*7 Гц.

На рис. 4 приведены графики нестационарных нагрузок (окружная сила, осевая сила и момент относительно центра тяжести лопатки), действующих на периферийный слой рабочей лопатки в течение 12 оборотов ротора (рис. 4 а, в, д) и их амплитудно-частотные спектры (рис. 4 б, г,е).

Как видно из графиков, основной вклад в нестационарные составляющие аэродинамических нагрузок вносят гармоники, с частотами кратными числам лопаток статора С0 VpOTxZco =200x44=8800 Гц и статора С1 VpOT xz С1=200 x33 = 6600 Гц.

входе в ротор Р0 = 100000 Па*

- полная температура в абсолютной системе координат на входе Т = 288°К*

- углы потока на входе в ступень в радиальном и окружном направлениях заданы*

- статическое давление на выходе за статором р2 = 102000 Па.

Колебания рабочих лопаток определялись с учетом первых пяти собственных форм, собственные частоты которых приведены в таблице.

Таблица

Собственные частоты

Номер моды 1 2 3 4 5

> 450 1350 1800 3150 4050

На первом этапе выполнен расчет нестационарного потока в компрессорной ступени с учетом окружной неравномерности потока, но без учета колебаний лопаток. Этот режим характеризуется массовым расходом 17.8 кг/сек и степенью повышения давления пк=1.2.

На рис. 3 приведены изолинии числа Маха для корневого, среднего и периферийного сечений рабочей лопатки. Как видно из рисунков, течение характеризуется дозвуковыми скоростями на входе (Миз ~ 0.62 * 0.68) и на выходе (Миз ~ 0.56 * 0.58).

С некоторого момента времени начинается решение связанной задачи, в которой учитываются деформации рабочих лопаток под действием мгновенных

0.05 -А/АО

Й

>,0.04 -|

< 0.03 -

■и

0.00 0.01 0.02 0.03 0.04 0.05 0.06

с

0 4000 8000 12000 16000 20000

Частота, Гц

а)

б)

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

А/АО -

0.007

М, Н*м -0.006

О 4000 8000 12000 16000 20000

Частота, Гц

0.00 0.01 0.02 0.03 0.04 0.05 0.06

0.00 0.01 0.02 0.03 0.04 0.05 0.06

с

в)

0.25 -г А/АО-0.20 --

ООО-----1-

А0=0.00504

г)

Д)

О 4000 8000 12000 16000 20000

Частота, Гц

е)

Рис. 4. Аэродинамическая нагрузка в периферийном слое рабочей лопатки а)-6) - окружная сила; в)-г) - осевая сила;

д)~е) - аэродинамический момент

>

0.005

3.10

0.004

0.003

-0.048 Ьу, мм

-0.049

0 016 АУА0 Е?0 012

АО =-0.05(14

0.00 0.01 0.02 0.03 0.04 0.05 0.06

(, с

0 1000 2000 3000 4000

Частота, Гц

а)

б)

в)

-0.0104

0.01 в • А/АО

А 0=- 0.047 2

Ф , град. -0.0108 -

0.03 -

А/АО

0.02 -

О 1000 2000 3000 4000

Частота. Гц

г)

0.00 0.01 0.02 0.03 0.04 0.05 0.06

I, с

д)

О 2000 4000 6000 8000 10000

Частота, Гц

е)

Рис. 5. Колебания периферийного сечения рабочей лопатки: а)-6) - в окружном направлении; в)-г) - в осевом

направлении; д)~е) - угол поворота

Максимальные значения амплитуд нестационарных сил имеют место в периферийном слое и составляют 2. 1 *2.3% от средних значений окружной и осевой сил, и -20% для аэродинамического момента. Следует обратить внимание на присутствие в спектре низкочастотных возмущений, вызванных колебаниями лопаток.

Колебания лопатки представлены графиками перемещения периферийного сечения рабочей лопатки (в окружном направлении (рис. 5а), осевом направлении (рис. 5 в) и угол поворота (рис. 5 д), а также их амплитудно-частотными спектрами в течение 12 оборотов ротора. Мы наблюдаем устойчивые автоколебания. Основной вклад в изгибные колебания вносит гармоника с частотой близкой к частоте 1-ой собственной формы ( -400 Гц), в то время как крутильные колебания включают гармоники с частотами -400 Гц (1-я собственная форма) и -1800 Гц (3-я собственная форма).

Следует подчеркнуть, что вынужденные высокочастотные колебания (6600 Гц, 8800 Гц) пренебре-

жимо малы по сравнению с амплитудами автоколебаний.

Выводы

Проведены численные исследования для полуторной ступени осевого компрессора (статор+рото-р+статор) с учетом взаимного аэродинамического взаимодействия лопаточных венцов и колебания лопаток.

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

Колебания лопаток ротора по всем формам являются устойчивыми.

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

Литература

1. Bakhle, M.A. Time domain flutter analysis of cascades using a full - potential solver [Текст] / M.A. Bakhle, T.S.R. Reddy and T.G.

Keith //J. AIAA.- 1992.-30, N1.- P.163-172.

2. Bolcs A. Aeroelasticity in Turbomachines - Comparison of Theoretical and Experimental Cascade Results [Текст] / A. Bolcs, T.H.

Fransson// Communication du LTAT-EPFL, Lausanne, Switzerland. - 1986. - 13. - 174p.

3. Moyroud F. A Modal Coupling for Fluid and Structure Analysis of Turbomachine Flutter. Application to a Fan Stage [Текст] / F. Moyroud, G. Jacquet-Richardet and T.H. Fransson // ASME Paper 96-GT-335.- 1996.-P. 1-19.

4. Chew, J.W. Part-speed flutter analysis of a wide - chord fan blade [Текст] / J.W. Chew, J.G. Marshall, M. Vandati and M. Imregun//

In: Fransson T.H. (Ed.). Unsteady Aerodynamics and Aeroelasticity of Turbomachines. Kluwer Academic Publishers, Dordrecht.-1998.-P.707-724.

5. Gnesin V.I. A coupled fluid-structure analysis for 3-D flutter in turbomachines [Текст] / V.I. Gnesin, R. Rzadkowski, and L. Kolodya-

zhnaya // ASME J. 2000- GT-380, International Gas Turbine and Aeroengine Congress, Munich,Germany.- 2000.- 8p.

6. Гнесин В.И. Численное анализ влияния статора и ротора на нестационарные нагрузки и режимы колебаний лопаток [Текст] /

В.И. Гнесин, Л.В. Колодяжная //Вестник НТУ ХПИ, Энергетические и теплотехнические процессы и оборудование.- 2009.-No3.- С. 23-32.

7. Gnesin V.I. A numerical modelling of stator- rotor interaction in a turbine stage with oscillating blades [Текст] / V.I. Gnesin, L.V

Kolodyazhnaya and R. Rzadkowski // J. of Fluid and Structure.-2004.- No 19, P. 1141-1153.

3

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