Разработайте задачу по небесной механике: определите параметры свободной возвратной траектории к Луне (free-return) для пилотируемой миссии и обсудите преимущества и ограничения такого профиля полёта
Задача (формулировка). Для пилотируемой миссии из круговой околоземной орбиты (LEO) высотой hLEO=200 kmh_{\rm LEO}=200\ \text{km}hLEO=200km над поверхностью Земли определить параметры классической свободной возвратной (free‑return) траектории к Луне в приближении «патч‑коник»: необходимые Δv\Delta vΔv для TLI, скорость на удалённости Луны (гиперболический избыток) v∞v_\inftyv∞, выбор перицентра облёта Луны \(r_{p,\Moon}\) (например, высота hp=100 kmh_{\rm p}=100\ \text{km}hp=100km над поверхностью Луны) и получаемый угол отклонения при пролёте δ\deltaδ. Оцените время перелёта и скорость при возвращении к околоземной орбите. Обсудите преимущества и ограничения free‑return профиля. Дано (числа для примера) - гравитационный параметр Земли μE=398600 km3/s2\mu_E=398600\ \text{km}^3/\text{s}^2μE=398600km3/s2, - радиус Земли RE=6378 kmR_E=6378\ \text{km}RE=6378km, - среднее расстояние Луны от центра Земли \(r_{\Moon}=384400\ \text{km}\), - гравитационный параметр Луны \(\mu_{\Moon}=4902.8\ \text{km}^3/\text{s}^2\), - радиус Луны \(R_{\Moon}=1737\ \text{km}\). Решение (патч‑коник, «эллипс с апогеем в точке Луны»). 1) Параметры исходной орбиты и трансфера: - радиус LEO rLEO=RE+hLEO=6578 kmr_{\rm LEO}=R_E+h_{\rm LEO}=6578\ \text{km}rLEO=RE+hLEO=6578km. - полуось трансферной эллиптической орбиты \(\displaystyle a=\frac{r_{\rm LEO}+r_{\Moon}}{2}\). 2) Скорости и Δv\Delta vΔv TLI: - орбитальная скорость в LEO: vcirc=μErLEO\displaystyle v_{\rm circ}=\sqrt{\frac{\mu_E}{r_{\rm LEO}}}vcirc=rLEOμE. - скорость на перигее трансферной эллипсы (в точке инжекции): vtrans(rLEO)=μE (2rLEO−1a)\displaystyle v_{\rm trans}(r_{\rm LEO})=\sqrt{\mu_E\!\left(\frac{2}{r_{\rm LEO}}-\frac{1}{a}\right)}vtrans(rLEO)=μE(rLEO2−a1). - требуемый импульс для TLI: ΔvTLI=vtrans(rLEO)−vcirc\displaystyle \Delta v_{\rm TLI}=v_{\rm trans}(r_{\rm LEO})-v_{\rm circ}ΔvTLI=vtrans(rLEO)−vcirc. Подстановка чисел: a=6578+3844002=195489 km,
a=\frac{6578+384400}{2}=195489\ \text{km}, a=26578+384400=195489km,vcirc=3986006578≈7.789 km/s,
v_{\rm circ}=\sqrt{\frac{398600}{6578}}\approx 7.789\ \text{km/s}, vcirc=6578398600≈7.789km/s,vtrans(rLEO)=398600(26578−1195489)≈10.917 km/s,
v_{\rm trans}(r_{\rm LEO})=\sqrt{398600\left(\frac{2}{6578}-\frac{1}{195489}\right)}\approx 10.917\ \text{km/s}, vtrans(rLEO)=398600(65782−1954891)≈10.917km/s,ΔvTLI≈10.917−7.789=3.13 km/s.
\Delta v_{\rm TLI}\approx 10.917-7.789=3.13\ \text{km/s}. ΔvTLI≈10.917−7.789=3.13km/s. 3) Скорость в точке апогея (на расстоянии Луны) и v∞v_\inftyv∞ относительно Луны: - скорость на апогее трансфера: \(\displaystyle v_{\rm trans}(r_{\Moon})=\sqrt{\mu_E\!\left(\frac{2}{r_{\Moon}}-\frac{1}{a}\right)}\). - орбитальная скорость Луны вокруг Земли: \(\displaystyle v_{\Moon}=\sqrt{\frac{\mu_E}{r_{\Moon}}}\). - гиперболический избыток относительно Луны: \(\displaystyle v_\infty=\big|v_{\rm trans}(r_{\Moon})-v_{\Moon}\big|\). Численно: \[ v_{\rm trans}(r_{\Moon})\approx 0.188\ \text{km/s},\quad v_{\Moon}\approx 1.018\ \text{km/s}, \] v∞≈0.831 km/s.
v_\infty\approx 0.831\ \text{km/s}. v∞≈0.831km/s. 4) Параметры облёта Луны и угол отклонения δ\deltaδ. Для заданного перицентра облёта (радиус центра Луны до перицентра) \(\displaystyle r_{p,\Moon}=R_{\Moon}+h_{\rm p}\) (при hp=100 kmh_{\rm p}=100\ \text{km}hp=100km получаем \(r_{p,\Moon}=1837\ \text{km}\)) можно оценить параметр столкновения bbb и угол отклонения: - связь bbb и rpr_prp: \(\displaystyle b=r_p\sqrt{1+\frac{2\mu_{\Moon}}{r_p v_\infty^2}}\). - угол отклонения асимптоты (гравитационный поворот): \(\displaystyle \delta=2\arctan\!\left(\frac{\mu_{\Moon}}{b\,v_\infty^2}\right)\). Численно (для hp=100 kmh_{\rm p}=100\ \text{km}hp=100km): \[ \frac{2\mu_{\Moon}}{r_p v_\infty^2}\approx 7.74,\quad b\approx 5434\ \text{km}, \] δ≈2arctan(4902.85434⋅0.8312)≈105∘.
\delta\approx 2\arctan\left(\frac{4902.8}{5434\cdot 0.831^2}\right)\approx 105^\circ. δ≈2arctan(5434⋅0.83124902.8)≈105∘. Такой поворот (>90°) достаточен, чтобы переориентировать траекторию назад к Земле — получается классический «circumlunar free‑return» (космический аппарат проходит за Луной и возвращается к Земле без дополнительного серьёзного Δv\Delta vΔv). 5) Время перелёта (полупериод трансферной эллипсы): TOF=πa3μE≈π(195489)3398600≈4.98 суток.
{\rm TOF}=\pi\sqrt{\frac{a^3}{\mu_E}}\approx \pi\sqrt{\frac{(195489)^3}{398600}}\approx 4.98\ \text{суток}. TOF=πμEa3≈π398600(195489)3≈4.98суток.
(Это значение для чистого «Hohmann»‑типа трансфера; на практике траектории короче — ~3 сут — возможны за счёт более энергичных инжекций.) 6) Скорость при возвращении к околоземной высоте Возвратная симметрия даёт приближённую скорость вблизи Земли равную скорости трансфера в перигее, т.е. ~10.9 km/s10.9\ \text{km/s}10.9km/s на высоте LEO (приближённо скорость при входе в атмосферу порядка 11 км/s). Краткое обсуждение преимуществ free‑return профиля - Отказоустойчивость: при неработающем двигателе после TLI аппарат сам вернётся к Земле без крупных манёвров. - Экономия топлива: минимальные дополнительные Δv\Delta vΔv для гарантированного возвращения. - Простота навигации: траектория геометрически надёжна и опробована (пример — миссии «Аполлон»). Ограничения и недостатки - Меньшая гибкость: жёсткие требования к времени инжекции и точности, ограниченные варианты облёта/посадки. - Длительность перелёта: «энергетически экономный» трансфер может быть длиннее по времени; при необходимости быстрее попасть к Луне нужен больший Δv\Delta vΔv. - Ограничения по параметрам облёта: чтобы обеспечить нужный угол поворота, требуется достаточно малый перицентр облёта (близко к Луне), что накладывает риски (атмосфера/твердость рельефа не актуальны для Луны, но навигационная точность критична). - Невозможность выполнить мягкую посадку на Луну в рамках простой free‑return без дополнительных манёвров (free‑return — про облёт, не посадку). Итог (рекомендация). Для стандартной пилотируемой миссии из LEO 200 km циркумлунарный free‑return с TLI Δv≈3.1 km/s\Delta v\approx 3.1\ \text{km/s}Δv≈3.1km/s, гиперболическим v∞≈0.83 km/sv_\infty\approx 0.83\ \text{km/s}v∞≈0.83km/s и перицентром облёта hp∼100 kmh_{\rm p}\sim 100\ \text{km}hp∼100km даёт угол отклонения ∼105∘\sim105^\circ∼105∘ и безопасный самовозврат. При проектировании следует учитывать требуемую точность инжекции, желаемое время полёта и требования к аварийным процедурам; при необходимости уменьшить время полёта — увеличить Δv\Delta vΔv и отказаться от чистого free‑return.
Дано (числа для примера)
- гравитационный параметр Земли μE=398600 km3/s2\mu_E=398600\ \text{km}^3/\text{s}^2μE =398600 km3/s2,
- радиус Земли RE=6378 kmR_E=6378\ \text{km}RE =6378 km,
- среднее расстояние Луны от центра Земли \(r_{\Moon}=384400\ \text{km}\),
- гравитационный параметр Луны \(\mu_{\Moon}=4902.8\ \text{km}^3/\text{s}^2\),
- радиус Луны \(R_{\Moon}=1737\ \text{km}\).
Решение (патч‑коник, «эллипс с апогеем в точке Луны»).
1) Параметры исходной орбиты и трансфера:
- радиус LEO rLEO=RE+hLEO=6578 kmr_{\rm LEO}=R_E+h_{\rm LEO}=6578\ \text{km}rLEO =RE +hLEO =6578 km.
- полуось трансферной эллиптической орбиты \(\displaystyle a=\frac{r_{\rm LEO}+r_{\Moon}}{2}\).
2) Скорости и Δv\Delta vΔv TLI:
- орбитальная скорость в LEO: vcirc=μErLEO\displaystyle v_{\rm circ}=\sqrt{\frac{\mu_E}{r_{\rm LEO}}}vcirc =rLEO μE .
- скорость на перигее трансферной эллипсы (в точке инжекции): vtrans(rLEO)=μE (2rLEO−1a)\displaystyle v_{\rm trans}(r_{\rm LEO})=\sqrt{\mu_E\!\left(\frac{2}{r_{\rm LEO}}-\frac{1}{a}\right)}vtrans (rLEO )=μE (rLEO 2 −a1 ) .
- требуемый импульс для TLI: ΔvTLI=vtrans(rLEO)−vcirc\displaystyle \Delta v_{\rm TLI}=v_{\rm trans}(r_{\rm LEO})-v_{\rm circ}ΔvTLI =vtrans (rLEO )−vcirc .
Подстановка чисел:
a=6578+3844002=195489 km, a=\frac{6578+384400}{2}=195489\ \text{km},
a=26578+384400 =195489 km, vcirc=3986006578≈7.789 km/s, v_{\rm circ}=\sqrt{\frac{398600}{6578}}\approx 7.789\ \text{km/s},
vcirc =6578398600 ≈7.789 km/s, vtrans(rLEO)=398600(26578−1195489)≈10.917 km/s, v_{\rm trans}(r_{\rm LEO})=\sqrt{398600\left(\frac{2}{6578}-\frac{1}{195489}\right)}\approx 10.917\ \text{km/s},
vtrans (rLEO )=398600(65782 −1954891 ) ≈10.917 km/s, ΔvTLI≈10.917−7.789=3.13 km/s. \Delta v_{\rm TLI}\approx 10.917-7.789=3.13\ \text{km/s}.
ΔvTLI ≈10.917−7.789=3.13 km/s.
3) Скорость в точке апогея (на расстоянии Луны) и v∞v_\inftyv∞ относительно Луны:
- скорость на апогее трансфера: \(\displaystyle v_{\rm trans}(r_{\Moon})=\sqrt{\mu_E\!\left(\frac{2}{r_{\Moon}}-\frac{1}{a}\right)}\).
- орбитальная скорость Луны вокруг Земли: \(\displaystyle v_{\Moon}=\sqrt{\frac{\mu_E}{r_{\Moon}}}\).
- гиперболический избыток относительно Луны: \(\displaystyle v_\infty=\big|v_{\rm trans}(r_{\Moon})-v_{\Moon}\big|\).
Численно:
\[
v_{\rm trans}(r_{\Moon})\approx 0.188\ \text{km/s},\quad v_{\Moon}\approx 1.018\ \text{km/s},
\]
v∞≈0.831 km/s. v_\infty\approx 0.831\ \text{km/s}.
v∞ ≈0.831 km/s.
4) Параметры облёта Луны и угол отклонения δ\deltaδ.
Для заданного перицентра облёта (радиус центра Луны до перицентра) \(\displaystyle r_{p,\Moon}=R_{\Moon}+h_{\rm p}\) (при hp=100 kmh_{\rm p}=100\ \text{km}hp =100 km получаем \(r_{p,\Moon}=1837\ \text{km}\)) можно оценить параметр столкновения bbb и угол отклонения:
- связь bbb и rpr_prp : \(\displaystyle b=r_p\sqrt{1+\frac{2\mu_{\Moon}}{r_p v_\infty^2}}\).
- угол отклонения асимптоты (гравитационный поворот): \(\displaystyle \delta=2\arctan\!\left(\frac{\mu_{\Moon}}{b\,v_\infty^2}\right)\).
Численно (для hp=100 kmh_{\rm p}=100\ \text{km}hp =100 km):
\[
\frac{2\mu_{\Moon}}{r_p v_\infty^2}\approx 7.74,\quad b\approx 5434\ \text{km},
\]
δ≈2arctan(4902.85434⋅0.8312)≈105∘. \delta\approx 2\arctan\left(\frac{4902.8}{5434\cdot 0.831^2}\right)\approx 105^\circ.
δ≈2arctan(5434⋅0.83124902.8 )≈105∘.
Такой поворот (>90°) достаточен, чтобы переориентировать траекторию назад к Земле — получается классический «circumlunar free‑return» (космический аппарат проходит за Луной и возвращается к Земле без дополнительного серьёзного Δv\Delta vΔv).
5) Время перелёта (полупериод трансферной эллипсы):
TOF=πa3μE≈π(195489)3398600≈4.98 суток. {\rm TOF}=\pi\sqrt{\frac{a^3}{\mu_E}}\approx \pi\sqrt{\frac{(195489)^3}{398600}}\approx 4.98\ \text{суток}.
TOF=πμE a3 ≈π398600(195489)3 ≈4.98 суток. (Это значение для чистого «Hohmann»‑типа трансфера; на практике траектории короче — ~3 сут — возможны за счёт более энергичных инжекций.)
6) Скорость при возвращении к околоземной высоте
Возвратная симметрия даёт приближённую скорость вблизи Земли равную скорости трансфера в перигее, т.е. ~10.9 km/s10.9\ \text{km/s}10.9 km/s на высоте LEO (приближённо скорость при входе в атмосферу порядка 11 км/s).
Краткое обсуждение преимуществ free‑return профиля
- Отказоустойчивость: при неработающем двигателе после TLI аппарат сам вернётся к Земле без крупных манёвров.
- Экономия топлива: минимальные дополнительные Δv\Delta vΔv для гарантированного возвращения.
- Простота навигации: траектория геометрически надёжна и опробована (пример — миссии «Аполлон»).
Ограничения и недостатки
- Меньшая гибкость: жёсткие требования к времени инжекции и точности, ограниченные варианты облёта/посадки.
- Длительность перелёта: «энергетически экономный» трансфер может быть длиннее по времени; при необходимости быстрее попасть к Луне нужен больший Δv\Delta vΔv.
- Ограничения по параметрам облёта: чтобы обеспечить нужный угол поворота, требуется достаточно малый перицентр облёта (близко к Луне), что накладывает риски (атмосфера/твердость рельефа не актуальны для Луны, но навигационная точность критична).
- Невозможность выполнить мягкую посадку на Луну в рамках простой free‑return без дополнительных манёвров (free‑return — про облёт, не посадку).
Итог (рекомендация). Для стандартной пилотируемой миссии из LEO 200 km циркумлунарный free‑return с TLI Δv≈3.1 km/s\Delta v\approx 3.1\ \text{km/s}Δv≈3.1 km/s, гиперболическим v∞≈0.83 km/sv_\infty\approx 0.83\ \text{km/s}v∞ ≈0.83 km/s и перицентром облёта hp∼100 kmh_{\rm p}\sim 100\ \text{km}hp ∼100 km даёт угол отклонения ∼105∘\sim105^\circ∼105∘ и безопасный самовозврат. При проектировании следует учитывать требуемую точность инжекции, желаемое время полёта и требования к аварийным процедурам; при необходимости уменьшить время полёта — увеличить Δv\Delta vΔv и отказаться от чистого free‑return.