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

книги из ГПНТБ / Жуковский М.И. Расчет обтекания решеток профилей турбомашин

.pdf
Скачиваний:
18
Добавлен:
30.10.2023
Размер:
15.1 Mб
Скачать

Продолжение табл. 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

X

У

т

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) принимает вид

С =

(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)

 

 

т

— z

 

 

 

После разделения выражения (47) на вещественную и мнимую

части, получаем:

 

 

 

 

 

 

Q

1п

У(х-гтУ + у\

 

 

 

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