
книги / Центробежные компрессоры
..pdfсгвителЬного Характера Движения газа в зазорах от имеющего место в опытах с диском, вращающимся в закрытом корпусе.
Для улучшения совпадения расчетных |
и |
опытных данных |
в [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
где индексом л отмечены параметры перед уплотнением, а давле ние за уплотнением принято равным давлению в сечении 0—0.
Коэффициент протечек получен делением (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 (см. данные в приложении), рассчитанные по измеренным пара