Астероидно-кометная опасность: вчера, сегодня, завтра

Таким образом, через 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.

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

^{Рис. 10.4. Результат длительного воздействия постоянной тяги dg_W, приложенной по оси W}

Обратимся к двум выбранным ранее примерам. В первом случае увода на расстояние, равное диаметру Земли, положив dn_max/R_э = 2, получим, что требуемое значение ускорения dg_W равно: dg_W /g₀ = (R_э/r₀) = 4,25 10^-5. Например, для Апофиса необходимая тяга составит~ 11,5 кН (килоньютон). Следовательно, такая тяга, приложенная к Апофису в течение полугода, даст изменение текущих координат на диаметр Земли. Во втором случае (увод из зоны резонансного возврата) тяга, прилагаемая также в течение полугода, оказывается в 1000 раз меньше, т. е. составит всего лишь около 11,5 Н.

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

Следовательно, изменения координат относительно старых значений описываются формулами