Difference between revisions of "Metnum03-Wildan Firdaus"

From ccitonlinewiki
Jump to: navigation, search
Line 349: Line 349:
 
[[File:K2D.png|500px|thumb|center]]
 
[[File:K2D.png|500px|thumb|center]]
  
[[File:K Lokal 2D1.png|500px|thumb|center]]
+
[[File:K Lokal 2D1.png|800px|thumb|center]]
  
[[File:K Lokal 2D2.png|500px|thumb|center]]
+
[[File:K Lokal 2D2.png|800px|thumb|center]]
  
 
Setelah didapatkan nilai K lokal, kita konversi menjadi K global pada setiap batang. Kemudian nilai K pada setiap batang ini dijumlahkan menjadi K global total. Disini saya menggunakan perhitungan manual menggunakan bantuan excel, sehingga menghasilkan sebagai berikut.
 
Setelah didapatkan nilai K lokal, kita konversi menjadi K global pada setiap batang. Kemudian nilai K pada setiap batang ini dijumlahkan menjadi K global total. Disini saya menggunakan perhitungan manual menggunakan bantuan excel, sehingga menghasilkan sebagai berikut.
Line 363: Line 363:
 
Sehingga untuk U dan R nya menjadi seperti berikut :
 
Sehingga untuk U dan R nya menjadi seperti berikut :
  
[[File:Gauss_Jordan2D.png|500px|thumb|center]]
+
[[File:Gauss_Jordan2D.png|800px|thumb|center]]
  
[[File:U2D.png|500px|thumb|center]]
+
[[File:U2D.png|800px|thumb|center]]
  
[[File:Reaksi 2D.png|500px|thumb|center]]
+
[[File:Reaksi 2D.png|800px|thumb|center]]
  
  
Line 383: Line 383:
  
  
[[File:K3D1.png|500px|thumb|center]]
+
[[File:K3D1.png|800px|thumb|center]]
  
[[File:K3D2.png|500px|thumb|center]]
+
[[File:K3D2.png|800px|thumb|center]]
  
[[File:K Lokal 3D.png|500px|thumb|center]]
+
[[File:K Lokal 3D.png|800px|thumb|center]]
  
  
Line 404: Line 404:
 
Sehingga untuk U dan R nya menjadi seperti berikut :
 
Sehingga untuk U dan R nya menjadi seperti berikut :
  
[[File:Gauss_Jordan3D.png|500px|thumb|center]]
+
[[File:Gauss_Jordan3D.png|800px|thumb|center]]
  
 
[[File:U3D.png|500px|thumb|center]]
 
[[File:U3D.png|500px|thumb|center]]
  
[[File:Reaksi 3D.png|500px|thumb|center]]
+
[[File:Reaksi 3D.png|800px|thumb|center]]

Revision as of 12:17, 7 December 2020

Assalamualaikum Warahmatullahi Wabarakatuh

BIODATA DIRI

Nama : Wildan Firdaus

NPM  : 1906435574

Fakultas/ Jurusan : Teknik/ Teknik Mesin

Kelas Metode Numerik 03

Pertemuan Metode Numerik 1 : 09 November 2020

Metode numerik merupakan teknik penyelesaian permsalahan matematis dengan menggunakan operasi hitungan yaitu operasi tambah, kurang, kali, dan bagi. Metode ini digunakan karena banyak permasalahan matematis tidak dapat diselesaikan menggunakan metode analitik. Jikapun terdapat penyelesaiannya secara analitik, proses penyelesaiaannya sering kali cukup rumit dan memakan banyak waktu sehingga tidak efisien. Sebelum UTS kami mempelajari berbagai materi pokok mengenai metode numerik yaitu mencari akar akar persamaan non linear, turunan numerik, regresi dan interpolasi.

1. Mencari akar akar persamaan non linear

Pada materi ini saya mempelajari cara cara menghitung akar persamaan non liner secara numerik dengan 2 metode yaitu close methods dan open methods. Closed Methods merupakan metode pencarian akar-akar dengan menggunakan batas atas dan batas bawah untuk mencari akar dan mengukur persentasi error yang didapat untuk mendapat tingkat akurasi dari iterasi. Metode yang saya pelajari ada 3 yaitu Graphical Methods, Bisection Methods, dan False-Position Methods. Sedangkan Open Methods merupakan metode pencarian yang hanya menggunakan 1 titik untuk menemukan akar-akar. Metode ini dapat menggunakan turunan suatu fungsi untuk menentukan titik pengujian baru yang semakin dekat dengan nilai akar-akar yang diinginakn. Metode yang saya pelajari ada 3 yaitu Fixed-Point Iteration, Newton-Rapshon, dan Secant Method

2. Turunan numerik

Turunan Numerik adalah menentukan nilai turunan fungsi f yang diberikan dalam bentuk tabel. Terdapat 3 pendekatan dalam menghitung turunan numerik:

Turunan Numerik.png

3. Regresi

Regresi linear adalah sebuah pendekatan untuk memodelkan hubungan antara variable terikat Y dan satu atau lebih variable bebas yang disebut X. Salah satu kegunaan dari regresi linear adalah untuk melakukan prediksi berdasarkan data-data yang telah dimiliki sebelumnya. Hubungan di antara variable-variabel tersebut disebut sebagai model regresi linear.

4. Interpolasi

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

PR 1 Belajar OpenModelica

Pada tugas minggu ini saya menggunakan materi interpolasi dalam penerapan metode numerik menggunakan aplikasi open modelica. Saya menggunakan interpolasi untuk mencari data properties udara berupa massa jenis, kapasitas panas, viskositas dinamik, viskositas kinematik, konduktivitas termal, difusivitas termal dan prandtl number untuk udara dengan temperatur 100 derajat celcius atau setara dengan 373 derajat kelvin. Data data yang dikaetahu adalah sebagai berikut :


Tabel A4.png


Dari data yang ada pada tabel tersebut kita akan menggunakan properties pada temperatur 350K dan 400K untuk melakukan interpolasi. Data tersebut kita input ke aplikasi modelica dengan tambahan kata "parameter real" di depan.


Parameter Real.png


Setelah itu kita input data yang akan di hitung menggunakan kata "real" didepan data yang akan kita cari, lalu masukan persamaan yang akan digunakan untuk menghitung interpolasi berikut. Persamaan nya berupa (y-y1)/(y2-y1) = (x-x1)/(x2-x1)


Data Yang akan Dicari.png


Setelah itu kita simulasikan untuk mendapatkan hasil dari persamaan tersebut , sehingga mendapatkan angka sebagai berikut


Hasil Interpolasi.png


Untuk penjelasan lebih lengkapnya saya cantumkan dalam video berikut :

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


Pertemuan Metode Numerik 2 : 16 November 2020

Pada pertemuan kedua kami semua diperintahkan untuh mencoba melakukan simulasi simpel yaitu menjumlahkan suatu variabel dengan angka 10 ( y = x + 10 ). Hal yang pertama dilakukan adalah menginput parameter real yaitu berupa nilai x = 2, setelah itu kita menginput real yaitu berupa variabel yang nilai nya akan ditentukan, kemudian dilanjutkan dengan menginput equation atau persamaan y = 10 + x.


Lat 1.png


Setelah itu kita melakukan simulasi dan dihasilkan hasil sebagai berikut :


Lat 1.3.png


Setelah disimulasi kita masih bisa mengubah variabel x dengan nilai berapapun, dan melalukan simulasi ulang / resimulate sehingga akan menghasilkan nilai yang berbeda. Contohnya saya akan mengubah variabel x = 15 dan saya tunjukan pada gambar berikut


Lat 1.2.png


PR 2 Aljabar Simultan dengan OpenModelica

Pada petemuan kedua pak Dai memberikan PR kepada mahasiswa yaitu membuat class dengan type function untuk menyelesaikan persamaan aljabar simultan (metoda gauss elimination, gauss seidel ataupun metoda lain dan sebuah class untuk menjalankan fungsi tersebut. Dari PR yang diberikan saya akan mencoba menggunakan metode gauss elimination untuk mendapatkan solusi dari persamaan berikut dengan menggunakan aplikasi openmodelica

   w + 2x - y + z = 6
  -w + x + 2y - z = 3
  2w - x + 2y + 2z= 14
   w + x - y  + 2z= 11

Dari persamaan berikut saya modelkan ke aplikasi openmodelica menggunakan bantuan fungsi solve


FunctionSolve.png


Setelah itu saya membuat class dengan dengan meng input koefisien persamaan diatas kedalam bentuk matrix


Class1.png


Setelah itu dilakukan simulasi sehingga mendapatkan hasil sebagai berikut


Hasil Eliminasi.png


Pertemuan Metode Numerik 3 : 23 November 2020

Pada pertemuan ketiga pak DAI memerintahkan semua mahasiswa untuk mencoba membuktikan hasil deltaX dari suatu sistem pegas yang ada didalam buku Steven C Chapra pada sub bab 12.4 dengan free body diagram sistem sebagai berikut :


Pegas 1 Metnum.png


Pers 1 2 dan 3.png


Karena sistem berada dalam kondisi statis maka m d^2x/dt^2 = 0. Sehingga menghasilkan persamaan sebagai berikut


Pers Total.png


Data yang diketahui dari soal adalah sebagai berikut

  k = 10N/m
  m1 = 2kg 
  m2 = 3kg
  m3 = 2.5

Dengan mensubstitusi semua besaran yang diketahui ke persamaan total maka akan didapatkan matrix sebagai berikut


K dan W.png


Dari persamaan berikut saya modelkan ke aplikasi openmodelica menggunakan bantuan fungsi solve dengan persamaan dasar W = K * x


Function pegas.png


Dari matrix yang telah didapat, input ke dalam class openmodelica


Class pegas.png


Setelah itu semua disimulasi akan mendapatkan nilai x1, x2 dan x3 sebagai berikut


Hasil pegas.png

Dan hasil dari perhitungan di buku adalah sebagai berikut

Hasil pegas dari buku.png


Dari kedua hasil perhitungan diatas dapat disimpulkan bahwa hasil dari perhitungan di buku dan open modelica mempunyai hasilyang sama.

PR 3 Menghitung Dispalcement dan Gaya Reaksi pada Balconi Truss dengan OpenModelica

Pada petemuan kedua pak Dai memberikan PR kepada mahasiswa yaitu menghitung dispalcement dan gaya reaksi pada balconi truss dengan OpenModelica.


Truss 0.png


Dari PR yang diberikan saya akan mencoba menggunakan aljabar simultan untuk mencari nilai displacement dan gaya reaksi disetiap join. Hal yang pertama dilakukan adalah menentukan node, sudut terhadap sumbu X positif dan element pada setiap batang.

Truss 2.png

Setelah itu kita menentukan kekakuan setiap batang dengan rumus k = AE/L. Untuk hal itu kita harus menentukan luas penampang dan panjang setiap batang. Dari soal sudah diketahui modulus elastisitas setiap batang E = luas penampang setiap batang sebesar 8inch persegi dan panjang batang 1, 3, 4 dan 6 adalah 3 feet atau sebesar 36 inch dan sudut antara batang 4 dengan batang 2 dan 5 adalah 45 derajat. Oleh karena itu kita dapat menghitung panjang batang 2 dan 5 dengan 36/cos45 = 50.9 inch. Sehingga k untuk batang 1, 3, 4 dan 6 adalah

Truss 3.1.png

dan K untuk batang 2 dan 5 adalah :

Truss 3.2.png

Setelah itu kita membuat persamaan K dalam bentuk matrix untuk setiap elemen.Dengan rumus sebagai berikut

Truss 4.png


Analisa Batang 1, 3 dan 6

K matrix untuk element digunakan dengan cara menginput setiap sudut yang terbentuk pada batng ke dalam marix. Karena batang 1,3 dan 6 pada mempunyai K = 4,22 x 10^5 dan sudut 0 terhadap sumbu X maka ketiga pada batang tersebut mempunyai matrix K sebagai berikut :


Truss 5.1.png


Lalu diaplikasikan ke batang 1, 3 dan 6


Truss 6.1.png


Lalu ketiga matrix tersebut diinput ke dalam matrix global sesuai dengan posisi U


Truss 7.1.png


Analisa Batang 4

Batang 4 mempunyai K = 4,22 x 10^5 dan sudut 90 derajat terhadap sumbu X maka ketiga pada batang tersebut mempunyai matrix K sebagai berikut :


Truss 8.png


Lalu matrix tersebut diinput ke dalam matrix global sesuai dengan posisi U


Truss 9.png


Analisa Batang 2

Batang 2 mempunyai K = 2,98 x 10^5 dan sudut 135 derajat terhadap sumbu X maka ketiga pada batang tersebut mempunyai matrix K sebagai berikut :


Truss 10.png
Truss 11.png


Lalu matrix tersebut diinput ke dalam matrix global sesuai dengan posisi U


Truss 12.png


Analisa Batang 5

Batang 5 mempunyai K = 2,98 x 10^5 dan sudut 45 derajat terhadap sumbu X maka ketiga pada batang tersebut mempunyai matrix K sebagai berikut :


Truss 13.png


Lalu matrix tersebut diinput ke dalam matrix global sesuai dengan posisi U


Truss 14.png


Kemudian kita menjumlahkan semua matrix K global tersebut menjadi matrix K global total


Truss 15.png

Sehingga akan menjadi hasil seperti berikut


Truss 16.png


Tentukan kondisi batas dan kondisi pembebanan pada struktur. Dari struktur kita dapat mengetahui bahwa U1x, U1y, U3x, U3y = 0 dan beban sebesar 500 kearah bawah pada titik tumpuan 4 dan 5 ke arah sumbu Y. Karena nilai U1x, U1y, U3x, U3y = 0 maka kita bisa mengeliminasi baris dan kolom 1, 2, 6 dan 7, sehingga nantinya akan didapatkan matrix berukuran 6 x 6.


Truss 17.png

Setelah itu kita akan menghitung nilai U menggunakan simulasi pada aplikasi OpenModelica. Saya menggunakan konsep eliminasi gauss jordan dengan menginput nilai K dan F untuk mendapatkan nilai U. Berikut adalah fungsi dan class nya


Truss 18.png
Truss 19.png
Truss 20.png


Hasil dari simulasi OpenModelica adalah

  U2x = -0.0035545 inch
  U2y = -0.0102659 inch
  U4x = 0.00118483 inch
   
  U4y = -0.0114507 inch
  U5x = 0.00236967 inch
  U5y = -0.0195458 inch

Karena nilai U1x, U1y, U3x, U3y = 0, komponen U jika disusun akan menjadi seperti berikut


Truss 21.png


Setelah nilai U didapatkan, kita dapat menghitung besarnya gaya reaksi disetiap titik tumpuan dengan rumus R = Kglobal*U - F dalam bentuk matrix


Truss 22.png

Setelah itu kita akan menyelesaikan persamaan matrix tersebut menggunakan Openmodelica menggunakan class sebagai berikut


Truss 23.png

Sehingga menghasilkan nilai R sebagai berikut

Truss 24.png


Pada hasil dibuku nilai nilai R tersebut mengalami pembulatan, sehingga hasilnya menjadi seperti berikut


Truss 25.png


Pertemuan Metode Numerik 4 : 30 November 2020 Quiz 1

Pada pertemuan kali ini kami semua diberikan soal quiz oleh pak Dai . Soalnya adalah seperti berikut

Quiz 1.jpeg


Quiz 2.jpeg

Pertama kami diminta untuk membuat flowchart untuk proses pengerjaan soal. Flowchart yang digunakan untuk kedua soal ini sama, yaitu:

Metnum Flowchart.jpg


Pertama yang harus dilakukan adalah membuat tabel, berisi nama element, node, dan sudut setiap batang terhadap sumbu X positif. Setelah itu menentukan nilai kekakuan pada setiap batang.

2D.jpeg


Setelah itu kita membuat K lokal untuk setiap batang menggunakan fungsi dan class pada OpenModelica


K2D.png
K Lokal 2D1.png
K Lokal 2D2.png

Setelah didapatkan nilai K lokal, kita konversi menjadi K global pada setiap batang. Kemudian nilai K pada setiap batang ini dijumlahkan menjadi K global total. Disini saya menggunakan perhitungan manual menggunakan bantuan excel, sehingga menghasilkan sebagai berikut.

Exc1.png

Setelah dimasukan kondisi batas Ux1, Uy1, Ux3, Uy3 = 0. Maka matrixnya akan menjadi berukuran 4 x 4

Exc2.png

Sehingga untuk U dan R nya menjadi seperti berikut :

Gauss Jordan2D.png
U2D.png
Reaksi 2D.png



NO.8

Pertama yang harus dilakukan adalah membuat tabel, berisi nama element, node, dan sudut setiap batang terhadap sumbu X positif. Setelah itu menentukan nilai kekakuan pada setiap batang.


3D.jpeg


Setelah itu kita membuat K lokal untuk setiap batang menggunakan fungsi dan class pada OpenModelica


K3D1.png
K3D2.png
K Lokal 3D.png


Setelah didapatkan nilai K lokal, kita konversi menjadi K global pada setiap batang. Kemudian nilai K pada setiap batang ini dijumlahkan menjadi K global total. Disini saya menggunakan perhitungan manual menggunakan bantuan excel, sehingga menghasilkan sebagai berikut.


Exc3.png


Setelah dimasukan kondisi batas UxB, UyB, UzB, UxC, UyC, UzC, UxD, UyD, UzD= 0. Maka matrixnya akan menjadi berukuran 3 x 3


Exc4.png


Sehingga untuk U dan R nya menjadi seperti berikut :

Gauss Jordan3D.png
U3D.png
Reaksi 3D.png