CC BY
6Электронный журнал «Труды МАИ». Выпуск № 66
www.mai.ru/science/trudy/
УДК 550.312; 528.27; 531.76; 656.7:658.012.011.56; 656.7:004; 681.5.015.44
О возможности построения бесплатформенного управляющего навигационно-гравиметрического комплекса беспилотного
летательного аппарата
Афонин А. А.,* Сулаков А. С.**, Ямашев Г. Г., Михайлин Д. А., Мирзоян Л. А.***,
Курмаков Д. В.
Московский авиационный институт (национальный исследовательский университет), МАИ, Волоколамское шоссе, 4, Москва, А-80, ГСП-3, 125993, Россия
*е-таИ: Al_aa@mail.ru **е-таИ: kaf305-mai@mail.ru ***е-таН: lola-mai@mail.ru
Аннотация
В статье представлен новый подход к построению интегрированного информационно-измерительного управляющего комплекса беспилотного летательного аппарата. Рассматривается вариант комплекса, в состав которого входят система планирования полета и система автоматического управления беспилотным летательным аппаратом, а также бесплатформенный гравиинерциальный навигационный комплекс -основной источник ориентационной, навигационной и гравиметрической информации, как для работы системы управления, так и для реализации целевой функции полета -проведения автоматизированной векторной аэрогравиметрической съемки. Представлены функциональные и структурные схемы, математические модели и алгоритмы работы основных узлов комплекса.
Ключевые слова:
бесплатформенный гравиинерциальный навигационный комплекс, функциональный алгоритм, гироскоп, акселерометр, спутниковая навигационная система, магнитометрическая курсовая система, аномалия, векторная гравиметрия, система планирования полета, система автоматического управления, нейросетевая система, оптимальная оценка, закон управления, автомат тяги
Введение
Сегодня по-прежнему важной и значимой остается проблема проведения высокоточных информативных оперативных детальных гравиметрических измерений в целях геофизического мониторинга природно-техногенной сферы, решения ряда важных задач геодезии и физики Земли, а также реализации возможностей высокоточной автономной навигации с использованием информации о геофизических аномалиях. Кроме того, все большую остроту приобретают вопросы гравиразведки залежей полезных ископаемых, включая природные энергоносители, в особенности углеводородные соединения. В настоящее время наиболее оперативным и производительным, универсальным и перспективным способом, позволяющим выполнять гравиметрические измерения в трехмерном пространстве практически над всеми территориями и акваториями земного шара, включая удаленные и труднодоступные, является аэрогравиметрическая съемка. Данная работа посвящена разработке новых подходов к построению автоматизированного комплекса технических средств, обеспечивающего качественные и количественные преимущества при ее проведении как в части собственно мобильных гравиметрических измерений, так и в части планирования, информационного обеспечения и управления движением ЛА над исследуемым регионом.
Современная технология аэрогравиметрии остается практически неизменной на
протяжении уже нескольких десятков лет и заключается в исследовании территорий и
акваторий, покрываемых ЛА сетью параллельных и секущих измерительных галсов [1].
Собственно измерения аномалий ускорения силы тяжести (УСТ) выполняются
мобильными гравиметрическими комплексами, расположенными на борту ЛА (самолетов,
вертолетов, дирижаблей и др.), включающих высокоточный гравиметрический датчик,
установленный на виброизолированной прецизионной гиростабилизированной платформе
и интегрированный с корректирующими подсистемами. Достоинствами современной
мобильной гравиметрической съемки являются достаточно высокий уровень аппаратной,
программной и технологической отработанности, а также точность (~ 0,1 - 1 мГал),
удовлетворительная для большинства поставленных перед аэроморской гравиметрией
задач, в том числе для поиска и разведки месторождений полезных ископаемых. Главным
недостатком современной технологии аэрогравиметрии является необходимость
установки массивного однокомпонентного гравиметрического датчика на
гиростабилизированную платформу для его выставки по местной вертикали, что в
большой степени увеличивает массогабаритные показатели комплекса (до более чем 1 м
и 100 кг), энергопотребление (до сотен Вт), стоимость (до сотен тыс. долл. США), это в
2
подавляющем большинстве случаев предопределяет использование в качестве носителей гравиметрического оборудования самолетов и вертолетов среднего класса, небольшой маневренности, имеющих ограниченные возможности для оперативных и детальных гравиметрических исследований. При этом, как правило, стоимость эксплуатации воздушного судна оказывается чрезвычайно высокой, а детальность и оперативность измерений - все еще недостаточными в условиях всевозрастающих требований современной инновационной экономики. Кроме того, сохраняется потенциальная опасность проведения таких работ в воздухе.
Таким образом, становятся очевидными актуальность и перспективность разработки комплекса технических средств, способного обеспечить решение задачи аэрогравиметрии на качественно новом техническом и экономическом уровне с привлечением современных достижений в области бесплатформенных и спутниковых навигационных технологий, средств электроники, информатики и вычислительной техники на базе современных малоразмерных, маневренных, недорогих, автоматически управляемых беспилотных летательных аппаратов (БЛА).
Новая технология аэрогравиметрии
Для разрешения круга обрисованных проблем современной аэрогравиметрии необходимо устранить ее главный недостаток - необходимость применения дорогого и массивного гиростабилизатора в конструкции гравиметрического комплекса для стабилизации высокоточного, но громоздкого однокомпонентного гравиметрического датчика. Это, в свою очередь, позволит существенно снизить массу, габариты, энергопотребление и стоимость комплекса, позволит перейти к его применению на борту автоматизированных аэроморских маневренных малоразмерных носителей, включая БЛА. Кроме того, в особенности в рамках малоразмерного автоматически управляемого борта, актуальным является вопрос аппаратной и алгоритмической интеграции систем планирования полета (СПП) и автоматического управления (САУ) аппарата-носителя с перспективным измерителем параметров ориентации, навигации и гравиметрии -бесплатформенным гравиинерциальным навигационным комплексом (БГНК) [2,3]. Итак, работа направлена на исследование возможности создания интегрированного высокоточного малогабаритного управляющего навигационно-гравиметрического комплекса (УНК), отличающегося лучшими массогабаритными, стоимостными характеристиками, большей информативностью, меньшим энергопотреблением по сравнению с гиростабилизированными комплексами при соизмеримой погрешности измерения, а также обладающим не уступающими или лучшими характеристиками
3
надежности и качества управления подвижным объектом-носителем по сравнению с существующими аналогами.
Возможность применения нового комплекса технических средств на борту перспективного класса носителей - малоразмерных БЛА - предопределяет возможность перехода к новой технологии бесплатформенной автоматизированной аэрогравиметрии. В ее основе лежат:
• бесплатформенный принцип построения БГНК, включающего следующие измерительные подсистемы: инерциальный измерительный блок (ИИБ) в составе прецизионных датчиков угловой скорости (гироскопов) и кажущегося ускорения (акселерометров); резервные микромеханический ИИБ и магнитометрическую курсовую систему, свойства и характеристики которых непрерывно уточняются с использованием алгоритмов оптимального оценивания с привлечением, например, нейросетевых методов в рамках совместной эксплуатации с основным ИИБ; навигационную аппаратуру потребителя (НАП) спутниковой навигационной системы (СНС), работающей в дифференциальном или фазовом режиме; систему воздушных сигналов (СВС); лазерный или радиовысотомер; доплеровский измеритель относительной скорости и др.;
• уникальная методика модернизации используемых в современной практике алгоритмов определения параметров ориентации, навигации и гравиметрии, основанная на рациональной модификации структуры оптимальной калмановской фильтрации с расширением вектора состояния системы;
• новая нейросетевая СПП, реализующая стратегический уровень оптимального планирования текущих целей полета и тактический уровень формирования оптимальной траектории полета БЛА с учетом особенностей зависимостей оцениваемости горизонтальных и вертикальной проекций аномалии УСТ от линейных и угловых параметров движения БЛА;
• специализированные алгоритмы управления движением аппаратов-носителей с высокоточной САУ, позволяющие реализовать оптимальную траекторию движения, повышающую общую точность работы комплекса;
• использование в САУ статистических алгоритмов оптимального оценивания неизмеряемых на борту параметров объекта, применение оригинальных законов управления на основании оценок неизмеряемых параметров;
• аппаратная реализация алгоритмов комплекса (ориентации, навигации, гравиметрии, оценивания параметров и коррекции), а также нейросетевой САУ, на
базе современных высокопроизводительных микроконтроллеров и программируемых логических интегральных схем (ПЛИС); • итоговая минимизация массогабаритных, энергопотребляющих, стоимостных характеристик комплекса, повышенная надежность, адаптация для работы в составе малоразмерного автоматического объекта.
В свою очередь, недорогие малоразмерные динамичные БЛА, оснащенные УНК, будут отличаться большей мобильностью и маневренностью, значительно меньшими затратами на эксплуатацию по сравнению с традиционными носителями гравиметрического оборудования, что обеспечит возможность аэрогравиметрических измерений в труднодоступных регионах, удаленных от оборудованных аэродромов. Важной видится возможность повышения производительности съемки посредством вывода группы автоматизированных БЛА в исследуемый регион. Данная технология обеспечит также улучшение детальности и оперативности пространственных и площадных гравиметрических съемок, геофизического мониторинга территорий и акваторий, решения задач информационного обеспечения автономной навигации по геофизическим полям и др. Существенной является принципиальная возможность проведения векторных гравиметрических измерений на борту динамичного объекта (определение модуля силы тяжести и уклонений отвесной линии), тогда как традиционные технологии позволяют определять лишь модуль силы тяжести. Существенно увеличивает надежность работы комплекса применение резервных МКС и интегрированного в бортовой вычислительный модуль недорогого и компактного микромеханического ИИБ в совокупности с возможностью текущей идентификации их параметров и компенсации погрешностей, а также адаптацией алгоритмов комплекса на использование МКС и микромеханического ИИБ по мере необходимости. Вместе с тем отсутствие экипажа и обслуживающего персонала на борту в совокупности с автоматизацией работ кардинально снижает общие затраты на проведение гравиметрической съемки и исключает риски для жизни и здоровья людей.
В настоящее время законченные зарубежные и отечественные аналоги УНК отсутствуют. Проведенные исследования показывают, что разрабатываемый комплекс по ряду параметров должен превзойти отечественный и зарубежный достижения в области аэрогравиметрии. Ожидаемые характеристики УНК: погрешность измерения на подвижном основании - порядка десятых долей мГал, при частоте обновления выходного сигнала в единицы Гц; масса - до 15 кг; объем - до 20 л; энергопотребление - до 100 Вт; прогнозируемая стоимость - до 100 тыс. долл. США.
Таким образом, разработка УНК является перспективным направлением в рамках развития новой технологии аэрогравиметрии. Ввиду ожидаемых технологического и экономического эффектов от внедрения УНК в соответствующие отрасли в совокупности с огромными (по данным международной рейтинговой компании Standard & Poor's - на уровне сотен млн. долл. США в год) общемировыми вложениями в геофизические исследования и гравиметрическое оборудование, в успешной реализации проекта создания УНК могли бы быть заинтересованы компании геофизического, геодезического, геологоразведочного профиля, в том числе занимающиеся поиском и разведкой полезных ископаемых, в особенности углеводородных соединений, а также организации, занимающиеся вопросами высокоточной навигации и геофизического, экологического и природоохранного мониторинга. Потенциальных инвесторов должен привлечь и тот факт, что проект сопряжен с низкими техническими рисками, так как разработка в целом строится на базе уже существующего и стремительно развивающегося аппаратного обеспечения (прецизионные навигационные акселерометры и датчики угловой скорости, магнитометры, НАП СНС, ПЛИС, микроконтроллеры и др. измерительные и электронные модули) и преимущественно ориентирована на глубокую модернизацию современных существующих функциональных алгоритмов гравиинерциальных комплексов, способов управления работой и обработки данных измерительных подсистем, а также алгоритмов СИП и САУ БЛА с их аппаратной реализацией на базе микроконтроллеров и ПЛИС.
Вариант УНК ограниченного состава
В зависимости от целей применения и условий эксплуатации УНК БЛА может существенно менять свой состав, однако общая структура и базовый набор основных узлов остаются неизменными. На рис. 1 в качестве примера представлена блок-схема рационального варианта УНК с ограниченным набором измерительных подсистем, к которым относятся: ИИБ, содержащий блок акселерометров (БА) и блок гироскопов (БГ), НАП СНС, МКС и СВС, включающий датчики воздушной скорости, углов атаки и скольжения. В рамках вычислительного модуля УНК реализуется программно-алгоритмическое обеспечение комплекса, в том числе алгоритмы БГНК и СПП, законы управления САУ БЛА. Результаты выполненных измерений и вычислений фиксируются внешним информационным накопителем.
В данном разделе основное внимание сосредоточено на демонстрации особенностей формирования алгоритмов определения параметров ориентации, навигации и гравиметрии, оптимального оценивания и коррекции. Основные принципы, положенные в основу построения СПП и САУ БЛА, изложены в следующем разделе.
6
Одной из важнейших является задача начальной выставки, решаемая комплексом. Выставка производится в два этапа. На первом относительно грубом этапе начальной выставки (см. рис. 1, пунктирные стрелки), выполняемой во время стоянки на неподвижном основании, можно полагать и0 = (0,0,0)т. Я может быть определен как из внешнего источника данных о местоположении БЛА, так и посредством собственных измерений НАП СНС, предпочтительно длительно выполняемых с текущей или последующей статистической обработкой (осреднение, низкочастотная фильтрация, обработка статистическим алгоритмом, например, методом наименьших квадратов). Определение ориентации относительно плоскости горизонта (тангаж, крен) осуществляется посредством алгоритма акселерометрического построителя вертикали (АПВ) [4,5] с использованием аналогичной статистической обработки показаний БА. Начальная выставка в азимуте (курс) также проводится по статистически и алгоритмически обработанным показаниям МКС [4,5] или БГ, во время стоянки измеряющего только и. На основе Я и у0, -Э0, у0, определяются А0/Е0 и Л0. Начальные условия остальных параметров функционального алгоритма можно принять нулевыми.
Рисунок 1 - Блок-схема УНК ограниченного состава На рисунке: БПО - блок предварительной обработки; БОПО - блок определения параметров ориентации; БО - блок определения; БСО - блок статистической обработки; ВИ - внешняя информация; НУиП - блок вычисления начальных условий и параметров; МНУСТ - модель нормального УСТ; БРФА - блок реализации функциональных
алгоритмов БГНК; БИ - блок измерений; ОФК - блок реализации алгоритма оптимального фильтра Калмана; БВДП - блок вычисления дополнительных параметров; БВМС - блок вычисления поправки на магнитное склонение Думс; Ao/E - матрица ориентации базовой правой Земной экваториальной системы координат (ЭСК) (начало в центре Земли, орты ^ расположены в плоскости экватора, где £ - на линии пересечения плоскости экватора с Гринвичским меридианом, а орт ^ направлен по вектору и угловой скорости вращения Земли) относительно связанной с объектом системы координат (ССК); индексом 0 обозначены начальные условия соответствующих параметров; БП - блок пересчета показаний акселерометров в ЭСК; у, у - курс, тангаж, крен; ю - вектор абсолютной угловой скорости и п - вектор кажущегося ускорения объекта, измеренные БГ и БА; индексом o обозначены векторы, представленные в проекциях на оси ССК (без данного индекса - на оси ЭСК); В - вектор индукции геомагнитного поля; Я -геоцентрический радиус-вектор местоположения объекта; и - вектор относительной скорости объекта; индексы С, м, а отвечают параметрам, вычисляемым СНС, МКС и АПВ; символ 5 обозначает погрешность соответствующей величины, Л - ее оценку; юШ -частота М. Шулера; ^гц, FR - проекции напряженности нормального поля тяготения Земли на оси сопровождающей геоцентрической системы координат (СК); Г - матрица динамики; О - матрица шумов системы; Н - матрица измерений; Р - матрица дисперсий шумов системы; Ь - матрица дисперсий шумов измерений; Р - апостериорная матрица ковариации; К - матрица коэффициентов усиления ОФК; г - вектор измерений; Л=(Хо,Х1,Х2,Хз)т - вектор параметров Родрига-Гамильтона, характеризующий поворот ССК относительно инерциальной СК; А - вектор аномалии УСТ; Аи - А, заранее известный с ограниченной степенью точности (при отсутствии исходных данных Аи=0); ц - вектор коэффициентов затухания корреляционных функций, выбранных для описания соответствующих проекций случайных составляющих аномалии УСТ, о - вектор СКО
проекций случайных составляющих аномалии УСТ; t - время; р = (р0 ■■■ рт)т - вектор
псевдодальностей, р = (р0 ■■■ рт)т - вектор псевдоскоростей (т+1) используемых
навигационных спутников; а, в - углы атаки и скольжения; V - вектор воздушной скорости; Да^ Дрw и Да, Др - приращения углов и систематических составляющих скоростей вертикального и бокового ветра соответственно.
Второй высокоточный этап начальной выставки, также проводимый на стоянке, и основной этап работы комплекса как на стоянке, так и в полете отличаются способом формирования вектора измерения. В зависимости от используемой слабосвязанной
(сильносвязанной) схемы коррекции БИ формирует разности соответствующих координат (псевдодальностей) и проекций относительной скорости (псевдоскоростей), вычисленных БРФА и СНС, а также угла курса, вычисленного БВДП по показаниям БРФА и вычисленного БИ по показаниям МКС или БГ. Применение слабосвязанной схемы целесообразно на стоянке и при движении по маршруту, в случае маневренного полета и в сложных условиях эксплуатации БЛА предпочтительно применение сильносвязанной схемы. Помимо этого только на втором этапе начальной выставки формируются разности угла тангажа (аналогично для крена), вычисленного БВДП по показаниям БРФА и вычисленного БИ в соответствии с алгоритмом АПВ по показаниям БА.
Для решения функционального алгоритма на втором этапе начальной выставки необходим ввод начальных условий из НУиП и параметров модели нормального УСТ из МНУСТ. При переходе к основному рабочему этапу начальные условия непосредственно выбираются в БРФА из результатов второго этапа. Решение функционального алгоритма и алгоритма ОФК происходит совместно, параллельно: по текущим входным данным и показаниям БРФА формируется г, на каждом шаге алгоритма обновляются используемые в ОФК матрицы Н, Ь, Р, Г, О и вычисляется К, используемая в произведении с г для коррекции БГНК по замкнуто-разомкнутой (смешанной) схеме. Итоговая информация
БРФА - вычисленные Я, и, Л, А0/е, А, ю0, п0, Даж, Д(3Ж, Да, Др, а также оценки
систематических погрешностей комплекса, использованные при их компенсации, применяются в БВДП для вычисления набора дополнительных параметров (высота, географическая широта, путевой угол, курс, тангаж, крен, проекции истинного УСТ и др.), а также в СПП и САУ для управления объектом. Данные измерений и вычислений сохраняются в памяти для последующего анализа результатов работы.
Собственно функциональный алгоритм БГНК, представленный в ЭСК, имеет вид:
(
л=0,5 м«,„ л; <0 О =«о -§ «оп -
ЬГ 0!
©: 0а
А =
АО/Е
С сов(и^) 8т(и^) 0^ - ът(М) соб^ ) 0 0 0 1
(Эи Ькс _
ух у уг
0й 0й 5Г
у гх гу г
«о +АОж К 0-2> г;
[(-1) Е+2 [ ккТ ]-2^0 Lx
и=-2ихи+А0Жп0 +А-К3-5>г; Я=и-К 0-8 г;
п О =ПО -§ ПОп -
Ьк" ©
0
0 Ь к" 0 ©" 0" Ьк
у гх гу г
Ьп,„=К9-11>г; =К12-1\; ЬАп=К15-17>г;
(30-32) п т^ (33-35
Оп "' Оп
Г к5А к ЬА коЬА ^
кЪА к ЬА (»-»0 )
кЬА У кЬА ГоЬА у
Д Ь А +К 18-20) 'г ч
кЬА=К21-29/г; Ьки=К30-32/г; Ькп =К33-35г; 0и=К(30-41>г; 0п =К<42-44>г; Ьум=К(45)г;
слабосвязанная схема
Ь =К 40-48 г;
Ь ис11=К<49-51 г;
^-Яс +Ь я сп^
г=
и-ис+Ь исп
М +Ь V;
Мп У
или
сильносвязанная схема
" 40-(40+т))
Ь Рсп =К
Ь Р сп =К
=
'г;
|(47+м )-(47+2хт )
г;
Ь т=к< (48+2хт)> г;
г=
Р-Рс +Ь Рсп
р-Р с +Ь р сп У-У м +Ь V Мп
этап начальной выставки
(
г=
Л
,У-У А
(1)
О
и
п
где индексы л, ^ и х, У, 2 соответствуют проекциям на оси ЭсК и ссК соответственно;
где индексы л, С и X, У, 2
Г 0 -а3
= а3 0
1-а2 а,
а
Л
0 у
М =
^ 0 -а -а2 -а3 ^
0 а -а
0 а,-
¿2 -а1 0 У
■ матрицы, составленные из
элементов некоторого вектора а=(а1,а2,аз)т; к™'" - масштабные коэффициенты гироскопов или акселерометров с измерительной осью г, кю,и - составленные из них вектор-столбцы; 0®'п -
малые углы отклонений измерительных осей гироскопов или акселерометров от идеальной оси г вдоль оси ' ссК, 0Ю,И - составленные из них вектор-столбцы; индекс <г' показывает, что в матрице K сохранены с г-й по у'-ю строки; и - модуль u; E - единичная матрица
соответствующей размерности; gТн - нормальный вектор УсТ; кЬА - линейные коэффициенты
модели погрешности аномалии в проекции на ось г, зависящие от координаты или параметра у, к5,4 - составленный из них вектор-столбец; Dц - диагональная матрица, составленная из элементов ц; индекс сл обозначает случайные составляющие величин, а индекс п отвечает их систематическим постоянным составляющим; 8к - геоцентрический радиус-вектор местоположения к-го спутника и \к - его относительная скорость, вещаемые снс;
рк =у/- Я )2 + (Б* - Я )2 + - Я^ )2 - псевдодальность к-ого спутника
• (Sk - R-) • (V-k - U-) + (Sf- R) • (Vf - Ц) + (Sj - R) • (Vk - U,) = —---- • 11 ' ' 1 1-1-— - его псевдоскорость,
V(Sf - R )2 + (Sf - R )2 + (Si - R )2
вычисляемые БГНК; 5т - погрешность часов НАП СНС.
Курс, вычисляемый по показаниям МКС, описывается соотношением
В sin у + B cosy
= -arctg (---) + Думс.
В cos & - (В cos Y - В sin Y) sin &
При этом могут использоваться как итоговые значения тангажа и крена вычисляемые в БВДП, так и значения, полученные в соответствии с алгоритмом АПВ:
&а = arctg( пх/ jny2 + пг2 ); Ya = -arctg( njny ),
где также могут использоваться проекции как по из БА, так и п0 из БРФА.
В функциональном алгоритме БГНК gTH описан соотношениями [6], являющимися наиболее рациональными для реализации алгоритма в ЭСК:
gTH = ^Лц + Fr 1r - u X (u X R); 1r = R(R • R)-0,5 = [l-, l,, l, ]т;
Егц = (Zx 1r)(1 -12)-0,5; Nгц = 1R XЕгц; R = (R• R)0,5;
Fm = g« (q-e2)(a/R)41,(1 -1,2)05 [ 1 + e2(0,5e2 -30q/14)(q-e2)-1 ]x x{1 + [30q /14 -1,5e2 +12 (3,5e2 - 5q)] (q - e2 )-1 (ea / R)2 J;
Fr = -g (a / R)2 {1 - 0,5e2 - 0,125e4 + q (1,5 - 15e2 / 28) + [0,5e2 - 0,25e4 + + q (-0,5 + 15e2 /14)-12 (1,5e2 - 0,75e4 + q (-1,5 + 45e2 /14))](a / R)2 +
x
+ [0,375e2 - 15q /28 +l2 (0,625e2 - 25q / 28) -
| (1 - ¡1) (4,375е2 - 6,25д)] е2 (а / Я)4 },
где Егц (Восток), Nгц (Север), 1д (геоцентрическая вертикаль) - орты сопровождающей геоцентрической СК; ^ - величина нормального УСТ на экваторе; q - отношение центробежной силы, возникающей вследствие вращения Земли, к силе тяжести на экваторе; е = (а2 — Ь2)0'5а- и а, Ь - первый эксцентриситет и величины большой и малой полуосей референц-эллипсоида. При необходимости функциональный алгоритм позволяет использование других моделей нормального поля тяготения, предполагающих зависимость от экваториальных, геоцентрических или географических координат. При этом может потребоваться пересчет определяемых комплексом экваториальных координат в геоцентрические и далее в географические на каждом такте вычислений [7,8].
Параметры ориентации объекта (у, 0, у) в БВДП вычисляются согласно унифицированному алгоритму [4]:
Z = (x X r)|x X r|—1; v = -arctg(Nrr • Z1(Err • ZJ—j 3 = arctg(x • r(1 - (x • r)2)-0,5); y = arctg(y • Zi(z • ZJ1), J
где x, y, z - орты ССК (представляют собой 1-й, 2-й и соответственно 3-й столбцы Ao/E); Err, Nrr, r - орты восточного, северного и вертикального направлений географической СК. Эти орты вычисляются по формулам:
Err = Eгц; Nrr = rXErr; r = 1r +(1 + 02/4)-1 (0х(1л + 0x 1R/2)),
где 0 = 2Err tg ((фгц — ф) /2) - соответствующий вектор конечного поворота; Фщ = arcsin(^) - rеоцентрическая широта; ф - rеоrрафическая широта, которая может
быть определена в БВД11 одним из известных способов - посредством аналитическоrо [7] или численноrо [8] решения алгебраических уравнений, либо дополнительно решая в БРФА дифференциальное уравнение
ф = ((RR, + RR) tgф — R^R2 + R) (1 — e2 sin2 ф) cos2 ф Дa e2cos^(1 — e2 sin2 ф) х x^R¡ + R2 — (R2 + R2) (1 — e2 sin2 ф) + ((r2 + R2) tgф — R^R + R ) e2 sin ф cos3 ф).
В тоже время на ограниченных высотах можно пользоваться упрощенной формулой
r = —g TH gTH ,
имеющей небольшую методическую ошибку в северном направлении. Ошибка практически равна нулю вблизи полюсов, экватора и поверхности референц-эллипсоида, достигает максимума в средних широтах и практически линейно зависит от высоты. При этом, например, на высоте ±10 км и широте 45° ошибка ориентации r составит всего лишь ±1,7 уга. с, что удовлетворяет требованиям векторной мобильной травиметрии.
Таким образом, описана общая структура, состав и основные алгоритмы части УНК, обеспечивающей функции БГНК. Основные подходы к построению функциональных алгоритмов оценивания неизмеряемых на борту параметров, реализуемых также в БРФА и используемых в САУ БЛА, приведены в следующем разделе.
Основные принципы построения СПП и САУ
Управление полетом БЛА (см рис. 2) можно представить как иерархический процесс, который условно можно разбить на три уровня:
• стратегический уровень (уровень планирования текущей цели);
• тактический уровень (уровень формирования траектории полета);
• исполнительный уровень.
Рисунок 2 - Трехуровневая концепция управления полетом БЛА
стратегический уровень, или уровень принятия решения о выборе поведения, состоит в оперативном планировании в реальном времени определения текущей цели. Задачей тактического уровня является формирование траектории полета по заданным координатам цели. Эти уровни управления реализуются сПП. Третий - исполнительный уровень, это уровень сАУ БЛА, в задачи которого входит поддержание или отработка с помощью управляющих органов задающих воздействий по таким сигналам, как желаемое направление, скорость, высота полета к цели.
В работах [9,10] предложена нейросетевая система планирования траектории, которая реализует стратегический уровень управления полетом БЛА. Такая нейросетевая система состоит из иерархических нейросетевых структур на базе управляющего нейроэлемента и формируется эволюционным методом с использованием моделирующей среды. Реализация нейросетевых структур планирования в конвейерном режиме с помощью малогабаритных ПЛИс позволяет значительно снизить стоимостные и массогабаритные характеристики системы планирования траектории в целом.
Нейросетевая система планирования траектории взаимодействует с другими управляющими системами БЛА следующим образом. Получая информацию БГНК о текущем положении БЛА, а также информацию о появлении новых объектов для наблюдения, нейросетевая система формирует сигнал в виде координат текущей цели
полета, на основании которого в дальнейшем формируется задающий сигнал для регуляторов, обеспечивающий полет к выбранной цели.
На рис. 3 представлена структурная схема контура управления полетом БЛА, содержащая:
• нейросетевую систему планирования траектории на ПЛИС (стратегический уровень управления);
• преобразователь координат цели в задающие сигналы - текущие высоту и угол рыскания, которые обеспечат полет к выбранной нейросетевой системой цели (тактический уровень управления);
• регуляторы рыскания, крена, тангажа, высоты полета и скорости (исполнительный уровень управления).
Таким образом, встроенная нейросетевая система планирования траектории формирует сигналы, передающиеся через задающие воздействия ко всем регуляторам САУ БЛА, т.е. находится в основном контуре управления и замыкает его.
Рисунок 3 - Структурная схема контура управления полетом БЛА На рисунке: X, У, Н и Х0, У0, Н0 - соответственно текущее координаты (например, широта, долгота и высота) БЛА и попадающих в его поле видимости объектов - индекс 0; индексы з и ц отвечают задающим и целевым значениям или сигналам; 5н, 5э, 5В, 5др - управляющие
сигналы по отклонению руля направления, элеронов, руля высоты, положению дроссельной заслонки двигателя соответственно; ур - угол рыскания.
На исполнительном уровне, для обеспечения автоматического управления малоразмерным БЛА при любых погодных условиях и различных комплексных возмущающих воздействиях, а также при ограничениях скоростной и позиционной характеристик исполнительного привода, целесообразно применять новые алгоритмы управления движением БЛА. Алгоритмы используют оптимальные калмановские оценки неизмеряемых координат [11] - скоростей и углов вертикального и бокового ветра и обеспечивают устойчивость и требуемое качество управления движением. Так как в процессе полета меняются характеристики объекта, то необходимо гарантировать адаптивность алгоритмов управления к изменяющимся условиям полета. В частности, при управлении скоростью БЛА предлагается использовать алгоритм, содержащий нейронную сеть, которая на основании текущих полетных данных в режиме реального времени осуществляет аппроксимацию сложной зависимости потребной мощности двигателя от текущей высоты, угла атаки и заданного угла наклона траектории. Для простоты изложения принципов делаем допущение, что в САУ используются быстродействующие приводы рулей высоты, направления, элеронов и их динамикой можно пренебречь.
Движение БЛА в продольной плоскости описано следующей системой уравнений:
_ТТ _g
ДН = ^-°(ДЗ- Да); Да = Дюz - 7ДДа + Даw) - Y/Д5в; 180
Дю z = Мга(Да + Да w) + М'Юz Дюг + Щ Да + M5/ Д5в;
л а л л • 1 А 180 . 180 : _
ДЗ = Дю z; Да w = --Да w + Да + Па; Д а = 0; T nT V0 nT V0
z =
ГДН 1
0 01 Да
0 0 Дюг
0 0 ДЗ
0 0, Да w
V Да ,
+
где ДН=Н-Н0; ДНз=Нз-Н0; Дх - приращение некоторой величины х; индекс 0 - здесь и ниже отвечает установившимся параметрам движения БЛА; £ - вектор шумов измерений; У"' - производная приведенной к произведению массы БЛА на скорость установившегося
—5„
полета аэродинамической подъемной силы крыла по а; У а - производная приведенной к произведению массы БЛА на скорость установившегося полета аэродинамической
подъемной силы по 5в; М^ - производная приведенного к моменту инерции относительно связанной оси Oz момента 0 по a; M^z - приведенный к моменту инерции относительно связанной оси Oz демпфирующий момент Mza - производная приведенного к моменту инерции относительно связанной оси Oz момента 0 по а ; М®в -производная приведенного к моменту инерции относительно связанной оси Oz момента 0 по 5в; Т - постоянная времени формирующего фильтра; - белый шум.
Закон управления БЛА в продольной плоскости выглядит следующим образом: 8в = ktâ (АН -АН) + к^н АН + каАЗ + к^ Аш. + kak^ Ad w + kakAa/i a,
где ka, kH, кй, к , k^, kAa - коэффициенты закона управления. В отличие от классического закона управления высотой ôB = kH(AH - АНз) + каА& + kWzAœz [12], где kH -
его коэффициент, в новом законе управления появляются составляющие от приращения вертикальной скорости, оценок приращений угла и скорости вертикального ветра. Движение БЛА в поперечной плоскости описывается системой уравнений: лК
AZ
180
(А^р -АР);
Ар = Zр (АР + А^ ) + sin а0 Ашх + cos а0 Аш + ^ cos Ау + Zô- Абн ;
- Ко —0
АШх = Мр (АР +Ар„ ) + Мхш Ашх + МШ Аш» + М»; Абэ + М^ Абн ;
АШ y = M р (АР +АР„ ) + МШ Аш y + M*» Ашу + Мъ; Абэ + Мъ; Абн,
Ау = Ашх - tg§oАшу ; АФр =
1
т
180
лТК
—1— Ашу ; cos у
180
0
лТК
Г1 000000 0Л 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0
Лр ; А р =
(AZ Л
АР
А®х
Аш» Ау +
А^р
АР w
А ^ Р у
:
z
где ДZ=Z-Z0; ДZз=Zз-Z0; Z - горизонтальная координата, отсчитываемая перпендикулярно
от линии заданного пути; ^р - белый шум; ZР - производная приведенной к произведению массы БЛА на скорость установившегося полета боковой аэродинамической силы по в; Z5н - производная приведенной к произведению массы
БЛА на скорость установившегося полета боковой аэродинамической силы по 5н; М *, М в - производные приведенных к соответствующему моменту инерции моментов у и ур по в;
-га* -гау
М* , Му - приведенные к соответствующему моменту инерции демпфирующие
-га^ -гах
моменты у и ур; Мх , Му - приведенные к соответствующему моменту инерции
-5Э -5Э -5Н -5Н
перекрестные моменты у и ур; М* , Му , М* , Му - производные приведенных к соответствующему моменту инерции моментов у и ур по 5э и 5н.
Управление БЛА в поперечной плоскости осуществляется с помощью элеронов и руля направления, при этом канал элеронов используется для управления траекторией движения БЛА, а канал руля направления - для устранения скольжения. Законы управления БЛА в поперечной плоскости выглядят следующим образом: 5э = к1 к^ (^ - ^) + к^Ау + куАу + ках Агах + к^ А|3„ + £у£дрА р;
5н = к£ дУ + кга уАга у ,
где ку, kZ, ку, к^, кга, кга, к^, кдр - коэффициенты законов управления. В отличие от классического закона управления элеронами 5э = kYkZ(АZ -АZз) + к^Ау + куАу + кгаАга* [13], где ^ - его коэффициент, в новом законе появляются составляющие оценок приращений угла и скорости бокового ветра.
Итак, ожидается и подтверждено проведенными оценками, что предложенные законы управления рулем высоты и элеронами по сравнению с классическими законами обеспечат лучшее качество управления БЛА при наличии возмущающих факторов.
Скорость БЛА управляется автоматом тяги (рис. 4), обеспечивающим формирование потребной тяги двигателя. На его выходе - сигнал Р, определяющий положение дроссельной заслонки двигателя, а также величина оборотов двигателя пдв. Входными сигналами являются Уз и V, Н, а и заданный угол наклона траектории 9з.
Рисунок 4 - Структурная схема автомата тяги На рисунке: ПР - привод дроссельной заслонки; СУ - силовая установка; kдV -коэффициент влияния невязки по скорости; ^ - коэффициент обратной связи по оборотам двигателя; - коэффициент масштабирования величины потребной мощности в величину отклонения дроссельной заслонки; #потр - величина потребной мощности двигателя; Uдр - сигнал управления дроссельной заслонкой.
Закон управления автомата тяги записывается следующим образом [13]: тт _ к N - к п
т др кыр1у потр кпп дв ,
где величина потребной мощности представлена как определенная зависимость N _ /(V + кдг(V - V),Н, 63,а), полученная в результате заранее проведенного
аэродинамического расчета и расчета характеристик силовой установки. Для аппроксимации зависимости используется искусственная нейронная сеть прямого распространения, содержащая два промежуточных слоя. Коэффициенты кир _ 1/ Nтах и к _ 1/ пттах определяются исходя из рабочих параметров силовой
установки - максимальной мощности Nmсa и максимальной величины оборотов двигателя пдв тса- Коэффициент kдV определяется экспериментально.
За счет использования нейронной сети в автомате тяги, предполагается добиться адаптивности к различным типам БЛА, а также повысить скорость вычислений алгоритмов автомата тяги, по сравнению с классическими алгоритмами. Такие возможности подтверждены проведенными оценками.
Таким образом, изложены основные принципы построения СПП и САУ в рамках УНК. В дальнейшем предполагается их углубленная интеграция с алгоритмами БГНК в рамках комплекса, имитационное и полунатурное моделирование их совместной работы.
Заключение
Итак, в работе изложен новый подход к построению интегрированного высокоточного малогабаритного управляющего навигационно-гравиметрического
комплекса БЛА. Обоснована актуальность и перспективность разработки. Изложены основные положения концепции новой технологии бесплатформенной аэрогравиметрии на базе малоразмерных БЛА. Рассмотрен вариант комплекса в составе системы планирования полета и системы автоматического управления БЛА, а также бесплатформенного гравиинерциального навигационного комплекса. Представлены структуры, состав, математические модели и алгоритмы работы основных узлов комплекса. В итоге можно сделать вывод, что интегрированный комплекс позволит обеспечить решение задач управления, ориентации и навигации БЛА, а также векторных гравиметрических измерений в соответствии с новой технологией аэрогравиметрии.
Работа выполнена при финансовой поддержке Минобрнауки РФ по проектам 14.B37.21.1904 от 04.10.2012г. и 14.B37.21.1545 от 20.09.2012г. ФЦП "Научные и научно-педагогические кадры инновационной России на 2009-2013 годы", а также РФФИ в форме гранта № 12-08-00028-а.
Библиографический список
1. Применение гравиинерциальных технологий в геофизике. Сборник статей и докладов. Под общей редакцией академика РАН В.Г. Пешехонова. -СПб.: ГНЦ РФ ЦНИИ «Электроприбор», 2002.
2. Афонин А.А., Сулаков А.С. Алгоритм оценивания параметров гравитационного поля, ориентации и навигации бесплатформенного гравиинерциального комплекса // Мехатроника, автоматизация, управление. -М.: 2011, № 9.
3. Афонин А.А., Сулаков А.С. О коррекции параметров ориентации, навигации и гравиметрии в замкнутой схеме включения фильтра Калмана навигационного комплекса // Авиакосмическое приборостроение. -М.: 2011, №8.
4. Тювин А.В., Афонин А.А., Черноморский А.И. Об одной концепции векторных гравиметрических измерений // Авиакосмическое приборостроение, 2005. -М.: №3.
5. Афонин А.А., Тювин А.В., Карачевцев М.В. Об оценке погрешностей в определении параметров ориентации в бесплатформенном гравиинерциальном навигационном комплексе // Авиакосмическое приборостроение. - М.: 2007, №3.
6. Андреев В.Д. Теория инерциальной навигации. Автономные системы. -М.: Наука, Главная редакция физико-математической литературы, 1966.
7. Афонин А.А., Сулаков А.С., Тювин А.В., Сачков Г.П. Об одной модификации алгоритма бесплатформенного гравиинерциального навигационного комплекса // Труды XXXVIII «Всероссийской конференции по управлению движением кораблей и специальных подводных аппаратов». -М.: ИПУ РАН, 2012.
8. ГОСТ Р 51794-2001. Аппаратура радионавигационная глобальной спутниковой системы и глобальной системы позиционирования. Системы координат. Методы преобразования координат определения точек. Госстандарт России. Москва.
9. Лебедев Г.Н., Мирзоян Л.А. Нейросетевое планирование действий по облету наземных объектов группой летательных аппаратов // Авиакосмическое приборостроение. -М.: 2005, №12.
10. Лебедев Г.Н., Мирзоян Л.А. Конвейерное нейросетевое управление полетом на базе обучающей моделирующей среды // Материалы Международной научной конференции «Современные проблемы прикладной математики и математического моделирования», -Воронеж: ВГТА, 2005.
11. Лебедев Г.Н., Михайлин Д.А., Ефимов А.В. Оценка вектора состояния беспилотного летательного аппарата при наличии в объекте управления нелинейных элементов // Вестник Московского авиационного института, 2012, т. 19, №1.
12. Михалев И.А., Окоемов Б.Н., Чикулаев М.С. Системы автоматического управления самолетом. -М: Машиностроение, 1971.
13. Горбачев Ю.В., Курмаков Д.В., Рыбников С.И. Интегрированная система автоматического управления продольным движением летательного аппарата в строю при регулировании тягой двигателя и рулем высоты // Электронный журнал «Труды МАИ», 2011, №49.
