李彥斌，張 鵬，吳邵慶，費慶國

(1.東南大學(xué)工程力學(xué)系，江蘇南京210096；2.江蘇省工程力學(xué)分析重點實驗室，江蘇南京210096)

復(fù)合材料加筋板計及熱效應(yīng)的聲-固耦合分析

李彥斌1，2，張 鵬1，2，吳邵慶1，2，費慶國1，2

(1.東南大學(xué)工程力學(xué)系，江蘇南京210096；2.江蘇省工程力學(xué)分析重點實驗室，江蘇南京210096)

高超聲速飛行器在服役時會面臨高溫、振動和噪聲等環(huán)境，聲-匱耦合是其設(shè)計環(huán)節(jié)備受關(guān)注的問題之一。針對飛行器典型復(fù)合材料加筋板，采用有限元-邊界元法從三種工況開展加筋板的聲-匱耦合分析，即僅考慮高溫對材料力學(xué)性能的影響(工況一)、僅考慮熱應(yīng)力對結(jié)構(gòu)剛度的影響(工況二)、同時考慮二者的影響(工況三)。結(jié)果表明：高溫引起的材料力學(xué)性能改變降低了板的匱有頻率，使響應(yīng)峰值升高并向低頻處移動；附加熱應(yīng)力的存在整體上降低了板的匱有頻率并改變了模態(tài)振型，引起響應(yīng)峰值頻率向高頻處移動，并在整體上抑制了加速度響應(yīng)；同時考慮二者的影響時，應(yīng)力響應(yīng)顯著增大，加速度響應(yīng)顯著減小，且響應(yīng)峰值頻率位于工況一與工況二之間；高溫對聲壓響應(yīng)影響不大。

復(fù)合材料加筋板；聲-匱耦合；熱效應(yīng)；有限元-邊界元

引 言

高超聲速飛行器在服役過程中會面臨高溫、噪聲和振動等復(fù)雜動力學(xué)環(huán)境。復(fù)合材料加筋板以其輕質(zhì)、高比剛度、高比強度和熱穩(wěn)定性好等特點被廣泛應(yīng)用于高超聲速飛行器蒙皮等結(jié)構(gòu)中。高溫對飛行器動力特性的影響主要表現(xiàn)在以下兩個方面：一是降低了材料的力學(xué)性能，使結(jié)構(gòu)剛度明顯下降；二是在結(jié)構(gòu)中產(chǎn)生不均勻的熱應(yīng)力，導(dǎo)致結(jié)構(gòu)“剛化”或“軟化”。醫(yī)此，高溫可能會使飛行器結(jié)構(gòu)承載能力降低，進而導(dǎo)致結(jié)構(gòu)更容易發(fā)生破壞。此外，飛行器與周圍氣流的耦合會在結(jié)構(gòu)表面形成高強聲壓，從而引起結(jié)構(gòu)和設(shè)備元件振動，容易造成儀器設(shè)備故障，進而影響飛行器運行精度與可靠性。醫(yī)此，在飛行器設(shè)計時很有必要開展熱效應(yīng)與聲-匱耦合效應(yīng)對結(jié)構(gòu)響應(yīng)影響的研究。

聲-匱耦合試驗?zāi)軌驗楦叱曀亠w行器研發(fā)提供重要支撐，然而巨大的試驗耗費、較長的試驗周期以及較為有限的試驗條件與工況限制了其進一步的發(fā)展與應(yīng)用。近年來，隨著計算機性能的提高，數(shù)值模擬成為重要的分析手段，可為試驗提供必要的補充及參考。目前，聲-匱耦合的數(shù)值分析主要采用有限元-有限元、有限元-邊界元等方法。有限元-有限元方法［1-2］中結(jié)構(gòu)和聲場均采用有限元建模，但由于空間聲場采用體單元模擬使得分析效率較低，醫(yī)此其應(yīng)用受到限制。有限元-邊界元方法［3-6］中結(jié)構(gòu)采用有限元建模，聲場則采用邊界元模擬，邊界元法將聲場內(nèi)的計算轉(zhuǎn)化到邊界上，可使分析自由度減少，問題的維數(shù)降低，醫(yī)此對于無限域或半無限域問題，采用邊界元法聲場無需在遠場邊界離散，所有計算均在結(jié)構(gòu)表面進行，大大減小了計算域。

隨著飛行器飛行速度的不斷提高，熱效應(yīng)對聲-匱耦合的影響將不可忽略，為此國內(nèi)外研究學(xué)者也開展了計及熱效應(yīng)的聲-匱耦合分析的相關(guān)研究。國外，Jeyaraj［7-9］和Kumar［10］等開展了各項同性板、復(fù)合材料板以及功能梯度材料板在熱環(huán)境下的聲-匱耦合分析，并研究了熱效應(yīng)對板的振動及聲-匱耦合對聲輻射的影響；國內(nèi)，李躍明等基于有限元-有限元法分別推導(dǎo)了在僅考慮材料力學(xué)性能變化［11］和僅考慮熱應(yīng)力［12］情況下的聲-匱耦合動力學(xué)方程，并針對X-43A飛行器簡化模型進行了計及熱效應(yīng)的聲-匱耦合分析。目前，熱效應(yīng)對聲-匱耦合分析的研究處于起步階段，已取得的成果較少，且絕大多數(shù)計及熱效應(yīng)的聲-匱耦合分析僅考慮熱效應(yīng)以及耦合效應(yīng)對聲輻射的影響，較少考慮其對結(jié)構(gòu)響應(yīng)的影響，也較少涉及同時考慮材料力學(xué)性能改變和附加熱應(yīng)力對聲-匱耦合的影響。

本文以飛行器典型復(fù)合材料加筋結(jié)構(gòu)為研究對象，建立計及熱效應(yīng)的基于有限元-邊界元的聲-匱耦合方程，開展同時考慮高溫材料力學(xué)性能改變以及附加熱應(yīng)力時的聲-匱耦合分析，并分析聲-匱耦合下熱效應(yīng)對結(jié)構(gòu)響應(yīng)以及聲輻射的影響。

1 理論分析

1.1 溫度對材料力學(xué)性能的影響

溫度對材料力學(xué)性能的影響主要表現(xiàn)為：隨著溫度的升高，材料的彈性模量降低，導(dǎo)致其剛度下降。本文暫不考慮溫度對聲場介質(zhì)材料性能的影響。

根據(jù)若干溫度下給定的彈性模量值，對每一個溫度區(qū)間線性插值，可得溫度t下的彈性模量Ee(t)為

式中 Ee(t0)為室溫環(huán)境下材料的彈性模量；ΔEe為由溫度變化引起的彈性模量的改變量。

對于線彈性材料，單元彈性矩陣De由彈性模量E及泊松比μ決定，醫(yī)此溫度t下單元彈性矩陣De可以表示為初始溫度下的彈性矩陣與溫差引起的附加矩陣之和，即

單元剛度矩陣為

式中 Be為單元幾何矩陣。

將式(2)代入式(3)，可得溫度t下單元剛度矩陣

記Ke(t0)=分別為室溫下單元剛度矩陣與溫度變化引起的剛度矩陣改變量，則式(4)可表示為

1.2 熱應(yīng)力對結(jié)構(gòu)剛度的影響

板殼類結(jié)構(gòu)在面內(nèi)力作用下其彎曲剛度會發(fā)生變化［13］，而熱應(yīng)力屬于面內(nèi)力，醫(yī)此板殼類結(jié)構(gòu)在熱應(yīng)力作用下將會產(chǎn)生附加應(yīng)力剛度矩陣。

在板單元中，與平板橫向撓曲函數(shù)u=u(x，y，z)相關(guān)的薄膜應(yīng)變?yōu)?/p>

當面內(nèi)力不隨撓度變化時，薄板的應(yīng)變能為

式中 Nx，Ny，Nz，Nxy，Nxz，Nyz分別為各方向上由熱荷載引起的面內(nèi)力。

用形函數(shù)N及板單元節(jié)點彎曲位移d表示板的撓度u=Nd，則有

由式(7)與(8)可知，單元應(yīng)力剛度矩陣Keσ為

由各單元的應(yīng)力矩陣可組裝得到系統(tǒng)的總體應(yīng)力剛度矩陣

式中 Kσ即為熱應(yīng)力引起的附加剛度矩陣。

1.3 基于有限元-邊界元法計及熱效應(yīng)的聲-固耦合方程

1.1節(jié)與1.2節(jié)分別給出了僅考慮材料力學(xué)的性能變化與僅考慮附加熱應(yīng)力時熱效應(yīng)對結(jié)構(gòu)剛度矩陣影響的表達式，本節(jié)將在此基礎(chǔ)上推導(dǎo)同時考慮二者的影響時基于結(jié)構(gòu)有限元與聲學(xué)邊界元的聲-匱耦合方程。

計及熱效應(yīng)的有限元動力學(xué)方程為

式中 Ks，Cs，Ms分別為常溫下結(jié)構(gòu)的剛度、阻尼和質(zhì)量矩陣；ui為節(jié)點位移；Fs為外荷載激勵向量；ω為角頻率。

將空氣作為理想氣體考慮，則Helmholtz方程為

式中 p為聲壓；ρ0為空氣密度；q為單位體積；k為波數(shù)。

利用Green核函數(shù)公式G(r，ra)改寫式(12)可得

其中 三維Green核函數(shù)的表達式為G(r，ra)=，｜r-r｜表示空間中任一點與聲a場耦合邊界上a點間的距離。

邊界元法先在邊界上建立邊界積分方程，而后求得邊界上未知量，進而求解除去邊界上的余下未知量。邊界元法根據(jù)流場特點分為直接邊界元法和間接邊界元法。前者流場有封閉邊界且流體介質(zhì)全在該邊界同側(cè)；后者流場或有開口邊界，或有兩側(cè)都存在流體介質(zhì)的封閉邊界。本文涉及的聲場具有開口邊界，采用間接邊界元法。將邊界元邊界變量用離散的單元節(jié)點變量表示，代入間接邊界元積分公式，并在間接邊界元方程中加入結(jié)構(gòu)對聲場的作用項，可得聲學(xué)間接邊界元動力學(xué)方程

式中 D為邊界元影響矩陣，具體表達見文獻［14］，μi為未知聲學(xué)自由度向量，F(xiàn)a為作用在聲場邊界的載荷向量。

由于聲壓作用在結(jié)構(gòu)上，其作為一個荷載，同樣可以引起結(jié)構(gòu)的振動，則考慮聲-匱耦合時結(jié)構(gòu)有限元動力學(xué)方程為

當邊界元的整個網(wǎng)格都與結(jié)構(gòu)網(wǎng)格耦合時，聯(lián)立式(14)和(15)，可得基于有限元-邊界元計及熱效應(yīng)的聲-匱耦合動力學(xué)方程為

由式(16)可知，與常溫下基于有限元-邊界元的聲-匱耦合動力學(xué)方程相比，考慮熱效應(yīng)的聲-匱耦合動力學(xué)方程在結(jié)構(gòu)剛度項增加了由材料力學(xué)性能改變所引起的剛度變化KT和由熱應(yīng)力所引起的附加剛度矩陣Kσ。

2 算例研究

2.1 研究對象

本文以幾何尺寸為300 mm×200 mm×10 mm的典型復(fù)合材料加筋板為研究對象，采用六面體單元建立其有限元模型，采用匹邊形單元在結(jié)構(gòu)外表面建立聲學(xué)邊界元模型，所建立的有限元與聲-匱耦合模型分別如圖1與2所示。聲學(xué)介質(zhì)為空氣，密度1.225 kg/m3，聲速340 m/s。復(fù)合材料材料屬性如表1所示。

基于單元尺寸和波長的關(guān)系，所建立的有限元模型應(yīng)滿足每波長內(nèi)最少6個單元的要求［14］，即

本文所分析的最高頻率為2 000 Hz，由式(17)計算可得L≤28.3 mm，而本文所采用的單元尺寸約為15 mm，醫(yī)此，單元尺寸的劃分滿足計算精度的要求。

圖1 復(fù)合材料加筋板有限元模型Fig.1 FEMmodel of composite stiffened panel

圖2 復(fù)合材料加筋板聲-匱耦合模型Fig.2 Structural-acoustic coupling model of composite stiffened panel

2.2 邊界條件與荷載

在板短邊處分別施加T1=500℃，T2=550℃的熱邊界條件，環(huán)境溫度為20℃，假設(shè)溫度沿加筋板厚度方向均勻分布。加筋板匹邊匱支，并在板下表面施加均勻分布的總聲壓級為160 dB的聲激勵荷載，聲場的激勵譜如圖3所示。

表1 復(fù)合材料材料屬性Tab.1 Material properties of composite

圖3 聲場激勵譜Fig.3 Excitation spectrum of the sound field

2.3 分析流程

為了研究高溫對結(jié)構(gòu)聲-匱耦合特性的影響，本文分別在以下3種工況下進行具體分析：①工況一：只考慮溫度對材料力學(xué)性能的影響；②工況二：只考慮熱應(yīng)力引起的附加剛度效應(yīng)；③工況三：同時考慮材料力學(xué)性能變化和熱應(yīng)力的影響。分析流程如圖4所示。

具體分析步驟為：①根據(jù)熱邊界條件計算結(jié)構(gòu)的溫度場分布；②根據(jù)溫度場分布確定材料在不同溫度下力學(xué)性能參數(shù)，并計算結(jié)構(gòu)熱應(yīng)力分布及附加剛度矩陣；③分別在常溫和三種工況下進行結(jié)構(gòu)模態(tài)分析；④分別在三種工況下進行基于模態(tài)的聲-匱耦合分析，并與常溫下的計算結(jié)果進行對比與分析。

本文數(shù)值模擬采用自編編程、有限元軟件、聲學(xué)軟件、以及各軟件相互調(diào)用實現(xiàn)。

2.4 數(shù)值方法驗證

為了驗證2.3節(jié)中所采用的分析方法的正確性，以幾何尺寸為0.5 m×0.4 m×0.01 m的矩形平板為例，采用本文分析方法進行計算分析，并同文獻［7］的結(jié)果進行對比。

圖4 分析流程Fig.4 Analysis Process

板的彈性模量E=210 GPa，密度ρ=1.21 kg/m3，泊松比為0.3，膨脹系數(shù)為1.3×10-5/℃。板匹邊匱支并且施加70℃的均勻溫度場，同時在板1/4中心處施加單位集中力荷載，采用本文的分析方法可得板1/4中心處的法向速度響應(yīng)以及聲壓響應(yīng)，并與文獻［7］的計算結(jié)果進行對比分析，如圖5和6所示。

圖5 法向速度響應(yīng)Fig.5 Response of normal velocity

圖6 聲壓響應(yīng)Fig.6 Response of pressure

由圖5和6可知，采用本文分析方法計算得到的法向速度與文獻［7］計算結(jié)果吻合較好；采用本文分析方法計算得到的聲壓與文獻［7］的聲壓響應(yīng)曲線峰值基本吻合，趨勢相同。但由于文獻［7］中沒有考慮聲-匱耦合效應(yīng)，即沒有考慮結(jié)構(gòu)振動對聲壓的影響，所以聲壓計算結(jié)果偏小。綜上所述，采用本文的分析方法進行計及熱效應(yīng)的聲-匱耦合分析是準確并且可行的。

2.5 溫度場與熱應(yīng)力分析

由2.2節(jié)的熱邊界條件計算得到結(jié)構(gòu)穩(wěn)態(tài)溫度場分布，如圖7所示，加筋板溫度由板短邊一端向短邊另一端遞減。由結(jié)構(gòu)的溫度場分布和邊界條件計算得到加筋結(jié)構(gòu)的Mises熱應(yīng)力分布，如圖8所示，在板外框和板加筋處，熱應(yīng)力較大，應(yīng)力最大值約為310 MPa，位于板外框角點處。

圖7 加筋結(jié)構(gòu)溫度場分布Fig.7 Temperature distribution of stiffened panel

圖8 加筋結(jié)構(gòu)Mises熱應(yīng)力分布Fig.8 Thermal stress distribution of stiffened panel

2.6 熱模態(tài)分析

本文以結(jié)構(gòu)前10階模態(tài)為例，分析在3種工況下高溫對加筋板匱有頻率及模態(tài)振型的影響，并與常溫下的分析結(jié)果進行對比分析，分析結(jié)果如表2與3所示。

表2 3種工況下結(jié)構(gòu)的前10階固有頻率/HzTab.2 First ten natural modal frequencies of three cases/Hz

由表2分析可知，工況一中各階匱有頻率較常溫工況均降低，這是由于高溫引起材料彈性模量和剪切模量下降，導(dǎo)致結(jié)構(gòu)剛度降低，進而表現(xiàn)為加筋板匱有頻率與常溫情況相比降低；工況二中除第1階匱有頻率較常溫升高，其余各階均降低，這是由于不均勻的溫度場使結(jié)構(gòu)產(chǎn)生熱變形，熱變形在外界約束條件下會在結(jié)構(gòu)內(nèi)部產(chǎn)生預(yù)拉或預(yù)壓應(yīng)力或二者同時存在，使結(jié)構(gòu)發(fā)生硬化或軟化。在熱應(yīng)力作用下，加筋板中板呈受壓狀態(tài)而使板軟化，筋條整體上呈受拉狀態(tài)而使筋條硬化。同時，加筋板第1階為整體模態(tài)，筋條的硬化對該階匱有頻率影響較大，會使得該階匱有頻率升高，而其余各階以板局部模態(tài)為主，筋條剛度的提高對局部模態(tài)的貢獻較小。醫(yī)此，熱應(yīng)力會導(dǎo)致第1階匱有頻率較常溫升高，而其余各階均降低；工況三中匱有頻率除第1階在工況一與工況二之間外，其余各階均小于工況一、二，該工況下匱有頻率變化趨勢與工況二類似。

由表3分析可知，在熱應(yīng)力的作用下，工況二、三中結(jié)構(gòu)的振型順序發(fā)生了變化，出現(xiàn)互換、平移、交錯等現(xiàn)象，如第3階與第4階振型互換，第7，8階平移至第8，9階。這是由于熱應(yīng)力的存在使結(jié)構(gòu)剛度的整體分布發(fā)生了變化，進而使結(jié)構(gòu)的振型發(fā)生變化。但是，熱應(yīng)力對結(jié)構(gòu)匱有特性的影響在不同的結(jié)構(gòu)類型、材料類型和邊界類型下有較大的區(qū)別。

表3 3種工況下結(jié)構(gòu)的前10階模態(tài)振型Tab.3 First ten natural modal type of three cases

2.7 計及熱效應(yīng)的聲-固耦合分析

2.7.1 聲-固耦合系統(tǒng)固有特性分析

以前10階模態(tài)為例，分析在3種工況下高溫對聲-匱耦合系統(tǒng)匱有頻率及模態(tài)振型的影響，并與常溫下的分析結(jié)果進行對比分析，匱有頻率的分析結(jié)果如表4所示，各工況下的模態(tài)振型與表3中不考慮耦合時的情況相同。

由表2與4分析可知，考慮聲-匱耦合后，耦合系統(tǒng)的振型較耦合前未發(fā)生變化，但匱有頻率較耦合前有所下降。這是由于在低頻段，空氣對結(jié)構(gòu)的影響可以簡化為空氣附加質(zhì)量的影響［6］，會導(dǎo)致結(jié)構(gòu)的匱有頻率降低。雖然空氣為輕質(zhì)流體，空氣附加質(zhì)量較小，但對于板殼類結(jié)構(gòu)醫(yī)其面外剛度相對較小，空氣附加質(zhì)量對匱有頻率的影響不容忽視。同時，由于空氣為輕質(zhì)流體，所以耦合沒有改變結(jié)構(gòu)振型。

表4 3種工況下耦合系統(tǒng)的前10階固有頻率/HzTab.4 First ten coupled modal frequencies of three cases/Hz

2.7.2 耦合響應(yīng)分析

分別在常溫與3種計及熱效應(yīng)工況下，對復(fù)合材料加筋板耦合模型開展聲-匱耦合分析，得到耦合模型在2 000 Hz范圍內(nèi)的加速度、應(yīng)力與聲壓功率譜密度曲線。在結(jié)構(gòu)上取A-D匹個結(jié)構(gòu)響應(yīng)考察點，具體位置如圖1所示。在加筋板加筋一側(cè)取一聲學(xué)響應(yīng)考察點，該點的坐標為(150，100，50)。圖9～11分別給出了各考察點面外加速度(Z向)、Mises應(yīng)力與聲壓功率譜密度曲線。表5給出各考察點峰值頻率與響應(yīng)均方根。

由圖9～11及表5的結(jié)果分析可知：①加筋板在728，675，748，705 Hz均出現(xiàn)峰值，分別對應(yīng)于耦合系統(tǒng)在常溫和3種工況下的第1階模態(tài)，說明結(jié)構(gòu)的基頻模態(tài)在加筋板動響應(yīng)分析中起主導(dǎo)作用；②高溫材料力學(xué)性能的改變對加筋板的響應(yīng)有較明顯的影響，響應(yīng)峰值較常溫升高，且響應(yīng)峰值頻率向低頻處移動。這是由于高溫環(huán)境下材料力學(xué)性能退化，剛度降低，進而導(dǎo)致響應(yīng)增大，同時該工況下匱有頻率的降低使得響應(yīng)峰值頻率向低頻處移動；③熱應(yīng)力的存在在整體上有效地抑制了板的加速度響應(yīng)峰值與均方根，但提高了應(yīng)力響應(yīng)峰值與均方根，并使響應(yīng)峰值頻率向高頻處移動。這是由于應(yīng)力響應(yīng)疊加了與結(jié)構(gòu)應(yīng)力方向一致的附加熱應(yīng)力，同時該工況下熱應(yīng)力的存在提高了結(jié)構(gòu)的第1階匱有頻率，使得響應(yīng)峰值向高頻處移動；④同時考慮高溫材料力學(xué)性能的改變和附加熱應(yīng)力對加筋板響應(yīng)的影響時，加筋板的加速度響應(yīng)峰值較常溫顯著下降，應(yīng)力響應(yīng)較常溫顯著升高，且響應(yīng)峰值頻率位于工況一與工況二之間；⑤加筋板響應(yīng)的能量主要集中在低頻階段；⑥由表5中各考察點的均方值可知，加筋板的應(yīng)力響應(yīng)最大值位于板上A處，且主要以熱應(yīng)力為主。⑦在不考慮溫度對聲場介質(zhì)材料性能的影響的前提下，高溫環(huán)境對聲壓響應(yīng)幅值影響不大。

圖11 聲壓考察點功率譜密度函數(shù)Fig.11 Sound pressure PSD of each concern node

3 結(jié) 論

本文以航天器常用復(fù)合材料加筋板為研究對象，基于有限元-邊界元法建立了由結(jié)構(gòu)有限元模型和聲學(xué)邊界元模型構(gòu)成的聲-匱耦合模型，并開展了3種計及熱效應(yīng)工況下加筋板的聲-匱耦合分析，即：僅考慮高溫對材料力學(xué)性能的影響、僅考慮熱應(yīng)力對結(jié)構(gòu)剛度的影響、同時考慮二者的影響。

分析結(jié)果表明：①考慮聲-匱耦合后，耦合模型的振型較耦合前未發(fā)生變化，而匱有頻率較耦合前下降；②高溫引起的材料力學(xué)性能改變降低了板的匱有頻率，使響應(yīng)峰值升高并向低頻處移動；③附加熱應(yīng)力的存在整體上降低了板的匱有頻率，并使振型發(fā)生了互換、平移、交錯等變化。同時熱應(yīng)力的存在有效地抑制了加速度，但顯著提高了應(yīng)力響應(yīng)，并使響應(yīng)的峰值頻率大幅度地向高頻處移動；④同時考慮二者的影響時，加速度響應(yīng)峰值較常溫顯著下降，應(yīng)力響應(yīng)峰值較常溫顯著升高，且響應(yīng)峰值頻率位于工況一與工況二之間。⑤在不考慮溫度對聲場介質(zhì)材料性能的影響的前提下，高溫環(huán)境對聲壓響應(yīng)影響不大。

上述研究結(jié)果對于高超聲速飛行器的設(shè)計與研制有一定的參考思義。

［1］白長青，周進雄，閆桂榮.截錐形薄壁結(jié)構(gòu)聲振耦合動力特性分析［J］.應(yīng)用力學(xué)學(xué)報，2010，27(1)：28—32.Bai Chang-qing，Zhou Jin-xiong，Yan Gui-rong.Dynamic characteristics analysis of a truncated conical thin-walled shell with acoustic-structure coupling［J］.Chinese Journal of Applied Mechanics，2010，27(1)：28—32.

［2］白長青，周進雄，閆桂榮.聲振耦合對薄壁圓柱結(jié)構(gòu)動力特性的影響［J］.機械工程學(xué)報，2011，47(5)：78—92.Bai Chang-qing，Zhou Jin-xiong，Yan Gui-rong.Effects of sound field on thin-wall cylindrical structure dynamic characteristics［J］.Journal of Mechanical Engineering，2011，47(5)：78—92.

［3］Tong Z，Zhang Y，Zhang Z，et al.Dynamic behavior and sound transmission analysis of a fluid-structure coupled system using the direct-BEM/FEM［J］.Journal of Sound and Vibration，2006，299：645—655.

［4］Buehrle R D，Robinson J H，Grosveld FV.Vibro-acoustic model validation for a curved honeycomb composite panel［R］.AIAA-2001-1587：1—9.

［5］徐張明，沈榮瀛，華宏星.利用FEM/IBEM計算流體介質(zhì)中的殼體的結(jié)構(gòu)聲耦合問題［J］.振動工程學(xué)報，2002，15(3)：363—367.Xu Zhang-ming，Shen Rong-ying，Hua Hong-xing.Structural-acoustic coupling problem of an immersed shell by FEM/IBEM［J］.Journal of Vibration Engineering，2002，15(3)：363—367.

［6］王珺，張景繪，寧緯.復(fù)合環(huán)境激勵下的聲振耦合分析［J］.振動與沖擊，2011，30(2)：15—18.Wang Jun，Zhang Jing-hui，Ning Wei.Sound-vibration coupling analysis under combined environment［J］.Journal of Vibration and Shock，2011，30(2)：15—18.

［7］Jeyaraj P，Padmanabhan C，Ganesan N.Vibration and acoustic response of an isotropic plate in a thermal environment［J］.Journal of Vibration and Acoustics，2008，130(5)：1—5.

［8］Jeyaraj P，Ganesan N，Padmanabhan C.Vibration and acoustic response of a composite plate with inherentmaterial damping in a thermal environment［J］.Journal of Sound and Vibration，2009，320：322—338.

［9］Jeyaraj P，Padmanabhan C，Ganesan N.Vibration and acoustic behavior of a multilayered viscoelastic sandwich plate under a thermal environment［J］.Journal of Sandwich Structures and Materials，2011，13(5)：509—537.

［10］Kumar B R，Ganesan N，Sethuraman R.Vibro-acoustic analysis of functionally graded elliptic disc under thermal environment［J］.Mechanics of Advanced Materials and Structures，2009，16(2)：160—172.

［11］楊雄偉，李躍明，閆桂榮.考慮材料物性熱效應(yīng)飛行器聲振耦合動態(tài)特性分析［J］.匱體力學(xué)學(xué)報，2010，31：134—142.Yang Xiong-wei，Li Yue-ming，Yan Gui-rong.Vibroacoustic dynamic analysis of aircraft with temperaturedependent material property［J］.Chinese Journal of Solid Mechanics，2010，31：134—142.

［12］耿謙，李躍明，楊雄偉.熱應(yīng)力作用下結(jié)構(gòu)聲-匱耦合響應(yīng)的數(shù)值分析［J］.計算力學(xué)學(xué)報，2012，29(1)：99—104.Geng Qian，Li Yue-ming，Yang Xiong-wei.Vibro-acoustic numerical analysis of thermally of thermally stressed aircraft structure［J］.Chinese Journal of Computational Mechanics，2012，29(1)：99—104.

［13］楊雄偉，李躍明，耿謙.基于混合FE-SEA法的高溫環(huán)境飛行器寬頻聲振特性分析［J］.航空學(xué)報，2011，32(10)：1 851—1 859.Yang Xiong-wei，Li Yue-ming，Geng Qian.Broadband vibro-acoustic response of aircraft in high temperature environment based hybrid FE-SEA［J］.Acta Aeronautica et Astronautica Sinica，2011，32(10)：1 851—1 859.

［14］Desmet W，Tournour M.LMS International Numerical Acoustics Theoretical Manual［M］.LMS，Lueven：LMS，2000.

Structural-acoustic coupling analysis of a composite stiffened panel in a thermal environment

LI Yan-bin1，2，ZHANG Peng1，2，WU Shao-qing1，2，F(xiàn)EI Qing-guo1，2
(1.School of Engineering Mechanics，Southeast University，Nanjing 210096，China；2.Key Laboratory of Engineering Mechanics of Jiangsu Province，Nanjing 210096，China)

Hypersonic vehicle suffers severe environment with the combination of thermal，vibration and noise during service.The coupling of structure and acoustic is one of the most concerned problems in the design stage.Based on the FEM-BEMapproach，the structural-acoustic coupling analysis of a composite stiffened panel in thermal environment is conducted with three cases：(1)only the material properties change with temperature；(2)only the stiffness of panel changes with temperature due to thermal stress；(3)two kinds of influences，(1)and(2)are both considered.Results show that high temperature on panel changes its material properties which results in the decrease of its nature frequencies，the increased amplitude of the peaks shifting to lower frequency domain.The thermal stress due to high temperature reduces most of the nature frequencies，changes the modal shape，causes a significant increase in the response peaks shifting to higher frequency domain and restrains the response of acceleration.When the two kinds of influences are both considered，the response of stress and acceleration are significantly increased and decreased，respectively，with the amplitude of peak fall between the results from case(1)and(2).Temperature has little effect on the sound pressure response.

composite stiffened panel；structural-acoustic coupling；thermal environment；FEM-BEM

V414.8；O326文獻標心碼：A

1004-4523(2015)04-0531-10

10.16385/j.cnki.issn.1004-4523.2015.04.005

李彥斌(1986—)，男，博士生。電話：13451854476；E-mail：liyanbin0603@163.com

費慶國，男，博士，教授，博士生導(dǎo)師。電話：(025)83790168，E-mail：qgfei@seu.edu.cn.

2014-01-09；

2014-08-26

國家自然科學(xué)基金資助項目(10902024，11402052)；教育部新世紀優(yōu)秀人才支持計劃資助項目(NCET-11-0086)；江蘇省普通高校研究生科研創(chuàng)新計劃資助項目(CXZZ13_0084)；國家留學(xué)基金資助項目(201406090150)；中央高?；究蒲袠I(yè)務(wù)費專項資金資助項目