尹志新，李端芳，唐 萌，賀琳丹

（廣西大學(xué)機(jī)械工程學(xué)院，廣西南寧530004）

汽車虛擬樣機(jī)技術(shù)，是計(jì)算機(jī)技術(shù)應(yīng)用到機(jī)械行業(yè)中來以后產(chǎn)生的一種新的汽車振動(dòng)性能實(shí)驗(yàn)方法。但虛擬樣機(jī)技術(shù)實(shí)驗(yàn)的結(jié)果，強(qiáng)烈地依賴于所建立的仿真模型的真實(shí)程度。目前大多數(shù)汽車振動(dòng)分析的激振信號(hào)，采用隨機(jī)信號(hào)或變頻正弦激勵(lì)信號(hào)，雖然激振頻帶能包含汽車固有頻率，但激振信號(hào)并不能代替實(shí)際路面振動(dòng)情況。在Matlab平臺(tái)下，利用濾波白噪聲生成法能重構(gòu)時(shí)域路面不平度仿真模型，并考慮行駛速度對(duì)汽車振動(dòng)的影響。

1 路面不平度功率譜密度和濾波白噪聲生成法

1.1 路面不平度功率譜密度

作為車輛振動(dòng)輸入的路面不平度，主要采用路面功率譜密度描述路面不平度統(tǒng)計(jì)特性。根據(jù)文獻(xiàn)[1]的建議，路面功率譜密度Gq(n)用下式作為擬合表達(dá)式：

式中，

n為空間頻率(m-1)；

n0為參考空間頻率，n0=0.1m-1；

Gq(n0)為參考空間頻率n0下的路面功率譜密度值，稱為路面不平度系數(shù)；

W為頻率指數(shù)（w=2），為雙對(duì)數(shù)坐標(biāo)上斜線的斜率，它決定路面功率譜密度的頻率結(jié)構(gòu)。

考慮到車速v對(duì)系統(tǒng)輸入激勵(lì)狀態(tài)的影響，引用文獻(xiàn)[1]的計(jì)算結(jié)果，時(shí)間頻率功率譜密度可表示為

由式(1)、式(2)可得時(shí)間頻率的不平度垂直速度q觶(t)=dq(t)/dt的功率譜密度Gq觶(f)（單位為m2/s）與位移功率譜密度Gq(f)的關(guān)系式，

1.2 濾波白噪聲生成法

濾波白噪聲生成法的基本思想，是將路面不平度的隨機(jī)擾動(dòng)抽象為滿足一定條件的白噪聲，然后經(jīng)濾波器進(jìn)行適當(dāng)變換而擬合輸出為具有指定譜特征的隨機(jī)過程。濾波白噪聲道路模擬生成法具有計(jì)算量小、速度快、模擬精度高的優(yōu)點(diǎn)[2]。在仿真過程中，將式(3)所表示的路面速度功率譜密度中“白噪聲”轉(zhuǎn)換為有色噪聲，從而得到時(shí)域路面不平度的狀態(tài)方程為

式中，

f0為濾波器的下限截止頻率；

q(t)為隨機(jī)路面不平度幅值；

w(t)為均值為0、強(qiáng)度為1的均勻分布單位白噪聲[3]。

2 時(shí)域路面不平度在MATLAB/Simulink中的實(shí)現(xiàn)

Simulink是一個(gè)集建模、仿真和分析動(dòng)態(tài)系統(tǒng)為一體的軟件，是一種基于MATLAB的圖形設(shè)計(jì)環(huán)境。利用Simulink，不但可以輕松地完成對(duì)工程問題的仿真，而且能實(shí)時(shí)地觀察可能出現(xiàn)的結(jié)果。根據(jù)濾波白噪聲生成法的基本思想，采用Simulink模塊庫(kù)中的模塊構(gòu)建路面不平度濾波白噪聲生成法仿真模型，其建立方法如下：

（1）在Simulink的sources模塊庫(kù)中選中 Band-limited White Noise模塊作為線性濾波器的輸入，替代式(4)中的w(t)函數(shù)。設(shè)置Band-limited White Noise模塊參數(shù)，使之成為單位“白噪聲”。具體參數(shù)設(shè)置為 ：Noise power:[1]、Sample time:0.000357、Seed:[23341]

（2）選擇Simulink的User-Defined Function模塊庫(kù)中的S-Function模塊，通過使用M函數(shù)文件編寫S函數(shù)線性濾波程序，實(shí)現(xiàn)對(duì)單位“白噪聲”信號(hào)的變換和擬合。在編寫M函數(shù)文件時(shí)，將路面不平度系數(shù)和車速設(shè)置為S-函數(shù)參數(shù)，以實(shí)現(xiàn)對(duì)不同路面等級(jí)和不同車速下的路面不平度模擬。

（3）數(shù)據(jù)處理單元完成后，將模擬結(jié)果以圖表的形式顯示，選擇Simulink的Commonly Used Blocks模塊庫(kù)中的Scope模塊作為系統(tǒng)的輸出。Xn用來保存仿真采樣數(shù)據(jù)，用于生成仿真時(shí)域路面不平度功率譜。

最終的濾波白噪聲生成法路面不平度的Simulink仿真模型如圖1所示：

圖1 濾波白噪聲生成法路面不平度的Simulink仿真模型

利用建立的Simulink仿真模型輸出路面不平度曲線，選擇車速v=20 m/s(72 km/h)，仿真輸出A級(jí)、B級(jí)、C級(jí)路面的路面不平度曲線，路面不平度系數(shù)Gq(n0)分別為：1.6e-5、6.4e-5、2.56e-4，模擬路長(zhǎng)20 m，故仿真時(shí)間t=1 s。仿真結(jié)果如圖2所示：

圖2 A級(jí)路面不平度曲線

圖3 B級(jí)路面不平度曲線

圖4 C級(jí)路面不平度曲線

從圖中可以看出，隨著路面等級(jí)的降低和車速的提高，路面不平度曲線幅值增加。以上仿真結(jié)果符合車輛實(shí)際道路行駛情況。

3 仿真誤差分析

以國(guó)際標(biāo)準(zhǔn)路面不平度功率譜為參照，對(duì)濾波白噪聲法MATLAB/Simulink仿真模型生成的B級(jí)路面不平度功率譜仿真誤差進(jìn)行了分析。

國(guó)際標(biāo)準(zhǔn)的路面功率譜密度(PSD)可由式(3)直接得到，代入已知參數(shù)(Gq(n0)=6.4×10-5m3、n0=0.1m-1、W=2)，得簡(jiǎn)化公式為

在對(duì)數(shù)坐標(biāo)下，國(guó)際標(biāo)準(zhǔn)路面不平度功率譜為一直線，如圖5所示。

模擬的時(shí)域路面不平度功率譜由直接法估算。直接法又稱周期圖法，把仿真采集數(shù)據(jù)x(n)的N個(gè)觀測(cè)數(shù)據(jù)視為一能量有限的序列，直接計(jì)算X(n)的離散傅立葉變換，得X(k)，然后再取其幅值的平方，并除以N，作為序列X(n)真實(shí)功率譜的估計(jì)。編寫M-file程序計(jì)算在車速v=20 m/s下，濾波白噪聲生成法模擬仿真的B級(jí)路面的路面不平度功率譜。將國(guó)際標(biāo)準(zhǔn)的位移功率譜密度曲線和仿真得出的時(shí)域B級(jí)路面不平度功率譜曲線繪制在同一圖中（圖5）：

圖5 路面不平度功率譜曲線

從圖中可以看出，濾波白噪聲生成法在MATLAB/Simulink平臺(tái)下模擬產(chǎn)生的時(shí)域路面不平度，能夠很好的逼進(jìn)目標(biāo)譜，滿足譜一致的準(zhǔn)則[4]。即模擬譜密度曲線下方的面積與功率譜密度曲線下方的面積近似相等。

4 結(jié)束語

在MATLAB平臺(tái)下利用濾波白噪聲生成法重構(gòu)時(shí)域路面不平度模型具有仿真精度高，可視化，計(jì)算速度快的優(yōu)點(diǎn)。所生成的時(shí)域路面不平度曲線可作為汽車振動(dòng)模擬分析實(shí)驗(yàn)的路面激勵(lì)輸入，以提高其分析精度。

[1]余志生.汽車?yán)碚摚ǖ?版）[M].北京：清華大學(xué)出版社，2002.

[2]彭 佳，何 杰，李旭杰，等.路面不平度隨機(jī)激勵(lì)時(shí)域模型的仿真比較與評(píng)價(jià)[J].解放軍理工大學(xué)學(xué)報(bào)，2009，10（1）:7-82.

[3]潘雙夏，陳助碧，馮培恩.M-File S-函數(shù)在時(shí)域路面不平度建模中的應(yīng)用[J].中國(guó)工程機(jī)械學(xué)報(bào).2006，4（4）:379-384.

[4]唐光武，賀學(xué)鋒，顏永福.路面不平度的數(shù)學(xué)模型及計(jì)算機(jī)模擬研究[J].中國(guó)公路學(xué)報(bào)，2000，13（1）:114-117.

[5]Zhang Yonglin，Zang Jiafan.Numerical simulation of stochastic road process using white noise filt ration[J].Mechanical Systems and Signal Processing.2006，20（2）:363-372.