
книги из ГПНТБ / Митрохин В.Т. Выбор параметров и расчет центростремительной турбины на стационарных и переходных режимах
.pdf
|
|
|
|
|
|
|
|
|
|
Таблица 2.2 |
|
|
|
|
|
|
|
|
|
|
|
|
ю |
|
15° |
62°36' |
62,6° |
69° |
68° |
126,5° |
142,5° |
147,5° |
151° |
128,7° |
128,7° |
|
15° |
62°36' |
63° |
69° |
71,1° |
141,4° |
145,7° |
147° |
148,1° |
128,7° |
128,7° |
|
15° |
62°3б' |
63,6° |
69° |
75,2° |
105,9° |
138,8° |
143,1° |
144,7° |
128,7° |
128,7° |
|
15° |
62°36' |
64,2° |
69,2° |
80,8° |
104,5° |
134,7° |
138,4° |
139° |
128,7° |
128,7° |
|
15° |
62°36' |
65, Г |
70,2° |
88,3° |
104,5° |
127,8° |
130° |
129,7° |
128,7° |
128,7° |
|
15° |
62°36' |
66° |
74° |
97° |
104,5° |
114,5° |
115° |
115° |
128,7° |
128,7° |
|
0° |
31° |
29° |
28° |
3,5° |
—8,5° |
—12° |
-14° |
-16° |
0° |
0° |
|
0° |
31° |
29° |
27° |
5° |
— 6,5° |
- 1 0 , 5 ° |
-13° |
-14° |
0° |
0° |
|
0° |
31° |
29° |
25° |
10° |
- 5 ° |
—10° |
-10° |
-10° |
0° |
0° |
|
0° |
31° |
28,5° |
22° |
18,5° |
—2,5° |
—7,5° |
10° |
- 7,5° |
0° |
0° |
|
0° |
31° |
28° |
18° |
2° |
- 6 ° |
—12° |
-13,5° |
-3° |
0° |
0° |
|
0° |
31° |
17,5° |
16,5° |
0° |
—12° |
- 2 5 ° |
-20,5° |
0° |
0° |
0° |
0 |
|
0,8986 |
0,8921 |
0,8854 |
0,879 |
0,8736 |
0,8694 |
0,8673 |
0,8675 |
|
|
1 |
1 |
0,89,86 |
0,8907 |
0,8825 |
0,8744 |
0,8665 |
0,8585 |
0,8545 |
0,8535 |
|
|
2 |
1 |
0,8986 |
0,8892 |
0,879 |
0,8678 |
0,8557 |
0,8426 |
0,8367 |
0,8362 |
|
|
3 |
1 |
0,8986 |
0,8872 |
0,875 |
0,8585 |
0,8405 |
0,815 |
0,8062 |
0,805 |
|
|
4 |
1 |
0,8986 |
0,8852 |
0,87 |
0,8468 |
0,8186 |
0,7717 |
0,7473 |
0,7505 |
|
|
5 |
|
0,8986 |
0,8835 |
0,8645 |
0,8325 |
0,783 |
0,6815 |
0,599 |
0,566 |
|
|
стеснения % в |
узлах |
исходного |
приближения |
приведены в |
||
табл. 2.2. |
|
|
|
|
|
|
Из табл. 2.2 видно, что лопатки имеют сложную |
пространст |
|||||
венную форму: лопатка на входе |
в колесо |
имеет |
нерадиальное |
|||
направление р' = 62°36'. На входе |
в колесо |
лопатки |
имеют наи |
|||
больший положительный угол с радиусом |
(б'=31°). |
Между ли |
||||
ниями U и h величина |
б' меняет знак. Наибольшие |
отрицатель |
||||
ные значения 6' = —25° на линии |
k и далее, так как Еыходная |
|||||
кромка лопатки |
расположена радиально, стремятся |
к нулю. Ф,и- |
Li 1
i !
|
Ч |
|
Ч |
I I I ; |
I I I |
x 1,0 |
0,5 |
|
W |
|
|
WA |
Щ |
wB |
|
2 |
|
|
|
Wj |
|
|
|
|
Г |
|
|
|
I |
\ |
|
|
1 11 |
\ |
|
|
\ |
|
|
|
1 |
|
|
|
1 |
чч |
|
|
1 |
v |
|
> |
, |
I I I I I I i I I I I |
|
7 v g - j ^ - r > j |
|||
|
0 |
0,5 |
1,0 x |
Рис. 2.32. Распределение скоростей на профиле
зическая толщина профиля выбрана одинаковой и равной А = = 3 мм, поэтому степень загромождения (коэффициент стеснения у) изменяется по мере уменьшения радиуса от 0,898 до 0,566. Профили рабочих лопаток на выходе изогнуты так, чтобы •обес печить следующие значения углов потока р2 :
впериферийном сечении — 23°; •в среднем сечении — 26°;
вкорневом сечении — 28°.
Расчет проводился при расходе воздуха G = 0,035 кг/с и уг ловой скорости вращения со = 6700 1/с.
В приведенном предварительном расчете по одномерной тео рии получены следующие значения скоростей на входе и выходе из рабочего колеса: W\ = 132,5 м/с, ш 2 = 150,5 м/с, т. е. осущест вляется разгон потока в относительном движении. Предваритель но протяженность меридианного профиля вдоль оси (см. рис. 2.31) была Еыбрана исходя из рекомендаций работы [20]. Затем
было рассчитано распределение |
скоростей |
по обводам профиля, |
|
считая весь поток, как одну струйку тока |
(рис. 2.32). Диффузор |
||
|
ов |
А. |
|
ныи участок |
распределения |
проверялся по формпарамет- |
|
ру /, который |
получился равным / = 0,624, т. е. в среднем на про- |
||
|
|
|
71 |
филе при числе лопаток рабочего колеса z=12 должно осуще ствляться безотрывное течение.
Таким образом, рассматриваемое колесо было спрофилиро вано в соответствии с известными рекомендациями, основанными на приближенном расчете потока.
S 4 S3 S2 S7 SQ
-1о
•Li
•I2
Расчет осесимметричного те чения по методике, изложенной в разд. 2.5.1, показал следующее.
На рис. 2.33 приведены линии токов последовательных прибли жений. Начиная с третьего при ближения, в расчет был введен коэффициент релаксации <х=0,5. Всего для необходимой сходимо сти (2—3% но скорости) потре-
Рис. 2.33. Линии тока последовательных приближений
72
бовалось выполнить 7 приближений. Начиная со второго при ближения, на внутреннем обзоде меридианного профиля в райо не линии k получались отрицательные значения меридианных составляющих скорости ws. При возникновении отрывного тече ния на первый взгляд кажется, что начальная скорость, удовлет воряющая заданному расходу и рассчитанному градиенту ско ростей, может быть выбрана произвольно. В зависимости от выбранного значения скорости будет получаться разная протя-
W,M/C
1 I
> b
-г О |
0,5 |
|
|
|
1 |
|
5) |
Рис. 2.34. Распределение скоростей вдоль линий /
женность зоны отрицательных скоростей (зоны отрыва). Однако это противоречит единственности решения. Очевидно, что наибо лее устойчивым, и значит, осуществимым будет такое решение, при котором зона отрыва минимальна и, следовательно, отрица тельное значение скорости минимальное. В таком случае безраз лично, на каком из обводов меридианного профиля задавать начальную скорость, ее выбор здесь диктуется условием полу чения минимальной отрицательной скорости. Определенная из этого условия зона отрыва потока показана на рис. 2.33. Гра ница зоны определена при условии ш8 = 0. Эта граница, естест венно, условна: в реальных условиях вверх по течению предвари тельно наступит отрыв пограничного слоя и граничная линия тока будет плавной (без угловых точек).
Распределение относительных скоростей конечного прибли
жения вдоль линий / приведено на рис. 2.34, а |
и б. За |
исключе |
||
нием линии k |
(предотрывная зона), |
распределение |
скоростей |
|
близко к линейному. |
|
|
|
|
Представляет интерес рассмотрение изменения скорости по |
||||
I в начальных |
участках меридианного |
профиля |
(линии h и <_). |
73
В рабочих колесах с радиально расположенными входными кромками вследствие малой кривизны меридианного профиля на
начальном участке скорости w практически |
не изменяются |
вдоль /. В рассматриваемом случае нерадиально |
расположенных |
входных кромок (р' = 62°30/ ) возникновение градиента скорости на начальном участке обусловлено влиянием воздействия лопа ток на поток ( б ' ^ О ) .
w, м/с
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1/ |
|
|
|
|
300 |
|
|
1 |
|
|
|
; |
|
|
|
1 |
|
и |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1 |
|
|
|
11 |
|
|
|
|
|
200 |
|
|
1 |
|
|
|
|
|
|
|
/ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
/ |
|
|
|
|
|
|
|
|
/ |
|
|
|
|
|
|
|
|
/ |
|
|
|
|
|
|
|
|
/ |
|
|
|
|
|
|
|
|
/ |
|
|
|
|
|
100 |
|
|
. |
|
|
|
|
|
|
|
|
|
I |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ч - |
|
|
|
О |
0,1 0,2 |
0,3 |
OA |
0,5 |
0,7 |
0,8 |
0,9 / /s |
|
|
||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
I |
Рис. 2.35. |
Распределение |
скоростей |
по обводам |
меридиан |
||||
|
|
|
ного |
профиля: |
|
|
||
|
/ — по внутреннему обводу; / / — |
по внешнему |
обводу; |
|||||
|
по методу |
последовательных приближений; |
|
|
||||
|
по |
формулам (2.61) |
|
|
|
|
|
Распределение скоростей по внутреннему и внешнему обво дам меридианного профиля показано на рис. 2.35 (сплошные линии). На этом же рисунке (штриховые линии) даны распреде ления скоростей, рассчитанные приближенным методом по мето дике, изложенной в разд. 2.5.2. По оси абсцисс отложено без размерное расстояние по контуру обводов в долях от длины дуги внутреннего обвода. Можно отметить только качественную схо димость точного и приближенного методов расчета. Поэтому при расчетах последовательных приближений в качестве исходного следует рекомендовать использовать данные приближенного рас-
74
чета. Естественно, в зонах больших углов б', что не учитывается приближенной методикой, расхождения наибольшие (начальный участок и зона отрыва). При этом важно, что приближенная ме тодика расчета не указывает возможность отрыва невязкой жид кости, поскольку величины скоростей по этому методу всюду со храняют положительные значения. В этом смысле можно гово рить и о качественном различии точного и приближенного мето дов при расчете течений с отрывом потока.
w
|
|
|
WA |
2 |
|
WB |
// |
|
|
|
|
|
|
||
|
|
|
|
fЛЛ |
|
||
|
|
|
|
1 |
\ЛV |
|
|
|
|
|
|
1 |
\ |
|
1 |
- |
|
|
|
|
|
|
1 |
|
|
\ \ ч |
|
|
> |
1 _ —-- |
|
|
|
к' ч |
|
|
ч_:>< |
||
|
|
Чч |
0-\ |
|
|
|
|
|
|
ч |
ч ч_,.л |
|
|
|
|
л 1,0 |
0,5 |
|
0,5 |
1,0 /. |
|||
Рис. 2.36. Распределение скоростей по обводам профиля: |
|||||||
— |
• — |
по данным расчета последовательными приближениями; |
|||||
|
|
по приближенной |
методике |
|
|
||
После |
того |
как были определены |
линии |
тока, по методике, |
изложенной в разд. 2.5.3, было рассчитано распределение скоро
стей в средней струйке тока |
по обводам |
профиля |
лопатки. |
|||||
|
На рис. 2.36 распределения скоростей |
по обводам профиля |
||||||
(штрих-пунктирные |
кривые) |
сопоставлены |
с |
данными |
расчета |
|||
по |
приближенной |
методике |
(штриховые), |
когда |
осесимметрич- |
|||
ное |
течение не рассчитывалось и весь |
поток |
считался |
одной |
||||
струйкой тока. В зоне отрывного течения в меридианной |
плоско |
сти (х>0,5) наблюдаются качественное и количественное рас
хождения в распределениях скоростей по |
обводам профиля. |
Если при расчете среднего течения кривая |
Wo |
плавно умень- |
шалась, то в уточненном расчете появился резкий минимум. Уве личился отрицательный градиент скорости и проверка по формпараметру показала, что течение по обводу профиля лопатки близко к отрывному (/ = 1,0).
Из рассмотрения этого примера можно сделать некоторые об щие выводы.
1. Известные приближенные методики расчета и профилиро вания рабочих колес не всегда пригодны для профилирования колес с загнутыми на входе лопатками. При профилировании
75
таких колес обязателен подробный расчет осесимметричного те чения.
2.Характер осесимметричного течения в значительной степе ни влияет на течение по осесимметричной поверхности тока (от лопатки к лопатке).
3.Приближенный расчет распределения скоростей по обво дам меридианного профиля качественно согласуется с более точ ным для рабочих колес сложной формы. Результаты этого рас чета могут служить лучшим исходным приближением для последующих уточнений, чем обычно применяемый расчет одно мерного течения, когда струйки тока разбиваются по принципу равных кольцевых площадей.
Помимо рассмотренного выше отрывного течения от обводов меридианного профиля, возможно появление отрицательных зна чений скоростей по контуру профиля. Заранее указать такие зо ны без расчета по методике, изложенной в разд. 2.5.1, нельзя, особенно для рабочих колес с загнутыми на входе лопатками. Не приводя примеров таких расчетов, отметим, что для ликви дации отрывных зон на профиле наиболее радикальным сред ством является увеличение числа лопаток. Причем, если по ус ловиям размещения в диске такое увеличение числа лопаток нежелательно, следует часть лопаток (половину) распространять не на всю длину меридианного профиля.
2.8. РАСЧЕТ КОЭФФИЦИЕНТОВ ПОТЕРЬ В РАБОЧЕМ КОЛЕСЕ РАДИАЛЬНОЙ ТУРБИНЫ
При расчете коэффициентов потерь в рабочем колесе прихо дится учитывать режимы работы радиальной турбины. На рас четном режиме, как правило, поток входит в рабочее колесо без угла атаки, и потери в рабочем колесе определяются в 'основном трением в пограничном слое. При возникновении отрыва угол входа потока на колесо отличается от расчетного, и потери в рабочем колесе в основном связаны со срывным течением.
Приводимая ниже методика приближенной оценки коэффи циентов потерь на расчетном режиме может быть использована для любых Т И П О Е колес центростремительных турбин, а также турбин турбомуфт при дозвуковой скорости течения.
Методика оценки потерь трения в рабочем колесе радиальной турбины, как и описанная в разд. 2.3 методика расчета потерь трения в сопловых аппаратах, базируется на известных величи нах коэффициента потерь трения в плоском потоке. Зная его, можно по формуле (2.34) рассчитать коэффициент потерь трения на профиле рабочего колеса радиальной турбины. В формулу (2.34) входила величина й;Сп.к средней скорости в круговой непо движной решетке. При расчете коэффициента потерь трения в рабочем колесе в эту формулу следует подставлять величину
76
и>сп.к Д л я вращающейся решетки (2.14) или для несжимаемой жидкости (2.17).
Если геометрические параметры (s, рг, pi) и числа Re в плос ком потоке и в потоке на криволинейной осесимметричной по верхности тока рабочего колеса одинаковы, то коэффициенты потерь трения на Профиле рабочего колеса £Тр.р и в плоском по токе Стр.п связаны так:
|
|
|
Ц>.р — ^ т р . п |
( |
zr^ |
) |
~ ^ т р . п |
( •••с " - ) . |
(2.8 г) |
|||
т. е. относятся |
между |
собой |
|
|
|
|
||||||
примерно |
как кубы |
средних |
|
|
|
|
||||||
скоростей |
или |
средних |
|
без |
|
|
|
|
||||
размерных |
скоростей |
|
для |
|
|
|
|
|||||
сжимаемой |
жидкости. |
|
По |
|
|
|
|
|||||
формуле |
|
(2.84) |
можно |
|
рас |
|
|
|
|
|||
считать |
|
потери трения |
на |
|
|
|
|
|||||
профиле |
|
радиальной |
решет |
|
|
|
|
|||||
ки при условии, что течение |
|
|
|
|
||||||||
в пограничном слое близко к |
|
|
|
1,75 ji |
||||||||
плоскому. |
|
На |
возможность |
|
|
|
||||||
|
|
|
|
|
||||||||
такого |
приближения |
было |
Рис. 2.37. Зависимость А,СПА сп.п ОТ |Л И О |
|||||||||
указано |
в разд. 2.6. |
|
|
|
||||||||
|
|
|
при |
^2 = 0,5; |
s^=3,65; Т=3,5; а2 =Р,=90°; |
|||||||
На рис. 2.37 в качестве |
||||||||||||
|
|
р 2 = з о ° |
|
|||||||||
примера |
|
приведена |
зависи |
|
|
|
\х и сте |
|||||
мость величин |
ЯспАсп.п от коэффициента |
радиальности |
пени уширения профиля б при характерных для решеток рабочих
колес центростремительных |
турбин |
параметрах: Я2 = 0,5; s = 3,65; |
|||
а 2 = 90°; Pi = 90° и |
р2 = 30°. При |
ц<\ |
и б > 1 график на |
рис. 2.37 |
|
соответствует центростремительным турбинам, при ц>1 |
и 6<1 —• |
||||
центробежным, а |
при |
и |
6—1 — осевым турбинам. Из |
рис. 2.37 видно, что коэффициент потерь трения в решетке рабо
чего колеса центростремительной турбины существенно |
больше, |
а в центробежной меньше, чем в осевой. |
|
Мы уже отмечали, что несмотря на большие величины коэф |
|
фициентов потерь и, следовательно, меньшие значения |
в цент |
ростремительной турбине, абсолютная величина потерь в рабо чем колесе центростремительной турбины меньше, чем в осевой. Это объясняется тем, что коэффициент потерь отнесен к энергии выходной скорости, составляющей в центростремительной тур бине малую долю теплоперепада, срабатываемого в рабочем ко лесе.
Коэффициент потерь |
|
|
C = l - f - |
2 С — 2 - ( " ! - » ! ) |
(2.85) |
|
77
Рассмотрим теперь, как влияет величина скорости вращения рабочего колеса на величины средних безразмерных скоростей, т. е. в конечном счете на потери в рабочем колесе.
На |
рис. 2.37 |
мы |
рассматривали зависимость А,СпАсп.п П Р И а 2 = |
|||||
= 90°; |
на рис. 2.38 приведена зависимость А,спДсп.п |
от ц и 6 при |
||||||
тех же |
|3ь р2 и |
s, что на рис. 2.37, но на остановленном |
колесе |
|||||
( а 2 = Р г ) . Из |
сопоставления |
данных, |
приведенных |
на рис. 2.37 и |
||||
2.38, видно, |
что |
отношения |
ЯопАсп.п |
и, следовательно, |
коэффи |
|||
циенты |
потерь |
при |
02 = 90° и а2 = р2 |
существенно |
различаются. |
Это свидетельствует о невозможности непосредственного исполь зования данных продувок неподвижного колеса в расчете потерь во вращающемся рабочем колесе.
|
Перейдем к |
последователь |
|||||
ности |
расчета |
коэффициентов |
|||||
потерь трения в рабочем коле |
|||||||
се |
радиальной |
турбины. |
|
|
|||
|
Для расчета |
должны |
быть |
||||
предварительно |
определены |
||||||
следующие |
величины: |
|
|
||||
|
— угол входа потока |
в |
ра |
||||
бочее |
колесо рь |
|
|
|
|||
|
— угол |
выхода |
потока |
из |
|||
рабочего колеса р2 ; |
радиально |
||||||
|
— коэффициент |
||||||
сти |
ц; |
|
|
|
|
|
|
|
— степень уширения |
профи |
|||||
ля |
проточной части |
6; |
|
z; |
|||
|
— число |
лопаток |
колеса |
||||
— длина |
скелетной |
линии |
|||||
профиля s. |
|
|
|
|
|
||
|
Ступень |
турбины должна |
|||||
быть |
приближенно |
рассчитана |
|||||
исходя из предварительно выбранного |
коэффициента |
потерь |
£. |
||||
В результате расчета величина £ уточняется |
и |
если |
различие |
между предварительно выбранным и полученным расчетом зна чениями £ велико, расчет ступени турбины проводится заново.
По данным, приведенным на рис. 2.7, исходя из заданных значений Pi и Рг, определяется величина коэффициента потерь
трения в плоском потоке Сгро = V ^ l — ?о • По формуле (2.37) вычисляется оптимальный относительный шаг и по среднему ша гу рабочего колеса i C p = (^1+^2)/2 — его средний относительный шаг. Если средний относительный шаг отличается от оптималь ного, величина фо корректируется по графику (см. рис. 2.9). Для решеток рабочих колес центростремительных турбин эта поправ ка, как правило, достаточно велика, так как при большой длине скелетной линии профиля s решетка получается очень густой.
78
|
В итоге получим коэффициент потерь трения в плоском пото |
|||||||||
ке |
С1.р.п = ]/ Л 1 — <?2 |
при |
таких же |
значениях густот профиля s |
||||||
и |
углов j3i и Рг. как |
и в |
решетке |
рабочего |
колеса |
радиальной |
||||
турбины. |
|
|
|
|
|
|
|
|
|
|
|
Если скорости потока |
малы (М„-2 <0,4), то по формуле |
(2.17) |
|||||||
вычисляется средняя |
относительная |
скорость в |
решетке |
колеса |
||||||
й7Сп- Если скорости потока велики |
( М № 2 > 0 , 4 ) , |
то |
вычисляется |
|||||||
средняя приведенная скорость Я с п |
по формуле |
(2.14). |
|
рассчиты |
||||||
|
Полагая в формуле (2.17) или |
(2.14) р,= 1 и 6 = 1 , |
||||||||
ваем средние скорости в плоском потоке при тех же |
|3i, |
р2 |
и 1 , что |
|||||||
и в решетке рабочего колеса радиальной турбины. |
|
|
|
|||||||
|
Если числа Re в плоском потоке и в потоке на |
криволинейной |
осесимметричной поверхности одинаковы, то коэффициент по терь трения в решетке рабочего колеса радиальной турбины £Тр.р рассчитывается по формуле (2.34) или (2.35).
2.9. ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ПОТЕРЬ
,В РАБОЧЕМ КОЛЕСЕ ПРИ УГЛАХ ВХОДА ПОТОКА,
ОТЛИЧАЮЩИХСЯ ОТ РАСЧЕТНОГО
По приведенной выше методике можно оценить коэффициент потерь на расчетном режиме работы ступени, когда не возникает, отрыва потока.
При расчетах коэффициента потерь на режимах, когда воз никают потери, связанные с отрывом потока от кромок рабочих лопаток, неполным заполнением межлопаточного канала и после
дующим выравниванием (смешением) |
потока можно |
получить |
[32] аналитическую зависимость для |
коэффициента |
так назы |
ваемых ударных потерь, который ранее определялся только эм пирически [2].
Рассмотрим задачу выравнивания плоского течения несжи маемой жидкости через плоскую решетку радиальных пластин, вращающихся с угловой скоростью Q.
Поток, входящий на решетку под углом, отличным от ради ального направления, отрывается на кромках пластин. Под влия нием внутреннего трения поток на некотором радиусе выравни вается; на этом радиусе направление потока совпадает с ра диальным, а полное давление становится постоянным. Трением маловязкой жидкости о пластины при отрыве потока можно пре небречь или учесть его влияние отдельно.
Для плоской прямой решетки, которую можно рассматри вать как предельный случай радиальной решетки в относитель ном движении при бесконечном увеличении радиуса и одновре менном увеличении числа пластин, аналогичные задачи (в более общем случае сжимаемой жидкости) были решены теоретически Г. И. Тагановым [49] и Г. Ю. Степановым [48] исходя из основных уравнений сохранения.
79