Metnum 03 - Muhammad Fairuz Daffa

Biodata

Nama: Muhammad Fairuz Daffa

NPM: 1806181716

Kelas: Metnum 03

Materi Sebelum UTS

Sebelum UTS, kelas metnum 03 diajar oleh Pak Engkos, sebelum uts saya sendiri belajar mengenai:

1. Turunan Numerik

Turunan numerik terbagi 3 yaitu:

1. Turunan Maju (Forward)

2. Turunan Mundur (Backward)

3. Turunan Pusat (Centered)

Dari ketiga metode turunan ini, metode turunan pusat atau centered memiliki kemungkinan error yang paling kecil.

2. Pencarian Akar - Akar Persamaan

Untuk mencari akar - akar persamaan dapat digunakan 2 cara cara yaitu:

1. Close method (Bracketing) : Terdiri dari graphical method, false position dan bisection

2. Open method : Terdiri dari newton raphson, secant, dan simple fix point

3. Interpolasi dan Regresi

Interpolasi adalah metode menghasilkan titik-titik data baru dalam suatu jangkauan dari suatu set diskret data-data yang diketahui.

Regresi adalah salah satu metode untuk menentukan hubungan sebab-akibat antara variabel dengan variabel lainnya.

Pertemuan tanggal 9 November 2020

Pada oertemuan ini Pak Dai menginstruksikan mahasiswa untuk mendownload software open modelica untuk pengaplikasian dari mata kuliah metode numerik, dan juga meminta kita untuk mempelajari software open modelica melalui video tutorial dari youtube, dan membuat video youtube berisikan kami mempelajari atau mengaplikasikan ilmu dari video tutorial yang sudah kita tonton.

Saya mempelajari software openmodelica dari tautan dibawah ini

 https://youtu.be/Dw66ODbMS2A

Pertemuan Kedua, Tanggal 16 November 2020

Pada pertemuan ini, pak DAI mengevaluasi tugas kita satu per satu. Lalu sambil menunggu satu per satu di interview, pak DAI memberikan arahan untuk membuat program sederhana di software openmodelica.

Pertama pak DAI mengajarkan untuk membuat program untuk meenghitung penambahan sepuluh terhadap suatu angka (x+10).

Kedua pak DAI menginstruksikan untuk membuat program yang bertujuan untuk menghitung nilai rata - rata dari beberapa kumpulan data. Berikut adalah simulasi yang saya coba

Tugas Pertemuan 2

Pada tugas ini pak DAI menginstruksikan untuk membuat program yang bisa mensolve permasalahan aljabar simultan. Pada tugas ini saya membuat program untuk menyelesaikan persamaan matriks dengan metode gauss jordan.

1. Soal yang saya gunakan adalah sebagai berikut.

2. Saya menggunakan fitur function dan class untuk menyelesaikan persamaan gauss jordan diatas.

3. Pada function saya menggunakan formula sebagai berikut:

4. Untuk class saya menggunakan formula sebagai berikut:

5.Setelah itu saya lakukan pengecekan terhadap formula tersebut, didapat hasil sebagai berikut

6. Lalu setelah pengecekan berhasil, saya melakukan simulasi, lalu didapat hasil sebagai berikut

7. Hasil yang didapat adalah x1 = 1; x2 = 4; x3 = -2

Pertemuan 3: Senin, 23 November 2020

Pada pertemuan ini kita disuruh mensolve soal dari buku sebagai berikut:

Tugas 3

---

Pada tugas ini kita disuruh mensolve soal dibuku sebagai berikut:

Pada tugas ini kami diminta untuk menghitung defleksi pada setiap batang dan gaya reaksinya.

Pertama saya mencari deltaL atau displacement dari setiap batang menggunakan open modelica. Dengan coding sebagai berikut:

Lalu didapat hasil sebagai berikut:

Kedua saya mencari nilai U dengan rumus RG x U = F, jika di translate ke coding maka akan menjadi

Lalu didapat hasil sebagai berikut

Ketiga saya mencari R menggunakan openmodelica dengan rumus R = (K x U) - F, jika di translate ke coding maka akan menjadi sebagai berikut:

Lalu hasil yang saya dapat adalah sebagai berikut:

Pertemuan 4: Senin, 30 November 2020

Pada pertemuan ini pak DAI menjelaskan perbedaan antara statika struktur dan dinamika struktur.

QUIZ

Berikut adalah flowchart saya dalam mengerjakan soal

Nomor 1

---

Class class no1Quiz // DECLARATIONS // Data for each member: [element #, node i, node j, theta, area, modulus, length] // SI Units: [integer, integer, integer, degrees, meters^2, pascals, meters] // Imperial Units: [integer, integer, integer, degrees, inches^2, lb/inches, inches] parameter Real [:,7] member = [1, 1, 2, 0, 8, 1.9e6, 36; 2, 2, 3, 135, 8, 1.9e6, 50.9; 3, 3, 4, 0, 8, 1.9e6, 36; 4, 2, 4, 90, 8, 1.9e6, 36; 5, 2, 5, 45, 8, 1.9e6, 50.9; 6, 4, 5, 0, 8, 1.9e6, 36]; // External loads for each node: [node #, FX, FY] parameter Real [:,3] node_load = [1, 0, 0; 2, 0, 0; 3, 0, 0; 4, 0, -500; 5, 0, -500]; // Vector for equivalent stiffness constant of each member Real [size(member,1)] k; // Array for Stiffness Matrix Real [size(member,1),4,4] Ke; equation // Calculate stiffness constants for each member element k = {(member[i,5] * member[i,6] / member[i,7]) for i in 1:size(member,1)}; // Find stiffness matrix for each element Ke = stiffnesselement(member); end no1Quiz;

Matriks Elemen function stiffnesselement input Real [:,7] beam_mat; output Real [size(beam_mat,1),4,4] Ke_mat; protected Real theta; Real [3] StiffTrig; Real [4,4] StiffTransformation; Real [size(beam_mat,1)] k_vec; Real float_error = 10e-10; algorithm k_vec := {(beam_mat[i,5] * beam_mat[i,6] / beam_mat[i,7]) for i in 1:size(beam_mat,1)}; // Finding stiffness matrix of each element member for i in 1:size(beam_mat,1) loop // Clearing the matrices StiffTrig := zeros(3); StiffTransformation := zeros(4,4); // Converting degrees to radians theta := Modelica.SIunits.Conversions.from_deg(beam_mat[i,4]); // {cos^2, sin^2, sincos} StiffTrig := {(Modelica.Math.cos(theta))^2, (Modelica.Math.sin(theta))^2, (Modelica.Math.sin(theta)*Modelica.Math.cos(theta))}; // Handle float error elements in StiffTrig for t in 1:size(StiffTrig,1) loop if abs(StiffTrig[t]) <= float_error then StiffTrig[t] := 0; end if; end for; // Construct stiffness transformation matrix StiffTransformation := [ StiffTrig[1], StiffTrig[3], -1*StiffTrig[1], -1*StiffTrig[3]; StiffTrig[3], StiffTrig[2], -1*StiffTrig[3], -1*StiffTrig[2]; -1*StiffTrig[1], -1*StiffTrig[3], StiffTrig[1], StiffTrig[3]; -1*StiffTrig[3], -1*StiffTrig[2], StiffTrig[3], StiffTrig[2]]; // Multiply in stiffness constant of element, add final stiffness matrix to Ke_mat for m in 1:4 loop for n in 1:4 loop Ke_mat[i,m,n] := k_vec[i] * StiffTransformation[m,n]; end for; end for; end for; end stiffnesselement;

Pertemuan 5: Senin, 7 Desember 2020

Pada pertemuan ini, pak DAI membahas tentang soal quiz yang telah diberikan. Teman saya, Ahmad Mohammad Fahmi menjelaskan tentang coding yang dia buat dari awal sampai akhir. Namun, coding dari Fahmi masih ada yang kurang, karena total gaya pada ketiga sumbu untuk soal nomor 8 masih belum sama dengan 0.

Pak DAI menginstruksikan kita untuk bertanya kepada fahmi terkait hal yang belum dimengerti. Lalu pak DAI memberikan soal lagi dan menginstruksikan untuk memakai coding dari fahmi untuk mengerjakan soal tersebut. Lalu kami semua mencoba mengerjakan. Lalu pak DAI menginstruksikan untuk memberikan feedback terhadap codingan yang dibuat fahmi, dan memperbaikinya jika bisa.

Tugas Besar

Asumsi yang digunakan untuk melakukan perhitungan ini antara lain:

Variasi Stiffness terikat dengan variabel area. Memvariasikan Elastisitas tergolong sulit karena setiap material memiliki range yang tidak teratur dan dalam satu material yang sejenis (struktur biaya tetap) tidak terjadi perubahan nilai elastisitas yang berbanding lurus dengan perubahan biaya.

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).

Metodologi

---

Untuk jenis material yang sama:

1. Mencari harga untuk 6 ukuran batang dengan material yang sama yaitu SS304.

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.

Untuk area penampang yang sama:

1. Mencari harga untuk 4 jenis material dengan area penampang yang sama yaitu 171 mm^2.

2. Mengitung nilai safety factor pada 4 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.