CC BY
105
АВИАЦИОННАЯ И РАКЕТНО-КОСМИЧЕСКАЯ ТЕХНИКА
УДК 629.735
И. А. Кривошеев, К. Е. Рожков
АНАЛИЗ РАБОЧИХ ПРОЦЕССОВ И ХАРАКТЕРИСТИК ЛОПАТОЧНЫХ МАШИН
При проектировании, испытаниях и эксплуатации лопаточных машин (компрессоров, турбин) решаются разнообразные задачи, где используются характеристики лопаточных венцов (рабочих колес и направляющих аппаратов), ступеней, многоступенчатых компрессоров (каскадов). Выделяются задачи определения потребной геометрии (конструкции) по заданным термодинамическим и кинематическим параметрам, и наоборот - при заданной геометрии и некоторых режимных параметрах (например, частота вращения ротора, расход или скорость потока на входе) определения остальных параметров. Авторами предлагается развитие экспериментально-теоретических методов получения и использования характеристик лопаточных машин. Авиационные двигатели; расчет характеристик; имитационное моделирование
Сложную картину пространственного (трехмерного) течения через лопаточный венец можно представить как совокупность двухмерных течений: осесимметричного, течения по поверхности вращения и вторичного.
Применительно к решетке компрессора часто используется допущение, что поверхность вращения представляет собой круговой цилиндр. Решетка профилей, полученная в результате развертки на плоскость, представлена на рис. 1.
Рис. 1. Схема и геометрические параметры компрессорной решетки
Важнейшим параметром решетки являются шаг (0 и густота решетки (Ь / t - отношение к шагу величины хорды), соединяющей точки пересечения средней линии профиля (см. рис. 1, пунктир) с контуром профиля. Помимо шага и густоты также важен угол установки профиля
Контактная информация: (347) 273-06-35
Работа поддержана грантом Президента Российской
Федерации НШ-65497.2010.9 для ведущих научных
в решетке 0 (угол между хордой и фронтом решетки). Углы потока на входе р1 и на выходе из решетки р2, определенные треугольником скоростей, отличаются от конструктивных углов лопаток р1л и р2л, определяемых между касательными к средней линии и фронтом решетки на угол атаки 1 по передней кромке и на угол отставания потока 5 на выходе из решетки. Угол поворота потока в решетке Др определяется как ДР = Р2-Рь а разность конструктивных углов лопаток определяет угол изгиба профиля 0 = = Р2л-Р1л. Вводя углы 1 и 5, получим Др = 0 + + 1 - 5. Параметрами определяющими пропускную способность решетки при больших скоростях набегавшего потока, является ширина узкого сечения (так называемого «горла») решетки (Аг) и отношение его к ширине струи на входе (Аг/А1). Помимо угла изгиба профиля 0 форма профиля характеризуется величинами х^~ расстоянием вдоль хорды от носика профиля до точки максимального прогиба, стах - максимальной толщиной профиля и хс - координатой положения максимальной толщины. Используются соответствующие безразмерные величины: с = стах/Ь; х~ = хс/Ь; = Xf/b. Задача аэро-
динамики решеток состоит в том, чтобы при заданном треугольнике скоростей определить потери в решетке, режим работы решетки с оптимальными потерями и обеспечить заданный угол поворота потока, а для этого надо определить оптимальный угол атаки и угол отставания потока. Ранее, в связи с диффузорным течением в решетках компрессора эти характеристики определялись экспериментально. В настоящее время технические возможности позволяют получить эти характеристики еще до реального появления лопаток в металле. Рассмотрим, прежде всего, обобщенные данные исследования плоских компрессорных решеток при малых скоростях набегающего потока. Типичные ре-
школ
зультаты продувки плоской компрессорной решетки - зависимости угла поворота потока Др
и коэффициента потерь ^ = ^^2 от угла атаки - представлены на рис. 2.
При отрицательных и небольших положительных углах атаки угол отставания потока 5 практически не изменяется и, следовательно, угол выхода потока из решетки р2 остается практически неизменным.
Рис. 2. Результаты продувок типовой компрессорной решетки
Поэтому при увеличении угла атаки (т. е. при уменьшении угла входа потока на решетку Рі) угол поворота потока Др увеличивается. При некотором положительном угле атаки /кр возникает срывное течение с выпуклой поверхности профиля и, несмотря на уменьшение угла потока р1, величина Др не увеличивается, так как на срывном режиме увеличивается угол отставания 5 и р2 уменьшается. Потери в решетке при этом примерно в два раза больше, чем на режиме минимальных потерь, характеризуемом углом /'хтт при небольших отрицательных углах атаки. На режиме /'хішп угол поворота потока в решетке Др мал и, следовательно, мал теоретический напор ступени. Наибольший к.п.д решетки достигается при некотором оптимальном (положительном) угле атаки /ор1. На этом режиме потери в решетке хотя и несколько больше минимальных, но повышенное значение угла поворота потока ДР и, следовательно, величины напора ступени по сравнению с режимом 7'хтт обеспечивают максимальный к. п. д. решетки. Картина течения в плоской компрессорной решетке находится в полном соответствии с картиной течения в плоском диффузоре. При течении
в плоском диффузоре также наблюдается, что режим максимума к.п.д. наступает ранее, чем режим возникновения срыва. Эти режимы определяются двумя параметрами: углом раскрытия диффузора 0 и величиной b/n1. Если сопоставить параметры компрессорной решетки с эквивалентным диффузором, то будем иметь 0 = = t/b(sinP2 - sinPi) и b/n1 = b/(tsinPi), т. е. режим работы решетки определяется углом поворота потока DP = p2 - Pi, углом Pi или p2 = Dp - Pi и густотой решетки b/t.
Экспериментальные исследования плоских компрессорных решеток показали, что перечисленные параметры оказывают определяющее влияние на режим работы решеток, а влияние таких параметров, как кривизна профиля, его относительная толщина оказывается значительно менее существенным. Поэтому обобщенные данные по плоским решеткам строятся на основе влияния величин b/t, р2 (или р1) на Др. Прежде чем привести такие обобщенные данные по результатам испытаний решеток при малых скоростях набегающего потока, отметим, что выбор режима работы решетки вблизи Дртах не представляется возможным не только потому, что при этом к.п.д. ступени падает, но и в связи с тем, что при возникновении срывного режима возникает неустойчивость течения в компрессоре, когда эксплуатация компрессора невозможна или сильно затруднена. Поэтому режим работы решетки и все обобщения строят для некоторого номинального угла поворота потока Др* =
0,8Дртах для решеток РК и Да* = 0,8Датах для решеток НА, обеспечивая необходимый запас относительно срывных режимов.
Характеристики решеток профилей могут дать информацию об их влиянии на различные параметры работы.
Выполненная авторами Са/Са0 экстраполяция коэффициентов Ольштейна Л. Е К =
= /(Q^Qio^wi) рис. з. и К2 = /(Q^Qio^wi) рис. 4. - функций, которые используются для построения соответственно кривых Нт =
= /(0°, *wi) и Н0 = /(°0°,_Awi) показала, что при °а = °а° Ki = К2 = Яту = °а0 ив эту же точку приходит касательная к кривой К2 = = / (0а/0аО<А wi), построенная в точке °°а = Са°, т. е. Са/Са0 = °. Однако это не согласуется с практикой - продувочные кривые решеток, линии теоретического напора зависят от густоты решеток b/t, от кривизны 0 и других геометрических параметров решеток.
Поэтому, в связи с разработкой в НИЛ САПР-Д УГАТУ систем имитационного моделирования компрессоров (СИМ KOMPR, STU-
РЕКУ и т. д.), с участием авторов проведен анализ экспериментальных данных, различных эмпирических и теоретических зависимостей, выведен ряд собственных зависимостей и предложена следующая методика: предлагается брать за расчетный режим именно номинальный, который отличается от оптимального, т. е. расчетный режим отдельных ступеней смещен по отношению к режиму с максимальным КПД и находится в области с более низким КПД, но с более высоким запасом устойчивости.
0,3 Кі
0 Са / С а 0
7 0 8 0 9 1 1 21,
-0,2 / ^ж1 Т
Рис. 3
0 3 К2
Са / Са0
7 0 8 0 90 1 1 1 2 1
/^м>1 Т
-0 , 3
Рис. 4
Для расчета выведен ряд дополнительных формул и предложена комбинация использования зависимостей, обнаруженных Л. Е. Оль-штейном и Хауэлом. Последняя дополнительно уточняется по зависимостям, приведенным К. В. Холщевниковым Для иллюстрации предложенной методики ниже приведен расчет характеристики первой ступени компрессора низкого давления (КНД) ГТП АЛ-31СТ.
В качестве исходных данных для данного расчета приняты следующие геометрические параметры 1ст. КНД АЛ-31СТ: р1л = 41,47, р2л =
= 52,13о, 0 = 10,66о, Ь / ґ = 1,1, X/ = 0,55 (сверхзвуковая решетка), угол выхода из ВНА аіл = = 55,5о.
Последовательность определения НТ0, са0 и П0:
1. Для получения параметров, характеризующих номинальный режим, необходимо построить зависимость Нт0 = /(са0) . Упрощенно это делается по формуле Хауэла 1,55
■С„
1 +1,55
Более точно это можно
сделать по зависимостям, приведенным в [5], для этого для данной густоты построены зависимости Н т / С а = / (р / С а) (1), после преобразований получаем искомую зависимость.
2. Учитывая, что с учетом закрутки за ВНА — 1
=------------, строим линию идеального
с^ а1эф + с^ 01
теоретического напора (при нулевом угле отставания на выходе из РК) ^тид = 1 —
—а1эф + ^ р2л) (в данном случае это луч, идущий из т. {1;0} в т. {0;0,6826}).
3. Определяем значения 50 и Р20 на номинальном режиме, где по определению
ДРо = 0,8ДРтах . Для определения 50 используется зависимость полученная авторами при анализе нескольких источников [1, 2, 3, 4] где предлагается следующая зависимость для определения угла отставания на номинальном режиме б = тб-^/Ь, где т = 0,92:; + 0,00202 [5] и т = 0,18 + 0,92;:; — 0,00202 [1]. При сравнении с экспериментальными продувками решеток профилей было установлено что зависимость [2] является более точной и при сравнении с экспериментальными продувками дает погрешность не более 1°-2°, и поэтому при получении формулы для оптимального угла отставания используется данная зависимость
2 " 0,18 + 0,92• ~хґ -0,002Д2д Ад-А >|Ь
1 - 0,002 • (Ад-Дд) л Ьь
^0 =
Для определения Др используется зависимость Др = 0,308 х р20 - 2,33, приведенная в [1], где р20 = р2л - 50. Данная зависимость пригодна для решеток профилей с густотой Ь / ^ = 1, для решеток профилей с другой густотой используется поправочный коэффициент / = ДР/ДР(ЬЛ = 1), где / = -0,109 х (Ь/04 + 0,628 х (Ь/03 - 1,48 х х (Ь/02 + 2,292хЬ/^ - 0,331.
В дальнейшем возможно уточнение номинального режима численным моделированием в САЕ-системе.
4. Проводим луч из т.{1;0} с учетом угла отставания на номинальном режиме - в точку на
1
оси абсцисс Сп =------------—. Авторами пока-
“ ^ аlэф+ctg 020
0
3
0
3
зано, что именно этот луч характеризует номинальный режим на кривой Нт = /(00). В рассматриваемом примере луч проходит на оси абсцисс через точку {0;0,605} и определяется на пересечении с лучом Хауэлла номинальный режим, где значения Нто = 0,24 и оа° =0,46 ^ ctg Рю = 1 / Са0 - ctg а1эф = 1,4866, Рю = 33,93о ^ /0 = 7,54о. Располагая графиками эмпирических зависимостей [5] изменения коэффициентов К1 и К2, получим значения °°то = 0,24 и 0000 =0,46.
По выражению н = к +1 I • Н можно поТ 1 і Са0 0
строить требуемую кривую °°т = /(0°°). Авторами показано, что аппроксимация эмпирической зависимости для К1 в зоне 0°0 < 1 сводится к лучу, идущему из точки на оси ординат К1 = = 0,75. Далее необходимо определить коэффициент изоэнтропического напора и КПД ступени на расчетном (номинальном) режиме (рис. 5).
Для этого воспользуемся зависимостями
К2 +
с„
\
С
а 0
• Н
- К2 +
* = Н=_2_
Пк Н „ с
-Н
*0
К 2 +:
Н
а0
-а0
с Н
К +^НТ0 К +£^^°
' ' %
Далее задаемся диапазоном изменения [0; 1]. Находим изоэнтропическую работу как
Н
Н= ^=- -18,32 • 288,15 • 1а и строим график
Са2
Н =_Д^1а) (рис. 6).
При этом самая верхняя изолиния с; = =const характеризует границу помпажа.
Полученный результат нетрудно представить в традиционных координатах п = =Хц] (рис. 7). На приведенных графиках не показаны режимы запирания, возникающие в «горле» РК или НА.
Для верификации разработанного метода расчета произведено численное 3Б САБ-/САЕ-моделирование.
Для построения математической модели была выбрана первая ступень осевого компрессора низкого давления (КНД). Компрессор состоит из четырех ступеней. Первая ступень состоит из входного направляющего аппарата (ВНА), рабочего колеса первой ступени (РК1) и направляющего аппарата первой ступени (НА1). ВНА состоит из 23 лопаток, РК1 - из 37 лопаток, НА1 - из 54 лопаток.
а0
'а 0 '10
с
а
а
Рис.5. Построение характеристики 1 ступени КНД АЛ-31СТ Н8 = /(с ) и т)! = /(са) на основе
геометрии решеток, с предложенным авторами методом использования эмпирических (обобщенных)
зависимостей
Рис. б. Характеристика первой ступени КНД АЛ31-СТ Hs = f(kia)
1,8 1,7 1,6 1,5 1,4 1,3 1,2 1,1 1
10 20 ЗО 40 50 60 70
Рис. 7. Характеристика первой ступени компрессора АЛ31 -СТ п = f[G; Уи]
Для расчета данная модель разбита на три составные части: входной направляющий аппарат, рабочее колесо (ротор) и направляющий аппарат. На каждой части модели необходимо построить сетку. Для этого необходимо при помощи CAE-системы CFX-Mesh задать параметры сетки и построить её для каждого элемента в отдельности.
На рис. 9 показан результат построения сетки ВНА. Далее необходимо запустить CFX-Pre, загрузить все три части модели с построенной сеткой (рис. 10).
Затем необходимо создать новый вращающийся домен Rotor, нажав на кнопку Create a Domain. В настройках домена указать: месторасположение (Location) - часть модели с рабочим колесом; рабочее тело (Fluids list) - Air Ideal Gas; ссылочное давление (Reference pres-
Частота вращения рабочего колеса (на номинальном режиме работы в составе ГПА) составляет 9,3 тыс. об./мин. Рабочей средой является воздух. Трехмерная модель первой ступени компрессора была построена в системе и^гарЫс8 КХ6 и представлена на рис. 8.
Рис. 8. Расчетная CAD модель 1 ступени КНД (модель ПЧ для газодинамического расчета)
sure) - 101325 Па; указать на то, что домен -вращающийся Domain Motion > Option - Rotating; частота вращения (Angular Velocity) 9300 об./мин.; ось вращения - Global X. Во вкладке Fluid Models поменять настройки модели тепло-переноса (Heat Transfer) на Total Energy, так как мы ожидаем, что в модели будут участки, где необходимо учитывать сжимаемость воздуха. В неподвижном домене (Default Domain) создать граничное условие Inlet, нажав на кнопку Create a Boundary Condition, с параметрами: тип граничного условия - Inlet, опция задания Mass and Momentum - Mass Flow Rate; Mass Flow Rate = 5,33 kg/s; Heat Transfer > Total Temperature = 288 K
Рис. 9. Сетка построенная на модели ВНА
Рис 10. CAE -модель с приложенными граничными условиями
Аналогичным образом задаем граничное условие на выходе Outlet: тип граничного условия - Opening; опция задания Mass and Momentum - Opening Pres. and Dirn.; Relative Pressure -
0 Па; Heat Transfer > Opening Temperature = = 288 K.
Так как расчётная модель представляет собой 1/12 всего компрессора, то на боковых гранях необходимо установить условие периодичности. Для этого необходимо нажать кнопку Create a Domain Interface, выбрать Interface Side
1 > Domain (Filter) - Default Domain; Region List - PerVNAl; выбрать Interface Side 2 > Domain (Filter) - Default DomainRegion List -PerVNA2; Interface Models > Option - Rotational Periodicity, Rotation Axis - Global X (указали, что взаимосвязь - вращающаяся периодичность с осью вращения ОХ). Подобным образом установить условие периодичности для рабочего колеса и направляющего аппарата.
Критерием сходимости является достижение величиной среднеквадратичной невязки уровня 10-4 или ниже в совокупности с установлением интегральных характеристик, что достигается приблизительно за 500 шагов по времени и соответствует прохождению потоком ~10 продольных размеров ступени. При проектировании и доводке можно использовать, как для получения характеристик, так и для оптимизации узла, Ansys CFX, но при этом не стоит забывать, что при расчете возможны погрешности, возможные причины погрешности - это небольшое количество элементов, на которое была разбита модель. При расчетах характеристик были использованы экспериментальные данные, данные характеристик показывают, что возможно, зная геометрию лопатки, получить характеристики компрессора без экспериментальных исследований по продувке компрессора. Разработанный алгоритм может лежать в основе компьютерных программ по расчету компрессоров и двигателя в целом. Такие системы моделирования (KOMPR, Stupeny,...) разрабатываются в НИЛ САПР-Д каф. АД УГАТУ. В отличие от них, программа, составленная по данному алгоритму, будет способна рассчитывать характеристики компрессора по его геометрии, во всем диапазоне режимов работы, в том числе в области авторотации, помпажа, зуда (рис. 11). При этом не будет необходимости (как это реализовано сейчас) вначале подробно подбирать параметры и геометрические размеры для расчетного режима.
Также разработанный метод построения упрощенных характеристик компрессоров может быть использован для планирования экспериментов, например для определения области границы помпажа, нахождение которой позволит сохранить двигатель при испытаниях. Кроме того, по упрощенным характеристикам можно найти область интересующей точки (точка, которая рассматривается в эксперименте), что позволит уменьшить количество экспериментов.
Предложенный алгоритм полезен для сглаживания экспериментально полученных характеристик, для их экстраполяции (в зону запуска, авторотации, запирания, обратных течений при помпаже), для выявления запомпажных частей напорных веток, для прогнозирования и уточнения протекания границы помпажа.
ЗАКЛЮЧЕНИЕ
Таким образом, в результате рассмотрения характеристик компрессоров, закономерностей их протекания и методов их получения предложен оригинальный графоаналитический метод построения упрощенных и полных характеристик ступеней, их суммирования и получения характеристики многоступенчатых компрессоров, метод обобщенного построения универсальной характеристики компрессора, метод планирования на этой основе эксперимента по получению характеристик и метод идентификации характеристик компрессоров по результатам испытаний двигателей.
СПИСОК ЛИТЕРАТУРЫ
1. Белоусов А. Н., Мусаткин Н. Ф.. Радь-
ко В. М. Теория и расчет авиационных лопаточных машин. Самара: Гос. аэрокосм. ун-т, 2003. 344 с.
2. Бунимович А. И., Орлова Г. С. Сборник аэродинамических характеристик плоских компрессорных решеток. М.: ЦИАМ, 1955. 170 с.
3. Федоров Р.М. Альбом характеристик компрессоров. М.: ВВИА им. Жуковского, 1981. 181 с.
4. Холщевников К. В., Емин О. Н., Митрохин В. Т. Теория и расчет авиационных лопаточных машин. М.: Машиностроение, 1986. 432 с.
5. Холщевников К. В. Теория и расчет авиационных лопаточных машин. М.: Машиностроение, 1970. 610 с.
6. Кривошеев И. А., Ахмедзянов Д. А., Киша-
лов А. Е. Получение характеристик осевых компрессоров. Уфа: УГАТУ.
ОБ АВТОРАХ
Кривошеев Игорь Александрович, декан факультета АД, проф., науч. рук. НИЛ САПР-Д. Дипл. инж.-мех. (УГАТУ, 1976). Д-р техн. наук по тепл. двигателям ЛА (2000). Иссл. в обл. инф. техн. в двигателестроении.
Рожков Кирилл Евгеньевич,
асп. той же каф. Дипл. инженер-механик (УГАТУ, 2009), Иссл. в обл. моделирования рабочих процессов и анализа характеристик лопаточных машин.
