Опорные траектории и управляющие функции для вывода ударного летательного аппарата в точку начала атаки Текст научной статьи по специальности «Механика и машиностроение»

УДК 533.6.013

ГРНТИ 30.17.53

ОПОРНЫЕ ТРАЕКТОРИИ И УПРАВЛЯЮЩИЕ ФУНКЦИИ ДЛЯ ВЫВОДА УДАРНОГО ЛЕТАТЕЛЬНОГО АППАРАТА В ТОЧКУ НАЧАЛА АТАКИ

Э.Г. КОКОТЕК

ВУНЦВВС «ВВА имени профессора Н.Е. Жуковского и ЮА. Гагарина» (г. Воронеж)

О.В. АБРАМОВ

ВУНЦ ВВС «ВВА имени профессора Н.Е. Жуковского и ЮА. Гагарина» (г. Воронеж)

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

Ключевые слова: ударный летательный аппарат; точка начала атаки; опорная траектория; управляющие функции.

SUPPORTING TRAJECTORIES AND CONTROL FUNCTIONS FOR OUTPUT

OF IMPACT AIRCRAFT IN THE ATTACK START POINT

E.G. KOKOTEK

MESC AF "N.E. Zhukovsky and YA. Gagarin Air Force Academy" (Voronezh)

O.V. ABRAMOV

MESC AF "N.E. Zhukovsky and YA. Gagarin Air Force Academy" (Voronezh)

The article suggests close to natural support trajectories and, using the known equations of motion of the center of mass of an aircraft, expressions for control functions are output, which makes it possible to create an autonomous guidance system for an impact aircraft at the final stage of a combat flight.

Keywords: impact aircraft, attack start point, reference path, control functions.

Под ударным будем понимать летательный аппарат (самолет, крылатая ракета, управляемая авиационная бомба, беспилотный летательный аппарат), предназначенный для поражения наземных объектов и живой силы противника.

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

Анализ характеристик и возможностей систем управления и наведения ЛА на цель показывает, что достоинства и количество применяемых радио- и оптико-электронных средств являются и их недостатками. Поскольку данные средства имеют свои приемные тракты полезных сигналов, они же становятся и потенциальными объектами не только огневого, но и радиоэлектронного поражения [1]. Поэтому вместе с совершенствованием систем наведения активно разрабатываются и различные способы

противодействия системам управления и наведения (как активные, так и пассивные). На повестку дня встает вопрос создания автономных, помехозащищенных систем наведения.

В работе [2] была предложена методика решения вариационных задач динамики полета с помощью ЭВМ, основанная на использовании прямого метода Ритца (в последующем Ритца-Тараненко), который позволяет решать вариационные задачи с заданными граничными условиями и ограничениями на фазовые координаты и управляющие функции. В последующем этот метод широко использовался (и используется) для решения задач динамики полета при атаке воздушных и наземных целей. Но для этого требуются компьютеры с большой оперативной памятью и быстродействием. Решать такие задачи на борту ЛА затруднительно.

Мы же ставим перед собой задачу - найти простые, но близкие к естественным, опорные траектории и получить управляющие функции, которые позволили бы решать эту задачу на борту ЛА и тем самым обеспечить возможность создания автономной системы наведения ударного ЛА на цель, а также обеспечить возможность пространственно-временного моделирования авиационного удара на этапе его подготовки.

Будем предполагать, что после вывода ударного ЛА в район цели известны его фазовые координаты и координаты ТНА. Расчетная схема этого этапа боевого полета показана на рисунке 1.

V

Уътв

Рисунок 1 - Расчетная схема заключительного этапа боевого полета

Имеем при t = 0:

- фазовые координаты ЛА - Xog, Yog = H, Zog, V = Vo, в = во, ф = фo;

- координаты ТНА (или цели) - XgK, YgK, ZgK.

Здесь V - скорость ЛА; в - угол наклона траектории; ф - угол поворота в горизонтальной плоскости (ф > 0 - при развороте влево).

Пусть траектория в вертикальной плоскости парабола (близка к баллистической кривой), траектория в горизонтальной плоскости также парабола (что вполне допустимо при разворотах на угол р < 90°), т.е. можем записать:

Y = а-(х -X )2 + b-(х -X )

g 1 V g og) \ g og>

+ c;

(1)

Z = а2 -(х -X )2.

g 2 V g og>

(2)

Прежде всего, необходимо определить коэффициенты а1, аг, Ь, с.

Начнем с коэффициента а2. Поскольку нам известны Xgк и ZgK, то

ZgK = а2 -(xgK - Xog)2,

откуда

Z

а2 =

gx

(Xgx - Xog)2 '

(3)

Теперь обратимся к уравнению (1). Если записать его для начальных условий, то получим

С = Yog = H.

(4)

Для определения коэффициентов а1 и Ь нам потребуются известные из динамики полета [3] кинематические уравнения связи:

X dt

= V - cos в cos р;

(5)

dt

dZg dt

= V - sine;

= -V-cosesinр.

(6) (7)

В уравнении (7) знак «-», т.к. ф > 0 при развороте влево. Из (l)

dY dX dX dX ( ч

g = 2aX—g + b—g = —— (2а1 Xg + b).

dt

g

dt dt dt

Из (2)

dZ dX

= 2a2 X —g dt g dt

С учетом (5) и (6)

V ■ sine = V ■ cos в cos ^(2^ Xg + b).

Запишем последнее выражение при t = 0:

Va ■ sin в0 = Va ■ cos в0 cos 0(2aiXog + b) или, сокращая на Vo и, приняв ф = = 0, X0g = 0, запишем

sin#o = cos 6jb.

Откуда

sin 0o

b =-f = tg^o

cose„

(8)

Неизвестным остается ai. Запишем уравнение (1) при t = tK

YK = ai ■ X2k + b ■ XgK + H

и с учетом (8), получим

a1 =

Y^ - XgKtg6>Q - H

x2„

(9)

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

Теперь, имея опорные траектории, необходимо заставить ЛА двигаться по этим траекториям, т.е. необходимо найти закон управления.

Имеем известные [3] уравнения движения центра масс ЛА:

^^ = g(nX - sine) dt ёУ X J

V— = g(nY cos у - cos в) dt

т/ Qdv

V cose— = -g ■ nY sin у dt

(10)

P - X Y

где -= nX - продольная перегрузка; — = щ - нормальная перегрузка; у - угол

G G

крена.

Рисунок 2 - Опорные траектории для различных начальных и конечных условий

На этапе атаки можно считать, что Р = X, тогда пх = 0 и выражения (10) приобретают вид:

■ п. — = - g э1пв.

аХ

(11)

V — = g(п¥ соэ у - соэ в);

ах

(12)

V соэв-= -g ■ пг Э1П у,

ах

(13)

здесь пг, у - управляющие функции.

Уравнение (11) позволяет находить скорость по траектории движения, а (12) и (13) - получить выражения для управляющих функций (потребных нормальной перегрузки и угла крена) в самом общем случае движения. Из(12)найдем пт

т/ав а

V — = g■ пг соэу-g■ соэв ах

э1

и

или

или

1ГсЮ

g■ пг cosу =¥--Vg■ cosO

dt

V--v g ■ соб 0

dt 6

g■соб у

(14)

Из (13) получим у

Из (12)

делим (12) на (13)

ту а йд

g ■ пг sin у = -V соб 0—.

dt

g ■ пг соб у =¥--V g ■ соб0 ;

dt

g ■ пт Бт у

V соб0

dt

g■ щ собу vd0+ g■ соб0 dt

Сокращая, получим

V соб0

tgу =

йд йх

й0 а V--v я ■ соб 0

Л

(15)

„ й0 „ йт

Пусть — = ов - угловая скорость в вертикальном плоскости, а — = ог - угловая dt Л

скорость в горизонтальной плоскости.

Тогда

У-Ъ

V ■юв + я ■ соб 0 я

+ соб0

Я■соб у

соб у

(16)

Щу = -

V соб0Ог V ■юв + я ■ соб 0

V соб0(

соб0(

V ав +

Я■соб 0

V

О +

Я■соб 0

V

(17)

щ =

щ =

как видно они зависят от V, в, ов, сог, т.е. от вида траектории. Ранее записали

dY dX * ч

—— = —- (2a1 X + b) dt dt V 1 g '

или с учетом кинематических уравнении связи

V ■ sine = V ■ cosecos©{la1 Xg + b). Далее, деля обе части уравнения на cos в и сокращая на V, получим

tge = cos ©(2a1 Xg + b). Найдем производную по времени

d (tge) = —■ — = (- sin©)- — (2a X g + b)+ cos© dC cos2в dt ' dty 1 g '

^ dX ^ 2a, —-

1 dt

v J

или

de = cos 2ecos©(2a1V ■ cos в cos©)- cos 2esin©- — (2a1 Xg + b) dt dt

g + b) = ®В

Ранее записали

dZg dX

—- = -V ■ cosesin© или 2a2 X —- = -V ■ cosesin ©

dt ' 2 g dt

или

2a2 XgV ■ cosecos© = -V ■ cosesin©.

Далее, сокращая, можно окончательно получить

tg© = -2a2 X- .

Найдем производную по времени

d

1 d©

^(tg©)= 2 ,

dt cos © dt

dX

= -2a2--- = -2a2 V ■ cos в cos ©

dt

или

d© dt

= -2a2V ■ cos в cos3 © = ^

(18)

Вернемся к с

с =

2a,y ■ cos30cos2 (p-{2a1 Xg + bj-cos20sin((-2a2V■ cos0cos3 (j = = 2axV■ cos30cos2 ( + (2a:Xg + bj-V2a2 ■ cos30sin( cos3 ( = = V ■ cos3o(2ax ■ cos 2 ( + (2a1Xg + bj■ a2 ■ sin 2( cos 2 (( = = V ■ cos30 cos2 ((2a1 + (2a1 Xg + bj a2 ■ sin 2( j.

Таким образом

с = V

V ■ cos30 cos2 ((2a1 + (2a1 Xg + b\ a2 ■ sin 2( j

(19)

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

СПИСОК ЛИТЕРАТУРЫ

1. Сырокваш С.М., Мехеда В.И. Системы управления и наведения крылатых ракет и перспективы противодействия // Наука и военная безопасность. 2008. № 3. С. 60-64.

2. Тараненко В.Т. Опыт применения методов Ритца, Пуанкаре и Ляпунова к решению задач динамики полета: Монография. М.: ВВИА им. Н.Е. Жуковского, 1968. 140 с.

3. Нелюбов А.И., Новад А.А. Динамика полета боевых летательных аппаратов: Учебник. М.: ВВИА им. Н.Е. Жуковского, 1992. 439 с.

REFERENCES

1. Syrokvash S.M., Meheda V.I. Control systems and guidance of cruise missiles and the prospects for counteraction // Science and military security. 2008. № 3. P. 60-64.

2. Taranenko V.T. Experience in applying the methods of Ritz, Poincare and Lyapunov to solving flight dynamics problems: Monograph. M.: N.E. Zhukovsky AFEA, 1968. 140 p.

3. Nelyubov A.I., Novad A.A. Dynamics of the flight of combat aircraft: Manual. М.: N.E. Zhukovsky AFEA, 1992. 439 p.

© Кокотек Э.Г., Абрамов О.В., 2017

Кокотек Эдуард Генрихович, кандидат военных наук, доцент, профессор кафедры общепрофессиональных дисциплин, Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж), Россия, 394064, г. Воронеж, ул. Старых Большевиков, 54А, [email protected]

Абрамов Олег Васильевич, доктор технических наук, доцент, профессор кафедры общепрофессиональных дисциплин, Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж), Россия, 394064, г. Воронеж, ул. Старых Большевиков, 54А, [email protected]