Abstract and keywords
Abstract (English):
The relevance is due to insufficient knowledge and development of movement patterns for unmanned aerial vehicles (UAVs), which include all possible drones and, in particular, various copters in mobile and turbulent atmosphere. Unmanned aerial vehicles are used both in civilian industries and in military affairs [1 – 4]. When operating UAVs in real conditions, their movement is affected by disturbances from air masses [5, 6], which are generally random in nature. At the same time, the operator is faced with inadequate UAV behaviour caused by disturbances, which makes it difficult to control. The aim of the article is to simulate the complex movement of UAVs, mainly quadrocopters, in space to understand the dynamics of UAV flight under the influence of several forces. The models are based on kinematics and dynamics of a material point, in the form of which the UAV is represented, involved in complex relative and translational movements. Kinematic models of drone movement patterns based on coordinate transformation are de-scribed. Planar and spatial trajectories of UAV movement are presented.

Keywords:
unmanned aerial vehicles, movement trajectory, coordinate transformation, movement pattern
Text
Text (PDF): Read Download

Ведение

 

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

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

В то же время в работе не рассматриваются возмущения в виде природных или искусственных аномалий, вызванных воздействием электромагнитного излучения (средств радиоэлектронной борьбы [8, 9]) или терриконов с высоким содержанием металлических примесей, способных изменить траекторию движения дрона, вплоть до его падения, вследствие отказа электронной части. Это вопрос другой темы.

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

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

 

Модели движения беспилотных летательных аппаратов

 

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

Рис. 1. Сложное движение БПЛА в вихревом облаке

Fig. 1. Complex drone movement in a swirling cloud

 

Движение дрона – сложное, его движение по окружности радиуса R относительное, а вращение его плоскости – переносное [12]. Абсолютная скорость дрона определяется выражением:

v=vотн+vпер,                                                         (1)

где vотн=dsdt=Rdφdt – относительная скорость,                                                                            (2)

а vпер=R1ω1 – переносная скорость дрона,                                                                             (3)

где R1 – радиус окружности, описываемой той точкой окружности радиуса R, с которой совпадает в данный момент т. М (расстояние от т. М до оси вращения воздушных масс).

R1=2R+R(1+cosβ)=R(1+cosβ);

а ω1 – модуль угловой скорости плоскости окружности радиуса R:

ω1=ω1=dϕdt и vпер=R(1+cosβ)∙dϕdt.                                    (4)

Модуль переносной скорости: vпер=R1ω1=R1dϕdt.                                      

Вектор vпер направлен по касательной к окружности в сторону вращения тела.

Абсолютная скорость т. М – геометрическая сумма относительной и переносной скоростей v=vотн+vпер. Векторы vотн и vпер взаимно перпендикулярны.

 

Модуль абсолютной скорости найдем как геометрическую сумму относительной и переносной скоростей:

v=vотн2+vпер2 ,                                                                 (5)

тогда подставив выражения (2) и (4) в (5), получим:

v=Rdφdt2+R(1+cosβ)∙dϕdt2=Rdφdt2+(1+cosβ)∙dϕdt2.

 

Угол β определяется угловой скоростью дронов, т.е. β=ωt.

 

Моделирование условий

 

При ω = 10, ω1 = 20, R = 2 и при t = 0, β = 0 зависимость абсолютной скорости дронов от времени и угла поворота приведены на рис. 2 и 3.

 

Рис. 2. Зависимость абсолютной скорости дронов от времени

Fig. 2. Dependence of the absolute speed of drones on time

 

 

Рис. 3. Зависимость абсолютной скорости дронов от угла поворота

Fig. 3. Dependence of the absolute speed of drones on the angle of rotation

 

Моделирование траектория движения

 

Траектория движения разведывательных дронов строится путем преобразования координат при движении плоскости вращения дрона относительно оси О1О2, приняв за неподвижную систему, начало которой совпадает с точкой касания плоскости вращения БПЛА и оси вращения воздушных масс.

Преобразование координат имеет вид [13]:

z'=zx'=xcosψ-ysinψy'=xsinψ+ycosψ – координаты в неподвижной системе отсчета при вращении плоскости дрона,

где z=Rsinφx=R(1+cosφ)y=0 – координаты в неподвижной системе отсчета при отсутствии поворота.

 

Здесь φ=ωt, а ψ=ω1t и

z'=z=Rsinωtx'=R(1+cosωt)∙cosω1∙t-ysinω1∙ty'=R(1+cosωt)∙sinω1∙t+ycosω1∙t.                                       (6)

 

При ω1=π, ω=2π и R=1 и начальном положении дрона в т. А (см. рис. 1) проекции его траекторий на плоскости XY, XZ и YZ приведены на рис. 4 на основе использования (6).

а)                                                                                                          б)

в)

Рис. 4. Проекция траектории дрона на плоскость:

аXZ, бYZ, вXY

Fig. 4. Projection of the drone trajectory onto the plane:

а – XZ, bYZ, cXY

 

Пространственная траектория дрона строится с использованием функции MATHCAD, программа расчета приведена ниже:

ω≔2π   ω1≔π

y0    r≔1

y1trsinω∙t

x1tr∙(1+cosω∙t)∙cosω1∙t-y∙sin⁡ω1∙t

y1tr∙(1+cosω∙t)∙sinω1∙t+y∙cos⁡ω1∙t

Rtx1ty1ty1t

MGreateSpace(R, 0, 2, 200),

где R(t) – вектор параметрических уравнений координат.

Здесь числа 0 и 2 отвечают за нижнюю и верхнюю границы переменной, tgrid = 200 – число линий сетки CreateSpace [14].

Траектории БПЛА в пространстве приведены на рис. 5.

 

Рис. 5. Траектории дрона в пространстве

Fig. 5. Drone trajectories in space

 

Перспективы исследований

 

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

Применение разработанных моделей в сочетании с искусственным интеллектом на безе нейронных сетей позволит избежать сложностей для оператора при управлении БПЛА [15], когда основную часть расчетов по корректировании траектории движения дрона при возмущениях будет производить нейронная сеть, программа функционирования которой должна быть заложена в контроллер управления беспилотников [16], а задача оператора заключалась бы только в управлении перемещением дрона по намеченному пути.

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

 

Выводы

 

1. Получены и изучены кинематические модели движения дронов в пространстве при перемещении воздушных масс с целью выяснения степени сложности траекторий их движений.

2. Для получения моделей сложного движения использованы математические методы пространственного преобразования координат. Представлены модели траектории движения БПЛА типа дронов.

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

References

1. Tuev A.V. Application of Robotics and UAVs in Agriculture. Youth and Science. 2020;4:30.

2. Akhmedov B., Akhmedova G. Using UAV in Collecting Geospatial Data for the Purpose of Solving Engineering Problems in Urban Construction. Danish Scientific Journal. 2021;46-1:41-45.

3. Kharichkin A.I., Rogov K.S., Dranitsyn A.V. The Use of Modern Methods of Aaerophotography Based on UAVs for the Project of Engineering Protection Measures Against Dangerous Geological Processes. Bulletin of Science and Research Centre of Construction. 2021;2(29):123-135.

4. Turik A.A., Miroshnikov V.I., Goncharov S.A. The Use of UAVs by the Parties in the Conduct of Hostilities in the SAR. Bulletin of Scientific Conferences. 2018;5-1(33):112-115.

5. Malkin V.A., Rozhkov I.V., Sanko A.A. Influence of Wind Disturbances and Nonlinearities of Servo Drive on the Contour of Stabilization of the Flight Height of Unmanned Aerial Vehicle. System Analysis and Applied Information Science. 2019;2:23-30.

6. Sanko A.A., Sheinikov A.A. Modelling Wind Disturbances Acting on the Sensitive Element Block of the Navigation System in the UAV Vertical Channel. Aviation Bulletin. 2022;6:68-76.

7. Shakina N.P. Orographic Turbulence: Generation Mechanisms and Forecasting. Hydrometeorological Research and Forecasts. 2019;1(371):25-47.

8. Burlakov D.O. Analysing Methods for Detecting and Controlling Military Drones. Modern Scientific Researches and Innovations. 2022;9(137).

9. Titkov O.S., Mukhametzhanova A.O. Modern Ways of Anti-Unmanned Aerial Vehicle System Development. Aviation Systems. 2022;4:44-53.

10. Wolf D.A., Aleksandrov V.A., Rezkov I.G. Automation of the UAV Pilot’s Behaviour Using a Native Microcontroller. Industrial ACS and Controllers. 2023;3:9-16.

11. Adnastarontsau AA, Adnastarontsava DA, Fiodortsev RV, et al. Algorithm for Control of Unmanned Aerial Vehicles in the Process of Visual Tracking of Objects With a Variable Movement’s Trajectory. Devices and Methods of Measurements. 2021;12(1):46-57.

12. Targ S.M. Theoretical Mechanics - A Short Course. Moscow: Vysshaya shkola; 2002.

13. Rotation Matrix in Three-Dimensional Space [Internet] [cited 2023 Apr 25]. Available from: https://steptosleep.ru/rotation-matrix-in-three-dimensional-simple/

14. Akishin B.A., Erkenov N.Kh. Applied Mathematical Packages. Part 1. MathCAD. Saint Petersburg: RadioSoft; 2009.

15. Kudryavchenko I.V., Yakubov K.S., Palamarchuk D.V. Features of Artificial Intelligence Usage for Organisation of a SUAS Swarm. Modern Issues in Radioelectronics and Telecommunications. 2020;3:173.

16. Adnastarontsau AA, Adnastarontsava DA, Fiodortsev RV, et al. Algorithm for Control of Unmanned Aerial Vehicles in the Process of Visual Tracking of Objects With a Variable Movement’s Trajectory. Devices and Methods of Measurements. 2021;12(1):46-57.

Login or Create
* Forgot password?