扭力梁懸架有限元模態(tài)分析中的橡膠襯套動剛度賦值方法*

石月奎　李洪亮　王海洋（中國汽車技術研究中心，天津 300300）

?

扭力梁懸架有限元模態(tài)分析中的橡膠襯套動剛度賦值方法*

石月奎李洪亮王海洋
（中國汽車技術研究中心，天津 300300）

【摘要】依托典型車型扭力梁懸架的有限元模態(tài)分析結果，在對懸架邊界條件合理簡化的前提下，針對懸架橡膠襯套單元模型在分析模型構建中的動剛度賦值問題開展研究，提出了一種“分頻迭代”的處理方法。與試驗模態(tài)分析結果對比，采用“分頻迭代”方法進行處理，可將模態(tài)頻率的計算誤差控制在10%以內(nèi)，且相應的模態(tài)振型一致性良好，頻率、振型的一致性均優(yōu)于傳統(tǒng)的靜剛度當量方法，表明了所提方法的有效性與合理性。

主題詞：橡膠襯套動剛度特性扭力梁懸架有限元分析

1　引言

扭力梁懸架系統(tǒng)動態(tài)特性與整車NVH性能密切相關，可通過計算或試驗分析手段揭示其模態(tài)分布規(guī)律［1，2］，扭力梁懸架的有限元模態(tài)分析面臨的突出問題是橡膠襯套的建模（簡化）原則并不十分清晰。具體而言，往往是掌握著詳實的橡膠襯套動剛度參數(shù)，卻無法將其合理賦值于有限元模型中加以有效應用，并由此對有限元分析精度帶來不利影響。

為解決上述問題，依托具體車型扭力梁懸架的有限元模態(tài)分析，在對懸架安裝約束邊界條件合理簡化的前提下，結合試驗與計算研究，提出一種用于懸架橡膠襯套動剛度合理賦值的“分頻迭代”方法。該方法以對襯套單元模型的靜態(tài)剛度賦值及對懸架模態(tài)分布的試算為起點和依據(jù)，對懸架模態(tài)頻率的取值區(qū)間進行預估，進而將預估區(qū)間上襯套動剛度的“有效值”作為賦值依據(jù)，并通過對預估區(qū)間的反復修正確保賦值的合理性。

2　　扭力梁懸架的試驗模態(tài)分析

為了提供扭力梁懸架計算模態(tài)分析的驗證依據(jù)，首先實施了相應的專項測試，包括自由模態(tài)測試與約束模態(tài)測試。

2.1自由模態(tài)測試

以橡膠繩懸掛扭力梁本體，使其處于自由懸掛狀態(tài)，如圖1所示，測試采用力錘激勵，單輸入、多輸出的方式進行。

圖2是試驗模態(tài)分析的結果。扭力梁本體自由模態(tài)分布：1階振型為整體X向彎曲，頻率為152.00 Hz；2階振型為整體Z向彎曲，頻率為192.00 Hz（此處坐標系參照整車坐標系定義，下同）。

2.2約束模態(tài)測試

排除輪胎這一復雜彈性系統(tǒng)的影響以簡化建模和分析，樣車裝配狀態(tài)下保留懸架彈簧和減振器，僅去除輪胎系統(tǒng)，懸架系統(tǒng)與地面間通過方木塊及千斤頂（左、右對稱）建立連接關系，通過千斤頂調(diào)節(jié)垂向高度以保證輪心位置與原狀態(tài)（未去除輪胎時）一致。試驗狀態(tài)見圖3。

試驗模態(tài)分析結果如圖4所示，扭力梁懸架約束模態(tài)分布：1階為斜向彎曲，頻率為85.73 Hz；2階為橫梁繞Y軸橫擺，頻率為116.56 Hz；3階為2階斜向彎曲，頻率為172.00 Hz。

3　　扭力梁懸架有限元模型的建立

3.1扭力梁本體有限元模型與自由模態(tài)分析

懸架橡膠襯套動剛度賦值方法研究的前提是保證其他部件計算模型的精度，故先建立精度可靠的扭力梁本體有限元模型。借助Hypermesh軟件建立扭力梁本體的有限元模型，其中網(wǎng)格單元大小為8 mm，網(wǎng)格類型為四邊形或三角形殼單元，焊縫以RBE2單元模擬，所得模型如圖5所示。

利用MSC.Nastran軟件即可求得扭力梁本體有限元模型的自由模態(tài)分布。將其與自由模態(tài)測試結果對比，如表1所示，扭力梁模型與實測值的質(zhì)量誤差絕對值為0.75%，自由模態(tài)計算值與試驗值誤差絕對值在6%以內(nèi)。振型分布見圖6、圖7，計算振型與試驗振型吻合。這表明扭力梁本體的有限元模型正確且精度較好，可直接用于后續(xù)對橡膠襯套動剛度賦值方法的研究。

表1　　扭力梁本體試驗值與計算值結果對比

3.2扭力梁懸架邊界條件的簡化處理

考慮到扭力梁懸架與車身均通過橡膠襯套或襯墊連接，且橡膠襯套或襯墊的剛度遠小于車身安裝點剛度，其質(zhì)量也遠小于扭力梁的質(zhì)量，對扭力梁懸架邊界條件作如下簡化假設：

a.懸架橡膠襯套基于特性建模，忽略其質(zhì)量；

b.車身安裝點作為剛體處理。

基于簡化假設，在有限元軟件中：扭力梁前端橡膠襯套以CBUSH單元模擬，CBUSH單元可實現(xiàn)襯套6向剛度特性（3個線剛度與3個角剛度）的賦值計算；因?qū)嵦卣鞣治鰰r，不考慮結構的阻尼［3］，將減振器簡化為圓柱副，以RJIONT單元模擬，其上、下端橡膠襯套以CBUSH單元模擬；彈簧以CELAS單元模擬。

4　橡膠襯套動剛度賦值方法及驗證

通過設計專項試驗獲取襯套剛度特性［4，5］。襯套剛度特性包含3個線剛度與3個角剛度，對于非中心對稱橡膠襯套，其方向定義如圖8所示。試驗包括靜剛度測試與動剛度測試，其中動剛度特性測試需保證測試裝夾時襯套預緊狀態(tài)與實車裝配狀態(tài)一致，以使測試所得動剛度值與模態(tài)測試時襯套所表現(xiàn)出的動態(tài)剛度特性吻合。測試過程如圖9所示。所得橡膠襯套剛度特性如圖10、圖11所示（以扭力梁本體前端與車身連接處橡膠襯套的3向線剛度為例）。

根據(jù)以上分析與試驗數(shù)據(jù)，定義橡膠襯套賦值方法。

a.方法1

有限元CBUSH單元中，以襯套靜剛度賦值。其取值為襯套靜剛度特性測試結果中（如圖10示例）線性段斜率，當量值見表2，其中X向、Y向、Z向為線剛度，RX向、RY向、RZ向為角剛度（以扭力梁本體前端與車身連接處橡膠襯套為例，下同）。

b.方法2

有限元CBUSH單元中，以襯套動剛度賦值。因模態(tài)分析中（以錘擊法試驗模態(tài)分析為例），力錘對測試對象的激勵是動態(tài)激勵，由此引起的懸架響應為動態(tài)響應，而懸架的響應又作為懸架橡膠襯套的激勵，使得橡膠襯套表現(xiàn)出動態(tài)特性，則橡膠襯套在某模態(tài)頻率點處被激發(fā)的動剛度特性可能與此階模態(tài)頻率相關。為驗證此分析是否合理，在已知試驗模態(tài)分布的情況下，將橡膠襯套動剛度特性曲線上試驗模態(tài)頻率對應的襯套動剛度賦值到CBUSH單元，其中模態(tài)頻率四舍五入取整數(shù)值，如表2所示。襯套動剛度特性明顯比靜剛度特性“強”，且襯套動剛度特性隨激勵頻率增加而增“強”。

表2懸架橡膠襯套剛度特性當量值　N∕mm

c.方法3

有限元CBUSH單元中，以襯套動剛度賦值。方法2中動剛度賦值的前提是已知試驗模態(tài)分布，若上述對于襯套動剛度取值與模態(tài)頻率內(nèi)在關系的分析合理，則對于未知試驗模態(tài)分布的情況，可先借助橡膠襯套靜剛度值（方法1）試算懸架的約束模態(tài)分布，再基于此預估模態(tài)頻率區(qū)間襯套動剛度當量的“有效值”，同時記錄下與之相對應的頻率值，然后將上述動剛度有效值賦予計算模型，再次求解出懸架的模態(tài)頻率，并與所記錄動剛度有效值的對應頻率相比較，如二者充分接近，則處理過程收斂并結束，否則修改模態(tài)頻率的預估區(qū)間并重新處理。

將上述3種方法的計算約束模態(tài)分布與試驗模態(tài)結果對比，如表3所示。方法1計算結果誤差絕對值最大為22.23%，方法2計算結果誤差絕對值最大為8.95%，方法3計算結果誤差絕對值可控制在9.13%，且此例中迭代2～3次便獲得了滿意效果。3種方法所得計算模態(tài)振型均與試驗振型一致，以方法3為例，其對比結果如圖12～圖14，方法2與方法3均能獲取較高精度的模態(tài)信息，兩者間偏差在1%左右。這表明橡膠襯套單元模型動剛度可由襯套測試動剛度曲線上模態(tài)頻率對應的襯套動剛度來賦值，方法3可達到與方法2相近的精度，且后2種方法均優(yōu)于方法1。

表3　　模態(tài)分析結果對比

參考文獻

1宋立新，侯之超，毛顯紅.汽車后扭力梁工作模態(tài)的試驗分析.第22屆全國結構工程學術會議，烏魯木齊，2013.

2羅明軍，趙永玲，宋立新，等.汽車后扭力梁結構的模態(tài)分析.南昌大學學報（工科版），2013（2）：172～175.

3曹樹謙，張文德，蕭龍翔.振動結構模態(tài)分析：理論、實踐與應用.天津：天津大學出版社，2001.

4陳茜.橡膠襯套動靜態(tài)力學行為研究：［學位論文］.重慶：重慶理工大學，2014.

5趙振東，雷雨成.橡膠元件在汽車懸架中的應用分析.汽車技術，2006（1）：19～22.

（責任編輯斛畔）

修改稿收到日期為2016年1月14日。

中圖分類號：U462.1

文獻標識碼：A

文章編號：1000-3703（2016）06-0015-03

*基金項目：天津市重點項目（14TXSYJC00457）。

Rubber Bushings Dynamic Stiffness Value Assignment Method in Finite Element Modal Analysis of Torsion Beam Suspension

Shi Yuekui，Li Hongliang，Wang Haiyang
（China Automotive Technology and Research Center，Tianjin 300300）

【Abstract】Based on finite element modal analysis of typical automobile torsion beam suspension，we investigate dynamic stiffness value assignment of suspension rubber bushing model in the construction of analysis model，under the prerequisite of rational simplification of suspension boundary conditions，and propose a“Crossover iteration”method. Compared with test modal analysis results，the“Cross iteration”method can control the error of modal frequency calculation within 10%，and the corresponding modal shape has good consistency，the consistency of both frequency and shape are better than the traditional static stiffness equivalent method，proving effectiveness and rationality of the proposed method.

Key words:Rubber bushing，Dynamic stiffness characteristic，Torsion beam suspension，F(xiàn)inite element analysis