基于有限元與能量方法結(jié)合的前圍板阻尼片布置優(yōu)化

常光寶，張吉健，羅德洋，白長安

(1.上汽通用五菱汽車股份有限公司技術(shù)中心，廣西 柳州 545007；2.武漢融聲奇科技有限公司，武漢 430014；3.西安交通大學(xué) 機械工程學(xué)院，西安 710049)

0 引言

阻尼片在汽車行業(yè)使用廣泛，這些阻尼材料通常由一些柔性材料制作而成(如合成橡膠、橡膠、軟木、高密度泡沫和層壓材料)，被貼裝在板件上，用以減小車身板件的振動。正因為這樣，阻尼片能夠改善金屬薄板的隔聲效果。優(yōu)化這些布放位置不僅能夠增強隔聲效果，還能契合成本控制、輕量化的發(fā)展需求。

數(shù)值計算已經(jīng)成為除理論分析和實驗測試外的常用研究方式[1]，通過有限元的數(shù)值計算方法研究結(jié)構(gòu)隔聲特性建立在經(jīng)典力學(xué)基礎(chǔ)之上，主要包括有限元法和邊界元法等。王智文[2]、張若平[3]等對汽車地板阻尼材料與車內(nèi)振動噪聲性能的聯(lián)系進行了研究，并對車身結(jié)構(gòu)的模態(tài)應(yīng)變能進行分析，優(yōu)化阻尼，提高NVH性能。趙建軒[4]、鄭玲[5]等進一步結(jié)合壁板貢獻量，降低振動噪聲。梅勝軍等[6]將模態(tài)應(yīng)變能、壁板貢獻量和模態(tài)貢獻量綜合考慮優(yōu)化阻尼。張一麟等[7]通過觀察車身模態(tài)并結(jié)合壁板貢獻度的方法優(yōu)化阻尼。除了上述的有限元建模外，統(tǒng)計能量方法也是常用的數(shù)值計算方法[8]，劉春蕾等[9]將60 km/h勻速行駛工況作為計算工況，對司機耳旁聲壓級進行預(yù)測，驗證了統(tǒng)計能量分析方法預(yù)測噪聲的有效性。統(tǒng)計性方法建立在統(tǒng)計力學(xué)基礎(chǔ)之上，其主要代表有統(tǒng)計能量法等。黃鼎友等[10]利用統(tǒng)計能量法對前圍板的隔聲性能進行了預(yù)測，并與實驗測試的結(jié)果進行了對比，兩者趨勢基本一致。但是統(tǒng)計能量法基于統(tǒng)計思想，很難考慮模型細節(jié)的影響，因此在計算精度方面有著天然的劣勢，主要用于趨勢的預(yù)測。

在統(tǒng)計能量分析中，模型分子系統(tǒng)或者分塊的概念成為近年來的分析熱點。相對于單純的基于網(wǎng)格的有限元方法，分塊方法可把若干網(wǎng)格組合成一塊，這樣整個系統(tǒng)分為若干塊，相對于原始網(wǎng)格來說，分塊的自由度會降低很多，分析效率會提升很多?；谀芰刻幚淼姆椒梢詣?chuàng)建更緊湊的能量模型[11-14]，這種方法基于模態(tài)，并可以自動將結(jié)構(gòu)劃分成多塊，分別對每一塊計算整體的分布矩陣。這種方法將原來的有限元模型縮減到對每一塊的計算上，可以快速地進行指定載荷下的能量評估。

本文采用有限元模型和能量后處理相結(jié)合的方法，通過建立更精簡的能量模型描述結(jié)構(gòu)部件的動力學(xué)響應(yīng)，并通過編程的方式實現(xiàn)自動計算每個單元塊對阻尼的敏感性響應(yīng)結(jié)果。這種方法基于現(xiàn)有的結(jié)構(gòu)模態(tài)和剛度質(zhì)量矩陣，模型中的各結(jié)構(gòu)部件自動劃分成很多單元片，并且每個單元片上進行“分布矩陣”組裝。該方法可以大幅降低模型計算量，有效提高數(shù)值計算效率；同時基于其高效性，可以把該方法集成到優(yōu)化工作流程中。

1 理論介紹

有限元方法與能量分析方法結(jié)合的流程中，基于結(jié)構(gòu)動力學(xué)的有限元方法不贅述。能量分析方法是基于有限元結(jié)構(gòu)模態(tài)與模型的剛度-質(zhì)量矩陣(K_M矩陣)來進行分塊處理和能量分析，該過程稱為能量后處理過程[15]。

能量后處理過程中，首先將有限元模型分解成一系列單元集合(Patch，下文稱“單元塊”)。每一個單元塊是一系列屬于相同有限元模型部件的連續(xù)的單元集合，且假設(shè)在每一個單元塊上能量級均勻分布。因此，單元塊是這種能量后處理方法的基本構(gòu)成。圖1展示了自動劃分的某示意板件的單元塊，每個小單元塊里面有若干有限元單元網(wǎng)格。

圖1 單元塊自動劃分示意圖

(1)

式中：Ki、Mi和Kip、Mip分別為模態(tài)坐標系和物理坐標系下的第i個單元塊上的剛度陣和質(zhì)量陣。分布式矩陣Ki和Mi只需計算1次，后續(xù)各種載荷作用下的響應(yīng)計算均調(diào)用。

計算單元塊的分布式矩陣優(yōu)勢在于可以在單元塊的子結(jié)構(gòu)上快速地進行能量計算。也就是說，在各種載荷工況下都可以計算出第i個單元片的動能、勢能和耗散功率：

(2)

(3)

(4)

式中：Y為模態(tài)坐標位移，Ti、Vi、Wi分別為單元的動能、勢能和耗散功率，ηi為阻尼損耗因子，ω為圓頻率。

當(dāng)只有遲滯阻尼決定阻尼片的能量耗散時，可以通過阻尼損耗因子ηi來定義第i個單元塊上的阻尼。在每個單元塊上定義不同的阻尼因子，就會形成不均勻的阻尼分布。經(jīng)過循環(huán)計算后，可以計算得到哪些單元塊對阻尼片較為敏感。

板件的隔聲特性仿真分析通常需要用擴散聲場作為載荷來模擬無規(guī)入射的聲能量。擴散聲場可以通過多個采樣的不相關(guān)的、在空間中沿著任意方向傳播的平面波疊加獲得，可以用來模擬一個弱穩(wěn)態(tài)隨機過程，描述這一理論模型的解析表達式如下：

(5)

(6)

式中：k為波數(shù)，c為聲速，d為聲場中2個觀察點之間的距離。

2 技術(shù)路線

采用數(shù)值分析軟件與自編程序相結(jié)合的方式進行前圍板阻尼片位置優(yōu)化方法研究，研究對象為含鈑金件和復(fù)合層的前圍板結(jié)構(gòu)，其中板件的有限元模型已經(jīng)確定，采用基于有限元方法與能量方法的敏感性分析流程確定阻尼片的布置位置，技術(shù)路線如圖2所示。

圖2 技術(shù)路線框圖

具體實施所用方法和內(nèi)容如下所示：

1)在有限元分析軟件Nastran中，采用有限元分析的方法進行車輛前圍板的有限元模態(tài)分析；

2)在聲學(xué)分析軟件Actran中，采用能量分析方法建立能量分析模型，在輸入模態(tài)結(jié)果和剛度質(zhì)量矩陣(KM矩陣)數(shù)據(jù)后，設(shè)定N個有限元網(wǎng)格組成一個單元塊子系統(tǒng)。分析模型中考慮擴散聲場激勵加載到結(jié)構(gòu)模態(tài)上，擴散聲場采用25個抽樣樣本構(gòu)成隨機的分布式載荷激勵；

3)使用Python語言編寫板件不同單元塊對阻尼的敏感性分析程序。因為結(jié)構(gòu)阻尼被認為在各個單元塊上是相同的，優(yōu)化流程最終會得到在各個單元塊上阻尼材料的分布量。通過監(jiān)測擴散聲場載荷下的板的法向均方速度來判定板件不同區(qū)域?qū)ψ枘岬拿舾行裕徽业絥個阻尼片的最優(yōu)分布方案：

① 對于每一個單元塊，應(yīng)用一個局部阻尼值；

② 計算對于每一種阻尼材料分布方案的整個板件法向速度的均方值；

③ 每一步都進行針對待優(yōu)化板件的白車身有限元模型的能量后處理。使用分布矩陣計算每一種阻尼優(yōu)化方案的殘余矩陣。

4)挑選出n個單元塊的法向振動速度均方值最小的位置添加阻尼。

在采用傳統(tǒng)直接頻響法進行每一步計算時，模型自由度由單元數(shù)量決定。而能量后處理法中，模型自由度由單元塊數(shù)量決定。文中優(yōu)化流程的每一步計算在采用能量后處理方法時，計算模型的自由度都被降低到最少，計算效率更高。在最終驗證計算時，采用直接頻響法，以確保優(yōu)化結(jié)果的準確性。

3 前圍板隔聲量性能研究及聲學(xué)優(yōu)化案例

3.1 前圍板隔聲量性能分析

前圍板隔聲量計算流程如圖3所示。隔聲量分析模型中，定義前圍板結(jié)構(gòu)各部分材料參數(shù)，并導(dǎo)入一側(cè)聲空間模型，定義為空氣材料屬性；模型中通過定義耦合面的方式把結(jié)構(gòu)模型和空腔模型耦合起來，保證信息可以正常傳遞。模型中定義擴散聲場激勵加載到結(jié)構(gòu)模型上，擴散聲場采用25個抽樣樣本構(gòu)成隨機的分布式載荷激勵。

圖3 前圍板隔聲量計算流程框圖

隔聲量計算模型中，前圍板結(jié)構(gòu)有限元模型與聲空間有限元模型組合成前圍板隔聲量計算聲振耦合分析模型，模型組成示意圖如圖4所示。

圖4 隔聲量計算-聲振耦合模型組成示意圖

圖4中，①為擴散聲場激勵，模擬無規(guī)律入射聲場激勵；②為APML邊界，模擬無反射的聲場邊界。前圍板鈑金件和阻尼材料參數(shù)如表1所示，阻尼材料的阻尼隨頻率變化。

表1 材料參數(shù)

隔聲量計算完畢，得到隔聲量曲線，為了驗證仿真計算精度和方法的有效性，取相同模型的實際樣件做隔聲量測試。前圍板隔聲測試試驗在同濟大學(xué)聲學(xué)研究所進行，測試方法按照GB/T 19889.10—2006小建筑構(gòu)件空氣聲隔聲的實驗室測試方法進行，隔聲室由兩間相鄰的混響室組成，如圖5所示。試驗過程中前圍板試驗件裝置在公共隔墻的開口上，測試溫度為16 ℃，相對濕度為80%。試驗采用B&K公司的PULSE 3560C作為聲源的信號發(fā)生器及數(shù)采系統(tǒng)，激勵信號為白噪聲信號。

圖5 前圍板隔聲量測試實驗示意圖

通過3次測試獲得的帶阻尼片前圍板模型的隔聲量數(shù)據(jù)如圖6中紅色實線所示；圖中藍色虛線為采用專業(yè)聲學(xué)軟件計算得到的帶阻尼片前圍板模型的隔聲量數(shù)值模擬數(shù)據(jù)。

圖6 帶阻尼片前圍板隔聲仿真與測試結(jié)果曲線

從圖6中對比情況來看，兩者曲線走向趨勢一致，吻合度較高，其中在80 Hz以下頻段，隔聲計算值與試驗值差距較大，這是由于低頻一般受邊界條件影響較大，而在仿真模型中難以施加與試驗狀態(tài)完全一致的邊界條件；在80～2 500 Hz，隔聲計算值與試驗結(jié)果趨勢基本一致，計算模型滿足工程計算的要求。曲線均在100 Hz處(曲線起始頻率)隔聲量最小，符合經(jīng)典板件隔聲量曲線中結(jié)構(gòu)共振區(qū)隔聲量最低的趨勢，100～350 Hz為隔聲量增大區(qū)域。350～450 Hz隔聲量逐漸降低，達到吻合頻率區(qū)，該區(qū)間內(nèi)入射聲波的波長在壁板上的投影與壁板中彎曲波的波長相近或相等時，發(fā)生吻合效應(yīng)。450 Hz之后的隔聲量曲線呈現(xiàn)逐漸增大的趨勢，且實驗與數(shù)值仿真吻合度較高。計算模型實現(xiàn)了2 500 Hz以內(nèi)頻段的隔聲效果，為阻尼片在地板高頻段的NVH性能表現(xiàn)提供了預(yù)測手段。

3.2 基于能量后處理的阻尼敏感性分析

前圍板隔聲量計算模型搭建完成且完成隔聲量對標后，在后續(xù)敏感性分析和位置優(yōu)化中作為初始輸入的基準模型。在聲學(xué)分析軟件Actran中，采用能量分析方法建立能量分析模型，分析模型基于3.1章節(jié)隔聲量計算模型，在輸入模態(tài)結(jié)果和剛度質(zhì)量矩陣(KM矩陣)數(shù)據(jù)后，設(shè)定N個有限元網(wǎng)格組成一個單元塊子系統(tǒng)。其中前圍板結(jié)構(gòu)模態(tài)結(jié)果如圖7所示，圖中展示前6階結(jié)構(gòu)模態(tài)的模態(tài)位移振型云圖。

圖7 前圍板結(jié)構(gòu)有限元模態(tài)位移振型云圖(前六階模態(tài)振型)

其中，第一、二、五階顯示局部模態(tài)振型；第三、四、六階顯示前圍板板件整體模態(tài)振型。在做能量分析時，可以設(shè)定N個有限元網(wǎng)格為一個單元塊，若有限元網(wǎng)格數(shù)為S個，最終可以生成的單元塊數(shù)即S/N個。前圍板模型中，有限元網(wǎng)格數(shù)為20 945個，如每個單元塊包含800個有限元網(wǎng)格，則最終生成25個單元塊；為了測試該方法的可靠性，研究了每個單元塊中包含200個有限元網(wǎng)格、400個有限元網(wǎng)格、800個有限元網(wǎng)格，共3種計算模型。加載擴散聲場激勵后，完成能量分析模型。在Python運行程序里設(shè)定阻尼循環(huán)加載腳本，最終獲得每個單元塊上的平均振速結(jié)果，生成每個單元塊對阻尼的敏感性云圖，如圖8所示。

圖8的計算結(jié)果中，(a)(b)(c)中左側(cè)圖片為安裝不同分組方式形成的不同網(wǎng)格數(shù)的單元塊數(shù)的分布圖，其中，(a)中每塊單元塊的網(wǎng)格數(shù)最少，所以單元塊數(shù)最多，(b)次之，(c)中單元塊數(shù)最少。3個分組圖中，右側(cè)圖片為加載擴散聲場激勵后獲得的敏感性分析云圖，圖中顏色代表每塊單元塊上加載大阻尼系數(shù)時整個板件的法向均方振動速度，藍色表示加載到該區(qū)域單元塊后，整個板件的振動較?。煌?，紅色區(qū)域表示該區(qū)域加阻尼系數(shù)后整個板件的振動較大。

圖8 有限元不同分塊形式云圖

最終在方案選型階段粘貼阻尼片時，應(yīng)主要粘貼在藍色區(qū)域，同上考慮到阻尼片質(zhì)量和工藝加工的要求，從藍色區(qū)域逐次過渡到淺藍色、綠色區(qū)域、褐色區(qū)域、淺紅區(qū)域，最后再選擇深紅區(qū)域。從圖8的3種分塊方案來看，結(jié)果穩(wěn)定性較好，均在右上區(qū)域呈現(xiàn)藍色且圍繞中間孔洞圓周分布，趨勢線較為一致?？紤]到分塊越少，計算效率越高，后續(xù)分析時，采用圖8(c)方案即可。

3.3 基于敏感性分析的局部阻尼片優(yōu)化

乘員艙內(nèi)含大量隔吸聲材料，且吸聲材料在中高頻段的降噪效果更為明顯，改進空間較大；相對而言，在中低頻范圍內(nèi)的隔聲性能，特別是如何提升100～500 Hz的隔聲量特性尤為關(guān)鍵。本案例中前圍板上原始阻尼片厚度均為0.002 m，總質(zhì)量為0.85 kg，設(shè)定阻尼片增質(zhì)不超過15%，優(yōu)化阻尼片位置及其厚度。位置優(yōu)化依托敏感性分析結(jié)果，從藍色區(qū)域開始，選取的區(qū)域內(nèi)按原始厚度0.002 m，其總質(zhì)量滿足原始阻尼片約90%的質(zhì)量0.78 kg，并把此作為后續(xù)阻尼厚度優(yōu)化的初始模型。

在阻尼厚度優(yōu)化時，調(diào)用聲學(xué)軟件中自帶的NLOPT優(yōu)化算法，以阻尼片厚度為優(yōu)化變量，以隔聲量的平均值為優(yōu)化目標，迭代20次，計算得到隔聲量最大的阻尼厚度方案。圖9為迭代后模型與原始模型的隔聲量曲線。

圖9 阻尼片厚度優(yōu)化前后隔聲量曲線

圖9中，藍色虛線為初始模型隔聲量曲線，紅色實線為迭代優(yōu)化后模型的隔聲量曲線；2條曲線整體走勢一致，優(yōu)化后模型在140～190 Hz隔聲量優(yōu)于原始模型隔聲量。在幾個特征峰值頻率處，如190、230、250、340和350 Hz等頻率處，優(yōu)化后模型隔聲量優(yōu)于原始模型隔聲量，驗證了該方法對于優(yōu)化阻尼片敷設(shè)方案的可行性，同時說明了合理的阻尼片敷設(shè)方案對于提升隔聲量的重要性。

4 結(jié)論

1)基于現(xiàn)有結(jié)構(gòu)模態(tài)和剛度質(zhì)量矩陣的能量分析模型，在每個單元塊上進行“分布矩陣”組裝，并針對具體工況載荷進行計算，將板件結(jié)構(gòu)進行分塊處理，并通過腳本實現(xiàn)在每一塊板件施加阻尼的計算方式，找到板件對阻尼最敏感的位置。在敏感區(qū)域通過自動厚度優(yōu)化計算，找到各區(qū)域最合適的厚度，能夠有效地降低阻尼材料的使用量，減重降本。

2)借助聲學(xué)軟件基于Python語言的腳本開發(fā)，調(diào)用能量后處理方法，對阻尼片有限元建模和實現(xiàn)優(yōu)化功能。僅需要已有的有限元結(jié)構(gòu)模型作為輸入，不需要額外的計算，可非常方便地融入工業(yè)設(shè)計流程。

3)對于車輛整車阻尼片優(yōu)化布置，可基于整車結(jié)構(gòu)模態(tài)和結(jié)構(gòu)剛度質(zhì)量矩陣形成能量分析模型，快速獲得阻尼敏感性分析結(jié)果，且由于考慮真實工況激勵如路噪軸頭力等，可以快速獲得基于真實激勵的整車阻尼優(yōu)化布置方案。對汽車鈑金件類零件的前期設(shè)計有較好的指導(dǎo)意義。