S彎進氣道課設(shè)

課程設(shè)計報告題目（s彎進氣道設(shè)計）班級：0207102學(xué)號：020710214姓名：唐鴻雁設(shè)計過程：1:計算參數(shù)：已知發(fā)動機進口直徑（進氣道出口直徑）為137.5mm，擴壓比為1.3，據(jù)此可以算的，進氣道進口半徑是60.3mm。2：選擇參數(shù)方程中性線方程： 進氣道面積變化規(guī)律：式中：y為擴壓器中心線的縱坐標(biāo)；為擴壓器的縱向偏距300mm；x為擴壓器中心線的橫坐標(biāo)；L為擴壓器的長度900mm。代入已知參數(shù)，上述公式變?yōu)椋簓=x/1350-x/2187000000;A=3447.43*[3*(x/900)-8*(x/900)+6*(x/900)]+11423.1;r=sqrt(A/pi)pi=3.14159263：畫出中心線思路：將中心線分成90段，即將x的值從0，每次增加10，一直加到900，再根據(jù)中心線方程，算出對應(yīng)的y坐標(biāo)，這一功能可以通過CPP編程來實現(xiàn)。CPP源代碼如下：#include<iostream.h>#include<fstream.h>voidmain(){ inta;fstreamout;char*f;f="1.dat";/*創(chuàng)建文件來保存坐標(biāo)值，以方便導(dǎo)入UG*/out.open(f,ios::out); doublex[92]={0},y[92]; for(a=0;a<91;a++) { x[a+1]=x[a]+10; y[a]=x[a]*x[a]/1350-x[a]*x[a]*x[a]*x[a]/2187000000; } for(a=0;a<91;a++) {out<<x[a]<<""<<y[a]<<""<<0<<'\n';/*將坐標(biāo)值輸入文件*/ }}中性線UG圖如下：4：畫出進氣道思路：中心線被分成了九十段，將每段看成一條直線，可以算出它的斜率，同樣可以得到與之垂直的面。可以在此面上做圓，圓的半徑即進氣道對應(yīng)截面的半徑，r=sqrt(A/pi)。從每一個圓上找出若干個點，如90個點，這樣，共有九十個圓，8100個點，導(dǎo)入UG，可以得到進氣道截面。CPP源代碼如下：#include<iostream.h>#include<fstream.h>#include<math.h>#definepi3.1415926voidmain(){ inta,b;fstreamout;char*f;f="2.dat";out.open(f,ios::out); doublex[92]={0},A[92]={11423.1},r[92],y[92],k[92],x1[92],y1[92],z1[92],jiao1[92]={0},jiao2[92]; for(a=0;a<91;a++) { jiao1[a+1]=jiao1[a]+4*pi/180; x[a+1]=x[a]+10; y[a]=x[a]*x[a]/1350-x[a]*x[a]*x[a]*x[a]/2187000000; } for(a=1;a<91;a++) { A[a]=3447.43*(3*(x[a]/900)*(x[a]/900)*(x[a]/900)*(x[a]/900)-8*(x[a]/900)*(x[a]/900)*(x[a]/900)+6*(x[a]/900)*(x[a]/900))+11423.1; r[a]=sqrt(A[a]/pi); } for(a=0;a<91;a++) { k[a+1]=(x[a]-x[a+1])/(y[a+1]-y[a]);/*對應(yīng)直線的斜率*/ jiao2[a+1]=atan(k[a+1])+pi/2; }for(b=0;b<90;b++)/*將九十個圓循環(huán)*/ { for(a=0;a<90;a++) { x1[a+1]=x[b+1]+r[b+1]*cos(jiao1[a+1])*sin(jiao2[b+1]); z1[a+1]=r[b+1]*sin(jiao1[a+1]); y1[a+1]=y[b+1]-r[b+1]*cos(jiao1[a+1])*cos(jiao2[b+1]); } for(a=1;a<91;a++) { out<<x1[a]<<""<<y1[a]<<""<<z1[a]<<'\n';/*將每個圓上的九十個點按順序輸入文件中*/ } }}5：畫出唇口只要將所給的數(shù)據(jù)輸入dat文件中，導(dǎo)入UG做旋轉(zhuǎn)即可得到。6：將唇口和進氣道裝配后縫合，得到UG圖如下所示7：添加彈身后UG圖如下：8：性能參數(shù)計算查表可知：在6000m高空處，p=0.47218pa，ρ=0.6601kg/m，K=1.4，當(dāng)馬赫數(shù)M=0.65時，p/p=0.7528。p=0.62723pa從流體力學(xué)中可知：p/p=[1+(k-1)M/2]由此可以得到進口馬赫數(shù)為：M1=0.65已知：總壓回復(fù)系數(shù)σ=0.98，則出口總壓p=0.61468pa則出口馬赫數(shù)為：M2=0.62.9：課設(shè)總結(jié)及心得本次課設(shè)的題目并不是很難，但是在此過程中也還是遇見不少的問題。如在cpp編完，得到的文件怎么導(dǎo)入Ug，獲得型面等。在cpp程序編好的時候，就在我以為可以一次性出圖的時候，結(jié)果在UG中只出現(xiàn)了一些相應(yīng)的點，卻沒有出現(xiàn)面，雖然大概的形狀是正確的，但是就是沒有辦法得到型面，這讓我非常著急。后來，檢查程序發(fā)現(xiàn)，由于循環(huán)的時候多輸入了一個圓，而那個圓的半徑是無窮大的，所以在導(dǎo)入的時候Ug才提醒“極點”的問題。最后程序一改，果然成功了??！本次課設(shè)，讓我重新溫故了CPP和UG兩種軟件，也讓我明白各種軟件并不是相互獨立的，將他們聯(lián)系起來，往往可以實現(xiàn)一些非常強大的功能！參考文獻1王珊珊主編C++程序設(shè)計教程機械工業(yè)出版社200