Difference between revisions of "Raihan Tsaqif A"

From ccitonlinewiki
Jump to: navigation, search
(TUGAS BESAR METODE NUMERIK)
(TUGAS BESAR METODE NUMERIK)
Line 50: Line 50:
  
 
model Trusses_3D_Tugas_Besar_Simplified2
 
model Trusses_3D_Tugas_Besar_Simplified2
 
 
//define initial variable
 
//define initial variable
 
parameter Integer Points=16; //Number of Points
 
parameter Integer Points=16; //Number of Points
Line 56: Line 55:
 
parameter Real Area=3777777.778; //Area
 
parameter Real Area=3777777.778; //Area
 
parameter Real Elas=1; //Elasticity (equals to one in order to determine the displacement limit)
 
parameter Real Elas=1; //Elasticity (equals to one in order to determine the displacement limit)
 
 
//define connection
 
//define connection
 
parameter Integer C[Trusses,2]=[1,5;  
 
parameter Integer C[Trusses,2]=[1,5;  
Line 133: Line 131:
 
   q1[j]:=P[C[i,1],j];
 
   q1[j]:=P[C[i,1],j];
 
   q2[j]:=P[C[i,2],j];
 
   q2[j]:=P[C[i,2],j];
end for;
+
end for;    
     
 
 
   //Solving Matrix
 
   //Solving Matrix
 
   L:=Modelica.Math.Vectors.length(q2-q1);
 
   L:=Modelica.Math.Vectors.length(q2-q1);
Line 143: Line 140:
 
                     cy*cx,cy^2,cy*cz;
 
                     cy*cx,cy^2,cy*cz;
 
                     cz*cx,cz*cy,cz^2];
 
                     cz*cx,cz*cy,cz^2];
 
 
   //Transforming to global matrix
 
   //Transforming to global matrix
 
   g:=zeros(N,N);  
 
   g:=zeros(N,N);  
Line 152: Line 148:
 
     g[3*(C[i,1]-1)+m,3*(C[i,2]-1)+n]:=-X[m,n];
 
     g[3*(C[i,1]-1)+m,3*(C[i,2]-1)+n]:=-X[m,n];
 
   end for;   
 
   end for;   
 
 
G_star:=G+g;
 
G_star:=G+g;
 
G:=G_star;
 
G:=G_star;
 
end for;
 
end for;
 
 
//Implementing boundary
 
//Implementing boundary
 
for i in boundary loop
 
for i in boundary loop
Line 163: Line 157:
 
end for;
 
end for;
 
end for;
 
end for;
 
 
//Solving displacement
 
//Solving displacement
 
displacement:=Modelica.Math.Matrices.solve(G,F);
 
displacement:=Modelica.Math.Matrices.solve(G,F);
 
 
//Solving reaction
 
//Solving reaction
 
reaction:=(G_star*displacement)-F;
 
reaction:=(G_star*displacement)-F;
 
 
//Eliminating float error
 
//Eliminating float error
 
for i in 1:N loop
 
for i in 1:N loop
Line 175: Line 166:
 
  displacement[i]:=if abs(displacement[i])<=err then 0 else displacement[i];
 
  displacement[i]:=if abs(displacement[i])<=err then 0 else displacement[i];
 
end for;  
 
end for;  
 
 
//Checking Force
 
//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[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[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)});
 
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
 
for i in 1:3 loop
 
  check[i] := if abs(check[i])<=ers then 0 else check[i];
 
  check[i] := if abs(check[i])<=ers then 0 else check[i];
 
end for;
 
end for;
 
 
end Trusses_3D_Tugas_Besar_Simplified2;
 
end Trusses_3D_Tugas_Besar_Simplified2;

Revision as of 14:53, 6 January 2021

Biodata

Nama : Raihan Tsaqif A

NPM : 1906379176

TTL : Semarang, 8 Oktober 2001

Hobi : Belajar

Saya adalah mahasiswa teknik mesin angkatan 2019, saya tertarik mengambil jurusan teknik mesin karena saya melihat bahwa di seluruh dunia ini semua bekerja dengan mesin. Mesin tidak hanya sesuatu yang melakukan pembakaran akan tetapi sebuah mekanisme pada suatu benda juga bisa disebut mesin sehingga dapat disimpulkan bahwa jurusan ini memiliki cakupan ilmu teknik yang sangat luas. Tujuan saya dalam mengikuti mata kuliah ini adalah untuk menjadi insan yang berguna untuk lingkungan sekitar dan dapat mendobrak perkembangan teknologi di Indonesia.

Tugas Metode Numerik

Pada tugas kali ini saya menggunakan aplikasi modelica dalam menentukan perubahan momentum tiap satuan waktu atau sebutan lainnya adalah gaya. Oleh karena itu saya lampirkan video sebagai berikut, sehingga harapannya dapat memberikan penjelasan juga kepada yang lainnya.



QUIZ METNUM

Flowchart

Pada quiz kali ini saya akan melampirkan langkah pengerjaan melalui flowchart


Flowchart tsaqif.jpg


TUGAS BESAR METODE NUMERIK

Dalam tugas besar ini kami diminta untuk mendesain rangka sesuai yang tertera dibawah dengan menggunakan rangka yang optimal dan cost yang minimum

Soaltubes2.jpg

Solaaa.jpg

Adapun Hal yang Harus Diperhitungkan sebagai Plotting

a. Harga material di dapat dari nilai optimum defleksi dan cross section

b. Material (Elastisitas properti)

c. Area Cross Section Truss (L profile/truss siku), luas cross section sebagai X

d. Defleksi sebagai Y

Trusses Modelling

model Trusses_3D_Tugas_Besar_Simplified2 //define initial variable parameter Integer Points=16; //Number of Points parameter Integer Trusses=24; //Number of Trusses parameter Real Area=3777777.778; //Area parameter Real Elas=1; //Elasticity (equals to one in order to determine the displacement limit) //define connection parameter Integer C[Trusses,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[Points,3]=[0.3,-0.375,0; //1

                           -0.3,-0.375,0;    //2
                           -0.3,0.375,0;     //3
                           0.3,0.375,0;      //4
                           0.3,-0.375,0.6;   //5
                           -0.3,-0.375,0.6;  //6
                           -0.3,0.375,0.6;   //7
                           0.3,0.375,0.6;    //8
                           0.3,-0.375,1.2;   //9
                           -0.3,-0.375,1.2;  //10  
                           -0.3,0.375,1.2;   //11
                           0.3,0.375,1.2;    //12
                           0.3,-0.375,1.8;   //13
                           -0.3,-0.375,1.8;  //14
                           -0.3,0.375,1.8;   //15
                           0.3,0.375,1.8];   //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};

//define boundary parameter Integer b[:]={1,2,3,4}; //solution Real displacement[N], reaction[N]; Real check[3]; parameter Integer N=3*Points; Integer boundary[3*size(b,1)]=cat(1,(3*b).-2,(3*b).-1,3*b); 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; Real 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 i in boundary loop for j in 1:N loop

 G[i,j]:=id[i,j];

end for; 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; end Trusses_3D_Tugas_Besar_Simplified2;