Точки либраций в небесной механике и космодинамике - Маркеев А.П.
Скачать (прямая ссылка):
§ 6. Оценка точности построенной теории движения КА
6.1. Общие замечания. Среди всего множества траекторий
движения КА в окрестности Ь2 наибольший практический интерес
представляют условно-периодические траектории. Наиболее существенным достоинством предложенного в настоящей главе аналитического метода расчета движения КА является возможность приближенного описания многообразия таких траекторий.
294
ДВИЖЕНИЕ КОСМИЧЕСКОГО АППАРАТА ВБЛИЗИ U
[ГЛ. 14
Но поскольку реализованная процедура аналитического метода расчета лишь приближенно описывает движение КА, то тем самым и многообразие условно-периодических траектррий определено также приближенно. В частности, в точном решении задачи о движении КА, определяемом начальными данными, соответствующими условно-периодическому движению приближенной задачи, неизбежно будут присутствовать экспоненциально возрастающие функции времени. Для оценки точности приближенного метода было проведено сравнение с результатами численного интегрирования строгих уравнений движения в декартовых координатах. Эти результаты рассматривались как эталонные. Вообще говоря, достаточно точное вычисление координат КА в окрестности неустойчивой особой точки с помощью численного интегрирования также является некоторой проблемой, так как методические ошибки аппроксимации и ошибки округления экспоненциально возрастают. Оценки показали, что их суммарная погрешность на интервале 10 сут. не превышает примерно 10 м, что существенно меньше ошибок приближенного метода. Поэтому для наших целей результаты численного интегрирования можно принять за эталон.
Эталонные расчеты проводились в двух вариантах.
1*) Движение КА определяется решением точной системы уравнений ограниченной задачи трех тел (Земля — Луна — КА). Движение Луны эллиптическое.
2) Движение КА определяется решением точной системы уравнений ограниченной задачи четырех тел (Земля — Луна — Солнце — КА). Геоцентрические координаты и компоненты вектора скорости Луны rL (t), у і (t) определяются численным интегрированием уравнений задачи трех тел (Земля — Луна — Солнце).
Первый из этих вариантов служит для выявления возможных грубых ошибок алгоритма и Оценки методических ошибок, связанных с учетом в приближенной теории предыдущего параграфа членов лишь до конечного порядка малости. Второй вариант оценивает степень пригодности построенной в предыдущем параграфе теории в случае эллиптической задачи трех тел для описания реального движения при наличии солнечных возмущений. В обоих варйантах сравнение результатов проводилось на траекториях, соответствующих приближенно условно-периодйческим траекториям = Р*М = 0, для различных начальных значений g*(l\ ро(1) (і = 2, 3).
6.2. Результаты численных экспериментов в эллиптической задаче. Сравнение результатов, полученных приближенным методом, с результатами эталонных расчетов проводилось для трех траекторий, определяемых такими начальными данными при * — to'-
дЧі) = рЧі) = 0, д*<2> = p*V> = дЧз) = рЧз> = 6,
ОЦЕНКА ТОЧНОСТИ ТЕОРИИ ДВИЖЕНИЯ КА:
295
где 6 — 0,005, 0,01 и 0,02, Эти значения 6 соответствуют орбитам, удаленным от точки L2 примерно на 2000, 4000 и 8000 км.
Обозначим через г' (t — t0) и тэ (t — to) геоцентрические радиусы-векторы КА, вычисленные соответственно в приближенном и эталонном расчетах. Ошибку приближенного метода расчета будем оценивать величинами компонент Ах, А у, A z вектора
Дг = г' (t — t0) —r'3(t — t0).
В результате вычислений с разлйчными значениями 6 были получены зависимости компонент вектора Аг от времени. Эти зависимости представлены на рис. 45 сплошными линиями. Видно, что
Рис. 45. Ошибки аналитической теории в случае эллиптической ограниченной задачи трех тел.
ошибка достигает величины 50 — 70 км за время 15, 10 и 5 сут соответственно для б = 0,005, 0,01 и 0,02. В последующие моменты времени ошибки растут очень быстро и при t — to — 20, 15 и 10. сут соответственно превосходят 500 км.
Если в трех рассмотренных вариантах зафиксировать интервал времени t — t0 = 7 сут (это примерно половина периода либрации КА в окрестности L% системы Земля — Луна), то максимальные ошибки составят 0,9, 11 и 210 'км соответственно. Из-за влияния неучтенных членов нормализованного гамильтониана можно ожидать ошибку, пропорциональную б4. С этим примерно согласу-
296
ДВИЖЕНИЕ КОСМИЧЕСКОГО АППАРАТА ВБЛИЗИ Х2
[ГЛ. 14
ется зависимость изохронной ошибки от б, выявленная в результате проведенных расчетов.
Дополнительный расчет с начальными данными, соответствующими точке либрации (6 = 0), показал, что в этом случае максимальная ошибка 50 км достигается только через 20 сут (см. рис. 45, г), а в конце семисуточного интервала ошибка не превосходит 60 м.
На рис. 45 пунктирными линиями показана зависимость ошибок от времени для случая, когда приближенные расчеты выполнены в рамках линейной теории, т. е. в уравнениях движения отброшены нелинейные члены относительно переменных
Из рисунков видно, что в этом случае с увеличением б вдвое изохронная ошибка увеличивается в четыре раза. Сравнение сплошных и пунктирных кривых показывает, как и следовало ожидать, существенную роль нелинейных членов.
