Difference between revisions of "Calvin rahmat pratama"

From ccitonlinewiki
Jump to: navigation, search
(Tugas Minggu 3)
Line 53: Line 53:
  
 
==Tugas Minggu 3 ==
 
==Tugas Minggu 3 ==
model Trusses
+
model Trusses
parameter Integer N=10; //Global matrice = 2*points connected
+
parameter Integer N=10; //Global matrice = 2*points connected
parameter Real A=8;
+
parameter Real A=8;
parameter Real E=1.9e6;
+
parameter Real E=1.9e6;
Real G[N,N]; //global
+
Real G[N,N]; //global
Real Ginitial[N,N]; //global
+
Real Ginitial[N,N]; //global
Real Sol[N]; //global dispplacement
+
Real Sol[N]; //global dispplacement
Real X[N]={0,0,0,0,0,0,0,-500,0,-500};
+
Real X[N]={0,0,0,0,0,0,0,-500,0,-500};
Real R[N]; //global reaction force
+
Real R[N]; //global reaction force
Real SolMat[N,1];
+
Real SolMat[N,1];
Real XMat[N,1];
+
Real XMat[N,1];
  
 
//boundary coundition
 
//boundary coundition
Integer b1=1;
+
Integer b1=1;
Integer b2=3;
+
Integer b2=3;
  
 
//truss 1
 
//truss 1
parameter Real X1=0; //degree between truss
+
parameter Real X1=0; //degree between truss
Real k1=A*E/36;
+
Real k1=A*E/36;
Real K1[4,4]; //stiffness matrice
+
Real K1[4,4]; //stiffness matrice
Integer p1a=1;
+
Integer p1a=1;
Integer p1b=2;
+
Integer p1b=2;
Real G1[N,N];
+
Real G1[N,N];
  
 
//truss 2
 
//truss 2
parameter Real X2=135; //degree between truss
+
parameter Real X2=135; //degree between truss
Real k2=A*E/50.912;
+
Real k2=A*E/50.912;
Real K2[4,4]; //stiffness matrice
+
Real K2[4,4]; //stiffness matrice
Integer p2a=2;
+
Integer p2a=2;
Integer p2b=3;
+
Integer p2b=3;
Real G2[N,N];
+
Real G2[N,N];
  
 
//truss 3
 
//truss 3
parameter Real X3=0; //degree between truss
+
parameter Real X3=0; //degree between truss
Real k3=A*E/36;
+
Real k3=A*E/36;
Real K3[4,4]; //stiffness matrice
+
Real K3[4,4]; //stiffness matrice
Integer p3a=3;
+
Integer p3a=3;
Integer p3b=4;
+
Integer p3b=4;
Real G3[N,N];
+
Real G3[N,N];
  
 
//truss 4
 
//truss 4
parameter Real X4=90; //degree between truss
+
parameter Real X4=90; //degree between truss
Real k4=A*E/36;
+
Real k4=A*E/36;
Real K4[4,4]; //stiffness matrice
+
Real K4[4,4]; //stiffness matrice
Integer p4a=2;
+
Integer p4a=2;
Integer p4b=4;
+
Integer p4b=4;
Real G4[N,N];
+
Real G4[N,N];
 
+
 
//truss 5
 
//truss 5
parameter Real X5=45; //degree between truss
+
parameter Real X5=45; //degree between truss
Real k5=A*E/50.912;
+
Real k5=A*E/50.912;
Real K5[4,4]; //stiffness matrice
+
Real K5[4,4]; //stiffness matrice
Integer p5a=2;
+
Integer p5a=2;
Integer p5b=5;
+
Integer p5b=5;
Real G5[N,N];
+
Real G5[N,N];
  
 
//truss 6
 
//truss 6
parameter Real X6=0; //degree between truss
+
parameter Real X6=0; //degree between truss
Real k6=A*E/36;
+
Real k6=A*E/36;
Real K6[4,4]; //stiffness matrice
+
Real K6[4,4]; //stiffness matrice
Integer p6a=4;
+
Integer p6a=4;
Integer p6b=5;
+
Integer p6b=5;
Real G6[N,N];
+
Real G6[N,N];
  
 
/*
 
/*
Line 124: Line 124:
  
 
//creating global matrice
 
//creating global matrice
K1:=Matrices(X1);
+
K1:=Matrices(X1);
G1:=k1*GlobalMatrices(K1,N,p1a,p1b);
+
G1:=k1*GlobalMatrices(K1,N,p1a,p1b);
  
K2:=Matrices(X2);
+
K2:=Matrices(X2);
G2:=k2*GlobalMatrices(K2,N,p2a,p2b);
+
G2:=k2*GlobalMatrices(K2,N,p2a,p2b);
  
K3:=Matrices(X3);
+
K3:=Matrices(X3);
G3:=k3*GlobalMatrices(K3,N,p3a,p3b);
+
G3:=k3*GlobalMatrices(K3,N,p3a,p3b);
  
K4:=Matrices(X4);
+
K4:=Matrices(X4);
G4:=k4*GlobalMatrices(K4,N,p4a,p4b);
+
G4:=k4*GlobalMatrices(K4,N,p4a,p4b);
  
K5:=Matrices(X5);
+
K5:=Matrices(X5);
G5:=k5*GlobalMatrices(K5,N,p5a,p5b);
+
G5:=k5*GlobalMatrices(K5,N,p5a,p5b);
  
K6:=Matrices(X6);
+
K6:=Matrices(X6);
G6:=k6*GlobalMatrices(K6,N,p6a,p6b);
+
G6:=k6*GlobalMatrices(K6,N,p6a,p6b);
  
G:=G1+G2+G3+G4+G5+G6;
+
G:=G1+G2+G3+G4+G5+G6;
Ginitial:=G;
+
Ginitial:=G;
  
 
//implementing boundary condition
 
//implementing boundary condition
Line 153: Line 153:
 
end for;
 
end for;
  
G[2*b1-1,2*b1-1]:=1;
+
G[2*b1-1,2*b1-1]:=1;
G[2*b1,2*b1]:=1;
+
G[2*b1,2*b1]:=1;
G[2*b2-1,2*b2-1]:=1;
+
G[2*b2-1,2*b2-1]:=1;
G[2*b2,2*b2]:=1;
+
G[2*b2,2*b2]:=1;
  
 
//solving displacement
 
//solving displacement
Sol:=Gauss_Jordan(N,G,X);
+
Sol:=Gauss_Jordan(N,G,X);
  
 
//solving reaction force
 
//solving reaction force
SolMat:=matrix(Sol);
+
SolMat:=matrix(Sol);
XMat:=matrix(X);
+
XMat:=matrix(X);
R:=Reaction_Matrices(N,Ginitial,SolMat,XMat);
+
R:=Reaction_Matrices(N,Ginitial,SolMat,XMat);
  
 
end Trusses;
 
end Trusses;

Revision as of 12:40, 2 December 2020

Biodata Diri

Nama: Calvin Rahmat Pratama

NPM: 1806200974

Tempat,tanggal lahir: Batam,14 Oktober 2000

Hobi: Gaming dan Golf


Saya adalah mahasiswa FTUI angkatan 2018 dari jurusan Teknik mesin dan saya adalah salah satu ciptaan terbaik dari Tuhan yang Maha Esa karena pada prinsipnya Tuhan yang Maha Esa itu mendesain manusia dengan sebaik baiknya makhluk.

Teknik Mesin merupakan program studi yang saya gemari dikarenakan ayah saya juga merupakan lulusan Teknik Mesin.

Tujuan Pembelajaran

1. Memahami konsep dengan dengan baik tentang dasar-dasar metode numerik.

2. Mampu menerapkan pemahaman konsep tersebut didalam permodelan numerik.

3. Mampu menerapkan metode numerik dipersoalan keteknikan.

4. Sebagai nilai tambah bagi saya agar menjadi manusia yang beradab seperti pada Sila ke-2 Pancasila.

Metode Numerika

Sejauh ini yang dapat saya ketahui dalam metode numerika adalah metode numerika dapat menyelesaikan permasalahan yang tidak dapat diselesaikan oleh metoda analisis saya juga mempelajari cara menggunakan excel untuk menyelesaikan beberapa permasalahan metode numerika. Seperti menggunakan psuedocode dan metode newton rhapson.

Tugas Minggu 1

Review Minggu 1

Untuk pertemuan pertama saya mempelajari cara menggunakan OpenModelica dimulai dari persamaan yang cukup sederhana, OpenModelica sendiri merupakan aplikasi penghitungan permodelan yang lumayan kompleks, sehingga OpenModelica itu sendiri sangat membantu saat menyelesaikan permodelan masalah mulai dari yang cukup mudah sampai ke cukup kompleks

Review Minggu 2

Alasan mengapa menggunakan openmodelica dibanding aplikasi lainnya

- Secara penggunanaan lebih digunakan untuk permodelan bukan untuk programming

- Didalam melakukan perhitungan kode yang diinput di modelica itu di terjemahkan ke Bahasa C dan kemudian diterjemahkan lagi ke Bahasa permesinan sehingga bisa melakukan perhitungan permodelan

- Paling penting adalah openmodelica itu gratis

Pada pertemuan kedua saya juga mempelajari menggunakan modelica untuk "memanggil" fungsi. Pada pertemuan kedua juga mempelajari mengenai Fungsi tambah yang mana terlebih dahulu dari permodelan masalah di translasi kan terlebih dahulu ke dalam bahasa C++ yang nantinya akan di proses menjadi data-data perhitungan.

Tugas Minggu 2

Review Minggu 3

Pada minggu ke tiga ini saya mencoba untuk mempelajari kembali untuk membuat fungsi modelica 9x12 dengan metode gauss atau menggunakan metode lain, kemudian mencoba mengimplementasikan psuedocode pada Figure 9.4 di modelica yang mana nantinya akan di test coding dengan example 9.5.

Tugas Minggu 3

model Trusses
parameter Integer N=10; //Global matrice = 2*points connected
parameter Real A=8;
parameter Real E=1.9e6;
Real G[N,N]; //global
Real Ginitial[N,N]; //global
Real Sol[N]; //global dispplacement
Real X[N]={0,0,0,0,0,0,0,-500,0,-500};
Real R[N]; //global reaction force
Real SolMat[N,1];
Real XMat[N,1];

//boundary coundition

Integer b1=1;
Integer b2=3;

//truss 1

parameter Real X1=0; //degree between truss
Real k1=A*E/36;
Real K1[4,4]; //stiffness matrice
Integer p1a=1;
Integer p1b=2;
Real G1[N,N];

//truss 2

parameter Real X2=135; //degree between truss
Real k2=A*E/50.912;
Real K2[4,4]; //stiffness matrice
Integer p2a=2;
Integer p2b=3;
Real G2[N,N];

//truss 3

parameter Real X3=0; //degree between truss
Real k3=A*E/36;
Real K3[4,4]; //stiffness matrice
Integer p3a=3;
Integer p3b=4;
Real G3[N,N];

//truss 4

parameter Real X4=90; //degree between truss
Real k4=A*E/36;
Real K4[4,4]; //stiffness matrice
Integer p4a=2;
Integer p4b=4;
Real G4[N,N];

//truss 5

parameter Real X5=45; //degree between truss
Real k5=A*E/50.912;
Real K5[4,4]; //stiffness matrice
Integer p5a=2;
Integer p5b=5;
Real G5[N,N];

//truss 6

parameter Real X6=0; //degree between truss
Real k6=A*E/36;
Real K6[4,4]; //stiffness matrice
Integer p6a=4;
Integer p6b=5;
Real G6[N,N];

/* for each truss, please ensure pXa is lower then pXb (X represents truss element number)

  • /

algorithm

//creating global matrice

K1:=Matrices(X1);
G1:=k1*GlobalMatrices(K1,N,p1a,p1b);
K2:=Matrices(X2);
G2:=k2*GlobalMatrices(K2,N,p2a,p2b);
K3:=Matrices(X3);
G3:=k3*GlobalMatrices(K3,N,p3a,p3b);
K4:=Matrices(X4);
G4:=k4*GlobalMatrices(K4,N,p4a,p4b);
K5:=Matrices(X5);
G5:=k5*GlobalMatrices(K5,N,p5a,p5b);
K6:=Matrices(X6);
G6:=k6*GlobalMatrices(K6,N,p6a,p6b);
G:=G1+G2+G3+G4+G5+G6;
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_Matrices(N,Ginitial,SolMat,XMat);

end Trusses;