Kepler yörüngesi

Kepler Orbitand eksentriklerini çeşitli biçimlerde gösteren diyagram. Mavi bir hiperbolik yörünge (e> 1) 'dir. Yeşil bir parabolik yörünge (e = 1) 'dir. Kırmızı eliptik yörünge (0 <e <1). Gri bir dairesel yörünge (e = 0) 'dir. (e = 0).

Gök mekaniği olarak, Kepler yörüngesi (veya Keplersel yörünge) üç boyutlu uzayda iki boyutlu bir yörünge düzlemi oluşturan bir elips, parabol, hiperbol benzeri bir yörünge cismininin hareketini açıklar. (Kepler yörüngesi aynı zamanda düz bir çizgi de çizebilir). Kepler yörüngesi yalnızca nokta iki cismin nokta benzeri yerçekimsel çekimlerini dikkate alır, atmosfer sürüklemesi, güneş radyasyonu baskısı, dairesel olmayan cisim merkezi ve bunun gibi bir takım şeylerin diğer cisimlerle girdiği çekim ilişkileri nedeniyle ihmal eder. Böylece Kepler problemi olarak bilinen iki-cisim probleminin, özel durumlara bir çözüm olarak atfedilir. Klasik mekaniğin bir teorisi olarak, aynı zamanda genel görelilik etkilerini dikkate almaz. Kepler yörüngeleri çeşitli şekillerde altı yörünge unsurları içine parametrize edilebilir.

Çoğu uygulamalarda, tüm sistemin kütle merkezi olarak kabul edilen büyük bir merkezi cisim ve kütle merkezi vardır. Ayrıştırma yapılınca, iki objenin benzer kütleleri, ikisinin ortak kütle merkezinin Kepler yörüngesi ya da barisentileri olarak açıklanabilir.

Giriş

Eski çağlardan itibaren 16. ve 17. yüzyıllara kadar gezegenlerin hareketlerinin Antik Yunan filozofları Aristo ve Batlamyus tarafından öğretildiği gibi mükemmel dairesel jeosentrik yolları takip ettiğine inanılırdı. Gezegenlerin hareketler değişimleri (epidaire) üst üste küçük dairesel yollar şeklinde açıklanmıştır.. Gezegenlerin ölçümleri giderek daha kesin olmaya başladığında, bu teoriye revizyonlar önerilmiştir. Hala gezegenlerin güneş merkezli mükemmel dairesel yörüngede gittiğine inanılmasına rağmen 1543 yılında Nicolaus Copernicus, güneş sisteminin bir güneş merkezli modelini yayınladı.

Yasaların geliştirilmesi

1601 yılında Johannes Kepler, Tycho Brahe tarafından yapılan gezegenlerin kapsamlı ve titiz gözlemlerini elde etti. Kepler, sonraki beş yılı Mars gezegeninin gözlemlerini çeşitli eğrilere uydurmaya çalışarak geçirecekti. 1609 yılında Kepler, gezegensel hareketin üç yasasından ilk ikisini yayınladı. Birinci yasa şöyledir:

"Güneş sistemindeki bütün gezegenler, odaklarının birinde Güneş olan elips şeklindeki bir yörüngede hareket ederler."

Daha genel anlatımla, Kepler hareketini yapan bir nesnenin yolu elips olduğu kadar, parabol veya hiperbol de olabilir. Bu bir grup şekil konik kesitler olarak bilinir. Matematiksel olarak, merkezi bir cisim ile yörüngedeki bir cisim arasındaki mesafe şu şekilde ifade edilebilir:

$r(\theta )=\frac{a(1-e^2)}{1+e\cos (\theta )}$

• $r$ mesafedir
• $a$ yörüngenin boyutunu tanımlayan yarı büyük eksendir
• $e$ yörüngenin şeklini tanımlayan dış merkezliktir
• $\theta$ yörüngedeki nesnenin mevcut konumu ile yörüngedeki merkez cisme en yakın olduğu (enberi noktası olarak adlandırılan) konum arasındaki açı olan gerçek ayrıklıkdır.

Alternatif olarak, denklem şu şekilde ifade edilebilir:

$r(\theta )=\frac{p}{1+e\cos (\theta )}$

Burada $p$, eğrinin yarı özkirişi olarak adlandırılır. Denklemin bu formu, özellikle yarı büyük ekseni sonsuz olan parabolik yörüngelerle uğraşırken kullanışlıdır.

Bu yasaları gözlemleriyle geliştiren Kepler, hiçbir zaman bu hareketleri açıklayacak bir teori geliştirememiştir.

Isaac Newton

Isaac Newton, 1665 ve 1666 yılları arasında hareket, kütleçekim ve diferansiyel hesap ile ilgili çeşitli kavramlar geliştirdi. Bununla birlikte bu kavramlar, 1687 yılında hareket yasalarını ve evrensel kütleçekim yasasını ana hatlarıyla belirttiği Principia'ya kadar yayınlanmadı. Üç hareket yasasından ikincisi şöyle der:

Bir cismin ivmesi, cisme etki eden net kuvvetle paralel ve doğru orantılıdır. Net kuvvet yönündedir ve cismin kütlesi ile ters orantılıdır:

$𝐅=m𝐚=m\frac{d^2𝐫}{d{t}^2}$

• $𝐅$ kuvvet vektörüdür
• $m$ kuvvetin etki ettiği cismin kütlesidir
• $𝐚$ ivme vektörüdür, $𝐫$ konum vektörünün ikinci zaman türevidir.

Kurallara bakılırsa, denklemin bu şekli yalnızca aşağıda yapılan basitleştirici varsayımlara dayalı olarak doğru olan sabit kütleli bir nesne için geçerlidir.

Newton'un evrensel kütleçekim yasasının mekanizmaları: Bir nokta kütle olan m₁ diğer bir nokta kütle olan m₂'yi, iki kütlenin çarpımı ile orantılı ve aralarındaki (r) uzaklığının karesi ile ters orantılı olan bir F₂ kuvveti ile çeker. Kütleler veya mesafeden bağımsız olarak, |F₁| ve |F₂|'nin büyüklükleri her zaman eşit olacaktır. G kütle çekimi sabitidir.

Newton'un kütleçekim yasası şöyle der:

Her nokta kütle, her iki noktayı kesen doğru boyunca işaret edilen bir kuvvetle diğer nokta kütleleri çeker. Kuvvet, iki kütlenin çarpımı ile orantılı ve nokta kütleler arasındaki uzaklığın karesi ile ters orantılıdır:

$F=G\frac{m_1{m}_2}{r^2}$

• $F$ iki nokta kütlesi arasındaki kütleçekim kuvvetinin büyüklüğüdür.
• $G$ kütle çekimi sabitidir
• $m_1$ birinci nokta kütlesinin kütlesidir
• $m_2$ ikinci nokta kütlesinin kütlesidir
• $r$ iki nokta kütlesi arasındaki mesafedir

Newton, hareket yasaları ve evrensel kütleçekim yasasından astronomide yörünge hareketine özgü Kepler yasalarını türetebildi. Kepler yasaları gözlem verileriyle iyi bir şekilde desteklendiğinden, bu tutarlılık Newton'un genelleştirilmiş teorisinin ve birleşik göksel ve sıradan mekaniğin geçerliliğine güçlü bir destek sağladı. Bu hareket yasaları, 20. yüzyılın başlarında Albert Einstein özel ve genel görelilik kavramlarını ortaya koyana kadar modern gök mekaniğinin temelini oluşturdu. Çoğu kullanımda Kepler hareketi, gezegenlerin ve uyduların hareketlerini nispeten yüksek doğruluk derecelerine yaklaştırır ve astronomi ve astrodinamikte yaygın olarak kullanılır.

Sadeleştirilmiş iki cisim problemi

İki cisimli bir sistemde bir cismin hareketini çözmek için iki sadeleştirici varsayım yapılabilir:

1. Cisimler küresel olarak simetriktir ve nokta kütleler olarak ele alınabilir.
2. Cisimlere karşılıklı kütleçekimi dışında etki eden iç ve dış kuvvetler yoktur.

Büyük gök cisimlerinin şekilleri kürelere yakındır. Simetri ile, bir kütle noktasını homojen bir küreye doğru çeken net kütleçekim kuvveti, merkeze doğru yönlenmelidir. Kabuk teoremi (Isaac Newton tarafından da kanıtlanmış olan), kürenin yoğunluğu derinliğe göre değişse bile (çoğu gök cisimleri için olduğu gibi) bu kuvvetin büyüklüğünün, kütlenin tamamı kürenin ortasında toplanmış gibi aynı olduğunu belirtir. Bundan hemen sonra iki homojen küre arasındaki çekim, her ikisinin de kütlesi merkezinde yoğunlaşmış gibidir.

Asteroitler veya uzay araçları gibi daha küçük nesneler genellikle bir küreden güçlü bir şekilde ayrılan bir şekle sahiptir. Ancak bu düzensizlikler tarafından üretilen kütleçekim kuvvetleri, merkezi cismin kütleçekimi ile karşılaştırıldığında genellikle önemsizdir. Düzensiz bir şekil ile mükemmel bir küre arasındaki fark da mesafelerle azalır ve yörünge mesafelerinin çoğu, yörüngedeki küçük bir cismin çapıyla karşılaştırıldığında çok büyüktür. Bu nedenle bazı durumlar için şekil düzensizliği, kesinlik üzerinde önemli bir etki olmaksızın ihmal edilebilir. Bu etki, özellikle alçak yörüngelerde bulunan yapay Dünya uyduları için oldukça belirgindir.

Gezegenler değişen hızlarda dönerler ve bu merkezkaç kuvveti nedeniyle hafif basık bir şekil alabilirler. Böyle yassı bir şekille kütle çekimi, homojen bir küreninkinden bir miktar sapacaktır. Daha büyük mesafelerde bu yassılığın etkisi ihmal edilebilir hale gelir. Güneş Sistemindeki gezegensel hareketler, nokta kütleler olarak ele alınırsa yeterli hassasiyetle hesaplanabilirler.

Kütleleri $m_1$ ve $m_2$ olan iki nokta kütle nesnesi ve bazı eylemsiz referans çerçevesine göreceli konum vektörleri ${𝐫}_1$ ve ${𝐫}_2$, kütleçekim kuvvetlerine maruz kalır:

$m_1{\ddot{𝐫}}_1=\frac{-G{m}_1{m}_2}{r^2}\widehat{𝐫}$

$m_2{\ddot{𝐫}}_2=\frac{G{m}_1{m}_2}{r^2}\widehat{𝐫}$

burada $𝐫$, kütle 1'in kütle 2'ye göre göreceli konum vektörüdür ve şu şekilde ifade edilir:

$𝐫={𝐫}_1-{𝐫}_2$

$\widehat{𝐫}$, o yöndeki birim vektördür ve $r$, bu vektörün uzunluğudur.

Kendi kütlelerine bölmek ve ikinci denklemi birinciden çıkarmak, birinci cismin ikinciye göre ivmesi için hareket denklemini verir:

$\alpha$, yerçekimi parametresidir ve eşittir

$\alpha =G(m_1+m_2)$

Birçok uygulamada, üçüncü bir sadeleştirici varsayım yapılabilir:

3. Merkezi cisim ile karşılaştırıldığında yörüngedeki cismin kütlesi önemsizdir. Matematiksel olarak, m₁ >> m₂, yani Şablon:Math'dir. Genellikle $\mu =GM$ olarak gösterilen bu tür standart kütleçekim parametreleri, yörüngedeki uydularından çok daha büyük kütlelere $M$ sahip olan Güneş, büyük gezegenler ve Ay için yaygın olarak kullanılır.

Bu varsayım, sadeleştirilmiş iki cisim problemini çözmek için gerekli değildir, ancak özellikle Dünya yörüngesindeki uydular ve Güneş'in yörüngesindeki gezegenlerle ilgili hesaplamaları basitleştirir. Jüpiter'in kütlesi bile Güneş'inkinden 1047 kat daha azdır, bu da α değerinde %0,096'lık bir hata oluşturur. Dikkate değer istisnalar arasında Dünya-Ay sistemi (kütle oranı 81,3), Plüton-Charon sistemi (kütle oranı 8,9) ve ikili yıldız sistemleri sayılabilir.

Bu varsayımlar altında iki cisim problemi için diferansiyel denklem matematiksel olarak tamamen çözülebilir ve Kepler'in gezegensel hareket yasalarını takip ederek sonuçta oluşan yörüngeye "Kepler yörüngesi" denir. Tüm gezegenlerin yörüngeleri yüksek doğrulukta Güneş etrafındaki Kepler yörüngeleridir. Küçük sapmalar, gezegenler arasındaki çok daha zayıf kütleçekim etkilerinden ve Merkür durumunda genel görelilikten kaynaklanmaktadır. Dünya etrafındaki yapay uyduların yörüngeleri, makul bir tahminle Güneş, Ay ve Dünya'nın yassılığından dolayı küçük düzensizliklere sahip Kepler yörüngeleridir. Hareket denkleminin bütün kütleçekim ve kütleçekim olmayan kuvvetlerin (güneş radyasyon basıncı ve atmosfer direnci gibi) hesaba katılmasıyla sayısal olarak entegre edilmesi gereken yüksek doğruluktaki kullanımlarda, Kepler yörünge kavramları çok önemlidir ve yoğun olarak kullanılır.

Kepler yörünge öğeleri

Kepler yörünge öğeleri

Kepler yörüngesinin altı parametre ile tanımlanabilirliğinden söz etmek gerekir. Üç boyutlu uzayda hareket eden bir nesnenin hareket bir konum vektörü ve bir hız vektörü ile karakterize edilir. Her vektörün üç bileşeni vardır, bu yüzden uzayda bir yörünge tanımlamak için gerekli değerlerin toplam sayısı altıdır. Bir yörünge (Kepler Elementleri olarak bilinen) ve pozisyon ve hızına bağlı olarak ölçülebilen altı elementten oluşur ki bunların üçü daha önce gösterilmiştir. Bu altı uygun elementten beş tanesi yörünge için değişmezdir. (iki sürekli değişen vektöre bir tezat niteliğinde olmak üzere). Kendi yörüngesinde içindeki bir nesnenin gelecekteki konumu tahmin edilebilir ve yeni konum ve hız kolayca yörünge elemanları ile elde edilebilir.

Bunlardan ikisi eksenin boyutunu ve şeklin belirler:

• Yarımajör eksen ($a$)
• eksantriklik ($e$)

Üç tanesi yörünge alanının oryantasyonunu belirler:

• eğim ($i$)yörünge düzlemi ile referans düzlemi arasındaki açıyı tanımlar.
• Yükselen nodun uzunluğu ($\mathrm{\Omega }$) Referans yönü ve referans düzlemi üzerinde yörüngeye yukarı kesişme (yükselen nod) arasındaki açıyı tanımlar.
• Periapsis argüman($\omega$) artan nod ve periapsis arasındaki açıyı belirler.

Ve son olarak:

• Gerçek anomali ($\nu$)periapsis ölçülen yörünge boyunca yörüngedeki cismin konumunu tanımlar. Çeşitli başka değerlerin yerine gerçek bir anomali kullanılabilir, bunlardan en sık kullanılanları mean anomoli olan $M$ ve periapsisden yana geçen süre olan $T$'dir.

$i$, $\mathrm{\Omega }$ ve $\omega$ sadece yörünge düzlemi içinde nesnenin hareketini tartışırken referans çerçevesinde yörünge yönünü tanımlayan açısal ölçümlerdir ve kesin olarak gerekli değildir. Burada tamamlanmışlık adına belitilmişlerdir ancak aşağıdaki ispatlar için gerekli değillerdir.

Yukarıdaki diferansiyel denklemin (Şablon:Denklem notu) matematiksel çözümü

Herhangi bir merkez kuvveti altında yani r bir kuvvet paraleli ve hareket için, belirli bağıl açısal momentumu $𝐇=𝐫\times \dot{𝐫}$ sabit kalır:
$\dot{𝐇}=\frac{d}{dt}(𝐫\times \dot{𝐫})=\dot{𝐫}\times \dot{𝐫}+𝐫\times \ddot{𝐫}=\mathrm{𝟎}+\mathrm{𝟎}=\mathrm{𝟎}$

Pozisyonun çapraz ürünü vektör ve onun hızı olduğu için sabit kaldığından, aynı düzelmde bulunmak zorundadırlar. ($𝐇$'ye ortogonal). Bu vektör, fonksiyonun, bir düzlem eğrisi olduğu anlamına gelir. Denklemin kökeni etrafında simetri vardır, çünkü kutupsal koordinatlarda çözmek kolaydır. Bununla birlikte, bu denklem dikkat etmek önemlidir Şablon:Denklem notu doğrusal ivmeye refere eder $\left(\ddot{𝐫}\right)$, açısalın karşıtı olarak $\left(\ddot{\theta }\right)$ yahut radyal $\left(\ddot{r}\right)$ ivmelenme. Bu nedenle, denklemi değiştiren kişinin dikkatli olması gerekir.Şimdi kartezyen koordinat sistemine $(\widehat{𝐱},\widehat{𝐲})$ ve kutupsal birim vektörlere $(\widehat{𝐫},\hat{\mathit{\theta }})$ düzlemi octoganelinde bakalım $𝐇$:

$\widehat{𝐫}=\cos (\theta )\widehat{𝐱}+\sin (\theta )\widehat{𝐲}$
$\hat{\mathit{\theta }}=-\sin (\theta )\widehat{𝐱}+\cos (\theta )\widehat{𝐲}$

Şimdi vektör $𝐫$ fonksiyonunu yeniden yazabiliriz ve diferansiyali şöyledir as:

$𝐫=r(\cos \theta \widehat{𝐱}+\sin \theta \widehat{𝐲})=r\widehat{𝐫}$

$\dot{𝐫}=\dot{r}\widehat{𝐫}+r\dot{\theta }\hat{\mathit{\theta }}$

$\ddot{𝐫}=(\ddot{r}-r{\dot{\theta }}^2)\widehat{𝐫}+(r\ddot{\theta }+2\dot{r}\dot{\theta })\hat{\mathit{\theta }}$

(Vektör calculusunu inceleyiz). Bunları yerine yazarsak şunu buluruz:
(Şablon:Denklem notu)

$(\ddot{r}-r{\dot{\theta }}^2)\widehat{𝐫}+(r\ddot{\theta }+2\dot{r}\dot{\theta })\hat{\mathit{\theta }}=(-\frac{\mu }{r^2})\widehat{𝐫}+(0)\hat{\mathit{\theta }}$

Bu sıradan olmayan polar diferansiyel denklem verir:

Bu denklemi çözmek için, öncelikle her zaman diferansiyel lineerini ortadan kaldırmak gerekir. Şuna ulaşırız:

$H=|𝐫\times \dot{𝐫}|=|(r\cos (\theta ),r\sin (\theta ),0)\times (\dot{r}\cos (\theta )-r\sin (\theta )\dot{\theta },\dot{r}\sin (\theta )+r\cos (\theta )\dot{\theta },0)|=|(0,0,r^2\dot{\theta })|=r^2\dot{\theta }$ Şablon:NumBlk

(Şablon:Denklem notu)'ün zaman diferansiyalinii aldığımız zaman, şuna ulaşırız:

(Şablon:Denklem notu) ve (Şablon:Denklem kaynağı}) denklemleri 'ın diferansiyellerini elememize izin verir.. In order to eliminate the time derivatives of $\theta$'ın zaman diferansiyellerini ellimine etmek için uygun kısaltmaları bulmak bunun için de zincir kuralı kullanmamız gerekir.

$\ddot{r}-r{\dot{\theta }}^2=-\frac{\mu }{r^2}$

$\frac{d^{2}r}{d{\theta }^2}\cdot {\dot{\theta }}^2+\frac{dr}{d\theta }\cdot \ddot{\theta }-r{\dot{\theta }}^2=-\frac{\mu }{r^2}$

$\frac{d^{2}r}{d{\theta }^2}\cdot {\left(\frac{H}{r^2}\right)}^2+\frac{dr}{d\theta }\cdot (-\frac{2\cdot H\cdot \dot{r}}{r^3})-r{\left(\frac{H}{r^2}\right)}^2=-\frac{\mu }{r^2}$ Şablon:NumBlk Diferansiyel denklem (Şablon:Denklem notu) değişken değişimi ile analitik çözülebilir Şablon:NumBlk Bu verilenleri kullanarak, in (Şablon:Denklem notu)'deki tüm zaman diferansiyelleri yoksayılabilir ve bizi $r\theta$.
$\ddot{r}-r{\dot{\theta }}^2=-\frac{\mu }{r^2}$ Şablon:NumBlk

(Şablon:Denklem notu) ve (Şablon:Denklem notu)'da belirtilenleri $\frac{d^{2}r}{d{\theta }^2}$ ve $\frac{dr}{d\theta }$ için kullanarak

Şablon:NumBlk genel bir sonuca ulaşılır. Şablon:NumBlk

e ve ${\theta }_0$'un başlangıçtaki s ve $\frac{ds}{d\theta }$'n başlangıç değerlerinin integral sabitleri olursa, integralinin sabitini kullanmak yerine birim vektörleri kullanıma alınır ve seçilen yörünge alanındaki koordinat sistemini belirler, böylece sıfır değerini alır ve e is pozitiftir. Bu şu anlama gelir: 'in maksimum olduğu yerde sıfırdır ve böylece minimumdur. p parametresini $\frac{H^2}{\mu }$ olarak tanımlarsak şu ortaya çıkar.

Alternatif türevi

Bu denklemi polar diferansiyel denklemini kullanmadan çözmenin bir diğer yolu şöyledir:
$𝐮$'yu bir birim vektör olarak tanımlayın, örneğin,$𝐫=r𝐮$ and $\ddot{𝐫}=-\frac{\mu }{r^2}𝐮$ gibi. Burdan yola çıkarak

$𝐇=𝐫\times \dot{𝐫}=r𝐮\times \frac{d}{dt}(r𝐮)=r𝐮\times (r\dot{𝐮}+\dot{r}𝐮)=r^2(𝐮\times \dot{𝐮})+r\dot{r}(𝐮\times 𝐮)=r^2𝐮\times \dot{𝐮}$

Şimdi şunu değerlendiriniz

$\ddot{𝐫}\times 𝐇=-\frac{\mu }{r^2}𝐮\times (r^2𝐮\times \dot{𝐮})=-\mu 𝐮\times (𝐮\times \dot{𝐮})=-\mu [(𝐮\cdot \dot{𝐮})𝐮-(𝐮\cdot 𝐮)\dot{𝐮}]$

(Üçlü vektör ürünü). Şunu dikkate alınız

$𝐮\cdot 𝐮=|𝐮{|}^2=1$

$𝐮\cdot \dot{𝐮}=\frac{1}{2}(𝐮\cdot \dot{𝐮}+\dot{𝐮}\cdot 𝐮)=\frac{1}{2}\frac{d}{dt}(𝐮\cdot 𝐮)=0$

Bu verileri bir önceki denkleme yerleştirdiğimizde

$\ddot{𝐫}\times 𝐇=\mu \dot{𝐮}$

İki tarafın da integralini alırsak

$\dot{𝐫}\times 𝐇=\mu 𝐮+𝐜$

Burada c sabit vektördür. Bunu r' ile birleştirmek ortaya ilginç bir sonuç çıkarır

$𝐫\cdot (\dot{𝐫}\times 𝐇)=𝐫\cdot (\mu 𝐮+𝐜)=\mu 𝐫\cdot 𝐮+𝐫\cdot 𝐜=\mu r(𝐮\cdot 𝐮)+rc\cos (\theta )=r(\mu +c\cos (\theta ))$

Burada $\theta$ $\overline{r}$ ve $\overline{c}$ arasındaki açıdır. r'ye göre çözersek

$r=\frac{𝐫\cdot (\dot{𝐫}\times 𝐇)}{\mu +c\cos (\theta )}=\frac{(𝐫\times \dot{𝐫})\cdot 𝐇}{\mu +c\cos (\theta )}=\frac{|𝐇{|}^2}{\mu +c\cos (\theta )}$

Dikkat ediniz ki $(r,\theta )$ vektör fonksiyonunun polar koordinatlarıdır. Verilenleri yerine koyarsak ve, şu denkleme ulaşırız

Bu fokal noktada merkezi olan bir konik bölmenin polar koordinatlarının denklemidir. $\theta$ Argümanına gerçek anomali denir.

Yörünge denkleminin özellikleri

$e=0$ ise bu çapı p olan bir dairedir.

$0<e<1$ ise bu aşağıdakine sahip bir elipstir.

$e=1$ ise fokal uzunluğu $\frac{p}{2}$ olan bir dairedir. $e>1$ ise bir aşağıdakine sahip bir hiperboldür.

Sonraki görsel, bir elipsi (kırmızı), bir paraboolü (yeşil) ve bir hiperbolü (mavi) gösterir.

Dış merkezi 0.7 olan bir elips Kepler yörüngesi, Parabol Kepler yörüngesi ve dış merkezi 1.3 olan bir Kepler yörüngesi. Fokal noktaya olan uzaklık, denklemle de gösterildiği gibi polar açının yatay çizgiye yakın açısının bir fonksiyonudur(Şablon:Denklem notu)

Yatay çizgi üzerinde nokta odak noktasından sağa doğru çıkıyorsa nokta $\theta =0$ odak noktasına olan minimum uzaklık $\frac{p}{1+e}$ değerini alır (pericentre) .Elips için de odak uzaklığı maksimum değerini aldığı bir apocentre vardır$\frac{p}{1-e}$. Hiperbol için $\theta$ 'nın aralığı:

$[-{\cos }^{-1}(-\frac{1}{e})<\theta <{\cos }^{-1}(-\frac{1}{e})]$

Bir parabol için aralık;

$[-\pi <\theta <\pi ]$

Diferansiyelin zinci kuralını kullanarak (Şablon:Denklem notu), (Şablon:Denklem notu) denklemi ve p’nin tanımı $\frac{H^2}{\mu }$ olmak üzere radyal hız elemanı şudur.

ve teğetsel bileşen olan (hız bileşeni $V_r$)’ye dik )

Polar argümanı $\theta$ ve zaman t arasındaki bağlantı elips ve hiperbolik yörüngeler için küçük değişiklikler gösterebilir.

Elips şeklindeki bir yörünge için, dış merkezli anomali olan E’ye geçeriz ve burdan

Şablon:NumBlk Şablon:NumBlk ve sonuç olarak. Şablon:NumBlk Şablon:NumBlk

açısal momentum H olmak üzere Şablon:NumBlk

Zamana göre entegrasyon olan t ile şu sonuca ulaşılır.

zaman varsayımını yaparak integral sabiti 0 olur. p nin tanımı gereği Şablon:NumBlk şu şekilde yazılabiir Şablon:NumBlk

Hiperbolik bir yörünge için parametre yaratmak için hiperbolik fonksiyon kullanılır.

bu durumda şuna ulaşılır

açısal momentum H olmak üzere

Zamana göre entegrasyon olan t ile şu sonuca ulaşılır.

Şablon:NumBlk örneğin Şablon:NumBlk

t’nin hangi zamanda belirli bir gerçek anomaliye ($\theta$) ne ulaştığını bulmak için, buna uyan parametre olan E’nin (Şablon:Denklem notu) ile zamana olan bağını hesaplaması, elips için (Şablon:Denklem notu) ile bağını ve hiperbolik yörünge için.

Dikkat ediniz ki (Şablon:Denklem notu) ve (Şablon:Denklem notu) bu aralıkta bir yer belirler. $[-\mathrm{\infty }<t<\mathrm{\infty }]⟷[-\mathrm{\infty }<E<\mathrm{\infty }]$

Bazı ilave formüller

Ayrıca bakınız: Merkez denklemi - analitik açılımlar

(Şablon:Denklem notu)ve (Şablon:Denklem notu)’den edilinen elips yörünge için Şablon:NumBlk ve buna göre Şablon:NumBlk (Şablon:Denklem kaynağı)’dan devam ettiğimizde şu sonuç çıkar

${\tan }^2\frac{\theta }{2}=\frac{1-\cos \theta }{1+\cos \theta }=\frac{1-\frac{\cos E-e}{1-e\cdot \cos E}}{1+\frac{\cos E-e}{1-e\cdot \cos E}}=\frac{1-e\cdot \cos E-\cos E+e}{1-e\cdot \cos E+\cos E-e}=\frac{1+e}{1-e}\cdot \frac{1-\cos E}{1+\cos E}=\frac{1+e}{1-e}\cdot {\tan }^2\frac{E}{2}$

Dış merkezli anomalinin geometrik yapısını belirlemek için, açıkça vektörler$(\cos E,\sin E)$ ve $(\cos \theta ,\sin \theta )$‘nin x aksisiyle aynı yanda olmaları gerekir. Buradan yola çıkarak $(\cos \frac{E}{2},\sin \frac{E}{2})$ ve $(\cos \frac{\theta }{2},\sin \frac{\theta }{2})$ aynı çeyrek daireye sahiptir. Burdan yola çıkarak

ve

$(x,y)$ vektörünün polar argümanlarınınve n "$\arg (x,y)$" örneğinde olduğu gibi seçilmiş olması gerekir.

$\arg (x,y)$nın sayısal hesaplaması için, standart ATAN2(y,x) fonksiyonunun, ( ya ikili DATAN2(y,x)), FORTRAN gibi programlama dillerinde kullanılması mümkündür.

Dikkate alınız ki bu aralıklar arasında alan vardır.

$[-\mathrm{\infty }<\theta <\mathrm{\infty }]⟷[-\mathrm{\infty }<E<\mathrm{\infty }]$

Hiperbolik bir yörünge için (Şablon:Denklem notu) ve (Şablon:Denklem notu)dan yola çıkarak Şablon:NumBlk and therefore that Şablon:NumBlk

ve sonuç olarak

${\tan }^2\frac{\theta }{2}=\frac{1-\cos \theta }{1+\cos \theta }=\frac{1-\frac{e-\cosh E}{e\cdot \cosh E-1}}{1+\frac{e-\cosh E}{e\cdot \cosh E-1}}=\frac{e\cdot \cosh E-e+\cosh E}{e\cdot \cosh E+e-\cosh E}=\frac{e+1}{e-1}\cdot \frac{\cosh E-1}{\cosh E+1}=\frac{e+1}{e-1}\cdot {\tanh }^2\frac{E}{2}$

$\tan \frac{\theta }{2}$ ve $\tanh \frac{E}{2}$ aşağıdaki ifadeyi sağlayan aynı işarete sahiptir. Şablon:NumBlk Bu ilişki gerçek anomali ve E parametresi arasında geçiş yapmak için uygundur. (burada E (Şablon:Denklem notu)’e zaman yönünden bağlıdır.) Aralık arasındaki alana dikkat ediniz.

$[-{\cos }^{-1}(-\frac{1}{e})<\theta <{\cos }^{-1}(-\frac{1}{e})]⟷[-\mathrm{\infty }<E<\mathrm{\infty }]$

ve aşağıdaki bağla açıklanabilir.

${\tanh }^{-1}x=\frac{1}{2}\ln \left(\frac{1+x}{1-x}\right)$

(Şablon:Denklem notu)’ye bakarak, P’nin elips bir yörünge için yörünge periyodu Şablon:NumBlk

Kuvvet alanına potansiyel enerjinin ilgisi (Şablon:Denklem notu)

$-\frac{\mu }{r}$

(Şablon:Denklem notu), (Şablon:Denklem notu), (Şablon:Denklem notu) ve (Şablon:Denklem notu)’ kadar kinetik ve potansiyel enerjinin özeti

$\frac{{V_r}^2+{V_t}^2}{2}-\frac{\mu }{r}$

Elips yörünge için

Şablon:NumBlk (Şablon:Denklem notu), (Şablon:Denklem notu), (Şablon:Denklem notu) ve (Şablon:Denklem notu)’den sonra kinetik ve potansiyel enerjinin hiperbol yörünge için

$\frac{{V_r}^2+{V_t}^2}{2}-\frac{\mu }{r}$

Eylemsizlik koordinat sistemiyle ilişkisi

Şablon:NumBlk and from (Şablon:Denklem notu), (Şablon:Denklem notu), olan bir yörünge düzleminde pericentre’ye doğru (Şablon:Denklem notu) ve (Şablon:Denklem notu)’dan yapılan hız elementleri çıkarımı Şablon:NumBlk

Belirli bir başlangıç durumuna karşılık Kepler yörüngesinin belirlenmesi

$\hat{x},\hat{y}$

Bu diferansiyel denklemi (Şablon:Denklem notu) için "başlangıç değeri problemi" dir ki 6 boyutlu bir durum vektründe ilk denklemdir ve şöyle yazılır

Başlangıç değeri vektörünün herhangi bir değeri için, bu başlangıç değeri probleminin çözümüne uygun Kepler Yörüngesi aşağıdaki algoritmayla bulunur

Ortogonal birim vektörleri $(\hat{r},\hat{t})$ tanımalamak gerekirse

Şablon:NumBlk Şablon:NumBlk ve $r>0$ and $V_t>0$

(Şablon:Denklem notu), (Şablon:Denklem notu), (Şablon:Denklem notu) ve (Şablon:Denklem notu)’dan sonra sırasını takip eder.

1. Şablon:NumBlk

$e\ge 0$ and $\theta$‘yı tanımlarken

Şablon:NumBlk Şablon:NumBlk

aşağıdakinin olduğu yerde

1. Şablon:NumBlk

aynı r, $V_r$ and $V_t$ değerlerine sahip olan gerçek anomaliye sahip bir Kepler yörüngesi bulunur ve bu (Şablon:Denklem notu) ve (Şablon:Denklem notu)’de gösterilmiştir.

Eğer bu Kepler yörüngesi aynı vektörlerine sahipse ve gerçek anomalisi (Şablon:Denklem kaynağı) ve (Şablon:Denklem kaynağı)de tanımladığı gibi ise, Kepler yörüngesinin $(\hat{r},\hat{t})$ vektörleri istenen değer olan $(\overline{r_0},\overline{v_0})$’yi gerçek anomalisi $\theta$ olmak üzere alır.

Standart atalet sabit koordinat sistemi $(\hat{x},\hat{y})$ yörüngesel düzlemde ($\hat{x}$ pericentre'nin homojen küresindin merkezinden gelecek şekilde) konik alanın oryantasyonunu belirliyorsa (elips, parabol, hiperbol), bu ilişkiyle hesaplanabilir.

1. Şablon:NumBlk
2. Şablon:NumBlk

(Şablon:Denklem notu) ve (Şablon:Denklem notu)’ün $V_r=0$ olduğu durumda bağlantılı olduğuna dikkat ediniz ve

$V_t=V_0=\sqrt{\frac{\mu }{p}}=\sqrt{\frac{\mu }{\frac{{(r\cdot V_t)}^2}{\mu }}}$

örneğin

1. Şablon:NumBlk

Diaresel bir yörüngeye uyan başlangıç durumunu $(\overline{r_0},\overline{v_0})$ gösterir.

Oskülatör kepler yörüngesi

$(\overline{r},\overline{v})$ gibi herhangi bir durum vektörü için, bu durumdaki Kepler yörüngesi aşağıda tanılmlanan algoritma ile hesaplanabilir. İlk parametreler $p,e,\theta$ ile belirlenir ve $\hat{x},\hat{y}$ yörüngesel düzelemde ortogonal birim vektörleri de (Şablon:Denklem notu) ve (Şablon:Denklem notu) arasındaki ilişkiyi kullanarak bulunur.

Şimdi, eğer hareket denklemi aşağıdaki gibi ise

1. Şablon:NumBlk

aşağıdakinin

$\overline{\mathrm{F}}(\overline{r},\dot{\overline{r}},t)$

bir fonksyion olup

$-\mu \cdot \frac{\hat{r}}{r^2}$

aşağıdaki parametlere sahip olmaması gerekir.

$p,e,\theta ,\hat{x},\hat{y}$

$\overline{r},\dot{\overline{r}}$ tarafından tanımlanmış, zaman içinde değişecek olan durumlar, Keppler Yörüngesinin tersine yalnızca $\theta$’da değişecektir.

T zamanında bu şekilde "hareket denkleminin" çözümü olarak aynı "durum vektör" olan (Şablon:Denklem kaynağı) 'da hesaplanan Kepler yörüngesi bu durumda "oskülatör" olduğu söylenebilir. Bu konsept, şu durumda yararlıdır.

$\overline{\mathrm{F}}(\overline{r},\dot{\overline{r}},t)=-\mu \cdot \frac{\hat{r}}{r^2}+\overline{\mathrm{f}}(\overline{r},\dot{\overline{r}},t)$

olmak üzere

$\overline{\mathrm{f}}(\overline{r},\dot{\overline{r}},t)$

küçük bir “bozucu güç”tür çünkü, bir örnekle açıklamak gerekirse diğer gök cisimlerinden sönük bir yer çekimi bu durumda örnek gösterilebilir. Oskülasyon halindeki Kepler yörüngesinin parametreleri, ancak o zaman yavaşça değişecek ve oskülasyon halindeki Kepler yörüngesi iyi bir tahminle oskülasyonun öncesinde ve sonrasında gerçek yörüngeye hatrı sayılır bir zaman diliminde girecektir. Bu konsept aynı zamanda bir roket uçuşu için de faydalı olabilir çünkü itme kapatıldığı zaman hangi Keppler yörüngesinde devam edebileceği belirlenmiş olur. “Daire olmaya yakın” bir yörünge konsepti için, dışmerkezli vektörü şeklinde tanımlamak yararlıdır. (Şablon:Denklem notu), (Şablon:Denklem notu) ve (Şablon:Denklem notu)’dan yola çıkarak:

Örneğin $(\overline{r},\overline{v})$ durum vektörlerinin diferansiyel fonskyonlarıdır $\overline{e}$, ayrıca bu durum dairesel bir yörüngeye karşılık geliyor ise de geçerlidir.

Ayrıca bakınız

• İki cisim problemi
• Kütleçekimsel iki cisim problemi
• Kepler problemi
• Kepler'in gezegensel hareket yasaları
• Eliptik yörünge
• Hiperbolik yörünge
• Kaçış yörüngesi

Kaynakça

• El'Yasberg "Theory of flight of artificial earth satellites", Israel program for Scientific Translations (1967)
• Şablon:Kitap kaynağı
• Şablon:Kaynak