原创力文档- 1、本文档共6页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
内燃机热力计算程序(c语言)
#includestdio.h
#includemath.h
#includestdlib.h
const double Pi=3.1416,e=18.2,D=0.093,S=0.102,L=0.152,r=0.335,n=3600,R=287.08,Hugb=3.5*0.00001,L0=14.4,mp=0.18,tr=350.0,o=1.0;
double V(double ca) /*体积函数 v*/
{
double v,x,y;
x=2.0/(e-1.0)+1.0-cos(ca*Pi/180.0);
y=(1.0/r)*(1.0-sqrt(1.0-r*r*sin(ca*Pi/180.0)*sin(ca*Pi/180.0)));
v=(D*D*S/8.0)*Pi*(x+y);
return (v);
}
double AG(double v,double m,double T,double de)
{
double ag;
ag=0.205*(1+0.09)*pow((m*R*T*0.000001/v),0.7)*pow((S*n/30.0),0.7)/pow(T,0.2)*pow(de,0.3);
return (ag);
}
double Qw_ca(double ag,double v,double T) /*气缸对外界的瞬时传热dQw/d?*/
{
double qw_ca,Tw1,Tw2,Tw3,Pe=18.4/Pi;
Tw1=120.0+30.0*Pe+273.0;
Tw2=100.0+7.0*Pe+273.0;
Tw3=100.0+4.0*Pe+273.0;
qw_ca=(ag/(6.0*n))*((Pi*D*D/4.0)*(Tw1-T)+(Pi*D*D/4.0)*(Tw2-T)+(4.0/D)*v*(Tw3-T));
return (qw_ca);
}
double V_ca(double ca,double v) /******气缸容积的瞬时变化率 dv/d?******/
{
double v_ca,x,y;
x=sin(ca*Pi/180.0);
y=r*sin(ca*Pi/180.0)*cos(ca*Pi/180.0)/sqrt(1.0-r*r*sin(ca*Pi/180.0)*sin(ca*Pi/180.0));
v_ca=(D*D*S/8.0)*Pi*(x+y)*Pi/180.0;
return (v_ca);
}
double U_T(double T,double a) /*?u/?T*/
{
double u_T,x,y;
x=-3.0*(0.0957+0.0485/pow(a,0.75))*(T-273.0)*0.000001;
y=2.0*(7.768+3.36/pow(a,0.8))*(T-273.0)*0.0001;
u_T=144.55*(x+y+(489.6+46.4/pow(a,0.93))*0.01);
return (u_T);
}
double U(double a,double T) /*工质的内能 u*/
{
double x,y,z,u;
x=-(0.0975+0.0485/pow(a,0.75))*(T-273.0)*(T-273.0)*(T-273.0)*0.000001;
y=(7.768+3.36/pow(a,0.8))*(T-273.0)*(T-273.0)*0.0001;
z=(489.6+46.4/pow(a,0.93))*(T-273.0)*0.01;
u=144.55*(x+y+z+1356.8);
return (u);
}
double QB_ca(double ca) /*燃料燃烧时的瞬时放热率 dQB/d?*/
{
double qb_ca,y;
y=(ca-tr)/60.0; //?
qb_ca=Hu*gb*6.908*pow(y,mp)*exp(-6.905*pow(y,(1.0+mp)))*((1.0+mp)/60.0);
return (qb_ca);
}
double U_a(double a,double T) /*?u/?入*/
{
double x,y,z,u_a;
x=0.75*(0.0485/pow(a,1.75))*(T-273.0)*(T-273.0)*(T-273.0)*0.000001;
y=-0.8*(3.36/pow(a,1.8))*(T-273.0)*(T-273.0)*0.00
文档评论（0）