Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.
Для эффективного применения космической техники, выполняющей задачи на поверхности других планет, необходимо осуществлять ее посадку с заданной точностью. Однако наличие различных возмущений в атмосфере планеты приводит к отклонению траектории его движения от расчетной. Для компенсации возмущений и реализации принципа управляемой посадки известен метод управления движением спускаемого аппарата путем смещения центра масс. В результате изменения положения центра масс, спускаемый аппарат будет двигаться под некоторым углом к набегающему потоку, из-за чего возникнут перегрузки перпендикулярные направлению движения. Используя эти перегрузки, можно компенсировать возмущения, то есть осуществить управление движением.
Цель настоящей статьи — провести исследование маневренных свойств спускаемого аппарата управляемого путем смещения центра масс.
В этой статье рассмотрен спускаемый аппарат, состоящий из надувного тормозного устройства, космического аппарат и поворотного устройства (рис. 1). Надувное тормозное устройство представляет из себя конус с углом полураствора 60° со сферическим затуплением. Геометрия надувного тормозного устройства и его размеры представлены на рис. 2.
Внутри надувного тормозного устройства установлен космический аппарат и поворотное устройство. Последнее предназначено для поворота космического аппарата вокруг точки их сопряжения. За счет поворота космического аппарата происходит смещение центра масс спускаемого аппарата.
Для исследования движения спускаемого аппарата вокруг центра масс при изменении положения космического аппарата составлена система дифференциальных уравнений.
На рис. 3 представлена схема движения спускаемого аппарата. Здесь введены следующие обозначения: — полная аэродинамическая сила; —аэродинамическая продольная сила; — аэродинамическая нормальная сила; — угол атаки; — угол поворота космического аппарата; НТУ — надувное тормозное устройство; КА — космический аппарат; СА — спускаемый аппарат.
В общем случае векторное уравнение вращательного движения твердого тела можно записать в следующем виде:
(1)
где — вектор кинетического момента; — тензор инерции; — вектор угловой скорости спускаемого аппарата.
Перепишем уравнение (1) в виде:
Тогда динамические уравнения вращательного движения можно представить в виде
где — сумма моментов внешних сил.
Так как из всех внешних сил только аэродинамические создают момент, то сумма моментов запишется следующим образом:
Здесь — аэродинамический момент крена; — аэродинамический момент рыскания; — аэродинамический момент тангажа; — характерный размер; — демпфирующий момент; — коэффициент аэродинамической продольной силы; — коэффициент аэродинамической нормальной силы; — скоростной напор; — площадь миделя; — коэффициент демпфирующего момента; — скорость набегающего потока.
Коэффициент демпфирующего момента рассчитаем по формуле [1, с. 36]:
где — угол полураствора конуса; — координата точки от носка конуса, относительно которой рассчитывается коэффициент демпфирующего момента; — коэффициент давления; — отношение радиуса сферического затупления к радиусу миделя; — диаметр миделя; — угол атаки.
В качестве кинематических уравнений вращательного движения взяты кинематические уравнения Эйлера:
Здесь — угол крена; — угол рыскания; — угол тангажа.
Тензор инерции спускаемого аппарата складывается из тензора инерции надувного тормозного устройства и тензора инерции космического аппарата , который может поворачиваться относительно осей системы координат, связанной с надувным тормозным устройством.
Тензор инерции спускаемого аппарата рассчитывается следующим образом:
Здесь — тензор инерции надувного тормозного устройства; — тензор инерции космического аппарата при ; — матрица перехода из системы координат связанной с космическим аппаратом в систему координат связанной с надувным тормозным устройством.
Движение спускаемого аппарата моделировалось в стационарном набегающем потоке со скоростным напором Па при скорости поворота космического аппарата 0,001 рад/с.
Решение системы дифференциальных уравнений проведено численно, методом Рунге — Кутта 4-го порядка. Результаты интегрирования представлены на рис. 4–6.
На рис. 4 приведена зависимость изменения угла атаки от угла поворота космического аппарата. Видно, что величина установившегося угла атаки превышает величину угла поворота космического аппарата. Полученный результат зависит от положения центра масс элементов спускаемого аппарата по отношению к центру давления.
На рис. 5 представлена зависимость нормальной перегрузки от величины угла атаки. Анализируя зависимость можно заметить, что зависимость имеет линейный характер при изменении угла атаки от 0 до 15°. В этом диапазоне изменения угла атаки нормальная перегрузка изменяется от 0 до 40 единиц.
На рис. 6 представлена зависимость нормальной перегрузки от угла поворота космического аппарата. Зависимость имеет более короткий линейный участок, чем зависимость нормальной перегрузки от угла атаки. Связано это с тем, что зависимость угла атаки от угла поворота КА не линейная.
Рассмотрим изменение величины угла атаки от времени, при повороте космического аппарата на угол равный 4°. Зависимость изменения угла атаки от времени представлена на рис. 7. Можно заметить, что после достижения значения установившегося угла атаки, спускаемый аппарат совершает колебания вокруг этого установившегося угла. Колебания препятствуют процессу управления движением спускаемого аппарата. Следовательно, для обеспечения хорошего качества управления необходимо применить устройства увеличивающие демпфирующий момент СА.
Применим динамическое демпфирование для гашения колебаний.
Для реализации динамического демпфирования воспользуемся следующим способом. При прохождении спускаемым аппаратом положения равновесия будем создавать момент, направление которого противоположно направлению вращательного движения. В результате создания такого момента угловая скорость вращательного движения будет уменьшаться.
Схема описанного метода гашения колебаний приведена на рис. 8.
На рис. 9 представлена зависимость изменения угла атаки от времени. Из приведенной зависимости видно, что при применении динамического демпфирования амплитуда колебаний спускаемого аппарата вокруг положения равновесия на шестой секунде не превышает 1°.
В результате выполнения работы составлена математическая модель вращательного движения спускаемого аппарата с изменяющимся центром масс.
По результатом интегрирования математической модели получено:
Для гашения автоколебаний возникающих при управлении движением спускаемого аппарата применен метод динамического гашения колебаний.
