推力室高度特性

1、專(zhuān)題YF-20推力室高度特性的研究作者陳志遠(yuǎn)學(xué)號(hào)12151175院（系）名稱(chēng)宇航學(xué)院 宇航推進(jìn)系火箭發(fā)動(dòng)機(jī)原理大作業(yè) 推力室基本性能計(jì)算2015年4月12日11目錄理論方法3公式推導(dǎo)4程序設(shè)計(jì)5原始數(shù)據(jù)7結(jié)果與分析8參考文獻(xiàn)11理論方法因?yàn)橘|(zhì)量流率守恒，由出口和喉部面積比確定q()，然后利用c語(yǔ)言編制程序，利用二分法或者枚舉法即可反解出速度因數(shù) 。如此一來(lái)，可得氣動(dòng)函數(shù) =PePc 的值。查資料，得到海拔高度和大氣壓強(qiáng)的近似函數(shù)關(guān)系式Pa(h)。利用for循環(huán)，將原始數(shù)據(jù)和以上求出的數(shù)據(jù)代入推力系數(shù)的公式，求出0-80km高度上的推力系數(shù)，代入推力公式即可求出推力。同樣，將特征速度c*與推力系

2、數(shù)一起代入比沖的計(jì)算式，得到每一高度上的比沖?？紤]到c語(yǔ)言畫(huà)圖功能的局限性，畫(huà)出的圖像只能用星號(hào)等散點(diǎn)組成，不美觀也不直接，故將數(shù)據(jù)導(dǎo)出，放入Excel中畫(huà)出推力和比沖的高度特性曲線。公式推導(dǎo)由連續(xù)方程Kp*T*Acr=Kp*qeT*Ae可得q=AtAe=k+121k-1(1-k-1k+12)1k-1用程序反求出e，代入e=PePc=(1-k-1k+12)kk-1得到PePc.查資料得到大氣壓與海拔的近似函數(shù)關(guān)系為：Pa=P0e-h7924將以上各式代入推力系數(shù)公式：CF=2kk-11-pepck-1k+AeAtpepc-papc然后再將推力系數(shù)CF代入比沖的計(jì)算式Is=c*CF即可得到相應(yīng)的

3、比沖.程序設(shè)計(jì)#include <stdio.h>#include "math.h"/預(yù)設(shè)常數(shù)double Pc=6.98E6,ea=12.69,Dt=278E-3,m=263.87,R=380,Tc=3200;float k=1.15;/大氣壓強(qiáng)Pa/double Pa(int h)float p,p0=1.01325e5;p=1.01325E5*pow(2.71828,0.0-h/7924.0);return p;/gama函數(shù)/double gama(float k)double temp;temp=sqrt(k)*pow(2.0/(k+1.0),(k+1

4、.0)/(2*k-2.0);return temp;/由q(lbt)反求lbt（）double q2lbt(double q)double q0;double j;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<q) break;return j;/main函數(shù)void main()double Cf,cstar,F81,Is81,lbt,pp; /pp是Pe/Pcint h,i;lbt=q2lbt(1.0/e

5、a);pp=pow(1-(k-1)*lbt*lbt/(k+1),k/(k-1);cstar=sqrt(R*Tc)/gama(k);printf("h F Isn"); /數(shù)據(jù)記錄+輸出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);Fh/1000=Cf*Pc*0.25*3.14159*Dt*Dt;Ish/1000=Cf*cstar;printf("%-5d%-15f%-15fn",h/1000,Fh/1000,

6、Ish/1000);system("pause");原始數(shù)據(jù)燃燒室壓力=6.98MPa, 噴管面積比=12.69, 噴管喉部直徑=278mm, 推進(jìn)劑質(zhì)量流量=263.87kg/s,燃燒溫度= 3200K, 燃?xì)獗葻岜?1.15, 燃?xì)鈿怏w常數(shù)=380J/kg.K.結(jié)果與分析l 數(shù)據(jù)H/kmF/kNIs/ms-1H/kmF/kNIs/ms-11696.86512840.0541765.21713118.6172705.02132873.29142765.26953118.833712.21052902.5943765.31573119.0184718.54742928.41

7、644765.35643119.1845724.13292951.1845765.39223119.3316729.05632971.24546765.42393119.4597733.39592988.93147765.45173119.5738737.2213004.5248765.47633119.6739740.59263018.26149765.4983119.76110743.56453030.37250765.5173119.83911746.1843041.04851765.53393119.90812748.4933050.45852765.54873119.96813750

8、.52823058.75353765.56183120.02114752.32213066.06454765.57333120.06815753.90333072.50855765.58343120.1116755.29713078.18856765.59243120.14617756.52563083.19557765.60033120.17818757.60843087.60858765.60723120.20719758.56293091.49859765.61343120.23220759.40423094.92760765.61883120.25421760.14583097.949

9、61765.62353120.27322760.79943100.61362765.62773120.2923761.37563102.96163765.63143120.30524761.88343105.0364765.63473120.31925762.3313106.85565765.63763120.3326762.72563108.46366765.64013120.34127763.07343109.8867765.64233120.3528763.37993111.12968765.64433120.35829763.65013112.23169765.6463120.3653

10、0763.88833113.20170765.64763120.37131764.09823114.05771765.64893120.37732764.28333114.81172765.65013120.38133764.44643115.47673765.65123120.38634764.59013116.06274765.65213120.38935764.71683116.57875765.65293120.39336764.82853117.03376765.65363120.39637764.9273117.43477765.65423120.39838765.01383117

11、.78878765.65483120.40139765.09033118.179765.65533120.40340765.15773118.37580765.65573120.404l 推力特性曲線隨著飛行高度的增加，YF-20的推力從696kN變化到765kN，推力隨高度不斷增加，最后趨于穩(wěn)定。根據(jù)發(fā)動(dòng)機(jī)推力公式：F=mue+(Pe-Pa)Ae隨著飛行高度的增加，大氣壓Pa隨之下降，所以推力F隨Pa的下降而增加。同時(shí)由于在我查到的公式中，大氣壓和高度是負(fù)指數(shù)的關(guān)系，因而Pa的下降速度隨高度的增加而減慢，因此推力曲線最后趨于平緩。若火箭完全脫離大氣層，則達(dá)到最大推力Fmax=mue+PeAel 比沖特性曲線隨著飛行高度從海平面到80km，YF-20的比沖從2840 m/s 增加到3120 m/s. 從比沖的計(jì)算式Is=c*CF 其中c*=RTf(k)可