S彎進(jìn)氣道課設(shè)

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