ANALISA MODE SIMPANGAN AYUNAN BOLA BERONGGA BERISI …

1

ANALISA MODE SIMPANGAN AYUNAN BOLA BERONGGA BERISI FLUIDA MELALUI METODE SIMULASI DAN EKSPERIMEN.

Oleh,

Yusack Antonio Talangas

NIM : 192007029

TUGAS AKHIR

Diajukan kepada Program Studi Pendidikan Fisika, Fakultas Sains dan Matematika guna memenuhi sebagian dari persyaratan untuk mencapai gelar Sarjana Pendidikan

Program Studi Pendidikan Fisika

Fakultas Sains dan Matematika

Universitas Kristen Satya Wacana

Salatiga

2013

2

3

4

5

MOTTO dan PERSEMBAHAN

MOTTO

PERSEMBAHAN

Mereka berkata bahwa setiap orang membutuhkan tiga hal yang akan membuat

mereka berbahagia di dunia ini, yaitu; seseorang untuk dicintai, sesuatu untuk

dilakukan, dan sesuatu untuk diharapkan.

Dipersembahkan kepada :

Allah Bapa di Surga, Tuhan Yesus Kristus juru slamatku, Bunda Maria

Keluargaku tercinta, Bapak, Mama, Adik yang selalu memberikan dukungan, motivasi dan Doa untuk ku

Para pembaca

6

KATA PENGANTAR

’’ Adil ka talino bacuramin ka saruga basengat ka jubata’’ Puji syukur penulis panjatkan kepada Tuhan Yesus Kristus atas berkat dan perlindungan-Nya,

penulis dapat menyelesaikan skripsi ini dengan baik. Skripsi ini ditulis dan disusun untuk memenuhi sebagian persyaratan memperoleh gelar sarjana pendidikan (S.Pd) Fisika di Universitas Kristen Satya Wacana Salatiga.

Dalam penyusunan skripsi ini, tidak terlepas dari bantuan dan dukungan berbagai pihak. Atas segala bantuan dan dukungan tersebut, pada kesempatan kali ini penulis mengucapkan terima kasih kepada :

1. Tuhan Yesus Kristus 2. Keluargaku tercinta, Bapakku Alexander. D, Mamaku Syarifah Farida, Adikku Tirsa Marsalino,

Chandra, dan ayuk kalian merupakan inspirator dan motivator terbesar bagiku untuk menyelesaikan semua ini. terima kasih atas doa dan dukungan kalian.

3. Dr.Suryasatriya Trihandaru, M.Sc.nat yang mana merupakan Pembimbing Utama dan Nur Aji Wibowo, S.Si, M.Si sebagai Pembimbing pendamping. Terimakasih sudah meluangkan waktu, tenaga dan pikiran untuk memberikan masukan, dorongan, semangat serta dengan penuh kesabaran membimbing dan menuntunku selama penelitian hingga tulisan ini selesai.

4. Dosen-dosen Fisika (Ibu Dra. Marmi Sudarmi, M.Si, Bapak Dr. Suryasatriya Trihandaru, M.Sc.nat, Bapak Nur Aji Wibowo, S.Si, M.Si, Ibu Made Rai Suci, M.Pd, Ibu Diane Noviandini, S.Pd, Bapak Adita Sutrisno, S.Si., M.Sc, Bapak Andreas Setiawan, S.Si, MT, Ibu Debora Natalia Sudjito, S.Pd, Bapak Prof. Liek Wilardjo, Bapak Dr. Ferdy S. Rondonuwu, S.Pd., M.Sc) terima kasih telah mengajar dan mendidik, memberi bekal ilmu Pengetahuan.

5. Dekan Fakultas Sains dan Matematika bapak Dr. Suryasatriya Trihandaru, M.Sc.nat beserta jajarannya.

6. Pemerintah Kabupaten Landak yang telah memberikan bantuan berupa beasiswa. 7. Mas Tri, Mas Sigit dan Pak Tafip selaku laboran Fisika UKSW. Terima kasih untuk bantuannya

selama ini. 8. Pacarku tercinta terima kasih banyak atas doa, dukungan dan motivasinya. 9. Teman seperjuangan dalam suka maupun duka, pendidikan Fisika 2007 Kabupaten Landak (Angi,

Lia, Putri, Ica, Kx Devi, Kx Monic, Thamrin, Deo, Carles, Dodi, Wilson, Rodi, Apri, Aloy, Supri, Ota, Agus, Brama, Aska, Marius, Suwardi, Hengki ) terima kasih atas bantuan, dukungan dan perjuangan bersama selama ini.

10. Adik-adik, kakak angkatan Fisika dan Pendidikan Fisika, terima kasih atas dukungan dan bantuannya.

11. Semua pihak yang penulis tidak sebutkan satu persatu namanya yang turut dan terlibat dalam penyusunan tulisan ini terima kasih semua.

Akhirnya semoga tulisan ini bermanfaat dan menjadi berkat bagi pembaca khususnya bagi

pihak-pihak yang berkepentingan.

Salatiga, 2013

Yusack Antonio Talangas

7

ANALISA MODE SIMPANGAN AYUNAN BOLA BERONGGA BERISI FLUIDA

MELALUI METODE SIMULASI DAN EKSPERIMEN

Yusack Antonio Talangas, Nur Aji Wibowo, Suryasatriya Trihandaru 1 Program Studi Pendidikan Fisika, Fakultas Sains dan Matematika

2 Program Studi Fisika, Fakultas Sains dan Matematika Universitas Kristen Satya Wacana

Abstrak

Telah dilakukan penelitian gerak mode bandul sederhana yang dimodifikasi menjadi kulit bola yang setengahnya diisi oleh cairan. Gerak persamaan model ini ditentukan oleh persamaan lagrange dengan mengabaikan adanya gesekan udara. Untuk memecahkan persamaan lagrange dan simulasi gerakan dari model, MATLAB 7.0 digunakan dengan menerapkan skema Runge kutta numerik dan juga dilakukan eksperimen. Mode gerak pendulum dan permukaan air untuk variasi panjang log dan variasi radius kulit bola diperoleh. Sudut simpangan antara simpangan bandul, simpangan air, laju simpangan bandul, dan laju simpangan air berubah secara harmonik dan laju perubahan amplitudo menirun secara linear. Kata kunci : Bandul sederhana, MATLAB 7.0

I. PENDAHULUAN

Osilasi merupakan gerak bolak balik suatu benda yang melewati titik setimbang pada lintasan yang sama secara periodik. Osilasi menjadi salah satu jenis gerak yang sangat penting didalam ilmu fisika. Banyak kejadian fisika di alam yang berhubungan dengan gerak ini, misalnya osilasi roda keseimbangan arloji, dawai biola yang bergetar, massa yang diikatkan pada pegas lalu disimpangkan, atom dalam molekul atau dalam kisi zat padat, molekul udara ketika ada gelombang bunyi dan sebagainya [1].

Banyak penerapan osilasi yang sudah dipelajari, salah satunya adalah model bandul sederhana. Pada penelitian ini model bandul sederhana dikembangkan dengan memodifikasi bandul menjadi bola berongga yang diisi dengan air setengah dari volume bola berongga tersebut. Bola berongga tersebut direkatkan dengan sebatang kayu bermassa yang digantungkan, kemudian disimpangkan dari posisi setimbangnya dan dilepaskan, sehingga bandul berosilasi dalam bidang vertikal karena pengaruh medan gravitasi. Untuk mensimulasikan gerak bandul tersebut digunakan bahasa pemrograman MATLAB 7.0. Digunakannya bahasa ini dikarenakan kemampuannya yang tinggi dalam komputasi teknis. Keunggulan lain dari MATLAB adalah kemampuannya dalam menggabungkan komputasi, visualisasi, dan pemrograman dalam satu kesatuan yang mudah digunakan dengan notasi matematik yang sudah dikenal [2][6]. Selanjutnya dilakukan eksperimen untuk membandingkan hasil yang didapat dari simulasi dengan hasil dari eksperimen. Tujuan dari penelitian ini adalah untuk memperoleh persamaan gerak air maupun bola dan untuk

8

mendapatkan mode simpangan ayunan bola berongga berisi air dengan mengabaikan keberadaan gesekan udara. Implikasi dari penelitian ini diharapkan bisa bermanfaat sebagai pengaya khasanah ilmu dalam pembelajaran mekanika lanjut.

II. METODOLOGI

Ayunan matematis didalam penelitian ini dimodelkan dalam bentuk sebuah bola transparan dengan jari-jari sebesar R yang didalamnya diisi air dengan volume setengah volume bola. Massa bola dan massa air berturut-turut adalah mb dan ma Bola tersebut direkatkan pada kayu bermassa mk sepanjang l sehingga bola tidak bergerak relatif terhadap kayu saat diayunkan, seperti pada Gambar 1.

Gambar 1. Model ayunan

A. Koordinat Posisi

Jika pada Gambar 1 titik (0,0) merupakan titik pusat (origin), maka koordinat pusat massa dari kayu, bola, dan air dapat dinyatakan sebagai berikut

sin21lx (1)

y x x1 x2 x3

y1

y2 y3 z

R

l

Kayu

(0,0)

9

cos21ly (2)

sin2 Rlx (3) cos2 Rly (4)

sinsin3 zRlx (5) coscos3 zRly (6)

dengan simpangan bandul dan simpangan air

B. Energi Kinetik T Dan Energi Potensial U Jika bola tersebut disimpangkan dengan sudut simpangan ( ≈ 0,3 rad) yang cukup kecil, sehingga air tidak bergerak secara turbulen, dan air dapat diperlakukan sebagai benda tegar. Energi kinetik T dari sistem diberikan oleh

airbolakayu TTTT (7)

Karena pada energi kinetik kayu Tkayu, energi kinetik bola Tbola dan energi kinetik air Tair mengalami translasi dan rotasi, maka energi kinetik sistem secara berturut-turut diberikan oleh

222

21

221

kkkayu IlmT

(8)

222

21

21

bbbola IRlmT

(9)

22222

21cos2

21

aaair IzRlzRlmT (10)

dimana adalah laju simpangan bandul dan adalah laju simpangan air. Dengan momen inersia kayu Ik , momen inersia bola Ib dan momen inersia air Ia berturut-turut sebesar

2

121 mlI k (11)

2

32 RmI bb (12)

22

52 zmRmI aaa (13)

Untuk energi potensial U dari sistem diberikan oleh airbolakayu UUUU (15)

dengan Ukayu diberikan oleh

Bola

Air

10

cos2lgmU kkayu (16)

dengan Ubola diberikan oleh

cos2 RlgmUbola (17) dengan Uair diberikan oleh

coscos zRlgmU aair (18) dan pada persamaan (10), (13), dan (18) dipakai defenisi (Lampiran A)

Rz83

(14)

C. Lagrange

Persamaan Lagrange L ditentukan atas dasar energi dari sistem baik energi kinetik, energi potensial maupun energi redaman, disamping gaya luar yang bekerja. Dengan menggunakan persamaan T dan U dari persamaan (7) dan (15) diperoleh Lagrangian L sebagai berikut [3][5].

UTL (19)

222222

31

21

61 RmRlmlmL bbk

2222222

21

51cos2

21

zmRmzRlzRlm aaa

coscoscoscos21 zRlgmRlgmglm abk (20)

D. Euler-Lagrange

Kemudian persamaan L disubstitusikan ke dua persamaan Euler-Lagrange sebagai berikut [4][7].

1

LL

dtd

(21)

2

LL

dtd (22)

Pada ruas kanan persamaan (21) dan (22) telah diperhitungkan gaya gesekan dengan asumsi bahwa besar gaya gesekan sebanding dengan kecepatan sudut dengan arah berlawanan terhadap arah simpangan. Didalam penelitian ini, diasumsikan bahwa gesekan dengan udara sangat kecil sehingga dapat diabaikan 1=0 dan gesekan antara air dan dinding bola tidak diabaikan, sehingga 1 ≠0 Dengan menyelesaikan persamaan (21) dan (22) didapatkan hasil sebagai berikut

11

21 hh (24) dan

21 gg (25) Dan pada persamaan (24) dan (25) dipakai defenisi (Lampiran B).

E. SKEMA NUMERIK RUNGE-KUTTA Skema numerik Runge-Kutta dipakai untuk memudahkan dalam membuat bahasa program pada MATLAB 7.0, adapun skema numerik untuk persamaan (24) dan (25) adalah sebagai berikut

,,,,1 tf

,,,,2 tf

Kemudian dibuat vektor U

dengan elemen sebagai berikut:

4321

UUUU

phidotphi

thetadottheta

U

Sehingga diperoleh

,,,,

,,,,

,,,,

2

1 tF

tf

tfdtUd

dtd

dtd

dtd

dtd

Kemudian membuat program menggunakan bahasa pemrograman MATLAB 7.0 menggunakan metode Runge-kutta. yang nantinya akan dicari untuk Mode simpangan 1. Variasi panjang kayu (l = 0,18 m, l = 0,28 m, dan l = 0,38 m) dengan jari-jari bola (R)

0,18 m. 2. Variasi jari-jari bola (R = 0,18 m, R = 0,28 m, dan R = 0,38 m) dengan panjang kayu

(l) 0,18 m Selanjutnya dilakukan eksperimen dengan cara memvideokan gerak ayunan bandul kemudian video tersebut diekstrak dan dicari pola grafik nya, selanjutnya grafik yang didapat akan digunakan untuk perbandingan dengan grafik hasil dari simulasi.

III. HASIL DAN ANALISA

12

Dengan menyelesaikan persamaan (21) dan (22), diperoleh mode ketergantungan waktu dari simpangan bandul dan laju simpangan bandul seperti ditunjukan pada Gambar 2 dan 3.

0 2 4 6

-0.2

-0.1

0

0.1

0.2

0.3

0 2 4 6

-0.2

-0.1

0

0.1

0.2

0.3

(a) (b) Gambar 2. Grafik t untuk (a). Variasi panjang kayu.

(b). Variasi jari-jari bola.

0 2 4 6

-1

0

1

0 2 4 6

-1

0

1

(a) (b) Gambar 3. Grafik t untuk (a). Variasi panjang kayu.

(b). Variasi jari-jari bola.

Gambar 2(a) mendeskripsikan mode simpangan untuk tiga panjang kayu yang berbeda, sedangkan Gambar 2(b) mendeskripsikan mode simpangan untuk tiga ukuran bola yang berbeda. Dari Gambar 2(a-b) tersebut terlihat bahwa kayu-bola bergerak secara harmonik untuk sudut simpangan yang kecil (θ ≈ 0,3 rad).

(r

adia

n)

(r

adia

n)

t (s) t (s)

t (s) t (s)

13

Dari hasil praktikum diperoleh juga mode ketergantungan waktu dari simpangan bandul dan laju simpangan bandul . Seperti ditunjukan pada Gambar 3 dan 4. Dimana mode tersebut hampir mendekati mode yang didapat pada teori.

0 1 2 3

-0.2

0

0.2

3 4 5 6-0.3

-0.2

-0.1

0

0.1

0 1 2 3

-0.3

-0.2

-0.1

0

0.1

0.2

a b c Gambar 4. Grafik t pada variasi panjang kayu untuk.

a. l=0.18 m. b. l=0.28 m. c. l=0.38 m

0 1 2 3

-0.2

0

0.2

0 1 2 3-0.4

-0.2

0

0.2

0.4

4 5 6 7-0.4

-0.3

-0.2

-0.1

0

0.1

a b c Gambar 5. Grafik t pada variasi jari-jari bola untuk.

a. l=0.055 m. b. l=0.065 m. c. l=0.075 m

Terdapat kekhususan dari mode gerak tersebut yaitu bahwa amplitudo simpangan cenderung menurun seiring dengan bertambahnya waktu seperti ditunjukan pada Gambar 6 Dan semakin besar panjang kayu maupun jari-jari bola diperoleh bahwa laju perubahan amplitudo menurun setiap saat.

2 4 6 80.1

0.15

0.2

0.25

0.3

2 4 6 80.1

0.15

0.2

0.25

0.3

max

(rad

ian)

t (s) t (s)

max

(rad

ian)

(r

adia

n)

(r

adia

n)

(r

adia

n)

t (s) t (s) t (s)

(r

adia

n)

(r

adia

n)

(r

adia

n)

t (s) t (s) t (s)

14

(a) (b) Gambar 6. Grafik max t untuk (a). Variasi jari-jari bola.

(b). Variasi panjang kayu.

Gambar 7 mendeskripsikan keterkaitan periode terhadap waktu. Gambar 7(ab) menunjukkan bahwa periode dari ketiga nilai panjang kayu dan jari-jari bola berubah setiap saat dengan tingkat perubahan maksimum 2.55%.

1 2 3 4 5 6 70.96

0.97

0.98

0.99

1 2 3 4 5 6 70.96

0.97

0.98

0.99

1

(a) (b) Gambar 7. Grafik T Osilasi ke-n untuk (a). Variasi jari-jari bola.

(b). Variasi panjang kayu.

Didalam gerak ayunan ini besarnya kecepatan simpangan maksimum dari kayu menurun seiring dengan bertambahnya waktu dan laju perubahan kecepatan maksimum ini juga sangat dipengaruhi oleh panjang kayu dan jari-jari bola seperti yang terlihat pada Gambar 8 untuk batang kayu yang semakin panjang dan diameter bola yang semakin besar, perubahan kecepatan tetap sama.

2 4 6 80.5

1

1.5

1 2 3 4 5 6 7

0.8

1

1.2

1.4

1.6

(a) (b)

T (s

)

T (s

)

Osilasi ke-n Osilasi ke-n

t (s) t (s)

15

Gambar 8. Grafik max t untuk (a). Variasi jari-jari bola. (b). Variasi panjang kayu.

Hal lain yang menarik untuk diteliti dari sistem ini adalah bagaimana gerakan permukaan air yang berada didalam kulit bola berubah terhadap waktu. Pola gerakan tersebut diperoleh dengan menyelesaikan persamaan (21) dan (22) sehingga didapat mode ketergantungan waktu dari simpangan air dan laju simpangan air seperti ditunjukkan pada Gambar 9 dan 10

0 2 4 6

-0.2

-0.1

0

0.1

0.2

0 2 4 6

-0.2

-0.1

0

0.1

0.2

(a) (b) Gamabr 9. Grafik t untuk (a). Variasi panjang kayu.

(b). Variasi jari-jari bola.

0 2 4 6

-1

0

1

0 2 4 6

-1

0

1

(a) (b) Gambar 10. Grafik t untuk (a). Variasi panjang kayu.

(b). Variasi jari-jari bola. Dari hasil praktikum didapat mode ketergantungan waktu dari simpangan air dan laju

simpangan air seperti ditunjukkan pada Gambar 11 dan 12.

t (s)

ϕ (r

adia

n)

ϕ (r

adia

n)

t (s) t (s)

t (s)

(rad

ian)

(rad

ian)

16

0 1 2 3

-0.5

0

0.5

3 4 5 6

-0.5

0

0.5

0 1 2 3

-0.5

0

0.5

a b c Gambar 11. Grafik t pada variasi panjang kayu untuk.

a. l=0.18 m. b. l=0.28 m. c. l=0.38 m

0 1 2 3

-0.5

0

0.5

0 1 2 3

-0.5

0

0.5

4 5 6 7

-0.5

0

0.5

a b c Gambar 12. Grafik t pada variasi jari-jari bola untuk.

a. R=0.055 m. b. R=0.065 m. c. R=0.075 m

Gambar 9(a) menunjukan mode simpangan untuk tiga bola yang berbeda jari-jarinya, sedangkan Gambar 9(b) menunjukan mode simpangan untuk tiga kayu yang berbeda panjangnya. Dari Gambar 9(a-b) tersebut terlihat bahwa air bergerak secara harmonik. Kekhususan yang terdapat dari mode gerak tersebut adalah bahwa amplitudo simpangan mengalami peningkatan diwaktu pertama dan setelah itu turun setiap saat seiring dengan bertambahnya waktu seperti pada Gambar 13 dengan semakin besar panjang kayu dan diameter bola diperoleh bahwa laju perubahan amplitudo meningkat secara linear.

2 4 6 80.1

0.15

0.2

2 4 6 80.1

0.15

0.2

(a) (b)

max

(rad

ian)

max

(rad

ian)

t (s) t (s)

(r

adia

n)

(r

adia

n)

(r

adia

n)

t (s) t (s) t (s)

(r

adia

n)

(r

adia

n)

(r

adia

n)

t (s) t (s) t (s)

17

Gambar 13. Grafik max t untuk (a). Variasi jari-jari bola. (b). Variasi panjang kayu.

Gambar 14 menunjukan keterkaitan periode setiap waktunya. Gambar 14(a-b) menunjukan

bahwa periode untuk perbedaan besarnya panjang kayu dan jari-jari bola berubah dengan tingkat perubahan maksimum 7,71% .

1 2 3 4 5 6 7

0.92

0.94

0.96

0.98

1

1 2 3 4 5 6 70.92

0.94

0.96

0.98

1

(a) (b) Gambar 14. Grafik T Osilasi ke-n untuk (a). Variasi jari-jari bola.

(b). Variasi panjang kayu

Didalam gerak ayunan ini besarnya kecepatan maksimum dari air juga meningkat disaat pertama dan selanjutnya menurun setiap saat dan laju perubahan kecepatan maksimum ini juga sangat dipengaruhi oleh panjang kayu dan jari-jari bola seperti terlihat pada Gambar 15.

2 4 6 8

0.8

1

1.2

1.4

1.6

2 4 6 8

0.8

1

1.2

1.4

1.6

(a) (b)

T (s

)

Osilasi ke-n

T (s

)

Osilasi ke-n

t (s)

t (s)

(rad

ian)

(rad

ian)

18

Gambar 15. Grafik max t untuk (a). Variasi jari-jari bola. (b). Variasi panjang kayu

Sedangkan keterkaitan periode terhadap waktu terdeskripsikan pada Gambar 16. menunjukan bahwa periode cenderung sama untuk osilasi ke 2 dan seterusnya.

1 2 3 4 5 6 7

0.8

0.9

1 2 3 4 5 6 70.8

0.9

1

(a) (b)

Figure 16. Graph T Osilasi ke-n untuk (a). Variasi jari-jari bola. (b). Variasi panjang kayu

IV. KESIMPULAN

Telah dilakukan studi tentang pengembangan dari model ayunan sederhana dengan memodifikasi bandul menjadi bola berongga yang di isi air dengan volume setengah dari volume bola tersebut. Dengan menyelesaikan persamaan gerak dari model bandul tersebut dan mensimulasikan geraknya dengan menggunakan program MATLAB 7.0 serta melakukan eksperimen untuk membuktikan teori yang ada, didapat mode gerak dari bandul dan permukaan air. mode simpangan θ, , ϕ, dan bergerak secara harmonik dan laju

perubahan amplitudonya menurun. Diperoleh hasil yang menarik untuk dan terdapat perubahan yang terjadi setiap saat

V. DAFTAR PUSTAKA

[1]. Resnick , Halliday. 1985. Fisika Jilid 1. Jakarta: Erlangga. [2]. MATLAB Bahasa Komputasi Teknis, Penerbit ANDI Yogyakarta 2000 [3]. M. Z. Rafat. 2006. Dynamics of a double pendulum with distributed mass. Tesis tidak diterbitkan.

School of Physics, University of Sydney, Australia. [4]. Single and Double plane pendulum oleh: Gabriela Gonz´alez [5]. Matthew, West. 2004. Variational Integrators. California Institute of Technology, Pasadena,

California. [6]. A MATLAB Tutorial. Ed Overman. Department of Mathematics, The Ohio State University

T (s

)

Osilasi ke-n Osilasi ke-n

T (s

)

19

[7]. The Euler-Lagrange Equation: When Necessary is Not Sufficient. Shawn D. Ryan. The Pennsylvania State University