Enačba Ciolkovskega

Iz testwiki
Pojdi na navigacijo Pojdi na iskanje

Predloga:Astrodinamika

Graf obratne vrednosti masnega razmerja rakete (1/μ1) v odvisnosti od njene končne hitrosti (Δv/ve) po enačbi Ciolkovskega.

Enačba Ciolkovskega (tudi formula ~,Predloga:Efn (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]

Δv=velnm0m1=veln1μ1,

kjer je:

Enačba se lahko zapiše tudi s pomočjo specifičnega sunka sile (impulza) namesto efektivne izpušne hitrosti s pomočjo formule ve=g0Isp:

Δv=g0Ispln1μ1,

kjer je:

Zgodovina

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:

m(t)d𝐯dt𝐮1(t)dm1dt+𝐮2(t)dm2dt𝐅=0,

oziroma, če ni zunanjih sil (𝐅=0):

m(t)d𝐯dt+𝐮(t)Φm=0,

kjer je:

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]

[(d2xdt2X)δx+(d2ydt2Y)δy+(d2zdt2Z)δz]dm
(Δv1δξ+Δv2δν+Δv3δζ)1dtdμ=0.

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:

Δvvmaks=veln(1+m2m1)=veln(1+Z)=veln(1ζ),

kjer je:

Količina Z=m2m1=ζ1ζ=1μ11 se imenuje število Ciolkovskega.[5][6][7]

Količine R, μ1, ζ in Z so si med seboj sorazmerne po razpredelnici (v mnogih virih imajo različne označbe):

količina definicija R μ1 ζ Z
R m0m1 1μ1 11ζ 1+Z
μ1 m1m0 1R 1ζ 11+Z
ζ m2m0 11R 1μ1 Z1+Z
Z m2m1 R1 1μ11 ζ1ζ

Izpeljava

Izpeljava z izrekom o gibalni količini

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 (Fi) s spremembo gibalne količine celotnega sistema (vključno z raketo in izpuhom), kot sledi:

iFi=limΔt0G2G1Δt,

kjer je G1 gibalna količina rakete v času t=0:

G1=(m+Δm)v

in G2 gibalna količina rakete in izpuščene mase v času t=Δt:

G2=m(v+Δv)+Δmv'e

in kjer je glede na opazovalca:

  • v – hitrost rakete v času t=0,
  • v+Δv – hitrost rakete v času t=Δt,
  • v'e – hitrost mase dodane izpuhu (in izgubljene za raketo) v času Δt,
  • m+Δm – masa rakete v času t=0,
  • m – masa rakete v času t=Δt.

Izpušna hitrost v'e v opazovalčevem koordinatnem sistemu je povezana z izpušno hitrostjo v koordinatnem sistemu rakete ve z zvezo (ker ima izpušna hitrost negativno smer):

v'e=vve.

Sledi:

G2G1=mΔvveΔm

in, če se vzame dm=Δm, ker izpuh pozitivnega dela mase Δm povzroči zmanjšanje mase rakete:

iFi=mdvdt+vedmdt.

Brez zunanjih sil

Če ni zunanjih sil, po izreku o gibalni količini velja iFi=0 in potem:

mdvdt=vedmdt.

Od tod se diferenciala dt lahko pokrajšata:

mdv=vedm.

Če se privzame, da je izpušna hitrost ve konstantna, se lahko enačba integrira z ločitvijo spremenljivk:

v0v1dv=vem0m1dmm,

kar da:

v1v0=ve(lnm1lnm0)=velnm1m0,

oziroma standardno obliko enačbe Ciolkovskega:

Δv=velnμ1

ali zapisano enakovredno:

μ1=eΔv/ve      ali      1μ1=eΔv/ve      ali      Z=eΔv/ve1,

kjer je:

Vrednost m2=m0m1 je celotna masa pogonskega goriva, ki je izgorelo, in zato:

ζ=1μ1=1eΔv/ve,

kjer je ζ razmerje pogonske mase (del začetne skupne mase, ki se porabi kot delovna masa).

Za ubežno hitrost z Zemlje Δvvu=11.200m/s s konstantno tipično hitrostjo rakete ve=4.000m/s je potrebno razmerje enako:

ζ=1e11200/4000=1e2,8=0,9392,

to je skoraj 94 % mase pri izstrelitvi mora biti pogonsko gorivo.

Δv (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

Če se upošteva konstantno gravitacijsko polje s standardnim težnim pospeškom g0, ki deluje v nasprotni smeri gibanja rakete, ima izrek o gibalni količini obliko:

mdvdt=vedmdtmg0.

Enačba se lahko integrira:

v0v1dv=vem0m1dmmg0t0=0t1dt,

kar da:[8]

v1v0=velnm1m0g0t1=velnm1m0g0m0m1Φm,

oziroma:

Δv=velnμ1g0m2Φm.

Tu je t1 celotni izgorevalni čas in Φm konstantni masni pretok izgorevanja pogonskega goriva, ki odraža zmanjševanje začetne mase m0 do končne mase m1 v tem času:

m1=m0Φmt1.

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 g0t1 pove, da daljši ko je izgorevalni čas, manjša bo sprememba hitrosti rakete. Število Ciolkovskega je v tem primeru enako:

Z=exp[Δv+g0m2Φmve1].

Posebna teorija relativnosti

V klasični mehaniki 2. Newton zakon velja le za zaprte sisteme s konstantno maso (m=konst., Φm=dm/dt=0). V takšnih sistemih ni razlike med obema oblikama zakona:

  1. F=ma=mdvdt,
  2. F=dGdt=ddt(mv)=mdvdt+v(dmdt)=0=mdvdt.

Obe obliki sta za sistem s konstantno maso invarianti za Galilejevo transformacijo:

x=xut;v=vu,

kjer je u 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:

F=F.

Če masa ni konstantna po času, sili v obeh sistemih nista enaki:

F=ddt(mv).

Časovni odvod desne strani da:

F=mdvdt+(vu)dmdt

in tako:

F=ddt(mv)udmdtF

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 Gp=γmv:

F=ddt(γmv),

kjer je γ Lorentzev faktor. Če se upošteva posebna teorija relativnosti, se lahko izpelje naslednja enačba za relativistično raketo,[9] kjer je Δv spet končna hitrost rakete (ko je porabila vso svojo reakcijsko maso in je preostala njena mirovna masa m1) v inercialnem opazovalnem sistemu, kjer je raketa začela polet v mirovanju (z mirovno maso vključno z začetno maso pogonskega goriva m0), in c0 hitrost svetlobe v vakuumu:

1μ1=[1+Δvc01Δvc0]c02ve.

Enačba se lahko preuredi kot:

Δvc0=(1μ1)2vec01(1μ1)2vec0+1.

Če se vzameta zvezi:

(1μ1)2vec0=exp[2vec0ln1μ1] (tukaj »exp« označuje eksponentno funkcijo; glej tudi naravni logaritem in tudi enakost »potenc« v logaritemskih enakostih) in:
thx=e2x1e2x+1 (glej hiperbolična funkcija),

je to enakovredno izrazu:

Δv=c0th(vec0ln1μ1).

Druge izpeljave

S sunkom sile

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:

Δv=t0t1|FT|m0tΔmdt,

kjer je:

  • FTpotisk
  • m0 – začetna (suha) masa,
  • Δm – začetna masa minus končna (suha) masa.

Integral rezultante sil po času je celotni sunek sile I, če se privzame, da deluje le potisk:

t0t1Fdt=I.

Integral je enak:

Ilnm0lnm1Δm.

Sunek sile po spremembi mase je enak sili po velikosti masnega toka pogonskega goriva (Φm), ki je sama enakovredna izpušni hitrosti:

IΔm=FΦm=ve.

Integral se lahko izenači z:

Δv=velnm0m1.

S pospeškom

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 Φm (kg/s) in z izpušno hitrostjo relativno na raketo ve (m/s). To povzroča konstantno potisno silo, ki poganja raketo, in, ki je enaka:

FT=veΦm.

Masa goriva, ki jo ima raketa na krovu na začetku, je enaka m2=m0m1. Velikost masnega toka je definirana kot razmerje med skupno mokro maso rakete in izgorevalnim časom rakete, tako da bo pretekel čas t1=(m0m1)/Φm, ko bo vse gorivo izgorelo. Na raketo tako deluje konstantna sila (mve), 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 t po zagonu motorjev in preden se porabi vso gorivo dan kot:

a(t)=veΦmm0(Φmt).

Ker je hitrost določeni integral pospeška in, ker je treba integrirati od časa zagona motorjev (t0=0) do časa, ko zadnje gorivo zapusti raketo (t1=(m0m1)/Φm), naslednji določeni integral da hitrost v času, ko se porabi vso gorivo:

Δv=v1v0=0t1veΦmm0(Φmt)dt=veln(m0(m0m1))+veln(m00)=velnm0velnm1=velnm0m1.

Izrazi v enačbi

Delta v

Predloga:Glavni

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

Predloga:Glavni

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

Predloga:Glavni

Efektivna izpušna hitrost je običajno določena kot specifični impulz (sunek sile). Obe količini sta povezani med seboj kot:

ve=g0Isp,

kjer je:

Uporabnost

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

Oberthova dvostopenjska raketa Modell-B

Masa gorivnega sredstva v odvisnosti od začetne celotne mase rakete sledi iz enačbe Ciolkovskega ali po definiciji:

m2=(1eΔv/ve)m0=ζm0.

Efektivna izpušna hitrost naj je ve= 4.500 m/s in delta v Δv= 9.700 m/s (Zemlja na NZO, ki vključuje Δv za preseganje gravitacije in aerodinamičnega upora).

  • enostopenjska raketa do orbite: m2=(1e9,7/4,5)m0=0,884m0. 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 Δv= 5.000 m/s; m2=(1e5,0/4,5)m0=0,671m0. 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 Δv= 4.700 m/s; m2=(1e4,7/4,5)m0=0,648m0. 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

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:

Δv =veln10010080=veln5=1,61ve.

S podobnimi zaporednimi manjšimi stopnjami z enako ve za vsako stopnjo je delta v enaka:

Δv =3veln5 =4,83ve,

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:

Δv =veln(100/11,2) =2,19ve.

Č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

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:

F=dpdt=mdvdt+vdmdt.

Če se uporabi ta formula in m(t) 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 ve sploh ne pojavi v tej fromuli.

Glej tudi

Opombe

Predloga:Notelist

Sklici (A)

Predloga:Sklici

Sklici

Predloga:Sklici

Viri

Predloga:Refbegin

Predloga:Refend

Zunanje povezave

Predloga:Portal

  1. 1,0 1,1 1,2 Predloga:Sktxt.
  2. Predloga:Sktxt.
  3. Predloga:Sktxt.
  4. 4,0 4,1 Predloga:Sktxt.
  5. Predloga:Sktxt.
  6. Predloga:Sktxt.
  7. Predloga:Sktxt.
  8. Predloga:Sktxt.
  9. Predloga:Sktxt.
  10. Predloga:Sktxt.
  11. Predloga:Sktxt. »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].«
  12. Predloga:Sktxt.
  13. Predloga:Sktxt.


Napaka pri navajanju: Za skupino »A« obstajajo značke <ref>, vendar ustrezne značke <references group="A"/> ni bilo mogoče najti.