Уравнение Vis-Viva

Уравнение для моделирования движения орбитальных тел

В астродинамике уравнение vis-viva , также называемое законом инвариантности орбитальной энергии или формулой Бургаса [1] [ нужен лучший источник ] , является одним из уравнений, моделирующих движение орбитальных тел . Это прямой результат принципа сохранения механической энергии , который применяется, когда единственной силой, действующей на объект, является его собственный вес, который является гравитационной силой, определяемой произведением массы объекта и напряженности окружающего гравитационного поля .

Vis viva (лат. «живая сила») — термин из истории механики, сохранившийся только в этом контексте. Он представляет собой принцип, согласно которому разница между общей работой ускоряющих сил системыи работой замедляющих сил равна половине vis viva , накопленной или потерянной в системе за время выполнения работы.

Уравнение

Для любой кеплеровской орбиты ( эллиптической , параболической , гиперболической или радиальной ) уравнение vis-viva [2] выглядит следующим образом: [3] где: в 2 = Г М ( 2 г 1 а ) {\displaystyle v^{2}=GM\left({2 \over r}-{1 \over a}\right)}

Произведение GM также можно выразить как стандартный гравитационный параметр, используя греческую букву μ .

Вывод для эллиптических орбит (0 ≤ эксцентриситет < 1)

В уравнении vis-viva масса m вращающегося тела (например, космического корабля) считается пренебрежимо малой по сравнению с массой M центрального тела (например, Земли). Центральное тело и вращающееся тело также часто называют первичным и частицей соответственно. В конкретных случаях эллиптической или круговой орбиты уравнение vis-viva может быть легко выведено из законов сохранения энергии и импульса.

Удельная полная энергия постоянна по всей орбите. Таким образом, используя индексы a и p для обозначения апогея (апогея) и перигея (перигея) соответственно, ε = в а 2 2 Г М г а = в п 2 2 Г М г п {\displaystyle \varepsilon ={\frac {v_{a}^{2}}{2}}-{\frac {GM}{r_{a}}}={\frac {v_{p}^{2}}{2}}-{\frac {GM}{r_{p}}}}

Перестановка, в а 2 2 в п 2 2 = Г М г а Г М г п {\displaystyle {\frac {v_{a}^{2}}{2}}-{\frac {v_{p}^{2}}{2}}={\frac {GM}{r_{a}}}-{\frac {GM}{r_{p}}}}

Вспоминая, что для эллиптической орбиты (и, следовательно, также круговой орбиты) векторы скорости и радиуса перпендикулярны в апоцентре и перицентре, сохранение момента импульса требует определенного момента импульса , таким образом : час = г п в п = г а в а = constant {\displaystyle h=r_{p}v_{p}=r_{a}v_{a}={\text{constant}}} v p = r a r p v a {\displaystyle v_{p}={\frac {r_{a}}{r_{p}}}v_{a}} 1 2 ( 1 r a 2 r p 2 ) v a 2 = G M r a G M r p {\displaystyle {\frac {1}{2}}\left(1-{\frac {r_{a}^{2}}{r_{p}^{2}}}\right)v_{a}^{2}={\frac {GM}{r_{a}}}-{\frac {GM}{r_{p}}}} 1 2 ( r p 2 r a 2 r p 2 ) v a 2 = G M r a G M r p {\displaystyle {\frac {1}{2}}\left({\frac {r_{p}^{2}-r_{a}^{2}}{r_{p}^{2}}}\right)v_{a}^{2}={\frac {GM}{r_{a}}}-{\frac {GM}{r_{p}}}}

Выделяя кинетическую энергию в апоцентре и упрощая, 1 2 v a 2 = ( G M r a G M r p ) r p 2 r p 2 r a 2 1 2 v a 2 = G M ( r p r a r a r p ) r p 2 r p 2 r a 2 1 2 v a 2 = G M r p r a ( r p + r a ) {\displaystyle {\begin{aligned}{\frac {1}{2}}v_{a}^{2}&=\left({\frac {GM}{r_{a}}}-{\frac {GM}{r_{p}}}\right)\cdot {\frac {r_{p}^{2}}{r_{p}^{2}-r_{a}^{2}}}\\{\frac {1}{2}}v_{a}^{2}&=GM\left({\frac {r_{p}-r_{a}}{r_{a}r_{p}}}\right){\frac {r_{p}^{2}}{r_{p}^{2}-r_{a}^{2}}}\\{\frac {1}{2}}v_{a}^{2}&=GM{\frac {r_{p}}{r_{a}(r_{p}+r_{a})}}\end{aligned}}}

Из геометрии эллипса, где a — длина большой полуоси. Таким образом, 2 a = r p + r a {\displaystyle 2a=r_{p}+r_{a}} 1 2 v a 2 = G M 2 a r a r a ( 2 a ) = G M ( 1 r a 1 2 a ) = G M r a G M 2 a {\displaystyle {\frac {1}{2}}v_{a}^{2}=GM{\frac {2a-r_{a}}{r_{a}(2a)}}=GM\left({\frac {1}{r_{a}}}-{\frac {1}{2a}}\right)={\frac {GM}{r_{a}}}-{\frac {GM}{2a}}}

Подставим это в наше исходное выражение для удельной орбитальной энергии: ε = v 2 2 G M r = v p 2 2 G M r p = v a 2 2 G M r a = G M 2 a {\displaystyle \varepsilon ={\frac {v^{2}}{2}}-{\frac {GM}{r}}={\frac {v_{p}^{2}}{2}}-{\frac {GM}{r_{p}}}={\frac {v_{a}^{2}}{2}}-{\frac {GM}{r_{a}}}=-{\frac {GM}{2a}}}

Таким образом, и уравнение vis-viva может быть записано или ε = G M 2 a {\displaystyle \varepsilon =-{\frac {GM}{2a}}} v 2 2 G M r = G M 2 a {\displaystyle {\frac {v^{2}}{2}}-{\frac {GM}{r}}=-{\frac {GM}{2a}}} v 2 = G M ( 2 r 1 a ) {\displaystyle v^{2}=GM\left({\frac {2}{r}}-{\frac {1}{a}}\right)}

Следовательно, сохраняющийся угловой момент L = mh можно вывести с помощью и , где aбольшая полуось , а bмалая полуось эллиптической орбиты, следующим образом: и поочередно, r a + r p = 2 a {\displaystyle r_{a}+r_{p}=2a} r a r p = b 2 {\displaystyle r_{a}r_{p}=b^{2}} v a 2 = G M ( 2 r a 1 a ) = G M a ( 2 a r a r a ) = G M a ( r p r a ) = G M a ( b r a ) 2 {\displaystyle v_{a}^{2}=GM\left({\frac {2}{r_{a}}}-{\frac {1}{a}}\right)={\frac {GM}{a}}\left({\frac {2a-r_{a}}{r_{a}}}\right)={\frac {GM}{a}}\left({\frac {r_{p}}{r_{a}}}\right)={\frac {GM}{a}}\left({\frac {b}{r_{a}}}\right)^{2}} v p 2 = G M ( 2 r p 1 a ) = G M a ( 2 a r p r p ) = G M a ( r a r p ) = G M a ( b r p ) 2 {\displaystyle v_{p}^{2}=GM\left({\frac {2}{r_{p}}}-{\frac {1}{a}}\right)={\frac {GM}{a}}\left({\frac {2a-r_{p}}{r_{p}}}\right)={\frac {GM}{a}}\left({\frac {r_{a}}{r_{p}}}\right)={\frac {GM}{a}}\left({\frac {b}{r_{p}}}\right)^{2}}

Следовательно, удельный момент импульса , и h = r p v p = r a v a = b G M a {\displaystyle h=r_{p}v_{p}=r_{a}v_{a}=b{\sqrt {\frac {GM}{a}}}}

Полный угловой момент L = m h = m b G M a {\displaystyle L=mh=mb{\sqrt {\frac {GM}{a}}}}

Практические применения

Учитывая общую массу и скаляры r и v в одной точке орбиты, можно вычислить:

  • r и v в любой другой точке орбиты; [примечания 1] и
  • удельная орбитальная энергия , позволяющая классифицировать объект, вращающийся по орбите вокруг более крупного объекта, как не имеющий достаточно энергии, чтобы оставаться на орбите, следовательно, являющийся « суборбитальным » (например, баллистическая ракета), имеющий достаточно энергии, чтобы быть «орбитальным», но без возможности завершить полный оборот по орбите в любом случае, поскольку он в конечном итоге сталкивается с другим телом, или имеющий достаточно энергии, чтобы прийти из бесконечности и/или уйти в нее (например, метеор). ε {\displaystyle \varepsilon \,\!}

Формулу для скорости убегания можно получить из уравнения Vis-viva, взяв предел как приближающийся : a {\displaystyle a} {\displaystyle \infty } v e 2 = G M ( 2 r 0 ) v e = 2 G M r {\displaystyle v_{e}^{2}=GM\left({\frac {2}{r}}-0\right)\rightarrow v_{e}={\sqrt {\frac {2GM}{r}}}}

Примечания

  1. ^ Для задачи трех тел едва ли существует сопоставимое уравнение vis-viva: сохранение энергии уменьшает большее число степеней свободы всего на одну.

Ссылки

  1. ^ Иванов, Стефан: XXV Национальная олимпиада по астрономии, Бургас, 06-08.05.2022, Полезни формулы и справочные данные (Полезные формулы и справочные данные)
  2. ^ Том Логсдон (1998). Орбитальная механика: теория и приложения. John Wiley & Sons. ISBN 978-0-471-14636-0.
  3. ^ Лиссауэр, Джек Дж.; де Патер, Имке (2019). Фундаментальные планетарные науки: физика, химия и обитаемость . Нью-Йорк, Нью-Йорк, США: Cambridge University Press. стр. 29–31. ISBN 9781108411981.
Retrieved from "https://en.wikipedia.org/w/index.php?title=Vis-viva_equation&oldid=1217771597"