Труды МАИ. 2021. № 121 Trudy MAI, 2021, no. 121
Научная статья УДК 519.711.2
DOI: 10.34759/trd-2021-121-25
ВЫБОР АРХИТЕКТУРЫ СИСТЕМЫ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ ПРЕОБРАЗУЕМОГО БЕСПИЛОТНОГО ЛЕТАТЕЛЬНОГО АППАРАТА - «КОНВЕРТОПЛАНА»
Аполлонов Дмитрий Вадимович1, Бибикова Кристина Игоревна2®,
Шибаев Владимир Михайлович3 Гаврилова Александра Владиславовна4
1,2,3,4Центральный аэрогидродинамический институт им. Н.Е. Жуковского, ЦАГИ,
Жуковский, Московская область, Россия
2ccfstd@tsagi.ru:
Аннотация. В настоящей статье представлена математическая модель преобразуемого беспилотного летательного аппарата («конвертоплана») и архитектура единой для всех режимов полета системы автоматического управления конвертоплана. Выполнено моделирование типовых режимов полета БЛА, а также проведен анализ возможности обеспечения управляемости в случае отказа органов управления.
Ключевые слова: математическая модель, конвертоплан, беспилотный летательный аппарат, система автоматического управления, отказ
Для цитирования: Аполлонов Д.В., Бибикова К.И., Шибаев В.М., Гаврилова А.В. Выбор архитектуры системы автоматического управления преобразуемого
беспилотного летательного аппарата - «конвертоплана» // Труды МАИ. 2021. № 121. DOI: 10.34759/trd-2021-121-25
CHOOSING THE ARCHITECTURE OF THE AUTOMATIC CONTROL SYSTEM OF A CONVERTIBLE UNMANNED AERIAL
VEHICLE - A TILTROTOR
Dmitry V. Apollonov1, Kristina I. Bibikova20, Vladimir M. Shibaev3, Alexandra V. Gavrilova4
1,2,3,4Central Aerohydrodynamic Institute named after N.E. Zhukovsky, TsAGI, Zhukovsky,
Moscow Region, Russia 2ccfstd@tsagi.com0
Abstract. Unmanned Aerial Vehicles (UAVs) has been the research subject of several
recent applications. This paper describes the flight dynamics simulation and automatic
control system architecture development of the tilt-rotor UAV. This tape of vehicle
combines the high-speed cruise flight capabilities of a conventional airplane with the
hovering capabilities of a helicopter by tilting their four rotors. The paper presents a
mathematical model of convertible unmanned aerial vehicle (tiltrotor) and the architecture
of a unified automatic control system for all flight modes of the tiltrotor. The authors
propose a control system structure that uses a limited set of «top-level» control parameters,
including the total engines thrust and three moments acting on the UAV. The principle of
forming the total thrust and moments by means of redundant set of «lower level» effectors,
including engines orientation angles and each engine thrust, as well as aerodynamic
2
control surfaces deflection, is determined. The possibility of using the proposed control principle is confirmed by the results of flight dynamics simulation of UAV executing a typical application scenario including vertical take-off and landing modes, flight along a given route and transient acceleration and deceleration modes. The paper considers the issues of assessing the safety of UAV flight in the event of possible individual control failures, taking into account the existing redundancy. A criterion for assessing controllability in the form of the maximum guaranteed angular acceleration is proposed which allows comparing the consequences of failures of each control effector. The paper shows that in order to guarantee UAV controllability we need to impose requirements on the choice of controller operational algorithms, both in the UAV angular motion control loop and in trajectory motion control loop.
Keywords: mathematical model, tiltrotor, unmanned aerial vehicle, automatic control system, failure
For citation: Apollonov D.V., Bibikova K.I., Shibaev V.M., Gavrilova A.V. Choosing the architecture of the automatic control system of a convertible unmanned aerial vehicle - a tiltrotor. Trudy MAI, 2021, no.121. DOI: 10.34759/trd-2021-121-25
1. Введение
Беспилотная авиация является одним из наиболее перспективных и быстро развивающихся направлений в современной авиации. С помощью беспилотных летательных аппаратов (БЛА) решают широкий круг задач, и расширение этого круга требует воплощения на практике принципиально новых характеристик БЛА.
Особый интерес представляют БЛА преобразуемого типа (иначе называемые «конвертопланами»), поскольку их можно использовать для реализации сценариев применения, при которых нет возможности обеспечить аэродромное базирование, но при этом необходимо, чтобы БЛА развивал относительно высокие скорости и преодолевал большие расстояния. [10,12,23] Главной особенностью конвертоплана является изменение аэродинамической конфигурации в процессе выполнения типовых режимов полета. Именно поэтому основными проблемами, с которыми сталкиваются разработчики системы управления конвертоплана, являются необходимость обеспечения устойчивости и управляемости БЛА на переходных режимах, а также разработка единой архитектуры системы управления для всех режимов полета. В настоящей статье предлагается подход, позволяющий реализовать архитектуру системы управления конвертоплана, учитывающую указанную выше особенность динамики полета такого БЛА.
В качестве примера реализации концепции преобразуемого беспилотного летательного аппарата рассмотрен БЛА с крылом тандемной схемы, имеющий четыре поворотных винтокольцевых движителя вентиляторного типа на концах консолей крыльев (рис. 1).
Рисунок 1 - Аэродинамическая схема БЛА
Принцип перехода из взлетного режима в режим крейсерского полета заключается в постепенном повороте всех двигателей (наклоне осей вращения вентиляторов в продольной плоскости) от вертикального положения до горизонтального с целью получения горизонтальной проекции суммарной тяги (рис.2). При этом, благодаря реконфигурации, улучшается аэродинамическое качество конвертоплана на крейсерском режиме полета и может быть существенно увеличена скорость перемещения в горизонтальной плоскости.
Переходный режки ^ ^ 91 Переходный режим
Рисунок 2 - Этапы полёта конвертоплана
5
Конфигурацию, при которой оси вращения двигателей расположены перпендикулярно продольной оси БЛА, будем называть «квадрокоптерным» режимом. В данном режиме БЛА совершает вертикально взлёт и посадку, а также полёт на малых поступательных скоростях. В свою очередь, конфигурацию, при которой оси вращений двигателей расположены вдоль продольной оси БЛА, будем называть «самолётным» режимом. БЛА в данном режиме совершает горизонтальный полёт на относительно высоких скоростях. Соответственно, режимы с промежуточной ориентацией двигателей будем называть «переходными».
1. Особенности объекта исследования 1.1 Этапы полёта БЛА
Основными этапами полёта БЛА, рассматриваемыми в данной работе, являются взлёт, вертикальный набор высоты, разгон, полет по заданному маршруту с изменением высоты и направления полета, торможение, вертикальное снижение и посадка (исходим из того, что БЛА всегда производит взлет и посадку вертикально, т.е., «по-вертолетному»).
В квадрокоптерном режиме оси вращения двигателей направлены вверх. Подъемная сила создаётся за счёт работы винтомоторной группы. Таким образом, в этом режиме БЛА имеет возможность совершать взлёт и посадку вертикально, а также перемещаться с небольшими поступательными скоростями в любом направлении в горизонтальной плоскости
В самолетном режиме оси вращения двигателей направлены вдоль продольной оси летательного аппарата. Подъемная сила в таком случае создается в основном за
счет воздействия аэродинамического набегающего воздушного потока на поверхность БЛА. В этом режиме летательный аппарат может осуществлять горизонтальный полет, набор высоты и снижение, изменение направления полета.
Переходный режим реализуется путём отклонения оси вращения двигателей на угол причём £ может изменяться от 0° до 90°. Этот режим реализуется в процессе полёта при переходе от вертолётного режима полёта к самолётному и обратно.
1.2 Органы управления БЛА в квадрокоптерном режиме полёта
В квадрокоптерном режиме полёта управление подъемной силой производится одновременным соразмерным изменением тяги всех двигателей с учетом сохранения равновесия относительно поперечной оси в зависимости от положения центра масс БЛА. Поперечное управление (по углу крена у) производится аналогично БЛА квадрокоптерной схемы - путём изменения соотношения суммарной тяги двигателей, установленных слева и справа относительно продольной оси БЛА (рис.3). Для БЛА небольшой взлетной массы запаздывание в канале управления тягой каждого двигателя мало по сравнению с постоянной времени в канале крена и не оказывает заметного влияния на характеристики устойчивости и управляемости БЛА. Для аппаратов большего взлетного веса это запаздывание должно быть учтено при формировании алгоритмов САУ. Аналогично, для продольного управления -по углу тангажа 0 - путём изменения соотношения суммарной тяги двигателей, установленных перед центром масс БЛА и позади него.
Рисунок 3 - Управление БЛА по углу крена в "квадрокоптерном" режиме
Нумерация двигателей, а в дальнейшем - и аэродинамических рулей на панелях крыльев, производится по часовой стрелке, начиная с левого переднего.
В отличие от БЛА квадрокоптерной схемы, разворот по углу рыскания у в квадракоптерном режиме выбранного БЛА осуществляется путём дифференциального изменения угла установки двигателей Л^, расположенных на концах консолей переднего крыла (рис.4). Положительным считается отклонение вектора тяги левого двигателя назад, а правого - вперед на угол Л^.
Рисунок 4 - Управление БЛА по углу рыскания в "квадрокоптерном" режиме
1.3 Органы управления БЛА в самолётном режиме полёта
В самолётном режиме синхронное изменение тяги двигателей служит для управления скоростью полета, при этом сохраняется возможность использовать управление угловым движением путем независимого изменения тяги двигателей и дифференциального отклонения передних двигателей, однако по мере роста приборной скорости увеличиваются эффективности аэродинамических органов управления.
Так управление по углу тангажа 0 и крена у осуществляется путём отклонения рулевых поверхностей, расположенных вдоль задних кромок переднего и заднего крыльев (рис.5). При формировании управляющих сигналов рассматривались два варианта.
Например, при управлении по крену в первом случае рулевые поверхности на левых консолях переднего крыла - флапперон - и заднего крыла - элевон -отклоняются синхронно на один и тот же угол в одну сторону, а на правых консолях - на тот же самый угол, но в противоположную сторону, тогда параметром, определяющим управляющее воздействие будет считаться одна переменная - «угол отклонения элерона» 6а. Для моделирования системы управления при такой конфигурации системы управления при продувках в АДТ получены зависимости поперечного момента Мх от 6а - синхронного отклонения всех элевонов и флапперонов в указанных направлениях.
Исследуется и другой подход, при котором отклонения всех четырех аэродинамических рулей рассматриваются как независимые органы управления
^/1, ^/2, ^ез, Для моделирования системы управления при такой конфигурации системы управления при продувках в АДТ получены зависимости поперечного момента Мх от отклонения отдельных элементов аэродинамических рулей $/1, ^/2, ^е3, ^е4. Считаем, что
Аналогично, при управлении по тангажу рулевые поверхности на переднем крыле - флаппероны - и на заднем крыле - элевоны - могут отклоняться синхронно на один и тот же угол, но в противоположные стороны, тогда параметром, определяющим управляющее воздействие будет считаться одна переменная -«отклонение руля высоты» 8е. Для моделирования системы управления при такой конфигурации системы управления при продувках в АДТ получены зависимости продольного момента Мг от 8е - синхронного отклонения всех элевонов и флапперонов.
Аналогично каналу крена, рассматривается и другой подход, при котором отклонения всех четырех аэродинамических рулей рассматриваются как
М8/1 + М8/2 + М8- - М8Х
Рисунок 5 - Управление БЛА по углу крена в "самолётном" режиме
независимые органы управления 8^±,8^2,8е3,8е4. Для моделирования системы управления при такой конфигурации системы управления при продувках в АДТ получены зависимости продольного момента Мг от отклонения отдельных элементов аэродинамических рулей 8^2,8е3,8е4. Считаем, что
М¡е = м8/1 + М8/2 - М¡е3 - М8/4
Ниже, в разделе 8, проведен анализ возможности обеспечения продольной и поперечной управляемости при отказе органов управления на крейсерском режиме в предположении, что имеется возможность управлять всеми четырьмя органами управления независимо.
Разворот по углу рыскания у осуществляется посредством отклонения руля направления 8Г, расположенного на задней кромке вертикального оперения в нижней части фюзеляжа БЛА (рис.6).
Рисунок 6 - Управление углом рыскания в "самолётном" режиме
1.4 Управление БЛА на переходных режимах
Переходный режим реализуется при разгоне в процессе перехода от квадрокоптерной конфигурации к самолетной, или при обратной реконфигурации в процессе торможения. Реализация переходного режима заключается в синхронном изменении углового положения осей вращения всех четырех винтокольцевых движителей вентиляторного типа.
На малых скоростях полета углы установки осей вращения двигателей относительно нормальной оси БЛА небольшие <30 град. При этом высотой полета и угловым положением БЛА можно управлять как в квадрокоптерном режиме.
При сравнительно больших скоростях полета углы установки осей вращения двигателей > 60 град. Управление угловым положением в таком случае осуществляется аналогично управлению, реализуемому в самолётном режиме.
В большинстве работ, посвященных системам управления конвертопланов, предлагается формировать независимо законы управления для вертолетной и самолетной конфигурации, при этом в процессе разгона суммарное отклонение органов управления рассчитывается как линейная комбинация «вертолетного» и «самолетного» управления с весовыми коэффициентами, линейно меняющимися от нуля до единицы в выбранном диапазоне скоростей полета или углов установки двигателей (рис. 7) [2,17,26,28,29].
Рисунок 7 - Структура системы управления конвертопланом
В настоящей работе предлагается вклад каждого из органов управления определять автоматически в зависимости от скорости полета и ориентации двигателей путем задания матрицы расчетных и экспериментально определенных эффективностей.
2. Математическая модель динамики полета конвертоплана В общем виде структурная схема моделирования динамики управляемого движения конвертоплана выглядит следующим образом:
Рисунок 8 - Структурная схема моделирования
Для обеспечения возможности моделирования динамики полета БЛА на разных режимах формирование управляющих сигналов в системе автоматического управления разделено на два этапа. На первом этапе формируются потребные значения управляющих сил и моментов, определяемые отклонениями текущих значений параметров движения и ориентации БЛА от заданных с помощью опорной («заданной») траектории. При таком подходе алгоритмы в блоке «Расчет потребных сил и моментов» незначительно зависят от конфигурации БЛА, в основном - за счет изменения располагаемых предельных значений создаваемых сил и моментов. На втором этапе, в модуле «Расчет отклонений органов управления», выполняется вычисление управляющих сигналов, непосредственно изменяющих значения тяг всех четырех двигателей (Т1, Т2, Т3 и Т4), дифференциального отклонения осей вращения передней пары двигателей для управления по рысканью на режимах малых скоростей (А£,), отклонений рулевых поверхностей в каналах тангажа
(условно - «флаппероны»), крена («элевоны») и рысканья («руль направления») - (8/,
§e, 8Г).
Рисунок 9 - Структурная схема модуля ЛА
2.1 Исходные данные, используемые при моделировании
m = 30 [кг], масса;
! = 1[м], размах крыла;
Lxl2 = 1[м], расстояние от точки приложения сил ВМГ первого крыла до центра масс вдоль оси OX;
Lxз4 = 1[м], расстояние от точки приложения сил ВМГ второго крыла до центра масс вдоль оси OX;
Ly12 = 0.036 [м], расстояние от точки приложения сил ВМГ первого крыла до центра масс вдоль оси OY;
Ly34 = 0,14 [м], расстояние от точки приложения сил ВМГ второго крыла до центра масс вдоль оси OY;
Lz12 = 1,75 ,расстояние от точки приложения сил ВМГ первого крыла до центра масс вдоль оси 07, [м];
Lz34 = 1,75 ,расстояние от точки приложения сил ВМГ второго крыла до центра масс вдоль оси 07, [м];
1х - Момент инерции вокруг оси ОХ, [кг*м2]; 1у - Момент инерции вокруг оси 0Y, [кг*м2]; \z - Момент инерции вокруг оси 0Z, [кг*м2]; g - Ускорение свободного падения, [м/с2]; Ьа - Средняя аэродинамическая хорда, [м].
Также в качестве исходных данных использовались результаты продувок в АДТ. На рисунке 10 приведена зависимость коэффициента подъемной силы от угла атаки для одной из конфигураций БЛА:
а,град.
-10
Суа
0,8 0,6 0,4 0,2 0 -0,2 -0,4 -0,6 -0,8
4 6 8 10
Рисунок 10 - Коэффициент подъемной силы в зависимости от угла атаки
2.2 Уравнения для описания пространственного движения БЛА
В качестве математической модели динамики полета БЛА на всех режимах полета использовались следующие уравнения:
Уравнения движения центра масс ЛА в проекциях на оси связанной системы координат OXYZ [1,6]:
т (^ + — = Их1 — Ха cos а cos р + Уа бта— 2а cos а sin р — mgsin Э
/йУу N
т + лм-гУх — ) = Яу1 + Ха бш а соз р + Уа соб а + 2а бш а бш р — mg соб Эсоб у,
0
m + wxVy — wyVx) = Rz1 — Xa sin P + Za cos p + mg cos -9sin y,
Yx, Vy, Vz- проекции вектора скорости центра масс ЛА на оси связанной системы координат; wx, wy, wz - проекции вектора угловой скорости вращения ЛА относительно его связанных осей OXYZ; X, Y, Z - проекции вектора аэродинамических сил на оси связанной системы координат, Rxl, Ry1, Rz1 -проекции вектора тяги на оси связанной системы координат.
Изменение дальности полета Хд, высоты полета H и боковой дальности Zg определяются проекциями скорости на оси нормальной земной системы координат: ^- = Vx cos тд cos -ф — Vy (cos у sin тд cos -ф — sin у sin -ф) + VZ (cos y sin ^ + sin y sin t9 cos -ф);
dH
— = Vxsind + Vy cos y cos д — Vzsiny cos д;
^ = —Vx cos-d sin-ф + (cos7sin^ sin^ + sinycos-ф) + Vz(cosy cos-ф — siny sin$ sin^),
где у, 0, у - углы крена, тангажа и рысканья соответственно.
Вычисление угловых скоростей в связанной СК выполняется в соответствии с уравнениями:
~át—lxy~dt+ (Iz — + = МХ + МХА,
— 1ху~ГГ + Ох — IZ)MXMZ — 1хушушг = Му + М.
dt У dt 1хушушх L'ly 1 L'lyg>
lz~It+ (Iy — 1х)ШУШх + 1ху(шу2—шх2) = Mz + MZA,
18
где Iх, I у, I z - осевые моменты инерции; Мх, Му, Mz - проекции вектора аэродинамического момента на оси связанной системы координат; Мхд, Муд, MZA -проекции моментов сил, развиваемые ВМГ; 1ху - центробежный момент инерции.
Кинематические соотношения для углов Эйлера:
^ = (Шу cos Y-Mz sin Y)tgd,
d— = mz cos y + шу sin y,
d-ф Ir . \
— =-(шу cos y — MZsm.Y).
dt cos-y у ' z lJ
2.3 Силы и моменты, действующие на БЛА 2.3.1 Аэродинамические силы и моменты
Ха = СхацБ - сила лобового сопротивления [Н],
Ya = СуацБ - подъемная сила [Н],
Za = CzaqS - боковая сила [Н],
где Сха, Суа, Cza - безразмерные коэффициенты соответствующих сил; S - площадь крыла [м2]; q - скоростной напор, зависящий от скорости БЛА:
pV2
q = —
4 2
Вектор аэродинамического момента:
м =
м. м.
у
м
где проекции вектора аэродинамического момента на оси связанной системы координат вычисляются следующим образом:
Мх = тхц51,
Му = mуqSl, Мг = т^Ьа,
тх, ту, тг - безразмерные аэродинамические коэффициенты соответствующих моментов; I - размах крыла [м]; Ьа- средняя аэродинамическая хорда [м].
Коэффициенты сил и моментов в зависимости от фазовых координат и отклонений органов управления рассчитывались на основании результатов испытаний модели БЛА в АДТ.
2.3.2 Силы и моменты, создаваемые винтомоторной группой
У
X
Рисунок 11 - Направление сил и моментов, создаваемых ВМГ
Значения сил тяги ВМГ рассчитывается по формуле: Т = даарп204,
где Т1 - сила тяги 1-го двигателя; I - номер двигателя ВМГ(г=1,2,3,4); В - диаметр винта; ай - коэффициент тяги винта; щ - скорость вращения ¿-го двигателя; д -ускорение свободного падения; р - плотность воздушной среды.[15,20]
Проекции суммарного вектора силы тяги ВМГ на оси связанной системы координат БЛА при этом будут выглядеть следующим образом:
КХ1 = (Т- + т2)бш%1 + (Т3 + Т^йп^ ПУ1 = (Т + т2)соб^1 + (Тз + Т^СОБ^, Яг1 = 0,
где ^ 2 - углы отклонения оси 0дУд, направления сил тяги пары двигателей от направления оси О У1.
Момент силы тяги ВМГ, развиваемый вокруг оси ОХ1 БЛА.:
МХд = 1212(Т1 - АО - Т2 + АО) - Ь234(Т3СО5%2 + Т^СОБ^У,
Момент силы тяги ВМГ, развиваемый вокруг оси ОУ1 БЛА:
Муд = 1212(-Т1 - АО + Т2 + АОу + ^(Т^п^ - Т^т&У
где Ьг12 - расстояние от приложения сил ВМГ первого крыла до центра масс вдоль оси 07; Ьг34 - расстояние от приложения сил ВМГ второго крыла до центра масс вдоль оси 07
Момент силы тяги ВМГ, развиваемый вокруг оси 02х БЛА определяется соотношениями:
М = М1Х + М1у
Мгд Мгд + мгд ,
= 1х12(т- соя(^1 - АО + Т2 соя(Ь + АО) - Ьх34(Т3СО5^ + Т4СО5^),
м% = Ьу12(Т.1 - АО + Т2 + АО) - Ьуз^тЪ + Т^т^У
где Ьх12, Ьу12 - расстояние от приложения сил ВМГ первого крыла до центра масс вдоль осей ОХ и ОУ соответственно; Ьх34, Ьу34 - расстояние от приложения сил ВМГ второго крыла до центра масс вдоль оси ОХ и ОУ.
3. Расчёт потребных управляющих отклонений органов управления
Для реализации универсального алгоритма расчёта управляющих сигналов необходимо воспользоваться уравнениями движения летательного аппарата. Рассмотрим алгоритм на примере «квадрокоптерного» режима - режима полёта на малых скоростях или стабилизации заданных координат вблизи заданной точки. В данном режиме преобразуемый беспилотный летательный аппарат может быть описан уравнениями движения, аналогичными используемым при исследовании динамики полета квадрокоптера [8]:
г.. ТЕ
Х = (бЫ ф бЫ у — соб ф бЫ д соб у) —
ш
те
У = — д + (соб д соб у) —
ш
те
2 = (соб ф бЫ у + бЫ ф бЫ д соб у) —
ш
Ьг — !уу , ¡тр , Мх шх=----шуш2 + —ш2П+ ■—
1ХХ 1ХХ 1ХХ
IXX — 122 , Му Шу Шу=----ШхШ2+---МуУШу
1УУ 1УУ
¡УУ — ¡ХХ ¡тр _ , Мж 0)2= — ---Шх^у--—ШхЫ+ -—
1 гг 1 гг 1 гг
В процессе моделирования использовались система нелинейных уравнений
движения БЛА. Однако для выбора параметров алгоритмов системы управления
уравнения движения были линеаризованы:
' Х = -дд
•• Те У = --д т
1 = У9
. мх
*хх
му
Ыу= —
1УУ
. Мг
В качестве компонентов вектора управления в таком представлении выступают суммарная тяга и три управляющих момента - моменты крена, тангажа и рыскания. [9,19]
Так как суммарная тяга и управляющие моменты определяются положениями органов управления и режимами работы двигателей, можно получить закон этой
зависимости в матричном виде:
Мх Му
Мт
1 1 1 1 0 0 0 0
М11 мтх2 мтх3 мр м? 0 М8Г 11х Мр
мр Мту2 мр мр М? 0 0 Мр
мр М™ мр мр М? мр 0 0
Тг Т2 Тз Т
Ц
8е 8Г
8-т
,тогда
Т£
Мх Му М*
= В
Тг Т2 Тз
Т
Ц 8е
8
7
где
В =
1 1 1 1 0 ^р8е 0 0
мТ1 мТ2 мХ3 мТ4 м? 0 м8г 11х м8хг
Щ1 мТ Щ3 мТ4 0 0 м8/
мТ1 мТ2 мТ3 мТ4 м? м88е 0 0
Элементы матрицы В представляют собой эффективности управления, т.е. производные суммарной силы и моментов по каждой компоненте из вектора управляющих параметров. В зависимости от режима полета и конфигурации БЛА значения каждого элемента матрицы эффективности определяются с помощью расчетов или путем использования результатов продувок в АДТ.
Можно заметить, что количество неизвестных получается больше количества уравнений. Довольно известным методом решения такой системы уравнений для расчета потребных значений параметров управления в таком случае является использование псевдообратной матрицы [13,25]. Тогда с помощью псевдообратной матрицы можно найти одно из решений уравнений.
п
Т2
Тз те
Т4 ^ = В+ мх му
мг
*т
8Г
, где В+ = ВТ(ВВТ) 1,В+ - матрица, псевдообратная матрице В.
Следовательно, из уравнений в таком виде можно получить тягу двигателей,
изменение угла дифференциального наклона осей вращения передней пары
двигателей и отклонения аэродинамических органов управления. Полученные
25
параметры подаются в качестве параметров управления в блок расчета сил и моментов модели динамики полета БЛА. Таким образом, управление БЛА происходит путём изменения моментов относительно связанных осей летательного аппарата и тяги двигателей.
Для того, чтобы управляющие параметры использовались в равной степени в пределах диапазона своего изменения, вводится диагональная нормирующая матрица зависящая от упомянутых диапазонов:
w =
1
Т1т 0
0
0
0
0
0
0
2
0
1
Т2 т 0
0
0
0
0
0
2
0
0 1
Тз т 0
0
0
0
0
2
0
0
0 1
Т
2
4т 0
0
0
0
0 1
^т2 00
0
0
0
0
0 1
5,
0 0
0 0
ет
0
0
0
0
0
0
0
1
6,
2
-7т 00
0
0
0
0
0
0
0 1
8
гт
Тогда выражение для псевдообратной матрицы примет вид:
В+ = wBT(BwB1)
Тл-1
Отметим, что если все параметры приведены к единому диапазону значений, из теории псевдообратных матриц следует, что норма решения, т.е. корень из суммы
2
2
квадратов нормированных значений всех полученных управляющих параметров, будет минимальной [8].
Стоит отметить, что матрица В состоит из эффективностей, зависящих от режима полета и конфигурации БЛА. Это означает, что используя данный алгоритм, можно рассчитать управляющие воздействия для любого из режимов полета, вне зависимости от того, в какой конфигурации были рассчитаны потребные значения сил и моментов.
4. Формирование управляющих сил и моментов
4.1 Режим стабилизации координат вблизи заданной точки
При полёте на малых поступательных скоростях преобразуемый БЛА может быть описан уравнениями движения квадрокоптера. Управление в вертикальной плоскости реализуется изменением тяги. В горизонтальной плоскости БЛА управляется изменением угла тангажа. При этом структура системы управления будет состоять из контуров управления траекторным и угловым движением. [3,16,18,21] В контуре управления траекторным движением можно управлять линейным положением БЛА и формировать суммарную тягу и углы тангажа и крена. Также можно управлять БЛА по углу рыскания. В контуре управления угловым движением формируются управляющие моменты крена, тангажа и рыскания. Полученные моменты и суммарную тягу можно использовать для расчета управляющих воздействий с помощью псевдообратной матрицы. На выходе блока расчёта формируются управляющие сигналы, распределяемые по органам
управления. В данном режиме учитываются четыре тяги двигателей и угол дифференциального отклонения передних двигателей.
Рисунок 12 - Структурная схема системы управления в "вертолётном" режиме
4.2 Режим полета по заданной траектории
В режиме полета по заданной траектории в соответствии с режимом, определяемым путем задания значений скорости и высоты полета, формируются заданные значения составляющих сил Рхд и Руд в продольной плоскости. Управление в «самолётном» режиме реализуется за счёт соразразмерного изменения угла тангажа и суммарной тяги. Поэтому с помощью уравнений равновесия для горизонтального полёта определим зависимость этих переменных от потребных значений горизонтальной и вертикальной составляющих сил:
Рисунок 13 - Схема сил, приложенных к центру масс БЛА в горизонтальном полете
Уравнения равновесия сил, в зависимости от изменения угла отклонения осей
вращения двигателем и от угла тангажа, имеют вид:
(Fxg=Tzsin(i;-id)-Xa \Fyg = Tz costf - д) + Ya
Линеаризуя данные уравнения, представим их в матричном виде:
(Fxq = Tzsinl;-Tzcos^d + Xad
Fyg = TE cos f + Tzsin% *д + Yad
или
F Fx g
F Fy g
sint; -mgcos% cos% mgsinf + Ya
TE д
= A
TE д
где A =
sinf -mgcos%
cos% mgsin% + Ya
Решая эти уравнения, можно сформировать заданные значения суммарной тяги двигателей и угла тангажа в зависимости от потребных значений горизонтальной и вертикальной составляющих сил, рассчитываемых с помощью алгоритмов блока управления траекторным движением БЛА:
ТЕзад изад
= А-1
р
гхд р
гуд
Управление в боковом канале выполняется путем формирования заданных значений углов крена и рыскания, которые, в свою очередь определяются на основании заданных координат промежуточных пунктов маршрута. [5]
В зависимости от боковых отклонений от заданной линии пути, соединяющей промежуточные пункты маршрута, формируется заданный угол крена:
У зад = К2Аг — KvzVz,
где Аг - боковое смещение относительно заданной линии пути, Vz.- скорость в направлении, перпендикулярном линии пути.
¥ / ППМ2 ч
/ /___ X V ч ч
ППМЗ
V
щ
МЩМ1
О 50 100 150 200 250 300 350 400 450 500
г
Рисунок 14 - Траектория полёта по заданной линии пути
Коэффициенты КРг и Кг выбираются с учетом быстродействия управления в канале крена и диапазона допустимых значений угла крена
Заданная трасса полета прокладывается в виде последовательности отрезков прямых линий, составляющих ломаную линию. Отсчет отклонений Аг от заданной линии пути производится относительно этой линии по перпендикуляру к её отдельным составляющим. Переход отсчета с одного отрезка на другой начинается с боковым упреждением разворота (БУР).
Заданные значения углов тангажа, крена и рыскания передаются в контур управления угловым положением, где происходит формирование потребных управляющих моментов относительно всех трех осей связанной системы координат. Моменты в совокупности с суммарной тягой поступают на вход блока формирования отклонений органов управления, при этом используется описанный выше алгоритм с применением псевдообратной матрицы и оценок эффективности
располагаемых органов управления в каждом канале. В данном режиме полета, помимо тяги всех четырёх двигателей и угла дифференциального отклонения передних двигателей, задействованы и аэродинамические рулевые поверхности -элевоны, флапероны и руль направления (рис.15).
Рисунок 15 - Структурная схема системы управления в "самолетном" режиме
5. Область допустимых углов наклона осей вращения двигателей в зависимости от скорости и тяги
Чтобы выяснить, в каком диапазоне могут изменяться углы наклона осей вращения двигателей £ в зависимости от скорости полета, необходимо рассмотреть уравнения равновесия для горизонтального полёта БЛА:
рУ2
Т1Бт(£ — -в) — Сх^-5 = 0
тхсоб(^ -тд) + Су^—Б — тд = 0
рУ2 ~2
pV2
Tz [sin f cosd — cos^sind] — Cx——S = 0
2
Tz [cosf costi — sin^simd] + Cy^—S — mg = 0
pV2 ~2
Считая, что угол тангажа мал, полагаем cosd ~ 1 и sind ~ д, тогда запишем уравнения в следующем виде:
pV2 >
pV2
Tzsinf — Tzcos% ti — Cx——S = 0
2
Tzcos% — Tzsin% ti — Cy——S — mg = 0
2
pV2 pV2
Tz — sin% Cx——S + cos% Cy——S — mgcos% = 0 22
Ру2 г л
Tz = ——S(Cxsin% — Cycos%) + mgcos% 2
Тогда при <f = 0:
pV2 pV2
Cy-S ti — mg ti — Cx-S = 0
22
При f = 90°:
pV2 pV2
Cx—S ti + Cy — S — mg = 0
Таким образом, из условия равновесия получим балансировочные значения угла тангажа БЛА на различных углах наклона двигателей в зависимости от скорости:
Рисунок 16 - Балансировочное значение угла тангажа БЛА в зависимости от
скорости полета БЛА
Балансировочные значения тяги БЛА на различных углах наклона двигателей в зависимости от скорости имеют вид:
Рисунок 17 - Балансировочное значение тяги в зависимости от скорости полета БЛА
Таким образом, выбирая на каждой скорости угол наклона двигателей, для которого балансировочный угол тангажа (равный для режима ГП углу атаки БЛА), а тяга минимальна или, по крайней мере, находятся в диапазоне располагаемых значений, можно определить рекомендуемый закон изменения углов наклона осей вращения двигателей в зависимости от скорости полёта:
Зависимость угла наклона двигателей от скорости БЛА
О
(и
(и
т
(С
|_
ш я
Ч а
(С и
X
о
к
(С
X
о
>
100 90 80 70 60 50 40 30 20 10 0
ООО
ОЖОООО
офоооо
офооооо
О о (1)18)18)18)18)®
0 50 100 150
скорость БЛА, м/с
О Допустимые углы установки двигателей в зависимости от скорости
Рекомендуемый закон изменения угла установки двигателей
Рисунок 18 - Область изменения угла установки двигателей
6. Примеры моделирование типовых режимов полёта
На основе разработанных алгоритмов системы управления конвертопланом было проведено моделирование основных этапов полёта. В частности - полёт на малых скоростях вблизи заданной точки и полёт по заданной траектории на высоких скоростях.
На рисунке проиллюстрирована траектория перемещения конвертоплана вблизи контрольной точки по заданной траектории, с последовательным изменением и последующей стабилизацией высоты и смещений вдоль оси ОХ и 07.
Рисунок 19 - Траектория полёта вблизи заданной точки
При этом текущие параметры траектории БЛА отслеживают заданные значения с требуемыми характеристиками качества, которые в данном случае определялись требуемым временем переходного процесса и отсутствием статической ошибки и перерегулирования:
Высота Н, м
1 08 0.6 :.-0.2 0
/ \
/ \
/ \
/ V
/
,,10-з Тангаж тЭ, ра/
0 -5 -10 -15
\
И
V
Перемещение по оси ОХ, м
1 03 0.6 0.2 0
0 10 20 30 40 50 60 70 60
Рисунок 20 - Изменение параметров БЛА при полёте вблизи заданной точки
На основе разработанных алгоритмов системы управления также был промоделирован полёт по заданной линии пути, определяемой промежуточными пунктами маршрута (ППМ). На рисунке 21 изображена траектория полёта. Летательный аппарат выполняет взлет вертикально «по-вертолётному» с вертикально ориентированными двигателями и по достижении высоты 10 м начинает увеличивать скорость в горизонтальной плоскости. По мере увеличения скорости двигатели отклоняются в продольной плоскости в соответствии с рекомендуемым законом. При достижении максимального значения скорости двигатели располагаются полностью горизонтально по отношению к продольной оси БЛА. Таким образом, полет по маршруту выполняется в квадрокоптерной конфигурации (начальный этап полета, включая разгон до 10 м/с и допорот по
направлению к ППМ1), затем в переходной конфигурации, включая этап выполнения бокового упреждения разворота при прохождении точки ППМ2 и доворота по курсу до направления к точке ППМ3 и, наконец, в самолетной:
Рисунок 21 - Траектория полёта БЛА по заданной линии пути
На всех этапах полета текущие параметры БЛА соответствуют заданным значениям. На графике видно постепенное изменение конфигурации БЛА за счет
изменения угла наклона двигателей в соответствии с изменением скорости полёта.
0.1 с -а*
■м
9< Ы
Высога Н. м
1
1
м
ЗА
ю о
г*
аз о ■02
Угол крема у. рал
- 'Л 1 \ -
- - А \ - Г - -
Скорость VI, м/с
_ _ _ а к,-_ _____
-
•
Сумма риая тяга Тзиш. Н
Рисунок 22 - Изменение параметров БЛА при полёте по заданной траектории
7. Анализ возможности обеспечения управляемости при отказе органов управления
Избыточность органов управления по сравнению с количеством контролируемых параметров позволяет поставить задачу анализа отказобезопасности системы управления, т.е., рассмотреть возможность обеспечить продолжение управляемого полета после возникновения отказа одного или
нескольких исполнительных элементов. Системы управления, обеспечивающие такую возможность, называют «устойчивыми к отказам» (fault-tolerant) [22,24].
Рассматривались два варианта парирования последствий отказа аэродинамических рулей - элеронов и флапперонов - на режиме крейсерского полета при управлении углами крена и тангажа. В этом случае в штатном случае мы имеем четыре рулевых поверхности S^1,6^2,6e3,6e4, с помощью отклонений которых формируются заданные значения моментов относительно продольной и поперечной осей, как было изложено в разделе
В первом случае («пассивное парирование отказа»), при постоянном законе формирования отклонений рулей с использованием информации об эффективности каждого из них:
D = В+М
зад,
где
В+ = wBT(BwBT)-1,
D - столбец отклонений рулевых поверхностей
Sn Sf2
8ез Se4
M - столбец проекций суммарного момента на продольную и поперечную оси B - матрица эффективностей
W - весовая матрица, учитывающая степень влияния каждой рулевой поверхности на формирование моментов.
При этом фактическое значение, вычисляемое как М = ВО
при отказе, например, рулевой поверхности №4 (6е4), будет отличаться от Мзад. Это может рассматриваться как внешнее возмущение, которое должно парироваться САУ (рис. 23):
Рисунок 23 - Пассивное парирование отказов
В этом случае вместо ПИД-регуляторов должны использоваться алгоритмы, устойчивые к изменению характеристик системы управления («робастные»)[4,7]. В качестве такого алгоритма был выбран регулятор, реализованный на основе концепции обратной задачи динамики (ОЗД-регулятор) [11,27]. Примеры моделирования работы системы автоматического управления в случае, когда все органы управления были задействованы, и в различных комбинациях отказов приведены на рис. 24 - рис. 32. В приведенных примерах для сравнения в канале крена контура угловой стабилизации использовался ПД-регулятор, в канале тангажа - ОЗД-регулятор.
На графиках показаны как переходные процессы по углам крена и тангажа в сравнении с их заданными значениями, так и командные сигналы, сформированные на уровне заданных моментов Мх и Мх и на уровне отклонений аэродинамических органов управления 6е3,6е4, рассчитанных с учетом матрицы
коэффициентов эффективности управления.
Для сравнения приведены заданные значения управляющих параметров (моментов и органов управления) и фактически реализуемых, с учетом располагаемых диапазонов отклонений и отказов.
0.025 0.02 "Л 0.015
е.
с" 0.01
X ф
0.005 О
-0.005 0.12
4
[[сек]
0.1
ов
И 0.
а
Ч>" 0.06
I 0.04
№
I-
0.02 О
— — —Заданное значение Текущее значение
---Зал анное значение /щее значение -
\
\
\
^ V
4
1[сек]
Рисунок 24 - Изменение углов крена и тангажа при правильно работающей системе
(без отказов)
В канале крена - ПД-регулятор, в канале тангажа - ОЗД-регулятор
Рисунок 25 - Отклонение рулевых поверхностей при правильно работающей системе
(без отказов)
Рисунок 26 - Изменение моментов крена и тангажа при правильно работающей
системе (без отказов)
25
20
¥ 15
а
(** 10
I
Я 5
0
-5
I
— — —Заданное значений Текущее значение
4
t [сек]
Ц 0 09
а.
" 0.06
I 0.04
I
Р 0.02
О
I — — —Заданное значение Текущее значение
4
t [сек]
Рисунок 27 - Изменение углов крена и тангажа при отказе левого флапперона и
левого элевона; в канале крена ПД-регулятор, в канале тангажа - ОЗД-регулятор
44
Рисунок 28 - Отклонение рулей при отказе левого флапперона и левого элевона; в канале крена ПД-регулятор, в канале тангажа - ОЗД-регулятор
Рисунок 29 - Изменение моментов крена и тангажа при отказе левого флапперона и левого элевона; в канале крена ПД-регулятор, в канале тангажа - ОЗД-регулятор
Во втором случае («активное парирование отказа»), при возникновении отказа происходит реконфигурация системы управления рулями.
Рисунок 30 - Активное парирование отказов
Например, при отказе руля № 4 расчет заданных значений отклонений будет выполняться с весовой матрицей W, у которой диагональный элемент, соответствующий отказавшему рулю, будет положен равным нулю.
Ю =
10 0 0 0 10 0 0 0 10 0 0 0 1
=> w =
1000 0100 0010 0 0 0 0
В этом случае для системы управления отказ фактически не существует (в тех
случаях, когда реконфигурация позволяет сохранить управляемость, в частности -
отказавший руль находится в нейтральном положении и эффективности оставшихся
рулей достаточно для создания потребных угловых ускорений), так как
формируемое управляющее воздействие Мзад без искажения передается на объект управления.
Рисунок 31 - Изменение углов крена и тангажа при отказе левого флапперона и левого элевона, реконфигурация; в канале крена ПД-регулятор, в канале тангажа
ОЗД-регулятор
* 1
щ 1
Ё
---<Ь_фл1 заданное -фактическое значение
4
I [сек]
^_____ ---<и -фа >л2 заданное тическое значение
/
. ® -г
Е 1
С)
с;
* 1
из 1
Е с;
О
-!Г
С
Сз "= -1
4
I [сек]
--- -фа г < заданное тическое значение
к"
4
I, [сек]
¡1 я4 заданное тическое значение
4
I [сек]
Рисунок 32 -Отклонение рулей при отказе левого флапперона и левого элевона, реконфигурация; в канале крена ПД-регулятор, в канале тангажа - ОЗД-регулятор,
незначительный выход за ограничения
Рисунок 33 - Изменение моментов крена и тангажа при отказе левого флапперона и левого элевона, реконфигурация; в канале крена ПД-регулятор, в канале тангажа -ОЗД-регулятор, незначительный выход за ограничения
Наличие ограничений на перемещения органов управления (в данном случае -аэродинамических рулей) приводит к тому, что диапазон воспроизводимых угловых ускорений существенно уменьшается, поэтому даже с использованием активного парирования отказов рекомендуется учесть в настройках контроллеров управления располагаемый диапазон управляющих моментов (рис.).
20
¥ 15
с.
10
I
£ 5
1 — — —Заданное значение -Текущее значение
4
I [сек]
- - — Зaj данное значение сущее значение
0.1
¥ о.
0.05
4
t [сек]
Рисунок 34 - Изменение углов крена и тангажа при отказе левого флапперона и правого элевона, реконфигурация; в канале крена ПД-регулятор, в канале тангажа ОЗД-регулятор, существенный выход управления за ограничения
---£|_фл1 заданной -фактической значений
| 10 с
0
>|
с
е -ю
4
% [сек]
4
t [сек]
4
^ [сек]
4
I, [сек]
^ — ^ _ ---<1 _<1 -фа л2заданное тическое значение -
Ч >
---с!_элЗ заданное фактическое значение
---<1_эл4 заданное -фактическое значение
/ *
\ .'
Рисунок 35 -Отклонение рулей при отказе левого флапперона и правого элевона, реконфигурация; в канале крена ПД-регулятор, в канале тангажа - ОЗД-регулятор,
существенный выход за ограничения
Рисунок 36 - Изменение моментов крена и тангажа при отказе левого флапперона и правого элевона, реконфигурация; в канале крена ПД-регулятор, в канале тангажа -ОЗД-регулятор, существенный выход за ограничения
На рисунках ниже приведены для сравнения границы располагаемых угловых ускорений в каналах крена и тангажа при всех работающих органах управления и
при отказе одного или двух элементов.
Рисунок 37 - Области располагаемых угловых ускорений в каналах крена и тангажа в случае всех работающих органов управления и при отказе одного элемента
-5 г
-2 ■
м
Муу
N О
-3
■ 1 - все рули работают
2 - отказ флапперонов
3 - отказ элевонов
Мх/1хх
Рисунок 38 - Области располагаемых угловых ускорений в каналах крена и тангажа в случае всех работающих органов управления и при симметричном отказе двух
элементов
Рисунок 39 - Области располагаемых угловых ускорений в каналах крена и тангажа в случае всех работающих органов управления и при несимметричном отказе двух
элементов
Для оценки управляемости БЛА при активном парировании отказов предлагается использовать в качестве критерия значение максимального гарантированного углового ускорения в любом направлении, т.е., максимальный радиус вписанной окружности в многоугольники, отображающие границы располагаемых угловых ускорений, приведенные на рисунках выше.
Следует отметить, что в данной работе приведены результаты для двух проекций момента, поскольку это упрощает отображение на рисунках. В общем случае анализ проводится для всех трех проекций моментов и с учетом возможности раздельного управления тягой двигателей.
Для рассмотренных случаев отказов при значениях эффективностей аэродинамических органов управления, соответствующих данным продувок модели БЛА в АДТ при скорости 50 м/с, получены следующие оценки управляемости:
Число Руль 1 Руль 2 Руль 3 Руль 4 Располагаемое
отказавших Левый Правый Правый Левый угловое
рулей флапперон флапперон элевон элевон ускорение*
рад/с2
0 + + + + 0.74
1 - + + + 0.46
1 + - + + 0.44
2 - - + + 0.31
2 + + - - 0.41
2 - + + - 0.28
2 - + - + 0.25
* значение максимального гарантированного углового ускорения в произвольном направлении, т.е., максимальный радиус вписанной окружности в многоугольники, отображающие границы располагаемых угловых ускорений
8. Выводы
1. На основании проведенного анализа аэродинамической компоновки разработана и программно реализована математическая модель динамики полета преобразуемого БЛА.
2. Предложена и реализована в составе математической модели динамики полета БЛА единая архитектура САУ, обеспечивающая управление преобразуемого БЛА при выполнении типовых сценариев применения, в том числе:
• На режимах малых поступательных скоростей полета
• На режимах вертикального взлёта и посадки
• В крейсерском полёте
• На переходных режимах
3. Проведён анализ и предложены алгоритмы обеспечения управляемости БЛА с применением разработанной системы управления в случае отказов органов управления. В качестве критерия управляемости после отказа предложено использовать значение максимального гарантированного углового ускорения в любом направлении.
Список источников
1. Баженов С.Г., Основы динамики полета. - М.: Физматлит, 2021. - 432 с.
2. Бюшгенс А.Г., Воронин А.Ю., Кувшинов В.М., Леонтьев В.А. Синтез алгоритмов системы управления беспилотного летательного аппарата типа конвертоплан // Ученые записки ЦАГИ. 2018. Т. 49. № 2. С. 39-61.
3. Корнилов В.А., Молодяков Д.С., Синявская Ю.А. Система управления мультикоптером // Труды МАИ. 2012. № 62. URL: http://trudymai.ru/published.php?ID=35543
4. Виноградов С.С. Синтез нечёткого навигационного регулятора для малоразмерного вертолёта «Раптор» // Труды МАИ. 2012. № 73. URL: http://trudymai.ru/published.php?ID=48562
5. Гуськов Ю.П., Загайнов Г.И. Управление полетом самолетов. - М.: Машиностроение, 1991. - 272 с.
6. Ефремов А.В., Захарченко В.Ф., Овчаренко В.Н. и др. Динамика полёта. - М.: Машиностроение-Полет, 2017. - 775 с.
7. Жуматаева Ж.Е. Исследование робастной устойчивости системы управления летательным аппаратом // Труды МАИ. 2012. № 53. URL: http://trudymai.ru/published.php?ID=29614
8. Канатников А.Н. Избранные лекции по алгебре. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2002. - 336 с.
9. Кожевников В.А. Автоматическая стабилизация вертолетов. - М.: Машиностроение, 1977. - 152 с.
10. Комков В.А., Гудзев В.А., Курсаков А.В., Харитонов С.В. Выбор параметров беспилотного летательного аппарата вертикального взлета-посадки дальности полета // Труды МАИ. 2015. № 81. URL: http://trudymai.ru/published.php?ID=57706
11. Крутько П.Д. Обратные задачи динамики в теории автоматического управления. - М: Машиностроение, 2004. - 576 с.
12. Тищенко М.Н., Артамонов Б.Л. Проблемы повышения крейсерской скорости полета вертолета и пути их решения // Труды МАИ. 2012. № 55. URL: http ://trudymai.ru/published.php?ID=30012
13. Dong Il Park et al. Analysis of The Weighting Matrix for Load Distribution Using Weighted Pseudoinverse in a Redundantly Actuated System // ASME 7th Biennial Conference on Engineering Systems Design and Analysis, 2004, DOI:10.1115/ESDA2004-58404
14. Gareth D. Padfield. Helicopter Flight Dynamics Including a Treatment of Tiltrotor Aircraft, Wiley, 2018, 856 p.
15. Jingxian Liao. Mathematical modelling and model predictive controller design of a quad tiltrotor UAV, 2020. DOI:10.1177/0954406220971330
16. Ke Lu et al. Flight Dynamics Modeling and Dynamic Stability Analysis of Tilt-Rotor Aircraft // International Journal of Aerospace Engineering, 2019, no. 2, pp. 1-15. DOI:10.1155/2019/5737212
17. Koji Muraoka et al. Transition Flight of Quad Tilt Wing VTOL UAV, 2012.
18. Kristi Marie Kleinhesselink. Stability and Control Modeling of Tiltrotor Aircraft, 2007.
19. Luis Rodolfo Garcia Carrillo et al. Quad Rotorcraft Control. Vision-Based Hovering and Navigation, 2013.
20. Mark E. Dreier. Introduction to Helicopter and Tiltrotor Simulation, American Institute of Aeronautics and Astronautics, 2007, 637 p.
21. Masayuki Sato. Flight Controller Design and Demonstration of Quad-Tilt-Wing Unmanned Aerial Vehicle // Journal of Guidance, Control, and Dynamics, 2014, no. 38(6), pp. 1-12. DOI:10.2514/1.G000263
22. Mogens Blanke, Jan Lunze, Marcel Staroswiecki. Diagnosis and Fault-Tolerant Control, 2016. DOI:10.1007/978-3-662-47943-8
23. Olena Matiychyk et al. Integration Of Tiltrotor Aircraft Into Modern Air Transport Systems // Conference: "Science - Future of Lithuania. Transport Engineering and Management", 2017. URL: http://jmk.transportas.vgtu.lt/index.php/tran2017/tran2017/paper/view/106
24. Praneet Vayalali et al. Redistributed Pseudoinverse Control Allocation for Actuator Failure on a Compound Helicopter // Conference: The Vertical Flight Society's 76th Annual Forum & Technology Display, 2020.
25. Raza S.J. Use of the Pseudo- Inverse for Design of a Reconfigurable Flight Control System, 1985. DOI:10.2514/6.1985-1900
26. Takahito Mikami et al. Design of Flight Control System for Quad Tilt-Wing UAV, 2015. DOI:10.1109/ICUAS.2015.7152364
27. Thomas Lombaerts et al. Nonlinear Dynamic Inversion based Attitude Control for a Hovering Quad Tiltrotor eVTOL Vehicle // AIAA Scitech 2019 Forum, 2019. DQL10.2514/6.2019-0134
28. Zaibin Chen et al. Design of Flight Control System for a Novel Tilt-Rotor UAV // Complexity, 2020, no.(1-4), pp. 1-14. DQL10.1155/2020/4757381
29. Zhiwei Kong et al. Mathematical Modeling and Modal Switching Control of a Novel Tiltrotor UAV // Journal of Robotics, 2018, no. 3, pp. 1-12. DOI:10.1155/2018/8641731
References
1. Bazhenov S.G. Osnovy dinamiki poleta (Fundametals of flight dynamic), Moscow, Fizmatlit, 2021, 432 p.
2. Byushgens A.G., Voronin A.Yu., Kuvshinov V.M., Leont'ev V.A. Uchenye zapiski TsAGI, 2018, vol. 49, no. 2, pp. 39-61.
3. Kornilov V.A., Molodyakov D.S., Sinyavskaya Yu.A. Trudy MAI, 2012, no. 62. URL: http://trudymai.ru/eng/published.php?ID=35543
4. Vinogradov S.S. Trudy MAI, 2012, no. 73. URL: http://trudymai.ru/eng/published.php?ID=48562
5. Gus'kov Yu.P., Zagainov G.I. Upravlenie poletom samoletov (Aircraft flight control). Moscow, Mashinostroenie, 1991, 272 p.
6. Efremov A.V., Zakharchenko V.F., Ovcharenko V.N. et al. Dinamika poleta (Flight dynamics), Moscow, Mashinostroenie-Polet, 2017, 775 p.
7. Zhumataeva Zh.E. Trudy MAI, 2012, no. 53. URL: http://trudymai.ru/eng/published.php?ID=29614
8. Kanatnikov A.N. Izbrannye lektsii po algebre (Lectures on Algebra), Moscow, Izd-vo MGTU im. N.E. Baumana, 2002, 336 p.
9. Kozhevnikov V.A. Avtomaticheskaya stabilizatsiya vertoletov (Automatic stabilization of helicopters), Moscow, Mashinostroenie, 1977, 152 p.
10. Komkov V.A., Gudzev V.A., Kursakov A.V., Kharitonov S.V. Trudy MAI, 2015, no. 81. URL: http://trudymai.ru/eng/published.php?ID=57706
11. Krut'ko P.D. Obratnye zadachi dinamiki v teorii avtomaticheskogo upravleniya (Inverse problems of dynamics in the theory of automatic control), Moscow, Mashinostroenie, 2004, 576 p.
12. Tishchenko M.N., Artamonov B.L. Trudy MAI, 2012, no. 55. URL: http ://trudymai.ru/eng/published.php?ID=3 0012
13. Dong Il Park et al. Analysis of The Weighting Matrix for Load Distribution Using Weighted Pseudoinverse in a Redundantly Actuated System, ASME 7th Biennial Conference on Engineering Systems Design and Analysis, 2004, D0I:10.1115/ESDA2004-58404
14. Gareth D. Padfield. Helicopter Flight Dynamics Including a Treatment of Tiltrotor Aircraft, Wiley, 2018, 856 p.
15. Jingxian Liao. Mathematical modelling and model predictive controller design of a quadtiltrotor UAV, 2020. D0I:10.1177/0954406220971330
16. Ke Lu et al. Flight Dynamics Modeling and Dynamic Stability Analysis of Tilt-Rotor Aircraft, International Journal of Aerospace Engineering, 2019, no. 2, pp. 1-15. D0I:10.1155/2019/5737212
17. Koji Muraoka et al. Transition Flight of Quad Tilt Wing VTOL UAV, 2012.
18. Kristi Marie Kleinhesselink. Stability and Control Modeling of Tiltrotor Aircraft, 2007.
19. Luis Rodolfo Garcia Carrillo et al. Quad Rotorcraft Control. Vision-Based Hovering and Navigation, 2013.
20. Mark E. Dreier. Introduction to Helicopter and Tiltrotor Simulation, American Institute of Aeronautics and Astronautics, 2007, 637 p.
21. Masayuki Sato. Flight Controller Design and Demonstration of Quad-Tilt-Wing Unmanned Aerial Vehicle, Journal of Guidance, Control, and Dynamics, 2014, no. 38(6), pp. 1-12. D0I:10.2514/1.G000263
22. Mogens Blanke, Jan Lunze, Marcel Staroswiecki. Diagnosis and Fault-Tolerant Control, 2016. DOI: 10.1007/978-3-662-47943-8
23. Olena Matiychyk et al. Integration Of Tiltrotor Aircraft Into Modern Air Transport Systems, Conference: "Science - Future of Lithuania. Transport Engineering and Management", 2017. URL: http://imk.transportas.vgtu.lt/index.php/tran2017/tran2017/paper/view/106
24. Praneet Vayalali et al. Redistributed Pseudoinverse Control Allocation for Actuator Failure on a Compound Helicopter, Conference: The Vertical Flight Society's 76th Annual Forum & Technology Display, 2020.
25. Raza S.J. Use of the Pseudo- Inverse for Design of a Reconfigurable Flight Control System, 1985. D0I:10.2514/6.1985-1900
26. Takahito Mikami et al. Design of Flight Control System for Quad Tilt-Wing UAV, 2015. D0I:10.1109/ICUAS.2015.7152364
27. Thomas Lombaerts et al. Nonlinear Dynamic Inversion based Attitude Control for a Hovering Quad Tiltrotor eVTOL Vehicle, AIAA Scitech 2019 Forum, 2019. D01:10.2514/6.2019-0134
28. Zaibin Chen et al. Design of Flight Control System for a Novel Tilt-Rotor UAV, Complexity, 2020, no. (1-4), pp. 1-14. D0I:10.1155/2020/4757381
29. Zhiwei Kong et al. Mathematical Modeling and Modal Switching Control of a Novel Tiltrotor UAV, Journal of Robotics, 2018, no. 3, pp. 1-12. D0I:10.1155/2018/8641731
Статья поступила в редакцию 22.11.2021; одобрена после рецензирования 29.11.2021; принята к публикации 21.12.2021.
The article was submitted on 22.11.2021; approved after reviewing on 29.11.2021; accepted for publication on 21.11.2021.