王巽，張紅業(yè)，楊權(quán)，何家興

廣州汽車集團(tuán)股份有限公司汽車工程研究院，廣東廣州 511434

0 引言

近年來，隨著汽車電動(dòng)化進(jìn)程加快，汽車的噪聲、振動(dòng)、聲振粗糙度（noise、vibration、harshness，NVH）性能作為越來越重要的評(píng)價(jià)指標(biāo)影響消費(fèi)者的購(gòu)買選擇，各個(gè)主機(jī)廠也均在該領(lǐng)域投入了越來越多的資源進(jìn)行優(yōu)化改善。

來自地面的隨機(jī)激勵(lì)通過輪胎系統(tǒng)的初步衰減傳遞至懸架系統(tǒng)，作為噪聲、振動(dòng)傳遞路徑中剛性乘員艙接受振動(dòng)前的最后一環(huán)，懸架系統(tǒng)的隔振表現(xiàn)關(guān)乎整車的NVH性能，而襯套作為懸架系統(tǒng)內(nèi)控制臂與軸節(jié)、控制臂與車架及車架與車身連接的多自由度柔性體單元，需在較寬的頻率段提供有效的降噪隔振功能。扭梁縱臂前點(diǎn)襯套作為懸架襯套的一種，其動(dòng)態(tài)特性會(huì)顯著影響整車NVH性能表現(xiàn)[1]。

目前，廣泛應(yīng)用的傳統(tǒng)空實(shí)心橡膠襯套具有結(jié)構(gòu)簡(jiǎn)單、性能可靠等優(yōu)點(diǎn)，但限于結(jié)構(gòu)及材料特性，提供的最大阻尼角較?。?0°左右），且基礎(chǔ)動(dòng)剛度較大，高頻段容易產(chǎn)生動(dòng)態(tài)硬化，不能有效兼顧各頻率段隔振降噪等要求。在空實(shí)心橡膠襯套的基礎(chǔ)架構(gòu)上封裝液體形成液壓型橡膠襯套，可有效改善上述問題。一個(gè)典型的液壓襯套由外管、內(nèi)管、橡膠主簧、封裝液體、液體流通通道、兩個(gè)儲(chǔ)液室及橡膠主簧內(nèi)的剛性骨架組成。其中，橡膠主簧的主要作用為支撐連接系統(tǒng)的靜載荷，并提供一定的基礎(chǔ)阻尼作用；橡膠主簧內(nèi)部有剛性骨架，其主要作用為在襯套受外力時(shí)保持液體流動(dòng)通道的整體輪廓。骨架的結(jié)構(gòu)設(shè)計(jì)對(duì)整個(gè)襯套的耐久性及剛度調(diào)校有重要影響。液體流動(dòng)通道可分為慣性通道式和孔口式，慣性通道式具有較長(zhǎng)的流動(dòng)通道，橫截面一般為矩形；孔口式流動(dòng)通道又分為節(jié)流孔式和液柱共振式，節(jié)流孔式通道孔徑較小，阻尼峰值較低，作用區(qū)間主要位于低頻率段，液柱共振式孔徑較大，阻尼作用頻率區(qū)間較寬。當(dāng)襯套受到外部的動(dòng)態(tài)載荷時(shí)，橡膠主簧產(chǎn)生變形，儲(chǔ)液室的內(nèi)部液體在壓力的作用下在流動(dòng)通道內(nèi)來回流動(dòng)，通過液體在該過程中的摩擦損失、慣性損失及局部損失來衰減振動(dòng)能量[2]。

目前，業(yè)內(nèi)對(duì)液壓型橡膠襯套的動(dòng)特性研究有3種方法：集總參數(shù)法、試驗(yàn)研究法、有限元仿真分析法。其中，集總參數(shù)法通過建立液壓型襯套較準(zhǔn)確的參數(shù)模型，從理論概念的角度對(duì)襯套的動(dòng)態(tài)性能表現(xiàn)進(jìn)行分析，其分析難度在于集總參數(shù)模型的建立及關(guān)鍵參數(shù)的精確識(shí)別；試驗(yàn)研究法是通過試驗(yàn)實(shí)測(cè)的方法對(duì)液壓型襯套進(jìn)行解析研究，試驗(yàn)結(jié)果對(duì)參數(shù)模型的建立及設(shè)計(jì)方案的驗(yàn)證有一定指導(dǎo)意義，但做件-測(cè)試迭代過程工序煩瑣且較為耗時(shí)，試驗(yàn)過程需嚴(yán)格控制其余變量，對(duì)試驗(yàn)設(shè)計(jì)要求較高；有限元仿真分析法是基于流體力學(xué)（computational fluid dynamics, CFD）和流-固耦合分析（fluid structure interaction, FSI）等方法對(duì)液體在通道內(nèi)的流動(dòng)過程進(jìn)行模擬計(jì)算，經(jīng)過標(biāo)定過的有限元參數(shù)模型可以預(yù)測(cè)優(yōu)化后產(chǎn)品的性能表現(xiàn)，一定程度上縮短了產(chǎn)品的正向開發(fā)周期。該方法中有限元模型的搭建和流-固耦合邊界的處理為首要難點(diǎn)[3-4]。

1 液壓型襯套的集總參數(shù)模型

以某車型扭力梁?jiǎn)瓮ǖ酪簤盒鸵r套為例，建立集總參數(shù)模型（圖1）。

圖1 單慣性通道液壓型扭力梁襯套集總參數(shù)模型

圖1中，Br、Kr分別代表橡膠主簧的阻尼系數(shù)，N·s/m及動(dòng)剛度，N/m。Ap1與Ap2分別代表兩個(gè)液室內(nèi)的等效活塞面積，m2，本例中Ap1=Ap2=Ap；P1（t）與P2（t）分別代表兩液室內(nèi)的壓力值，N/m2。Ki1與Ki2代表兩個(gè)液室的體積剛度，N/m5。C1和C2為兩個(gè)液室的體積柔度，代表兩個(gè)液室受到壓力下的膨脹特性。Qi（t）為液體在流通通道內(nèi)的動(dòng)態(tài)流量，m3/s。li為流體通道的長(zhǎng)度值，m。Ai為流體通道的橫截面積，m2。液柱在通道內(nèi)的相對(duì)位移用Xi（t）表示，m。襯套的支反力為FT（t）。Ii為流體通道內(nèi)液體的慣性系數(shù)，kg/m4，其解析式為Ii=mi/Ai2=ρli/Ai，其中mi為流通通道內(nèi)液柱的質(zhì)量；ρ為流體通道內(nèi)液體的密度。R1為流通通道內(nèi)液體流動(dòng)的線性阻尼系數(shù)，R2為流通通道內(nèi)液體流動(dòng)的非線性阻尼系數(shù)，其解析式分別為R1=ε1/Ai2、R2=ε2/Ai5，其中ε1、ε2為與流動(dòng)通道內(nèi)的液體流動(dòng)速度正相關(guān)的阻尼系數(shù)。令KV=Ki1+Ki2，BV=B1+B2，其中，B1、B2為兩個(gè)儲(chǔ)液室的體積阻尼，N·s/m5。鑒于本文主要考察液壓襯套在定幅值下動(dòng)剛度及阻尼角的性能表現(xiàn)，為降低模型復(fù)雜程度，可忽略集總參數(shù)表達(dá)式中慣性通道內(nèi)流體的非線性阻尼，即R2=0。

建立單慣性通道液壓型襯套集總參數(shù)模型，得到動(dòng)剛度的計(jì)算解析式為[5-6]：

阻尼角的計(jì)算解析式為：

式中，K1為儲(chǔ)能剛度，其表達(dá)式為：

K2為耗能剛度，其表達(dá)式為：

經(jīng)數(shù)據(jù)擬合，本文中各參數(shù)的取值為Kr=2.2×105N/m，Br=220 N·s/m，AP=2.3×10-3m2，Ki1、Ki2=1.5×1010N/m5，B1、B2=1.55×105N·s/m5，Ii=2.6×106N·s2/m5，R1=9.6×106N·s/m5，將上述參數(shù)輸入模型中進(jìn)行計(jì)算，解析值與振幅0.05 mm工況下襯套的實(shí)測(cè)數(shù)據(jù)對(duì)比，如圖2所示。

圖2 單通道扭力梁液壓襯套動(dòng)態(tài)特性（振幅0.05 mm）

圖3 具有n個(gè)相同流體通道液壓襯套的集總參數(shù)模型

由圖2可知，動(dòng)剛度、阻尼角的解析值及實(shí)測(cè)值變化趨勢(shì)與峰值頻率基本一致，在0～150 Hz頻率段的相對(duì)誤差均在可接受范圍內(nèi)，證明了模型的可靠性及參數(shù)的準(zhǔn)確性。

假設(shè)襯套兩個(gè)液室之間具有n個(gè)尺寸相同的流體通道，搭建集總參數(shù)模型，得到該型襯套的動(dòng)剛度及阻尼角的解析表達(dá)式[7]。

式中，τX為無量綱的修正系數(shù)，與外部激勵(lì)振幅相關(guān)，本文中τX=0.89。

該車型在樣車驗(yàn)證階段存在路噪大、方向盤振動(dòng)明顯等問題，通過對(duì)噪聲、振動(dòng)傳遞路徑的分析，確認(rèn)扭梁液壓襯套為關(guān)鍵影響因素，經(jīng)與同級(jí)別標(biāo)桿車的零件性能對(duì)比確認(rèn)，扭梁液壓襯套應(yīng)滿足以下關(guān)鍵動(dòng)態(tài)性能指標(biāo)：在預(yù)載0 N、振幅0.05 mm、40 Hz測(cè)試工況下，動(dòng)剛度≤1600 N/mm，阻尼角≥42°?，F(xiàn)有液壓型襯套的兩個(gè)考核值為1810 N/mm、40°，不滿足以上要求，需進(jìn)行性能優(yōu)化。

2 扭梁液壓襯套動(dòng)態(tài)特性正交試驗(yàn)優(yōu)化

2.1 正交試驗(yàn)

由襯套集總參數(shù)表達(dá)式挑選出影響襯套動(dòng)態(tài)性能的6項(xiàng)主要因素：液室的等效活塞面積（Ap）、體積剛度（KV）、流體通道長(zhǎng)度（Ii）、流量阻尼系數(shù)（R1）、流道液體質(zhì)量的慣性系數(shù)（Ii）及流道個(gè)數(shù)（n）等。本文以現(xiàn)有襯套結(jié)構(gòu)為基準(zhǔn)設(shè)置3個(gè)因素水平：等效活塞面積Ap∈[0.8Ap，1.2Ap]，體積剛度KV∈[0.8KV，1.2KV]，流體通道長(zhǎng)度li∈[0.8li，1.2li]，流量阻尼系數(shù)R1∈[0.8R1，1.2R1]，流道液體質(zhì)量的慣性系數(shù)Ii∈[0.8Ii，1.2Ii]，流道個(gè)數(shù)n∈[1，3]，假設(shè)這6個(gè)因素之間不存在明顯的交互影響，設(shè)計(jì)6因素3水平的正交試驗(yàn)，見表1。

表1 因素水平

選擇6因素3水平正交表設(shè)計(jì)方法，設(shè)計(jì)出27組試驗(yàn)，采用上文中搭建的集總參數(shù)模型計(jì)算襯套在40 Hz處的動(dòng)剛度及40 Hz處的阻尼角的解析值，正交試驗(yàn)方案及結(jié)果見表2。

表2 正交試驗(yàn)方案及結(jié)果

參考正交試驗(yàn)結(jié)果，首先根據(jù)單一指標(biāo)的分析方法分別對(duì)兩組優(yōu)化目標(biāo)值進(jìn)行極差分析，優(yōu)選出二者各自的最優(yōu)組合，再對(duì)兩組優(yōu)選結(jié)果進(jìn)行綜合分析，得到最終的因素水平優(yōu)化組合。

2.2 試驗(yàn)數(shù)據(jù)分析

根據(jù)表2中各組合方案的解析結(jié)果，計(jì)算得到各個(gè)因素的極差，結(jié)果見表3，表中ki代表該因素采用第i水平時(shí)試驗(yàn)優(yōu)化指標(biāo)的平均值，極差（R）代表各因素中每一列ki的最大值與最小值的差值。根據(jù)統(tǒng)計(jì)學(xué)理論，極差的大小反映了該因素對(duì)目標(biāo)結(jié)果的影響程度，極差值越大，表示該因素對(duì)試驗(yàn)?zāi)繕?biāo)結(jié)果的影響越顯著。根據(jù)分析的極差結(jié)果繪制出各因素對(duì)目標(biāo)結(jié)果的影響趨勢(shì)，如圖4所示。

表3 極差分析結(jié)果

圖4 各因素對(duì)40 Hz處動(dòng)剛度和阻尼角的影響趨勢(shì)

對(duì)兩個(gè)優(yōu)化指標(biāo)分別進(jìn)行極差計(jì)算發(fā)現(xiàn)，對(duì)于40 Hz處動(dòng)剛度，影響因素的主次順序均為F＞A＞B=E＞D＞C，即流道個(gè)數(shù)＞等效活塞面積＞體積剛度=流道液體質(zhì)量的慣性系數(shù)＞流量阻尼系數(shù)＞流體通道長(zhǎng)度；對(duì)于40Hz處的阻尼角，影響因素的主次順序均為F＞B=D＞C＞A＞E，即流道個(gè)數(shù)＞體積剛度=流量阻尼系數(shù)＞流體通道長(zhǎng)度＞等效活塞面積＞流道液體質(zhì)量的慣性系數(shù)。

由圖4可知，以40 Hz處動(dòng)剛度為主要優(yōu)化指標(biāo)時(shí)，最優(yōu)的結(jié)構(gòu)參數(shù)組合為A2B3C2（C3）D1E1F2；以40 Hz處阻尼角為主要優(yōu)化指標(biāo)時(shí)，最優(yōu)的結(jié)構(gòu)參數(shù)組合為A3B2C2D3E1F1。

由上述結(jié)果可知，兩個(gè)最優(yōu)結(jié)構(gòu)參數(shù)組合并不完全相同，下一步需對(duì)試驗(yàn)的相關(guān)數(shù)據(jù)進(jìn)行方差分析，以確定試驗(yàn)結(jié)果存在波動(dòng)的主要原因，方差分析結(jié)果見表4。

表4 方差分析結(jié)果

由表中各因素對(duì)考察指標(biāo)影響的顯著性可以得出，因素A對(duì)40 Hz處動(dòng)剛度有顯著性影響，故選取A2為優(yōu)解；因素B對(duì)40 Hz處動(dòng)剛度和阻尼角均有顯著性影響，綜合考慮F比后，選取B3為較優(yōu)解；因素C對(duì)40 Hz處阻尼角有顯著性影響，選取C2為優(yōu)解；因素D對(duì)40 Hz處阻尼角有顯著性影響，選取D3為優(yōu)解；因素E對(duì)40 Hz處動(dòng)剛度有顯著性影響，選取E1為優(yōu)解；因素F對(duì)40 Hz處動(dòng)剛度和阻尼角均有顯著性影響，考慮到降低動(dòng)剛度對(duì)整車平順性的提升更加明顯，故選取F2為較優(yōu)解。綜合以上，得出最優(yōu)的結(jié)構(gòu)參數(shù)組合為A2B3C2D3E1F2。

3 優(yōu)化方案驗(yàn)證

在現(xiàn)有襯套結(jié)構(gòu)參數(shù)組合的基礎(chǔ)上對(duì)上述正交試驗(yàn)所得最優(yōu)結(jié)構(gòu)參數(shù)組合進(jìn)行樣件試制，優(yōu)化前后的液體流道模型如圖5所示。

圖5 優(yōu)化前后的液體流道模型

測(cè)試優(yōu)化后襯套在40 Hz處動(dòng)剛度及阻尼角，相關(guān)結(jié)果見表5。由表5可以看出，優(yōu)化后結(jié)構(gòu)40 Hz處動(dòng)剛度由1810 N/mm降低至1580 N/mm；40 Hz處的阻尼角由40°提高至42°，采用優(yōu)化后的襯套結(jié)構(gòu)能滿足既定的動(dòng)態(tài)性能目標(biāo)要求。

表5 優(yōu)化前后襯套的動(dòng)態(tài)性能表現(xiàn)

優(yōu)化結(jié)構(gòu)參數(shù)后，將襯套搭載在實(shí)車上，測(cè)試特征路面（壞路、粗糙路）車內(nèi)路噪及方向盤振動(dòng)的客觀表現(xiàn)，測(cè)試結(jié)果與原襯套客觀數(shù)據(jù)的對(duì)比如圖6所示。由圖6可得，優(yōu)化后襯套能顯著降低車內(nèi)路噪及方向盤振動(dòng)，證明了扭梁液壓襯套的結(jié)構(gòu)優(yōu)化對(duì)整車NVH性能表現(xiàn)有較大的正向作用。

圖6 實(shí)車搭載優(yōu)化前后扭梁襯套在特征路面的客觀測(cè)試數(shù)據(jù)

4 結(jié)論

（1）以40 Hz處動(dòng)剛度為指標(biāo)，得出最優(yōu)結(jié)構(gòu)因素組合為：A2B3C2（C3）D1E1F2，各因素對(duì)40 Hz處動(dòng)剛度的影響程度為：流道個(gè)數(shù)＞等效活塞面積＞體積剛度=流道液體質(zhì)量的慣性系數(shù)＞流量阻尼系數(shù)＞流體通道長(zhǎng)度；以40 Hz處的阻尼角為指標(biāo)，得出最優(yōu)結(jié)構(gòu)因素組合為：A3B2C2D3E1F1，各因素對(duì)40 Hz處的阻尼角的影響程度為：流道個(gè)數(shù)＞體積剛度=流量阻尼系數(shù)＞流體通道長(zhǎng)度＞等效活塞面積＞流道液體質(zhì)量的慣性系數(shù)。

（2）綜合考慮各因素對(duì)兩個(gè)考核指標(biāo)的影響，因降低動(dòng)剛度對(duì)整車平順性的提升更加明顯，最終確定最優(yōu)結(jié)構(gòu)參數(shù)組合為A2B3C2D3E1F2。

（3）根據(jù)最優(yōu)結(jié)構(gòu)參數(shù)組合進(jìn)行液壓襯套改善樣件試制并進(jìn)行動(dòng)態(tài)特性測(cè)試，結(jié)果表明，優(yōu)化后結(jié)構(gòu)40 Hz處動(dòng)剛度由1810 N/mm降低至1580 N/mm；40 Hz處的阻尼角由40°提高至42°，采用優(yōu)化后的襯套結(jié)構(gòu)能滿足既定的動(dòng)態(tài)性能目標(biāo)要求。搭載實(shí)車測(cè)試特征路面客觀數(shù)據(jù)可得，優(yōu)化后襯套能顯著降低車內(nèi)路噪及方向盤振動(dòng)，證明了扭梁液壓襯套的結(jié)構(gòu)優(yōu)化對(duì)整車NVH性能表現(xiàn)有較大的正向作用。