推力室高度特性.doc

《推力室高度特性.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、.专题YF-20推力室高度特性的研究作者陈志远学号12151175院（系）名称宇航学院宇航推进系火箭发动机原理大作业——推力室基本性能计算2015年4月12日..目录理论方法3公式推导3程序设计3原始数据3结果与分析3参考文献3..理论方法因为质量流率守恒，由出口和喉部面积比确定D_Dd______利用二分法或者枚举法即可反解出速度因数()_/_D_Dd的值。查资料，得到海拔高度和大气压强的近似函数关系式利用for循环，将原始数据和以上求出的数据代入推力系数的公式，求出0-80km高度上的推力系数，代入推力公式即可

2、求出推力。同样，将特征速度与推力系数一起代入比冲的计算式，得到每一高度上的比冲。考虑到c语言画图功能的局限性，画出的图像只能用星号等散点组成，不美观也不直接，故将数据导出，放入Excel中画出推力和比冲的高度特性曲线。..公式推导由连续方程可得用程序反求出，代入得到.查资料得到大气压与海拔的近似函数关系为：将以上各式代入推力系数公式：然后再将推力系数代入比冲的计算式..即可得到相应的比冲...程序设计#include#include"math.h"//预设常数doublePc=6.98E6,ea=

3、12.69,Dt=278E-3,m=263.87,R=380,Tc=3200;floatk=1.15;//大气压强Pa//doublePa(inth){floatp,p0=1.01325e5;p=1.01325E5*pow(2.71828,0.0-h/7924.0);returnp;}//gama函数//doublegama(floatk){doubletemp;temp=sqrt(k)*pow(2.0/(k+1.0),(k+1.0)/(2*k-2.0));returntemp;}//由q(lbt)反求lbt（do

4、ubleq2lbt(doubleq)..{doubleq0;doublej;for(j=1.000;j<=2.659;j=j+0.001){q0=pow((k+1.0)/2.0,1.0/(k-1.0))*j*pow((1.0-((k-1.0)*j*j/(k+1.0))),1.0/(k-1.0));if(q0

5、);pp=pow(1-(k-1)*lbt*lbt/(k+1),k/(k-1));cstar=sqrt(R*Tc)/gama(k);printf("hFIs");//数据记录+输出for(h=1000;h<=80000;h=h+1000){Cf=gama(k)*sqrt(2*k/(k-1)*(1-pow(pp,(k-1)/k)))+ea*(pp-Pa(h)/Pc);F[h/1000]=Cf*Pc*0.25*3.14159*Dt*Dt;Is[h/1000]=Cf*cstar;..printf("%-5d%-15f%

6、-15f",h/1000,F[h/1000],Is[h/1000]);}system("pause");}原始数据燃烧室压力=6.98MPa,喷管面积比=12.69,喷管喉部直径=278mm,推进剂质量流量=263.87kg/s,燃烧温度=3200K,燃气比热比=1.15,燃气气体常数=380J/kg.K...结果与分析l数据H/kmF/kNIs/ms-1H/kmF/kNIs/ms-11696.86512840.0541765.21713118.6172705.02132873.29142765.2695311

7、8.833712.21052902.5943765.31573119.0184718.54742928.41644765.35643119.1845724.13292951.1845765.39223119.3316729.05632971.24546765.42393119.4597733.39592988.93147765.45173119.5738737.2213004.5248765.47633119.6739740.59263018.26149765.4983119.76110743.56453030.3

8、7250765.5173119.83911746.1843041.04851765.53393119.90812748.4933050.45852765.54873119.96813750.52823058.75353765.56183120.02114752.32213066.06454765.57333120.06815753.90333072.5085