Проектирование газотурбинного двигателя для многоцелевого легкого боевого самолета

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Введение

Для современной авиации характерно применение различных типов газотурбинных двигателей. Это объясняется разнообразием типов самих летательных аппаратов и специфическим требованиям, предъявляемым каждым типом летательного аппарата к его силовой установке.

Наиболее простым и по этой причине первым, получившим широкое применение в авиации, газотурбинным двигателем был турбореактивный двигатель. Такие двигатели используются на самолетах с околозвуковыми максимальными скоростями полета, но при высокой температуре газа перед турбиной их применяют при скоростях, соответствующих числу M полета не более 2.0.

На самолетах с большой потребной тяговооруженностью и большими скоростями полета (М=2.0 - 3.5) широко используются турбореактивные двигатели с форсажной камерой (ТРДФ), в которой к потоку газов, покидающих турбину, подводится дополнительное тепло, образующееся при сжигании топлива. Это позволяет значительно повысить температуру газа перед реактивным соплом ТРДФ и, соответственно, существенно увеличить скорость истечения газов из двигателя и получить более высокую тягу, чем у ТРД.

Для этого вида двигателей характерно наличие турбокомпрессора - агрегата, состоящего из компрессора, камеры сгорания, турбины и форсажной камеры. В современных ТРД преимущественно применяются осевые компрессоры и турбины, хотя имеются двигатели (в основном маломощные), в которых используются центробежные или диагональные компрессоры и радиальные турбины.

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

Переход к более высоким температурам требует одновременного повышения степени повышения давления. В связи с этим широкое применение нашла схема двигателя с двухвальным газогенератором (ДТРД), основной задачей которой является повышение газодинамической устойчивости двигателя. Но и на этом этап совершенствования не может быть завершен, т.к. с развитием авиации растет потребность к более экономичным двигателям. Для этого наиболее выгодной является схема двухконтурного двигателя с двухвальным газогенератором (ДТРДД).

Целью проекта является проектирование газотурбинного двигателя для многоцелевого легкого боевого самолета. В качестве прототипа был взят двухвальный турбореактивный двухконтурный двигатель с форсажной камерой АИ-222К-25Ф, разработки «Ивченко-Прогресс». Особенности влияния расчетных значений степени повышения давлений, и температуры на удельную тягу и удельный расход топлива двигателя прослеживаются при выполнении данного проекта.

1. Профилирование лопатки первой ступени компрессора высокого давления

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

При учебном проектировании расчет решеток рабочего колеса и их лопаток проводят на трех характерных радиусах.

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

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

1.1 Исходные данные

В качестве исходных данных профилирования примем параметры потока и размеры проточной части РК дозвуковой ступени КВД. Основные исходные данные занесем в таблицу 1.1.

Таблица 1.1 Основные исходные данные для расчета

Параметры	Размерность	Сечение
		вт	ср	к
D	м	0,3178	0,3747	0,424
r/rк	-	0,7495	0,883726	1
Uср	м/с	-	374,57	-
U	м/с	322,07	374,57	430
С1аср	м/с	-	160	-
С2аср	м/с	-	158	-
С1uср	м/с	-	139,6	-
С2uср	м/с	-	224,5	-
скср	-	-	0,497197	-
Нтср	Дж/кг	-	16596,65	-
Т1*	К	418,2	418,2	418,2
Т2*	К	455,1	455,1	455,1
Нтср	Дж*с^2/кг*м^2	0,16	0,118292	0,08976
Принятые значения и допущение
Параметры	Размерность	Величина
k	-	1,4
R	Дж/кг*К	287
При расчете Дб принимаем б3=б1

Для последующего компьютерного расчета, необходимо рассчитать:

­ .

1.2 Выбор закона закрутки

Выбираем закон изменения параметров по радиусу (законы закрутки потока). Критерием выбора оптимального закона закрутки по радиусу является обеспечение дозвуковых скоростей и приемлемых углов потока (в частности, Mw1 и Mc2 0,85…0,90, в1 25o на периферии, в2 90о на втулке). Для расчета первой ступени проектируемого компрессора примем закон «твердое тело» (на входе) при заданном Нт(r).

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

1.3 Компьютерный расчет лопатки компрессора

Файл исходных данных OCK.DAT:

Таблица 1.2

08 12 07 1 3 1.400 287.00 ( дата, M, Ks, kг, Rг )

1.304 430.000 .612 .892 .985 .986 .749 .765

1.000 160.000 418.200 334000.0 139.600 0.424 .964

1.050 .950

_ _ _ _

Пi* Uк Hтвт КПД* Sна D21 d1вт d2вт

m С1асp Т1* P1* С1uсp D1к Кн

W1к/W1сp W1вт/W1сp

Тип ступени : 1 - дозвуковая ступень;

( M ) 2 - свеpхзвуковая ступень .

Закон кpутки : 1 - пеpвая ступень без ВНА ( С1u=0., А=В=D=0.);

( Ks ) 2 - " свободный вихpь "(на входе) пpи заданном Нт(r);

3 - " твеpдое тело " (на входе) пpи заданном Нт(r);

4 - Rок=const пpи заданном Нт(r) ;

5 - по значениям W1к/W1сp и W1вт/W1сp.

Расчет проводим с помощью программы OCK.EXE.

Файл полученных результатов OCK.REZ:

Таблица 1.3

ГДР СТ.ОК ДАТА 8.12. 7

ИСХОДНЫЕ ДАННЫЕ

M= 1 KR= 3 КГ= 1.400 RГ= 287.00

1.304 430.00 .612 .892 .985 .986 .749 .765

1.000 160.00 418.20 334000. 139.60 .424 .964

1.050 .950 А= .367 B= .000 D= .000

ГАЗОДИНАМИЧЕСКИЙ РАСЧЕТ СТ. ОК

PI1=1.304 PI2=1.304 HZ1=37105. HZ2=37105. T01=455.14 T02=455.14

P01= 435536. P02= 435536.

( GB= 23.346 ROK= .4992 HTO= .2082 WC= 19368.8 )

Таблица 1

N U CU CA T0 T P0 P

ROTH RO C W LC LW AL BE

11 430.00 158.02 120.99 418.20 398.48 334000. 282062.

1.0000 2.4663 199.02 297.68 .5318 .7733 37.441 23.982

12 401.35 147.49 145.16 418.20 396.88 334000. 278115.

.9334 2.4416 206.94 292.43 .5530 .7624 44.546 29.762

13 374.57 137.65 163.35 418.20 395.49 334000. 274704.

.8711 2.4202 213.61 287.78 .5708 .7526 49.882 34.585

14 348.38 128.02 178.32 418.20 394.21 334000. 271624.

.8102 2.4008 219.51 283.47 .5866 .7434 54.323 38.980

15 322.07 118.35 191.21 418.20 393.03 334000. 268774.

.7490 2.3828 224.88 279.40 .6009 .7347 58.244 43.186

* * * * * * * *

21 423.98 251.04 93.86 455.14 419.38 442169. 332067.

.9860 2.7589 268.01 196.76 .6865 .5134 20.499 28.490

22 394.16 247.83 139.66 455.14 414.86 442169. 319691.

.9166 2.6850 284.47 202.28 .7287 .5299 29.404 43.665

23 369.49 243.71 171.29 455.14 410.97 442169. 309329.

.8593 2.6226 297.88 212.50 .7631 .5578 35.101 53.711

24 346.73 239.64 196.80 455.14 407.27 442169. 299703.

.8064 2.5640 310.09 224.05 .7943 .5889 39.393 61.447

25 324.73 235.91 218.89 455.14 403.59 442169. 290311.

.7552 2.5064 321.82 236.22 .8244 .6216 42.856 67.914

* * * * * * * *

81 423.98 155.80 126.60 418.20 398.14 334000. 281207.

.9860 2.4610 200.75 296.56 .5365 .7709 39.096 25.271

82 394.16 144.84 150.39 418.20 396.50 334000. 277173.

.9166 2.4357 208.80 291.16 .5580 .7597 46.076 31.099

83 369.49 135.78 166.45 418.20 395.23 334000. 274087.

.8593 2.4163 214.80 286.92 .5740 .7508 50.795 35.459

84 346.73 127.42 179.18 418.20 394.14 334000. 271437.

.8064 2.3996 219.87 283.21 .5876 .7429 54.584 39.249

85 324.73 119.33 189.99 418.20 393.14 334000. 269052.

.7552 2.3845 224.36 279.80 .5996 .7356 57.867 42.768

Рис. 1.1 График изменения и по высоте лопатки на входе и выходе из рабочего колеса



Рис. 1.2 График изменения Lc и Lw по высоте лопатки на входе и выходе из рабочего колеса

Рис. 1.3 График изменения и по высоте лопатки на входе и выходе из рабочего колеса

Рис. 1.4 Треугольники скоростей в 1-м сечении лопатки рабочего колеса



Рис. 1.5 Треугольники скоростей во 2-м сечении лопатки рабочего колеса

Рис. 1.6 Треугольники скоростей в 3-м сечении лопатки рабочего колеса

Рис. 1.7 Треугольники скоростей в 4-м сечении лопатки рабочего колеса



Рис. 1.8 Треугольники скоростей в 5-м сечении лопатки рабочего колеса продолжение таблицы 1.3

Таблица 2

Профилирование лопатки РК по радиусу

----------------------------------------------------------

Паpаметp | Сечение по высоте лопатки

| 1(пеp) 2 3(сp) 4 5(вт)

----------------------------------------------------------

ro .9860 .9166 .8593 .8064 .7552

b 26.55 26.55 26.55 26.55 26.55

t 24.78 23.04 21.60 20.27 18.98

b/t 1.071 1.152 1.229 1.310 1.399

Cm .0300 .0400 .0500 .0600 .0700

i .0000 .0000 .0000 .0000 .0000

del 1.666 5.400 6.850 7.432 7.588

dbe 3.219 12.57 18.25 22.20 25.15

tet 4.885 17.97 25.10 29.63 32.73

be1l 25.27 31.10 35.46 39.25 42.77

be2l 30.16 49.07 60.56 68.88 75.50

Число рабочих лопаток - 53. шт.



Рис. 1.9 Профиль лопатки на r=0,986

Рис. 1.10 Профиль лопатки на r=0,917

Рис. 1.11 Профиль лопатки на r=0,859



Рис. 1.12 Профиль лопатки на r=0,806

Рис. 1.13 Профиль лопатки на r=0,755

1.4 Расчет и профилирование лопатки

Расчет производим для первой ступени осевого компрессора высокого давления.

· В таблице 1.4 - расчет Са и Сu при использовании выбранного закона, а также некоторых параметров планов скоростей (скоростей С и W, углов б, в, Дв) и чисел М.

Таблица 1.4

	Расчет Са и Сu(закон закрутки "твердое тело" при заданной Hт(r))
	Параметры	Размерность	Сечение
			вт	ср	к
	С1а	м/с	191,21	163,35	120,99
	С2а	м/с	218,89	171,29	93,86
	ск	-	0,450026	0,490936	0,524349
Нт=Нтср	Нт	Дж/кг	16596,65	16596,65	16596,65
	C1u	м/с	118,35	137,65	158,02
	C2u	м/с	235,91	243,71	251,04

Таблица 1.5

	Расчет некоторых параметров планов скоростей
	Параметры	Размерность	Сечение
			вт	ср	к
	W1	м/с	279,3977	287,7748	297,6772
	C1	м/с	224,8733	213,6135	199,0198
	a1	м/с	427,5678	427,5704	427,5738
	Mw1	-	0,653458	0,673046	0,696201
	W2	м/с	235,2369	215,5565	202,0801
	C2	м/с	321,8173	297,8839	268,0127
	a2	м/с	402,6735	406,34	410,479
	Mc2	-	0,799202	0,73309	0,652927
	б1	град	58,27398	49,9055	37,45892
	б2	град	42,87856	35,11902	20,5104
	в1	град	43,20758	34,60277	23,994
	в2	град	68,54904	52,64804	27,68991
	Дб	град	15,39543	14,78647	16,94852
	Дв	град	25,34146	18,04527	3,695912

При расчете считалось, что б3 ? б1.

Полученные выше кинематические параметры являются исходными для расчета профилей лопаток и решеток РК дозвукового осевого компрессора. Но в виду получения более правильных профилей данными для профилирования будем брать из компьютерного расчета лопатки (п. 1.3, таблица 1.3).

· Расчет густоты решеток профилей

Густота решетки b/t (b - хорда пера лопатки; t - шаг решетки) является параметром, в значительной степени определяющим аэродинамическую нагруженность лопаточного венца. Уменьшение значений густоты решетки в сравнении с оптимальными означает недогрузку ступени, а увеличение густоты отрицательно сказывается на КПД ступени.

Определение густоты решетки РК производится на так называемом номинальном режиме, который характеризуется бессрывным обтеканием решетки при отсутствии резкого роста потерь.

Отношение рассчитанного угла поворота потока к углу на номинальном режиме * определяет запас по срыву компрессорной решетки и выбирается в зависимости от положения ступени в компрессоре. Величину /* для первой ступени КВД принимаем равной 0,9.

· Выбор углов атаки лопаточного венца на номинальном режиме

Расчетные углы атаки i по радиусу лопатки выбирают на номинальном режиме. Отклонение от номинальных углов атаки может привести к существенному снижению КПД и сокращению диапазона устойчивой работы компрессора. В современной практике углы атаки рассчитывают по эмпирическим формулам, учитывающим накопленный статический материал.

В первом приближении допустимо принимать углы атаки i=0 по всей высоте лопатки.

· Расчет углов изгиба профиля пера

Из геометрических соображений угол изгиба профиля пера лопатки равен:

,

где - угол отставания потока в лопаточном венце.

,

- относительный прогиб средней линии профиля. При дуге окружности =0.5.

Исходя из полуэмпирической зависимости для угла отставания потока на номинальном режиме можно определить:

Расчет параметров решетки на среднем радиусе:

Таблица 1.6

Расчет параметров решетки на среднем радиусе
Параметры	Размерность	Величина
Dk	м	0,424
Dcp	м	0,3747
Dвт	м	0,3178
h=( Dk - Dвт)/2	м	0,0531
	-	2
	м	0,02655
(Дв)cp	град	18,045271
(Дв/Дв*)cp	-	0,7256
	град	24,8694474
в2cp	град	52,6480371
	град	11,99
	-	2,07418243
	-	1,196
	м	0,022199
	шт	53,0005035
z	шт	53
	м	0,02219921
	м	0,02655025
	-	1,999981

· Выбор дуги средней линии профиля

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

Принимаем форму средней линии профиля пера лопатки - дуга окружности. Тогда для средней линии по дуге окружности целесообразно использовать соотношение для радиуса дуги:

.

Расчет углов средней линии профиля, углов наклона кромок пера и угла установки профиля в решетке

Углы средней линии профиля лопатки на входе и на выходе:

Углы наклона входной и выходной кромок пера лопатки:

где k=0,5 при средней линии по дуге окружности.

Углы установки (выноса) профиля в решетке:

· Расчет осевого размера лопаточного венца

Осевой размер лопаточного венца на соответствующем радиусе:

.

· Выбор относительной толщины профиля

При выборе относительной максимальной толщины профиля пера по хорде стремятся получить достаточно тонкие профили, необходимые для обеспечения высокого КПД решетки (М<Mкр), но при этом учитывают прочностные характеристики лопатки, особенно в призамковых областях рабочих лопаток. Для коротких рабочих лопаток () можно применять по всей высоте лопатки (порядка 0,05…0,06). Принимаем .

· Построение профилей лопаток и решеток профилей

Построение профилей лопаток состоит из этапов построения средней линии и самого аэродинамического профиля. Из одинаковых профилей, расположенных с заданной густотой b/t под углами установки г к фронту решетки, составляют решетки профилей. Существуют аналитические и графоаналитические методы построения профилей лопаток и решеток профилей. Аналитические методы являются более точными и широко применяются в практике современного турбомашиностроения. Графоаналитические методы целесообразнее для учебных целей в связи с наглядностью и меньшей трудоемкостью вычислений. В данном курсовом проекте будем использовать графоаналитический метод построения профилей лопаток. Расчет параметров лопаток и решеток профилей по радиусу:

Таблица 1.7

Расчет параметров лопаток и решеток профилей по радиусу
Параметры	Размерность	Сечение
		вт	ср	к
b=bcp=const	м	0,02655	0,02655	0,02655
D	м	0,3178	0,3747	0,424
	м	0,018828	0,022199208	0,02512
b/t	-	1,410136	1,196	1,056937
i	град	0	0	0
в2	град	68,54904	52,64803712	27,68991
Xf	-	0,5	0,5	0,5
	-	0,27	0,30	0,35
Дв	град	25,34146	18,04527101	3,695912
(Дв/Дв*)cp	-	0,9	0,9	0,9
Дв*	град	28,15718	20,05030112	4,106569
	град	35,26054	26,40814001	4,742393
	град	8,103358	7,357838881	1,635824
в1	град	43,20758	34,60276611	23,994
	град	43,20758	34,60276611	23,994
	град	76,6524	60,005876	29,32574
k	-	0,5	0,5	0,5
	град	17,63027	13,20407	2,371197
	град	17,63027	13,20407	2,371197
	м	0,083589	0,1132209	0,641499
	м	0,043852	0,058146127	0,321024
	м	0,026982	0,026795017	0,026566
	град	52,02272	41,20480111	25,1796
	м	0,020921	0,017482761	0,011291
C	-	0,08	0,06	0,04
Cbmax	м	0,002124	0,001593015	0,001062
Ar	м	0,01208	0,012375	0,01017
	м	0,012885	0,012600968	0,01021
Ar/A1	-	0,937501	0,982067385	0,996087
M1kp	-	0,66	0,723	0,746
M1max	-	0,746	0,823	0,846
	-	0,703	0,773	0,796
Mw1	-	0,653458	0,673046464	0,696201

· Построение средней линии профиля

При выбранной дуге средней линии пера профиля в виде дуги окружности разбиваем хорду на равное число участков (через 10% всей длины хорды, совпадающей с осью абсцисс). Ординаты средней линии вычисляют по приближенной зависимости:

.

Таблица 1.8

Результаты расчета средней линии
Хср.л,%b	0	10	20	30	40	50	60	70	80	90	100
Хср.л,мм	0	2,6550	5,3101	7,9651	10,6201	13,2751	15,8302	18,4852	21,1402	23,7952	26,550
yср.л вт,мм	0	0,759	1,349	1,771	2,0239	2,10828	2,03018	1,78354	1,3682	0,78427	0
yср.л ср,мм	0	0,56	0,996	1,3075	1,4942	1,55651	1,49885	1,31676	1,0101	0,57901	0
yср.л к,мм	0	0,099	0,176	0,2308	0,2637	0,27471	0,26454	0,2324	0,1783	0,10219	0

Рис 1.12 Средняя линия профиля (втулка)

Рис 1.13 Средняя линия профиля (среднее)

Рис 1.14 Средняя линия профиля (периферия)

· Построение аэродинамического профиля и решетки профилей

В качестве исходного аэродинамического профиля в проекте используем симметричный (yB=yH) профиль А-40 с расположением относительной максимальной толщины профиля равной =0,1, на расстоянии 40% длины хорды от входной кромки профиля ().

Координаты исходного аэродинамического профиля А-40 ( в процентах от длины хорды b, - в процентах от величины максимальной толщины профиля Сbmax) представлены в табл. 1.9;

Таблица 1.9

		Координаты рассчитанного профиля
Х	У	Х,м	Сечение
			Втулочное	Среднее	Периферийное
			Ув=Ун, мм
0	0	0,00	0	0	0
0,01	11,4	0,27	0,242	0,182	0,121
0,015	14,3	0,40	0,304	0,228	0,152
0,025	18,5	0,66	0,393	0,295	0,196
0,05	25,5	1,33	0,542	0,406	0,271
0,075	30,9	1,99	0,656	0,492	0,328
0,1	35,25	2,66	0,749	0,562	0,374
0,15	41,6	3,98	0,884	0,663	0,442
0,2	45,5	5,31	0,966	0,725	0,483
0,25	47,88	6,64	1,017	0,763	0,508
0,3	49,27	7,97	1,047	0,785	0,523
0,35	49,86	9,29	1,059	0,794	0,530
0,4	50	10,62	1,062	0,797	0,531
0,5	48,58	13,28	1,032	0,774	0,516
0,6	44,42	15,93	0,943	0,708	0,472
0,7	37,83	18,59	0,804	0,603	0,402
0,8	28,5	21,24	0,605	0,454	0,303
0,9	17,22	23,90	0,366	0,274	0,183
0,95	10,03	25,22	0,213	0,160	0,107
1	0	26,55	0	0	0

Таблица 1.10

Расчетные радиусы скругления
	Сечение
	Втулочное	Среднее	Периферийное
R1	0,1168211	0,0876158	0,058410555
R2	0,0106201	0,0079651	0,00531005

Вывод

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

В качестве исходного аэродинамического профиля использовали симметричный профиль А-40. При профилировании использовался закон «твердого тела» на входе при заданном Нr.

Рассчитанная и построенная решетка профилей первой ступени осевого компрессора высокого давления удовлетворяет требованиям и сможет обеспечить требуемые параметры. Исходные данные взяты из газодинамического расчета осевого компрессора высокого давления. Полученные профили решетки профилей изображены на рис. 1.4- 1.8.

Полученные числа Мw1 потока во всех сечениях находятся на допустимом диапазоне, т.е. Mw1<M1* - т.е. не происходит запирания решеток профилей.

2. Профилирование лопатки первой ступени турбины высокого давления

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

При учебном проектировании расчет решеток рабочего колеса и их лопаток проводят на трех характерных радиусах.

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

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

2.1 Выбор закона закрутки

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

Примем закон закрутки и . В нашем случае, у этого закона есть ряд преимуществ:

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

2. Высокий КПД

3. Применение этого закона значительно упрощает технологию изготовления лопаток СА и РК, позволяет создать хорошую конструктивную базу для их монтажа.

4. При , лопатки СА и РК первой ступени турбины являются некручеными и имеют почти постоянный профиль по высоте, что способствует организации внутреннего охлаждения.

2.2 Расчет и профилирование лопатки

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

Исходные данные:

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

При законе закрутки и найдем следующие параметры:

1) Радиусы струек тока в сечениях 1-1 и 2-2:

2) Окружные скорости на входе и выходе из рабочего колеса:

3) Окружную составляющую абсолютной скорости на выходе из лопаток соплового аппарата:

4) Расходную составляющую абсолютной скорости на выходе из лопаток СА:

5) Угол потока в абсолютном движении на выходе из лопаток СА:

6) Абсолютную скорость на выходе из лопаток СА:

7) Приведенную скорость перед РК в абсолютном движении:

8) Угол потока в относительном движении на входе в лопатки РК:

9) Относительную скорость на входе в лопатки РК:

10) Окружную составляющую абсолютной скорости на выходе из лопаток РК:

где работа на окружности колеса

11) Окружную составляющую относительной скорости на выходе из лопаток РК:

12) Расходную составляющую абсолютной скорости на выходе из лопаток РК:

13) Угол потока в относительном движении на выходе из лопаток РК:

14) Угол потока в абсолютном движении на выходе из лопаток РК:

15) Относительную скорость потока на выходе из лопаток РК:

16) Абсолютную скорость потока на выходе из лопаток РК:

17) Угол поворота потока в относительном движении:

18) Температуру газа за РК по заторможенным параметрам потока в относительном движении:



19) Приведенную скорость за РК в относительном движении:

20) Термодинамическую степень реактивности:

Результаты расчетов параметров потока в межвенцовых зазорах ступени во втулочном и периферийном сечениях вносим в табл. 2.1:

Таблица 2.1

№ п/п	Параметр	Размерность	Значения параметров при законах закрутки потока и в сечениях
			втулочном	среднем	периферийном
1		м	0,2	0,2165	0,233
2		м	0,2001	0,2225	0,2449
3		м/с	405,54	439	472,46
4		м/с	405,596	451	496,4
5		м/с	696,134	652	613,61
6		м/с	208,85	196	184,09
7		град	16,7	16,7	16,7
8		м/с	726,78	680	640,42
9		-	1,0288	0,657	0,916
10		град	35,71	42,7	52,52
11		м/с	357,81	289,456	237,58
12		м/с	-116,84	-121	-79,498
13		м/с	522,436	572	575,898
14		м/с	182,67	200	201,36
15		град	19,2	19,2	19,2
16		град	57,4	58,8	68,46
17		м/с	533,45	605,96	610,09
18		м/с	216,84	233	216,49
19		град	125,09	118,1	108,28
20		К	1255,88	1281	1266,8
21		-	0,3	0,899	0,342
22		-	0,292	0,38	0,45

Расчет треугольников скоростей в межвенцовых зазорах по высоте можно считать законченными, т.к. полученные параметры во втулочном сечении удовлетворяют условиям , , но на втулке не выполняется условие , что может привести к дополнительным гидравлическим потерям и увеличению нагрузки на лопатку.

По результатам расчетов в масштабе строим треугольники скоростей на трех радиусах (рис 2.1-2.3).



На следующем этапе проектирования определяем геометрические параметры решеток профилей. При этом расчет параметров и построение профилей выполняем для обеспечения закрутки потока по закону и .

1) Радиусы расположения среднего втулочного и периферийного сечений проточной части:

2) Шаг решетки:

3) Хорда профиля лопатки:



4) Угол установки профиля в решетке:

5) Геометрический угол решетки на входе:

, , .

6) Геометрический угол решетки на выходе, принимаем равным углу потока, т.к. :

.

7) Угол отгиба выходной кромки, выбираемый с учетом характера его изменения по высоте:

.

8) Ширина горла межлопаточного канала:

9) Относительная толщина профиля лопатки:

10) Абсолютная толщина профиля лопатки:

11) Относительное удаление максимальной толщины профиля:

12) Абсолютное удаление максимальной толщины профиля:

13) Радиус скругления входной кромки:

13) Радиус скругления выходной кромки:

14) Угол заострения входной кромки:

16) Угол заострения выходной кромки:

17) Угол, образованный лучом, проходящим через центры окружностей радиусами R1 и R2, и фронтом решетки:

Результаты расчетов геометрических параметров заносим в таблицу 2.2:

Таблица 2.2

№ п/п	Параметр	Размерность	Значения параметров в сечениях
			втулочном	среднем	периферийном
1		м	0,20005	0,2195	0,23895
2		м	0,01367	0,015	0,01633
3		м	0,0213	0,0213	0,0213
4		град	71,16	60,6	56,73
5		град	44	52	58
6		град	19,2	19,2	19,2
7		град	18	18	18
8		м	0,00445	0,00493	0,00537
9		-	0,264	0,24	0,216
10		м	0,0056	0,0051	0,0046
11		-	0,275	0,278	0,2788
12		м	0,00586	0,00592	0,00594
13		м	0,001456	0,001326	0,001196
14		м	0,0005	0,0005	0,0005
15		град	44,413	41,007	37,885
16		град	10	10	10
17		град	68,59	58,38	54,86

2.3 Компьютерный расчет лопатки турбины

Файл исходных данных OCT.DAT:

Таблица 2.3.1

18 11 07 1 2 1.300 290. Дата, nr, kz, kг, Rг

.433 .445 .033 .045 D1c,D2c,h1,h2

.949 .955 .957 .380 fi,psi,Л1,Roтc

196.00 200.00 625.00 -121.00 C1ac,C2ac,C1uc,C2uc

16.70 42.70 19.20 20.50 21.00 alf1c,be1c,be2c,G1,G2

19369.6 1260. n,T2*

Лопатка СА - nr=0, лопатка РК - nr=1.

Закон кpутки: 0 - C1u*r=const, C2u*r=const;

( kz ) 1 - alf1(r)=const, L(r)=const;

2 - alf1(r)=const, be2(r)=const.

Для профилирования используем программу OCT.exe.

Файл полученных результатов OCT.REZ:

Таблица 2.3.2

Дата 18.11. 7 KZ= 2 Кг= 1.300 Rг= 290.0

D1ср= .4330 D2ср= .4450 h1= .0330 h2= .0450

C1aср=196.00 C2aср=200.00 C1uср=625.00 C2uср=******

alf1с= 16.70 be1ср= 42.70 be2ср= 19.20

Л1= .957 Фи= .949 Пси= .955 Rтс= .380

n =19369.6 T2*= 1260.0

Таблица 1

Изменение параметров потока по радиусу

----------------------------------------------------------

Паpаметp | Сечение по высоте лопатки

| 1(пеp) 2 3(сp) 4 5(вт)

----------------------------------------------------------

r .2390 .2292 .2195 .2097 .2000

ro 1.000 .9592 .9184 .8776 .8368

U 484.8 465.0 445.2 425.5 405.7

C1u 576.0 596.1 617.9 641.6 667.3

C1a 184.1 189.8 196.0 202.7 210.1

alf1 16.70 16.70 16.70 16.70 16.70

C1 604.7 625.6 648.3 672.8 699.6

be1 63.65 55.36 48.62 43.17 38.76

C2u -103.1 -111.7 -121.2 -131.8 -143.6

W2u 587.9 576.7 566.4 557.3 549.3

C2a 207.6 203.6 200.0 196.8 194.0

be2 19.45 19.45 19.45 19.45 19.45

Л1 .8927 .9236 .9570 .9933 1.033

Rт .4606 .4226 .3800 .3321 .2780

T2w 1393. 1387. 1382. 1377. 1372.

Л2w .9224 .9068 .8925 .8797 .8687

Л1w .3040 .3421 .3881 .4411 .5003

Л2 .3606 .3613 .3639 .3685 .3755

dbe 96.90 105.2 111.9 117.4 121.8

alf2 63.59 61.26 58.78 56.18 53.48

( Данные по пpофилиpованию записаны в файл gfrt.dat )

Рис. 2.4 График изменения и по высоте лопатки на входе и выходе из рабочего колеса

Рис. 2.5 График изменения Lc и Lw по высоте лопатки на входе и выходе из рабочего колеса.



Рис. 2.6 График изменения и по высоте лопатки на входе и выходе из рабочего колеса

лопатка двигатель сопло турбина

Рис. 2.7 Треугольники скоростей в 1,2 и 3 сечениях лопатки рабочего колеса



Рис. 2.8 Треугольники скоростей в 4 и 5 сечениях лопатки рабочего колеса

Рис 2.9 Решетка профилей на периферии

Рис 2.10 Решетка профилей на r=0.959



Рис 2.11 Решетка профилей на r=0.918

Рис 2.12 Решетка профилей на r=0.878

Рис 2.13 Решетка профилей на r=0.837

Вывод

В этой части курсового проекта были получены решетки профилей первой ступени рабочего колеса турбины высокого давления на трех радиусах: втулочном, среднем и периферийном. Результаты ручного счета и счета на ЭВМ незначительно отличаются. Это можно объяснить принятыми допущениями при ручном счете. Полученные профили лопаток имеют среднюю относительную толщину (=0,264). Это связано с тем, что лопатка охлаждаемая, так как работает при температурах (Т*г=1585 К). Полученные параметры удовлетворяют требованиям: W2>W1, углы >55 ( 57,4, 58,8, 68,46 градусов н втулочном, среднем и периферийном радиусах). Форма межлопаточного канала решетки профилей - конфузорная.

3. Проектирование камеры сгорания

Камеры сгорания газотурбинных двигателей и установок представляют собой наиболее сложный узел, в котором одновременно протекают различные процессы, течения, физико-химические процессы горения, тепловые процессы, связанные с тепловыми потоками и термическими нагрузками деталей. Большинство из этих процессов плохо поддаются расчетам, поэтому при создании КС требуется большой объем экспериментальных и доводочных работ. Проектировочный расчет является первым приближением в создании КС новых двигателях с одновременным использованием предыдущего опыта каждой конкретной двигателестроительной фирмы. Особое внимание при создании новых двигателей в последнее время уделяется образованию вредных веществ в КС, исходя и удовлетворения экологических нормам. Анализ аварийных ситуаций при эксплуатации авиадвигателей свидетельствует о случаях, причинами которых были колебательные процессы, возникшие в КС и связанные с ее акустическими нормами.

3.1 Формирование исходных данных

Рис 3.1 Геометрия камеры сгорания

1. Температура воздуха по заторможенным параметрам на выходе из компрессора (сечение К-К) =717,5 К.

2. Температура газа по заторможенным параметрам на выходе из КС

=1585 К.

3. Полное давление воздуха на выходе из компрессора =1805000 Па.

4.Коэффициент восстановления полного давления КС, обусловленный гидравлическими потерями =0,96.

5. Коэффициент восстановления полного давления КС, обусловленный тепловыми потерями =0,98.

6. Коэффициент полноты сгорания КС =0,995.

7. Коэффициент полноты сгорания в зоне горения =0,85.

8. Стехиометрическое количество воздуха для используемого топлива (керосин) = 14,8 .

9. Теплотворная способность топлива =43000 .

10. Коэффициент избытка воздуха КС =2,893.

11. Коэффициент избытка воздуха на выходе из фронтового устройства =0,5

12. Коэффициент избытка воздуха в конце зоны горения =1,6

13. Диаметр компрессора на выходе из НА =0,412 м.

14. Относительный внутренний диаметр компрессора на выходе .

15. Диаметр турбины на входе в СА =0,466 м.

16. Относительный внутренний диаметр турбины на входе в СА

 .

17.Относительный диаметр КС прототипа

.

18. Относительный внутренний диаметр КС прототипа

.

19. Относительный диаметр жаровой трубы

.

20. Относительный внутренний диаметр жаровой трубы

.

21. Относительный диаметр фронтового устройства

.

22. Относительный внутренний диаметр фронтового устройства

.

23. Относительная длина диффузора КС

.

24. Относительная длина жаровой трубы

.

25. Относительная длина головки жаровой трубы

.

26. Относительная длина зоны горения

.

27. Относительная длина газосборника жаровой трубы

.

3.2 Компьютерный расчет кольцевой камеры сгорания

Файл исходных данных GDR KS.DAT:

Таблица 3.1

1

23.1 717.5 1585.0 1805000 0.960 0.980 0.99 0.85

14.8 43e6 2.893 0.5 1.6

0.412 0.903 0.4660 0.858 1.36 0.466

0.917 0.584 0.5690 0.300 1.00

0.755 1.830 0.6100 0.351 0.35

1.000 1.000 1.000

Расчет выполняем с помощью программы GDR KS.EXE.

Файл полученных результатов GDR KS.REZ:

Таблица 3.2

Геометрический расчет кольцевой камеры сгорания

-------------------------------------------------------------------------

Исходные данные

Gв Tк* Tг* Pк* б гидр б тепл КПДг КПДзг

23.100 717.500 1585.000 1.805E+06 0.960 0.980 0.990 0.85

ALFA0 Hu ALFA ALFAф ALFAзг

14.8 4.300E+07 2.893 0.500 1.600

_ _ _ _

Dк dк Dт dт Dкс dкс

0.412 0.903 0.466 0.858 1.360 0.466

_ _ _ _

Dж dж Dф dф E

0.917 0.584 0.569 0.300 1.000

_ _ _ _ _

lд lж lг lгс lзг

0.755 1.830 0.610 0.351 0.350

_

Kн Kвн lц

1.000 1.000 1.000

-------------------------------------------------------------------------

Результаты расчета:

Gт Gф Gзг Gв см Gг

0.540 3.99 12.78 10.32 23.64

Tзг* Cp зг [RO]зг [RO]к Dкс ср

1790.3 1225.9 3.50 8.77 0.437

Dк Dк вн Dкс h кс Dкс вн

0.412 0.372 0.560 0.150 0.261

Dт Dт вн Dж h ж Dж вн

0.466 0.400 0.514 0.107 0.300

Fкс Fж Vж тр h в h н

0.1930 0.1366 0.0267 0.0214 0.0214

l кс l д l ж тр l г l зг

0.276 0.081 0.196 0.065 0.037

l з см l гс Z Dф Dф вн

0.158 0.038 12 0.061 0.018

Qv Cж тр Cсм

1714200 26.71 20.88

Fкс опт Fж опт

0.0747 0.0462

EJco EJcн EJno

27.2 3.7 13.2



Рис. 3.2 Эскиз камеры сгорания.

Вывод

В этой части курсового проекта был проведен расчет кольцевой камеры сгорания. Теплонапряженность камеры равна QV=1714200Дж/К2Па4. Умеренные значения теплонапряженности позволяют увеличить ресурс. Топливом служит авиационный керосин, низшая теплотворная способность которого Hu=43000 кДж/кг.

Расчет камеры сгорания проводился на основании геометрических размеров камеры сгорания двигателя-прототипа и в результате расчета было получено отклонение геометрии от оптимальных значений, в виду изменения параметров цикла двигателя.

4. Проектирование выходного устройства

Исходя из выбранного прототипа проектируем эжекторное сопло, внутренняя часть которого является реактивным насадком двигателя, а наружная - элементом летательного аппарата. Этот вид сопел отличается от сопел Лаваля тем, что у них сверхзвуковой контур полностью или частично заменен границей свободной сверхзвуковой струи, а для улучшения характеристик осуществляется подача в сопло вторичного (эжектируемого, пассивного) воздуха. Сопла этого вида называются также соплами с разрывом сверхзвукового контура.

Сопло с разрывом сверхзвукового контура проще регулируется, чем сопло Лаваля. Это объясняется тем, что площадь критического сечения первичного сопла может изменятся за счет использования обычных регулируемых створок. Дополнительная возможность регулирования обеспечивается подачей вторичного воздуха, благодаря чему сама граница свободной струи может в той или иной степени изменять свою форму, приспосабливаясь к изменению положения регулируемых створок и перепада давлений в сопле. Подача вторичного воздуха способствует при этом снижения внутренних потерь в сопле, а также обеспечивает интенсивное охлаждение элементов конструкций.

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

4.1 Газодинамический расчет сопла

4.1.1 Формирование исходных данных

Исходными данными для расчета эжекторного сопла являются известные параметры активного и пассивного потоков:

­ Расход активного потока ;

­ Доля пассивного потока ;

­ Давление в сечении А-А ;

­ Давление в сечении П-П ;

­ Температура в сечении А-А ;

­ Температура в сечении П-П ;

­ Константа для газа ;

­ Константа для воздуха ;

­ Показатель изоэнтропы ;

­ Универсальная газовая постоянная ;

­ Показатель изоэнтропы ;

­ Универсальная газовая постоянная ;

­ Температура воздуха ;

­ Расход воздуха ;

­ Расход газа ;

­ Давление окружающей среды .

Где некоторые величины были рассчитаны следующим образом:

1.

2. ;

3.

4.1.2 Компьютерный расчет эжекторного сопла

Таблица 4.1

Газодинамический расчет ЭЖЕКТОРНОГО СОПЛА

-------------------------------------------

Исходные Данные

Gа = 48.49 кг/с

y = 0.05

Pа° = 226094 Па

Pп° = 98306 Па

Tа° = 2000 К

Tп° = 288.15 К

mг = 0.03918

mв = 0.0404

Kг = 1.30

Rг = 290.0 Дж/кг·К

Kв = 1.4

Rв = 287.0 Дж/кг·К

Tв° = 288.15 К

Gв = 48.49 кг/с

Gг = 48.49 кг/с

Pн = 101325 Па

Результаты Расчета

Fа = 0.2448 мм

rа = 0.2791 м

Pа` = 51935 Па

Gп = 2.425 кг/с

Fп` = 0.0104 мм

п(Ла`)= 0.230

Ла` = 1.4855

q(Ла`)= 0.763

Fа` = 0.3206 мм

F` = 0.3310 мм

rа` = 0.3195 м

r` = 0.3246 м

Cpг = 1257 Дж/кг·К

Cpв = 1005 Дж/кг·К

Tс° = 1240 К

Pc° = 220009 Па

п(Лс) = 0.461

Лс = 1.1207

q(Лс) = 0.984

Fc = 0.4027 мм

rc = 0.3580 м

Рис. 4.1 Эскиз выходного устройства

4.2 Профилирование сопла

4.2.1 Формирование исходных данных

Исходные данные для профилирования сопла выбираются исходя из его газодинамического расчета, а при использовании программы soplo.exe, они вбиваются в раздел «профилирование», автоматически после газодинамического просчета сопла.

­ Радиус в сечении А'- А' ;

­ Радиус в сечении А- А ;

­ Угол ;

­ Радиус в сечении С- С ;

­ Радиус ;

­ Угол ;

­ Радиус скругления ;

­ Угол ;

­ Угол .

4.2.2 Компьютерное профилирование эжекторного сопла

Таблица 4.2

Профилирование ЭЖЕКТОРНОГО СОПЛА

------------- Исходные Данные ----------------

ra` = 0.3195 м

ra = 0.2791 м

J` = 14.0 °

rc = 0.3580 м

r` = 0.3246 м

Jc = 9.0 °

R` = 0.3246 м

Jg = 30.0 °

Jн = 4.0 °

----------- Результаты Расчета -------------

Lg = 0.1620 м (длина дозвуковой части сопла)

Lc` = 0.2109 м (длина сверхзвуковой части сопла)

Lc = 0.2364 м (уточненная длина)

Рис. 4.2 Чертеж выходного устройства

Вывод

В этой части курсового проекта проведен расчет эжекторного сопла. Для внутреннего части: площадь сопла равна 0,2448 м2, длина 0,162 м; для наружного части: площадь сопла равна 0,4027 м2, длина 0,2109 м. Следует также отметить указанные величины соответствуют режиму работы «полный форсаж» при

Заключение

В данном курсовом проекте было выполнено профилирование рабочего колеса ступени компрессора высокого давления и газовой турбины, расчёт кольцевой камеры сгорания, расчёт эжекторного сопла на базе данных, полученных в результате термогазодинамического просчета двигателя: к*=18 , Тг*=1585К, G=50,28 кг/с.

На этапе расчета и профилирования первой ступени КВД были определены компоненты треугольников скоростей, а также геометрические параметры решеток профилей на трех радиусах. Полученные дозвуковые (Мw1=0,7263) профили. При профилировании применялся закон крутки, «твердое тело» (на входе) при заданном Нт(r).

На этапе профилирования ступени газовой турбины были определены компоненты треугольников скоростей в межвенцовых зазорах, а также геометрические параметры решёток профилей на трёх радиусах. На втулке газовый поток дозвуковой (л1=0,8927). Применение законов 1=const и 2=const значительно упрощает технологию изготовления лопаток СА и РК, а также позволяет создать хорошую конструктивную базу для их монтажа в статоре и роторе.

В разделе “Расчёт камеры сгорания” была просчитана кольцевая камера сгорания прототипа АИ-222К-25Ф, но с другими параметрами цикла. Теплонапряженность камеры получилась нормальная (QV=1714200Дж/К2Па4.), что сокращает ресурс двигателя в целом, но скорости потоков в камере получились малы.

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

Список использованной литературы

1. Незым В.Ю. Расчет и построение решеток профилей дозвукового осевого компрессора. Харьков, ХАИ 1988г.

2. Коваль В.А. Профилирование лопаток авиационных турбин. Харьков, ХАИ 1986 г.

3. А. Н. Анютин Расчет и профилирование на ЭВМ лопаток осевой газовой турбины. Харьков, ХАИ 1991 г.

4. Герасименко В.П., Никишов А.А. Проектирование камер сгорания газотурбинных двигателей. Харьков, ХАИ 1999 г.

5. Анютин А. Н., Дегтярев О. Д., Ковалев В. И. , Незым В. Ю.

Профилирование выходных устройств ГТУ. Харьков, ХАИ 2001 г.

Размещено на Allbest.ru