Квант

43

Л. Баканина

Задачи о спутниках

4 октября 1957 года по правау считают началом нового века в освенении Космosa. Запуск первого искусственного спутника Земли открыл невообразимые ранее просторы для научных изысканий. Стали возникными исследованиями космических изучений, не возможности атмосферой. Человек получил возможность по-новому, so стороны взглянуть на свою родную планету: на ее поверхности, атмосферу и даже глубины.

Искусственные спутники породили множество проблем, связанных с их собственным движением, и вызвали новую волну интересса к вопросам, которые занимались еще Кеплер, Галилей, Ньюотон.

Точный расчет движения спутника — задачарезвычайно сложная. Однако в первом приближении она совпадет к задаче взаимодействия двух тел. Именно такие задачи и будут рассмотрены в этой статье.

Вспомним основные физические законы, которые необходимо знать при решении задач о движении спутников.

Прежде всего это закон всемирного тяотения — любые две материалы точки прятягивают друг друга с силой F⃗\vec{F}, направленной по линии, их соединяющей, а по абсолотной величины прямого пропорциональной их массам m1m_1 и m2m_2 и обратно пропорциональной квадрату расстояния rr между ними:

∣F⃗∣=Gm1m2r2|\vec{F}| = G \frac{m_1 m_2}{r^2}

(G=6,67⋅10−11 m3/(кг⋅с2)G = 6,67 \cdot 10^{-11} \text{ m}^3 / (\text{кг} \cdot \text{с}^2) — гравиитционная постоянная).

По поводу этого закона в третьем номере «Кванта» за 1980 год была опубlikована статья В. Можаева. В ней было показано, в частности, что потенциальная энергия гравиитонного взаимодействия тел равна

Ep=−Gm1m2r2E_p = -G \frac{m_1 m_2}{r^2}

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

ΔEp=mgΔh\Delta E_p = mg\Delta h

является частьим случае приведенной выше обшей формулы. Она справедлива, если рассматривается взаимодействие тела массой mm с Землей, причем тело находится в близи поверхности Земли и перемещается на высоту Δh\Delta h, малую по сравнению с радиусом Земли. В этом случае

ΔEp=(−GmM3R3+Δh−(−GmM3R3))=GM3mΔh(R3+Δh)R3≈mGM3R32Δh.\Delta E_p = \left( -G \frac{m M_3}{R_3 + \Delta h} - \left( -G \frac{m M_3}{R_3} \right) \right) = \frac{G M_3 m \Delta h}{(R_3 + \Delta h) R_3} \approx m \frac{G M_3}{R_3^2} \Delta h.

Но, как следует из закона всемирного тяотения, ускорение свободного падения на поверхности Земли равно

g=GM3/R32.g = G M_3 / R_3^2.

Таким образом, действительно

ΔEp=mgΔh.\Delta E_p = mg\Delta h.

Кроме закона всемирного тяотения, следует вспомнить законы Кеплер —

  1. каждая планета движется по эллипсу, в одном из фокусов которого находится Солнце;
  2. радиус-вектор планеты за одинаковые промежутки времени описывается равные площади;
  3. квадраты перiodов обращения планет относится как кубы больших полуосей их орбит.

44

Эти законы, полученные при обработке результатов многоолетних наблюдений датского астронама Тихо Браге, справедливы и для спутников. Так, например, для спутника, движущегося вокруг Земли, Земля выступает в роли Солнца, а спутник — в роли планеты.

Рассмотрим теперь несколько конкретных задач.

З адача 1. Покажите, что перiod спутника, обращающийся воркре пленеты в неосредственной близости от ее поверхности, зависит только от средней плотности пленеты. Вычислите перiod такого спутника для нейтронной звезды, считая, что плотность нейтронной звезды порядка плотности атомных ядер Q=1017Q = 10^{17} ке/м3^3.

Сила притяжения спутника к планете сообщает спутнику центростремильноеускорение:

GmMr2=mω2r,G \frac{mM}{r^2} = m\omega^2 r,

где mm — масса спутника, rr — его расстояние до центра планеты, примерное равное радиусу планеты, M=4/3πr3M = 4/3\pi r^3 — масса планеты. Отсюда можно найти угловую скорость ω\omega, а значит, и период TT обращения спутника:

T=2πω=3πGω.T = \frac{2\pi}{\omega} = \sqrt{\frac{3\pi}{G\omega}}.

Для спутника нейтронной звезды T≈1,2⋅10−3T \approx 1,2 \cdot 10^{-3} с.

З адача 2. Определите отношение массы Марса к массе Земли по параметрам орбиты советской автоматической станции «Марс-2»: максимальное удаление от поверхности (в апоцентре) a=25000a = 25000 км, минимальное (в апоцентре) p=1380p = 1380 км, период обращения T=18T = 18 и 00 мин. Радуус Марса rM=3400r_M = 3400 км, радуус Земли r3=6400r_3 = 6400 км.

Из условия задачи известны период TT и длина большой полуоси aM=2rM+a+p=33180a_M = 2r_M + a + p = 33180 км эллипса, по которым движется автоматическая станция «Марс-2». Легко вычислить периоды T3T_3 и TMT_M спутников, врашаюшихся по круговым орбитам

вблизи поверхности Земли и Марса, и сравнить все три периода.

Для спутника Земли можно написать

GmcM3r32=mc4π2T32r3,G \frac{m_c M_3}{r_3^2} = m_c \frac{4\pi^2}{T_3^2} r_3,

откуда

T32(2r3)3=π22GM3.\frac{T_3^2}{(2r_3)^3} = \frac{\pi^2}{2GM_3}.

Для аналогичного спутника Марса

TM2(2rM)3=π22GMM.\frac{T_M^2}{(2r_M)^3} = \frac{\pi^2}{2GM_M}.

По третьему закону Кеплера

TM2T2=(2rM)3aM3.\frac{T_M^2}{T^2} = \frac{(2r_M)^3}{a_M^3}.

Таким образом, из последних трех равенств получаем

4π2T32r3=g,\frac{4\pi^2}{T_3^2} r_3 = g, и T32=4π2r3gT_3^2 = \frac{4\pi^2 r_3}{g}.

Итак, окончательно

MMM3=π2aM32r32gT2=0,11.\frac{M_M}{M_3} = \frac{\pi^2 a_M^3}{2r_3^2 g T^2} = 0,11.

З адача 3. Искусственный спутник Земли запущен с экватора и врачается по круговой орбите в направлении врачия Земли. Найдите отношение радиуса орбиты спутника к радиусу Земли, при котором спутник перiodически проходит над точкой запускаровно через суток.

Так как угловая скорость спутника ωc\omega_c может быть и больше, и меньше угловой скорости Земли ω3\omega_3, возможно два решения задачи.

  1. Пусть ωc>ω3\omega_c > \omega_3, тогда

(ωc−ω3)2T3=2π,(\omega_c - \omega_3) 2T_3 = 2\pi,

ω3=2πT3 и ωc=3πT3\omega_3 = \frac{2\pi}{T_3} \text{ и } \omega_c = \frac{3\pi}{T_3}

(здесь T3T_3 — один сутки — период

45

образения Земли вокруг своей оси).
Согласно второму закону Ньютона,
GmcM3r2=mcωc2r.G \frac{m_c M_3}{r^2} = m_c \omega_c^2 r.

Подставив сюда выражение для ωc\omega_c и зная, что GM3=gr32GM_3 = gr_3^2, найдем отношение радиуса орбиты спутника rr к радиусу Земли r3r_3:

(rr3)3=gT329π2r3≈130,rr3≈5.\left( \frac{r}{r_3} \right)^3 = \frac{g T_3^2}{9 \pi^2 r_3} \approx 130, \frac{r}{r_3} \approx 5.

  1. Пусть ωc<ω3\omega_c < \omega_3, тогда
    (ω3−ωc)2T3=2π,ωc=πT3,(\omega_3 - \omega_c) 2T_3 = 2\pi, \omega_c = \frac{\pi}{T_3},
    (rr3)3=gT32π2r3≈1170,rr3≈10,5.\left( \frac{r}{r_3} \right)^3 = \frac{g T_3^2}{\pi^2 r_3} \approx 1170, \frac{r}{r_3} \approx 10,5.

Задача 4. Спутник, запущенный на крутовую орбиту высотой H=500H = 500 km над поверхностью Земли, тормозится в верхих слоях атмосферы. Уеловое ускорение спутника равно β=−3⋅10−13\beta = -3 \cdot 10^{-13} рад/с². На какой высоте оказывается спутник через месяц?

Обозначим через T0T_0 и r0r_0 начальные периоды обращения и радиуса орбиты спутника, а через TT и rr — период и радиус орбиты спутня месяц после запуска. По третьему закону Кеплера
T02T2=r03r3;\frac{T_0^2}{T^2} = \frac{r_0^3}{r^3};

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

Пусть период увеличился на ΔT\Delta T, а радиус орбиты — на Δr\Delta r. Поскольку уловое ускорение β\beta мало, можно считать, что ΔT≪T\Delta T \ll T и Δr≪r\Delta r \ll r. Из равенства
T02(r0+Δr)3=(T0+ΔT)2r03,T_0^2 (r_0 + \Delta r)^3 = (T_0 + \Delta T)^2 r_0^3,
пренебрега всеми степенями ΔT/T\Delta T / T и Δr/r\Delta r / r, кроме первого, получим
2ΔTT0=3Δrr0.2 \frac{\Delta T}{T_0} = 3 \frac{\Delta r}{r_0}.

Перiod обращения TT связаан с угловой скоростью ω\omega соответствием T=2π/ωT = 2\pi / \omega. Продиференцируем его: T′(ω)=−2π/ω2T'(\omega) = -2\pi / \omega^2. Ввиду малости ΔT\Delta T и Δω\Delta \omega будем считать, что ΔTΔω=T′(ω)\frac{\Delta T}{\Delta \omega} = T'(\omega). Тогда
ΔTT0=−Δωω0=−βtω0.\frac{\Delta T}{T_0} = -\frac{\Delta \omega}{\omega_0} = -\frac{\beta t}{\omega_0}.

Начальную угловую скорость ω0\omega_0 найдем из условия, что на высоте H≪r3H \ll r_3 центростремительное ускорение спутника приблизительно равно gg:
ω02r0=g,ω0=gr0=gr3+H.\omega_0^2 r_0 = g, \omega_0 = \sqrt{\frac{g}{r_0}} = \sqrt{\frac{g}{r_3 + H}}.
Окончательно
Δr=23ΔTT0r0=−23βt(r3+H)r3+Hg≈3 км.\Delta r = \frac{2}{3} \frac{\Delta T}{T_0} r_0 = -\frac{2}{3} \beta t (r^3 + H) \sqrt{\frac{r_3 + H}{g}} \approx 3 \text{ км}.
Итак, через месяц спутник будет находиться на высоте H′=503H' = 503 км.

Упражнения

  1. Искусственный спутник Земли запущен с экватора и врашается по круговой орбите в плоскости экватора в направлении осевого врашения Земли. Радиус орбиты спутника в α=3\alpha = 3 раза больше радиуса Земли r3=6400r_3 = 6400 км. Через какое время спутник в первой раз пройдёт над точкой запуска?
  2. Вычислите первую космическую скорость при старте с поверхности Юпитера, если известно, что спутник Юпитера Ганимед врашается по почти круговой орбите радиусом r=1⋅106r = 1 \cdot 10^6 км с периодом T=7,15T = 7,15 земных суток. Радиус Юпитера r0=70⋅103r_0 = 70 \cdot 10^3 км.
  3. Известно, что в настоящие время Луна удаляется от Земли со скоростью v⃗(∣v⃗∣=3,3\vec{v} (|\vec{v}| = 3,3 см/год). Найдите угловое ускорение Луны. Среднее расстояние от Луны до Земли R=3,84⋅105R = 3,84 \cdot 10^5 км, угловая скорость врашения Луны вокруг Земли ω=2,56⋅10−6\omega = 2,56 \cdot 10^{-6} рад/с.