基于遺傳算法的圓柱空腔吸聲覆蓋層參數優化研究

陶 猛， 趙 陽， 王廣瑋

(1.貴州大學 機械工程學院，貴陽 550025；2.上海交通大學 機械系統與振動國家重點實驗室, 上海 200240)

在水下航行器表面敷設吸聲覆蓋層既可減小主動聲吶探測回波，亦可降低自身噪聲的外輻射，作為提高水下航行器隱身性能的技術手段被廣泛采用。經半個多世紀發展，吸聲覆蓋層研究已取得豐富成果，如共振空腔吸聲覆蓋層[1-3]、微粒復合吸聲覆蓋層[4-6]、夾芯復合吸聲覆蓋層[7]等。圓柱空腔吸聲覆蓋層為共振吸聲覆蓋層的一種，通常由內部排列的周期性空腔粘彈性介質形成，空腔的基本形狀為圓柱、圓臺或喇叭形等截面且隨厚度變化。對該結構的吸聲覆蓋層研究方法有：① 基于粘彈性介質中波傳播特性的解析計算方法[3,8-9]，其區別在于所得波傳播特性精確[3]或等效[8-9]；② 基于有限元的數值計算方法[1-2,10]。數值計算方法可計算相對復雜結構，但物理意義不明確；解析方法雖物理意義明確，便于分析其本質，但無法計算復雜的聲學結構。

吸聲覆蓋層的機理分析及聲學特性計算雖有較多研究成果，但對吸聲覆蓋層優化設計研究較少。吸聲覆蓋層設計涵蓋材料、結構、工藝等多方面，主要技術為聲學設計，即對吸聲覆蓋層基底材料及空腔結構設計、優化。吸聲覆蓋層優化設計即通過合適的優化算法合理設計吸聲覆蓋層結構尺寸與材料參數，使吸聲覆蓋層達到最優降低反射系數要求。但目前研究成果多集中于單個結構參數或單個材料參數[3,6,10]對吸聲性能影響規律，而多參數同時變化進而改變吸聲性能研究成果較缺乏。由于吸聲覆蓋層聲學性能與材料屬性、結構參數呈非線性復雜關系，不適用經典的搜索算法。隨計算機性能提高及現代優化算法如遺傳算法、模擬退火算法等發展，使建立合理有效的吸聲覆蓋層材料及結構參數綜合優化方法成為可能。

本文對圓柱空腔吸聲覆蓋層的二維理論模型進行簡化，以獲得能快速準確計算吸聲覆蓋層聲學性能的簡化模型。并結合遺傳算法，通過對吸聲覆蓋層基底介質材料參數及空腔結構的優化分析，建立圓柱空腔吸聲覆蓋層材料與結構參數的綜合優化設計方法。

1 吸聲覆蓋層聲學特性簡化模型

圓柱空腔吸聲覆蓋層結構沿厚度方向形成空腔，其基本形狀為圓柱型或圓臺形，見圖1(a)，按正三角形周期性排列相同空腔。因對稱性只取其中一個單元分析。由于建立六面棱柱體理論模型較困難，故用圓柱代替簡化模型的復雜性，見圖1(b)。

圖1 典型圓柱空腔吸聲覆蓋層及單元結構示意圖

湯渭霖等[3,11]建立的圓柱空腔吸聲覆蓋層二維理論將吸聲覆蓋層模型簡化為截面連續變化的多個粘彈性圓柱管組合，其聲學特性由波在圓柱管中傳播、損耗特性決定。但該模型因考慮高階軸對稱波對聲學特性的貢獻，使計算分析過程復雜繁瑣。因此，需基于二維理論尋找合理快速的近似方法，該方法主要基于兩個近似：①平面波垂直入射時，吸聲覆蓋層低頻性能由粘彈性圓柱管中最低階軸對稱波傳播特性決定[11]，進而粘彈性圓柱管等效阻抗由最低階軸對稱波傳播波數確定。②對圓柱-圓臺空腔，近似為截面階梯變化的多個圓柱管組合，見圖1(c)，每段圓柱管有等效阻抗，整個圓柱管可視為多個阻抗變化的傳輸線組合。

平面波垂直入射時，吸聲覆蓋層中軸對稱波傳播特性可通過求解軸對稱波的特征方程獲得，即：

(1)

其中：

(2)

式(1)的詳細建立過程見文獻[3]。對給定的結構參數及材料屬性，求解式(1)可得一組復根Km(m=1,2,…)，并由此確定kl,m，kt,m。此求解在復平面進行，K的實部為傳播常數，虛部反映傳播衰減。

二維理論給出吸聲覆蓋層單元-粘彈性圓柱管中位移、應力表達式[3]：

(3)

其中：RΦ,m為聲波在吸聲覆蓋層后端面(z=0)即背襯反射系數;下標m表示第m階軸對稱波傳播模式，j表示圓柱管截面第j圓環區域，m=j為能計算吸聲覆蓋層反射系數或吸聲系數的必要條件。式(3)中ζmj(V)，ζmj(W)為關于吸聲覆蓋層結構、材料的復雜函數，即：

(4)

其中：αA,m,αB,m,αC,m為式(1)的代數余子式。參數V,W可表示為：

(5)

當僅考慮最低階軸對稱波(m=j=1)時，粘彈性圓柱管的等效阻抗可寫為：

(6)

經推導、簡化，式(6)可寫成簡單形式：

(7)

至此，粘彈性圓柱管前端面(z=H，上標f表示)、后端面(z=0，上標b表示)處位移及應力可通過傳遞矩陣T相聯系，即：

(8)

總傳遞矩陣可寫成每段圓柱管傳遞矩陣的連乘形式：

(9)

當粘彈性圓柱管后端面(背襯)的阻抗Zb已知時，前端面表面阻抗可由式(8)獲得：

(10)

吸聲覆蓋層反射系數可寫成：

(11)

式中：Zw為水介質特性阻抗。

2 基于遺傳算法的吸聲覆蓋層參數優化

由以上推導看出，吸聲覆蓋層反射系數或吸聲系數為關于基底介質材料屬性及結構參數的復雜函數，單純優化吸聲覆蓋層結構參數或材料屬性，難在寬頻范圍內具備較高吸聲系數，而二者同時優化則可能實現寬頻吸聲目標。吸聲覆蓋層結構參數包括穿孔率、空腔高度等；材料參數包括基底介質密度、泊松比、楊氏模量及損耗因子等。而楊氏模量及損耗因子大小與頻率有關，較難用傳統的搜索算法解決。因此本文用遺傳算法進行吸聲覆蓋層多目標參數優化研究。

圖2 遺傳算法流程圖

圖3 圓柱-圓臺空腔吸聲覆蓋層單元結構示意圖

與傳統優化方法相比，遺傳算法、模擬退火算法等在解決大規模非線性優化問題時優勢明顯。其中遺傳算法為借鑒生物界自然選擇、自然遺傳機制的隨機搜索算法，較適用于處理傳統搜索算法難以解決的復雜、非線性優化問題。遺傳算法流程圖見圖2。遺傳算法由隨機產生的初始解開始搜索，通過選擇、交叉、變異操作逐步迭代產生新解。該算法優點即將參數編碼成染色體后進行優化，而不針對參數本身，因而不受函數約束條件限制；搜索過程由解的一個集合開始，而非單個個體，具有隱含并行搜索特性，可大大減少陷入局部最小的可能性，且計算時不依賴于梯度信息，不要求目標函數連續可導，較適合求解傳統搜索算法難以解決的大規模、非線性組合優化問題。

用遺傳算法對圖3吸聲覆蓋層進行優化。該吸聲覆蓋層單元間按正三角形方式排列，據上述簡化方式，每個單元可簡化為半徑a的圓柱體，且每個單元中包含圓柱、圓臺兩段空腔結構，圓柱空腔半徑及圓臺空腔上表面半徑均為b1，圓臺空腔下表面半徑為b2，圓柱空腔及圓臺空腔高度分別為h1，h2，空腔兩端均有厚度為hc的封口薄層，整個吸聲覆蓋層高度為H=h1+h2+2hc。平面波由圓柱空腔一側垂直入射到吸聲覆蓋層，吸聲覆蓋層背襯條件為絕對軟背襯。

用遺傳算法進行優化時，目標函數選擇為：

(12)

據遺傳算法理論，在MATLAB軟件中編寫實現基本遺傳算法程序以尋找目標函數最優解。遺傳算法參數設置：種群規模100，進化次數100，交叉概率0.8，變異概率0.2。

2.1 算例1(材料屬性優化)

只優化吸聲覆蓋層基底介質的粘彈性動態力學參數，即材料楊氏模量及損耗因子。由于粘彈性材料的動態力學參數非常數，隨頻率變化呈色散關系，故假設粘彈性材料的楊氏模量及損耗因子隨頻率成二次函數關系[6]，即：

(13)

式中：f為頻率；E為楊氏模量(Pa)。若只考慮吸聲覆蓋層由一種基底介質構成，需優化6個參數p0～p2，q0～q2。由此可計算不同頻率處的楊氏模量及損耗因子，進而用簡化模型計算反射系數。

經優化的6個參數及設定的覆蓋層單元各結構參數見表1?！?”表示未優化參數(下同)，在該算例中為單元各結構參數。吸聲覆蓋層厚度H=50 mm。

表1 吸聲覆蓋層優化后材料屬性及未優化結構參數

吸聲覆蓋層基底介質的楊氏模量及損耗因子經優化的吸聲系數見圖4，經優化的楊氏模量及損耗因子見圖5。由圖4看出，吸聲覆蓋層最低吸聲頻率約2 kHz，在2～8 kHz頻率范圍內，平均吸聲系數約0.86。

圖4 吸聲覆蓋層材料參數優化后吸聲系數

圖5 優化后楊氏模量及損耗因子

2.2 算例2(材料屬性、結構參數綜合優化)

在算例1基礎上進行材料屬性、結構參數綜合優化。待優化結構參數包括圓柱空腔半徑b1、圓臺下表面空腔半徑b2、圓柱空腔高度h1，圓臺空腔高度可據h2=H-h1-2hc計算獲得，而吸聲覆蓋層厚度H與封口薄層厚度hc設定不變。因此，待優化參數包括與基底材料楊氏模量與損耗因子有關的6個系數、3個結構參數計9個參數。

表2 吸聲覆蓋層優化后材料屬性與結構參數

圖6 吸聲覆蓋層材料及結構參數優化后吸聲系數

圖7 優化后楊氏模量及損耗因子

經優化后9個參數及設定的其它結構參數見表2，吸聲覆蓋層材料及結構參數經綜合優化后吸聲系數見圖6，經綜合優化后基底介質楊氏模量及損耗因子見圖7。由圖6看出，吸聲覆蓋層最低吸聲頻率仍為2 kHz左右，而在 2～8 kHz頻率范圍內，平均吸聲系數約0.89。與算例1比較知，若同時優化吸聲覆蓋層材料及結構參數，可獲得吸聲性能更好的材料屬性與結構參數組合。

2.3 算例3(兩層材料屬性及結構參數綜合優化)

考慮更復雜情況，設吸聲覆蓋層由兩種基底介質構成：包含圓柱空腔部分為第一種介質，包含圓臺空腔部分為第二種介質，即兩種介質以圓柱空腔及圓臺空腔的交界處為邊界。需優化的參數包括與第一種材料楊氏模量及損耗因子相關的6個系數(以上標“(1)”表示)、與第二種材料楊氏模量及損耗因子相關的6個系數(以上標“(2)”表示)、3個結構參數計15個待優化參數。

表3 吸聲覆蓋層優化后兩層材料屬性及優化后結構參數

優化后15個參數、設定的其它結構參數見表3，吸聲覆蓋層,兩層材料屬性及結構參數經綜合優化后吸聲系數見圖8，優化后第一層介質楊氏模量及損耗因子見圖9，優化后第二層材料楊氏模量及損耗因子見圖10。由圖8看出，吸聲覆蓋層最低吸聲頻率略高于2 kHz, 在2～8 kHz頻率范圍內，平均吸聲系數約0.93。與前兩種優化模型相比，平均吸聲系數明顯提高，原因為考慮吸聲覆蓋層分層后，吸聲覆蓋層實現阻抗漸變形式，使波在多層結構中傳播時遭遇阻抗變化并發生多次反射、透射，模式轉換可產生更多的能量損耗，從而提高寬頻吸聲能力。而阻抗漸變形式不僅可通過結構分層實現，也可通過材料分層實現。

圖8 吸聲覆蓋層材料及結構參數優化后吸聲系數

3 結 論

(1) 圓柱空腔吸聲覆蓋層材料屬性與結構參數存在一組使吸聲覆蓋層在工作頻段內具有高吸聲系數的最佳值，該值可由參數優化方法獲得。

(2) 較多需優化參數尤其粘彈性材料動態力學參數(楊氏模量、損耗因子)為頻率函數，常規搜索算法無法解決此類大規模非線性優化問題。

(3) 本文建立的圓柱空腔吸聲覆蓋層性能分析簡化模型，通過用遺傳算法進行吸聲覆蓋層材料屬性及結構參數優化分析，仿真算例表明，綜合優化材料屬性及結構參數可獲得更佳的寬頻吸聲性能。

(4) 本文所建優化模型可進行更多參數優化分析，亦為吸聲覆蓋層優化設計提供可行計算方法。

參 考 文 獻

[1]Hennion A C, Decarpigny J N. Analysis of the scattering of a plane acoustic wave by a doubly periodic structure using the finite element method: application to alberich anechoic coatings [J]. Journal of the Acoustical Society of America, 1991, 90(6):3356-3367.

[2]Easwaran V, Munjal M L. Analysis of reflection characteristics of a normal incidence plane wave on resonant sound absorbers: a finite element approach [J]. Journal of the Acoustical Society of America, 1993, 93(3):1308-1318.

[3]湯渭霖, 何世平, 范 軍. 含圓柱形空腔吸聲覆蓋層的二維理論[J]. 聲學學報, 2005, 30(4):289-295.

TANG Wei-lin, HE Shi-ping, FAN Jun. Two-dimensional model for acoustic absorption of viscoelastic coating containing cylindrical holes [J]. Acta Acustica, 2005, 30(4): 289-295.

[4]Zhao H G, Liu Y Z, Wen J H, et al. Tri-component phononic crystals for underwater anechoic coatings[J]. Physics Letters A, 2007, 367(3): 224-232.

[5]Ivansson S M. Numerical design of alberich anechoic coatings with superellipsoidal cavities of mixed sizes[J]. Journal of the Acoustical Society of America, 2008, 124(4): 1974-1984.

[6]王仁乾, 馬黎黎. 吸聲材料的物理參數對消聲瓦吸聲性能的影響[J]. 哈爾濱工程大學學報, 2004, 25(3): 288-294.

WANG Ren-qian, MA Li-li. Effects of physical parameters of the absorption material on absorption capability of anechoic tiles[J]. Journal of Harbin Engineering University, 2004, 25(3): 288-294.

[7]朱 錫, 羅 忠, 周 欣,等. 斜入射下水中隱身夾芯復合材料殼板結構聲學設計[J]. 船舶力學, 2009, 13(4): 646-656.

ZHU Xi, LUO Zhong, ZHOU Xin, et al. Acoustic stealth design on the underwater sandwich composite shell structure in oblique incidence [J]. Journal of Ship Mechanics, 2009, 13(4): 646-656.

[8]何祚鏞, 王 曼. 水下非均勻復合層結構吸聲的理論研究[J]. 應用聲學, 1996, 15(5):12-19.

HE Zuo-yong, WANG Man. Investigation of the sound absorption of non-homogeneous composite multiple-layer structures in water [J]. Applied Acoustics, 1996, 15(5): 12-19.

[9]朱蓓麗, 任克明. 等效參數法研究帶圓柱通道橡膠體的聲學性能[J]. 上海交通大學學報, 1997, 31(7): 20-24.

ZHU Bei-li, REN Ke-ming. Acoustic properties analysis of elastomer with cylindrical bores using equivalent parametric method [J]. Journal of Shanghai Jiaotong University, 1997, 31(7): 20-24.

[10]Panigrahi S N, Jog C S, Munjal M L. Multi-focus design of underwater noise control linings based on finite element analysis[J]. Applied Acoustics, 2008, 69(12): 1141-1153.

[11]何世平, 湯渭霖, 何 琳,等. 變截面圓柱形空腔覆蓋層吸聲系數的二維近似解[J]. 船舶力學, 2006, 10(1): 120-127.

HE Shi-ping, TANG Wei-lin, HE Lin, et al. Analysis of acoustic characteristics of anechoic coating containing varying sectional cylindrical cavity [J]. Journal of Ship Mechanics, 2006, 10(1): 120-127.