Файл: Дракин, И. И. Основы проектирования беспилотных летательных аппаратов с учетом экономической эффективности.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