王長明 曹正林 李響 劉再生 趙晉

（中國第一汽車股份有限公司技術(shù)中心汽車振動(dòng)噪聲和安全控制綜合技術(shù)國家重點(diǎn)實(shí)驗(yàn)室）

基于實(shí)測路譜的比利時(shí)路面載荷預(yù)測與試驗(yàn)對比研究

王長明 曹正林 李響 劉再生 趙晉

（中國第一汽車股份有限公司技術(shù)中心汽車振動(dòng)噪聲和安全控制綜合技術(shù)國家重點(diǎn)實(shí)驗(yàn)室）

基于某C級轎車實(shí)測數(shù)據(jù)在ADAMS/car中建立了整車多體動(dòng)力學(xué)模型和虛擬2D比利時(shí)試驗(yàn)場路面模型，并以該路面耐久性行駛工況為輸入條件，結(jié)合多體仿真、載荷虛擬迭代和道路載荷數(shù)據(jù)采集等技術(shù)獲取了路面車輛的工作載荷。試驗(yàn)結(jié)果表明，比利時(shí)路面載荷預(yù)測仿真結(jié)果與測試結(jié)果在時(shí)域、幅值頻次趨勢基本一致，該方法可有效預(yù)測比利時(shí)路面車輛的動(dòng)態(tài)載荷，同時(shí)可為整車及零部件的疲勞壽命分析和結(jié)構(gòu)優(yōu)化提供載荷輸入條件。

1 前言

目前，國內(nèi)外汽車廠家為達(dá)到減少道路試驗(yàn)及臺(tái)架試驗(yàn)次數(shù)、縮短產(chǎn)品研發(fā)周期及降低研發(fā)成本，通常采用多體仿真的方法為整車及零部件的疲勞分析提供載荷輸入，利用該方法計(jì)算載荷的過程稱之為載荷預(yù)測。載荷預(yù)測技術(shù)包括道路載荷數(shù)據(jù)采集（RLDA）、多體仿真（MBS）、載荷虛擬迭代（LVIT）和虛擬道路載荷數(shù)據(jù)采集（VRLDA）等。

開發(fā)設(shè)計(jì)階段獲取載荷數(shù)據(jù)的途徑之一是CAE仿真，主要采用半解析法和完全解析法[1～5]。半解析法與試驗(yàn)結(jié)果吻合較好，分析結(jié)果精度較高；完全解析法需要建立高精度整車模型和數(shù)字化試驗(yàn)路面模型，因此其在汽車產(chǎn)品設(shè)計(jì)開發(fā)中還處于初步探索階段。

本文以大量實(shí)車樣件測試數(shù)據(jù)和比利時(shí)路面高程數(shù)據(jù)為參數(shù)輸入，在ADAMS/CAR中建立了某C級轎車整車多體動(dòng)力學(xué)模型和虛擬2D比利時(shí)試驗(yàn)場路面模型，用上述2種方法進(jìn)行了后懸架模型動(dòng)態(tài)特性臺(tái)架試驗(yàn)和整車模型比利時(shí)路面的仿真計(jì)算，并分別與臺(tái)架試驗(yàn)結(jié)果和試驗(yàn)場比利時(shí)路面的測試結(jié)果進(jìn)行了對比，驗(yàn)證了該計(jì)算模型的仿真精度。

2 比利時(shí)路面譜和試驗(yàn)車載荷譜采集

2.1 比利時(shí)路面譜采集與數(shù)據(jù)處理

在一汽農(nóng)安汽車試驗(yàn)場耐久試驗(yàn)路的比利時(shí)路段進(jìn)行路面譜采集[6，7]。由于所采集的路面高程數(shù)據(jù)存在零漂、野點(diǎn)和趨勢項(xiàng)等問題，所以必須對數(shù)據(jù)進(jìn)行預(yù)處理。

根據(jù)GB7031—86《車輛振動(dòng)輸入路面不平度表示方法》的建議進(jìn)行數(shù)據(jù)預(yù)處理，GB7031—86中規(guī)定采用垂向位移單邊功率譜Gq（n）和時(shí)間功率譜Gq（f）來描述道路不平度的統(tǒng)計(jì)特性：

式中，n為空間頻率；n0為參考空間頻率，為0.1 m-1；Gq（n0）為參考空間頻率n0時(shí)的路面功率譜密度，即路面不平度系數(shù)；ω為頻率指數(shù)；f為時(shí)間頻率；v為車速。

圖1為數(shù)據(jù)預(yù)處理前、后比利時(shí)路面信號(hào)對比結(jié)果，可見處理后的路面譜數(shù)據(jù)能夠真實(shí)反映實(shí)際道路的路形情況[8～10]。

2.2 試驗(yàn)車載荷譜采集

在一汽農(nóng)安汽車試驗(yàn)場的比利時(shí)路段上進(jìn)行試驗(yàn)車道路載荷數(shù)據(jù)采集（RLDA），主要測試參數(shù)見表1。采集完成后，對載荷譜數(shù)據(jù)進(jìn)行預(yù)處理使其可直接用于模型驗(yàn)證分析。

表1 測試參數(shù)

3 虛擬試驗(yàn)場道路載荷數(shù)據(jù)采集

利用完全解析法實(shí)現(xiàn)虛擬試驗(yàn)場道路載荷數(shù)據(jù)采集（VRLDA）仿真分析，采集流程見圖2。

3.1 整車多體動(dòng)力學(xué)模型

3.1.1 模型參數(shù)獲取

整車多體動(dòng)力學(xué)模型參數(shù)主要來自試驗(yàn)數(shù)據(jù)和整車設(shè)計(jì)目標(biāo)值，為此，對影響仿真精度的主要參數(shù)進(jìn)行了實(shí)車及底盤彈性件的臺(tái)架試驗(yàn)測試，圖3～圖5為部分試驗(yàn)測試裝置和部分試驗(yàn)結(jié)果。

3.1.2 整車仿真模型的建立

根據(jù)獲取的整車多體動(dòng)力學(xué)模型參數(shù)，在ADAMS/CAR中建立了試驗(yàn)車整車多體動(dòng)力學(xué)模型，如圖6所示。整車模型中包含車身、前懸架、后懸架、動(dòng)力總成、轉(zhuǎn)向系統(tǒng)和輪胎模型等，表2為滿載狀態(tài)下整車試驗(yàn)輪荷與計(jì)算模型輪荷對比。

表2 滿載狀態(tài)下整車試驗(yàn)輪荷與模型計(jì)算輪荷對比kN

3.2 數(shù)字化試驗(yàn)路模型

數(shù)字化試驗(yàn)路的輪廓以實(shí)際采集比利時(shí)路段高程信號(hào)（即路面不平度）的預(yù)處理數(shù)據(jù)為輸入，在ADAMS/CAR中建立了350 m長的2D數(shù)字化比利時(shí)路面模型。

3.3 駕駛員模型

為保證數(shù)據(jù)采集過程中車輛通過試驗(yàn)場各典型路段時(shí)均為勻速行駛，根據(jù)試驗(yàn)數(shù)據(jù)采集特點(diǎn)，采用驅(qū)動(dòng)控制文件的方式控制整車模型的運(yùn)動(dòng)行為。驅(qū)動(dòng)控制文件通過駕駛員模型控制轉(zhuǎn)向、油門、制動(dòng)以及變速器和離合器模型參數(shù)來實(shí)現(xiàn)駕駛目標(biāo)事件的整車模型運(yùn)動(dòng)行為[11]。

駕駛員模型的控制目標(biāo)是保證整車模型沿比利時(shí)路面勻速直線行駛，目標(biāo)車速為實(shí)際比利時(shí)道路試驗(yàn)車速的平均值。

3.4 仿真計(jì)算

根據(jù)所建立的整車多體動(dòng)力學(xué)模型、數(shù)字化比利時(shí)路面模型和駕駛員模型對虛擬道路行駛工況進(jìn)行仿真計(jì)算，計(jì)算完成后將每個(gè)鉸接點(diǎn)的載荷提取出來，得到虛擬比利時(shí)道路載荷譜數(shù)據(jù)。

4 仿真與試驗(yàn)的對比驗(yàn)證

4.1 整車模型的動(dòng)態(tài)特性驗(yàn)證

對整車多體動(dòng)力學(xué)模型進(jìn)行了除輪胎模型外的對比試驗(yàn)驗(yàn)證，鑒于篇幅僅以后懸架總成的模型驗(yàn)證為例進(jìn)行說明。

對測試六分力數(shù)據(jù)進(jìn)行編輯得到多軸臺(tái)架試驗(yàn)驅(qū)動(dòng)譜[12]，以該驅(qū)動(dòng)譜為輸入，通過半解析法中的載荷虛擬迭代技術(shù)實(shí)現(xiàn)對后懸架總成的虛擬臺(tái)架試驗(yàn)的驅(qū)動(dòng)；再通過仿真分析結(jié)果與臺(tái)架動(dòng)態(tài)試驗(yàn)測試結(jié)果的對比來驗(yàn)證仿真計(jì)算模型的精度。圖7和圖8分別為臺(tái)架試驗(yàn)和仿真計(jì)算的前下控制臂連桿軸向載荷的時(shí)間歷程信號(hào)及其穿級計(jì)數(shù)對比，圖9為臺(tái)架試驗(yàn)和仿真計(jì)算的3個(gè)連桿軸向載荷的相對損傷對比。

由圖7～圖9可看出，臺(tái)架試驗(yàn)結(jié)果與仿真計(jì)算結(jié)果在時(shí)域、幅值頻次和相對損傷等方面一致性較好，說明由實(shí)車測試參數(shù)建立的后懸架總成多體動(dòng)力學(xué)模型的精度較高，可依據(jù)此進(jìn)行比利時(shí)路面載荷預(yù)測。

4.2 比利時(shí)道路仿真與試驗(yàn)結(jié)果對比

表3為實(shí)際道路測試的車輪垂向力、軸頭加速度和螺旋彈簧位移與比利時(shí)道路仿真計(jì)算結(jié)果的標(biāo)準(zhǔn)差值對比，圖10為左前輪垂向力時(shí)間歷程信號(hào)的測試與計(jì)算結(jié)果對比，圖11為穿級計(jì)數(shù)的測試與計(jì)算結(jié)果對比，圖12為測試與計(jì)算的左前輪垂向力、軸頭加速度和螺旋彈簧位移的功率譜密度的對比，圖13為測試與計(jì)算的車輪垂向力、軸頭加速度和螺旋彈簧位移數(shù)據(jù)結(jié)果的相對損傷值比較。

由表3和圖10～圖13可知，仿真與測試結(jié)果在標(biāo)準(zhǔn)差值、幅值頻次和相對損傷等方面一致性較好，說明由實(shí)測路面譜建立的2D虛擬試驗(yàn)場路面可較好地體現(xiàn)實(shí)際路面不平度的統(tǒng)計(jì)特性。但仿真和測試結(jié)果在時(shí)域和頻域內(nèi)還存在一定的誤差，引起誤差的原因主要是整車模型存在大量剛體簡化，同時(shí)未考慮襯套的動(dòng)態(tài)特性與輪胎模型的動(dòng)力學(xué)特性誤差較大的問題。

表3 測試與計(jì)算結(jié)果的標(biāo)準(zhǔn)差對比結(jié)果

5 結(jié)束語

基于實(shí)車實(shí)測數(shù)據(jù)在ADAMS/car中建立了某C級轎車多體模型和虛擬2D比利時(shí)試驗(yàn)場路面模型，并以該路面耐久性行駛工況為輸入條件，結(jié)合多體仿真、載荷虛擬迭代和道路載荷數(shù)據(jù)采集等技術(shù)，進(jìn)行了路面譜、車輛載荷譜和整車模型參數(shù)的測試與試驗(yàn)，并對其進(jìn)行仿真模擬計(jì)算。通過仿真結(jié)果與測試結(jié)果的對比，表明該方法可有效預(yù)測比利時(shí)路面車輛的動(dòng)態(tài)載荷，通過該方法得到的載荷可作為產(chǎn)品設(shè)計(jì)初期整車與零部件的疲勞壽命分析和結(jié)構(gòu)優(yōu)化的載荷輸入條件。

1李飛，郭孔輝，等.汽車耐久性分析底盤載荷預(yù)測方法研究綜述.科學(xué)技術(shù)與工程，2010，10（24）：5960～5964.

2 Joselito Menezes da Cruz.A Semi-Analytical Method to Generate Load Cases For CAE Durability Using Virtual Vehicle Prototypes.SAE，2003-01-3667.

3徐剛，周鋐等.基于虛擬試驗(yàn)臺(tái)的疲勞壽命預(yù)測研究.同濟(jì)大學(xué)學(xué)報(bào)，2009，37（1）：97～100.

4 Daeoh Kang，Seungjin Heo and Hoiyoung Kim.Virtual Road Profile Modeling Using Equivalent Damage Method For VPG Simulation.SAE，2009-01-0814.

5邵建，董益亮，等.基于多體模型仿真的載荷譜虛擬迭代技術(shù)分析.重慶理工大學(xué)學(xué)報(bào)，2010，24（10）：84～87.

6韓毅.中國典型汽車道路譜數(shù)據(jù)采集與處理研究：[學(xué)位論文].北京：清華大學(xué)，2010.

7劉再生，霍福祥，等.基于路譜輸入的汽車臺(tái)架耐久性試驗(yàn)方法研究.汽車技術(shù)，2010（9）：47～50.

8石峰，段虎明，等.路面不平度的測量.武漢理工大學(xué)學(xué)報(bào)，2010（3）：406～410.

9韓毅，楊殿閣，等.基于路形測量的車載道路譜檢測儀研究與實(shí)現(xiàn).中國汽車工程學(xué)會(huì)年會(huì)論文集，2009: 770～774.

10段虎明，石峰，等.道路路面的刨面的測量研究與實(shí)踐.振動(dòng)與沖擊，2011，30（3）：161～164.

11 Adams/SmartDriver help

12劉再生，霍福祥，等.轎車懸架臺(tái)架多軸疲勞試驗(yàn)載荷開發(fā).汽車技術(shù)，2012（6）：47～50.

（責(zé)任編輯文楫）

修改稿收到日期為2014年3月4日。

Research on Belgian Road Load Prediction and Experimental Study Based on M easured Road Profile

Wang Changming，Cao Zhenglin，Li Xiang，Liu Zaisheng，Zhao Jin
（ChinaFAWCo.，Ltd R&DCenter，StateKey LaboratoryofTechnologyon AutomobileVibrationand Noise&Safetyand Control）

A multi-body dynamics model of a C-class car and a 2D virtual Belgian road model are built in Adams/ car based on the measured data，and the durability driving conditions are used as input condition.The workload of the road vehicle is obtained by the combination of technologies like multi-body simulation，load virtual iteration and the road load data acquisition and etc.The test results indicate that the Belgian road load predication&simulation results and test results are basically consistently in time domain，amplitude frequency trend，this method can effectively predict the dynamic loads of the vehicle on the Belgian road.Meanwhile，this method can also be used in fatigue life analysis and structure optimization of the vehicle and components as input load condition.

Digital Road,Load Prediction,MBS

數(shù)字化路面載荷預(yù)測多體仿真

U461.5+2

A

1000-3703（2014）08-0046-04