




版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
航空發(fā)動(dòng)機(jī)強(qiáng)度振動(dòng)上機(jī)作業(yè)題題目三班級(jí):140411姓名:蘇雨學(xué)號(hào):14041032一:題目要求3-1某轉(zhuǎn)子葉片根部固定,其材料密度ρ=2850kg/m3,彈性模量E=71.54GPa,葉片長(zhǎng)0.1m,各截面位置、面積、慣性矩列于下表,試求其前3階固有靜頻。二:分析公式如題目所示,已知材料密度,彈性模量,各截面位置、面積、慣性矩,需要求解3階靜頻率。求解靜頻公式如下:1.彈性線歸一化2.振型逼近法3.二階振型和固有頻率求解4.三階振型和固有頻率求解三:編程計(jì)算程序使用c語(yǔ)言編寫,源代碼如下:#include<stdio.h>#include<math.h>#include<stdlib.h>intmain(void){floatrou=2850;floatE=71540000000;floatX[11]={0.0,0.01,0.02,0.03,0.04,0.05,0.06,0.07,0.08,0.09,0.10};floatA[11]={0.00017,0.000146,0.000126,0.000109,0.000096,0.000086,0.000077,0.000073,0.00007,0.000068,0.000068};floatI[11]={0.000000000279,0.000000000212,0.000000000157,0.000000000108,0.000000000084,0.000000000061,0.000000000045,0.000000000037,0.000000000032,0.000000000030,0.000000000030};floatAb[10];floatIb[10];inti=0;while(i<=9){ Ab[i]=((A[i+1]+A[i])/2); i=i+1;}i=0;while(i<=9){ Ib[i]=((I[i+1]+I[i])/2); i=i+1;}floatY0js[10]={0,0.1,0.2,0.4,0.5,0.6,0.7,0.8,0.9,1};floatY0sj[10]={0,0,0,0,0,0,0,0,0,0};floatY0xxx[10]={0,0,0,0,0,0,0,0,0,0};floatone[10]={0,0,0,0,0,0,0,0,0,0};floattwo[10]={0,0,0,0,0,0,0,0,0,0};floatthree[10]={0,0,0,0,0,0,0,0,0,0};floatfour[10]={0,0,0,0,0,0,0,0,0,0};floatwucha=10;floatxz=2.58;inta=0;intb=0;intc=0;intd=0;inte=0;intf=0;intg=0;inth=0;intj=0;intk=0;intp=0;intq=0;while(wucha>=0.000001)//給定擬合精確度{ q=0; while(9>=q){ one[q]=0; two[q]=0; three[q]=0; four[q]=0; q=q+1;}//這步是給數(shù)組清零,千萬(wàn)不能忘! a=0;b=9;while(9>=a){ b=a; while(9>=b+1){ one[a]=one[a]+(Ab[b+1]*Y0js[b+1]*0.01); b=b+1; }a=a+1; }//第一重循環(huán)c=0;d=0;while(9>=c){ d=c; while(9>=d+1){ two[c]=two[c]+(one[d+1]*0.01); d=d+1;}c=c+1;}//第二重循環(huán)e=0;f=0;while(9>=e){ f=0; while(e>=f){ three[e]=three[e]+two[f]*0.01*(1/Ib[f]); f=f+1;}e=e+1;}//第三重循環(huán)g=0;h=0;while(9>=g){ h=0; while(g>=h){ four[g]=four[g]+three[h]*0.01; h=h+1;}g=g+1;}//第四重循環(huán)k=0;while(9>=k){ Y0sj[k]=four[k]/four[9]; Y0xxx[k]=four[k];k=k+1;}//求出實(shí)際y0,以便和假設(shè)yo對(duì)比迭代wucha=0;j=0;while(j<=9){ wucha=wucha+fabs(Y0js[j]-Y0sj[j]); j=j+1;}//假設(shè)的y0與求出y0之間的誤差p=0;while(9>=p){ Y0js[p]=Y0sj[p];p=p+1;}//令實(shí)際值等于假設(shè)值,再次迭代運(yùn)算}floatomega=0;omega=5000*sqrt((1/four[9]));printf("一階固有靜頻為:%.5fHZ\n",omega/(2*3.1415926));floatY2js[10]={0,-0.1,-0.3,-0.5,-0.4,-0.2,-0.1,0.5,0.8,1};floatY2[10]={0,0,0,0,0,0,0,0,0,0};floatxiuzhen=4;floatY2sj[10]={0,0,0,0,0,0,0,0,0,0};floatb11=0;floatC21=0;floata21=0;a21=10;while(fabs(a21)>=0.00000001)//給定擬合精確度{ b11=0; C21=0; a21=0; i=0;while(i<=9){ b11=b11+2850*Ab[i]*Y0xxx[i]*Y0xxx[i]*0.01; i++;}i=0;while(i<=9){ C21=C21+2850*Ab[i]*Y0xxx[i]*Y2js[i]*0.01; i++;}a21=C21/b11;i=0;while(i<=9){ Y2[i]=Y2js[i]-a21*Y0xxx[i]; i++;}p=0;while(9>=p){ Y2js[p]=Y2[p];p=p+1;}//令實(shí)際值等于假設(shè)值,再次迭代運(yùn)算}i=0;while(i<=9){Y2sj[i]=Y2[i]/Y2[9];i++; }//歸一化q=0; while(9>=q){ one[q]=0; two[q]=0; three[q]=0; four[q]=0; q=q+1;}//這步是給數(shù)組清零,千萬(wàn)不能忘!a=0;b=0;while(9>=a){ b=a; while(9>=b+2){ one[a]=one[a]+(Ab[b+2]*Y2sj[b+2]*0.01); b=b+1;}a=a+1;}//第一重循環(huán)c=0;d=0;while(9>=c){ d=c; while(9>=d+1){ two[c]=two[c]+(one[d+1]*0.01); d=d+1;}c=c+1;}//第二重循環(huán)e=0;f=0;while(9>=e){ f=0; while(e>=f){ three[e]=three[e]+two[f]*0.01*(1/Ib[f]); f=f+1;}e=e+1;}//第三重循環(huán)g=0;h=0;while(9>=g){ h=0; while(g>=h){ four[g]=four[g]+three[h]*0.01; h=h+1;}g=g+1;}//第四重循環(huán)omega=sqrt((E/rou)*(1/four[9]));omega=omega*xz;printf("二階固有靜頻為:%.5fHZ\n",omega/(2*3.1415926));floatY3js[10]={0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1};floatY3[10]={0,0,0,0,0,0,0,0,0,0};floatY3sj[10]={0,0,0,0,0,0,0,0,0,0};b11=0;floatb12=0;floatb21=0;floatb22=0;floatC31=0;floatC32=0;floata31=0;floata32=0;a32=10;floatshoulian=10;while(shoulian>=0.0001)//給定擬合精確度{ while(i<=9){ Y3js[i]=Y3[i]; i++;} b12=0;b21=0;b22=0;C31=0;C32=0;a31=0;a32=0;i=0;while(i<=9){ b11=b11+2850*Ab[i]*Y0xxx[i]*Y0xxx[i]*0.01; i++;}i=0;while(i<=9){ b12=b12+2850*Ab[i]*Y0xxx[i]*Y2[i]*0.01; i++;}i=0;while(i<=9){ b21=b21+2850*Ab[i]*Y0xxx[i]*Y2[i]*0.01; i++;}i=0;while(i<=9){ b22=b22+2850*Ab[i]*Y2[b]*Y2[i]*0.01; i++;}i=0;while(i<=9){ C31=C31+2850*Ab[i]*Y0xxx[i]*Y3js[i]*0.01; i++;}i=0;while(i<=9){ C32=C32+2850*Ab[i]*Y2[i]*Y3js[i]*0.01; i++;}a31=C31/b11;a32=C32/b22;i=0;while(i<=9){ Y3[i]=Y3js[i]-a31*Y0xxx[i]-a32*Y2[i]; i++;}shoulian=fabs(Y3js[0]-Y3[0]);}q=0; while(9>=q){ one[q]=0; two[q]=0; three[q]=0; four[q]=0; q=q+1;}//這步是給數(shù)組清零,千萬(wàn)不能忘!a=0;b=0;while(9>=a){ b=a; while(9>=b+2){ one[a]=one[a]+(Ab[b+2]*Y3[b+2]*0.01); b=b+1;}a=a+1;}//第一重循環(huán)c=0;d=0;while(9>=c){ d=c; while(9>=d+1){ two[c]=two[c]+(one[d+1]*0.01); d=d+1;}c=c+1;}//第二重循環(huán)e=0;f=0;while(9>=e){ f=0; while(e>=f){ three[e]=three[e]+two[f]*0.01*(1/Ib[f]); f=f+1;}e=e+1;}//第三重循環(huán)g=0;h=0;while(9>=g)
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 【正版授權(quán)】 ISO 7436:2024 EN Fasteners - Slotted set screws with cup point
- 2025年度二手房買賣交易合同(附帶房屋抵押權(quán)解除及貸款還款計(jì)劃)
- 2025年度新能源項(xiàng)目安全生產(chǎn)責(zé)任書(shū)匯編
- 2025年度科技創(chuàng)新項(xiàng)目資金擔(dān)保合同
- 2025年高科技車間承包服務(wù)協(xié)議
- 2025年度社區(qū)配套車位代理銷售服務(wù)合同
- 傳統(tǒng)藝術(shù)與現(xiàn)代美術(shù)課程融合計(jì)劃
- 如何進(jìn)行有效的課堂觀察計(jì)劃
- 提升自我管理能力的策略計(jì)劃
- 關(guān)注員工個(gè)人發(fā)展的措施計(jì)劃
- 部編版人教版語(yǔ)文八年級(jí)下冊(cè)全冊(cè)課件
- 2024年太倉(cāng)高新控股有限公司招聘筆試沖刺題(帶答案解析)
- 人教版七年級(jí)地理下冊(cè)《全冊(cè)完整》
- 2024年02月中央軍委后勤保障部2024年公開(kāi)招考專業(yè)技能崗位文職人員筆試參考題庫(kù)附帶答案詳解
- 小學(xué)數(shù)學(xué)三年級(jí)下冊(cè)第八單元《數(shù)學(xué)廣角-搭配(二)》大單元集體備課整體設(shè)計(jì)
- 10kv高壓送電專項(xiàng)方案
- (高清版)TDT 1031.6-2011 土地復(fù)墾方案編制規(guī)程 第6部分:建設(shè)項(xiàng)目
- 煤炭供應(yīng)鏈管理與協(xié)同創(chuàng)新
- 2024年江蘇省高中學(xué)業(yè)水平測(cè)試生物試卷
- 露天采場(chǎng)危險(xiǎn)有害因素辨識(shí)
- 食品感官評(píng)價(jià)員培訓(xùn)方案
評(píng)論
0/150
提交評(píng)論