藺朝莉，彭 勇，陳志軍，呂向飛

（1.重慶電子工程職業(yè)學(xué)院，重慶 400054； 2.重慶長安汽車股份有限公司，重慶 400102；3.重慶大學(xué) 機械工程學(xué)院機械傳動國家重點實驗室，重慶 400044）

懸架液壓襯套是一種新型車用隔振元件，由于其優(yōu)異的性能，正逐步在汽車上推廣應(yīng)用。液壓襯套是典型的非線性彈性體和非線性流體的耦合件。依靠橡膠彈性體的彈性作用和流體的阻尼作用的綜合組合，使得液壓懸置具有傳統(tǒng)橡膠襯套所不具備的良好的隔振和抗沖擊能力的雙重優(yōu)勢。傳統(tǒng)的橡膠襯套往往無法兼顧平順性和操縱性，而液壓襯套可以在保證懸架系統(tǒng)具備較好操縱穩(wěn)定性的前提下，盡可能提高懸架系統(tǒng)對路面沖擊載荷的衰減作用，減小路面的振動沖擊通過懸架連接路徑向車身的傳遞[1-3]。

國內(nèi)外學(xué)者和工程師們對液壓襯套進行了大量的仿真與試驗研究，積累了大量的設(shè)計經(jīng)驗。目前，液壓襯套的研究方法有流固耦合有限元法和集總參數(shù)模型方法。流固耦合有限元法的分析精度較高，但計算效率低。為了獲得液壓襯套的頻變動特性，需要對流固耦合有限元模型在所分析的頻率范圍內(nèi)，每隔一定的頻率區(qū)間分別進行加載計算，并提取響應(yīng)，通過輸入位移激勵和輸出力響應(yīng)之間的關(guān)系推算出其動特性。而集總參數(shù)方法利用力學(xué)等效方法將液壓襯套的復(fù)雜性能參數(shù)歸結(jié)為主簧剛度、等效泵壓面積、上液室體積剛度、液體慣性系數(shù)和液體阻尼系數(shù)等幾個有限的集總參數(shù)，基于力學(xué)模型對液壓襯套的動特性進行仿真[4-7]。集總參數(shù)模型的仿真精度取決于其性能參數(shù)的精確度，常用的集總參數(shù)模型獲取方法有試驗法、經(jīng)驗公式法和智能算法識別法等[8] -[15]，而這些方法大多依賴于樣件試驗，且無法在方案設(shè)計階段就能夠?qū)σ簤阂r套的幾何結(jié)構(gòu)參數(shù)對性能的影響作出準(zhǔn)確定量分析。因此，將流固耦合有限元模型和集總參數(shù)模型結(jié)合的模型方法是當(dāng)前的研究熱點，這種組合模型能夠充分發(fā)揮2種模型的優(yōu)勢，取長補短，在提高分析精度的同時，保證計算效率。

本文在分析液壓襯套結(jié)構(gòu)和工作原理的基礎(chǔ)上，采用AMESim建立液壓襯套的仿真模型，通過流固耦合有限元技術(shù)，提取了集總參數(shù)中的各模型參數(shù)，將參數(shù)代入AMESim模型，建立流固耦合-集總參數(shù)組合模型，并通過試驗驗證了仿真模型的精度。

1 液壓襯套結(jié)構(gòu)和工作原理

圖1所示是液壓襯套的實際結(jié)構(gòu)，液壓襯套主要由橡膠主簧、流道和上下液室組成。

圖1 液壓襯套實際結(jié)構(gòu)

橡膠主簧是主要的承載和受力元件，為液壓襯套系統(tǒng)提供剛度，而液體通過流道在上下液室的來回震蕩運動為液壓襯套系統(tǒng)提供液體阻尼。

圖2所示是由液壓襯套的結(jié)構(gòu)模型抽象簡化而來的力學(xué)模型。

圖2 液壓襯套力學(xué)模型

圖2中Kr、Cr為襯套橡膠主簧的動剛度和阻尼系數(shù)，Ap為襯套的等效泵壓面積，K1和K2為上下液室的體積剛度，其中K1隨著空氣腔的開閉而改變，P1、P2為上下液室內(nèi)的壓強，Q為流道內(nèi)的液體流量，I、R為流道內(nèi)液體的慣性系數(shù)和流量阻尼系數(shù)，X為作用于襯套上端的位移激勵，F(xiàn)為液壓襯套的傳遞力。

根據(jù)流體力學(xué)的原理可得以下公式

路面振動位移激勵經(jīng)襯套傳遞到車身上的力為

對式（1）至式（4）進行拉氏變換，合并可得液壓襯套的復(fù)剛度表達式如下

其中：s=jω為復(fù)變量。

按照定義，液壓襯套的動剛度為復(fù)剛度表達式的實部，液壓襯套的阻尼角為復(fù)剛度表達式的虛部與實部比值的正切值即

2 液壓襯套的AMESim模型

根據(jù)上文對液壓襯套工作原理的分析以及所推導(dǎo)的力學(xué)模型，采用AMESim建立液壓襯套的動特性仿真分析模型。AMESim建模過程一般分為4步：草圖模式下搭建草圖模型；子模型模式下選擇合理的子模型；參數(shù)模式下設(shè)置正確的參數(shù)；運行模式下進行仿真運算以及分析。圖3所示的模型即為液壓襯套的AMESim模型。

圖3 AMESim液壓襯套模型

3 液壓襯套參數(shù)的獲取

3.1 液壓襯套主簧參數(shù)

為了計算橡膠主簧靜剛度，對液壓襯套的橡膠主簧進行靜載荷分析。橡膠材料選用Mooney-Rivlin模型，模型共包含單元32658個，其中橡膠單元20948個，金屬單元8068個，其余為塑料單元。將載荷施加于有限元模型后可獲得橡膠主簧的垂向位移云圖。由仿真的載荷-位移曲線，就可以提取出主簧剛度。

3.2 液固耦合有限元參數(shù)提取

在液壓襯套三維數(shù)模的基礎(chǔ)上，將各部分結(jié)構(gòu)導(dǎo)入Hypermesh中，進行幾何清理，并畫出網(wǎng)格。分別將畫好的固體網(wǎng)格和液體網(wǎng)格模型導(dǎo)入ADINA中，進行流固耦合模型的設(shè)置。圖4所示為把固體模型和液體模型結(jié)合起來的液壓襯套液固耦合模型。

圖4 液固耦合模型

在建立了流固耦合有限元模型的基礎(chǔ)上，進行加載，從而獲取集總參數(shù)模型中各個物理參數(shù)。

3.2.1 上下液室體積剛度提取

上液室體積剛度定義為上液室壓強的變化量與體積變化量之比。仿真計算時，在上液室液體下表面施加一垂向速度v，經(jīng)過一定時間t，可以得到流入上液室的液體體積和上液室的均布壓強P，由此計算上液室體積剛度

根據(jù)仿真獲得的壓力即可通過式（8）的求解獲得上下液室的體積剛度。

3.2.2 等效泵壓面積提取

為提取液壓襯套的等效泵壓面積，在主簧中心位置施加一個向下的力，使之推動液體向下運動，并帶動活塞向下移動。通過計算活塞掃過的體積與主簧的位移之比，便可以獲得液壓襯套的等效泵壓面積Ap。即等效泵壓面積的計算公式為

式中：S為活塞位移，D為活塞下端直徑，X為主簧位移。

3.2.3 液體慣性系數(shù)和阻尼系數(shù)提取

由式（3）知，只要提取流道兩端的壓強差ΔP=P1-P2，并計算流道內(nèi)液體的流量Q及流量的變化率˙，便可以獲得流道內(nèi)液體的慣性系數(shù)I和流量阻尼系數(shù)R。

在時刻i，流道兩端的壓差、流量和流量的變化率分別用 ΔPi、Qi和表示 ，取n個時刻點，由式（3）得

其中

對于式（10），方程中有2個未知數(shù)I和R，引入最小二乘法，通過最小化誤差的平方和獲得方程組的解，即

將上下液室的仿真數(shù)據(jù)代入式（10）至式（12）求解，可得流道慣性系數(shù)和流道阻尼系數(shù)。

4 液壓襯套動特性仿真與試驗

采用液固耦合有限元方法提取得到的液壓襯套物理參數(shù)如表1所示。

表1 基于流固耦合有限元的參數(shù)提取值

將各個參數(shù)代入AMESim模型中便可以得到液壓襯套的動態(tài)特性。通過MTS 831試驗臺對液壓襯套進行掃頻試驗，獲取不同頻率下的動剛度和阻尼角，如圖5所示。

圖5 仿真與試驗對比

由圖5所示的試驗和仿真數(shù)據(jù)的對比可知，流固耦合-集總參數(shù)組合模型可以對液壓襯套的動特性進行較高精度的仿真。動剛度的最大仿真誤差為5.2%，而阻尼角的最大仿真誤差為12%，仿真的誤差均控制在15%以內(nèi)，滿足工程要求。

5 模型參數(shù)影響分析

載荷、材料參數(shù)、網(wǎng)格精度、幾何簡化程度和分析頻率范圍等因素都會對流固耦合所提取出的模型參數(shù)產(chǎn)生影響。在實際建模中，有必要根據(jù)實際測試數(shù)據(jù)對流固耦合有限元模型進行相應(yīng)的修正，提高其分析精度。為分析建模的主要參數(shù)在變載變頻下提取時的非線性變化對建模影響程度，在建立的流固耦合-集總參數(shù)組合模型的基礎(chǔ)上，分析表1中列出5個主要的模型參數(shù)對動剛度和阻尼角曲線的影響，本小節(jié)的分析同時可以為液壓襯套的性能改進和參數(shù)優(yōu)化提供依據(jù)。

圖6所示是等效泵壓面積變化對動剛度和阻尼角的影響規(guī)律圖。

圖6 等效泵壓面積的影響

由圖6中的曲線可知，隨著等效泵壓面積的增加，動剛度和阻尼角都增大。

圖7所示是上液室體積剛度變化對動剛度和阻尼角的影響規(guī)律圖。

由圖7中的曲線可知，隨著上液室體積剛度的增加，動剛度和阻尼角都增大。

圖8所示是橡膠主簧剛度變化對動剛度和阻尼角的影響規(guī)律圖。

由圖8中的曲線可知，隨著橡膠主簧剛度的增加，動剛度增大，而阻尼角減小。這是由于橡膠剛度增大削弱了液體阻尼效應(yīng)的作用。

圖9所示是流道慣性系數(shù)變化對動剛度和阻尼角的影響規(guī)律圖。

由圖9中的曲線可知，隨著流道慣性系數(shù)的增加，動剛度和阻尼角曲線向左移動，說明液柱的共振頻率減小。

圖7 上液室體積剛度的影響

圖8 橡膠主簧剛度的影響

圖10所示是流道阻尼系數(shù)變化對動剛度和阻尼角的影響規(guī)律圖。

圖9 流道慣性系數(shù)的影響

圖10 流道阻尼系數(shù)的影響

由圖10中的曲線可知，隨著流道阻尼系數(shù)的增加，動剛度和阻尼角減小，說明增加流道阻尼系數(shù)削弱了液體的共振能力，減小了動剛度和阻尼角的共振峰幅值大小。

6 結(jié)語

（1）通過集總參數(shù)數(shù)學(xué)模型和AMESim軟件建立了液壓襯套的動特性仿真模型，可以方便地進行動特性仿真和參數(shù)影響分析，為液壓襯套的設(shè)計分析提供了一種實用便捷的技術(shù)手段；

（2）根據(jù)流固耦合有限元模型，在ADINA軟件中提取了液壓襯套集總參數(shù)模型中的物理參數(shù)，在保證集總參數(shù)模型參數(shù)精度的同時，減小了建模計算量，提高了分析效率；

（3）根據(jù)液壓襯套動特性的臺架試驗，對液壓襯套動特性仿真的流固耦合-集總參數(shù)組合模型的仿真精度進行了驗證。

（4）在流固耦合-集總參數(shù)組合模型的基礎(chǔ)上，分析了主要模型參數(shù)對動特性的頻變特性的影響，可以為液壓襯套的性能改進和參數(shù)優(yōu)化提供參考。