CC BY
139
Доклады БГУИР
Doklady BGUIR
2017, № 3 (105) 2017, No. 3 (105)
УДК 629.13.014-506.4
РОБАСТНЫЙ АВТОПИЛОТ КАНАЛА ТАНГАЖА БЕСПИЛОТНОГО ЛЕТАТЕЛЬНОГО АППАРАТА
Ю.В. ГРИДНЕВ, А.Г. ИВАНОВ
Физико-технический институт Национальной академии наук Беларуси,
Республика Беларусь
Объединенный институт проблем информатики Национальной академии наук Беларуси,
Республика Беларусь
Поступила в редакцию 17 марта 2017
Аннотация. Рассматривается система автоматического управления по тангажу при полете беспилотных летательных аппаратов (БЛА) с различной аэродинамической схемой и конфигурацией планера. Для такой системы предлагается робастная схема автопилота, которая позволяет при изменении параметров канала управления планера самолета подстраивать параметры автопилота так, чтобы динамические и переходные ошибки управления были постоянными при смене планера БЛА. Компьютерное моделирование робастного автопилота подтвердило правильность теории робастных систем.
Ключевые слова: робастный автопилот канала тангажа, математическая модель, моделирование.
Abstract. Automatic pitch control system of unmanned aircraft vehicles (UAVs) with different aerodynamic configuration and airframe construction is reflected here. Robust autopilot scheme is proposed for such control system. While adjusting the parameters of UAV's control channel, robust autopilot scheme allows presetting autopilot the way the dynamic and transient errors of control system stay invariant under UAV changing. Computer modeling of robust autopilot confirms theory of the robust systems.
Keywords: robust autopilot of the pitch channel, mathematical model, modeling.
Doklady BGUIR. 2017, Vol. 105, No. 3, pp. 40-44
Robust autopilot of the pitch channel of unmanned aircraft vehicles
Yu.V. Gridnev, A.G. Ivanov
Современный этап эволюционного развития авиации характеризуется стремительным развитием беспилотных авиационных комплексов (БАК), в состав которых входят беспилотный летательный аппарат (БЛА), наземный пункт управления (НПУ) и радиоканал связи между НПУ и БЛА. Интерес к БАК обусловлен их применением для решения гражданских и военных целевых задач, а также новыми достижениями микросистемной авионики.
Управление полетом БЛА осуществляется пятью контурами систем автоматического управления (САУ), основными элементами которых являются последовательно включенные автопилот (АП) и планер БЛА. Количественные оценки управления полетом самолета и его контуров САУ можно определить, если известны передаточные характеристики АП и планера БЛА. Математическая модель пленера БЛА как объекта управления описывается его движением центра масс и вращением вокруг него. Каждая из этих двух составляющих обладает тремя степенями свободы и управляется силами тяги и рулями БЛА. Уравнение продольного движения БЛА в вертикальной плоскости описывается силами и моментами (рис. 1)
Рис. 1. Силы и моменты, действующие на планер БЛА
К центру масс БЛА (ЦМ) приложены сила тяжести О, сила тяги двигателя Р и полная аэродинамическая сила Я, возникающая при взаимодействии БЛА с набегающим потоком воздуха. Управление полетом БЛА при его продольном движении осуществляется с помощью изменения силы тяги двигателя Р, совпадающей по направлению с продольной осью БЛА и положением руля высоты (угла ЪРВ). Изменение угла ЪРВ поворачивает самолет в вертикальной плоскости вокруг оси 02 за счет появления вращающего момента. Кроме того, на БЛА действуют сила тяжести О, направленная вертикально, и полная аэродинамическая
сила Я = У + Q, как равнодействующая всех аэродинамических сил, которая определяется
подъемной силой У, и силой лобового сопротивления Q . Подъемная сила, как составляющая полной аэродинамической силы, перпендикулярна вектору скорости БЛА, и определяется по формуле
Р V2
Р V2
= с V^ = с/
-Бс1а = У аа .
(1)
где Су = Суаа - коэффициент подъемной силы, определяющийся через его производную Суа и
угол атаки а ; д =
рИ 2
скоростной напор воздуха с учетом плотности воздуха и скорости
полета; ^ - площадь поверхности крыла и части фюзеляжа БЛА, которые и создают подъемную силу.
Управляющая сила руля высоты приложена к фокусу крепления руля, расположенного на расстоянии хрв позади ЦМ, и определяется как
У = С
Р V2
$РВ = Су
Р V2
ирв ирв у
(2)
2 у 2
Учитывая расстояние х^ между ЦМ и ЦД в процессе полета БЛА возникает общий
аэродинамический момент М, который проецируется на оси 0ХУ2 связанной системы координат:
М =Мх +Му +М-. (3)
Момент тангажа М2 вращает БЛА вокруг оси 02 в вертикальной плоскости и зависит от угла поворота руля высоты Ърв, угла атаки а и угловой скорости вращения БЛА в вертикальной плоскости гаг. С учетом направления вращения БЛА вокруг оси 02, момент тангажа М может быть представлен тремя слагаемыми:
м = м62ъРВ -м> -М;ю ,
(4)
где Мг,Ма,М; - частные производные от момента тангажа по соответствующим аргументам.
Первое слагаемое выражения (4) является управляющим моментом и зависит от угла поворота руля 8РВ и расстояния хрв от центра масс до фокуса руля:
М2упр = м%в = С/ Р—2$РВ • хрв • 5РВ . (5)
Второе слагаемое выражения (6) является стабилизирующим моментом и зависит от угла атаки а и расстояния хР от ЦД до ЦМ:
Р —2
Мгстаб = М> = хр -а . (6)
Третье слагаемое является демпфирующим моментом и представляет собой момент сопротивления, который возникает при вращении БЛА в потоке воздуха:
М = Мш• ш = С ш• Р—- Б—ш . (7)
гдемп 2 2 у ^ у 2
Пространственные движения БЛА определяются уравнением сил
=тЖу «7аа + Ра±758га, уравнением моментов ./ГЭ = АС Ъ,,,, -М^и-М'^ок и
уравнениями кинематики полета (уравнением связи углов & = 9 + а и скоростным
Гхд =КсО8 0 уравнением <
[уя=Г* щ9
Для получения передаточной функции канала тангажа планера БЛА по углу атаки а необходимо получить зависимость а = /(8РВ) с учетом динамики вращения БЛА вокруг оси а( р)
OZ, т.е. К(р) =-. Для этого используем уравнения сил, моментов и уравнений кинематики
5 (р)
■ 7а + р
с учетом предварительных преобразований Ш = V ■ 9 , и 1<\ = тУ9 = (¥"' + р)а,9 =-а ,
тУ
где Тэ --аэродинамическая постоянная времени БЛА.
С учетом сказанного выше основное уравнение динамики вращения БЛА вокруг оси 02 можно записать в виде
^ф + а)=М*5-М?а-М?ф + а), (8)
которое позволяет, используя преобразования Лапласа, получить передаточная характеристика канала тангажа САУ по углу атаки а:
Ка (Р) =аР =-. (9)
0( Р) п а2 — 2 2\
а (1+ р + р )
а а1
Выражение (9) показывает, что передаточная характеристика канала тангажа планера БЛА соответствует колебательному звену:
Ка(Р) = л^\ Т2 2 . (Ю)
1 + 2с,Тр + Т р
„ М0 М0 М\Т, М0 ^х. где К = —- =-2-=-—-» —- = 0 0 - коэффициент передачи канала,
а а Ма + М1_ М^Т, + МШг Ма БаХр фф р ,
который определяется отношением площади руля к площади крыла с учетом их плеч
относительно ЦМ; Т = I— =
V а
1 1
— — —
-2-« —— =-2--постоянная времени;
М а + МШ_ V Ма V су ^ Хр у
2 Т
Е « ^ + Р \—— коэффициент демпфирования. 2тУ
В итоге общая передаточная характеристика планера БЛА канала тангажа с учетом связи углов определяется выражением:
К (1+рт■)
к& (р) =
5(р) р(1 + 2%Тр + Т2 р2)' К,
(11)
3 р, = Кр АЗ
где К (р)= —~ - коэффициент передачи канала тангажа.
Тэ
Второй подсистемой САУ является автопилот, управляющее воздействие которого в канале тангажа формируется на основе сигнала ошибки, представляющего собой разность АЗ = Зз — & между задаваемым углом тангажа & и измеренным углом 0. После преобразования сигнала АЗ в различных устройствах АП на его выходе формируется сигнал управления рулевого привода
^АЗ — КЛ, (12)
р
где р - оператор Лапласа; - угловая скорость ЛА вокруг оси г; Кр,Ки Кш - коэффициент передачи автопилота по положению, по скорости полета и по угловой скорости.
Для конкретного типа планера БЛА с его заданными передаточными характеристиками выбирают в канале тангажа оптимальные коэффициенты Кр, К, Кш автопилота, позволяющие повысить точностные характеристики канала тангажа САУ. Изменение параметров планера БЛА (при его замене) приводит к ухудшению показателей качества характеристик контура САУ, т.к. данный АП в измененном контуре САУ не будет оптимальным.
Разрешить проблему сохранения оптимальности САУ возможно за счет применения робастного автопилота, который является универсальным для любых планеров БЛА с весом от 0,5 кг до 50 кг. Синтез структуры робастного автопилота канала тангажа производится в соответствии с концепцией обратных задач динамики [1], которая позволяет определить управляющее воздействие на входе сервопривода автопилота. Для робастного АП выходной сигнал у^) описывается дифференциальным уравнением третьего порядка [1]:
а3у + а2у + а1у + а0у = а0хвк, (13)
где №>, а\. а2. аз - коэффициенты уравнения; у,у,у - третья, вторая и первая производные
сигнала на выходе АП; хвх - сигнал на входе АП.
Коэффициенты а0, а, а выбираются исходя из требуемых характеристик переходного процесса. Согласно обратной задачи динамики, сигнал управления на входе сервопривода определяется видом
м(/) = к «01(кгхВК - аха - а2а - а
(14)
Математическая модель предлагаемого робастного автопилота показана на рис. 2.
Рис. 2. Робастный АП в САУ канала тангажа
Для проверки работоспособности и характеристик робастного автопилота канала тангажа разработана компьютерная модель, результаты моделирования которой показаны на рис. 3 и 4.
а б
Рис. 3. Переходные процессы при К = 4: а - классический АП; б - робастный АП
а б
Рис. 4. Переходные процессы при К = 8: а - классический АП; б - робастный АП
Из рисунков видно, что при увеличении коэффициента К в 2 раза, классический АП на
грани возбуждения, а робастный АП сохранил характеристики канала тангажа.
Список литературы
1. Крутько П.Д. Обратные задачи динамики в теории автоматического управления. Цикл лекций. М: Машиностроение, 2004. 576 с.
2. Малкин В.А., Гриднев Ю.В., Пальцев А.Н., Цанава А.А. Робастный автопилот канала тангажа / Патент РБ № 8404.
3. Малкин В.А., Гриднев Ю.В., Пальцев А.Н., Цанава А.А. Робастный автопилот канала тангажа летательного аппарата / Патент № 18251.
References
1. Krut'ko P.D. Obratnye zadachi dinamiki v teorii avtomaticheskogo upravlenija. Cikl lekcij. M: Mashinostroenie, 2004. 576 s. (in Russ.)
2. Malkin V.A., Gridnev Ju.V., Pal'cev A.N., Canava A.A. Robastnyj avtopilot kanala tangazha / Patent RB № 8404. (in Russ.)
3. Malkin V.A., Gridnev Ju.V., Pal'cev A.N., Canava A.A. Robastnyj avtopilot kanala tangazha letatel'nogo apparata / Patent № 18251. (in Russ.)
Сведения об авторах
Information about the authors
Гриднев Ю.В., к.т.н., профессор, ведущий научный сотрудник ФТИ НАН Беларуси.
Gridnev Yu.V., PhD, professor, leading researcher of PTI NAS Belarus.
Иванов А.Г., аспирант ОИПИ НАН Беларуси.
Ivanov A.G., PG student of UIIP NAS Belarus.
Адрес для корреспонденции
220012, Республика Беларусь г. Минск, ул. Сурганова, д. 6, ОИПИ НАН Беларуси тел. +375-17-284-20-86; e-mail: [email protected] Иванов Александр Геннадьевич
Address for correspondence
220012, Republic of Belarus, Minsk, Surganova st., 6, UIIP NAS Belarus tel. +375-17-284-20-86; e-mail: [email protected] Ivanov Alexander Gennad'evich
