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

книги из ГПНТБ / Дракин, И. И. Основы проектирования беспилотных летательных аппаратов с учетом экономической эффективности

.pdf
Скачиваний:
37
Добавлен:
21.10.2023
Размер:
8.98 Mб
Скачать

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

а

bd

G f

Go

 

где ~ — знак пропорциональности, d — диаметр бака, при этом

d ~ Q f .

Критическое напряжение сжатия в цилиндрической оболочке, как известно,

акр

d

O f '

 

 

 

следовательно, учитывая, что ст = сгКр, получим

8 — O f .

Вес бака пропорционален толщине (приведенной) стенок б и его поверхности

S —Q2'3 .

следовательно,

^ Н) = ^ - = %.несO f -

Коэффициент масштабного эффекта для несущих баков низкого давления

г’ б.нес —__

-f Гі

^б.нес^

£

 

J 0U 0

— “ g1" /о^„Гн',б .н ес’

 

ИЛИ

Ce.Hec = - f /o^afäec-

(2.10)

 

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

.зависит от веса ЛА. Относительный вес ненесущих (н. н) баков низкого давления

_ і _

^бін.н^^б.н.н^о .

Следовательно,

чб.н,н

------ / > гаін) .

Г

__

£ J0ѵ б.Н.Н

( 2 . 1 0 ' )

80

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

Для определения коэффициента аэродинамического масштаб­ ного эффекта воспользуемся формулами (1.42) и (1.43) и пере­ пишем выражение (1. 61) в виде

Учитывая формулу (1. 60),

ди0

 

 

(2.11)

°0ПДІ- = - - Д - !\ а ,

-

 

следовательно,

 

3

(2.12)

^а-=------ ^о:Ата ( 1+ «)-

 

Следует заметить, что величина рта в формуле (2. 12)

должна

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

так как

последнее не имеет «масштабного эффекта».

 

 

Общее значение коэффициента «масштабного эффекта»

С =

Са + С„.п + :б,

 

(2.13)

где £а, £н.п, £б определяются по формулам (2.9), (2. 10),,(2. 11), (2. 12). В выражении (2. 13) не учитывается масштабный эффект двигателей, который для РДТТ и Ж РД практически равен нулю. Для ВРД (ТРД, РПД и ПВРД) масштабный эффект близок к нулю. Для ТРД теоретически £дВ~ Р 1/4, однако практика разви­ тия ТРД последних лет за рубежом не подтверждает этого.

Пример.

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

ступени зенит­

ной ракеты,

рассмотренной в

гл.

I, разд. 7.

Принимаем:

f0

= 3,73; рф =0,3;

Цкр = 0,06; [Лоп =

0,01; «о = 0,1;

а =

0,18; рітa =

0,1795; рт =

0,4252.

По формуле

(2. 12)

 

 

 

 

 

Са = — — -3,73-0,1795-1,18 = — 0,263,

3

по формуле (2. 9)

С„.п = — -3,73-0,07 =,- 0,087,

3

по формуле (2. 10)

Иб нес = — -3,73-0,4252-0,1 =0,0264,

6

следовательно, по формуле (2. 13)

С = - 0 ,1 5

81

По формуле (2. 5)

/ рр

3,73

— = 2 ,5 8 ,

1,3 + 0,15

по формуле (2. 7)

АСо = 2,58-1,3-ДОп.„ = 3,35 ДСіп.

4.ВЛИЯНИЕ ИЗМЕНЕНИЙ ВЕСА КОНСТРУКЦИИ

ИЛОБОВОГО СОПРОТИВЛЕНИЯ Л А НА ПОЛЕТНЫЙ ВЕС

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

dO0= d O K-\-dOT,

полезная нагрузка считается постоянной.

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

ва будет

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

рт,

вызванного изменением

лобового

сопротивления,

а также

за счет факторов, приведенных в разд. 3.

Обозначаем

частные

дифференциалы

рк

и

рт,

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

изменения схе­

мы конструкции, через

?

 

и 6?cxpT, тогда

 

 

 

 

 

 

 

 

 

 

йсхцк

Gn 5G0

4-

 

 

-j-

 

 

 

 

 

 

 

*

 

 

 

 

 

dG0= G0âCK\i.K

\>-gdG0

 

 

äGn

аdGT

 

 

 

 

 

 

 

 

 

 

 

 

ф т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

GodGn dG0\

 

 

 

 

 

 

+ GA*!** + M Go +dG T,туг1

 

 

 

 

 

 

так как последние три члена равны

 

то

 

 

 

)J ^ - l r f G 0 +

 

dG0

G0

[<icxiAK

(1 -ф а) <icxuT] -j-G,

Фк'

 

 

 

 

dGn

 

 

 

 

 

•1 • «

 

отсюдаdGn

 

 

 

 

 

+

“ )

rf'Go,

dG0

J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1—

 

 

___ ФхМ-к +

(1 "Г а)

^CxM'T

 

 

 

 

 

 

 

 

Gn

Фк'

 

 

Ф т

 

 

 

 

 

 

 

 

 

 

dGn

+ (1 + а) dGn

 

— [+r + О + афт]

 

 

Учитывая выражения

(1.6),

(2.4), (2. 5), находим

 

(2:14)

 

dGn

 

 

1

Рф — £)

 

 

 

 

 

 

 

 

«о2

G , і.н (1 +

[ФхІ-’-к “Ь ( 1 +

а ) ^схі1т.

82

Интегрируем левую часть

 

О

Q А

G

о dOp

GQ_

XGp/Gp

 

 

 

 

 

 

 

 

 

fJ

 

(jQ2 _

1 +

AGQ/GQ

 

 

 

 

G 0

 

 

 

 

 

 

 

Интеграл правой части с учетом формулы ('2. 5) будет

 

 

 

 

7Г" / р [АсхР-к + (1 +

а .^ с х !Лт] •

 

 

Обрзначим

 

ао

 

 

 

 

 

(2.

15)

SG 0— / р [ДСХ!Ѵ “Ы

1 -j- а ) Дсхрт],

тогда интеграл выражения (2. 14) получит вид:

 

 

 

 

ДОр __

 

5G0

 

(2. 16)

 

 

 

 

Gp

 

1—

йГ0

 

При малых значениях öG0:

 

 

 

 

 

 

 

 

 

 

AG0

8G0.

 

 

 

 

 

 

 

Gp

 

 

 

 

В этих случаях при Агхрт = 0 из формулы (2. 15)

(2.17)

 

 

 

 

AOp~/pAGK.

 

 

Из последней формулы особенно ясно, что величина fp оправды­ вает свое название коэффициента роста. Выражение (2. 17) яв­ ляется точным при AGK— И).

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

ного веса будет

 

AG0=

д„ HG0+

ДсхО0.

(2. 18)

Для повышения

точности следует вначале определить ACXG0 по

формуле

(2. 16), Jприр

этом' J р

нужно найти новое значение

 

 

 

£ г ,— , ^ Go + A CXG Q

£ \ \ -f

 

 

 

 

 

г

 

— У р Т

А / р і

затем по формуле (2. 7) определяется значение

 

 

 

AII. Ä

— fp

(1 “ЬРф) А^и.н-

Возможен и обратный процесс, т. е. вначале определяется Ап. HG0,

а затем

fp

и ACXG0.

 

 

 

 

 

 

 

 

 

 

 

 

 

Значение Дсхрт, входящее в формулу (2. 15), не является пер­ вичной, заданной величиной. Первичной величиной, вытекающей из тех или иных конструктивных изменений является величина

83

ACT. Найдем связь между АсхРт и Дст. Представим выражение (1.61) с учетом формулы (1. 43) в виде

о

Частный дифференциал этого выражения по Сто будет

X

^<4Ча= (^3о) j"

о

и,следовательно,

(2.19)

°0

Это выражение будет справедливым при соблюдении следую­ щих положений:

1)если при исследуемых конструктивных изменениях закон изменения скоростей не изменяется;

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

3)изменение баллистического коэффициента Дст по числам М подчиняется тому же закону, что и ст, т. е.

л’= Ч л+1)'

Первые два положения выполняются сами собой и вполне реали­ зуемы. Третье положение является допущением. Оно справедли­ во или почти справедливо, если удлинение носовой частя корпу­ са А „> 2 и ACT соответствует хорошо обтекаемому телу. Если третье положение выполняется в недостаточной степени, то тог­ да вместо формулы (2. 191 следует применить формулу, вытека­ ющую из выражения для рта

 

 

 

 

X

 

 

 

 

 

Acxl4=A.HTa=A°o j* у ( Лі + ")< 7 dt’

(2.20)

где

Ai

и

В 1

о

 

 

 

должны соответствовать закономерности изменения

ACT по М.

4.1. Выбор формы носовой части корпуса

Произведем на численном примере сравнение носовых частей корпуса с конической и параболической образующими. В каче­ стве исходных данных принимаем: диаметр корпуса d = 0,6 м, удлинение носовой части Ан = 3,0; /„=1,8 м. Полетный вес, бал­ листический коэффициент и значение рТа берем полученные в

84

гл. I разд. 7, коэффициент роста берем, полученный в примере разд. 3.2, т. е.

G0= 1865 кгс, з0= 0,667 ■ ІО-4, ата=0,1769

(без учета индуктивного сопротивления), /р = 2,58.

Сравнение носовых частей с различными образующими дол­ жно производиться при компоновочной или эксплуатационной равноценности сравниваемых вариантов. В данном случае про­ изводим сравнение при одинаковой длине и одинаковом объеме сравниваемых вариантов носовых частей. На рис. 2. 1 приведены

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

Равенство объемов и длин приводит к тому, что площади бо­ ковых поверхностей рассматриваемых носовых частей отличают­ ся на 1,3%. Поэтому разностью в сопротивлении трения и весе можно пренебречь. Учитываем только влияние волнового сопро­ тивления.

Коэффициент волнового сопротивления для конической носо­ вой части определяем по формуле Краснова Н. Ф. [30]

с ,в.к = (0,0016 + 0,002М-2)[ф7, •

где ßK — угол полураствора конуса в градусах, который в рас­ сматриваемом случае равен 13,42°. Значение числа М опреде­ ляем, исходя из средней скорости 1/ср = 1071 м/с и средней высо­

ты полета /г<щ=30300 м (см. разд.

7 гл. I), для

которых М

ср =

= 3,52. Следовательно, схв.к = 0,144.

Коэффициент

волнового

со­

противления для параболической носовой части определяем по рис. 122 из работы Краснова Н. Ф. [30], согласно которой при

7Н. ч = 3 и М = 3,52

сх

в. п = 0,086.

 

Разница в баллистических коэффициентах

Дз

0,144 — о ,086_ _

0,0877- ю - 4.

 

1865

4

 

 

85

Переходя с помощью формулы (1.43) от М ср = 3,52 к М0=2,5, находим

д3()=0,1025-10-4.

По формуле (2. 19)

Следовательно, по формуле

(2. 15), принимая как в разд. 7 гл. 1

а = 0,15,

5G0= 2,58-1,15-0,0271 =

0,0803,

 

по формуле (2. 16)

= 0,0872.

 

 

 

Go

 

параболи­

Это'значительная величина — 8,7%

от G0. Применяя

ческую

носовую часть, стартовый

вес

анализируемой

ступени

БЛА будет на 163 кгс меньше, чем при конической носовой части. Следует заметить, что если бы сравнение делалось без усло­ вия постоянства объема носовых частей, и начало конической носовой части совпадало бы с началом параболической носовой части, то при конической носовой части стартовый вес был бы на 2% меньше, чем при параболической носовой части. Однако при этом объем конической части был бы на 38% меньше объема

параболической носовой части.

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

нения головок

самонаведения

форма носовой части корпуса

в значительной

степени будет

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

свойствами оболочки.

5. ВЛИЯНИЕ УДЕЛЬНОГО ИМПУЛЬСА ДВИГАТЕЛЯ

НА ПОЛЕТНЫЙ ВЕС

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

Для выявления влияния величины A I на относительный вес топлива рт воспользуемся выражением для цт (1.56). Полагаем, что в полете значение /= /cp = const. Для ракетных двигателей это допустимо, для ВРД это допущение может дать заметную ошибку. ОбЬзначаем

86

(см. выражения (1.52), (1.53) и (1.54)], тогда формулу (1.56) можно представить в виде

—а jAT= 1 — е 1.

Дифференцируя это выражение, находим

і+м

PfXT

 

Л Л = (

d l =

 

~дГ

Делаем подстановку х = 1//, находим

І + М - J

ае

/2

d l.

А/рт= е ~ ах° — е

или

 

Z+а Д / ^ е I

 

L

_ е /

'

1 + Л / / Д

 

Учитывая,что

 

 

а

Г

а

I

М/І

\

 

 

 

 

 

 

 

 

 

 

е

1= 1 — рт,

 

 

 

 

 

получаем

Д /Рт

(1

 

Н-Т

'І +ММII

1

 

(2. 21)

Формулу

(2. 15)

теперь

можно

 

 

Ц

 

 

 

 

представить в расширенном

виде

ЬО0

 

 

а)(ДаРта +

Д/Рт)]-

(2- 22)

 

—/ р [ДсхІ^к + ІІ +

Входящие в это выражение величины определяются по форму­ лам (2.20) и (2.21). Величина AG0/G0 определяется по форму­ ле (2. 16).

\

5.1. Оптимальное давление в камере сгорания РДТТ

Для иллюстрации применения формул (2.21) и (2.22) при­ водим методику и численный расчет определения оптимального давления в камере сгорания РДТТ на БЛА низких высот. Ис­ ходные данные при давлении в камере сгорания рк = 60 кгс/см2 принимаем следующие: рт = 0,3; адв = 0,25. Удельный импульс определяем по работе {27]

при рк = 60 кгс/с /= 219,5.

Вес реактивного сопла, согласно работе [3], Gc = 2,5- 10~4/s= 2,5- 10~4GT/,

где / s суммарный импульс реактивной силы.

87

Следовательно,

относительный вес сопла

 

при рк = 60 кгс/см2

 

 

G ,

=2,5- ІО-4-219,5 =

0,055.

 

 

 

G T

 

Вес камеры сгорания

Щ. = а„ — а =0,195.

 

БЛА низких высот,

Ввиду относительной

длительности полета

а следовательно, сравнительно

небольшой

тяговооруженности,

 

 

 

 

 

влиянием

свободного простран­

 

 

 

 

 

ства двигателя на изменение веса

 

 

 

 

 

двигателя

можно

пренебречь.

 

 

 

 

 

Так как вес камеры сгорания про­

 

 

 

 

 

порционален

номинальному дав­

 

 

 

 

 

лению в камере сгорания, то

 

 

 

 

 

 

 

 

0,195

Рі +

2,5-10- V ,

20 W

60

 

80рК)кгс/смг

 

 

60

 

Рис.

2 .2.

 

32,5/ѵ-

или, учитывая выражение для /,

« д в =

 

1730

578

•10-4.

 

 

 

 

 

 

Р

рк —

 

 

 

 

Принимая за исходное значение к

60 кгс/см2,

 

При этом

А аДв=^32, 5Д/7К +

28,8

•ІО-4.

 

 

 

Acxftc"

AG.

 

:A < V T

 

 

 

По формуле (2. 21)

 

 

 

 

 

 

 

 

(

 

— “ І'Л

 

 

 

 

 

A;!XT= 0 ,7 U - 0 ,7 1+а/// ) .

 

будет при

Оптимальное значение рк по весовому

критерию

(Арк. т)тіп, где

 

 

Аіхк.т= А схРк+

Г,25Д/іхт.

 

 

\

Соответствующие вычисления приведены в табл. 2. 4, график за­ висимости Дрк. т от рк приведены на рис. 2.2. Как видно из гра­ фика, оптимальное значение

Ак.опт^Зб кгс/см2.

88

Т а б л и ц а 2. 4

Р к , кгс/см*

20

40

60

80

Ар

 

- 4 0

— 20

0

20

А«яв

 

— 0,136

— 0,067

0

0,066

(Д///)/(1 +

 

196,4

213,7

219,5

222,3

А///)

— 0,1178

— 0,0271

0

0,0126

/

А/щ

 

0,0290

0,00658

0

— 0,00342

 

 

 

 

 

(1 + а ) Д/р.т

0,03625

0,00822

0

— 0,00428

АСх^к

 

— 0,0408

— 0,0201

0

0,0198

Д[АКТ

 

— 0,0045

— 0,0119

0

0,0155

Заметим, что значение рк, 0Пт будет повышаться приуменьшении относительного веса конструкции и времени работы двигателя. На Рк. опт существенное влияние может оказывать рецептура топлива и высота полета.

6. СРАВНИТЕЛЬНЫ Й ЭКОНОМИЧЕСКИЙ АНАЛИЗ

КОНСТРУКТИВНО-ТЕХНОЛОГИЧЕСКИХ ВАРИАНТОВ

АГРЕГАТОВ И ДЕТАЛЕЙ

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

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

ксущественной разнице в затратах на элементы стартового и

эксплуатационного комплексов.

Представляем производственные затраты на БЛА без полез­ ной нагрузки в соответствии с формулой (2. 1)

Q0= M ° 0—'° я— GmH)v“ + a A 4

где Qо— производственные затраты на ЛА без полезной нагруз­ ки; 0Д — вес анализируемой детали или агрегата; Gn. н— вес полезной нагрузки.

В приведенном выражении первое слагаемое правой части —• затраты на ЛА без анализируемой детали и полезной нагрузки.

Производственные затраты на ЛА будут изменяться вслед­ ствие влияния следующих факторов: изменения веса анализи­ руемой детали, изменения полетного веса, вызванного изменени-

89

Соседние файлы в папке книги из ГПНТБ