| 查看: 963 | 回复: 16 | |||
| 当前主题已经存档。 | |||
li343至尊木虫 (职业作家)
|
[交流]
【求助】哪位英雄帮我看一下这个小程序【已完结】
|
||
|
程序编译正常，有结果（不知道是否正确），但发现v一直为1，无论温度T怎么给，v都不变，这就是个问题了，请高手帮看一下是啥原因。 先谢谢了！ #include #include void main() { int i,j; double t,a1,a2,a3,P,Pc,Tc,T,Tr; double c0,c1,c2,c3,Mole,a_T,a1_T,v,R,X0,w; double m,ac,b,Afa_T; double h1,h2,H; a1=6.6400537; a2=8.3412027; a3=27.498677; if(-100 c0=3.089329; c1=0.1020401; c2=-2.400100e-4; c3=4.126521e-7; } if(298<=T<1500) { c0=-2.207131; c1=0.1002983; c2=-0.5530711e-4; c3=0.1187821e-7; } X0=249.4758; Tc=408.16; //临界温度/K Pc=3684e+3; //临界压力/kPa Mole=58.12; //分子量 w=0.1756; //偏心因子 R=8.314; //气体常数 T=338.15; Tr=T/Tc; t=1-Tr; m=0.37464+1.54226*w-0.22992*w*w; ac=0.45724*pow((R*Tc),2)/Pc; b=0.07780*R*Tc/Pc; Afa_T=pow((1+m*(1-pow(Tr,0.5))),2); a_T=ac*Afa_T; a1_T=-m*a_T*pow(Tr,0.5)*pow(Afa_T,-0.5)/T; double V[1000]; V[0]=0.7; for(i=0;i<=1000;i++) { V[i+1]=pow(((R*T-P*b)*pow(V,2)+(2*R*T*b+3*P*b-a_T)*V-(P*pow(b,3)+R*T*pow(b,2)-a_T*b))/P,(1/3)); v=V; } printf("%f\n",m); printf("%f\n",ac); printf("%f\n",b); printf("%f\n",Afa_T); printf("%f\n",a_T); printf("%f\n",a1_T); printf("\n"  ;P=Pc*exp((a1+a2*pow(t,1.78)+a3*pow(t,6.06))*log(Tr)); //饱和蒸汽压,Pa h1=c0*T+c1*pow(T,2)/2+c2*pow(T,3)/3+c3*pow(T,4)/4; h2=(a_T-T*a1_T)*log((v-0.4142*b)/(v+2.4142*b))/(2*sqrt(2)*a1_T*b); H=4.1868*h1/Mole+h2+P*v-R*T+X0; //焓值,J/(kg*K) printf("%lf\n",P/1000); printf("%lf\n",v); printf("%lf\n",H/1000); } [ Last edited by 余泽成 on 2009-11-26 at 20:25 ] |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
sitonyl
铁杆木虫 (文坛精英)
- 应助: 33 (小学生)
- 金币: 12935.9
- 红花: 3
- 帖子: 10448
- 在线: 338.8小时
- 虫号: 722261
- 注册: 2009-03-14
- 性别: GG
- 专业: 导航、制导与传感技术
2楼2009-07-30 10:48:09
li343
至尊木虫 (职业作家)
- 应助: 32 (小学生)
- 金币: 16345.4
- 散金: 138
- 红花: 11
- 沙发: 1
- 帖子: 3033
- 在线: 998.3小时
- 虫号: 579854
- 注册: 2008-07-11
- 性别: GG
- 专业: 工程热物理相关交叉领域
3楼2009-07-30 11:32:18
4楼2009-07-30 12:35:46
5楼2009-07-30 12:45:29
6楼2009-07-30 13:02:43
7楼2009-07-30 21:14:03
tjpm
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1025.7
- 散金: 10
- 红花: 1
- 帖子: 367
- 在线: 19.3小时
- 虫号: 717334
- 注册: 2009-03-07
- 专业: 凝聚态物性 II ：电子结构
8楼2009-07-30 22:31:35
9楼2009-08-07 09:59:42
10楼2009-08-09 21:10:26