Основы метода отслеживания струй Текст научной статьи по специальности «Механика и машиностроение»

!55Ы_1992-6502_(Рпп^_

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

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

УДК 621.452.3

Основы метода отслеживания струй

12 3

И. А. Кривошеев , К. Е. Рожков , Н. Б. Симонов

1 krivosh777@mail.ru, 2 rke@mail.ru, 3 sventigo@yandex.ru 1-3ФГБОУ ВПО «Уфимский государственный авиационный технический университет» (УГАТУ)

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

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

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

ВВЕДЕНИЕ

В современных методиках двумерного расчета лопаточных машин (в 2Б-расчетах, проводимых на ранних стадиях проектирования компрессоров и турбин), рассматриваются так называемые «треугольники» скоростей, которые отображают векторное сложение скоростей на входе и выходе из элементарных решеток профилей и ступеней (выделяемых по методу Жуковского Н. Е. в слоях переменной толщины вдоль поверхностей тока). Аналогичным образом векторные соотношения изображают для интегральных (осредненных по расходу или по площади) параметров на входе и выходе из лопаточных венцов. В соответствии с принципом Даламбера треугольник скоростей - это векторное уравнение, состоящее из скоростей: относительной, окружной (переносной) и абсолютной. В общем виде такое уравнение выглядит следующим образом:

с = и + #, (1)

где и - окружная скорость; - относительная скорость; с - абсолютная скорость рабочего тела.

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

ПОСТАНОВКА ЗАДАЧИ

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

ПРЕДЛАГАЕМОЕ РЕШЕНИЕ

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

по гипотезе Жуковского Н. Е. слоя вдоль поверхности тока). Эта струйка проходит через рабочее колесо, осевой зазор и направляющий аппарат в общем случае. Так происходит от ступени к ступени (через весь каскад и даже через весь двигатель). В настоящее время в литературе этот метод получил название МОС (метод отслеживания струй)

В лопаточных машинах [2] при продвижении частицы рабочего тела вдоль проточной части (ПЧ) происходит изменение термогазодинамических параметров (статические и полные значения давления, температуры в относительном и абсолютном движениях). Эти изменения связаны с изменениями направления потока. Анализируя кривизну траекторий и линий тока, (в относительном и абсолютном движении), можно определить соответствующее изменение полных давлений и температуры. Таким же образом вдоль каждой траектории и линии тока на соответствующих участках по ширине лопаточных венцов можно определить степень повышения давления, коэффициент потерь, в т.ч. полного давления для рабочих колес (РК), ступеней и многоступенчатых каскадов (компрессоров и турбин). В выделенном эквивалентном канале можно определить расход воздуха, коэффициента полезного действия и другие параметры.

Для анализа эффективности использования метода МОС построена ЗБ-САБ/САЕ модель исследуемого объекта (ступени компрессора), произведен расчет в газодинамическом пакете (Ап8у8 СБХ), далее результаты трехмерного расчета (ЗБ-САБ/САЕ) сопоставлены с результатами, полученными по 2Б-методике на основе МОС, описанной ниже.

В основе современных теоретических и экспериментальных исследований течений в тур-бомашинах лежит гипотеза Жуковского Н. Е. о том, что течение в слоях переменной толщины, выделяемой вдоль поверхностей тока мало зависит от соседних слоев. Это позволяет использовать 2Б-моделирование течений в выделяемых вдоль поверхностей тока (на различной высоте ПЧ) элементарных решетках, элементарных ступеней и элементарном многоступенчатом каскаде (компрессора, турбины). При этом используется метод экспериментального исследования течений в плоских бесконечных решетках эквивалентных профилей, предложенный Жуковским Н. Е. в 1889 г. Для получения такой плоской решетки рассечем кольцевую решетку двумя соосными цилиндрическими сечениями с радиусами г и . Развернем на плоскость полученную кольцевую решетку и увеличим количество и высоту лопаток беспре-

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

Рис. 1 Выделение в ПЧ ступени компрессора слоя переменной толщины по методу Жуковского Н.Е.: а - без учета искривления поверхностей тока; б - с учетом искривления и условий радиального равновесия

На рис. 2, а показана выделенная элементарная решетка, в которой происходит изменение параметров: изображены прямые неизогнутые профили (пластины). Их толщину можно считать бесконечно малой, по сравнению с другими размерными величинами. Это соответствует гипотезе Жуковского Н. Е. На рис. 2, б выделена плоская решетка уже в реальном компрессоре, между двумя, близко расположенными поверхностями тока.

Рис. 2 Общий вид плоской элементарной решетки (а) для выделенного слоя переменной толщины (б)

Проведены 3D-CAD/CAE расчеты ступеней компрессоров с ВНА и без ВНА. На рис. 3, а, б показаны линии и поверхности тока, непрерывные траектории движения частиц в абсолютном и относительном движении в ступени компрессора, на рис. 3, в, г показано течение в ступени компрессора без ВНА в относительном и абсолютном движении, на рис. 3, д, е показано течение в ступени компрессора с ВНА в относительном и абсолютном движении.

Для проведения дальнейшего расчета получены результаты 3D-CAD/CAE-моделирования компрессора низкого давления, получены картины течения в абсолютном и относительном движении по длине компрессора вблизи втулки (рис. 4).

В настоящее время при 3D-CAD/CAE-моделировании [3] в середине осевого зазора (ОЗ) между РК и НА применяется скачкообразное переключение «перескок» из относительной в абсолютную систему координат, задавая метод перескока (Frozen rotor или Stage). При этом полагают, что окружная скорость продлевается за срез решетки до середины осевого за-

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

Для ликвидации этих недостатков и этих вопросов авторами предлагается метод отслеживания струй (МОС). Суть метода - это одновременное рассмотрение траекторий и линий тока для течения рабочего тела вдоль лопаточной машины в рамках выделенного слоя переменной толщины. Однако в развитие метода Жуковского Н. Е. предложено прослеживать непрерывное движение одновременно в абсолютной и в относительной системе координат, сочетая методы Эйлера и Лагранжа.

Рассмотрим предложенный метод на примере. На развертке слоя переменной толщины изобразим траектории движения частицы рабочего тела одновременно в относительной и абсолютной системах координат. Для этого в цилиндрической системе координат вдоль оси двигателя выделяем координату £ - аппликату (для осевых элементарных решеток, ступеней и каскадов - вдоль оси двигателя, для диагональных удобнее использовать координату вдоль поверхности тока в меридиональной плоскости). Выделяем также х„ - тангенциальная (окружная) координата в относительном движении, х - тангенциальная (окружная) координата в абсолютном движении (рис. 5).

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

В НА

О'- РК < эз НА

1

ш

д

Рис. 3 Линии и поверхности тока в ступени компрессора: а - без ВНА; б - с ВНА; в - течение в ступени компрессора без ВНА в относительном движении; г - течение в ступени компрессора без ВНА в абсолютном движении; д - течение в ступени компрессора с ВНА в относительном движении; е - течение в ступени компрессора с ВНА

в абсолютном движении

в

г

е

Рис. 4 Течение (траектории и линии тока) в компрессоре низкого давления

(вдоль втулки)

В традиционных методах расчета - с использованием треугольников скоростей (в 2Б-моделях) все параметры записывались в конечных разностях типа :

Шг = Д(и • си). (2)

Рис. 5 Система координат, лопаточные венцы, и траектория движения частицы

В отличие от этого, в методе МОС предлагается перейти к дифференциальным уравнениям и интегральным соотношениям типа

ДЯг| 2 = г с • с1и +С2 и • Сс .

г с и «к,,. и

(4)

1^1 '"Ч ** ""и

Рассмотрение изменения са(Б) позволяет от фактора времени . перейти к анализу кривизны линий тока и траекторий, связать их с изменением Г*, Г,Д Р*, Рк*, окружных сил, действующих на газ (на частицу) и на лопатки. Действительно:

с£

Са =

ск

, СБ С. =—.

Са

С учетом этого:

Сх Сх

Си т Са ,

ск йБ л л Сх Схл

ас = ас--+ са-й (—),

и а йБ а СБ

СНГ 2 2 Сх Сг с Сх Сс

-г = г а (------ +---а--

СБ СБ СБ г СБ СБ

_ сС 2 х + ^),

(5)

(6)

(7)

(8)

СНг = а (и • си ) = си

- и • Сс„

(3)

т-т г 2 г,dx dr ca dx dca

HT = Jr2-w (------ +---- +

dS dS r dS dS

d x.

(10)

+ c

a dS2

)dS.

Для определения параметров вдоль линии тока воспользуемся дифференциальным соотношениям для определения величины теоретического напора лопаточной машины:

dHT = cu ■ du + u - dcu =

dx dx dx (11)

= du - c---+ u[dc---+ c ■ d(—)].

a dS a dS a dS

В случае использования допущения

C _ const получаем упрощенные соотношения

* *

C dT _ dI dHT p dS dS dS

d2 x

du dx dca dx

= ca----+ u[—---+ ca---

dS dS dS dS dS

dT r2-w2

(12)

dS Cp

(dx dr ~ca + dx dca + _ d2 x)

dS dS r dS dS a dS2 Cp '

(13)

T =J

i 2 2ч

(r -W )

Cp

d2 x.

(14)

dx dr c dx dc „ „

- (-----^+---^+ca—2)dS.

dS dS r dS dS a dS2

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

P * т* _к_

ид _ (L-\k-1 P* -\T*> '

r0 10

T* —

p * _ p*(-)k 1

''ИД — r0\^* J ■

T

(15)

(16)

Для получения на основе кривых идеального давления Р*ид (.5) реальных значений давления Р* (.5) необходимо при расчете учесть потери полного давления, связанные с геометрией лопаточной машины и режимом течения:

Р* = Рид * "О. (17)

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

О = П О,. =п (ОрК „'ООз'ОНА ), . (18)

п=1 п=1

При этом предлагается в осевом зазоре вычислить потери полного давления в относительном и абсолютном движениях и результат ос-реднить. В рабочем колесе вычисляем потери полного давления в относительном движении, в направляющем аппарате - в абсолютном движении. Можно использовать эмпирические результаты экспериментальных исследований, продуть или в ЗВ-СЛБ/СЛБ смоделировать

движение в эквивалентных каналах.

Для построения траектории в абсолютном и относительном движении воспользуемся зависимостью:

, , , dS dS

dx + dx„, = u-dt = u--= =-,

c c

a a

x (S) + xw (S) = J dS.

(19)

(20)

Для значений в относительном движении справедливо:

dxw

Wu =-Т , (21)

dt

dxw dxw

wu = —- = ca---

u dt a dS

dH

dS

Tw 2 2 dx dr c dx dc = r - о (-

a , w a

dS dS r dS dS

d2 xw

(22)

(23)

a dS2

. 2 2 dxw dr ca dxw dca

ЯГ,-2 w a . w a

Tw = J r -w (------1-----

dS dS r dS dS

- d^

"a dS2

-)dS,

C = const,

m' 2 2

dT w r -W

dS C

d2 x

(+dx^.dc^+c -- -w)

dS dS r dS dS a dS2 dT dI dHT

C _w__w = Tw =

p dS dS dS

(24)

(25)

(26)

du dxw dca dxw = ca----- + u[—^---

dS dS dS dS

d2 x2 V~dS~

a

г* =, (г2 •а ) Сх, Сг са

№ Ср СБ ' СБ ' г

сСх„, Сса С х

(28)

, а

+ с

СБ СБ а СБ2

)СБ.

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

• в окружном направлении:

- Сси

К =-Т, (29)

Ж

где К - удельная сила, действующая на частицу в окружном направлении;

• полная (интегральная) сила действующая на лопатку, рассчитывается как:

К = и2 • с (

,Сх Сс„ _ а хч

■ + са--2) =

СБ СБ а СБ2

= с

сС2 х

а1Б'

(30)

• в осевом направлении:

— йса Сса

а 1 са ,

йг йБ

• сила, действующая на участке СБ

, 2 _ Сх Сса _ а2хч , ,

= [и •са — + са ■ -¿Г^РСа^ аБ аБ аБ

* (са2 С7Xx')Рcadгds, СБ

(31)

(32)

• интегральная сила, действующая по высо-

те:

2 _ сСх Сса _ С х.

К (г) = Л [и 2 • са + •са -Т)]рсаСг \ ds I СБ СБ СБ2 I

• I 2 С х I

• полная интегральная сила:

, ,г 2 - ,Сх Сса _ С2х

К = Ц [и (---- + са--2)]рсаСгСя '.

т аКСБ СБ а СБ2 а

(33)

!П(с2 С^х)РСа СБ

(34)

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

Ссг 2 = Рс2

' С2 г >

ч СБ2 у

(35)

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

Ссг Сг

с2 СР + р— = — ^

г Сг

С2 г

СБ2

1

Сх СБ

Сг

(36)

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

* *

р - Р2 1-о

(37)

£ =-

Р

1 - л(Х)

Отсюда искомое значение коэффициента восстановления полного давления

и = 1 - ^[1 - л(Л)], где значение Я должно подставляться для входа в компрессорную решетку (и для выхода - из турбинной решетки) профилей. При этом для расчета местных сопротивлений необходимо учитывать диффузорность каналов. Диффузорность каналов можно рассчитывать по-разному, в т.ч. по Либляйну с помощью формулы из [2]:

В = 1 + 0,5 хЬ / г х

А,

". (38)

К1 "1 В соответствии с расчетными данными по графику, представленному в работе [5] можно определить величину потерь в соответствии с конфигурацией межлопаточного канала. Кроме того, величину потерь можно определить так, как описано в [3, 5], при продувке всей ступени по готовой геометрии. Такой метод определения потерь будет наиболее удобным, т.к. не требует никаких математических описаний. Но для моделирования течения в ступени или каскаде необходимо иметь все геометрические параметры ступени или каскада.

ВЫВОДЫ

Таким образом, предложен достаточно простой и эффективный метод моделирования и проектирования лопаточных машин (компрессоров и турбин) в составе ГТД. Предлагаемый оригинальный метод отслеживания струй (МОС), в отличие от использования треугольников скоростей, позволяет проследить непре-

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

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

1. Кривошеев И. А. Использование имитационного моделирования и поддержки принятия решений при системном проектировании авиационного двигателя в составе летательного аппарата. Авиационно-космическая техника и технология. Харьков: «ХАИ», 2007. № 6(42) . С. 63 -81. [I. A. Krivosheev. "Using simulation and decision maling support during system designing aircraft engine as part of the aircraft". (in Russian) in Aerospace technic and technology. Kharkiv: «KhAI», 2007, № 6(42) . pp. 63-81.]

2. Холщевников К. В., Емин О. Н., Митрохин В. Т. Теория и расчет авиационных лопаточных машин. М.: Машиностроение. 1986. 432 c. [K. V. Holschevnikov, O. N. Emin, and V. T. Mitrokhin, Theory and Design Aviation Turbomachinery, (in Russian). Moscow: Mashinostroenie, 1986.]

3. Щелковский М. Ю. Верификация программного комплекса ANSYS CFX для численного анализа трехмерного вязкого течения в компрессоре. Восточно-Европейский журнал передовых технологий. 3/10(57) 2012. С. 60 - 65. [M. Y. Shchelkovsky. "Verification software package ANSYS CFX for numerical analysis of three-dimensional viscous flow in the compressor".(in Russian) in Easter-European Journal of Enterprise Technologies. 3/10(57) 2012. pp. 60-65.]

4. Щелковский М. Ю. Параметрическое исследование газодинамических характеристик компрессорных решеток. Авиационно-космическая техника и технология, Харьков: «ХАИ», 2012, № 7(94) . С. 95 - 106. [M. Y. Shchelkovsky. "Parametric study of gas dynamic characteristics of compressor grids".(in Russian) in Aerospace technic and technology. Kharkiv: «KhAI», 2012, № 7(94) . pp. 95-106.]

5. Кривошеев И. А., Симонов Н. Б. Сравнение результатов моделирования с экспериментальными данными о продувках решеток плоских неизогнутых пластин. Молодежный Вестник УГАТУ. № 4 (9), 2013. С. 17-23. [I. A. Krivosheev, N. B. Simonov. "Comparison of simulation results with experimental data of flat undeflected plates lattices expulsions".(in Russian) in Youth Vestnik UGATU, № 4 (9) / 2013. pp. 17-23.]

ОБ АВТОРАХ

КРИВОШЕЕВ Игорь Александрович, проф. каф. авиац. двигателей, науч. рук. НИЛ САПР-Д. Дипл. инж.-мех. (УАИ, 1976). Д-р техн. наук по тепл. двиг. ЛА (СГАУ, 2000).

Иссл. в обл. инф. технологий в двигателестроении. РОЖКОВ Кирилл Евгеньевич, асп. каф. авиац. двигателей. Дипл. инж.-мех. (УГАТУ, 2009). Готовит дис. в обл. моде-лир. раб. процессов лопаточных машин. СИМОНОВ Николай Борисович, асп. каф. авиац. двигателей. Маг.-инж. (УГАТУ, 2014). Готовит дис. в обл. моделир. раб. процессов лопаточных машин.

METADATA

Title: Streamline tracking method basics.

Authors: I. A. Krivosheev1, K. E. Rojkov2, N. B. Simonov3.

Affiliation:

1-3Ufa State Aviation Technical University (UGATU), Russia. Email: 1krivosh777@mail.ru, 2rke@mail.ru,

3sventigo@yandex.ru. Language: Russian.

Source: Vestnik UGATU, vol. 19, no. 3 (69), pp. 103-110, 2015.

ISSN 2225-2789 (Online), ISSN 1992-6502 (Print). Abstract: Turbomachinery (compressors and turbines) modeling and designing have been considered as part of GTE. Original streamline tracking method (STO) have been proposed. Unlike using of velocity triangles, STM allows to trace continuous (along the flow surfaces) changing in kinematic thermodynamic parameters of working fluid, changing in circumferential forces and torques, acting on profiles and working fluid particles, under conditions of shock-free inlet and outlet of the flow through profiles lattices to create the flow part shape, interscapular channels and taking it into account to profile blades Key words: turbomachinery, compressors, turbines, profiles lattices, flow surface, streamline, trajectory, temperature, pressure, velocity, working fluid. About authors:

KRIVOSHEEV, Igor Aleksandrovich, prof. Aviation Engines Dept. Dipl. engineer (USATU 1976). Dr. Tech. Sciences (USATU, 2000).

ROJKOV, Kirill Evgenevich, Dipl. engineer (USATU, 2009). Graduate student.

SIMONOV, Nikolai Borisovich, Mag.-eng. (USATU, 2014). Graduate student.