原创力文档- 1、本文档共3页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
油库设计与管理上机作业中国石油大学(华东)
油气储运工程10-1班 岳韬 #includestdio.h
#includemath.h
void main()
{
double tys0,rout,tbi0,Fjr1,Fjr2,tbi1=-0.01,g=9.8,a1bi,D=15.25,Bta,Gr,Pr,A,ty,ipsil,n,e=2,t,d4t,miut,lamidat,cd,tj=-3.0,tys1=-0.01,tyz=40.0,H=11.74,h,Fbi,Fdi,Fding,F,tzl=-2,v=3.0,miuqi,Re,C,a3bi,a2bi,Kbi,K,cp,Kdi,Kding,deta=2,routy,cty,Q1,Q3,Q,tjrbi1=100,tjrbi0,ct,a2,K0,Gz1,Gz2,lamidaqi;
h=H*0.85;
Fbi=3.1415926*D*0.85*H;
Fdi=0.25*3.14*D*D;
Fding=(H-h)*3.1415926*D+0.25*3.1415926*D*D;
//总传热面积
F=Fbi+Fdi+Fding;
//油品加热的初始温度
while(deta=0.01tys1=40.0)
{ tbi1=-0.01;
e=2;
tys1=tys1+0.005;
if((tyz-tj)/(tys1-tj)2) ty=tj+(tyz-tys1)/(log((tyz-tj)/(tys1-tj)));
else ty=(tys1+tyz)/2.0;//油品平均温度
//管壁传热系数
while(e=1tbi1=40.0)
{
tbi1=tbi1+0.005;
t=0.5*(ty+tbi1);
//温度t时相对密度
d4t=0.89-0.00062*(t-15);
//粘温特性
miut=(pow(10,pow(10,8.921-3.439*log10(t+273)))-0.6)/1000000;
lamidat=117.5/890*(1-0.00054*t);
//比热容
cd=2.018+0.00322*(t-100);
Bta=(2.73355-2.25355*d4t)/1000;
Gr=(h*h*h*Bta*(ty-tbi1)*g)/(miut*miut);
Pr=1000000*miut*d4t*cd/lamidat;
A=Pr*Gr;
if(A=0.001A500){ipsil=1.18;n=0.125;}
if(A=500{ipsil=0.54;n=0.25;}
if({ipsil=0.135;n=0.3333;}
a1bi=ipsil*lamidat/h*pow(A,n);
lamidaqi=(2.361+(2.442-2.361)/10.0*(tzl+10))/100;
miuqi=(12.43+(13.28-12.43)/10*(tzl+10))/1000000;
Re=D*v/miuqi;
if(Re=5Re=30){C=0.81;n=0.4;}
if(Re80Re=5000){C=0.625;n=0.46;}
if(Re5000Re=50000){C=0.197;n=0.60;}
if(Re50000){C=0.023;n=0.8;}
a2bi=C*lamidaqi/D*pow(Re,n);
a3bi=0.96*5.67*(pow((tbi1+273)/100,4)-pow((tzl+273)/100,4))/(tbi1-tj);
Kbi=1/(1/a1bi+1/(a2bi+a3bi));
tbi0=-Kbi/a1bi*(ty-tj)+ty;
e=fabs(tbi0-tbi1);
}
Kding=0.7;
Kdi=0.12;
K=(Fdi*Kdi+Kbi*Fbi+Fding*Kding)/F;
cp=2.018+0.00322*(ty-100);
routy=(0.89-0.00062*(ty-15))*1000;
tys0=tj+(tyz-tj)*pow(2.718,-14*3600*24*K*F/(0.25*3.14*D*D*h*routy)/cp/1000);
deta=fabs(tys0-tys1);
}
printf(起始加热温度tys=%f\n罐壁温tbi1=%f\n罐传热系数K=%f\n
文档评论（0）