Gemi Hidrostatiği Eğriler Ödevi

Burada eğrilerin değerleri Python'da numpy kütüphanesiyle nümerik olarak hesaplanıyor. Bizi verilen seri 60 için offset tablosunu matris olarak içeriye aktarıyorum. Offset için oluşturduğum bu matrisin satırı her bir su hattı için yarı genişlik değerini veriyor. Bu benim yapacak nümerik hesapta dikey eksende yani z olarak düşüneceğim. Matrisin sütunu ise her posta için yarı genişlik değerini veriyor. Bu da benim yapacağım nümerik hesapta yatay eksende yani x olarak düşüneceğim. wl diye tanımladığım liste dikey yönünde su hattı konumunu veriyor. Bunu eğer z 'ye göre integral alırsam yani trapez yaparsam dz olarak kullanacağım veya eğer trapezin içinde gelmiyorsa dikeyde moment kolu olarakta kullanırım. pn diye tanımladığım liste ise yatay yönünde posta konumlarını veriyor. Bunu eğer x 'a göre integral alırsam yani trapez yaparsam dx olarak kullanacağım veya eğer trapezin içinde gelmiyorsa yatayda moment kolu olarak kullanırım. Sonuc tablosunu ise pandas kütüphanesiyle önce tablo dönüştürüyorum, sonra ise bunu dataframe_image kütüphanesiyle resim olarak çıktısını aldım. EKSİK OLAN EĞRİ İÇİN ISLAK YÜZEY ALAN EĞRİSİ def arc_length(x, y): npts = len(x) arc = np.sqrt((x[0] - x[1])**2 + (y[0] - y[1])**2) for k in range(1, npts): arc = arc + np.sqrt((x[k - 1] - x[k])**2 + (y[k - 1] - y[k])**2) return arc/2 l = np.zeros((row, col)) for i in range(row): for j in range(1, col): l[i, j] = round(arc_length(offset[i - 1, j + 1 : j], wl[j + 1 : j - 1]), 3) S = np.zeros(col) for i in range(1, col): S[i] = S[i-1] + 2*np.trapz(l[:, i], pn) Burada sadece değerleri verir videoyu izleyerek onuda çizime ekleyebilirsiniz. Hesaplar yapılırken kullanılan kaynak Prof. Dr. Fahri Çelik hocamızın avesis'teki ders notlarındandır. https://drive.google.com/file/d/1xTcT... Kaynak kodu indirmek için bu linki kullanabilirsiniz https://github.com/nediry/ship-hydros... https://github.com/nediry/ship-hydros...