astronomi: Juli 2011

Bab 2

Orbit dalam ruang

Pada paragraf yang lalu telah diuraikan bahwa suatu lintasan di dalam ruang ditentukan oleh bentuk orbit dan orientasinya. Bentuk geometri suatu orbit dicerminkan oleh oleh dua unsur yaitu elemen geometri dan elemen orientasi. Elemen orientasi adalah sudut simpul naik, W, argumen perihelium w dan inklinasi i . Sedangkan elemen geometri ialah setengah sumbu panjang elip, a , eksentrisitas, e .

Gb. 2-1 Orbit anggota Tata Surya relatif terhadap bidang ekliptika dengan Matahari sebagai salah satu titik api lintasan berbentuk Elip

Periode orbit, P dan saat terakhir melewati titik terdekat dengan titik fokus lintasannya yang berbentuk elip,T, disebut elemen dinamik. Seandainya kala edar P diketahui maka masalah yang harus dipecahkan adalah bagaimana menyatakan koordinat polar benda langit sebagai fungsi waktu. Dari pengetahuan ini kita akan dapat menentukan posisi benda langit tersebut dalam koordinat ekuatorial, asensio rekta, a dan deklinasi, d. Untuk keperluan ini tinjaulah ilustrasi yang diragakan dalam Gb. 2-2

Gb.2-2 Ilustrasi orbit elip dan lintasan bantu Kepler

(lingkaran putus-putus dengan jejari a)

Andaikan m adalah Satelit yang bergerak mengorbit Bumi, m₁ dan misalkan pula koordinat polar titik massa m pada saat t adalah (r,f). Dalam hal ini r, menyatakan jarak m terhadap m₁ sedangkan f, adalah sudut yang dibentuk oleh radius vektor r terhadap sumbu referensi yang kita pilih. Selanjutnya definisikan besaran berikut;

a) Anomali benar (true anomaly) f, adalah sudut yang diukur searah dengan gerak titik perige terhadap garis vektor yang menghubungkan m dengan m₁

b) Anomali eksentrik (eccentric anomaly) E, yaitu sudut pada pusat lingkaran yang diukur dari perige dalam arah yang sama seperti halnya f

c) Anomali rata-rata (mean anomaly) M, dinyatakan sebagai sudut yang ditempuh oleh radius vektor r, rata-rata selama satu satuan waktu sejak meninggalkan titik perige.

(2-1)

Harga n dapat ditentukan dari kaedah hukum Kepler III yaitu ;

(2-2)

Dalam hal ini : T = saat terakhir melewati perige

k = konstanta Gauss

m = dinyatakan dalam massa matahari

n = dalam radian persatuan waktu

2.1 Pernyataan persamaan lintasan

Untuk membahas persamaan lintasan akan digunakan bantuan geometri seperti yang diperlihatkan dalam Gb. 2-2.

Kita lihat bahwa;

S₂K : LK = b : a atau r Sin f : a Sin E = b : a (2-3)

Terlihat pula bahwa;

S₁K = r Cos f = ae – a Cos E atau r Cos f = a Cos E –ae (2-4)

dan S₂K = r Sin f = a

(2-5)

Dari kedua pernyatan ini dapat dihitung bahwa;

(r Cos f)² +(r Sin f)² = a( 1 – e² ) Sin²E + ( aCos E – ae )² (2-6)

Atau r = a ( 1 – eCos E ) (2-7)

Dengan mengingat hubungan goniometri, Cos f = 1- 2 Sin² (f/2) dengan cara yang sama kita peroleh;

Atau

(2-8)

Berdasarkan hukum Kepler III dapat juga diturunkan dengan cara berikut; yaitu pada saat T, m ada di titik terdekat dengan massa m₁ selanjutnya terlihat pula bahwa;

Luas S₁S₂P =

Disamping itu luas S₁S₂P dapat juga dihitung dengan cara yang lain yakni;

Luas S₁S₂P = Luas KPS₂ + Luas S₁S₂K

, atau dapat ditulis

M = E –e Sin E (2-9)

Persamaan ini disebut dengan persamaan Kepler. Nilai E dapat dihitung dari persamaan Kepler bila M dan e diketahui. Bila eksentrisitas, e cukup kecil, dalam hal ini e < 0,2 seperti halnya orbit planet dan mayoritas satelit buatan, persamaan ini dapat diuraikan dalam deret Fourier yang bentuknya dinyatakan oleh persamaan berikut;

(2-10)

Dalam hal ini J_k adalah fungsi Bessel, contoh untuk k= 2 adalah;

(2-11)

Untuk keperluan praktis akan lebih mudah kalau persamaan Kepler diselesaikan dengan metoda numerik Newton-Raphson. Caranya diberikan dalam algoritma dan flowchart berikut;

Algoritma Newton-Raphson(f(e),f’(E),E₀,e, M dan E)

1. Berikan nilai pemula E₀ untuk harga E

2. Hitung

3. Test apakah, | E-E₀ | £ e | E₀ | bila ya proses dihentikan dan E adalah nilai yang memenuhi. Bila tidak ambil E₀ = E kembali ke langkah 2

Simbol pada algoritma diatas adalah f(E) persamaan Kepler dan f’(E) turunan pertamanya. Sedangkan e adalah presesi yang kita inginkan dan E₀ harga pendekatan awal anomali eksentrik yang kita ambil. Bila E telah dapat ditentukan maka r dan f dapat kita hitung dari persamaan (2-7) dan (2-6).

Sewaktu menggunakan algoritma ini perlu diperhatikan nilai pemula E₀ . Perlu dihindari titik stasioner yaitu titik dimana turunan pertama fungsi Kepler, f’(E₀) = 0 dan titik belok, yaitu titik dimana terjadi peralihan dari cekung ke atas ke cekung ke bawah atau sebaliknya. Titik belok memenuhi syarat f“(E₀) = 0, iterasi tidak akan pernah konvergen pada kedua titik ini. Untuk itu dalam program komputer perlu dibuat subroutine guna menghindari kedua kasus ini. Flowchart pada Gb 2-.3 tidak meninjau kasus seperti ini.

Flowchart Metoda Newton-Raphson

Gb.2-3 Flowchart solusi persamaan Kepler. Dalam hal proses tidak konvergen

ulangi proses dengan mengambil harga E₀ yang berbeda

2.2 Contoh Kasus

Sebuah titik massa m, dilempar dari planet Bumi dengan tujuan Mars. Andaikan dalam perjalanannya ke planet Mars, titik massa itu hanya dipengaruhi oleh gaya gravitasi Matahari. Tentukankanlah kordinat (r,f) titik massa tersebut bila diketahui jarak Mars dari Matahari pada saat titik massa m dilemparkan adalah 1,38 SA .

Penyelesaian

Deskripsi persoalan ini dijelaskan dengan diagram bantu seperti yang diragakan dalam Gb.2-4

Gb.2-4 Diagram lintasan Mars, titik massa m dan orbit Bumi

Perhatikan gambar diatas, untuk titik m berlaku;

Jarak Bumi-Mars pada saat itu merupakan sumbu panjang lintasan elip yang akan ditempuhnya, jadi 2a = 1 + 1,38 = 2,38 SA dengn demikian a = 1,19 SA
Periode P dapat dicari dari hukum Kepler III;

= 0,0172(1,19) ^-1,5 = 0,0132 rad/hari

sehinggga P = 2p/n = 476 hari

Menentukan eksentrisitas e;

Posisi perihelion r_p = a(1- e ) = 1 SA

Posisi aphelion r_a = a (1 + e) = 1,38 SA

Dengan demikian kita peroleh eksentrisitas e = 0,16 sehingga radius vektornya dapat dihitung dari ;

Untuk menentukan radius vektor r dan anomali benar f pada saat tertentu, kita harus mengetahui lebih dahulu posisi Bumi dan Mars pada saat yang bersangkutan konfigurasi umum ini diragakan pada gambar 2-15

Gb.2-5 Konfigurasi planet Mars dan Bumi. Jarak Mars dari Bumi dapat

dihitung dengan rumus kosinus ;

Sebagai contoh, misalkan kita ingin mendaratkan wahana pada saat jarak Bumi dan Mars minimum yaitu sekitar tangggal 20 Agustus 1961. Jadi titik massa m harus kita luncurkan P/2 atau 238 hari sebelum 20 Agustus 1961 jadi tanggal 26 Desember 1960. Anomali benar dan jarak wahana dari Matahari untuk berbagai tanggal diberikan dalam Tabel 2-3 berikut

Tabel 2-1 Jarak wahana dan anomali benar untuk berbagai saat pengamatan

No	Tanggal [1961]	t-T (hari)	M (Radian)	f (derajad)	r (SA)
4.	Februari, 1	35	0,463	36,4	1,03
5.	Maret ,1	64	0,847	64,6	1,09
6.	April, 1	95	1,255	90,4	1,16
7.	May, 1	125	1,655	112,6	1,24
8.	Juni, 1	156	2,062	132,6	1,30
9.	Juli, 1	186	2,460	150,8	1,35
10.	Agustus, 1	217	2,872	168,6	1,37

Jika posisi wahana diketahui maka jarak wahana dari Bumi bisa dihitung, untuk itu paragraf berikut menjelaskan cara menghitung koordinat ekuatorial wahana tersebut

Gb.2-6 Posisi m dalam sistem kartesis XYZ. m₁ menyatakan matahari dan m, menunjukkan wahana

Dari gambar diatas kita lihat bahwa ;

x = r cos b cos l

y = r cos b sin l (2-12)

z = r sin b

atau dalam sistem koordinat yang baru dimana sumbu x’ dambil sebagai garis nodal, maka dapat dilihat bahwa;

x’ = r Cos b Cos ( l - W ) = r Cos ( f + w )

y’ = r Cos b Sin ( l - W ) = r Sin ( f + w ) Cos I (2-13)

z’ = r Sin b = r Sin ( f + w ) Sin i

Oleh sebab itu jika r,b,l dan W diketahui maka x’,y’ dan z’ bisa dihitung. Selanjutnya dari pernyataan ini dapat diturunkan beberapa hal yaitu; menentukan hubungan koordinat ekuatorial heliosentrik dan elemen posisi wahana

Gb.2-7 Lintasan titik massa m dalam ruang. Sumbu x mengarah

pada titik vernal ekuinok.(posisi matahari terbit tanggal 21 Maret)

Dari pernyataan (2-13) dapat dilakukan beberapa kombinasi bila ini dilakukan maka dari pernyataan diatas kita lihat bahwa;

Tan ( l - W ) = Tan ( f + w ) Cos i (2-14)

Tan b = Sin ( l - w ) Tan i (2-15)

Pernyataan ini menunjukkan bahwa bila inklinasi, i = p/2 maka Tan ( l - W ) =0 atau l = W , sedangkan b tidak dapat didefinisikan, demikian pula halnya apabila Sin ( l - w ) = 0 maka berakibat b =0.

Jika m₁ menyatakan Matahari dan kita ingin menentukan posisi m dalam tata koordinat ekuatorial, maka kedudukan m dengan koordinat (b, l) bila dilihat dari Bumi merupakan posisi m dalam koordinat ekuatorial heliosentrik. Untuk menentukan (a,d) bila dilihat dari Bumi dapat dilakukan dengan melakukan transformasi koordinat ekuatorial heliosentrik ke koordinat ekuatorial geosentrik;

Gb.2-8 Konversi posisi ekuatorial heliosentrik ke

tata koordinat ekuatorial geosentrik

Dalam gambar diatas p menyatakan planet atau benda langit lainnya sedangkan koordinat ekuatorial, ( a, d ) = ( l, l ) menyatakan kedudukan titik massa m bila dilihat dari Matahari. Selanjutnya E, menyatakan Bumi sebagai titik asal koordinat jadi posisinya adalah (0,0,0) dan S menyatakan Matahari dengan koordinat ( X, Y, Z ) dapat dilihat pada Nautical Almanac atau dihitung dengan menggunakan algoritma Meeus (1997). Andaikan bidang x-h adalah bidang ekuator Bumi, maka Matahari akan mempunyai koordinat (X,Y,Z). Kedudukan relatif P terhadap S adalah;

x = x’ + X = r Cos d Cos a

h = h’ + Y = r Cos d Sin a (2-16)

z = z’ + Z = r Sin d

Akibatnya kita mempunyai ;

Tan a =

dan Sin d = x ( x² + h²+ z² ) ^-1/2 (2-17)

Dengan demikian a dan d dapat kita tentukan.

2-3 Menentukan Elemen Orbit

Menghitung orbit benda langit identik dengan menentukan elemen orbitnya yaitu: a, e, i, W, w dan T. Karena ada enam konstanta yang harus dihitung maka paling sedikit harus ada tiga pasang data pengamatan mengenai a dan d sebagai fungsi waktu. Misalkan ( a₁ , d₁) , ( a₂ , d₂) menyatakan posisi ekuatorial geosentrik planet tersebut pada saat t₁ dan t₂ dan ( l, b ) adalah longitude dan latitude planet tersebut. Jarak planet ke Bumi dinyatakan sebagai r maka dari pengetahuan tentang transformasi koordinat yang telah dipelajari pada Astronomi Bola, dapat ditunjukkan bahwa;

Sin b = Sin d Cos e - Cos d Sin e Sin a

Cos b Sin l = Sin d Sin e - Cos d Cos e Sin a (2-18)

Cos b Cos l = Cos d Cos a

Dalam hal ini untuk perhitungan yang tidak memerlukan ketelitian tinggi dapat diambil nilai e = 23⁰27’ dengan demikian dari pernyataan diatas dapat dihitung ( l₁, b₁ ) dan ( l₂, b₂ ), yaitu nilai ( l, b ) pada saat t₁ dan t₂. Maka koordinat siku-siku ekliptika geosentrik adalah;

Untuk waktu t₁

x₁ = r₁ Cos b₁ Cos l₁

y₁ = r₁ Cos b₁ Sin l₁

z₁ = r₁ Sin b₁

Untuk waktu t₂

x₂ = r₂ Cos b₂ Cos l₂

y₂ = r₂ Cos b₂ Sin l₂ (2-19)

z₂ = r₂ Sin b₂

Karena bidang orbit Bumi identik dengan bidang ekuator Matahari maka dapat dianggap latitude Matahari , B @ 0

Untuk waktu t₁

X₁ = R₁ Cos L₁

Y₁ = R₁ Sin L₁

Z₁ = 0

Untuk waktu t₂

X₂ = R₂ Cos L₂

Y₂ = R₂ Sin L₂ (2-20)

Z₂ = 0

Dalam hal ini R dan L masing-masing adalah jarak Bumi-Matahari dan longitude geosentrik Matahari. Koordinat ini (X_,Y, Z ) dapat dilihat pada Nautical Almanac untuk setiap waktu t. Untuk lebih jelas perhatikanlah gambar 2-8 dengan S,P,E dan r masing-masing menyatakan Matahari, Planet, Bumi dan jarak matahari ke Planet.

Gb.2-9 Konversi posisi ekliptika helionsentrik ke sistem ekliptika geosentrik

Dalam sistem baru ini l dan b adalah longitud planet P dan ini adalah koordinat heliosentrik P. Kemudian jika X₀,Y₀,Z₀ menyatakan koordinat kartesis P didalam sistem heliosentrik, maka kita mempunyai;

x₀ = x – X = r Cos b Cos l

y₀ = y – Y = r Cos b Sin l (2-21)

z₀ = z – Z = r Sin b

Pada saat t₁ kita dapatkan;

x₀(t₁) = x₁ – X₁ = r₁ Cos b₁ Cos l₁

y₀(t₁) = y₁ – Y₁ = r₁ Cos b₁ Sin l₁ (2-22)

z₀(t₁) = z₁ – Z₁ = r₁ Sin b₁

Sedangkan untuk t₂ kita peroleh;

x₀(t₂) = x₂ – X₂ = r₂ Cos b₂ Cos l₂

y₀(t₂) = y₂ – Y₂ = r₂ Cos b₂ Sin l₂ (2-23)

z₀(t₂) = z₂ – Z₂ = r₂ Sin b₂

Substitusi persamaan (2-20) dan (2-21) pada pernyataan diatas maka kita peroleh;

r₁ Cos b₁ Cos l₁ = r₁ Cos b₁ Cos l₁ – R₁ Cos L₁

r₁ Cos b₁ Sin l₁ = r₁ Cos b₁ Sin l₁ – R₁ Sin L₁ (2-24)

r₁ Sin l₁ = r₁ Sin b₁

dan ;

r₂ Cos b₂ Cos l₂ = r₂ Cos b₂ Cos l₂ – R₂ Cos L₂

r₂ Cos b₂ Sin l₂ = r₂ Cos b₂ Sin l₂ – R₂ Sin L₂ (2-25)

r₂ Sin b₂ = r₂ Sin b₂

dalam pernyataan ini harga R₁ , R₂ , L₁, L₂ , atau X₁ , Y₁, X₂ , Y₂ dapat dilihat pada Nautical Almanac untuk saat t₁ dan t₂ . Selanjutnya nilai l dan b dapat kita hitung. Andaikan lintasan planet mengelilingi Matahari dalam bentuk lingkaran, dengan perkataan lain r₁ = r₂ jadi hanya satu besaran r yang perlu ditentukan. Selanjutnya r, l dan b untuk t₁ dan t₂ dapat kita tentukan. Dengan mengambil kuadrat pernyataan r₂ Cos b₂ Sin l₂ dan r₂ Cos b₂ Cos l₂ dari persamaan (2-25) dan kemudian menjumlahkannya diperoleh;

r₂² = r₂² + R₂² -2 r₂ R₂Cos b₂ Cos ( l₂ –l₂) (2-26)

tetapi r₁ = r₂ akibatnya jika r₁ diketahui maka r₂ dapat dihitung dengan demikian pernyataan (2-23) dapat digunakan untk mencari l₂ dan b₂ . demikian pula jika r₁ dapat ditaksir, r₁ dapat ditentukan dengan begitu l₁ dan b₁ dapat dihitung. Sekarang kita harus melihat bagaimana besaran ini dapat dipergunakan untuk menentukan elemen orbit. Dalam gambar 2-9, misalkan P₁ dan P₂ menyatakan posisi planet pada saat t₁ dan t₂ .

Gb.2-10 Segitiga bola dan bidang ekliptika

Dengan menggunakan hukum kosinus untuk segitiga bola P₁NP₂ kita mempunyai hubungan;

Cos A = Sin b₁ Sin b₂ + Cos b₁ Cos b₂ Cos (l₂ – l₁) (2-27)

Dalam hal ini ;

A adalah busur lingkaran yang ditempuh planet dalam interval waktu (t₂ – t₁). Jika koordinat (l₁ , b₁ ) dan ( l₂ , b₂ ) diketahui maka A dapat dihitung dari hukum Kepler III;

, P dinyatakan dalam satuan hari

Massa planet dapat diabaikan karena ia jauh lebih kecil dari massa Matahari

maka ;

(2-28)

Busur A ditempuh dalam waktu;

(2-29)

Jadi dengan perkataan lain luas busur tempuhannya adalah;

(2-30)

Nilai A yang dihitung dengan persamaan (2-29 ) haruslah sesuai dengan pernyataan (2-30) dan ini hanya berlaku bila ( l₁, b₁ ) dan ( l₂, b₂ ) menunjukkan hasil yang benar. Jadi (2-29) dan (2-30) dapat kita gunakan untuk menentukan ( l₁, b₁ ) dan ( l₂, b₂ ) dengan cara iterasi numerik. Prosedurnya sebagai berikut;

Algoritma ( r₀ , t_i , l_I , b_I, R_i , L_i ) i= 1,2

Berikan harga r₀ sembarang pada r₁
Tentukan ( r₁ , l₁ , b₁ ) dari pernyataan (2-24)
Dari pernyataan (2-26) hitung r₂ dalam hal ini r₂ = r₁
Gunakan (2-25) untuk menghitung l₂ dan b₂

Untuk menentukan i, w, dan W perhatikanlah segitiga bola yang diragakan pada Gb.2-12 berikut.

Gb. 2-11 Aplikasi rumus Napier dalam segitiga bola untuk menghitung elemen orbit dan analoginya pada hubungan i, w, dan W suatu lintasan pada segitiga bola

Dengan kaedah Napier untuk saat t₁ kita memperoleh;

Sin ( l₁ - W ) = tan (90⁰ – I ) Tan b₁ (2-21)

Hal yang sama untuk t₂;

Sin ( l₂ - W ) = tan (90⁰ – I ) Tan b₂ (2-22)

Selanjutnya gunakan cara berikut;

( l₂ - W ) = ( l₁ - W )+ ( l₂ – l₁ ), maka pernyataan (2-230) dapat ditulis sebagai;

Tan i { Sin ( l₁ - W ) Cos ( l₂ – l₁ ) + Cos ( l₁ - W ) Sin( l₂ – l₁ )} = Tan b₂ (2-23)

Gabungkan (2-23) pada (2-22) maka kita peroleh hasilnya;

Tan b₁ Cos ( l₁ - W ) + Tan I Cos( l₁ - W ) Sin( l₂ – l₁ )} = Tan b₂ (2-24)

Substitusi (2-21) pada (2-25) maka kita peroleh hasil sebagai berikut;

Tan b₁Cos( l₂ – l₁ ) +

Tan I Cos ( l₁ - W ) Sin( l₂ – l₁ )} = Tan b₂ atau dapat juga ditulis sebagai;

(2-25)

Dengan menggunakan kembali pada pernyataan (2-21) ataupun pada (2-22 ) nilai i dapat kita hitung. Harga W dapat dicari dari pernyataan ( lihat gambar 2-21 )

Sin W Cos i = Sin( l₂ – W) Cos b₁– Cos ( l₁ – W) Sin b₁ Cos 90⁰

= Sin ( l₁ – W) Cos b₁

Dan

Cos w = Cos( l₁ – W) Cos b₁ + Sin ( l₁ – W) Sin b₁ Cos 90⁰

= Cos ( l₁ – W) Cos b₁

Dari kedua persamaan diatas diperoleh;

(2-26)

Perlu diiingat bahwa w dan ( l₁ – W) berada dalam kuadran yang sama. Jadi dengan proses diatas bila P diketahui maka i, W dan w dapat ditentukan. Selanjutnya tinjaulah kasus jika radius vektor sebuah objek diketahui pada tiga posisi di langit untuk watu yang berbeda. Maka elemen lintasan dapat kita tentukan dengan cara berikut. Misalkan benda langit bergerak mengitari Matahari dengan lintasan elip, periodenya P andaikan pula posisi koordinat heliosentrik ekliptika diberikan oleh pernyataan;

pada saat t₁

pada saat t₂ (2-27)

pada saat t₃

Selanjutnya kita andaikan bahwa koordinat polar (l,b,r) pada tiga saat tersebut diketahui dengan demikian koordinat kartesis

dan

pada saat itu dapat ditentukan. Misalkan N₁, N₂ dan N₃ vektor satuan seperti yang diperlihatkan pada gambar 2-9

vektor satuan pada garis nodal dengan arah ke titik simpul naik

vektor satuan yang tegak lurus pada

dan terletak pada bidang orbit

vektor satuan yang tegak lurus

dan

jadi

Dapat dilihat bahwa;

= ( Cos W )

+( Sin W )

= ( Cos i Sin W )

+( Cos i Cos W )

+( Sin i )

(2-28)

= ( Sin i Sin W )

- (Cos W Sin i)

+( Cos i )

Dari vektor

dan

cari vektor satuan yang tegak lurus

dan

dengan cara sebagai berikut;

A₁

+A₂

+A₃

(2-29)

Vektor ini tegak lurus terhadap bidang orbit dan identik dengan

oleh sebab itu dapat ditulis;

A₁ = ( Sin i Sin W )

A₂ = - (Cos W Sin i) (2-30)

A₃ = ( Cos i )

Dari sini kita peroleh;

i = Arc Cos (A₃ ) dan W = Arc Tg

(2-31)

Jadi bila

dan

diketahui maka i, W dapat kita hitung dan hanya berlaku bila

dan

non-collinear, persamaan lintasan dapat ditulis kembali sebagai;

(2-32)

Dalam hal ini f adalah anmali benar dan u sudut yang dibentuk dari titik simpul naik ke radius vektor pada bidang orbit, selanjutnya kita lihat bahwa bila dinyatakan dalam besaran skalar maka;

e Cos w(r₁ Cos u₁ – r₂ Cos u₂ ) + e Sin w(r₁ Sin u₁ – r₂ Sin u₂ ) = r₂ – r₁

e Cos w(r₁ Cos u₁ – r₃ Cos u₃ ) + e Sin w(r₁ Sin u₁ – r₃ Sin u₃ ) = r₃ – r₁ (2-33)

Tetapi r Cos u =

dan r Sin u =

(2-34)

Oleh sebab itu persamaan (2-33 ) dapat ditulis sebagai;

(2-35)

Karena W dan i diketahui maka sistem persamaan linier ini dapat diselesaikan bila bentuk ( e Cos w ) dan ( e Sin w ) telah ditentukan.

Misalkan diketahui ( e Cos w ) = a₁ dan ( e Sin w ) = a₂ maka kita peroleh;

dan

(2-36)

Karena e > 0, kuadran w ditentukan oleh tanda aljabar dari besaran ( e Cos w ) dan besaran ( e Sin w ). Dari persamaan (2-240 ) setengah sumbu panjang elip dapat kita tentukan, yaitu;

(2-37)

Setiap harga r yang dipergunakan harus memberikan hal yang sama, karena tadi kita andaikan periode P diketahui maka dengan menggunakan kaedah hukum Kepler III, a dapat ditentukan, demikian pula sebaliknya bila P tidak diketahui maka a harus dihitung lebih dahulu. Saat melewati perihelion dapat dicari dengan bantuan pernyataan; r = a(1 – eCosE) dan persamaan Kepler; M=E-eSinE, dalam hal ini M dapat ditentukan pada setiap saat pengamatan. Harga T dapat diperoleh dari bentuk;

(2-38)

Untuk mencari elemen orbit a, e, i, W, w dan T sebenarnya hanya diperlukan dua posisi dalam koordinat polar, dengan menggunakan konstanta luas dan persamaan (2- ) nilai a,e dan w dapat diturunkan. Berikut diberikan sebuah contoh cara menentukan elemen orbit dari suatu benda langit.

2-4 Ilustrasi

Koordinat heliosentrik sebuah objek yang bergerak diberikan oleh tabel 2-1 dibawah ini. Dari pengamatan diketahui gerak harian rata-rata objek adalah, 4⁰,0923/hari. Pertanyaannya tentukanlah elemen orbit benda langit tersebut

Tabel 2-2 Posisi polar objek pada tahun 1960

No	Tanggal	l	b	r(SA)
1	Juni, 1, 0^h UT	142⁰ 45’ 40”	6⁰ 58’43”	0,34200
2	Juni 6, o^h UT	166⁰ 37’ 56”	6⁰ 08’47”	0,37022
3	Juni 11, 0^h UT	186⁰ 58’ 43”	4⁰ 35’51”	0,39867

Dari pernyataan (2-22 ) transformasi koordinat polar ke koordinat kartesis kita peroleh harga x,y dan z untuk ketiga data pengamatan tersebut;

Tabel 2-3 Posisi kartesis objek pada tahun 1960

No	Tanggal	R(SA)	x	y	z
1	Juni, 1, 0^h UT	0,34200	-0,27025	0,20542	0,0415
2	Juni 6, o^h UT	0,37022	-0,35812	0,08510	0,03964
3	Juni 11, 0^h UT	0,39867	-0,39444	-0,04828	0,03195

Jadi;

Dengan demikian hasil kali vektor;

Jadi i = Cos^-1 (0,99254)=7⁰ 0’ (2-39)

(2-40)

Oleh sebab itu;

(2-41)

Kita lihat bahwa N₃ identik pernyataan momentum sudut, sehingga:

(2-42)

Dari persaaman (2-41 ) dan (2-42 )

dan

(2-43)

dan

Oleh sebab itu dengan melihat persamaan ( ) kita peroleh;

0,14815 e Cos w + 0,01586 e Sin w = 0,02822 (2-44)

0,27135 e Cos w + 0,07904 e Sin w = 0,05667

Dengan menyelesaikan persamaan ini diperoleh;

e Sin w = 0,09966 dan e Cos w = 0,17981, sehingga didapat;

e = 0,2056 dan w = 28⁰ 59’,8

Selanjutnya, gunakanlah data ini untuk menentukan r, E₁ dan T . Setengah sumbupanjang a, dapat dihitung dari pernyataan (2-246).

atau a= 0,38702

karena nilai e sudah diketahui maka, untuk tanggal Juni 10,00

r₁ = a (1-e Cos E₁ ) atau 0,34200=0,38702(1-0,2056 Cos E₁ )

atau

Cos E₁ = 0,56576 atau E₁ = 55⁰ 32’,7

Dari persamaan Kepler diperoleh M₁ (ingat E₁ dinyatakan dalam radian). Jadi;

M₁ = E₁ – e Sin E₁ = 0,96944- (0,2056)(0,82457) = 0,79991 radian

Oleh sebab itu dari persamaan diatas, kita peroleh

atau 0,79991 = 0,07142(t₁ – T ) dengan demikian (t₁ – T )= 11,200

Dan karena t= Juni 1,00 = “May 32,00 “ maka kita peroleh T= 1960 May 20,800 .

Kesimpulan akhir diragakan dalam tabel berikut;

Tabel 2-4 Eleman orbit objek

No	Elemen orbit	Data
1.	Saat terakhir lewat perihelion, T	1960 May 20,800
2.	Setengah sumbu-panjang elip, a	0,38702 SA
3.	Eksentrisitas, e	0,2056
4.	Inklinasi, i	7⁰ 0’
5.	Sudut simpul naik, W	47⁰ 47’,6
6.	Argumen perihelium, w	28⁰ 59’,8

Untuk memeriksa apakah harga a yang kita peroleh sudah benar, dapat digunakan hukum Kepler III;

Karena m << 1 maka

atau a =

= 0,38712 SA

Dalam hal ini tampak sampai desimal ketiga hasil ini cukup signifikan dengan nilai setengah sumbu panjang a, yang diragakan dalam tabel diatas.

2-5 Orbit parabolik

Dari pembahasan terdahulu. Bila kita mempunyai suatu sistem orbit yang berbentuk elip, maka pada lintasan tersebut berlaku;

Konstanta Kepler
Persamaan energi sistem dengan M = m₁ + m₂ dan

m₂ << m₁

Persamaan lintasan dan

Apabila lintasan berubah menjadi parabola maka E = 0 atau dengan perkataan lain

, kemudian nyatakan v sebagai rh dan untuk saat t= T misalkan r = q maka kita peroleh

, dengan demikian untuk mencari persamaan lintasan yang berbentuk parabola dapat dilakukan dengan mengganti h pada pernyataan elip, kita peroleh;

(2-45)

Dalam hal ini f disebut anomali benar, diukur dari perihelion, untuk lebih jelas perhatikan gambar berikut ini

Gb.2-12 Ilustrasi komet yang melintasi Matahari dalam orbit parabola

Selanjutnya dari konstanta luas (setelah mengganti q dengan f) , kita peroleh;

oleh sebab itu dapat ditulis;

(2-46)

Andaikan;

Pada saat T komet ada di perihelion (f = 0)
Pada saat t, komet ada di tempat lain (f¹0)

Maka bila persamaan diatas diintegrasikan dari saat T ke t , ruas kanan harus kita integrasikan dari 0 sampai f, bila diselesaikan diperoleh;

(2-47)

Untuk menyederhanakan persamaan ini misalkanlah;

akibatnya,

Substitusikan persamaan ini pada (2-47) dan ambillah M sebagai satuan maka kita peroleh;

(2-48)

jadi

atau

(2-49)

Oleh sebab itu untuk menentukan f harus diselesaikan lebih dahulu tiga persamaan berikut secara berurutan.

Harap diingat, dalam hal massa matahari M diambil sebagai satuan, maka nilai k adalah konstanta Gauss dan q dalam satuan astronomi

Studi Kasus 1. Komet dalam orbit parabola

Bongkahan es raksasa datang dari awan Oort dan bergerak mengelilingi Matahari, ketika mendekati Jupiter orbitnya diganggu sehingga menjadi parabola. Sesaat setelah melewati Matahari bongkahan es tadi mencair dan kandungan gas beku menguap ke dalam ruang antar planet, selanjutnya benda terlihat sebagai komet

Persoalannya

(a) hitung kecepatan lingkaran benda tersebut V₀ pada jarak 2 tahun cahaya dari Matahari

(b) hitung kecepatan komet di titik A(lihat gambar), ketika ia berjarak r_A=2 SA dan ketika ia berada di perihelion r_P= 1SA

Penyelesaian

(a) untuk orbit lingkaran;

jarak komet ke matahari R = 2 tahun cahaya= 2x9,5 10¹⁵ = 19´10¹⁵ meter

= 84 m/det

(b)

Gb.2-13 Lintasan parabola sebuah komet, P titik perihelion sedangkan A titik sembarang pada orbit, p menyatakan lotus rectum, q jarak perihelion dan hubungannya adalah p=2q

Persamaan energi yang berlaku adalah;

Karena

1) r₀ = 2 ly merupakan jarak yang jauh lebih besar dibandingkan denagn r_A= 2 SA

2) V₀ << V_A

maka dapat ditulis

km/det

Untuk lintasan parabola berlaku, e=1, jadi diperoleh bentuk;

Atau secara singkat ;

Dalam hal ini q titik terdekat komet(perihelium)

Studi Kasus 2. Menentukan massa bintang ganda visual

Pendekatan two-body problem dapat digunakan untuk menentukan massa bintang ganda visual, bila magnitudo bolometrik(magnitudo untuk seluruh panjang gelombang)diketahui. Untuk itu dalam mempelajari dinamika system bintang berdua visual ada beberapa pernyataan yang dapat digunakan untuk menghitung jarak dan massa bintang

1. Paralak dinamik

Tinjau hukum harmonik;

Untuk bintang ganda visual M₁ dan M₂ hampir sama besarnya, massa bintang yang satu tidak bisa diabaikan terhadap massa yang lain, selain itu setengah sumbu panjang orbit, a dinyatakan dalam detik busur dan jarak dinyatakan dalam parsek sedangkan massa dalam satuan massa Matahari. Karena paralak p=1/d dalam detik busur, maka pernyataan diatas menjadi;

Jika P dalam tahun, massa dalam satuan massa matahari maka bentuk pernyataan diatas ini menjadi

Pernyataan diatas, dikenal sebagai paralak dinamik dalam hal ini;

p- paralak dalam detik busur dan setengah sumbu panjang, a dalam detik busur

P-periode revolusi dinyatakan dalam tahun

M_i massa bintang ke- i dalam satuan massa matahari

2. Magnitude bolometric versus paralak

M_b – magnitude absolute bolometric

m_b – magnitude semu bolometrik

p –paralak

3. Hubungan massa-luminositas

Log M = 0,1 (4,6 - M_b ) bila 0 < M_b < 7,5

Log M = 0,145 (5,2 - M_b ) bila 7,5 < M_b < 11

Dalam hal ini M – massa bintang

Sebagai contoh akan dihitung massa dan jarak bintang ganda visual ADS 1733 dengan elemen orbit sebagai berikut;

Informasi tentang bintang ganda visual ADS 1733

Elemen Orbit	Luminositas
a = 1”,673	m_b1=8,1 magnitude semu bolometrik bintang primer
e = 0.426.	m_b2=9,1 magnitude semu bolometrik bintang sekunder
P = 168.303 tahun.

Dalam nomenklatur simbol setengah sumbu panjang orbit elips, a umumnya diganti dengan a. Hasil iterasi diperlihatkan dalam tabel berikut;

Iterasi	M₁ + M₂	p	¹M_b	²M_b	M₁	M₂
0	2	0,04356	6,305457	7,305457	0,691367	0,495118
1	1,593243	0,04699	6,470037	7,470037	0,654402	0,468646
2	1,358145	0,049559	6,585597	7,585597	0,629634	0,450909
3	1,219344	0,051372	6,66363	7,66363	0,613442	0,439313
4	1,13605	0,052598	6,714845	7,714845	0,603042	0,431865
5	1,085478	0,053402	6,747805	7,747805	0,596442	0,427138
6	1,054529	0,05392	6,768743	7,768743	0,592287	0,424163
7	1,035489	0,054248	6,781931	7,781931	0,589685	0,422299
8	1,023737	0,054455	6,790193	7,790193	0,58806	0,421136
9	1,016466	0,054585	6,795352	7,795352	0,587049	0,420411
10	1,011963	0,054665	6,798566	7,798566	0,586419	0,41996

Dalam tabel diatas sebagai nilai awal diambil M₁ + M₂ = 1

Kesimpulan

Paralak dinamik bintang ganda tersebut adalah 0”,054 atau jaraknya d=18,52 parsek
Massa dari bintang ganda tersebut adalah M₁ = 0,58 M₀ dan M₂ = 0,42 M₀
Magnitudo absolut bolometrik bintang tersebut adalah ¹M_b = 6,79 dan ²M_b = 7,79

Studi Kasus 3. Menentukan periode dari luas daerah yang disapu

Diketahui sebuah planet bergerak dalam orbit elips, dengan F adalah posisi Matahari seperti gambar berikut ini, busur BPB’ ditempuh dalam waktu 2T₁. Sedangkan untuk busur B’AB, diperlukan waktu 2T₂

Pertanyaannya, buktikan bahwa

Bukti

Tinjau lintasan setengah elips BPB’

Menurut hukum Kepler : Dua kali luas daerah yang disapu persatuan waktu adalah tetap yaitu sebesar h(momentum sudut) dengan;

Luas D BFB’ =

Luas daerah PBFB’ adalah: Luas BPB’ – Luas DBFB’ =

= hT₁

Luas daerah BFB’A = Sisa luas daerah =

= hT₂

Rasio luas kedua daerah tersebut(PBFB’/BFB’A) adalah ;

Oleh sebab itu jika T₁ dan T₂ diketahui maka setengah periode orbit,T dapat dicari, yaitu T= T₁ + T₂ atau periode P=2T

Studi Kasus 4. Menentukan definisi 1 satuan astronomi pada saat asteroid mendekati Bumi

Beberapa dekade yang lalu Eros mendekati Bumi, banyak informasi yang dapat dipelajari tatkala ada benda langit yang mendekati Bumi. Pada saat oposisi dilakukan pengamatan Eros dari dua observatorium A dan B yang terpisah sejauh 1519 kilometer, masing-masing observatorium mengamati Eros dan bintang standard yang sama (lihat gambar). Sudut diantara dua objek tadi adalah SAE =6² sedangkan sudut EBS = 8². Ketika pengamatan dilakukan Eros dan Bumi sedang berada diperihelium. Andaikan Bumi dan Eros adalah co-planar hitunglah paralak harian Eros. Selain itu definisi satuan astronomi juga bisa direvisi kembali dengan datangnya Eros. Jika Eros mempunyai periode P= 642 hari dan eksentrisitasnya, e = 0,223 tentukanlah jarak Bumi ke Matahari pada saat Eros diamati dalam satuan kilometer dan bandingkan hasilnya dengan data sebelumnya

Penyelesaian

a) Jarak Eros;

Untuk mengukur jarak Eros ditentukan sudut SAE dan sudut SBE dengan satu bintang standar,S bintang terlihat sejajar baik dari titik A maupun titik B

Pendekatan yang dilakukan

jarak AB bisa diabaikan terhadap jarak Eros-Bumi
bintang yang sama terlihat sejajar dari A dan B

dengan demikian;

Tg AEB' =Tg AEB=AB'/ B'E=AB/BE dalam hal ini BE adalah jarak eros ke Bumi, d dengan demikian, sudut AEB= 6"+8"=14"

kilometer=22,38 ´10⁶ km

Besaran ini merupakan jarak minimum Eros ke Bumi, sehingga paralaknya menjadi maksimum yaitu;

(b) dari hukum Kepler:

Pada saat Eros dan Bumi di perihelion berlaku;

SE=a(1-e)=1,457´(1-0,223)=1,47´0,777=1,132 SA

SB=a(1-e)=1´ (1-0,0167)=0,983=0,983 SA

Jarak Eros ke Bumi = SE-SB®22,38 ´10⁶ km= 0,149 SA

Sehingga ; definisi 1 SA = 22,38x10⁶ / 0,149 = 150,2 ´10⁶ km » 150 ´10⁶ kilometer

Pesan dari soal ini adalah, ternyata dengan mengamati asteroid orang bisa merevisi kembali definisi satu satuan astronomi.

Gambar 7 Geometri posisi Bumi dan Eros pada saat pengamatan

dalam hal ini S-matahari, B-Bumi dan E- Eros

Studi Kasus 5. Menentukan paralak trigonometri dari dua tempat di Bumi

Pada tanggal 13 April 2029, sebuah asteroid 99942-Apophis mendekati Bumi, pada saat itu jaraknya adalah 0,10 LD(lunar distance = jarak rerata Bumi-Bulan). Sekelompok astronom akan mengukur paralak asteroid tersebut dari Observatoire de Paris dan Naval Observatory Washington, secara simultan. Posisi geografi kedua observatorium tersebut adalah;

Observatoire de Paris(France):

l₁= 2^o20’14² Timur = - 2^o20’14² dan j₁= 48^o50’11² Utara = + 48^o50’11²

Naval Observatory Washington(USA);

l₂= 77^o03’56² Barat = +77^o03’56² dan j₂= 38^o55’17² Utara = +38^o55’17²

Pertanyaannya:Hitunglah jarak linier kedua observatorium tersebut dan berapakah paralak asteroid tersebut bila dihitung ?

Penyelesaian

Jika koordinat geografi, longitude(bujur), l dan latitude(lintang), j dua titik di permukaan Bumi maka jarak sudut keduanya d, dapat dihitung dari;

Jarak liniernya dapat dihitung dari;

R, jari-jari Bumi yaitu 6371 kilometer dan d dalam derajad maka S dalam kilometer

Jika d dalam radian maka

dalam hal ini S dalam kilometer

Dengan memasukkan data diatas diperoleh jarak Washington-Paris adalah, S= 6181,6 km

Jarak asteroid d = 0,1 LD = 0,1´384400 km = 38440 km

Efek projeksi kedudukan asteroid relatif terhadap bintang latar belakang

Paralak asteroid a, dapat dihitung dari;

= 9^o,21

jadi paralaknya adalah a = 9^o,21

Ragam Soal Latihan

1) Planet Mars mempunyai elongasi j₁= 60⁰ dan pada saat bersamaan sebuah asteroid tampak dengan sudut phase b = 30⁰ dan elongasi j₂ = 45⁰. Jika jarak Mars-Matahari 1,5 Satuan Astronomi dan kamu sekarang berumur 17 tahun, namun sejak lahir kamu tinggal di Mars berapakah umurmu sekarang dalam penanggalan Mars ?. Selanjutnya hitunglah.

a) jarak asteroid itu dari Matahari (r₁ )

b) jarak Mars dari Bumi (r₂)

c) jarak asteroid dari Mars (r₃ )

d) tempo yang diperlukan asteroid untuk kembali ke posisi semula relatif terhadap bintang latar belakang

2) Orbit Parabola

a. Energi Total, E_t, sistem dua benda yang bergerak mengitari pusat massanya dapat dinyatakan dalam pernyataan:

E_k + E_p = E_t

Dalam hal ini E_k dan E_p masing-masing menyatakan energi kinetis dan energi potensial. Energi total dapat berharga negatif, nol dan positif. Deskripsikan kriteria energi sistem untuk lintasan elips, parabola dan hiperbola! Uraikan jawab saudara!

b. Komet periode panjang dianggap mempunyai lintasan parabola. Anomali benar f komet tersebut dapat dicari dari persamaan Baker:

sedangkan jaraknya ke Matahari dihitung dari pernyataan:

Ingat, jika massa Matahari, M, diambil sebagai satuan, jarak dalam satuan astronomi (SA), waktu dalam tahun (year), dan konstanta Gauss, k, maka konstanta ini nilainya adalah k = 0,01720209895

Pertanyaannya: Seandainya komet X berada di perihelium pada tanggal 6 Januari 2009 pada pukul 24:00 dengan jarak q =1,2 SA, berapakah jaraknya ke Matahari dan anomali benar komet tersebut pada tanggal 8 Januari 2009 pada jam yang sama? Untuk menjawab pertanyaan ini, gunakan metoda numerik! Lengkapi prosedur perhitungan dengan diagram alir (flowchart)!

3) Sebuah asteroid bergerak mengelilingi Matahari dengan periode 4,5 tahun dan mempunyai setengah sumbu pendek lintasannya yang berbentuk elip 3,2 SA. Bila eksentrisitas asteroid itu e = 0,4 . Berapakah luas daerah yang telah disapu oleh asteroid itu selama 1,5 tahun. Seandainya asteroid itu beroposisi pada tanggal 31 Juli 2006, tanggal berapakah ia akan beroposisi kembali?

4) Sebuah asteroid bergerak dengan orbit elips, jika eksentrisitasnya adalah e buktikan bahwa rasio kecepatan kuadrat di aphelion terhadap kecepatan kuadrat di pherihelion adalah;

5) Sekelompok peneliti cuaca meluncurkan roket dari titik A yang terletak di ekuator menuju pulau kecil B(30⁰LU) yang ada diatasnya. Jika tempo yang dibutuhkan untuk tiba di pulau B tersebut adalah 10 menit. Pertanyaannya, apakah roket itu akan jatuh di pulau B itu ?, jika tidak dimanakah ia jatuh ?. Ambil untuk Bumi, jejari, R = 6378 km dan periode rotasi = 24 jam

6) Berikut diberikan data bintang ganda visual;

Tabel informasi tentang 3 bintang ganda visual

Star	Visual		P [tahun]	a [²]	Spectrum		p [²]
Star	m₁	m₂	P [tahun]	a [²]	S₁	S₂	p [²]
a Centauri	-0.04	1.38	79.9	17.58	G2	K5	0. ²751
h Cas	3.47	7.22	480.0	11.99	G0	M0	0. ²176
e Hyd	3.7	4.8	15.0	0.21	G0	?	0. ²014

Ingat dalam fotometri, jika kita ingin menggunakan hubungan massa-luminositas magnitude visual harus dinyatakan dulu dalam magnitude bolometrik. Pertanyaannya;

a. Jika kita menggunakan iterasi perlukah dilakukan koreksi terhadap magnitude semu ?.

b. Hitunglah massa masing-masing bintang dengan 2 cara;

1. sebagai tebakan awal ambil M₁+M₂= 1

2. sebagai nilai awal paralak dinamik ambil p = 0.1

c. Hitunglah galat relatif paralak dinamik tiap bintang jika sebagi acuan diambil data paralak yang dianut orang selama ini,

a Centauri = 0.”751

h Cas = 0.”176

e Hyd = 0.”014

astronomi

Jumat, 15 Juli 2011

orbit dalam ruang

2) Orbit Parabola