基于IWO算法的軸系-基座-殼體系統聲輻射優化研究

付衛華，蔡耀全，劉金林，賴國軍

(1. 廣州航海學院 船舶與海洋工程學院，廣東 廣州 510000；2. 中國人民解放軍92067部隊，廣東 湛江 524000；3. 海軍工程大學 動力工程學院，湖北 武漢 430033)

0 引 言

螺旋槳推進軸系引起的船體尾部振動和噪聲輻射，一直是人們重點關注和研究的問題。對于復雜連續結構的問題，目前研究的主要方法有有限元法、功率流有限元法和統計能量法等。在中低頻段，采用功率流有限元法是分析復雜結構振動能量密度和強度的有效方法，該方法采用有限元法進行數值計算，從能量的角度來分析問題。因此，它不僅對各節點的能量進行分析，還能夠彌補傳統有限元法只適用于低頻的缺陷。

關于功率流有限元的研究主要有：伍先俊等[1-2]從功率流理論著手，結合Ansys和Isight軟件，對組件功率流的計算方法和隔振系統進行研究，取得良好優化效果；賀云南等[3]將功率流有限元法和邊界元結合，開發了以功率流有限元法（ PFFEM）為基礎的聲輻射預測程序，對水下復雜結構—點力激勵下潛艇的振動和聲輻射特性進行預測，獲得了良好效果；趙群等[4]基于振動的功率流理論和概率攝動法，提出了頻域內振動傳遞路徑的路徑功率流傳遞度的新概念和方法，以此為基礎研究功率流傳遞度對工程中不確定因素的敏感程度，從而指導系統結構參數的優化設計；陳爐云等[5]將功率流有限元法和聲學邊界元方程結合，將法向功率流定義為目標函數進行結構噪聲問題研究，用遺傳算法對一加肋板結構進行優化計算分析，計算結果驗證了該方法的可行性和有效性。楊德慶等[6]給出基于有限元法的各種力學構件如隔振器、梁、板及其組合結構有限元功率流計算公式，研究了有限元功率流落差與振級落差這2種評價方法的優缺點。肖功煜等[7]基于功率流有限元法研究了板、殼結構的振動功率流特性，分析結構阻尼、阻尼器、加強筋等參數對結構功率流的影響。

從公開的文獻資料來看，以功率流有限元法為基礎對簡單的梁、板、殼等結構進行優化研究較多，對復雜系統的結構優化研究較少。因此本文結合功率流有限元法和聲學邊界元法，對軸系-基座-殼體這一復雜系統進行研究，給出相應的減振降噪措施，利用入侵性野草算法（IWO）[8]，以流經傳遞路徑總功率流為目標函數，對相關設計參數進行優化，最后通過場點聲壓來分析比較該優化方法對系統聲輻射優化的效果。

1 入侵性野草算法

野草算法（Invasive Weed Optimization，IWO）是近年來提出的一種簡單、有效的基于種群的優化算法，最早由Mehrabian和Lucus為解決數值優化問題而提出。該算法自提出以來，其較強的魯棒性、自適應性和隨機性使得其得到了廣泛關注，并應用于解決實際多參數連續函數的優化問題中[9-10]。該算法啟發于野草繁殖這一生物現象，其核心思想是每棵野草對環境具有不同的適應度，適應度高的野草會繁殖出更多的種子，從而具備該野草基因的種群存活率越高。因此該算法的目標找到更高適應度的野草。

該算法可分為4步：1）初始化；2）繁殖；3）空間分布；4）競爭性生存。設野草種群中野草的初始數量Pinit，最大種群規模為Pmax，最大迭代次數Imax、野草最大和最小可生成種子數Smax和Smin、非線性指數n，種子散布的初始步長σinit和最終值步長σfinal。野草種群中每根野草可根據自身適應度及種群中所有個體的最小適應度和最大適應度確定產生種子的數目。這樣適應度高的可以產生較多的種子，適應度低產生較少的種子，采用這種機制可保證適應度高的野草能夠最終生存下來。確定種子數量公式如下：

式中：Fcur，Fmax和Fmin分別為當前種群中當前野草適應度，最大適應度和最小適應度。

根據實際問題若每個野草個體變量維數D，某個父輩野草為并隨機產生Pinit個初始解，分配到D維空間中。新種子變量值由產生它的父輩加上某個數值Z。該值在D維空間中服從均值為0，標準方差為σcur的正態分布（即σcur值隨著迭代次數的增加而減小。假設則新種子的變量值為：

其中σcur計算公式如下：

式中：σini，σfin分別為標準方差初始值和最終設定值；Imax，I分別為最大迭代次數和當前迭代次數；n為非線性指數，可控制標準偏差衰減速度。根據式（3）可知，當迭代次數增加時σcur逐漸減小，這樣適應度高的種群逐漸聚集，算法由全局搜索向局部搜索轉移，從而優化野草種群的散布空間，最終目標值趨向收斂。

2 系統減振優化框架構建

2.1 系統模型建立

圖 1 系統尾部模型Fig. 1 The system stern model

以研究的推進軸系試驗平臺軸系為基礎，建立軸系-基座-殼體系統模型。圖1為系統尾部模型，殼體采用shell181單元，后尾軸承采用5個COMBI214單元支撐，前尾軸承采用4個COMBI214單元支撐，推力軸承支撐采用1個MTRI27單元（集成縱向推力塊油膜剛度）和3個COMBI214單元支撐，電機支撐采用2個COMBI214單元支撐在電機殼體上。對于該系統其振動噪聲主要來源激振力引起的軸系和尾部縱向耦合振動以及回旋耦合振動，因此本文提出的系統減振措施如下：在推力軸承和半聯軸器間設1個縱向減振器；推力軸承右端加5個縱向吸振器；電機外殼到殼體之間采用4個隔振器支撐，電機外殼橫向對稱布置6個吸振器；前尾軸承處內殼外設4個吸振器上下左右均布，內殼內采用4個吸振器均布。縱向減振器采用MATRI27單元模擬，便于模擬軸系轉動，所有吸振器和隔振器采用COMBI14單元，質量單元為MASS21。所優化的參數有縱向減振器質量X1、剛度X2和阻尼X3，推力軸承縱向吸振器質量X4、剛度X5和阻尼X6，電機外殼吸振器質量X7、剛度X8和阻尼X9，電機隔振器質量X10、剛度X11和阻尼X12，前尾軸承處內殼外吸振器質量X13、剛度X14和阻尼X15以及內殼內吸振器質量X16、剛度X17和阻尼X18。對于上述參數均給出上下限。因此，在采用IWO算法時，每個野草個體變量維數為D=18，則

2.2 適應度函數建立

本文研究的系統結構-聲輻射優化往往是以尾部聲場點的聲壓為目標函數，而場點聲壓一般是通過聲邊界元方程求解得來，若能求得結構邊界表面上的聲壓分布即可對聲場內任意場點的聲壓進行求解，在結構表面處，離散形式的聲壓[Pf]由下式計算：

式中，[A]和[B]為K×K階復系數方陣邊界元法影響系數矩陣，為激勵頻率的函數，與結構表面尺寸形狀和插值函數有關，vn為結構表面節點法向速度。在結構邊界表面處，由于聲場的變化會引起流體對該結構處節點動壓力[Fp]，此時結構邊界處節點力與流體動壓力是等值方向的，因而流體動壓力可表示為：

式中，[G]為法矢量方向余弦轉換陣，[N]為邊界結構形狀函數矩陣。則有：

對于軸系-基座-殼體系統，主要是研究系統在穩態激勵作用下軸系傳遞到殼體上的功率流特性。因此穩態激勵下，功率流計算式如下：

則結構表面處第i個節點功率流法向分量可表達為：

式中，為結構表面節點法向位移共軛矩陣。設結構表面有s個節點，則結構表面總功率流

對于本文所研究的系統，從軸系到殼體及結構表面其結構是確定的，因此[G][S]是一個常量，流入到結構表面的法向功率流主要是由螺旋槳激振力、電機激振力以及半聯軸器和高彈不對中產生的激振力引起，傳遞的主要路徑為后尾軸承、前尾軸承、推力軸承以及電機隔振器。假設軸系上功率流Ps傳遞到軸承基座時功率流傳遞為η1，軸承基座上功率流Pb傳遞到殼體表面時功率流傳遞率為η2，則有：

由式（10）可知，當軸承基座到殼體表面結構確定時，則傳遞率η2確定，則減少流經傳遞路徑的功率流在很大程度上就能減小結構表面的功率流。因此對于該系統結構-聲輻射優化的目標函數先從尾部場點聲壓轉化到結構表面功率流，最終轉換成減少流經傳遞路徑的總功率流，即最終優化目標是流經傳遞路徑的總功率流最小。這樣處理具有許多優點：將聲壓（矢量）參數求解轉化成功率流（標量）求解，對求解結果評估時可不用考慮場點位置等因素；大大減少計算量，因為該系統由于尾部結構表面劃分成許多單元和節點，編號無序，這就增加了提取節點信息的難度，并且需要計算每個節點的功率流，計算量較大；避免了同時采用有限元和邊界元之間轉化時造成的誤差等問題。

對于第R個傳遞路徑有I個評價點，則流經該路徑的總功率流評價公式為：

則適應度函數為：

式中：R為傳遞路徑數量R=4，對于后尾軸承和電機支撐軸承分別為I=5和I=4。對于前尾軸和推力軸承，考慮到兩者上還需安裝動力吸振器等因素，因此不能從軸承支撐處來評價，需從兩者與殼體連接處的功率流來進行評價，對于前尾軸承，其連接處如圖2和圖3所示。對于這2處的功率流計算方法參考文獻[6]。前尾軸承處共60個單元，評價節點共計I=224，推力軸承連接處為24個單元，評價節點共計I=75。由適應度函數可知，當流入傳遞路徑的總功率流越小，則適應度越高，種群就會聚集下來，最終找到最合適的種子。

圖 2 前尾軸承連接處Fig. 2 Former stern bearing connected

圖 3 推力軸承連接處Fig. 3 Thrust bearing bearing connected

2.3 參數設置

初始參數的設置對優化結果有一定的影響，特別是野草種群中野草數量的最大值Smax、非線性指數n，種子散布的標準差初始值σinit和標準差最終值σfinal。Smax并不是越大越好，太大反而影響IWO算法的效率，一般選取10～20之間的整數最佳，這里設定Smax=20；非線性指數n一般設置為3；標準偏差初始值應盡量大一些，以保證算法在迭代初期的搜索能力，一般設置為每個變量的解空間長度的1%～5%。標準偏差最終值影響算法的局部尋優能力，減小該數值可提高尋優的精度，但該值過小，會影響算法的收斂速度，且對精度提高不大，因此設置應當適中。表1給出了IWO算法參數的詳細設置。

表 1 IWO算法的參數設置Tab. 1 Parameters of IWO algorithm

2.4 優化過程

本文所用程序均在Matlab中進行，在適應度函數中需要調用Ansys進行頻率響應計算，將傳遞路徑上相關節點的力和位移信息導出，用于功率流計算。進行頻率響應分析時，模型計算參數如下：螺旋槳、半聯軸器、高彈聯軸器及電機轉子的極轉動慣量分別為3 250 kg·m2，125 kg·m2，365 kg·m2和 2 550 kg·m2；軸系密度、泊松比、彈性模量和切變模量分別為7 850 kg/m3，0.26、2.1E11 Pa和0.769E11 Pa；高彈聯軸器扭轉剛度kn=2.6E6 N·m/rad，徑向剛度為7.8E6 N/m；殼體分雙層殼體；軸系轉速為70 r/min。相對應于該工況，分別計算支撐軸承的油膜動力特性系數和推力塊上油膜動力特性系數[11]，代入系統模型進行計算；計算所施加外載荷有螺旋槳激振力6個分量、電機驅動力矩、高彈和半聯軸器不對中產生的激振力。由于軸系不對中和螺旋槳重力造成軸系呈彎曲狀態，因而軸系在旋轉時產生陀螺效應，為此本文先進行靜態計算，再對軸系運轉進行有預應力的頻率響應分析，這樣就可以將軸系變形時產生的彎曲應力和位移考慮到動態結果中，得到軸系在彎曲狀態下的結果。優化過程見圖4。

3 系統結構-聲輻射優化分析

經迭代優化計算，得到相應結果。圖5給出了適應度評價函數值隨迭代次數的變化曲線。算法在約第172代基本達到收斂，且在進化初期具有較快的上升速度，能夠避免陷入局部收斂，表明算法的收斂性能良好，可以尋求全局最優。圖6為迭代過程與傳遞路徑總功率流的變化曲線，反應了迭代過程中每個種子生成的總功率流。表2為最合適種子，即參數優化結果，共計18個參數。

將優化后的參數代入到系統模型進行計算，得到相應結果。圖7為流經傳遞路徑的總功率流與殼體尾部某場點聲壓頻率響應對比圖。由圖可知：優化后流經傳遞路徑的總功率流與場點聲壓頻率響應特性從整個頻率段來看大體一致，但在中高頻處某些點響應規律也不盡相同，這是由于流經傳遞路徑的總功率流不僅僅是傳遞到殼體尾部，而是傳遞到整個殼體，且在這些中高頻處，軸系與殼體的耦合振動往往是以整個殼體與軸系的共振響應為主，因此響應規律略有不同；而在中低頻處，系統振動主要是尾部結構和軸系的耦合振動，場點也是選在尾部附近，因此兩者響應規律基本一致。由此可見，通過流經傳遞路的總功率流來間接反映場點聲壓的響應規律可行。由圖7還可知，優化后流經傳遞路徑的總功率流有所減少，尤其是在頻率19 Hz，27 Hz和34 Hz處，減少最為明顯，但是，在頻率為9 Hz，23 Hz和49 Hz處有所增大，對于中高頻段其總功率流在優化后略有減小。

圖 4 IWO算法優化過程Fig. 4 Optimization procedure of IWO algorithm

圖 5 適應度評價曲線Fig. 5 Fitness value evaluation curve

圖 6 迭代過程Fig. 6 Iterative process

表 2 設計參數優化值Tab. 2 The value of design variable

圖 7 功率流和聲輻射對比圖Fig. 7 Compare power flow and acoustic radiation

為驗證圖7中優化前后功率流的結果，采用聲學邊界元法，計算優化前后尾部某場點聲壓頻率響應進行對比。由圖8可知，在頻率19 Hz，27 Hz和34 Hz處，場點聲壓減小最為明顯，尤其是對主要聲壓做貢獻的頻率點19 Hz聲壓減小幅度最大。而在頻率為9 Hz，23 Hz和49 Hz處場點聲壓有所增大，在23 Hz處增幅最大，對于中高頻段其總功率流在優化后略有減小。對于這些聲壓增大的頻率點，主要原因是系統尾部結構復雜，使得其共振頻率點較多，頻帶較寬，而對于減振器、吸振器和隔振器本身減振頻帶較窄，加上本文對同一類吸振器和隔振器均采用相同頻率吸振，這就造成整體減振頻帶不寬，會導致某些頻率響應值不減反而增大。但是對于主要頻率段還是能起到良好減振效果。通過計算對比分析，以流經傳遞路徑的總功率流為目標函數，通過本文所采取的優化措施，能夠有效地對軸系-基座-殼體這一復雜系統進行聲輻射優化。若能采用多個不同參數的吸振器和隔振器進行優化，則能增加減振降噪的效果。

圖 8 場點聲壓頻率響應對比圖Fig. 8 Compare the frequency response of filed point acoustic pressure

為更清楚了解本文采用的優化方法對系統減振降噪的效果，將頻率為19 Hz時優化前場點聲壓云圖9和優化后場點聲壓云圖10。由圖9可見，聲壓最高處分布在殼體尾部和殼體中部范圍，峰值接近200 dB，在殼體中部水平方向上聲輻射范圍最大最廣。由圖10可見，經過優化后，殼體尾部和殼體中部聲輻射明顯減弱，峰值降到187 dB以內，尤其是殼體中部聲輻射減弱最為明顯。由此可見本文采用的優化計算方法有效可行。

圖 9 19 Hz優化前場點聲壓云圖Fig. 9 19 Hz contours of filed point acoustic pressure before optimization

圖 10 19 Hz優化后場點聲壓云圖Fig. 10 19 Hz contours of filed point acoustic pressure after optimization

4 結 語

本文以流經傳遞路徑的總功率流為目標函數，采用IWO算法對縱向減振器、吸振器和隔振器的參數為設計變量，對軸系-基座-殼體系統進行聲輻射優化，通過計算分析得到如下結論：

1）入侵性野草算是一種簡單有效的優化方法，具有較好的適應性和隨機性，能夠很好的全局尋優。

2）減少流經傳遞路徑的總功率流能夠有效地減少殼體表面聲輻射。

3）以流經傳遞路徑的總功率流為目標函數能夠大大減少計算量和簡化問題，結合功率流有限元法和聲學邊界元法，從而為軸系-基座-殼體提供有效的優化和計算方法。

4）對于減振元件的參數設計應該多樣化以增加減振頻帶，從而更有效地對軸系-基座-殼體減振降噪。

[1] 伍先俊, 程廣利, 朱石堅. 最小振動功率流隔振系統ANSYS優化設計[J]. 武漢理工大學學報(交通科學與工程版), 2005,29 (2): 186–189.

[2] 伍先俊, 朱石堅. 組件功率流計算法和iSIGHT 環境下隔振系統優化設計[J]. 船舶力學, 2006, 10 (2): 138–145.WU Xian-jun, ZHU Shi-jian. Vibration power flow calculation based on component modal technique and isolation system optimization using iSIGHT[J]. Journal of Ship Mechanics,2006, 10 (2): 138–145.

[3] 賀云南, 何琳, 呂志強, 等. 功率流有限元法結果分析輻射噪聲[J]. 船舶力學, 2006, 10 (5): 150–154.HE Yun-nan, HE Lin, LV Zhi-qiang, et al. Radiation noise analysis using the results of power flow finite element method[J]. Journal of Ship Mechanics, 2006, 10(5): 150–154.

[4] 趙群, 張義民, 趙晉芳. 振動傳遞路徑的功率流傳遞度靈敏度分析[J]. 振動與沖擊, 2009, 28(7): 183–186.ZHAO Qun, ZHANG Yi-min, ZHAO Jin-fang. Sensitivity analysis of powerflow transfer probability for a vibration transfer path[J]. Journal of Vibration and Shock, 2009, 28(7):183–186.

[5] 陳爐云, 張裕芳. 基于功率流分析的結構聲優化研究[J]. 振動與沖擊, 2010, 29(10): 191–194.

[6] 楊德慶, 羅放, 陳靜. 有限元功率流落差計算方法研究[J]. 噪聲振動與控制, 2009, 6: 127–131.YANG De-qing, LUO fang, CHEN Jing. Power flow level difference finite element analysis[J]. Noise and Vibration Control, 2009, 6: 127–131.

[7] 肖功煜, 劉微, 朱翔. 船舶板殼結構的振動功率流分析[J]. 船海工程, 2011, 40(6): 45–48.XIAO Gong-yu, LIU Wei, ZHU Xiang. Structure vibration and energy flow characteristics of plates and shells of ships[J]. Ship& Ocean Engineering, 2011, 40(6): 45–48.

[8] FOUDAZI A, A R MALLAHZADEH, Pattern synthesis for multi-feed reflector antennas using invasive weed optimisation[J]. Microwaves, Antennas & Propagation, IET,2012. 6(14): 1583–1589.

[9] MEHRABIAN A R, LUCAS C. A novel numerical optimization algorithm inspired from weed colonization[J].Ecological Informatics, 2006, 1(3): 355–366.

[10] DASTRANJ A, H ABIRI, A MALLAHZADEH. Design of a broadband cosecant squared pattern reflector antenna using IWO algorithm[C]// Antennas and Propagation, IEEE Transactions on, 2013. 61(7): 3895–3900.

[11] 石磊. 計入支承系統特性的船舶推進軸系動態校中研究[D].大連: 大連理工大學, 2010.