Enačba Ciolkovskega
Del niza o področju |
Astrodinamika |
---|
Enačba Ciolkovskega (tudi formula ~,[a] (osnovna) raketna ~, enačba idealne rakete in Oberthova enačba) v astrodinamiki in raketodinamiki opisuje gibanje vesoljskega plovila, ki sledi osnovnemu načelu rakete, naprave, ki lahko s pomočjo pospeška nase s potiskom in izpuščanjem dela svoje mase z veliko hitrostjo, ter se tako giblje zaradi ohranitve gibalne količine.
Enačba povezuje delto v (največjo spremembo hitrosti) rakete, če nanjo ne deluje nobena zunanja sila) z efektivno izpušno hitrostjo in začetno in končno maso rakete ali drugi reaktivni motor.
Za vsak tak manever (ali polet, ki vsebuje niz takšnih manevrov) je enačba Ciolkovskega enaka:[1]
kjer je:
- – delta v, največja sprememba hitrosti plovila (brez delovanja zunanjih sil),
- – začetna skupna masa vklučno s pogonskim gorivom (mokra masa),
- – končna skupna masa brez pogonskega goriva (suha masa),[1]
- – celotni čas izgorevanja pogonskega goriva (čas delovanja raketnega motorja)[1]
- razmerje končne in začetne skupne mase,[b]
- – efektivna hitrost izhajanja izpušnega plina v okvirju plovila,
- – funkcija naravnega logaritma.
Enačba se lahko zapiše tudi s pomočjo specifičnega sunka sile (impulza) namesto efektivne izpušne hitrosti s pomočjo formule :
kjer je:
- – standardni težni pospešek = 9,80665 m/s2,
- – specifični sunek sile (impulz), izražen kot časovna perioda v sekundah.
Zgodovina
[uredi | uredi kodo]Enačba se imenuje po Konstantinu Edvardoviču Ciolkovskem, ki jo je neodvisno izpeljal in objavil leta 1903 v svojem delu Raziskovanje vesoljskih prostorov z reaktivnimi pripravami (Изслѣдованіе мировыхъ пространствъ реактивными приборами).[A 1] Uvedel jo je 10. (22.) maja 1897 v rokopisu Raketa (Ракета).[A 2] Enačbo je pred tem izpeljal angleški matematik William Moore leta 1813,[A 3] Peter Guthrie Tait in William John Steele sta jo izpeljala leta 1856 v delu A treatise on the dynamics of a particle.[A 4] Škotski znanstvenik in minister William Leitch je izpeljal osnove raketne tehnike leta 1861.
Čeprav je enačba preprosta vaja v infinitezimalnem računu, jo je Ciolkovski prvi uporabil pri vprašanju ali lahko rakete dosežejo potrebne hitrosti za vesoljske polete. Pri tem morajo biti veliko bolje energetsko grajene z razliko od dotedanjih artilerijskih raket gnanih s smodnikom. Enačba se lahko izpelje z integriranjem diferencialne enačbe Meščerskega sistema točkastih teles s spremenljivimi masami:
oziroma, če ni zunanjih sil ():
kjer je:
- – masni pretok.
Enačbo je leta 1897 izpeljal Ivan Vsevolodovič Meščerski.[2] Jiří Buquoy je pred njim raziskal sisteme s spremenljivo maso. Leta 1812 je eksplicitno formuliral pravilno dinamično enačbo gibanja za primer kadar se masa gibajočega telesa spreminja.[A 5][3] Njegovo delo so leta 1815 predstavili v Francoski akademiji znanosti. Razen kratkega Poissonovega članka Sur le Mouvement d’un Système de Corps, en supposant le Masses Variables iz leta 1819 njegove zamisli niso vzbudile posebnega zanimanja in so se sčasoma pozabile.[A 6] Kakor je poudaril Meščerski, je enačbo gibanja planeta s spremenljivo maso pri stalnih trkih in ločitvah izpeljal von Seeliger leta 1890.[A 7][4] Cayley je leta 1857 obravnaval sorodni razred dinamičnih problemov trkov sistemov z delci infinitezimalne mase in vpeljal variacijsko načelo:[A 8][4]
Cayley je posebej obravnaval verigo padajočo z mize. Delo Meščerskega na zahodu ni bilo znano do leta 1949, ko so izdali njegova zbrana dela.
Robert Esnault-Pelterie, Robert Hutchings Goddard, Hermann Oberth (1924) in mnogi drugi so neodvisno od Ciolkovskega kasneje zapisali enačbo.
Ciolkovski je enačbo zapisal v obliki:
kjer je:
- – masa pogonskega goriva,
- – razmerje pogonske mase.
Količina se imenuje število Ciolkovskega.[5][6][7]
Količine , , in so si med seboj sorazmerne po razpredelnici (v mnogih virih imajo različne označbe):
količina | definicija | ||||
---|---|---|---|---|---|
– | |||||
– | |||||
– | |||||
– |
Izpeljava
[uredi | uredi kodo]Izpeljava z izrekom o gibalni količini
[uredi | uredi kodo]Obravnava se naslednji sistem s spremenljivo maso:
V tej izpeljavi »raketa« pomeni »raketa in vso njeno neizgorelo pogonsko gorivo.«
2. Newtonov zakon povezuje zunanje sile () s spremembo gibalne količine celotnega sistema (vključno z raketo in izpuhom), kot sledi:
kjer je gibalna količina rakete v času :
in gibalna količina rakete in izpuščene mase v času :
in kjer je glede na opazovalca:
- – hitrost rakete v času ,
- – hitrost rakete v času ,
- – hitrost mase dodane izpuhu (in izgubljene za raketo) v času ,
- – masa rakete v času ,
- – masa rakete v času .
Izpušna hitrost v opazovalčevem koordinatnem sistemu je povezana z izpušno hitrostjo v koordinatnem sistemu rakete z zvezo (ker ima izpušna hitrost negativno smer):
Sledi:
in, če se vzame , ker izpuh pozitivnega dela mase povzroči zmanjšanje mase rakete:
Brez zunanjih sil
[uredi | uredi kodo]Če ni zunanjih sil, po izreku o gibalni količini velja in potem:
Od tod se diferenciala lahko pokrajšata:
Če se privzame, da je izpušna hitrost konstantna, se lahko enačba integrira z ločitvijo spremenljivk:
kar da:
oziroma standardno obliko enačbe Ciolkovskega:
ali zapisano enakovredno:
ali ali
kjer je:
- – začetna skupna masa vključno s pogonskim gorivom,
- – končna skupna masa,
- – hitrost raketnega izpuha glede na raketo (specifični sunek sile ali, če se računa v času, pomnožen s standardnim težnim pospeškom na Zemlji).
Vrednost je celotna masa pogonskega goriva, ki je izgorelo, in zato:
kjer je razmerje pogonske mase (del začetne skupne mase, ki se porabi kot delovna masa).
Za ubežno hitrost z Zemlje s konstantno tipično hitrostjo rakete je potrebno razmerje enako:
to je skoraj 94 % mase pri izstrelitvi mora biti pogonsko gorivo.
(delta v) je integral po času velikosti pospeška, ki ga proizvede raketni motor (kar bi bil dejanski pospešek, če ne bi bilo zunanjih sil). V prostem prostoru za primer pospeška v smeri hitrosti je to povečanje hitrosti. V primeru pospeška v nasprotni smeri (zaviranje) je to zmanjšanje hitrosti. Tudi gravitacija in aerodinamični upor pospešujeta ali zavirata plovilo, in se lahko prištejeta ali odštejeta od spremembe hitrosti, ki jo plovilo doživi. Zato delta v po navadi ni dejanska sprememba hitrosti ali hitrost plovila.
Vpliv gravitacije
[uredi | uredi kodo]Če se upošteva konstantno gravitacijsko polje s standardnim težnim pospeškom , ki deluje v nasprotni smeri gibanja rakete, ima izrek o gibalni količini obliko:
Enačba se lahko integrira:
kar da:[8]
oziroma:
Tu je celotni izgorevalni čas in konstantni masni pretok izgorevanja pogonskega goriva, ki odraža zmanjševanje začetne mase do končne mase v tem času:
V primeru odprtega prostora brez gravitacijskih vplivov je vseeno kako dolgo izgoreva gorivo. Sprememba hitrosti ne vpliva na izgorevalni čas. V prisotnosti gravitacije pa ima velik vpliv. Člen pove, da daljši ko je izgorevalni čas, manjša bo sprememba hitrosti rakete. Število Ciolkovskega je v tem primeru enako:
Posebna teorija relativnosti
[uredi | uredi kodo]V klasični mehaniki 2. Newton zakon velja le za zaprte sisteme s konstantno maso (, ). V takšnih sistemih ni razlike med obema oblikama zakona:
Obe obliki sta za sistem s konstantno maso invarianti za Galilejevo transformacijo:
kjer je hitrost v inercialnem opazovalnem sistemu relativnem krajevnemu. Če se Galilejeva transformacija zapiše za obe obliki, sledi, da sta tudi sili v obeh sistemih enaki:
Če masa ni konstantna po času, sili v obeh sistemih nista enaki:
Časovni odvod desne strani da:
in tako:
Druga oblika zakona velja tako za zaprte sisteme s konstantno maso in odprte sisteme s spremenljivo maso. Tako tudi v posebni teoriji relativnosti druga oblika zakona velja, prva pa ne, in je potrebna relativistična definicija gibalne količine :
kjer je Lorentzev faktor. Če se upošteva posebna teorija relativnosti, se lahko izpelje naslednja enačba za relativistično raketo,[9] kjer je spet končna hitrost rakete (ko je porabila vso svojo reakcijsko maso in je preostala njena mirovna masa ) v inercialnem opazovalnem sistemu, kjer je raketa začela polet v mirovanju (z mirovno maso vključno z začetno maso pogonskega goriva ), in hitrost svetlobe v vakuumu:
Enačba se lahko preuredi kot:
Če se vzameta zvezi:
- (tukaj »exp« označuje eksponentno funkcijo; glej tudi naravni logaritem in tudi enakost »potenc« v logaritemskih enakostih) in:
- (glej hiperbolična funkcija),
je to enakovredno izrazu:
Druge izpeljave
[uredi | uredi kodo]S sunkom sile
[uredi | uredi kodo]Enačba Ciolkovskega se lahko izpelje tudi iz osnovnega integrala pospeška v obliki sile (potiska) po masi. Enačbe delte v se zapiše v obliki:
kjer je:
- – potisk
- – začetna (suha) masa,
- – začetna masa minus končna (suha) masa.
Integral rezultante sil po času je celotni sunek sile , če se privzame, da deluje le potisk:
Integral je enak:
Sunek sile po spremembi mase je enak sili po velikosti masnega toka pogonskega goriva (), ki je sama enakovredna izpušni hitrosti:
Integral se lahko izenači z:
S pospeškom
[uredi | uredi kodo]Naj raketa v prostoru miruje brez delovanja zunanjih sil (1. Newtonov zakon). Od trenutka, ko se vžgejo njeni motorji (ura se nastavi na 0), raketa izpušča plinsko maso s konstantno velikostjo masnega toka (kg/s) in z izpušno hitrostjo relativno na raketo (m/s). To povzroča konstantno potisno silo, ki poganja raketo, in, ki je enaka:
Masa goriva, ki jo ima raketa na krovu na začetku, je enaka . Velikost masnega toka je definirana kot razmerje med skupno mokro maso rakete in izgorevalnim časom rakete, tako da bo pretekel čas , ko bo vse gorivo izgorelo. Na raketo tako deluje konstantna sila (), istočasno pa se njena skupna masa stalno zmanjšuje, ker izpušča plin. Po 2. Newtonovem zakonu ima to lahko le eno posledico – njen pospešek se stalno povečuje. Za določitev pospeška je treba pogonsko silo deliti s skupno maso rakete. Tako bo pospešek v poljubnem času po zagonu motorjev in preden se porabi vso gorivo dan kot:
Ker je hitrost določeni integral pospeška in, ker je treba integrirati od časa zagona motorjev () do časa, ko zadnje gorivo zapusti raketo (), naslednji določeni integral da hitrost v času, ko se porabi vso gorivo:
Izrazi v enačbi
[uredi | uredi kodo]Delta v
[uredi | uredi kodo]Delta v (dobesedno »sprememba hitrosti«), označena kot Δv, izgovorjena kot delta-ve, kakor se rabi v dinamiki leta vesoljskih plovil, je mera sunka sile, ki je potreben za izvedbo manevra, kot je izstrelitev ali pristanek na planetu ali naravnemu satelitu, ali v vesoljskem prostoru kot orbitalni manever. Je skalarna količina in ima enoto hitrosti. Če se rabi v tem kontekstu, ni enaka kot fizikalna sprememba hitrosti plovila.
Delto v proizvajajo reakcijski motroji, kot so na prier raketni motorji, in je sorazmerna s potiskom na enoto mase in izgorevalnim časom. Rabi se za določevanje mase raketnega goriva, potrebnega za dani manever prek enačbe Ciolkovskega.
Za mnogokratne manevre se delte v seštevajo linearno.
Za medplanetarne odprave se delta v običajno prikaže s konturnim grafom svinjske zarebrnice, ki prikazuje zahtevano delto v odprave kot funkcijo datuma izstrelitve.
Masno razmerje
[uredi | uredi kodo]V zrakoplovno-vesoljski tehniki je masno razmerje goriva del mase plovila, ki ne doseže cilja, in se po navadi rabi kot mera učinkovitosti plovila. Masno razmerje goriva je razmerje med maso goriva in začetno maso plovila. V vesoljskem plovilu je cilj po navadi tir, za letalo pa njegov kraj pristanka. Večje masno razmerje predstavlja manjšo težo konstrukcije. Druga pomembna mera je tovorno razmerje, ki je del začetne teže kot tovora.
Efektivna izpušna hitrost
[uredi | uredi kodo]Efektivna izpušna hitrost je običajno določena kot specifični impulz (sunek sile). Obe količini sta povezani med seboj kot:
kjer je:
- – specifični impulz (sunek sile), merjen v m/s, kar je enako efektivni izpušni hitrosti, merjeni v m/s,
- – standardni težni pospešek = 9,80665 m/s2,
- – specifični sunek sile, izražen kot časovna perioda v sekundah.
Uporabnost
[uredi | uredi kodo]Enačba Ciolkovskega vsebuje osnove fizike poleta rakete v eni kratki enačbi. Velja tudi za plovila s pogoni podobnim raketnim, kadar je efektivna izpušna hitrost konstantna in se lahko sešteva ali integrira, kadar se efektivna izpušna hitrost spreminja. Enačba Ciolkovskega obravnava le reakcijsko silo iz raketnega motorja – pri tem ne upošteva drugih sil, ki bi lahko delovale na raketo, kot so na primer aerodinamične ali gravitacijske sile. Če se rabi za izračun potrebnega raketnega goriva za izstrelitev s planeta (ali pristanek na) planetu z atmosfero, je treba vplive teh sil upoštevati v potrebni delti v (glej zglede spodaj). V »trinoštvu raketne enačbe«, kakor se imenuje, obstaja meja količine tovora, ki ga raketa lahko nosi s seboj, saj večje količine raketnega goriva povečujejo skupno maso in tako povečujejo tudi porabo goriva.[10] Enačba ne velja za neraketne sisteme, kot so: aerozaviranje, vesoljski topovi, vesoljska dvigala, izstrelitvene zanke ali vesoljski povodci.
Enačba Ciolkovskega se lahko rabi za orbitalne manevre za izračun koliko raketnega goriva je potrebno za spremembo na določeno novo orbito ali za poiskanje nove orbite kot rezultat določenega izgoretja raketnega goriva. Kadar se uporabi pri orbitalnih manevrih, se privzame manever sunka sile v katerem se raketno gorivo porabi, delta v pa se uporabi trenutno. To privzetje je relativno točno za kratkočasovna izgorevanja, kot so popravki s srednjim dosegom in manevri orbitalne vložitve. Ko se čas izgorevanja povečuje, je rezultat manj točen zaradi vpliva gravitacije na plovilo med samim manevrom. Za pogone z dolgim trajanjem in majhnim potiskom, kot je na primer električni pogon, se za napovedovanje orbitalnega gibanja uporabi bolj zamotana analiza na podlagi širjenja vektorja stanja vesoljskega plovila in integracije potiska.
Zgledi
[uredi | uredi kodo]Masa gorivnega sredstva v odvisnosti od začetne celotne mase rakete sledi iz enačbe Ciolkovskega ali po definiciji:
Efektivna izpušna hitrost naj je 4.500 m/s in delta v 9.700 m/s (Zemlja na NZO, ki vključuje za preseganje gravitacije in aerodinamičnega upora).
- enostopenjska raketa do orbite: . Tako mora biti vsaj 88,4 % začetne celotne mase raketno gorivo. Preostalih 11,6 % je namenjeno konstrukciji, motorjem, rezervoarju in tovoru.
- dvostopenjska raketa do orbite: naj prva stopnja proizvede 5.000 m/s; . Tako mora biti vsaj 67,1 % začetne celotne mase raketno gorivo za prvo stopnjo. Preostale mase je 32,9 %. Ko se prva stopnja odvrže, preostane ta masa enaka 32,9 % minus masa rezervoarja in motorjev prve stopnje. Naj je te mase 8 % začetne celotne mase, in tako preostane 24,9 % mase. Druga stopnja naj zagotovi 4.700 m/s; . Tako mora biti vsaj 64,8 % preostale mase raketno gorivo, kar je 16,2 % izvirne celotne mase, in preostane 8,7 % mase za rezervoar in motorje druge stopnje, tovor, ter v primeru raketoplana tudi orbiter. Tako je skupaj 32,9 − 16,2 = 16,7 % izvirne izstrelitvene mase na razpolago za vse motorje, rezervoarje in tovor.
Stopnje
[uredi | uredi kodo]V primeru zaporedno potiskanih raketnih stopenj enačba Ciolkovskega velja za vsako stopnjo posebej, kjer je za vsako stopnjo začetna masa v enačbi skupna masa rakete po odvržbi predhodnje stopnje, končna masa v enačbi pa je skupna masa rakete pred izvržbo obravnavane stopnje. Za vsako stopnjo je specifični sunek sile lahko različen.
Če je na primer 80 % mase rakete raketno gorivo prve stopnje, 10 % suha masa prve stopnje in 10 % preostale mase rakete, je delta v enaka:
S podobnimi zaporednimi manjšimi stopnjami z enako za vsako stopnjo je delta v enaka:
tovor pa je 10 · 10 · 10 = 0,1 % začetne mase.
Primerljiva enostopenjska raketa do orbite s tovorom tudi 0,1 % začetne mase ima lahko maso 11,1 % za gorivne rezervoarje in motorje, ter 88,8 % za gorivo. To da delto v:
Če se motor nove stopnje prižge preden se je odvrgla predhodnja stopnja in, če imajo sočasno delujoči motorji različni specifični sunek sile, kakor je večkrat v primeru trdogorivnih raketnih potisnikih in tekočegorivne stopnje, so razmere bolj zapletene.
Pogostno napačno pojmovanje
[uredi | uredi kodo]Kadar se raketa obravnava kot sistem s spremenljivo maso, se ne more neposredno analizirati z 2. Newtonovim zakonom, saj zakon velja le za sisteme s konstantno maso.[11][12][13] Lahko pride do zmede, ker je enačba Ciolkovskega podobna relativistični enačbi za silo:
Če se uporabi ta formula in kot spremenljiva masa rakete, se zdi, da je enačba Ciolkovskega na ta način izpeljana. Vendar izpeljava ni pravilna. Pri tem se efektivna izpušna hitrost sploh ne pojavi v tej fromuli.
Glej tudi
[uredi | uredi kodo]- zaloga delte v
- masno razmerje
- Oberthov pojav
- relativistična raketa
- pogon vesoljskega plovila
- sistem s spremenljivo maso
- delovna masa
Opombe
[uredi | uredi kodo]- ↑ Največkrat v ruskih virih (формула Циолковского).
- ↑ Masno razmerje je največkrat definirano kot razmerje med začetno in končno skupno maso, označeno kot R – v tem zapisu zato obratna vrednost .
Sklici (A)
[uredi | uredi kodo]- ↑ Ciolkovski (1903).
- ↑ Архив Российской академии наук (АРАН). Ф. 555. Оп. 1. Д. 32. Лл. 1—2, 5, 11, 20.
- ↑ Moore (1813).
- ↑ Tait; Steele (1856).
- ↑ Buquoy (1812).
- ↑ Poisson (1819).
- ↑ Seeliger (1890).
- ↑ Cayley (1857).
Sklici
[uredi | uredi kodo]- ↑ 1,0 1,1 1,2 Strnad (1977), Izreki gibanja, Izrek o gibalni količini, Nasprotna sila curka, str. 72.
- ↑ Kosmodemjanski (1949), str. 13.
- ↑ Pesce; Casetta (2007).
- ↑ 4,0 4,1 Irschik (2012), str. 718.
- ↑ Borodovski (1998).
- ↑ Kosmodemjanski (1969), str. 116.
- ↑ Alemasov (1962).
- ↑ Peraire; Widnall (2008), str. 6.
- ↑ Forward (1995).
- ↑ Pettit (2012).
- ↑ Plastino; Muzzio (1992). »Lahko zaključimo s poudarkom, da 2. Newtonov zakon velja le za konstantno maso. Kadar se masa spreminja zaradi akrecije ali ablacije, je treba uporabiti [alternativno enačbo, ki eksplicitno upošteva spreminjajočo se maso].«
- ↑ Halliday; Resnick (1978), str. 199.
- ↑ Kleppner; Kolenkow (1973), str. 133–134.
Viri
[uredi | uredi kodo]- Alemasov, Vjačeslav Jevgenjevič (1962), Теория ракетных двигателей, Ripol Klassik, ISBN 9785458258111
- Awrejcewicz, Jan (2011), Classical mechanics: dynamics (PDF), doi:10.1007/978-1-4614-3740-6, ISBN 978-1-4614-3739-0, arhivirano iz prvotnega spletišča (PDF) dne 24. junija 2018, pridobljeno 23. junija 2018
- Borodovski, Vitalij Nikolajevič (1998), Отечественные ракеты. История и будущее, IMPB RAN
- Bors, Dorota; Stańczy, Robert (2012), Some nonlinear second order equation modelling, arXiv:1212.5826
- Buquoy, Jiří (1812), Analytische Bestimmung des Gesetzes der Virtuellen Geschwindigkeiten in Mechanischer und Statischer Hinsicht, Leipzig: Breitkopf und Härtel, doi:10.3931/e-rara-14843
- Cayley, Arthur (1857), »On a class of dynamical problems«, Proceedings of the Royal Society of London, VIII: 506–511
- Ciolkovski, Konstantin Edvardovič (1903), Изслѣдованіе мировыхъ пространствъ реактивными приборами (razpoložljivo na spletu tukaj Arhivirano 2011-08-15 na Wayback Machine. v formatu RAR PDF)
- Forward, Robert Lull (1995), A Transparent Derivation of the Relativistic Rocket Equation (PDF), Ameriški inštitut za aeronavtiko in astronavtiko (AIAA), arhivirano iz prvotnega spletišča (PDF) dne 6. septembra 2018, pridobljeno 12. junija 2018 (Glej desno stran enačbe 15 na zadnji strani, kjer je z R označeno masno razmerje – razmerje med začetno in končno maso (obratna vrednost razmerja μ1), z w pa izpušna hitrost, kar odgovarja ve v zapisu tega članka)
- Halliday, David; Resnick, Robert (1978), Physics, zv. 1, str. 199, ISBN 0-471-03710-9,
Pomembno je poudariti, da se ne da izpeljati splošni izraz za 2. Newtonov zakon za sisteme s spreminjajočo maso, če se v izrazu F = dP/dt = d(Mv) masa obravnava kot spremenljivka. [...] Lahko se uporabi izraz F = dP/dt za analizo sistemov s spremenljivo maso le, če se ga uporabi za celotni sistem s konstantno maso z deli, med katerimi se mase izmenjujejo.
[Poudarki kot v izvirniku] - Irschik, Hans (2012), »The Cayley variational principle for continuous-impact problems: A continuum mechanics based version in the presence of a singular surface« (PDF), Journal of Theoretical and Applied Mechanics, varšava, 50 (3): 717–727
- Kleppner, Daniel; Kolenkow, Robert (1973), An Introduction to Mechanics, McGraw-Hill, str. 133–134, ISBN 0-07-035048-5,
Spomnimo, da je izraz F = dP/dt izpeljan za sistem sestavljen iz določene množice delcev[. ... B]istveno je obravnavati isto množico delcev čez cel časovni interval[. ...] Posledično se masa sistema ne more spreminjati v tem obravnavanem časovnem intervalu.
- Kosmodemjanski, Arkadij Aleksandrovič (1949), »Научная деятельность Ивана Всеволодовича Мещерского, str. 9–25«, v Meščerski, Ivan Vsevolodovič, Работы по механике тел переменной массы (1. izd.), Moskva: GITTL
- Kosmodemjanski, Arkadij Aleksandrovič (1969), »First works by K. E. Tsiolkovsky and I. V. Meshchersky on rocket dynamics«, v Essays on the History of Rocketry and Astronautics, str. 115–124
- Levantovski, Lev Vladimirovič (1980), Механика космического полета в элементарном изложении, Moskva: Nauka
- Lundin, Lukas (2016), The Rocket Equation (PDF)
- Matvejev, Aleksej Nikolajevič (1976), Механика и теория относительности, Moskva: Visšaja škola, COBISS 13653081 (angleški prevod Mechanics and theory of relativity (1989))
- Moore, William; (Kraljeva vojaška akademija, Woolwich) (1813), A Treatise on the Motion of Rockets. To which is added, An Essay on Naval Gunnery, London: G. and S. Robinson
- Peraire, Jaime; Widnall, Sheila (2008), Variable Mass Systems: The Rocket Equation (PDF), MIT OpenCourseWare, MIT
- Pesce, Celso Pupo; Casetta, Leonardo (2007), »Variable mass systems dynamics in engineering mechanics education« (PDF), Proceedings of COBEM 2007
- Pettit, Don (5. januar 2012), »NASA – The Tyranny of the Rocket Equation«, www.nasa.gov (v angleščini), arhivirano iz prvotnega spletišča dne 13. aprila 2016, pridobljeno 18. aprila 2016
- Plastino, Angel Ricardo; Muzzio, Juan Carlos (1992), »On the use and abuse of Newton's second law for variable mass problems«, Celestial Mechanics and Dynamical Astronomy, Nizozemska: Kluwer Academic Publishers, 53 (3): 227–232, Bibcode:1992CeMDA..53..227P, doi:10.1007/BF00052611, ISSN 0923-2958, pridobljeno 30. decembra 2011
- Poisson, Siméon-Denis (1819), »Sur le Mouvement d'un Système de Corps, en supposant le Masses Variables«, Bull. Sci. Soc. Philomat., Pariz (Avril): 60–62
- Rodrigues, Hilário Antônio; Panza, Nelson; Portes, Dirceu Jr; Soares, Alexandre de Souza (2014a), »A model of oscillator with variable mass« (PDF), Revista Mexicana de Física, 60: 31–38
- Rodrigues, Hilário Antônio; Panza, Nelson; Portes, Dirceu Jr; Soares, Alexandre de Souza (2014b), »The motion of a leaking oscillator: a study for the physics class«, Physics Education, 49 (5): 557–563, arXiv:1407.5200, doi:10.1088/0031-9120/49/5/557, ISSN 1361-6552
- Romano, Giovanni; Barretta, Raffaele; Diaco, Marina (2017), »Solid–fluid interaction: a continuum mechanics assessment« (PDF), Acta Mechanica, 228: 851–869, doi:10.1007/s00707-016-1738-7, ISSN 1619-6937
- Seeliger, Hugo (1890), »Über Zusammenstösse und Theilungen planetarischer Massen«, Abh. der Königl. Bayer. Akademie der Wiss., Cl. II, XVII (Abth. II): 459–490
- Strnad, Janez (1977), Fizika, 1. del, Mehanika, toplota, Ljubljana: Državna založba Slovenije, COBISS 4171521
- Tait, Peter Guthrie; Steele, William John (1856), A treatise on the dynamics of a particle