книги из ГПНТБ / Дракин, И. И. Основы проектирования беспилотных летательных аппаратов с учетом экономической эффективности
.pdfрактеристик, пропорциональны полетному весу 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