【正文】 版社, 2001[17] 馬建敏，肖燦章，李農(nóng)勝. 高頻噪聲場(chǎng)中薄壁結(jié)構(gòu)響應(yīng)的預(yù)估 [ J ].航空學(xué)報(bào), 1992, 13 (5) : 317~321[18] Marlana N. Behnke, Anurag Sharma, Adam Przekop, Stephen A. Rizzi，ThermalAcoustic Analysis of a Metallic Integrated Thermal Protection System Structure, AIAA, 2010[19] R. Vaicaitis, ACOUSTIC FATIGUEA MONTE CARLO APPROACH, AIAA, 1987[20] R. Vaicaitis, MONTE CARLO SIMULATION FOR APPLICATION TO AEROSPACE STRUCTURES, DTIC, 1990[21] Adam Przekop, Stephen A. Rizzi，An Investigation of HighCycle Fatigue Models for Metallic Structures Exhibiting SnapThrough Response, AIAA, 2007[22] M. Shinozuka, Y. K. Wen, Monte Carlo Solution of Nonlinear Vibrations, AIAA, 1972[23] Roussos L A, Heitman K E, Rucker C E. Predicted and Measured Strain Response of Rectangular Panels Due to Acous2tic Loading [ J ], AIAA2 862 1931: 1—10 .[24] Prasad C B, Chuh Mei Multi p le2 mode large deflection random response of beams with nonlinear damping subjected t o acoustic excitation[R ]. AIAA2 87—2712.[25] James Locke, A Finite Element Formulation for The Large Deflection Random Response of Thermally Buckled Beans AIAA 891354CP[26] Clarkson, B. L., Review of sonic fatigue technology NASA CR4587, April NASA Langley Research Center, Hampton, VA USA, 1994[27] K. AbdelMotagaly. B. Duan, Nonlinear Response of Composite Panels under Combined Acoustic Excitation and Aerodynamic Pressure NASA/CP199920913 6/PT[28] M. Shinozuka, Y. K. Wen, Monte Carlo Solution of Nonlinear Vibrations, AIAA, 1972[29] Stephen A. Rizz, Experimental Research Activities in Dynamic Response and Sonic Fatigue of Hypersonic Vehicle Structure at NASA Langley Research center AIAA930383[30] Jack D. Leatherwood, Acoustic Testing of High Temperature Panels AIAA903939CP…………致 謝本文是在指導(dǎo)老師沙云東教授的悉心指導(dǎo)下完成的，從論文的選題立意、理論分析、具體工作以及最后論文審稿、定稿等各方面，沙老師都給予了細(xì)致入微的指導(dǎo)，促使我順利完成畢業(yè)設(shè)計(jì)內(nèi)容及論文的撰寫(xiě)工作。沙老師以其淵博的理論知識(shí)，嚴(yán)謹(jǐn)?shù)闹螌W(xué)態(tài)度，對(duì)學(xué)問(wèn)孜孜不倦的探索精神和對(duì)學(xué)生的循循善誘的指導(dǎo)，深深的打動(dòng)了我，從中學(xué)到了很多知識(shí)和為人處事的道理，這將使我終身受益。值此論文完成之際，謹(jǐn)向沙老師表示我最崇高的敬意和最真摯的感謝。在課題的完成過(guò)程中，感謝我的學(xué)長(zhǎng)高志軍給予的無(wú)私幫助，還有同學(xué)李朋賓、梁桄瑋、郭峰給予我的幫助和建議。在此，向他們表示誠(chéng)摯的感謝。對(duì)所有曾給予我關(guān)心和幫助的老師、同學(xué)及朋友表示感謝。在此論文完成之際，我還要深深感謝我的親人，正是他們的支持和鼓勵(lì)才讓我有了這次學(xué)習(xí)和研究的機(jī)會(huì)。祝愿他們平平安安、身體健康。附錄Ⅰ程序代碼1批量計(jì)算跳變響應(yīng)程序代碼*DIM,TEMP_FILE,CHAR,7,1,1,*DIM,TEMP_LOAD,ARRAY,7,1,1,*DIM,SP_LOAD,CHAR,7,1,1*SET,TEMP_FILE(1),39。39。,39。39。,39。39。,39。39。,39。39。,39。39。,39。39。*SET,TEMP_LOAD(1),0,10,20,30,40,50,60*SET,SP_LOAD(1),39。15539。,39。16039。,39。16539。*DIM,TP,TABLE,1999,1,TIME,TEMP*DIM,SP,TABLE,18384,1,TIME,PRESSURE,*DO,I,1,7PARSAV,ALL,MEMO,/CLEAR,STARTPARRES,NEW,MEMO,*TREAD,TP,TEMP_FILE(I),39。TXT39。,39。DATA\39。,*TREAD,SP,SP%SP_LOAD(2)%,39。TXT39。,39。DATA\39。,/CONFIG,NRES,100000/FILNAME,PreT%TEMP_LOAD(I)%SPL%SP_LOAD(2)%F1500/PREP7ET,1,SHELL181 R,1, MP,DENS,1,2700 MP,EX,1,73E9MP,PRXY,1, MP,ALPX,1,MP,REFT,1,0 RECTNG,0,0, LESIZE,ALL, AMESH,ALL LSEL,S,LINE,1,4,1 NSLL,S,1D,ALL,ALL ALLSEL ANTYPE,4BF,ALL,TEMP,%TP%SFE,ALL,PRES,%SP%NLGEOM,ONSSTIF,OFFKBC,1TIME,DELTIM,1E4AUTOTS,OFFOUTRES,ERASEOUTRES,NSOL,ALLOUTRES,STRS,ALLLSWRITE,1BF,ALL,TEMP,%TEMP_LOAD(I)%SFE,ALL,PRES,%SP%NLGEOM,ONSSTIF,OFFKBC,1TIME,DELTIM,1E4AUTOTS,OFFOUTRES,ERASEOUTRES,NSOL,ALLOUTRES,STRS,ALLLSWRITE,2/SOLLSSOLVE,1,2,1SAVEFINISH*ENDDO2[ANSYS]固支板屈曲后模態(tài)頻率/CONFIG,NRES,100000/VIEW,1,2,3/FILNAME,BeamBuckling/TITLE,EigenBuckling/PREP7ET,1,SHELL181R,1,MP,DENS,1,2700MP,EX,1,73E9MP,PRXY,1,MP,ALPX,1,MP,REFT,1,0RECTNG,0,0,LESIZE,ALL,AMESH,ALLLSEL,S,LINE,2,4,2NSLL,S,1D,ALL,ALLALLSELBF,ALL,TEMP,1/SOLUANTYPE,STATICPSTRES,ONSOLVEFINISH/SOLUANTYPE,BUCKLEPSTRES,ONBUCOPT,LANB,10SOLVEFINISH3功率譜密度生成程序clc。clear。format long [Name Path]=uigetfile(39。*.txt39。,39。Select File39。,39。File Name39。)。InputPath=[Path Name]。 w=load(InputPath)。Fs=1/(w(2,1)w(1,1))。[Pxx,f]=pwelch(w(:,2),[],[],[],Fs)。figure(1)semilogy(f,Pxx,39。color39。,39。red39。,39。linewidth39。,)xlabel(39。Frequency [Hz]39。,39。fontsize39。,11)。ylabel(39。UZ Power Spectral Density [m^2/Hz]39。,39。fontsize39。,11)。 axis([0 1500 1E6 1E16])。 % axis([0 1e3 1e15 1e5])。 % set(gca,39。position39。,[ ],39。fontsize39。,10)。% print(i,39。djpeg39。,[FileName{i},39。_DisPSD39。])。[Pxx,f]=pwelch(w(:,3),[],[],[],Fs)。figure(2)semilogy(f,Pxx,39。color39。,39。red39。,39。linewidth39。,)xlabel(39。Frequency [Hz]39。,39。fontsize39。,11)。ylabel(39。SX Power Spectral Density [m^2/Hz]39。,39。fontsize39。,11)。 %axis([0 2000 1E+4 1E+10])。 [Pxx,f]=pwelch(w(:,4),[],[],[],Fs)。figure(3)semilogy(f,Pxx,39。color39。,39。red39。,39。linewidth39。,)xlabel(39。Frequency [Hz]39。,39。fontsize39。,11)。ylabel(39。SY Power Spectral Density [m^2/Hz]39。,39。fontsize39。,11)。 %axis([0 2000 1E+4 1E+10])。 [Pxx,f]=pwelch(w(:,5),[],[],[],Fs)。figure(4)semilogy(f,Pxx,39。color39。,39。red39。,39。linewidth39。,)xlabel(39。Frequency [Hz]39。,39。fontsize39。,11)。ylabel(39。SXY Power Spectral Density [m^2/Hz]39。,39。fontsize39。,11)。 %axis([0 2000 1E+4 1E+10])。 [Pxx,f]=pwelch(w(:,6),[],[],[],Fs)。figure(5)semilogy(f,Pxx,39。color39。,39。red39。,39。linewidth39。,)xlabel(39。Frequency [Hz]39。,39。fontsize39。,11)。ylabel(39。SEQV Power Spectral Density [m^2/Hz]39。,39。fontsize39。,11)。 axis([0 1500 1E+14 1E+6])。附錄Ⅱ位移及應(yīng)力特性分布圖肋板160dB聲載荷下