Équation de Tsiolkovski

L'équation de Tsiolkovski est l'équation fondamentale de l'astronautique, reliant l'accroissement de vitesse au cours d'une phase de propulsion d'un astronef doté d'un moteur à réaction au rapport de sa masse initiale à sa masse finale.

On la doit à Constantin Tsiolkovski et, indépendamment, à Hermann Oberth.

Histoire

L'équation de Tsiolkovski est considérée comme l'équation fondamentale de l'astronautique. Son éponyme est Constantin Tsiolkovski (1857-1935), qui l'a déduite puis publiée en 1903,,.

Bien que cette équation soit souvent attribuée au grand penseur autodidacte Tsiolkovski, une forme de l'équation figure déjà dans un traité du mathématicien britannique William Moore (en) (fl. c.1806-1823) paru en 1813 ,,, puis dans un article du général-major belge Casimir-Érasme Coquilhat (1811-1890) paru en 1873 , ,,.

L'expérience de la barque, de Tsiolkovski

Expérience de la barque, de Tsiolkovski.

Afin de faire comprendre le principe de la propulsion à réaction, Constantin Tsiolkovski a proposé une fameuse expérience « de la barque ». Une personne se trouve démunie d'avirons dans une barque à l'écart de la rive. Elle veut rejoindre cette rive. Elle remarque que la barque est chargée d'une certaine quantité de pierres et a l'idée de lancer, une à une et le plus vivement possible, ces pierres dans la direction opposée à la rive. Effectivement, à la quantité de mouvement des pierres jetées dans un sens correspond une quantité de mouvement égale pour la barque dans l'autre sens.

Énoncé

L'équation de Tsiolkovski s'écrit :

Δ v → = v → f − v → i = − v → e ln ⁡ m i m f {\displaystyle \Delta {\vec {v}}={\vec {v}}_{\mathrm {f} }-{\vec {v}}_{\mathrm {i} }=-{\vec {v}}_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}}

\Delta {\vec {v}}={\vec {v}}_{\mathrm {f} }-{\vec {v}}_{\mathrm {i} }=-{\vec {v}}_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}

, vectoriellement, ou Δ v = v e ln ⁡ m i m f {\displaystyle \Delta v=v_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}}

\Delta v=v_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}

, numériquement

où :

Δ v {\displaystyle \Delta v} $\Delta v$ est la variation de vitesse de l'astronef entre la fin et le début de la phase propulsée considérée ;
v e {\displaystyle v_{\mathrm {e} }} $v_{\mathrm {e} }$ est la vitesse d'éjection des gaz ;
m i {\displaystyle m_{\mathrm {i} }} $m_{\mathrm {i} }$ est la masse totale de l'astronef au début de la phase propulsée (indice i pour initial) ;
m f {\displaystyle m_{\mathrm {f} }} $m_{\mathrm {f} }$ est la masse totale de l'astronef à l'issue de la phase propulsée (indice f pour final), exprimée dans la même unité que m i {\displaystyle m_{\mathrm {i} }} $m_{\mathrm {i} }$ ;
ln {\displaystyle \ln } $\ln$ est la fonction logarithme népérien.

Établissement

Démonstration

Cette équation est établie en intégrant l'équation de conservation de la quantité de mouvement entre le début et la fin de la phase propulsée sous les hypothèses suivantes :

l'étude du mouvement est faite dans un référentiel d'inertie ;
l'astronef n'est soumis qu'à la force de poussée fournie par ses moteurs, aucune autre action extérieure (gravité, efforts aérodynamiques) n'est prise en compte (voir à la fin de l'article pour la prise en compte de la gravité) ;
la vitesse d'éjection des gaz est constante.

À un instant donné, lorsque le vaisseau de masse m {\displaystyle m} $m$ se déplaçant à la vitesse v → {\displaystyle {\vec {v}}} ${\vec {v}}$ éjecte une petite quantité d'ergol à la vitesse v e → {\displaystyle {\vec {v_{\mathrm {e} }}}} ${\vec {v_{\mathrm {e} }}}$ , on note d m {\displaystyle \mathrm {d} m} $\mathrm {d} m$ sa variation de masse et d v → {\displaystyle \mathrm {d} {\vec {v}}} $\mathrm {d} {\vec {v}}$ sa variation de vitesse. La variation de quantité de mouvement du système isolé (vaisseau + ergol éjecté) est nécessairement nulle, il vient donc :

m ⋅ d v → + ( − d m ) ⋅ v e → = 0 → {\displaystyle m\cdot \mathrm {d} {\vec {v}}+(-\mathrm {d} m)\cdot {\vec {v_{\mathrm {e} }}}={\vec {0}}}

m\cdot \mathrm {d} {\vec {v}}+(-\mathrm {d} m)\cdot {\vec {v_{\mathrm {e} }}}={\vec {0}}

Pour obtenir la variation de vitesse Δ v → {\displaystyle \Delta {\vec {v}}} $\Delta {\vec {v}}$ de l'astronef quand sa masse passe de m i {\displaystyle m_{\mathrm {i} }} $m_{\mathrm {i} }$ à m f {\displaystyle m_{\mathrm {f} }} $m_{\mathrm {f} }$ , on peut intégrer cette petite variation de vitesse :

Δ v → = ∫ v → = v i → v f → d v → = ∫ m = m i m f v e → m d m = v e → ∫ m i m f d m m = v e → ln ⁡ m f m i = − v e → ln ⁡ m i m f {\displaystyle \Delta {\vec {v}}=\int _{{\vec {v}}={\vec {v_{\mathrm {i} }}}}^{\vec {v_{\mathrm {f} }}}\mathrm {d} {\vec {v}}=\int _{m=m_{\mathrm {i} }}^{m_{\mathrm {f} }}{\frac {\vec {v_{\mathrm {e} }}}{m}}\mathrm {d} m={\vec {v_{\mathrm {e} }}}\int _{m_{\mathrm {i} }}^{m_{\mathrm {f} }}{\frac {\mathrm {d} m}{m}}={\vec {v_{\mathrm {e} }}}\ln {\frac {m_{\mathrm {f} }}{m_{\mathrm {i} }}}=-{\vec {v_{\mathrm {e} }}}\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}}

\Delta {\vec {v}}=\int _{{\vec {v}}={\vec {v_{\mathrm {i} }}}}^{\vec {v_{\mathrm {f} }}}\mathrm {d} {\vec {v}}=\int _{m=m_{\mathrm {i} }}^{m_{\mathrm {f} }}{\frac {\vec {v_{\mathrm {e} }}}{m}}\mathrm {d} m={\vec {v_{\mathrm {e} }}}\int _{m_{\mathrm {i} }}^{m_{\mathrm {f} }}{\frac {\mathrm {d} m}{m}}={\vec {v_{\mathrm {e} }}}\ln {\frac {m_{\mathrm {f} }}{m_{\mathrm {i} }}}=-{\vec {v_{\mathrm {e} }}}\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}

Puisque m i > m f {\displaystyle m_{\mathrm {i} }>m_{\mathrm {f} }} $m_{\mathrm {i} }>m_{\mathrm {f} }$ , ln ⁡ ( m i / m f ) > 0 {\displaystyle \ln(m_{\mathrm {i} }/m_{\mathrm {f} })>0} $\ln(m_{\mathrm {i} }/m_{\mathrm {f} })>0$ , la variation de vitesse du vaisseau a donc comme prévu la même direction que la projection des ergols et en sens opposé.

Commentaires

Il est souvent dit que pour trouver cette équation, il faut que le débit massique d'ergol soit constant pendant la phase de propulsion, mais ce n'est pas obligatoire, même si cela simplifie le travail d'intégration dans un premier temps.

L'équation est valable aussi bien lors d'une phase d'accélération (la poussée est dans la direction de la vitesse, Δ v {\displaystyle \Delta v} $\Delta v$ est positif : c'est un accroissement de vitesse) que de décélération (la poussée est de direction opposée à la vitesse, Δ v {\displaystyle \Delta v} $\Delta v$ est négatif : c'est une réduction de vitesse).

La différence entre la masse initiale m i {\displaystyle m_{\mathrm {i} }} $m_{\mathrm {i} }$ et la masse finale m f {\displaystyle m_{\mathrm {f} }} $m_{\mathrm {f} }$ correspond à la masse que la fusée a éjectée durant sa propulsion ; on l'appelle masse d'appui (« d'appui » parce que c'est la masse sur laquelle la fusée s'est appuyée pour se propulser).

Pour les fusées thermo-chimiques (Ariane, Soyouz, navette, etc.), la masse d'appui est la masse des ergols (poudre, ou dioxygène et dihydrogène), laquelle est également source d'énergie chimique : c'est donc la masse d'appui elle-même qui contient l'énergie qui servira à sa propre éjection. Ce n'est plus le cas pour les moteurs ioniques (qui représentent sans doute l'avenir de la conquête spatiale). Ceux-ci sont régis tout pareillement par l'équation de Tsiolkovski, mais leur masse d'appui est constituée d'un gaz neutre (du xénon) ; c'est la très forte vitesse d'éjection de cette masse d'appui qui rend ces moteurs très économes en masse d'appui (il leur faut cependant une source d'énergie pour réaliser l'éjection). À ce titre, le fonctionnement des moteurs ioniques est comparable à celui des fusées à eau, dans lesquelles l'eau n'est utilisée que pour sa masse (l'énergie résidant dans l'air comprimé).

Dans le cas où la phase propulsée est réalisée au moyen de plusieurs étages fonctionnant successivement, la même équation de Tsiolkovski peut être utilisée pour le vol de chaque étage. On peut ainsi montrer l'intérêt de telles fusées à plusieurs étages. Voir l'exemple dans la section suivante.

Malgré l'apparente simplicité de cette équation et des hypothèses qui la sous-tendent, elle constitue une approximation utile au calcul des manœuvres de changement d'orbite, ces manœuvres étant qualifiées d'impulsionnelles, c’est-à-dire effectuées en un temps suffisamment bref pour que les hypothèses de l'équation de Tsiolkovski restent approximativement valables.

Temps nécessaire

Si le vaisseau utilise un débit massique d'ergol constant q {\displaystyle q} $q$ on peut écrire :

Δ t = Δ m q {\displaystyle \Delta t={\frac {\Delta m}{q}}}

\Delta t={\frac {\Delta m}{q}}

Or l'équation de Tsiolkovski peut s'écrire :

Δ v v e = ln ⁡ m i m i − Δ m {\displaystyle {\frac {\Delta v}{v_{\mathrm {e} }}}=\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {i} }-\Delta m}}}

{\frac {\Delta v}{v_{\mathrm {e} }}}=\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {i} }-\Delta m}}

C'est-à-dire, en changeant de signe et en passant à l'exponentielle :

m i − Δ m m i = e − Δ v / v e {\displaystyle {\frac {m_{\mathrm {i} }-\Delta m}{m_{\mathrm {i} }}}=\mathrm {e} ^{-\Delta v/v_{\mathrm {e} }}}

{\frac {m_{\mathrm {i} }-\Delta m}{m_{\mathrm {i} }}}=\mathrm {e} ^{-\Delta v/v_{\mathrm {e} }}

On en tire l'expression de Δ m {\displaystyle \Delta m} $\Delta m$ qu'on reporte dans celle de Δ t {\displaystyle \Delta t} $\Delta t$ :

Δ t = m i q ( 1 − e − Δ v / v e ) {\displaystyle \Delta t={\frac {m_{\mathrm {i} }}{q}}\,\left(1-\mathrm {e} ^{-\Delta v/v_{\mathrm {e} }}\right)}

\Delta t={\frac {m_{\mathrm {i} }}{q}}\,\left(1-\mathrm {e} ^{-\Delta v/v_{\mathrm {e} }}\right)

Exemple

L'exemple qui suit a pour objet de montrer l'intérêt des fusées à plusieurs étages.

Soit une fusée à deux étages ayant les caractéristiques suivantes :

la masse d'ergols embarqués par chaque étage (premier étage : 100 tonnes ; deuxième étage : 20 tonnes) représente 10 fois sa masse à vide ;
la vitesse d'éjection des gaz est de 4 000 m/s ;

et supposons qu'elle emporte une charge utile de 2 t. Résumons ces données dans un tableau :

Étage	Masse d'ergols (t)	Masse à vide (t)	Masse totale (t)	Vitesse d'éjection des gaz (m/s)
Premier étage	m e 1 = 100 {\displaystyle m_{\mathrm {e} {\mathfrak {1}}}=100} $m_{\mathrm {e} {\mathfrak {1}}}=100$	m v 1 = 10 {\displaystyle m_{\mathrm {v} {\mathfrak {1}}}=10} $m_{\mathrm {v} {\mathfrak {1}}}=10$	m t 1 = 110 {\displaystyle m_{\mathrm {t} {\mathfrak {1}}}=110} $m_{\mathrm {t} {\mathfrak {1}}}=110$	v e = 4 000 {\displaystyle v_{\mathrm {e} }=4\,000} $v_{\mathrm {e} }=4\,000$
Deuxième étage	m e 2 = 20 {\displaystyle m_{\mathrm {e} {\mathfrak {2}}}=20} $m_{\mathrm {e} {\mathfrak {2}}}=20$	m v 2 = 2 {\displaystyle m_{\mathrm {v} {\mathfrak {2}}}=2} $m_{\mathrm {v} {\mathfrak {2}}}=2$	m t 2 = 22 {\displaystyle m_{\mathrm {t} {\mathfrak {2}}}=22} $m_{\mathrm {t} {\mathfrak {2}}}=22$	v e = 4 000 {\displaystyle v_{\mathrm {e} }=4\,000} $v_{\mathrm {e} }=4\,000$
Charge utile			m c u = 2 {\displaystyle m_{\mathrm {cu} }=2} $m_{\mathrm {cu} }=2$
Total fusée	m e = 120 {\displaystyle m_{\mathrm {e} {\mathfrak {}}}=120} $m_{\mathrm {e} {\mathfrak {}}}=120$	m v = 12 {\displaystyle m_{\mathrm {v} {\mathfrak {}}}=12} $m_{\mathrm {v} {\mathfrak {}}}=12$	m t = 134 {\displaystyle m_{\mathrm {t} {\mathfrak {}}}=134} $m_{\mathrm {t} {\mathfrak {}}}=134$

On peut alors mener les calculs d'incréments de vitesse, comme suit, en employant deux fois l'équation de Tsiolkovski, aux étapes 3 et 6 :

Étape de calcul		Formule	Masse (t)	Vitesse (m/s)
1	Masse à l'allumage du premier étage	m i 1 = m t {\displaystyle m_{\mathrm {i} {\mathfrak {1}}}=m_{\mathrm {t} {\mathfrak {}}}} $m_{\mathrm {i} {\mathfrak {1}}}=m_{\mathrm {t} {\mathfrak {}}}$	134 {\displaystyle 134} $134$
2	Masse à l'extinction du premier étage	m f 1 = m i 1 − m e 1 {\displaystyle m_{\mathrm {f} {\mathfrak {1}}}=m_{\mathrm {i} {\mathfrak {1}}}-m_{\mathrm {e} {\mathfrak {1}}}} $m_{\mathrm {f} {\mathfrak {1}}}=m_{\mathrm {i} {\mathfrak {1}}}-m_{\mathrm {e} {\mathfrak {1}}}$	34 {\displaystyle 34} $34$
3	Incrément de vitesse du premier étage	Δ v 1 = v e ln ⁡ m i 1 m f 1 {\displaystyle \Delta v_{1}=v_{\mathrm {e} }\ln {\tfrac {m_{\mathrm {i} {\mathfrak {1}}}}{m_{\mathrm {f} {\mathfrak {1}}}}}} $\Delta v_{1}=v_{\mathrm {e} }\ln {\tfrac {m_{\mathrm {i} {\mathfrak {1}}}}{m_{\mathrm {f} {\mathfrak {1}}}}}$		5 486 {\displaystyle 5\,486} $5\,486$
4	Masse à l'allumage du second étage	m i 2 = m f 1 − m v 1 {\displaystyle m_{\mathrm {i} {\mathfrak {2}}}=m_{\mathrm {f} {\mathfrak {1}}}-m_{v{\mathfrak {1}}}} $m_{\mathrm {i} {\mathfrak {2}}}=m_{\mathrm {f} {\mathfrak {1}}}-m_{v{\mathfrak {1}}}$	24 {\displaystyle 24} $24$
5	Masse à l'extinction du second étage	m f 2 = m i 2 − m e 2 {\displaystyle m_{\mathrm {f} {\mathfrak {2}}}=m_{\mathrm {i} {\mathfrak {2}}}-m_{e{\mathfrak {2}}}} $m_{\mathrm {f} {\mathfrak {2}}}=m_{\mathrm {i} {\mathfrak {2}}}-m_{e{\mathfrak {2}}}$	4 {\displaystyle 4} $4$
6	Incrément de vitesse du deuxième étage	Δ v 2 = v e ln ⁡ m i 2 m f 2 {\displaystyle \Delta v_{2}=v_{\mathrm {e} }\ln {\tfrac {m_{\mathrm {i} {\mathfrak {2}}}}{m_{\mathrm {f} {\mathfrak {2}}}}}} $\Delta v_{2}=v_{\mathrm {e} }\ln {\tfrac {m_{\mathrm {i} {\mathfrak {2}}}}{m_{\mathrm {f} {\mathfrak {2}}}}}$		7 167 {\displaystyle 7\,167} $7\,167$
7	Vitesse finale	Δ v = Δ v 1 + Δ v 2 {\displaystyle \!\,\Delta v=\Delta v_{1}+\Delta v_{2}} $\!\,\Delta v=\Delta v_{1}+\Delta v_{2}$		12 653 {\displaystyle 12\,653} $12\,653$

Par comparaison, une fusée comportant un seul étage avec la même quantité totale d'ergols (120 t) et la même masse à vide totale (12 t) imprimerait à une charge utile de même masse (2 t) une vitesse environ 30 % inférieure :

Étape de calcul		Formule	Masse (t)	Vitesse (m/s)
1	Masse à l'allumage de l'étage (unique)	m i = m t {\displaystyle m_{\mathrm {i} {\mathfrak {}}}=m_{\mathrm {t} {\mathfrak {}}}} $m_{\mathrm {i} {\mathfrak {}}}=m_{\mathrm {t} {\mathfrak {}}}$	134 {\displaystyle 134} $134$
2	Masse à l'extinction de l'étage	m f = m i − m e = m v + m c u {\displaystyle m_{\mathrm {f} {\mathfrak {}}}=m_{\mathrm {i} {\mathfrak {}}}-m_{\mathrm {e} {\mathfrak {}}}=m_{\mathrm {v} {\mathfrak {}}}+m_{cu{\mathfrak {}}}} $m_{\mathrm {f} {\mathfrak {}}}=m_{\mathrm {i} {\mathfrak {}}}-m_{\mathrm {e} {\mathfrak {}}}=m_{\mathrm {v} {\mathfrak {}}}+m_{cu{\mathfrak {}}}$	14 {\displaystyle 14} $14$
3	Vitesse finale	Δ v = v e ln ⁡ m i m f {\displaystyle \Delta v=v_{\mathrm {e} }\ln {\tfrac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}} $\Delta v=v_{\mathrm {e} }\ln {\tfrac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}$		9 034 {\displaystyle 9\,034} $9\,034$

Pertes par pesanteur

Cette section ne cite pas suffisamment ses sources (août 2023). Pour l'améliorer, ajoutez des références de qualité et vérifiables (comment faire ?) ou le modèle {{Référence nécessaire}} sur les passages nécessitant une source.

Les calculs ont été effectués dans l'hypothèse d'une absence de pesanteur (manœuvres en orbite). Lorsque cette gravité agit, un terme simple doit être ajouté à l'équation de Tsiolkovski. Celle-ci devient :

Δ v → = − v → e ln ⁡ m i m f + √ ( 2 g → Δ h ) {\displaystyle \Delta {\vec {v}}=-{\vec {v}}_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}+\surd {\bigl (}2{\vec {g}}\,\Delta h{\bigr )}}

\Delta {\vec {v}}=-{\vec {v}}_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}+\surd {\bigl (}2{\vec {g}}\,\Delta h{\bigr )}

, vectoriellement, ou Δ v = v e ln ⁡ m i m f − √ ( 2 g Δ h ) {\displaystyle \Delta v=\ v_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}-\surd {\bigl (}2g\Delta h{\bigr )}}

\Delta v=\ v_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}-\surd {\bigl (}2g\Delta h{\bigr )}

, si l'on projette l'équation radialement,

g → {\displaystyle {\vec {g}}} ${\vec {g}}$ étant l'accélération locale de la pesanteur et Δ h {\displaystyle \Delta h} $\Delta h$ la différence entre l'altitude de satellisation et celle du point de départ de la fusée. On fait l'hypothèse que g {\displaystyle g} $g$ est constant pendant la propulsion (alors que g {\displaystyle g} $g$ diminue très légèrement avec l'altitude).

Notes et références

Notes

L'équation de Tsiolkovski est aussi connue comme la formule de Tsiolkovski ou la loi de Tsiolkovski.
À partir de la vitesse d'éjection des gaz v e {\displaystyle v_{\mathrm {e} }} $v_{\mathrm {e} }$ , on définit une grandeur dérivée, souvent usitée en astronautique, appelée l’impulsion spécifique du moteur, par la relation I s p = v e g 0 {\displaystyle I_{sp}={\frac {v_{\mathrm {e} }}{g_{0}}}} $I_{sp}={\frac {v_{\mathrm {e} }}{g_{0}}}$ , où g 0 {\displaystyle g_{0}} $g_{0}$ est l’accélération de la pesanteur (9,81 m/s2). L’impulsion spécifique a la dimension d’un temps et s’exprime donc en secondes.
Naturellement, d m < 0 {\displaystyle \mathrm {d} m<0} $\mathrm {d} m<0$ . La masse d'ergol éjecté est donc − d m {\displaystyle -\mathrm {d} m} $-\mathrm {d} m$ .

Références

Bonnal et al. 2014, § 8.2.3, p. 168, col. 1.
Bouchez 2010, § 1.2, p. 3, col. 2.
Rax 2005, § 1.2.3, p. 29.
Macdonald, Norris et Spencer 2014, § 1.1.1, p. 2, col. 1.
Macdonald, Norris et Spencer 2014, réf., p. 24, 1.
Tsiolkovski 1903.
Voir une mise en application de cette formule dans Journal des sciences militaires des armées de terre et de mer, tome II, 1826,
Macdonald, Norris et Spencer 2014, § 1.1.1, p. 2, col. 1-2.
Macdonald, Norris et Spencer 2014, réf., p. 24, 2.
Moore 1813.
Mémoires de la Société royale des sciences de Liège, deuxième série, tome V, p. 315, .
voir le lien ESA :
Macdonald, Norris et Spencer 2014, § 1.1.1, p. 2, col. 2.
Macdonald, Norris et Spencer 2014, réf., p. 24, 3.
Coquilhat 1873.

Voir aussi

Publications originales

(en) Williams Moore, A treatise on the motion of rockets : to which is added an essay on navel gunnery, in theory and practice ; designed for the use of the army and navy, and in all places of military, naval and scientific instruction, Londres, 1813.
Casimir-Érasme Coquilhat, « Trajectoires des fusées volantes dans le vide », Mémoires de la Société royale des sciences de Liège, 2e série, t. V,‎ nov. 1873, art. no 5, −33 p. (OCLC 248870034, lire en ligne), réimpr. :
- Casimir-Érasme Coquilhat, Trajectoires des fusées volantes dans le vide : 1873, Dijon, Nielrow, janv. 2018, 1re éd., 1 vol., 51 (ISBN 978-2-9559619-5-7, EAN 9782955961957, OCLC 1043726396, BNF 45456511).
(ru) Konstantin Èduardovič Ciolkovskij, « Исследовние мировых пространств реактивными приборами » , Научное обозрение, vol. 7, no 5,‎ 1903, p. 45-75 (résumé) — publication originale de l'article de Tsiolkovski énonçant l'équation fondamentale de l'astronautique.

Bibliographie

: document utilisé comme source pour la rédaction de cet article.

Marc Bouchez, « Propulsion aérospatiale : lois simplifiées pour le dimensionnement et exercices d'application », Techniques de l'ingénieur,‎ 10 juillet 2010, art. BM 3 003, § 1.2 (« Formule de Tsiolkovski : cas d'une accélération constante en palier »), p. 3, col. 2 (résumé, lire en ligne).
(en) V. I. Feodosiev et G. B. Siniarev (trad. du russe par S. N. Samburof), Introduction to rocket technology , New York et Londres, Academic Press, janvier 1959 (réimpr. mai 2014), 1re éd., X-344 p., in-8o (15,2 × 23,5 cm) (ISBN 978-1-4832-3201-0, EAN 9781483232010, OCLC 840470516, BNF 44667202, DOI 10.1016/C2013-0-12448-8, présentation en ligne, lire en ligne), chap. Ier (OCLC 7333041790, DOI 10.1016/B978-1-4832-3201-0.50006-1, résumé), section B (« Tsiolkovski equation for the ideal velocity of a rocket ») , p. 16-22.
(en) Malcolm Macdonald et Viorel Badescu (éditeur), The international handbook of space technology , Heidelberg et Chichester, Springer et Praxis, coll. « Praxis books / Astronautical engineering », février 2014, 1re éd., XV-731 p., 21,6 × 27,9 cm (ISBN 978-3-642-41100-7 et 978-3-662-50608-0, OCLC 958469721, BNF 44707981, DOI 10.1007/978-3-642-41101-4, SUDOC 199104573, lire en ligne) :
- (en) Malcolm Macdonald, Pat Norris et David B. Spencer, « Introduction », dans op. cit., chap. 1er , p. 1-24 (OCLC 7327569974, DOI 10.1007/978-3-642-41101-4_1).
- (en) Christophe Bonnal, Alessandro Ciucci, Michael H. Obersteiner et Oskar Haidn, « Launch systems », dans op. cit., chap. 8, p. 165-195 (OCLC 7326880462, DOI 10.1007/978-3-642-41101-4_8.
Jean-Marcel Rax (préf. de Bernard Bigot), Physique des plasmas : cours et applications, Paris, Dunod, coll. « Sciences sup », août 2005, 1re éd., XII-426 p., 17 × 24 cm (ISBN 978-2-10-007250-7, EAN 9782100072507, OCLC 77230441, BNF 40025642, SUDOC 090320409, présentation en ligne, lire en ligne).
(en) Jean-Jacques Serra, Philippe Jung et Théo Pirard, « Casimir Coquilhat's theory on rocket motion : the rocket equation established in 1871 ! », dans John Harlow (éditeur), History of rocketry and astronautics (actes du 42e symposium d'histoire de l'Académie internationale d'astronautique, tenu à Glasgow en 2008), San Diego, Univelt, coll. « AAS history series » (no 39) et coll. « IAA history symposia » (no 28), 2013, 1re éd., XII-345 p., 25 cm (ISBN 978-0-87703-589-3 et 978-0-87703-590-9, EAN 9780877035893, OCLC 840430825, présentation en ligne), 2e partie, chap. 7, p. 121-134 (résumé).