張建華，楚武利，楊曉彤，張晶輝

(1.西安航空學(xué)院飛行器學(xué)院，陜西西安710129；2.西北工業(yè)大學(xué)動(dòng)力與能源學(xué)院，陜西西安710072；3.先進(jìn)航空發(fā)動(dòng)機(jī)協(xié)同創(chuàng)新中心，北京100191；4.中國(guó)航發(fā)西安航空發(fā)動(dòng)機(jī)有限公司，陜西西安710021)

0 引 言

船用離心風(fēng)機(jī)一般在高溫、高濕且較為密閉的環(huán)境下運(yùn)行。根據(jù)GB/T 11865《船用離心通風(fēng)機(jī)》[1]的國(guó)家標(biāo)準(zhǔn)，船用離心風(fēng)機(jī)不僅要滿(mǎn)足普通風(fēng)機(jī)的安裝和運(yùn)行要求，還要滿(mǎn)足低振動(dòng)烈度、噪聲輻射運(yùn)行要求。當(dāng)風(fēng)機(jī)用于艙室內(nèi)通風(fēng)換氣時(shí)，風(fēng)機(jī)的出口連接長(zhǎng)管道，此時(shí)風(fēng)機(jī)的氣動(dòng)噪聲不能直接傳遞到外界，流體激勵(lì)風(fēng)機(jī)機(jī)殼引發(fā)的殼體表面振動(dòng)，振動(dòng)激發(fā)產(chǎn)生的結(jié)構(gòu)振動(dòng)噪聲則成為主要噪聲。事實(shí)上風(fēng)機(jī)的非定常流動(dòng)誘發(fā)的噪聲屬于流固干涉噪聲(流體激勵(lì)力對(duì)殼體振動(dòng)引發(fā)的在空氣中傳播的振動(dòng)噪聲)，葉輪和蝸殼是彈性體結(jié)構(gòu)，尤其在大型風(fēng)機(jī)中蝸殼的振動(dòng)不可忽略[2]。離心風(fēng)機(jī)的葉輪內(nèi)部流動(dòng)誘發(fā)殼體振動(dòng)的振動(dòng)噪聲鮮有學(xué)者研究，但這類(lèi)噪聲問(wèn)題在大型風(fēng)機(jī)系統(tǒng)中也顯得尤為突出，目前葉輪內(nèi)部流動(dòng)誘發(fā)蝸殼殼體振動(dòng)噪聲的研究?jī)H限于計(jì)算方法。20世紀(jì)末，Koopmann等[2]首次提出了一種基于邊界元(Boundary Element Method, BEM)計(jì)算和實(shí)驗(yàn)測(cè)量相結(jié)合的預(yù)測(cè)方法。此方法將分離氣動(dòng)噪聲，單獨(dú)計(jì)算大型離心風(fēng)機(jī)葉輪內(nèi)部的非定常流動(dòng)誘發(fā)的蝸殼殼體振動(dòng)噪聲，振動(dòng)噪聲計(jì)算所需的壓力脈動(dòng)通過(guò)實(shí)驗(yàn)獲得；隨后，Hwang 等[3]采用和文獻(xiàn)[2]相同的方法成功預(yù)測(cè)了冰箱用壓縮機(jī)的振動(dòng)聲輻射；21世紀(jì)初，蔡建成等[4]、Lu等[5]采用數(shù)值計(jì)算和實(shí)驗(yàn)相結(jié)合的方法，預(yù)測(cè)了工業(yè)用T9-19 No.4A離心風(fēng)機(jī)蝸殼的振動(dòng)聲輻射。振動(dòng)噪聲研究的目的是獲得振動(dòng)噪聲產(chǎn)生機(jī)理，提出有效的降噪方法。目前，結(jié)構(gòu)減振方法非常引人注目，它實(shí)際上是通過(guò)修改受控對(duì)象的動(dòng)力學(xué)特性參數(shù)使振動(dòng)滿(mǎn)足特定要求，不需要添加任何子系統(tǒng)的控制方法。而結(jié)構(gòu)減振主要以結(jié)構(gòu)優(yōu)化為主，離心風(fēng)機(jī)的蝸殼屬于薄壁結(jié)構(gòu)，而薄壁結(jié)構(gòu)的振動(dòng)輻射聲功率是關(guān)于結(jié)構(gòu)振動(dòng)速度的二次型函數(shù)[6]，通過(guò)優(yōu)化設(shè)計(jì)使得結(jié)構(gòu)振動(dòng)速度降低，則必能使其輻射的聲功率在一定范圍內(nèi)降低。周正等[7]、Lu等[5]通過(guò)優(yōu)化降低結(jié)構(gòu)表面振動(dòng)速度方法獲得了降噪效果。然而，上述學(xué)者在預(yù)測(cè)殼體振動(dòng)聲輻射的計(jì)算中，將計(jì)算分為兩個(gè)過(guò)程：(1) 以結(jié)構(gòu)有限元法獲得結(jié)構(gòu)的振動(dòng)速度；(2) 以聲學(xué)邊界元法獲得結(jié)構(gòu)的振動(dòng)聲輻射，計(jì)算過(guò)程較為繁瑣，容易出錯(cuò)。因此，有必要采用振動(dòng)、噪聲同步計(jì)算的方法，減小計(jì)算過(guò)程中可能存在的誤差。本文采用數(shù)值計(jì)算方法，給出一種基于振動(dòng)噪聲同步計(jì)算的振動(dòng)噪聲一體化優(yōu)化方法，在保持殼體質(zhì)量不變的前提下，有效地降低殼體振動(dòng)聲輻射。

1 數(shù)值計(jì)算模型

1.1 離心風(fēng)機(jī)模型參數(shù)

以某船用離心風(fēng)機(jī)為研究對(duì)象，該風(fēng)機(jī)依據(jù)使用環(huán)境不同，分為兩種安裝狀態(tài)：進(jìn)出口開(kāi)口和進(jìn)出口連接有封閉長(zhǎng)管道，本文所研究的風(fēng)機(jī)處在第二種安裝環(huán)境下，該風(fēng)機(jī)具體參數(shù)見(jiàn)表1。

1.2 非定常流場(chǎng)計(jì)算

1.2.1 流場(chǎng)網(wǎng)格

葉輪出口不穩(wěn)定流動(dòng)對(duì)殼體壁面的周期性沖擊是殼體振動(dòng)的主要激勵(lì)源。為了獲得較為精確的振動(dòng)激勵(lì)源，風(fēng)機(jī)內(nèi)部的流動(dòng)結(jié)構(gòu)全部采用高質(zhì)量的六面體結(jié)構(gòu)化網(wǎng)格處理。由于風(fēng)機(jī)計(jì)算模型較復(fù)雜，將風(fēng)機(jī)整個(gè)計(jì)算域分成四個(gè)主要部分： 過(guò)渡段和進(jìn)口管道計(jì)算域(包含深入葉輪內(nèi)的軸向進(jìn)氣空腔區(qū)域)，翼型葉輪網(wǎng)格計(jì)算域、帶有前后空腔的蝸殼整體網(wǎng)格計(jì)算域，節(jié)流閥和出口延伸段計(jì)算域。風(fēng)機(jī)計(jì)算域的網(wǎng)格數(shù)如表2所示。

表1 風(fēng)機(jī)參數(shù)Table 1 Parameters of centrifugal fan

表2 風(fēng)機(jī)計(jì)算域的網(wǎng)格數(shù)Table 2 Grid numbers in the calculational domain of fan system

各計(jì)算域之間通過(guò)交界面(interface)相互連接，在轉(zhuǎn)動(dòng)域和靜止域之間設(shè)定兩個(gè)重要的轉(zhuǎn)/靜交界面：其一，與葉輪相連接的進(jìn)口空腔和葉輪進(jìn)口(進(jìn)口-葉輪)交界面；其二，葉輪出口和蝸殼進(jìn)口(葉輪-蝸殼)交界面。為了減小計(jì)算模型和樣機(jī)的誤差，網(wǎng)格劃分過(guò)程中考慮了葉輪、過(guò)渡段厚度(2 mm)以及葉輪蓋板厚度(輪盤(pán)厚度為3 mm，輪蓋厚度為6 mm)、過(guò)渡段與葉輪進(jìn)口徑向間隙1mm，詳細(xì)參數(shù)如圖1所示。圖2顯示了風(fēng)機(jī)的總壓升系數(shù)和氣動(dòng)效率隨流量的變化曲線(xiàn)，網(wǎng)格數(shù)增加一倍后，總壓升基本保持不變，此曲線(xiàn)進(jìn)一步證明了網(wǎng)格數(shù)超過(guò)2.8×106后，網(wǎng)格達(dá)到無(wú)關(guān)解(網(wǎng)格數(shù)量對(duì)計(jì)算結(jié)果沒(méi)有影響)。

1.2.2 流場(chǎng)計(jì)算模型和邊界條件

圖1 風(fēng)機(jī)網(wǎng)格截面圖Fig.1 Surface mesh of blade and hub

圖2網(wǎng)格無(wú)關(guān)性驗(yàn)證Fig.2 Grid independence validation

基于商業(yè)軟件ANSYS CFX求解連續(xù)方程、動(dòng)量方程和標(biāo)準(zhǔn)k-ε湍流控制方程。各控制方程均采用有限元體積法離散。控制方程的空間離散采用高精度的高階格式差分，非定常計(jì)算的時(shí)間推進(jìn)格式為雙時(shí)間步長(zhǎng)全隱式格式，時(shí)間項(xiàng)使用二階后向歐拉差分離散，一個(gè)耦合求解器[8]依照壓力的耦合求解算法耦合求解連續(xù)方程和動(dòng)量方程。定常流場(chǎng)以穩(wěn)態(tài)解為初場(chǎng)，模擬設(shè)計(jì)點(diǎn)下最佳效率點(diǎn)(Best Efficiency Point, BEP)風(fēng)機(jī)內(nèi)部的非定常流動(dòng)狀況(本文的振動(dòng)噪聲的計(jì)算流量點(diǎn))。非定常時(shí)間步設(shè)定：葉輪有12個(gè)葉片通道，每個(gè)通道給定30個(gè)時(shí)間步，葉輪轉(zhuǎn)動(dòng)一周分成360個(gè)時(shí)間步，即葉輪每轉(zhuǎn)動(dòng)1°消耗一個(gè)時(shí)間步長(zhǎng)，每個(gè)時(shí)間步長(zhǎng)為5.7089×10-5s，這個(gè)時(shí)間步長(zhǎng)足夠進(jìn)行動(dòng)態(tài)壓力信號(hào)采集。計(jì)算結(jié)果顯示，當(dāng)葉輪轉(zhuǎn)動(dòng)5 400個(gè)時(shí)間步長(zhǎng)，即轉(zhuǎn)動(dòng)15周后，設(shè)定的監(jiān)測(cè)點(diǎn)達(dá)到穩(wěn)定的周期性波動(dòng)，此時(shí)判定計(jì)算收斂。

圖3 風(fēng)機(jī)氣動(dòng)性能曲線(xiàn)對(duì)比Fig.3 Comparison between numerical and experimental curves of fan aerodynamic performance

圖4 監(jiān)測(cè)截面1上三個(gè)軸向測(cè)量點(diǎn)的蝸殼壁P1截面上壓力脈動(dòng)頻譜Fig.4 Power spectra of the pressure fluctuation on P1 section of volute wall measured at three axial measurement points

圖3給出了風(fēng)機(jī)氣動(dòng)性能數(shù)值和實(shí)驗(yàn)的對(duì)比曲線(xiàn)。在設(shè)計(jì)轉(zhuǎn)速下(2 920 r·min-1)的最佳效率點(diǎn)(BEP)對(duì)應(yīng)的流量為Q=3.361 kg·s-1(φ=0.166)和總壓升PT=3182 Pa(ψ=0.420)。圖4給出了蝸殼壁面P1截面葉輪出口點(diǎn)位壓力脈動(dòng)頻譜圖。測(cè)點(diǎn)分布和描述詳見(jiàn)作者已發(fā)表文獻(xiàn)[8]。在P01點(diǎn)位的所有的數(shù)值和實(shí)驗(yàn)測(cè)量頻譜中，葉輪出口范圍內(nèi)(Z/B=0.07～0.36，Z為蝸殼到后板的距離，B為蝸殼軸向?qū)挾?，葉片通過(guò)頻率(Blade Passing Frequency, BPF)(對(duì)應(yīng)基頻584 Hz)處觀察到明顯的峰值。且在基頻上，計(jì)算和實(shí)驗(yàn)吻合較好，這表明對(duì)于文中最為關(guān)注的基頻噪聲，噪聲源計(jì)算模型是可靠的，詳細(xì)實(shí)驗(yàn)和計(jì)算分析過(guò)程參考文獻(xiàn)[8]。

1.3 振動(dòng)噪聲計(jì)算模型

1.3.1 振動(dòng)噪聲數(shù)學(xué)模型

定義振動(dòng)表面聲輻射的基爾霍夫(Kirchhoff)公式為[9]

式中：p表示聲壓，nv表示振動(dòng)的法向速度，H(f)為赫維賽德(Heaviside)函數(shù)，而δ(f)為狄拉克(Dirac)函數(shù)，方程右端第一項(xiàng)是由于振動(dòng)表面排開(kāi)流體體積產(chǎn)生的聲音，是單極子源；方程右端第二項(xiàng)是結(jié)構(gòu)處于其振動(dòng)產(chǎn)生的聲場(chǎng)中對(duì)聲場(chǎng)的散射影響，為偶極子源。式(1)經(jīng)過(guò)快速傅里葉變換后得到頻域結(jié)構(gòu)表面聲輻射控制方程即基爾霍夫-亥姆霍茲(Kirchhoff- Helmholtz)方程為

式中：W表示計(jì)權(quán)函數(shù)，q為聲源項(xiàng)(式(1)中右端的所有項(xiàng))。根據(jù)有限元理論獲得結(jié)構(gòu)振動(dòng)表面離散點(diǎn)(結(jié)構(gòu)網(wǎng)格節(jié)點(diǎn))的聲壓，通過(guò)插值函數(shù)逼近得到其他結(jié)構(gòu)節(jié)點(diǎn)聲壓，隨后求得聲功率值。

1.3.2 振動(dòng)噪聲數(shù)值計(jì)算過(guò)程

風(fēng)機(jī)蝸殼殼體有限元模型采用Shell63單元構(gòu)建。按照風(fēng)機(jī)殼體厚度不同分為三個(gè)板塊有：前板厚度TF、后板的厚度TB均為6 mm，殼體側(cè)板的厚度TS為5 mm，如圖5所示。基于表面四邊形網(wǎng)格劃分46 182個(gè)shell63單元。模型材料為鋼，密度為ρ=7 800 kg·m-3，彈性模量E=2.06×1011 Pa，泊松比γ=0.3。振動(dòng)邊界約束條件為：蝸殼前板設(shè)定10個(gè)均勻分布的固定螺栓，蝸殼后板設(shè)定4個(gè)均勻分布的固定螺栓，所有螺栓的三個(gè)方向平動(dòng)自由度為0，有限元網(wǎng)格詳細(xì)描述參考文獻(xiàn)[10]。針對(duì)殼體振動(dòng)噪聲，文中給出了具體的計(jì)算方法，計(jì)算流程如圖6所示。根據(jù)振動(dòng)激勵(lì)源計(jì)算、振動(dòng)計(jì)算、振動(dòng)噪聲數(shù)值計(jì)算將整個(gè)計(jì)算過(guò)程分為三個(gè)主要步驟：(1) 基于非定常流場(chǎng)計(jì)算獲得振動(dòng)激勵(lì)源；(2) 基于構(gòu)建的殼體結(jié)構(gòu)有限元模型計(jì)算模態(tài)參與因子；(3)加載模態(tài)參與因子和振動(dòng)激勵(lì)源，基于振動(dòng)噪聲單向計(jì)算方法計(jì)算殼體振動(dòng)聲輻射，詳細(xì)分析過(guò)程見(jiàn)參考文獻(xiàn)[11]。

圖5 蝸殼有限元模型Fig.5 FEM model of volute casing

圖6 振動(dòng)噪聲計(jì)算流程示意圖Fig.6 The flow chart of numerical evaluation method for casing vibro-acoustic calculations

1.3.3 振動(dòng)噪聲計(jì)算數(shù)值驗(yàn)證

風(fēng)機(jī)運(yùn)行時(shí)，內(nèi)部非定常流動(dòng)激勵(lì)機(jī)殼振動(dòng)，機(jī)殼產(chǎn)生微小變形，然而其表面變形的位移量級(jí)遠(yuǎn)小于流體壁面邊界層厚度，因此對(duì)于此類(lèi)設(shè)備結(jié)構(gòu)的振動(dòng)噪聲問(wèn)題一般采用考慮流體脈動(dòng)壓力單向作用于薄殼體壁面，計(jì)算殼體振動(dòng)在空氣中輻射噪聲時(shí)，忽略空氣的反作用。Jiang等[12]對(duì)離心泵內(nèi)部的非定常流動(dòng)誘發(fā)振動(dòng)問(wèn)題的計(jì)算和實(shí)驗(yàn)結(jié)果驗(yàn)證了此方法的合理性。

定義總振級(jí)的公式為

式中：aref=1×10-6m·s-2表示參考加速度值，afei表示頻譜內(nèi)單個(gè)頻率的振動(dòng)加速度有效值。

圖7給出了各振動(dòng)測(cè)點(diǎn)在20～3 000 Hz頻率范圍內(nèi)總振級(jí)的計(jì)算和實(shí)驗(yàn)對(duì)比。從圖7中可以看出絕大部分振動(dòng)測(cè)點(diǎn)的計(jì)算和實(shí)驗(yàn)吻合良好，雖然少數(shù)測(cè)點(diǎn)(測(cè)點(diǎn)1、4、13)有限元計(jì)算結(jié)果和實(shí)驗(yàn)值誤差較大，最大誤差為15 dB，但是總體來(lái)說(shuō)計(jì)算和實(shí)驗(yàn)誤差控制在一個(gè)量級(jí)范圍(20 dB)內(nèi)，這與文獻(xiàn)[5]中對(duì)T9-19 No.4A離心風(fēng)機(jī)蝸殼振動(dòng)的計(jì)算和實(shí)驗(yàn)所得的誤差分析一致。這表明，本文所采用的振動(dòng)噪聲計(jì)算方法是合理且有效的。

圖7 振動(dòng)加速度的計(jì)算值和實(shí)驗(yàn)值對(duì)比Fig.7 The comparison between calculated and tested amplitudes of vibration acceleration

2 振動(dòng)噪聲優(yōu)化過(guò)程

2.1 振動(dòng)噪聲優(yōu)化數(shù)學(xué)模型

對(duì)于開(kāi)口域的噪聲控制問(wèn)題，目標(biāo)函數(shù)一般選取外部輻射的聲功率，此種方法已經(jīng)被諸多學(xué)者所驗(yàn)證[13]。結(jié)構(gòu)表面輻射聲功率和外部激振力Fi、模態(tài)振型φi、頻率放大因子βi之間的關(guān)系定義為[14]

式中：Wo,active表示結(jié)構(gòu)表面有功輻射聲功率，Re(In)表示聲強(qiáng)實(shí)部。當(dāng)外加激振力Fi確定的情況下(流動(dòng)激勵(lì)力不變)，減小結(jié)構(gòu)振型φi可以達(dá)到降噪的效果。結(jié)構(gòu)的振型和固有頻率取決于結(jié)構(gòu)的幾何形狀、厚度分布、結(jié)構(gòu)剛度以約束位置，當(dāng)結(jié)構(gòu)幾何形狀、結(jié)構(gòu)剛度和約束位置一定時(shí)，可以通過(guò)控制結(jié)構(gòu)厚度分布來(lái)改變結(jié)構(gòu)的振型。由于風(fēng)機(jī)蝸殼的厚度一般小于10 mm，我們將蝸殼板厚限制為4～10 mm。本部分優(yōu)化的數(shù)學(xué)模型如下：

定義蝸殼表面輻射聲功率(Kirchhoff Sound Power of Wall, Kirchhoff SPW)為目標(biāo)函數(shù)Ws：

設(shè)計(jì)變量TF、TS、TB，約束條件為

聲功率級(jí)定義為

式中：Wref表示參考聲功率值，Wref=1×10-12W。

2.2 優(yōu)化過(guò)程及分析

為了提高優(yōu)化效率，文中采用了試驗(yàn)設(shè)計(jì)方法(Design of Experiment, DOE)和近似模型(Approximation Models, AM)相結(jié)合的方法來(lái)搜尋結(jié)構(gòu)設(shè)計(jì)變量和目標(biāo)函數(shù)之間的確定關(guān)系，為后面的優(yōu)化提供基礎(chǔ)數(shù)據(jù)和模型。蝸殼結(jié)構(gòu)的振動(dòng)噪聲優(yōu)化涉及多個(gè)程序的順序啟動(dòng)，因此本文將多個(gè)不同的計(jì)算軟件(UG、ANSA、Nastran、LMS)整合到多學(xué)科仿真優(yōu)化平臺(tái)Isight中，基于多學(xué)科多目標(biāo)優(yōu)化平臺(tái)Isight完成風(fēng)機(jī)蝸殼殼體振動(dòng)噪聲的優(yōu)化過(guò)程。圖8給出了Isight的集成示意圖，從圖中可以看出，優(yōu)化過(guò)程主要分為3個(gè)步驟：(1) 基于DOE創(chuàng)建樣本點(diǎn)設(shè)計(jì)空間(Design space)；此部分耗時(shí)最長(zhǎng)，也是計(jì)算的中心環(huán)節(jié)。針對(duì)本文所研究的船用離心風(fēng)機(jī)蝸殼，以各蝸殼板厚(前板TF、后板TB、側(cè)板TS)為設(shè)計(jì)變量，每個(gè)設(shè)計(jì)變量給定5個(gè)水平，表3給出了各設(shè)計(jì)變量的水平分布。由于設(shè)計(jì)變量較少，在三個(gè)變量所構(gòu)建的設(shè)計(jì)空間中采用全因子法進(jìn)行樣本點(diǎn)的采集，共采集125個(gè)樣本點(diǎn)。對(duì)于每個(gè)樣本點(diǎn)采用前文的振動(dòng)噪聲單向數(shù)值計(jì)算方法進(jìn)行數(shù)值計(jì)算，獲得各變量組合蝸殼結(jié)構(gòu)表面的輻射聲功率。(2) 基于步驟(1)中完成的樣本點(diǎn)通過(guò)插值方法(RBF 近似模型，對(duì)于非線(xiàn)性問(wèn)題最為有效)獲取優(yōu)化的近似數(shù)學(xué)函數(shù)。圖9給出了近似模型的擬合精度，從圖中可以看出近似擬合計(jì)算值和數(shù)值計(jì)算值基本重合，蝸殼表面輻射聲功率(Kirchhoff SPW)的R2無(wú)限接近1，因而所構(gòu)建的近似模型完全能替代實(shí)際的蝸殼表面輻射聲功率、模型，用于后面的優(yōu)化設(shè)計(jì)工作。(3) 基于步驟(2)得到的擬合函數(shù)，采用自適應(yīng)模擬退火算法全局優(yōu)化(Adaptive Simulated Annealing, ASA)獲得優(yōu)化解，以全局優(yōu)化解為初值，通過(guò)混合整型序列二次規(guī)劃局部尋優(yōu)獲得最優(yōu)解。

圖8 Isight優(yōu)化整合示意圖Fig.8 The schematic diagram of Isight optimizing integration

表3 蝸殼各板厚設(shè)計(jì)參數(shù)的水平分布Table 3 The horizontal distributions of the design parameters of volute wall thickness (TF, TS and TB)

圖9 蝸殼表面聲輻射值擬合誤差分析Fig.9 Analysis of fitting error

圖10給出了設(shè)計(jì)變量(TF、TS、TB)和響應(yīng)目標(biāo)(蝸殼結(jié)構(gòu)表面輻射聲功率Kirchhoff SPW)之間的相關(guān)性系數(shù)分布。圖10中的正值表示目標(biāo)響應(yīng)和設(shè)計(jì)變量之間的正比例關(guān)系，反之，為反比例關(guān)系；系數(shù)絕對(duì)值越接近1，相關(guān)性程度越高。從圖10中可以看出，TS對(duì)響應(yīng)相關(guān)性程度最高，其次是TF，再次是TB，TS、TB與輻射聲功率成反比例關(guān)系，這表明TS、TB越大，輻射聲功率越小。然而，圖11蝸殼各板塊的主效應(yīng)圖則顯示，隨著TF厚度增大蝸殼表面輻射聲功率呈現(xiàn)先減小后增大趨勢(shì)，當(dāng)TF值為6.3 mm時(shí)，蝸殼輻射聲功率最??；TB相對(duì)于TF則呈現(xiàn)相反的變化趨勢(shì)，TB隨著厚度的增加，蝸殼表面輻射聲功率表現(xiàn)出先增大后減小的趨勢(shì)；而TS隨著厚度的增加，蝸殼表面輻射聲功率值呈現(xiàn)直線(xiàn)下降趨勢(shì)，但是當(dāng)TS大于8.8 mm之后，蝸殼表面輻射聲功率幅值基本上保持不變。以上分析表明，TF、TS、TB均存在最佳合理范圍，這需要優(yōu)化來(lái)獲得最合理的TB，TS、TB值滿(mǎn)足蝸殼表面輻射聲功率值最小。圖12給出了此類(lèi)單目標(biāo)優(yōu)化的流程圖。單目標(biāo)優(yōu)化采用模擬退火算法在樣本空間內(nèi)進(jìn)行全局搜索，以前文的近似模型最優(yōu)結(jié)果為初值，全局優(yōu)化迭代了10 000步，耗時(shí)12 min；隨后以ASA最優(yōu)結(jié)果為初值采用混合整型序列二次規(guī)劃(Mixed-Integer Squential Quadratic Programming,MISQP)進(jìn)行局部尋優(yōu)，此過(guò)程迭代了12步，耗時(shí)僅為數(shù)秒。

圖10 壁厚設(shè)計(jì)參數(shù)和蝸殼結(jié)構(gòu)表面輻射聲功率的相關(guān)性分析Fig.10 The correlation analysis between the design parameters of volute wall and the acoustic power radiated from the surface of volute structure (Kirchhoff SPW)

圖11 蝸殼壁厚和輻射聲功率關(guān)系Fig.11 The relationships between the volute wall thickness and Kirchhoff SPW

圖12 單目標(biāo)優(yōu)化流程圖Fig.12 The flow chart of the single-objective optimization(Kirchhoff SPW)

表4給出了保持質(zhì)量不變的優(yōu)化結(jié)果，若保持質(zhì)量不變，蝸殼結(jié)構(gòu)表面的輻射聲功率也有較大程度的減弱，平均減弱了6.3 dB。由于工程上常用尺寸圓整或是近似厚度的鋼板，因此設(shè)定蝸殼3個(gè)板厚度分別為T(mén)F=4.5 mm，TS=7.5 mm，TB=4 mm，經(jīng)過(guò)數(shù)值計(jì)算得到蝸殼表面的輻射聲功率為59.42 dB，厚度誤差為1.26%，在許可范圍內(nèi)。圖13的蝸殼結(jié)構(gòu)表面的輻射聲功率頻譜也顯示出，基頻聲功率幅值最為突出，優(yōu)化明顯地改善了蝸殼表面的基頻輻射聲功率，降低了6.23 dB。為了驗(yàn)證頻率步長(zhǎng)對(duì)聲功率頻譜尖峰處的數(shù)值影響，圖14給出了以1 Hz為頻率步長(zhǎng)，頻率波段在573～593 Hz之間的聲功率頻譜曲線(xiàn)。從圖14中可以看出，基頻處的聲功率依然最為突出，優(yōu)化后基頻處的聲功率值降低了6.03 dB。從圖15和圖16中基頻下蝸殼表面以及不同截面的振動(dòng)聲輻射分布對(duì)比中可以分析出：優(yōu)化使得蝸殼側(cè)板靠近蝸舌附近的振動(dòng)輻射聲壓大幅度減小，蝸殼其他區(qū)域的振動(dòng)輻射聲壓也有不同程度的減弱，優(yōu)化也改變了聲輻射的指向性分布，在蝸殼后板一側(cè)產(chǎn)生了非常明顯的指向性。

表4 蝸殼總質(zhì)量不變的單目標(biāo)優(yōu)化結(jié)果Table 4 The single-objective optimization results while the volute total mass keeping constant

圖13 蝸殼結(jié)構(gòu)表面輻射聲功率譜Fig.13 The radiated sound power spectrum onthe volute structural surface

圖14 蝸殼結(jié)構(gòu)表面的輻射聲功率譜Fig.14 The radiated sound power spectrum on the surface of volute

圖15 基頻下原型風(fēng)機(jī)蝸殼振動(dòng)聲輻射分布Fig.15 The sound radiation distribution of the original volute vibration at the fundamental frequency

圖16 基頻下優(yōu)化的風(fēng)機(jī)蝸殼振動(dòng)聲輻射分布Fig.16 The sound radiation distribution of the optimized volute vibration at the fundamental frequency

3 結(jié) 論

(1) 對(duì)于大型帶有長(zhǎng)管道的船用風(fēng)機(jī)系統(tǒng)，研究表明，風(fēng)機(jī)的振動(dòng)噪聲最為突出，為了降低風(fēng)機(jī)振動(dòng)噪聲，本文基于DOE提出了一種考慮振動(dòng)噪聲干涉作用(流體脈動(dòng)壓力單向作用于蝸殼壁面，在計(jì)算蝸殼振動(dòng)在空氣中輻射噪聲時(shí)忽略空氣的反作用)的優(yōu)化設(shè)計(jì)方法，優(yōu)化后風(fēng)機(jī)殼體結(jié)構(gòu)表面輻射聲功率均有不同程度降低，基頻處的噪聲降幅最為顯著，最大降低了6.23 dB。

(2) 優(yōu)化后，蝸殼表面的基頻(584 Hz)輻射聲功率均有較大程度降低，尤其是風(fēng)機(jī)蝸舌附近聲輻射功率降幅最為顯著，且優(yōu)化也改善了振動(dòng)噪聲的指向性分布。

本文的優(yōu)化方法可供此類(lèi)風(fēng)機(jī)振動(dòng)噪聲控制參考和借鑒應(yīng)用，獲得預(yù)期的社會(huì)經(jīng)濟(jì)價(jià)值。