Difference between revisions of "Raihan Tsaqif A"
(→TUGAS BESAR METODE NUMERIK) |
(→TUGAS BESAR METODE NUMERIK) |
||
| Line 81: | Line 81: | ||
14,15;//3rd floor | 14,15;//3rd floor | ||
15,16;//3rd floor | 15,16;//3rd floor | ||
| − | 13,16];//3rd floor | + | 13,16];//3rd floor |
| − | |||
//define coordinates (please put orderly) | //define coordinates (please put orderly) | ||
parameter Real P[Points,3]=[0.3,-0.375,0; //1 | parameter Real P[Points,3]=[0.3,-0.375,0; //1 | ||
| Line 99: | Line 98: | ||
-0.3,-0.375,1.8; //14 | -0.3,-0.375,1.8; //14 | ||
-0.3,0.375,1.8; //15 | -0.3,0.375,1.8; //15 | ||
| − | 0.3,0.375,1.8]; //16 | + | 0.3,0.375,1.8]; //16 |
| − | |||
//define external force (please put orderly) | //define external force (please put orderly) | ||
parameter Real F[Points*3]={0,0,0, | parameter Real F[Points*3]={0,0,0, | ||
| Line 118: | Line 116: | ||
0,0,-1000, | 0,0,-1000, | ||
0,0,-500}; | 0,0,-500}; | ||
| − | |||
//define boundary | //define boundary | ||
parameter Integer b[:]={1,2,3,4}; | parameter Integer b[:]={1,2,3,4}; | ||
| − | |||
//solution | //solution | ||
Real displacement[N], reaction[N]; | Real displacement[N], reaction[N]; | ||
Real check[3]; | Real check[3]; | ||
| − | |||
parameter Integer N=3*Points; | parameter Integer N=3*Points; | ||
Integer boundary[3*size(b,1)]=cat(1,(3*b).-2,(3*b).-1,3*b); | Integer boundary[3*size(b,1)]=cat(1,(3*b).-2,(3*b).-1,3*b); | ||
| Line 131: | Line 126: | ||
Real err=10e-10; | Real err=10e-10; | ||
Real ers=10e-4; | Real ers=10e-4; | ||
| − | |||
algorithm | algorithm | ||
//Creating Global Matrix | //Creating Global Matrix | ||
Revision as of 14:52, 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
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
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;
