Transcript Document

Gravitace
Newtonův univerzální zákon: Každá částice ve vesmíru přitahuje každou jinou
částici silou, která je úměrná součinu jejich hmotností a nepřímo
úměrná druhé mocnině jejich vzdálenosti.
„jako v nebi tak i na Zemi“
Vzorcem: označíme jejich hmotnosti m1 a m2 jejich vzdálenost r
m1m2
Fg  G 2
r
G=6,67310-11Nm2kg-2
m2
Splňují zákon akce a reakce
Fg
Fg
m1
K tomu: síly od různých částic se sčítají (vektorově)…princip superpozice
+vliv síly na pohyb dá Newtonův základní zákon dynamiky:
dp
F
dt
● Vše, co potřebujeme znát pro určení pohybu těles vlivem gravitace.
● Matematický popis místo mechanického stroje.
● Vše ostatní jsou důsledky.
Hned tady na Zemi
Uvidíme odkud se bere zrychlení g:
Z minula:
Fg  mg
Nyní:
Fg  G
m MZ
rZ2
Porovnáním:
mg G
Odtud:
g G
Hrubý odhad:
MZ=5.97361024kg
rZ=6378km
m MZ
rZ2
MZ
rZ2
g  6 1011
6 1024
6 10 
6 2
ms-2  10ms-2
Navíc…skoro konstantní při změně r menší než poloměr Země
Keplerovy zákony pohybu planet
Logicky…důsledky Newtonova zákona
Historicky…Newton naopak z nich přišel na svůj zákon
1. Všechny planety se pohybují po elipsách
se Sluncem v jednom ohnisku.
2. Průvodič od Slunce k planetě pokryje
stejné plochy za stejné časové intervaly.
3. Druhá mocnina doby oběhu planety je úměrná
třetí mocnině hlavní poloosy elipsy.
1.zákon
a je hlavní poloosa
b je vedlejší poloosa
Pro vzdálenost ohniska c platí c2=a2-b2
Excentricita e = c/a
● 1. zákon je přímý důsledek závislosti síly na vzdálenosti jako 1/r2 (což nedokážeme).
● Trajektorie planet jsou skoro kružnice (nejvíc excentrický je Merkur…0,4%).
● Naopak trajektorie komet jsou velmi excentrické.
2. zákon
Za dobu dt urazí planeta dráhu vdt
Vektorový součin má velikost plochy
rovnoběžníka z obou vektorů.
Trojúhelník má poloviční plochu.
dA 
1
1
1
1
rdr  r vdt 
r  mv d t 
L dt
2
2
2m
2m
dA
1

L
d t 2m
Tímto jsme zavedli moment hybnosti L (více o něm za chvíli).
2. Zákon říká, že se L nemění…uvidíme brzy proč.
3. zákon
Ukážeme pro kruhovou dráhu, kdy hlavní poloosa je rovná poloměru r.
Doba oběhu T dá oběhovou rychlost v = 2r / T
Dostředivé zrychlení
(z minula)
v 2 2r / T 
4 2 r
a

 2
r
r
T
2. Newtonův zákon
Fg  M P a
2
M SM P
4 2 r
G
 MP 2
2
r
T
4 2 3
T 
r
GM S
2
Hmotnost planety se zkrátila mezi oběma stranami rovnice.
Gravitační pole
Gravitační síla působí mezi libovolnými dvěma objekty na libovolné vzdálenosti.
Možná interpretace: jeden „zdrojový“ objekt vytváří gravitační pole,
které působí na druhý „testovací“ objekt.
Pole působí ve stejném místě.
(místo kritizovaného působení na dálku od zdrojového objektu)
Působící síla je úměrná hmotnosti testovacího objektu. Síla vydělená hmotností
testovacího objektu = intenzita g. Charakterizuje pole samotné.
g
Fg
m

ma
a
m
Tedy intenzita = gravitační zrychlení (už jsme viděli blízko povrchu Země)
Pro pole od bodové částice:
GM
g 2
r
směřuje k zdrojové částici
Gravitační pole je konzervativní…už jsme se dozvěděli minule a ukázali
blízko povrchu Země, kde je gravitační zrychlení konstantní.
Teď ukážeme pro libovolnou vzdálenost:
Síla je radiální
(ve směru průvodiče)
 1 1
dr
W   F  d r   F r  d r  GMm  2  GMm   
r
 rf ri 
ri
ri
ri
rf
rf
rf
Takže práce závisí jen na počáteční a koncové vzdálenosti od zdroje, a tudíž
nezávisí na dráze.
Můžeme proto zavést potenciální energii vůči zvolenému referenčnímu bodu r0:
 1 1
d~
r
~
~
~
~
Epot r; r0    Fr   d r   F r  d r  GMm ~ 2  GMm  
r
 r0 r 
r
r
r
r0
r0
r0
Blízko povrchu Země
Ve výšce h je r = rZ + h
1
1 1
1 1

1 
1 
  GMm 
  GMm  1  h / rZ 
Epot r; r0   GMm 
 r0 rZ  h 
 r0 rZ 1  h / rZ 
 r0 rZ

Poslední úprava…
…rozvoj do řady v malém parametru h / rZ přičemž jsme nechali jen první člen.
Takže:
1 1
GM
Epot r; r0   GMm    m 2 h
rZ
 r0 rZ 
První člen = potenciální energie místa na povrchu Země
Druhý člen = potenciální energie vůči povrchu Země = mgh jako minule.
Vzorec pro Epot se zjednoduší, když jako referenční bod r0 zvolíme nekonečno.
Pro jednoduchost označení ho pak nebudeme psát:
Epot r   
GMm
r
Pak potenciální energie je vždycky záporná.
Celková energie…přidáme kinetickou, která je kladná
Pro kruhovou dráhu studovanou výše:
E  Ekin  Epot
1 2
mM 1 mM
mM
1 m M Epot
 mv G
 G
G
 G

0
2
r
2
r
r
2
r
2
Takže celková energie je záporná.
Volba referenčního bodu pro potenciální energii v nekonečnu pak znamená,
že objekt vázaný v gravitačním poli má celkovou energii zápornou.
Volný objekt má energii kladnou.
Energie rovná nule…na rozhraní mezi volným a vázaným
…objekt akorát tak unikne do nekonečna
0E
Musí mít únikovou rychlost
Hrubý odhad pro Zemi:
1 2
mM
mv G
2
r
v
2GM
r
v
2  6  1011  6  1024 1
7
1
1
ms

12

10
ms

10
kms
6  106
Únikové rychlosti pro různé planety
Důsledky pro atmosféru na planetě (pro nás):
● planety lehčí než Země…neudrží žádnou atmosféru.
● planety těžší než Země…v atmosféře převládá H2, He.
● u nás tak akorát…při pokojové teplotě převládá N2, O2.
Doposud jsme uvažovali jen bodové zdroje gravitačního pole i tam, kde to nebyla
pravda:
Gravitační zrychlení na povrchu Země jsme dostali, jako kdyby veškerá
hmota Země byla soustředěná v jejím středu.
Důvodem jsou slupkové teorémy o gravitačním poli tenké sférické slupky:
(1) Gravitační pole vně slupky je stejné, jako kdyby veškerá její hmota byla soustředěna
v jejím středu.
(2) Gravitační pole uvnitř slupky je nulové (gravitační Faradayova klec).
Kulatý zdroj (Slunce, Země) gravitačního pole si můžeme představit
jako složený ze sférických slupek.
Proto gravitační pole vně kulatého zdroje pole je stejné, jako kdyby
celá jeho hmotnost byla ve středu.
To věděl už Newton…teorém 30 a 31 v Principiích
Důkaz slupkových teorémů:
Slupkové teorémy plynou z obecnější vlastnosti gravitačního pole
…z Gaussova zákona: tok intenzity gravitačního pole ven z uzavřené
plochy je roven -4G hmotnost uvnitř plochy.
Znaménko mínus kvůli tomu, že gravitační intenzita směřuje ke zdroji
Nejprve proto dokážeme Gaussův zákon:
Začneme gravitačním polem od jednoho hmotného bodu hmotnosti M. Toto pole klesá
se vzdáleností od zdroje jako 1 / r2 (jako pokles intenzity záření od bodového zdroje).
Nejlépe je Gaussův zákon vidět pro sférickou plochu,
kde tok je roven součinu plochy a intenzity:
2větší poloměr…4větší plocha sféry
…ale též 4nižší intenzita, takže součin
plochy a intenzity se nezmění.
2
Tok = g 4r  
GM
4r 2  4GM
2
r
Obecnou plochu rozdělíme na kousky
Tok gravitační intenzity kouskem plochy…je potřeba vzít kolmý průmět plochy dA
k intenzitě g:
dA
r
dA
Pak:
d A  r 2 d Ω
d
M
kde d je malý prostorový úhel
Proto:
g r  d A  GM
1 2
r dΩ
2
r
r2 se zkrátí a integrál po celé ploše dá celkový prostorový úhel 4
 g r d A

 GM  d Ω  4GM
Kvůli superpozici tohle platí, když hmotu M libovolně rozložíme v prostoru
(to rozložení si můžeme představit jako složené z hmotných bodů
…též v druhé půlce přednášky)
Tím jsme dokázali Gaussův zákon. (potkáme jej taky v elektrostatice)
g(r)
Důkaz slupkových teorémů z Gaussova zákona:
Pro sférickou slupku hmotnosti M zvolíme jako uzavřenou plochu
sféru se středem ve středu slupky. Tato sféra může být větší nebo menší než slupka.
Ze symetrie má intenzita všude na sféře stejnou velikost a směr do středu
2


g
r
dA


4

r
g r 


g
● Sféra větší než slupka obsahuje uvnitř celkovou hmotnost
slupky M, takže
takže
 4r 2 g r   4GM
GM
g r   2
r
r
M
Jako pro pole od bodové částice
● Sféra menší než slupka nemá uvnitř žádnou hmotu, takže intenzita je nulová.
Tím jsme dokázali oba slupkové teorémy z Gaussova zákona.
Soustava hmotných bodů
nebo těleso, když nemůžeme zanedbat jeho rozměry (pak nekonečně mnoho bodů)
Uvidíme, že:
● Pohyb těžiště je stejný jako pohyb hmotného bodu.
● Rotační pohyb vůči těžišti se řídí podobnými zákony.
Nejprve zavedeme těžiště
Bod, ve kterém je v průměru veškerá hmota objektu.
Proto je těžiště přímo pod jakýmkoliv bodem zavěšení:
Matematicky:
rT 
m1r1  m2r2   m1r1  m2r2  

m1  m2  
M
● M je celková hmotnost soustavy
● Vážený průměr; váha každého bodu je jeho hmotnost.
● Pro nekonečně mnoho bodů v tělese přejde součet na integrál.
Druhou interpretaci těžiště dostaneme derivováním podle času:
d
d
d
Mv T  M
rT  m1 r1  m2 r2    p1  p 2    p celk
dt
dt
dt
Slovy: Rychlost těžiště  celková hmotnost je celková hybnost.
Tedy těžiště je bod, do kterého můžeme soustředit
celkovou hmotnost a celkový pohyb tělesa.
Sečteme všechny vnější síly. Podle Newtonova základního zákona dynamiky
dají celkovou změnu hybnosti.
Dosazením do rovnice nahoře:
d
d
Fvnejs, celk  p celk  M
v T  Ma T
dt
dt
1. věta impulsová
Fvnejs, celk  MaT
Těžiště se pohybuje po parabole, jako hmotný bod, v němž je
soustředěna veškerá hmotnost.
Co vnitřní síly mezi různými body soustavy?
Uvážíme nejprve jednu dvojici bodů. Podle 3. Newtonova zákona síly mezi nimi
tvoří pár akce a reakce. Proto změna hybnosti jednoho z obou bodů je stejně velká
a opačného směru vůči změně hybnosti druhého z bodů (viz kladivo a hřebík minule).
Příspěvek ke změně celkové hybnosti je proto nula.
Toto platí pro všechny dvojice bodů, takže vnitřní síly nemají vliv na pohyb těžiště.
Vnitřní síly nemají vliv na pohyb těžiště, ale mohou mít vliv na pohyb částí.
Nejlépe vidět ve srážkách (rozpadech) a reaktivním pohybu raket.
Srážka:
Těžiště si srážky ani nevšimne,
ale každá z koulí ano.
Hybnost se zachová při srážce vždy.
Pokud se zachová i mechanická energie, srážka je pružná. Jinak je nepružná.
Příklad: balistické kyvadlo
● Zařízení na měření rychlosti střely
● Nejprve srážka, pak vychýlení kyvadla
● Z výšky h můžeme určit počáteční rychlost.
Srážka: ze zákona zachování hybnosti dostaneme rychlost po srážce
m1v1,i  m1  m2 vf  vf 
m1v1,i
m1  m2
Energie:
2 2
m
1
1
1
m1
1
1 v1,i
2
2
Ef  m1  m2 vf 
 m1v1,i
 m1v12,i  Ei
2
2 m1  m2 2
m1  m2 2
Mechanická energie se nezachovává při srážce…srážka je nepružná.
Energie se zachová při následném pohybu dřeva se střelou
Ef  m1  m2 gh
m1  m2
v1,i 
2 gh
m1
Číselně: pro m1=5g, m2=1kg, h = 5cm je v1,i 
1.005kg
2 10ms- 2  0.05m  200ms-1
0.005kg
Reaktivní pohyb rakety
v+ v
Raketa se odstrkuje vypouštěním spálených plynů jako puštěný balonek.
M+m
v
● Rychlost výtoku paliva vůči lodi ve
● Změna hybnosti paliva o hmotnosti m je m ve
proti směru pohybu rakety.
● Změna hybnosti rakety o hmotnosti M je stejně velká
ve směru pohybu rakety.
● Změna hybnosti rakety je dána změnou rychlosti rakety
v jako M v.
● Vytékající palivo zmenšuje hmotnost rakety: m = - M.
● V limitě, kdy změny jdou k nule má zákon zachování
hybnosti tvar:
M
m
ve
M d v  ve d M
Toto je diferenciální rovnice pro rychlost v rakety jako funkci její hmotnosti M.
Rovnici řešíme separací proměnných a integrací:
M d v  ve d M
dM
d v   ve
M
vf
Mf
 d v  v 
e
vi
Mi
dM
M
Mi
vf  vi  ve ln
Mf
Potřebujeme co největší rychlost vypouštěných plynů a co největší část hmotnosti
lodě ve formě paliva.
Rotační pohyb vůči těžišti
Krasobruslař dá ruce k sobě a zvýší obrátky.
Množství rotačního pohybu dává moment hybnosti:
● Pro jeden hmotný bod:
L  r p
● Pro soustavu hm. bodů:
L celk  r1  p1  r2  p 2  
● Pro  bodů v tělese: součet přejde na integrál.
Vlastnosti L:
● Zvětší se při zvětšení hybnosti nebo vzdálenosti.
● Díky vlastnostem vektorového součinu se sčítá orientace
hybnosti vůči těžišti, ne hybnost.
Součet hybností vůči těžišti by dal nulu!
Časová změna pro jeden bod:
d
d
dr
dp
L  r  p  
p  r
 rF  τ
dt
dt
dt
dt
● První člen v derivaci součinu je nula, protože vektory jsou rovnoběžné.
● Pro druhý člen jsme použili Newtonova základního zákona.
● Poslední rovnost definuje moment síly .
Moment radiální síly je nulový…moment hybnosti se nemění…2.Keplerův zákon.
Planeta se chová jako ruce krasobruslaře.
Pro soustavu bodů postupujeme podobně jako u translačního pohybu:
● Sečteme všechny vnější momenty síly a momenty hybnosti.
● Páry akce-reakce vnitřních sil
nepřispívají ke změně celkového momentu hybnosti.
d
τ vnejs, celk  L celk
dt
2. věta impulsová
Tuhé těleso: vzájemné vzdálenosti
bodů se nemění
● Poloha každého bodu vůči těžišti daná úhlem otáčení kolem osy
● Časová derivace úhlu dá úhlovou rychlost
dθ
ω
dt
● Jako vektory…ve směru osy rotace podle pravidla pravé ruky (tady směrem k nám)
● Ještě jedna derivace podle času…úhlové zrychlení
dω
α
dt
● Rychlost daného místa v tělese daná úhlovou rychlostí pomocí vektorového součinu
vi  ω  ri
Celkový moment hybnosti soustavy hmotných bodů:
L   Li   ri  mi vi   miri ω  ri 
i
i
i
Pokud osa rotace je osou symetrie, dá se  vytknout a platí
L  ω mi ri 2  Jω
i
Veličina J se nazývá moment setrvačnosti.
Časová derivace:
dL
dω
J
 Jα
dt
dt
takže 2. Impulsová věta dostane tvar:
τ vnejs, celk  Jα
Kinetická energie:
1
1
1
1 2
2
2 2
Ekin   mi vi   mi ω  ri     mi ri   J
2 i
2
i 2
i 2

2
Moment setrvačnosti se chová jako hmotnost pro posuvný pohyb.
Příklad: moment setrvačnosti válce
Poloměr R, délka L, hustota 
Odtud hmotnost
M  R 2 L
Rozdělíme na slupky sestávající z bodů
stejné vzdálenosti r od osy a sečteme jejich
momenty setrvačnosti:
R
R
4
2
R
R
1
J   d m r 2  L  2r d r  r 2  2L  r 3 d r 2L
 R 2 L
 MR2
4
2 2
0
0
Podobnosti translačního a
rotačního pohybu
dr
dθ
v
ω
dt
dt
dv
dω
a
α
dt
dt
p  mv
L  Jω
dp
dL
F
 ma τ 
 Jα
dt
dt
rf
W   Fdr
ri
1 2
Ekin  m v
2
θf
W   τ dθ
θi
1
Ekin  J 2
2
Příklad: válec na podložce
…translační a rotační pohyb současně
Válec o hmotnosti m = 1kg a poloměru r = 10cm je roztočen kolem své podélné
osy úhlovou rychlostí 0 = 120 s-1 a položen na drsnou vodorovnou podložku.
Účinkem smykového tření, jehož koeficient má hodnotu  = 0,2, se válec
uvede do zrychleného posuvného pohybu a současně začne být brzděn ve
svém pohybu otáčivém. Za jakou dobu válec přejde do čistého valivého pohybu
bez prokluzování?
0
m
T

(t)>0
T

r
v(t)
F=mg
Čas = 0
F=mg
Čas = t
Rovnice pro posuvný pohyb těžiště:
Odtud zrychlení:
Integrace dá:
d vt 
F m
dt
d vt  F mg
 
 g
dt
m
m
vt    g d t  g  1d t  gt  C1
Integrační konstantu C1 určíme z počáteční podmínky v(t = 0)=0
g 0  C1  0  C1  0
Takže rychlost těžiště jako funkce času je:
vt   gt
Rovnice pro otáčivý pohyb kolem těžiště:
d  t 
dt
d  t 
F r
m gr
2g

 1 2 
dt
J
r
2 mr
Odtud úhlové zrychlení:
Integrace dá:
F  r  J
2g
 2g 
 2g 
t  C2
dt  
  1d t  
r 
r 
r


 t     
Integrační konstantu C2 určíme z počáteční podmínky (t = 0)=  0
2g
0  
0  C 2  C 2  0
r
Takže úhlová rychlost jako funkce času je:
2 g
 t   0 
t
r
Podmínka pro valivý pohyb bez prokluzování:
Dosazení za v(t) a (t) dá:
Odtud doba:
Číselně:


gt  r  0 
r0
t
3g
101 m 120s 1
t
 2s
2
3  0.2 10ms
vt   rt 
2g 
t
r 
I minule (volný pád) i nyní jsme vyřešili diferenciální rovnici z Newtonova zákona
pro konstantní sílu nebo konstantní moment síly.
Příště…oscilace…síla závislá na poloze
Navíc oscilace jsou všude, kde je stabilní rovnováha (viz minule).
Když se oscilace mohou šířit v prostoru, vzniká vlnění.