АЛГОРИТМЫ ОПРЕДЕЛЕНИЯ КУРСОВОГО УГЛА ДЛЯ ВЫЧИСЛЕНИЯ ФАЗОВОГО ЦЕНТРА АНТЕННЫ КОСМИЧЕСКОГО АППАРАТА «ГЛОНАСС» НА УЧАСТКАХ УПРЕЖДАЮЩЕГО РАЗВОРОТА

Полный текст

Полная орбитальная группировка ГЛОНАСС содержит 24 КА на круговых орбитах с наклонением 64,8° в трех орбитальных плоскостях по восемь КА в каждой. Долгота восходящих узлов трех орбитальных плоскостей различается номинально на 120°. В каждой орбитальной плоскости восемь КА разнесены по аргументу широты через 45°. Такая конфигурация подсистем КА позволяет обеспечить непрерывное и глобальное покрытие земной поверхности и околоземного пространства навигационным полем [1]. Орбита навигационного КА системы ГЛОНАСС имеет особые точки (при углах Солнце-объект-Земля (СОЗ) близких к 0° - малые углы СОЗ и близких к 180° - большие углы СОЗ), при прохождении которых КА совершает разворот вокруг оси, ориентированной на Землю (ось минус ОХ). С целью уменьшения влияния сил солнечного давления на движение центра масс КА и для минимизации погрешности отслеживания курсового угла, КА совершает симметричный упреждающий разворот вокруг оси минус ОХ. При прохождении больших углов СОЗ при угле между плоскостью орбиты и направлением на Солнце, меньше заданного значения, формируется упреждающий разворот относительно оси, направленной на Землю, до момента прохождения максимального угла СОЗ. При угле между осью минус OX и направлением на Солнце, меньше заданного значения определяются: - угол в плоскости орбиты между текущим положением КА и положением КА в момент прохождения максимального угла СОЗ; - угол разворота относительно оси минус OX в процессе движения по орбите от текущего положения КА до точки прохождения максимального угла СОЗ; - время разворота по оси минус OX на указанный выше угол при заданной скорости разворота; - время прохождения по орбите угла от текущего положения до точки прохождения максимального угла СОЗ. Как только время, необходимое для разворота вокруг оси минус ОХ на вычисленный угол разворота, становится больше времени прохождения по орбите угла от текущего положения КА до точки прохождения максимального угла СОЗ, формируется признак больших углов СОЗ и признак на упреждающий разворот КА. Выключение упреждающего разворота осуществляется после прохождения больших углов СОЗ при угле между плоскостью XOY и плоскостью СОЗ меньше заданного значения. Малые углы СОЗ характеризуются наличием тени Земли. Приближение теневого участка определяется путем сравнения угла между осью минус ОХ и направлением на Солнце с угловым размером Земли, если угол между плоскостью орбиты и направлением на Солнце меньше углового размера Земли. При входе в тень отслеживание плоскостью направления на Солнце осуществляется по разности углов: угла, полученного интегрированием углов скорости относительно оси минус ОХ после входа КА в тень Земли, и угла, равного разности текущего угла рыскания и угла рыскания на момент входа КА в тень Земли, рассчитанных программным комплексом баллистических задач (ПК БЗ). При прохождении в тени Земли малых углов СОЗ, при которых скорость вращения плоскости СОЗ превышает поисковую скорость относительно оси минус ОХ, формируется упреждающий разворот относительно этой оси до момента прохождения минимального угла СОЗ (середины теневого участка). При угле СОЗ, меньше заданного значения, с использованием информации ПК БЗ определяются: - угол в плоскости орбиты между текущим положением КА и положением КА в момент прохождения минимального угла СОЗ; - угол разворота относительно оси минус OX от текущего положения КА до точки прохождения минимального угла СОЗ; 199 Вестник СибГАУ. N 4(50). 2013 - время разворота (Ti) относительно оси минус OX на угол, определенный в предыдущем пункте, при заданной скорости разворота; - время прохождения (T2) по орбите угла от текущего положения до точки прохождения минимального угла СОЗ. При T1 > T2 формируется признак упреждающего разворота. Скорость разворота ~0,25 °/с, время разворота на угол 180° не превышает 15 минут. Угол между плоскостью орбиты и направлением на Солнце, при котором формируется упреждающий разворот, не превышает 2°. Угол СОЗ, при котором разрешается формирование упреждающего разворота, не более 5°. Выключение упреждающего разворота осуществляется после прохождения малых углов СОЗ при угле между плоскостью XOY и плоскостью СОЗ, меньше заданного значения. Для формирования признака упреждающего разворота и знака разворота используются следующие баллистические параметры: - текущий угол СОЗ; - текущий угол между плоскостью орбиты и направлением на Солнце. Положение связанных осей КА относительно солнечно-земной и орбитальной систем координат перед формированием упреждающего разворота показано на рис. 1 [4]. Ниже представлена математическая модель вычисления курсового угла при упреждающем развороте: Ψ x = Ψ н Ψ1 + “n ■ t + - ,0 < t < L ■ t, ^разг < t < ^разг + 2tn Ψ1 +ψ 2 + “n *( - ^разг - 2tn )- “ ■ (t ^разг 2tn ) t + 2t < t < 2t + 2t ^разг ' ^Ln — -^разг ' ^Łn ω ■t. где Ψ1 = Ψнач + “нач ■ *разг +' разг Ψ2 =“n ■ (разг + 2tn ). В модели упреждающий разворот разделен на три взаимосвязанных участка: - разгон (участок набора поисковой скорости) -этот участок описывает движение КА вокруг оси минус ОХ (ось направленная на Землю) при наборе поисковой скорости (максимально допустимая скорость вращения аппарата). Набор поисковой скорости ха-рактрезуется наличием постоянного углового ускорения; - разворот на поисковой скорости - описывает движение КА вокруг оси минус ОХ с постоянной максимально допустимой скоростью вращения КА; - торможения - описывает движение КА вокруг оси минус ОХ при уменьшениии угловой скорости КА с постоянным отрицательным ускорением. На рис. 2 показан принцип работы алгоритма, позволяющего заблаговременно сформировать начальные данные, необходимые для выполнения упреждающего разворота. Χ,Χο,Χζ Рис. 1. Положение связанной системы координат (OXYZ) относительно солнечно-земной (OXzYzZz) и орбитальной (OX0Y0Z0) систем координат: as - угол между плоскостью орбиты и направлением на Солнце; СОЗ - текущий угол Солнце-объект-Земля; E - угол по орбите от текущей точки до точки минимального угла СОЗ; ψ^ - значение курсового угла перед началом упреждающего разворота; ψ^ -угол разворота относительно оси ОХ в процессе движения по орбите от текущего положения КА до точки прохождения максимального угла СОЗ, в которой ось OY должна быть перпендикулярна плоскости орбиты Сравнение текущего угла между плоскостью орбиты и направлением на Солнце с пороговым значением Сравнение текущего угла СОЗ с пороговым значением Расчет скорости изменения угла между плоскостью орбиты и направлением на Солнце Расчет угла между плоскостью орбиты и направлением на Солнце на момент упреждающего разворота_ Расчет времени начала и конца упреждающего разворота (^Конец программьГ^) Рис. 2. Блок-схема алгоритма прогнозирования начальных данных упреждающего разворота 200 Авиационная и ракетно-космическая техника Участок торможения Отслеживание модельного угла разворота и модельной скорости разворота Выдача управляющего сигнала на двигатели-маховики Сравнение текущего угла разворота с углом окончания разворота Получение информации о времени начала упреждающего _разворота_ Обработка начальных данных для упреждающего разворота, полученных с KA Сравнение текущего времени и времени начала разворота Участок разгона Расчет текущего угла разворота Участок разворота на поисковой скорости Расчет текущего угла разворота Вычисление угла разворота для участка торможения Участок торможения Расчет текущего угла разворота и текущего значения скорости разворота Сравнение текущего угла разворота с углом окончания разворота Расчет поправок к координатам фазового центра антенны Сравнение текущего времени со временем необходимым для на поисковой скорости Рис. 3. Блок-схема алгоритма вычисления курсового угла при упреждающем развороте на борту космического аппарата Рис. 4. Блок-схема алгоритма вычисления в аппаратуре потребителя поправки положения фазового центра антенны Сформированные начальные данные передаются в аппаратуру потребителя через навигационный кадр. Эталонные модели упреждающего разворота на борту и в аппаратуре потребителя одинаковые, следовательно, время начала упреждающего разворота, фактический угол разворота КА по каналу рыскания и расчетный угол разворота КА по каналу рыскания на борту и в аппаратуре потребителя будут также одинаковы. На КА управляющие двигатели-маховики (УДМ) имеют разброс по крутизне моментной характеристики. Момент сопротивления на валу УДМ также иска жает моментную характеристику УДМ. Момент инерции КА относительно оси ОХ также имеет разброс [2]. • Mx Исходя из этого, величина CD = —— не имеет по- Jx стоянного значения. Поэтому при проведении упреждающего разворота будут возникать существенные погрешности по углу разворота КА. Для минимизации погрешностей, обусловленных разбросом характеристик УДМ и разбросом момента инерции КА относительно оси ОХ, рассматривается 201 Вестник СибГАУ. № 4(50). 2013 алгоритм формирования управляющих воздействий по каналу ОХ, с дополнительным формированием управляющих воздействий на двигатель-маховик по оси ОХ при рассогласовании по углу и скорости фактического и модельного движения КА. На рис. 3 представлен принцип работы алгоритма вычисления курсового угла при упреждающем развороте на борту навигационного КА. Потребитель проводит измерения положения фазового центра антенны, а должен знать положение центра масс КА. В связи с этим при прохождении особых точек орбиты в аппаратуре потребителя необходимо вычислять поправки на положение фазового центра антенны. На рис. 4 представлен принцип работы алгоритма вычисления в аппаратуре потребителя поправки положения фазового центра антенны. Алгоритмы, описанные в данной статье, являются универсальными для всех типов космических аппаратов навигационной группировки ГЛОНАСС. Погрешность определения ориентации по каналу рыскания с использованием разработанных алгоритмов при прохождении упреждающего разворота не превышает 15', при этом погрешность определения положения фазового центра антенн в плоскости OY0Z0 относительно номинального положения не бу дет превышать 4,5 мм при расстоянии от оси минус ОХ до центра антенны 1 м.

Об авторах

Алексей Владимирович Фатеев

ОАО «Информационные спутниковые системы» имени академика М. Ф. Решетнева»

Email: mr.fateev@inbox.ru
инженер 2 категории Российская Федерация, 662972, г. Железногорск Красноярского края, ул. Ленина, 52

Данил Витальевич Емельянов

ОАО «Информационные спутниковые системы» имени академика М. Ф. Решетнева»

инженер-конструктор 1 категории Российская Федерация, 662972, г. Железногорск Красноярского края, ул. Ленина, 52

Юрий Александрович Тентилов

ОАО «Информационные спутниковые системы» имени академика М. Ф. Решетнева»

начальник сектора систем ориентации и стабилизации (9352) Российская Федерация, 662972, г. Железногорск Красноярского края, ул. Ленина, 52

Андрей Викторович Овчинников

ОАО «Информационные спутниковые системы» имени академика М. Ф. Решетнева»

начальник отдела систем ориентации и стабилизации (935) Российская Федерация, 662972, г. Железногорск Красноярского края, ул. Ленина, 52

Михаил Васильевич Лукьяненко

Сибирский государственный аэрокосмический университет имени академика М. Ф. Решетнева

Email: sibgau-sau@mail.ru
профессор, заведующий кафедрой систем автоматического управления Российская Федерация, 660014, Красноярск, просп. им. газ. «Красноярский рабочий», 31

Список литературы

Дополнительные файлы