qsv1z1.m

画一种新型气固分离器气相流场的matlab程序

clear;x1=-150;y1=-300:0.1:0;%颗粒出口左侧直边
t2=pi*(-50:1:50)/50;x2=150*cos(t2);y2=abs(150*sin(t2));%上部半圆形壳体
x3=150;y3=y1;%入口右侧直边
x4=-150:0.1:-90;y4=-300;%颗粒出口横边
x5=-90:80;y5=-100;%排气管下横边
x6=80;y6=-300:0.1:-100;%入口左侧直边
x7=-90;y7=y6;%颗粒出口右侧直边
x8=80:150;y8=y4;%入口横边
x9=80*cos(t2);y9=80*sin(t2);%排气管
plot(x1,y1,'b',x2,y2,'b',x3,y3,'b',x4,y4,'b',x5,y5,'b',x6,y6,'b',x7,y7,'b',x8,y8,'b',x9,y9,'b','LineWidth',1),axis image,axis off,grid off,box off,hold on;
vx=[-3.6,-3.5,-3.5,-3.4,-3.3,0.9,4.1,-3.5,4.1,-5.8,-7.1,-5.2,-7.0,-10.1,-13.6,-17.2,-16.9,-8.8,-10.0,-12.9,-13.8,-13.9,-12.6,-15.7,-3.1,-8.2,-6.2,-3.2,-1.4,-1.2,0.3,1.5,1.5,2.7,4.5,6.4,8.8,11.3,14.7];
vy=[10.2,10.2,11.1,10.7,10.3,4.3,11.3,9.2,11.4,15.5,14.3,12.2,12.9,12.3,10.4,7.5,0.6,7.7,7.5,5.3,1.9,-1.3,-3.6,-10.2,-12.0,-3.2,-6.2,-7.1,-6.1,-6.2,-6.2,-5.2,-4.3,-3.7,-0.4,0.5,0.9,2.7,6.3];
x=[140,127.3,114.5,101.8,89.1,76.4,63.6,140,123.3,106.7,90,140,122.6,104,84.3,63.6,41.8,121.2,100.,77.5,53.6,29.7,5.4,-18.2,-41.8,0,-37.2,-59.2,-73.8,-83.8,-88.5,-90,-140,-116.7,-93.3,-70,-46.7,-23.4,0];
y=[-230.4,-157.3,-128.1,-106.8,-89.8,-75.6,-63.6,-100,-36.6,-14.6,0,0,18.4,35.7,51.7,66.3,79.7,70,82.1,90.8,95.2,97.5,95.4,89.6,79.7,140,116.7,93.3,70,46.7,23.3,0,0,-36.3,-58.3,-73,-82.6,-88.2,-90];
quiver(x,y,vx,vy,0.4)
xa=[115.9,-119.4];ya=[-269.6,-208.8];u=[0,0];v=[60.8,-60.8];
quiver(xa,ya,u,v,0.4)
xf1=linspace(-51.8,0,30);
yf1=-61/51.8*xf1;
xf2=linspace(-20.8,0);
yf2=-77.3/20.8*xf2;
xf3=-xf1;
yf3=yf1;
xf4=-xf2;
yf4=yf2;
plot(xf1,yf1,'m',xf2,yf2,'m',xf3,yf3,'m',xf4,yf4,'m')