Difference between revisions of "Kelompok 4"
(→TUGAS 2: Gauss Jordan dalam Penyelesaian Statika Struktur) |
(→TUGAS 2: Gauss Jordan dalam Penyelesaian Statika Struktur) |
||
Line 30: | Line 30: | ||
Pada tugas ini kami menyelesaikan suatu persamaan Statika Struktur dengan metode Gauss Jordan menggunakan python. | Pada tugas ini kami menyelesaikan suatu persamaan Statika Struktur dengan metode Gauss Jordan menggunakan python. | ||
− | Kami mengambil contoh soal | + | Kami mengambil contoh soal sebagai berikut: |
− | |||
+ | Sebuah batang panjang 12 m diberi beban 1500N ditengah batang. | ||
+ | Lalu kami menghitung support reaction yang terjadi pada dua ujung batang. | ||
import numpy as np | import numpy as np | ||
Line 53: | Line 54: | ||
a[i][f]=a[i][f]-op*a[i-1][f] | a[i][f]=a[i][f]-op*a[i-1][f] | ||
− | |||
if a[n-1][n-1]!=0: | if a[n-1][n-1]!=0: | ||
c[n-1]=b[n-1][0]/a[n-1][n-1] | c[n-1]=b[n-1][0]/a[n-1][n-1] |
Revision as of 21:44, 15 October 2019
TUGAS KELOMPOK 4
Bismillahirrahmanirrahim.
PERKENALAN
Perkenalkan kami dari kelompok 4 yang beranggotakan:
1. Satria Nur Muhammad
2. Stefanus Nandiwardhana
3. Rafiq Ali Abdillah
TUGAS 1: Gauss Jordan dalam Penyelesaian Persamaan Aljabar Linear
Pada tugas ini kami menyelesaikan suatu persamaan Aljabar Linear dengan metode Gauss Jordan menggunakan python.
TUGAS 2: Gauss Jordan dalam Penyelesaian Statika Struktur
Pada tugas ini kami menyelesaikan suatu persamaan Statika Struktur dengan metode Gauss Jordan menggunakan python.
Kami mengambil contoh soal sebagai berikut:
Sebuah batang panjang 12 m diberi beban 1500N ditengah batang. Lalu kami menghitung support reaction yang terjadi pada dua ujung batang.
import numpy as np a=np.array ([[1,1], [7,-5]]) b=np.array ([[1500], [0]]) c=[1,1,1,1,1,1,1,1] n=len(b) for k in range(0,n-1): for i in range(n-1,k,-1): if a[i][k]!=0.0: op=a[i][k]/a[i-1][k] b[i][0]=b[i][0]-op*b[i-1][0] for f in range (0,n): a[i][f]=a[i][f]-op*a[i-1][f] if a[n-1][n-1]!=0: c[n-1]=b[n-1][0]/a[n-1][n-1] else: c[n-1]=0 for i in range(n-2,-1,-1): sigma=0 for k in range(0,n,1): y=a[i][k]*c[k] sigma=sigma+y sigma=sigma-a[i][i] if a[i][i]!=0: c[i]=(b[i][0]-sigma)/a[i][i] else: a[i][i]=0 for v in range(0,n,1): print(a[v]) print(c[0:n])