УДК 533.6.013
РЕАЛИЗАЦИЯ ТРАЕКТОРНОГО УПРАВЛЕНИЯ В ГОРИЗОНТАЛЬНОЙ ПЛОСКОСТИ С АДАПТИВНЫМ АЛГОРИТМОМ УПРАВЛЕНИЯ ПОПЕРЕЧНЫМ ДВИЖЕНИЕМ САМОЛЕТА
Д.В. ВЕРЕЩИКОВ, П.С. КОСТИН, Ю.О. ВЕРЕЩАГИН
Статья представлена доктором технических наук, профессором Левицким С.В.
В статье представлены результаты исследований по оценке эффективности применения адаптивного алгоритма управления маневренным самолетом, реализованного в поперечном канале управления, на основе идентификации ряда аэродинамических характеристик самолета и возмущений, вызываемых образованием несимметричной конфигурации размещения внешних подвесок. Путем моделирования движения самолета в ходе выполнения S-образных маневров с использованием программно-моделирующего комплекса и пилотажного стенда продемонстрировано существенное улучшение характеристик управляемости несимметричного самолета и, таким образом, повышение эффективности его применения.
Ключевые слова: устойчивость, управляемость, адаптация, идентификация, маневрирование.
Одной из главных задач при создании и модернизации современных маневренных самолетов является обеспечение заданных пилотажных характеристик или характеристик устойчивости и управляемости (ХУУ) во всем эксплуатационном диапазоне полета. Эта задача в последнее время стала более трудоемкой в силу следующих причин [1]: сложной аэродинамической компоновки современных маневренных самолетов; значительного расширения эксплуатационного диапазона для обеспечения многофункциональности (что является одним из основных требований к современному самолету); создания большого количества серий и модификаций базового самолета; необходимости обеспечения заданного уровня безопасности полетов [2], в тот момент, когда необходимо добиться компромисса между максимальным использованием возможностей самолета и наличием эксплуатационных ограничений. Все это заставляет авиационных специалистов разрабатывать и применять новые концептуальные решения для обеспечения заданных ХУУ.
Новым решением при синтезе систем управления (СУ) современных маневренных самолетов является реализация адаптивных алгоритмов на основе идентификации нескольких аэродинамических параметров в режиме реального времени - «on-line». Такие алгоритмы в настоящее время применяются на самолетах Су-30МКИ [3], Су-30СМ, Як-130, Су-35, Т-50. Необходимо заметить, что реализация этих алгоритмов на борту самолета стала возможна благодаря широкому использованию в системах управления цифровых вычислителей, которые обеспечивают необходимый объем вычислений с приемлемой скоростью.
Анализ СУ современных самолетов показал, что адаптивные алгоритмы на основе идентификации ряда аэродинамических характеристик применяются только для обеспечения заданной динамики продольного движения. При этом в боковом канале управления обеспечение заданной динамики при действии возмущений различной природы не предусмотрено. Решение этой задачи является весьма актуальным для современных маневренных самолетов, у которых при энергичном пространственном маневрировании возможно проявление существенной взаимосвязи продольного и бокового движений, приводящей к неблагоприятным последствиям [4; 5]. Ситуация еще более осложняется, когда на самолет при этом действуют различные возмущения, например от несимметричной внешней подвески, атмосферных возмущений, отказов и неисправностей СУ.
Авторами было принято решение провести исследования по обеспечению заданных динамики и устойчивости поперечного движения самолета путем использования адаптивного алго-
ритма управления на основе идентификации аэродинамических характеристик поперечного движения. Для этого был создан программно-моделирующий комплекс, основными элементами которого являются пилотажный стенд маневренного самолета [6] и математическая модель динамики пространственного движения самолета с комплексной системой управления, близкого по своим аэродинамическим, массовым и геометрическим характеристикам самолету Як-130.
В качестве алгоритма идентификации был принят рекуррентный метод наименьших квадратов. При этом рассматривается модель поперечного движения в следующем виде
ЬЩг+^ - 1у ¡уГ + 1у (¡2Г -1^ - ^ШГ- ) = + (1)
+М ¡Г + М Г Г + М М X, где 1Х, 1у, Iг - осевые моменты инерции; I , I хг, 1гу - центробежные моменты инерции; Г, ¡у, ¡г - скорости крена, рыскания, тангажа соответственно; М3Х, М3Х - производные
момента крена, характеризующие отклонение рулевых поверхностей; Мх , М°Г, М¡у -производные момента крена, характеризующие поперечную статическую устойчивость, демпфирование и перекрестную связь; 8Э, 5Н - углы отклонения элеронов и руля направления соответственно; в - угол скольжения.
Инерционные перекрестные связи, отражающие влияние на динамику вращательного движения самолета относительно центра масс инерционных моментов представляются следующим образом
М кр = (1г - 1у ¡¡у - 1у (Г2-¡¡) +1 +I хгЩГ - Г (I уГг -1 ГУ). (2)
В целях увеличения сходимости оценок параметров при идентификации выражение (1) с учетом выражения (2) преобразуется к виду (3) и, таким образом, исключается инерционная составляющая поперечного движения
Г = Мкр_ = 1 (МЧ + МГг + МГуг + М вв+МХ). (3)
ш м1х1х
В обобщенном (векторном) виде рассматриваемая модель выглядит следующим образом
х = А(х, г )х + В(х, г )и + (4)
где х - п -мерный вектор состояния; и - т -мерный вектор управления; £ - вектор внешних
неконтролируемых возмущений; А, В - матрицы параметров объекта управления соответствующей размерности с элементами, являющимися функцией времени и состояния объекта; г - текущее время.
Оценка (идентификация) неизвестных параметров матриц А и В осуществляется с использованием следующего алгоритма [7; 8]:
1) е(г) = Г (г) - г (г -1), где е(г) - ошибка на текущее время; Р(г )и (г +1) Л + иТ (г +1) Р(г )и (г +1)
ния; Р(г) - обратная корреляционная матрица; и(г) - вектор входных сигналов; Л - коэффициент забывания;
2) к (г +1) = -—Т/ '
—-— , где к (г) - вектор столбец коэффициентов усиле-
Р(г +1) =1 [р(г) - к (г + 1)иТ (г +1) Р(г)];
3) Р(г +1) = —1Ч"Т'
4) у (г +1) = у (г) + К (г + 1)е(г), где у (г) - вектор оценок неизвестных параметров на г -такте просчета.
Адаптация поперечного движения самолета осуществляется путем определения расчетного (потребного) угла отклонения элеронов (3Э)расчет для обеспечения модельной скорости крена
а>хмод, которая может задаваться различными способами
, ч р-(й*-Р-^)++ мвв+м:*5н)
(3Э) =-,-, (5)
V э / расчет -о 5 ^ '
М 3■
" x
где M®x, Mx, Mв, Mх", AMx , M хэ - полученные оценки аэродинамических производных крена, отнесенные к осевому моменту инерции Ix .
После определения (8Э )расчет его значение сравниватеся с текущим 8Э, и, если они не равны, формируется корректирующий сигнал, направленный на устранение этой разницы
u = S -(S ) (6)
кор э \ э/расчет ' \ /
В процессе проведения исследований работоспособности представленного адаптивного алгоритма моделировался полет с несимметричной внешней подвеской, которая вносит определенные коррективы в динамику полета самолета, при этом изменяются моменты инерции и дополнительно создаются моменты от весовой несбалансированности и интерференционного взаимодействия крыла и подвески. При наличии несимметричной внешней подвески выражения (1)-
(3), (5) изменятся, в них появляются дополнительные члены ^(р-подеески rz_ )ny cosy - момент от
i
весовой несбалансированности, AMx - интерференционный момент от аэродинамического взаимодействия подвески и крыла самолета:
ГТ т \ т г 2 2S т day т daz
+ (Iz - Iy )VyVz + Izy (®z -Юу) - Iy—jy - ^ + (Iyx^z - Izxay
, (7)
= MSs S + МЮЮ + + Мв + M?SH +X (GnodeecKUrZi )ny cosy + AMx_p
i
где GmdeecKU - вес несимметричной i -й подвески; rz - расстояние по оси OZ связанной системы координат до центра масс несимметричной i -й подвески; у - угол крена; ny - нормальная
перегрузка;
M =Y (G е r )n cos y - (I -1 )a a -1 (a2-a2) +1 +1 -a (I a -1 a ); (8)
кп /i\ подеескиг xx' y / \ z y/ z y zy\ z y/ xy ^ xz ^ x \ xy z xz y s
da' da M Kp 1
(MSS S + M a®x + M a ay + Me в + M SS"SH + AM ); (9)
1. 1. T T \ x э x x x y x ' x
dt dt Ix Ix
a - (ax д)+M У ay + Mi в в + Ms"SH + AMx + M р)
J. X x \ xмод! X y x Г x н xu кр)
(S) =-dt-,-. (10)
V э J расчет s v '
м 3
Структурная схема адаптивного алгоритма на основе идентификации аэродинамических параметром в полете с несимметрично расположенными внешними подвесками представлена на рис. 1.
Рис. 1. Структурная схема адаптивного алгоритма управления поперечным движением самолета
Для оценки эффективности применения адаптивного алгоритма выполнено моделирование Б - образного маневра, применяемого для устранения рассогласования по линии заданного пути (ЛЗП) при решении задач прицеливания, дозаправки в воздухе, полета строем, движения при заходе на посадку. При этом рассматривалось движение симметричного с точки зрения размещения внешних подвесок самолета и при наличии несимметричной внешней подвески. С методической точки зрения моделировалось движение самолета как с идеальным регулятором [9], так и с учетом запаздывания, вносимого летчиком. Закон работы идеального регулятора описывается выражением
к /
у-
V
к; (г - г Зад)+к^¥-¥зад)+к; —
+к?юх, (11)
где к; (; - ;Зад ) + к^(¥-¥зад ) + ку~ = У зад - заданный крен; ; , ;зад - текущая и заданная
Ж;
Ж;
Ж
координаты по оси О;нормальной земной системы координат; у , Щзад - текущий и заданный курсы самолета; кгэ , к2; , к^ , к; , к- коэффициенты усиления при соответствующих сигналах.
Учет запаздывания, вносимого в управление летчиком, выполнялся путем внесения в контур управления передаточной функции
Жл (р) = к ле , (12)
где Т - эффективное запаздывание летчика; кл - коэффициент усиление летчика [10].
Результаты моделирования представлены на рис. 2-4. Моделирование выполнялось при Т = 0,2 , так как при Т < 0,2 летчик оценивает объект как плохой в управлении, а при Т > 0,25 появляется раскачка при точном выдерживании одной из фазовых координат движения (в данном случае координаты ; ). Высота и скорость движения самолета выдерживались постоянными Н = 2000м, V = 400 км/ч .
апшс/рнчнип салюлег
жшмефй^ш самолет
чиый самолет ' ал/орнгао.и
с алагггмэшм а.']
сникефнчнмн самолет
нешд1е1ричныйв||1вш
несд>шсгаи*шыЯ самолет с йЛ$ПШВШМ алгоритмом
а
б
Рис. 2. Выход на ЛЗП при А; = 50м : а - идеальный регулятор; б - регулятор с учетом свойств летчика
симметричный самолет
несимметричный самолет •
несшшетршшшг еа.иолет : с ддапшв/шм алгоритмом
симметричный самолет
несиммстричный самолет
несимметричный самолет с аддитивным алгоритмом
б
Рис. 3. Выход на ЛЗП при А; = 25м : а - идеальный регулятор; б - регулятор с учетом свойств летчика
20
а ' б
Рис. 4. Выход на ЛЗП при А; = 10м : а - идеальный регулятор; б - регулятор с учетом свойств летчика
Анализ результатов моделирования пространственного движения маневренного дозвукового самолета при наличии несимметричной внешней подвески дает основания утверждать, что без адаптивного алгоритма и идеальный регулятор, и тем более летчик тратят значительно больше времени на вывод самолета на ЛЗП. При этом имеют место значительные колебания по углу крена. Это вызвано разной реакцией самолета на отклонение ручки управления и, как следствие, элеронов в разные стороны.
Применение адаптивного алгоритма как с идеальным регулятором, так и в случае запаздывания, вносимого летчиком, в значительной степени улучшает качество управления как по времени выхода на ЛЗП, так и по характеру изменения угла крена.
Таким образом, представленный алгоритм автоматизации управления обеспечивает улучшение характеристик управляемости самолета в поперечном канале управления при наличии возмущений от несимметричных внешних подвесок, приближая его с точки зрения управляемости к симметричной конфигурации. Это позволяет сократить время на реализацию Б - образных маневров и повысить точность выдерживания траектории, что в свою очередь приводит к повышению эффективности применения самолета.
ЛИТЕРАТУРА
1. Балык О.А., Москалев П.Б. Реализация принципов адаптивного управления самолетом // Научно- технический сборник 929 ГЛИЦ им. В.П. Чкалова. - 2006. - № 4. - С. 67-74.
2. Кротов С.А. К вопросу о контроле отказобезопасности функциональных систем воздушных судов в процессе эксплуатации // Научный Вестник МГТУ ГА. - 2013. - № 197. - С. 79-84.
3. Москалев П.Б. Применение самонастраивающихся алгоритмов в системе управления маневренных самолетов: сб. трудов междунар. научн.-технич. симпозиума в рамках международного авиасалона МАКС - ЦАГИ. Жуковский, 2001.
4. Левицкий С.В., Свиридов Н.А. Динамика полета. - М.: ВВИА им. проф. Н.Е. Жуковского, 2008.
5. Верещиков Д.В., Салтыков С.Н. Особенности аэродинамики и динамики полета самолета с внешними подвесками при сверхзвуковых скоростях полета. - Воронеж: ВАИУ, 2009.
6. Верещагин Ю.О., Костин П.С., Подкуйко Т.А., Верещиков Д.В. Применение полунатурного моделирования для исследования пилотажных характеристик самолетов на кафедре авиационных комплексов и конструкции летательных аппаратов // Вестник Воронежского государственного технического университета. - 2011. - № 11.2. - Т. 7. - С. 12-14.
7. Александров А.Г. Оптимальные и адаптивные системы. - М.: Высшая школа, 1989.
8. Сергиенко А.Б. Алгоритмы адаптивной фильтрации: особенности реализации в МАТЬАБ // Среда разработки. - 2003. - № 1. - С. 18-28.
9. Вавилов Ю.А. Системы автоматического управления полетом. - М.: ВВИА им. проф. Н.Е. Жуковского, 2009.
10. Лысенко Н.М. Системы управления и бортовые цифровые вычислительные комплексы летательных аппаратов / под ред. Н.М. Лысенко. - М.: ВВИА им. проф. Н.Е. Жуковского, 1990.
REALISATION TRAJECTORY CONTROL IN A HORIZONTAL WITH ADAPTIVE ALGORITHM OF MANAGEMENT OF THWART MOTION OF THE PLANE
Vereshikov D.V., Kostin P.S., Vereshagin Yu.O.
In article results of researches are presented according to efficiency of application of adaptive algorithm of management by the maneuverable plane realised in a cross-section control path, on the basis of identification of some aerodynamic characteristics of the plane and the indignations caused by formation by an asymmetrical configuration of placing of external suspension. By modelling of movement of the plane during performance S - figurative maneuvers with use of a pro-grammno-modelling complex and the flight stand substantial improvement of characteristics of controllability of the asymmetrical plane and, thus, increase of efficiency of its application is shown.
Keywords: stability, controllability, adaptation, identification, manoeuvring.
Сведения об авторах
Верещиков Дмитрий Викторович, 1970 г.р., окончил ВАИУ (г. Иркутск) (1993), кандидат технических наук, доцент, начальник кафедры авиационных комплексов и конструкции летательных аппаратов ВУНЦ ВВС «ВВА» (г. Воронеж), автор более 150 научных работ, область научных интересов - динамика полета, моделирование авиационных комплексов.
Костин Павел Сергеевич, 1985 г.р., окончил ВАИУ (г. Воронеж) (2012), адъюнкт кафедры авиационных комплексов и конструкции летательных аппаратов ВУНЦ ВВС «ВВА» (г. Воронеж), автор 12 научных работ, область научных интересов - системы управления летательных аппаратов.
Верещагин Юрий Олегович, 1989 г.р., окончил ВАИУ (г. Воронеж) (2011), адъюнкт кафедры авиационных комплексов и конструкции летательных аппаратов ВУНЦ ВВС «ВВА» (г. Воронеж), автор 16 научных работ, область научных интересов - системы управления летательных аппаратов, моделирование авиационных комплексов.