Динамика пассивного космического объекта при его бесконтактной транспортировке ионным потоком

Содержание

Слайд 2

Доклады на конференциях V.S. Aslanov, A.S. Ledkov, "Chaotic motion of a

Доклады на конференциях

V.S. Aslanov, A.S. Ledkov, "Chaotic motion of a

passive space object during its contactless transportation by ion beam", ITNT-2020, Samara
V.S. Aslanov, A.S. Ledkov, M. Konstantinov "Influence of space debris attitude motion on ion beam assisted removal mission costs", IAC-2020, Dubai
G. Popov, V.S. Aslanov, V. Petukhov, A.S. Ledkov, "Influence of space debris attitude motion on ion beam assisted removal mission costs", IAC-2020, Dubai
Слайд 3

Синдром Кеслера – лавинообразное увеличение числа обломков КМ 2007 – Разрушение

Синдром Кеслера – лавинообразное увеличение числа обломков КМ
2007 – Разрушение Fengyun-1С
2009

- Столкновение Iridium 33 и Cosmos-2251

Проблема космического мусора

34000 объектов>10 cm
5400 объектов >1m
2000 активных спутников
US Space Surveillance Catalogue

Подходы к активному уделению космического мусора:
стыковка или жесткий захват
захват и тросовая буксировка
бесконтактная транспортировки

Donald J. Kessler

Слайд 4

Бесконтактное воздействие на объект Бесконтактное взаимодействие: Электростатическое Гравитационное Ионный поток Магнитное

Бесконтактное воздействие на объект

Бесконтактное взаимодействие:
Электростатическое
Гравитационное
Ионный поток
Магнитное
Лазер

Преимущества
Отсутствие механического контакта
Возможность транспортировки быстро вращающихся

объектов
Слайд 5

Схема транспортировки ионным потоком Авторы идея: C. Bombardelli и J. Pelaez

Схема транспортировки ионным потоком

Авторы идея:
C. Bombardelli и J. Pelaez (Ion Beam Shepherd)
S.

Kitamura (Ion Beam Irradiation Reorbiter)

Conceptual scheme

[1] Bombardelli C., Pelaez J. Sistema de modificación de la posición y actitud de cuerpos en órbita por medio de satélites guía, Patent No. P201030354, filed 11 March 2010
[2] Kitamura, S., Large Space Debris reorbiter using ion beam irradiation, 61st International Astronautical Congress, International Astronautical Federation, Paris, France, 2010.

Слайд 6

Вычисление ионных сил и моментов [3] A.P. Alpatov, S.V. Khoroshylov, A.I.

Вычисление ионных сил и моментов

[3] A.P. Alpatov, S.V. Khoroshylov, A.I. Maslova,

Contactless de-orbiting of space debris by the ion beam. Dynamics and Control. — Кyiv: Akademperiodyka, 2019. — 170 p. DOI: 10.15407/akademperiodyka.383.170. (Chapter 3)
[4] V.S. Aslanov, A.S. Ledkov, Attitude motion of cylindrical space debris during its removal by ion beam, Mathematical Problems in Engineering. (2017) Article ID 1986374. DOI: 10.1155/2017/1986374.

d=15m β=0

Слайд 7

Задачи исследования Исследовать особенности динамики пассивного цилиндрического объекта при его бесконтактной

Задачи исследования

Исследовать особенности динамики пассивного цилиндрического объекта при его бесконтактной транспортировке

в плоском и пространственном случае.
Найти закон управления ионным потоком, обеспечивающий стабилизацию транспортируемого объекта.
Оценить выигрыш в затратах топлива при управляемой транспортировке.
Слайд 8

Математическая модель. Плоское движение Допущения Активный космический аппарат рассматривается как материальная

Математическая модель. Плоское движение

Допущения
Активный космический аппарат рассматривается как материальная точка
Пассивный объект

рассматривается как цилиндр
Учитывается влияние только гравитационного поля Земли, сил тяги и передаваемых ионным потоком сил и моментов
Слайд 9

Математическая модель Уравнения Лагранжа 2-го рода Обобщенные координаты r - расстояние

Математическая модель

Уравнения Лагранжа 2-го рода

Обобщенные координаты
r - расстояние до центра Земли
ν

- угол истинной аномалии
θ - угол отклонения оси объекта
x - относительная координата КА
y – относительная координата КА

(1)

Обобщенные силы

(2)

Слайд 10

Математическая модель Функция Лагранжа (3) Обобщенные силы (4)

Математическая модель

Функция Лагранжа

(3)

Обобщенные силы

(4)

Слайд 11

Уравнения движения (5)

Уравнения движения

(5)

Слайд 12

Упрощенная математическая модель Допущения: x=const, y=const, 1/r (6) Центр масс движется

Упрощенная математическая модель

Допущения: x=const, y=const, 1/r<<1

(6)

Центр масс движется по Кеплеровой орбите

(7)

(8)

(9)

В

случае малого эксцентриситета e<<1
Слайд 13

Движение по орбите с малым эксцентриситетом (10) Момент ионного потока (11) (12) Невозмущенная система (e=0)

Движение по орбите с малым эксцентриситетом

(10)

Момент ионного потока

(11)

(12)

Невозмущенная система (e=0)

Слайд 14

Параметры моделирования Параметры системы Параметры ионного потока

Параметры моделирования

Параметры системы

Параметры ионного потока

Слайд 15

Рассматриваемые случаи

Рассматриваемые случаи

Слайд 16

Момент ионного потока [3] A.P. Alpatov, S.V. Khoroshylov, A.I. Maslova, Contactless

Момент ионного потока

[3] A.P. Alpatov, S.V. Khoroshylov, A.I. Maslova, Contactless de-orbiting

of space debris by the ion beam. Dynamics and Control. — Кyiv: Akademperiodyka, 2019. — 170 p. DOI: 10.15407/akademperiodyka.383.170. (Chapter 3)
[4] V.S. Aslanov, A.S. Ledkov, Attitude motion of cylindrical space debris during its removal by ion beam, Mathematical Problems in Engineering. (2017) Article ID 1986374. DOI: 10.1155/2017/1986374.
Слайд 17

Бифуркационные диаграммы Case 1 Case 2 Case 3

Бифуркационные диаграммы

Case 1

Case 2

Case 3

Слайд 18

Фазовый портрет невозмущенной системы (Case 3)

Фазовый портрет невозмущенной системы (Case 3)

Слайд 19

Хаотическое движение (Случай 3) e=0.001 e=0.01 Спектр показателей Ляпунова: Λ=[0.1, 0.001, -0.201] (10)

Хаотическое движение (Случай 3)

e=0.001

e=0.01

Спектр показателей Ляпунова: Λ=[0.1, 0.001, -0.201]

(10)

Слайд 20

Вывод Топология фазового пространства может сильно изменяться в процессе спуска. С

Вывод

Топология фазового пространства может сильно изменяться в процессе спуска. С увеличением

эксцентриситета толщина хаотического слоя увеличивается.
Слайд 21

Energy integral Исследование невозмущенного движения

Energy integral

Исследование невозмущенного движения

Слайд 22

Осредненная сила ионного сопротивления

Осредненная сила ионного сопротивления

Слайд 23

Относительное положение КА Направление ионного потока Закон управления двигателями КА

Относительное положение КА

Направление ионного потока

Закон управления двигателями КА

Слайд 24

Управление колебаниями объекта

Управление колебаниями объекта

Слайд 25

Результаты численного моделирования

Результаты численного моделирования

Слайд 26

Результаты численного моделирования

Результаты численного моделирования

Слайд 27

1 2 Результаты численного моделирования

1

2

Результаты численного моделирования

Слайд 28

Управление угловым движением объекта при его транспортировке позволяет существенно снизить затраты

Управление угловым движением объекта при его транспортировке позволяет существенно снизить затраты

топлива
Даже без приема во внимание влияния угловых колебаний объекта, рассматриваемая технология позволяет решать задачу бесконтактной транспортировки

Вывод

Слайд 29

Математическая модель. Системы координат - Инерциальная СК - Орбитальная СК Обобщенные координаты - Связанная СК

Математическая модель. Системы координат

- Инерциальная СК

- Орбитальная СК

Обобщенные координаты

- Связанная СК

Слайд 30

Движение центра масс Гравитационная сила Ионная сила Уравнения движения (1) (2) (3) (4)

Движение центра масс

Гравитационная сила

Ионная сила

Уравнения движения

(1)

(2)

(3)

(4)

Слайд 31

Движение относительно центра масс Момент количества движения Тензор инерции (5) Угловая

Движение относительно центра масс

Момент количества движения

Тензор инерции

(5)

Угловая скорость

Гравитационный момент

Ионный момент

(6)

Слайд 32

Рассмотрим движение на GEO (7) (8) (9) (10) (11) (12)

Рассмотрим движение на GEO

(7)

(8)

(9)

(10)

(11)

(12)

Слайд 33

Уравнения движения ЦМ на GEO Control parameter (13) (14) (15) (16)

Уравнения движения ЦМ на GEO

Control parameter

(13)

(14)

(15)

(16)

Слайд 34

Положение равновесия (8) (11) (17)

Положение равновесия

(8)

(11)

(17)

Слайд 35

Параметры системы

Параметры системы

Слайд 36

Момент ионного потока [1] A.P. Alpatov, S.V. Khoroshylov, A.I. Maslova, Contactless

Момент ионного потока

[1] A.P. Alpatov, S.V. Khoroshylov, A.I. Maslova, Contactless de-orbiting

of space debris by the ion beam. Dynamics and Control. — Кyiv: Akademperiodyka, 2019. — 170 p. DOI: 10.15407/akademperiodyka.383.170. (Chapter 3)
[2] V.S. Aslanov, A.S. Ledkov, Attitude motion of cylindrical space debris during its removal by ion beam, Mathematical Problems in Engineering. (2017) Article ID 1986374. DOI: 10.1155/2017/1986374.
Слайд 37

Невозмущенное движение

Невозмущенное движение

Слайд 38

Закон управления тягой ионного двигателя Закон управления Функция Хевисайда

Закон управления тягой ионного двигателя

Закон управления

Функция Хевисайда

Слайд 39

Результаты численного моделирования Начальные условия

Результаты численного моделирования

Начальные условия

Слайд 40

Case a Case b Результаты численного моделирования

Case a

Case b

Результаты численного моделирования

Слайд 41

Результаты и выводы Исследовано невозмущенное движение системы Предложен закон управления тягой

Результаты и выводы

Исследовано невозмущенное движение системы
Предложен закон управления тягой ионного двигателя.
Эффективность

закона подтверждаете результатами численного моделирования.