高压下空气密度的编程计算.doc

高压下空气密度的编程计算 (2010-11-16 11:35:13)转载这几天做输气的大作业，要用到高压下（20MPa）空气的密度值，翻书、上网找了半天也没有结果，最多压力就到2.5MPa，于是乎，我痛苦地决定，自己动手，丰衣足食！利用课上讲过的PR方程改进式PRSV2，虽然用BWRS方程解会更精确但没有氧气的有关数据。C语言编程源代码：#include #include #include #define R 8.3143#define e 10e-6voidmain()start: double x2=0.21,0.79;double Tc2=154.77,126.2;double Pc2=5090000.0,3400000.0;double w2=0.02128,0.03726;double k2=0.0,k02=0.0,k12=0.01512,0.01996,k22=-0.009,0.3162,k32=0.49,0.535;double a3=0.0,ac3=0.0,alpha2=0.0,b3=0.0,Tr2=0.0;double P=0.0,T=0.0;double v50=0.0,F50=0.0;int i=0,j=0;char choice=y;printf(nn输入系统压力值（单位kPa）：);scanf(%lf,&P);printf(输入系统温度（单位K）：);scanf(%lf,&T);for(i=0;i=1;i+)Tri=T/Tci;aci=0.45724*pow(R,2)*pow(Tci,2)/Pci;bi=0.0778*R*Tci/Pci;b2+=bi;k0i=0.378893+1.4897153*wi-0.17131748*pow(wi,2)+0.0196544*pow(wi,3);ki=k0i+(k1i+k2i*(k3i-Tri)*(1.0-sqrt(Tri)*(1.0+sqrt(Tri)*(0.7-Tri);alphai=pow(1.0+ki*(1.0-sqrt(Tri),2);ai=aci*alphai;for(i=0;i=1;i+)for(j=0;j=1;j+)a2=xi*xj*sqrt(ai)*sqrt(aj);v0=0.0;F0=R*T/(-b2)+a2/pow(b2,2);v1=R*T/P;for(i=1,j=1;i+,j+)Fi=R*T/(vi-b2)-a2/(vi*(vi+b2)+b2*(vi-b2)-P;vi+1=(vi-1*Fi-vi*Fi-1)/(Fi-Fi-1);if(fabs(vi+1-vi)=e)break;printf(%.0lfkPa、%.2lfK下，空气密度:%.4lfkg/m3,P,T,1/vi*28.84);printf(nn是否继续？（y/n）);scanf( %c,&choice);if(tolower(choice) = y)goto start;程序运行截图：（看起来还是能用滴，随便找个编译器，把上面的代码复制进去运行就OK了，虽然没啥技术含量，精度也不算太高，但还是希望能帮到像我当初一样苦苦寻找高压下空气密度数据的盆友们）这几天因为课程需要编了很多程序，有的运行的很顺利，有的卡了我好几天，你觉得哪儿都对，但就是得不到想要的结果，让你恨不得把电脑砸了，但是自己的东西肯定心疼，想砸又不忍心，真纠结，这是我这段时间明白的第一件事对“纠结”这个词有了更加深刻的感性认识。每个人都有自己不知道错在哪里的时候，正所谓“不识庐山真面目，只缘身在此山中”嘛，后来大家坐下来一起讨论、互相“挑刺