Векторная величина в небесной механике
В небесной механике удельный относительный момент импульса (часто обозначаемый или ) тела — это момент импульса этого тела, деленный на его массу. [1] В случае двух вращающихся тел это векторное произведение их относительного положения и относительного линейного импульса , деленное на массу рассматриваемого тела. час → {\displaystyle {\vec {h}}} час {\displaystyle \mathbf {ч} }
Удельный относительный угловой момент играет ключевую роль в анализе задачи двух тел , поскольку он остается постоянным для заданной орбиты в идеальных условиях. « Удельный » в этом контексте означает угловой момент на единицу массы. Единицей СИ для удельного относительного углового момента является квадратный метр в секунду.
Определение Удельный относительный момент импульса определяется как векторное произведение вектора относительного положения и вектора относительной скорости . г {\displaystyle \mathbf {r} } в {\displaystyle \mathbf {v} } час = г × в = Л м {\displaystyle \mathbf {h} =\mathbf {r} \times \mathbf {v} ={\frac {\mathbf {L} }{m}}}
где — вектор момента импульса, определяемый как . Л {\displaystyle \mathbf {L} } г × м в {\displaystyle \mathbf {r} \times m\mathbf {v} }
Вектор всегда перпендикулярен мгновенной оскулирующей орбитальной плоскости , которая совпадает с мгновенной возмущенной орбитой . Он не обязательно перпендикулярен средней орбитальной плоскости во времени. час {\displaystyle \mathbf {ч} }
Доказательство постоянства в случае двух тел Вектор расстояния , вектор скорости , истинная аномалия и угол траектории полета на орбите вокруг . Также изображены наиболее важные меры эллипса ( среди которых, обратите внимание, что истинная аномалия обозначена как ). г {\displaystyle \mathbf {r} } в {\displaystyle \mathbf {v} } θ {\displaystyle \тета} ϕ {\displaystyle \фи} м 2 {\displaystyle m_{2}} м 1 {\displaystyle m_{1}} θ {\displaystyle \тета} ν {\displaystyle \nu} При определенных условиях можно доказать, что удельный момент импульса постоянен. Условия этого доказательства включают:
Масса одного объекта намного больше массы другого. ( ) м 1 ≫ м 2 {\displaystyle m_{1}\gg m_{2}} Система координат инерциальная . Каждый объект можно рассматривать как сферически симметричную точечную массу . На систему не действуют никакие другие силы, кроме силы тяготения, которая соединяет два тела.
Доказательство Доказательство начинается с уравнения движения двух тел , выведенного из закона всемирного тяготения Ньютона :
г ¨ + Г м 1 г 2 г г = 0 {\displaystyle {\ddot {\mathbf {r} }}+{\frac {Gm_{1}}{r^{2}}}{\frac {\mathbf {r} }{r}}=0}
где:
г {\displaystyle \mathbf {r} } — радиус-вектор от до со скалярной величиной . м 1 {\displaystyle m_{1}} м 2 {\displaystyle m_{2}} г {\displaystyle r} г ¨ {\displaystyle {\ddot {\mathbf {r} }}} является второй производной по времени от . ( ускорение ) г {\displaystyle \mathbf {r} } Г {\displaystyle G} Гравитационная постоянная .Перекрестное произведение вектора положения и уравнения движения равно:
г × г ¨ + г × Г м 1 г 2 г г = 0 {\displaystyle \mathbf {r} \times {\ddot {\mathbf {r} }}+\mathbf {r} \times {\frac {Gm_{1}}{r^{2}}}{\frac {\mathbf {r} }{r}}=0}
Поскольку второй член исчезает: г × г = 0 {\displaystyle \mathbf {r} \times \mathbf {r} =0}
г × г ¨ = 0 {\displaystyle \mathbf {r} \times {\ddot {\mathbf {r} }}=0}
Также можно сделать вывод, что: г г т ( г × г ˙ ) = г ˙ × г ˙ + г × г ¨ = г × г ¨ {\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\left(\mathbf {r} \times {\dot {\mathbf {r} }}\right)={\dot {\mathbf {r} }}\times {\dot {\mathbf {r} }}+\mathbf {r} \times {\ddot {\mathbf {r} }}=\mathbf {r} \times {\ddot {\mathbf {r} }}}
Объединение этих двух уравнений дает: d d t ( r × r ˙ ) = 0 {\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\left(\mathbf {r} \times {\dot {\mathbf {r} }}\right)=0}
Поскольку производная по времени равна нулю, величина постоянна. Используя вектор скорости вместо скорости изменения положения, а для удельного момента импульса: является постоянной. r × r ˙ {\displaystyle \mathbf {r} \times {\dot {\mathbf {r} }}} v {\displaystyle \mathbf {v} } h {\displaystyle \mathbf {h} } h = r × v {\displaystyle \mathbf {h} =\mathbf {r} \times \mathbf {v} }
Это отличается от обычного построения импульса, поскольку не включает массу рассматриваемого объекта. r × p {\displaystyle \mathbf {r} \times \mathbf {p} }
Законы движения планет Кеплера Законы движения планет Кеплера можно доказать почти напрямую с помощью приведенных выше соотношений.
Первый закон Доказательство снова начинается с уравнения задачи двух тел. На этот раз векторное произведение умножается на удельный относительный угловой момент r ¨ × h = − μ r 2 r r × h {\displaystyle {\ddot {\mathbf {r} }}\times \mathbf {h} =-{\frac {\mu }{r^{2}}}{\frac {\mathbf {r} }{r}}\times \mathbf {h} }
Левая часть равна производной, поскольку момент импульса постоянен. d d t ( r ˙ × h ) {\textstyle {\frac {\mathrm {d} }{\mathrm {d} t}}\left({\dot {\mathbf {r} }}\times \mathbf {h} \right)}
После некоторых шагов (включая использование тройного векторного произведения и определение скаляра как радиальной скорости , а не нормы вектора ) правая часть становится: r ˙ {\displaystyle {\dot {r}}} r ˙ {\displaystyle {\dot {\mathbf {r} }}} − μ r 3 ( r × h ) = − μ r 3 ( ( r ⋅ v ) r − r 2 v ) = − ( μ r 2 r ˙ r − μ r v ) = μ d d t ( r r ) {\displaystyle -{\frac {\mu }{r^{3}}}\left(\mathbf {r} \times \mathbf {h} \right)=-{\frac {\mu }{r^{3}}}\left(\left(\mathbf {r} \cdot \mathbf {v} \right)\mathbf {r} -r^{2}\mathbf {v} \right)=-\left({\frac {\mu }{r^{2}}}{\dot {r}}\mathbf {r} -{\frac {\mu }{r}}\mathbf {v} \right)=\mu {\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\mathbf {r} }{r}}\right)}
Приравнивая эти два выражения и интегрируя по времени, получаем (с постоянной интегрирования ) C {\displaystyle \mathbf {C} } r ˙ × h = μ r r + C {\displaystyle {\dot {\mathbf {r} }}\times \mathbf {h} =\mu {\frac {\mathbf {r} }{r}}+\mathbf {C} }
Теперь это уравнение умножается ( скалярное произведение ) на и переставляется r {\displaystyle \mathbf {r} } r ⋅ ( r ˙ × h ) = r ⋅ ( μ r r + C ) ⇒ ( r × r ˙ ) ⋅ h = μ r + r C cos θ ⇒ h 2 = μ r + r C cos θ {\displaystyle {\begin{aligned}\mathbf {r} \cdot \left({\dot {\mathbf {r} }}\times \mathbf {h} \right)&=\mathbf {r} \cdot \left(\mu {\frac {\mathbf {r} }{r}}+\mathbf {C} \right)\\\Rightarrow \left(\mathbf {r} \times {\dot {\mathbf {r} }}\right)\cdot \mathbf {h} &=\mu r+rC\cos \theta \\\Rightarrow h^{2}&=\mu r+rC\cos \theta \end{aligned}}}
В итоге получаем уравнение орбиты [1] r = h 2 μ 1 + C μ cos θ {\displaystyle r={\frac {\frac {h^{2}}{\mu }}{1+{\frac {C}{\mu }}\cos \theta }}}
которое является уравнением конического сечения в полярных координатах с полушириной и эксцентриситетом . p = h 2 μ {\textstyle p={\frac {h^{2}}{\mu }}} e = C μ {\textstyle e={\frac {C}{\mu }}}
Второй закон Второй закон немедленно следует из второго из трех уравнений для расчета абсолютного значения удельного относительного момента импульса. [1]
Если связать эту форму уравнения с соотношением для площади сектора с бесконечно малым углом (треугольник с одной очень малой стороной), то уравнение d t = r 2 h d θ {\textstyle \mathrm {d} t={\frac {r^{2}}{h}}\,\mathrm {d} \theta } d A = r 2 2 d θ {\textstyle \mathrm {d} A={\frac {r^{2}}{2}}\,\mathrm {d} \theta } d θ {\displaystyle \mathrm {d} \theta } d t = 2 h d A {\displaystyle \mathrm {d} t={\frac {2}{h}}\,\mathrm {d} A}
Третий закон Третий закон Кеплера является прямым следствием второго закона. Интегрирование по одному обороту дает орбитальный период [1] T = 2 π a b h {\displaystyle T={\frac {2\pi ab}{h}}}
для площади эллипса. Заменив малую полуось на и удельный относительный момент импульса на, получим π a b {\displaystyle \pi ab} b = a p {\displaystyle b={\sqrt {ap}}} h = μ p {\displaystyle h={\sqrt {\mu p}}} T = 2 π a 3 μ {\displaystyle T=2\pi {\sqrt {\frac {a^{3}}{\mu }}}}
Таким образом, между большой полуосью и периодом обращения спутника существует связь, которую можно свести к константе центрального тела.
Смотрите также
Ссылки
^ abcd Валладо, Дэвид А. (2001). Основы астродинамики и приложения (2-е изд.). Дордрехт: Kluwer Academic Publishers. стр. 20–30. ISBN 0-7923-6903-3 .