330 likes | 500 Views
Numeriska ber äkningar i Naturvetenskap och Teknik. Dagens ämne:. Lite celest mekanik. F. Numeriska ber äkningar i Naturvetenskap och Teknik. Koordinatsystem. Kartesiska koordinater. Enhetsvektorerna är ortogonala och normerade. Cylinderkoordinater.
E N D
Numeriska beräkningar i Naturvetenskap och Teknik Dagens ämne: Lite celest mekanik F
Numeriska beräkningar i Naturvetenskap och Teknik Koordinatsystem Kartesiska koordinater Enhetsvektorerna är ortogonala och normerade
Cylinderkoordinater Numeriska beräkningar i Naturvetenskap och Teknik
Vektor- och skalärprodukt i cylinderkoordinater Numeriska beräkningar i Naturvetenskap och Teknik Ortogonala Högersystem
Sfäriska koordinater Numeriska beräkningar i Naturvetenskap och Teknik
Lite inledande mekanik Numeriska beräkningar i Naturvetenskap och Teknik Kraftlagen Momentet Rörelsemängdsmomentet ger:
Numeriska beräkningar i Naturvetenskap och Teknik Rörelsemängdsmomentet är konstant...
Numeriska beräkningar i Naturvetenskap och Teknik r x p är vinkelrät mot r, dvs r är vinkelrät mot L som är konstant. Centralkraft 1. Rörelsemängdsmomentet är en rörelsekonstant 2. Rörelsen sker i ett plan
Numeriska beräkningar i Naturvetenskap och Teknik För att sätta upp rörelseekvationerna behöver vi känna accelerationen i cylinderkoordinater.
Hastigheten i cylindriska koordinater Numeriska beräkningar i Naturvetenskap och Teknik Rörelse i planet givet av centralkraften Radiell hastighet vinkelhastighet
Accelerationen i cylindriska koordinater Numeriska beräkningar i Naturvetenskap och Teknik
Accelerationen i cylindriska koordinater Numeriska beräkningar i Naturvetenskap och Teknik
Accelerationen i cylindriska koordinater Numeriska beräkningar i Naturvetenskap och Teknik med ins. enl. ovan
Numeriska beräkningar i Naturvetenskap och Teknik Rörelsekvationerna i centralkraftsystemet med accelerationen i planet kan detta också skrivas:
Numeriska beräkningar i Naturvetenskap och Teknik Rörelsekvationerna i planet i cylinderkoordinater Beror av kraftens form Kan integreras utan att kraften specifieras Man utnyttjar nu följande trick... dvs vilket ger
Sektorhastigheten Numeriska beräkningar i Naturvetenskap och Teknik Keplers andra lag
Numeriska beräkningar i Naturvetenskap och Teknik Rörelsekvationerna i planet i cylinderkoordinater nu används men
Numeriska beräkningar i Naturvetenskap och Teknik Energin är en andra rörelsekonstant...
En andra rörelsekonstant Numeriska beräkningar i Naturvetenskap och Teknik För en konservativ kraft, dvs en kraft som har en potential Nytt trick... multiplicera med Detta är lika med
Fortsätt med att titta på Numeriska beräkningar i Naturvetenskap och Teknik vänsterledet i ekv nedan v.l. kan skrivas Vi har nu tidsderivator på båda sidor av denna ekvation! dvs
Numeriska beräkningar i Naturvetenskap och Teknik Hastigheten är Från L konstant har vi (fortfarande)
Lösningen till rörelsekvationerna Numeriska beräkningar i Naturvetenskap och Teknik Man kan nu antingen välja att försöka integrera lösningen i tidsvariabeln eller söka en lösning som funktion av vinkeln. Vi börjar med det senare fallet:
Lösningen till rörelsekvationerna Numeriska beräkningar i Naturvetenskap och Teknik I detta läge har man således men Binet!
Lösningen till rörelsekvationerna Numeriska beräkningar i Naturvetenskap och Teknik Binets ekvation för keplerfallet (1/r2 ) Andra ordningens diff ekv. (löses med den sekulära ekvationen!)
Olika typer av banor Numeriska beräkningar i Naturvetenskap och Teknik Referensriktning då α lika med noll
Numeriska beräkningar i Naturvetenskap och Teknik Olika typer av banor Undersöks på egen hand i projektet!
Banrörelse ρ(t) Numeriska beräkningar i Naturvetenskap och Teknik
Banrörelse ρ(t) Numeriska beräkningar i Naturvetenskap och Teknik Denna integral kan i princip lösas för t(ρ) men är inverteringen ρ(t) är inte möjlig i enkla funktioner. Samma sak gäller för vinkeln som funktion av tiden. Vad kan man göra?
Numeriska beräkningar i Naturvetenskap och Teknik Ytterligare ett variabel byte... Halva storaxeln Eccentriska anomalin Genomsnitts anomalin
Numeriska beräkningar i Naturvetenskap och Teknik efter detta variabelbyte... Keplers 3e lag (kan också fås genom geometrisk betraktelse)
Generellt vid tiden t Numeriska beräkningar i Naturvetenskap och Teknik Keplers ekvation Hur få ρ(t)? Endast numerisk lösning ger sedan ρ (detta var vår substitution)!
Numeriska beräkningar i Naturvetenskap och Teknik Tvåkropparsproblemet För två kroppar under ömsesidig vxv ersättes m med reducerade massan ovan: Trekropparsproblemet... Har lett till många försök till lösning (Poincare mfl). Det existerar serieutvecklingslösningar...Läs gärna själv historien bakom inkluderande ex.vis Mittag-Lefflers pris.
Notera att volymelementet i cylinderkoordinater är: Numeriska beräkningar i Naturvetenskap och Teknik