Bugün Sistem Analizi dersinin final sınavı vardı. Final soruları aşağıda:

Tüm sorular bunlar. Aslında sayısal işlemin olmadığı kolay bir sınavdı. Bakalım sonuçlar ne olacak.

Dün Matlab üzerinde eşanjör projelendirmesine başladım. Programı yazmam yaklaşık 1.5 günümü aldı. Şuan için akış yönü farketmeden, her türlü akışkan için Kd, hi ve hd değerlerini buldurabilmek mümkün. Kd’yi başlangıçta rastgele atan program, iterasyon yöntemi ile Kd’yi hi ve hd ile arasındaki bağıntıya eşitleyerek, size doğru Kd değerini ve bu değere bağlı olarak çıkan boru sayısı ve h değerlerini veriyor. Programın başından bir kesit:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21

clc
clear
disp (' ');
disp (' ');
disp ('ESANJOR PROJELENDIRME PROGRAMI v1.1');
disp (' ');
disp (' ');
disp ('Programlama: Yusuf ÖZTÜRK');
disp ('Yıldız Teknik Üniversitesi Makina Mühendesliği Bölümü');
disp ('Termodinamik ve Isı Tekniği Anabilim Dalı');
disp (' ');
disp (' ');
disp ('Hesaplamaya baslamak icin lutfen bir tusa basiniz.. ');
pause
clc
Tog = input ('Soguk akiskan giris sicakligini (Tog) giriniz: ');
Toc = input ('Soguk akiskan cikis sicakligini (Toc) giriniz: ');
Tsg = input ('Sicak akiskan giris sicakligini (Tsg) giriniz: ');
Tsc = input ('Sicak akiskan cikis sicakligini (Tsc) giriniz: ');
sogukakiskankonum = input ('Soguk akiskan boru icinden gidecekse 1, boru disindan gidecekse 2: ');
Debi = input ('Sicak akiskan debisi (Ms) icin 1, Soguk akiskan debisi (Mo) icin 2: ');

İlerleyen akademik hayatımda kullanabilmem amacıyla kodlarımı paylaşmıyorum fakat Matlab üzerinde biraz zaman geçirerek, güzel bir eşanjör projelendirme programı yapabilirsiniz. Cidden matlab ile yapabileceklerinizin sınırı yok. Kolay gelsin.

Sistem analizi dersinde Tamer Kepçeler hocamızın verdiği değerler doğrultusunda, dişli çark üzerinde hacim hesabı yapılmıştır. Bulunan hacim, sınır değerlerini aşmayacak şekilde optimize edilerek ortaya koyulmaktadır. Matlab komutlarını aşağıda veriyorum.

Öncelikle ilk function dosyamızı yazalım.

1
2

function f=GearObjFunct(x)
f=(x(6)*(x(3)+x(1)*x(3))/2+x(7)*(x(4)+x(2)*x(4))/2+40)*5*x(5)*(x(2)*x(4)*x(7)+50);

Yukardaki function dosyasını GearObjFunct.m olarak kaydedelim. Diğer function dosyamız,

1
2
3
4
5
6
7

function [C,Ceq]=GearNonLinConstr(x)
C(1)=(400/(x(3)*x(5)))-x(6)^2;
C(2)=(400*x(1)/(x(4)*x(5)))-x(7)^7;
C(3)=40-x(3)*x(6);
C(4)=60-x(4)*x(7);
C(5)=12-x(1)*x(2);
Ceq=[];

Yukardaki function dosyasını ise GearNonLinConstr.m olarak kaydedelim. Şimdi gelelim hacmi verecek olan komutlara,

1
2
3
4
5
6

x0=[2,1.5,14,17,20,0.5,1];
LBnd=[2,1.5,14,17,20,0.5,1];
UBnd=[4,5,20,23,60,2.5,3];
A=zeros(4,7);
b=[0,0,0,0];
[x,f]=fmincon('GearObjFunct',x0,A,b,[],[],LBnd,UBnd,'GearNonLinConstr')

Peki X olarak verdiğimiz girdiler neler? Hemen onları vereyim.

GearObjFunct içersinde, dişli kutusunun hacmini verecek olan bağıntıyı yazdık. İkinci function dosyasında ise parametreler arasındaki bağıntıları belirttik. Belirttik çünkü parametreler birbirleriyle bağıntılı olduğu için, bir parametrenin değişmesi, diğer parametrelerinden de değişmesini sağlamalı. Yoksa herhangi bir şey belirtmemiş olsak, optimizasyon işlemi sırasında tek parametre değiştirilir, diğerleri ise sabit kalırdı. İşlemde ise LBnd ve UBnd ile X[Y] değerlerinin olması gereken sınırları belirttik. Örneğin X[1] değerinin 2 ile 4 arasında olmasını isteriz. Böylece optimizasyon da bu sınırlar arasında kalmış oldu. İyi çalışmalar..

Sistem Analizi dersinde işlenilen çelik çatı konstrüksiyonlarında kullanılan profiller üzerindeki kesitlerin atalet ve mukavemet momentlerini arıyorsanız doğru yerdesiniz. Aşağıdaki döküman ile bir I profil ya da U profil üzerinde nasıl eğilme momenti vs. bulabileceğinizi görebilirsiniz. Bu konu da Tamer Kepçeler hocamın talep ettiği bir ödevdir ve üzerinde son derece önemle durulması gereken bir konudur. Kolay gelsin.

Dökümanı hazırlayarak, biz öğrencilerin kullanımına sunan Uludağ Üniversitesine teşekkürler..

Ansys üzerinde akışkanlar mekaniği ile ilgili analizler yaptırabilirsiniz. Flotran ile bu dökümantasyonda 2 boyutlu dış akış problemleri ve 2 boyutlu laminar akış problemlerine giriş yapacağız. Akışkanlar mekaniği temeliniz iyi değilse, acilen İsmail Teke hocamızın ders notlarını internette aramanızı tavsiye ederim.

Flotran modülü dökümantasyonunu aşağıdaki bağlantıdan indirebilirsiniz:

PDF dökümanını hazırladıkları için Dokuz Eylül Üniversitesine tekrar teşekkürler..

Daha önceki sistem analiziyle ilgili yazılarımda, Ansys dökümantasyonları yayınlama başlamıştım. Önceki yazıları okuyanlar bilecektir, Tamer Kepçeler hocamızın tavsiyesi üzerine Ansys programını kullanmaya başlamıştım. Ansys programı ile bir çok sistemin analizi çok kolay bir şekilde yapılabiliyordu. Şimdiki eğitim ise Ansys Workbench üzerine. Bu dökümantasyonu hazırlayan Dokuz Eylül Üniversitesi’ne teşekkürlerimi iletiyorum.

PDF formatındaki bu dökümanı, aşağıdaki bağlantıdan indirebilirsiniz:

Tamer Kepçeler hocamız ile uzay çatı sistemlerine giriş yapmıştık. Kuleli Askeri Lisesi’ne Uzay çatı sistemi kurulurken az çok bilgilenmiştim bu alanda. Uzay çatıların hesaplanması oldukça zor, tabi Ansys gibi bir program kullanmıyorsanız. Bu arada belirtmek istiyorum, uzay çatı dediğimiz şey, çok geniş alanların herhangi bir kolon kullanılmadan kapatılmasıdır. Havalanları, alışveriş merkezleri gibi mekanlarda uzay çatı sistemleri kullanılmaktadır.

Yukardaki dökümanda Ansys programı ile uygulamalı bir uzay çatı analizi bulabilirsiniz. Bu dökümanı hazırladıkları için Firat Üniversitesi’nden Aydın Turgut ve Yavuz Solmaz hocalarımıza teşekkürlerimi iletiyorum.