压缩机部件MATLAB仿真

admin

[mw_shl_code=applescript,true]L=0.054;  
R=0.034;  
e=0.005;  
n=50;  
Tc=327.55;  
Te=280.15;  
Tsuc=Te+5;%输入初始参数  
Pc=refpropm('P','T',Tc,'Q',1,'R22');%已知Q=1和Tc，求Pc  
Pe=refpropm('P','T',Te,'Q',1,'R22');%已知Q=1和Te，求Pe  
k=refpropm('K','T',Te,'Q',1,'R22');%求进口处k  
Td=Tsuc*(Pc/Pe)^((k-1)/k);%排气温度计算  
m0=refpropm('D','T',Tsuc,'P',Pe,'R22');%求进气密度  
hsuc=refpropm('H','T',Tsuc,'P',Pe,'R22')/1000;%求进气比焓kJ  
hend=refpropm('H','T',Td,'P',Pc,'R22')/1000;%中止比焓  
h1=refpropm('H','T',Te,'Q',1,'R22')/1000;%入口工质比焓  
h2=refpropm('H','P',Pc,'Q',0,'R22')/1000;%饱和液体比焓，以上为利用refpropm函数调用部分  
e0=e/R;%偏心率  
vsuc=1/m0;%比容  
Vth=60*n*3.1415926*R*R*L*(2-e0);%理论输气量
dPc=25*(Tc-273.15)^(-1.01)*10^(-0.15*e0);%排气压损  
Iv=1-0.015*(((Pc-dPc)/Pe)^(1/k)-1);  
It=0.00257*Tc+0.00106*(Tsuc-Te);  
Id=0.87;  
I=Iv*It*Id*1;  
Mcom=(I*Vth)/vsuc;%以上为流量计算  
dpe=0;  
Li=(It*Id)/(1-(1.5*(dpe+dPc*(Pc/Pe)^(1/k)))/((hend-hsuc)/vsuc));  
Lt=It;  
Ll=Id;  
Lm=0.7;  
Lmo=0.8;  
Lel=Li*Lt*Ll*Lm*Lmo;  
Nth=I*Vth*Pe*k/(k-1)*((Pc/Pe)^(1-1/k)-1);  
Nin=Nth/Lel;%以上为功率计算  
Q=Mcom*(h1-h2);%制冷量  
COP=Q/Nin;%COP计算  
[/mw_shl_code]