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

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

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

потоку вдоль обводов профилей. В области безотрывных режимов хорошо отработанные компрессорные решетки профилей указанных типов имеют малые потери энергии.

Вопрос об оптимальных эпюрах скоростей для компрессорных решеток разработан слабо, известны лишь отдельные исследова­ ния [97].

Степень диффузорности эпюры скоростей на выходном участке выпуклой стороны профиля зависит главным образом от отношения й1/&2, характеризующего диффузорность канала (фиг. 128), и от кри­ визны обвода профиля. Указанное отношение bjb2 в свою очередь определяется треугольником скоростей.

Представляет интерес вопрос о типе диффузорной эпюры скоро­ стей. Как показано в работе [97], выпуклая эпюра скоростей имеет некоторые преимущества по сравнению с вогнутой (фиг. 129), в осо­ бенности при больших положительных углах атаки. Выбор того или иного типа эпюры зависит от числа М и от треугольника скоро­ стей, реализуемого решеткой.

Рассмотренные типы эпюр скоростей для турбинных и компрес­

сорных решеток профилей рекомендуются для проектирования реше­ ток, обтекаемых, в основном, с дозвуковыми скоростями.

33. ПОСТРОЕНИЕ ГРУППЫ РЕШЕТОК ПРОФИЛЕЙ ТУРБОМАШИН

Основными параметрами, определяющими при всех прочих рав­ ных условиях тип профиля, являются расчетные углы натекания

и выхода потока j32. С целью получения достаточно простого инже­ нерного .способа проектирования решеток выделим сначала влияние основного параметра 1 2. Проектирование решетки в соответствии с заданным расчетным углом натекания осуществим затем путем

последующей вариации входной части профиля. Такой подход

к разработке профилей позволит избежать сложного интерполиро­ вания функций, определяющих обтекание решетки, по двум пере­ менным.

Поставим задачу о построении некоторой конструктивной схемы,

которая позволяла бы на основе обобщения способа построения одного или двух профилей проектировать профили на заданные промежу­ точные условия. Проектируемые профили вместе с исходными будем называть группой профилей [87].

В целях систематизации профилей проектируемые таким способом обводы должны быть в конструктивном отношении однотипными.

Поэтому основное условие, которому должно удовлетворять данное построение, заключается в абсолютной «жесткости» конструктивной схемы. При этом обеспечивается возможность получения профилей единственным образом в зависимости от заданного значения 2, что позволяет путем интерполирования решения задачи об обтекании

1 Точность, с которой проектируемые решетки должны обеспечивать угол выхода потока, составляет 15 ч- 30'; в отношении же угла натекания допустимы во многих случаях отклонения до нескольких градусов.

238

решеток как потенциальным потоком, так и потоком реальной жидкости достаточно быстро определять эпюры скоростей и экспери­ ментальные характеристики, если они известны для нескольких

исходных решеток данной группы.

Вследствие периодичности функций, определяющих решение задачи об обтекании решетки, все построения осуществляются при фиксированном значении шага решетки. Изложение метода прово­

дится на примере построения группы турбинных решеток профилей,

однако метод полностью применим и для проектирования решеток других типов.

Примем для разработки группы в качестве исходного профиль ТЗА. Построим серию однотипных профилей так, чтобы точки сопря­ жения дуг АВ и ВС, ВС и CD (фиг. 130) всех профилей лежали на некоторых прямых, уравнения которых заданы. Центры дуг АВ,

ВС и CD для всех профилей поместим также на некоторых заданных прямых х. Далее, диаметры входной и выходной кромок DE и AG выберем для всех профилей одинаковыми. Точка А сопряжения дуг АВ и AG является точкой сопряжения соответствующих дуг для

1 Выбор прямых линий для геометрического места характерных точек данного построения сделан с целью упрощения метода.

239

всех профилей; центры окружностей входных кромок поместим на вертикали, проходящей через точку, соответствующую центру дудуги DE для профиля ТЗА.

Задавая различные значения р2ак, получим выпуклые стороны

для нескольких профилей, форма которых определится единствен­ ным образом и будет зависеть только от угла |32вк (фиг. 106), харак­ теризующего при фиксированном шаге определенный угол выхода потока.

Построение вогнутой части профиля выполняется аналогично.

Точки сопряжения дуг EF и FG находятся на одной прямой с цен­ трами дуг FG.

Выбор направлений прямых, проходящих через точки В и С,

а также прямых, соединяющих центры всех окружностей, опреде­

ляется типом профилей, их толщинами и т. п. Построение группы профилей осуществляется особенно просто, если профили, выбран­ ные в качестве исходных, являются однотипными.

Получаемые по данному способу решетки профилей имеют опти­ мальный угол натекания, зависящий также от угла выхода потока 2.

При помощи дальнейшей вариации формы входной части профиля по способу, изложенному в п. 34, удается без проведения дополни­ тельных расчетов получать набор решеток для оптимальных углов натекания, изменяемых в достаточно широких пределах. Решетки,

проектируемые по методам пп. 33 и 34, обеспечивают с достаточной для практики точностью заданный треугольник скоростей и имеют благоприятное обтекание, близкое к обтеканию исходных решеток.

Рассмотрим подробнее построение группы решеток профилей на численном примере. Ширину профиля ТЗА возьмем равной 50 мм.

Начало координат выберем в точке А. Относительный шаг положим равным 0,6; при этом в принятом масштабе абсолютное значение

шага равно 30,828 мм. В соответствии со сказанным выше зададим

следующие уравнения для прямых, проходящих через точки В и С.

Прямую, проходящую через точку В, проведем под углом 135°

коси х; ее уравнение имеет вид 1

у= — х + 34,075.

Прямая, проходящая через точку С, проводится под углом 30е

коси х\ уравнение этой прямой запишется так

у= 0,577 х + 0,797.

Центры для дуг АВ лежат на прямой, параллельной оси х; ее уравнение

у = —22,440.

Абсциссы центров и радиусы находятся из формул (фиг. 130)

а! = 22,440 ctg 2gK;

п _ 22,440

1 ~ sin р2вк •

1 Все величины х, у и свободные члены уравнений даются в миллиметрах.

240

Центры для дуг ВС находятся на прямой, составляющей с осью х угол + Ю°); уравнение этой прямой имеет вид

у = — 5,671х + 168,366.

Центры дуг CD лежат на прямой

у = — 0,577х + 12,946,

составляющей с осью х угол

+ 60°^.

Центры дуг, оформляющих

входные кромки, лежат на верти­

кали

х = 49,056.

Уравнение линии, на которой лежат центры дуг EF и FG, имеет

вид

у = — 1,732х + 36,852.

Указанная прямая составляет с осью х угол (^ + 30°). Радиус Re,

дуги FG и координаты центра а6, Ь6 находятся из условия получения определенной толщины профиля или угла выхода р2-

В данном случае для радиуса дуги FG принимаем формулу

R= 391’^°—375,719.

ьcos 2вх

Наконец, координаты центров аь и Ь5 и радиусы R5 дуг EF находятся

при помощи формул:

а _ *6 + j/б — а4 — -I- — 2 (Уб — bi) (Уб — kx6) + 2/?.,х6 V 1

+ /г2

2[х6 — at + (ye — bi)k-\

1 -{ ^21

*5 = Уб — Ue — as)

^6 = IUe-«5) (/bbF|;

k = tg<p6,

здесь заданными величинами являются радиус 7?6 дуги FG, коор­ динаты центра этой дуги а6, Ь6, радиус дуги входной кромки

и координаты центра этой дуги а4, by, кроме того, известными яв­ ляются координаты точки F (х6, у(>) и <р6 — угол наклона прямей FO6.

Приведенное построение является примерным, и оно может изме­ няться в зависимости от задаваемых условий. Так, для получения более тонких профилей, уклон прямой, соединяющей центры дуг FG, следует уменьшить. При необходимости получения профилей с дру­

гой кривизной дуги АВ следует изменить уклон прямой 1—1. При необходимости проектирования профилей с выходной кромкой, отличной от принятой для этой группы, используется метод, изло­

женный выше (п. 31).

Пусть обтекания двух решеток ТЗА и Т4, принятых для построе­ ния данной группы, известны (табл. 27). Основные размеры, опре­ деляющие обводы этих профилей, даны в табл. 28. Значения вели­ чин р2, 2вк> и 7 для этих решеток даны в табл. 29.

16 М. И. Жуковский

700

241

 

 

 

Результаты расчета обтекания решеток

 

 

Профиль ТЗА

 

 

 

7 = 0,6;

q = 0,846;

= 27°30';

р2 23°50'

 

X

У

ds

S

 

 

 

 

 

 

 

 

Wz

0

—0,059

0,201

0,5189

0,048

—2,690

10

—0,095

0,290

0,5799

0,070

—2,712

20

—0,137

0,387

0,6228

0,093

—2,702

30

—0,183

0,486

0,6319

0,117

—2,768

40

—0,235

0,586

0,6581

0,142

—2,725

50

—0,294

0,687

0,6790

0,168

—2,749

60

—0,358

0,787

0,6945

0,195

—2,931

70

—0,443

0,884

0,7986

0,224

—2,949

80

—0,567

0,975

0,9711

0,259

—2,832

90

—0,734

1,042

1,0799

0,299

—2,727

100

—0,920

1,061

1,0382

0,340

—2,642

ПО

— 1,085

1,034

0,8819

0,377

—2,642

120

—1,214

0,978

0,7392

0,409

—2,680

130

—1,314

0,910

0,6493

0,436

—2,731

140

—1,391

0,834

0,6043

0,460

—2,750

150

— 1,456

0,753

0,5851

0,483

—2,707

160

—1,511

0,671

0,5421

0,505

—2,754

170

—1,553

0,591

0,5015

0,525

—2,721

180

—1,586

0,513

0,4593

0,544

—2,579

190

—1,612

0,446

0,3452

0,560

—2,738

200

—1,622

0,399

0,2233

0,571

—2,874

210

—1,613

0,368

0,1587

0,578

—1,680

220

—1,589

0,360

0,2086

0,585

0,898

230

—1,535

0,393

0,5221

0,600

1,325

240

—1,420

0,466

1,0105

0,630

1,351

250

—1,217

0,548

1,5054

0,679

1,361

260

—0,921

0,586

1,8935

0,746

1,398

270

—0,587

0,523

1,9535

0,820

1,484

280

—0,310

0,369

1,6359

0,890

1,614

290

—0,142

0,198

1,1166

0,944

1,813

300

—0,061

0,072

0,6096

0,977

2,149

310

—0,025

0,010

0,2382

0,994

2,681

320

—0,007

—0,003

0,0825

0

0,676

330

0

0,014

0,1609

0,005

—2,682

340

—0,010

0,054

0,3242

0,014

—2,567

350

—0,032

0,120

0,4549

0,029

—2,531

242

Таблица 2

'.профилей ТЗА и Т4

Профиль Т4

7 = 0,595; q = 0,781; = 45°30'; ?2 = 32°20'

й°

 

ds

 

S

 

 

X

У

 

</b

 

L

 

 

 

 

 

 

0

—0,057

0,133

 

0,4529

 

0,038

 

10

—0,097

0.209

 

0,5235

 

0,059

 

20

—0.144

0,292

 

0,5674

 

0,082

 

30

—0.198

0,378

0,5956

 

0,101

 

40

-0.260

0,465

 

0,6412

j

0,133

 

 

—0,332

 

 

0,162

 

50

0,556

 

0,6780

 

 

60

- 0,412

0.646

 

0,6985

 

0,191

 

70

—0,504

0,730

 

0,7417

 

0,222

 

80

—0,618

0,803 '

 

0,8146

 

0,255

 

90

—0,757

0.854

 

0,8718

'

0,291

!

100

- 0,907

0,876

 

0,8597

 

0,328

 

 

 

 

ПО

1,051

0,866

 

0,7912

 

0,363

;

120

— 1,177

0,831

 

0,7100

,

0,395

 

130

— 1,283

0,781

 

0,6374

 

0,424

 

140

-- 1,372

0,723

 

0.5892

 

0,450

 

150

- 1.448

0,659

 

0,5377

 

0,474

 

160

— 1,510

0,596

 

0,4753

 

0,496

 

170

—-1,559

0,535

 

0,4282

 

0,515

 

180

— 1.598

0,479

 

0,3354

 

0,532

 

190

— 1.621.

0,440

 

0,1922

 

0,543

 

200

— 1.620

0,413

 

0,1451

 

0,550

 

210

— 1.600

0,398

 

0,1604

 

0,557

 

220

--1,559

0,409

 

0,3813

 

0,568

 

230

- 1,470

0,449

 

0,7365

 

0,592

 

240

— 1,319

0,500

 

1.0926

 

0,631

 

250

-1,102

0,534

 

1,4242

 

0,685

 

260

—0,833

0.518

 

1,6355

 

0,750

 

270

—0,561 "

0.434

 

1,5924

 

0,819

 

280

—0,340

0.304

 

1,3251

 

0,882

 

290

-0,185

0,175

 

0,9842

 

0,931

1

300

—0,088

0,072

0,6416

 

0,966

 

310

- -0,035

0,007

 

0,3219

 

0,986

 

320

—0.010

—0,017

 

0.1061

 

0,996

 

330

0

—0,010

 

0,1023

 

0

 

340

—0,006

0,018

 

0,2415

 

0,007

 

350

-О',027

0.068

 

0,3680

 

0,020

 

W

wz

—2,032

—2,074

—2,137 i

—2,197 |

—2,166 | i

—2,178 }

—2,287 2,373

—2,360

—2,282

—2,209

—2,150 !

—2,097 |

—2,054 j

—1,974 |

— 1,918 j

--1,886

1,742

1,707

1,938

—0,986 0,745 1,097 1,027 ; 1.032 i

1,052 |

1,097 i

1,196

1,337

1.464

1,588

1,829

1,813

1,540

1,900

1.936

И?

243

 

Координаты обводов профилей ТЗА и Т4

Таблица 28

 

 

Номер точки

Координаты точек сопря­

Координаты центров

Радиусы

сопряжения

жения дуг

 

 

 

 

дуг

Хп

Уп

 

 

Ьп

Rn

 

 

 

 

 

Профи. I ь ТЗА

 

 

 

1

0

0

79,444

- 22,440

82,552

9

10,743

23,332

27,548

12,136

20,193

3

42.720

25,461

17,179

 

3,028

33,994

4

49,984

11Л 38

49,056

11,686

0,961

5

48,456

10,936

28,998

— 13,373

31,138

6

0,748

—0,277

0,385

—0,109

0,400

 

 

Профи л ь Т4

 

 

 

1

0

0

56,968

—22,140

61,228

2

13,447

20,627

28,714

 

5519

21,479

3

39,919

23,844

24,579

— 1,244

29,406

4

49,883

13,737

49,056

 

13,247

0,961

5

48,601

12,400

31,981

— 18,540

35,123

6

14,420

11,876

36,970

—27,181

45,100

7

0,694

—0,384

0,372

—0,147

0,400

 

 

 

 

 

 

Таблица 29

Профиль

 

 

,

мм

1

t

 

 

г,

:

 

 

 

 

 

1

0.6

ТЗА

23°50'

15°46'

30,828

|

Т4

32°20'

21°30'

30,828

|

0,595

 

Координаты обвода

профиля Т

 

 

Таблица 30

 

 

 

 

Номер точки

Координаты точек со­

Координаты

центров

Радиусы

j сопряжения

пряжения дуг

 

 

 

 

Дуг

хп

Уп

ап

 

Ьп

Rn

 

 

1

0

0

72,269

—22,440

75,673

2

11,466

22,609

27,840

 

10,477

20,378

3

42,056

25,078

19,361

 

1,769

32,533

4

49,960

12,816

49,056

 

12,490

0,961

5

48,512

11,697

30,045

—15,188

32,616

6

13,737

13,059

30,868

— 16,612

34,262

7

0,734

—0,309

0,382

—0,119

0,400

244

Спроектируем при помощи данного построения (фиг. 130) решетку

профилей Т с углом выхода р2 = 26°. Интерполируя зависимость р2вк

от р2 для решеток ТЗА и Т4,

найдем,

что для угла

2 равного 26°,

угол р2вк составляет 17°15'.

Обвод

профиля,

соответствующий

этому значению 2вк, показан

штрих-пунктиром

на

фиг. 130. Зна­

чения координат центров дуг, описывающих профиль Т, и их радиусы,

приведены в табл. 30. Распределение скоростей, определенное рас-

Фиг. 131. Эпюры скоростей для решеток:---------------------

ТЗА;

--------------------Т;----------------------Т4; .............

—интерполя­

ция для решетки Т.

 

где нанесены также величины скоростей для решеток ТЗА и Т4.

Угол выхода оказался равным 26°10', что отличается от задания

на 10'. Сравнение профилей дано на фиг. 130. Результаты расчета обтекания решетки Т приведены в табл. 31. Скорости, вычисленные при помощи интерполяции даны в табл. 32.

Рассчитаем распределение скоростей обтекания данной решетки интерполяционным путем при помощи эпюр скоростей для реше­ ток ТЗА и Т4. В данном случае пользуемся линейной интерполяцией, так как известны лишь два распределения скоростей для решеток ТЗА и Т4 \ При помощи табл. 29 находим 21

—** *

17°15' — 15°46'

-------------------= O.ZDD.

2вк2вк

21°30' — 15°46'

Далее вычисляем значения густоты q решетки кругов

— 0,255

1 Если известны эпюры для нескольких исходных решеток, то можно интерпо­ лировать с учетом более высоких разностей.

2 Величины, относящиеся к решеткам ТЗА и Т4А, отмечены соответственно значками * и **.

245

 

Результаты

расчета обтекания

решетки Т

 

9 = 0,828;

рх — 33°;

2 = 26°101;

г = 0,597

0

X

У

ds

S

 

т

 

 

 

 

0

—0,059

0,171

0,4984

0,047

10

—0,094

0,255

0,5529

0,068

20

—0,138

0,347

0,6088

0,091

30

—0 186

0,444

0,6261

0,116

40

—0,242

0,540

0,6576

0,142

50

—0,304

0,640

0,6779

0,168

60

—0,373

0,737

0,6920

0,196

70

—0,460

0,830

0,7858

0,225

80

—0,580

0,917

0,9110

0,259

90

—0,734

0,980

0,9935

0,298

100

—0,907

1,001

0,9810

0,337

110

—1,066

0,982

0,8514

0,374

120

—1,194

0,934

0,7342

0,405

130

—1,297

0,870

0,6570

0,433

140

—1,380

0,800

0,5970

0,458

150

— 1,448

0,723

0,5739

0,482

160

— 1,506

0,645

0,5312

0,504

170

—1,552

0,571

0,4800

0,524

180

—1,587

0,499

0,4296

0,542

190

—1,615

0,440

0,3000

0,557

200

—1,624

0,402

0,1809

0,566

210

— 1,610

0,376

0,1597

О',573

220

-1,579

0,375

0,2599

0,582

230

— 1,515

0,414

0,6328

0,600

240

— 1,383

0,482

1,0728

0,634

250

— 1,167

0,545

1,4955

0,685

260

—0,878

0,558

1,8032

0,751

270

—0,567

0,482

1,8237

0,824

280

—0,313

0,333

1,5166

0,891

290

—0,152

0,175

1,0593

0,942

300

—0,067

0,057

0,6176

0,976

310

—0,027

—0,005

0,2521

0,993

320

—0,008

—0,019

0,0904

0

330

—0,002

—0,005

0,1345

0,004

340

—0,008

0,030

0,2999

0,013

350

—0,032

0,092

0,4473

0,028

Таблица 31

N 1 а

—2,480

—2,547

—2,509

—2,564

—2,531

—2,573

—2,743

—2,748

—2,705

—2,625

—2,498

—2,490

—2,472

—2,448

—2,470

—2,385

—2,364

—2,323

—2,177

—2,364

—2,393

—0,624 1,139 1,207 1,212

1,250

1,305

1,399

1,543

1,728

1,972

2,490

1,180

—2,497

—2,348

—2,242

246

 

 

 

 

 

 

 

 

 

 

 

Таблица 32

 

 

 

 

Интерполяция скоростей для решетки Т

 

 

 

 

 

 

 

 

q = 0,829;

= 33°;

= 26°;

1 = 0,597

 

 

 

 

е

ds

s

w

e

ds

s

w

9

ds

s

w

 

~dh

~L

Wz

 

Иь

~~L

wz

 

db

Wz

0

0,5021

GJMl

—2,479

120

0,7318

0,406

—2,487

240

1,0314

0,632

1,268

10

0,5655

0,068

—2,507

130

0,6463

0,434

—2,493

250

1,4847

0,682

1,267

20

0,6087

0,091

—2,524

140

0,6005

0,459

—2,458

260

1,8277

0,748

1,296

30

0,6227

0,116

—2,591

150

0,5730

0,482

—2,391

270

1,8615

0,821

1,380

40

0,6538

0,141

—2,556

160

0,5251

0,504

—2,392

280

1,5567

0,889

1,513

50

0,6788

0,168

—2,580

170

0,4828

0,524

—2,310

290

1,0829

0,942

1,699

60

0,6955

0,195

—2,739

180

0,4277

0,542

—2,185

300

0,6177

0,976

1,979

70

0,7841

0,225

—2,767

190

0,3062

0,557

—2,312

310

0,2595

0,993

2,420

80

0,9312

0,259

—2,662

200

0,2034

0,567

—2,121

320

0,0885

0

1,171

90

1,0269

0,298

—2,556

210

0,1591

0,574

—0,613

330

0,1460

0,005

—2,335

100

0,9927

0,338

—2,483

220

0,2526

0,582

1,183

340

0,3031

0,014

—2,347

НО

0,8583

0,375

—2,481

230

0,5991

0,599

1,282

350

0,4327

0,028

—2,335

и производной

Эпюра скоростей, вычисленная таким путем, дана на фиг. 131.

Для сравнения нанесены величины скоростей, полученные точным расчетом. Соответствие эпюр является вполне удовлетворительным.

Результаты вычислений скоростей w (s) для решетки Т при помощи интерполяции приведены в табл. 32.

Аналогично при помощи интерполяции можно определять и экс­ периментальные характеристики, если они известны для исходных решеток.

Указанный в данном параграфе метод профилирования может

быть использован в некоторых случаях для проектирования только выходной части профиля в соответствии с заданным углом выхода потока. Проектирование входной части может вестись при этом в зависимости от заданной толщины профиля путем использования формулы (4) для обеспечения оптимального угла натекания.

Следует отметить, что при решении задачи о сведении сложной зависимости формы турбинных профилей к одной переменной 2 (или р2вк) не Делалось никаких предположений о характере потока. Следовательно, указанный метод принципиально применим для про­ филирования решеток, используемых при больших скоростях потока, а также для профилирования решеток других типов. Таким образом,

задача о проектировании решеток профилей различных типов сво­ дится к построению нескольких групп профилей.

Для построения группы в качестве исходных профилей могут

быть взяты совершенные в аэродинамическом отношении решетки, разработанные при помощи методов, изложенных в гл. II и III,

или решетки, отработанные экспериментальным путем. Две исход­ ные решетки на минимальный и максимальный углы выхода потока следует выбрать с благоприятным обтеканием. При этом обтекание промежуточных решеток будет также благоприятным, и уклон эпюры скоростей для промежуточных профилей на соответственных участках обвода будет ограничен сверху и снизу величинами уклонов эпюр для исходных профилей, т. е. 1

dw*

dw

.. dw**

' '

ds

ds

ds ’

Выполнение этого неравенства проиллюстрируем для выходной части выпуклой стороны профиля, обтекание которой представляет наибольший интерес.

1 Здесь значками * и ** отмечены величины производных для исходных про­ филей, причем значки * и ** относятся соответственно к профилям с наибольшим и наименьшим уклонами.

248