Metnum03-Iza Azmar Aminudin

                                                               بِسْمِ اللّهِ الرَّحْمَنِ الرَّحِيْمِ

Assalamu'alaikum Wr. Wb,

Salam Sejahtera untuk kita semua,

Biodata Diri

Perkenalkan saya:

Nama  : Iza Azamr Aminudin

NPM  : 1806233316

Kelas  : Metnum-03

Melalui laman web ini, insya Allah saya akan berbagi kepada teman-teman sekalian terkait pembelajaran mata kuliah Metode Numerik yang saya pelajari bersama pa DAI dan teman-teman di kelas Metnum 03. Semoga yang saya tulis dalam laman ini dapat memberi hikmah bagi saya dan teman-teman pembaca, aamiiin.

Pendahuluan | Filosofi dalam Belajar Metode Numerik

Salah satu kewajiban hakiki untuk manusia adalah untuk belajar memetik hikmah dari semua pembelajaran yang akan terus ada selama hidup ini. Selain itu, Nabi Muhammad SAW sebagai idola saya pernah berkata bahwa sebaik-baiknya manusia adalah manusia yang paling bermanfaat bagi manusia (H.R Ahmad). Kedua hal ini insya Allah terus memotivasi saya untuk terus belajar dan memanfaatkan ilmu yang saya pelajari baik bagi kehidupan saya sendiri hingga untuk kepentingan orang banyak lainnya. Melalui mata kuliah metode numerik ini saya berharap untuk mendapatkan ilmu baik yang berupa materi pembelajaran seperti yang diajarkan oleh Pa DAI, dari buku, maupun melalui pengajaran yang dilakukan oleh asisten dosen dan teman-teman saya. Selain itu, saya berharap dengan menekuni pembelajaran metode numerik ini dapat menambah pengetahuan dan "ketidaktahuan" tentang hal-hal yang belum saya pelajari maupun mengerti sehingga terus menambah motivasi saya dalam belajar dan juga menjaga saya dari sikap riya dan takabur yang dapat membuat saya menjadi orang yang berpikiran tertutup dan menututp keran ilmu pengetahuan,

Semangat teman-teman semua dalam menimba ilmu dalam mata kuliah ini.

Dari ‘Abdullah bin ‘Amr dan Anas bin Malik radhiyallahu ‘anhuma, Rasulullah shallallahu ‘alaihi wa sallam bersabda,

قيِّدُوا العِلمَ بالكِتابِ

“Jagalah ilmu dengan menulis.” (Shahih Al-Jami’, no.4434. Syaikh Al-Albani mengatakan bahwa hadits ini sahih).

Pertemuan 1: Review Materi Pra-UTS dan Pengantar Software Open-Modelica | 9 November 2020

Review Materi Sebelum UTS

Pembelajaran mata kuliah Metode Numerik-03 sebelumnya diampu oleh Bapak Dr. Ir. Engkos A. Kosasih, M.T. Selama pembelajaran sebelum UTS, kami mempelajari beberapa konsep yang digunakan dalam menyelesaikan permasalahan dengan metode numerik. Diantara materi-materi tersebut diantanya adalah sebagai berikut.

1. Deret Maclaurin dan Deret Taylor

Teorema Ketunggalan Andaikan f memenuhi uraian berikut,

                                                          f(x) = c0 + c1(x - a) + c2(x - a)^2  + c3(x - a)^3 .  .  .

untuk semua x dalam selang di sekitar a, maka cn = Jadi suatu fungsi tidak dapat digambarkan oleh dua deret pangkat dari (x - a).

2. Gauss Elimination

Tutorial Open-Modelica

https://www.youtube.com/watch?v=l7I8YlQBdug

Pertemuan 2: Review Tugas Modelica dan Simulasi Penyelesaian Metode Numerik | 16 November 2020

Pada kesempatan kali ini, kami mempelajari cara membuat program sederhana didalam modelica untuk menjumlahkan angka 10 kepada sebuah variable dan mencari rata-rata nilai dengan menggunakan software modelica.

1. Membuat program sederhana didalam modelica untuk menjumlahkan angka 10 kepada sebuah variable software modelica.

Memasukan Coding untuk permasalahan

Menguji logic coding

Melakukan Simulasi

2. Mencari rata-rata nilai dengan menggunakan software modelica.

Memasukan Coding untuk permasalahan

Menguji logic coding

Melakukan Simulasi

3. Tugas 02

Melakukan Eliminasi Gauss untuk menyelesaikan persamaan linear tiga variabel. Tugas kali ini yaitu Pr2 metnum ... Buat class dengan type function untuk menyelesaikan persamaa2n aljabar simultan (metoda gauss elim, gauss seidel ataupun metoda lain) ... Dan sebuah class untuk menjalankan fungsi tersebut. Waktu 1 minggu dan dikumpulkan di wikipage air.eng.

Berikut adalah persamaan yang akan saya uji

Persamaan Aljabar Linear

Kemudian saya membuat class sebagai berikut.

Class Gauss Elimination

Dan membuat function sebagai berikut

Function Gauss Elimination

Setelah function dan class dibuat, barulah saya melakukan simulasi untuk mendapatkan hasilnya. Berikut adalah hasil yang saya dapatkan.

Grafik Hasil Function

Pertemuan 3: Menyelesaikan Spring-mass System| 23 November 2020

Menyelesaikan persamaan Spring-mass System dengan menggunakan software modelica

Soal Spring-mass System

Penjabaran Soal Spring-mass System

Persamaan Dari Problem di Soal

Persamaan Hk. Hooke dari Problem di Soal

Matriks Persamaan Hk. Hooke

Nilai Delta x dapat dicari dari matriks tersebut dengan eliminasi Gauss menggunakan software modelica.

Coding di Class Modelica

Coding di Function Modelica

Sehingga didapatlah hasil sebagai berikut.

Hasil nilai delta x dan grafiknya dalam modelica

Tugas 3 Dalam tugas 3 kami diberikan suatu permasalahan metode numerik dalam suatu truses

Soal

Memisalkan trus dalam node-node

Pemisalan truss menjadi node node

Berikut adalah persamaan

penyelesaian dengan rumus

Persamaan solusi dari soal

Matrix penyelesaian dari soal

Matrix solusi dari trusses untuk diselesaikan dengan openmodelica

Coding di Class Modelica

Coding di Class Modelica

Coding di Class Modelica

hasil grafik dan perhitungan dari modelica

Coding di Class Modelica

Pertemuan 4: Membahas Penyelesaian Trusses dan Quis 1| 30 November 2020

Soal Quis 1

Soal Quiz no. 4

Soal Quiz 8

Langkah-langkah dalam menyelesaikan soal

Flow chart

Menyelesaikan Soal no. 4

Proses mencari persamaan matriks untuk menentukan nilai konstanta kekakuan global

Solusi

Solusi

Setelah mendapatkan persamaan dari matriks untuk mencari kekakuan global dari masing-masing titik, langkah selanjutnya adalah dengan menggunakan perangkat lunak Open Modelica untuk mencari kekakuan global dari masing-masing titik tersebut

Berikut adalah Tampilan Coding di Software Modelica Beserta rincian Codingnya

Coding Untuk Matriks Mencari nilai Konstanta Kekakuan Global

Coding

Tampilan Modelica

Coding Untuk Matriks Mencari nilai Konstanta Kekakuan Global

Selanjtnya dilakukan pengecekan Coding dan dilakukanlah simulasi

Solusi

Berikut adalah hasil yang didapatkan dari simulasi menggunakan Open Modelica untuk mencari Kekakuan Lokal

Solusi

Setelah itu, langkah selanjutnya adalah mencari displacement (U) dan reaction (R). Untuk mencari displacement, perlu diperhatikan boundary dan external load yang ada dan dimasukan kedalam matriks. Karena titik 1 dan 3 merupakan pinned, maka tidak ada displacement pada titik 1 dan 3(U1x=U1Y=U3X=U3Y=0). External load berada pada titik 2 dan 4 (pada arah x sebesar -4000Cos(15) sedangkan pada arah y sebesar 4000Cos(15)). Pada Open Modelica, menggunakan persamaan A*U=F dan R=(KG*U)-R untuk mencari displacement dan reaction setiap titik:

Berikut adalah coding yang digunakan untuk mencari displacement dan reaction. Dalam simulasinya, harus menggunakan NaiveGauss agar simulasi dapat dilakukan

Tampilan Modelica Untuk Coding Model U dan R serta function NaiveGauss

Coding untuk mencari nilai U dan R

Coding function NaiveGauss untuk mencari nilai U dan R

Berikut adalah hasil yang didapatkan

Solusi untuk U dan R

Pertemuan 5: Menyelesaikan Persamaan Truss dengan Open Modelica untuk 2 dimensi dan 3 dimensi| 7 Desember 2020

Terdapat beberapa cara yang dapat dilakukan untuk menyelesaikan suatu persoalan Truss dengan menggunakan Open Modelica. Pada pertemuan ini kami membahas tentang soal quiz yang telah diberikan pada pertemuan pekan 4. Jadi fungsi mempelajari soal ini adalah agar kita dapat lebih memahami pengaplikasian dari software openmodelica pada contoh kasus yang real. Pertemuan ini menjelaskan bagaimana alur pengerjaan soal quiz yang diberikan agar dapat membuat semua persamaan yang dapat digunakan pada kasus yang lebih umum.

Pada pertemuan ini, Pak Dai menunjuk ahmad muhammad fahmi untuk menjelaskan kepada teman-teman di kelas tentang fungsi-fungsi yang saya gunakan untuk mengerjakan soal quiz. Beberapa hal yang saya jelaskan antara lain:

1. Cara melakukan looping menggunakan fungsi for.

2. Cara membuat data dalam bentuk array.

3. Cara menggunakan fungsi if.

4. Apa itu floating number dan bagaimana cara menghilangkannya.

5. Apa itu protected variable.

6. Fungsi dari modelica class dengan specialization class dan function.

Tugas 5

soal trusses 3D

Pertemuan 6: Muhasabah diri | 14 Desember 2020

Pada pertemua ini, kami diminta untuk menilai pemahaman diri kami sendiri untuk mata kuliah metode numerik. Kami diminta untuk mempresentasikan seberapa dalam pemahaman kami terkait materi metode numerik. Kami diminta untuk menjelaskan mengenai konsep fisika dari soal trusses yang selama ini kami garap. Lalu setelah itu, kami menjelaskan bagaimana cara untuk mensimulasikan soal trusses tersebut dalam perangkat openmodelica. Dalam soal kami mencari beberapa parameter seperti displacement lokal maupun global, konstanta stiffness lokal maupun global, serta reaction force yang dialami oleh tumpuan truss.

Dari kelas hari ini, saya mendapatkan beberapa bahan evaluasi yang perlu untuk saya perbaiki kedepannya. Pertama, jika dibandingkan dengan rekan-rekan sayan seperti ahmad muhammad fahmi, saya masih harus lebih banyak belajar lagi mengenai penggunaan openmodelica untuk mensimulasikan soal trusses ini. Kedepannya saya akan banyak berlatih dan mengeksplor lagi penggunaan openmodelica.

Pertemuan 7: aplikasi metode numerik dalam kasus optimisasi | 21 Desember 2020

Pada pertemuan kali ini kami mempelajari simulasi optimalisasi One-Dimensional Unconstrained menggunakan software openmodelica bersama dengan pa Dai dan ibu Candra. Kami mempelajari bagaimana cara mengaplikasikan perangkat openmodelica untuk mengoptimalkan suatu sistem trusses. Pada pertemuan ini kami mempelajari metode golden ratio untuk melakukan simulasi dalam suatu sistem.

Sinopsis Tugas Besar: Mengoptimalkan Sistem Trusses untuk Mencapai Best Efficiency Point (BEP)

Tugas Besar

Pada tugas besar kali ini, kami diberi tugas untuk mengaplikasikan ilmu metode numerik untuk mengoptimasi suatu sistem trusses dengam menggunakan perangkat lunak openmodelica. Berikut adalah geometri rangka beserta parameter-parameternya.

Tugas Besar Metnum

Tugas kami adalah mencari Best Efficiency point untuk sistem trusses ini mengacu pada harga dan juga kemampuan fisik dari trusses itu sendiri. Dalam tugas besar kali ini, kami mengerjakan pembuatan coding tersebut secara bersama-sama dari 2 kelas metode numerik 02 dan 03 dengan dipimpin oleh Ahmad Muhammad fahmi, Josiah Enrico, dan juga Christoper.

Flow Chart

Berikut adalah flow chart pengerjaan tugas besar dan pembuatan codinng untuk tugas besar

Flowchart Tugas Besar Metnum

Flowchart membuat coding tugas besar

Pengerjaan Tugas Besar: Optimasi Sistem Trusses untuk Mencapai Best Efficiency Point (BEP)

1. Memahami karakteristik fisik trusses

Kasus kali ini, kami diberikan suatu bentuk trusses dengan 16 komponen penyusunnya (truss) dan terdapat 24 point dalam trusses tersebut. Trusses tersebut tersebut disusun secara hotizontal dan vertikal sehingga sudut yang dihasilkan dalam trusses hanya berbentuk 0 dan 90 derajat menurut pada sumbu x negatif. Kami diminta untuk membuat optimasi dari trusses ini berdasarkan pada harga serta kemampuan fisiknya. Governing Equation yang digunakan dalam sistem ini adalah dari formulasi finite element untuk defleksi sebagai berikut.

                                                                     F=(AE/L)∆L

Dengan F adalah gaya yang dialami oleh trusses, A adalah luas penampang trusses, E adalah elastisitas bahan, L adalah panjang trusses dan ∆L adalah pertambahan panjang.

Selain itu, untuk mengevaluasi perhitungan , kami menggunakan hukum hooke untuk mengetahui hubungan antara stiffness trusses dengan gaya-gaya yang terjadi dalam trusses. Persamaan dasarnya adalah sebagai berikut.

                                                                       k=AE/L

Dengan k adalah konstanta kekakuan.

Nantinya, trusses ini akan diselesaikan dengan metode numerik dengan bentuk matriks yang memiliki persamaan dasar untuk dislacement sebagai berikut sebagai berikut.

                                                                     {U}=[T]{u}

dengan U adalah displacement dari perspektif global, dan u adalah displacement untuk perspektif lokal serta T adalah matriks transformasi.

Sementara hubungan antara nilai force yang dialam trusses antara lokal dan global dihubungkan dengan persamaan berikut.

                                                                     {F}=[T]{f}

dengan F adalah force yang terjadi pada node dari perspektif global, dan f adalah force yang dialami node untuk perspektif lokal serta T adalah matriks transformasi.

Persamaan tersebut dielaborasikan menjadi persamaan sebagai berikut

                                                                   {F}=[T][K][T]⁻¹{U}

Tugas Besar Metnum

2. Mendefinisikan Asumsi dalam Perhitungan

Asumsi yang digunakan untuk melakukan perhitungan ini antara lain:

• Beban akan terdistribusi hanya pada node (karena bersifat trusses).
• Safety factor minimal bernilai 2.
• Batas displacement 0,001 m sebelum buckling (pada truss paling atas).
• Ketinggian trusses pada tiap lantai sama yaitu 0,6 m.
• tidak ada bending karena bersifat truss
• Beban terdistribusi pada node

Metodologi

A. Untuk jenis material yang sama:

1. Mencari harga untuk 6 ukuran batang dengan material yang sama yaitu ASTM A36. 2. Mengitung nilai safety factor pada 6 ukuran batang dengan coding awal. 3. Membuat rasio antara safety factor dengan harga total. 4. Membuat persamaan antara rasio (safe/harga) dengan area dengan metode curve fitting. 5. Melakukan optimasi menggunakan metode golden section.

B. Untuk area penampang yang sama:

1. Mencari harga untuk 3 jenis material dengan area penampang yang sama yaitu 171 mm^2. 2. Mengitung nilai safety factor pada 3 jenis batang dengan coding awal. 3. Membuat rasio antara safety factor dengan harga total. 4. Membuat persamaan antara rasio (safe/harga) dengan jenis material dengan metode curve fitting. 5. Melakukan optimasi menggunakan metode golden section.

3. Mendefinisikan bagian-bagian dalam trusses

Asumsi node dan trusses

Constraint

- Node 1,2,3,4 (lantai dasar) fixed

- Beban F1 dan F2 terdistribusi ke node sekitaranya, sehingga:

1. Node 13 & 16 = 1000N

2. Node 14 & 15 = 500N