BAB 2 TINJAUAN PUSTAKA

Bab ini membahas

penyelesaian masalah yang akan disimulasi, dan algoritma

program yang akan digunakan.

Bab IV

Hasil dan Pembahasan Bab ini memberikan hasil uji coba simulasi dengan membandingkan pola lintasan yang didasarkan pada sudut variasi sudut lintang. Titik yang diambil yaitu pada kutub bumi (90o ), khatulistiwa (0o ), dan sembarang tempat antara kutub bumi dan khatulistiwa; Medan (3o ) dan Jakarta (6o ).

Bab V

Kesimpulan dan saran Bab ini memberikan kesimpulan dari hasil perancangan program yang telah dilakukan dan juga memberikan saran-saran untuk penelitian selanjutnya.

BAB 2

TINJAUAN PUSTAKA

2.1. Rotasi Bumi

Rotasi adalah

perputaran

benda

pada

suatu sumbu yang

tetap,

misalnya

perputaran gasing dan perputaran planet Bumi pada sumbunya. Pada rotasi bumi, gerak rotasi ini terjadi pada garis/poros/sumbu utara-selatan (garis tegak dan sedikit miring ke kanan). Adapun kecepatan rotasi bumi dapat dihitung sesuai dengan persamaan 2.1. v = 2?rBumi cos ?/24

(2.1)

Dengan ? ?adalah besarnya derajat lintang suatu tempat di permukaan bumi. Faktor cos ? muncul pada persamaan 2.1 karena tiap titik di permukaan bumi berputar terhadap sumbu

Universitas Sumatera Utara

putar bumi, bukan terhadap titik pusat bumi. Seperti yang telah diketahui bahwa rBumi pada khatulistiwa sama dengan 6375 km, sehingga kecepatan bumi berdasarkan persamaan 2.1 yaitu sebesar 1669 cos ? km/jam. Jika diambil posisi pada garis khatulistiwa, atau pada ? = 0 maka akan diperoleh kecepatan bumi sebesar 463 meter per detik yang ternyata sudah melampaui besar kecepatan suara di udara. Pengaruh dari rotasi bumi ini secara garis besar, yaitu: a. Pergantian Siang dan malam b. Perbedaan waktu c. Perbedaan percepatan gravitasi bumi d. Pembelokan arah angin e. Pembelokan arus laut f.

Peredaran semu harian benda-benda langit

Salah satu cara untuk menganalisa rotasi bumi beserta pengaruhnya adalah dengan menggunakan pendulum Foucault (akan dijelaskan lebih lanjut pada subbab 2.3). (Pikatan, 2009)

r Ω

Gambar 2.1. Rotasi Bumi

2.2. Gaya Coriolis

Universitas Sumatera Utara

Gaya Coriolis adalah suatu proses alam yang dinamakan sesuai dengan penemunya yaitu Gaspard Gustave Coriolis (1844) dari Prancis. Ia menemukan bahwa rotasi bumi selalu menyebabkan simpangan terhadap setiap gerakan yang terjadi pada permukaan bumi. Besar kecilnya penyimpangan tergantung dari lintang geografi, dimana gerakan itu terjadi. Pada belahan bumi utara setiap gerak di belokkan ke kanan dan pada belahan bumi selatan di belokkan ke kiri dari arah kecepatan. Pembelokkan atau perubahan arah ini hanya dapat di terangkan dengan menganggap atau menambahkan adanya gaya tidak nyata yang bekerja pada suatu benda. Gaya khayal atau gaya fiktif ini dinamakan gaya coriolis. Gaya semacam ini yang bekerja pada satu satuan massa atau percepatan yang ditimbulkan disebut percepatan coriolis.

Besarnya percepatan coriolis secara matematis dapat dinyatakan sebagai berikut.

a c = v 2Ω sinφ

(2.2)

v adalah kecepatan benda yang bergerak, ? kecepatan sudut rotasi bumi yang besarnya sama dengan 7, 27x10− 5 rad / s , ? besarnya derajat lintang. Dari persamaan tersebut di atas dapat dilihat bahwa untuk kecepatan yang sama pengaruh coriolis terbesar terdapat di kutub. Pengaruh ini makin berkurang dengan berkurangnya lintang tempat dan sama dengan nol di khatulistiwa. Sering pula persamaan di atas ditulis:

a c = vf

dengan f = 2Ωsinφ

(2.3)

f dinamakan parameter coriolis. Arah percepatan coriolis ialah tegak lurus pada kecepatan, ke kanan di belahan bumi utara dan ke kiri di belahan bumi selatan. (Susilo Prawirowardoyo,1996).

2.3. Pendulum Foucault

Universitas Sumatera Utara

Pendulum Foucault adalah suatu alat yang berguna untuk menunjukkan arah rotasi bumi. Alat ini ditemukan oleh Jean Bernard Léon Foucault. Alat eksperimen ini terdiri atas bandul panjang yang bebas bergerak kesana kemari pada latar vertikal. Diagram dari pendulum Foucault diberikan pada gambar 2.2.

Gambar 2.2 pendulum Foucault

Pertunjukan pertama bandul Foucault kepada khalayak terjadi pada bulan Februari 1851 di Ruang Meridian yang ada di Observatorium Paris. Beberapa minggu kemudian, Léon Foucault membuat bandul terkenalnya ketika ia menggantung potongan logam seberat 28 kg dengan kabel sepanjang 67 meter dari kubah Panthéon di Paris. Hal itu untuk membuktikan bahwa bumi berputar pada porosnya. Pesawat osilasi dari pendulum Foucault berputar sepanjang hari sebagai akibat dari rotasi bumi. Pesawat osilasi menyelesaikan seluruh lingkaran dalam interval T, yang tergantung pada lintang geografis.

ini selanjutnya digunakan untuk menganalisis rotasi bumi pada porosnya, dengan bandul yang cukup berat yang digantungkan pada tali panjang, agar dapat bertahan berayunayun cukup lama. Oleh rotasi bumi, tentunya bandul akan berosilasi karena seolah-olah tertinggal oleh rotasi bumi. Untuk menyelidiki gerakan bandul yang memperlihatkan adannya rotasi itu, dapat tinjau bandul yang dipasang di tempatnya misalnya ?.(peter,2000)

Universitas Sumatera Utara

Gambar 2.3 Pendulum yang menunjukkan rotasi bumi.

pendulum yang terdiri dari partikel P dengan massa dan kawat (diperkirakan tak bermassa ) panjang, dan yang menyangga pada titik Q dekat lintang permukaan bumi seperti yang ditunjukkan pada gambar 2.3 untuk mempermudah, jarak vertikal h antara titik O ',yang mana pada permukaan bum, penyangga pada Q adalah diambil menjadi h=L. Skala dari gambar 2.3 terlalu dibesarkan dengan diasumsikan panjang L dari pendulum terlalu kecil dari radius RE di bumi, dengan bumi yang homogen sebagai massa M E . Hal ini juga diasumsikan bahwa gaya pada P adalah tegangan tali QP, besarnya N dan gaya gravitasi

bumi.

Besarnya gaya gravitasi yang

bekerja pada P ini diperkirakan sebagai GM E m / R 2 E = mg , dimana g = 9,81 m/s2 adalah percepatan gravitasi di permukaan bumi. Dan juga karena L sangat kecil dibandingkan dengan RE , gaya gravitasi yang bekerja pada P dinyatakan dalam arah − zˆ dimana zˆ adalah sebuah vektor satuan yang diarahkan ke atas sepanjang gatis vertikal OO’. Koordinar rectangular x, y, dan z ditunjukkan pada gambar 2.3 digunakan untuk menganalisis efek ? sebagai gaya pendulum yang bergerak relatif terhadap bumi. Istilah partikel yang diam di bumi dinyatakan dalam vektor satuan

( xˆ, yˆ, zˆ = xˆ × yˆ ) yang

r

ditunjukkan pada 2.3 dengan resultan gaya F dengan persamaan:

Universitas Sumatera Utara

r QP  xxˆ + yyˆ + ( z − L )zˆ  F ≈ − mgzˆ − N = −mgzˆ − N   L   QP

(2.4)

Dimana x, y dan z adalah koordinat rectangular dengan P adalah titik axis pada permukaan bumi dan titik tengah di O’. perhatikan bahwa, panjang pendulum L adalah konstan, variabel x, y dan z dinyatakan dengan persamaan:

x 2 + y 2 + ( z − L ) = L2 2

v

Dengan ? menjadi lokasi lintang titik O’ di bumi, dengan kecepatan angular ω dari

{xˆ, yˆ, zˆ} yang ditunjukkan pada gambar 2.3

r ω ≈ Ω (zˆ sin λ − xˆ cos λ )

(2.5)

dimana ? = 2p rad/hari ≈ 7,27.10-5 rad/s. Pada persamaan 2.5 rotasi garis dari matahari ke pusat O di bumi. Dengan kecepatan angular 2p rad/tahun. Diabaikan karena lebih kecil ? di sekitar titik pusat bumi. Dengan posisi vektor bumi dari O’ ke P dinyatakan sebagai

O' P = xxˆ + yyˆ + zzˆ dan di pusat O di bumi sekarang pada kerangka inersia. Posisi vektor sebenarnya dari P adalah:

r r p = RE zˆ + xxˆ + yyˆ + zzˆ

Universitas Sumatera Utara

r

r

Oleh karena itu, kecepatan υ p dan percepatan a p di P adalah:

r drp r r r υp = = x& xˆ + y& yˆ + z&zˆ + ω × rp dt

= ( x& − Ωy sin λ )xˆ + [ y& + Ωx sin λ + Ω(R E + z ) cos λ ]yˆ + ( z& − Ωy cos λ )zˆ

(2.6)

r r dυ p  dυ p  r r v  ap = =  + ω ×υ p  dt  dt  {xˆ , ˆy , ˆz}

[ ] ( + [&z& − 2Ωy& cos λ − Ω x (sin λ ) cos λ − Ω ( R + z ) cos λ ]zˆ

)

= &x& − 2Ωy& sin λ − Ω 2 ( RE + z )(sin λ ) cos λ xˆ + y&& + 2Ω x& sin λ + 2Ωz& cos λ − Ω 2 y yˆ 2

2

2

(2.7)

E

r

r

Dengan F dan a p sebenarnya ditunjukkan pada persamaan 2.4 dan 2.7 di bawah ini adalah

v

r

persamaan skalar yang diperoleh langsung dari F = ma p

[

]

) x m &x& − 2Ωy& sin λ − Ω 2 x sin 2 λ − Ω 2 ( RE + z )(sin λ ) cos λ = − N L

(

)

m &y& + 2Ωx& sin λ + 2Ωz& cos λ − Ω 2 y = −

[

(2.8)

y N L

(2.9)

]

m &z& − 2Ωy& cos λ − Ω 2 x (sin λ ) cos λ − Ω 2 (RE + z ) cos 2 λ = − mg −

z−L N L

(2.10)

Universitas Sumatera Utara

Persamaan 2.8 sampai 2.10 di subsitusikan ke persamaan 2.5 yang mempengaruhi gaya dari pendulu m relatif terhadap perputaran bumi. Efek dari ? sangat kecil pada gerakan pendulum dan sekarang dibahas secara rinci. dimana ? = 2p rad/hari ≈ 7,27.10-5 rad/s sangat kecil, perkiraan sekarang dibuat untuk mengurangi kompleksitas pada persamaan 2.8 sampai 2.10 tapi masih mempertahankan efek ? tentunya. Karena ayunan cukup kecil, gerakan pendulum boleh dikatakan berada pada bidang datar, maksudnya komponen gerakan ke atas dan ke bawahnya boleh diabaikan yang berarti z& dan &z& adalah 0. Maka persamaan pendulum Foucault adalah ditunjukkan pada persamaan 2.11 dan 2.12.

&x& − 2Ω sin λy& +

g x=0 L

(2.11)

&y& + 2Ω sin λx& +

g y =0 L

(2.12)

Persamaan 2.11 dan 2.12 diatas menghubungkan gerakan-gerakan sepanjang sumbu x dan sepanjang sumbu y yang berarti menentukan bentuk lintasan pendulum.(Marcelo da Silva,2004)

2.4. Metode Runge-Kutta

Salah satu metode numerik yang digunakan dalam penyelesaian persamaan differesial adalah metode Runge-Kutta. Metode ini mencapai ketelitian suatu pendekatan deret Taylor tanpa memerlukan kalkulasi turunan yang lebih tinggi. Banyak perubahan terjadi, tetapi semuanya dapat ditampung dalam bentuk umum dari persamaan 2.13 yi+1 = yi + f (xi , y i , h) h

(2.13)

dimana f (xi , yi , h) disebut suatu fungsi yang dapat diinterpretasikan sebagai sebuah slope rata-rata sepanjang interval. Fungsi tersebut dapat ditulis dalam bentuk umum dalam persamaan 2.14.

Universitas Sumatera Utara

f = a 1 k 1 + a 2 k 2 + … + an kn

(2.14)

dimana setiap a adalah konstanta dan setiap k besarnya adalah persamaan-persamaan 2.15.

k1 = f(xi , yi ) k2 = f(xi + p 1h, yi + q 11 k1 h) k3 = f(xi + p 2h, yi + q 21 k1 h + q 22k2h)

(2.15)

M kn = f(xi + p n-1h, yi + q n-1,1 k1 h + q n-1,2 k2 h + ...+ q n-1,n -1 kn-1h)

Semua harga k berhubungan secara rekurensi. Artinya k1 muncul dalam persamaan untuk k2, yang muncul lagi dalam persamaan untuk k3, dan seterusnya. Rekurensi ini membuat metode RK efisien untuk kalkulasi oleh komputer (Raymond et al, 1991).

Berbagai jenis metode Runge-Kutta dapat direncanakan dengan melaksanakan jumlah suku-suku yang berbeda pada fungsi tersebut seperti dinyatakan oleh n. untuk n = 1 atau RK orde pertama ternyata adalah metode Euler, yaitu persamaan 2.16. y1 = y0 + h f(x 0 ,y0 )

(2.16)

Dalam deret Taylor didapatkan persamaan 2.17.

y 0 = y(x 0 + h) = y 0 + h f(x0 , y 0 ) +

h2 f ' ( x0 , y 0 ) + ... 2!

(2.17)

Untuk metode RK orde kedua diberikan oleh persamaan-persamaan 2.18.

Universitas Sumatera Utara

k 1 = hf(x , y) 1 1   k 2 = hf  x + h, y + k 1  2 2   ∆y = k 2 , dengan h = ∆x

(2.18)

Metode RK orde tiga diberikan oleh persamaan-persamaan 2.19.

k 1 = hf(x , y) 1 1   k 2 = hf  x + h, y + k 1  2 2   k 3 = hf ( x + h, y + 2 k 2 − k1 ) 1 ∆y = ( k1 + 4k 2 + k 3 ) 6

(2.19)

Metode RK orde empat diberikan oleh persamaan-persamaan 2.20.

k 1 = hf(x , y) 1 1   k 2 = hf  x + h, y + k1  2 2   1 1   k 3 = hf  x + h, y + k 2  2 2   k 4 = hf ( x + h, y + k 3 ) 1 ∆y = ( k1 + 2k 2 + 2k 3 + k 4 ) 6 y ( x + h) = y ( x) + ∆y

(2.20)

Sedangkan untuk menyelesaikan persamaan differensial orde dua digunakan metode RK orde empat dengan terlebih dahulu membuat permisalan. Ditinjau persamaan differensial orde dua seperti pada persamaan 2.21.

d2y dy = f (x , y , ) 2 dx dx

(2.21)

Universitas Sumatera Utara

Dengan y(x 0 ) = y0, dan y’(x0 )= y0 ’ . Persamaan 2.19. dibuat permisalan sehingga diperoleh persamaan-persamaan 2.22.

dy = y′ = z dx dz = z ′ = y ′′ = f ( x, y , y ′) = f ( x , y , z ) dx

(2.22)

Persamaan-persamaan 2.16. merupakan persamaan-persamaan simultan yang dapat juga dituliskan sebagai f1 (x,y,z)=z dan f 2 (x,y,z)=f(x,y,z). Berdasarkan persamaan-persamaan 2.16 tersebut, persamaan differensial orde tersebut diselesaikan dengan mengikuti aturan metode RK orde empat pada persamaan 2.14 (Kandasamy et al,1997).

Contoh Runge kutta:

y " + xy ' + y = 0 dengan y (0 ) = 1; y ' (0 ) = 1; y (0,1) = y1 y " = − xy' − y = z ' y' = z dy = z = f 1 ( x, y, z ) dx

d2y = − xz − y = f 2 (x , y , z ) dx 2 Dengan memberikan nilai awal y 0 = 1, z 0 = y ' 0 = 0 maka:

k 1 = hf1 ( x0 , y 0 , z 0 ) = 0,1 f1 (0,1,0 ) = 0,1(0 ) = 0 l1 = hf 2 ( x0 , y 0 , z 0 ) = 0,1 f 2 (0,1,0 ) = 0,1(− (0)(0 ) − (1)) = 0,1(− 1) = −0,1 1 1 1   k 2 = hf 1  x 0 + h, y0 + k1 , z 0 + l1  2 2 2  

Universitas Sumatera Utara

1 1 1   = (0,1) f 1  0 + (0,1),1 + (0),0 + (− 0,1)  2 2 2   = 0,1 f 1 (0,05,1,0,005) k 2 = − 0,005

1 1 1   l2 = hf 2  x0 + h, y0 + k 1 , z 0 + l1  2 2 2   = 0,1 f 2 (0.05,1,− 0.005) l 2 = − 0, 09975

1 1 1   k 3 = hf 1  x 0 + h, y 0 + k 2 , z 0 + l 2  2 2 2   = 0,1 f1 (0.05,0.9975,− 0.0499) = -0.00499

1 1 1   l3 = hf 2  x 0 + h, y 0 + k 2 , z 0 + l2  2 2 2   = 0,1 f 2 (0.05,0.9975,− 0.0499 ) = -0.09950

k 4 = hf1 ( x0 + h, y 0 +k 3 , z 0 + l3 ) = 0,1 f1 (0.1,0.99511,−0.0995 ) = -0.00995

l4 = hf 2 (x 0 + h, y 0 + k3 , z 0 + l3 ) = 0,1 f 2 (0.1,0.99511, −0.0995 ) = -0.0985

Universitas Sumatera Utara

∴ y1 = y 0 + ∆y = 1+

1 [0 + 2(− 0.005 ) + 2(− 0.00499 ) − 0.00995] 2

= 0.9950

2.5 Mathematica 6

Mathematica adalah pemograman komputer yang dapat di kerjakan dengan matematika. Mathematica sering digunakan untuk instruksi, pekerjaan rumah, riset, dan penulisan. Mathematica cocok digunakan untuk numerik dan kemampuan dalam proses pemrograman sangat baik. Mathematica dapat juga dipakai sebagai pemodelan, perhitungan intensif, proyek disertasi, dan lain sebagainya. Elemen dasar :

A. Membuka paket Mathematica 6 a. Cari ikon Mathematica 6, kemudian klik dengan cepat dua kali, jendela kerja Mathematica 6 muncul.

b. Perintah (Command) dari Mathematica 6 di tuliskan dengan memakai symbol Matematika.

c. Setelah selesai menggunakan Mathematica 6, ketik quit dan [enter], atau klik file/EXIT

B. Operasi Mathematica 6

Universitas Sumatera Utara

Tabel 2.1 Operator Mathematica

Operator

Simbol

Contoh

Penambahan

+

2+3

Pengurangan

-

5-4

Perkalian

*

3*2

C. Variabel

Variabel pada Mathematica 6 harus diberi nama.Nama variabel harus dimulai dengan huruf, dan bisa

Pembagian

/

6/3

Pemangkatan

^

3^2

diikuti dengan huruf lain atau angka maksimum 31 karakter. Nama varibel dengan huruf besar (kapital)

dianggap berbeda dengan nama variabel yang ditulis dengan huruf kecil.

D. konstanta/tetapan

Beberapa tetapan yang berlaku pada Mathematica adalah sebagai berikut:

pi nilai p = 3,14452… eps nilai epsilon, bilangan natural e = 2.71828 inf nilai tak berhingga ~

E. Tanda baca

% Digunakan untuk mengawali komentar (command)

Universitas Sumatera Utara

, digunakan untuk memisahkan dua pernyataaan dalam sebaris ; digunakan untuk memisahkan dua pernyataan tanpa echo … digunakan untuk melanjutkan statemen ke baris berikutnya.

F. Kontrol Program

Mathematica 6 menyediakan beberapa instruksi yang memugkinkan pengguna membuat program atau fungsi, antara lain instruksi pemilihan (seleksi) dan instruksi perulangan (loop)

Instruksi seleksi: a. Pemilihan bersyarat: if (syarat-1) instruksi-1

else if instruksi-2 else instruksi-3 end.

Pemilihan diatas digunakan untuk memilih satu diantara beberapa instruksi sesuai dengan syarat yang dipenuhi. Bila syarat 1 dipenuhi maka laksanakan instruksi 1 , bila syarat 2 dipenuhi, maka laksanakan instruksi 2 bila tidak ada syarat yang dipenuhi maka laksanakan instruksi-3.

b. Pemilihan kasus switch variabel case 1 {nilai-1} instruksi-1 case 2 {nilai-2} instruksi-2 case 3 {nilai-3} instruksi-3 … Otherwise instruksi n End.

Universitas Sumatera Utara

Instruksi seleksi ini akan memilih satu instruksi berdasarkan nilai yang diberikan pada variabel. Bila nilainya adalah nilai -1 maka instruksi 1 dilaksanakan. Bila nilainya adalah 2 maka instruksi 2 yang dilaksanakan. Bila nilainya adalah 2 maka instruksi 2 yang dilaksanakan. Instruksi perulangan a. Perulangan dengan for for var = n1:n2:n3

instruksi- instruksi

end contoh: For[i=0,i<4,i++,Print[i]] 0 1 2 3 …………………………………………………………………………………… For[i=1;t=x,i^2<10,i++,t=t^2+i;Print[t]] 1+x2 2+(1+x^2)2 3+(2+(1+x^2)^2)2 Perulangan yang dibatasi oleh nilai var, mulai dari n1 hingga n3 dengan perubahan nilai sebesar n2 pada setiap putaran. Apakah n2=1 maka n2 tidak perlu ditulis, sehingga bentuknya menjadi: for var = n1: n3 instruksi- instruksi

Universitas Sumatera Utara

end b. Perulangan denganWhile while (syarat) insruksi- instruksi end Perulangan yang ditentukan oleh suatu syarat. Selama syarat terpenuhi maka perulangan akan berlangsung. (Wolfram,1991).

BAB 3

ANALISIS MASALAH DAN PERANCANGAN PROGRAM

3.1. Analisis Masalah

3.1.1. Persamaan Gerak Pendulum Foucault

Pada subbab ini akan dijelaskan persamaan pendulum Foucault yang akan diselesaikan yaitu persamaan 3.1 dan 3.2.

&x& − 2Ω sin λy& +

g x=0 L

(3.1)

&y& + 2Ω sin λx& +

g y =0 L

(3.2)

Dimana: x,y = koordinat titik pendulum bob dilihat dari permukaan bumi ?

= kecepatan sudut bumi (rad/s)

Universitas Sumatera Utara