压缩机示功图编程MATLAB.doc

压缩机示功图测试实验报告 表1 试验压缩机性能参数 压缩机的型号名称 压缩机额定排气压力（kgf/cm2） 压缩机额定排气量（m3/min） 电动机型号名称 电动机装置功率（kw） 电动机与压缩机的传动方式 压缩机活塞行程L (mm) 60 压缩机连杆长度L1 (mm) 50 压缩机活塞直径D(mm) 80 压缩机额定转速N（rpm） 电机效率ηe (%) 0.9 传动效率ηz(%) 0.97 表2 压缩机示功图技术数据 1. 试验日期 2. 大气压力（pa） 101325 3. 压缩机试验排气压力Pd (kgf/cm2) 604.860Kpa 4. 压缩机试验转速n (rpm) 939.830 5. 行程长度L (mm) 60 6. 吸入线长度L1 (mm) 42.266 7. 排气损失面积Fp（mm*100Kpa） 4.4942 8. 吸气损失面积Fx（mm*100Kpa） 19.773 9. 示功图面积Fi（mm*100Kpa） 135.7 10. 示功图名义面积Fh（mm*100Kpa） 111.43 11. 平均指示压力Pi=Fi/L（100Kpa） 2.2616 12. 指示功率Ni (kw) 1.0684 13. 电机输入功率Ne (kw ). 6 14. 压缩机轴功率Nz=Ne *ηe *ηz(kw) 5.076 15. 压缩机机械效率ηm=3 N i / Nz (%) 63.144 16. 容积效率λv= L1 /L (%) 70.443 17. 气阀功率损失系数 △C = 100 (Fx + Fp ) / Fh（%） 21.778 附录1：数据处理示例程序片段 MATLAB代码 global A Py Pq Pt ND N Xy Xt G A=importdata(yasuoji.data); Py=A.data(1); Pq=A.data(2); Pt=A.data(3); ND=A.data(4); N=A.data(5); Xy=A.data(6)-101.325; Xt=A.data(7); G=A.data(8); L=60; L1=50; D=80; dx=0.9; cx=0.97; global j global y global x global p for i=1:G p(i)=A.data(i+8)*100; j(i)=2*(i-1); t=j(i)/180*pi; y(i)=(L/2/L1)*sin(t); x(i)=L/2+L1-L/2*cos(t)-L1*cos(atan((y(i)/sqrt(1-y(i)^2)/180*pi))); end x(G+1)=x(1); p(G+1)=p(1); plot(x,p); X=[0 70]; P1=[0 0]; P2=[Pq Pq]; P3=[Xy Xy]; line(X,P1); line(X,P2); line(X,P3); xlabel(DD3ì??ò?x/mm) ylabel(?1??p/Kpa) title(?1???úê?1|í?) global l1 l2 l3 l4 l5 l6 for i=1:G a=p(i); b=p(i+1); if a=XybXy n1=i; l1=x(i); end if a=XybXy n2=i; l2=x(i); end l3=l2-l1; if a=PqbPq n3=i; l4=x(i); end if a=PqbPq n4=i; l5=x(i); end l6=l4-l5; end disp(l3); disp(l6); s=polyarea(x,p)/100; x1=x(n1:n2); x2=x(n3:n4); p0=p(n1:n2); p3=p(n3:n4); p1=abs(p(n1:n2)-Xy); p2=abs(p(n3:n4)-Pq); p0(1)=Xy; p0(length(x1))=Xy; p3(1)=Pq; p3(length(x2))=Pq; s1=abs(trapz(x1,p1))/100; s2=abs(trapz(x2,p2))/100; patch(x1,p0,b) patch(x2,p3,g) Ni=s*(pi*(D^2)/4)* N/60*10^-7; Pi=s/L; di