Panduan BEM: Perhitungan, Coding, dan Interpretasi

Halo semuanya! Kali ini kita akan membedah Metode Elemen Batas atau BEM. Kita mengacu pada paper klasik C.A. Brebbia tahun seribu sembilan ratus tujuh puluh delapan. Penerapannya adalah untuk melihat distribusi panas pada sebuah dinding yang kita bagi menjadi titik-titik elemen batas di permukaannya.

Langkah perhitungannya dimulai dari Identitas Green Kedua, yaitu c i dikali u i ditambah integral pada batas gamma dari u dikali q bintang de gamma sama dengan integral q dikali u bintang de gamma. Kemudian, kita melakukan diskritisasi elemen menjadi bentuk sigma H topi i j dikali u j sama dengan sigma G i j dikali q j.

Selanjutnya, kita membentuk matriks H dikali vektor U sama dengan matriks G dikali vektor Q. Di sini, matriks G berasal dari integrasi solusi fundamental u bintang, sedangkan matriks H dari turunan solusi q bintang. Vektor U adalah temperatur dan Q adalah fluks panas yang akan kita cari.

Untuk struktur kodingnya, kita menggunakan perulangan bertingkat. Untuk setiap node sumber i dan setiap elemen batas j, program menghitung nilai integral G i j dan H i j. Setelah semua komponen matriks terisi, kita menyusunnya ke dalam sistem persamaan linear A x sama dengan B untuk diselesaikan.

Mari kita baca output grafiknya. Sumbu x adalah posisi dan sumbu y adalah temperatur T. Kita bisa melihat nilai batas di tiap node, temperatur di titik internal dalam dinding, serta kemiringan garis yang merepresentasikan gradien suhu atau fluks panas secara detail.

Kesimpulannya, analisis BEM ini sangat efisien karena hanya memerlukan data pada batas dinding saja. Akurasinya sangat baik untuk masalah tak hingga dan sangat cocok untuk implementasi desain termal pada konstruksi bangunan.

Itulah penjelasan rinci mengenai penerapan Metode Elemen Batas mulai dari teori paper hingga hasil simulasinya. Semoga penjelasan ini membantu kalian memahami distribusi panas dengan lebih mendalam. Sampai jumpa di pembahasan berikutnya!