|
|
| (8 intermediate revisions by the same user not shown) |
| Line 1: |
Line 1: |
| − | {| class="wikitable" style="float:right; margin-left: 10px;"
| |
| − | !colspan="2"|Profile
| |
| − | |-
| |
| − | |colspan="2"|[[File:Rezky.jpg|250px|Center]]
| |
| − | |-
| |
| − | |'''Nama Lengkap'''
| |
| − | |Rezky Alfian Fatra
| |
| − | |-
| |
| − | |'''NPM'''
| |
| − | |1906379043
| |
| − | |-
| |
| − | |style="vertical-align:top;"|'''Jurusan'''
| |
| − | |Teknik Mesin
| |
| − | |-
| |
| − | |colspan="2" style="text-align:center;width:200px;"|
| |
| − | |}
| |
| | | | |
| − |
| |
| − | '''Rezky Alfian Fatra''' (lahir di Jakarta, 13 November 2000; umur 20 tahun) adalah seorang mahasiswa yang mengambil jurusan Teknik Mesin di Universitas Indonesia angkatan 2019. Selain menjalankan tugas utama menjadi seorang mahasiswa, '''Rezky''' juga dikenal aktif berkontribusi di Lembaga Dakwah Kampus,khususnya di IMM yang memberikan wawasan berupa ilmu keagamaan kepada mahasiswa Departemen Teknik Mesin.Salah satu hal yang dilakukan oleh Rezky adalah dengan membuat program kerja yang menjadi wadah bagi mahasiswa Departemen Teknik Mesin untuk meningkatkan wawasan maupun kemampuan keagamaanya.
| |
| − |
| |
| − | __TOC__
| |
| − |
| |
| − | == Aktivitas ==
| |
| − |
| |
| − | Rezky Alfian Fatra atau panggilan akrabnya adalah Ekky adalah seorang mahasiswa yang aktif dalam melakukan kontribusi semenjak SMA, khususnya dalam bidang dakwah agar bisa memberikan manfaat bagi orang banyak.Ini terbukti dengan banyaknya organisasi rohani islam yang pernah ia jalani
| |
| − |
| |
| − | Sangat besar harapan saya agar terus bisa berkontribusi memberikan kebaikan yang bisa bermanfaat dan berguna bagi orang banyak, seperti yang pernah disabdakan oleh Nabi Muhammad SAW : "Sebaik-baiknya manusia adalah manusia yang bermanfaat untuk orang lain."
| |
| − |
| |
| − | Dalam aktivitasnya menjadi seorang pendakwah yang memberikan manfaat bagi orang banyak, Rezky juga tidak lupa untuk terus menambah ilmunya dengan mengikuti kajian keagamaan yang diselenggarakan oleh Manazil Ilmi Insani (Meninggikan Ilmu Manusia)
| |
| − | *https://www.youtube.com/channel/UCv_8XF7jxSQ27uhbAo_039w
| |
| − |
| |
| − |
| |
| − |
| |
| − | == Metode Numerik ==
| |
| − |
| |
| − | Metode numerik merupakan sebuah metode yang digunakan untuk menyelesaikan persoalan dimana perhitungan yang dilakukan secara analitik sudah sulit atau bahkan tidak bisa untuk dilakukan.Metode numerik ini berawal dari sebuah pemikiran bahwa permasalahan dapat diselesaikan dengan menggunakan pendekatan-pendekatan yang dapat dipertanggungjawabkan dengan metode analitik.Dengan menggunakan metode numerik, kita bisa mendapatkan hasil yang cepat dan tepat (meskipun dengan menggunakan pendekatan).
| |
| − |
| |
| − | Pada kelas Metode Numerik-02 semester ganjil 2020/2021, dari awal perkuliahan sampai UTS diisi oleh Bapak Dr. Ir. Engkos A. Kosasih, M.T.
| |
| − | Dimana pada pertemuan yang telah dilakukan telah membahas materi seperti:
| |
| − |
| |
| − | *Pseudocode
| |
| − | *Sistem Persamaan
| |
| − | *Turunan Numerik
| |
| − | *Metode Newton Rhapson, Secant, dan Bisection
| |
| − | *Regresi Linear
| |
| − | *Interpolasi
| |
| − |
| |
| − | Setelah UTS pembelajaran diisi oleh Bapak Dr. Ir. Ahmad Indra Siswantara
| |
| − |
| |
| − | == Pertemuan Perkuliahan ==
| |
| − |
| |
| − | === Pertemuan Pertama (Rabu, 11 November 2020) ===
| |
| − |
| |
| − | Pada pertemuan pertama yang diisi oleh Bapak Dr. Ir. Ahmad Indra Siswantara. Mahasiswa dijelaskan mengenai esensi atau tujuan pembelajaran metode numerik, yaitu:
| |
| − |
| |
| − | 1. Memahami konsep-konsep dan prinsip-prinsip dasar dalam metode numerik seperti:Persamaan algoritma, aljabar, pencocokan kurva, persamaan diferensial parsial, dan lainnya
| |
| − |
| |
| − | 2. Mengerti dan mampu menerapkan pemahaman atau aplikasi terhadap konsep metode numerik
| |
| − |
| |
| − | 3. Mampu menerapkan metode numerik dalam hal terkait persoalan keteknikan
| |
| − |
| |
| − | 4. Mendapat added value (nilai tambah)/adab sehingga mahasiswa menjadi orang yang memiliki budi pekerti yang baik
| |
| − |
| |
| − | Selain empat poin yang telah disebutkan, Mahasiswa juga mendapatkan pelajaran terkait penggunaan aplikasi OpenModelica. Dimana Mahasiswa dijelaskan dengan sangat sistematis mengenai penggunaan sederhana OpenModelica.Bahkan Bapak Dr. Ir. Ahmad Indra Siswantara juga menjelaskan bahwa komputer memiliki limit atau batasan yang tidak dapat dilampauinya.Pada misalnya saja komputer tidak mampu menentukan nilai 1/0 yang hanya Tuhan yang tahu nilainya. Pelajaran atau hikmah yang dapat diambil dari kasus ini adalah manusia haruslah beriman dengan Rabb-Nya dan jangan pernah melupakan Tuhan yang mengetahui dan berkuasa atas segala sesuatu.
| |
| − |
| |
| − | Tugas minggu pertama
| |
| − | [[File:Rezky Alfian Fatra Video1.mp4|720px|center]]
| |
| − |
| |
| − | === Pertemuan Kedua (Rabu, 18 November 2020) ===
| |
| − | Pada pertemuan minggu ke-2, beberapa mahasiswa dimintai untuk mempresentasikan hasil tugas yang diberikan pada minggu pertama.Selain itu, pada pertemuan kali ini, Bapak Dr. Ir. Ahmad Indra Siswantara juga mengajari mahasiswa mengenai fitur class dan function.Alasan Bapak Dr. Ir. Ahmad Indra Siswantara memilih untuk menggunakan open modelica, yaitu karena aplikasi ini lebih mudah untuk digunakan (tidak butuh spesifikasi yang tinggi) dan bisa digunakan secara gratis oleh mahasiswa.
| |
| − |
| |
| − | Untuk tugas minggu ke dua, mahasiswa diminta untuk mempelajari lagi open modelica dan melakukan sistem panggil class dan function.Berikut tugas saya dengan menggunakan contoh soal Aljabar Linier dengan ordo 4x4
| |
| − |
| |
| − | [[File:Rezky Alfian Fatra Video2.mp4|720px|center]]
| |
| − |
| |
| − | === Pertemuan Ketiga (Rabu, 25 November 2020) ===
| |
| − |
| |
| − | '''Tugas 3'''
| |
| − | *Berikut adalah soal untuk tugas kali ini:
| |
| − | *Penyelesaian:
| |
| − | [[File:Soal Trusses 2 Rezky.jpg|480px|thumb|Soal Trusses]]
| |
| − |
| |
| − | '''Persamaan'''
| |
| − |
| |
| − | class Trusses_Tugas3_RezkyAlfianFatra
| |
| − |
| |
| − | parameter Integer N=8; //Global matrice = 2*points connected
| |
| − | parameter Real A=0.001; //Area m2
| |
| − | parameter Real E=200e9; //Pa
| |
| − | Real G[N,N]; //global
| |
| − | Real Ginitial[N,N]; //global
| |
| − | Real Sol[N]; //global dispplacement
| |
| − | Real X[N]={0,0,-1035.2762,-3863.7033,0,0,-1035.2762,-3863.7033};
| |
| − | Real R[N]; //global reaction force
| |
| − | Real SolMat[N,1];
| |
| − | Real XMat[N,1];
| |
| − |
| |
| − | //boundary condition
| |
| − | Integer b1=1;
| |
| − | Integer b2=3;
| |
| − |
| |
| − | //truss 1
| |
| − | parameter Real X1=0; //degree between truss
| |
| − | Real k1=A*E/1;
| |
| − | Real K1[4,4]; //stiffness matrice
| |
| − | Integer p1a=1;
| |
| − | Integer p1b=2;
| |
| − | Real G1[N,N];
| |
| − |
| |
| − | //truss 2
| |
| − | parameter Real X2=0; //degree between truss
| |
| − | Real k2=A*E/1;
| |
| − | Real K2[4,4]; //stiffness matrice
| |
| − | Integer p2a=2;
| |
| − | Integer p2b=3;
| |
| − | Real G2[N,N];
| |
| − |
| |
| − | //truss 3
| |
| − | parameter Real X3=90; //degree between truss
| |
| − | Real k3=A*E/1.25;
| |
| − | Real K3[4,4]; //stiffness matrice
| |
| − | Integer p3a=2;
| |
| − | Integer p3b=4;
| |
| − | Real G3[N,N];
| |
| − |
| |
| − | //truss 4
| |
| − | parameter Real X4=90+38.6598; //degree between truss
| |
| − | Real k4=A*E/1.6;
| |
| − | Real K4[4,4]; //stiffness matrice
| |
| − | Integer p4a=1;
| |
| − | Integer p4b=4;
| |
| − | Real G4[N,N];
| |
| − |
| |
| − | //truss 5
| |
| − | parameter Real X5=90-38.6598; //degree between truss
| |
| − | Real k5=A*E/1.6;
| |
| − | Real K5[4,4]; //stiffness matrice
| |
| − | Integer p5a=3;
| |
| − | Integer p5b=4;
| |
| − | Real G5[N,N];
| |
| − |
| |
| − | /*
| |
| − | for each truss, please ensure pXa is lower then pXb (X represents truss element number)
| |
| − | */
| |
| − |
| |
| − | algorithm
| |
| − |
| |
| − | //creating global matrice
| |
| − | K1:=Stiffness_Matrices(X1);
| |
| − | G1:=k1*Local_Global(K1,N,p1a,p1b);
| |
| − |
| |
| − | K2:=Stiffness_Matrices(X2);
| |
| − | G2:=k2*Local_Global(K2,N,p2a,p2b);
| |
| − |
| |
| − | K3:=Stiffness_Matrices(X3);
| |
| − | G3:=k3*Local_Global(K3,N,p3a,p3b);
| |
| − |
| |
| − | K4:=Stiffness_Matrices(X4);
| |
| − | G4:=k4*Local_Global(K4,N,p4a,p4b);
| |
| − |
| |
| − | K5:=Stiffness_Matrices(X5);
| |
| − | G5:=k5*Local_Global(K5,N,p5a,p5b);
| |
| − |
| |
| − | G:=G1+G2+G3+G4+G5;
| |
| − | Ginitial:=G;
| |
| − |
| |
| − | //implementing boundary condition
| |
| − | for i in 1:N loop
| |
| − | G[2*b1-1,i]:=0;
| |
| − | G[2*b1,i]:=0;
| |
| − | G[2*b2-1,i]:=0;
| |
| − | G[2*b2,i]:=0;
| |
| − | end for;
| |
| − |
| |
| − | G[2*b1-1,2*b1-1]:=1;
| |
| − | G[2*b1,2*b1]:=1;
| |
| − | G[2*b2-1,2*b2-1]:=1;
| |
| − | G[2*b2,2*b2]:=1;
| |
| − |
| |
| − | //solving displacement
| |
| − | Sol:=Gauss_Jordan(N,G,X);
| |
| − |
| |
| − | //solving reaction force
| |
| − | SolMat:=matrix(Sol);
| |
| − | XMat:=matrix(X);
| |
| − | R:=Reaction_Trusses(N,Ginitial,SolMat,XMat);
| |
| − |
| |
| − | end Trusses_Tugas3_RezkyAlfianFatra;
| |
| − |
| |
| − | '''Grafik Displacement'''
| |
| − | [[File:Grafik Displacement Rezky.jpg|480px|Center]]
| |
| − |
| |
| − |
| |
| − | '''Grafik Reaction Forces'''
| |
| − | [[File:Trusses 2 Reactiion Rezky.jpg|480px|Center]]
| |
| − |
| |
| − |
| |
| − |
| |
| − | '''Fungsi Panggil'''
| |
| − | {| class="wikitable"
| |
| − | |-
| |
| − | | style='border-style: none none solid solid;' |
| |
| − | ''Matrice Transformation''
| |
| − |
| |
| − | function Stiffness_Matrices
| |
| − | input Real A;
| |
| − | Real Y;
| |
| − | output Real X[4,4];
| |
| − | Real float_error = 10e-10;
| |
| − |
| |
| − | final constant Real pi=2*Modelica.Math.asin(1.0);
| |
| − |
| |
| − | algorithm
| |
| − |
| |
| − | Y:=A/180*pi;
| |
| − |
| |
| − | X:=[(Modelica.Math.cos(Y))^2,Modelica.Math.cos(Y)*Modelica.Math.sin(Y),-(Modelica.Math.cos(Y))^2,-Modelica.Math.cos(Y)*Modelica.Math.sin(Y);
| |
| − |
| |
| − | Modelica.Math.cos(Y)*Modelica.Math.sin(Y),(Modelica.Math.sin(Y))^2,-Modelica.Math.cos(Y)*Modelica.Math.sin(Y),-(Modelica.Math.sin(Y))^2;
| |
| − |
| |
| − | -(Modelica.Math.cos(Y))^2,-Modelica.Math.cos(Y)*Modelica.Math.sin(Y),(Modelica.Math.cos(Y))^2,Modelica.Math.cos(Y)*Modelica.Math.sin(Y);
| |
| − |
| |
| − | -Modelica.Math.cos(Y)*Modelica.Math.sin(Y),-(Modelica.Math.sin(Y))^2,Modelica.Math.cos(Y)*Modelica.Math.sin(Y),(Modelica.Math.sin(Y))^2];
| |
| − |
| |
| − | for i in 1:4 loop
| |
| − | for j in 1:4 loop
| |
| − | if abs(X[i,j]) <= float_error then
| |
| − | X[i,j] := 0;
| |
| − | end if;
| |
| − | end for;
| |
| − | end for;
| |
| − |
| |
| − | end Stiffness_Matrices;
| |
| − |
| |
| − | | style="width: 20cm;"|
| |
| − | ''Global Element Matrice''
| |
| − |
| |
| − | function Local_Global
| |
| − | input Real Y[4,4];
| |
| − | input Integer B;
| |
| − | input Integer p1;
| |
| − | input Integer p2;
| |
| − | output Real G[B,B];
| |
| − |
| |
| − | algorithm
| |
| − |
| |
| − | for i in 1:B loop
| |
| − | for j in 1:B loop
| |
| − | G[i,j]:=0;
| |
| − | end for;
| |
| − | end for;
| |
| − |
| |
| − | G[2*p1,2*p1]:=Y[2,2];
| |
| − | G[2*p1-1,2*p1-1]:=Y[1,1];
| |
| − | G[2*p1,2*p1-1]:=Y[2,1];
| |
| − | G[2*p1-1,2*p1]:=Y[1,2];
| |
| − |
| |
| − | G[2*p2,2*p2]:=Y[4,4];
| |
| − | G[2*p2-1,2*p2-1]:=Y[3,3];
| |
| − | G[2*p2,2*p2-1]:=Y[4,3];
| |
| − | G[2*p2-1,2*p2]:=Y[3,4];
| |
| − |
| |
| − | G[2*p2,2*p1]:=Y[4,2];
| |
| − | G[2*p2-1,2*p1-1]:=Y[3,1];
| |
| − | G[2*p2,2*p1-1]:=Y[4,1];
| |
| − | G[2*p2-1,2*p1]:=Y[3,2];
| |
| − |
| |
| − | G[2*p1,2*p2]:=Y[2,4];
| |
| − | G[2*p1-1,2*p2-1]:=Y[1,3];
| |
| − | G[2*p1,2*p2-1]:=Y[2,3];
| |
| − | G[2*p1-1,2*p2]:=Y[1,4];
| |
| − |
| |
| − | end Local_Global;
| |
| − |
| |
| − | | style="width: 20cm;"|
| |
| − | ''Reaction Matrice Equation''
| |
| − |
| |
| − | function Reaction_Trusses
| |
| − | input Integer N;
| |
| − | input Real A[N,N];
| |
| − | input Real B[N,1];
| |
| − | input Real C[N,1];
| |
| − | Real X[N,1];
| |
| − | output Real Sol[N];
| |
| − | Real float_error = 10e-10;
| |
| − |
| |
| − | algorithm
| |
| − | X:=A*B-C;
| |
| − |
| |
| − | for i in 1:N loop
| |
| − | if abs(X[i,1]) <= float_error then
| |
| − | X[i,1] := 0;
| |
| − | end if;
| |
| − | end for;
| |
| − |
| |
| − | for i in 1:N loop
| |
| − | Sol[i]:=X[i,1];
| |
| − | end for;
| |
| − |
| |
| − | end Reaction_Trusses;
| |
| − |
| |
| − | |}
| |
| − |
| |
| − | '''Gauss Jordan'''
| |
| − | {| class="wikitable"
| |
| − | |-
| |
| − | | style='border-style: none none solid solid;' |
| |
| − |
| |
| − | function Gauss_Jordan
| |
| − | input Integer N;
| |
| − | input Real A[N,N];
| |
| − | input Real B[N];
| |
| − | output Real X[N];
| |
| − | Real float_error = 10e-10;
| |
| − |
| |
| − | algorithm
| |
| − | X:=Modelica.Math.Matrices.solve(A,B);
| |
| − |
| |
| − | for i in 1:N loop
| |
| − | if abs(X[i]) <= float_error then
| |
| − | X[i] := 0;
| |
| − | end if;
| |
| − | end for;
| |
| − |
| |
| − | end Gauss_Jordan;
| |
| − |
| |
| − | |}
| |
| − |
| |
| − | === Pertemuan Keempat (Rabu, 02 Desember 2020) ===
| |
| − | [[File:Flowchart Rezky.jpg|480px|center]]
| |
| − |
| |
| − | ===Tugas 4===
| |
| − | Pr No 8
| |
| − | kami diminta untuk menyelesaikan permasalahan 3D dengan membuat class diagram, flowchart, dan penyelesaiannya dengan openmodelica
| |
| − | [[File:Prno8.png|720px|center]]
| |
| − | Pertama saya membuat class diagram
| |
| − | [[File:PR8 RezkyA.jpeg|400px|center]]
| |
| − | Kemudian saya membuat flowchart
| |
| − | [[File:Pr8bag2 Rezky.jpeg|400px|center]]
| |
| − |
| |
| − | Setelah itu, saya membuat Free Body Diagram (FBD) untuk mencari data-data yang diperlukan pada persoalan ini.Data tersbut berupa node, modulus elastisitas, luas, sudut, dan panjang dari setiap batang.Karena pada kasus ini merupakan permasalahan 3D (3 Dimensi), maka sudut akan berada pada arah X, Y, dan Z.Rumus untuk mencari panjang batang, bisa kita dapatkan di buku.Untuk gaya eksternal yang ada, kita dapat melihatnya dari gambar yang hanya berada di node 1 arah sumbu y.Kemudian setelah itu, data-data yang ada kita ubah menjadi bentuk matriks agar bisa diselesaikan dengan menggunakan openmodelica
| |
| − | [[File:Pr8bag3 Rezky.jpeg|400px|center]]
| |
| − |
| |
| − | Untuk mendapatkan nilai K dapat menggunakan rumus K=AE/L.Setelah mendapatkan nilai K pada setiap batang, nilai tersebut digunakan untuk membuat matriks K elemental (Ke) pada setiap batang
| |
| − | [[File:MatriksKekakuanElemen Rezky.png|600px|center]]
| |
| − |
| |
| − | Kemudian matriks 6x6 tersebut akan dibuat menjadi matrks 12x12 sesuai dengan jumlah batang yang ada (3)*jumlah node (4).Ini disebut dengan matriks kekakuan global (Kg)
| |
| − | [[File:MatriksKekakuanGlobal Rezky.png|600px|center]]
| |
| − |
| |
| − | Selanjutnya matriks kekakuan pada setiap batang dijumlahkan agar mendapatkan matriks kekakuan global total
| |
| − | [[File:JumlahMatriksKekakuanGlobal Rezky.png|600px|center]]
| |
| − |
| |
| − | Setelah itu, boundary diterapkan pada matriks kekakuan global total.Pada kasus yang ini, node 2, node 3, dan node 4 merupakan pinned. Sehingga U2X=U2Y=U2Z=U3X=U3Y=U3Z=U4X=U4Y=U4Z=0.
| |
| − | [[File:BatasKekakuanMatriksGlobal Rezky.png|600px|center]]
| |
| − |
| |
| − | Untuk mencari nilai U atau nilai displacement lainnya menggunakan Eliminasi Gauss Jordan dengan menggunakan persamaan XU=F. Dengan X adalah matriks K global yang sudah menerapkan boundary (KgB) dan F didefinisikan sebagai eksternal load
| |
| − | [[File:GaussJordan Rezky.png|600px|center]]
| |
| − |
| |
| − | Sedangkan untuk mencari reaction force dapat menggunakan persamaan R=Kg*U-F
| |
| − | [[File:ReactionForce Rezky.png|600px|center]]
| |
| − |
| |
| − | Function-function tersebut akan dipanggil ke dalam sebuah class untuk penyelesaiannya. Pada class juga akan diinput data-data yang telah dicari pada langkah pertama
| |
| − | [[File:Tugas8 RezkyAlfian.png|600px|center]]
| |
| − |
| |
| − | Sehingga, setelah disimulate, kita akan mendapatkan nilai displacement (U)
| |
| − | [[File:NilaiU Rezky.png|200px|center]]
| |
| − |
| |
| − | dan Nilai R
| |
| − | [[File:NilaiR Rezky.png|200px|center]]
| |
| − |
| |
| − | Berikut video penjelasannya
| |
| − | [[File:Video RezkyAlfian No8.mp4|720px|center]]
| |
| − |
| |
| − | === Pertemuan Keenam ===
| |
| − | Sistem optimasi
| |
| − | dalam sebuah sistem fluida ada yang berupa zat cair seperti aliran darah, dan juga ada yang berupa gas seperti sistem pernapasan.Yang ditampilkan pada grafik dibawah ini adalah sebuah karakteristik dari sistem dan komponen.
| |
| − | sistem fluida tidak hanya perpipaan, tapi juga mesin kerja berupa pompa dan tenaga, ada valve, dan juga ada pipa belok.Sistem fluida, aliran air dimasukan dari yang tediri pompa, pipa, pipa belok (T,L) valve.Pada setiap komponen memiliki tujuan yang sama, yaitu mengalirkan air dari sebuah tempat ke tempat yang kita tuju.Sistem tersebut adalah pompa, dengan menggunakan sistem pendukung, memberikan tenaga pada pompa
| |
| − |
| |
| − | [[File:Grafikfluid.jpeg|720px|center]]
| |
| − |
| |
| − | *Yang melengkung pada grafik tersebut adalah efisiensi
| |
| − | *Sumbu datar adalah jumlah aliran yang diantarkan
| |
| − | *Sumbu y adalah efisiensi dari pompa
| |
| − | *Garis hijau merupakan tenaga air yang keluar dari pompa
| |
| − | *Pada saat jumlah air yang dikeluarkan tekanan akan turun dan breakhouse power akan naik
| |
| − | *Pada efisiensi, makin tinggi kapasitasnya maka efisiensinya akan naik. Terdapat juga titik maksimum dan kemudan setelah itu efisiensi akan menurun
| |
| − | *Untuk efisiensi pada sistem, mempertimbangkan heatloss pada kondisi tertentu
| |
| − |
| |
| − | Pembahasan berikutnya mengenai metode Bracket Optimization Using Golden Ratio Method yang dijelaskan oleh bu Chandra.Jika ingin menyelesaikan suatu permasalahan dengan menggunakan cara optimasi, kita harus paham terlebih dahulu mengenai materinya
| |
| − |
| |
| − | [[File:Bracktergraph.png|500px|center]]
| |
| − | *Menentukan dan mengetahui batasan titik awal terendah (xL) dan tertinggi (xH)
| |
| − | *Mendapatkan nilai x1 dan x2 dari golden ratio (d)
| |
| − | *Menentukan xL dan xH berdasarkan nilai batasan maksimum dan minimal baru
| |
| − | *Mendapatkan nilai x1 dan x2 baru
| |
| − |
| |
| − | Jika dalam soal diketahui untuk menyelesaikan persamaan y=2sinx-x^2/10,kemudian fungsi ini dimasukkan ke dalam openmodelica, kemudian fungsi tersebut akan dipanggil
| |
| − | [[File:Rezky Metnum.png|720px|center]]
| |
| − |
| |
| − | Setelah di cek dan tidak terjadi masalah dalam simulate openmodelica, dapat dibuat model yang akan memanggil fungsi tersebut
| |
| − |
| |
| − | [[File:Brackopt.png|720px|center]]
| |
| − |
| |
| − | == Progress Tugas Besar ==
| |
| − | Tugas besar yang diberikan kemarin untuk kelas metnum 02 dan 03 adalah menghitung dengan menggunnakan optimasi untuk rangka batang:
| |
| − | [[File:600px-Soaltubes.jpeg|720px|center]]
| |
| − | suatu rangka batang yang telah dibentuk seperti gambar diatas, dengan spesifikasi sebagai berikut:
| |
| − | [[File:Soaltubes2.jpeg|720px|center]]
| |
| − | Diharapkan dari persoalan di atas dapat mendesain rangka batang dengan harga yang terjangkau dan optimal.
| |
| − |
| |
| − | Faktor atau variabel bebasnya meliputi : Harga, Material yang digunakan, dan cross section
| |
| − |
| |
| − | Dari soal yang diberikan, mahasiswa diharuskan mendesign rangka pada batang dengan harga yang terjangkau dan optimal
| |
| − |
| |
| − | Pada tugas besar ini, kami diberikan sebuah tugas untuk melakukan perhitungan pemilihan material dan luas penampang pada trusses yang akan digunakan untuk membuat sebuah konstruksi pada gambar di atas :
| |
| − |
| |
| − | 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 terjadinya buckling (pada truss paling atas).
| |
| − | *Ketinggian trusses untuk setiap lantai sama yaitu 0,6 m.
| |
| − |
| |
| − | Untuk jenis material yang sama:
| |
| − | #Mencari harga untuk 6 ukuran batang dengan material yang sama yaitu dengan menggunakan SS304.
| |
| − | #Mengitung nilai safety factor yang digunakan pada 6 ukuran batang.
| |
| − | #Membuat rasio antara safety factor dengan harga total batang.
| |
| − | #Membuat persamaan antara rasio (safe/harga) dengan area dengan menggunakan metode curve fitting.
| |
| − | #Melakukan optimasi menggunakan metode golden section.
| |
| − |
| |
| − | Untuk area penampang yang sama:
| |
| − | #Mencari harga untuk 4 jenis material dengan area penampang yang sama, yaitu 224 mm^2.
| |
| − | #Mengitung nilai safety factor pada 4 jenis batang.
| |
| − | #Membuat rasio antara safety factor dengan harga total.
| |
| − | #Membuat persamaan antara rasio (safe/harga) dengan jenis material dengan metode curve fitting.
| |
| − | #Melakukan optimasi menggunakan metode golden section.
| |
| − |
| |
| − | '''Coding yang digunakan'''
| |
| − | ----
| |
| − | *Untuk perhitungan displacement, reaction force, stress, dan safety factor:
| |
| − | model Trusses_3D_Tugas_Besar_Safety
| |
| − |
| |
| − | //define initial variable
| |
| − | parameter Integer Points=size(P,1); //Number of Points
| |
| − | parameter Integer Trusses=size(C,1); //Number of Trusses
| |
| − | parameter Real Yield=215e6; //Yield Strength (Pa)
| |
| − | parameter Real Area=0.000224; //Area L Profile (Dimension=0.03, Thickness=0,004) (m2)
| |
| − | parameter Real Elas=193e6; //Elasticity SS 304 (Pa)
| |
| − |
| |
| − | //define connection
| |
| − | parameter Integer C[:,2]=[1,5;
| |
| − | 2,6;
| |
| − | 3,7;
| |
| − | 4,8;
| |
| − | 5,6; //1st floor
| |
| − | 6,7; //1st floor
| |
| − | 7,8; //1st floor
| |
| − | 5,8; //1st floor
| |
| − | 5,9;
| |
| − | 6,10;
| |
| − | 7,11;
| |
| − | 8,12;
| |
| − | 9,10; //2nd floor
| |
| − | 10,11;//2nd floor
| |
| − | 11,12;//2nd floor
| |
| − | 9,12; //2nd floor
| |
| − | 9,13;
| |
| − | 10,14;
| |
| − | 11,15;
| |
| − | 12,16;
| |
| − | 13,14;//3rd floor
| |
| − | 14,15;//3rd floor
| |
| − | 15,16;//3rd floor
| |
| − | 13,16];//3rd floor
| |
| − |
| |
| − | //define coordinates (please put orderly)
| |
| − | parameter Real P[:,6]=[0.3,-0.375,0,1,1,1; //1
| |
| − | -0.3,-0.375,0,1,1,1; //2
| |
| − | -0.3,0.375,0,1,1,1; //3
| |
| − | 0.3,0.375,0,1,1,1; //4
| |
| − |
| |
| − | 0.3,-0.375,0.6,0,0,0; //5
| |
| − | -0.3,-0.375,0.6,0,0,0; //6
| |
| − | -0.3,0.375,0.6,0,0,0; //7
| |
| − | 0.3,0.375,0.6,0,0,0; //8
| |
| − |
| |
| − | 0.3,-0.375,1.2,0,0,0; //9
| |
| − | -0.3,-0.375,1.2,0,0,0; //10
| |
| − | -0.3,0.375,1.2,0,0,0; //11
| |
| − | 0.3,0.375,1.2,0,0,0; //12
| |
| − |
| |
| − | 0.3,-0.375,1.8,0,0,0; //13
| |
| − | -0.3,-0.375,1.8,0,0,0; //14
| |
| − | -0.3,0.375,1.8,0,0,0; //15
| |
| − | 0.3,0.375,1.8,0,0,0]; //16
| |
| − |
| |
| − | //define external force (please put orderly)
| |
| − | parameter Real F[Points*3]={0,0,0,
| |
| − | 0,0,0,
| |
| − | 0,0,0,
| |
| − | 0,0,0,
| |
| − | 0,0,0,
| |
| − | 0,0,0,
| |
| − | 0,0,0,
| |
| − | 0,0,0,
| |
| − | 0,0,0,
| |
| − | 0,0,0,
| |
| − | 0,0,0,
| |
| − | 0,0,0,
| |
| − | 0,0,-500,
| |
| − | 0,0,-1000,
| |
| − | 0,0,-1000,
| |
| − | 0,0,-500};
| |
| − |
| |
| − | //solution
| |
| − | Real displacement[N], reaction[N];
| |
| − | Real check[3];
| |
| − |
| |
| − | Real stress1[Trusses];
| |
| − | Real safety[Trusses];
| |
| − | Real dis[3];
| |
| − | Real Str[3];
| |
| − |
| |
| − | protected
| |
| − | parameter Integer N=3*Points;
| |
| − | Real q1[3], q2[3], g[N,N], G[N,N], G_star[N,N], id[N,N]=identity(N), cx, cy, cz, L, X[3,3];
| |
| − | Real err=10e-10, ers=10e-4;
| |
| − |
| |
| − | algorithm
| |
| − | //Creating Global Matrix
| |
| − | G:=id;
| |
| − | for i in 1:Trusses loop
| |
| − | for j in 1:3 loop
| |
| − | q1[j]:=P[C[i,1],j];
| |
| − | q2[j]:=P[C[i,2],j];
| |
| − | end for;
| |
| − |
| |
| − | //Solving Matrix
| |
| − | L:=Modelica.Math.Vectors.length(q2-q1);
| |
| − | cx:=(q2[1]-q1[1])/L;
| |
| − | cy:=(q2[2]-q1[2])/L;
| |
| − | cz:=(q2[3]-q1[3])/L;
| |
| − | X:=(Area*Elas/L)*[cx^2,cx*cy,cx*cz;
| |
| − | cy*cx,cy^2,cy*cz;
| |
| − | cz*cx,cz*cy,cz^2];
| |
| − |
| |
| − | //Transforming to global matrix
| |
| − | g:=zeros(N,N);
| |
| − | for m,n in 1:3 loop
| |
| − | g[3*(C[i,1]-1)+m,3*(C[i,1]-1)+n]:=X[m,n];
| |
| − | g[3*(C[i,2]-1)+m,3*(C[i,2]-1)+n]:=X[m,n];
| |
| − | g[3*(C[i,2]-1)+m,3*(C[i,1]-1)+n]:=-X[m,n];
| |
| − | g[3*(C[i,1]-1)+m,3*(C[i,2]-1)+n]:=-X[m,n];
| |
| − | end for;
| |
| − |
| |
| − | G_star:=G+g;
| |
| − | G:=G_star;
| |
| − | end for;
| |
| − |
| |
| − | //Implementing boundary
| |
| − | for x in 1:Points loop
| |
| − | if P[x,4] <> 0 then
| |
| − | for a in 1:Points*3 loop
| |
| − | G[(x*3)-2,a]:=0;
| |
| − | G[(x*3)-2,(x*3)-2]:=1;
| |
| − | end for;
| |
| − | end if;
| |
| − | if P[x,5] <> 0 then
| |
| − | for a in 1:Points*3 loop
| |
| − | G[(x*3)-1,a]:=0;
| |
| − | G[(x*3)-1,(x*3)-1]:=1;
| |
| − | end for;
| |
| − | end if;
| |
| − | if P[x,6] <> 0 then
| |
| − | for a in 1:Points*3 loop
| |
| − | G[x*3,a]:=0;
| |
| − | G[x*3,x*3]:=1;
| |
| − | end for;
| |
| − | end if;
| |
| − | end for;
| |
| − |
| |
| − | //Solving displacement
| |
| − | displacement:=Modelica.Math.Matrices.solve(G,F);
| |
| − |
| |
| − | //Solving reaction
| |
| − | reaction:=(G_star*displacement)-F;
| |
| − |
| |
| − | //Eliminating float error
| |
| − | for i in 1:N loop
| |
| − | reaction[i]:=if abs(reaction[i])<=err then 0 else reaction[i];
| |
| − | displacement[i]:=if abs(displacement[i])<=err then 0 else displacement[i];
| |
| − | end for;
| |
| − |
| |
| − | //Checking Force
| |
| − | check[1]:=sum({reaction[i] for i in (1:3:(N-2))})+sum({F[i] for i in (1:3:(N-2))});
| |
| − | check[2]:=sum({reaction[i] for i in (2:3:(N-1))})+sum({F[i] for i in (2:3:(N-1))});
| |
| − | check[3]:=sum({reaction[i] for i in (3:3:N)})+sum({F[i] for i in (3:3:N)});
| |
| − |
| |
| − | for i in 1:3 loop
| |
| − | check[i] := if abs(check[i])<=ers then 0 else check[i];
| |
| − | end for;
| |
| − |
| |
| − | //Calculating stress in each truss
| |
| − | for i in 1:Trusses loop
| |
| − | for j in 1:3 loop
| |
| − | q1[j]:=P[C[i,1],j];
| |
| − | q2[j]:=P[C[i,2],j];
| |
| − | dis[j]:=abs(displacement[3*(C[i,1]-1)+j]-displacement[3*(C[i,2]-1)+j]);
| |
| − | end for;
| |
| − |
| |
| − | //Solving Matrix
| |
| − | L:=Modelica.Math.Vectors.length(q2-q1);
| |
| − | cx:=(q2[1]-q1[1])/L;
| |
| − | cy:=(q2[2]-q1[2])/L;
| |
| − | cz:=(q2[3]-q1[3])/L;
| |
| − | X:=(Elas/L)*[cx^2,cx*cy,cx*cz;
| |
| − | cy*cx,cy^2,cy*cz;
| |
| − | cz*cx,cz*cy,cz^2];
| |
| − |
| |
| − | Str:=(X*dis);
| |
| − | stress1[i]:=Modelica.Math.Vectors.length(Str);
| |
| − | end for;
| |
| − |
| |
| − | //Safety factor
| |
| − | for i in 1:Trusses loop
| |
| − | if stress1[i]>0 then
| |
| − | safety[i]:=Yield/stress1[i];
| |
| − | else
| |
| − | safety[i]:=0;
| |
| − | end if;
| |
| − | end for;
| |
| − |
| |
| − | end Trusses_3D_Tugas_Besar_Safety;
| |
| − |
| |
| − | Berikut ini adalah variasi data untuk material yang sama, yaitu SS304 dan menggunakan area factor yang sama, yaitu 0.000224
| |
| − | [[File:TB Metnum1.png|500px|center]]
| |
| − |
| |
| − | Setelah melakukan simulate, maka kita akan mendapatkan nilai dari safety factor, kemudian kita mencari Rasio dengan melakukan perbandingan Safety factor dengan total harga, maka akan didapatkan :
| |
| − | [[File:TB Metnum2.png|500px|center]]
| |
| − |
| |
| − | Kemudian, berdasarkan hasil data yang telah didapat, kita bisa mendapatkan Luas area dan material yang baik digunakan pada rangka yang akan dibuat dengan menggunakan optimasi dengan curve fitting dan golden bracket.
| |
| − | Fungsi untuk Curve Fittingnya adalah sebagai berikut
| |
| − |
| |
| − | {| class="wikitable"
| |
| − | |-
| |
| − | | style='border-style: none none solid solid;' |
| |
| − |
| |
| − | function Curve_Fitting
| |
| − |
| |
| − | input Real X[:];
| |
| − | input Real Y[size(X,1)];
| |
| − | input Integer order=2;
| |
| − | output Real Coe[order+1];
| |
| − |
| |
| − | protected
| |
| − | Real Z[size(X,1),order+1];
| |
| − | Real ZTr[order+1,size(X,1)];
| |
| − | Real A[order+1,order+1];
| |
| − | Real B[order+1];
| |
| − |
| |
| − | algorithm
| |
| − |
| |
| − | for i in 1:size(X,1) loop
| |
| − | for j in 1:(order+1) loop
| |
| − | Z[i,j]:=X[i]^(order+1-j);
| |
| − | end for;
| |
| − | end for;
| |
| − | ZTr:=transpose(Z);
| |
| − |
| |
| − | A:=ZTr*Z;
| |
| − | B:=ZTr*Y;
| |
| − | Coe:=Modelica.Math.Matrices.solve(A,B);
| |
| − | //Coe:=fill(2,size(Coe,1));
| |
| − |
| |
| − | end Curve_Fitting;
| |
| − |
| |
| − | |}
| |
| − |
| |
| − | saya melakukan optimasi golden ratio untuk data material tetap
| |
| − |
| |
| − | {| class="wikitable"
| |
| − | |-
| |
| − | | style='border-style: none none solid solid;' |
| |
| − |
| |
| − | model Opt_Gold
| |
| − |
| |
| − | parameter Real xd[:]={0.000111,0.000141,0.000171,0.000231,0.000304,0.000375};
| |
| − | parameter Real yd[size(xd,1)]={678148,829776,985330,1308215,1722235,2147211};
| |
| − | parameter Real xlo=193e9;
| |
| − | parameter Real xhi=197e9;
| |
| − | parameter Integer N=10; // maximum iteration
| |
| − | parameter Real es=0.0001; // maximum error
| |
| − |
| |
| − | Real f1[N], f2[N], x1[N], x2[N], ea[N], y[3];
| |
| − | Real xopt, fx;
| |
| − | protected
| |
| − | Real d, xl, xu, xint, R=(5^(1/2)-1)/2;
| |
| − |
| |
| − | algorithm
| |
| − | xl := xlo;
| |
| − | xu := xhi;
| |
| − | y := Curve_Fitting(xd,yd);
| |
| − |
| |
| − | for i in 1:N loop
| |
| − | d:= R*(xu-xl);
| |
| − | x1[i]:=xl+d;
| |
| − | x2[i]:=xu-d;
| |
| − | f1[i]:=y[1]*x1[i]^2+y[2]*x1[i]+y[3];
| |
| − | f2[i]:=y[1]*x2[i]^2+y[2]*x2[i]+y[3];
| |
| − | xint:=xu-xl;
| |
| − |
| |
| − | if f1[i]>f2[i] then
| |
| − | xl:=x2[i];
| |
| − | xopt:=x1[i];
| |
| − | fx:=f1[i];
| |
| − | else
| |
| − | xu:=x1[i];
| |
| − | xopt:=x2[i];
| |
| − | fx:=f2[i];
| |
| − | end if;
| |
| − |
| |
| − | ea[i]:=(1-R)*abs((xint)/xopt);
| |
| − | if ea[i]<es then
| |
| − | break;
| |
| − | end if;
| |
| − | end for;
| |
| − |
| |
| − | end Opt_Gold;
| |
| − |
| |
| − | |}
| |
| − |
| |
| − | Setelah dilakukan stimulate akan didapatkan hasil optimum pada material yang tetap
| |
| − | [[File:TB_Metnum3.png|500px|center]]
| |
| − |
| |
| − | Maka dapat disimpulkan bahwa area yang paling optimum untuk material yang tetap adalah 0.000373673 mm^2
| |
| − |
| |
| − | Untuk area yang tetap, saya juga melakukan optimasi golden ratio
| |
| − |
| |
| − | {| class="wikitable"
| |
| − | |-
| |
| − | | style='border-style: none none solid solid;' |
| |
| − |
| |
| − | model Opt_Gold
| |
| − |
| |
| − | parameter Real xd[:]={193000000000,195000000000,197000000000};
| |
| − | parameter Real yd[size(xd,1)]={1109361,732878,531151};
| |
| − | parameter Real xlo=111e-193000000000;
| |
| − | parameter Real xhi=197000000000;
| |
| − | parameter Integer N=10; // maximum iteration
| |
| − | parameter Real es=0.0001; // maximum error
| |
| − |
| |
| − | Real f1[N], f2[N], x1[N], x2[N], ea[N], y[3];
| |
| − | Real xopt, fx;
| |
| − | protected
| |
| − | Real d, xl, xu, xint, R=(5^(1/2)-1)/2;
| |
| − |
| |
| − | algorithm
| |
| − | xl := xlo;
| |
| − | xu := xhi;
| |
| − | y := Curve_Fitting(xd,yd);
| |
| − |
| |
| − | for i in 1:N loop
| |
| − | d:= R*(xu-xl);
| |
| − | x1[i]:=xl+d;
| |
| − | x2[i]:=xu-d;
| |
| − | f1[i]:=y[1]*x1[i]^2+y[2]*x1[i]+y[3];
| |
| − | f2[i]:=y[1]*x2[i]^2+y[2]*x2[i]+y[3];
| |
| − | xint:=xu-xl;
| |
| − |
| |
| − | if f1[i]>f2[i] then
| |
| − | xl:=x2[i];
| |
| − | xopt:=x1[i];
| |
| − | fx:=f1[i];
| |
| − | else
| |
| − | xu:=x1[i];
| |
| − | xopt:=x2[i];
| |
| − | fx:=f2[i];
| |
| − | end if;
| |
| − |
| |
| − | ea[i]:=(1-R)*abs((xint)/xopt);
| |
| − | if ea[i]<es then
| |
| − | break;
| |
| − | end if;
| |
| − | end for;
| |
| − |
| |
| − | end Opt_Gold;
| |
| − |
| |
| − | |}
| |
| − |
| |
| − |
| |
| − | [[File:QQ-3.png|500px|center]]
| |
| − |
| |
| − |
| |
| − | Untuk area tetap, saya mendapatkan material dengan elastisitas 1.93e11 N/m^2, yang mana itu merupakan material ss304
| |
