Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.
Как правило, орбиту любого движущегося тела в космосе можно охарактеризовать шестью основными математическими элементами, которые однозначно определяют конкретную орбиту. В качестве альтернативы для описания орбиты можно использовать три дифференциальных уравнения второго порядка относительно векторов положения и скорости точки. Однако такой вариант допускает аналитические решения в замкнутой форме только для модели Земли в виде однородного шара. Учет возмущений обычно усложняет систему дифференциальных уравнений так, что получить точные аналитические решения становится невозможным.
С целью учета отличия гравитационного поля Земли от поля однородного шара, гравитационный потенциал Земли может быть построен на основе модели с использованием зональных сферических гармоник, из которых J2 имеет наибольшее значение по причине ее сравнительно большой величины [1]. Соответствующая J2 потенциальная функцияϕJ2 имеет вид
где Re — радиус Земли; r — радиальное расстояние; z — полярная компонента радиус-вектора спутника r.
Ускорение из-за возмущения J2 получается путем взятия градиента потенциальной функции ϕJ2. Тогда ускорение в инерциальной системе отсчета в декартовых координатах определяется как вектор
Уравнения относительного движения ОКМ под действием силы тяги, магнитного и гравитационного возмущений записываются в геоцентрической инерциальной системе отсчета. Уравнения выводятся с использованием декартовых координат. При этом
Здесь вектор ag — ускорение, обусловленное гравитацией и полученное на основе потенциала однородного шара; вектор u — ускорение, обусловленное силой тяги; вектор aL — ускорение, вызванное силой Лоренца, действующей на ОКМ, движущийся в магнитном поле коллектора [2].
Выражение для силы Лоренца, действующей на проводящий сферический объект в локальной системе LVLH O (r, α, h), записывается как
где параметр — проводимость ОКМ; — средний радиус проводящей поверхности ОКМ; — плотность ОКМ; — индукция магнитного поля коллектора.
Векторы, заданные в локальной системе координат LVLH, а именно ускорение из-за силы Лоренца и вектор силы тяги, могут быть преобразованы в геоцентрическую инерциальную систему отсчета.
Для анализа движения системы “коллектор + ОКМ” находится ее полная энергия, которая в данном случае эквивалентна Гамильтониану . Учитывая, что потенциал ϕJ2 не нарушает консервативности гравитационного поля, Гамильтониан должен оставаться постоянным [3]. Векторы относительного смещения и скорости рассчитываются как разности между соответствующими векторами для ОКМ и коллектора. Анализ построенной нелинейной математической модели производится с использованием численного интегрирования дифференциальных уравнений движения коллектора и ОКМ с различными наборами начальных условий. Результирующие векторы могут быть преобразованы в переменные состояния относительного движения, обычно выражаемые в системе отсчета LVLH, центр которой связан с коллектором.
Новизна работы заключается в использовании нового бесконтактного электромагнитного способа захвата и удержания ОКМ с целью его последующего удаления из околоземного космического пространства. Применение электромагнита, питаемого от бортовых аккумуляторов коллектора, в сочетании с солнечной энергией для создания достаточно сильного магнитного поля, способного захватить и удерживать металлические ОКМ в космосе, не обсуждалось в предыдущих работах, связанных с удалением космического мусора [4]. Результаты численного моделирования показывают эффективность предложенного способа удаления ОКМ. Учет возмущающего воздействия второй зональной гармоники геопотенциала не нарушает работоспособности предложенного способа. Представленный способ выгодно отличается от некоторых известных электродинамических способов тем, что обеспечивает контроль над траекторией удаляемого ОКМ.
