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

2015. Т. 19, № 3 (69). С. 65-69

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

УДК 621.51.226.2.53

Особенности изменения

радиального зазора и износа прирабатываемого покрытия в осевом компрессоре при флаттере лопаток

А. А. Хориков

[email protected]

ФГУП «Центральный институт авиационного моторостроения имени П. И. Баранова» (ФГУП «ЦИАМ им. П. И. Баранова»)

Поступила в редакцию 1 июля 2015 г.

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

Ключевые слова: радиальный зазор; флаттер лопаток; истираемое покрытие; износ покрытия; корпус компрессора, диаметральные формы.

ВВЕДЕНИЕ

Для обеспечения высоких параметров компрессоров разработчики всегда стремятся к тому, чтобы радиальный зазор между рабочими лопатками внутренней поверхности корпуса компрессора был минимальным. С целью избежания жесткого касания лопаток о корпус, на его внутреннюю поверхность наносится прирабатываемое (истираемое) покрытие. Вследствие закрученности лопаток или конусности корпуса при колебаниях лопатки происходит заметное изменение зазора [1, 2]. Особенно сильное изменение зазора происходит при флаттере лопаток, т.к. этому виду колебаний лопаток свойственны высокие амплитуды колебаний.

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

ступенчатого вентилятора. Видно, что при флаттере имеет место существенно большая неоднородность сигнала, как по амплитуде, так и по времени. В частности, при отсутствии флаттера временной интервал между импульсами сигналов одинаков и равен А( = 0,003 с, а при флаттере изменяется от t = 0,0025 с до t = 0,0035 с. При этом на прирабатываемом покрытии были обнаружены характерные особенности его износа.

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

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

Яял ■ с • --с

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

Рис. 1. Осциллограммы с датчиков зазора: а - нет флаттера; б - флаттер

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ИЗМЕНЕНИЯ ЗАЗОРА ПРИ ФЛАТТЕРЕ ЛОПАТОК

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

датливом диске. Соответствующая деформационная картина венца и сигналов с датчиков пульсаций представлена на рис. 2.

Рис. 2. Деформационные картины лопаточного венца и соответствующие им сигналы с датчиков изменения зазоров: а - неколеблющаяся решетка профилей; б - решетка профилей на податливом диске

при бегущих волнах деформации Предположим в первом приближении при отсутствии колебаний лопаток гармонический характер сигнала с датчика изменения зазора Л

л((,п) = Л0[ + а cos(ч/-фои)], (1)

где Шр = 2лfp, фоп = 2^n/N- угловая координата; fp - частота вращения; N - количество лопаток в колесе; а - амплитуда переменной составляющей сигнала с датчика зазора; n - порядковый номер лопатки. Начало отсчета может быть выбрано на любой лопатке.

При возникновении флаттера лопаток бегущая по колесу волна деформации с числом волн z угловую координату фоп в формуле (1)

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

Фп = ФСП + Л sin (t - 2яп / N), (2)

где 90 - величина, пропорциональная амплитуде крутильной составляющей колебаний периферийного сечения лопатки, а Ш2 = Ш л + zQp в соответствии с обобщенным эффектом

Доплера; Ш л - частота колебаний во вращающейся системе координат, z = 0,1...N -1.

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

А0 =Л0 (t,n)= A [l + S0 cos (ш/ - 2nn/N)], (3)

где 80 - величина, пропорциональная амплитуде колебаний лопатки в осевом направлении.

Между 9 и 5 справедливо соотношение

9 = d5 / dx, (4)

вследствие чего между ними существует сдвиг фаз 90°.

Подставляя выражения (3) и (2) в формулу (1) и принимая начало отсчета сигнала с лопатки с номером n = 0, получим

Л(() = A (l + 50 cos шг t) х

x[l + а cos (ш t + 90 sin шzt)].

Преобразуем выражение (5), используя формулу косинуса разности двух аргументов:

(5)

A(t) = A(l + 50 cosшzt)х

1 + а cos ш t • cos(90 sin шzt) -- а sinшt • sin (90 sin шzt) Поскольку справедливы разложения

cos (9o • sin шzt) = 10 (9q ) +

да

+ 2X 12k (90 )cos 2k Шzt '

k =1

sin(90 • sinшzt) =

да

= 2X12k-1 (90 )sin (2k - 1) Шzt'

k=1

(6)

(7)

где через In (&0) обозначена функция Бесселя первого рода n-го порядка от аргумента 90, то, применив эти разложения к формуле (6), после обычных алгебраических преобразований и раскрытий произведений функций будем иметь

Л (t) = A (l + S0 cos cozt) [l + а 10 () cos ct +

+ -

+

(- 1k E !k 9 )(cos Шt - Ыzt)-

k=1

да

XIk (90)cos (шt + Ыzt) ].

(8)

k=1

Из полученного выражения следует, что при возникновении флаттера датчик зазора, расположенный над лопатками компрессорного колеса, будет регистрировать следующие диагностические частоты флаттера:

coz = сл + zap, - частоту, обусловленную эффектом Доплера; а> = юр • N - частоту следования лопаток;

= c — kcz - нижнюю боковую полосу

Ш

z-k

частот;

Ш

z+k

= ф + - верхнюю боковую полосу частот.

Подставляя выражение (3) и (2) в формулу (1), будем иметь для любого номера п

A(t,n) = A

1 -50 cos I шЛ -

2лп

(

1- а cos

шрt-90 sin I шzt-

N 2жп

N

(9)

Формулу (9) можно применить для определения величины зазора по окружности колеса,

2жп

если угловую дискретную координату — за-

х

х

2л

менить на непрерывную координату -х, где

Л

х - расстояние по дуге, а Л - длина волны деформации. Тогда получим

х)=А

1-д0 соъI с А-

2пх

(

1- а cos

сорt-$0 ътI сС -

Л

2пх ^

Т

(10)

Минимальный зазор по координате х бу-

дет иметь место при выполнении двух условий к = 0,1,2...

2пх .

со А--=2пк,

2 Л

юр1-30 ът I с А -

2жх

Т

= 2як, к = 0,1,2. (11)

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

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

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

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

Поставленная задача решается тем, что реализуют способ диагностики флаттера рабочего колеса в составе осевой турбомашины на заданном рабочем режиме [3]. Способ характеризуется тем, что на внутреннюю поверхность корпуса турбомашины в зоне вращения рабочего колеса наносят истираемое покрытие, толщину которого устанавливают исходя из отсутствия радиального зазора между торцами лопаток рабочего колеса и поверхностью покрытия на максимальной частоте вращения. Выводят турбомашину на заданную частоту вращения. Обеспечивают подачу газового потока на входе в турбомашину, устанавливают заданные значения параметров температуры и давления газового потока и дросселируют турбомашину при заданных значениях параметров температуры и давления газового потока на заданной частоте вращения до заданной точки на напорной характеристике турбомашины. Диагностирование наличия или отсутствия флаттера на заданном рабочем режиме осуществляют после останова турбомашины по результату анализа характера износа истираемого покрытия.

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

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

х

X

Рис. 3. Схема экспериментального компрессора

Рис. 4. Износ истираемого покрытия при флаттере лопаток

Из рис. 4 следует, что в данном случае реализовался флаттер с тремя волнами деформации, т.е. реализовалась третья диаметральная форма колебаний. Диаметры колебаний показаны пунктирными линиями. Характер износа истираемого покрытия показан в виде круговой диаграммы толстой волнистой линией.

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

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

1. Beaupain O. Numerical simulation of blade/casing rub interaction. European Workshop on New Aero Engine Concepts. Munich, 2010, pp. 5-16.

2. Хориков А. А. Исследование колебаний лопаток

компрессоров датчиками измерения радиального зазора // Авиационно-космическая техника и технология.№ 8/55. 2008. С. 77-82. [A. A. Khorikov, "Research compressor blades oscillations with sensors measuring radial clearance", (in Russian), Aerospace technics and technology, no. 8/55, pp. 77-82, 2008.].

3.

Хориков А. А. Способ диагностики флаттера лопа-

ток рабочего колеса в составе осевой турбомашины. Патент на изобретение № 2525061. [A. A. Khorikov, The diagnostic method of rotor blades flutter in the composition of axial turbomachine. Patent 2525061 Russian Federation].

METADATA

Title: Radial clearance changes features and burned coatings depreciation in the axial compressor at blades flutter.

Authors: A. A. Khorikov.

Affiliation:

Central Institute of Aviation Motors after P.I. Baranov.

Email: [email protected].

Language: Russian

Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), -vol. 19, no. 3 (69), pp. 65-69, 2015.

ISSN 2225-2789 (Online), ISSN 1992-6502 (print).

Abstract: Mathematical model of change of radial clearance, and the specific characteristics of the coating wear burnished on the body in the event of axial compressor rotor blades of flutter was proposed. It is shown that in the event of blade flutter traveling wave deformation variation of the radial clearance can be represented by a set of Bessel functions of the first kind. A method for diagnosing a flutter on the rotor blades of the compressor characteristic features burnished wear coating has been proposed. Examples where, based on analyzing wear burnished coatings can infer the presence or absence of flutter in this mode of operation of the compressor, and to determine the number diametral waveform at which flutter have been implemented are shown.

Key words: Radial clearance, flutter blades, abradable coating, the wear coating, compressor housing, diametric form.

About authors:

KHORIKOV, Anatoliy Alekseevich, Doctor of Science, pro-

fessor, head of department. Strength, dynamics, calculation and experimental researches.

ОБ АВТОРАХ

ХОРИКОВ Анатолий Алексеевич, д-р техн. наук, проф, нач. отд., ФГУП «ЦИАМ им. П. И. Баранова».