Главная Случайная страница



Полезное:

Как сделать разговор полезным и приятным Как сделать объемную звезду своими руками Как сделать то, что делать не хочется? Как сделать погремушку Как сделать неотразимый комплимент Как сделать так чтобы женщины сами знакомились с вами Как сделать идею коммерческой Как сделать хорошую растяжку ног? Как сделать наш разум здоровым? Как сделать, чтобы люди обманывали меньше Вопрос 4. Как сделать так, чтобы вас уважали и ценили? Как сделать лучше себе и другим людям Как сделать свидание интересным?


Категории:

АрхитектураАстрономияБиологияГеографияГеологияИнформатикаИскусствоИсторияКулинарияКультураМаркетингМатематикаМедицинаМенеджментОхрана трудаПравоПроизводствоПсихологияРелигияСоциологияСпортТехникаФизикаФилософияХимияЭкологияЭкономикаЭлектроника







Расчет теоретического значения тяги без учета входного импульса





PjR = GN · GCN - FCN · PH ·106 , (2.6.10.49)

действительной тяги, но без учета входного импульса

Р/OCT = РOCT - Pj · PjR - FCN · PH ·106. (2.6.10.50)

расчет относительного коэффициента тяги осуществляется по формуле

 

(2.6.10.51)

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

PjR = GZ · GCZ - FCZ · PC ·106 , (2.6.10.52)

РOCT = Р/OCT + Рj · PjR - FCZ · PH ·106. (2.6.10.53)

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

Р/OCT = РOCT - RBX , (2.6.10.54)

после чего

Р/OCT = РOCT .

 

 

2.6.11. Модуль узла типа теплообменник или воздухо-воздушный радиатор

Подпрограмма TEPLOO

 

Алгоритм модуля обеспечивает расчет теплообменников регенеративного и рекуперативного типов, а также теплообменников типа воздухо-воздушного радиатора, как на установившихся, так и на переходных режимах.

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

Алгоритм модуля позволяет проводить расчеты как с использованием характеристик (см. раздел 2.2), так и без характеристик. При этом должны быть заданы степень регенерации m, отношение водяных эквивалентов R, коэффициенты восстановления полного давления с “холодной” и “горячей” сторон теплообменника sX, sГ.

Характеристики задаются только в аппроксимированном виде. Коэффициенты полиномов получаются с помощью специальной программы “APRTO”, включенной в сервисное обслуживание программного комплекса ГРАД.

В алгоритме предусмотрена возможность учета утечек из “холодной” части и перетечек из “холодной” части в “горячую” часть теплообменника. Алгоритм модуля

позволяет включать в схему ГТД его отдельные части “холодную” или “горячую” с соблюдением теплового баланса.

Входные данные модуля узла типа теплообменник делятся на три части.

1. Текущие значения параметров воздуха (газа) на входе в “холодную” (“горячую”) части теплообменника. В зависимости от номера контура, в котором находится “холодная” или “горячая” часть (эти номера записаны в разрядах "В" и "С" условного номера узла), считываются из соответствующего канала передачи информации типа "Контур", т.е. из групп “21000”...”25000” массива “А”. В число этих параметров входят:



- температура и давление торможения,

G1 - расход воздуха (газа),

qТ1 - относительный расход топлива ,

- энтальпия и энтропия торможения воздуха (газа).

Кроме того, из группы “26000” при расчете переходных процессов считываются шаг интегрирования Dt, значения температур торможения на входе и выходе в “холодную” часть на предыдущем шаге интегрирования , .

2. Параметры узла. Передающиеся в подпрограмму через массив “А” в составе информационной подгруппы модуля “9ВС100”.

Характеристики узла передаются в подпрограмму через массив “А” в составе информационной подгруппы “9ВС300”.

 

Описание алгоритма

 

При наличии характеристик рассчитывается значение по полиному

(2.6.11.1)

и определяется

.

 

Рассчитывается температура на выходе из “холодной” части теплообменника по формуле

 

Т/ = mТ* + (1-m) Т* (2.6.11.2)

 

рассчитывается прибавка температуры DТХ2 по формуле

 

(2.6.11.3)

где Dti - шаг интегрирования по времени, с ;

DT1Xi* - изменение температуры на входе в “холодную” часть за время Dti ;

rt0 - постоянная времени теплообменника

 

Приближенно значение rt0 можно вычислить по формуле

(2.6.11.4)

 

 

Рассчитывается значение приведенного значения расхода

(2.6.11.5)

 

Рассчитывается относительный приведенный расход

и рассчитывается по полиному коэффициент восстановления давления

(2.6.11.6)

Определяется давление и расход на выходе из “холодной” части теплообменника

 

G2X = G1X(1 - qГ - q), (2.6.11.7)

 

Р2X* = Р1X* ·s. (2.6.11.8)

 

Рассчитываются тепловой поток, средний температурный перепад, площадь теплообменной поверхности при заданном коэффициенте температуропроводности

 

QX = GX2(H1X* - H2X*) ·1000, (2.6.11.9)

 

CP = (DТГ + DТX) / 2 (2.6.11.10)

 

FX = QX / (кDТCP), (2.6.11.11)

 

где DТCP - средний температурный напор

QX - тепловой поток с “холодной” стороны

FX - площадь поверхности теплообменника с “холодной” стороны с заданным коэффициентом температуропроводности и тепловым потоком

 

Рассчитывается смешение газа и воздуха (при наличии утечек) по формуле

(2.6.11.13)

Затем уточняются по НСН* значения Т* и значение qТ

qТ = GТ / GВ .

 

рассчитываются теплоемкости СРХ, СРГ и параметр R = (СРХ ·GХ) / (СРГ ·GГ).

Рассчитывается относительное значение степени регенерации

и рассчитывается параметр

(2.6.11.14)

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



Т* = mТ* ·R + (1- mR) T1Г* (2.6.11.15)

 

 

Рассчитывается приведенное значение расхода с “горячей” стороны

(2.6.11.16)

и уточняется разность расходов

DGГ = GГ - GХ (2.6.11.17)

Рассчитываются

, .

Рассчитывается коэффициент восстановления давления с “горячей” стороны теплообменника

(2.6.11.18)

 

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

(2.6.11.19)

 

 

Рассчитывается тепловой поток от “горячей” стороны в кВт

QГ = GГ* - Н*) (2.6.11.20)

 

2.6.12. Модуль узла типа Агрегаты

Подпрограмма AGREGT

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

Предусмотрена возможность задания сразу же нескольких характеристик. Их подключение в процессе расчета осуществляется программой управления при занесении соответствующего признака в элемент “10ВС104” подгруппы “10ВС100”.

Также предусмотрена возможность учета КПД агрегатов, которое задается в подгруппе “10ВС100” (элемент “10ВС106”). При задании в элементе “10ВС106” нуля, можно использовать характеристики в виде

МКР = f(x) .

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

Входные данные модуля узла типа Агрегаты на три части.

1. Текущее значение частоты вращения вала, на котором установлен узел n ,об/мин., и суммарная мощность на валу от которого приводится данный узел, SNe ,кВт. Данные параметры передаются через группу “1000” и группу “2000”

2. Параметры узла, передающиеся в подпрограмму через массив “А” в составе информационной подгруппы модуля “10ВС100”.

3. Характеристики узла передаются в подпрограмму через массив “А” в составе информационной подгруппы “10ВС300”.

Выходные данные делятся на две части.

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

NПОТ - потребная мощность для привода агрегатов , которая суммируется в соответствующем условном адресе группы “2000”.

2. Результаты расчета данного узла переписываются в подгруппу результатов модуля “10ВС200”.

 

Описание алгоритма

 

рассчитывается мощность агрегата

(2.6.12.1)

Если характеристика задана в координатах то определяется .

Затем вычисляется абсолютное значение мощности

(2.6.12.2)

 

Отбор мощности принимается постоянным

 

NeАГР = NeН.

 

Рассчитывается потребная мощность на агрегаты с учетом КПД.

NeПОТ = NeАГР / hРЕД, кВт (2.6.12.3)

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

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

, кН·м (2.6.12.4)

 

 

2.6.13. Модуль узла типа Воздушный винт

Подпрограмма VINT

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

Входные данные для модуля типа Воздушный винт делятся на четыре части.

1. Текущие значения параметров воздуха на входе в воздушный винт с учетом скорости и высоты полета. В число этих параметров входят:

ТH* - температура торможения;

РH* - давление торможения;

МП - число Маха полета;

NeВВ - мощность на выводном валу;

nВВ - частота вращения выводного вала. При заданном значении iРЕД считывается из группы "1000" и рассчитывается с учетом iРЕД. При i = 0считывается из группы "2000", куда значение nВВ заносит модуль "Редуктор".

2. Параметры узла.

Они передаются в подпрограмму через массив "А" в составе информационной подгруппы модуля "10ВС100".

3. Характеристики узла.

Передается в подпрограмму через массив "А" в составе информационной подгруппы "10ВС300".

Выходные данные делятся на две части:

1. Значение текущего параметра, который передается в группу "2000".

NeВ - мощность воздушного винта, которая суммируется в соответствующий элемент группы "2000" с учетом знака.

2. Результаты расчета данного узла переписываются в подгруппу результатов "10ВС200".

 

Описание алгоритма

рассчитываются

РН = РН* / (1 + 0,2Мп2)3,5 , (2.6.13.1)

ТН = ТН* / (1 + 0,2Мп2) . (2.6.13.2)

 

Рассчитывается скорость полета

VП = МП аН , м/с (2.6.13.3)

плотность воздуха

, кг/м3 (2.6.13.4)

где m - молярная масса воздуха на уровне моря, 28,96442 кг/кмоль;

R* - универсальная газовая постоянная, 8314,32 Дж/кмоль;

динамическая вязкость воздуха

,кг·с/м2 (2.6.13.5)

где br =1,458*10-6 кг/с·м·к 0,5 и S = 110,4 К - эмпирические коэффициенты Сатерленда

кинематическая вязкость воздуха

, м2/с (2.6.13.6)

 

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

NeВ = NeВВ .

Рассчитывается коэффициент мощности

, (2.6.13.7)

где nS - частота вращения винта в об/с,

DВ - диаметр винта в м.

 

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

hВМГ = hВ · КRe · КМ · КГ (2.6.13.8)

и полезная мощность NeПОЛ = NeВ · hВМГ .

Если скорость полета VП ¹ 0, то задается через входные данные узла коэффициент качества винта (1...1,5) и рассчитывается тяга винта

РВ = КМ ·NeВ, кН (2.6.13.9)

и коэффициент мощности

.

Если VП = 0, то рассчитывается полезная тяга

 

РПОЛ = NeПОЛ / VН, кН (2.6.13.10)

 

и качество винта

 

RВ = РПОЛ / NeВ, кН/кВт (2.6.13.11)

 

 

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

. (2.6.13.12)

при заданных значениях jВ и lВ определяются b и a.

Рассчитывается мощность винта

, кВт (2.6.13.13)

тяга винта

, кН (2.6.13.14)

кпд винта

. (2.6.13.15)

 

С учетом hВМГ рассчитывается полезная мощность

 

NeПОЛ = NeВ · hВМГ , (2.6.13.16)

 

полезная тяга

 

РПОЛ = NeПОЛ / VП , кН (2.6.13.17)

 

и качество винта

 

RВ = РВ / NeВ , кН/кВт . (2.6.13.18)

 

Расчет несущего винта. Определяется ометаемая винтом площадь

, м2 (2.6.13.19)

и радиус винта

 

RВ = DВ / 2 , м . (2.6.13.20)

 

Рассчитывается коэффициент заполнения

 

КZ = (ZЛ · FЛ) / FВ , где (2.6.13.21)

 

ZЛ - число лопастей,

FЛ - площадь лопасти, м2 ,

FВ - ометаемая винтами площадь, м2.

 

Если КZ £ 0, то КZ = 0,7854 · DВ2 .

Все конструктивные размеры сводятся к параметру

 

АК = FВ · КZ · R3 · 2 · 29676 ·10-3 . (2.6.13.22)

 

Обращаемся к подпрограмме "PSP" и определяем

 

МКР = f (jВ).

 

Рассчитывается мощность создаваемая винтом

, кВт (2.6.13.23)

тяга винта

, (2.6.13.24)

скорость концов винта

, м/с (2.6.13.25)

число Маха на концах винта

 

МК = WK / аH, (2.6.13.26)

 

число Рейнольдса на концах винта

 

Re = WK · BK / nK, (2.6.13.27)

 

где BK - хорда на конце винта, м.

 

Расчет мощности гидротормоза

 

NeГТ = (АГТ ·nS3) / 1000, кВт (2.6.13.28)

 

 

2.6.14.. Модуль узла типа Редуктор

Подпрограмма REDUKR

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

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

1. Текущее значение частоты вращения вала, на котором установлен узел п., об/мин, и суммарная мощность на валу от которого приводится данный узел SN , кВт.

Данные параметры передаются через группы "1000" и "2000".

2. Параметры узла, передающиеся в подпрограмму через массив "А" в составе информационной подгруппы модуля "10ВС100".

3. Характеристики узла передаются в подпрограмму через массив "А" в составе информационной подгруппы "10ВС300".

Выходные данные делятся на две части.

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

NeВВ - мощность на выходном валу, кВт;

nВВ - частота вращения выходного вала, об/мин;

МКР - крутящий момент на выходном валу, кН м.

2. Результаты расчета данного узла.

Переписываются в подгруппу результатов модуля "10ВС200".

Описание алгоритма модуля

 

Расчет без учета характеристик на барботаж и с учетом КПД зацепления

 

, кВт (2.6.14.1)

(2.6.14.3)

 

Рассчитывается потерянная мощность

 

NeПОТ = NeЗАЦ + NeБАР, кВт (2.6.14.4)

 

крутящий момент на выходном валу

, кН·м (2.6.14.5)

и частота вращения выходного вала

 

МКР = nТС / iРЕД, об/мин (2.6.14.6)

 

 

2.6.15. Модуль узла типа "Двухпозиционный переключающий клапан"

Подпрограмма KLAPAN

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

Входные данные делятся на три части.

1. Значения текущих параметров потоков воздуха (газа) считываются из групп текущих параметров массива "А" в соответствии с их расположением по контурам двигателя. К этим параметрам относятся:

- температура и давление торможения во входных сечениях клапана;

G, G - расходы воздуха (газа) во входных сечениях клапана;

qТ1В , qТ1С - относительные расходы топлива во входных сечениях клапана;

- энтальпии и энтропии торможения во входных сечениях клапана.

 

2. Параметры узла передаются в подпрограмму "KLAPAN" через массив "А" с составе информационной подгруппы "11ВС100".

3. Характеристики узла передаются в подпрограмму "KLAPAN" через массив "А" в составе информационной подгруппы "11ВС300".

Выходные данные делятся на две части.

1. Две группы значений текущих параметров потоков воздуха (газа) после прохождения переключающего клапана, передаваемые последующим модулям по каналам типа "Контур" (т.е. через группы 21000...25000 массива "А"). В их число входят:

- температура и давление торможения на выходе из узла (контуров "В" и "С");

G, G - расходы воздуха (газа) на выходе из узла;

qТ2В , qТ2С - относительные расходы топлива на выходе из узла;

- энтальпии и энтропии торможения выходящих потоков на выходе из узла.

 

2. Результаты расчета данного узла заносятся в информационную подгруппу результатов 11ВС200.

Описание алгоритма

 

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

этого контура не определяются. Алгоритм блока 2 реализован с помощью дважды выполняемого цикла. Внутри цикла происходит считывание текущих параметров потока соответствующего контура, последовательное обращение к подпрограммам "FUNKZI" и "BSP", вычисление значений "l" и "М". Найденные значения заносятся в подгруппу результатов. Выполняется расчет коэффициента потерь полного давления в контуре В. При неравенстве нулю номера используемой характеристики происходит обращение к подпрограмме PSP, в которой определяется значение параметра . Умножением на номинальное значение определяется искомая величина. Полученное значение sВ заносится в подгруппу результатов. Выполняется расчет коэффициента потерь полного давления в контуре "С". Корректируются текущие значения давлений в соответствующих контурах - РВ* и РС*. В зависимости от значения признака «Пр» происходит (Пр = 1) или не происходит (Пр = 0) замена текущих значений параметров потока в контурах "В" и "С". Замена параметров осуществляется с помощью одной переменной (буфера). Занесение оставшихся параметров в подгруппу результатов расчета.

 

 

2.6.16. Модуль узла типа Эжектор

Устройство перепуска или камера смешения 2

Подпрограмма EGECT

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

В алгоритме модуля предусмотрено использование характеристики вида

s = f(x), позволяющей учесть дополнительные потери полного давления с учетом их изменения по режимам работы.

Входные данные модуля делятся на три части.

1. Значения текущих параметров газовых потоков на входе в эти устройства, которые считываются в зависимости от значения разрядов "С" и "D" условного номера узла из группы "21000"..."25000" или из группы "5000" массива "А" (при расчете эжектора). В число этих параметров входят:

- температура и давление торможения на входе в узел;

Gi - расход газа (воздуха) на входе в узел;

qTi - относительный расход топлива на входе в узел;

- энтальпия и энтропия торможения на входе в узел.

Индекс i определяет номер контура и принимает значение в зависимости от значений разрядов "С" и "D" - условного номера узла.

2. Параметры узла. Передаются через массив "А" в составе информационной подгруппы "12ВС100", где "В" при расчете эжектора может принимать значения 0, 1, 2, при расчете устройства перепуска 3, 4, 5, 6, 7, а при расчете камеры смешения 8, 9

3. Характеристики узла. Передаются в подпрограмму через массив «А» в составе подгруппы «12ВС300».

Выходные данные этого модуля делятся на две части:

1. Текущие значения параметров газа после смешения на выходе из данного узла и передаваемые последующему модулю по каналам передачи информации типа Контур (т.е. через группы «21000»…»25000» массива «А»). Запись значений этих параметров в одну из этих групп осуществляется в зависимости от значений разрядов "С" и "D" условного номера узла. В число этих параметров входят:

- температура и давление торможения на выходе из узла;

G2 - расход газа на выходе из узла;

qT2 - относительный расход топлива на выходе из узла;

- энтальпия и энтропия торможения на выходе из узла.

 

2. Результаты расчета данного узла переписываются в информационную подгруппу результатов модуля «12ВС200».

 

Описание алгоритма

Поскольку в данном алгоритме реализована возможность расчета всех трех устройств, то выбор последовательности расчета по алгоритму определяется заданием соответствующих значений разрядов "В", "C" и "D" в условном номере узла NY = ABCD.

При расчете эжектора разряд "С" - номенклатура двигателя, через который проводится активный (эжектирующий) поток газа, а D = 0 - означает, что пассивный (эжектируемый) поток подводится из атмосферы. При расчете устройства перепуска разряд "С" указывает на номер контура двигателя, в который осуществляется перепуск газа из контура с номером "D". При расчете камеры смешения разряд "С" указывает на номер контура основного потока и одновременно номер контура, по которому будет направлена смесь газов. Разряд в этом случае указывает на номер контура, по которому к камере смешения подводится подмешиваемый поток газа.

 

Определяются площади проходных сечений на входе в узел.

FC = F2 / (1 + a), (2.6.16.1)

FD = F2 - FC, (2.6.16.2)

 

Рассчитывается устройство перепуска (2 < В £ 7). Переменная iTiP принимает значение 2. Корректируется значение расхода GD.

GD = GD · n ,

WPAR (3, 2) = GD/. (2.6.16.3)

 

Корректируется значение параметра

n = WPAR (3, 2) / WPAR (3, 1) . (2.6.16.4)

 

Рассчитывается эжектор (В2). Переменная принимает значение 1. Определяется значение расхода

 

GD = WPAR (3, 1) · n, (2.6.16.5)

 

Если пред активным соплом (контур "С") осуществляется впрыск охлаждающей жидкости, о чем свидетельствует dТ > 0, то уточняется значение температуры ТC*.

 

ТC*/ = ТC* - dТ , (2.6.16.6)

уточняются значения Нj* и Sj*

Определяется суммарный подошедший расход газа.

 

G1S = GC + GD , (2.6.16.7)

где GC = WPAR (3, 1), GD = WPAR (3, 2).

Номер активного контура iA принимает значение 1, т.е. контур С, а пассивного iPA=2

 

Проверяется условие

D*< PC*) Ú iT iP = 1. (2.6.16.8)

 

Если условие выполняется, переход к блоку 10, иначе к блоку 9, где iA принимает значение 2, а iPA=1, т.е. активным будет считаться контур "D", после чего также переход к блоку 10.

 

определяются критическая температура

 

ТKPA = f (TA* , qTA) ,

 

энтальпия

 

НKPA = f (TKPA, qTA),

 

газовая постоянная RA = f(qTA) ,

 

где ТA* = WPAR (2, iA), qTA = WPAR (4, iA).

 

Критическая скорость истечения из активного сопла

(2.6.16.9)

 

Расчет параметров пассивного потока.

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

 

RП = f (qТП), RП ® WPAR (7, iPA),

 

где qТП = WPAR (4, iPA).

Затем с помощью подпрограммы BSP по заданным значениям ТП* , РП*, GП, FП, qТП, RП определяются статические и критические параметры

 

За статическое давление активного потока принимается статическое давление пассивного потока, т.е.

 

РА = Р = WPAR (10, iPA). (2.6.16.10)

 

Давление по заторможенным параметрам

 

РА* = WPAR (2, iA).

 

Энтропия по заторможенным параметрам

 

SА* = WPAR (6, iA).

 

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

 

RА = WPAR (7, iA).

 

Площадь проходного сечения

 

FА = WPAR (8, iA).

 

Затем определяется энтропия по статическим параметрам

 

SА = SА* - RА · ln (PА* / PА). (2.6.16.11)

 

После чего с помощью подпрограммы FUNKZI определяются:

 

Статическая температура ТА = f (SA, qTA) ,

статическая энтальпия НА = f (TA, qTA) ,

 

скорость истечения

, (2.6.16.12)

 

Уточнение параметров активного сопла, работающего в критическом режиме.

 

СA = СKPA ,

ТA = ТKPA .

 

определяется энтропия:

 

SA = f (TKPA, qTA).

 

По уравнению адиабаты уточняется давление

. (2.6.16.13)

 

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

. (2.6.16.14)

 

 

Рассчитывается камера смешения. Определяется значение n.

 

n = WPAR (3, iPA) / WPAR (3, iA) (2.6.16.15)

 

Расчет параметров смеси.

Энтальпия по заторможенным параметрам

. (2.6.16.16)

 

Относительный расход топлива

(2.6.16.17)

определяется газовая постоянная

 

RСМ = f (qТСМ),

 

температура торможения

 

ТСМ* = f (НСМ*, qТСМ),

 

энтропия по заторможенным параметрам

 

SСМ* = f (TСМ* , qТСМ),

 

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

 

ТКРСМ = f (TСМ*, qТСМ),

 

критическая энтальпия

 

НКРСМ = f (TКРСМ, qТСМ),

 

Критическая скорость

, (2.6.16.18)

скорость смеси

 

ССМ = l2 · СКРСМ , (2.6.16.19)

 

энтальпия по статическим параметрам

, (2.6.16.20)

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

 

ТCM = f (НCM, qTCM),

 

энтропия по статическим параметрам

 

SCM = f (TCM, qTCM),

 

суммарный подошедший расход газа

 

GCM = GA + GП , (2.6.16.21)

где GА = WPAR (3, iA), GП = WPAR (3, iРA),

вспомогательный параметр

, (2.6.16.22)

суммарный импульс смеси газов

 

JS = GСМ ·ССМ + PF · 106 . (2.6.16.23)

 

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

(2.6.16.24)

Расчет при запирании сопла пассивного потока.

 

Р*СМ = Р*А ,

ССМ = lСМ · САКР , (2.6.16.25)

, (2.6.16.26)

RСМ = RА, qТСМ = qТА,

ТСМ = f (НСМ , qТСМ),

SСМ = f (ТСМ , qТСМ),

, (2.6.16.27)

. (2.6.16.28)

Расчет при критическом истечении из активного сопла

 

PMAX = 0,999 · PП* . (2.6.16.29)

 

Если РA > РMAX (ранее было принято РA = РП), то

 

РA = 0,9 · РMAX(2.6.16.30)

 

Определение критического давления в пассивном потоке

, (2.6.16.31)

где HП* = WPAR (5, iPA), CКРП = WPAR (12, iPA).

 

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

ТКРП = f (НКРП , qТП),

Энтропия

SКРП = f (ТКРП , qТП),

 

критическое давление (уравнение адиабаты)

, (2.6.16.32)

РMIN = РКРП . (2.6.16.33)

 

Решение системы уравнений при расчете статического давления в сечении выравнивания статических давлений двух потоков. В основе метода решения лежит комбинация методов Ньютона, дихотомии и штрафных функций. В процессе решения уточняются статические давления двух потоков и значение суммарной потребной площади F. В качестве исходного статического давления принимается статическое давление активного потока РА, в пределах границ его изменения от РMIN до РMAX, определяемых по формулам (2.6.16.29 и 2.6.16.33).

 

Определяется статическое давление смеси газов на выходе из камеры смешения

РCM = РF / F2P . (2.6.16.34)

 

давление заторможенного потока смеси

(2.6.16.35)

 








Date: 2015-05-04; view: 398; Нарушение авторских прав



mydocx.ru - 2015-2021 year. (0.208 sec.) Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав - Пожаловаться на публикацию