范 璐，陳偉波，劉立剛

（泛亞汽車技術(shù)中心有限公司，上海 201201）

前言

結(jié)構(gòu)耐久性能是車輛最重要的性能之一，也是整車開發(fā)中最關(guān)鍵的環(huán)節(jié)之一。為了快速考核懸架的結(jié)構(gòu)耐久性能，開發(fā)出了臺(tái)架道路模擬試驗(yàn)方法，其試驗(yàn)設(shè)備以MTS329多軸道路模擬試驗(yàn)機(jī)為代表。過(guò)去幾十年行業(yè)內(nèi)普遍將試車場(chǎng)道路數(shù)據(jù)采集的載荷作為臺(tái)架試驗(yàn)的載荷源，然而這種方法由于引入了試車場(chǎng)道路數(shù)據(jù)采集，所以臺(tái)架道路模擬試驗(yàn)不得不推遲到工程樣車制造完成之后進(jìn)行。若臺(tái)架道路模擬試驗(yàn)的載荷來(lái)源能夠脫離試車場(chǎng)道路數(shù)據(jù)采集，則臺(tái)架道路模擬試驗(yàn)可至少提前半年進(jìn)行，能有效縮短系統(tǒng)級(jí)和整車級(jí)臺(tái)架結(jié)構(gòu)耐久試驗(yàn)的周期。

隨著計(jì)算機(jī)技術(shù)和虛擬仿真技術(shù)的發(fā)展，多體動(dòng)力學(xué)軟件已能準(zhǔn)確模擬車輛在試車場(chǎng)所受的載荷，因此，可將該載荷作為臺(tái)架道路模擬試驗(yàn)的載荷源［1］。美國(guó)通用汽車公司在多年前就開始掃描試車場(chǎng)3D路面［2］，并將數(shù)字路面導(dǎo)入多體動(dòng)力學(xué)軟件中，生成用于臺(tái)架道路模擬試驗(yàn)的虛擬載荷信號(hào)，在騾子車階段完全取代了試車場(chǎng)道路數(shù)據(jù)采集［3］。歐洲車企也能利用虛擬載荷進(jìn)行前期零部件級(jí)試驗(yàn)，但對(duì)懸架或整車級(jí)多軸道路模擬試驗(yàn)運(yùn)用較少。MTS公司結(jié)合自身的技術(shù)優(yōu)勢(shì)，開發(fā)出了一套HSRC（hybrid system response convergence）混合試驗(yàn)方法，采用虛擬輪胎模型和路面模型與實(shí)際臺(tái)架和車輛相結(jié)合的模式，也無(wú)需試車場(chǎng)道路數(shù)據(jù)采集［4］。國(guó)內(nèi)絕大多數(shù)整車廠目前仍然采用試車場(chǎng)道路數(shù)據(jù)采集的載荷作為臺(tái)架道路模擬試驗(yàn)的載荷源。一些學(xué)者專家嘗試通過(guò)建立虛擬試驗(yàn)臺(tái)預(yù)測(cè)零部件的疲勞壽命［5-6］，但這些方法仍然無(wú)法擺脫試車場(chǎng)道路數(shù)據(jù)采集，如果虛擬試驗(yàn)臺(tái)的載荷能來(lái)自于虛擬仿真分析，則效率會(huì)顯著提高。鑒于此，國(guó)內(nèi)各整車廠陸續(xù)開始掃描試車場(chǎng)路面，生成虛擬路面模型用于仿真分析。榮兵等［7］采用插值方法編程實(shí)現(xiàn)了三維虛擬路面重構(gòu)，并研究發(fā)現(xiàn)曲線規(guī)則格柵（CRG）格式的路面具有更高的計(jì)算效率。孫成智等［8］通過(guò)試車場(chǎng)掃描，建立了3D數(shù)字路面模型，導(dǎo)入整車模型后能有效進(jìn)行載荷預(yù)測(cè)。

本文中詳細(xì)論述了基于虛擬載荷的前懸架和后懸架臺(tái)架道路模擬試驗(yàn)方法，介紹了虛擬載荷生成方法及其數(shù)據(jù)處理流程，并在國(guó)內(nèi)首次將虛擬載荷成功運(yùn)用于物理臺(tái)架道路模擬試驗(yàn)，通過(guò)對(duì)比物理試驗(yàn)的失效結(jié)果驗(yàn)證了此試驗(yàn)方法的有效性。

1 基于虛擬載荷的試驗(yàn)方法概述

傳統(tǒng)基于試車場(chǎng)道路載荷的臺(tái)架道路模擬試驗(yàn)方法如圖1所示，在試車場(chǎng)道路數(shù)據(jù)采集之前，須對(duì)物理工程樣車的關(guān)鍵零部件貼應(yīng)變片，必要時(shí)還須在車身和軸頭安裝加速度傳感器、在輪心安裝六分力測(cè)量輪等。工程樣車在試車場(chǎng)采譜時(shí)，受天氣、路面、駕駛員和車輛自身狀態(tài)等條件的限制，往往需要耗費(fèi)很長(zhǎng)時(shí)間。采譜過(guò)程全部結(jié)束之后，才能進(jìn)行數(shù)據(jù)處理和臺(tái)架道路模擬試驗(yàn)。

圖1 基于試車場(chǎng)道路載荷的試驗(yàn)方法

基于虛擬載荷的臺(tái)架道路模擬試驗(yàn)完全省去了試車場(chǎng)道路數(shù)據(jù)采集這一環(huán)節(jié)，試驗(yàn)載荷來(lái)自于虛擬仿真分析。圖2示出此試驗(yàn)方法的基本流程。模型對(duì)標(biāo)包括靜態(tài)對(duì)標(biāo)和動(dòng)態(tài)對(duì)標(biāo)，其中靜態(tài)對(duì)標(biāo)是在相同激勵(lì)下將K&C試驗(yàn)臺(tái)的響應(yīng)信號(hào)與初始車輛模型響應(yīng)信號(hào)進(jìn)行對(duì)標(biāo)，而動(dòng)態(tài)對(duì)標(biāo)則是基于多軸道路模擬試驗(yàn)臺(tái)。若初始車輛模型的響應(yīng)信號(hào)和臺(tái)架的響應(yīng)信號(hào)差別較大，則須進(jìn)一步修正模型，若初始車輛模型的響應(yīng)信號(hào)與臺(tái)架響應(yīng)信號(hào)一致，則該模型可用來(lái)生成虛擬載荷。結(jié)合駕駛員模型、該車型對(duì)應(yīng)的輪胎模型和試車場(chǎng)3D數(shù)字路面模型，利用多體動(dòng)力學(xué)軟件可輸出虛擬載荷時(shí)間歷程。經(jīng)過(guò)數(shù)據(jù)處理將虛擬載荷轉(zhuǎn)化成迭代目標(biāo)信號(hào)，可進(jìn)行零部件級(jí)試驗(yàn)、系統(tǒng)級(jí)試驗(yàn)和整車級(jí)試驗(yàn)。由于篇幅限制，本文中僅討論懸架子系統(tǒng)的試驗(yàn)方法。

圖2 基于虛擬載荷的試驗(yàn)方法

2 虛擬載荷生成

虛擬載荷由多體動(dòng)力學(xué)軟件Motion View生成，虛擬模型包括駕駛員模型、車輛模型、輪胎模型和路面模型，構(gòu)成“人-車-路”閉環(huán)系統(tǒng)。

駕駛員模型控制轉(zhuǎn)向盤轉(zhuǎn)角，保證車輛模型按照預(yù)定軌跡和車速在路面上正確行駛，其控制方程為

式中：β為轉(zhuǎn)向盤轉(zhuǎn)角；Sr為傳動(dòng)比；δ為前輪轉(zhuǎn)角。

根據(jù)2自由度車輛模型，可得前輪轉(zhuǎn)角與軌道誤差之間的關(guān)系［9］為

式中：a為整車質(zhì)心至前軸的距離；b為整車質(zhì)心至后軸的距離；m為整車質(zhì)量；v為車輛行駛速度；Cf為前軸側(cè)偏剛度；Cr為后軸側(cè)偏剛度；d為前視距離；Δy為即時(shí)狀態(tài)橫向位移誤差，即軌道誤差。

車輛模型基于硬點(diǎn)數(shù)模坐標(biāo)、襯套剛度曲線、緩沖塊剛度曲線、減振器阻尼特性曲線等參數(shù)建立，其中前緩沖塊剛度曲線如圖3所示，前減振器阻尼特性曲線如圖4所示。

輪胎模型采用基于柔性環(huán)假設(shè)的3D非線性FTire模型，仿真頻率高達(dá)200 Hz，且能處理極其粗糙的路面特征。通過(guò)多種類型的輪胎特性試驗(yàn)得到其力學(xué)特性數(shù)據(jù)，采用配套辨識(shí)工具進(jìn)行輪胎模型參數(shù)辨識(shí)。表1詳細(xì)列出了需要進(jìn)行的輪胎特性試驗(yàn)種類。

圖3 前緩沖塊剛度曲線

圖4 前減振器阻尼特性曲線

表1 輪胎特性試驗(yàn)種類

虛擬路面模型源自國(guó)內(nèi)某試車場(chǎng)，建模方式主要有3種，包括根據(jù)施工圖紙的路面幾何尺寸進(jìn)行建模，以及根據(jù)2D或3D路面高程掃描數(shù)據(jù)進(jìn)行建模。根據(jù)施工圖紙建模適用于橫截面形狀規(guī)則的路面，如方坑路等。對(duì)于橫截面形狀不規(guī)則的路面，無(wú)法根據(jù)施工圖紙直接建模，因此需要對(duì)路面高程進(jìn)行掃描，2D掃描通過(guò)激光位移傳感器對(duì)路面進(jìn)行測(cè)量，最終獲取路面高程信息，生成2D路面模型。3D掃描采用2個(gè)高分辨率激光掃描探頭結(jié)合高精度慣性補(bǔ)償系統(tǒng)，對(duì)路面進(jìn)行360°環(huán)繞掃描，獲取指定分辨率的點(diǎn)云數(shù)據(jù)，后處理生成3D路面模型。常用的路面模型文件有等效容積（RDF）格式、規(guī)則柵格（RGR）格式和曲線規(guī)則柵格（CRG）格式。鑒于計(jì)算效率，RDF格式僅適用于橫截面規(guī)則的路面，而對(duì)于橫截面不規(guī)則的路面，需要采用RGR格式或者CRG格式。RGR格式路面通過(guò)一系列離散矩形節(jié)點(diǎn)的高度來(lái)模擬路面高程，這些節(jié)點(diǎn)在水平面均為等間距分布。由于柵格的規(guī)則性，路面文件中每個(gè)節(jié)點(diǎn)的水平面坐標(biāo)值都無(wú)需存儲(chǔ)，只須儲(chǔ)存水平面坐標(biāo)增量、柵格原點(diǎn)坐標(biāo)、柵格旋轉(zhuǎn)角度和節(jié)點(diǎn)高度值。如果實(shí)際道路為曲線路徑，RGR文件數(shù)據(jù)可被帶有路徑軌跡的3D樣條路面文件引用，這樣規(guī)則柵格就可沿著路徑軌跡鋪設(shè)。CRG格式路面類似于帶路徑中心線的RGR格式路面，不同點(diǎn)在于節(jié)點(diǎn)在路面橫向的間距值可任意指定，從而能更靈活地定義各種路面。表2匯總了此次針對(duì)國(guó)內(nèi)某試車場(chǎng)的各種路面建模方法。

表2 虛擬路面建模方法

將駕駛員模型、輪胎特性文件、虛擬路面模型與車輛模型集成，構(gòu)成完整的虛擬模型，進(jìn)行仿真分析，生成虛擬載荷。圖5展示了Motion View中的虛擬模型。

3 數(shù)據(jù)處理與載荷相關(guān)性分析

圖5 Motion View虛擬模型

對(duì)于由多體動(dòng)力學(xué)軟件生成的虛擬載荷，須先進(jìn)行數(shù)據(jù)處理，才能得到迭代目標(biāo)信號(hào)，用于臺(tái)架試驗(yàn)迭代。表3和表4分別列出了前、后懸架系統(tǒng)的虛擬載荷通道。數(shù)據(jù)處理的流程如圖6所示，首先修改虛擬載荷信號(hào)的頭文件，使時(shí)域分辨率、頻域分辨率、通道名稱和通道順序等滿足試驗(yàn)要求。對(duì)于模型與臺(tái)架單位定義不同的通道，須對(duì)其進(jìn)行單位轉(zhuǎn)換。由于試驗(yàn)臺(tái)架繞Z軸的轉(zhuǎn)角很小，無(wú)法全面模擬車輛在試車場(chǎng)的轉(zhuǎn)彎運(yùn)動(dòng)，試驗(yàn)時(shí)將車輛轉(zhuǎn)向機(jī)鎖死在中間位置，只模擬繞Z軸的力矩，并不模擬繞Z軸的運(yùn)動(dòng)，因此須對(duì)虛擬載荷信號(hào)進(jìn)行轉(zhuǎn)向修正，修正公式為

式中：Fx，F(xiàn)y分別為修正前的輪心縱向力和側(cè)向力；FxE，F(xiàn)yE分別為修正后的輪心縱向力和側(cè)向力。

表3 前懸架系統(tǒng)虛擬載荷通道

表4 后懸架系統(tǒng)虛擬載荷通道

圖6 數(shù)據(jù)處理流程

為滿足臺(tái)架試驗(yàn)快速驗(yàn)證的需求，往往須對(duì)信號(hào)進(jìn)行壓縮，刪除對(duì)損傷貢獻(xiàn)量小的信號(hào)，以加速試驗(yàn)周期。采用基于 Miner法則的偽損傷計(jì)算方法［10］，以保留一定百分比的原始信號(hào)偽損傷為標(biāo)準(zhǔn)，對(duì)信號(hào)進(jìn)行時(shí)域壓縮，其中偽損傷保留的百分比由企業(yè)標(biāo)準(zhǔn)確定。

最后按要求提取通道，通道順序需與試驗(yàn)站臺(tái)設(shè)置文件中的通道順序一致。

為驗(yàn)證基于多體動(dòng)力學(xué)軟件生成的虛擬載荷的準(zhǔn)確性，對(duì)虛擬載荷與實(shí)際道路載荷進(jìn)行相關(guān)性分析。由于載荷通道數(shù)較多，此處僅以比利時(shí)路面左前輪心縱向力為例進(jìn)行說(shuō)明。表3列出了兩種載荷源下左前輪心縱向力統(tǒng)計(jì)學(xué)特性的對(duì)比，圖7～圖9分別展示了雨流計(jì)數(shù)、穿級(jí)計(jì)數(shù)和功率譜密度的對(duì)比?？梢?jiàn)，虛擬載荷與實(shí)際道路載荷相關(guān)性很好。

圖7 雨流計(jì)數(shù)對(duì)比

圖8 穿級(jí)計(jì)數(shù)對(duì)比

圖9 功率譜密度對(duì)比

4 實(shí)車試驗(yàn)驗(yàn)證

實(shí)車試驗(yàn)驗(yàn)證的總體思路為：首先根據(jù)上述多體動(dòng)力學(xué)軟件生成的虛擬載荷信號(hào)，進(jìn)行一輪前、后懸架系統(tǒng)耐久試驗(yàn)；其次對(duì)相同狀態(tài)的零件，通過(guò)傳統(tǒng)試車場(chǎng)數(shù)據(jù)采集，得到試車場(chǎng)的真實(shí)道路載荷，再進(jìn)行一輪前、后懸架系統(tǒng)耐久試驗(yàn)；最后對(duì)比分析兩輪耐久試驗(yàn)的零件失效位置和發(fā)生失效的時(shí)間，若一致，則說(shuō)明基于虛擬載荷的懸架道路模擬試驗(yàn)方法是有效的。

4.1 夾具設(shè)計(jì)

試驗(yàn)車輛的前懸架為麥弗遜結(jié)構(gòu)，后懸架為扭力梁結(jié)構(gòu)。對(duì)于虛擬模型里輸出的懸架與車身分界面間的作用力，在實(shí)車試驗(yàn)中可用力傳感器進(jìn)行測(cè)量。在傳統(tǒng)的懸架試驗(yàn)夾具基礎(chǔ)上進(jìn)行改進(jìn)，在懸架與夾具的分界面上加裝力傳感器，其所測(cè)得的力信號(hào)能準(zhǔn)確表征懸架與車身分界面的力。圖10和圖11分別示出前懸架和后懸架試驗(yàn)中力傳感器的安裝位置，其中前懸架有2個(gè)力傳感器安裝點(diǎn)，后懸架有6個(gè)力傳感器安裝點(diǎn)。

圖10 前懸架試驗(yàn)中力傳感器安裝位置

圖11 后懸架試驗(yàn)中力傳感器安裝位置

4.2 模型辨識(shí)

對(duì)系統(tǒng)施加特定白噪聲激勵(lì)信號(hào)，得到此激勵(lì)下的響應(yīng)信號(hào)，采用頻率響應(yīng)函數(shù)法進(jìn)行模型辨識(shí)，求解公式為

式中：H（f）為系統(tǒng)模型；Gyx（f）為響應(yīng)信號(hào)和激勵(lì)信號(hào)的互功率譜；Gxx（f）為激勵(lì)信號(hào)的自功率譜。

4.3 迭代

迭代的目的是得到試驗(yàn)臺(tái)架的驅(qū)動(dòng)信號(hào)，使被試件在此驅(qū)動(dòng)下的響應(yīng)信號(hào)與迭代目標(biāo)信號(hào)一致，其基本流程如圖12所示。對(duì)于懸架系統(tǒng)試驗(yàn)，迭代質(zhì)量評(píng)價(jià)指標(biāo)主要包括各通道均方根值誤差、最終響應(yīng)信號(hào)與迭代目標(biāo)信號(hào)的偽損傷值對(duì)比、幅值對(duì)比和功率譜密度對(duì)比等。

圖12 迭代流程圖

傳統(tǒng)的懸架道路模擬試驗(yàn)迭代通道通常包括由測(cè)量輪測(cè)出的輪心六分力信號(hào)和由關(guān)鍵零部件應(yīng)變片測(cè)出的應(yīng)變信號(hào)，但由于在虛擬模型中計(jì)算應(yīng)變耗時(shí)長(zhǎng)、效率低，傳統(tǒng)的迭代方法和迭代通道選擇需要改變。懸架與夾具連接面上的力傳感器信號(hào)能直接反映零件的真實(shí)受力狀態(tài)，因此選用此力信號(hào)作為首要迭代控制目標(biāo)，測(cè)量輪的六分力信號(hào)作為次要迭代控制通道。

應(yīng)用MTS公司的 RPC軟件［11］對(duì)臺(tái)架進(jìn)行迭代，迭代路面包括比利時(shí)路、方坑路、扭曲路等。由于篇幅有限，文中僅展示損傷占比最大的比利時(shí)路面的迭代效果。表5列出了比利時(shí)路面各通道最終響應(yīng)信號(hào)與迭代目標(biāo)信號(hào)的偽損傷值之比。由表可見(jiàn)，各通道偽損傷比值均在60%～180%之間，滿足企業(yè)規(guī)范要求。圖13截取了部分左前輪心縱向力時(shí)域信號(hào)，可見(jiàn)時(shí)域上最終響應(yīng)信號(hào)與迭代目標(biāo)信號(hào)一致。圖14展示了左前輪心縱向力功率譜密度對(duì)比，可見(jiàn)頻域上最終響應(yīng)信號(hào)與迭代目標(biāo)信號(hào)一致。其余31個(gè)通道的迭代結(jié)果與左前輪心縱向力類似。分析表明迭代效果很好，臺(tái)架最終響應(yīng)信號(hào)復(fù)現(xiàn)了迭代目標(biāo)信號(hào)。

表5 迭代各通道偽損傷對(duì)比 %

4.4 耐久試驗(yàn)

圖13 左前輪心縱向力時(shí)域信號(hào)對(duì)比

圖14 左前輪心縱向力功率譜密度對(duì)比

提取各路面迭代的最終驅(qū)動(dòng)信號(hào)，按照企業(yè)標(biāo)準(zhǔn)中規(guī)定的各路面循環(huán)次數(shù)進(jìn)行一輪基于虛擬載荷的耐久試驗(yàn)。同時(shí)，采用相同狀態(tài)的懸架系統(tǒng)零件，基于傳統(tǒng)的試車場(chǎng)道路載荷再進(jìn)行一輪臺(tái)架耐久試驗(yàn)。試驗(yàn)臺(tái)架如圖15所示。兩輪試驗(yàn)均發(fā)現(xiàn)左右兩側(cè)前控制臂襯套、后橋左右兩側(cè)襯套和后橋本體焊縫等處開裂，其中控制臂襯套開裂如圖16所示。表6列出了兩輪耐久試驗(yàn)中3處開裂發(fā)生的時(shí)間。由表可見(jiàn)，基于虛擬載荷的懸架試驗(yàn)不僅能完全復(fù)現(xiàn)失效位置，且失效發(fā)生的時(shí)間也非常接近，相對(duì)誤差＜7%。

圖15 試驗(yàn)臺(tái)架

圖16 控制臂襯套開裂

表6 耐久試驗(yàn)失效發(fā)生時(shí)間對(duì)比

為進(jìn)一步揭示虛擬載荷與試車場(chǎng)真實(shí)道路載荷的當(dāng)量關(guān)系，在后橋7處高應(yīng)力區(qū)貼應(yīng)變片（見(jiàn)圖17），分別用基于虛擬載荷的臺(tái)架驅(qū)動(dòng)和基于試車場(chǎng)載荷的臺(tái)架驅(qū)動(dòng)作為激勵(lì)，采集應(yīng)變片的時(shí)域信號(hào)，并進(jìn)行偽損傷分析計(jì)算。以基于試車場(chǎng)載荷的偽損傷值為基準(zhǔn)，表7列出了后橋7處應(yīng)變信號(hào)偽損傷值對(duì)比，可見(jiàn)虛擬載荷可用來(lái)表征試車場(chǎng)真實(shí)道路載荷，臺(tái)架試驗(yàn)結(jié)果能真實(shí)反映試車場(chǎng)試驗(yàn)結(jié)果。

圖17 后橋貼片位置

表7 后橋應(yīng)變信號(hào)偽損傷對(duì)比 %

5 結(jié)論

基于虛擬載荷的懸架道路模擬試驗(yàn)方法能準(zhǔn)確復(fù)現(xiàn)被試件的耐久失效情況，將此試驗(yàn)方法加入整車開發(fā)流程中，可有效縮短整車開發(fā)周期，且由于不需要用于試車場(chǎng)道路數(shù)據(jù)采集的工程樣車，能為項(xiàng)目節(jié)省至少幾百萬(wàn)元人民幣的制造費(fèi)用。