Астероидно-кометная опасность: вчера, сегодня, завтра - Шустов Борис (читать книги бесплатно TXT) 📗

где P — невозмущенный период обращения астероида по орбите, а амплитуда колебаний равна

^{Рис. 10.1. Влияние импульса скорости, приложенного по оси W, на орбиту астероида}

Выберем в качестве удобного масштаба текущих отклонений небесного тела в линейных единицах экваториальный радиус Земли R_э = 6378 км. Тогда можно получить выражения для величины максимального смещения dn_max по оси Z и относительного приращения скорости dV_Z/V₀, необходимого для обеспечения смещения dn_max:

Отметим, что фаза колебания останется привязанной к точке коррекции орбиты небесного тела.

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

В первом из них примем, что орбита проходит через центр Земли, и для избежания удара необходимо получить расстояние увода, равное величине 2R_э (коэффициент 2 берется для гарантии). Нетрудно рассчитать, что для выполнения условия избежания удара dn_max/R_э ≥ 2 необходимо обеспечить величину относительного приращения скорости, равную dV_Z/V₀ = 85 10^-6. Принимая для оценки требуемого изменения скорости значение V₀ = 30 км/с, получим минимальное значение требуемого приращения скорости, равное dV_Z = 2,6 м/с.

Во втором случае примем, что необходимо получить отклонение угрожающего тела на ∼ 12 км, т. е. на величину порядка 0,001R_э. Это значение соответствует уходу от конкретной зоны резонансного возврата протяженностью ∼ 1 км, реализуемому с большим запасом. Очевидно, что теперь величина необходимого импульса коррекции скорости уменьшится в тысячу раз и составит величину всего лишь порядка нескольких миллиметров в секунду (dV_Z = 2,6 мм/с).

Рассматривая результат коррекции орбиты по нормали, можно сразу видеть основной недостаток такой коррекции — ее периодичность и результативность лишь в относительно небольшие интервалы времени, как это следует из характера изменения величины dZ_W. Действительно, максимальное значение увода существует практически в течение ∼ 1/6 периода обращения тела, и оно наступит лишь через 1/4 этого периода.

Далее рассмотрим результат приложения импульса скорости по радиусвектору поражающего тела (рис. 10.2).

^{Рис. 10.2. Влияние импульса скорости, приложенного по оси S, на орбиту астероида}

Обращаясь к производным, приведенным в работе [Эльясберг, 1965], сразу получаем, что небесное тело будет совершать периодические колебания dr, dl относительно текущего невозмущенного положения по радиус-вектору r (т. е. по координате S) и вдоль орбиты l (практически по координате T):

Можно видеть, что и здесь результат коррекции периодичен, но колебания небесного тела возникают уже по двум координатам. Рассматривая отклонение тела в подвижной относительной системе координат, привязанной к текущему первоначальному движению поражающего астероида, можно видеть, что результат коррекции представляет собой эллипс с полуосями dV_r/V₀ и 2(dV_r/V₀) (справа внизу на рис. 10.2).

Таким образом, через 1/4 витка расстояние увода от начальной точки составит √5 (dV_r/V₀), причем в дальнейшем, в течение полупериода обращения, эта величина увода будет увеличиваться. Ее максимальное значение будет иметь место через полпериода обращения тела, отсчитываемого от момента коррекции, и составит 4(dV_r/V₀). При положительном направлении импульса dV_r (от центра орбиты) отклонение по оси T будет направлено в сторону, противоположную движению астероида, т. е. он будет отставать от точки коррекции.

Сравнивая эти результаты с предыдущими, видим, что эффективность коррекции, определяемая по максимальному расстоянию увода, повысилась в 4 раза. Однако и в этом случае периодический характер увода сохраняется. Таким образом, создание импульса коррекции по осям W или S в общем случае полностью задачу увода не решает.

Рассмотрим результаты коррекции скорости, проводимой по оси Т. И в этом случае результатом коррекции является изменение тех же двух координат в плоскости орбиты, но имеющее уже существенно другой характер. Приложение импульса по оси T изменяет период обращения угрожающего тела на постоянную величину.

Как следствие, по координате T (вдоль орбиты), кроме периодических колебаний, появляется дополнительный уход, линейно нарастающий со временем. Характер изменения текущих координат показан на рис. 10.3, а полные выражения для компонент увода имеют следующий вид:

^{Рис. 10.3. Влияние импульса скорости, приложенного по оси T, на орбиту астероида}

Можно видеть, что радиальная компонента увода осталась прежней (правда, с другим знаком), зато у другой компоненты появился линейно нарастающий член, а колебания стали вдвое больше по амплитуде. Тогда средний увод по координате T за один виток будет составлять: dl/R_э = 6π(r₀/R_э)(dV_T /V₀), а сам увод станет пропорционален количеству витков. Любопытно отметить, что при величине отношения dV_T /V₀, обеспечивающей целочисленное отношение (1/3)(V₀/dV_T) = N_в, астероид через N_в витков вернется в ту же точку. Можно видеть, что условием возврата является то, что отношение V₀/dV_T должно быть кратно трем.

Этот случай представляет собой своеобразный вариант резонансного возвращения. Правда, практически это произойдет не скоро (как минимум, через десятки тысяч витков орбиты), и по прошествии значительного времени, вследствие возмущений реальной орбиты, возврата в ту же самую точку почти наверняка не произойдет.

Обращаясь к двум последним примерам, рассмотренным выше, получим, что в первом случае требуемый увод можно получить за один виток, создав приращение скорости, равное всего лишь 0,14 м/с. Во втором случае требуемое приращение будет выражаться малыми долями миллиметра в секунду и составит всего лишь 0,14 мм/с.

10.4. Эффективность непрерывного воздействия на орбиту астероида

Сначала рассмотрим результат длительного воздействия постоянной тяги dg_W, приложенной по оси W (рис. 10.4).

Можно видеть, что результатом приложения тяги вдоль оси W является поворот плоскости орбиты относительно начальной точки приложения непрерывного воздействия. Заметим, что время приложения тяги принималось равным периоду обращения небесного тела.

Решение линеаризованных уравнений имеет вид:

Здесь g₀ — ускорение астероида, вызванное притяжением Солнца: g₀ = µ_c/r₀²,

где µ_c — гравитационный параметр Солнца, а r₀ — радиус орбиты астероида. Величина dg_W представляет собой ускорение тела массой M, создаваемое тягой F. Например, по имеющимся оценкам, масса Апофиса составляет M = 4,6 10¹⁰ кг. Схема изменения орбиты под воздействием тяги dg_W показана на рис. 10.4.