基于Abaqus軟件的輪胎噪聲分析

黃 偉，曹金鳳

（青島理工大學(xué) 理學(xué)院，山東 青島 266520）

隨著汽車(chē)工業(yè)的不斷發(fā)展，人們對(duì)汽車(chē)駕乘舒適性的要求越來(lái)越高。噪聲是汽車(chē)舒適性的重要方面，而輪胎是汽車(chē)重要的噪聲源之一[1]。輪胎噪聲在車(chē)輛高速行駛時(shí)是整車(chē)噪聲的主要部分。研究輪胎噪聲的產(chǎn)生原因并提出降噪措施具有重要的工程意義。

輪胎噪聲的產(chǎn)生機(jī)理十分復(fù)雜，包括空氣動(dòng)力學(xué)機(jī)理、振動(dòng)機(jī)理、增加或減弱機(jī)理等[2]。與空氣動(dòng)力學(xué)機(jī)理相關(guān)的噪聲包括泵浦噪聲、空氣柱共鳴噪聲、亥姆霍茲共振噪聲和空氣湍流噪聲；與振動(dòng)機(jī)理相關(guān)的噪聲包括胎面花紋塊振動(dòng)噪聲、花紋塊粘著振動(dòng)和滑移噪聲；與增強(qiáng)或減弱機(jī)理相關(guān)的噪聲包括喇叭口效應(yīng)、輪胎共振、機(jī)械阻抗效應(yīng)和聲學(xué)阻抗效應(yīng)。目前業(yè)界對(duì)各種噪聲機(jī)制對(duì)輪胎噪聲的貢獻(xiàn)尚未有明確結(jié)論。馮希金等[3]提出了一種泵浦噪聲源識(shí)別的新方法，將試驗(yàn)與數(shù)值分析結(jié)合，通過(guò)對(duì)花紋輪胎噪聲測(cè)試、輪胎聲場(chǎng)阻抗分析和聲源辨識(shí)，利用反演法進(jìn)行試驗(yàn)與仿真，以驗(yàn)證方法的有效性。包秀圖等[4]利用有限元分析得到輪胎表面節(jié)點(diǎn)的振動(dòng)速度，并將其轉(zhuǎn)化為聲學(xué)計(jì)算的頻域邊界條件，再利用聲學(xué)邊界元理論計(jì)算輪胎的低頻振動(dòng)噪聲。

本工作基于非線性有限元分析軟件Abaqus，以385/65R22.5輪胎為例，結(jié)合試驗(yàn)進(jìn)行對(duì)比分析，為低噪聲輪胎的設(shè)計(jì)與開(kāi)發(fā)提供一種新思路。

1 材料模型

輪胎的力學(xué)行為與橡膠材料的性能密切相關(guān)，而輪胎有限元仿真結(jié)果的精度取決于材料參數(shù)的真實(shí)有效性。為了獲得準(zhǔn)確的橡膠材料參數(shù)，在數(shù)值模擬前進(jìn)行了單向拉伸、純剪切和等雙軸拉伸試驗(yàn)。試驗(yàn)在易瑞博科技（北京）有限公司的WDW-10型微控電子萬(wàn)能試驗(yàn)機(jī)上進(jìn)行。為了模擬準(zhǔn)靜態(tài)的拉伸狀態(tài)，試樣加載時(shí)拉伸需十分緩慢，應(yīng)變速率控制在1%·s-1。為了減少橡膠材料的軟化效應(yīng)帶來(lái)的不利影響，3項(xiàng)試驗(yàn)均需要在試驗(yàn)前做5次預(yù)拉伸循環(huán)，以便在試驗(yàn)過(guò)程中獲得穩(wěn)定的拉伸曲線。將3種拉伸曲線作為穩(wěn)定的超彈性本構(gòu)關(guān)系曲線[5]，經(jīng)數(shù)據(jù)分析后得到超彈性本構(gòu)關(guān)系曲線，如圖1所示。

圖1 橡膠的本構(gòu)關(guān)系曲線

2 輪胎的滾動(dòng)仿真分析

2.1 輪胎滾動(dòng)的數(shù)值模擬

輪胎軸對(duì)稱(chēng)有限元模型如圖2所示。首先在Abaqus軟件中建立光面輪胎的二維有限元模型，利用關(guān)鍵詞“symmetric model generation，revolve”將二維模型轉(zhuǎn)換成三維模型；再將三維光面輪胎模型與三維花紋有限元模型進(jìn)行組裝，形成完整的帶花紋輪胎的有限元模型，如圖3所示。

圖2 輪胎二維軸對(duì)稱(chēng)有限元模型

圖3 帶花紋的三維輪胎有限元模型

在Abaqus/Standard求解器中模擬輪胎裝配、充氣和垂直預(yù)加載工況，充氣壓力為600 kPa，垂直預(yù)加載力（Fz）為2 700 kN；通過(guò)import命令實(shí)現(xiàn)將分析結(jié)果由Standard求解器向Explicit求解器的轉(zhuǎn)換，以模擬輪胎的加速及勻速滾動(dòng)，輪胎的滾動(dòng)速度為70 km·h-1，并保證勻速轉(zhuǎn)動(dòng)距離不少于3周，以獲得輪胎滾動(dòng)狀態(tài)的穩(wěn)定分析結(jié)果。

2.2 流-固耦合仿真分析

基于Abaqus軟件的強(qiáng)大求解器功能，聯(lián)合其他軟件（FlowVision和VirtualLab）進(jìn)行仿真，不僅可以縮短分析時(shí)間，而且可以發(fā)揮兩種軟件的優(yōu)勢(shì)。對(duì)輪胎進(jìn)行流-固耦合分析時(shí)，首先在Abaqus軟件中建立物理場(chǎng)的仿真環(huán)境，流-固耦合作用面上的單元網(wǎng)格和節(jié)點(diǎn)信息從輪胎滾動(dòng)數(shù)值模擬結(jié)果文件odb中提取，將流-固耦合面單元信息提取并寫(xiě)入重啟動(dòng)文件；在流-固耦合重啟動(dòng)文件中，使用Abaqus軟件中的“Co-Simulation”關(guān)鍵詞調(diào)用外部流場(chǎng)仿真軟件，實(shí)現(xiàn)聯(lián)合仿真。在滾動(dòng)輪胎在接地過(guò)程中，花紋溝的泵浦效應(yīng)以及輪胎表面的轉(zhuǎn)和振動(dòng)將導(dǎo)致周?chē)諝猱a(chǎn)生不連續(xù)的壓力波動(dòng)，所發(fā)聲音都屬于氣動(dòng)噪聲[6-7]。

本工作主要研究輪胎花紋對(duì)噪聲的影響，以花紋橫溝寬度為自變量，將流-固耦合仿真后的壓力波動(dòng)結(jié)果導(dǎo)入VirtualLab中進(jìn)行聲學(xué)仿真建模[8-9]。由于在流-固耦合分析時(shí)輪胎花紋與胎體同時(shí)參與分析，聲學(xué)仿真時(shí)應(yīng)同時(shí)考慮輪胎內(nèi)表面及其中的流場(chǎng)，因此需封閉胎側(cè)部分，以防止其對(duì)輪胎外部噪聲產(chǎn)生影響。輪胎結(jié)構(gòu)圖的網(wǎng)格處理如圖4所示，基本場(chǎng)點(diǎn)網(wǎng)格的設(shè)置如圖5所示。

圖4 輪胎結(jié)構(gòu)圖的網(wǎng)格處理

圖5 聲學(xué)場(chǎng)點(diǎn)網(wǎng)格

3 試驗(yàn)結(jié)果分析

3.1 試驗(yàn)方案設(shè)計(jì)

本研究涉及的輪胎試驗(yàn)在半消聲室進(jìn)行，以排除外界環(huán)境因素的干擾。為提供自由場(chǎng)或半自由場(chǎng)的試驗(yàn)測(cè)試環(huán)境[9]，半消聲室背景噪聲設(shè)置為25 dB。將試驗(yàn)輪胎安裝在載重汽車(chē)右前輪上，通過(guò)轉(zhuǎn)鼓驅(qū)動(dòng)輪胎，實(shí)現(xiàn)輪胎定速轉(zhuǎn)動(dòng)。輪胎近場(chǎng)噪聲試驗(yàn)傳感器布置如圖6所示。為測(cè)試不同角度近場(chǎng)噪聲，在輪胎不同方位布置9個(gè)傳感器，同時(shí)采集噪聲數(shù)據(jù)，每個(gè)傳感器距輪胎中心位置均為1 m。

圖6 輪胎近場(chǎng)噪聲試驗(yàn)傳感器布置

3.2 聲學(xué)仿真結(jié)果與試驗(yàn)結(jié)果對(duì)比

輪胎滾動(dòng)仿真是利用轉(zhuǎn)鼓的驅(qū)動(dòng)進(jìn)行，輪胎中心保持不變。為研究仿真分析結(jié)果與試驗(yàn)工況的吻合度，選取1，3，5，7和9共5個(gè)傳感器的數(shù)據(jù)進(jìn)行對(duì)比分析，將仿真結(jié)果與試驗(yàn)結(jié)果取均值，以1號(hào)傳感器為0°方位依次類(lèi)推，將數(shù)據(jù)繪制于極坐標(biāo)下，如圖7所示。可見(jiàn)輪胎近場(chǎng)噪聲試驗(yàn)與仿真結(jié)果都具有明顯的指向性[9]，輪胎接地前后端的聲壓級(jí)最大，且接地后端的聲壓級(jí)大于接地前端，90°方向（5號(hào)）的聲壓級(jí)最小，總體呈“M”字形分布。

圖7 聲壓級(jí)隨測(cè)點(diǎn)方位角度的變化曲線

為了進(jìn)一步研究輪胎噪聲產(chǎn)生的機(jī)理[10]，結(jié)合圖8和9對(duì)90°方向聲場(chǎng)測(cè)點(diǎn)的試驗(yàn)和仿真頻譜特性進(jìn)行對(duì)比分析。結(jié)果表明，試驗(yàn)與仿真結(jié)果在主要的峰值頻率（400，800，1 000和1 500 Hz）上具有較好的吻合度，仿真結(jié)果的峰值頻率主要分布在360，610，1 090，1 420和1 610 Hz，試驗(yàn)結(jié)果峰值頻率主要分布在356，617，1 078，1 424和1 617 Hz。分析比較可知，頻率在300～1 500 Hz之間時(shí)，仿真分析結(jié)果與試驗(yàn)結(jié)果的峰值頻率的誤差范圍為4～12 Hz。

圖8 不同花紋橫溝寬度的輪胎在90°測(cè)點(diǎn)的仿真結(jié)果

圖9 不同花紋橫溝寬度的輪胎在90°測(cè)點(diǎn)的試驗(yàn)結(jié)果

節(jié)距噪聲（泵浦噪聲）的基本頻率（f）計(jì)算如下：

式中，Ve為輪胎滾動(dòng)速度，rt為輪胎滾動(dòng)半徑，np為花紋周節(jié)數(shù)。輪胎的運(yùn)行工況為：Ve＝70 km·h-1，rt＝0.52 m，np＝54，計(jì)算得到f為321 Hz。仿真結(jié)果體現(xiàn)了節(jié)距噪聲的特性，試驗(yàn)結(jié)果也在二階、三階、四階和五階節(jié)距噪聲頻率附近出現(xiàn)峰值，尤其是第5階節(jié)距噪聲頻率與理論結(jié)果吻合較好。

通過(guò)仿真分析結(jié)果與試驗(yàn)聲壓頻譜對(duì)比發(fā)現(xiàn)，花紋橫溝寬度為8和12 mm的輪胎峰值趨勢(shì)較為一致，花紋橫溝寬度為8 mm的輪胎在中高頻段的聲壓級(jí)較為平緩，衰減速度相對(duì)較快；花紋橫溝寬度為12 mm的輪胎聲壓級(jí)較陡，衰減速度較慢。

4 結(jié)論

（1）輪胎近場(chǎng)噪聲主要分布在輪胎接地前端和后端，試驗(yàn)與仿真分析結(jié)果均表明其具有明顯的指向性，在0°和180°方向的聲壓級(jí)最大，在3 000 Hz左右接地后端的聲壓級(jí)明顯大于接地前端，90°測(cè)點(diǎn)的聲壓級(jí)最小，試驗(yàn)與仿真聲壓級(jí)在極坐標(biāo)上呈“M”字形分布。

（2）通過(guò)改變花紋橫溝寬度，發(fā)現(xiàn)輪胎近場(chǎng)噪聲的聲壓級(jí)在中高頻段衰減速度與花紋橫溝寬度呈負(fù)相關(guān)，花紋橫溝寬度為8 mm的輪胎在中高頻段的聲壓級(jí)相對(duì)平緩，而花紋橫溝寬度為12 mm的輪胎的聲壓級(jí)則表現(xiàn)相對(duì)陡勢(shì)。節(jié)距噪聲基本頻率與花紋橫溝寬度無(wú)關(guān)。

（3）通過(guò)對(duì)比90°測(cè)點(diǎn)的噪聲頻譜特性，仿真結(jié)果與試驗(yàn)結(jié)果非常接近，誤差范圍為4～12 Hz。