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

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

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

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

горения.

 

 

 

 

 

Блок JL L L (2 8 9 0 -3 1 8 0 )

- «блок управления работой

циклов пе­

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

раметров

и рк . В этом

же блоке осуществляется

организация

повторного счета оптимальных

вариантов РДТТ для

каждого ти­

па заряда

с выдачей на печать

полной информации

(3 0 0 0 -

3 0 4 0 ) .

 

 

 

 

 

2 .2 . Оптимизация однокамерного двухрежимного

РДТТ

Упрощенная блок-схема программы оптимизации однокамерно­

го двухрежимного РДТТ

0PTRJ32 приведена на рис.

а лис­

тинг фортран-программы -

в раздаточном материале.Структура про­

граммы и принципы ее построения в основном аналогичны програм­ ме 0PTRD1 рассмотренной выше. Блоки, обозначенные на схе­ ме рис. 2 .4 одноименными цифрами, имеют ту же структуру и назначение, что и в схеме рис. 2 .1 . Поэтому остановимся лишь на некоторых наиболее важных особенностях программы 0P T R B 2. Блоки I и П описания переменных и ввода исходных данных от­ личаются лишь количеством описываемых переменных и вводимых величин, необходимых для расчета РДТТ рассматриваемых типов. В частности, в ибХодных данных дополнительно задаются: при­

знак

типа оптимизируемого

двигателя - с регулируемым

( J1 IV =

= 4 )

или нерегулируемым

(1DY = 3 ) соплом; признаки

типов

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

Поскольку оптимизируемыми параметрами в программе

0 P T R B 2

являются рабочее давление рк и степень расширения

сопла J

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

мерного РДТТ с регулируемым соплом задается

также

отноше­

ние рабочих давлений на стартовом и маршевом режимах

р отн =

= р к / р К

, которое может задаваться в принципе

произвольно,

но в процессе оптимизации считается неизменным. Соответствен­ но несколько больший объем вычислений производится и в бло­ ке Ш.

В программе 0PTRJD2 организовано два цикла перебора ти­ пов и вариантов зарядов - внешний цикл перебора стартовых за - 31

|Д Ввод исходных данни ч— I Описание, дераменише и массивов

ШВ госдад » кеиотмт^гМ

'>М‘ 0>л*двоп='<?,

г«кЯ’ 7 Г

■ f c g g g p H g g S H © !

*&Шлк&тох<&та)4&**л(£г/гн

11

 

 

 

 

 

 

Конец

 

 

 

 

 

 

программа

 

 

 

'QVAD г

 

jy j Анализ результатов

 

т

 

 

 

 

и

 

 

. -■ -------,

 

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

 

 

^ ^ в т Г р д а т

й'коррввдия

 

 

 

 

 

I

______ ТЧТТГПГП • •

« Л Т Ч П А Т Л П Г Ш !

 

I

 

^

2 * а н Г с Т \

^Шй# 1|хода вычисления

 

 

|ШЕ

3

I

 

УПс[ Анализ резтльта^

I

I

 

сяммо(ояяа-я )

 

L

[V/w

 

 

 

тов расчета старто­

 

 

 

 

 

вого заряда и коррек-[

 

 

дх ан см )

 

foO£T*Q)

цжя хода вычислений

1П

I

 

 

 

 

 

 

 

 

|

 

V'*4x#fica)(a*/9#*x) |-----ЧУШм)Анализ резуль^

 

топлива

 

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

 

Тотов расчета мар-

 

 

Я

\*нву&(расчеттм )

 

кор^кция^хода*вы-

 

 

1

 

3-------J

 

числений

 

 

заряда

 

Печать'результатов|-

 

 

 

 

 

 

 

маршевоговариантов

 

|П Сравнение вариант

 

 

 

 

 

 

 

 

 

да

[нет

 

 

 

 

 

Перебор

 

 

 

 

 

 

 

 

 

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

Рис. 2 .4

рядов (начало цикла отмечено меткой 1 1 0 ) и вложенный в не­

го цикл перебора

маршевых

зарядов (начало - метка 1 0 ),

кото­

рые работают

точно так

же,

как и в программе 0PTRD 1

Со­

ответственные изменения внесены и в блок

ХП. Циклы перебо­

ра величин р к

и

построены так же, как и в программе

0PTRJJ1.

В качестве

переменны* , задающих в циклах перебора зарядов

типы стартовых и маршевых зарядов, используются

переменные

с идентификаторами К1

и К

соответственно,

а для

обозначения

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

отличающихся числом щелевых пропилов и т.п.,

- переменные

N1

и N

Принятая последовательность циклов пе­

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

Газодинамический расчет двухрежимного РДТТ в блоке ГУ. осуществляется с помощью СпП GDRD2.

Блок У, схема которого показана на рис. 2 .5 , расширен по сравнению с аналогичным блоком программы 0 P T R JJ1 , по­ скольку при расчете однокамерных двухрежимных двигателей мо­ жет возникнуть больше причин, по которым он не реализуется.

Если .газодинамический расчет в блоке 1ST завершился нор­

мально -

J

О ST= О,

расчет продолжается в блоке "ЗЙС. В слу­

чае

если

J

OST = 7

или J O S Т =10, коррекция вычислений про­

исходит так же, как в

программе

0P T R JJ1

В случае JO ST =

= 6

(ситуация) которая возникает

при расчете

РДТТ с регулируе­

мым соплом), если потребная площадь критического сечения соп­ ла на стартовом режиме превышает площадь среза сопла на мар-

шевом режиме, расчет производится заново при следующем уве­

личенном значении/^

или р к (если

увеличивать нельзя).

В случае JO S T -

8 £давление на стартовом режиме меньше

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

ме, как и в случае

JO S T = 7 .

В случае J O S T

= 9 (степень расширения сопла на марше­

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

отсечку части раструба сопла на маршевом режиме

(признак

JO T

= 1 ) , расчет продолжается в блоке 1У , иначе

расчет пре­

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

программы продолжается с

конца цикла перебора р к

(блок

1 2 ).

В

программе ОРТ ЯД 2

расчет внутреннего диаметра

каме­

ры для стартового и маршевого зарядов производится соответ­ ственно в блоках "У1 с и Y I M , а расчет зарядов - в блоках УЛс и У11м соответственно. Эти блоки аналогичны соответствующим блокам программы 0РТЯЪ1.

Блок УПСс - блок анализа результатов расчета стартовых

за­

рядов и коррекции хода вычислений отличается от блока YELU

 

программы

0Р Т Я Д 1

(см . рис. 2 .3 ) лишь тем,

что из

него

 

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

при

 

сT O S T -

5 ,

поскольку

в программе ОРТRЯ 2

расчет

старто­

вого заряда сигаретного типа не предусмотрен.

0 Р Т Я Я 1

 

Блок YIIIM аналогичен рлоку Y1IL программы

Во

всех случаях

J O S T Ф О в блоках УПТс и УШм на печать выво­

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

Блок I X . Обращение к СпП МЯУ2 , в которой производится расчет массы и основных размеров двухрежимного однокамерно­ го РДТТ.

Блок X аналогичен

соответствующему

блоку программы

0 P T R . E 1 но вместо

СпП

P R IN 1 используется СпП P R IN 2 .

Блок X L

такой же,

как

в программе

0Р Т .:Я 1

Блок X II

отличается

от такого

же блока программы 0РТЯЯ1

тем, что

дополнительно управляет работой цикла перебора стартовых за­ рядов.

2 .3 . Оптимизация двухкамерного РДТТ

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

стем влияют на условия заряжания стартовой камеры двигателя. Вместе с тем длина такого газовода, толщина его наружной

изоляции, а следовательно, и м ассе. маршевой ^амеры зависят от варьируемых параметров стартовой камеры. Однако совместное варьирование параметров двухкамерных РДТТ, число которых в двое больше, чем для однорежимного РДТТ, существенно увели­ чивает трудоемкость и время решения задачи. Вместе с тем взаимное влияние параметров разных камер друг на друга не очень сильно. Поэтому в настоящем пособии, учитывая его учеб­ ную направленность, для автоматизированной оптимизации двух­ камерных РДТТ, конструктивные схемы которых приведены на рис. 0 .1 , разработана программа 0PTRH 2H , в которой по критерию минимума суммарной массы камеры последовательно оптимизируются вначале маршевая, а затем стартовая камер» РДТТ и их заряды. Урощенная блок-схема ее приведена на рис. 2 .6 . В этой программе оптимизация каждой камеры произ­ водится по такой же расчетной схеме, как и для однорежимного

двигателя, с

помощью разработанной на

основе программы

0 P T R JJ1

специальной подпрограммы

O PK RJJ , учитывающей

конструктивные особенности двухкамерного двигателя. Расчет массы конструкции каждого варианта камеры в этой СпП осуще­ ствляется с помощью СпП M BV 2K При этом масса внутренне­ го газовода и промежуточного днища относится к массе перед­

ней камеры

двигателя.

_

 

 

Оптимизация параметров

( рк п /а )

и вариантов

топливного

заряда для

каждой камеры

в программе

0P T R U 2 K

осуществля­

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

метры заряд стартовой

камеры. В процессе оптимизации

каж­

дой камеры сокращенная

информация о результатах расчета

каж­

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

 

Рис.

2 .6

При этом, во-первых, если

это

необходимо, уточняются длина

и масса внутреннего газовода

с учетом оптимальных параметров

и размеров стартовой камеры

и, во-вторых, выводится на печать

с помощью СпП P R IN Z K полная информация о расчетных харак­

теристиках оптимального варианта РДТТ, достаточная для его

конструктивной проработки. В

случае если двухкамерный РДТТ

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

дится на печать.

В заключение заметим, что идентификаторы параметров стар­

товой камеры и заряда отличаются от маршевых дополнительным индексом " 1 " или "с"

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

РЕЖИМЕ ОБРАБОТКИ 3 .1 . Подготовка задания для оптимизации параметров РДТТ

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

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

цесса используется диск (том) с серийным именем A P S 6 0 2 , на котором указанные выше программы и СпП, имеющие свойимена, записаны в виде исходных модулей в одноименных раз­

делах библиотеки

S L IB и в

виде загрузочных модулей под те­

ми же ^менами в

библиотеке

Z A G L IB 2 .

Для решения задачи оптимизации некоторого типа РДТТ с использованием пакетного варианта С А П Р -602 пользователь должен подготовить пакет перфокарт (п/к), состоящий из двух частей:

1 ) комплекта с заданием на запуск и обработку необходимой программы и задаваемой пользователем исходной информацией

на языке

управления заданием;

 

 

 

2 ) комплекта п/к с числовыми исходными данными для вво­

да их с помощью оператора

REA11

 

 

Ниже в качестве примера приведен текст

задания на запуск

программы 0PTRJ32

из библиотеки S L IB

с распечаткой

тек­

ста программы на АЦПУ:

 

 

 

 

fISf6 6 0 14 / LJ JO B и ( E lP L ,

К6 0 2 , 1 4 0 0 ), MS&LEVEL=(2,0),

10

// u u u u u REGION = 100К, TIME = ( 3 ,0 )

 

20

// W u u EXECu FORTGCLG, PARM. FORT = SOURCE

JO

// FORT. SYSIN t_i EE uuESN=SLIB (OPTRE2 ), E IS P =SH R,

4 0

/ / u u u u

UNIT = SYSEA , VOL = SER = A P S 602

 

SO

I I I КEE. SYSLIB LJ ЕЕ и ESN = SYS1. FORTLIB, EISP = SHR

60

//u u ЕЕ и ESN=LOAE,EISP =SHR, UNIT=SYSEA, VOL=SER PS602

70

//uuEE^ESN=ZA&LIB2,SISP=SHR^ UNIT=SYSEA,V0L=SER*APS602

80

IIGO. F T 0 6 F 0 0 1 u E E u SYSOUT=A, ECB = BLKSIZE =132 ,

90

//«uuw SPACE = (CYL, ( 4 ,0 ,0 ) ) , UNIT = SYSMA

100

// GO. SYSIN и Е Е u * комплект

п/к с исходной

110

 

 

числовой

информацией

 

 

/ *

 

 

 

 

 

 

//

 

 

 

 

 

 

П/к N° 1 0

и комплект

п/к с исходной числовой информацией яв­

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

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

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

п/к

]sfe

3 0 - 1 1 0

нужно заменить перфокартами следующего содержа­

ния:

 

 

 

 

// J/UJU

EXEC и PGM = OPTED2

30

 

HSTEPLIB^DDU UNIT^SYSDA, VOL^SER =APS602yDISP~SHRy 40

 

ff u u и u DSN = ZAGLIB2

 

SO

 

//FT06 FQ01кл DDXJ SYSOUT=A,DCB = BLRS1ZE=732,

60

 

/ / o u L iu

SPACE =(CYL, (4,0,0)), UNIT=SYSDA

70

 

// FT05F001 L J D D LJ *

 

80

 

Далее следует комплект п/к с исходной числовой информаци­

ей

и т.д.

 

 

 

 

Обратим внимание на п/к №

8 0 , 9 0 , которые обеспечивают

печать информации в строках с

количеством символов больше

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

Правила комплектования п/к с исходной информацией приведе­ ны ниже в п. 3 .2 *“ 3 .4 .

3 .2 . Исходные данные для оптимиздпии однорежимного РДТТ и топливного заряда (программа 0PTRD 1)

Числовые значения

вводимых в

программу вегичин набивают­

ся на

1 2 перфокартах в

 

порядке и формате, указанном в табл.3.1.*

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 3 .1

Номер

 

 

 

 

 

 

 

 

 

 

перфо­

 

Вводимые

величины31

 

 

Формат3

карты

 

 

 

 

 

 

 

 

 

 

1

Р т

о п

 

t

cr у

t

min, *

t

max

F 1 0 .4

 

у AS у г /f» Гн max

>

 

 

 

2

у

* f pA у ^о$ > ^8t об Р об » ^ое

F10.ii

 

 

3

d те *d CL м ах, Ч о P c iV n Vz

 

 

F 1 0 .4

 

 

 

 

4

($ р ) , ЛИ2 , р мз ,

 

 

>Р бр у ^ адг ’ Р а д г

FW .4

х Наименование и размерность вводимых величин указаны в

списке

основных обозначении

 

(см. приложение 1 ).

 

* * Вместо

спецификации

преобразования F 1 0 .4

мбжно ис­

пользовать спецификацию

 

Ё 1 0 .4 , аналогичную по действию при

вводе.

 

 

 

 

 

 

 

 

 

 

5

L

d Т

2

 

F 10.9

 

 

г у “'г »*к >

 

 

6

Рзд-сТ) ^ н -с т у

Рм.допг

 

F I0 .9

7

 

fa та

»^ а уpKini >Ркт а »

^доп

F1Q.9

 

 

 

8

MZ(t), MZ(2),MZ(3),MZ(4), MZ(S),MZ(6), MZ(7), MZ(8)

11

9

MS(1), MS(2), MS(3), MSW

 

11

10

MW), ML(2), ML(3), ML(9), Ml(5), M W ), ML(7), ML(8),ML(9)

12

11

MH(1), MH(2), MH(3),MH(9),MN(5), MN(6), MH(7),

22

 

MH(8), MH(9), MH(10)

 

12

ic , J U ,

J C T

 

 

 

 

 

 

 

12,11,11

Гак,

например,

если Zr =

0 ,5 ;

dr =

2;

Тк =

2 9 5 0 ; Л^=0 ,8 см/с,

то н а“5-й перфокарте

в

первых

40 колонках должны

быть на­

биты

следующие символы

( и и

— символ

пробела):

 

u UUU 0»5 UJLJLJL-ILJUJLJUJ 2 фLJLJUJLJLJ

ZSSO.LJUJLJULJLJLJU 0,8 UJLJ UJ

1

 

 

11

 

 

 

21

 

 

31

W

Если

ic = 2 ,

JU = 0 ,

JC T =

1 ,т о

на

12-й

перфокарте

должны

быть

набиты в

первых

4 колонках символы

и 2 0 1 .

 

Массивом

AfZ, набиваемым на 8

п/к, в

программу вводят­

ся признаки

типов зарядов

( JZ A R

),

которые пользователь за­

дает в качестве альтернативных при оптимизации РДТТ. Число­

вые

значения признаков JZ A R (см.

с. 6 , 7 )

набиваются в

пер­

вых

8 колонках перфокарты подряд

в порядке

возрастания,

на­

чиная с первой колонки. Если число альтернативных типов заря­ дов меньше 8 , то в оставшихся колонках вплоть до 8-й набива­

ются символы

0 .Например,

если альтернативные

типы зарядов -

щелевой заряд

типа

SZC

 

( JZ A R

= 3 ) , заряд со звездообраз­

ным каналом типа

ZV ( JZ A R = 4 ) ,

одношашечный или много­

шашечный заряд типа TZN

( JZ A R

= 8 ),

тогда

на 8-й п/к на­

биваются символы 3 4 8 0 0

0 0 0

На этой п/к должен быть на­

бит хотя бы один символ, отличный

от

и не превышающий 8 .

Массив MS вводит в

программу количество щелевых пропилов

{ п щ )» задаваемое

пользователем

для различных

вариантов ще­

левых зарядов

типа SZJJ

или SZC

 

В программе

предусмотрено

задание до 4 вариантов с количеством щелевых пропилов от 3 до 6. Числа набиваются в первых четырех колонках п/к N? 9 подряд в порядке возрастания, начиная с первой колонки.

Массив ML вводит

в программу количество лучей звездооб­

разного канала ( )

для альтернативных вариантов заряда ти-

 

3 9

па ZV , задаваемых пользователем. Предусмотрена возможность

расчета зарядов с числом лучей от 4

до

1 2 . Числа набиваются

на 10-й перфокарте в колонках с 1

до

18-й подряд в поряд­

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

Массив МН вводит в программу количество шашек тгш для различных альтернативных вариантов трубчатых зарядов типа TZN. Предусмотрена возможность расчета зарядов с количеством ша­

шек

1 ,3 ,4 ,5 ,6 ,7 ,1 2 ,1 3 ,1 8 ,1 9 . Числа шашек для

альтернативных

вариантов заряда типа TSN набиваются на п/к N°

1 1

в

первых

20 колонках (на каждое число две колонки) подряд

в

порядке

возрастания

 

 

 

Если число вариантов, задаваемых массивами MS

, ML или

МН

, меньше максимального предусмотренного размерностью

массива, то в оставшихся колонках соответствующей перфокарты

пробиваются символы 0

Так, например, если для заряда типа

TSN

заданы варианты

заряда

с числом шашек 5 ,6 ,1 2 ,1 8 , то

на п/к № 11 в первых

20 колонках должны быть набиты симво­

лы

и 5 и 6 1 2 1 8 0 0 0 0 0 0 0 0 0 0 0 0

 

Если массивом MZ не предусмотрен расчет зарядов

типа SZJ)9

S Z C

, Z V или TZN ,

то элементам соответствующих

массивов

MS

, ML или МН могут быть

заданы произвольные значения, но

в нужных форматах, например

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

бы один элемент (первый) соответствующего массива должен от­ личаться от 0 .

Программой предусмотрено задание величины Л д оп (п/к № 7) в двух формах. Если пользователь знает, например, по прототипу или предыдущим расчетам допустимую величину относительной

скорости потока в каналах заряда, то

Л доп задается в

относи­

тельной форме и должна быть меньше

1 (практически

0 ,8 ).

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

ной скорости газового

потока

в каналах заряда

(например, поро­

говая скорость эрозионного горения), тогда

Л доп

задается,

как

эта скорость, с размерностью

м/с,и Л#оп >

1. В

этом случае

от­

носительная

величина

Л$оп

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

в

программе.

 

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

ную печать

под заголовком "Исходные данные"

в порядке ввода

в форме: <

идентификатор входной величины > =

< значение >

 

Массивы выводятся в

форме:

< идентификатор массива >= < зна­

чения элементов, разделенные

запятыми >

 

 

 

 

Во избежание аварийных остановов вычислительной системы или получения абсурдных результатов некоторые вводимые вели­ чины должны удовлетворять следующим условиям:

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