BAB 2 TINJAUAN PUSTAKA

BAB 2 TINJAUAN PUSTAKA

2.1.Dasar Fluida Dalam buku yang berjudul “Fundamental of Fluid Mechanics” karya Bruce R. Munson, Donald F. Young, Theodore H. Okiishi, dan Wade W. Huebsch, fluida didefinisikan sebagai zat yang berdeformasi terus-menerus selama dipengaruhi suatu tegangan geser. Sebuah tengangan (gaya per satuan luas) geser terbentuk apabila sebuah gaya tangensial bekerja pada sebuah permukaan. Apabila bendabenda padat biasa seperti baja atau logam-logam lainnya dikenai oleh suatu tegangan geser, mula-mula benda ini akan berdeformasi (biasanya sangat kecil), tetapi tidak akan terus-menerus berdeformasi (mengalir). Namun, cairan yang biasa seperti air, minyak, dan udara memenuhi definisi dari sebuah fluida artinya, zat-zat tersebut akan mengalir apabila padanya bekerja sebuah tegangan geser. Beberapa bahan, seperti lumpur, aspal, dempul, odol dan lain sebagainya tidak mudah untuk diklasifikasikan karena bahan-bahan tersebut akan berperilaku seperti benda padat jika tegangan geser yang bekerja kecil, tetapi jika tegangan tersebut melampaui suatu nilai kritis tertentu, zat-zat tersebut akan mengalir. Ilmu yang mempelajari bahan-bahan tersebut disebut rheology dan tidak termasuk dalam cakupan mekanika fluida klasik.

2.2 Bilangan Reynolds

Joseph et al. (1996) dalam bukunya menjelaskan tentang bilangan Reynolds dimana jika diperhatikan gerak dinamis dari aliran kental dengan skala kecepatan 𝑈 dan skala panjang 𝐿. Dua parameter cairan yang paling penting yang mempengaruhi gerak adalah 𝜌 kepadatan dan viskositas 𝜇. Empat parameter ini (𝑈, 𝐿, 𝜌, 𝜇) dapat dikombinasikan ke dalam kelompok berdimensi tunggal yang disebut bilangan Reynolds atau sering dituliskan 𝑅𝑒 (Osborne Reynolds (1883))

6

𝑅𝑒 =

𝜌. 𝐿. 𝑈 𝐿. 𝑈 = 𝜇 𝑣

(1)

𝜇 di mana 𝑣 = ⁄𝜌 adalah rasio nyaman yang disebut viskositas kinematik fluida. Bilangan Reynolds adalah parameter dominan yang mempengaruhi hampir semua arus kental.

2.3 Laminar dan Turbulent Flow

Pentingnya jumlah Reynolds dengan indah digambarkan dalam percobaan klasik oleh Reynolds sendiri, menggunakan zat warna yang beruntun untuk memvisualisasikan aliran melalui pipa halus, seperti pada Gambar. 2.1. Jika jumlah Reynolds rendah, zat warna yang beruntun tetap lurus dan halus [Gambar. 2.1 (sebagai kondisi yang disebut laminar atau merampingkan aliran. Dalam Reynolds berbagai jumlah menengah [Gambar. 2.1 (b)], zat warna yang beruntun memperlihatkan perilaku yang tidak menentu, dan pengukuran titik, katakanlah, kecepatan terhadap waktu menunjukkan tidak beraturan "semburan" aktivitas . rentang peralihan ini disebut aliran transisi.

Pada nomor Reynolds tinggi [Gambar. 2.1 (c)], istirahat pewarna beruntun dan campuran pada tingkat intens, mengisi tabung dengan warna. Pengukuran kecepatan titik menunjukkan fluktuasi acak kontinu disebut turbulensi, dan arus sesaat menjalin seperti spaghetti. Ini adalah aliran turbulen, dan memiliki gesekan dan panas transfer cukup karakter yang berbeda dibandingkan dengan aliran laminar.

7

Gambar 2.1. Visualisasi zat warna yang beruntun dan pengukuran kecepatan aliran saluran (setelah percobaan terkenal oleh Osborne Reynolds pada tahun 1883): (a) aliran laminar, Re rendah, (b) aliran transisi, Re moderat, dan (c) aliran turbulen, Re besar.

2.4.Persamaan Navier-Stokes

Persamaan Navier-Stokes adalah dasar persamaan differensial parsial yang menguraikan aliran fluida yang tak dapat dimampatkan. Dengan menggunakan tingkat tekanan dan tingkat tegangan tensor. Hal ini, dapat ditunjukkan dari persamaan Fj sebagai bagian komponen kekuatan merekat dari F pada suatu wadah yang tak berputar yaitu sebagai berikut 𝐹𝑖 𝜕 𝜕𝑢𝑖 𝜕𝑢𝑗 = [𝜂 ( + + 𝜆𝛿𝑖𝑗 ∇. 𝒖)] 𝑉 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑖

(2)

dimana 𝜂 adalah kecepatan dinamik, 𝜆 adalah koefisien kecepata kedua, 𝛿𝑖𝑗 adalah Kronecker delta, ∇. 𝒖 adalah divergen (Tritton 1988, Faber 1995).

8 Dalam buku yang berjudul “Fundamental of Fluid Mechanics” karya Bruce R. Munson, Donald F. Young, Theodore H. Okiishi, dan Wade W. Huebsch, persamaan Navier-Stokes adalah persamaan diferensial dasar yang menggambarkan aliran fluida Newtonian. Suatu persamaan tegangan dapat disubitusikan terhadap persamaan differensial untuk benda yang bergerak yakni : 𝜌𝑔𝑥 +

𝜕𝜎𝑥𝑥 𝜕𝜏𝑦𝑥 𝜕𝜏𝑧𝑥 𝜕𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑢 + + = 𝜌( +𝑢 +𝑣 +𝑤 ) 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

(3)

𝜌𝑔𝑦 +

𝜕𝜎𝑥𝑦 𝜕𝜏𝑦𝑦 𝜕𝜏𝑧𝑦 𝜕𝑣 𝜕𝑣 𝜕𝑣 𝜕𝑣 + + = 𝜌( +𝑢 +𝑣 +𝑤 ) 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

(4)

𝜌𝑔𝑧 +

𝜕𝜎𝑥𝑧 𝜕𝜏𝑦𝑧 𝜕𝜏𝑧𝑧 𝜕𝑤 𝜕𝑤 𝜕𝑤 𝜕𝑤 + + = 𝜌( +𝑢 +𝑣 +𝑤 ) 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

(5)

dan disederhanakan menggunakan persamaan kontinuitas 𝜕𝑢 𝜕𝑥

𝜕𝑣

+ 𝜕𝑦 = 0

(6)

sehingga diperoleh : (terhadap x) 𝜕𝑢

𝜕𝑢

𝜕𝑢

𝜕𝑢

𝜕𝑝

𝜕2 𝑢

𝜕2 𝑢

𝜕2 𝑢

𝜕𝑢

𝜕𝑢

𝜕𝑝

𝜕2 𝑣

𝜕2 𝑣

𝜕2 𝑣

𝜕𝑢

𝜕𝑢

𝜕𝑝

𝜕2 𝑤

𝜕2 𝑤

𝜕2 𝑤

𝜌 ( 𝜕𝑡 + 𝑢 𝜕𝑥 + 𝑣 𝜕𝑦 + 𝑤 𝜕𝑧 ) = − 𝜕𝑥 + 𝜌𝑔𝑥 + 𝜇 (𝜕𝑥 2 + 𝜕𝑦 2 + 𝜕𝑧 2 )

(7)

(terhadap y) 𝜕𝑢

𝜕𝑢

𝜌 ( 𝜕𝑡 + 𝑢 𝜕𝑥 + 𝑣 𝜕𝑦 + 𝑤 𝜕𝑧 ) = − 𝜕𝑦 + 𝜌𝑔𝑦 + 𝜇 (𝜕𝑥 2 + 𝜕𝑦 2 + 𝜕𝑧 2 )

(8)

(terhadap z) 𝜕𝑢

𝜕𝑢

𝜌 ( 𝜕𝑡 + 𝑢 𝜕𝑥 + 𝑣 𝜕𝑦 + 𝑤 𝜕𝑧 ) = − 𝜕𝑧 + 𝜌𝑔𝑧 + 𝜇 ( 𝜕𝑥 2 + 𝜕𝑦 2 +

𝜕𝑧 2

)

(9)

dimana u, v dan w adalah komponen-komponen kecepatan dari x, y dan z. Dapat dilihat telah disusun kembali persamaan tersebut di mana terminologi percepatan ditunjukkan pada sisi sebelah kiri dan terminologi ketegangan sebelah kanan. Persamaan-persamaan tersebut secara umum disebut dengan persamaan Navier-Stokes yang diambil dari nama ahli matematika Francis L. M. H. Navier (1785-1836) dan mekanik Inggris Bapak G. G. Stokes (18191903). Ketiga persamaan tersebut ketika dikombinasikan dengan persamaan kekekalan massa (persamaan kontinuitas), memperlihatkan uraian matematika yang lengkap dari suatu aliran fluida Newtonian tak termampatkan. Diperoleh empat persamaan dan empat tak diketahui (u, v, w dan p) dan oleh karena itu

9 masalahnya adalah mana “yang baik diambil” pada sifat-sifat matematika. Sayangnya, karena kompleksitas umum dari persamaan Navier-Stokes (yaitu nonlinier, tingkat-kedua, persamaan differensial parsial) kompleksitas tersebut tidak dapat dikerjakan dengan penyelesaian yang sangat baik kecuali pada beberapa permasalahan. Namun, pada beberapa permasalahan yang solusinya telah diperoleh dan dibandingkan dengan hasil eksperimen, ternyata hasilnya hampir dapat diterima. Oleh karena itu, persamaan Navier-Stokes dibuat sebagai pendekatan persamaan differensial untuk fluida Newtonian tak termampatkan.

Dari sisi koordinat polar silinder (tabung), persamaan Navier-Stokes dapat ditulis sebagai : (terhadap r) 𝜕𝑣

𝜌 ( 𝜕𝑡𝑟 + 𝑣𝑟 1 𝜕2 𝑣𝑟 𝑟 2 𝜕𝜃2

𝜕𝑣𝑟 𝜕𝑟

2 𝜕𝑣𝜃

− 𝑟2

𝑣𝜃 𝜕𝑣𝑟

+

𝑟 𝜕𝜃 𝜕2 𝑣𝑟

+

𝜕𝜃

𝜕𝑧 2

−

𝑣𝜃2 𝑟

+ 𝑣𝑧

𝜕𝑣𝑟 𝜕𝑧

𝜕𝑝

1 𝜕

) = − 𝜕𝜃 + 𝜌𝑔𝑟 + 𝜇 [𝑟 𝜕𝑟 (𝑟

𝜕𝑣𝑟 𝜕𝑟

𝑣

) − 𝑟𝑟2 +

]

(10)

(terhadap 𝜃) 𝜕𝑣

𝜕𝑣𝜃

𝜌 ( 𝜕𝑡𝜃 + 𝑣𝑟

𝜕𝑟

1 𝜕2 𝑣𝜃

𝑣𝜃

+ 𝑟2 𝑟2

+

𝑣𝜃 𝜕𝑣𝜃 𝑟 𝜕𝜃

2 𝜕𝑣𝑟

𝑣𝑟 𝑣𝜃

𝜕2 𝑣𝜃

+ 𝑟2

𝜕𝜃

𝜕𝑣𝑧

𝑣𝜃 𝜕𝑣𝑧

𝜕𝜃2

+

+

𝜕𝑧 2

𝑟

+ 𝑣𝑧

𝜕𝑣𝜃

1 𝜕𝑝

1 𝜕

) = − 𝑟 𝜕𝑥 + 𝜌𝑔𝜃 + 𝜇 [𝑟 𝜕𝑟 (𝑟 𝜕𝑧

]

𝜕𝑣𝜃 𝜕𝑟

)−

(11)

(terhadap z) 𝜕𝑣

𝜌 ( 𝜕𝑡𝑧 + 𝑣𝑟 2 𝜕𝑣𝜃 𝑟 2 𝜕𝜃

+

𝜕𝑟

𝜕2 𝑣𝑧 𝜕𝑧 2

+

𝑟 𝜕𝜃

+ 𝑣𝑧

]

𝜕𝑣𝑧

𝜕𝑝

1 𝜕

) = − 𝜕𝑧 + 𝜌𝑔𝑧 + 𝜇 [𝑟 𝜕𝑟 (𝑟 𝜕𝑧

𝜕𝑣𝑧

1 𝜕2 𝑣𝑧

) + 𝑟2 𝜕𝑟

𝜕𝜃2

−

(12)

2.5.Metode Elemen Hingga Hidayat (2005) dalam bukunya yang berjudul “Teori dan Penerapan Metode Elemen Hingga”, disampaikan bahwa Metode Elemen Hingga merupakan prosedur numerik yang diterima secara luas untuk menyelesaikan persamaan differensial dalam teknik dan fisika. Metode ini menjadi dasar komputasional

10

dari system computer untuk perancangan. MEH mempunyai kemampuan yang sangat baik untuk menyelesaikan persoalan transien satu dimensi dan dua dimensi.

Metode Elemen Hingga (MEH) merupakan prosedur numerik untuk menyelesaikan permasalahan fisik yang diatur dengan persamaan diferensial atau teorema energi. Karakteristik MEH yang membedakan dengan prosedur numerik lainnya adalah : 1. MEH menggunakan penyelesaian integral untuk menghasilkan sistem persamaan aljabar 2. MEH menggunakan fungsi-fungsi kontinu sebagian (continuous piecewise smooth functions) untuk mendeteksi kuantitas atau beberapa kuantitas yang tidak diketahui

Secara umum MEH terdiri dari lima langkah dasar : 1. Mendiskritisasikan daerah-daerah yang meliputi langkah-langkah penempatan titik-titik nodal, penomoroan titik-titik nodal dan penentuaan koordinatnya. 2. Menentukan derajat atau orde persamaan pendekatan : linear atau kuadratik. Persamaan harus dinyatakan sebagai fungsi nodal. Persamaan ditentukan untuk tiap elemen. 3. Menyusun system persamaan-persamaan. 4. Menyelesaikan system persamaan-persamaan. 5. Menghitung kuantitas yang dicari. Kuantitas dapat merupakan komponen tegangan, heat flow, fluid velocities, dan lain-lain.

Persamaan dalam MEH biasanya berbentuk : [𝑘]{𝑢} = {𝐹}

(13)

dengan [𝑘] merupakan matrik bujur sangkar yang disebut matrik kekakuan, {𝑢} merupakan vector kolom dengan komponen matrik berupa nilai nodal yang tidak diketahui. Nilai nodal dapat berupa simpangan atau temperature,

11 sedangkan {𝐹} berupa matrik kolom yaitu ‘gaya’ yang bekerja pada nodal. Gaya dapat berupa 𝐹 (gaya) atau 𝑄 (kalor).

Dalam menyelesaikan masalah fisik yang berhubungan dengan persamaan differensial, cara terbaiknya adalah : 1. Mencari solusi analitisnya. Pada banyak kondisi, solusi analitis sulit diperoleh, sehingga digunakan metode numerik untuk mencari solusi pendekatannya. 2. Beberapa prosedur untuk mendapatkan penyelesaiann persamaan differensial dengan metode numerik adalah : a. Metode beda hingga b. Metode varisional c. Metode Residual Berat Dari ketiga metode tersebut, akan menggunakan metode residual berat yaitu metode Galerkin.

2.6. Elemen Linier 1 Dimensi

Pada bagian ini akan dibahas pembagian daerah satu dimensi menjadi elemenelemen linier dan mengembangkan persamaan untuk satu elemen. Persamaan elemen ini digeneralisasi untuk memperoleh persamaan kontinu sebagian untuk daerah satu dimensi tersebut. Daerah satu dimensi merupakan segmen garis atau suatu garis. Pembagian segmen garis menjadi elemen-elemen yang lebih kecil dengan menggunakan nodal. Ketentuan untuk elemen dan nodal adalah : 1. Nomor nodal dengan urutan dari kiri ke kanan 2. Nomor elemen dengan urutan dari kiri ke kanan; di dalam tanda kurung ( -- ) Sedangkan ketentuan penempatan nodal : 1. Tempatkan nodal-nodal dengan lebih rapat pada daerah di mana parameter yang tidak diketahui berubah dengan cepat dan tempatkan

12

nodal-nodal secara berjauhan jika unknown parameter nya konstan atau relative konstan. 2. Tempatkan nodal di manapun terdapat perubahan nilai koefisien D dan Q. 3. Tempatkan nodal di manapun jika diinginkan mengetahui nilai 𝜙

Elemen linier 1 dimensi adalah garis dengan panjang L dengan nodal pada ujung-ujungnya. Nodal dinyatakan dengan I dan j dan nilai nodal dengan 𝜙i dan 𝜙j . Elemen linier 1 dimensi ditunjukkan pada Gambar 2.2.

𝜙 = 𝑎1 + 𝑎2 𝑥

𝜙2 𝜙1

L 𝑥𝑖 𝑥𝑗

Gambar 2.2. Elemen Linier Parameter 𝜙 berubah secara linier antara nodal i dan j. Persamaan 𝜙 adalah : 𝜙 = 𝑎1 + 𝑎2 𝑥

(14)

Koefisien a1 dan a2 ditentukan dari nilai kondisi nodal : 𝜙 = Φi

di

x = Xi

𝜙 = Φj

di

x = Xj sehingga diperoleh

Φi = a1 + a2Xi

dan Φj = a1 + a2Xj

(15)

Eliminasi persamaan (15), maka dapat diperoleh 𝑎1 =

Φ𝑖 𝑋𝑗 −Φ𝑗 X𝑖 𝑋𝑗 −𝑋𝑖

dan

𝑎2 =

Φ𝑗 −Φ𝑖 𝑋𝑗 −𝑋𝑖

(16)

13

Substitusi persamaan (16) ke (14) diperoleh : 𝑋𝑗 −𝑥

𝜙=(

𝐿

𝑥−𝑋𝑖

) Φ𝑖 + (

𝐿

) Φ𝑗

(17)

Dengan L = Xj - Xi Persamaan (17) adalah bentuk fungsi interpolasi elemen hingga standar. Fungsi linear x pada persamaan (17) adalah fungsi bentuk yang dinyatakan dengan N dan tanda indeks yang sesuai dengan nodalnya. Fungsi bentuk pada persamaan (17) dinyatakan dalam Ni dan Nj sebagai berikut :

𝑁𝑖 =

𝑋𝑗 − 𝑥 𝐿

dan

𝑁𝑗 =

𝑥− 𝑋𝑖 𝐿

Sehingga dapat ditulis : 𝜙 = 𝑁𝑖 Φ𝑖 + 𝑁𝑗 Φ𝑗 dan dinyatakan dalam bentuk persamaan matrik sebagai : 𝜙 = [𝑁]{Φ} dengan [𝑁] = [𝑁𝑖 𝑁𝑗 ] merupakan vector baris fungsi bentuk dan

{Φ} =

Φi {Φ } merupakan vector kolom yang memuat nilai-nilai nodal elemen. 𝑗

Fungsi bentuk mempunyai karakteristik sebagai berikut : 1. Fungsi bentuk bernilai 1 (Φ = 1) pada nodalnya dan bernilai nol (0) di nodal yang lain. 2. Jumlah 2 fungsi bentuk = 1, untuk kasus elemen linear 1 dimensi. 3. Fungsi bentuk merupakan polynomial dengan bentuk yang sama dengan persamaan interpolasi awal. 4. Turunan fungsi bentuk terhadap x = 0 untuk elemen linear 1 dimensi Berikut ini gambar fungsi bentuk linear 𝑁𝑖 dan 𝑁𝑗 : 𝑁𝑗 (𝑥)

𝑁𝑖 (𝑥)

X 𝑥𝑖

𝑥𝑗

X 𝑥𝑖

Gambar 2.3 Fungsi Bentuk Linear

𝑥𝑗

14

Contoh Ilustrasi : Elemen 1 dimensi digunakan untuk mendekati distribusi temperatur pada sirip. T pada nodal i dan j adalah 120oC dan 90oC. Tentukan T pada titik yang berjarak 4 cm dari titik asal dan gradient T dalam elemen tersebut. Koordinat nodal i dan j masing-masing adalah 1,5 dan 6 cm dari titik asal Penyelesaian :

𝜙𝑖 = 120𝑜 𝐶 𝜙𝑗 = 90𝑜 𝐶

i

j

1,5

6

Gambar 2.4. Elemen Satu Dimensi untuk Pendekatan Distribusi Temperatur Temperatur 𝜙 dalam elemen ditentukan dengna persamaan (17) : 𝑋𝑗 − 𝑥 x − Xi 𝜙= ( ) Φi + ( ) Φj 𝐿 L Data elemen : Xi = 1,5 cm

Xj = 6,0 cm

Φi = 120oC

Φj = 90oC

x = 4,0 cm

L = 4,5 cm

Diperoleh : 6−4

4−1,5

𝜙 = ( 4,5 ) 120 + ( = 103,3oC

4,5

) 90

15

Gradien temperature adalah turunan Φ terhadap x Φ𝑗 − Φ𝑖 𝑑Φ = 𝑑𝑥 𝐿 Diperoleh

𝑑Φ 𝑑𝑥

=(

90−120 4,5

) = −𝟔, 𝟔𝟕𝒐 𝑪 /𝒄𝒎

Persamaan kontinu sebagian untuk 1 dimensi disusun dengan menghubungkan beberapa persamaan linear. Persamaan linear tersebut dapat ditulis sebagai berikut:

𝜙 𝑒) = 𝑁𝑖 (𝑒) Φ𝑖 + 𝑁𝑗 (𝑒) Φ𝑗

(18)

Dengan : 𝑁𝑖 (𝑒) =

𝑋𝑗 − 𝑥 𝑋𝑗 −𝑋𝑖

𝑑𝑎𝑛 𝑁𝑗 (𝑒) =

𝑥− 𝑋𝑖

(19)

𝑋𝑗 −𝑋𝑖

Indeks (e) menunjukkan elemen. Nilai i,j dan e ditentukan dari grid elemen hingga. Misalkan batang termal seperti pada Gambar 2.4 Persamaan untuk tiap elemen : 𝜙 (1) = 𝑁1(1) Φ1 + 𝑁1(1) Φ2 𝜙 (2) = 𝑁2(2) Φ2 + 𝑁3(2) Φ3 𝜙 (3) = 𝑁3(3) Φ3 + 𝑁4(3) Φ4 𝜙 (4) = 𝑁4(4) Φ4 + 𝑁5(4) Φ5 temperatur

(1)

(2)

(3)

(4)

Gambar 2.5 Batang Termal dengan Beberapa Elemen

16

e

i

j

(1)

1

2

(2)

2

3

(3)

3

4

(4)

4

5

Tabel 2.1. Data Elemen Perhatikan bahwa N2(1) dan N2(2) adalah persamaan yang berbeda. 𝑁2(1) =

𝑥−𝑋1 𝑋2 − 𝑋1

𝑑𝑎𝑛 𝑁2(2) =

𝑋3 −𝑥 𝑋3 − 𝑋2

Masing-masing persamaan pada persamaan (2.7) berlaku untuk elemen yang sesuai dan tidak dapat dipakai di luar elemen yang bersangkutan. Untuk selanjutnya, jika persamaan dalam bentuk 𝜙 (𝑒) = N𝑖 Φ𝑖 + 𝑁𝑗 Φ𝑗 maka Ni dan Nj yang dimaksud adalah Ni(e) dan Nj(e) sedangkan Φ𝑖 dan Φ𝑗 menyatakan nilainilai nodal elemen (e).

2.7. Elemen Hingga

Persamaan elemen hingga diperoleh dari perumusan Galerkin. Penyelesaian integral residual barat (weighted residual integral) menghasilkan

satu

persamaan nodal yang dipakai secara berulang-ulang untuk menghasilkan system persamaan-persamaan linear. Suatu sistem persamaan linear diperoleh dari penyelesaian integral residual berat : 𝐻

𝑑2 𝜙

− ∫0 𝑊(𝑥) (𝐷 𝑑𝑥 2 + 𝑄) 𝑑𝑥 = 0

(20)

Dengan fungsi berat yang disusun menggunakan fungsi bentuk Ni dan Nj. Metode elemen hingga dengan fungsi berat Galerkin menentukan fungsi berat untuk nodal s, Ws, terdiri dari fungsi-fungsi bentuk untuk nodal s. Misalkan fungsi berat untuk nodal 3 pada grid linear, seperti pada Gambar 2.6, terdiri dari fungsi-fungsi bentuk untuk nodal 3 :

𝑁 (2) 𝑊3 (𝑥) = { 3 𝑁3(3)

𝑋2 ≤𝑥 ≤ 𝑋3 𝑋3 ≤𝑥 ≤ 𝑋4

(21)

17

Atau secara umum untuk fungsi berat Ws :

𝑁 (𝑒) 𝑊𝑠 (𝑥) = { 𝑠 𝑁𝑠(𝑒+1)

𝑁𝑠(𝑒) =

𝑥− 𝑋𝑟 𝑋𝑠− 𝑋𝑟

𝑋𝑟 ≤𝑥 ≤ 𝑋𝑠 𝑋𝑠 ≤𝑥 ≤ 𝑋𝑡

dan

𝑁𝑠(𝑒+1) =

(22) 𝑋𝑡−𝑥 𝑋𝑡− 𝑋𝑠

Gambar 2.6.. Fungsi Berat untuk Nodal 3

18

Gambar 2.7. Fungsi-fungsi Berat untuk (a) Nodal Pertama ,(b) Nodal Bagian Dalam (c) Nodal Terakhir dalam Grid 1 Dimensi

Fungsi berat untuk nodal pertama : W1(x) = N1(1) dan untuk nodal terakhir : Wp(x) = Np(p-1)

(23)

Selanjutnya selesaikan integral residual berat dengan menggunakan urutan nodal r,s dan t. Persamaan (3.1) menjadi : Rs = Rs(e) + Rs(e+1) 𝑑2 𝜙

𝑋

𝑅𝑠 = − ∫𝑋 𝑠 [𝑁𝑠 (𝐷 𝑑𝑥 2 + 𝑄)] 𝑟

(𝑒)

𝑋

𝑑2 𝜙

𝑑𝑥 − ∫𝑋 𝑡 [𝑁𝑠 (𝐷 𝑑𝑥 2 + 𝑄)]

(𝑒+1)

𝑠

𝑑𝑥 = 0 (24)

Karena fungsi berat Ws = 0 untuk x < Xr dan x > Xt maka Ws (x) terdiri dari 2 persamaan terpisah dalam interval Xr ≤ x ≤ Xt. Rs(e) dan Rs(e+1) adalah kontribusi elemen (e) dan (e+1) kepada persamaan residual Rs pada nodal s.

Perhatikan persamaan integral (24) dan persamaan turunan sebagai berikut : 𝑑

𝑑𝜙

𝑑2𝜙

𝑑𝑁 𝑑𝜙

𝑑𝜙

𝑑𝑁𝑠 𝑑𝜙

= 𝑁𝑠 2 + 2𝑠 (𝑁 ) 𝑠 𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑁𝑠

𝑑2𝜙 𝑑𝑥 2

=

𝑑

(𝑁𝑠 𝑑𝑥 ) − 𝑑𝑥

𝑑𝑥 𝑑𝑥

(25) (26)

19

Substitusi ke persamaan (24) diperoleh : (𝑒+1) 𝑋𝑠 𝑑2𝜙 − ∫𝑋 (𝑁𝑠 𝐷 2 ) 𝑑𝑥 𝑑𝑥 𝑟 𝑥

𝑠 ∫𝑥 (𝐷 𝑟

𝑑𝑁𝑠 𝑑𝜙 (𝑒) 𝑑𝑥 𝑑𝑥

)

= − (𝐷𝑁𝑠

𝑑𝜙 (𝑒) 𝑥𝑠 ) |𝑥 𝑑𝑥 𝑟

+

𝑑𝑥 (27)

Untuk elemen (e) sedangkan untuk elemen (e+1) : (𝑒+1) 𝑋𝑠 𝑑2𝜙 − ∫𝑋 (𝑁𝑠 𝐷 2 ) 𝑑𝑥 𝑑𝑥 𝑟 𝑥

𝑡 ∫𝑥 (𝐷 𝑠

𝑑𝑁𝑠 𝑑𝜙 (𝑒+1) 𝑑𝑥 𝑑𝑥

)

= − (𝐷𝑁𝑠

𝑑𝜙 (𝑒) 𝑥𝑡 ) |𝑥 𝑑𝑥 𝑠

+

𝑑𝑥 (28)

Telah diketahui sebelumnya bahwa 𝑁𝑠 𝑑𝜙 (𝑒)

(𝐷𝑁𝑠 𝑑𝑥 )

(𝑒)

=

𝑥− 𝑋𝑟 𝑋𝑠 − 𝑋𝑟

, dan

𝑥𝑠 𝑋 − 𝑋 𝑑𝜙 𝑋 − 𝑋 𝑑𝜙 𝑑𝜙 |𝑥 = (𝐷 𝑠 𝑟 ) − (𝐷 𝑟 𝑟 ) = 𝐷 |𝑥 = 𝑋 𝑋 − 𝑋 𝑑𝑥 𝑋 − 𝑋 𝑑𝑥 𝑑𝑥 𝑟

𝑠

𝑟

𝑠

𝑠

𝑟

Persamaan residual menjadi : 𝐻

𝑅𝑠 = 𝑅𝑠 (𝑒) + 𝑅𝑠 (𝑒+1) = − ∫0 𝑊𝑠 (𝐷

𝑑2𝜙 𝑑𝑥 2

+ 𝑄) 𝑑𝑥

𝑅𝑠 = 𝑅𝑠 (𝑒) + 𝑅𝑠 (𝑒+1) = − (𝐷 (𝑒)

𝑁𝑠 𝑄)

𝑑𝜙 (𝑒)

) 𝑑𝑥

𝑑𝑥 + (𝐷

𝑋

𝑑𝑁𝑠 𝑑𝜙

𝑟

𝑑𝑥 𝑑𝑥

|𝑥= 𝑋 + ∫𝑋 𝑠 (𝐷 𝑠 𝑑𝜙 (𝑒+1)

) 𝑑𝑥

|𝑥=𝑋 + 𝑠

− 𝑋

𝑑𝑁𝑠 𝑑𝜙

𝑠

𝑑𝑥 𝑑𝑥

𝑡 ∫𝑋 (𝐷

−

(𝑒+1)

𝑁𝑠 𝑄)

(29) Penyelesaian persamaan integral dalam persamaan (29) : Dimulai dari elemen (e) .

𝜙 (𝑒) = 𝑁𝑟 Φ𝑟 + 𝑁𝑠 Φ𝑠 𝜙 (𝑒) = dengan :

𝑋𝑠−𝑥 𝐿

Φ𝑟 +

𝑥− 𝑋𝑟 𝐿

Φ𝑠

(30)

20

𝑁𝑠 (𝑒) =

𝑥− 𝑋𝑟 𝐿

𝑑𝑁𝑠

,

𝑑𝑥

=

1

(31)

𝐿

dan 𝑑𝜙 (𝑒) 𝑑𝑥

1

=

𝐿

(−Φ𝑟 + Φ𝑠 )

(32)

Substitusi dan penyelesaian integral memberikan : 𝑋

𝑑𝑁𝑠 𝑑𝜙

𝑟

𝑑𝑥 𝑑𝑥

𝑠 ∫𝑋 𝐷

𝐷

𝑑𝑥 =

𝐿

(−Φ𝑟 + Φ𝑠 )

(33) dan 𝑋

𝑄𝐿

𝑟

2

𝑠 ∫𝑋 𝑄𝑁𝑠 𝑑𝑥 =

(34)

Maka untuk elemen (e) diperoleh :

𝑅𝑠 (𝑒) = − (𝐷

𝑑𝜙 (𝑒)

) 𝑑𝑥

|𝑥=𝑋 + 𝑠

𝐷 𝐿

(−Φ𝑟 + Φ𝑠 ) −

𝑄𝐿 2

(35) Untuk elemen (e+1)

𝜙 (𝑒+1) = 𝑁𝑠 Φ𝑠 + 𝑁𝑡 Φ𝑡 𝜙 (𝑒+1) = ( 𝑁𝑠

dengan :

𝑋𝑡−𝑥 𝐿

(𝑒+1)

𝑥− 𝑋𝑠

) Φ𝑠 + 𝑋𝑡 −𝑥

= (

𝐿

𝐿

Φ𝑡

(36)

𝑑𝑁𝑠 (𝑒+1)

);

𝑑𝑥

= −

1 𝐿

(37) dan 𝑑𝜙 (𝑒+1) 𝑑𝑥

=

1 𝐿

(−Φ𝑠 + Φ𝑡 )

(38)

Penyelesaian integral menghasilkan : 𝑋

𝑑𝑁𝑠 𝑑𝜙

𝑠

𝑑𝑥 𝑑𝑥

𝑡 ∫𝑋 𝐷

𝑑𝑥 =

𝑋

𝑄𝐿

𝑠

2

𝑡 ∫𝑋 𝑄𝑁𝑠 𝑑𝑥 =

𝐷 𝐿

(Φ𝑠 + Φ𝑡 )

(39) (40)

Kontribusi elemen (e+1) terhadap persamaan residual :

𝑅𝑠 (𝑒+1) = 𝐷

𝑑𝜙 𝑑𝑥

|𝑥=𝑋 + 𝑠

𝐷 𝐿

(Φ𝑠 + Φ𝑡 ) −

𝑄𝐿 2

(41)

21

Persamaan residual untuk nodal s : 𝑑𝜙 (𝑒+1)

𝑅𝑠 = (𝐷

) 𝑑𝑥

D (e+1)

(L)

𝑑𝜙 (𝑒)

|𝑥=𝑋 − (𝐷 ) 𝑑𝑥 𝑠

D (e+1)

] Φs − ( L )

𝐷 (𝑒)

|𝑥=𝑋 − ( ) 𝐿 𝑠

𝑄𝐿 (𝑒)

Φt − ( ) 2

𝑄𝐿 (𝑒+1)

− ( ) 2

D (e)

Φr [( ) L

+

=0 (42)

D dan Q adalah konstanta yang sama seperti ditentukan pada persamaan :

𝐷

𝑑2𝜙 𝑑𝑥 2

+ 𝑄=0

Suku ERROR pada persamaan (42) : 𝑑𝜙(𝑒+1)

𝐷 𝑑𝑥

𝑑𝜙(𝑒)

|𝑥=𝑋 − 𝐷 𝑑𝑥 𝑠

|𝑥=𝑋

(43)

𝑠

Adanya suku ini menunjukkan bahwa metode elemen hingga merupakan pendekatan. Jika suku error dihilangkan, maka persamaan residual untuk nodal s adalah : 𝐷 (𝑠−1)

𝑅𝑠 = − ( 𝐿 ) 𝑄𝐿 (𝑠−1)

(2)

𝐷 (𝑠−1)

Φs−1 + [( 𝐿 )

𝑄𝐿 (𝑠)

−( 2 )

𝐷 (𝑠)

𝐷 (𝑠)

+ ( 𝐿 ) ] Φs − ( 𝐿 )

Φ𝑠+1 −

=0 (44)

Contoh penerapan persamaan (44) pada analisis batang tumpuan sederhana dengan momen terkonsentrasi pada ujung-ujungnya. Persamaan differensial 𝑑2 𝜙

pengatur untuk semua defleksi pada batang adalah : 𝐸𝐼 𝑑𝑥 2 − 𝑀(𝑥) = 0 E

D

Q

L

1

2,4 x 1010

- 106

200

2

4,0 x 1010

- 106

200

3

4,0 x 1010

- 106

200

4

2,4 x 1010

- 106

200

Tabel 2.2. Data Elemen

22

Bentuk persamaan (44) dengan Q dan L konstan adalah :

𝑅𝑠 =

−𝐷 (𝑠−1) 𝑌𝑠−1 +(𝐷 (𝑠−1) + 𝐷 (𝑠) )𝑌𝑠 − 𝐷(𝑠) 𝑌𝑠+1 𝐿

− 𝑄𝐿 = 0

Y = nilai defleksi nodal (𝜙)

Gambar 2.8. Batang Tumpuan Sederhana Persamaan residual untuk nodal 2,3 dan 4 adalah : R2 = - 1,2 Y1 + 3,2 Y2 – 2,0 Y3 + 2 = 0 R3 = - 2,0 Y2 + 4,0 Y3 – 2,0 Y4 + 2 = 0 R4 = - 2,0 Y3 + 3,2 Y4 – 1,2 Y5 + 2 = 0 (untuk 3 persamaan ini 108 ‘dihilangkan’) Tumpuan pada kedua ujung batang menunjukkan Y (0) = Y (800 cm) = 0 sehingga kondisi batas Y1 = 0 dan Y5 = 0, selanjutnya diperoleh set persamaan : R2 = 3,2 Y2 – 2,0 Y3

=-2

R3 = -2,0 Y2 + 4,0 Y3 - 2,0 Y4 = - 2 R4 =

- 2,0 Y3 + 3,2 Y4 = - 2

23

Diselesaikan dan diperoleh : Y2 = -2,50 cm

Y3 = -3,0 cm

Y4 = -2,5 cm

a). Perhitungan defleksi di x = 300 cm, berada pada elemen (2) Y(2) = N2(2)Y2 + N3(2)Y3 =(

𝑋3 − 𝑥 𝑋3 − 𝑋2

𝑥− 𝑋2

) 𝑌2 + (𝑋

3 − 𝑋2

Diketahui X2 = 200 cm ;

) 𝑌3

X3 = 400 cm

Maka nilai simpangan di x = 300 cm : =(

400 – 300

300−200

) (−2,5) + (300−200) (−3,0) 400−200

= - ½ (2,5 + 3,0) = - 2,75 cm b). Perhitungan slope di elemen (1) : 𝑑𝑦 (1) 𝑑𝑥

=

1 𝐿

(−𝑌1 + 𝑌2 ) =

−2,5−0 200

= - 0,0125 cm/ cm

Sistem persamaan-persamaan linear pada contoh di atas dapat dinyatakan dalam notasi matrik : 𝑅2 𝑌2 3,2 −2 0 −2 0 {𝑅3 } = [−2 4 −2] {𝑌3 } − {−2} = {0} 𝑅4 0 −2 3,2 𝑌4 −2 0 atau dalam bentuk persamaan matrik {R} = [K] {Y} – {F} = {0} dengan [K] menyatakan matrik system, {Y} menyatakan vektor simpangan, {F} menyatakan vektor gaya luar dan {R} menyatakan vektor residu untuk tiap elemen.

2.8. Formula Weak

Sebelum menerapkan Metode Elemen Hingga untuk memecahkan persamaan dengan kondisi batas, perlu untuk mengubah persamaan menjadi bentuk yang lebih cocok. Untuk melakukan itu ada dua alternatif: 1. Turunan satu dapat memperoleh masalah minimalisasi setara, yang memiliki tepat solusi sama dengan persamaan diferensial. 2. Turunan satu dapat memperoleh apa yang disebut formulasi lemah. Kedua metode akhirnya mengarah kepada hasil yang sama persis, namun, karena untuk persamaan umum untuk diperlakukan tidak ada masalah

24

minimisasi yang setara , maka akan membatasi diri untuk metode kedua. Awalnya formulasi weak atau lemah telah diperkenalkan oleh matematika murni untuk menyelidiki perilaku solusi dari persamaan diferensial parsial, dan untuk membuktikan keberadaan dan keunikan dari solusi. Kemudian skema numerik telah didasarkan pada formulasi ini yang menyebabkan solusi perkiraan dengan cara yang konstruktif.

Dapat dilihat bahwa kondisi batas penting secara otomatis menunjukkan bahwa fungsi tes yang sesuai adalah sama dengan nol, sedangkan kondisi batas natural tidak memaksakan kondisi apapun baik dengan tidak diketahui atau fungsi tes. Hal ini tidak segera jelas apakah kondisi batas penting atau alami, kecuali dalam kasus di mana terdapat masalah minimisasi sesuai. Secara umum, bagaimanapun, dapat dikatakan bahwa untuk persamaan diferensial orde kedua, semua kondisi batas yang mengandung turunan pertama yang alami, dan fungsi yang diberikan pada batas sangat penting.

Dalam masalah rangka keempat situasinya lebih kompleks. Namun, untuk masalah fisik, secara umum, dapat dinyatakan bahwa jika kondisi batas mengandung turunan kedua atau ketiga mereka natural, sedangkan kondisi batas yang hanya berisi fungsi atau urutan pertama turunan sangat penting. Cara termudah untuk memeriksa apakah kondisi batas penting atau natural adalah untuk mempertimbangkan integral batas. Jika dalam beberapa cara syarat batas dapat diganti, kondisi batas wajar. Jika tidak kondisi ini penting dan fungsi pengujian harus dipilih sedemikian rupa sehingga integral batas lenyap.

2.9. Metode Galerkin

Titik awal adalah yang disebut sebagai metode Galerkin. Dalam metode ini solusinya c didekati oleh kombinasi linear dari fungsi ekspansi yang disebut fungsi dasar: 𝑛

𝑐 𝑛 (𝑥) = ∑ 𝑐𝑗 𝜑𝑗 (𝑥) + 𝑐0 (𝑥) 𝑗=1

(45)

25

di mana parameter 𝑐𝑗 harus ditentukan. Fungsi dasar 𝜑𝑗 (𝑥) harus independen linear. Selain itu harus sedemikian rupa sehingga fungsi sewenang-wenang dalam ruang solusi dapat didekati dengan akurasi yang sewenang-wenang, tersedia dalam jumlah yang memadai fungsi dasar yang digunakan dalam kombinasi linear (45). Fungsi 𝑐0 (𝑥) harus dipilih sedemikian rupa sehingga 𝑐 𝑛 (𝑥) memenuhi kondisi batas penting. Secara umum ini berarti bahwa 𝜑𝑗 (𝑥) = 0 dan 𝑐0 (𝑥) = 𝑔.Dalam rangka untuk menentukan parameter 𝑐𝑗 (𝑗 = 1,2, … , 𝑛) fungsi tes 𝑣 dipilih dalam ruang yang direntang oleh fungsi dasar 𝜑1 (𝑥)untuk 𝜑𝑛 (𝑥). Hal ini cukup subtitusi 𝑣(𝑥) = 𝜑1 (𝑥)

(46)

ke persamaan yang akan dicari yang telah ditentukan formula weaknya terlebih dahulu. Setelah itu subtitusi persamaan (45) dan (46) ke formula weak yang diperoleh, disebut sebagai formula Galerkin.