14.4: Lagrange's Planetary Equations
- Page ID
- 6878
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
We now go to Equation 14.2.8 to obtain Lagrange’s Planetary Equations, which will enable us to calculate the rates of change of the orbital elements if we know the form of the perturbing function:
\[ \begin{align} \dot{a} &= - \frac{2a^2}{GMm} \frac{\partial R}{\partial T} , \label{14.4.1} \\[5pt] \dot{e} &= - \frac{a(1-e^2)}{GMme} \frac{\partial R}{\partial T} , \label{14.4.2} \\[5pt] i &= - \frac{1}{\sqrt{GMm^2 a (1-e^2) \sin i}} \frac{\partial R}{\partial \Omega} - \frac{1}{me} \sqrt{\frac{1 - e^2}{GMa}} \frac{\partial R}{\partial ω} , \label{14.4.3} \\[5pt] \dot{ω} &= \frac{1}{me} \sqrt{\frac{1 - e^2}{GMa}} \frac{\partial R}{\partial e} - \frac{1}{\sqrt{GMm^2 a (1 - e^2)} \tan i } \frac{\partial R}{\partial i} , \label{14.4.4} \\[5pt] \dot{Ω} &= \frac{1}{\sqrt{GMm^2 (1 - e^2) \sin i}} \frac{\partial R}{\partial i} , \label{14.4.5} \\[5pt] \dot{T} &= \frac{2a^2}{GMm} \frac{\partial R}{\partial a} + \frac{a(1 - e^2)}{GMme} \frac{\partial R}{\partial e}. \label{14.4.6} \end{align}\]