湯超群 趙雨雷 吳兵 梁濤

摘? 要：本文首先通過對(duì)減振器測(cè)試所得的頻響函數(shù)進(jìn)行數(shù)據(jù)處理與擬合，得到了減振器垂直方向的等效剛度、阻尼曲線;所得曲線可作為減振器有限元模型的數(shù)據(jù)輸入。其次通過對(duì)減振器進(jìn)行模態(tài)試驗(yàn)、仿真結(jié)果分析對(duì)比，調(diào)整了減振器有限元模型CBUSH單元的徑向剛度，使仿真、試驗(yàn)結(jié)果趨于一致，保證了所建減振器有限元模型的有效性。通過上述試驗(yàn)、仿真分析工作，提供了一種可用于整車路噪仿真分析的減振器有限元建模方法。

關(guān)鍵詞：減振器;頻響函數(shù);有限元;整車路噪

中圖分類號(hào)：U463.3? ? ?文獻(xiàn)標(biāo)識(shí)碼：A? ? 文章編號(hào)：1005-2550（2021）02-0042-05

Abstract： Firstly， through data processing and fitting of the frequency response function obtained from the vibration absorber test， the equivalent stiffness and damping curve of the vertical direction of the shock absorber are obtained. The curve obtained can be used as the data input of finite element model of shock absorber. Secondly， by modal test and simulation analysis， The radial stiffness of CBUSH element of finite element model of shock absorber is adjusted ， So that the simulation and test results are consistent. Then， further assurance the validity of finite element model of shock absorber. Through the above experiment and simulation， the finite element modeling method of shock absorber， which can be used for road noise simulation analysis of the whole vehicle， is standardized.

Key Words： Shock Absorber; Frequency Response Function; Finite Element; Road Noise

引? ? 言

作為用戶抱怨最多，消費(fèi)者購車重要參考指標(biāo)的路面振動(dòng)噪聲問題，越來越得到主機(jī)廠的重視。通過CAE仿真分析技術(shù)，對(duì)整車路噪進(jìn)行仿真分析研究，優(yōu)化整車結(jié)構(gòu)，把問題暴露在產(chǎn)品開發(fā)前期，減少后期問題整改，是整車路噪研究的重要手段[2]。

懸架作為路噪的第二級(jí)隔振系統(tǒng)，包括彈性元件、減震器和襯套。作為整車主要的阻尼部件，減震器的非線性特性尤為明顯，而整車路噪仿真分析模型為線性模型，難以準(zhǔn)確模擬減振器剛度、阻尼特性。為此，在整車路噪仿真分析中，考慮減振器剛度、阻尼因素，準(zhǔn)確建立減震器仿真分析模型尤為重要[3]。

為了規(guī)范用于整車路噪仿真分析的減振器有限元模型建模方法，本文首先通過減振器的頻率響應(yīng)函數(shù)測(cè)試及數(shù)據(jù)擬合結(jié)果獲取了減振器的等效動(dòng)剛度、阻尼參數(shù)。其次，用Hypermesh軟件完成減振器的網(wǎng)格劃分的同時(shí)，基于測(cè)試所得等效動(dòng)剛度、阻尼參數(shù)，用CBUSH單元模擬減振器垂直方向的剛度、阻尼特性。減振器中的橡膠元件的剛度特性也用CBUSH模擬。最后，通過減振器的模態(tài)測(cè)試與仿真結(jié)果分析對(duì)比，修正了減振器橡膠剛度參數(shù)，驗(yàn)證了所建減振器有限元模型的有效性。

1? ? 理論背景

整車在平整路面激勵(lì)下，減振器拉桿處于小滑移運(yùn)動(dòng)狀態(tài)。此時(shí)，減振器可近似為單自由度的質(zhì)量-彈簧-阻尼系統(tǒng)，如下圖1所示。該單自由度系統(tǒng)的移動(dòng)質(zhì)量m由減震器活塞桿、活塞組成，為已知量。簡(jiǎn)化為單自由度系統(tǒng)的減振器等效剛度k、等效阻尼c可通過測(cè)試得到的柔度系數(shù)（） 計(jì)算得到。

由單自由度質(zhì)量-剛度-阻尼系統(tǒng)柔度公式（1）：

2? ? 頻響函數(shù)測(cè)試

頻響函數(shù)測(cè)試中，需保證減振器測(cè)試狀態(tài)與整車安裝狀態(tài)一致。如下圖2，可通過調(diào)整彈性繩拉力，保證減振器測(cè)試狀態(tài)壓縮量與整車安裝狀態(tài)一致。

使用AB膠將減振器底座固定在實(shí)心平臺(tái)上。激振器用彈性繩懸吊起來，由上往下激勵(lì)減振器拉桿頂端。為保證激勵(lì)過程不產(chǎn)生晃動(dòng)，試驗(yàn)中需選用質(zhì)量較大的激振器。將兩相傳感器（可同時(shí)測(cè)量激勵(lì)力與響應(yīng)加速度）安裝在激振桿頂端。調(diào)整激振桿位置或激振器高度，使激振桿與拉桿接觸但無預(yù)載。用AB膠將傳感器與拉桿固定。為確保全頻范圍內(nèi)良好激勵(lì)，同時(shí)去除激振桿與儲(chǔ)液缸之間的靜摩擦，獲得良好的激振效果，需選擇合適的激勵(lì)信號(hào)。這里選用2Hz正弦信號(hào)與300Hz隨機(jī)信號(hào)的組合信號(hào)作為激勵(lì)信號(hào)。正弦信號(hào)幅值為100N，隨機(jī)信號(hào)幅值為50N。分析帶寬設(shè)為512Hz，頻率分辨率設(shè)為1Hz，采集模式選擇Free run。測(cè)試得到的減振器拉桿頂端頻響函數(shù)如下圖3所示：

將計(jì)算得到的柔度函數(shù)用實(shí)部（Real）、虛部（Imag）的形式表達(dá)，并導(dǎo)出到Excel中。按公式（5）編輯公式，計(jì)算得到減振器垂直方向的等效剛度、等效阻尼曲線，分別如下圖4、5所示：

進(jìn)一步對(duì)計(jì)算得到的剛度、阻尼曲線進(jìn)行線性擬合。根據(jù)曲線走勢(shì)，將曲線分為0-120Hz、120Hz-250Hz兩段進(jìn)行線性擬合，擬合所得的線性剛度、阻尼如公式（6）、公式（7）：

3? ? 減振器有限元模型的建立

對(duì)導(dǎo)出的減振器3D數(shù)模進(jìn)行簡(jiǎn)化處理后，所建減振器有限元模型如圖6所示。模型中用RB2單元模擬螺栓、焊縫連接。減振器中的橡膠元件用CBUSH單元模擬。減振器拉桿用CBAR單元模擬。因簡(jiǎn)化了減振器中的橡膠，減振液、螺母、焊縫等，需用CONM2單元對(duì)模型進(jìn)行配重。

如圖6，減振器上部連接位置，即圖中1號(hào)位置建立有CBar單元、1個(gè)CBUSH單元、1個(gè)集中質(zhì)量CONM2單元和RBE2連接單元。減振器中部連接位置，即圖中的2號(hào)位置，建立有CBar單元、1個(gè)CBUSH單元、1個(gè)集中質(zhì)量CONM2單元和RBE2連接單元。減振器活塞底部連接位置，即圖中的3號(hào)位置，建立有CBar單元、2個(gè)CBUSH單元、1個(gè)集中質(zhì)量CONM2單元和RBE2連接單元。其中1#位置處的CBUSH單元用來對(duì)該處的橡膠特性進(jìn)行模擬;2#位置處的CBUSH單元用來模擬減振器垂直方向的等效剛度、阻尼特性。減振器的底部連接位置，即圖中的4號(hào)位置，建立有1個(gè)集中質(zhì)量單元和RBE2連接單元。因需保持所建減振器模型與減振器試驗(yàn)狀態(tài)一致，對(duì)4號(hào)位置RBE2單元中心節(jié)點(diǎn)添加全局約束。模型中的CBUSH單元參考坐標(biāo)系以減振器儲(chǔ)液筒軸線方向?yàn)閦軸，減振器儲(chǔ)液筒半徑方向?yàn)閤、y方向。

按減振器實(shí)際運(yùn)動(dòng)特性及內(nèi)部橡膠作用，賦予模型中CBUSH單元?jiǎng)偠?、阻尼特性。即，只考慮上圖中1號(hào)、2號(hào)位置CBUSH，以及3號(hào)位置1#CBUSH單元的徑向剛度特性，以及3號(hào)位置2#CBUSH單元的垂向剛度、阻尼特性。按經(jīng)驗(yàn)賦給1號(hào)、2號(hào)位置CBUSH單元，以及3號(hào)位置1#CBUSH單元12000N/mm的徑向剛度初始值。3號(hào)位置2#CBUSH單元的垂向剛度、阻尼依據(jù)公式（6）、（7）擬合所得的減振器等效線性剛度、阻尼曲線賦予，如下圖7所示。圖中初始剛度、阻尼值為50Hz時(shí)線性等效剛度、阻尼曲線所對(duì)應(yīng)的值。

4? ? 試驗(yàn)?zāi)B(tài)與仿真模態(tài)對(duì)比

為驗(yàn)證所建減振器有限元模型的有效性，以及進(jìn)一步修正模型參數(shù)，分別對(duì)所建減振器有限元模型進(jìn)行約束模態(tài)仿真分析和相同狀態(tài)下的減振器模態(tài)測(cè)試。

減振器約束模態(tài)試驗(yàn)結(jié)果如下圖8所示：

減振器約束模態(tài)仿真結(jié)果如下圖9所示：

對(duì)比減振器模態(tài)試驗(yàn)、仿真結(jié)果發(fā)現(xiàn)：第一階振型均為X向的整體彎曲模態(tài);仿真頻率為122.6Hz，試驗(yàn)頻率為118.9Hz。第二階振型均為Y向整體彎曲模態(tài);仿真結(jié)果為137.2Hz，試驗(yàn)結(jié)果為131.0Hz。第三階振型，試驗(yàn)為活塞桿的X向彎曲模態(tài)，仿真為活塞桿的Y向彎曲模態(tài);試驗(yàn)頻率為306.0Hz，仿真頻率為316.0Hz。第四階振型，試驗(yàn)為活塞桿的Y向彎曲模態(tài)，仿真為活塞桿的X向彎曲模態(tài);試驗(yàn)頻率值327.0Hz，仿真頻率值為324.5Hz。減振器模態(tài)試驗(yàn)、仿真結(jié)果存在一定的差異。

因減振器仿真各階模態(tài)頻率主要由1號(hào)、2號(hào)位置CBUSH，以及3號(hào)位置1#CBUSH單元的徑向剛度決定。通過調(diào)整上述各位置處CBUSH單元的徑向剛度，可使仿真與試驗(yàn)?zāi)B(tài)頻率結(jié)果一致。調(diào)整后的各CBUSH單元徑向剛度數(shù)值如下表1所示：

調(diào)整后，減振器模態(tài)振型及對(duì)應(yīng)頻率仿真、試驗(yàn)結(jié)果對(duì)比如下表2所示：

由上表2可知，調(diào)整CBUSH單元?jiǎng)偠群蟮臏p振器模態(tài)仿真分析結(jié)果與試驗(yàn)一致，各振型對(duì)應(yīng)模態(tài)頻率誤差不超過5%。所建減振器有限元模型可裝配到整車路噪模型中用于整車路噪仿真分析。

5? ? 結(jié)論

（1）本文首先對(duì)減振器進(jìn)行了頻響函數(shù)測(cè)試。通過對(duì)測(cè)試所得數(shù)據(jù)的處理與擬合，得到了減振器線性等效剛度、阻尼曲線。所得曲線可作為減振器有限元模型的數(shù)據(jù)輸入。

（2）其次對(duì)減振器進(jìn)行了模態(tài)試驗(yàn)、仿真對(duì)比分析。通過調(diào)整減振器有限元模型CBUSH單元徑向剛度，使模態(tài)仿真、試驗(yàn)結(jié)果一致，保證了所建減振器有限元模型的有效性。

（3）本文通過上述試驗(yàn)、仿真分析工作，提出了一種可用于整車路噪仿真分析的減振器有限元模型建模方法。

參考文獻(xiàn)：

[1]龐劍，諶剛，何華. 汽車噪聲與振動(dòng)-理論與應(yīng)用[M]. 北京：北京理工大學(xué)出版社，2006.

[2]陳釗，霍俊焱，吳星晨，鄧江華，王海洋.整車路噪仿真及優(yōu)化技術(shù)應(yīng)用[J].噪聲與振動(dòng)控制，2019，39（03）：147-151.

[3]沈陽，黃元毅，梁靜強(qiáng)，常光寶，呂兆平.路面譜激勵(lì)整車噪聲分析流程開發(fā)研究[J].重慶交通大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，38（07）：119-125.

[4]劉德柱，趙唐雷，劉艷華，王海.基于AMESim某車型減振器仿真方法研究[J].汽車實(shí)用技術(shù)，2019（10）：127-130.

[5]劉曉東，李俊杰，曹逢源，姜洪雷.硅橡膠減振器的有限元建模與動(dòng)態(tài)特性分析[J].噪聲與振動(dòng)控制，2018，38（06）：194-198.