CC BY
194
№ 5 (53) 2014
А. В. Людомир, техник Государственного научно-исследовательского института авиационных систем, г. Москва, gluek@bk.ru
В. С. Орлов, начальник сектора Государственного научно-исследовательского института авиационных систем, г. Москва, orlov@gosniias.ru
имитационное моделирование динамической воздушной обстановки в управляемом воздушном пространстве
Рассмотрено имитационное моделирование потоков воздушных судов в системе управления воздушным движением. Приведена структура модели воздушного судна, указаны соотношения учета ЛТХ в соответствии с базой данных, публикуемой организацией Евроконтроль, произведено обоснование модели ошибок самолетовождения, базирующейся на концепции «навигации, основанной на характеристиках» PBN (Performance Based Navigation). Рассмотрен пример самолетовождения с обработкой цифровых команд от модели диспетчера управления воздушным движением.
Ключевые слова: моделирование, управление воздушным движением, PBN, CPDLC.
введение
Для решения многих задач исследования и оценки безопасности полетов применяется подход на основе имитационного моделирования потоков воздушных судов (ВС) в сегрегированном воздушном пространстве под контролем системы управления воздушным движением (УВД). Движение потока ВС реализуется в составе модели динамической воздушной обстановки (МДВО). Условия, в которых проходят полеты, динамически изменяются непосредственно в процессе моделирования. Это касается прежде всего состояния метеорологической обстановки и состояния структуры воздушного пространства (динамически изменяющихся ограничений использования воздушного пространства (ИВП), изменения структуры секторов и пр.).
МДВО применяется в программных средствах, где требуется имитация воздушного движения. С использованием МДВО решаются задачи оценки безопасности полетов, проводится моделирование взаимодейст-
вия между наземными службами УВД и бортом, проводится разработка и анализ бортовых алгоритмов, в том числе обнаружения и предотвращения опасных сближений в воздухе.
В зависимости от целей моделирования в составе МДВО имитируются модели связи и наблюдения, модель ошибок навигации, реализующая современную концепцию навигации, основанную на характеристиках PBN, модель взаимодействия с диспетчером УВД по линии цифровой связи CPDLC (Controller-Pilot Data Link Communications).
В основу представленной в работе модели положены принципы системной динамики и дискретно-событийного моделирования.
описание модели системы Увд
Модель системы УВД рассматривается как сложная система, ее обобщенная структура представлена на рис. 1. Модель имитирует функции диспетчерского управления воздушным движением и формирует текущее состояние воздушных судов, выпол-
No. 5 (53) 2014
Модель метео
Система планирования
Структура ВП
Ветер, ОИВП
Планы полетов
Модель диспетчера
Команды, разрешения
CPDLC, голо(
Запросы
Модель АС УВД
Модель ВС
Модель ДВО
Воздушная обстановка
Рис. 1. Структура модели системы УВД
няющих под контролем УВД полеты по заданным маршрутам в заданных временных интервалах.
В качестве подсистем в модели рассматриваются:
— ВС (собственно ВС, экипаж, бортовая система управления), выполняющие рейсы в соответствии с планами полетов и командами управления и составляющие в совокупности управляемый поток воздушного движения (ВД);
— наземный сегмент УВД (наземные измерительные, радиотехнические и другие средства, диспетчерские пункты управления, которые могут рассматриваться в качестве измерительных и регулирующих подсистем системы управления);
— ограничения ИВП (метео- и административные), ветер.
Функции системы УВД реализованы в модели автоматизированных систем управления воздушным движением (АС УВД). Модель имитирует соответствующие действия диспетчеров районного центра (РЦ), подхода, аэродромной зоны во всей области
моделируемого воздушного пространства (ВП). В АС УВД моделируется функциональное взаимодействие наземной системы УВД и ВС. Модель имитирует действия системы УВД по управлению ВС в целом, обеспечивая контроль за ВС на всех этапах его движения.
Входными данными для модуля диспетчерского управления являются структура ВП, планы полетов, данные наблюдения в своем секторе ответственности, сообщения от диспетчеров соседних секторов, донесения от бортов, метеосводки.
Выходными данными этого модуля являются команды, диспетчерские разрешения, информация о погодных и административных ограничениях.
Модель динамической воздушной обстановки
В состав модели управляемого полета ВС (модель ВС) входят модель бортовой системы самолетовождения и собственно модель движения ВС.
№ 5 (53) 2014
Модель бортовой системы самолетовождения отражает функции бортовых средств: информационных, связи, а также обеспечения выработки команд управления. Бортовая система самолетовождения при выработке команд управления использует плановые данные о маршруте и режимах полета, аэронавигационные данные, метеоданные, данные о параметрах движения и намерениях других ВС, находящихся в некоторой окрестности этого ВС, команды органов УВД.
Модель движения ВС описывает процесс перемещения ВС во времени и пространстве в соответствии с командами управления, поступающими на вход системы управления.
Задача моделирования процессов выполнения управляемых потоков ВД характеризуется большой размерностью, и как следствие, вычислительной сложностью. Так, при исследовании ВД РФ необходимо моделировать процессы регулирования и управления воздушным движением, в котором участвует несколько сотен ВС. Количество диспетчерских пунктов, обеспечивающих полеты в ВП РФ, исчисляется не одним десятком. Это накладывает ограничения на степень подробности описания процессов управления и движения, особенно при ускоренном моделировании. Исходя из этого был принят ряд допущений.
1. Моделируется прямолинейное движение центра масс ВС, модель вращательного движения ЛА необходима только в объеме, влияющем на прямолинейное движение.
2. Система управления реализуется посредством трех каналов:
— курса (управление по направлению);
— вертикального (управление по высоте);
— продольного (управление по скорости).
3. Фазы вылета и посадки описываются упрощенно, наиболее детально моделируется маршрутная часть полета. Это объясняется целью построения модели, направленной на исследование процесса выполнения полетов в верхнем воздушном пространстве.
Тем не менее принятые структурные решения обеспечивают возможность расширения динамической модели с более подробным описанием полетов в приаэродромном пространстве и движения по аэродрому.
Очевидна необходимость разработки различных вариантов моделей ВС, отличающихся уровнем детализации описания динамических процессов и процессов выработки команд управления. Например, при моделировании эксплуатационных процедур разрешения конфликтов в воздухе, прогнозируемых за 3-5 минут и более до минимального сближения, допустимо описывать процессы управления по курсу без учета запаздывания в контуре стабилизации. При моделировании же экстренного уклонения от столкновения (за несколько секунд до встречи) учет такого запаздывания необходим.
Другим примером может служить задача разрешения конфликтных ситуаций путем маневрирования в горизонтальной плоскости, где необходимо описывать движение в канале курса более подробно, а движение в вертикальном и продольном каналах — минимально подробно.
Адаптация моделей требуется также при исследовании различных вариантов и методов бортовых алгоритмов управления ВС (алгоритмов выработки команд в системе управления) в частности:
— точное исполнение плана полета (например, четырехмерных траекторий);
— учет случайных факторов, влияющих на точность исполнения плана полета (внешних, включающих метеоусловия, ограничения ВП, и внутренних, включающих аппаратные, погрешностей, отказы);
— учет возможных мер регулирования УВД, приводящих к изменениям плана полета (задержки, управление интервалами по дистанции и в точках слияния маршрутов);
— исполнение полета по свободным маршрутам (концепция Free Flight).
Также предусматривается возможность моделирования различных методов и алго-
No. 5 (53) 2014
ритмов управления маневрированием (например, при разрешении конфликта).
Такой модульный подход обеспечивает в рамках единой структуры модели настройку на необходимую детализацию описания конкретных схем и алгоритмов управления. В программной реализации это обеспечивается применением объектно ориентированного подхода, а также выбором структуры модели управляемого полета ВС, входных и выходных сигналов отдельных подсистем этой модели.
структура модели воздушного судна
В составе МДВО реализуется модель ВС, структура которой представлена на рис. 2. Модель состоит из следующих компонентов: бортовой системы самолетовождения, системы управления, модели движения ВС и информационной подсистемы.
состав компонентов имитационной модели и их описание
Бортовая система самолетовождения
Эта модель наиболее подвержена изменениям.
В составе входных параметров следующие:
— оценки текущих параметров собственного движения;
— собственный план полета;
— оценки характеристик воздушного движения в заданной окрестности (при наличии таких данных);
— команды диспетчера;
— метеоданные.
Оценки собственного состояния в виде измеренного вектора состояния известны внутри модели ВС в полном объеме; параметры окружающих ВС в виде измеренных векторов состояния поступают от моделей окружающих ВС.
План полета состоит из участков маршрута вылета SID (Standard Instrument Departure), крейсерского полета и маршрута прилета STAR (Standard Instrument Arrival). Все участки условно разделяются на прямолинейные, скругления траектории и смены высоты. План полета представляет собой последовательность поименованных точек воздушного пространства с указанием координат, высоты, скорости полета и радиуса разворота на участке после соответствующей точки.
В качестве команд диспетчера в модель поступают указания по изменению времени и схемы вылета, времени пролета точки, изменению эшелона, схемы захода на посадку, указания на выполнение маневра ожидания.
Метеоданные задают географические зоны опасных метеоявлений и схему ветров. Зоны метеоявлений представляются в виде многоугольных призм с указанием контуров и степени опасности. Ветер представляется в виде эквивалентного ветра на участках воздушных трасс.
Рис. 2. Структура модели ВС
92 ,
№ 5 (53) 2014
Выходные параметры включают:
— назначаемый угол курса;
— назначаемую высоту;
— назначаемую скорость.
Архитектура системы самолетовождения
такова, что методы выработки команд управления при одних и тех же входных и выходных параметрах могут различаться. Модель реализована в объектно ориентированной среде программирования Delphi, что позволяет использовать различные реализации методов для объекта управления. Введен класс ВС, для которого описан ряд свойств и методов, позволяющих настроить каждый экземпляр класса под конкретные методы работы информационных и управляющих систем. В отношении выработки команд управления определен абстрактный класс подсистемы управления, для которого описан ряд наследников, реализующих различные методы управления.
Реализованы следующие варианты самолетовождения:
— выдерживание плановой траектории;
— полет с выдерживанием правил выполнения полетов;
— отработка команд, получаемых от модели диспетчера УВД;
— предотвращение опасных сближений в автоматическом режиме;
— различные комбинации представленных вариантов.
Планы полета
Вариант системы самолетовождения с обработкой команд от модели диспетчера УВД представлен на рис. 3.
Команды диспетчера УВД передаются в систему самолетовождения в виде команд CPDLC. Состав условных команд, составленных в соответствии с [1], приведен в табл. 1.
Команды хранятся и обрабатываются в стеке команд, откуда в расшифрованном виде поступают в блок расчета плановой траектории. Блок расчета плановой траектории предназначен для расчета и поддержания в актуальном состоянии подробных данных (подробного плана полета ВС) по программной траектории полета ВС, выдерживание которой является целью работы системы управления полетом ВС. Далее этот подробный план назовем путевой таблицей (ПТ). Путевая таблица представляет собой набор данных по всей последовательности элементарных участков траектории, описывающих программную полную траекторию от взлета до посадки.
Элементарный участок траектории — это часть полной траектории, отвечающая следующим условиям:
— в горизонтальной проекции участок является частью ортодромии или дугой окружности;
— в вертикальной плоскости на протяжении всего участка идет либо горизонталь-
Стек Блок расчета Путевая таблица Формирование команд
плановой траектории управления
Бортовая система самолетовождения
Рис. 3. Вариант системы самолетовождения при взаимодействии с моделью диспетчера
93
No. 5 (53) 2014
Таблица 1
Ситуация Фразеология Значение
УКАЗАНИЯ, СВЯЗАННЫЕ С ВЫЛЕТОМ CLD (time) РАЗРЕШЕНИЕ ВЫЛЕТА (время)
ВЫДЕРЖИВАНИЕ УКАЗАННЫХ УРОВНЕЙ MAINTAIN (level) ВЫДЕРЖИВАЙТЕ (уровень)
УКАЗАНИЯ ПО ЭШЕЛОНИРОВАНИЮ CROSS (significantpoint), AT (time) СЛЕДУЙТЕ ЧЕРЕЗ (основная точка) В (время)
УКАЗАНИЯ, СВЯЗАННЫЕ С ПОЛЕТОМ ПО ЛИНИИ ПУТИ (СМЕЩЕННОЙ), ПАРАЛЛЕЛЬНОЙ РАЗРЕШЕННОМУ МАРШРУТУ PROCEED OFFSET (distance) RIGHT/LEFT OF (route) СЛЕДУЙТЕ СО СМЕЩЕНИЕМ (расстояние) ВПРАВО/ВЛЕВО ОТ (маршрут)
CANCEL OFFSET ПРЕКРАТИТЕ ПОЛЕТ СО СМЕЩЕНИЕМ
РАЗРЕШЕНИЯ, СВЯЗАННЫЕ С ПОЛЕТОМ В ЗОНЕ ОЖИДАНИЯ CLEARED TO (significant point) [MAINTAIN (or CLIMB or DESCEND TO) (level)] СЛЕДУЙТЕ ДО (основная точка) [ВЫДЕРЖИВАЙТЕ (или НАБИРАЙТЕ ВЫСОТУ, или СНИЖАЙТЕСЬ ДО (уровень)]
ный полет, либо набор высоты, либо снижение.
Расчет подробного плана полета ведется на основе плана полета и является его детализацией. Рассчитанный подробный план полета должен представлять траекторию, выполнение которой возможно с точки зрения ЛТХ ВС. При управлении в соответствии с путевой таблицей блок формирования команд управления формирует команды управления в виде требуемых команд управления: курс, высота (эшелон), скорость.
Система управления
Система управления описывает динамику процессов отработки команд управления и изменения сил и моментов, действующих на ВС.
Описание динамики процессов стабилизации в виде динамического эквивалента можно дать с различной степенью подробности. Для большого круга исследовательских задач допустимо ограничиться простейшей моделью динамики в форме апериодического звена с заданной добротностью (для каждого канала отдельно), при-
ближенно описывающей запаздывание в отработке команд управления.
Динамическое звено для каждого из трех каналов описывается дифференциальным уравнением первого порядка
x = Dx (x * -x),
где x — текущее значение переменной состояния (курс у, высота H, скорость V), x* — команда управления (требуемые у*, H*, V*), DX — добротность в соответствующем канале управления.
Летно-технические характеристики для различных типов ВС описаны в базе данных BADA (Base of Aircraft Data) Eurocontrol [2]. База содержит коэффициенты летных качеств и эксплуатационных процедур для различных типов ВС. Коэффициенты используются для расчетов тяги, силы сопротивления и расхода топлива, а также для указания скоростей номинального крейсерского режима, режимов набора и снижения.
В канале курса накладывается ограничение на максимально возможную угловую скорость изменения курса при координированном развороте
№ 5 (53) 2014
V max
g ■ g(Ym
V
где угп — максимальным угол крена при координированном развороте определен в БД BADA Eurocontrol, V — текущая воздушная скорость ВС. Ограничения в канале курса:
еСЛИ К >Vmax,
если К <-уmax,
Т0 V = Vmax , Т0 V =-Vmax.
В канале высоты ограничения по максимальным скоростям набора и снижения Hmax и Hmin определяются из БД BADA Eurocontol в зависимости от текущей высоты полета H. Ограничения по высоте:
если Н > /"Hmax, Т0 H = Hmax,
если " < Hmin , Т0 H = ""min.
На канал горизонтальной скорости накладываются ограничения в зависимости от высоты полета в соответствии с БД BADA Eurocontrol. Максимальная воздушная скорость ограничивается значением максимальной индикаторной скорости VM,. Минимальная скорость для ВС является функцией от скорости сваливания ВС Vsta¡¡:
VT = О/mln X Vsta„, для Фейсерской фазы
полета коэффициент CVmin = 1,3.
Для пересчета индикаторной скорости в путевую используется модель стандартной атмосферы (Atm), которая описывает зависимость плотности атмосферы и температуры от высоты полета.
Vmin = Atm/mn,, H), Vmax = AtmVO, H).
Ограничения по горизонтальной скорости накладываются на команду управления:
еСЛИ V* > Vmax, ТО V* = ^
если V* < Vm¡„, ТО V* = Vm¡„.
Модель движения ВС
Модель движения ВС рассчитывает текущие траекторные параметры (долготу ф
и широту X, высоту H, путевую V и вертикальную VH скорости, путевой угол f) и другие характеристики эффективности выполнения полета. Движение ВС имитируется путем изменения во времени значений параметров текущего состояния ВС, которые получаются на основе численного интегрирования дифференциальных уравнений, описывающих изменения. Численное интегрирование дифференциальных уравнений производится методом Эйлера первого порядка.
dH / dt = Hi, dV / dt = V d у / dt = \j/,
¿ф / dt = [VB cos f + VwcosVJ / R + H),
dX / dt = [Vb sin f + Vwsin fw) / / [ (R3 + H) cos B),
где ф, X — широта и долгота текущего положения ВС, Vв — текущая воздушная скорость, Vw — величина ветра, fw — курсовой угол ветра, Яз — радиус Земли.
Учет ветра осуществляется в виде эквивалентного ветра на текущем участке для заданной высоты полета.
Интегральная модель погрешностей навигации
Комитет по будущим аэронавигационным системам (FANS) разработал концепцию «навигации, основанную на характеристиках» PBN [3]. Показателем точности навигации, реализуемой посредством применения концепции PBN, принят уровень требуемых навигационных характеристик (RNP — Required Navigation Performance). Параметр RNP k указывает на точность боковой навигации в морских милях, которая должна выдерживаться по крайней мере в течение 95% полетного времени всеми воздушными судами.
В МДВО при моделировании ошибок выдерживания плановой траектории ВС разработана модель случайного процесса оши-
95
No. 5 (53) 2014
бок самолетовождения в соответствии с уровнем RNP к, так как это рекомендуется в [4]. Отклонение ВС, выполняющего полет по RNP к, разыгрывается в соответствии с функцией плотности вероятности f (у):
1 к при Л = —
1 л
Г(у) = ^ ,
ln(0,05)
Параметр X выбирается из условия
а
J f (у )dy = 0,95.
Спектральная плотность случайного процесса распределения ошибок самолетовождения определяется полосой пропускания автопилота. Принимается, что гармоники, соответствующие этой полосе возмущающего воздействия, пропускаются, а более высокочастотные фильтруются. Выходом модели погрешностей навигации является ошибка бокового отклонения ВС от плановой траектории.
Информационная подсистема
На вход информационной системы поступает текущий вектор состояния ВС от модели движения и разыгранная ошибка бокового отклонения от модели погрешности навигации. Зашумленный вектор состояния ВС формируется добавлением разыгранных ошибок боковых отклонений к текущему истинному положению ВС.
Таким образом, бортовая система самолетовождения формирует параметры управ-
ления в соответствии с измеренным состоянием ВС, по которому модель движения формирует параметры истинного положения ВС. Реализованные траектории потока ВС на участке ЛУТЕР — УСТОК приведены на рис. 4, траектории движения показаны сплошными черными линиями.
Заключение
Имитационная модель динамической воздушной обстановки реализована в объектно ориентированной среде Delphi с применением СУБД Oracle. Модель воспроизводит движение воздушных судов в структурированном воздушном пространстве под управлением системы УВД. В созданной модели нашли отражение подходы на основе системной динамики и дискретно-событийного моделирования. Это позволяет представлять конкретные схемы и алгоритмы управления на различном уровне абстракции, описывать систему управления движением ВС и методы выработки команд самолетовождения с разной степенью детализации.
С использованием модели проведены исследования по оценке безопасности полетов при изменении структуры ВП РФ (создание укрупненных РЦ). Проведены исследования по внедрению перспективных бортовых алгоритмов, в том числе обнаружения и предотвращения опасных сближений в воздухе.
корук
Рис. 4. Моделирование полета ВС
№ 5 (53) 2014
Список литературы
1. Правила аэронавигационного обслуживания; Организация воздушного движения; Международная организация гражданской авиации. Doc 4444. Изд. пятнадцатое. ИКАО, 2007.
2. User manual for the base of aircraft data (BADA) revision 3.8, EUROCONTROL, 2010.
3. Руководство по навигации, основанной на характеристиках (PBN); Международная организация гражданской авиации. Doc 9613. Изд. третье. ИКАО, 2009.
4. Единые принципы моделирования риска столкновения в обоснование Руководства по методике планирования воздушного пространства для определения минимумов эшелонирования. Doc 9689. ИКАО, 2009.
5. Емельянов В. В., Ясиновский С. И. Имитационное моделирование систем. М.: Изд. МГТУ им. Н. Э. Баумана, 2009. — 584 с.
References
1. Air traffic management, Procedures for air navigation services. Doc 4444. Fifteenth edition. ICAO, 2007.
2. User manual for the base of aircraft data (BADA) revision 3.8. EUROCONTROL, 2010.
3. Performance-based Navigation (PBN) Manual. Doc 9613. Third Edition. ICAO, 2008.
4. A Unified Framework for Collision Risk Modelling in Support of the Manual on Airspace Planning (Doc 9689). Cir 319. ICAO, 2009.
5. Emelyanov V. V., Yasinovskiy S. I. Imitatcionnoe modelirovanie sistem, MGTU im. N. E. Baumana, 2009. — 584 p.
A. Lyudomir, Technician of GosNIIAS, Moscow, gluek@bk.ru
V. Orlov, Head of Research Sector of GosNIIAS, Moscow, orlov@gosniias.ru
Simulation of dynamical environment in controlled airspace
Paper deals with the solutions of dynamical airspace environment development. The aircraft simulator models generate airspace environment. Each aircraft model simulates an aircraft flight according to its flight plan. The article presents the air traffic control system structure and determines the aircraft simulation model, which includes air navigation, and aircraft movement submodels. Object-oriented implementation allows air navigation submodel to perform several different modes of aircraft control. There are modes of aircraft control that are implemented in simulation model: aircraft flight along its flight path, flight under flight rules, flight with reactions on the air traffic control dispatcher commands, self-separation between aircraft and dangerous meteo areas during flight and various combinations of these modes. Aircraft movement submodel simulates aircraft flight as mass point and imply restrictions and characteristics according to aircraft type. It applies flight dynamics specified by Eurocontrol Base of Aircraft Data. The article describes equations taking into account flight characteristics. The aircraft model simulates navigation deviations based on performance-based navigation concept. The paper presents an example of air navigation under Air Traffic Control with aircraft reactions to Controller — Pilot Data Link Communications commands. Simulation of dynamical environment model implemented in Delphi object-oriented programming environment with Oracle database tools. There are researches, which were conducted using this model: flight safety assessment for new structure of Russian Federation airspace, implementation of airborne ADS-B applications, including conflict detection and resolution application. Keywords: simulation, Air traffic control, PBN, CPDLC.
4.....97
