
книги из ГПНТБ / Митрохин В.Т. Выбор параметров и расчет центростремительной турбины на стационарных и переходных режимах
.pdfстоятельство должно быть дополнительно учтено при выборе рас четной точки турбины.
Для того чтобы узнать, насколько изменяется величина ци при отклонении от оптимальных параметров, в частности, при переходе к значению а2 = 90°, проведем расчет г\и в зависимости от «! и р при |х = 0,5 и ai = р 2 = 15° (см. рис. 1.10). Одновременно этим расчетом проверим достаточность условий (1.21), (1.22) и
(1.23), как условий максимума т)„. |
|
|
|
|||||
Пи |
|
|
|
Согласно выражениям |
(1.21)г |
|||
9=0А |
|
|
(1.22) и (1.23) максимальное зна |
|||||
0,9 |
|
|
чение т]м при принятых |
JLI, |
а] и (32 |
|||
|
|
10° |
достигается при |
= 0,73, р = 0,47 |
||||
|
L |
|
||||||
|
|
|
и 02=143°. Рис. 1.10 подтвержда |
|||||
0,85 |
|
|
|
|||||
1 |
\ 9 - о^п |
ет, что условия |
(1.21), |
(1.22) и |
||||
|
Ч),6 |
(1.23) действительно являются до |
||||||
0,8 |
/ |
|
|
статочными условиями |
максиму |
|||
|
|
|
ма Г]и. |
|
|
|
||
0,75 |
|
|
Далее из графика следует, что |
|||||
|
Q,h 0,6/0,8 |
|
|
разница между значениями к. п. д. |
||||
0,7 |
|
|
в точке его максимума и при угле |
|||||
1,0 |
1,2 1,<+ и, выхода потока из турбины a2 = 90a |
|||||||
Рис. 1.10. Зависимость т)и |
от щ и р |
не очень велика. Поэтому, |
исходя |
|||||
из эффективной |
работы |
затурбин- |
||||||
|
|
|
|
ного устройства, а также по проч ностным соображениям иногда можно рекомендовать не стре миться выдерживать оптимальные соотношения по й\ и р.
Проведенный анализ показывает, что, определив кюпт и pomno формулам (1.23) и (1.21), в расчете центростремительной тур бины можно применять значения
h рас ~ |
м 1 опт |
0,1 |
(1.26) |
|
Ррас ~ |
Ропт |
' 0> 1 |
||
|
Выбирать значения й\ и р большими оптимальных не реко мендуется в основном по прочностным соображениям. Однако в некоторых случаях турбины работают на правой ветви зависи
мости к. п. д. от |
M I (например, когда на приводимом турбиной |
агрегате нельзя |
уменьшать число оборотов, а выбирать малый |
диаметр турбины нельзя по конструктивным соображениям). Из
рис. |
1.10 видно, что зависимость к. п. д. от й\ обрывается |
на пра |
||
вой |
ветви при некотором |
значении Гц тем меньшем, |
чем мень |
|
ше р. При малых р ( р < 0 , 4 ) |
зависимость ч\ = ц(й\) обрывается на |
|||
левой ветви, не достигая максимального значения. |
|
|
||
В этих точках величина относительной скорости |
на |
выходе |
из колеса Й72 =0. Режимы, при которых г£2 = 0, характерны лишь для центростремительной турбины. Условие й>2 = 0 накладывает
20

определенные ограничения на выбор параметров ступени. Оста новимся на этом вопросе подробнее.
Назовем |
степень реактивности, при которой W2 = 0, |
минималь |
||||||||
ной и определим ее. |
|
|
|
|
|
|
||||
Несложные преобразования приводят к соотношению |
||||||||||
|
|
- ^ c o s a ^ ± |
V |
4<f>2 coS 2 а х |
г4 — 4 (<р2 _ 1) ( l f-Щ) |
1.27) |
||||
Pmln = 1 |
|
|
|
|
|
_ _ |
_ _ |
|||
Из |
выражения |
(1.27) |
видно, |
|
|
|||||
что величина pm m зависит от й\, ц, |
|
|
||||||||
а\ и ф. На рис. 1.11 |
приведена эта |
|
|
|||||||
зависимость, |
рассчитанная |
по |
|
|
||||||
формуле |
(1-27). |
Зависимость |
|
|
||||||
(1.27) |
особенно упрощается |
при |
|
|
||||||
Ф = 1 . |
В |
этом |
случае |
после |
рас |
|
|
|||
крытия |
неопределенности |
|
|
|
|
|||||
( P m l n = 1 — |
1 — \y?U2 . 2 |
(1.28) |
|
|
||||||
^2 COS ОЦИ! |
|
|
||||||||
|
|
|
|
|
|
|
|
|||
В случае |
осевой турбины |
(|х = |
|
|
||||||
= 1) принципиально также возни |
|
|
||||||||
кают |
ограничения |
по |
выбору р, |
|
|
|||||
но диапазон |
возможных значений |
|
|
рнастолько широк (например,
при Ц) = 0,4 и ai = 15° р т т = 0 , 8 ) ,
что нет смысла говорить о мини мальной реактивности в осевой турбине и, тем более, — в центро бежной.
В заключение отметим, что вопросам выбора оптимальных параметров ступени уделено вни мание во многих работах [4, 8, 17, 18, 27, 31, 36, 39, 41]. Большинство авторов при расчете оптимальных параметров исходят из зависимо сти к. п. д. от двух переменных, считая остальные параметры пос тоянными. В частности, в работе [18] получены примерные соотно шения между оптимальными зна
чениями р И H i При ф = 1|)=1 и а г = 90°.
Формула связи р и й\ по данным работы [18] следующая:
и. |
1 |
(1.29) |
|
||
2 COS _! ^ 1 — р |
|
|
|
|
21
тах [17] и [31]. |
|
|
|
|
|
|
- |
|
|||
На рис. 1.12 приведено сравнение |
оптимальных |
р и и 1, полу- |
|||||||||
ченных по формулам (1.21) и (1.23) и по данным рга б о т ы |
[18]. |
||||||||||
Как и |
следовало |
ожидать, |
|
|
|
|
|
||||
по работе |
[18], поскольку |
в ней |
|
|
|
|
|
||||
принято |
|
аг = 90°, |
при |
данной |
|
|
|
|
|
||
степени |
реактивности |
р |
опти |
|
|
|
|
|
|||
мальные |
й~1 ниже, чем по фор |
0,7 |
|
|
|
|
|||||
муле (1.21). Кроме того, отме |
|
|
|
|
|
||||||
тим, что в этой работе |
дается |
|
|
|
|
|
|||||
только связь |
р и |
M I , а не сами |
0,6 |
|
|
|
|
||||
их значения, в то время как по |
|
|
|
|
|
||||||
изложенной |
выше |
методике |
|
|
ч |
|
|||||
можно вычислить |
оптимальные |
0,5 |
|
|
|||||||
величины |
|
р и нь |
причем для |
|
|
||||||
расчета |
использована |
зависи |
|
|
|
|
|
||||
мость т)и |
|
от трех |
переменных. |
°>t |
|
|
|
|
|||
Так же, как в работе [18], из |
9 |
|
|
|
|
||||||
зависимости |
к. п. д. от |
двух |
|
|
|
|
|||||
переменных исходили и в рабо |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
0,2 |
|
|
\ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
20 |
30 |
k0 |
рг |
Рис. 1.12. Сравнение оптимальных значе
ний щ и р при ccj =20°:
/ — по формулам (1.23) и (1.21); 2 — по данным работы [18]
Рис. 1.13. Сравнение оптимальных значений и\ и р при ц = 0,5; <р = 0,95; •ф = 0,85:
а — по формулам |
(1.23) и |
(1.21); б — по |
данным |
работы |
[17]; |
——— |
1, =15°; |
а, =20°; а, =30'
Так, в работе [17] из условия минимума пот'ерь выведены сле дующие зависимости для оптимальных р и и\1с\\
ЛЬ, |
f< Р(1 |
2) ф2 cos2p2 |
|||||
v- V |
(—1 -ф |
|
cos |
|
рг) |
||
|
2 |
2 |
|
2 |
|
2 |
|
р = 1 - |
|
|
|
|
1 |
|
|
2? |
|
cos ctj + |
|
|
Хотя в отличие от формул в работе [18] эти формулы позволя ют независимо одна от другой определить «i/ci и р, в работе [17] использована зависимость к. п. д. только от двух переменных ве-
22
личин и получены заниженные значения щ и р по сравнению с определенными по формулам (1.23) и (1.21).
Это сравнение величин и\ и р приведено на рис. 1.13. Согласно работе [36] для расчета оптимальных параметров
ступени приходится решать сложные уравнения шестой степени, что очень неудобно в практике расчета турбины.
Зависимость к. п. д. от трех переменных использует для ана лиза оптимальных условий работы центростремительной турби ны М. С. Приходько [39]. Однако, как уже обмечалось, Б работе
dc
[39] ошибочно привлечено лишнее условие — = 0 , в связи с чем, dm
по сути дела, проведен анализ оптимальных условий работы сту пени турбины для двух независимых переменных.
В рассмотренных выше работах оптимальные параметры сту пени выбирались исходя из одномерной мод'ели течения. Инте ресна попытка, предпринятая в работе [31], дополнить общепри нятую методику расчета ступени по одномерной теории учетом
параметров потока по радиусу |
Е Ы Х О Д Н О Г О |
сечения. Однако, |
как |
|
мы видели, принятый при этом в качестве основного условия |
осе |
|||
вой выход потока из ступени |
не является |
обязательным. |
При |
|
определении оптимальных параметров ступени не |
использова |
|||
лось уравнение неразрывности |
и, следовательно, не |
накладыва |
лось условие на производительность турбины. В практике расче та турбин расход рабочего тела, как правило, достаточно точно определен. Рассматриваемая методика может быть использова на, когда расход рабочего тела через турбину задан. Оптималь ные параметры ступени т, й\ и р зависят от ai, Рг и ц. При про извольном варьировании значениями ai, Р2 и \х могут получиться явно неконструктивные соотношения размеров проточной части. Поэтому выбор параметров и основных размеров ступени по предложенной методике должен проводиться в два этапа: сна чала определяются оптимальные значения И\ и р при произволь но выбранных ai, Р2 и д., затем рассчитываются основные разме ры ступени при заданном расходе рабочего тела. Если при этом получаются неудовлетворительные соотношения размеров про точной части, то определяются такие значения ai, Р2 и ц, при ко торых проточная часть конструктивно и технологически выпол нима и обеспечивает высокую экономичность. Для определенных таким образом си, Рг и д. заново рассчитываются оптимальные значения т, и\ и р.
Для облегчения расчетов на первом этапе приведем наиболее распространенные в практике применения радиальных турбин в турбокомпрессорах и турбодетандерах значения некоторых пара метров ступеней:
•—угол |
выхода |
потока |
из соплового |
аппарата ai = |
124-30; |
— угол |
выхода |
потока |
из рабочего |
колеса Рг = 20 |
-М5; |
— степень радиальности ц = 0,3-^0,5;
23
— относительный диаметр втулки колеса dBT/di |
= 0,15ч-0,4; |
||||||
— относительная |
длина |
сопловой |
лопатки |
A/di = 0,02-f-0,14; |
|||
— относительная |
длина |
рабочей лопатки |
на |
выходе У,й\ = |
|||
= 0,1ч-0,3; |
|
|
|
|
|
|
|
— коэффициент |
скорости |
в сопловом |
аппарате |
<р = 0,95-^0,97; |
|||
— коэффициент |
скорости |
в рабочем |
колесе |
ф = 0,8ч-0,9. |
|||
При выбранных |
параметрах ступени |
оптимальные значения |
|||||
Ропт, «опт и йюпт располагаются в диапазонах: |
|
|
|
||||
« о п т = 1,15-^-1,25; ропт = 0,2-^0,5; H i o n |
T = 0,6-f-0,75. |
||||||
Если параметры |
ступени |
выбраны |
близкими |
к |
оптимальным, |
а профилирование отдельных элементов проточной части выпол нено так, что обеспечивается безотрывность потока на расчетном режиме, то можно ожидать получить значение к. п. д. ци в пре делах 0,90—0,93, а значение мощностного к. п. д. г]т = 0,86-^-0,9.
Глава II
СТ А Ц И О Н А Р Н О Е Т Е Ч Е Н И Е В С О П Л О В Ы Х АППАРАТАХ
ИРАБОЧИХ КОЛЕСАХ
Впредыдущей главе были определены основные параметры ступени. При этом предполагалось, что потери в сопловом аппа рате и в рабочем колесе известны и могут быть получены путем
соответствующего профилирования решеток соплового аппарата и рабочего колеса.
Задачей настоящей главы является изучение течения газа че рез сопловые и рабочие решетки с тем, чтобы на этой основе:
а) суметь рассчитать действительные параметры потока при
известной (выбранной) форме и размерах |
проточной части; |
||
б) определить формы и размеры, обеспечивающие наиболь |
|||
шую экономичность |
процесса. |
|
|
Расчет течения |
в |
сопловом аппарате |
центростремительной |
турбины не является |
существенно новой задачей по сравнению |
с расчетом течения в сопловом аппарате осевой турбины, но те чение в колесе центростремительной турбины значительно слож нее. Обычно при расчете осевых турбин используют данные про дувки плоских турбинных решеток, т. е. считают, что поверхности тока не отличаются «ли мало отличаются от поверхности соосных круговых цилиндров. Теория расчета таких течений доста точно развита и накоплен громадный опыт проектирования осе вых турбин. В частности, созданы методики расчета, основанные на большом количестве экспериментальных данных продувки плоских решеток.
Потери энергии в решетках осевой и радиальной турбин скла дываются из следующих видов потерь:
а) потери трения на профиле; б) кромочные потери; в) вторичные потери;
г) волновые потери, связанные с образованием зон сверхзву ковых скоростей и появлением скачков уплотнения.
Сумму потерь трения на профиле и кромочных обычно на зывают профильными потерями.
Кратко рассмотрим, следуя Г. Ю. Степанову и В. Л. Эпштейну [47], как рассчитываются профильные и вторичные потери.
25
Потери трения на профиле при безотрывном обтекании не сжимаемой жидкостью выражаются известной формулой
*t sin (32
где 02 — потери импульса соответственно на выпуклой (сп) и вогнутой (вог) сторонах профиля в выходном сечении решетки.
Величина потери импульса определяется следующим выра жением:
/ |
7 |
_ |
_ \ 0 ' 8 |
8**=0,0361 Re - ° - 2 w - 3 - 2 8 |
j |
wz>™ds |
В выходном сечении решетки будем иметь
/ Т _ |
|
_\°>8 |
|
|
8r=0,0361Re-0 -2 М w |
3 - 8 |
6 ^ |
, |
(2.2) |
где й> и s—-относительные величины (скорости отнесены к ско рости w2, а длины отнесены к хорде профиля).
Если в формулу (2.2) подставить величины средних скоростей на профиле по вогнутой и выпуклой сторонам, то формула (2.1) примет следующий вид:
|
j. |
_ |
0,072 |
/ |
s \о,8 |
3,09 |
|
/ да„ |
(2.3) |
|
|
|
~ R e ° > 2 sin Э 2 к |
I |
t |
w2 I |
|
\ |
w2 |
||
|
T P |
|
|
|||||||
где wca и |
|
|
— средние скорости на профиле. |
обычно |
пренебре |
|||||
Формула |
|
(2.3), где второй член в скобках |
||||||||
жимо мал по сравнению с первым, может |
быть непосредственно |
|||||||||
применена для расчета |
потерь |
на трение |
в |
решетке |
соплового |
|||||
аппарата |
радиальной турбины, |
если будет |
определена средняя |
|||||||
скорость |
на |
|
профиле. |
|
|
|
|
|
|
|
Применительно к расчету потерь на трение на профиле ра бочего колеса радиальной турбины формула (2.3) требует не которых оговорок. Дело в том, что в пограничном слое рабочего колеса действуют центробежные и кориолисовы силы инерции, не учтенные формулой (2.2). Однако практика расчета и экспе римента показала, что учет этих сил приводит к незначительно му изменению величины 6**. Так, К- М. Каминская [1], сравни вая полученные расчетным путем величины потерь трения на профиле неподвижной и вращающейся пластин, показала, что напряжение трения в ламинарном пограничном слое в каждом сечении вращающейся пластины очень близко к напряжению трения на продольно обтекаемой неподвижной пластине.
Экспериментальное исследование гидравлического сопротив ления вращающихся труб при турбулентном режиме течения
26
[30] показало, что их сопротивление приблизительно равно со противлению неподвижных труб. Ниже будут приведены еще не которые экспериментальные обоснования, позволяющие считать возможным использование теории плоского пограничного слоя.
Это позволяет приближенно принимать выражение (2.2) и при расчете потерь трения вращающейся решетки, а отсюда и фор мулу (2.3) для расчета £Тр-
Отметим также, что формула (2.3), выведенная для несжи маемой жидкости, может быть применена при вычислении сред них скоростей на профиле и для сжимаемой жидкости [47].
Зная потери на трение, можно вычислить профильные по тери, поскольку зависимости для коэффициентов кромочных по терь известны [47]. Вторичные потери вычисляют в зависимости от профильных потерь в решетке.
Таким образом видно, что в расчете потерь энергии основу составляют потери трения, поскольку все остальные виды потерь можно затем рассчитать в зависимости от потерь трения. Для расчета этих потерь необходимо определить средние скорости на профиле решетки радиальной турбины.
2.1. РАСЧЕТ СРЕДНИХ СКОРОСТЕЙ
Рассмотрим двухмерное неплоское движение идеальной сжи маемой жидкости через вращающуюся решетку турбины с про извольной формой меридианного профиля проточной части. Для определения средних скоростгй на профиле применим уравнения
неразрывности, моментов количест ва движения, энергии, состояния и условие равенства нулю циркуля ции абсолютной скорости [32].
Для контура abcdefsgha |
(рис. 2.1) |
|
|
|||
условие равенства нулю |
циркуляции |
|
|
|||
абсолютной |
скорости имеет следую- » |
|
|
|||
щий вид: |
|
|
|
|
|
|
Рис. 2.1. К |
расчету средних |
ско |
|
|
||
|
ростей |
|
|
|
|
|
|
{Wl |
COS |
- j - |
щ) tx - f - (w2 COS (32 |
— U2) t2 |
= |
|
|
scn |
|
sbot |
|
|
|
= |
f (®cn + |
« , ) r f s - j (wB0T + |
us)ds, |
(2.4) |
|
|
|
b |
|
о |
|
|
где us — проекции окружной скорости на направление дуги про филя s, а параметры газа перед и за решеткой (на линиях а — b и с — d) принимаются равномерными (средними).
27
Уравнение моментов количества движения для того же кон тура
[(те»! cos р! + щ) гх + (w2 cos р2 — и2) r2 ] G =
|
|
^вог |
|
|
|
scn |
|
|
|
|
||
|
= |
J |
pB0Tlr |
cos aB0Tds— |
|
§ pcJr |
cos acnds, |
(2.5) |
||||
|
|
о |
|
|
|
|
0 |
|
|
|
|
|
где |
G = p2uy2 /2 |
4 — расход |
газа. |
|
|
|
|
|
|
|||
|
Уравнение энергии |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
2 |
_ * |
|
2 |
|
|
|
|
|
|
|
|
|
И, |
|
Ыч |
|
|
|
|
|
|
|
S 7 ^ . |
|
|
£ - = С / « . |
|
^ - |
|
( 2 - 6 ) |
||
|
Циркуляция Г окружной скорости по профилю равна потоку |
|||||||||||
вихря окружной скорости через площадь сечения профиля |
|
|||||||||||
|
|
J u/ls = |
^ |
Rot udf |
= |
j^ |
rQ |
sin |
ydf, |
|
||
где |
й — угловая |
скорость |
вращения; |
|
|
|
|
|||||
|
Y — угол |
наклона |
осесимметричной |
поверхности тока |
к оси. |
|||||||
|
Окончательно |
|
|
|
|
|
|
|
|
|
||
|
|
|
usds |
2Ы/ |
sm Ycp ~ |
|
• |
|
|
|
||
|
Тогда уравнение (2.4) |
можно переписать так: |
|
|||||||||
|
|
(w1 |
cos ^ + |
И]) ^ + (w2 |
cos р2 — |
и2 ) г?2 = |
|
|||||
|
|
|
j |
|
|
* _ |
|
|
|
|
|
|
|
|
|
, |
, |
|
, |
, |
i S M j / slir rl |
i T c p |
|
||
|
|
|
|
|
|
Г |
, |
i |
^ " 1 / |
|
T c p |
|
|
|
|
|
?— |
j |
|
таВ0^Н |
|
• |
(2-7) |
||
|
|
|
о |
|
|
0 |
|
|
|
|
|
|
Введем средние величины
|
s |
|
|
w=— |
^wds; |
h=— |
^hds; rj = — J r r f s . |
|
о |
б |
о |
Из уравнения |
(2.6) |
получим |
|
где
''г |
1 |
«2 |
/"1
28
Перепишем уравнение (2.7). в безразмерном виде:
|
|
|
|
2tr- |
|
г |
I . |
|
|
|
|
|
|
1 |
|||
|
|
|
|
|
'ср |
|
|
X |
|
|
|
|
1 |
|
|
V |
|
|
|
|
|
|
1 |
тср |
||
X \ Y |
|
|
|
|
|
|||
г |
|
|
|
|
|
|
|
|
|
|
COS й -2-+ |
X2 cos32 + XUa |
( - 1 - 1 |
||||
|
|
2 Х Ц | / |
sin |
7 с р |
|
(2.8) |
||
|
|
|
|
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
'ср = |
|
|
|
|
|
Тер |
Т«>! |
х - |
|
1 |
Х 2 |
/ 1 |
|
|
2 ( х + |
1) |
|
|
|
||||
|
|
|
|
|
|
|||
Обозн ачая s/^cp |
— s и полагая |
5 с п ~ £ в о г ~ 5 , получим оконча |
||||||
тельно |
|
|
|
|
|
|
|
|
/ |
Тп |
^-сп — Кот2 |
|
|
ср |
• Х2 cos8 -4-Х ( 1 |
И 2 X " - / S I N Т С Р |
где
f
•+
(2.9)
r\t.ср
Относительная площадь профиля / в рабочих колесах ра диальной турбины мала, так как профили лопаток этих колес обычно имеют малую толщину. Поэтому в дальнейших выклад-
2Х |
Ц |
/ sin TfCp |
ках член |
|
—=—• опустим. |
|
|
S |
Выразим давления р, входящие в правую часть уравнения (2.5), формулой
и разложим ЩХ) в ряд по степеням X относительно среднего значения Я
п(X) = пfx)+±=± гг(X) + ( Х - 1 ) 2 п - ( Х ) + . . . .
29