TELAAH ORBIT SATELIT LAPAN-TUBSAT Nizam Ahmad dan Thomas Djamaluddin Pusat Pemanfaaian Sains Aniariksa LAPAN Email:
[email protected] ABSTRACT The study of orbit characteristics of LAPAN-TUBSAT can be done by u s i n g some microsatellites which identical in t h e case of mission a n d its orbit. This knowledge is useful to reduce t h e failure of mission. This characteristic can be s e e n from the change of orbital elements c a u s e d by perturbation forces. LAPAN-TUBSAT is placed in Low Earth Orbit (LEO) at altitude about 630 km which h a s orbit n e a r polar ( i * 98°). At this altitude, t h e perturbations come from non-gravitational force s u c h as atmospheric drag a n d gravitation force from e a r t h oblatness. From t h e simulation can be predicted t h a t in the early y e a r s of its operation, the variation in altitude a n d semi major axis is relatively small. It m e a n s t h a t this satellite could have t h e life time more t h a n 50 years. The effect of earth oblatness will cause t h e regression of the n o d e s of about l ° / d a y a n d the rotation of t h e line of apsides of about -3°/day. These changes a r e not too critical for n e a r polar orbit which m e a n s t h a t satellite keeps conducted t h e mission goal s u c h as satellite imaging. ABSTRAK Kajian karakteristik orbit satelit LAPAN-TUBSAT dapat dilakukan melalui beberapa satelit mikro yang identik dalam hal mist dan orbitnya. Pengetahuan a k a n karakteristik ini sangat perlu u n t u k memperkecil kegagalan misi satelit. Karakteristik ini terlihat dari s e b e r a p a b e s a r p e r u b a h a n elemen orbit akibat gangguan t e r h a d a p misi satelit. Satelit LAPAN-TUBSAT yang ditempatkan p a d a ketinggian orbit r e n d a h bumi (LEO), dengan ketinggian sekitar 630 km, memiliki orbit hampir polar (s* 98°). Pada ketinggian ini, gaya gangguan non-gravitasional p a d a orbit berasal dari h a m b a t a n atmosfer d a n gaya gravitasi akibat efek Bumi pepat. Dari simulasi, dapat diperkirakan b a h w a p a d a t a h u n p e r t a m a satelit ini beroperasi, variasi ketinggian d a n s u m b u semi mayor relatif kecil. Satelit diperkirakan dapat bertahan dengan kala hidup orbit di a t a s 50 t a h u n . Efek Bumi pepat menyebabkan satelit mengalami presesi nodal sekitar i V h a r i d a n rotasi apsida sekitar -3°/hari. P e r u b a h a n ini tidak terlalu beresiko t e r h a d a p satelit dengan orbit hampir polar, malah akan memberikan k e u n t u n g a n bagi misi pencitraan satelit. Kata k u n c i : LAPAN-TUBSAT, Orbit, Lingkungan antariksa
1
1
PENDAHULUAN
Satelit mikro p e r t a m a Indonesia dibuat a t a s d a s a r kerjasama a n t a r a LAPAN dan Technical University of Berlin (TV Berlin), k a r e n a n y a d i n a m a k a n LAPAN-TUBSAT. Satelit dengan berat sekitar 57 kg dan u k u r a n kurang dari 1 m ini bermisi pengawasan (Surveillance). Satelit ini memiliki m u a t a n d u a kamera penginderaan j a u h dengan f B 1000 mm d a n f • 50 mm, k o m u n i k a s i d a t a simpan d a n a n t a r [store and forward), serta kontrol sikap tiga s u m b u [3-axis attitude control}. Satelit ini telah diluncurkan secara piggy back (menumpang) dengan roket PSLV (Polar Satellite Launch Vehicle) milik India pada tanggal 10 J a n u a r i 2007. Dalam p e r e n c a n a a n awal a d a tiga alternatif orbit yang digunakan oleh satelit ini (Hardienata, 2005) yaitu orbit sinkron-matahari (Sun-Synchronous Orbit, SSO) dengan inklinasi, i • 97,9° d a n ketinggian H = 630 ]ga, orbit ekuatorial dengan inklinasi : = 2,5° d a n ketinggian H = 550 km, serta orbit semi-ekuatorial dengan inklinasi i - 37° dan ketinggian H = 6 0 0 km. Pemilihan orbit ini sangat vital bagi optimasi aplikasi dan operasional satelit. Hasil simulasi awal memperlihatkan b a h w a p e n e m p a t a n satelit di ketiga orbit tersebut memiliki karakteristik tersendiri y a i t u orbit SSO dengan inklinasi 97,9° d a n 37° a k a n memberikan hasil yang lebih baik u n t u k aplikasi satelit k a r e n a m e m u n g k i n k a n n y a meliput s e l u r u h wilayah Indonesia d a n orbit ekuator dengan inklinasi 2,5° a k a n memberikan hasil yang lebih baik bagi operasional satelit k a r e n a waktu lintas ulang (revisit) yang lebih sering (Hardienata, 2005). Untuk p e m a n t a u a n d a n pengambilan data, LAPAN mengoperasikan d u a s t a s i u n Bumi (SB) di Rumpin (Bogor) dan Biak. Lokasi ini dipilih b e r d a s a r k a n k e m a m p u a n m e l a k u k a n liputan t e r h a d a p wilayah k e p u l a u a n Indonesia. SB u t a m a berlokasi di Rumpin, Bogor yang memiliki koordinat dengan lintang 06°22'27" S dan bujur 106°37'87" T. SB ini a k a n m e m a n t a u satelit di bagian b a r a t wilayah Indonesia. SB yang k e d u a yakni Biak memiliki koordinat dengan lintang 01°10'36" S d a n bujur 136°06'04" T. SB ini a k a n m e m a n t a u satelit di bagian timur wilayah Indonesia. Satelit LAPAN-TUBSAT ditempatkan p a d a ketinggian orbit r e n d a h Bumi sekitar 6 3 0 km. Mengingat kondisi lingkungan antariksa diketinggian orbit satelit tidak lepas dari p e n g a r u h gaya gangguan, m a k a perlu a d a n y a kajian t e r h a d a p karakteristik orbit akibat p e n g a r u h gaya gangguan ini (Hasting, 1996). Hal ini selain mengurangi resiko kerugian anggaran j u g a agar k e b u t u h a n misi satelit ini nantinya dapat terpenuhi. G a m b a r a n u m u m karakteristik orbit LAPAN-TUBSAT ini diperoleh melalui kajian t e r h a d a p beberapa satelit mikro yang identik dalam hal misi m a u p u n orbitnya (Ahmad, 2006). P e r m a s a l a h a n y a n g k e r a p terjadi p a d a g a n g g u a n orbit adalah b e r u b a h n y a ketinggian satelit akibat gaya h a m b a t a n atmosfer (Atmospheric drag) d a n terjadinya presesi nodal serta rotasi apsida akibat efek b u m i p e p a t 2
{Earth oblatness)., Makalah ini a k a n m e m b a h a s hasil analisis karakteristik orbit satelit mikro sejenis LAPAN-TUBSAT d a n prakiraan gangguan d a n kala hidup LAPAN-TUBSAT. 2
DATA DAN PENGOLAHAN
Secara s t r u k t u r , LAPAN-TUBSAT d a p a t dibagi menjadi d u a bagian (rak), rak bawah d a n rak a t a s (Hardienata, 2005). Rak bawah berisi sistem kontrol sikap, sistem telemetry and telecommand (TTC), m u a t a n k a m e r a dengan fokus lensa 1000 mm d a n sistem S-band.
STAR SENSOR
G a m b a r 2 - 1 : Rak bawah LAPAN-TUBSAT Rak a t a s terdiri dari baterai, kontrol tenaga, sistem data handling (PCDH), m u a t a n k a m e r a dengan fokus lensa 50 mm d a n air coil.
i
G a m b a r 2-2 : Rak a t a s LAPAN-TUBSAT 3
Secara keseluruhan, satelit ini memiliki u k u r a n sekitar 45 x 45 x 27 cm, berbentuk kotak dengan m a s s a sekitar 57 kg. Lebih detailnya, d a t a - d a t a satelit ini dapat dilihat p a d a Tabel 2 - 1 .
Nation
T a b e l 2 - 1 : DATA SATELIT LAPAN-TUBSAT Indonesia, Germany
Type/App
Earth Observing
Operator
LAPAN/TU Berlin
Contractors
TU-Berlin
Equipment Configuration Propulsion
1 x 3CCD Color-Camera with 6 m GSD, 1 CCD Color-Camera 45 x 45 x 27 cm b u s None
Lifetime
1 year (design)
Mass
57 kg
Orbit
635 km x 635 km, 98°
Telaah awal orbit satelit LAPAN-TUBSAT ini dilakukan melalui satelit mikro lain yang telah diluncurkan dan memiliki k e s a m a a n dalam h a l misi m a u p u n orbitnya. Melalui kajian karakteristik orbit satelit mikro (Nizam, 2006a), a d a b e b e r a p a satelit yang dapat digunakan u n t u k memperkirakan orbit LAPAN-TUBSAT ini, diantaranya adalah Alsat 1 (Aljazira), S u n s a t (Afrika Selatan), Bilsat 1 (Turki), Fedsat (Australia) d a n Rocsat 1 (China). B e r d a s a r k a n inklinasi orbit, satelit Alsat 1, Sunsat, Bilsat 1 d a n Fedsat m e r u p a k a n satelit Near Polar d a n satelit Rocsat 1 adalah satelit Ekuatorial. D a t a - d a t a orbit kelima satelit ini dapat dilihat p a d a tabel 2-2 (Spacetrack, com) Tabel 2-2 : DATA - DATA ELEMEN ORBIT MIKROSAT Name Alsat-1 Sunsat Bilsat-1 Fedsat Rocsat 1
Period 98.44 99.77 98.44 100.85 96.1
Incl, (deg) 98.1 96.48 98.1 98.53 34.98
Apogee {km) 691 854 692 806 594
Perigee (km} 677 641 676 792 549
:
Metodologi yang digunakan dalam telaah awal orbit satelit LAPANTUBSAT ini adalah, p e r t a m a memanfatkan d a t a tiga alternatif orbit tersebut. J a d i b e r d a s a r k a n inklinasinya, LAPAN-TUBSAT telah memiliki alternatif orbit near polar d a n equatorial. Berikutnya adalah mencari satelit mikro yang memiliki misi d a n orbit yang identik dengan satelit LAPAN-TUBSAT d a n ini 4
dilakukan t e r h a d a p kelima mikrosat p a d a Tabel 2-2 di a t a s . Terakhir meaganalisis b e b e r a p a elemen orbit yang secara langsung dipengaruhi oleh gangguan (evolusi ketinggian, presesi nodal d a n rotasi apsida) u n t u k memperoleh karakteristik orbitnya. Hal ini s e c a r a tidak langsung memberikan kesamaan g a m b a r a n p a d a karakteristik LAPAN-TUBSAT nantinya. Untuk simulasi lintasaii digunakan perangkat lunak STSplus dengan menggunakan d a t a - d a t a orbit (TLE) mikrosat pada Tabel 2-2 di a t a s dengan sedikit p e n y e s u a i a n epoch d a n propagasi orbit melalui SGP4. 3
HASIL
Prediksi evolusi ketinggian satelit LAPAN-TUBSAT dilakukan dengan menggunakan model propagasi SGP4. Karena satelit ini b e b e r a p a kali mengalami p e n u n d a a n peluncuran, m a k a r a n c a n g a n orbitnya j u g a mengalami p e r u b a h a n . Namun, berdasarkan alternatif orbit tersebut, diambil jadwal p e l u n c u r a n [Launch Opportunity) p a d a tanggal 10 J a n u a r i 2007. Bila mengacu p a d a d a t a TLE awal orbit mikrosat, m a k a dengan m e n g g u n a k a n data TLE awal mikrosat ini dilakukan analogi sebagai berikut, A. Diluncurkan dengan orbit near polar (= 10 Januari 2007) 1. LAPAN-TUBSAT mengacu pada satelit Alsat 1 1 27559U 02054A 0 7 0 1 0 . 0 0 0 0 0 0 0 0 - . 0 0 0 0 0 0 4 6 00000-0 00000+0 0 2 27559 098.0000 217.0176 0043370 072.5437 288.0513 14.54233099
10 03
2. LAPAN-TUBSAT m e n g a c u p a d a s a t e l i t Sunsat 1 25636U 99008C 07010.00000000-.00000032 00000-0 00000+0 0 2 25636 098.0000 009.5999 0150608 252.6338 105.8740 14.40957591
11 24
3. LAPAN-TUBSAT m e n g a c u pada satelit Bilsat 1 1 27943U 03042E 0 7 0 1 0 . 0 0 0 0 0 0 0 0 - . 0 0 0 0 0 0 4 8 00000-0 00000+0 0 2 27943 098.0000 157.4488 0011745 216.7135 143.3214 14.62490622
14 111
4. LAPAN-TUBSAT m e n g a c u pada satelit Fedsat 1 27598U 02056B 07010.00000000-.00000044 00000-0 00000+0 0 2 27598 098.0000 057.7273 0009377 032.3043 327.8701 14.27685195
11 57
B. Diluncurkan dengan orbit equatorial (* 10 Januari 2007) LAPAN-TUBSAT m e n g a c u pada satelit Rocsat 1 1 25616U 99002A 07010.00000000-.00000027 00000-0 00000+0 0 2 25616 02.50000 284.7377 0009437 177.6235 182.4516 14.90907854
11 00
Dengan m e n g g u n a k a n d a t a - d a t a TLB di a t a s sebagai input, m a k a diperoleh perkiraan evolusi ketinggiannya selama 5 t a h u n sejak mengorbit yang dapat dilihat p a d a Gambar 3 - 1 . 5
Alsat 1
Sunsat
Fedsat
Bilsat 1
Rocsat 1
G a m b a r 3 - 1 : Prediksi evolusi h LAPAN-TUBSAT G a m b a r 3-1 m e r u p a k a n profil ketinggian orbit satelit LAPAN-TUBSAT hasil propagasi dengan m e n g g u n a k a n SGP4, yaitu propagasi orbit yang memperhitungkan pengaruh gaya gangguan yang berasal dari h a m b a t a n atmosfer d a n efek Bumi pepat dalam kalkulasi p e r u b a h a n elemen orbitnya. Prakiraan evolusi ketinggian orbit LAPAN-TUBSAT ini dilakukan dengan mengadopsi d a t a TLE awal satelit mikro p a d a Tabel 2-2 yang disesuaikan p a d a epoch p e l u n c u r a n satelit, yaitu tanggal 10 J a n u a r i 2 0 0 7 . Dari sini terlihat bahwa evolusi ketinggian orbit satelit selama 5 t a h u n memperlihatkan ketinggian satelit masih di a t a s 500 km.
Gambar 3-2 : Prediksi evolusi a LAPAN-TUBSAT 6
Evolusi pada elemen s u m b u semi mayor dapat dilihat pada Gambar 3-2. Dari profil s u m b u semi-mayor (a) p a d a Gambar 3-2 tersebut d a p a t dilihat bahwa r a t a - r a t a a selama evolusi orbit 5 taKun m a s i h berkisar sekitar 7000an km. B e r d a s a r k a n kajian karakteristik orbit beberapa mikrosat, LAPANTUBSAT diperkirakan a k a n mengalami p e r u b a h a n p a d a empat elemen orbit, yakni ketinggian (h) a t a u s u m b u semi mayor (a), Asensio rekta titik tanjak (/?) dan Argumen perigee (co) (Nizam, 2006a). LAPAN-TUBSAT diperkirakan a k a n memiliki pola evolusi ketinggian yang s a m a dengan mikrosat p a d a Tabel 2-2. Bila orbitnya ekuatorial, p e n u r u n a n s u m b u semi-mayor a k a n lebih cepat dibandingkan orbit near polar. P a d a ketinggian di a t a s 600 km, kala h i d u p {Life Time) orbit diperkirakan bisa b e r t a h a n di a t a s 10 t a h u n . (Wertz, 2001). Selain itu, satelit ini a k a n mengalami regresi nodal r a t a - r a t a sebesar 1 p u t a r a n / t a h u n , bila orbitnya near polar d a n bila orbitnya equatorial, presesi rata-ratanya berkisar antara 5-6 p u t a r a n / t a h u n . D a m p a k berikutnya t e r h a d a p orbit satelit ini adalah rotasi apsida yang m e n y e b a b k a n elemen co berlawanan a r a h dengan gerak satelit u n t u k orbit near polar. Akibatnya satelit a k a n m e n c a p a i perigee lebih l a m b a t u n t u k setiap revolusinya. Untuk orbit equatorial, elemen co a k a n bergerak s e a r a h dengan gerak satelit. Akibatnya satelit a k a n mencapai perigee lebih cepat untuk setiap revolusinya (Vallado, 2001). Simulasi lintasan satelit LAPAN-TUBSAT d e n g a n m e n g g u n a k a n d a t a TLE di a t a s d a p a t d i l a k u k a n dengan m e n g g u n a k a n p e r a n g k a t l u n a k STSplus yang didalamnya telah m e n g g u n a k a n model propagasi SGP4. Simulasinya dapat dilihat p a d a Gambar 3-1 hingga Gambar 3-7. 1. LAPAN-TUBSAT menga cu p a d a satelit Alsat 1
Gambar 3-3 : Prediksi lintasan LAPAN-TUBSAT 1
7
Satelit a k a n b e r a d a p a d a lintang 72,22° U dan bujur 109,16° B melintasi k a w a s a n sekitar S a n Francisco d a n New York dengan ketinggian orbit sekitar 753 km. 2. LAPAN-TUBSAT m e n g a c u p a d a satelit S u n s a t
G a m b a r 3-4 : Prediksi lintasan LAPAN-TUBSAT 2 Satelit a k a n b e r a d a p a d a lintang 39,24° U d a n bujur 170,18° B serta berada di sekitar daerah San Francisco d a n Tokyo p a d a ketinggian orbit sekitar 689 km. 3. LAPAN-TUBSAT mengacu p a d a satelit Bilsat 1
Gambar 3-5 : Prediksi lintasan LAPAN-TUBSAT 3 Satelit a k a n b e r a d a p a d a lintang 79,32° S d a n bujur 33,64° T setelah wilayah timur Indonesia menuju daerah South Atlantic Anomaly (SAA) dengan ketinggian orbit sekitar 718 km. 8
4. LAPAN-TUBSAT mengacu p a d a satelit Fedsat
G a m b a r 3-6 : Prediksi lintasan LAPAN-TUBSAT 4 Satelit a k a n berada p a d a lintang 19,19° S d a n bujur 114,03° B sekitar wilayah Rio De Janeiro menuju kota San Francisco dengan ketinggian orbit sekitar 809 km. 5. LAPAN-TUBSAT mengacu p a d a satelit Rocsat 1
G a m b a r 3-7 : Prediksi lintasan LAPAN-TUBSAT 5 Satelit a k a n b e r a d a p a d a lintang 2,01° S d a n bujur 28,76° B di sekitar daerah ekuator Bumi dengan ketinggian orbit sekitar 589 km. Lintasan p a d a Gambar 3-7 di a t a s memperlihatkan satelit ini a k a n melewati wilayah Indonesia di sekitar ekuator Bumi. 9
4
PEMBAHASAN
D a t a - d a t a TLE yang digunakan di atas m e r u p a k a n d a t a - d a t a awal satelit mikro ketika mengorbit awal {insertion). Analogi yang dilakukan p a d a LAPAN-TUBSAT adalah dengan modifikasi epoch yang sesuai d e n g a n Launch Oppurtinity-nya. yaitu p a d a tanggal 10 J a n u a r i 2007 serta p e n y e s u a i a n inklinasi b e r d a s a r k a n alternatif orbit awal LAPAN-TUBSAT. Merujuk p a d a hasil Djamaluddin (2005), dengan menganalisis orbit satelit mikro diperkirakan evolusi ketinggian j a n g k a panjang satelit LAPAN-TUBSAT a k a n mengalami p e n u r u n a n signifikan secara berkala setiap p u n c a k aktivitas m a t a h a r i (Gambar 4-1).
G a m b a r 4 - 1 : Perkiraan evolusi ketinggian LAPAN TUBSAT B e r d a s a r k a n analisis kondisi fisik antariksanya, satelit diperkirakan a k a n mengalami s u h u ekstrem atmosfer a t a s sekitar 700 - 2000 K bila b e r t a h a n sampai aktivitas matahari maksimum. Pada t a h u n pertama t a h u n 2007 s a a t m a t a h a r i minimum, rentang temperatur berkisar 600 - 1500 K. Gangguan rutin akibat peningkatan flux proton d a n elektron yang mungkin dialami adalah saat mengorbit di a t a s SAA. Pada t a h u n - t a h u n p e r t a m a kemungkinan gangguan akibat aktivitas matahari relatif sedikit dan gaya h a m b a t a n atmosfer tidak a k a n bernilai ekstrem selama periode ini sehingga satelit a m a n dari a n c a m a n p e l u r u h a n ketinggian (decay) a t a u p u n reentry. Hal ini j u g a terlihat melalui profil ketinggian d a n s u m b u semi-mayor p a d a Gambar 3-1 d a n 3-2 di a t a s . Dalam periode 5 t a h u n , variasi h d a n a tidak berbeda j a u h u n t u k setiap epoch. Dari sini dapat diperkirakan seberapa besar koefisien balistik yang dimiliki oleh LAPAN-TUBSAT yang nantinya berpengaruh t e r h a d a p kala hidup orbit satelit.
10
ZAMMi mam* $**!*)
Gambar 4-2 : Prediksi LT orbit satelit (Wertz, 2001) Dari Gambar 4-2, p a d a tingkat aktivitas m a t a h a r i m a k s i m u m a k a n menyebabkan LT satelit berkurang d a n sebaliknya, tergantung dari nilai koefisien balistik satelit yang m e n y a t a k a n seberapa besar p e n g a r u h gaya h a m b a t a n atmosfer t e r h a d a p satelit. Nilai koefisien ini j u g a a k a n berbedabeda u n t u k setiap ketinggian satelit. Semakin r e n d a h ketinggian, m a k a nilai koefisien ini a k a n semakin tinggi yang berarti satelit a k a n semakin cepat mengalami p e l u r u h a n ketinggian (decay). Kala hidup orbit satelit LAPAN-TUBSAT dapat dihitung dengan mengetahui d a t a - d a t a fisik satelit yang b e r u k u r a n 45 x 45 x 27 cm d a n berbentuk kotak dengan m a s s a (m) sebesar 57 kg serta l u a s p e n a m p a n g satelit (A) sebesar 0,1225 m?. Untuk b e n t u k kotak C D = 2 - 4 (Wertz & Larson, 1990, p. 207). Sehingga B = C D A/m = 0.004339 - 0.008679 m V k g atau koefisien balistik 1/B = 115.2263 - 230.4527 kg/m2 diperkirakan satelit LAPAN-TUBSAT a k a n b e r t a h a n lebih dari 50 t a h u n . Tetapi, kala hidup instrumennya diperkirakan h a n y a b e r t a h a n sekitar 1 t a h u n , s e s u a i prakiraan perancangannya. 5
KESIMPULAN
Kajian orbit t e r h a d a p beberapa satelit mikro dapat digunakan u n t u k memperkirakan karakteristik orbit satelit LAPAN-TUBSAT. Melalui d a t a - d a t a TLE s a a t satelit mengorbit, dapat dilakukan propagasi u n t u k b e b e r a p a waktu kedepan. Dengan mengacu p a d a orbit beberapa mikrosat di a t a s , dapat diperkirakan b a h w a p a d a t a h u n - t a h u n awal, orbit satelit a k a n berevolusi dengan tingkat gangguan relatif kecil. Variasi elemen orbit t i d a k j a u h berbeda selama 5 t a h u n sejak peluncurannya. Gangguan u t a m a masih berasal dari h a m b a t a n atmosfer d a n efek Bumi pepat. Dengan mengambil analogi kondisi lingkungan a n t a r i k s a 11
LAPAN-TUBSAT yang identik dengan kondisi lingkungan antariksa mikrosat di a t a s , dapat diperkirakan evolusi ketinggian orbitnya setelah 5 t a h u n m a s i h berada di a t a s 500 km. Tetapi dalam j a n g k a panjang ketinggiannya a k a n mengalami p e n u r u n a n signifikan setiap m a s a aktif m a t a h a r i N a m u n diperkirakan kala h i d u p orbitnya bisa m e n c a p a i lebih dari 50 t a h u n , walau p u n kala hidup i n s t r u m e n n y a h a n y a sekitar 1 t a h u n . UCAPAN TERIMAKASIH Penulis m e n g u c a p k a n terima kasih k e p a d a Abdul R a c h m a n yang telah m e m b a n t u dalam pemograman serta k e p a d a Neflia yang banyak m e m b a n t u dalam pencarian d a t a - d a t a y a n g diperlukan. Terakhir u c a p a n terima kasih j u g a ditujukan k e p a d a Dr. Suwarto Hardienata dan Dr. T.S. Kelso a t a s konsultasinya melalui e-mail DATTAR RUJUKAN Ahmad, N., 2007. Karakterisuk Orbit Satetit Mikro di Ketinggian LEO, J u r n a l S a i n s Dirgantara, Vol. 5 No. 1, Desember 2007 (nomor ini). Djamaluddin, T., 2005. Pengaruh Aktivitas Matahari Pada Kala Hidup Satelit, J u r n a l S a i n s Dirgantara, Vol. 3, No 1. Desember 2005, him. 65. Hardhienata, S., 2005. Modelling and Simulation of LAPAN-TUBSAT Micro Satellite Orbit, International Conference on Instrumentation, Communication a n d Information Technology (ICICI) 2005 Proc., August 3 r d -5 t h , 2005, Bandung, Indonesia. Hasting D.; Garret H. 1996. Spacecraft Environment Interaction, Cambridge University Press. Ransom, D. 2006. STS Orbit Plus Satellite Tracking and Orbit Visualization Program, h t t p : / / w w w . d r a n s o m . c o m / Oktober 2006. Space-track 2006. Catalog Number, http://www.space_track.org/_Oktober 2006. Vallado, D. A., 2 0 0 1 . Fundamental of Astrodynamics ans Applications, Kluwer Academic Publishers. Wertz, J. R., 2 0 0 1 . Afission Geometry : Orbit and Constellation Design and Management, Kluwer Academic Publishers. Wertz, J. R.; a n d Larson, W. J. 1999. Space Mission Analysis And Design, Kluwer Academic Publishers.
'
12