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

А. В. Бараненко. Холодильные машины

.pdf
Скачиваний:
522
Добавлен:
10.02.2015
Размер:
32.49 Mб
Скачать

гидравлического кпд (9.30), который в этом случае представля­

ют так:

~

 

lэф

 

] .2

(

)

l1ПOJ1· (9.207)

11" =Т = [

 

<Р1и

1 + ~пр + ~тp

 

а

<Р2и(п2СР/ Цср) -

и1СР

 

 

 

Определять по ре3ультатам специальных исследований потери

на протечки и трение для осевых компрессоров достаточно труд­

но, особенно если принять во внимание. что в подавляющем боль­ шинстве компрессоров РК не имеют бандажа. Это. приводит

К тому, что разделить потери на протечки и трение практически

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

потерь на протечки путем введения дополнительных коэффици­

еитов к потерям на трение.

Мощность Трения ротора осевой ступени дисковой конструк­ ции (в киловаттах) находят по формуле Стодолы [74. 89]

(9.208)

Здесь D, ивт - диаметр втулки РК, м, и окружная скорость на

ее поверхности, м/с; р = О,5(Р! + Р2) - средняя плотность газа

по обе стороны диска, кг/м3; Ктр - опытный коэффициент; для

воздуха Ктр = 1,07.

Косвенный учет протечек в практике расчетов осуществляют

введением дополнительного коэффициента Кор =1,3 + 2,0, и тог­

да суммариые потери на протечки и трение определяют по фор­

муле, полученной на основе зависимости (9.198),

R.

+ R.

=

КщlVтр

= КпрNтр (9209)

I-'пр

I-'тp

[(

)] 2

Gl··

 

 

G Ч>2и

D2cp/D1cp - <Р1и и1ср

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

эффициент Кпр.тр'; КпрКтр' который изменяется в пределах

Кпртр = 1,4 + 2,2.

Потери на протечки и трение в раСчетах авиационных осевых

компрессоров [89] учитывают с помощью коэффициентов, по форме

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

(9.210)

Здесь la - эйлерова работа ступени (9.202);Zop = (!1Gпр/G)Ia -

удельная работа, дополнительно затрачиваемая на сжатие пере­

текающего газа !1Gx; ~ = N тp/G - удельная работа, дополнитель­

но затрачиваемая на преодоление сил трения диска о газ; 1'\" -

гидравлический КПД (9.207); 1'\пр - коэффициент, учитываю­ щий потери на протечки; 1'\тр - коэффициент, учитывающий поте­

ри на трение. (В работе [89] 1'\пр обозначен как 1'\3 и 1'\тр - как 1'\(.) Обычно при расчетах осевых компрессоров принимают сразу

произведение, которое представляет собой

(9.211)

 

- величи~у, обратную 1 + J3ПР + Ртр' где f3пр и ~TP имеют тот же

смысл и определяются также по формулам (9.19) и (9.198). Под­ ставив зависимость (9.211) в (9.210), приходим к выражению

для гидравлического КПД (9.207).

у ступеней осевых компрессоров 1 + ~ПР + Ртр = 1,01 + 1,04, что

соответствует llпрllтр = 0,99 + 0,96. Меньшие значенияl + Рпр + ~TP

(БОльшие з~ачения 1'\пр1'\тJ соответствуют первым ступеням с ма­

лыми Vl' БОльшие - последним ступеням ~ большими v1При

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

линейное изменение 1 + ~op + ~тp (или 1'\пр1'\тJ по ступеням.

Теоретическую (эйлерову) работу элементарной ступени за­

писывают с учетом того, чтоиi = и2 = u = idem в таком виде:

(9.212)

Она совпадает с внутренней работой 1з = 1, , так как в элемен:

тарной ступени нет потерь на протечки и трение РПР = Ртр = О

По этой же причине внутренний (политропный) КПД элементар­

ной ступени в соответствии с уравнением (9.207) совпадает с гид­

р8В:1Iическим кпд. ЭТО обстоятельство уже было отражено в за­ писи зависимости (9.185) для коэффиnиента эффективной работы элементарной ступени '1'.

Коэффициент реактивности ступени осевого компрессора оп­

ределяют так же, как и ДЛЯ центробежного:

n =i2 - i1

=i2 -11

 

1-2

1; -1;

i; -t; .

(9.218)

Так как

616

611

то

а1-2

=1-

c~ -cf

 

 

2

(.*

.*) =

 

 

 

 

Z2 -

Zl

 

 

=1 _

222

 

2

 

с2и

+ C2z - С1и -

C1z

 

2 (1 + Рпр + Ртр) [<I>2и(D2ср/D1сР)..... <l>1и] ufcp

 

2

2

 

2

 

2

 

= 1 _ <l>2и + <l>2z

- <l>1и

-

<l>lz

(9.214)

 

 

 

 

 

 

для 3JlUU!HmapнoucmyneHU, УI«YЮрой п1 = п2 = D, Рпр + Ртр = О

И C1z = C2z = cz ' коэффициент реактивности записывают в более

простом виде

a~~2 = 1- С22и -

С12и

=1- <l>22 и

-

<l>12и

(9.215)

2(С2и -

с) и

2 (<I>2u

-

<l>1и)

 

Из треугольников скоростей элементарной ступени (рис. 9.32, б,

9.33) следует, что

(9.216)

и

(9.217)

Используя эти соотношения, можно представить коэффициент

реактивности элементарной ступени в таком виде:

а(К)

= 1- С2и + с= 1- <l>2и + <l>1и

1-2

2

(9.218)

Как следует из полученных зависимостей, коэффициент реак­

тивности такой элементарной ступени полностью определяется ее треугольниками скоростей. Поэтому его называют "инематичес­

"им "оэффициентом реа"тивности. Мы будем обозначать его

с верхним 'индексом .к•.

Различают элементарные ступени с коэффициентами реактив-

носtи a~~2 = О; 0,5 и 1,0 (см. рис. 9.33).

,

В ступени ca~~2 =,0 (рис. 9.33, а) в РК относительная ско­

рость изменяется только по направлению, оставаясь неизменной

618

Рис. 9.33. Схемы лопаточных аппаратов и треуroJIЬRИКИ скоростей эле­

мевтарных осевых ступев:ей при различных зв:ачев:иях кив:ематического

коэффициеl;lТа реактивв:ости: а - a~~2 = О; б - a~~2" 0,5; в - a~~2" 1

по значеНИЮW2 = w 1 + f).w и W 1 = W 2 Из-за этого В РК измене­

ния энтальпии газа не происходит и t2 = t1 ' что В соответствии

с формулой (9.213) и обусловливает равенство нулю коэффициен­

та реактивности. Весь теплоперепад, т. е. все увеличение энталь­ пии газа, срабатывается в НА за счет уменьшения в нем абсолют-

619

ной скорости газаС1 = ~2 + L\c и С1 < С2 • Отметим, что у элемен­

тарной ступени, как это видно из рис. 9.32, всегда

.6.с = L\C/J, =-L\w = -L\wu'

(9.219)

В РК такой ступени вся работа, подводимая к газу отлопаточ­

ного аппарата, преобразуется в кинетическую энергию выходя-

щего потока. Поэтому ступени сa~~2 = О ~азывают .активны- ми•. Уровень абсолютныхскоростей в НАступени сa~~2 = О весь­

ма высок. это ограничи,вает работу, подведенную к газу в ступе­

ни, так как при высоких скоростях ка входе и связанных с ними

числах Маха растут и потери в НА. Особенно значительно поте-

ри увеличиваются, если MI:

~ М ' гдеМкр = 0.9 + 0,95 - кри­

2

кр

тическое значение ЧИСJlа Маха.

Вступени cQi~2 = 0,5 (рис. 9.33, б)~ДHaполовинатеплопере­

пада срабатывается в РК, а другая - в НА. В РК относительная скорость уменьшается по значению 102 < 101' В НА настолько же

уменьшается абсолютная скорость С3(I) < С2 • Треугольники ско-

ростей ступени сa~~2 = 0,5 симметричны, так что W1 = С2'

W 2 =Сl'Рl =(12 И 132 =(11'Вследствие этого скорости потока в РК и НА изменяются в одинаковой степени W 2 /W1 = С12 ,а значит,

будет t2 -it ='з(1) - ~ =0,5 (t; - tn . .

Такое распределение теплоперепадов между РК и НА приводит

К тому, ЧТО скорости при входе в них одинаковы W 1 = С2 ' а это,

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

каждый лопаточный аппарат, подвести в ступени к газу наиболь­

шую работу и, значит, получить наибольшие отношения давле­

ний при сохранении достаточно высокого КПД.

. СТупень ca~~2 = 1,0 (рис. 9.33,8) отличается тем, что весь

теплоперепад срабатывается TOJILKO в РК. Абсолютная скорость в НА изменяется только по направлению С1 = С2 + L\c , оставаясь

одинаковой по значению С1 = С2• Это значит, что перепад энталь-

пий в НА t3(l) - tz = О, а следовательно, i2 - i1 =i; -i; и в соот­

ветствии с формулой (9.213) будет a~~2 =1,0.

В дозвуковых стационарных осевых компрессорах применяют

ступени с a~~2 = 0,5 + 1,0. ПО КПД они примерно равноценны. но

наибольшие отношения давления могут быть получены при

a~~2 = 0,5. В дозвуковых И трансзвуковых транспортных (авиа-

ционных) осевых компрессорах используют преимущественно

высоконапорные ступени с a~~2 = 0,5 + 0,6 и лишь у последних

ступеней a~~2 может возрасти до -0,7. В сверхзвуковых транс­

портных осевых компрессорах могут применяться ступени

с a~~2 = 0,3 + 0,5 .

Недостатками ступеней сa~~2 = 0,5 являются крутая характе­

ристика и более узкий по расходу/диапазон работы, чем у ступе-

ней с a~~2 =1,0 .

Отметим, что ступени сQi~2 = 0,7 + 0,8 позволяют полу­

чить(хl == (хз == 900 (рис. 9.34). это Означает, что поток при входе

в РК и выходе из НА имеет осевое направление. Применение таких ступеней в качестве первой и последней в многоступенча­

том осевом компрессоре дает возможность упростить его конструк­

цию, отказавшись от установки входиого направляющего перед

первой и спрямляющего аппаратов после последней ступени[74].

В разделе, посвященном профилированию лопаток осевого ком­

прессора по высоте, будут рассматриваться элементарные ступе-

ни~ у которых иl = и2 = и , но c1z ~ C2z' В зависимость (~~.286)

для коэффициента реактивности таких ступеней кроме а1-2 вхо­

дит член, содержащий коэффициенты расхода Ч'lz и Ч'2z' Однако и в этом случае связь между окружными составляющими скоростей определяется кинематическим коэффнциентом реактивности. Это

своЙство a~~2 будет использовано при определении углов потока

в ступенях с разными коэффициентами расхода при IJходе и вы­

ходе из колеса.

Лопаточные решетки осевых компрессоров. В прикладной теории И инженерных расчетах установившихся режимов работы

лопаточных машин, течение

газа в которых всегда нестаци­

онарно, используют параметры

осредненного потока, не зави­

сящие от времени. Такое допу­

щение позволяет рассматри­

вать течение газа как устан0- вившееся и стационарное, что

намного упрощает задачу ис­ следования.

Линия то"а представляет

 

и

 

собой векторную линию, в каж­

 

 

 

 

дой точке которой вектор ско­

Рис. 9.34. ТреyroJIЬВИlCИ скоростей эле-

рости газа направлен по каса-

тельной к ней. Трае"тория

меитариой ступени с a~~2 =0.785

621

620

движения частиц газа - это геометрическое место точек последо­ вательllых положений частицы в пространстве в следующие друг за другом моменты времени. В установившемся стационарном те­ чении линии тока и траектории совпадают [33]. Линии тока,

пересекающие входное сечение круговой решетки на одном ра­

диусе, образуют поверхность то"а, форма которой при дальней­

шем движении потока после входного сечения в общем случае

может изменяться. Векторы скорости газа lJаправлены по каса-

тельной к поверхности тока.

.

Обобщение и анализ опытных данных позволили установить,

что в круговых лопаточных решетках, применяющихся в осевых

машинах, форма осредненцых поверхностей тока близка к ци­

линдрической, причем цилиндры практически соосны с решет­

кой. Это дает основание расчленить ступень осевой машины

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

тарных ступеней. Как уже отмечалось выше, развертка на плос­

кость цилинДрического сечения элементарной ступени, что при­

мерно соответствует развертке осредненных поверхностей' тока,

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

решеток профилеЙ.

Рассмотрим геометрические параметры плоской решетки и про­

филей, ее составл~ющих (рис. 9.35). Профили, образующие ре_

Шe'fку, расположены на одинаковом расстоянии, под углом к пря-

Рис. 9.35. ГеометричесICИе и режиllDlЫе параметры лопаточкой реmетICИ колеса

осевоro JCOашрессора

мой линии ММ, называемой фронтом решет"и. ГеОl\fетрическое

место центров окружностей, вписанных в контур профиля, пока­

занное на рисунке штрихпунктирной линией, образует среднюю

линию профиля, которая делит пополам расстояние между его

выпуклым (спинкой) и вогнутым (корытом) участками. Отрезок, стягивающий среднюю линию профиля, называется хордой nро­

филя.

Профиль характеризуется следующими геометрическими пара­

метрами: длиной хорды Ь; Щlксимальной толщиной Стах' которая равна наибольшему диаметру окружности, вписанной в контур профиля; максимальной вогнутостью f тах профиля, отсчитывае­

мой от хорды до средней линии; расстояниями от носика до сече­

ний профиля с максимальной толщиной x~ и вогнутостью а, от­

считываемых вдоль хорды; углами Хl и Х2 изгиба входной и вы­

ходной кромок профиля, образуемыми хордой и касательными к средней линии на этих участках; входным ~lл И выходным ~2л

лопаточными углами, между касательными к средней линии при

входе и выходе и фронтом решетки. Для решетки неподвижного НА эти углы обозначаются <Хи <ХЗл соответственно, углом е кри-

визны или изогнутости профиля е =Хl + Х2 = Р2л - Рlл; углом ус-

. тановки профиля \) (углом выноса) между хордой и фронтом ре­

шетки.

Безразмерные линейные параметры оп~еляются в долях дли­ ны хорды Ь и обозначаются теми же символами с чертой сверху.

Режимные параметры профиля в решетке следующие: углы потока ~1 при входе и ~2 при выходе из решетки, угол атаки

(натекания) i2 =Рlл - Рl ; угол отставания потока при выходе из

решетки о = Р2л - Р2

Решетка характеризуется такими геометрическими парамет­

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

параллельного фронту решетки; относительным шагом решетки

t/b Или густотой решетки b/t; углом отклонения потока в решет­

ке ~p = Р2 - Рl или ~a =аз - а2 между векторами скорости W 2

и W 1 В РК или СЗ и С2 В НА. В компрессорной решетке с диффузор­

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

в пределах L\~ =L\<x =15 + 250 •

Силы, действующие на профиль в решетке. Прежде чем ана­ лизировать силы, действующие на профиль в решетке, заметим,

что в идеальной несжимаемой жидкости, не имеющей вязкости

и связанных с нею потерь на трение, машины динамического

принципа действия работать не могут. Это непосредственно сле­ дует из установленного Даламбером парадокса, справедливого для

тел, которые ограничены замкнутой поверхностью: при безвихре­

вом стационарном обте"ании тела "оне'!ного размера идеаль­ ной несжu.м.аемоЙ жид"остью и отсутствии во"руг тела ис-

622

623

Рие.9.36. Бесциркуляциоввое (а) и ЦИРКУJJJlЦИОВВое (6)

обтекание профИЛJl лопатки

. .

точн-иков ши сmoков zл.aвн-ы.й вектор CIlA давАен-ия патока на

теАОравен- Н-УАЮ[33]. Обтекание профилей в этом случае будет бесциркуляционным (рис. 9.36, а), а положение задней точки В

схода потока - произвольным, зависящим 'rQлько от положения

профиля в потоке. Однако в реальной вязкой жидкости при без­ отрывном течении точка В располагается на задней кромке про­

филя, очерченной в выполненных конструкциях окружностью

малого радиуса. Кроме того, теоретическое исследование крыло­

вых профилей различного типа показало, что если задняя кром­

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

кривизны, то при таком характере обтекания, как показано на

рис. 9.36, а, скорость потока в этой точке должна быть беско­ нечно большой. Это приводит К физически невозможным беско­ нечно большим отрицательным давлениям у задней кромки. 0606- щая эти факты, Н. Е.ЖуковскиЙ и с. А.Чаплыгин сформулиро­

вали такой постулат: среди бескон-ечного ЧUCАа теоретически

возможны.х обтеканий nрофuл.я с угАовой точкой на задней кром­

ке в действитеАьности осущесmвАяется nл.авное обтекание с кон-ечной скоростью 8 этой точке [33]. Это означает, что точ­ ка схода потока В при безотрывном течении всегда будет совпа­ дать с угловой точкой на задней кромке профиля (рис. 9.36,6). Чтобы обеспечить такой характер обтекания, на течение идеаль­

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

Жуковского-Чаплыгина. (Заметим, что при течении вязкой жид­

кости циркуляция возникает самопроизвольно вследствие дейст­

вия сил трения, обусловленных вязкостью).

ЦиРКУАяцией скорости w по контуру L называют интеграл

г = Jw_ds = Jwcosyds,

(9.220)

L

L

 

в котором у - угол между вектором скорости w и касательной

к элементу контура ds; ш, - проекция'Скорости w на касатель-

ную к контуру. Циркуляция ПО замкнутому контуру Г = §w,dS.

L

Будем рассматривать безвихревое или потенциальное устано­

вившееся плоское течение жидкости в решетке, отличающееся

тем, что для него существует некоторая функция ср(и, z), назы­

ваемая потенциалом скорости, частные пр()изводные которой по координатам x~ у равны проекциям скорости потока на соответ-

ствующие оси: 'Шх =d<pjdX; Wz = d<p/dz. В гидродинамике дока­

зывается [33], что при потенциальном, безвихревом течении цир­

куляция ПО произвол~ному контуру, охватывающему профиль, равна циркуляции по поверхности самого профиля, если внутри

выбранного контура нет источников и стоков. Используя это по­

ложение, можно найти циркуляv;ию следующим образом.

Выделим в решетке охватывающий один из профилей контур abcd (рис. 9.37) таким образом, что произвольно проведенные,­

но эквидистантные участки аЬ и cd будут смещены на расстоя­ ние, равное шагу решетки t, а участки ьс и da, параллельные

фронту решетки, будут находиться на таком удалении от нее,

чтобы поток можно было считать невозмущенным. Это о.значает, что вдоль участков Ьс и da векторы скоростей W 1 и W 2 будут постоянны по значению и направлению. Определим циркуляцию по замкнутому контуру abcd, применив формулу (9.220) последо-

Рие. 9.37. СИлы, действующие на ПрофИJIЬ в осевой решетке

624

40

П/р л. с. Тимофеевекого

625

 

 

вательно ко всем четырем отрезкам контура. Примем во внима­ ние, что на выбранных нами эквидистантных участках аЬ и cd,

длинаи характер распределения скоростей вдоль которых одина­

ковы, циркуляции равны по величине и противоположны по зна-

ку ГаЬ = -r("d. (Здесь принято, что положительное направление обхода контура - против часовой стрелки.) Тогда

4

Г = raocd = I,r, = гос + ГМ = -twzu + tw1u

=

'=1

.

 

=t(w1u -W2u) = t(C2u ).

(9.221)

Полученная зависимость справедлива как для РК, так и для НА.

Теорему Н.Е.Жу"овс"ого о подъемной силе, установленную

как для единичного профиля, так и для решетки профилей [33, 74,89], мы приводим без доказательства, которое желающие мо­ гут найти в цитируемых работах. Ее можно сформулировать так: подъемная си-лаR, действующая на профиль в nлос"оЙреш.ет"е единичной высоты от обте"ающего его nото"а идеальной не­ сжимаемой жид"ости, равна произведению плотности жид"ос­ ти р, цир"уляции с"орости во"руг профиля Г и средней ве"тор­ ной с"орости nото"а W m; направление ве"тора подъемной силы

R определяется поворотом ве"тора сре,дней с"орости на 1[/2

в сторону, противоположную направлению цир"уляции. (На­ правление циркуляции обычно совпадает с направлением движе­

ния потока на выпуклой стороне профиля и противоположно ему на вогнутой; см. рис. 9.28.)

Математическая запись теоремы Н. Е. Жуковского такова:

Rид = prwm

(9.222)

Высота решетки l, равная единице, здесь явно не записана,

что необходимо учитывать при определении размерности подъем­

ной силы: ньютон на метр (Н/м).

Здесь циркуляцию Г находят по зависимости (9.221), а сред­

няя скорость W т равна

W

= W 1 + W2 И

W1 + W2

(9.223)

m

2

W m =--2-·

 

Угол наклона средней скорости к фронту решетки определяют

по соотношению (рис. 9.37)

ctg~m = ~tg~l + ctg~2 .

(9.224)

2

 

Плотность потока при малых дозвуковых скоростях газа МШ1 ~ 0,3 можно принить постоянной И равной плотности тор­ можения при входе в решетку. При более высоких скоростях

0,3 < Мш1 < 0,8 + 0,9 необходимо учитывать влияние сжимаемости

626

и принима'l'Ь среднее арифметическое или среднегармоническое значение плотности [74,89]:

р = Р1 + Р2 или Р = 2Р1Р2. (9.225); (9.226)

2

Р1 + Р2

Для неподвижных решеток вместо относительных скоростей W

записываются абсолютные скорости с и вместо углов ~ COOTBeТCT~ вующие им углы а в абсолютном движении.

Проекцию подъемной силы R на направление переносной скорос­

ти и для участка плоской решетки единичной высоты и длиной t

находят из уравнения импульсов в проекциях на направление и

Rиид = -рtioz(W2U - W1u ) = -prwz

(9.227)

Проекция R на осевое направление Оси z в общем случае равна алгебраической сумме сил давления, действующих по обе сторо­

ны решетки, и силы, возникающей из-за изменения импульса

в проекциях на направление

(9.228)

Для элементарной ступени и для плоской решетки единичной

высоты при р = const

будет W2z = w1z

= Wz

, и тогда зависимость

упрощается

 

 

..

(9.229)

Заметив, что при отсутствии потерь и течении несжимаемой

жидкости из уравнения (7.10) после интегрирования следует

р/р + w 2 /Z= const,

и записав, используя это выражение, разность давлений, найдем

(9.230)

Отношение

(9.231)

указывает на ортогональность векторов R и W т ' что непосредст­

венно видно из силовых треугольников, показанных на рис. 9.37,

если рассматривать их с учетом этого соотношения.

Подъемная сила профиля при течении идеальной, невязкой

жидкости определяется зависимостью

w 2

(9.232)

Rид =R/lИ)!, =С/lидрЬf,

627

где C/I И,II - коэффициент подъемной силы при обтекании идеаль­

ной жидкостью; Ь - хорда профиля.

При обтекании профиля вязкой жидкостью подъемная сила

уменьшится вследствие того. что часть энергии потока будет за­

трачиваться на преодоление сопротивления от сил трения. В ре­

зультате появится сила профильного сопротивления

w2

Rж = СжрЬf. (9.233)

гдеСж ---: коэффициен;профИЛЬНОroсопротивления.изопределения

которогос.ледуеточевиднОесоотношение ежид =О,таккак Rжид =О.

Рассмотрим обтекание одной и той же решетки потоками иде­

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

ростей перед и за решеткой. На этом режиме в обоих случаях

будут одинаковыми составляющие скорости wz ' W 1u И w2u ' а зна­

чит, как следует из уравнения (9.227), и окружные составляю-

щие подъемной силы Rжи,ll = Ru (рис. 9.37), так как уравнение

импульсов справедливо для любого течения и не зависит от по­

терь и характера течения между рассматриваемыми сечениями.

Видно, что при обтекании вязкой жидкостью из-за появления

силы профильного сопротивления Rж подъемная сила R/I стала

меньше. а вектор суммарной силы R сместился относительно силы

RИJI на угол '1', определяемый зависимостью

-

tg'l'= Rж/R/I = Сж/Су = Jl = l/К.

(9.234)

.В этом выражении

 

 

 

R

= С

W 2

(9.235)

pb---1!!..

/1

11

2

 

- подъемная сила ПрофИJlЯ; С - коэффициент подъемной силы

при обтекании профиля вязкоl жидкостью;

(9.236)

-качество профиля и Jl - обратное качество профиля.

Хара"теристи"и nлос"их решето" определяют опытным пу­

.тем на специальных газодинамических стендах. в которых ре­

шетку обтекает поток газа. Ковструкция стенда позволяет изме­

нять угол атаки и скорость потока при входе в широких преде­

лах [75. 89]. В процессе исп~таний по результатам измерен~

опредеЛяют скорости W p W 2' углы потока 131,132' давлеНИЯР1' Р1 и Р2' р; статические и заторможенного потока при входе и выхо­

де из решетки.

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

ют. а аэродинамические коэффициенты рассчитывают по полу­

ченным данным.

Типичная характеристика плоской решетки при малых чис­

лах MWj s ~3 + 0,.4 представлена на рис. 9.38. Углы отклонения

потока дЛЯ РК и НА обозначены одинаково и определяются так:

для РК е =АР = tl2 - tl1 И дЛЯ НА е =Аа =аа - а2' Исследован­

ная решетка составлена из профилей с близким к предельноМу для осевых компрессоров углом изогиутости профиля е = 400, вход-

ным и выходным лопаточными углами Р1л = 42.5 о, 132л =82.5·.

с максимальным прогибом средней линии, расположенным посе­

редине хорды а/Ь =o~ и густотой решетки b/t = 1,06 . На основа­

нии накопленного опыта установлено, что при одних и тех же

значениях чисел Маха М и Рейнольдса R, углы отклонения пото-

ка &. коэффициенты C/I и Сж однозначно определяются углом ата­

ки i при входе в решетку. (Здесь и в дальнейшем в скобках указа­

ны величины, относящиеся к НА, так как приводимые зависи­

мости имеют один и тот же вид как для рабочей, так и для направляющей решеток.) Из характеристик видно, что по мере увеличения угла атаки коэффициенты e/l и С/lИ,ll возрастают, од­ нако при больших значениях i их увеличение прекращается. Уroл

поворота потока в решетке & с ростом i также увеличивается,

достигая своего максимального значения ефах при положитель­

ных углах атаки, но по мере их дальнейшего увеличения & начи­

нает уменьшаться. Это объясняется тем, что при превышении

некоторого критического угла атаки i KP ' соответствующего ешах ,

пограничны'й слой отрывается от стенки профиля. Угол поворота

потока в решетке начинает уменьшаться, наблюдается резкое уве­

личение потерь, а значит, и коэф-

фициента еж решетки. Как пока-

зывает опыт, при t = tKP коэффи-

циеНТПрофильногосопротивления

Сж возрастает пример но вдвое по

сравнению с минимальным значе-

нием, соответствующим углу ата-

ки iопт• Отрыв потока возникает по­ степенно. При этом угол поворота

потока еще продолжает увеличи­ ваться, хотя и медленнее, чем рань­

ше, а потери в решетке и коэффи­

циеН'l' начинают возрас.тать. Режи­

му начала отрыва потока от сте­

нок профиля соответствует угол

атаки iио

НОJttин.альн.ым.. илирасчетн.ым.. режи.моJttработы решет"и приня­

то считать режим, при котором

угол поворота потока в ней состав-

&;о

 

 

 

 

 

 

 

 

 

~

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;0

 

 

 

 

 

 

 

~

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С!/..S;

 

~

~

 

 

.!:-

 

 

 

 

 

 

 

 

 

~

 

0,5

 

 

J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s!.~-

r -

V

1--

 

 

 

 

 

 

 

 

 

zo

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

о

 

DV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0;0

 

О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'-.... cJ(

 

 

 

 

 

.....

Г!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

!

 

 

 

li

 

 

 

 

 

-20 -15 -10J -

 

 

 

 

 

 

 

 

 

1.. 1.0

51

10

 

Рис. 9.38. Газодинамические ХА­

рактеристики плоской решетки,

полученные при ее продувке на спеЦИAJIЬНОМ стенде

628

629

 

ляет 80% от максимального. Параметры номинального режима

обозначаются индексом .*.: е*, L\f3*, L\a*, t* .. е;, е; и т. п.

Анализ большого числа опытных данных, полученных при

небольших числах М на различных решетках, позволил А. Хоу­

эллу установить, что возможно обобщение опытных характерис­

тик, если построить их в относительных, безразмерных коорди-

*

 

t*)

 

"

 

 

 

~

(рис. 9.39); Обычно номиналь-

натах в виде е/е, еж = f (i -

 

ный угол атаки находится в пределах i*

= ±5 0

Видно, что при

увеличении i от отрицательных значений до

i*

угол поворота

потока в решетке & возрастает с увеличением угла атаки практи­

чески линейно. Это свидетельствует о том, что угол отставания

потока от лопаток 8 (см. рис. 9.35) в этой области характеристи­

ки примерно постоянен. При увеличении i свыше i* из-за воз­

никновения отрывного течения угол отставания довол~но резко уве­

личивается, что приводит к уменьшению &. Штриховой линией на

рис. 9.39 показан режим с постоянным углом отставания 8.

Некоторые параметры решеток на режимах, близких к номи­

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

зависят от относительной густоты решетки b/t и одного из ха­ рактерных углов: В, f3m{a m) или f32{аз) . Так, номинальный угол поворота потока е* в решетке может быть представлен в виде

 

 

 

 

 

 

 

 

S~O

"

 

 

 

 

 

/ ....... -.....

50

 

 

 

l/

к§.

 

АД

~o

 

 

 

 

 

 

е;"

 

 

 

 

 

 

 

 

 

СК

48

 

 

 

1/

 

 

 

 

 

 

0,10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;f

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,08

 

0.6

 

 

 

 

 

 

 

I

 

0.06

 

2\

 

 

 

 

 

 

I

 

0,0+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i'{--

 

 

J

 

 

о,ог

 

 

 

 

 

 

 

 

 

 

 

 

 

..

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

0,0 -

4~ -

О)

о

 

 

Рис. 9.40. Зависимость иоми·

 

Рис. 9.39. Обобщеиная харак·

 

теристика ПJIоских решеток

В8JIЬИого угла поворота потока

 

 

от угла выхода потока и густо­

ты решетки при &* = 0,8&...,.; t =±5'

е* = f(b/t, f32(аз)) (рис. 9.(0). Видно, что с увеличением густоты

решетки растет и е* , так как расстояние между профилями умень­

шается и благодаря их взаимному влиянию отрыв потока насТу­ пает при более высоких углах атаки.

Обратное качество профиля может быть обобщено в виде зависи-

мости J.1 = f(b/t, ctgf3m(am)) (рис. 9.(1) [75, 89]. Значения котан­

генсов средних углов определяют по формулам (9.258) и (9.259), полученным ниже. Зависимость обратного качества от КОэффици­ ентов расхода и реактивности (рис. 9.(2) для решетки с густотой

b/t =1,0, найденная пересчетом соответствующей характеристи­

ки на рис. 9.41 [89], показывает, что с увеличением нагрузки на

решетку J.1 возрастает. На рис. 9.41 и 9.42 дЛЯ РК а(К) = a~~2

и для НА а(К) = 1-a~~2'

Влияние числа Мшt (для решетки НА Мс2 ) на характеристи­

ку решетки неОДН03начно и в значительной .степени определяется

режимом ее работы.

Минимальная площадь межлопаточного канала единичной высоты f 1е определяется наименьшим диаметром окружности, вписанной в канал (рис. 9.43, а). Площад~.потока'1 при входе

в канал определяется как 11 =tsinf31 ' где t - шаг решетки.

Межлопаточный канал можно условно разделить на два участ­ ка. Первый - 'косой срез на входе и второй - канал, располо­

женный после сечения ''.Режим работы косого среза переме­

нен и зависит от отношения пл~щадей 1/11 При 'высоких 131'

Рис.9.41. Зависимость обратного качества профиля от ryc-

тоты решетки и средиего угла потока (для РК

n(К) = n(К) ;

дЛЯ НА n(К) = 1- n~~2 )

1-2

631

630

0,7 0,9

Рис.9.42. Зависимость обра'I'JЮГО качества про­

фили от коэффициеита расхода и кинемати­

ческого КОr..t!~еита реl1К'l'ИВROCти (b/t ",1,0)

(да PIC а(к = al~2; д.!UI НА .а(к) = 1- a~~2)

000ТIIe'IC'l'IJбольшим

расходам газа, '1е < ' 1 И

поток В косом срезе ус-

коряется: наблюдается

конфузорное течение

газа. При малых ~1'

когда расходы меньше

ноМинального, '> ' 1 И

поток В косом срезе за­

медляется. Этот режим

с точки зрения уравне-

"ия скоростей и газоди­

намических потерь наиболее благоприятен, так как соответст­

вует общему характеру диффузорного течения в элеМентах компрес­

сора, однако только до момента возникновения отрывного течения.

Ма1tси,м,а.яьны.м ЧUCJЮ,м, Мш1тах при входе в решеm1tу называ~т

такое, при котором в минимальном сечении канала устанавли­

вается скорость, равная скорости звука: Мш1е = 1. Массовый

расход газа на этом режиме достигает своего максимума. Теоре­ тически при течении без потерь изменение параметров газа опре­

.деляется уравнением изоэнтропы, и соответствующие максималь-

ные числа Маха обозначены MwlmaxT (рис. 9.43, бj. В этом

случае при f1e /f1 =1 и M w1 = 1 будет и MwlmaxT = 1. При

действительном характере течения с потерями и том же отно-

шении площадей Мш1е =1 устанавливается при Мш1тах == 0,73.

При отношениях f1e /f1 > 1 Mw1max возрастает до - 0,98 и за­

тем из-за увеличения потерь начинает уменьшаться.

Криmичес1tи,м, число,м, Мш1кр при входе в решеm1tу называют

такое, при котором скорость потока достигает скорости звука

в какой-либо области течения, прилегающей к профилю. Обыч­ но это происходит у стенки профиля вблизи от входного сече­ ния канала 1е .Увеличение МШ1 свыше Мш1кр приводит к воз­

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

сверх- к дозвуковому течению образуется система скачков уп­

лотнения, вызывающая отрыв потока и повышенные газодина­

мические (волновые) потери. Чтобы сохранить высокий КПД

c;zпени, следует работать на таких режимах, при которых wl < Mw1кp

Наибольшего значения 0,65-0,70 число Mw1KP достигает при f1e /f1 = 1,15 . На режимах M w1 < Mw1кp минимальные значения СХ и соответсТВУЮЩJlе им Су практически неизменны, но при

отклонении режима обтекания от оптимального коэффициент про-

Рис. 9.43. Влияние режима иатекаиии JЮТока на решетку иа пре. дельвое число Маха при входе: а - измеиение IШощади потока

Dpи JlXоде в решетку; (j - зависимости максимальвых и критичес­

кого чисел Маха от отиоmеиии ПJIощадей f1./f1

фильного сопротивления Сх увеличивается (рис~ 9.(4). Коэффи­

lI.иент подъемной силы Су остается практич~ки таким же, но его

наибольшее значение с ростом МШ1 уменьшается.

Для определения минимального коэффициента профильного сопротивления рекомендуется эмпирическая формула, получен­

ная в результате обобщения опытных данных [89]

-

ь

(9.237)

CxDliв =0.012+0,042fmax +0,0023-,

 

t

 

где fmп =fmaх/Ь.

Этазависимость справедлива при Mw1 < Mw1кp для pemeroкс про­

филями относительной толщины 0,05 < ё < 0,010 при углах ус­

тановки 45· < '\)< 100·.

Угол изогнутости профиля, как следует из рис. 9.35, опреде­ ляют с помощью такой зависимости:

е=Е+~-t. (9.238)

Для решеток осевых компрессоров номинальный угол поворо­

та потока находится в пределах &* = 10 + 30·, а угол атаки реко­

мендуется для номинального режима выбирать в диапацоне t* =о + 5 •. Чтобы при профилировании решетки найти е необхо­

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

632

633

 

 

нального режима 8 можно опре­

 

 

делить из такой эмпирической за­

 

С)(

висимости [12,74,75,89]:

 

'~--~-~~-4-т-+-г~o,20

8* = тe.fiiЬ '

(9.239)

 

 

~-1~-+--~~~~o,~ где

 

.~--f--+--~-t-+---iО,10

т = 0,92 (а/Ь),2 + 0,002~;, (9.240)

 

 

причем у осевых компрессоров

 

 

а/Ь =0,40 + 0,50 .

 

 

 

Совместное решение уравнений

-5

 

(9.238) и (9.239) дает такое выра­

жение:

Рис. 9.44. Влияние числа Маха при

8* = тfi7Ь(&* -i*)

 

входе в решетку ка ее характерис-

(9.241)

тики

 

 

1- тfi7Ь

 

Для решеток НА в формулу (9.240) подставляют угол а;.

Связь между rидраВJIическим КПД элементарной ступени и основными параметрами плоской решетки. Из определения гидравлического КПД, в~раженного зависимостью (9.207), для

ступени осевого компрессора можно получить такие выражения:

11" = lэф

=

-41-3 =1_41-3 =

 

.

 

 

za

I

 

41-2

'r2-3

= 1- L\11h1-2 -

 

(9.242)

=1- - ,- - - ,-

 

L\11h2-3'

ээ

в этих выражениях 41-2 ='rl-2 + 'r2-3 - работа, затраченная

на преодоление сопротивлений в ступени; 'rl-2' 'r2-3 - работа,

затраченная на преодоление сопротивлений в РК и НА соответст­

венно; 'э' lэф - теоретическая (эйлерова) и эффективная работа ступени; L\11hl-2' L\11h2-3 - снищение гидравлического КПД из­

за потерь в РК и НА соответственно.

Мощность, потерянная в РК элементарной ступени, отнесен­

ная к единице высоты лопаток, выраженная в ваттах на метр

(Вт/м), затрачивается на преодоление силы профильного сопро­

тивления R XP K (9.233)

w3

 

.

 

N r1 - 2 =Z1 Rxp.KW m =Z1 CXP.KPbp K f.

(9.243)

Массовый расход газа, также отнесенный к единице высоты

лопаток, кг/(с,м),

G =PWzZ1t1 =PCZZ2t2'

(9.244)

в этих формулах Z1'

t1 ,

z2'

t2 -

числа лопаток и их шаг по

фронту решеток РК и НА соответственно,

Ьрк - длина хорды

профиля РК.

 

 

 

 

 

Удельная работа, потерянная в РК, с учетом того, что W z =Cz '

будет равна

 

 

 

 

 

 

 

 

 

 

(9.245)

Заменив размерные величины безразмерными, получим

 

- с

Ь

-3

2

l,. 1-2

Р·К

q> m 1

-

Хр.К -t-T-U

(9.246)

 

 

 

1

q>z

 

Подобным же образом найдем мощность и удельную работу,

потерянную в НА,

 

 

 

 

 

 

 

 

 

 

Nr2-З

= Z2Rхи.аСm = Z2Схи.аРЬ

С3

;

(9.247)

 

;

l,. .

= Nr2-З

= С

Ьи.а

C~.!. = с

Ьн.а

<p~ .!..и2•

(9.248)

2-3

G

ХИ.а t

2

C

ХИ.а

t

2

2 q>z

 

 

 

 

2

 

z

 

 

 

 

В этих выражениях Z2' t2

-

число лооаток НА и их шаг

соответственно; Ьна -

хорда профиля НА.

 

 

 

Теоретическая (эйлерова) работа ступени l'yIожет быть выраже­ на как через параметры решетки РК, так и через параметры решетки НА. Для этого необходимо отметить, что в соответствии

с уравнением равновесия на установившемся режиме окружные

составляющие суммарных сил, действующих от потока на профи­ ли решеток РК и НА, равны по значению и противоположны по

направлению Ru р.К =- Ru Н.а

Теоретическая (эйлерова) мощность ступени с учетом этого об-

стоятельства

-

 

 

Nэ = ZI Rup.Ku = Z2Ruи.аU

(9.249)

Окружные усилия на профилях РК и НА (см. рис. 9.37):

RUP.K = Rp.kCOs[90-(~m + "'р.К)] = Rp.Ksin(~m + "'р.К); (9.250)

RuH•a = Rн.acos[90- (ат + 'l'и.~)] = Rн.asin(am + "'и.а). (9.251)

Заметив, что в обоих случаях справедлива зависимость

RII /R = COS""

найдем из этих выражений:

 

sin (~т + '" р.К)

 

R UPK = RIIP.K-~---=--"':"

 

 

COS"'P.K

 

Ruи•а = RlIи•а

sin (ат + '" и.а)

.

COS \11

= RlIи.аSlП ат(1 + ctg amtg "'и.а)·(9.253)

ТИ.а

635

634