Файл: Дракин, И. И. Основы проектирования беспилотных летательных аппаратов с учетом экономической эффективности.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 20.10.2024
Просмотров: 117
Скачиваний: 2
3.3. Пример оптимизации параметров корпуса
Вкачестве основных исходных данных для примера исполь зуем произвольные характеристики БЛА класса ЗВ, близкие к характеристикам, приведенным в статье [14] для франко-англий
ского БЛА «Кроталь».
Стартовый |
вес в к г с ..................................... |
км |
75 |
8 |
|
|
|||
Дальность полета у земли в |
. |
600 |
|
|
|||||
Средняя скорость |
полета в м/с . |
|
|
||||||
Конечная скорость полета в м/с . |
. |
500 |
|
|
|||||
Двигатель |
РДТТ, |
удельный3 |
импульс |
|
240 |
|
|
||
[56] в |
с ............................................................................... |
|
|
|
|
|
|
|
|
Объем корпуса в см .......................................... |
|
м |
|
0,06 |
|
|
|||
Объем отсеков с |
оборудованием в |
|
0,032 |
. |
2.1 0 |
||||
|
|
|
|
|
3 |
|
|
||
Удельный вес конструкции отсеков оборудования в кгс/м |
|||||||||
Удельный вес конструкции двигателя диаметром |
1 м в кгс/м 50 |
||||||||
Относительные |
веса конструкций: рНп = |
0,05, |
7[хуп= |
0,02; |
ß,j,= |
||||
= 0,25, адв=0,4, соотношение удельных стоимостей топлива и кон |
|||||||||
струкции принимается |
(по разд. 7. 1, гл. II), |
<т/< к =0,0835 |
|
||||||
|
|
|
|
|
|
7 |
|
|
|
Оптимальное удлинение носовой части
Определяем входящие в формулу (3.82) величины. Ввиду ко роткой длины полета, причем на постоянной высоте (h m 0), оп ределяем величину шм н по формуле (3.79), ограничиваясь од ним членом суммы.
По формулам (1.64) и (1.67)
|
Ф—ехр |
500 — 600 |
1,043, |
|
|
w М |
9,81-240 |
|
8000 |
1830. |
|
следовательно, 1,043(1 + |
1,25/3,12)0,125-360000 |
|
|||
|
н |
2-240 |
|
600 |
|
Ввиду слабого влияния длины на коэффициент аэродинами ческого трения определяем длину корпуса по статистической формуле для зарубежных маневренных БЛА
і_ і_ /=0,650(5 = 0,65-75s =2 ,74 м.
Среднее число Рейнольдса
Rert |
2,74-600 |
1,1310s. |
|
1,46-10-5 |
|
Определяем Cf, принимая пограничный слой турбулентный, по работе [33]
с/ф= с /тілі=0,0021 -0,77 = 0,00162.
143
По формуле (3.72), ограничиваясь одним членом суммы, анало гично величине wM n, находим
Принимая длину |
да/ = 2,12. |
|
Іф, |
находим |
носовой части |
О |
|
приведенную длину носовой части для случая турбулентного пограничного слоя, согласно работе [21]:
Ф'
Следовательно,
ReH= 1,89 ■ ІО7, Су = 0,0027, г/н= 0,0021.
По формуле (3. 78)
wfH= 2 ,7 5 .
Оптимизацию удлинения носовой части проводим по эконо мическому критерию [см. формулу (3.6')]. Определяем вначале значения величин г и и, входящих в формулу (3. 6')- Для этого предварительно определяем рта по формуле (1.75). Учитывая, что для данного БЛА удлинение носовой части должно быть больше среднестатистического (низкая высота полета), прини маем
0 = 7 -ІО-4, 4 = 0,5, В = 1,25 (см. разд. 4.1 гл. I).
Дальнейшие вычисления приведены в табл. 3. 3.
В табл. 3. 3 получено при экономическом критерии
При весовом критерии
Ѵв= 1 + а = 1,4.
Следовательно,
^= 5 . 6 5 ,
т.е. на 12% больше, чем при экономическом критерии.
Определение диаметра корпуса |
носовой части |
Определяем вначале сх ф при Л4ср = 1,77. Для |
|
при %j= 5,03 согласно формуле (3.77) волновое |
сопротивление |
сх в = 0,0449
Донное сопротивление не учитываем, так как у РДТТ будет сравнительно большая тяговооруженность (рср = 8,15), поэтому диаметр выходного сечения сопла будет близок к диаметру кор пуса. Сопротивление трения определяем по статистическому уд линению. Принимаем Лф= 15.
144
Формулы и замечания
|
J cp ■ |
|
D |
|
(А + В/Мср), |
|
Л1ср = 1 ,7 7 |
|
|
||||||||||||
|
|
|
|
|
|
|
|||||||||||||||
|
пга = |
^ |
- ^ Ѵ срг2 |
, |
|
|
1.043, |
|
|
|
|||||||||||
|
(Хт у |
|
^ср =: 600 |
м/с, |
т = |
13,32 |
с |
|
|
|
|
||||||||||
|
|
|
Ѵк - Ѵ |
|
п |
|
|
К к = 500 м/с, |
|
/ = |
240 с |
|
|||||||||
|
|
|
|
gl |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
[іт= |
1— е |
1’’тѴГ+ |
(Атв |
|
|
|
|
|
|
||||||||
|
|
|
/о = |
|
|
|
|
|
1 + |
Рф |
а) fJ-т] |
|
|
|
|
|
|||||
|
|
|
1— [(*£■ + |
(1 + |
|
|
|
|
|
||||||||||||
|
|
|
|
p-g- = |
рнп1-Т р.уп — 0,07 |
|
|
|
|
|
|
||||||||||
|
|
|
|
3= - |
|
g Лн-Т й(1+ а) |
|
|
|
|
|||||||||||
|
|
|
|
: CaСщі -- |
g |
/орнп |
|
|
|
О |
|
|
|
|
|||||||
|
|
|
|
|
|
|
+ |
|
Снп + |
Сб > |
Сб = |
|
|
|
|
|
|||||
|
|
|
|
|
|
д О п |
|
1 + |
/о1’ М-г а |
|
|
|
|
|
|
||||||
|
|
|
|
|
/ р |
= |
|
|
|
Р ф - С |
|
|
|
|
|
|
|
|
|||
|
|
K |
|
,. / |
|
,,, |
|
, |
|
ѵт9т |
^ |
|
^рт \ .. |
|
|
|
|||||
|
|
|
|
, |
|
Фт |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Ф-к |
|
G,О |
г0 |
|
+, |
a G^ 0 |
^Рт — |
1 |
|
Цнп — |
Q |
а Рт |
t |
|||||||
|
|
^рнп, ~ |
|
|
1 |
|
|||||||||||||||
|
öGn |
|
|
/р |
ö° |
|
|
|
+ |
|
G0 |
|
|
ѵ к, |
|
|
|||||
|
G ? |
|
: |
OOQ |
|
|
|
|
|
dt/Q |
|
«3 |
|
|
|
о |
|
|
|||
|
|
|
r> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
V |
|
|
д |
О |
0 |
|
ѵ к? к |
|
dG0 J |
|
|
|
|||||||
|
|
|
ѵ к = |
0,7, |
ѵ т = |
1 |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
I |
Рк + |
v.tqT |
рт |
|
I Ѵк, |
|
|
|
|||||
|
|
|
и !Ч к — |
/ р |
ѵк?к |
|
|
|
|
||||||||||||
Yi = |
|
Рк |
|
|
|
/ |
|
|
|
|
Р'ДВ |
|
|
|
|
||||||
|
|
|
рнп T “ Рун + р-ф |
|
|
|
|
|
|
||||||||||||
|
г |
|
|
|
|
|
|
|
|
Г + И 1 / / |
|
Г + И \ |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
J / V |
1 + |
ѵк? к / |
|
|||||
|
L Ѵ |
- + |
а + |
( 1 |
+ |
а ) -------- |
|
/ |
- ^ — \ |
||||||||||||
|
|
|
к^к |
|
|
|
|
|
|
|
v Kq K |
|
|
|
|
|
|
|
|||
|
|
|
|
^.опт = |
0,656 X |
|
|
|
|
|
|
|
|||||||||
X |
|
|
|
|
|
|
|
|
|
|
0,37 |
|
|||||||||
|
|
|
|
|
|
|
|
|
W M н |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
g'o6 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
,1——45 W)1j, Yi— — |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
____T)U7 = |
0,6 |
|
|
|
|
|
|
|
|
|
Т а б л и ц а 3.3
Номер формулы Величина
(1.43) |
2-10-4 |
|
(1.75) |
0,261 |
|
(1.62) |
0,213 |
|
(1.60) |
0,453 |
|
( 1.6) |
4,25 |
|
(2. 12) |
-0,518 |
|
(2.9) |
0,0709 |
|
(2.13) |
-0,447 |
|
(2.5) |
2,51 |
|
(2.12') |
-0,087 |
|
(2.9), |
-0,0181 |
|
(2. 11) |
-0,05 |
|
(3 .4'), |
||
разд. |
2.2 |
|
гл. |
II |
0,637 |
(3.5') |
||
( З .б ') |
0,92 |
|
(3.82) |
5,03 |
145
Так как поверхность параболической носовой части
Учитывая найденные выше значения Х*= 5,03, с/ф= 0 ,00162,
получаем
схг?= 4 (^-5,03 + 9,97] 0,00162 = 0,0852.
Следовательно,
схф= 0 ,0449 + 0,0863 = 0,132. По формуле (3. 71), аналогично Wf, находим
®Ф= 172.
Пользуясь формулой (3.75), исходными данными, найденными величинами и учитывая, что
|
«дн = |
2. |
Яб=1> |
|
|
|
находим при экономическом критерии, при -у = 0,92, |
||||||
d = yrb,77- |
ІО-3 — 0,332аГ4. |
|
|
|||
Задаваясь вначале |
|
|
d = |
0,2, последова |
||
произвольным значением |
|
|||||
тельными приближениями находим |
|
|
||||
|
d a |
=0,185 м. |
|
|
||
При весовом критерии, при |
у = |
1,4, учитывая |
|
|||
получим |
Хң.СЛГТ |
5,65, |
|
|
||
й?о п т = 0 - 1 8 3 М - |
|
|
||||
Удлинение корпуса находим по формуле (3. 69) |
|
|||||
X |
4-0,06 |
5,03 |
1---- ^ 0 ,6 |
=15,5. |
||
'ф-опт |
я-0,1853 |
|
|
15 |
|
|
О П Т И М И ЗА Ц И Я |
К И НГ л а в а IV |
Е Т И Ч Е СК И Х Х А Р А К Т ЕР И СТ И К |
ИР А СХ О Д А ТО П Л И В А
1.ОПТИМИЗАЦИЯ ТРАЕКТОРИЙ ПАРАМЕТРИЧЕСКИМ МЕТОДОМ
Траектория полета является одной из основных кинематиче ских характеристик ЛА. В реальном полете траектория опреде ляется методом наведения и законом скоростей. При проектиро вании ЛА необходимо вначале найти желательную траекторию, а затем определять параметры системы управления, которые могли бы обеспечить заданную траекторию, хотя бы прибли женно.
Следует заметить, что точное решение задачи определения траектории должно решаться с учетом метода наведения. Одна ко, существующая аппаратура управления и методы наведения настолько гибки, что позволяют осуществлять, хотя бы прибли женно, возможные для ЛА траектории. Поэтому в данном раз деле излагаются методы оптимизации траекторий без учета мето дов наведения.
Оптимизация траекторий может производиться на основе раз личных критериев. Наиболее распространенными критериями для оптимизации траекторий являются: стартовый вес ЛА, потреб ный вес топлива для полета, ошибки наведения, время полета. При заданной системе управления и определенной компоновоч ной схеме стартовый вес будет определяться в основном расхо дом топлива. При этом предполагается, что закон скоростей по времени задан (он определяется требованиями к начальной, средней и конечной скоростям, а также требованиями к продоль ным и поперечным перегрузкам).
Ошибки наведения зависят от конечной части траектории, поэтому могут быть учтены в виде ограничений на поперечные перегрузки в конце полета. Время полета определяется средней скоростью и дальностью полета, и может быть опре делено из ограничений. Время полета как критерий может пред ставлять интерес при заданном расходе топлива, т. е. при экс плуатационных задачах.
Исходя из приведенных соображений, в качестве критерия оптимальности траектории принимаем расход топлива. Заметим,
147
что этот критерий в данном случае почти эквивалентен экономи ческому критерию.
В процессе эксплуатации летательных аппаратов какого-ни будь типа могут встречаться различные исходные условия для по-,
лета: дальность полета, высота пуска Л А, |
высота в конце поле |
та и др. Очевидно, что при оптимизации |
траектории следует |
выбирать из общего комплекса возможных условий для полета такие условия, которые ведут к максимальному расходу топли ва. Например, задаются дальности полета на большой и на ма лой высоте, причем последняя дальность меньше. В этих слу чаях необходимо путем приближенных сравнительных расчетов выявить расчетные условия для оптимизации.
Оптимизация траектории является вариационной задачей. Вариационным задачам динамики полета посвящено много ли тературы, например, работы [40], [44], [46], [49], [60]. Однако да же простейшие вариационные задачи динамики полета требуют значительной вычислительной работы, так как приводятся к не скольким дифференциальным уравнениям, которые, как прави ло, необходимо решать с помощью ЭЦВМ . В связи с этим развит ряд методов численного решения таких задач.
Применение ЭЦВМ позволило использовать для оптимиза ции кинематических характеристик статистические методы поиска [51], которые задачу оптимизации функции позволяют сводить к оптимизации параметров.
В данном параграфе оптимизация траектории проводится с целью правильного определения потребного запаса топлива, а также с целью разработки метода наведения БЛА. Такая оп тимизация проводится на этапе предэскизного проектирования, когда не требуется высокой точности расчета, так как исходные данные — приближенные. Действительно, на этапе, предшеству ющем оптимизации траектории, конструктивный и динамический облик ЛА еще не выявлен в достаточной мере, следовательно, величина аэродинамического сопротивления может быть опре делена или приближенно или по статистике. Также приближен но могут быть определены весовые характеристики. Следует учи тывать и такие факторы, как нестабильность атмосферных усло вий: плотности воздуха, его температуры, скорости ветра. Поэто му заданные летные характеристики следует рассматривать как среднестатистические.
Вследствие указанных причин, следует по существу опреде лять не оптимальную траекторию, а зону оптимальных траекто рий. В кач'естве примера на рис. 4. 1 для произвольного ЛА при водится оптимальная траектория А и зона оптимальных траек торий между кривыми В я С. Оптимальная траектория соответ ствует номинальным атмосферным условиям и статистическому значению баллистического коэффициента. Кривые В и С соответ ствуют траекториям с расходом топлива на 3% больше, чем при полете по траектории А.
148
Более простыми являются параметрические методы оптими зации траекторий. Эти методы базируются на том, что траекто рии обладают следующими свойствами: непрерывность, отсутст вие особых точек, сравнительно большие радиусы кривизны тра екторий. Одним из вариантов параметрического метода оптимизации траекторий является аппроксимация траектории ломаной кривой. Разбивая горизонтальную ось х на. п участков и задаваясь для каждого значения ду величиной высоты полета
hi, |
задачу оптимизации траек |
h,m\ |
|
|
|
||||
|
|
|
|
hi. |
зо |
|
А |
|
|
тории в вертикальной плоско- |
|
|
і |
|
|||||
сти |
можно |
свести к |
задаче |
|
|
|
1 |
|
|
оптимизации |
параметров |
|
20 |
|
|
|
|
||
При этом для каждой |
комби |
м |
|
|
|
|
|||
|
|
|
п, |
|
|
|
|
|
|
нации параметров определяет- |
|
|
|
|
|
||||
ся расход топлива. |
|
тем |
|
—-А> |
|
|
|||
Чем больше участков |
|
|
|
|
|
|
|||
точнее будет |
оптимизация, но |
|
|
10 20 30 |
00 50 |
х,км |
|||
тем |
больше будет параметров |
|
|
|
|
|
|||
и, следовательно, тем больший |
|
|
Рис. 4. |
I. |
|
||||
объем расчетной работы. Прак |
|
|
|
||||||
тически этот |
параметрический |
|
|
|
|
|
|||
метод оптимизации траекторий |
|
|
|
случае |
осо |
||||
может осуществляться с помощью ЭЦВМ . В этом |
бенно удобны методы статистического поиска, допускающие при приемлемом времени расчета большое количество параметров и большое количество ограничений.
Значительное сокращение объема вычислительной работы мо жет быть достигнуто, если траекторию аппроксимировать подхо
дящим видом уравнения |
|
(4. і; |
|
^ f |
(•*"> |
Хі> /2 ’ • • • > Xm’ ?1> ^2’ ?л)> |
|
|
|
где величины %и %2, . . . , %т определяются на основе исходных данных, величины | ь |г, ■ ■ •, E« — являются варьируемыми пара метрами. Уравнение (4. I) соответствует траектории, располо женной в вертикальной плоскости.
Для траекторий с наклонным (не вертикальным) стартом под ходящим видом уравнения траектории может быть
|
|
^—^0+ l l X |
Ъ * 2 + ••■ + Х т * " 1 + І і Х т |
+ |
Х + |
|
|
|
+ |
Ъ#*+2+ . . . + ІпЛт +11. |
|
(4.2) |
|
Следует заметить, что на этапеп |
предэскизного |
проектирования |
||||
можно ограничиться значением |
= 2 . |
|
|
|||
В случае вертикального старта, когда радиус кривизны в точ |
||||||
ке старта |
г |
= —оо, первый член в уравнении (4.2) |
следует иметь |
|||
|
|
|
|
|
|
1 |
в дробной степени; подходящим значением может являться Хіх2 ■
149