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

книги / Оптимизация параметров РДТТ и топливных зарядов при автоматизированном проектировании летательных аппаратов

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
2.39 Mб
Скачать

ся от воздействия газов топливным зарядом (типа ТСН и ТНС), то он снаружи не теплоизолируется. Абсолютный диаметр газово-

да определяется по формуле dr = dr d кр

, где d Кр -

критичес­

кий диаметр сопла РД ТТ, рассчитываемый

в программе

оптими­

зации РДТТ. Для телескопической схемы РДТТ диаметр внутрен­ него газовода dr3 принимается равным диаметру внутренней ка­

меры

двигателя DKc> М асса

стыковых узлов газовода принимает>-

ся пропорциональной квадрату ^диаметра газовода

+П'СТГ~

Ю <£*•

Таким образом, масса внешнего газовода рассчитывается

по фор­

муле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'>5 i

3 f u3}

,0 d r ,

(1 .1 3 )

а внутреннего газовода

- по формуле

 

 

 

 

 

газ

\(d - S ) S р M d - S

-1,58

)1 ,S $

р

} + W d*

(1 .1 4 )

гз Lv гз гз/ гз'об 1

гз гз

7 из d

из i J из 3

г3

 

 

1 .2 . Краткое

описание

СпП расчета

массы РДТТ

На основе схемы расчета массы РДТТ, рассмотренной выше, разработаны три специальные подпрограммы расчета массы РДТТ и основных геометрических параметров камеры: MJJV1 для однорежимных однокамерных РДТТ; MUV2 - для двухрежим­ ных однокамерных РДТТ и MDV2K- Для двухрежимных двухка­

мерных РДТТ. В этих СпП рассчитываются массовые и геометри­ ческие характеристики РДТТ и основных его элементов, такие, как fnAB , (ХАВ , L гз , 1 С, LAy , L A B и др: Входные переменные для работы этих и других СпП, используемых в программах опти­ мизации двигателей, которые необходимо задавать в исходных данных, приведены ниже в п. 3 .2 - 3 .4 .

Остановимся на некоторых особенностях использования СпП МЛУ2 К При однократном обращении к этой СпП рассчиты­ ваются параметры только одной камеры двухкамерного двигате­ ля - стартовой или маршевой в зависимости от состава входных величин. Поэтому для получения полных данных РДТТ необходи­ мо двукратное обращение к СпП и суммирование результатов, ко­

торое осуществляется в

программе О РТRВ 2 К .

При расчете

маршевой камеры JTZ =

2 , а

для стартовой

J T Z = 1.

Расчет массы дополнительных элементов в вышеприведенных

СпП осуществляется также в специальной подпрограмме

SUBROUTINE

МП OP

,

обращение к

которой

осуществля­

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

подпрограмм.

 

 

Вычисление переменной Л

, KOTopasi используется в приведен­

ных выше СпП, а также при расчете зарядов, производится в

SUBROUTINE

J3KBH.

 

 

 

 

Г л а в а 2 . ПРОГРАММЫ ОПТИМИЗАЦИИ ПАРАМЕТРОВ РДТТ И ТОПЛИВНЫХ ЗАРЯДОВ

2 .1 . Оптимизация однорежимного РДТТ

Упрощенная блок-схема прораммы оптимизации однорежимно­ го РДТТ - 0PTRJ11 показана на рже. 2 .1 . (см. примечание

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

В числе исходных данных задаются: признаки типа и конструк­ тивных особенностей оптимизируемого РДТТ; признаки типа и вариантов топливных зарядов, которые пользователь программы считает целесообразным просмотреть при оптимизации из числа типовых; требуемые внешние характеристики РДТТ и условия его работы; характеристики конструкционных материалов и твер>- дого топлива; диапазон и шаг изменения варьируемых значений рабочего давления в РДТТ р к и степени расширения сопла Более подробный перечень входных величин приводится ниже при

описании правил использования программы

0 P T R B 1

Для кон­

троля все вводимые величины выводятся на печать.

 

В блоке Ш ( 6 2 0 - 1 2 6 0 )*' вычисляются

некоторые

константы

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

граммы: PML, РМА

} &0РТ7 J0P }JJP1 и др. В этом

же блоке

анализируются исходные

данные, и если в них имеются

ошибки,

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

В прюграмме

0PTRB1

в

качестве переменной,

задающей

тип

заряда, в цикле

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

J

идентификатором

К

» которой

в оператора

с меткой 1 0 при­

 

сваивается значение

признака

типа заряда

JZ A R

очередного

 

* В скобках указаны номера строк в листинге программы.

рГ Описание переменных и массивоЁГ]

[П|Ввод исходных данных |—»-[шВычисление констант;/*/

р * ------------------------------

S -—< ^

 

7 >

-

 

 

г

J*ei

%9%/nc/ii<&/naxda/na>d&ncndLni

2 0 V

---------

 

 

 

r € > - /7ъ&ит 0

vLnin

 

 

 

L.

ЕН

Ен

Е

&

•о

04

«

I

I

s s

©

н

о

04

о

х>

&

е

R

GDMl

I------ KJ У1Анализ результатов L.»_

1—

----- --- ■

J

газодинамического

I

 

 

 

Расчета РДТТ и кор-

 

 

 

ш------- рекция хода вычисле­

 

----------1------------ -

ний

__________

 

У1111 Анализ

резулъта-

 

 

atZMxeeg**; |-

 

 

 

 

тов расчета зарядов

 

 

 

 

и коррекция хода вы­

 

ТгТMDVJCoacvtvn/7?#)^ ^ Га{Н числений

 

 

3

т

Полная печать

результатов

.Сокращенная печать результатов

да

Л^дбо/г~У^&6id&mindL

Л / / \

dast£i*&Jii -~{8ду~

1

- U -

[нет

Перебор вариантов топливного заряда

Рис. 2 .1

альтернативного варианта из числа заданных пользователем в массиве M Z . Переменной, задающей вариант заряда, отличающий­ ся числом щелевых пропилов, количеством лучей звездообразного канала или трубчатых шашек, является переменная с идентифика­ тором N , которой в операторах ( 1 3 1 0 , 1 3 3 0 ^ 1 3 5 0 ) при­ сваивается очередное значение из числа заданных пользователем в массивах MS , ML или МН соответственно.

Цикл перебора типа зарядов организован

с помощью операто­

ра присваивания и начинается

оператором с

меткой 1 0 ( 1 2 7 0 )

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

меткой 9 5 и последующим услов­

ным оператором ( 3 1 5 0 - 3 1 7 0 ) , которые осуществляют переход к расчету следующего типа заряда или выход из цикла. Если для очередного варьируемого типа заряда, например щелевого, звез­ дообразного или многошашечного трубчатого, пользователем про­ граммы предусмотрен расчет различных его вариантов, отличаю­ щихся количеством щелевых пропилов, лучей звездообразного ка­ нала или шашек, то перебор этих вариантов организуется допол­ нительными циклами, начинающимися вычисляемым оператором

перехода

&0 ТО . . . ( 1 3

0 0 ) и операторами с метками

 

1 2 ,

1 3 ,

1 4 -

и заканчивающимися

операторами сметками

9 1

,

9 2 ,

9 3 ,

9 4

и 9 5

( 3 0 6 0 - 3 1 6 0 ) ,

в

которых осуществляется

переход

к

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

дельно допустимый

свод ‘горения меньше, а коэффициент попе­

речного заполнения больше при прочих равных условиях.

 

Цикл перебора значений рабочего давления в

РДТТ рк

начи­

нается оператором с

меткой

2 3 ( 1 5 2 0 ) и заканчивается

опера­

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

3 6

и 8 7 ,

в которых переменная цикла

Р К

увеличивается на величину заданного шага U P К

или скорректиро­

ванного шага U P f

, ^ и условным оператором с

меткой 8 8

( 2 9 8 0 ) , организующим возврат к началу цикла или* выход из него, если перебор закончен. Этот цикл при расчете каждого варианта заряда начинается с минимального заданного пользова­ телем значения давления в РД ТТ, которое затем увеличивается с заданным в исходных данных шагом до тех пор, пока не будет превышено максимальное заданное в исходных данных значение рабочего давления. В процессе работы программы минимальное давление в РШ 'Г, при котором проводятся расчеты, может быть скорректировано в блоке анализа и коррекции, о чем будет ска­ зано ниже.

2 4

Последний внутренний цикл перебора значений степени расши­

рения сопла

,

организованный так же, как и все внешние, на­

чинается пустым

оператором с

меткой

2 5 ( 1 5 7 0 ) и заканчива­

ется оператором

с меткой 8 0

( 2 8 9 0 )

и последующими тремя

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

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

Блок ТЕ ( 1 5 8 0 - 1 6 8 0 ) - обращение к СпП GJDRB1 , в ко­ торой производится газодинамический расчет однорежимного РДТТ и определяются основные параметры двигателя и топливн^ го заряда [4 ] . При некоторых сочетаниях варьируемых i. парамет­ ров практическая реализация двигателя может оказаться невоз­ можной, например, если варьируемое значение рк задано меньше минимально допустимого из условий полного сгорания топлива

р км , которое вычисляется в этой

СпП, или если

при заданных

значениях

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

безотрывное

течение газа

всопле или диаметр среза сопла превышает допустимое значение

^ат ах 9 заданное в исходных данных. В этом случае работа под­ программы GJ7RJJ1 прекращается, служебной переменной JQ ST

присваиваете^ значение 7 ,

в случае если

меньше минимально

допустимого p kMiили 1 0 ,

если не обеспечивается безотрывное

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

 

При нормальном завершении работы СпП GJDRJD1

переменная

J 0 S T

=■ 0

После завершения работы СпП переменным

JF A

и J Р К

присваивается значение 1

для исключения

повторяющихся

расчетов в

СпП GURJ11

[ 4

J , а при нормальном

завершении это

же

значение присваивается переменной JM Z

блоке Е Л ).

 

 

Блок У

( 1 6 9 0 - 2 0 8 0 )

анализа

результатов

газодинамическо­

го

расчета

и коррекции хода

вычислений

(упрощенная

схема

бло­

ка

показана

на рис. 2 .2 ) . В

случае, если

перемейная

J 0 S T - 0 ,

расчет

продолжается в блоке

У-L,

иначе на печать выводятся со­

общение о том, что РДТТ не реализуется и значения варьируе­

мых параметров р к и

,

служебной переменной

J ОST , а так­

же минимально допустимое

рабочее давление в камере р Км или

максимальное реализуемое

значение степени расширения сопла

f m ax И кРаткая информация о коррекции хода

вычислений.

Нели JO S T = 7 , наименьшее варьируемое значение р к коррек­ тируется и принимается равным минимально допустимому значе­

нию Ркм (1 9 9 0 ,

1 4 9 0 ) , рассчитанному в

СпП & BRB1

По­

скольку диапазон

варьируемых значений р к

задается пользова­

телем произвольно, может оказаться, что скорректированное в блоке У минимальное значение р км станет больше заданного в исходных данных максимального варьируемого значения В этом случае работа программы оптимизации прекращается

(2 0 0 0 ) и на печать выводится соответствующая информация. В другом случае корректируется первый шаг перебора давления

d p (2 0 4 0 ) , чтобы сохранить заданные пользователем варьи­ руемые значения рк , превышающие минимально допустимое зна­ чение, и далее расчет вновь производится с начала цикла пере­ бора р к (переход к метке 2 0 ) .

Рис. 2 .2

В случае J О ST - 1 0 дальнейший ход вычислений зависит от предельного реализуемого значения степени расширения соп­ ла , рассчитываемого в СпП G-BRB1 Может оказать­ ся, что эта величина меньше минимального варьируемого значе­

ния j а

, заданного в исходных данных,

так

как последняя

может

задаваться пользователем произвольно.

Из

физических со ­

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

допустимое значение f fn ax

было

не меньше

минимального за­

данного значения ( 2 0 7 0 ) .

Однако

указанная

коррекция опять

может привести к тому, что скорректированное значение р

• К fPilft,

как и в случае JO S T = 7 , станет больше максимального варьи­ руемого значения, заданного в исходных данных. Поэтому даль­

нейший ход вычислений в этом случае аналогичен случаю

JO S T =

= 7 ( &0 ТО 3 6 ) . Если предельная реализуемая величина

f

больше минимального варьируемого ее значения, заданного в ис­

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

и зарядов при

следующем большем значении варьируемого параметра р

(пере­

ход к оператору с

меткой 8 5

(2 9 3 0 )) .

 

 

Б л о к а х ( 2 0 9

0 - 2 1 2 0 ) . В

блоке рассчитывается толщина

стенки камеры, а

также с помощью СпП В К В И

внутренний

диаметр камеры и толщина изоляционного покрытия. Заметим,

что в тех

случаях, когда пользователь использует

программу

0 P T R B 1

только для выбора наиболее рационального типа или

варианта

заряда при одних и тех же значениях

ъ рк

, рас­

смотренная выше последовательность работы программы реализу­ ется только для первого заряда, в последующем блоки 33? и У обходятся, так к$к вычисляемые в них величины не меняются при

неизменных J

и р к

Для этого используется служебная перемен­

ная JA IZ в операторе

( 1 4 7

0 ) .

Блок УН,

( 2 1 6 0 - 2

1 8 0 ) .

Обращение к СпП Z A K JJO 9 в ко­

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

ным причинам. В этом случае служебной переменной

JO S T

при-•

сваивается определенное условное числовое значение, соответст­

вующее причине, по которой реализация заряда невозможна (см.

п. 3 .5 ) , и осуществляется выход из СпП Z A R JJ) В

случае

 

нормального завершения работы этой СпП переменная

JO S T

=0.

Блок УШ - блок анализа результатов расчета топливных за­ рядов и коррекции хода вычислений ( 2 1 9 0 - 2 4 5 0 ) . Упрощенная схема блока приведена на рис. 2 .3 . При нормальном завершении

работы СпП Z A R j p ( J 0 5 T = i 1)

работа программы

0P T R D 1

продолжаете:: в блоке I X . Если

реализация заряда

невозможна,

на печать выдается сообщение об этом и информация о типе и

варианте заряда,

значениях варьируемых параметров рк ъ

f a

и признака JO ST

, а также расчетные значения внутреннего

диаметра камеры

РД ТТ, толщины обечайки камеры, свода

горе­

ния заряда и относительной скорости газового потока в каналах заряда, если она рассчитывалась, а для сигаретного заряда еще и необходимое для реализации заряда увеличение скорости горе­ ния топлива [АТС]. Эта информация может использоваться при необходимости для коррекции исходных данных.

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

ка JO S T

В случае JO S T

= 1 (рассматриваемый вариант за­

ряда при заданном на этом

шаге значении р не реализуется по

своду горения), если для рассматриваемого заряда при других сочетаниях варьируемых параметров р к и /а , рассмотренных ранее, двигатель мог быть реализованх , то расчет повторяется для оптимального сочетания этих параметров из числа рассмо'г-

ренных

ранее с- выводом полной информации на печать ( 2 3 1 0 -

2 3 2 0 ,

3 0 0 0 - 3 0 4 0 , 2 5 0 0 - 2 5 3 0 ) ; далее расчет продолжается

для другого типа или варианта заряда.

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

сочетании ранее рассмотренных параметров рк

и f

реализация

оказалась невозможной ( Л = О ), то в этом

случае

дальнейшая

работа программы заЬисит^ от тогЬ типа заряда, который про­ сматривается на текущем шаге.

Напомним, что заряды в программе оптимизации перебирают­ ся в порядке уменьшения предельно допустимого свода горения. Поэтому если рассматриваемый заряд имеет признак JZ A R у

не превышающий 5 , и не

может быть реализован по своду горе­

ния даже при минимальном варьируемом значении р к

, то заря­

ды с большим значением

признака JZ A R и подавно не

могут

быть реализованы. Если рассматриваемый заряд является трех­

шашечным

(тип TCH

J Z A R

= 6)

или двухшашечным телеско­

пическим

(типа ТНН,

J Z A R -

7 ),

то имеет смысл просмотреть,

если это предусмотрено пользователем, еще одношашечный труб­ чатый заряд типа TZN %у которого предельно допустимый свод

х Признаком возможности реализации некоторого варианта

двигателя и заряда служит служебная переменная Л , которой в этом случае присваивается значение 1 ( 2 8 0 0 ) .

горения больше; в противном случае работа программы заверша­ ется. Рассмотренные логические операции реализуются операто­

рами

( 2 3 3 0 - 2 3 5 0 ) .

 

 

В

случае J O S T = 2

или J O S T - 4 (заряд не

реализуется

по предельно допустимой

скорости газового потока

в каналах

топливного заряда) расчет продолжается при большем значении степени расширения сопла или большем давлении в камере РДТТ,

если

увеличение

недопустимо

(

2 3 6 0 ,

2 3 9 0 ,

2 4 0 0 ) ; и то и

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

скорость

газового

потока

в канале заряда.

В

случае J

0 6 Т = 3

(заряд

не реализуется

по требуемому

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

ров с

выдачей на печать полной информации ( £ 3 6 0 - 2 3

8 0 , 3 0 0 0 -

3 0 4 0

, 2 5 0 0 - 2 5 3 0 )

и последующим переходом к рассмотрению

другого типа или варианта заряда; если при других сочетаниях

значений рк и р е а л и з а ц и я была невозможной, осуществляется

переход к рассмотрению другого типа или

варианта

заряда

( 2 3 7 0 ) . В случае

JO S T = 5 (потребная площадь горения сига­

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

камеры

РДТТ)

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

требную площадь горения

заряда и з-за увеличения скорости го­

рения.

JO S T = 1 3

 

 

В случае

(заряд

со звездообразным каналом

после коррекции формы канала в СпП ZV не обеспечивает требуе­

мого постоянства периметра горения), если при других ранее

рассмотренных

сочетаниях

рк и f a

реализация была возможна,

расчет повторяется для оптимального сочетания этих параметров (из ранее рассмотренных) с выдачей на печать полной информа­ ции и последующим переходом к расчету нового типа или вариан­

та заряда

так же, как и в случае JO S T = 3 ; иначе расчет про­

должается

сразу для варианта с большим числом лучей звезды,

если это предусмотрено исходными данными, или для другого ти­ па заряда.

Блок IX ( 2 4 6 0 - 2 4 9 0 ) . Обращение к СпП /Ш7К/ , в кото­ рой производится расчет массы и основных размеров однорежим­ ного РДТТ.

Блок X . Блок печати ( 2 5 0 0 - 2 7 1 0 ) . В этом блоке осуществля­ ется выдача на печать полной или сокращенной информации о ре­ зультатах расчета. Полная информация о конструктивно-баллисти­ ческих параметрах РДТТ й топливного заряда выводится на пе­ чать с помощью СпП P R IN 1 ( 2 5 1 0 - 2 5 3 0 ) для каждого из

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

тания варьируемых параметров р к и

. Для остальных вариан­

тов на печать

выводится сокращенная информация.

Блок 33 . (

2 7 2 0 - 2 8 8 0 ) . В этом

блоке осуществляется срав­

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

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

и 1сьоп

( 2 7 8 0 - 2 7 9 0 ) .

В этом же блоке осуществляется коррекция минимального

варьируемого значения fc u m in

9 с котоРого начинается пере­

бор при различных значениях параметра

р к . Как уже отмеча­

лось выше, для каждого варианта заряда перебор значений сте­ пени расширения сопла при минимальном давлении начинается с минимального заданного в исходных данных значения / ЛАА1 Пе­ ребор может осуществляться в диапазоне, не превышающем за­ данный в исходных данных, но, как правило, в процессе счета этот диапазон постоянно корректируется так, чтобы исключить расчеты с явно неоптимальными значениями f a

Как известно, оптимальная по весовому критерию степень рас­ ширения сопла РДТТ несколько меньше степени расширения так называемого расчетного сопла, у которого давление на срезе равно атмосферному, и это оптимальное значение зависит при прочих равных условиях от давления в РДТТ. С увеличением давления оптимальная степень расширения сопла растет. В соо'г-

ветствии с этим в программе

0 P T R B 1 перебор значений

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

тех пор, по­

ка очередной шаг не приведет

к увеличению массы

РДТТ по

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

мом значении J

не обеспечивается безотрывное течение газов

в сопле или диаметр среза

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

В этом случае в

качестве

оптимального значения f принимает­

ся последнее варьируемое значение, при котором реализация дви­

гателя была возможна.

Оптимальное значение у

, полученное при

очередном значении р~к

, используется как минимальное при пере­

б о р е н а

следующем шаге

по р'

, так

как меньшие значения

будут явно хуже во

всех

отношениях.

Однако для щелевых

зарядов перебор варьируемых значений^

всегда

начинается

с

заданного

минимального значения

.

так как щелевой

за -

ряд,

реализуемый при малом

значении f

, может

оказаться

не­

реализуемым при большем. Дело в том, что. с увеличением сте­ пени расширения сопла и ростом при этом удельного импульса гяги потребная площадь горения заряда уменьшается. Это при

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