Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Осевые и центробежные компрессоры двигателей летательных аппаратов. Теория, конструкция и расчет

.pdf
Скачиваний:
42
Добавлен:
12.11.2023
Размер:
12.67 Mб
Скачать

В книге К.В. Холщевникова* дается вывод уравнения для случая конической ступени компрессора с криволинейными поверхностями тока (при dcr=0). Здесь в учебных целях рассматривается упрощен­ ная модель течения в ступени: течение воздуха осуществляется по цилиндрическим линиям тока, т.е. параллельно оси компрессора; равновесие рассматривается в осевом зазоре, и, следовательно, отсут­ ствуют силы от воздействия лопаток на воздух; напор и потери оста­ ются постоянными по радиусу проточной части ступени. На рис. 5.1,а приводится принятая схема течения воздуха в ступени. Выделим в по­ токе воздуха в пределах осевого зазора элементарный объем (рис. 5.1,6) и рассмотрим схему воздействия на элементарную частицу воз­ духа сил гидродинамических давлений и центробежной силы на про­ извольном радиусе г (рис. 5.1,в). Размеры частицы: толщина dr, протя­ женность по окружности rd(p, протяженность по оси ступени 1, масса частицы dm = р • г • dip dr • 1, где р — плотность воздуха.

На нижнюю грань частицы действует гидродинамическое давле­ ние р и соответствующая сила dP\ = pdf\, где df\ = г • rfcp • 1 — площадь

нижней грани частицы. На верхнюю грань действует гидродинамичес-

Рис. 5.1. Расчетная схема течения воздуха в ступени (а, схема (в) и диаграмма (г) сил, действующих на элементарный объем воздуха в осевом зазоре

*См. сноску на стр. 38.

кое давление p +dp и соответствующая сила dP2 = (p +dp)df2 , где df2 = (г + dr) dip • 1 — площадь верхней грани частицы. На боковые гра-

 

dp

ни действует гидродинамическое давление Р + ~тг и соответствующие

 

z

силы dP ^-dP ^-

\

dfo, где dfy = df4 = dr \ — площадь боковых

V J

граней частицы. Кроме того, из-за наличия окружной составляющей скорости си на частицу действует центробежная сила de.

На рис. 5.1,г показана диаграмма сил, участвующих в рассматри­ ваемом явлении. Очевидно, для того чтобы выполнилось требование сг= 0, сумма проекций всех сил на радиальное направление должна

быть равной нулю. Имеем

de + dP\ - dP2+ 2dP3sin dip

= 0 .

(5.1)

2

 

 

Это и есть уравнение радиального равновесия в общем виде.

Учитывая, что de = dm

си

уцс, где уцс = — — центростремительное ус­

корение, подставляя значение сил в (5.1), отбрасывая бесконечно ма-

d(D

самим углом

do

лые высших порядков и заменяя величину sin1-г12

^

радианах), получаем окончательно

 

 

 

f - p - = o -

 

(5.2)

dr r г

 

 

 

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

5.3. УРАВНЕНИЕ СВЯЗИ ОКРУЖНЫХ Си И ОСЕВЫХ Са СКОРОСТЕЙ В РАДИАЛЬНО УРАВНОВЕШЕННОМ ПОТОКЕ

В общем случае на любом радиусе ступени имеем

с2=с% +с£ +с?.

При cr =О

С2 = Сд + Сц.

(5.3)

Изменение по радиусу ступени окружной скорости си предопреде­

ляется уравнением (5.2). Возникает вопрос: а как при этом должна из­ меняться по радиусу ступени осевая скорость са . Ответ на этот вопрос

дает искомое уравнение связи.

Обобщенное уравнение Бернулли (см. разд. 1.3) в дифференци­

альной форме имеет вид

 

\

 

( л

 

dH, = d

с_

+& +dLR .

\

2

J Р

*

Строго говоря, оно справедливо только вдоль струи тока или вдоль поверхности тока. Однако если допустить, что изменение состо­ яния воздуха вдоль всех линий тока протекает с постоянным показа­ телем политропы п =const, т.е. Я2(г) = const, то уравнение Бернулли

можно продифференцировать по радиусу. Не учитывая силы трения и принимая во внимание соотношение (5.3), получаем

dHz

\ dp

1

dcfu dcj

(5.4)

dr

р dr

2

dr + dr

 

Для простоты будем считать, что напор Hz постоянен по высоте

dHz

проточной части. Тогда— =0, и из соамостного решении (5.2) и (5.4)

окончательно имеем

dc„ dc„

(5.5)

d r +Ca d r ~ ° -

dHz

Для случая - г - * 0 вместо (5.5) получается уравнение

аг

dHz

i

1 d(cu r)2

d £

(5.6)

dr

~ 2

j2 dr

+ dr

 

Уравнения связи (5.5) и (5.6) позволяют рассчитать распределение осевых скоростей са по радиусу проточной части в осевых зазорах ци­ линдрической ступени компрессора при радиально уравновешенном потоке.

Таким образом, уравнение радиального равновесия (5.2) совместно с уравнениями (5.5) или (5.6) служат базой для расчета окружных и осевых скоростей, а значит, и всех других параметров ступени: Н2,

Рст » » К 2» ^Р» b/t

5.4.СПОСОБЫ ЗАКРУТКИ ЛОПАТОК

ВСТУПЕНИ ОСЕВОГО КОМПРЕССОРА

Анализ уравнений (5.5) и (5.6) показывает, что если задаться зако­ ном изменения cju и по радиусу проточной части ступени компрес­

сора, то соответственно получим Cjfl=/i(r) и

А так как ок­

ружная скорость зависит от радиуса (и = (йг)9 то однозначно опреде­ ляются и все остальные параметры треугольников скоростей, а зна­ чит, и законы изменения по радиусу таких величин, как Hz , рст, ДР, b /t и др.*

Закрутка лопаток по закону постоянства циркуляции. Для по­ лучения небольших потерь в ступени целесообразно, чтобы слои воз­ духа по высоте проточной части в осевом направлении не «проскаль­ зывали» друг относительно друга, т.е. чтобы ca(r) =const. Для этого

случая

dca

= 0, и из уравнения (5.5) следует, что

4dc„

-+ cu - ^ = ° ,

ИЛИ

 

 

г

си = 0 .

 

Решая последнее уравнение, имеем

 

lnr

ctt = const,

 

или окончательно

 

 

r c u =const.

(5.7)

*

Eckert В. Axialkompressoren and radialkompressoren. — Berlin: Springer-Verlag, 1953.

Уравнение (5.7) определяет способ закрутки лопаток, именуемый законом постоянства циркуляции.

Применяя его к контрольному сечению 1-1 перед рабочим коле­ сом, получаем С\а(г) = const, г -с 1м = const. В процессе выводов принято

допущенйе, что Ят(г)* const. Учитывая, что при цилиндрических ли­ ниях тока в контрольных сечениях и\ - и2 и гх= г2, а Ят = ©г (с^ - с1и), при гх • с1ы= const ДОЛЖНО быть и 7*2 *С2и = Const.

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

Из уравнения (5.7) легко видеть, что при этом законе С\и и уменьшаются по сравнению с с1мср и С2нср к периферии (при г>гср) и увеличиваются к втулке (при г<гср) (рис. 5.2).

Рассмотрим изменение других параметров по радиусу ступени при этом способе закрутки. Как отмечалось в разд. 3.5, степень реактив­

ности ступени

 

 

 

 

 

c j - c f

(5.8)

 

Рст 2Нт

2Нт

 

 

 

 

Для случая с\а =const и с^ - const

 

вместо последнего равенства можно за­

 

писать

 

 

 

 

Рст~ 1

c j j - c f u

 

 

 

2Нт

 

 

 

Отсюда следует, что при законе закрут­

 

ки г • си = const с учетом Ят(г) = const

 

 

Рис. 5.2. Изменение а и и съл по

 

 

 

 

высоте проточной части

 

Рст~ ^

const

(5.9)

 

?

 

 

 

В соответствии с (5.9) Рст>Рст.ср при г> Гср и рст<рсг.ср при г< г^ (рис. 5.3).

Такое изменение степени реактивности рст(г) приводит к возраста­ нию относительной скорости w ^w ^w ^p) в верхней половине рабо­ чих лопаток и абсолютной скорости С222ср) в нижней половине. Это, в свою очередь, при относительно длинных лопатках может при­ вести к MWi > кр в зоне 1 и MCz > Мс2кР в зоне 2 (рис. 5.4). Последнее недопустимо из-за возможности срывного характера течения в соот­

ветствующих сечениях ступени. Особенно сильно меняется рст по ра­ диусу и соответственно W\ и с2 в ступенях с малыми rfBT, т.е. при от-

Лл

носительно длинных лопатках ——. ^ср

Рст.ср

Рст

 

Рис. 5.3. Изменение pet по радиу­

Рис. 5.4. Схема срывного тече­

су г ступени

 

ния в ступени:

 

 

1 — зона Мщ >Мцг| кр j

2 — зона Мс2 >Мс2 Ч»

Поэтому из-за того, что в первых ступенях осевых компрессоров обычно dBT = 0,35—0,45 и длина лопаток большая, способ закрутки rcu = const, са(г) = const здесь не применим. Такой способ целесообраз­ но применять в средних и последних ступенях компрессора, где dBT> —0,6.

Способ закрутки с постоянной реактивностью по радиусу. При выборе способа закрутки лопаток удобно задаваться изменением сте­ пени реактивности по радиусу в форме

В

(5.10)

Рст” 1 г т + 1 *

где В и т — постоянные. Отсюда, например, при Ш= 1 получается

_В_

Рст “ 1 9 »

т.е. способ закрутки rctt = const. При m = - 1 из (5.10) следует, что

Рст = 1 - ^ = 1 - 5 = const.

(5.11)

Такой способ закрутки принято называть способом «с постоянной по радиусу реактивностью».

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

целей уравнением (3.13), справедливым не только для гср, но и для любого текущего радиуса г. Разрешая его относительно с\и , получаем

= 1 - Р с т ) - § .

(5.12)

Учитывая, что по условию Ят = Ят ср=и(с2й - с 1и) = const и, следова­

тельно,

Ят

с 2и ~ с \и + ~ »

найдем также, что

С2и = ы(1-Рст) + ^ -

(5.13)

Значения осевых составляющих скоростей С\а(г) и с^ г) при закрут-

ке рст= const не будут постоянными по радиусу. Для их определения необходимо воспользоваться уравнением связи, например, в форме (5.5), если Ят(г) = const.

dc\и

Для контрольного сечения 1-1 С\и и

находятся по уравнению

(5.12), и дифференциальное уравнение (5.5) для определения с\а(г) принимает вид

d(c2la)

Нт ц2

,

(5.14)

- ^ -

= 2 (1 -р ст) - ^ - 4 ^ (

1 - р ст)2

После его интегрирования с учетом того, что и = саг, a константа интегрирования находится по известному значению с1а(гср), получим

с 1а = У <?1аср- 2(1 - Рст)V - М2р) + 2(1 - рст) Ятср In j -

(5.15)

гср

Таким же образом для осевой составляющей скорости за рабочим колесом (в сечении 2-2) найдем

С2а= УС2а ср- 2(1 - Рст)2^ 2- Ы2р) - 2(1 - рст) Ятср In -f-

(5.16)

гср

 

Анализ формул (5.15) и (5.16) показывает, что при закрутке лопа­ ток по способу pCT = const, # T = constc увеличением радиуса осевые скорости С\а и С2а уменьшаются, притом тем сильнее, чем длиннее ло­

патки. При dBT < 0,5 и с а{г) Ф const появляются дополнительные потери энергии трения из-за проскальзывания слоев по радиусу. Вместе с тем для рассматриваемого способа закрутки лопаток рст= const характерны малые изменения по радиусу чисел Маха MWj и МСг, что гарантирует

отсутствие критических чисел Маха по высоте проточной части, если на среднем диаметре Mw cp<Mw ,ф и МС2Ср<МС2Кр. Таким образом,

способ закрутки рст= const особо подходит для использования в пер­ вых и средних ступенях компрессора, где dBT< 0,5—0,6.

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

cunrm = const,

(5.17)

где

с1ц+с2ц

2

Здесь т — тот же самый показатель степени, что и в формуле (5.10). При разных его численных значениях в пределах от + 1 до - 1 получа­ ют различные способы закрутки лопаток по радиусу: так, при т = 1,0

— ступень с постоянной по радиусу циркуляцией, т.е. rcu =const, при т = - 1,0 — ступень с постоянной по радиусу степенью реактивности рст= const.

Как показывает анализ результатов расчетов и опыт проектирова­ ния, в первых ступенях каскада низкого давления целесообразно ис­ пользовать промежуточный закон закрутки, когда показателю степени т в формулах (5.17) и (5.10) придают промежуточные значения 1,0>т > - 1,0. В этом случае осевые составляющие абсолютной скоро­ сти потока воздуха на входе в рабочее колесо в различных сечениях по радиусу определяются по формуле

c la = V H ^ Cle c p

m - 1 Дт.ср ..

m +1 с,2

1

‘'la ср

 

i ^

Рст.ср)

 

Р

-

т

с 2

- 5

'с р

 

 

с 1а ср

Г

(5.18)

 

 

 

 

 

 

- ( m +

1) "

 

PcTcpJ

Г . - f i

l

_

 

 

 

 

 

 

 

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

- & —

У , + ”Lz 1 0

 

Рст.ср)

 

1

-

с 2а ср

 

тт

c l ср

_

/

. т t

^т.ср

ч

 

г

 

1 - ( т

 

1 - ——

 

 

m + lc ,2

(1_Рстср)

 

_

u.v

 

 

 

с 1а ср

 

 

 

 

 

 

 

- 2 т

 

г ср +

с р ,

_

(5.19)

+ 1)"

_

- Г

-

гср

где г = 7 и гср= 7 ^ — относительные радиусы расчетных сечении; г —

' к

 

г к

текущий расчетный радиус.

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

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

минимальные потери и допустимые числа Маха

и МС2.

Отметим также, что условие Ят(г) = const, принятое при выводе рас­

четных формул для длинных лопаток (dBT<0,35), что характерно в случае вентиляторных ступеней, вносит заметные погрешности в рас­ четы пространственного потока. Последнее объясняется трудностью соблюдения этого условия в периферийных и корневых сечениях про­ точной части. Поэтому при расчете подобных ртупеней задают харак­ тер изменения напора по радиусу НТ(г)Ф const. При этом в лопатках вентилятора напорность у втулки существенно уменьшается.

5.5. РАСЧЕТ НА ЭВМ КИНЕМАТИЧЕСКИХ ПАРАМЕТРОВ ПОТОКА ПО ВЫСОТЕ ПРОТОЧНОЙ ЧАСТИ СТУПЕНИ КОМПРЕССОРА

Поясним, как практически пользуются законами закрутки лопаток. Ис­ ходными данными являются результаты численного расчета осевого компрес­ сора на_среднем диаметре (см. разд. 4.5). Поэтому считается, что параметры

^т.ср» ^т.ср»

ср » Рст.ср» ^li/cp» ^lacp» ^2аср» ^2иср» P lcp 1 fecp» (^Оср H®

среднем диаметре ступени известны.

Взадачу расчета входит определение кинематических параметров потока

внескольких (обычно 5—7) сечениях по высоте проточной части. Для приме­ ра расчет проведем в трех сечениях (рис. 5.5).

Ш

^ г3 щ0,5935(f3-0,979}

гг ~гср 9б55(гг =0,768)

Г, =0,2827(г , - 0,955}

Рис. 5.5. Схема расчетных сечений по радиусу лопатки

Расчет начинается с выделения расчетных сечений по высоте проточной части, т.е. радиусов Г |, г2 , ...» гп . Расчетные сечения в корневой и перифе­

рийной частях (см. гj и г3 на рис. 5.5) выбирают на расстоянии 2—4 мм соот­

ветственно от втулочного и наружного радиусов лопатки на выходе из решет­ ки. В соответствии с высказанными в разд. 5.4 рекомендациями выбирается закон закрутки.

Впрактике проектирования компрессоров расчеты параметров потока по радиусу проводятся для всех ступеней. С учетом большого числа ступеней и расчетных сечений по высоте каждой ступени расчеты выполняются на ЭВМ. Ниже приводится алгоритм расчета параметров потока по радиусу проточной части ступени.

Втабл. 5.1 дан пример расчета параметров потока по радиусу первых сту­ пеней КНД и КВД (значения в скобках). При этом для ступени КНД был вы­

бран промежуточный закон закрутки лопатки с показателем ступени т = - 0 ,5 , а для ступени КВД закон закрутки рст= const (m = - 1). Исходные данные на

среднем радиусе для расчета закрутки лопаток по высоте взяты из примера расчета компрессора, приведенного в разд. 8 .2 . Относительные радиусы (г] =0,466, Г2 = гср = 0,768; г3 = 0,9790 определяются как отношение радиуса в

расчетном сечении (rj = 0,2827 м; г2 = гср = 0,4655 м; г3 = 0,5935 м) к периферий­ ному радиусу на входе в рабочее колесо (гк = 0,6065 м) (см. рис. 5 .5 ).