Difference between revisions of "Kelompok 12"

From ccitonlinewiki
Jump to: navigation, search
(Tugas Finite Element Analysis dengan Eliminasi Gauss)
Line 24: Line 24:
 
== Tugas Finite Element Analysis dengan Eliminasi Gauss ==
 
== Tugas Finite Element Analysis dengan Eliminasi Gauss ==
  
 +
Dengan kode yang telah dibuat untuk menyelesaikan masalah pada sistem persamaan linear, kode tersebut juga dapat diaplikasikan pada simulasi Finite Element Analysis
 +
Contoh Soal
  
Soal Statika Sturktur dengan Element Analysis menggunakan Eliminasi Gauss
+
[[File:Soal Gauss Simulasi.png|700px]]
  
Pertimbangkan tiga elemen sistem di mana Node 1 dilampirkan ke fix support, menghasilkan kendala perpindahan U1 = 0, k1 = 50 lb / in, k2 = 75 lb / in, dan k3= 100 lb/in.  F1=0, F2=75, F3=50, F4=25 untuk kondisi ini tentukan perpindahan nodal U2, U3, dan U4
+
Seperti pada gambar diatas, suatu struktur diberikan gaya P keatas pada jarak 1 m dari lantai dan gaya F kebawah pada ujung struktur, hitunglah nodal displacement yang terjadi.
  
 +
Penyelesaian
  
import numpy as np
+
[[File:Element.png|500px]]
  
  A = np.array ([[50,-50, 0, 0],[-50, 125, -75, 0],[0, -75, 175, -100], [0, 0, -100, 100]])
+
* untuk menyelesaikan masalah ini, kita dapat membagikan struktur tersebut menjadi 3 element dengan tinggi 1 meter dengan luas penampang rata-rata berturut-turut dari lantai 950 m2, 850 m2, dan 750 m2
  b = np.array ([[0],[75],[50],[25]])
+
* dengan rumus stiffness k = A*E/L kita dapat menghitung k1, k2, k3 berturut-turut 9500 N/m, 8500 N/m, dan 7500 N/m
 +
* kemudian kita masukkan nilai-nilai stiffness '''k''' dan gaya '''F''' dan '''P''' beserta constraint yang ada ke matriks untuk penyelesaian
  
 
+
[[File:Matrix Simulasi.png|500px]]
def gaussElimin(A,b):
 
  n = len(b)
 
  for k in range (0,n-1):
 
  for i in range(k+1,n):
 
  if A[i,k] != 0.0:
 
  lam = A[i,k]/A[k,k]
 
  A[i,k+1:n] = A[i,k+1:n] - lam*A[k,k+1:n]
 
  b[i] = b[i] - lam*b[k]
 
  for k in range(n-1,-1,-1):
 
  b[k] = (b[k] - np.dot(A[k,k+1:n],b[k+1:n]))/A[k,k]
 
  return b
 
  x = gaussElimin(A,b)
 
  print (x)
 

Revision as of 02:00, 8 October 2019

Perkenalkan, kami kelompok 12 dengan anggota kelompok sebagai berikut:

1. Elvin

2. Khairun Naziri Batubara

3. Adam Ilham Maulana


Tugas Metode Eliminasi Gauss

Video diatas merupakan penjelasan program eliminasi Gauss yang dibuat pada python. Rangkaian kode dari video diatas ditampilkan sebagai berikut:

Code Gauss Elimination Kel.12.png

Console Gauss Elimination Kel.12.png


Tugas Finite Element Analysis dengan Eliminasi Gauss

Dengan kode yang telah dibuat untuk menyelesaikan masalah pada sistem persamaan linear, kode tersebut juga dapat diaplikasikan pada simulasi Finite Element Analysis Contoh Soal

Soal Gauss Simulasi.png

Seperti pada gambar diatas, suatu struktur diberikan gaya P keatas pada jarak 1 m dari lantai dan gaya F kebawah pada ujung struktur, hitunglah nodal displacement yang terjadi.

Penyelesaian

Element.png

  • untuk menyelesaikan masalah ini, kita dapat membagikan struktur tersebut menjadi 3 element dengan tinggi 1 meter dengan luas penampang rata-rata berturut-turut dari lantai 950 m2, 850 m2, dan 750 m2
  • dengan rumus stiffness k = A*E/L kita dapat menghitung k1, k2, k3 berturut-turut 9500 N/m, 8500 N/m, dan 7500 N/m
  • kemudian kita masukkan nilai-nilai stiffness k dan gaya F dan P beserta constraint yang ada ke matriks untuk penyelesaian

Matrix Simulasi.png