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

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

.pdf
Скачиваний:
8
Добавлен:
12.11.2023
Размер:
510.95 Кб
Скачать

Затем по рис. 3 приложения 1 пособия [1] с графиков можно снять для принятой температуры газа искомое значение коэффициента избытка воздуха. В дальнейшем это значение используется для определения (для смеси продуктов сгорания топлива и избыточного воздуха) показателя адиабаты kг и изобарной теплоёмкости cpг (рис. 6 и 7 прил. 1 или табл. 1 и 2 прил. 2 [1]

для авиационного керосина как топлива; табл. 3 и 4 прил. 2 [1] для природного газа в качестве топлива) по предварительно найденной средней температуре газа в турбине:

T *

T * 800 2.

(2)

г ср

Г

 

Показатель адиабаты kв и изобарная теплоёмкость c p в воздуха определяются по средней температуре воздуха в компрессоре:

T *

T * T 2 .

(3)

к ср

К н

 

2. Свободная энергия в двигателе. Поскольку рассматривается расчёт наземного ТВаД, постольку расчёт свободной удельной работы цикла несколько упрощается относительно реактивных ГТД [1]:

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

cp вTн

 

 

* kв 1

 

 

 

 

 

 

 

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

L

c T

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kв

1

 

св

 

p г

Г

 

 

 

 

 

 

 

 

 

kг 1

 

 

* *

 

 

к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

*

 

*

 

*

 

kг

 

к т

 

 

 

 

 

 

 

 

 

 

 

 

 

ВЗ

 

КС

ПК

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Удельную эффективную мощность находят по формуле

 

 

2

 

 

 

 

cС

 

*

Ne уд

2

ТС тр.

Lсв

 

 

 

2 ВУ

 

(4)

(5)

Здесь коэффициент скорости ВУ выходного устройства использован для

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

4. Удельный эффективный расход топлива определяется в пересчёте на единицу мощности, измеряемой в кДж/кг:

 

 

3600

 

ce

 

 

 

 

.

(6)

L

N

e уд

10 3

0

 

 

 

 

Результаты расчёта рекомендуется свести в таблицу, состоящую из представленных ниже блоков. Для удобства анализа полученных результатов и выбора расчётной точки необходимо построить графики, аналогичные показанным на рис. 7.1 пособия [1].

11

 

 

Результаты предварительного расчёта

 

 

 

 

 

 

 

 

 

 

 

 

T * , К

cС , м/с

*

α

Lсв, Дж/кг

Neуд , Дж/кг

ce,

кг

 

 

 

Г

 

к

 

 

 

 

кВт ч

 

 

 

 

 

 

 

 

 

 

190

6,8

 

 

 

 

 

 

1250

 

 

 

 

 

 

 

 

 

 

 

 

 

6,8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Однако повышенная удельная мощность может быть получена не только за счёт уменьшения скорости истечения cС , но и температуры газа TГ*.

Рост последней нежелателен вследствие снижения ресурса и необходимости применять мероприятия по его восстановлению.

Для повышения экономичности следует выбирать как можно большие значения степени повышения давления *к. Но в этом случае может умень-

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

1.4.Расчёт термогазодинамических параметров

вхарактерных сечениях

1.4.1.Сечение «Вх» перед компрессором

Так как наземные установки работают в стендовом режиме (скорость полёта Mн = 0), то температуру торможения необходимо принять равной атмосферной для условий МСА:

T *

T

288 К.

(7)

Вх

н

 

 

Давление торможения определяется по атмосферному давлению с учётом потерь в воздухозаборнике (входном устройстве):

pВх* pн *ВЗ 101325

*ВЗ .

(8)

Статическая температура воздуха:

 

 

 

 

 

T *

c2

 

 

T

Вх

,

 

(9)

 

 

Вх

Вх

2cp в

 

 

 

 

 

 

12

где изобарную теплоёмкость воздуха cpв принимать равной 1005 Дж/(кг К)

для температуры торможения на входе. Более точное значение, найденное для средней температуры TВх* ср TВх* TВх 2 , не требуется по причине не-

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

адиабаты:

kв

pВх pВх* TTВх* kв 1 , (10)Вх

где показатель адиабаты воздуха kв принимать равным 1,4 для тех же условий, что и для изобарной теплоёмкости, используемой в расчёте статической температуры (9).

Плотность воздуха по уравнению состояния:

Вх

pВх

,

(11)

 

 

RвTВх

 

где принять газовую постоянную Rв = 287 Дж/(кг К).

1.4.2. Сечение «1» перед первым рабочим колесом компрессора

Температура торможения в ВНА не изменяется: T1* TВх* . Давление торможения уменьшается на величину потерь в ВНА:

p* p*

*

,

(12)

1

Вх

НА

 

 

где коэффициент сохранения полного давления в ВНА может принимать зна-

чения 0,98…0,985.

 

T T *

 

 

c2

 

Статическая температура:

 

 

1

 

.

 

(13)

 

 

 

1

Вх

 

 

2cp в

 

 

 

 

 

 

 

Статическое давление:

p1 p1* T1

TВх*

kв

 

kв 1

.

(14)

Изобарную теплоёмкость воздуха cpв и показатель адиабаты воздуха kв принимать такими же, как и в сечении «Вх».

Статическая плотность:

 

 

p1

.

(15)

 

 

1

 

RвT1

 

 

 

 

 

13

1.4.3. Сечение «К» за компрессором

Потребная удельная работа компрессора:

*

 

*

kв 1

 

 

*

(16)

kв

 

 

Lк сp вTВх

к

 

 

1

к .

 

 

 

 

 

 

 

 

Здесь в первом приближении можно принимать значения изобарной теплоёмкости и показателя адиабаты воздуха такими же, как и для сечения «Вх». После нахождения температуры торможения за компрессором их необходи-

мо уточнить для средней температуры Tк*ср TК* Tн 2 по рис. 6 и 7 приложения 1 пособия [1].

Температура торможения за компрессором:

 

 

 

T * T *

 

Lк

.

 

(17)

 

 

 

 

К

Вх

 

cp в

 

 

 

 

 

 

 

 

Затем по температуре T *

уточнить значения kв и c

pв

и повторить вычисле-

к ср

 

 

 

 

 

 

ния (16) и (17). Обычно достаточно одного-двух уточнений.

Давление торможения на выходе из компрессора:

 

pК* pВх* *к .

 

(18)

Статические параметры состояния воздуха за компрессором TК , pК

и К рассчитываются по

формулам,

аналогичным

(9)…(11). Разумеется,

необходимо использовать индексы, соответствующие сечению «К». Величины cp в , kв находить по указанным выше приложениям [1] для температуры за

компрессором TК* . Газовую постоянную для расчёта плотности искать из термодинамического соотношения

R

kв 1

c

p в

.

(19)

 

в

kв

 

 

 

 

 

 

1.4.4. Сечение «Г» за камерой сгорания

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

Для начала расчёта относительного расхода топлива принять теплоёмкость газа по средней температуре в камере сгорания Tг*ср TК* TГ* 2 для

14

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

Относительный расход топлива, приходящийся на 1 кг воздуха:

 

 

 

с

p г

T *

T *

 

 

 

 

 

 

 

Г

К

 

qт

H

u

 

г

с

p г

T * T * .

(20)

 

 

 

 

 

Г К

 

Уточнение коэффициента избытка воздуха:

 

1

 

,

(21)

q

L

 

 

 

 

т

 

0

 

 

где теоретически необходимое количество воздуха L0 для сжигания 1 кг топлива принимать согласно используемому виду топлива (авиационный керосин или природный газ).

По уточнённому значению α и температуре Tг*ср скорректировать изо-

барную теплоёмкость газа в камере сгорания и повторить расчёт по (20) и (21). Повторение прекратить при достижении точности 0,1…0,2% (до 4-й значащей цифры) по относительному расходу топлива.

Давление торможения газа за камерой сгорания:

pГ* pК* *КС .

(22)

Для температуры TГ* и уточнённого коэффициента α по данным [1] определить изобарную теплоёмкость сpг и показатель адиабаты kг газа.

Вычислить газовую постоянную

R

kг 1

c

p г

.

(23)

 

г

kг

 

 

 

 

 

 

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

Статические параметры состояния газа в сечении за камерой сгорания TГ, pГ и Г рассчитываются по формулам, аналогичным (9)…(11). Исполь-

зовать индексы, соответствующие состоянию газа в сечении «Г». Теплоёмкость, показатель адиабаты и газовую постоянную можно принимать по температуре торможения, то есть без коррекции на температуру, среднюю из полной и статической температур.

15

1.4.5. Сечение «ТК» за турбиной компрессора

 

 

 

 

 

 

 

 

Определение параметров

 

 

 

 

 

 

 

за первым сопловым аппаратом

 

 

 

 

 

Удельную работу турбины компрессора предлагается рассчитывать по

балансу мощностей компрессора и турбины. С учётом того, что в современ-

ных двигателях параметры цикла велики и требуют, как правило, интенсив-

ного охлаждения первого соплового аппарата (СА1), в баланс мощностей

необходимо ввести возврат из СА1 охлаждающего воздуха перед первым ра-

бочим колесом:

 

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

отб

 

охлСА1

 

 

 

GввхLк Gг Gотб GохлСА1 LТК м Gг 1

G

 

 

G

LТК м .

 

 

 

 

 

г

 

 

 

г

 

 

 

 

Используя известные соотношения Gг 1 qт Gввх

и

1 Gотб

Gг , вводя

поправку на возврат в проточную часть охлаждающего первый СА воздуха

СА1 GохлСА1 Gг ,

окончательно

формулу

для

определения

потребной

удельной работы турбины компрессора можно записать так:

 

 

 

 

 

Lк

 

 

 

Lк

 

 

 

 

 

 

 

LТК 1 qт СА1 м

1 qт ТК м .

 

 

 

(24)

 

 

 

Для определения коэффициента

 

 

возврата воздуха СА1 в первом

 

 

сопловом аппарате необходимо за-

 

 

дать

рабочую

 

температуру

лопаток

 

 

СА1 (в диапазоне от 1000 до 1150 К в

 

 

зависимости от требуемого ресурса).

 

 

Затем найти потребную эффектив-

 

 

ность охлаждения:

 

 

 

 

 

 

 

 

 

 

 

 

 

T *

T

 

 

(25)

 

 

 

 

 

 

Г

 

лСА1 .

 

 

 

 

 

 

 

 

T

*

T

*

 

 

 

 

 

 

 

 

 

 

Г

К

 

 

 

 

С графиков на рис. 2 снять значение

Рис. 2. Эффективность различных схем

доли охлаждающего воздуха, необхо-

димой

для достижения

этой

эффек-

воздушного охлаждения: 1 – конвектив-

тивности.

 

 

 

 

 

 

 

 

 

ное охлаждение с продольным течени-

 

Поскольку при больших темпе-

ем охлаждающего воздуха; 2 – то же

 

ратурах газа за камерой сгорания (бо-

с поперечным течением охладителя;

3 – конвективно-пленочное (комбини-

лее

чем

примерно

1380…1400 К)

рованное) охлаждение; 4 – пористое

необходимо

 

охлаждать

и

полки

(эффузионное) охлаждение

сопловой

лопатки,

 

то

полученное

 

 

значение

необходимо умножить на

 

 

поправочный коэффициент. Тогда

16

СА1 kп

 

 

(26)

Gв .

Для температуры газа 1450 К и выше можно принимать kп 3 . Для не-

охлаждаемых полок СА1 при сниженных температурах газа поправочный коэффициент принимать равным единице.

Для дальнейших расчётов потребуется найти температуру за первым сопловым аппаратом. Приближённо это можно вычислить так:

 

 

с

T * с

p в

T *

СА1

 

 

T *

 

 

p г Г

К

.

(27)

 

 

 

 

 

СА

 

 

сp г сp в

СА1

 

 

 

 

 

Здесь изобарную теплоёмкость газа принимать согласно данным [1] по температуре TГ* и уточнённому в расчёте сечения «Г» коэффициенту избытка воздуха в камере сгорания, изобарную теплоёмкость воздуха находить по температуре за компрессором TК* . Обе теплоёмкости определяются без ите-

рационных уточнений по средней температуре в компрессоре и турбине. Полное давление за первым СА необходимо определять с учётом по-

терь от выдува охлаждающего воздуха (потерь смешения) по приближённой формуле:

p

*

p

*

0,985

 

 

СА1

 

,

(28)

 

 

 

 

 

СА

 

Г

 

 

3

 

 

 

где 0,985 – среднее значение коэффициента сохранения полного давления в неохлаждаемом сопловом аппарате, а вычитаемое моделирует потери смешения, пропорциональные количеству возвращаемого в проточную часть охлаждающего первый СА воздуха. При необходимости или желании полное давление за СА можно рассчитать по методике, аналогичной расчёту камеры смешения двухконтурного турбореактивного двигателя со смешением потоков (ТРДДсм).

Определение параметров за турбиной компрессора

Приняв предварительно изобарную теплоёмкость и показатель адиабаты такими же, как и для расчёта температуры TСА* , найти в первом приближении температуру за турбиной компрессора:

T *

T *

 

LТК

.

(29)

 

ТК

СА

 

cp г

 

 

 

 

 

По средней температуре TТК* ср TСА* TТК* 2 уточнить теплоёмкость и пока-

затель адиабаты. Пересчитать температуру за турбиной компрессора. Уточнение можно прекратить при достижении точности нахождения TТК* не более 1…2 К (примерно 0,1 %).

17

Степень понижения полного давления в турбине компрессора:

*

 

 

 

 

 

1

 

 

 

.

(30)

 

 

 

 

 

 

 

 

ТК

 

 

 

 

 

 

 

 

kг

 

 

 

 

 

 

LТК

kг 1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

c

 

T *

*

 

 

 

 

 

 

 

 

p г

 

 

 

 

 

 

 

 

 

 

СА

ТК

 

 

 

 

Полное давление за турбиной компрессора:

pТК* pСА*

*ТК .

(31)

Для расчёта статических параметров задать число Маха MТК. Критерием адекватности принятого значения MТК можно считать приемлемую форму проточной части турбины компрессора (определяется при расчёте диаметральных размеров).

Статическая температура за турбиной компрессора:

T

T

*

 

 

k

г

1

M

2

 

1

(32)

 

1

 

2

 

 

.

ТК

ТК

 

 

 

 

ТК

 

 

Здесь kг определять для температуры T * .

Для этой же температуры найти

 

 

 

 

 

 

 

ТК

 

 

 

 

 

и теплоёмкость cp г .

Статическое давление найти по уравнению изоэнтропы:

 

pТК pТК*

TТК TТК*

 

kг

 

 

 

.

(33)

 

kг 1

Газовая постоянная:

R kг

1c

p г

.

 

 

(34)

 

г

kг

 

 

 

 

 

 

 

 

 

 

 

Плотность газа:

ТК pТК

RгTТК .

(35)

Скорость потока:

cТК MТК

kгRгTТК .

(36)

1.4.6. Сечение «ТК′» («Д»)

Осевую скорость потока в данном сечении принимать с учётом наличия или отсутствия переходного канала. При наличии такого канала скорость назначается независимо от скорости за турбиной компрессора. При отсутствии – равной скорости за турбиной компрессора cТК cТК .

В данном сечении температуру торможения принимать равной температуре торможения за турбиной компрессора TТ*К TТК* . При наличии пере-

ходного диффузорного канала давление торможения рассчитывать с учётом потерь в данном канале:

18

pТ*

К pТК* *ПК ,

(37)

иначе pТ*К pТК* . Коэффициент

сохранения

полного давления принимать

равным *ПК 0,985 0,995 . Чем больше снижение скорости в переходном

канале, тем меньшее значение *ПК надо выбирать. Статическая температура перед силовой турбиной:

 

T *

 

c2

 

T

 

ТК

.

(38)

 

ТК

ТК

 

2c p г

 

Здесь теплоёмкость cpг и показатель kг адиабаты принимать такими же, как в сечении «ТК». Скорость газа (108) предварительно определяется при нахождении размеров данного сечения исходя из потребной площади поперечного сечения «ТК′».

Статическое давление найти по уравнению изоэнтропы:

 

 

p*

T

 

 

T *

 

kг

 

p

ТК

 

 

kг 1

.

(39)

 

ТК

ТК

 

ТК

 

 

 

 

Плотность газа:

ТК pТК

RгTТК .

(40)

1.4.7. Сечения «Т» и «С»

Определение потребного полного давления за силовой турбиной

Предварительное значение удельной работы силовой турбины:

 

2

 

 

 

 

cС

 

 

*

 

 

ТС .

LТС Lсв

2

2

 

С

 

Предварительная температура торможения за силовой турбиной:

TТ* TТ*К LТС .

cp г

(41)

(42)

Здесь предварительно принимать теплоёмкость такой же, как и в сечении «ТК». В дальнейшем при уточнениях переопределять по средней в турбине

температуре торможения TТС* ср TТ*К TТ* 2 .

Критическая скорость определяется температурой торможения за силовой турбиной, так как TС* TТ* :

a

крС

a

крТ

2

kг

R T * .

(43)

 

 

 

 

kг 1

г Т

 

 

 

 

 

 

 

 

19

Принятые для критической скорости показатель адиабаты и газовая постоянная требуют уточнения после нахождения более точного значения TТ* .

Статическое давление за ВУ pс = 101325 Па. Полное давление на выходе из ВУ:

*

 

 

 

 

pС

 

 

 

 

,

(44)

pС

 

 

 

 

 

 

 

 

 

 

k

 

1

 

 

kг

 

 

г

2

kг 1

 

 

 

 

 

 

 

 

 

 

 

 

1

kг

1

С

 

 

 

 

 

 

 

где С cСaкрС . Скорость cс истечения из выходного устройства для ТВаД

наземного применения можно принимать не на выходе из него (скорость выхлопа), а на выходе из силовой турбины cС cТ . Это необходимо для учёта

потерь в выхлопной шахте (скорость за турбиной больше скорости выхлопа в атмосферу) и для согласования методики расчёта ТВаД с методикой расчёта ТВД.

Потребное полное давление за силовой турбиной

pТ* pС*

*ВУ ,

(45)

где *ВУ *РС – коэффициент сохранения полного давления в выходном устройстве (реактивном сопле).

Расчёт параметров силовой турбины

Располагаемая степень понижения полного давления:

*

p*

 

p* .

(46)

ТС

ТК

 

Т

 

Удельная работа силовой турбины:

 

 

1

 

 

 

*

 

 

*

 

LТС cp гTТК 1

 

 

 

 

ТС .

(47)

*

kг 1

 

 

kг

 

 

 

 

 

ТС

 

 

 

Далее необходимо повторить расчёт, начиная с уточнения температуры (42) за турбиной. Обычно достаточно одного-двух уточнений.

Для расчёта статических параметров задать число Маха MТ. Критерием адекватности принятого значения MТ можно считать приемлемую форму проточной части силовой турбины (определяется при расчёте диаметральных размеров).

Статическая температура за силовой турбиной:

 

 

k

г

1

M2

 

1

(48)

T T * 1

 

2

 

.

Т Т

 

 

 

Т

 

 

20

Соседние файлы в папке книги