Difference between revisions of "Raihan Tsaqif A"
(→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 55: | Line 56: | ||
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 79: | 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 96: | Line 99: | ||
-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 114: | Line 118: | ||
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 124: | Line 131: | ||
Real err=10e-10; | Real err=10e-10; | ||
Real ers=10e-4; | Real ers=10e-4; | ||
| + | |||
algorithm | algorithm | ||
//Creating Global Matrix | //Creating Global Matrix | ||
| Line 131: | Line 139: | ||
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 140: | Line 149: | ||
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 148: | Line 158: | ||
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 157: | Line 169: | ||
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 166: | Line 181: | ||
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:56, 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;
