推力室高度特性.docxVIP

北京航咳航夭大爭 B E I H A N G UNIVERSITY 火箭发动机原理大作业 ——推力室基本性能计算 YF?20推力室高度特性的研究陈志远院（系）名称宇航学院宇航推进系 院（系）名称 2015年4月12日 目录 TOC \o 1-5 \h \z \o Current Document 理论方法 3 \o Current Document 公式推导 4 \o Current Document 程序设计 5 \o Current Document 原始数据 7 \o Current Document 结果与分析 8 \o Current Document 参考文献 11 理论方法 因为质量流率守恒，由出口和喉部而积比呦定q(A),然后利用c语言编制程 序，利用二分法或者枚举法即可反解出速度因数入。如此一來，可得气动函数 = 2的值。 査资料，得到海拔高度和大气压强的近似函数关系式Pa(h)o 利用for循环，将原始数据和以上求出的数据代入推力系数的公式，求出 0-80km高度上的推力系数，代入推力公式即可求出推力。 同样，将特征速度c?与推力系数一起代入比冲的计算式，得到每一高度上的 比冲。 考虑到c语言画图功能的局限性，画出的图像只能用星号等散点组成，不美 观也不直接，故将数据导出，放入Excel中倆岀推力和比冲的高度特性曲线。 公式推导 PP(A)yfT7 PP(A)yfT7 *心=K 由连续方程 可得 用程序反求出儿，代入 Pe k-1 0 fc 咫)飞y — e*尸 得到去 查资料得到大气圧与海拔的近似函数关系为: 将以上各式代入推力系数公式: 然后再将推力系数6代入比冲的计算式 即可得到相应的比冲. 程序设计 include stdio.h include Mmath.hH 〃预设常数 double Pc=6?98E6,=12?69,Dt=278E-为01=263.87,R=?8o,Tc=32OO； float k=i.i5； 〃大气压强Pa〃 double Pa(int h) { float p,po=i?(n 72525; p=i.oi325E5*pow(2.71828,0.041/7924.0); return p; } //gam a 函数// doublegama(float k) { double temp; temp=sqrt(k)*pow(2.o/(k+i.o)Xk+i.o)/( *k-.c)))； return temp; } 〃山q(lbt)反求Ibt (A) double q2lbt(double q) double qo; double j; for(j=iooo;j=2.6594=j+o.ooi) qo=pow((k+i.o)/2.o,i.o/(k-i.o))*j*pow((i o-((k-d.o)*j*j/(k+i.o))),i.o/(k-i.o)); if(qoq)break; } return j; } //main函数 void main() { doubleCf,cstar,F[8i],ls[ i]」bt,pp;//pp是Pe/Pc inthj； Ibt=q2lbt( .o/ea); pp=pow(i-(k-i)*lbt*lbt/(k+i),k/(k-i)); cstar=sqrt(R*Tc)/gama(k); printf(Mh F ls\nH);// 数据记录 +输出 for(h=iooo; h=8oooo; h=h+iooo) { Cf=gama(k)*sqrt( *k/(k-i)*( -pow(pp,(k-)/k)))+ea*(pp-Pa(h)/Pc); F[h/iooo]=Cf*Pc*o.25*s.i4i59*Dt*Dt; ls[h/iooo]=Cf*cstar； printf( ,X-5d%-i5fX-i5f\n,,fh/woo/F[h/iooo],ls[h/ioo(D； } system(MpauseH); 原始数据 原始数据 燃烧室床力pc=6.98MPa/ 喷管面积比6 2.69, 喷管喉部直径4 =278mm/ 推进剂质量流量lii =263.87kg/s/ 燃烧温度Tc = 3200K, 燃气比热比k =1.15, 燃气气体常数R =380J/kg.K. 结果与分析 数据 H/km FAN Is/ms-1 H/km F/kN Is/msT 1 696. 8651 2840. 05 41 765. 2171 3118.617 2 705. 0213 2873. 291 42 765. 2695 3118. 83 3 712. 2105 2902. 59 43 765. 3157 3119. 018