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

Tamer Hoca’nın derslerini düzenli olarak takip ediyorsanız, zorlanmadan yapacağınız bir sınavdı. Notlar serbest olduğu için stres yapmanıza gerek yok. Aldığım notu da buraya yazmayı düşünüyorum :)

Bu arada klinker üretimi ile birlikte çimento üretimini de anlatmış, çok güzel bir animasyon daha buldum.

Animasyonun üçüncü sayfasında klinker üretimini görebilirsiniz.

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.

Böylece istediğiniz moment ve hıza sahip servo motorun ne kadarlık bir güce ihtiyacı olduğunu bulabilirsiniz. Bu arada unutmadan, az bir güçle çok bir güç elde edilmez :)

Ö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..

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

Muhtemelen pek bir şey anlamadınız yukardaki işlemlerden. Bu yüzden projenin tamamına göz atsanız daha iyi olacak. Autocad’de yapılmış bir proje ile sadece hesapsal olarak ayarlanmış bir proje ekliyorum aşağıya. Biraz düşünerek okuduğunuzda, projenin nasıl yapıldığını az çok anlayacaksınız. Kolay gelsin.

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..

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