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

УДК 629.7.05

СИНТЕЗ И АНАЛИЗ АЛГОРИТМОВ УПРАВЛЕНИЯ БОКОВЫМ ДВИЖЕНИЕМ БЕСПИЛОТНОГО ЛЕТАТЕЛЬНОГО АППАРАТА С ДВУМЯ УПРАВЛЯЮЩИМИ ПОВЕРХНОСТЯМИ

В.М. Глумов, A.M. Пучков, А.Е. Селезнев

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

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

ВВЕДЕНИЕ

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

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

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

Траектория полета БПЛА после отделения его от носителя в верхних слоях атмосферы имеет вид пространственной кривой с большим перепадом высот и содержит криволинейные участки маневрирования в горизонтальной плоскости. При движении по заданной траектории система автоматического управления (САУ) БПЛА должна обеспечивать высокие требования на динамические характеристики движения и точность реализации траектории. Следовательно, отсутствие руля направления требует дополнительного анализа возможностей САУ и модернизации традиционных алгоритмов управления [1, 2].

Частично задача построения подобного автопилота стабилизации углового положения ЛА в

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

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

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

1. ПОСТАНОВКА ЗАДАЧИ

Линеаризованные уравнения, описывающие изменения угловых координат в боковом движении БПЛА при полете на постоянной высоте при малом угле атаки, имеют вид [2]:

= % - Я1„[Р0) + РЛО] + ^(0, ^ = -С21*(0 - *2„[Р(0 + РЛО] -

- Ь21у(?) - ^„Уэ^

- ( -) =

^ = -Ь31%(0 - *3„[Р(0 + М>] -

- С31х(?) - ¿¡А^

где в — угол скольжения, р^ — приращение угла скольжения, возникающее при действии ветра на корпус аппарата, у — угол крена, у — угол рыскания; ах и юу — проекции вектора угловой скорости на оси, соответственно, Ох и Оу связанной с корпусом системы координат; а1„, с1„, с21, а2„, Ь21, ^2„,

31, а3„, с31

с3, и Л„ — динамические коэффициенты,

Ь

характеризующие эффективность управляющих органов, аэродинамическое демпфирование, флюгерную устойчивость, поперечную статическую устойчивость (их значения зависят от аэродинамических характеристик БПЛА, его геометрических, центровочных параметров и массы, от скорости и высоты полета [2, 6]); 8э — воздействие, управляющее изменением угла крена. В БПЛА с двумя управляющими поверхностями воздействие 8э формируется посредством их отклонения в соответствии с выражением 8э = 0,5(8п - 8Л), где 8п и 8Л — отклонения правой и левой поверхностей, согласно которому сигнал управления углом крена подается на приводы.

В рассматриваемом варианте БПЛА управление углом рыскания осуществляется исключительно через управление углом крена. В общем случае традиционный алгоритм управления углом крена ЛА с математической моделью (1) имеет вид [1, 2]

= кДУупр« - У«] + кл^

где к' и к(йх — передаточные коэффициенты; уупр — сигнал управления по крену, формируемый на основе информации о требуемом движении по углу рыскания.

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

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

— процессов стабилизации угла рыскания при действии возмущений р^, вызывающих изменение угла скольжения;

— переходных процессов при координированном управлении БПЛА.

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

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

Для качественного аналитического исследования математической модели (1) приняты допущения:

— ^21 = 0, а20 = 0, т. е. ненулевой угол скольжения и изменения угла рыскания не влияют на изменения угла крена;

— ^30 = 0, т. е. управляющее воздействие 8э не влияет непосредственно на изменение угла рыскания.

При принятых допущениях справедливо утверждение, что с помощью воздействия 8э осуществляется обособленное управление изменением угла крена, который, в свою очередь, полностью определяет изменение углов скольжения и рыскания. В этом случае при анализе и синтезе системы стабилизации угла рыскания оказывается возможным рассматривать математическую модель бокового движения БПЛА в виде последовательного соединения системы стабилизации угла крена, реализующей алгоритм (2), и динамического звена, которое описывает динамику изменения угла рыскания при изменении угла крена. Тогда математическая модель движения БПЛА по координате у описывается передаточной функцией Ж,(5) = у(5)/уупр(^), которая получается из системы уравнений (1) при алгоритме (2) в виде

— (s) B„( s) s—( s):

где K = d20kY — коэффициент передачи в системе стабилизации угла крена; B (s) = s2 + (c21 + kffixd20)s + + K — характеристический полином системы стабилизации угла крена; A^(s) = c31s2 + a10c31s + a30c10;

Bv(s) = ^ + (b31 + a10)s + a30 + b31a10-

При условиях, когда справедливы уравнения (1), всегда выполняются неравенства c10 > 0 и

c31 < 0, так как c31 = —k31 , k31 = const > 0, mz —

спиральный момент рыскания, который традиционно положителен [2]. Для рассматриваемого случая a30 > 0, так как БПЛА обладает флюгерной устойчивостью, и a10 = —k10 Св, k10 = const > 0,

сЦ < 0 [2]. Следовательно, справедливо неравенство а30с10/с31 < 0, которое определяет два действительных нуля разных знаков в передаточной функции № (5). Если учесть, что обычно |а30с10/с31| . |а10|,

то нули функции №^(5) с достаточной точностью вычисляются по приближенной формуле

N

1, 2

±л/-а30С10/c31 .

При большой степени флюгерной устойчивости справедливы неравенства а30 . Ь31, а30 > а10, что позволяет считать а30 + Ь31а10 « а30. Тогда корни уравнения В (5) = 0 определяются выражением

_ b31 + a

10

1, 4

± La

30 .

Корни уравнения В^(^) = 0 определяются параметрами системы стабилизации угла крена:

Л,4 ^ С21 + ± 2 Л С21 + ^30 )2 - 4К .

Коэффициенты к и кшх выбираются из условия, чтобы корни Р3 4 были действительными, отрицательными и различными.

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

3. АНАЛИЗ АЛГОРИТМОВ СТАБИЛИЗАЦИИ УГЛА РЫСКАНИЯ

Наличие нулевого полюса в передаточной функции (3) объекта управления позволяет обеспечить астатизм первого порядка в САУ и делает целесообразным применение пропорционального или пропорционально-дифференциального алгоритма стабилизации угла рыскания.

3.1. Пропорциональный алгоритм стабилизации

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

Yynp ^упр):*

где 1п — коэффициент передачи, у — управляющий сигнал по углу рыскания.

В динамической системе, которая описывается передаточной функцией (3) и алгоритмом (7), из-

Рис. 1. Корневой годограф при пропорциональном алгоритме

менение коэффициента 1п определяет изменение корней ее характеристического уравнения

4 — 2

п (s + P) + in - П (s + N) = 0,

l = 0 c31 j = 1

где Р„ = 0; Р; (/ = 1, 4) — полюсы и Л-, 7 = 1, 2 — нули передаточной функции (3), вычисляемые по формулам (4)—(6). Траектории изменения корней уравнения (8) на комплексной плоскости корней целесообразно и удобно проанализировать с помощью метода корневого годографа [8].

Возможный вид корневого годографа, иллюстрирующий изменение корней уравнения (8) при изменении 1п от 0 до да, приведен на рис. 1.

Отметим, что при 1п = 0 корни характеристического уравнения (8) совпадают с полюсами Р^

I = 0, 4, а при 1п ^ да они совпадают с нулями Л-, 7 = 1, 2. В рассматриваемом случае три траектории корневого годографа уравнения (8) стремятся в бесконечность, а две траектории — к нулям Л-, 7 = 1, 2. Представляет интерес поведение пары комплексных полюсов Р1 2, расположенных близко к мнимой оси. Для определения угла выхода траектории корневого годографа из полюса Р1 составим уравнение фазовых углов для точек, близких к полюсу Р1:

+ 60 ^2 03 ^4— 0,

где 9° 2 — фазовые углы векторов, проведенных к полюсу Р1 из нулей Л-, 7 = 1, 2; 9„, 9г, I = 2, 4, — фазовые углы векторов, проведенных к полюсу Р1 из

полюсов Р„ и Рг, I = 2, 4 . Если учесть условия, при

которых справедливы формулы (4)—(6), то можно считать, что 910 + 92° « 0; 90 « 90°; 92 « 90°; 94 « 0. Тогда 91 « —93, и поскольку |93| < 90°, то и |91| < 90°. Из этого следует, что при увеличении i комплексно сопряженные полюса Р1 2 устремляются в сторону правой полуплоскости к положительному нулю N2, что обусловливает большую колебательность переходных процессов в динамической системе с передаточной функцией (3) и алгоритмом (7) даже при тех значениях in, при которых еще имеет место устойчивость.

3.2. Пропорционально-дифференциальный алгоритм стабилизации

В соответствии с пропорционально-дифференциальным алгоритмом стабилизации угла рыскания сигнал управления по крену формируется в виде

Y — i (ш — ш ) + рв ,

'упр nVT тупр7 rn y>

где рп — коэффициент передачи по угловой скорости рыскания.

Для анализа динамической системы с передаточной функцией (3) с алгоритмом (10) рассмотрим ее характеристическое уравнение в виде

4

I

I = 0

I Bs1 + (PnS + in) I Ajsj — 0,

4

nj j = 0 j

2

где X В/ = (*); I А s- = КД(*).

1 = „ - = „

Уравнение (12) преобразуем к виду

В/ + В/ + (В2 + A2pn)s3 + (В1 + А2/п + АЛ)/ +

+ (В„ + АА + A„Pn)s + А„'п =

Рассмотрим особенности области устойчивости для коэффициентов 1п и рп, применяя метод ^-раз-биения [9]. Заменой s на у'ю получаем параметрическое уравнение границы ^-разбиения:

1

Д

(-A2 в + A0) -A1

2

А1Ю

-A2 в2 + A0

-B4 в4 + B22 - B

42 -B3 В + B1B

2 2 2 где Д(ю) = (-А2ю2 + А„)(-А2ю2 + А„) - А1(-А1ю2).

Отметим, что в данном случае на плоскости (/п, рп) имеется особая прямая 1п = 0. Для определения точек пересечения границы области устой-

2

х

n

X

чивости, определяемой уравнением (12), с осью рп положим /п = 0. Тогда для определения соответствующих значений рп получим из уравнения (12) систему уравнений

В44 - В2Ю + В0 = (А2Ю - ¿0)Pn, В3 Ю2 - В1 = А1Рп.

Система уравнений (13) преобразуется в квадратное уравнение относительно коэффициента рп. Если учесть условия, при которых определялись нули и полюсы передаточной функции (3), то это квадратное уравнение имеет два действительных корня разных знаков.

Для определения точек пересечения с осью /п положим рп = 0. Тогда для определения соответствующих значений получим систему уравнений

В4Ю4 - В2Ю2 + В0

¿1 »п,

В3 Ю4 - В1Ю2 = -(-¿2Ю2 + ¿0) »п.

При решении системы уравнений (14) получается одно кубическое уравнение относительно ю2. Поэтому граница области устойчивости с осью /п может иметь одну или три точки пересечения. Отметим, что на практике в диапазоне частот, наиболее интересном в смысле построения области устойчивости, имеет место одно действительное решение.

Из решений уравнений (13) и (14) следует, что область устойчивости на плоскости (/п, рп) расположена справа от оси рп в левой нижней части I квадранта и в левой верхней части IV квадранта. Таким образом, в алгоритме (10) могут быть использованы как положительные, так и нулевое или отрицательные значения рп.

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

Алгоритм (10) добавляет в разомкнутой системе третий нуль N3 = —/п рп1. Рассмотрим корневой годограф при изменении /п от 0 до при условии фиксированного отношения /п/рп, т. е. при неподвижном нуле N3. Если рп > 0, то нуль N3 лежит в левой полуплоскости, и в уравнении фазовых углов (9) добавится член 93° е [0°; 90°]. Очевидно, 91 е [—90°; 90°]. Это означает, что при увеличении /п ветви корневого годографа из полюсов Р1 2 будут

направлены в сторону правой полуплоскости. По этой причине добиться требуемого качества стабилизации при законе (10) с рп > 0 не представляется возможным.

В случае, если рп < 0, уравнение фазовых углов будет иметь вид

01 + 00 + 03 00 01 02 03 04 = П.

При этом 03 е [90°; 180°]. Рассуждая так же, как и в случае уравнения (9) при пропорциональном алгоритме (7), получаем: 01 « —03 + 0^ — п. Таким образом, чем ближе мы будем приближать нуль N3 к началу координат, тем на больший угол повернутся в отрицательном направлении ветви корневого годографа при выходе из полюсов Р1 2. Это означает, что при увеличении /п при фиксированном отношении /п/рп, обеспечивающем близкое к началу координат расположение нуля N3, комплексно сопряженные полюса Р1 2 будут двигаться вглубь левой полуплоскости.

Отметим, что при фиксированном отношении »п/Рп значение коэффициента /п должна выбираться из условия 0 < /п < »пР , где — значение /п, при котором система (3), (10) находится на границе устойчивости. Значение определяется из совместного решения уравнения границы области устойчивости (12) и уравнения рп = — N3^ /п.

Таким образом, при рп < 0 существует возможность уменьшить степень колебательности системы (3), (10) путем сдвига комплексно сопряженных полюсов влево от мнимой оси.

3.3. Выбор постоянной времени устройства согласования каналов рыскания и крена

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

Для варианта системы с пропорциональным алгоритмом стабилизации (7) и устройством со-

гласования сигнал управления по крену принимает вид

1

упр

T„s + 1

Уупр)'

где Тп — постоянная времени вводимого инерционного звена.

Введение апериодического

звена добавляет дополнительный полюс Р5 = — Тп-1 в разомкнутой системе. Угол выхода корневого годографа из полюса Р1 определяется из уравнения фазовых углов, которое получается из уравнения (9) при добавлении в него фазового угла 95 полюса Р5:

Э^ + 00 00 01 02 03 04 05 = 0.

При условиях, справедливых для формул (4)—(6), как и для уравнения (9), искомый угол выхода: 01 « —03 — 05. Таким образом, введение апериодического звена поворачивает ветви корневого годографа при выходе из полюсов в отрицательном направлении относительно мнимой оси. Выбрав достаточно большую постоянную времени Тп, можно добиться, чтобы при возрастании 1п пара комплексно сопряженных полюсов передаточной функции стала двигаться не в сторону правой полуплоскости, а вглубь левой полуплоскости. Относительный коэффициент затухания, соответствующий этой паре комплексно-сопряженных корней, при росте 1п в формуле (15) будет увеличиваться, что приводит к уменьшению колебательности переходных процессов.

При выборе постоянной времени Тп, чтобы не ухудшить динамические свойства системы, необходимо учитывать, что сопрягающая частота вводимого апериодического фильтра (исопр = 1/Тп) должна в 4+8 раз превышать частоту среза передаточной функции разомкнутой системы (3), которую приближенно можно оценить по формуле: иср = а30с10/(а30 + Ь31«10), откуда получаем условие для выбора значения Тп:

Ри <

a 30 + b 3i а a30c10

10

Для варианта системы с пропорционально-дифференциальным алгоритмом стабилизации (10) и устройством согласования сигнал управления по крену принимает вид

у = —-—[/

упр Т 5 + 1 тупр7 г^у! V /

Из полученных ранее результатов следует, что при алгоритме (17) (для рп < 0) имеет место соотношение: 01 « 0° — 03 — 03 — п. Таким образом, ветвь корневого годографа, выходящая из полюса Р1, будет повернута относительно своего положения на углы —05 (что соответствует введению апериодического звена) и 03 — п (что соответствует введению угловой скорости в алгоритме (17)). При соответствующем выборе данных углов можно обеспечить требуемое качество переходных процессов в канале управления углом рыскания БПЛА рассматриваемого типа.

Отметим, что для случаев, когда —05 = 0° — п (т. е. 1/Тп = — ¿п/рп), корневой годограф будет такой же, как и при алгоритме стабилизации вида

= -P„s + 1 .

Yynp а а + 1 гл(у Уупр^

Динамическое звено с передаточной функцией Wc(s) = (— + 1)/(Tns + 1) представляет собой аппроксимацию первого порядка ряда Паде для звена чистого запаздывания (время запаздывания т = 2Tn) [10]. Следовательно, если технически возможно и целесообразно, то введение чистого запаздывания в пропорциональный алгоритм (7) может положительно сказаться на динамике переходных процессов. Однако это требует проведения дополнительных исследований.

На основе алгоритмов (2) и (17) рассмотрена задача координированного управления БПЛА с двумя аэродинамическими рулями [3].

4. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ АЛГОРИТМОВ СТАБИЛИЗАЦИИ

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

Для оценки влияния полной математической модели на результаты качественного анализа алгоритмов стабилизации угла рыскания рассмотрим систему стабилизации углового положения при боковом движении самолета «Lockheed Jetstar» [11] при отказе руля направления. В этом случае боковое угловое движение самолета описывается математической моделью (1). Предполагается, что в результате данного отказа автопилот стабилизации

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

Инерционно-массовые характеристики самолета в рассматриваемой конфигурации: т =

= 17 330 кг, / = 161 030 кг-м2, / = 330 140 кг-м2,

х у

I, = 184 210 кг-м2, 1ху = 6861 кг-м2. Характеристика рассматриваемого режима: установившийся полет на уровне моря с постоянной скоростью V = 178 м/с с балансировочным углом атаки абал = 2,7°. Динамические коэффициенты ЛА на рассматриваемом режиме: а10 = 0,241 с

-1

c10 = 0,055 с

-1

c21 = 1,799 с

-1

a20 = 9,2 с 2, b21 = 0,178 с 1,

a30 = 6,859 с

-2

¿20 = 5,694 с 2, Ь31 = 0,374 с 1, с31 = -0,203 с-1, ¿30 = 0,088 с-2.

В таблице приведены значения полюсов разомкнутой системы, которая описывается уравнениями (1) и (2) при к = 5, ктх = 2, для трех вариантов расчета:

1) по приближенным формулам (5) и (6) при а20 = 0, = 0, ¿31 = 0;

2) на основе характеристического уравнения при а20 = 0, Ь21 = 0, ^30 = 0;

3) на основе характеристического уравнения при а20 ^ 0, Ь21 ^ 0, ^30 ^ 0.

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

Моделирование системы стабилизации угла рыскания, содержащей объект (1) и канал управления углом крена, работающий по алгоритму (2), проводилось с одним из рассмотренных ранее четырех алгоритмов стабилизации:

— пропорциональный (7) с /п = 3, которому соответствуют полюсы замкнутой системы Р0 = = -0,181; Р12 = -0,118 ± 2,369/; Р3 = -2,369; Р4 = -10,4; '

Расчетные значения полюсов разомкнутой системы

Вариант P1,2 P3 P4

1 0,3075 ± 2,619i -2,72 10,467

2 -0,3075 ± 2,618i -2,72 10,467

3 -0,331 ± 2,616i -2,663 -10,471

Рис. 2. Области устойчивости: - — закон управления без

инерционного звена;-----закон управления с инерционным

звеном

— пропорционально-дифференциальный (10) с /п = 3; рп = -3, которому соответствуют полюсы замкнутой системы Р0 = -0,227; Р12 = -0,417 ± 1,73/; Р3 = -4,7; Р4 = -9,355;

— пропорционально-дифференциальный (17) с /п = 3; рп = -3, в котором значение Тп = 1 с выбрано с учетом условия (16); для рассматриваемого объекта Т < 2,5 с;

п

— пропорциональный с чистым запаздыванием

Y = 3е °'8s(w

'упр v т

^упр^'

Параметры алгоритмов стабилизации i и pn выбраны внутри соответствующих областей устойчивости, которые представлены на рис. 2. Границы областей устойчивости рассчитаны на основе выражений (12)—(14) для принятых значений коэффициентов уравнений (1) и Tn = 1 с. Кроме границ областей устойчивости, на рис. 2 показана прямая линия, соответствующая уравнению pn = — N3-1 in, где N3 = 1 — принятое при моделировании фиксированное положение нуля N3 (см. п. 3.2). Видно, что введение инерционного звена с достаточно большой постоянной времени значительно расширяет область устойчивости. Это позволяет для алгоритма (17) выбрать большие значения коэффициента in, благодаря чему становится возможным повысить быстродействие системы стабилизации. Однако, как видно из рис. 2, диапазон увеличения коэффициента in в алгоритме (17) при сохранении фиксированного значения N3 существенно сокращается по сравнению с вариантом, когда pn = const. Отметим, что максимальная погрешность расчета области устойчивости при a2° = 0, b21 = 0, d3° = 0 составляет около 10 % (например,

Рис. 3. Корневые годографы для доминирующих полюсов замкнутой системы при различных алгоритмах стабилизации бокового движения:

а — алгоритм

при рп = —3 граничное значение 1п = 7 при наличии всех коэффициентов в модели (1) и 1п = 7,8 при а20 = 0, Ь21 = 0, V).

На рис. 3 представлены корневые годографы для полюсов, близких к мнимой оси (доминирующих полюсов), моделируемой системы с указанными ранее четырьмя исследуемыми алгоритмами стабилизации. На траекториях корней черными квадратиками отмечены полюсы замкнутой системы при 1п = 3; рп = —3. Корневой годограф на рис. 3, а подтверждает отмеченный в п. 3.1 факт, что при увеличении коэффициента 1п в пропорциональном алгоритме (7) пара комплексных полюсов Р1 2 замкнутой системы перемещается в сторону правой полуплоскости.

Корневой годограф на рис. 3, б построен при изменении коэффициента 1п с одновременным изменением рп при условии рп/п = —1 (направление изменения 1п и рп показано на рисунке стрелкой) и показывает, что введение в алгоритм (10) сигнала угловой скорости при отрицательном рп и фиксированном положении нуля N изменяет первоначальное направление движения комплексно сопряженных полюсов Р1 2.

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

РА = —1).

При введении в пропорциональный закон чистого запаздывания в виде (18) при увеличении коэффициента 1п комплексно-сопряженные полюса Р1 2 перемещаются вглубь правой полуплоскости (рис. 3, г).

На рис. 4 и 5 показаны осциллограммы реакции системы стабилизации на ступенчатое управляющее воздействие шупр = 1° [?] и на ступенчатый порыв ветра, которому соответствует эквивалентный скачок угла скольжения р^ = 1° [?] соответственно. Из приведенных результатов моделирования видно, что при пропорциональном законе управления наблюдается большая колебательность переходных процессов. Применение алгоритмов стабилизации (10) и (17) позволяет существенно улуч-

V, град,

О 5 10 15 20 25 30

Время, с

Рис. 4. Отработка ступенчатой команды управления при различных алгоритмах стабилизации угла рыскания: —— пропорциональный закон управления; —@--закон управления по рассогласованию и угловой скорости; —д— — закон управления по рассогласованию и угловой скорости с инерционным звеном; — —  — пропорциональный закон управления с запаздыванием

Рис. 5. Отработка ступенчатого порыва ветра при различных алгоритмах стабилизации угла рыскания: —— пропорциональный закон управления; —— закон управления по рассогласованию и угловой скорости; —д— — закон управления по рассогласованию и угловой скорости с инерционным гнездом; — —  — пропорциональный закон с запаздыванием

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

Отметим, что введение в пропорциональный алгоритм чистого запаздывания, равного 0,8 с, позволило существенно улучшить качество переходного процесса. Этот факт, безусловно, интересен. Однако при непосредственной реализации применение алгоритма (17) представляется оправданным и простым в реализации по сравнению с реализацией чистого запаздывания.

ЗАКЛЮЧЕНИЕ

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

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

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

улучшения динамики переходного процесса при введении в пропорциональный алгоритм управления вида (7) чистого запаздывания.

ЛИТЕРАТУРА

1. Сыров А. С., Пучков А.М., Рутковский В.Ю., Глумов В.М. Задачи управления движением многорежимных беспилотных летательных аппаратов // Проблемы управления. — 2014. — № 4. — С. 45—52.

2. Боднер В.А. Теория автоматического управления полетом. — М.: Наука, 1964. — 608 с.

3. Сыров А.С., Пучков А.М., Селезнев А.Е., Глумов В.М. Алгоритмы модернизированного координированного управления беспилотным летательным аппаратом // Тр. XII Все-рос. совещания по проблемам управления (16—19 июня 2014 г., Москва, ИПУ РАН). — М., 2014. — С. 3407—3416.

4. Пат. № 2532720 РФ, кл. G05D. Двухканальное устройство координированного управления / А.С. Сыров, А.М. Пучков // Опубл. 10.11.2014, Бюл. № 31.

5. Основы формирования облика систем управления авиационного ракетного вооружения / Э.М. Абадеев, Ю.П. Балы-ко, В.В. Ляпунов и др.; под ред. В.В. Трусова. — М.: Изд.-торговая корпорация «Дашков и К°», 2012. — 176 с.

6. Колесников A.A. Новые нелинейные методы управления полетом. — М.: Физматлит, 2013. — 196 с.

7. Боднер В.А., Козлов М.С. Стабилизация летательных аппаратов и автопилоты. — М.: Оборонгиз, 1961. — 508 с.

8. Удерман Э.Г. Метод корневого годографа в теории автоматических систем. — М.: Наука, 1972. — 448 с.

9. Красовский A.A., Поспелов Г.С. Основы автоматики и технической кибернетики. — М.—Л.: ГЭИ, 1962. — 968 с.

10. Титов НИ, Успенский В.К. Моделирование систем с запаздыванием. — Л.: Энергия, 1969. — 97 с.

11. Heffley R., Jewell W. Aircraft handling qualities data / NASA Contractor Report-2144. December 1972. — P. 166—193.

Статья представлена к публикации членом редколлегии Б.В. Павловым.

Глумов Виктор Михайлович — д-р техн. наук, вед. науч. сотрудник, Институт проблем управления им. В.А. Трапезникова РАН, г. Москва,  (495) 334-87-79, И vglum@ipu.ru,

Пучков Александр Михайлович — д-р техн. наук, нач. лаборатории,

Московское опытно-конструкторское бюро «Марс»,  (495) 684-44-56, И office@mokb-mars.ru,

Селезнев Антон Евгеньевич — аспирант, инженер 1 кат., Московское опытно-конструкторское бюро «Марс»,  (495) 684-44-56, И office@mokb-mars.ru.

W

Не забудьте подписаться!

Подписку на журнал «Проблемы управления» можно оформить в любом почтовом отделении (подписной индекс 81708 в каталоге Роспечати или 38006 в объединенном каталоге «Пресса России»), а также через редакцию с любого месяца, при этом почтовые расходы редакция берет на себя. Отдельные номера редакция высылает по первому требованию.