Файл: Дракин, И. И. Основы проектирования беспилотных летательных аппаратов с учетом экономической эффективности.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 20.10.2024

Просмотров: 116

Скачиваний: 2

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Представление траектории в виде алгебраической функции позволяет при вычислениях ограничиться малым количеством рассчитываемых точек траектории (4—6). При выражении тра­ ектории комплексом координат ее точек необходимо иметь малые расстояния между точками, так как в процессе определения аэро­ динамических характеристик необходимо определять первую (dhjclx) и вторую (d2hjdx2) производные координат. Это требует значительного увеличения количества рассчитываемых точек. Малое количество параметров и малое количество расчетных то­ чек весьма существенно уменьшают объем расчетной работы при оптимизации траектории параметрическим методом.

Если принять один варьируемый параметр, то для определе­ ния оптимальной траектории следует определить расход топли­ ва для 4 траекторий. В случае принятия в уравнении (4. 2) двух варьируемых параметров, придется определить расход топлива уже для 16 траекторий. Как видим, объем расчетов увеличи­ вается в 4 раза.

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

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

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

h%

при втором

приближении

можно представить в виде /?2 —

 

 

(4.3)

где

hi

 

 

из расчета первогоприближения

— высота, полученнаяі

 

 

 

 

 

(4.4)

 

 

 

/ii = fLo+'^7Ux‘ + t1x?,+1.

Величина

Ah

<=і

 

 

параметра |2,

 

является функцией координаты х и

 

 

 

Ah=

ср(£2,

л-),

(4.5)

 

 

 

 

 

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

ном и конечном значениях

х

 

=0.

(4.6)

 

хк)

Д/г0 = <р(І , 0) =0, ДЛК= ф(|2,

 

2

 

 

 

 

Дальнейшая задача состоит в подборе подходящей функции для Ah, удовлетворяющей условиям (4. 6) и позволяющей воз­ можно в большей степени учесть факторы, влияющие на опти-

150


мальность траектории. Одной из таких функций хможет являться функция

где

s> m \

Л/г = 53 (лгк ■— x)sx ,

(4.7)

 

увеличение s понижает перегрузки

в конце полета,

т соответствует (4. 4). Эта функция учитывает тот факт, что для уменьшения расхода топлива, идущего на преодоление аэроди­ намического сопротивления, следует возможно быстрее перехо­ дить на большие высоты полета, т. е. делать более крутой старт. Конечно, при этом возникает и отрицательный фактор — увели­ чение поперечных перегрузок, что приводит к возрастанию ин­ дуктивного сопротивления. Функция (4. 7) мала по величине в. конце полета, где оптимальная форма траектории в достаточной степени определяется в результате первой оптимизации или соот­ ветствующими ограничениями по перегрузкам.

Возможны для величины Ah и другие функции, например,

р

Л/г = £2(хк- х ) ( 1 - е _ ^ ) .

(4.8)

От величины

зависит положение

Ahmax,

например, при р=10„

 

величина ЛЛтах будет при х/хк«0,22.

1.1. Пример определения оптимальной траектории

Для иллюстрации оптимизации траектории параметрическим методом приводим пример расчета в сокращенном виде. В ка­ честве исходных данных принимаем следующие произвольные значения: х = 0, хк = 50 км, А0= 0, /гк = 25 км; Ѵо=600 м/с, Ѵк= = 800 м/с, т= 80 с, ао=10“ 4 м2/кгс, /ср=>240, ограничений на тра­ екторию, кроме указанных не накладывается.

Так как траектория ограничена только положением началь­ ной и конечной точек полета, то для оптимизации первого при­

ближения, учитывая,

что вследствие того, что /г0= 0, значение

Хг= 0, уравнение (4.2)

принимает вид

где / 1

 

h = y„iXJr ’iI*2,

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

(хк =

50 км, /гк= 2 5 км),

ул = 0,5 — 50ф.

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

h = (0,5 — 50ф) х ф- фх2,

здесь значение х и h в км.

Задаемся значениями h: 0, —0,01, —0,02, —0,03. Кривые, со­ ответствующие получающимся при этом уравнениям, приведены на рис. 4. 2. Пользуясь методами определения расхода топлива, закона изменения скоростей и аэродинамического сопротивле­

151


ния,.изложенными в гл. 1, найдены относительные веса топлива при полете по каждой траектории с указанными значениями На рис. 4. 3 приведен график зависимости относительного веса

топлива рт от величины параметра

кривая

А.

При определе­

 

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

Как видно из графика на рис. 4.3, оптимальная траектория соответствует | 4= —0,015, при этом значение рт достигает ми­ нимального значения. Уравнение этой траектории

А = 1,25л:— 0,015л2.

По сравнению с прямолинейной траекторией, при полете по оп­ тимальной траектории расход топлива на 13% меньше.

Переходим к последующей оптимизации траектории. Учиты­ вая, что в найденном уравнении траектории величина т = 1, при­ нимаем в выражении (4. 7) s = 3; некоторое увеличение s по срав­ нению с (т + 1) сделано в целях предотвращения роста перегру­ зок в конце полета. Тогда

ДЛ = $2(50 — x f x .

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

h = 1,25х — 0,015л2-|-^2 (50 — х)3х.

Задаемся значениями £2 в интервале 0—16-10—6. Заметим, что £г = 0 соответствует найденному выше уравнению оптимальной траектории первого приближения. Для новых траекторий опре-

152

деляем относительный вес топлива и затем находим графически значение £2, соответствующее минимуму [іт-

Зависимость величины рт от |г приведена на рис. 4.3, см. кри­ вую В. Как видно из графика, оптимальное значение

 

 

L n„T= 6 , 5 ■ ІО -6,

 

 

 

h.KM

 

 

при этом

 

 

 

 

 

 

 

 

30

е = = Г —

-

 

 

IK min =

0,424.

 

 

 

 

 

 

Эта

величина

на

 

4%

меньше

20

L

 

 

Ю

 

 

значения

рт,

полученного

в

 

I

 

результате

оптимизации

пер­

 

 

 

вого приближения.

 

 

 

 

 

 

10 20 30 U0

50 Х,КМ

А

На рис. 4.4 приведены опти­

 

 

 

мальные траектории:

кривая

 

 

 

 

— траектория первого

при­

 

Рис. 4. 4.

 

ближения, кривая

В

— траекто­

 

 

 

 

 

рия

второго

приближения.

 

 

 

Дальнейшее

уточнение оптимизации не имеет смысла, так как

возможная/ ,

экономия топлива будет меньше возможных ошибок

к

величине рт, вызванных

неточностями исходных данных (оо,

 

Т,

Q и др).

 

 

 

 

 

 

 

 

 

 

2. ОПТИМАЛЬНАЯ КОНЕЧНАЯ СКОРОСТЬ МАНЕВРЕННЫХ БЛА

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

где /ѵ

 

‘■Ук'

су

йѴІ

 

 

а к

Рк

 

 

■ удельная нагрузка на крыло в конце полета;

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

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

Увеличение конечной скорости полета ведет к увеличению расхода топлива [см., например, формулы (1.62), (1.69)], сле­ довательно, при этом увеличивается и полетный вес. Однако при увеличении конечной скорости уменьшается потребная площадь крыльев, что ведет к уменьшению их веса и аэродинамического сопротивления во всем полете, а следовательно, к уменьшению пблетного веса. Очевидно, в каждом случае существует опти­

153


малыше значение конечной скорости, при котором полетный вес будет минимальным.

•Варьируя конечную скорость полета, предполагается, что как схема всего БЛА, так и его агрегатов и, в частности, крыльев и двигательной установки, остаются неизменными. Оптимизацию конечной скорости проводим вначале на базе весового крите­ рия — GoПоэтому условием оптимальности Ск будет

- ^ - = 0 .

(4.9)

dVк

Для многоступенчатого БЛА под весом G0 в данном случае следует понимать вес последней маневрирующей на высоте сту­ пени. Действительно, согласно формуле (1.2), вес предшествую­ щих ступеней будет практически пропорционален весу послед­ ней ступени, поэтому вполне возможно оптимизировать величи­ ну Ѵк по весу последней ступени.

Учитывая, что в данном случае G0 зависит как от ркр, так и

от рт, то

ÖGQ

d\L

Ö G Q

d V K

dGg

âu.l:

p

кр

Ф-Г

d V K

 

d V K

 

 

d[XT

Из формулы (1.1) следует

 

 

 

10дО

 

 

 

 

 

Go1 )

Од

[1—(Ак+ш)]2

I

(fXK-f[AT)

о

 

О „я

 

°°('+

или учитывая, что

 

(‘ +ѣ)

<?!Х

[1—0*к+М-т)]2

 

1—0%+ш)

то

dGp _

 

<Ѵг

адв)

 

GQ(1 +

 

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

дщ

 

. 1— (ftc + щ)

^Н-кр

 

(1

d\

0.

 

 

 

лт

 

lv~K

 

 

dv7

 

 

 

 

 

 

Сопоставляя это уравнение с уравнением (3. 11), а также с выражением (3. 10), нетрудно видеть, что полученное уравне­ ние можно применить для оптимизации Ѵк и при экономическом

154