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

книги / Центробежные компрессоры

..pdf
Скачиваний:
14
Добавлен:
20.11.2023
Размер:
14.98 Mб
Скачать

При выборе размеров входа следует учитывать нежелатель­ ность слишком малых втулочных отношений. Если rmJr0 < < 0,4 -4-0,45, ухудшаются условия обтекания корневых сечений лопаток. Минимальные числа Мш на входе достигаются в том случае, если по формуле (2.11) относительная скорость миними­ зируется на fi. Скорость — средняя по высоте лопатки — при этом больше, чем в случае минимизации w[ на некотором сред­ нем между г0 и гвт радиусе rv Рекомендуется с помощью урав­

Уодmax

нения

(2.11)

определить г0

двумя

способами: минимизи­

1,0

руя входную

относительную

 

 

скорость на

го и на сред-

£ РКmin

0,015

0,0Z

0,025

0,05 0,055bz/Dz

 

 

 

Рис. 2.54. Максимальные к. п. д.

Рис. 2.55. Зависимость

 

и

оптимальные

коэффициенты

минимального

коэф­

 

напора для колес (-------- ) и сту­

фициента потерь колес

 

пеней РК +

БЛД (-----------) в за­

«радиальная звезда» от

 

 

висимости от bjD^:

 

отношения w

 

1 — г0 =

0,468;

2 — Го = 0,526;

 

 

 

 

 

3 — г,

- 0,61

 

 

 

 

неквадратичном

радиусе

гг =

[0,5 (r0 + гвт)]0’5.

Выбор

г0опт,

лежащего

между

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

произвести

с

учетом

Ми и влияния

на уровень

скоростей гё)1,

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

Рл по

высоте

лопатки.

 

 

 

(величину

b2/D2)

следует назначать

исходя

Размеры выхода

из тех же соображений, что и для РК ПЦК с учетом нежелатель­

ности

w < 0,7 4-0,75. При выборе

меридиональных

обводов

нужно

учесть как конструктивные и

технологические

соображе­

ния, так и их сильное влияние на характер W. По возможности следует обеспечить быстрое замедление потока в начале и уско­ рение в конце межлопаточных каналов, чтобы получить ха­

рактер W, подобный показанному на рис. 2.51 и рис. 2.52 (сплош­ ные линии).

Средняя нагрузка на лопатки в связи с их значительной про­ тяженностью обычно несколько ниже, чем у РК ПЦК, и обычно составляет величину порядка 0,25. У РК с развитой осевой частью на средней ОПТ Ддоср ^ 0,2 в осевой части и примерно 0,3—0,35 — в радиальной. Зависимости рл = / (1т ) на осесимметричных

Щ

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

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

зон

Фр ^

0,06 —0,09 для РК

«радиальная

звезда». Очевидно,

они

могут быть

использованы

и при Фр >

0,09,

однако

следует

иметь в

виду,

что особенно

высокорасходные

ступени

обычно

Рис. 2.56. Измеренные распределения скоростей в осерадиальном пространственном РК ЛПИ с Ф0ит ~ 0,11^0,12 [43]:

— # — — лопатки первою яруса, —А — — лопатки второго яр>са

применяют для крупных ПЦК, экономичность которых весьма важна. Очевидно, что более высокий к. п. д. могут обеспечить осерадиальные РК закрытого типа с лопатками, спрофилирован­ ными по всей длине и коэффициентами напора ниже, чем у РК «радиальная звезда». Опыта использования этих колес отечествен­ ное компрессоростроение пока почти не имеет. Ступени с такими РК при умеренной расходности (Фр ^ 0,080) разработаны СКБ-К. Высокорасходные пространственные РК изучались в ЛПИ [43]. Показано, что использование методик, основанных на идеях пп. 2.3 и 2.8, позволяет получить желаемые распределения ско­ ростей на ОПТ и в МП.

для

На рис. 2.56 показаны измеренные распределения скоростей

колеса 703

[43] на периферийной и втулочной ОПТ. Осред­

нение по

0 измеренных w показало, что заданный закон ст =

— / (б) =

const

(«равноскоростное»

колесо)

примерно

выдержан

за

счет выбора

соответствующего

«навала»

лопаток

(см. п. 2.3).

На основе изложенных выше рекомендаций авторами совместно с А. М. Симоновым и М. Р. Полесом спроектирована и испытана ступень с высокорасходным РК 705.

Для изучения влияния меридионального профилирования и измерений неподвижно установленными в корпусе датчиками давления колесо на первом этапе было сделано полуоткрытым. Несмотря на это, эффективность колеса и ступени (рис. 2.57) высока, величина Фр ^ 0,16 превышает в среднем вдвое расходность колес первых ступеней выпускаемых промышленностью компрессоров. Конструкторские проработки показали возмож­ ность примерно двукратного уменьшения металло- и трудоем­ кости компрессоров за счет при­ менения этого и подобных ра­ бочих колес.

2.8. Приближенный метод

 

 

 

 

 

 

 

 

 

решения обратной задачи

0,09

0,11

0,13

0,15

0,17

0,19

0,21

Ф

 

 

 

 

 

 

Под

решением

обратной за­

Рис. 2.57. Характеристики колеса 705

дачи

понимается

определение

и ступени

РК +

БЛД

при

 

b2fD2 =

формы

лопаточной

решетки,

=

0,08,

г3 =

1,9,

Ми =

0,6

обеспечивающей то

или

иное

в области

течения. Как

пока­

распределение параметров

потока

зано

выше, распределение

скоростей невязкого потока по поверх­

ности лопаток на ОПТ дает возможность достаточно глубоко проанализировать свойства решетки и может быть положено в основу определения ее оптимальной формы. Наиболее важным является характер распределения относительной скорости по задней поверхности лопаток w3 = f (I).

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

определяется W на ОПТ по одному из точных методов. В случае

неудовлетворительного характера W в форму решетки вносятся коррективы и расчет обтекания повторяется. Описанные ниже приемы проектирования сводят необходимые корректировки к ми­ нимуму.

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

ников скоростей, формы профилей лопаток, средней нагрузки Дй)^р и желательному w3 = f (I) посвящены пп. 2.2, 2.5, 2.6 и 2.7. При изложении метода будем исходить из того, что отно­

сительные скорости

w-L и й)2» толщина слоя

6 = f

(?) и собственно

форма ОПТ известны, выбраны средняя

нагрузка и желательное

w3 = / (0*

Выбранной

меридиональной

форме,

величине

w3 =

= f (I)

и числу лопаток

z

будет соответствовать

некоторое зна­

чение

коэффициента

теоретического

напора,

которое

может

не совпасть

с заданным.

Поэтому

метод

решения обратной

задачи должен включать этап оценки

и

уточнения w3 =

f (/), 2

или б = f (/).

 

 

 

 

 

 

 

 

Определению формы лопаток на ОПТ предшествует определе­ ние линии тока осредненного по шагу потока Р = /(/), которое в соответствии с теоремой моментов и уравнением неразрывности обеспечивает заданное или возможно близкое к нему wà = /(/). Далее тем или иным способом определяется рл = / (0> т. е- форма средней линии лопатки на ОПТ, обеспечивающая необходимую траекторию потока. Форма профиля лопаток полагается выбран­ ной с учетом конструктивно-технологических соображений.

Основные уравнения

для расчета:

 

 

 

 

Ф

 

Ф .

/0

W

 

2г%гБ sin Р ’

ф

2ет гЪ

(2 -7 4 )

й>3 =.

-[-

0,5 Ай>;

wn =

w — 0,5 tiw.

(2.75)

Здесь w — среднерасходная скорость, поэтому уравнения (2.75) строго справедливы только в случае линейного изменения w = = /(?), что не вполне соответствует действительности (см. п. 2.3).

Поскольку ^форма лопатки определяется в зависимости от распределения скоростей невязкого потока, расчет степени сжа­ тия ведется по уравнению адиабатного процесса

е = [1 + ( k - 1)МЦг2- ш 2)Гт_1).

(2.76)

Коэффициент стеснения определяется по действительной тол щине лопаток в данном сечении

т = 1

zb

(2.77)

2 nr sin рл *

 

 

Систему уравнений приходится решать итерациями, в первых шагах которых при расчете .коэффициента стеснения т можно принять рл = р. Функция 6 = f (I) определяется выбранной формой профиля лопатки. У наиболее распространенных лопаток постоянной толщины с заострениями на входе и выходе для удоб­ ства решения можно принять закон изменения толщины в пре­ делах входного и выходного участков лопатки 1Ъ, например в виде степенной функции

ô = ô0CH(//^)tt,

(2.78)

где при значениях п = 2 для входного и п = 1,5 для выходного участков получаются формы, близкие к опробованным в практике.

Для построения

конструктивного чертежа лопатки уравне­

ние (2.78) в полной

степени не подходит, так как в точке I =

получается излом, который следует сгладить [это в малой степени влияет на т = / (I) и может в процессе профилирования в расчет не приниматься].

Изменение толщины слоев b в общем случае определяется расчетом осесимметричного меридионального течения. Для РК ПЦК за расчетную ОПТ принимается радиальная плоскость, за высоту слоя — высота лопатки b в направлении г. Тогда при рекомендуемых прямолинейных образующих основного и напра­ вляющего дисков

 

6 =

 

 

 

 

(2 79)

Нагрузка Aw может

быть

найдена

из уравнения

моментов

 

d (с„г) =

г Apbrdr/m,

 

где т = 2nrxpwr

 

(0

для

невязкого

потока при р = / (t) =

Из условия р* = f

= const

 

 

 

 

 

 

 

 

Рп -

Рг = Р А® (а»п +

о»э)/2.

 

При линейной

зависимости

до = f

(t)

справедливо

равенство

до = (доп + до3)/2.

Тогда

в

безразмерном

виде

 

Ддо =

 

z

Так как до3 = до +

0,5Ддо

= I(си — г)2 + ф2]0-5,

 

- de

 

 

(2 80)

Г~Т-*- + с 1 -SB-

drг

J

W '

 

[см. уравнение (2.75)]

и до =

ф [ГЧ Г + С“~ ]

т

(2.81)

®8 = [(Си — Г)* + ф2]0’5

z

Wu — г? + ф2]0,5

 

Численное решение уравнения (2.81) при известных соотноше­

ниях между входящими в него параметрами и аргументом f

или I

учетом dl = drlsin Р)

дает зависимость

си =

f (г) или

си =

=

/ (0»

что позволяет определить р = / (/)

или

р = / (г) из

оче­

видного

равенства

 

 

 

 

 

 

P =

arctg [ср/(г — си)].

 

(2.82)

Рассмотрим вопрос об ограничениях, налагаемых на вид функции w3= f (I) при проектировании лопаток. Очевидно, что, будучи связана со среднерасходной скоростью w, при поло­ жительной нагрузке Aw скорость на задней стороне может вы­ бираться только в некоторых пределах. Нас интересует прежде всего минимально достижимая скорость в начальной части ло-

Наток, величина которой определяет общий уровень скоростей в решетке и максимальное число М — оба параметра определяют эффективность.

Относительная скорость

w = (ш2 -)- wl)0,5 =

(wj -)- и2 — 2cuu +

+ cl)0-5 или

в безразмерном

виде

 

 

 

 

w = [ф2 +

г 2 — 2си (г +

с„)]°.б =

[ w \ ep — 2си (г + сй)]0-5, (2.83)

где

 

 

 

 

 

 

 

 

ЙУнеР = (Ф2 +

г2) =

Ф2

2 +

Г2 .

 

 

 

 

f

JL f

,

 

 

 

 

 

b2

р2

 

Последнее соотношение дает величину относительной скорости при «нерабочих», не отклоняющих поток лопатках. Так как при

этом Aw = 0, то w3 = wn =

ш. Относительная скорость с увели­

чением г быстро возрастает, и при

г = 1 численное

значение

w

больше единицы.

 

удержать w3 в

минимально

Таким образом, для того чтобы

возможных пределах при

Aw > 0,

надо вводить нагрузку,

но

ограниченно, чтобы снижение wy сопутствующее росту на­ грузки, т. е. росту си [см. уравнение (2.83)], не привело согласно

(2.81) к

росту w3.

Aw,

Если

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

то согласно (2.81) можно получить любую зависимость w3 =

/(/);

задача решается изменением знака при dcjdf. При некоторых значениях w3 на большей или меньшей части лопатки в области входа при гг ^ г ^ гвх обтекание будет носить «турбинный» характер, как при обтекании с /1пс < 0. Выше было показано, что для РК ПЦК режим небольшого отрицательного угла атаки, улучшающий обтекание задней и несущественно ухудшающий обтекание передней поверхности лопаток, обеспечивает близкую к минимально возможной величину £РК. При этом w3m3X ^ w[. Исходя из этого, при профилировании РК ПЦК в качестве пер­ вого шага можно выбрать условие w3 = w{ = const для основной части лопатки. В пределах от некоторого радиуса гвых до г = 1 следует обеспечить плавную разгрузку лопатки, например, в со­

ответствии

со степенной

зависимостью

 

 

 

 

 

W3— Ш3<ВЬ1Х

К (Г

Гвых)

 

 

При г =

1, т. е. на выходной кромке лопаток,

должна

быть

обеспечена

полная разгрузка,

т. е.

w3 = w2.

Тогда

К =

= (^з.вых — ®а)/(1 — гВых)п, и поэтому

г — 'вых у

 

 

 

вых

(0У3. вых

 

(2 84)

 

1

^вых /

 

 

 

 

 

 

 

У высокоэффективных РК ПЦК численное значение радиуса гвых> на котором начинается разгрузка, соответствует примерно

0,8

от радиальной

протяженности

лопатки,

т. е. гвых =

4-

+

0,8 (1 — Гх), а

показатель степени п = 2.

значения си = f (г)

 

Определенные

по

выбранному w3

— f (г)

[см. уравнение (2.81)] позволяют определить распределение ско­

ростей по передней стороне и

коэффициент

напора

фт = си2.

Следует проверить соответствие

полученного

напора

заданному

и приемлемость

wn =

f (f). Как уже указывалось выше, нежела­

тельны значения

wn <

0,25, а также излишне большой пик wa max

в случае отрицательной нагрузки

на входе wnmax > 1,1 w[. Если

этот пик больше, при повторном расчете следует принять w3l >w[. Увеличить wn можно за счет увеличения отношения w или

принятия

вместо

w3 = const

распределения

w3 <

1

(что

неже­

лательно

и

может

быть

принято, если ^

—-ч-----

 

7

увеличение w приводит к неприемлемому

'

 

уменьшению

б2

и Ù).

При

си2 < Фт. P

J

 

1

1

1 1

/ 1

следует

либо

 

увеличить

число

лопаток

°>ô

 

 

 

 

Т

и уменьшить

d>,

либо

принять

w3i>w[.

Q6__________ /

'

Два последних пути приемлемы, если по-

J

 

 

 

 

лученное

при

первом расчете значение wu

ДО

 

 

 

 

не слишком мало. В противном случае

 

 

 

 

 

 

может

потребоваться

принять

w3\ >

w[

ДО

 

 

 

 

и w3 <

1

одновременно.

 

 

 

 

 

 

 

 

 

 

 

Последний

 

этап

профилирования

 

0

0,2

0,4

0,6

0Л г-г,

определение

углов

лопаток

Рл = f

(г),

 

 

 

 

 

rz~Tf

которые

 

соответствуют

рассчитанному

Рис

2 58.

Эксперимен­

P =

f (г). Расчетные данные Стейница [96]

тальная

зависимость для

по

связи

формы

лопаток

с траекторией

углов отставания

потока

относительного

потока,

полученные

для

в межлопаточных каналах

средненапорных РК ПЦК

решетки «радиальная звезда», не вполне

с /Сдо ~

0,2

и Jinc. р ^ 0

соответствуют

 

 

условиям

 

течения

в

 

 

 

 

 

 

РК

ПЦК. Для

колес с

W,

близким

к рекомендуемому

в ме­

тоде ЛПИ, можно использовать экспериментальную зависимость (см. рис. 2.58), полученную В. П. Митрофановым по данным изме­ рений в относительном движении В. В. Козлова и В. И. Зараева. Величина 02 определяется соотношениями, приведенными в п. 2.4.

Для решеток с произвольным W может быть использован подход, предложенный А. М. Симоновым и Б. Н. Савиным при­ менительно к решению прямой задачи. Для нахождения связи относительной траектории и формы решетки используется тот факт, что при увеличении нагрузки по длине канала I струйки тока утоньшаются у задней и утолщаются у передней стороны. Поток как бы отжимается в области увеличивающихся скоро­ стей и поворачивает к задней поверхности лопаток. Это харак­ терно для входной части канала, а на выходе, в областях умень­ шения нагрузки лопаток, формируется отставание потока. Тор­ можение вдоль задней и разгон потока у передней стороны с соот­ ветствующим изменением толщин струек тока отклоняют поток от направления лопаток в сторону, противоположную вращению.

Для вывода расчетных соотношений в рассмотрение прини­

мается вектор относительной

скорости в

виде составляющих

w =

W|3 +

Aw„,

(2.85)

где wp — проекция относительной

скорости

на направление |3Л.

Составляющая Awu определяет несоответствие направления по­ тока и лопаток.

Для определения составляющей Awu рассмотрим уравнение

неразрывности

при

течении на радиальной ОПТ для

элемента

с размерами

6,

dr, dQ

 

 

 

 

 

 

 

 

 

 

 

sin $аг9Ь) +

pb

=

О,

(2.86)

где 0 — относительная угловая координата,

0 =

г 0/(2лт).

Делением

(2.86)

на

массовый расход

газа

через

РК т =

= 2яг6ршр ср sin рлт

получено

 

 

 

 

 

 

д

/ Awu \

.

2ягт2 sin рл

д /

1\ _ ^

(2.87)

 

дв VЩ ср

/

 

г

 

дг \

ср т } ~

 

 

 

 

Величина

Аш„ = /(0)

может

быть

найдена

интегрированием

(2.87) при

известной

функции

Шр =

/ (0). Принимая_ допущение

о линейном характере

изменения скорости

= / (0),

получим

 

 

 

 

 

=

1

Aw

(0 -0 ,5 ).

 

(2.88)

 

 

 

 

СР

Й&СР

 

С использованием (2.88), в пренебрежении производной от

коэффициента стеснения т по

радиусу,

интегрирование

(2.87)

по угловой координате

приводит к зависимости

 

 

 

 

nfTsinpnd (Дй/юр СР)

(02 -

0).

(2.89)

®3 ср

 

г

dr

 

 

 

 

 

Найдем осредненное по расходу значение Awucp

 

 

1

 

 

 

 

 

 

 

|

&wuw$dB

 

 

 

 

 

Дйиср _

о________ ягт sin рл d

/

Аш

\

/2 90)

®р ср

-1

.

dr

Vгг>р СР У‘

*■ ' '

 

J йрй0йрСр

 

 

 

 

 

 

о

 

 

 

 

 

 

Полученное уравнение, дополненное уравнением связи компо­ нент скорости с учетом а>р ср sin рл = <р, ctg рл = ctg р — Дшы/ф, позволяет найти угол лопатки Рл при известном законе измене­

ния

нагрузки и известной функции Р = / (г).

В ЛПИ описанный метод профилирования использован в блоке

программ, которые по заданным Фр, фт. р,

М„ и конструктив­

ным

параметрам (см. п. 2.2) осуществляют

меридиональное и

радиальное профилирование РК ПЦК, расчет обтекания методом

интегральных уравнений в заданном диапазоне Ф и определение характеристик ступени РК + БЛД по математической модели потерь и напора (см. гл. 6). В состав ступени могут быть вклю­ чены также и другие неподвижные элементы с использованием эмпирических зависимостей для потерь из гл. 5. Перебор зада­ ваемых распределений скоростей, параметров, определяющих меридиональные обводы, и другие факторы позволяют произ­ вести аналитическую оптимизацию ПЧ,

Глава 3 БЕЗЛОПАТОЧНЫЙ ДИФФУЗОР

3.1. Особенности течения, одномерный анализ

Безлопаточные диффузоры имеют простейшую геометриче­ скую форму. У промышленных центробежных компрессоров это обычно кольцевой канал, образованный плоскими радиальными стенками корпуса. Иногда эти стенки имеют криволинейную образующую или прямую, наклоненную по отношению к радиаль­ ному направлению. Выбор того или иного закона изменения ши­ рины диффузора b = f (lm) называют меридиональным профили­ рованием БЛД. Простота формы определяет простоту и нагляд­ ность одномерного анализа, расчета с использованием эмпириче­ ских коэффициентов.

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

соответственно:

 

или

(3.1)

где /Стр < 1 учитывает уменьшение момента количества движе­ ния из-за трения о стенки.

Отсюда

Для невязкого несжимаемого потока в БЛД постоянной ши­

рины

(/Стр =

1, е = const, В = const):

 

 

 

 

си =

фт/г;

сг = ф2;г;

с = сг/г; а =

а 2 = const.

(3.3^

Траектория потока в таком БЛД, так называемая логариф­

мическая спираль, — линия,

касательная

к

которой

в любой

точке

наклонена

под одинаковым углом

(90 — а)0 к

радиусу.

Из (3.2) следует, что увеличение плотности газа по радиусу умень­ шает угол а. Так же действует расширение БЛД в меридиональ­ ной плоскости. В обоих случаях сг быстрее уменьшается по ра­ диусу, чем в случае течения несжимаемого газа при b = f (г) = = const. Наоборот, торможение составляющей си касательными силами вязкости на стенках БЛД (/СтР < 1) и сужение диффу­ зора в меридиональной плоскости способствуют увеличению а по мере продвижения газа вглубь диффузора (речь идет об осредненном по b потоке, в действительности картина сложнее, так как силы вязкости действуют на частицы газа в разной степени в зависимости от их расстояния от стенок).

При поэлементных расчетах ПЧ потери в БЛД иногда реко­ мендуют определять как потери трения, причем коэффициенты трения К подбираются таким образом, чтобы расчетные потери соответствовали экспериментально определенным. Для оценки минимально возможных потерь в зависимости от параметров по­ тока и диффузора величина к может браться из опытов с трубами. В силу переменности скорости и гидравлического диаметра по длине БЛД потери трения определяются интегралом

 

 

(3.4)

где

 

 

dl =

dr!sin a;

dr = 4//U = 2b.

Наиболее просты

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

БЛД с b = const в предположении а = const. Последнее часто принимается при расчете ступеней ПЦК, так как противополож­

ное влияние

вязкости и

сжимаемости на направление потока

в

БЛД примерно

взаимно

компенсируется. В этом случае с =

=

с/с2г=

и из

(3.4)

 

 

 

 

 

(3.5)

Обратно пропорциональная зависимость £тр от sin а и б очевидна. Интересно, что с ростом г3 нарастание потерь постепенно замед­ ляется. При гг ОО

m