книги из ГПНТБ / Жуковский М.И. Расчет обтекания решеток профилей турбомашин
.pdfПродолжение табл. 11
|
Решетка № 7; |
q = 0,702; |
|
|
Решетка № 8; q = 0,514 |
|
||
|
7= 0,644 |
|
|
7= 1,111 |
|
|||
в° |
X |
У |
S |
Ь° |
X |
У |
S |
|
т~ |
"Г |
|||||||
|
|
|
|
|||||
0 |
1,478 |
0,197 |
0,032 |
0 |
0,411 |
0,312 |
0,125 |
|
10 |
1,435 |
0,255 |
0,053 |
15 |
0,374 |
0,412 |
0,183 |
|
20 |
1,383 |
0,321 |
0,078 |
30 |
0,334 |
0,516 |
0,244 |
|
30 |
1,321 |
0,393 |
0,106 |
45 |
0,292 |
0,618 |
0,304 |
|
40 |
1,253 |
0,469 |
0,136 |
60 |
0,250 |
0,710 |
0,360 |
|
50 |
1.176 |
0,545 |
0,167 |
75 |
0,206 |
0,790 |
0,409 |
|
60 |
1,092 |
0,621 |
0,201 |
90 |
0,167 |
0,856 |
0,451 |
|
70 |
0,998 |
0,693 |
0,236 |
105 |
0,136 |
0,900 |
0,480 |
|
80 |
0,893 |
0,759 |
0,272 |
120 |
0,112 |
0,921 |
0,498 |
|
90 |
0,779 |
0,815 |
0,309 |
135 |
0,100 |
0,916 |
0,508 |
|
100 |
0.663 |
0,852 |
0,345 |
150 |
0,106 |
0,888 |
0,525 |
|
ПО |
0.548 |
0,869 |
0,379 |
165 |
0,129 |
0,833 |
0,558 |
|
120 |
0,439 |
0,869 |
0,411 |
180 |
0,167 |
0,753 |
0,606 |
|
130 |
0,344 |
0,856 |
0,439 |
195 |
0,213 |
0,658 |
0,664 |
|
140 |
0,262 |
0,833 |
0,465 |
210 |
0,265 |
0,553 |
0,728 |
|
150 |
0,192 |
0,806 |
0,487 |
225 |
0,318 |
0,445 |
0,794 |
|
160 |
0,140 |
0,780 |
0,504 |
240 |
0,367 |
0,341 |
0,851 |
|
170 |
0,109 |
0,754 |
0,515 |
255 |
0,410 |
0,248 |
0,912 |
|
180 |
0,100 |
0,724 |
0,525 |
270 |
0,442 |
0,174 |
0,957 |
|
190 |
0,106 |
0,693 |
0,535 |
285 |
0,466 |
0,122 |
0,987 |
|
200 |
0,125 |
0,669 |
0,544 |
300 |
0,476 |
0,100 |
0 |
|
210 |
0,168 |
0,659 |
0,558 |
315 |
0,473 |
0,112 |
0,009 |
|
220 |
0,246 |
0,661 |
' 0,581 |
330 |
0,461 |
0,155 |
0,034 |
|
230 |
0,367 |
0,659 |
0,616 |
345 |
0,440 |
0,224 |
0,074 |
|
240 |
0,523 |
0,638 |
0,662 |
|
|
|
|
|
250 |
0,700 |
0,591 |
0,716 |
|
|
|
|
|
260 |
0,882 |
0,521 |
0,773 |
|
|
|
|
|
270 |
1,056 |
0,430 |
0,831 |
|
|
|
|
|
280 |
1,207 |
0,333 |
0,883 |
|
|
|
|
|
290 |
1,327 |
0,244 |
0,927 |
|
|
|
|
|
300 |
1,415 |
0,169 |
0,961 |
|
|
|
|
|
310 |
1,469 |
0,119 |
0,983 |
|
|
|
|
|
320 |
1,499 |
0,100 |
0,994 |
|
|
|
|
|
330 |
1,518 |
0,107 |
0 |
|
|
|
|
|
340 |
1,525 |
0,124 |
0,006 |
|
|
|
|
|
350 |
1,510 |
0,152 |
0,016 |
|
|
|
|
112
Продолжение табл. 11
|
Решетка № 9; |
q = 0,725; |
|
Решетка № 10; |
q = 0,462; |
|
|
|
t = 0,619 |
|
|
t = 1,226 |
|
||
|
|
|
S |
|
|
|
S |
6° |
X |
У |
т |
0° |
X |
У |
|
|
|
||||||
|
|
|
|
||||
0 |
1,418 |
1,047 |
0,954 |
» |
0,417 |
0,674 |
0,900 |
15 |
1,423 |
1,052 |
0,965 |
15 |
0,457 |
0,743 |
0,948 |
30 |
1,364 |
1,052 |
0,983 |
30 |
0,482 |
0,791 |
0,981 |
45 |
1,253 |
1,042 |
0 |
45 |
0,492 |
0,816 |
0,996 |
60 |
1,091 |
1,014 |
0,033 |
60 |
0,486 |
0,813 |
0,007 |
75 |
0,901 |
0,970 |
0,092 |
75 |
0,461 |
0,785 |
0,029 |
90 |
0,714 |
0,895 |
0,152 |
90 |
0,421 |
0,738 |
0,066 |
105 |
0,549 |
0,802 |
0,200 |
• 105 |
0,370 |
0,678 |
0,113 |
120 |
0,419 |
0,704 |
0,240 |
120 |
0,316 |
0,608 |
0,167 |
135 |
0,314 |
0,601 |
0,284 |
135 |
0,263 |
0,531 |
0,223 |
150 |
0,231 |
0,495 |
0,325 |
150 |
0,215 |
0,449 |
0,281 |
165 |
0,167 |
0,386 |
0,364 |
165 |
0,173 |
0,364 |
0,337 |
180 |
0,122 |
0,284 |
0,398 |
180 |
0,139 |
0,283 |
0,390 |
195 |
0,100 |
0,198 |
0,425 |
195 |
0,115 |
0,211 |
0,435 |
210 |
0,103 |
0,132 |
0,445 |
210 |
0,103 |
0,154 |
0,471 |
225 |
0,129 |
0,100 |
0,458 |
225 |
0,100 |
0,115 |
0,494 |
240 |
0,177 |
0,118 |
0,476 |
240 |
0,105 |
0,100 |
0,506 |
255 |
0,257 |
0,204 |
0,512 |
255 |
0,120 |
0,113 |
0,519 |
270 |
0,381 |
0,360 |
0,573 |
270 |
0,145 |
0,153 |
0,547 |
285 |
0,563 |
0,559 |
0,653 |
285 |
0,179 |
0,215 |
0,589 |
300 |
0,784 |
0,744 |
0,740 |
300 |
0,218 |
0,296 |
0,644 |
315 |
1,004 |
0,888 |
0,819 |
315 |
0,265 |
0,391 |
0,707 |
330 |
1,201 |
0,982 |
0,885 |
330 |
0,317 |
0,491 |
0,775 |
345 |
1,345 |
1,028 |
0,930 |
345 |
0,369 |
0,587 |
0,841 |
8 М И. Жуковский 700 |
113 |
и |
Приведем два примера расчета решетки |
при |
вариации |
шага |
|||
угла |
установки. |
|
|
|
|
|
|
|
Изменение шага решетки. Пусть необходимо рассчитать обте |
||||||
кание |
решетки направляющих |
лопаток с |
относительным |
шагом |
|||
t |
= 0,72, если известен расчет этой решетки при |
~ 0,56. |
Угол |
||||
установки профилей |
— 46°35'. |
Профили при t — 0,56 и 7 = 0,72 |
|
Фиг. |
47. |
Расчет |
близкой |
решетки при |
вариации |
шага: |
|
|
-------------- |
|
профиль |
близкой решетки;---------------------- |
профиль |
задан |
||
|
ной решетки; —о—о—о — профиль решетки, рассчитанной со смеще |
|||||||
|
|
|
|
нием 6 (1-е приближение). |
|
|
||
показаны |
на |
фиг. 47, а. |
Используя график функции &о = /( о»%) |
|||||
(фиг. |
22), |
получим % = 333° для рассчитанной решетки, имеющей |
||||||
шаг 1 |
7о = 0,56 (густота |
b — 1,78) и 90 = 329° для рассчитываемой |
||||||
решетки с шагом 71 = 0,72 |
(Ь |
1,39). Поворот |
круга равен, следо |
|||||
вательно, дб = 4°. |
|
для шага t2 = 0,56, осуществляем смеще |
||||||
Построив график х (9) |
ние по 6, как это показано на фиг. 47, б. Снятые с графика смещенные
1 Для приближенных расчетов смещения по 6 можно принять величину Ьп равной хорде профиля Ь, а угол р0 равным углу установки профиля.
114
на дб = 4° значения х*(0) |
умножением на , после чего снимаем |
_ 62 |
_ _ |
спрофиля величины у*. Приближенные значения х*, у* для решетки
сti = 0,72 принимаются в качестве нулевого приближения. Первое
Фиг. 48. Расчет близкой решетки при вариации угла установки
-------------------профиль близкой решетки;----------------------- |
|
про |
|
филь заданной решетки; —о—о—о------ |
профиль решетки, |
рас |
|
считанной со смещением 6 (1-е приближение). |
|
|
|
приближение, вычисленное по методу, изложенному |
в |
п. 14, нане |
|
сено на эту же фиг. 47, а. Дальнейший ход расчета |
осуществляется |
путем переноса точек на заданный профиль при помощи формулы (28) и вычисления последующих приближений.
Изменение угла установки. Необходимо рассчитать обтекание
решетки направляющих лопаток, повернутых в сторону уменьшения выходного угла на — 5° (фиг. 48, а). Обтекание исходной решетки
известно. Шаг решетки t = 0,59. Используя график фиг. 22, имеем
8* |
115 |
t = 0,59; ~ = 1,70; |
= 44°06'; |
% |
= 333°48' и |
7 - 0,59; |
|
-y-= 1,70; |
= 49°06'; |
90 = 330°08'. |
Смещение no 6 равно, |
||
3°. На фиг. |
48,6 нанесено |
значение |
x (6) |
для исходного |
профиля. |
При помощи снятых с этого графика величин х (9) с учетом смещения
Д9 = 3° определяются у (9). Указанные значения х и у пересчиты ваются для профиля, повернутого по известным формулам поворота координат
х* — х cos Дрй + у sin дрй;
у* = у cos дрй — xsin дрй,
после чего по смещенным значениям х* с контура повернутого про филя снимаются величины у*. Эти величины х* и у* используются
в качестве нулевого приближения. Результаты вычисления х(1), z/<» первого приближения показаны на фиг. 48, а. Полученный про филь уже мало отличается от заданного.
16. КРУГОВЫЕ ПЛОСКИЕ РЕШЕТКИ
Задача об определении потока в круговой решетке рассматрива лась Шпанхаке [42], П. А. Вальтером, П. В. Мелентьевым [33],
Д. А. Войташевским [39, 43], Г. И. Костычевым [44], Л. А. Дор фманом [45], Г. П. Таушкановым и др.
Профили в плоской круговой решетке расположены симметрично
относительно центра со сдвигом, равным где п — число профилей
(фиг. 2). Течение в таких решетках схематизирует движение в направ ляющих и рабочих лопатках радиальных турбомашин, а также обте кание лопаток в коническом слое постоянной толщины.
Задача об определении обтекания неподвижной круговой решетки может быть сведена к расчету прямой решетки. Известны также спо собы [44], позволяющие решать задачу без перехода к прямой решетке.
Конформное отображение круговой решетки на прямую выпол няется следующим образом.
Пусть прямая решетка профилей расположена в плоскости z. Функция, осуществляющая конформное отображение плоскости z
на плоскость С круговой решетки с центром в начале координат,
представляется |
в |
виде |
|
|
|
|
|
С = Aekz, |
(37) |
где |
А — |
— постоянная. |
|
|
|
Перемещению на период — в плоскости С соответствует смещение |
|||
на |
it в плоскости |
z. Следовательно, |
|
116
Таким образом, равенство (37) принимает вид
С = |
2к |
(38) |
Разрешая (38) относительно z, получаем следующую формулу, выполняющую отображение круговой решетки на прямую
Расчет обтекания прямой решетки в плоскости z выполняется затем по методу, изложенному в п. 14. В случае редкой круговой решетки в качестве вспомогательной области может быть взята
также внешность или внутренность круга с симметрично располо женными вихреисточниками (фиг. 10). Решетчатая область исполь зуется в расчетах круговых решеток в работах П. В. Мелентьева [33]
иЛ. А. Дорфмана [45].
Вслучае вращающейся круговой решетки течение в относительном потоке является вихревым. Граничным условием при решении задачи в абсолютном движении является равенство между нормаль ной составляющей скорости vn абсолютного движения и нормальной составляющей скорости движения профиля ип
vn(s) = un(s). |
(40) |
При этом должно выполняться условие
[ v„ (s) ds = 0,
L
где интегрирование ведется по контуру профиля.
Задача об обтекании вращающейся круговой решетки решается по следующей общепринятой схеме [42]. Сначала рассчитывается обтекание круговой решетки, создаваемое вихреисточником и вихрем.
Затем определяется так называемый поток вытеснения, учитывающий
вращение круговой решетки в неподвижной жидкости.
Применим изложенный в пп. 11 и 14 метод для решения данной
задачи в случае круговых решеток большой густоты. Выполнив преобразование (39), получим в плоскости z прямую решетку профи лей, которую затем отобразим на решетку кругов (п. 14). Таким обра
зом определяется соответствие z (С) между точками профиля и окруж ности.
Вследствие потенциальности абсолютного движения поток вытес
нения может быть определен следующим образом. На основании оче видной формулы
(wx — iwv)e^ = ws — iwn,
( — угол касательной к профилю, ws и wn — проекции скорости
соответственно на касательную и нормаль к контуру профиля)
117
аналитическая функция vs — ivn |
может быть представлена в виде |
|
ряда |
|
|
п=0 |
|
|
ИЛИ |
|
|
vs -ivn = C+ % C_S~n + 2 СЛп . |
(42) |
|
n=l |
л=0 |
|
где |
|
|
С = С + iC".
Первые два коэффициента С и Ci разложения (41) или (42) опреде
ляются из граничных условий задачи. Так как при С — ±°° ско'
рость потока вытеснения равна нулю, то
С = C_j = 0.
Тогда разложение (42) может быть представлено в виде
vs~ivn= |
|
n=0 |
(44) |
|
|
n=2 |
|
||
ИЛИ |
|
|
|
|
vs = У I A_n |
cos nfi + 3] (B_„ — B„] sin«9 4- |
|
||
n=2 |
|
n=2 |
|
|
|
4- A i cos 9 — Bi sin 9 4~ An, |
(45) |
||
v„ = У (B-n 4- Bn) cos «9— У (A_n — An) sin n9 4~ |
|
|||
n=2 |
|
n=2 |
|
|
|
4- Bi cos f) |
4- 4j sin 9 4- Bo, |
(46) |
|
где |
|
|
|
|
Cffl = Am 4- iBm |
(m = 0, |
x 1, r2,. ..) |
|
Определение коэффициентов разложения (45) производится сле
дующим образом. Величины vn (s) определяются на окружности в решетке кругов по найденному из расчета прямой решетки отобра жению z (С). Функция vn (9, q) разлагается затем в ряд, представ ляющий собой*мнимую часть (44). В результате определяются коэф фициенты
Рп = -4—п + Ап и Qn = В_„ — В„.
Далее составляются и решаются системы уравнений, аналогичные уравнениям (65)—(68). гл. II. После определения всех коэффициентов
А—п, В_п величины А„ и В„ находятся по формулам, аналогичным (23),
118
и при помощи ряда (45) вычисляются проекции скорости потока вытеснения на касательную к контуру профиля. Проекции абсолют ной скорости на касательную к профилю результирующего движения
находятся |
затем путем суммирования скорости потока вытеснения |
и скорости |
поступательно-циркуляционного течения. |
17.ПРИБЛИЖЕННЫЕ МЕТОДЫ РАСЧЕТА ОБТЕКАНИЯ ГУСТЫХ ТУРБИННЫХ РЕШЕТОК ПРОФИЛЕЙ
Приближенный расчет распределения скоростей по профилю в густой решетке может быть выполнен путем использования решения задачи о течении в канале. В случае густых решеток обтекание выход
ной |
части |
профиля |
практически не зависит от условий течения |
на |
входном |
участке |
профиля. |
Сравнительно слабым в этом случае является влияние сосед них профилей на течение
всредней части канала. Поэтому
вприближенной постановке те чение в средней части межпро фильного канала можно рас сматривать как течение в еди ничном канале, пренебрегая взаимным влиянием профилей.
Известно несколько спосо бов решения задачи об опреде
лении потока в канале [32, |
|
46, |
47]. |
Дальнейшее развитие метода |
|
||
расчета потока в |
канале дано Фиг. 49. К расчету канала R2 > Ri. |
||
в работах Г. Ю. Степанова [48], |
O1O2 — Е • |
||
Ю. И. Швеца [49], |
Г. С. Самой |
||
|
ловича [50], А. Н. Шерстюка [51], М. Е. Дейча, М. И. Жуковского и др.
В работе А. Н. Шерстюка 151 ] предложены приближенные формулы для расчета течения в каналах, удовлетворяющих усло вию 7?1 < 7?о —/ [7?i и 7?о —радиусы кривизны соответственно выпуклой и вогнутой стенок канала, I — ширина канала (фиг. 49)].
Автором разработан простой метод расчета потока несжимаемой
жидкости в каналах с границами, описанными дугами окружностей 1 с произвольным эксцентриситетом. Этот метод, основанный на исполь
зовании точного решения |
для течений, создаваемых источником |
|
и |
стоком, а также парой |
вихрей, позволяет рассчитывать течение |
в |
каналах как при 7?i < |
так и ПРИ |
Рассмотрим возможные типы каналов, которые образуются между
линиями тока для течений, создаваемых источником и стоком и парой
вихрей.
1 Турбинные профили, обводы которых описаны дугами окружностей, широко применяются в турбостроении.
119
Источник и сток. Поместим на оси |
ох в точках х = +т источ |
||||
ник и сток (фиг. 50). |
|
|
|
|
|
Характеристическая функция этого течения имеет вид |
|
||||
|
|
J_ln |
|
(47) |
|
|
|
2л |
т |
— z |
|
|
|
|
|||
После разделения выражения (47) на вещественную и мнимую |
|||||
части, получаем: |
|
|
|
|
|
|
Q |
1п |
У(х-гтУ + у\ |
|
|
|
2я |
|
X(х — т)2 4- у2 ’ |
|
|
0 = -Д- [arctg-^------ arctg—-—] . |
(49) |
||||
• |
2л ( |
+ т |
х—т\ |
' |
Фиг. 50. Картина течения, образуемого источником и стоком: а — величины т- d и Л; б —- линии ? и ф.
Линиями тока t|> = С2 являются окружности, проходящие через точки х — + т с центрами, находящимися на оси оу. Положения
центров h и радиусы |
г этих |
окружностей находятся из формул |
|||||
|
|
|
|
/i |
= mctg^>; |
(50) |
|
|
|
|
|
R = -^. |
(51) |
||
|
|
|
|
|
|
sin ф |
|
Здесь |
|
2л0 |
|
|
|
~ 2л<р |
2лС1 % |
— |
= |
2лС2 |
= С, и |
||||
|
|
|
|
= —= С,. |
|||
Уравнение окружностей |
имеет вид |
|
|||||
|
|
|
|
|
|
|
\2 |
120
Изопотенциальные линии представляют собой семейство окруж
ностей с центрами, расположенными на оси ох. |
Положение центров d |
|
и радиусы R этих окружностей определяются при помощи формул |
||
(фиг. 50, б) |
|
|
d — т cth |
|
(52) |
R = sh <р . |
|
(53) |
Уравнение окружностей ср — Ci имеет вид |
|
|
(х — meth св)2 + у2 — /—. |
(54) |
|
\ sh <р |
/ |
|
При изменении h от — со до 0 эксцентриситет |
е двух окружно |
стей <|> = С2 радиусов Ri и R2 определяется по формуле
s = yrR2i — m2 — R2— m2 ,
откуда
tn — |
4е2 |
(55) |
Из формулы (55) следует, что |
|
|
7?2<^7?1<С7?24*£- |
|
|
Предельное значение R (пг |
< R <оо), при |
котором выполняются |
соотношения (55) и (56), равно т, при этом h = 0. Окружность радиуса R = т изображена на фиг. 50, б жирной линией. Для поло
жительных величин |
h (0< h < со) |
радиусы окружностей возра |
||
стают с увеличением |
h. |
|
|
находящихся |
Далее, для двух каких-либо дуг окружностей, |
||||
вне основного круга, эксцентриситет определяется по формуле |
||||
е — |
/?2 — tn2, — |
R2 — т2 , |
|
|
откуда |
|
|
|
|
|
|
|
|
(57) |
Из равенства (57) |
для случая Ri < R2 получаем |
|
||
|
|
/?2>^1>^2-S- |
(58) |
|
Неравенствами (56) |
и (58) определяются, таким |
образом, два |
типа каналов, соответственно имеющих
> Т?2
и
< R%-
121