趙 華，樂 晨，谷小軍，朱繼宏，張衛(wèi)紅

(1. 西北工業(yè)大學機電學院，西安 710072；2. 北京宇航系統(tǒng)工程研究所，北京 100076；3. 西北工業(yè)大學無人系統(tǒng)研究院，西安 10072)

0 引言

在火箭發(fā)射過程中，往往會伴隨著嚴重的振動和噪聲，這不僅會引起火箭整流罩內(nèi)儀器設(shè)備出現(xiàn)故障，甚至可能導致發(fā)射任務的失敗。因此，研究火箭發(fā)射過程中整流罩內(nèi)的噪聲響應和降噪方案將變得非常重要。

運載火箭發(fā)射階段整流罩內(nèi)的噪聲來源于很多方面，但主要包含兩部分，一是火箭發(fā)動機的排氣噪聲以及整流罩與空氣摩擦產(chǎn)生的氣動噪聲,二是火箭結(jié)構(gòu)振動傳遞到整流罩的振動噪聲。通過將兩者疊加，可以得到一個總的整流罩外部噪聲源。

由于整流罩噪聲分析的頻譜很寬，一般被劃分為低頻和中高頻單獨分析。對于低頻，系統(tǒng)的各動力學參數(shù)比較確定，結(jié)構(gòu)和聲腔的模態(tài)密度也較小，因此一般采用有限元方法進行分析；但是對于中高頻，結(jié)構(gòu)系統(tǒng)較為復雜，系統(tǒng)的各動力學參數(shù)具有隨機性，并且結(jié)構(gòu)和聲腔的模態(tài)密度比較大，因此常采用統(tǒng)計能量分析(SEA)方法進行分析。如Hwang等采用SEA方法分析了在發(fā)射階段聲激勵下火星航天器的聲振響應，在高頻段其仿真與實驗測試結(jié)果的聲壓級基本吻合。Larko等運用SEA方法仿真分析了航天器上蜂窩夾層板的高頻聲振響應，仿真與實驗測試結(jié)果誤差較小。姚德源等通過SEA方法預測了飛行器圓柱殼儀器艙段聲振環(huán)境。孫目等運用Auto SEA2軟件預測了某儀器艙在高頻段的聲振響應以及衛(wèi)星整流罩內(nèi)聲場環(huán)境，仿真與試驗數(shù)據(jù)基本吻合。另外，王懷志等采用能量有限元分析(EFEA)方法預示了雙星整流罩中頻聲振環(huán)境，并與SEA方法和實驗測試結(jié)果進行了對比。曾耀祥等基于SEA方法分析了某火箭儀器艙的聲振環(huán)境，并與實驗數(shù)據(jù)進行了對比。在本文中由于結(jié)構(gòu)系統(tǒng)較為復雜且主要關(guān)注整流罩中高頻段內(nèi)的噪聲響應，因此通過SEA方法對其進行分析。

在以往整流罩的噪聲響應分析和降噪方案研究中，沒有考慮整流罩內(nèi)有效負載等結(jié)構(gòu)占用聲腔的影響，這勢必會影響整流罩內(nèi)噪聲響應分析的準確性。本文通過將衛(wèi)星負載等簡化為圓柱和立方體兩種包絡結(jié)構(gòu)分別在直徑和長度變化下對整流罩內(nèi)的噪聲響應進行分析。并研究不同降噪材料在不同參數(shù)下的降噪效果，給以后的整流罩降噪設(shè)計提供一些參考。

1 基于SEA方法的整流罩分析模型

SEA方法是利用統(tǒng)計學的概念研究各子系統(tǒng)之間的能量傳遞和損耗，通過系統(tǒng)振動能量與模態(tài)之間的內(nèi)在聯(lián)系建立聲學和結(jié)構(gòu)振動的動力學方程。SEA方法模型的參數(shù)主要有：模態(tài)密度、輸入功率、內(nèi)損耗因子和耦合損耗因子。對于有個子系統(tǒng)的SEA模型，當整個系統(tǒng)處于穩(wěn)態(tài)時，可以建立如下能量平衡方程

(1)

式中

[]=

(2)

表示分析帶寬內(nèi)的中心頻率,表示子系統(tǒng)的內(nèi)損耗因子，表示子系統(tǒng)和間的耦合損耗因子，表示子系統(tǒng)的模態(tài)密度，表示子系統(tǒng)的輸入能量，,表示平衡狀態(tài)下子系統(tǒng)存儲的能量。

本文采用VAOne軟件進行火箭整流罩噪聲響應的建模和分析，整流罩結(jié)構(gòu)簡化為板殼結(jié)構(gòu)，劃分各子系統(tǒng)時不僅要考慮到實際各子系統(tǒng)的連接，同時為了保證能夠采用SEA方法進行求解，各子系統(tǒng)在分析帶寬內(nèi)的振型數(shù)要大于1。整流罩結(jié)構(gòu)主要為蜂窩夾層板，由于無法在VAOne軟件中直接定義，因此采用三明治夾芯板理論將蜂窩夾層板的夾心結(jié)構(gòu)等效為高度不變的正交各向異性材料板，上下面板保持原有屬性。圖1是整流罩結(jié)構(gòu)在考慮衛(wèi)星負載下的SEA方法模型(以圓柱形衛(wèi)星包絡結(jié)構(gòu)為例)，為方便后面討論，對本文關(guān)注的圓柱聲腔和部分圓錐聲腔編號1～4。圖2是作用于整流罩外部的聲壓載荷，表1和表2分別是蜂窩結(jié)構(gòu)的尺寸參數(shù)和材料參數(shù)。

(a) 結(jié)構(gòu)子系統(tǒng)

(b) 聲腔子系統(tǒng)圖1 整流罩SEA模型Fig.1 Statistical energy analysis model for fairing

圖2 整流罩外部的聲壓載荷Fig.2 Sound pressure loads acting on the fairing

表1 蜂窩結(jié)構(gòu)的尺寸參數(shù)Tab.1 Dimensional parameters of honeycomb structures

表2 蜂窩結(jié)構(gòu)的材料參數(shù)Tab.2 Material parameters of honeycomb structures

2 火箭整流罩的噪聲響應分析

為使仿真結(jié)果更加準確，因此需要考慮整流罩內(nèi)衛(wèi)星等有效負載結(jié)構(gòu)對噪聲響應的影響。整流罩內(nèi)衛(wèi)星結(jié)構(gòu)形式非常復雜，給建模帶來不小困難，因此，在VAOne軟件建模中將衛(wèi)星結(jié)構(gòu)簡化為圓柱或立方體的包絡形式。對兩種簡化形式分別進行仿真分析，討論隨包絡結(jié)構(gòu)尺寸改變的聲壓級變化規(guī)律，以便為后續(xù)整流罩內(nèi)衛(wèi)星的排布提供一些參考。

2.1 圓柱形衛(wèi)星包絡結(jié)構(gòu)的噪聲響應分析

針對圓柱形衛(wèi)星包絡結(jié)構(gòu)，分析不同直徑和高度下的噪聲響應結(jié)果。其中，直徑分別為4,5,6 m，高度分別為5,10,13 m。

2.1.1 直徑變化對圓柱(錐)聲腔聲壓級的影響

在原始模型的基礎(chǔ)上僅改變圓柱形衛(wèi)星包絡結(jié)構(gòu)的直徑，分別在高度為5,10,13 m這3種情況下分析直徑的變化對圓柱(錐)聲腔聲壓級的影響。分析結(jié)果如圖3所示，在相同高度下，圓柱(錐)聲腔的聲壓級隨著直徑的增大而減小，表明直徑的增加有利于圓柱(錐)聲腔聲壓級的降低。

(a) 高度5 m

(b) 高度10 m

(c) 高度13 m圖3 相同高度下圓柱(錐)聲腔聲壓級隨直徑的變化規(guī)律Fig.3 Variation of sound pressure level with diameter for a cylindrical (conical) cavity at the same height

2.1.2 高度變化對圓柱(錐)聲腔聲壓級的影響

同理，在原始模型的基礎(chǔ)上僅改變圓柱形衛(wèi)星包絡結(jié)構(gòu)的高度，在直徑分別為5 m和6 m的情況下分析高度變化對圓柱(錐)聲腔聲壓級的影響。分析結(jié)果如圖4所示，在相同直徑下，各圓柱(錐)聲腔的總聲壓級隨著高度的增加大部分呈降低的趨勢，表明高度的增加有助于降低各聲腔的總聲壓級。

(a) 直徑5 m

(b) 直徑6 m圖4 相同直徑下圓柱(錐)聲腔聲壓級隨高度的變化規(guī)律Fig.4 Variation of sound pressure level with height for a cylindrical (conical) cavity at the same diameter

2.2 立方體形衛(wèi)星包絡結(jié)構(gòu)下的噪聲響應分析

針對立方體形衛(wèi)星包絡結(jié)構(gòu)，分析立方體不同底面對角線長度和高度下的噪聲響應結(jié)果。其中，底面對角線長度分別為4，5,6 m，高度分別為5,10,13 m。

2.2.1 底面對角線長度變化對圓柱(錐)聲腔聲壓級的影響

與圓柱形衛(wèi)星包絡結(jié)構(gòu)下的分析方式一樣，通過控制變量法，在高度分別為5,10，13 m這3種情況下分析底面對角線長度的變化對圓柱(錐)聲腔聲壓級的影響。分析結(jié)果如圖5所示，在相同高度下，圓柱(錐)聲腔的聲壓級隨著底面對角線長度變化沒有明顯改變。

(a) 高度5 m

(b) 高度10 m

(c) 高度13 m圖5 相同高度下圓柱(錐)聲腔聲壓級隨底面對角線長度的變化規(guī)律Fig.5 Variation of sound pressure level with the diagonal length of the bottom surface for a cylindrical (conical) cavity at the same height

2.2.2 高度變化對圓柱(錐)聲腔聲壓級的影響

針對立方體衛(wèi)星包絡結(jié)構(gòu)，在底面對角線長度分別為5 m和6 m的情況下分析高度變化對圓柱(錐)聲腔聲壓級的影響。分析結(jié)果如圖6所示，在相同底面對角線長度下，隨著高度的增加，各圓柱(錐)聲腔的總聲壓級呈降低的趨勢。

(a) 底面對角線長度5 m

(b) 底面對角線長度6 m圖6 相同底部對角線長度下圓柱(錐)聲腔聲壓級隨高度的變化規(guī)律Fig.6 Variation of sound pressure level with height for a cylindrical (conical) cavity at the same diagonal length of the bottom surface

3 降噪材料對整流罩噪聲響應的影響

對于整流罩的減振降噪，主動控制或者半主動控制對于低頻噪聲有比較好的降噪效果，但由于成本和結(jié)構(gòu)的復雜性，實際實施起來有些困難。因此，對于中高頻段噪聲可以采用較為簡便的被動吸聲的方法進行減振降噪，其中泡沫和纖維材料在實際應用當中的效果比較突出。

以直徑和高度均為5 m的圓柱形衛(wèi)星包絡結(jié)構(gòu)為分析模型，下面針對3種不同材料的降噪效果進行分析討論。圖7表示3種材料在不同鋪設(shè)厚度下的整流罩聲壓級響應值，并且與未鋪設(shè)降噪材料下的整流罩聲壓級響應值進行對比，分析3種不同降噪材料的降噪效果，圖8對3種材料不同鋪設(shè)厚度下的總聲壓級和質(zhì)量進行了對比。

從圖7和圖8結(jié)果中可以發(fā)現(xiàn)：1)在較低頻率時，降噪材料對整流罩內(nèi)聲壓級響應的影響十分微弱，僅在1 000 Hz以上時聲壓級響應隨厚度變化有較為明顯的改變，且隨著厚度的增加而愈發(fā)明顯。2)在相同厚度下，毛氈和玻璃纖維的質(zhì)量相同且相比于三聚氰胺泡沫更重，但同時發(fā)現(xiàn)毛氈和玻璃纖維的降噪效果也更為明顯。因此，若要從玻璃纖維、毛氈和三聚氰胺泡沫中選擇最優(yōu)的一種，則需要比較三者單位質(zhì)量下的降噪值(降噪值指鋪設(shè)降噪材料與未鋪設(shè)降噪材料下的總聲壓級的差值)。

(a) 玻璃纖維

(b) 三聚氰胺泡沫

(c) 毛氈圖7 降噪材料在不同鋪設(shè)厚度下整流罩的聲壓級響應隨頻率的變化規(guī)律Fig.7 Variation of sound pressure level response of fairing with frequency under different laying thickness of noise reduction material

(a) 不同厚度下的總聲壓級

(b) 不同厚度下的質(zhì)量圖8 降噪材料在不同鋪設(shè)厚度下總聲壓級響應和質(zhì)量對比Fig.8 Comparison of total sound pressure level response and quality of noise reduction materials at different laying thicknesses

表3～5分別是3種不同材料在不同鋪設(shè)厚度下的降噪值、質(zhì)量和單位質(zhì)量下的降噪值，圖9是3種材料在不同厚度下的降噪值及單位質(zhì)量降噪值的曲線圖。從中可以發(fā)現(xiàn)：1)在相同厚度下，玻璃纖維、毛氈、三聚氰胺泡沫的降噪值依次降低，并且玻璃纖維和毛氈降噪值隨厚度增加更為明顯，而三聚氰胺泡沫降噪值隨厚度的增加變化則較為平緩。2)三聚氰胺泡沫的單位質(zhì)量降噪能力明顯高于玻璃纖維和毛氈，但在相同厚度下的降噪能力則明顯弱于玻璃纖維和毛氈。因此，在有降噪值和鋪設(shè)降噪材料厚度要求時則需要綜合考慮降噪值、質(zhì)量和單位質(zhì)量降噪值，以獲得最優(yōu)的降噪材料鋪設(shè)方案。

表3 單位質(zhì)量玻璃纖維的降噪值Tab.3 Noise reduction value of glass fiber with per unit mass

表4 單位質(zhì)量三聚氰胺泡沫的降噪值Tab.4 Noise reduction value of melamine foam with per unit mass

表5 單位質(zhì)量毛氈的降噪值Tab.5 Noise reduction value of felt with per unit mass

(a) 3種材料不同厚度下的降噪值

(b) 3種材料單位質(zhì)量降噪值對比圖9 降噪材料在不同鋪設(shè)厚度下的降噪值及單位質(zhì)量降噪值Fig.9 Noise reduction values and noise reduction values with per unit mass for different laying thicknesses of noise reduction materials

4 結(jié)論

本文基于SEA方法建立了考慮衛(wèi)星結(jié)構(gòu)的重型火箭整流罩噪聲響應分析模型，并通過在整流罩內(nèi)壁鋪設(shè)降噪材料，分析各種降噪材料在不同參數(shù)下的降噪效果，具體工作和創(chuàng)新點如下：

1)相比于前人工作，本文提出了考慮衛(wèi)星結(jié)構(gòu)等負載對整流罩內(nèi)噪聲響應的影響?？紤]到衛(wèi)星結(jié)構(gòu)形式的復雜性，給建模帶來一定困難，因此將衛(wèi)星結(jié)構(gòu)簡化為圓柱或立方體的包絡形式，由此建立了重型火箭整流罩的SEA模型。

2)針對圓柱和立方體兩種不同的衛(wèi)星包絡形式，分別采用控制變量法分析了衛(wèi)星包絡結(jié)構(gòu)高度和直徑(底面對角線長度)的變化對圓柱(錐)聲腔聲壓級的影響。圓柱形衛(wèi)星包絡結(jié)構(gòu)直徑和高度的增加有助于降低各圓柱(錐)聲腔的聲壓級，另外，隨著立方體形衛(wèi)星包絡結(jié)構(gòu)高度的增加，各圓柱(錐)聲腔的聲壓級呈降低的趨勢，但底面對角線長度變化對圓柱(錐)聲腔的聲壓級沒有明顯影響，以上結(jié)論可為后續(xù)整流罩內(nèi)衛(wèi)星的排布形式提供一些參考。

3)通過仿真分析玻璃纖維、毛氈和三聚氰胺泡沫在不同鋪設(shè)厚度下的降噪值和單位質(zhì)量降噪值，可以發(fā)現(xiàn)在相同鋪設(shè)厚度下，三聚氰胺泡沫的降噪能力明顯弱于玻璃纖維和毛氈，但若考慮到降噪材料的質(zhì)量約束時，三聚氰胺泡沫的單位質(zhì)量降噪能力則明顯優(yōu)于另外兩者。因此，在實際應用中不僅要滿足降噪要求，同時也要考慮降噪材料鋪設(shè)質(zhì)量和厚度約束，以此通過鋪設(shè)單一降噪材料或多種降噪材料混合的方式獲得最優(yōu)的鋪設(shè)方案。