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

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

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

сгвителЬного Характера Движения газа в зазорах от имеющего место в опытах с диском, вращающимся в закрытом корпусе.

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

и

опытных данных

в [65] дается рекомендация принимать 0,1 IV

=

0,11-10"3. Однако

вцелом формула (1.62) не может служить основой для определения потерь дискового трения для произвольного РК, так как ее струк­ тура не соответствует физической картине движения газа в зазорах.

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

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

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

максимального расхода,

когда при

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

за ступенью статическое

давление

в РК может не возрастать,

а уменьшаться). У основного диска возможны такое же (для кон­ цевой ступени) и противоположное (для промежуточной ступени) направления течения. Поэтому при расчете потерь трения дисков необходимо учесть работу, получаемую газом от РК и вносимую в зазор в случае движения от периферии к центру. При противо­ положном направлении течения следовало бы учесть работу, полу­ чаемую газом в зазоре от диска, а затем возвращаемую основному потоку. В первом случае происходит возврат части энергии рабо­ чему колесу, что уменьшает мощность трения по сравнению со случаем отсутствия радиального течения. При си < и газ отстает от вращающейся поверхности и возникающая сила трения, про­ порциональная квадрату относительной скорости wu = си — и, тормозит рабочее колесо; при си > и газ как бы обгоняет рабочее колесо, и в этом случае сила трения создает момент, совпадающий с направлением вращения ротора, т. е. мощность трения отрица­ тельна; в зазоре газ вращает колесо за счет энергии, полученной от него же ранее. Подчеркнем, что сила трения на элементе вра­ щающейся поверхности пропорциональна квадрату скорости дви­ жения газа относительно этого диска, т. е. скорости wU9 а не пере­ носной скорости и, как это принято в (1.60). Так как мощность и сила трения должны быть положительны при wu < 0 (за поло­ жительное принято направление и), в формуле для элементарной силы трения введен знак «минус», а для учета возможного изме­ нения знака wu — ее модуль

dNTP= и dPTр = —Àjp0,51wuI wuu dF,

(1.63)

где — коэффициент трения, связывающий силу трения со сред­ ней по ширине зазора кинетической энергией газа, вращающегося относительно РК; dF — площадь элемента поверхности (dF — = 2nrdr — для радиальной поверхности).

Для вывода расчетных уравнений мощности трения рассмотрим элемент газа высотой dr, вращающийся на радиусе г в зазоре ме­ жду рабочим колесом и корпусом с некоторой скоростью си. На боковые поверхности выделенного элемента действуют силы трения о неподвижную и подвижную поверхности зазора. Силы сопроти­ вления поверхностей по аналогии с (1.63) таковы:

dP„ — 0,5Àip | wu I wu dF\ dPB —

dF.

По теореме об изменении момента количества

движения эле­

мента массы dm — р2nrdrB (В — зазор между корпусом и диском)

d m dJÿ ~ = d M г = ~ d P »r + d P nF

где dm/dT = mnp — массовый расход протечек через уплотнение. Подставляя значения сил трения и приводя предыдущее выра­

жение к безразмерному виду, получим

dcu =

or [—cl + \cu — r\(cu — r)]d r ---- cf- - dr,

(1.64)

где a = \г /Ф пр

комплекс безразмерных величин,

входящих

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

/ (г); Фпр = Ф|Зпр.

 

Уравнение (1.64) — дифференциальное уравнение первого по­ рядка, решением которого на ЭВМ методом Рунге — Кутта полу­ чены зависимости cu = f(r, à) при задании граничных условий: си — °u2 = Фт (ел = 0) при г — 1 — движение газа от периферии; си = 0 при г = гп — движение газа к периферии в зазоре у основ­ ного диска.

Решения са — f (г, фт, а) и си / (г, а, гя) получены для диапа­ зонов фт = 0,2 ч-1 и а = 2-н20 при течении к центру и при гл= = 0,3 и а = 0,5^-100 при течении к периферии. Графические зависимости представлены на рис. 1.9 и 1.10.

Зависимости cu — f (г, фт, а) подтверждают предположение о возможности возврата энергии рабочему колесу. Прямая линия

на графике рис. 1.9

изображает изменение окружной

скорости,

В зависимости о т а в

пределах г = 0,3 -н0,65 при фт =

0,6 и при

больших значениях г с возрастанием фт систановится

больше й,

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

обработать

более тщательно (что

эквивалентно уменьшению а),

то начало

возникновения эффекта

будет происходить раньше,

а значит, и большая доля энергии затем будет возвращена колесу. В той части диска, где си > и, его поверхность можно сделать более

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

Интегрируя выражение (1.63) в предположении а = const, найдем коэффициент потерь трения радиальной поверхности од­ ного диска

 

гл

 

Ртр1==

I r — cu\ ( r ~ c u)P dr.

(1.65)

Рис. 1.9. Изменение окружной составляющей скорости в зазоре у покрывающего диска

Введем в рассмотрение условный коэффициент потерь трения

ВТР= Р T P l

Ф'фт

Ртрхфт

r — cu\(r — cu) P d f.

(1.66)

^1®

РпрЯ

 

 

 

 

Как видно из (1 .66), ВтР зависит от и от распределения

скорости

си = / (г),

которое определяется значениями а и i|\r (при

течении

от периферии).

 

 

Полученные в результате интегрирования (1.66) зависимости

Втр. п =

/ (?л. Фт>а) Для покрывающего диска (течение от перифе­

рии) графически представлены на рис. 1.11. Зависимости 5 TPl0 — = f (/‘л, а) (течение к периферии) показаны на рис. 1 .12 . Обращает внимание существенное влияние гл. Это связано с тем, что в усло­ виях радиального течения газа в зазорах вид зависимости си =

=f (г) не совпадает с неявно принятым при выводе формул (1.60)

и(1.62) соотношением си = 0,5ы = 0,5со/\ Отметим, что при Рпр = 2

2 Селезнев К. П.

33

=

0

= оо) это условие должно выполняться;

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

( 1 .6 6 )

дает одинаковые значения 5 тр. 0 =

£ тр. „ =

0 ,0 5 (при гл=

=

0). Видно,

что при

рпр Ф 0 потери трения

на

основном диске

при течении

к центру

несколько больше, чем 0 , 0 5 ,

причем они

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

росте протечек

(малые

а),

что связано

с

подводом

энергии

к

протекающему

газу

(как

указывалось

в п. 1 .3 , соответствующий напор не^является полностью потерян­ ным). Зато при типичных а «=* 5 - ь 1 0 при г л != » 0 ,6 и грт = 0 ,6 - 4 - 0 ,8 значения Втр. п я» 0 ,0 3 - т -0 ,0 1 , т. е. в среднем в 1 ,5 — 5 раз меньше,

Рис. 1.10. Изменение окружной составляющей скорости в зазоре у основного диска — течение к периферии

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

ного диска (концевая ступень) значения ВтР

0,3 -4-0,4, т. е.

меньше 0,5.

Значение ВтР позволяет определить коэффициент дискового трения РК по соотношению, вытекающему из (1.61), (1.65) и (1 .66):

Ртр = 1,1(Втр. 0 + ВтР. п) ^ г>

(1.67)

где коэффициент 1,1 по аналогии с (1.62) приближенно учитывает

трение цилиндрических участков дисков на г2, а е

е2 в соответ­

ствии с обычно принимаемым условием р <=* р2 [60; 65].

Для использования формулы (1.67) необходимо

знать вели­

чины Хх и рпр, от которых, в частности, зависит коэффициент а, определяющий Втр. Для расчета Рпр рекомендуется соотношение (1.79). Впредь до получения исчерпывающих экспериментальных

данных для

^ рекомендуется

следующее ориентировочное соот­

ношение:

К =

4V,

(1.68)

 

где X' соответствует значениям /Стр =

лХ’102 по данным Цумбуша,

приводимым

в работе [65).

 

 

Соотношение (1 .68) вытекает из равенства сил трения, опре­ деляемых через X и %' при отсутствии радиального течения, когда w„ = си = 0,5мг:

О.брА,!^ dF = 0,5pV«2dF.

Для удобства использования соотношения (1.68) эксперимен­ тальные данные Цумбуша, приводимые обычно в графическом

виде [60; 65], представлены ниже

в виде аппроксимирующих

алгебраических

выражений:

 

R e „ ^ 6 -105

 

 

 

 

 

 

 

для

гладких

дисков при

 

 

 

 

 

 

 

 

X’ =

(10'2/л) (2,12 -

0,231 lg Re*)

при

B/D2=

0,2;

 

 

 

X' =

(10-2/я) (1,96 -

0,215 lg Re„)

при

B/D2 =

0,1;

( 1 .6 9 )

 

 

X’ = (10~2/я) (1,51 — 0,154 lg Re„) при

B/D2 = 0,015;

 

 

для

гладких

дисков при

Re„ < 6 -1 0 8

 

 

 

 

 

 

 

 

X' =

(10-2/я) (27,5/lg Re„ — 4)

при

B/D2 — 0,2 ;

 

 

 

 

X' =

(10~2/я) (24/lg Reu — 3,46)

при

B/D2 =

0,1;

 

(1.69а)

 

 

X’ = (10-2M)(20,4/lgReu - 2 ,9 ) при

fi/D2 =

0,015; ,

 

для

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

при Re„ s? 6-10®

 

 

 

 

 

 

X' = (10~2/я) (6,65/lg Rea -

0,35)

при

B/D2= 0,155;

 

 

 

X' =

(10_2/JT) (5,54/lg Re„ — 0,3)

при

B/D2 = 0,01.

(1.696)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Косвенная оценка значений

при условиях

радиального те­

чения от периферии с варьированием

/гёпр и си2 была

проведена

 

____________________

 

по

распределению давлений

в зазо-

Tpov

I

Г

I

 

рах у вращающихся дисков РК мо­

0,06

 

 

 

 

 

дельных

ступеней

[(Р.

 

Набавани,

 

 

 

 

 

ЛПИ) и на специальной

эксперимен­

 

 

 

 

г„-0,2

 

тальной установке (С. С. Евгеньев,

0,050

 

 

 

/ ,0,3

 

СКВ-К).

Измеренные Ар = р2 рл

 

 

 

 

П ь ч

 

позволяют,

используя

 

выражение

 

 

 

 

 

 

(1.78), определить зависимость си =

0,05

 

40

60

80

 

=

/ (г),

определяемую

величиной

го

 

а, т. е. величиной 5^, при прочих

Рис.

1.12. Условный коэффици­

 

заданных параметрах. При характер­

ент

трения диска

при течении

 

ных для РК ПЦК значениях

Фдр =

 

 

к периферии

 

=

Рпрф

 

-*-2) 10~4 численное зна-

 

 

 

 

 

 

чение*

оказалось

практически не

зависящим от этого параметра. Зато, по предварительным дан­

ным, значение фт на величину

Х1

влияло

очень

сильно: диапа­

зону фт = 0 , 5

-г -0 ,9 (М„ = 0 , 6

и

0 ,7 5 ) соответствовали

пример­

ные значения

^ от 4V до 0,5А/.

 

Xt ^

f (фт)

рекомен­

Впредь до

подтверждения

зависимости

дуется соотношение (1.68), возможно завышающее ртр для высоко­ напорных РК.

Соотношение между потерями на трение дисков по предлагае­ мой методике (1.67) и по обычно рекомендуемому соотношению (1.62)

Ртр/Р;р = 2 ,5 ( В тр. п + ВгР. о) ( W ,

(1 .7 0 )

где как ртр обозначена величина, определяемая по (1.62).

Для расчета Ртр величины ВтР. 0 и ВтР. п могут быть определены интерполяцией по графикам на рис. 1.11 и 1.12. Аналитически значения Втр при характерных величинах гл «=* 0,6 для покры­ вающего и 0,3 для основного дисков даются следующими зависимостями, аппроксимирующими результаты расчетов по уравнению (1 .66):

 

 

Втр.п =

Ю' 4(—956 +

147а - 7,2а2 +

0,12а3) +

 

 

 

 

 

 

+ ( 1IVT?) (0,22/У 5 - 0 ,1 3 /а -

0,03);

 

(1.71)

 

 

 

 

Втр. о =

0,043+ 0,086/|/5.

 

 

 

 

 

 

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

аппроксимирующие зависимости:

 

 

 

 

 

 

 

Втр. п =

0,09 -

0,ЗД /5 + гл (0,4/а2 -

0,0075) +

 

 

 

 

+ 1/К ф Д 0,16/У 5 - 0 , 0 2 - г л (0,00035-0 ,3 5 /а 3)];

 

 

 

Втр. 0=

10- 1

[0,442 — 0,071гл + (0,857 +

0,052гл) а-0,5].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.72)

2

Полиномы (1.72)

аппроксимируют функцию

Втр. п в областях

а ^

20;

0,2 ^

фт ^

1; 0,2 < гл ^ 0,7,

а

5 TP. 0

при

20 ^

^

а ^

100;

0,2

гл ^

0,5 с погрешностью

не более 5 %.

 

 

 

К достоинствам описанного метода расчета ртр следует отнести

количественный

учет различий в условиях течения

газа

у

по­

верхностей дисков. Например, можно рассчитать уменьшение

ртр

при износе уплотнений

(рост

|Зпр — уменьшение а), которое упо­

мянуто в работе

[60].

подход

к расчету течения газа

в зазорах

 

Изложенный

выше

между дисками и корпусом позволяет также уточнить перепад давления Арл при расчете протечек в лабиринтном уплотнении. Численное значение этих протечек определяется по известной формуле Стодола [60]

^ пр = [ж2лгл0дрСР J / 2 Д/7л/(2рсР)>

(1.73)

где р — коэффициент расхода, зависящий от конструкции лаби­ ринтного уплотнения; 0Л — радиальный зазор в уплотнении.

Из сравнения (1.73) с выводом той же формулы в книге [65] следует, что точное значение средней плотности газа в процессе дросселирования через лабиринты таково:

Рср ^ (Рд/^Рд) (Рд Н“ Ро)у

(1.74)

37

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

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

 

Рпр = (2fxrЛбл/Ф) J/2 eCPA p Jzt

(1.75)

где Арл =

ApJ(poul) — относительный перепад

давлений в ла­

биринтном

уплотнении.

 

В работе [60] изложена общепринятая методика определения перепада давлений в уплотнении, в которой предполагается, что газ в щели вращается с угловой скоростью, равной половине ско­ рости вращения рабочего колеса, прирост давления в РК опреде­ ляется только разностью и2 и иъ а плотности газа в РК, лаби­ ринте и зазоре одинаковы. В результате выкладок предложена формула для расчета коэффициента протечек, которая дана в при­ нятых в данной работе обозначениях:

Рпр = (2цгл0л/ф) 8СРY l 3/(42)] (1 — rf).

(1.76)

Принятые допущения, сделанные для исключения Арл из рас­ четных соотношений, должны приводить к заметным ошибкам, поскольку обычно w Ф 1 для РК ПЦК, что приводит к дополни­

тельному

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

выше,

си =

= / (г) в зазоре не соответствует и = 0,5сог.

 

 

Уточним перепад давлений в уплотнении следующим образом:

 

Д/>л = Рг — Ро — Ар,

 

(1-77)

где Д/7 =

— рл — падение давления при течении

газа

в за­

зоре от входа в него до уплотнения.

 

 

Интегрируя уравнение равновесия элемента в зазоре, можно

получить перепад давлений в

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

 

 

 

Ар = [Др/(Рср. кul)] =

J

(р/рср. к) {cl/r) dr =

 

 

 

 

1

 

 

 

 

— (рср/рср.к) j

{du/г) dr,

 

(1.78)

 

 

1

 

 

 

где рср. к — средняя плотность газа в рабочем колесе; рср — сред­ няя плотность газа в зазоре.

При определении Ар по (1.78) потеря давления от сопротив­ ления протеканию т пр через зазор считается пренебрежимо ма­ лой, а отношение плотностей можно принять равным единице.

Перепад давлений в рабочем колесе определим через давления по полным параметрам и скоростям в сечениях 0—0 и 2— 2

Р2 Ро ==zР2 Ро — рср. к0,5 (^2 — Со) ^ рср. к п — ^ д ) •

Используя понятие теоретической степени реактивности ра­ бочего колеса QT, определим динамический напор

0,5 ( с2 -

со) «

(1 - Пт) К

 

 

и, так как h*n = hTу\1у

 

 

 

 

{р%— Ро)/(рср. К^г) ^

Фт (^2 1 “Ь От).

 

Тогда перепад давлений на лабиринтном уплотнении

 

ДРл = Фт(*)2 -

1 + Q i) — Ар-

 

 

Формула для определения коэффициента протечек (1.75) при­

мет вид

 

 

 

 

Рпр — (2ргЛ0Л/Ф\ j / 2еср. у8Ср. к/^ЕфтСПг — 1

От) — ДрЬ (1 -79)

где р ^ 0,7 — для обычных ступенчатых уплотнений

[65];

еср. убСр. к = Рср. урср. к/(ро)2æ

1,05 -т- 1,15 — для

обычных

РК

ПЦК (произведение средних коэффициентов сжатия газа в уплот­

нении и РК); 1 — От ^ 0,5фт.

 

/ (г) можно использовать

Полученные выше зависимости^ =

при интегрировании (1.78). Результаты

расчетов Ар = / (гл,

фт, а) представлены на рис.

1.13 и 1.14, а в диапазоне параметров

гл =

0,2 ч-0,8; а =

2 4-20;

фт =

0,2 4-1

при течении к центру и

для

гл = 0,2 4-0,5;

а = 20 4-100

при

течении к периферии

ап­

проксимированы аналитическими зависимостями.

для

Для типичных

значений гл = 0,6

и

0,3 соответственно

покрывающего и основного дисков получены относительно простые соотношения соответственно для течения от периферии и от центра:

Ар = г|4'65 (0,91а -

0,2/а2 - 0,02) +

а/(61,7

39 /

а + 16,9а); )

 

 

Ар = 0,0048 (10 + У а).

 

( 1.80)

 

 

 

 

 

 

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

Ар = i]4-66 {ri'5/[а (0,18 + 1,56гл3)] -

гл/(9,11 +

Ю7,88г3/ ) +

 

+ a_2 (0,848r“ 1’3 -

1,84)}+

а/а,

(1.81)

где

 

 

 

 

 

 

а =

102[—0,64 + 8,69гл -

26,55г2 +

25,93г2 +

+

(0,19 -

4,72гл + 14,15г2 -

13, 16г2) у а+

 

-Ь (-0,034 -|- 1,12гл — 2,87?\ + 2,61 г») а];

Ар = 0,048 - 0,031 ri + у а (0,0068 -

0,0066гл).

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

Так как входящие в формулы (1.80) и (1.81) значения а завйсят от РПр, расчет по предложенной методике ведется итера­ циями (не более двух-трех). Расчеты^Рпр по (1.79) показывают,

что для типичных РК ПЦК протечки у покрывающих дисков со­ ставляют 0,7—0,8 от величины, даваемой соотношением (1.76).

Уточнение расчетов Рпр и ртр довольно заметно влияет на оценку эффективности РК и на точность расчета проточной части. Для примера на рис. 1.15 сравниваются значения коэффициента потерь РК унифицированной ступени ЛПИ типа 1а-0,25-0,06 (см. данные в приложении), рассчитанные по измеренным пара­