程穩(wěn)正 曹正林 趙晉

（中國第一汽車股份有限公司技術中心 汽車振動噪聲和安全控制綜合技術國家重點試驗室）

乘用車整車誤用試驗軸頭載荷仿真研究

程穩(wěn)正 曹正林 趙晉

（中國第一汽車股份有限公司技術中心 汽車振動噪聲和安全控制綜合技術國家重點試驗室）

以某C級轎車為研究對象，建立考慮強非線性特征的整車極限誤用工況全柔性仿真模型，按規(guī)范進行兩種誤用試驗的仿真研究，得到軸頭載荷，并與同工況實車測試進行對比研究。模型計算精度達到設計要求，表明這種在車型開發(fā)早期獲取誤用試驗軸頭載荷的仿真方法有效可行。

1 前言

為確保誤用工況下車輛的可靠性，準確獲得軸頭力是強度設計的關鍵。文獻[1]、[2]提出的軸頭力經驗公式均基于耐久性路面測試數(shù)據，不適用于誤用工況。分析不同車型誤用試驗數(shù)據后發(fā)現(xiàn)，在相同的試驗方法下軸頭G載荷相差很大，數(shù)據沿用性差。在數(shù)據積累有限的情況下，為降低誤用試驗成本、提高試驗通過率，迫切需要在車型開發(fā)早期能夠仿真誤用試驗，獲得載荷水平及失效模式。

誤用時存在強非線性，采用傳遞函數(shù)迭代軸頭力的方法存在較大誤差[3]，因此對于整車誤用試驗的仿真一直是個難點。Ray Hanley、Ch Oertel、Manfred Baecker等分別研究了誤用試驗的輪胎行為，驗證了仿真輪胎模型[4～6]；寶馬公司的Micro Hauke開展了整車誤用行為仿真研究，對部分連接點載荷進行了比較[7]；方杰等進行了輪胎機械特性虛擬試驗場研究[8]；曹正林等建立了虛擬試車場動力學輪胎模型[9]，開展了虛擬試車場的轎車懸架耐久性強化試驗仿真研究[10]。

本文以某C級轎車為研究對象，建立考慮強非線性特征的整車模型，根據試驗規(guī)范仿真了兩種誤用試驗，并與測試結果進行對比研究。

2 非線性處理方法

2.1 輪胎模型

以實車輪胎為研究對象，進行輪胎徑向剛度、扭轉特性和側偏特性的仿真分析，將計算結果與相同工況下的輪胎試驗臺測試結果進行對比，以驗證所建立的輪胎模型剛度特性的有效性。圖1為輪胎扭轉剛度仿真模型。圖2為試驗與仿真比較分析結果。具體建模及驗證過程見文獻[9]。

2.2 橡膠襯套

橡膠襯套在誤用工況下的特性影響與之相連接結構的過載行為，需要考慮其在非線性下的力學行為。橡膠材料采用Mooney-Rivlin本構模型，相關的本構參數(shù)可根據相關測試數(shù)據擬合獲取。

擬合過程如下：

a.按試驗方法進行襯套主要受力方向（徑向、軸向）剛度有限元分析，以試驗測試剛度為目標，調整材料參數(shù)，直至達到合理的誤差范圍內；

b.以第1步得到的材料參數(shù)進行非主要受力方向（扭轉方向）剛度有限元分析，驗證設置參數(shù)的有效性；

c.重復迭代，直至3個方向的剛度分析結果均在合理的誤差范圍內；

d.擬合得到的本構參數(shù)作為橡膠材料的屬性定義。

圖3為襯套徑向剛度有限元分析模型。圖4為襯套有限元分析的變形情況。

表1為剛度分析結果與試驗結果的比較?？芍瑒偠日`差在可接受的范圍內，因此得到的材料本構參數(shù)可作為橡膠襯套建模的材料屬性輸入。

表1 襯套剛度

對于缺少試驗數(shù)據的橡膠襯套材料，可根據式（1）～式（4）計算得到［11］。

式中，E為彈性模量；μ為泊松比，對于橡膠材料μ≈0.5；G′為剪切模量；C1、C2為Mooney-Rivlin材料參數(shù)；HA為橡膠的邵爾硬度。

2.3 滑柱總成

前、后懸架的滑柱總成是懸架建模的關鍵，其包括減振器、彈簧和緩沖塊。減振器部分簡化為一個Cylin?der運動副和一個減振器單元，減振器單元的屬性來自于減振器的外特性曲線，圖5為前、后減振器懸架調校的速度-載荷曲線。但圖5中曲線速度只到±1 m/s，而誤用工況減振器的相對速度能達到7.5 m/s［7］，因此需要將該曲線進行延伸以考慮極限沖擊情況。彈簧部分簡化為一個線性剛度的彈簧單元。緩沖塊部分簡化為一個非線性的彈簧單元。整個滑柱總成模型見圖6。

彈簧的預載可以通過設定彈簧單元的初始偏移值來實現(xiàn)，初始偏移值計算方法如下：

式中，F(xiàn)為彈簧預載；k為彈簧剛度；L0為彈簧自由長度；l為模型中定義彈簧單元的兩個節(jié)點距離；l0為彈簧單元的初始偏移值。

3 整車模型與仿真

3.1 底盤模型

將副車架、控制臂、轉向節(jié)、滑柱、輪胎等零件組裝到一起構成底盤模型，圖7為前、后軸總成模型。各個零件之間用球頭、橡膠襯套連接，球頭采用球頭單元，橡膠襯套總成采用實體網格建模，襯套總成的內、外套管與各自相連接零件定義連接關系。

3.2 整車模型組裝及配重

整車的仿真模型主要包括車身、底盤、動力總成等。依據車型各主要子系統(tǒng)質量進行模型配重，模型不包括內飾、電子電器、被動安全等部分，這些質量用質量單元將其附加到相關位置的節(jié)點上。滿載狀態(tài)5個乘員按每人75 kg計算，后備箱貨物按60 kg計算。配重后的模型質心位置與實車位置比較見表2。

表2 質心位置比較 mm

3.3 路面模型

依據試驗規(guī)范對路面的要求建立路面模型，共仿真兩種路面，其分別為橫向溝槽與橫木障礙。路面模型如圖8所示。

3.4 誤用試驗仿真

將整車模型及路面模型組裝在一起，定義相關的載荷條件與控制參數(shù)，建立整車仿真模型。整個計算模型單元為442 299個，節(jié)點數(shù)為337 982個。為了節(jié)省運算時間，將部分零件及總成定義為剛體材料，整車誤用仿真模型見圖9。

施加整車G載荷，定義行駛速度，預留一定的穩(wěn)定G載荷的時間，載荷結果輸出時間間隔0.005 s，圖10為橫向溝槽試驗過程仿真。

4 實車試驗

在試驗車底盤零件上粘貼電阻應變片，模擬實車受力狀態(tài)對載荷進行標定，在此基礎上進行實車誤用試驗輪心載荷測試。根據該轎車前、后懸架的部件特點，對表征車輪載荷的相應部件測點進行篩選和布置。主要測試載荷為前輪縱向載荷、垂向載荷和后輪垂向載荷。其余方向由于部件測點應變與車輪載荷相對關系不明顯而未進行測試。前輪縱向載荷測點位置為前下控制臂中間部位，垂向載荷表征點為前下控制臂球頭附近；后輪垂向載荷表征點為彈簧臂中間部位。

根據誤用工況試驗規(guī)范要求，實車以指定的行駛車速通過橫木障礙和橫向溝槽路面，每個工況進行兩次試驗，并進行數(shù)據采集。

5 仿真與試驗結果比較

5.1 載荷歷程

圖11為橫向溝槽誤用試驗后左輪垂向軸頭力仿真與測試結果的比較。由于測試信號為動態(tài)載荷，為便于比較將其數(shù)值加上靜態(tài)輪荷得到絕對載荷。從載荷峰值看，仿真結果與測試結果相差不大。由于減振器在相對速度較大時的特性由經驗給出，與實際性能存在差別，載荷衰減階段仿真的軸頭載荷歷程存在振蕩，因此與測試結果存在差別。從時間軸看，從輪胎再次接觸路面到出現(xiàn)極大峰值，仿真結果與測試結果存在0.03 s的時間差，表明輪胎變形響應時間和滑柱總成下跳所處的位置與實際存在差別。仿真模型較準確的獲得了軸頭力峰值，但要想準確模擬整個動態(tài)過程，模型仍需進一步調試。

5.2 相對輪荷

表3統(tǒng)計了載荷峰值時仿真和試驗的G載荷（兩次試驗的平均值）的比值?？芍瑱M向溝槽工況仿真結果與試驗結果吻合良好，而橫木障礙工況下前輪精度較高，后輪垂向載荷仿真結果比試驗結果大1.35倍?？傮w看，載荷峰值仿真結果具有較高的精度。

表4為兩種試驗仿真G載荷的比值?？芍瑱M向溝槽工況均明顯大于橫木障礙工況。對于懸架系統(tǒng)強度設計，應重點關注橫向溝槽工況，結構設計時以該工況的軸頭力載荷進行零部件的強度校核。

表3 仿真與試驗G載荷比值

表4 兩種試驗仿真G載荷比值

需要指出的是，誤用工況車輛通過障礙時產生的軸頭力受到輪胎、懸架參數(shù)、軸距、質心等多個影響因素的影響，不同級別的車型、不同類型的懸架、不同大小的輪胎均會影響軸頭載荷的大小，因此需要不斷完善仿真模型及不斷的數(shù)據積累才能找到較為可靠的誤用工況軸頭載荷。

6 結束語

建立了某C級轎車的整車誤用工況虛擬試驗模型，進行了兩種誤用試驗的仿真分析，獲得了軸頭載荷值，并同試驗結果進行了比較，結果表明：

a.橫向溝槽試驗工況前、后輪軸頭載荷具有較高的仿真精度；橫木障礙工況前輪軸頭載荷精度較高，后輪垂向載荷大于試驗結果。

b.對比兩種試驗的分析結果，橫向溝槽工況存在較大的垂向軸頭力和軸頭縱向力，對于懸架系統(tǒng)強度設計，應重點關注該工況。

c.通過建立整車誤用工況虛擬試驗，仿真結果達到預期精度要求，該方法可行。

1 Cheng Chao，Sudhakar Medepalli.Statistical Relationship between Corner Weight and Spindle Load.SAE papers 2007-01-1346.

2 Jacob Szczudlak，Mark Villaire，Rajesh Rao.Developing Ge?neric Load Cases by Defining Maximum Spindle Loads as a Function of Corner Weight&Tire Sidewall Height.SAE pa?pers 2013-01-1435.

3 靳暢，周鋐，慕樂.基于頻響函數(shù)法的路面激勵下車輪軸頭力的估計.汽車工程，2014,36（4）:469～474.

4 Ray Hanley，David Crolla，Micro Hauke.Tire Modeling for Misuse Situations.SAE papers 2001-01-0148.

5 Manfred Baecker，Axel Gallrein，Hitoshi Haga.A Tire Mod?el for Very Large Tire Deformations and its Application in Very Severe Events.SAE papers 2010-01-0373.

6 Ch Oertel，A Fandre.Tire Model RMOD-K 7 and Misuse Load Cases.SAE papers 2009-01-0582.

7 Micro Hauke.Simulation of full vehicle misuse behavior.SAE papers 2004-01-0192.

8 方杰，吳光強.輪胎機械特性虛擬試驗場.計算機仿真，2007,6:243～247.

9 曹正林，李駿，郭孔輝.虛擬試車場動力學輪胎模型技術研究.汽車技術，2012（4）:13～16.

10 曹正林，李駿，郭孔輝.基于虛擬試車場的轎車懸架耐久性強化試驗仿真研究.機械工程學報，2012（5）：122～127.

11 戴永謙，宋希庚，薛冬新，等.發(fā)動機懸置軟墊初始斷裂位置及斷裂面預判.農業(yè)機械學報，2005（10）：23～25.

（責任編輯晨 曦）

修改稿收到日期為2015年6月1日。

6 結束語

a.建立某農用車發(fā)動機GT-POWER模型并進行試驗驗證，最大功率計算值與試驗值誤差為2.48%，模型較為準確。

b.建立整車冷卻系統(tǒng)GT-COOL模型并進行適當簡化，將散熱器模型作多孔介質區(qū)域處理，求取其芯部的流動特性參數(shù)和換熱特性參數(shù)；基于GT-POWER模型，求取冷卻系統(tǒng)中發(fā)動機模型的熱邊界條件和散熱器換熱特性參數(shù)。仿真求得散熱器出水口溫度為362.7 K，試驗測得散熱器出水口溫度為362 K，誤差在0.2%以內。

c.將散熱器模型和發(fā)動機模型在STAR-CCM+軟件中耦合，分析水泵流量對系統(tǒng)特性參數(shù)的影響，得出水泵流量在0.8 kg/s時該散熱器性能得到充分發(fā)揮且燃燒室溫度比較合適。

參考文獻

1 劉鵬，劉志博.楊小松.汽車冷卻系統(tǒng)匹配性探討.重型汽車,2013（2）：16～17.

2 楊小松.汽車冷卻系統(tǒng)匹配性探討.汽車研究與開發(fā),1999（2）：24～26.

3 徐立平.基于AMESim發(fā)動機冷卻系統(tǒng)匹配仿真分析.機械研究與應用,2014,129(27):42～43.

4 高思遠，趙長祿，李云龍，等.冷卻液溫度對柴油機熱功轉換效率的影響.農業(yè)機械學報,2012,43(2):29～32.

5 韓愷,趙長祿,張付軍，等.基于一維CFD仿真技術的裝甲車輛冷卻空氣系統(tǒng)設計.北京理工大學學報,2010,30（4）:416～420.

6 郭學勤,劉瑞祥.基于GT-COOL的對發(fā)動機冷卻系散熱能力影響研究.農業(yè)裝備與車輛工程,2010,227（6）: 30～32.

7 Jae Su Kim,Allan Wang.Transient Conjugate CFD Simula?tion of the Radiator Thermal Cycle.SAE Paper,2006-01-1157.

8 黎蘇，鄭清平，張坤鵬,等.天然氣發(fā)動機氣缸蓋流固耦合傳熱研究.農業(yè)機械學報，2012,43（7）:24～26.

9 童正明,陳丹,梁淑君,等.管帶式汽車散熱器試驗及數(shù)值模擬研究.汽車技術,2013,（10）:51～54.

10 郭良平,張衛(wèi)正,王長圓,等.柴油機氣缸蓋傳熱規(guī)律研究.北京理工大學學報,2011,31（3）:278～282.

11 白敏麗,丁鐵新,呂繼祖.活塞組-氣缸套耦合傳熱模擬.內燃機學報,2005,23（2）:169～175.

12 愈小莉，鄭飛，嚴兆大．內燃機氣缸體內表面穩(wěn)態(tài)傳熱邊界條件的研究．內燃機學報，1987，5（4）：324～332.

13《STAR CCM+User Guide》.CD-adapco,2011.

（責任編輯晨 曦）

修改稿收到日期為2015年8月1日。

Simulation Research on Spindle Loads of Passenger Car in Misuse Tests

Cheng Wenzheng,Cao Zhenglin,Zhao Jin
（State Key Laboratory of Technology on Automobile NVH&Safety Control,China FAW Co.,Ltd R&D Center）

In this study,a C-class car is taken as a research object,a full flexible vehicle model for misused tests simulation is established with consideration of the extremely nonlinear characteristic.Then,simulation research is conducted for two misuse tests according to specifications,and the spindle load is obtained,which is compared with that obtained from real vehicle test.The results show that the model calculation precision satisfies design requirement, indicating that the method to obtain the spindle loads in misuse test during the early development phase is feasible.

Misuse test，Spindle load，Virtual test，Passenger car

誤用試驗 軸頭載荷 虛擬試驗 乘用車

U467

A

1000-3703（2015）11-0034-05