工热上机.doc

下载文档 降价啦

5
0
约5.19千字
约 7页
2018-05-28 发布于河南
举报
版权申诉
保障服务

工热上机.doc

1、本文档共7页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

工热上机

郎肯循环 #includestdio.h #includemath.h void main() { double p1,p2,t1,t2,h1,h2,s1,s2,h3,sy,sx,hy,hx,nt,x,nex; printf(Input nex,p1(Mpa),p2(Mpa)和t1(oC):\n); scanf(%lf%lf%lf%lf,nex,p1,p2,t1); t1=t1+273; p1=p1*10; h1=2018.24+1.693*t1+0.0002721*t1*t1+(3.643-594.5*p1)*pow(t1/100,-3.1)+(0.6156-2.696*pow(10,6)*pow(p1,3))*pow(t1/100,-13.5); s1=1.693*log(t1)-0.4795*log(p1)-2.9347+0.0005442*t1+(0.02747-4.495*p1)*pow(t1/100,-4.1)+(0.005731-2.51*10000*pow(p1,3))*pow(t1/100,-14.5); sy=0.2776*log(p2)+1.9308; sx=-0.3445*log(p2)+6.5627; hy=94.298*log(p2)+627.9; hx=38.619*log(p2)+2762.3; t2=22.502*log(p2)+149.88; s2=s1; h3=hy; if(s2sx) {x=1; h2=2018.24+1.693*t2+0.0002721*t2*t2+(3.643-594.5*p2)*pow(t2/100,-3.1)+(0.6156-2.696*pow(10,6)*pow(p2,3))*pow(t2/100,-13.5); } else if(s2sx) {x=(s2-sy)/(sx-sy); h2=(1-x)*hy+x*hx; } nt=nex*(h1-h2)/(h1-h3); printf(\n朗肯循环热效率为%lf 。\n,nt); } 第一题：n=0.390885 第二题：见截图第三题：n=0.332252 再热循环 #includestdio.h #includestdlib.h #includemath.h float s(float p,float t) { float sx; sx=1.693*log(t)-0.4795*log(p)-2.9347+0.0005442*t+(0.02747-4.495*p)*pow((t/100),-4.1)+(0.005731-2.51*pow(10,4)*pow(p,3))*pow((t/100),-14.5); return (sx); } float h(float p,float t) { float hx; hx=2018.24+1.693*t+0.0002721*pow(t,2)+(3.634-594.5*p)*pow((t/100),-3.1)+(0.6156-2.696*pow(10,6)*pow(p,3))*pow((t/100),-13.5); return (hx); } float t(float p,float sm) { float t=273.149,s; do { t=t+0.001; s=1.693*log(t)-0.4795*log(p)-2.9347+0.0005442*t+(0.02747-4.495*p)*pow((t/100),-4.1)+(0.005731-2.51*pow(10,4)*pow(p,3))*pow((t/100),-14.5); }while(abs(10000*(s-sm))/10000.00.0001); return (t); } void cal() { float n,x,p1,p2,t1,h1,h2,h3,h5,s1,s2,s3,s5,nex; float pa,pb,ta,tb,t2,ha,hb,sa,s6,s7,h6,h7; printf(nex p1(MPa) t1(℃) p2(MPa) pa(MPa) ta(℃) \n); scanf(%f%f%f%f%f%f,nex,p1,t1,p2,pa,ta); t1=t1+273.15; ta=ta+