6-Amaliy mashg’ulot Mavzu: Diffuziyali modelni Ishdan maqsad

Ishni bajarish buyicha metodik ko`rsatma.

Bu ishimizda diffuziyali model bu - mahsulot bo’lib quyidagi parametrlar berilgan bolib, quyidagilar berilgan:

u=0.001m/s – mahsulotning o’rtacha tezligi
p=1200 kg/m3 – mahsulot zichligi
cm=800 Dj/(kg.grad) – mahsulotning o’rtacha issiqlik sig’imi
l=0,2 m – issiqlik muhitining uzunligi
tbx=315 K – mahsulotning boshlang’ich temperaturasi
tc=415 K – mahsulotning chiqishdagi temperaturasi
d1= 0,2 m – silindrning tashqi diametri
d2=0,17 m – silidrning ichki diametri
n=10 – elementar zonalar soni
t=300 sek – vaqt.
Berilganlardan foydalanib mahsulot temperaturasing vaqt bo’yicha o’zgarish sxemasini komp’yuterda qurish talab etiladi.


Ishni bajarish tartibi.

1. 
   MatLAB dasturi ishga tushiriladi.
   
2. 
   File-function for dextrud nomli fayl yaratilib quyidagilar kiritiladi
   

function dtm=dextrud(t,tm,u,l,rom,cm,d1,d2,kt,tbx,tc);
n=length(tm);
k1=u*n/l;
k2=kt*d1*l/((d1^2-d2^2)*rom*cm*n);
dtm(1,1)=k1*tbx-(k1+k2)*tm(1,1)+k2*tc;
for i=2:n
dtm(i,1)=k1*tm(i-1)-(k1+k2)*tm(i,1)+k2*tc;
end;

3. 
   Skript-file for dextrud nomli ikkinchi fayl yaratilib, berilgan parametr qiymatlari quyidagi ketma-ketlikda kiritiladi
   

u=0.001;n=10;rom=1200;cm=800;l=0.2;
d1=0.2;d2=0.17;kt=650*n/l;tbx=315;tc=415;
tfin=300;
tspan=[0:tfin];
for i=1:n
tm0(i)=315;
end;
[t,tm]=ode15s(@dextrud,tspan,tm0,[],u,l,rom,cm,d1,d2,kt,tbx,tc);

4. 
   Temperaturaning vaqt bo’yicha o’zgarish grafigini chizish uchun quyidagilar kiritiladi
   

plot(t,tm);grid on;
xlabel('Vaqt');ylabel('Temperatura');

5. 
   Temperaturaning vaqt bo’yicha o’zgarish grafigini bosmaga chiqariladi.
   

6. 
   Grafikdagi 10 ta zona bo’yicha olingan natijalar jadvalga yoziladi.