Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

EXTRAPOLATION OF THE EXPERIMENTAL OBTAINED COMPRESSOR SPEED CURVE ON COMPRESSOR´S MAP

Krivosheev I.A. 1 Rozhkov K.E. 1 Simonov N.B. 1
1 Ufa State Aviation Technical University
Experimental obtained compressor speed curve extrapolation method has been described. This method can be applied to calculated speed curve. Theoretical head rise nominal and maximal point determining method has been proposed using nomogram formation method. Graphic explanations extrapolation method has been demonstrated in general and particular cases for specific compressors. Method verification has been carried out for example low pressure compressors (already have been using for more than 5 years) at different rotating speed. Accomplished by authors analysis showed that stability loss of the compressor system (GTE or a compressor with a receiver and a nozzle at the stand) is corresponded with point on speed curve where isentropic head rise maximum is achieved. Results comparison has been carried out with known numerical and experimental engine´s characteristics.
extrapolation
speed curve
compressor´s map
compressor

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

При экспериментальном получении характеристик компрессоров (рис.1), особенно если испытания проводятся в составе ГТД, возникает проблема – достаточно трудно получить напорные ветки в «нижней части» характеристики (при малых частотах вращения ротора, начиная с нулевой), а такая характеристика нужна для моделирования запуска двигателя. Кроме того, практически невозможно точно определить точки на границе помпажа, поскольку вводить двигатель в режим помпажных автоколебаний нельзя – возникает заброс температуры в камере сгорания, ухудшается охлаждение лопаток турбины и возможен их обрыв, из-за колебаний возможен обрыв лопаток компрессора. Сложно находить точки в зоне запирания (направляющих аппаратов и рабочих колес, при малых частотах вращения – последних ступеней, при больших частотах вращения – первых ступеней). Тем более не удается получать запомпажную часть напорных веток (а такая характеристика требуется для моделирования помпажных колебаний в двигателе, настройки и отладки системы защиты от помпажа в составе САУ ГТД). Обычно характеристику компрессора получают расчетным или экспериментальным путем в виде и .

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

В данной статье предложен способ экстраполяции экспериментально получаемых напорных веток, на которых несколько точек (в области устойчивой работы компрессора) определены экспериментально. Способ позволяет выделить на экстраполированных напорных ветках точки протекания границы помпажа и границы запирания, а также точки номинального и максимального значения теоретического напора. Это позволяет с высокой степенью достоверности нанести на характеристике компрессора границу помпажа и линию запираний. Кроме того, способ позволяет на основе полученных напорных веток построить новые – с требуемым шагом по частоте вращения, вплоть до «нулевой» (начало раскрутки ротора при запуске). Предлагается определенная последовательность действий над компрессором и определенная последовательность обработки результатов испытаний. При этом в отличие от известных способов, для построения границы помпажа не требуется при испытаниях вводить компрессор в помпаж и фиксировать помпажные точки на напорных ветках характеристики. Достаточно экспериментально зафиксировать на каждой напорной ветке по три точки (желательно одна – вблизи номинального режима, а две другие выше и ниже его по напору, но до запирания), в обычном диапазоне устойчивой работы компрессора. Далее с использованием результатов этих измерений производятся расчеты и построения на основе установленных авторами закономерностей и выведенных (в развитие метода Ольштейна Л.Е.) формул. По результатам строится граница помпажа и экстраполируются напорные ветки (а также строятся новые), что и составляет искомую характеристику компрессора. Как известно, помпажные автоколебания являются недопустимым режимом работы ГТД (в т.ч. компрессора в его составе) и возможность получения границы помпажа по результатам испытаний в зоне режимов устойчивой работы важна для практики.

Описание предлагаемой методики экстраполяции напорных веток и характеристик компрессоров

В предлагаемой методике, как и в известных, при испытаниях компрессора авиационного газотурбинного двигателя на стенде или в составе двигателя, измерения производят на каждой напорной ветке при фиксированной частоте вращения, при этом различные точки получают за счет изменения степени дросселирования проходного сечения за компрессором или перед ним. Однако, в отличие от известных способов [1-5], для определения границы газодинамической устойчивости и построения напорных веток предлагается следующая последовательность действий. Выполняются испытания компрессора – либо автономно (как показано на рис. 2), либо в составе двигателя (на стенде или на самолете, в т.ч. в эксплуатации – при контроле параметров). Испытания начинаются с наименьшей частоты, при которой двигатель или компрессор на стенде работают устойчиво. При каждом испытании компрессора измеряются и фиксируются: давление и температура окружающей среды, параметры в трех точках напорной ветки (при заданной частоте вращения) вблизи номинального режима: частота вращения n, статическое давление в мерном сечении лемнискаты на входе, заторможенные давление и температура за компрессором. По этим параметрам стандартным образом определяется расход воздуха G, газодинамическая функция плотности тока q(λ1a) и приведенная скорость λ1a на входе в компрессор. В каждой точке определяется степень повышения давления в компрессоре и изоэнтропический и теоретический напор по формулам:

(1)

(2)

(3)

Строится фрагмент характеристики компрессора (рис. 3) в координатах и . Как показано на рис.1, каждая напорная ветвь (при фиксированной приведенной частоте вращения nпр) в координатах (λ1a) изображается дугой, ограниченной справа вертикальной линией режимов запирания, и расположенным над отрезком прямой . Как показал анализ, на практике линейность зависимости особенно четко обеспечивается при низких частотах вращения компрессора (0,3-0,4 от ). Такая линейность прежде всего наблюдается для элементарной (по Н.Е.Жуковскому) решетки профилей или элементарной ступени, но согласно публикациям Ольштейна Л.Е., при изображении характеристики ступени в виде такая линейность также имеет место вблизи номинального режима.

Рис. 1. Графическая интерпретация отрезков напорных веток, построенных по результатам испытаний компрессора (в области устойчивой работы) для теоретического (λ1a) и изоэнтропического напоров (λ1a) при фиксированной частоте вращения n

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

Рис. 2. Линейная экстраполяция отрезка напорной ветки (λ1a) до оси ординат, определение значения , построение касательной к напорной ветке (λ1a), определение в точке касания номинальной точки «0» и характерной точки «xs» на оси абсцисс, определение значений λ1a0 и λ1axs для данной частоты вращения n

Выполненный авторами анализ показал, что потеря устойчивости системой с компрессором (ГТД или компрессор с ресивером и соплом на стенде) соответствует точке на напорной ветке, где значение изоэнтропического напора достигает максимума . Для определения положения точки "m", соответствующей максимуму теоретического напора, предложено использовать номограммы в виде семейства изолиний в пространстве координат при условии отсутствия закрутки на входе . Для удобства предложено шаг выбирать переменным, так, чтобы каждые две соседние изолинии соответствовали – одна номинальному значению , а следующая – максимальному значению . Для этого соблюдается условие =0,8 Как известно, номограммой является графическое представление функции от нескольких переменных, позволяющее с помощью простых геометрических операций исследовать функциональные зависимости без выполнения вычислений. Для данной методики использованы номограммы, которые можно описать следующим соотношением:

(4)

(5)

(6)

Номограммы строятся таким образом, что на напорной ветке, имеющей точку касания с номограммой =const, в точке касания ("m") достигается значение , которому соответствует максимальное значение , а в точке "0" пересечения с номограммой , достигается значение , которой соответствует номинальное значение поворота потока в решетке , как показано рисунке 3.

Рис. 3. Номограммы и напорная ветка в координатах для конкретной решетки профилей (решетка №15) из альбома продувок решеток профилей [4]

Выполненный анализ показал, что известный классический метод задания универсальной кривой в общем случае дает большую погрешность. Вместо этого для определения координат т. "m" предложено соотношение и как еще более общий способ – нахождение для заранее определенной т. "о" наклона линейной части на биссектрисе между этой линией и касательной, построенной из т."о" к номограмме =const. Таким образом, определяется т."m". Далее, в отличие от метода Ольштейна Л.Е. и в соответствии с атласом продувок решеток профилей [1] в окрестности т. "m" предлагается аппроксимировать напорную ветку в координатах функцией Безье, для которой в т. "о" задан луч линейной части характеристики , а в т."m" задана касательная к номограмме =const. Для этого через 2 точки (максимальную т. "m", номинальную т. "о") проводим дугу окружности в координатах и в координатах достраиваем эту дугу ветку до номограммы=const. Из этой точки "f" проводится участок вертикально вниз до точки "k", где луч из т. "x" касается одной из номограмм (для последней определяется соответствующее значение =const). Точка "х" - где линейная часть пересекается с осью абсцисс, =0. Таким образом, получается характеристика решетки профилей с учетом точки срыва и включая участок левее зоны срыва. Аналогичным образом, для конкретной частоты вращения n может быть получена экстраполированная напорная ветка в виде вместе с участком за пределами границы помпажа. Данный метод построения характеристик решеток профилей с использованием предложенных номограмм сопоставлен авторами с результатам продувок решеток профилей, опубликованных в атласе [1] и с результатами 3D CAD/CAE-моделирования, выполненного в системах UG и Ansys CFX. Степень и область адекватности признаны достаточными для использования при проектировании, доводке и диагностике состояния компрессоров в эксплуатации. Поэтому авторами решено использовать данную методику для определения номинального и максимального режима работы решеток профилей, для расчета их характеристик при создании СИМ (системы имитационного моделирования) VENEC. Такая система создается для моделирования и автоматизированного проектирования компрессоров (всех типов, включая диагональные и центробежные), а в дальнейшем – и турбин. Совместно с методами 3D CAD/CAE-моделирования ([5]), в сочетании с использованием обобщения экспериментальных данных, такие системы успешно применяются при проектировании и доводке ГТД.

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

В приведенных ниже примерах показано, как предложенные номограммы применяются для построения характеристики компрессора в целом. Для иллюстрации разработанной методики воспользуемся экспериментальными характеристиками двух компрессоров: трехступенчатого компрессора низкого давления (КНД двигателя АЛ-55) и четырехступенчатого компрессора с входным направляющим аппаратом (ВНА). На рисунке 4 показаны отрезки напорных веток и , соответствующих 30 % максимальной частоты вращения КНД АЛ-55. На рисунке 5 показана экстраполяция построенной по результатам испытаний при n=30 % от максимальной частоты вращения линейной части напорной ветки , определение точки на оси ординат, поворот вокруг этой точки касательной к напорной ветке изоэнтропического напора трехступенчатого КНД, определение в точке касания номинальной точки «0s» и на отрезке соответственно точки «0» на линии теоретического напора. На рисунке 6 показаны номограммы и напорная ветка для трехступенчатого КНД при n=30 % от максимальной частоты вращения, полученная по данной методике, в координатах .

Рис. 4. Отрезки напорных веток в виде и трехступенчатого компрессора низкого давления (КНД АЛ-55), построенных по экспериментальным результатам испытания при nпр=30 % от максимальной частоты вращения

Рис. 5. Экстраполяция для построенных по результатам испытаний при n=30 % от максимальной частоты вращения линейной части напорной ветки , определение точки на оси ординат, поворот вокруг этой точки касательной к напорной ветке изоэнтропического напора трехступенчатого КНД, определение в точке касания номинальной точки «0s» и на отрезке соответственно точки «0»

На рис. 6–8 показаны результаты построения напорных веток для трехступенчатого КНД при n=30 %, 90 %, 100 % от максимальной частоты вращения, полученные по данной методике, в координатах .

Рис. 6. Напорная ветка трехступенчатого КНД при n=30 % от максимальной частоты вращения, полученная по данной методике, в координатах и

Рис. 7. Напорная ветка трехступенчатого КНД при 90 % максимальной частоты вращения, полученная по данной методике, в координатах и

Рис. 8. Напорная ветка трехступенчатого КНД при n=100 % максимальной частоты вращения, полученная по методике, в координатах

Анализ показал, что данная методика может быть использована и для компрессоров с поворотным ВНА (входным направляющим аппаратом) и при повороте вместе с ним других НА (направляющих аппаратов) в составе компрессора. В качестве примера на рисунках 9–10 показаны номограммы и напорные ветки, построенные для четырехступенчатого компрессора с ВНА.

Рис. 9. Напорная ветка четырехступенчатого компрессора с ВНА при n=100 % максимальной частоты вращения, полученная по методике, в координатах

Рис. 10. Напорная ветка четырехступенчатого компрессора с ВНА при n=95 % максимальной частоты вращения, полученная по методике, в координатах

Выводы

Показано, что традиционные методы универсального представления характеристик решеток профилей в общем случае дают большую погрешность. Предложен новый метод, позволяющий при заданной геометрии решетки профилей построить ее характеристику, включая точку срыва и левую часть (за точкой срыва). Предложен способ определения номинального режима для каскада (компрессора). Показан способ определения на напорной ветке теоретического напора точек, соответствующих номинальному режиму и режиму максимального поворота потока (в применении к осредненной ступени и осредненной решетки профилей). Описана методика экстраполяции экспериментально полученных напорных веток на характеристике компрессора. Она позволяет по экспериментальным точкам, полученным в режиме устойчивой работы компрессора, построить границу помпажа, получить запомпажные части напорных веток, построить (в дополнение к имеющимся) дополнительные напорные ветки (вплоть до «нулевой», соответствующей началу раскрутки ротора при запуске). Приведены графические пояснения методики в общих и частных случаях. Проверка методики выполнена на примере компрессоров низкого давления ряда современных двигателей (уже находящихся в эксплуатации более пяти лет – трехступенчатого КНД АЛ-55 и четырехступенчатого компрессора с ВНА). Проведено сопоставление результатов с известными расчетными и экспериментальными характеристиками двигателей и с опубликованными результатами продувок решеток профилей, а также с результатами 3D CAD/CAE-моделирования в системах UG и Ansys CFX. Предложенная методика экстраполяции экспериментально полученных напорных веток на характеристике при малом числе испытаний предполагает определенную последовательность действий над компрессором и порядок последующей обработки результатов испытаний либо контроля параметров в эксплуатации. Она позволяет упростить и объединить экспериментальные исследования компрессоров и экстраполяцию напорных веток на характеристиках с помощью SCADA систем.

Работа выполнена при финансовой поддержке Министерства образования и науки РФ (Госзадание).

Рецензенты:

Ахмедзянов Д.А., д.т.н., профессор кафедры авиационных двигателей ФГБОУ ВПО «Уфимский государственный авиационный технический университет», г. Уфа;

Кузьмичев В.С., д.т.н., профессор кафедры теории двигателей летательных аппаратов ФГАОУ ВО «Самарский государственный аэрокосмический университет имени академика С.П. Королева (национальный исследовательский университет)», г. Самара.