黃志強，鄭旺輝

(北京機(jī)械設(shè)備研究所，北京 100854)

0 引言

隨機(jī)路面中激勵點的高程值變化可以用另一個詞匯來表達(dá)：路面不平度。路面不平度[1]定義為道路表面相對于理想平面的偏離，它描述了隨機(jī)路面下車輛的振動輸入。同時，相較于車輛發(fā)動機(jī)等激勵源，路面不平度是車輛振動系統(tǒng)的主要振源，在建立車輛系統(tǒng)整車數(shù)學(xué)或者仿真模型進(jìn)行平順性分析時，可以忽略發(fā)動機(jī)等激勵源，而將路面不平度激勵作為唯一的激勵源施加在車上。

路面不平度通常用來描述路面的起伏程度[1]。對路面譜的研究首先在于獲取路面譜，最直接的方法就是測量，為此，研究者們發(fā)明了多種路面譜測量設(shè)備。近年來，隨著傳感器技術(shù)、計算機(jī)技術(shù)和信號處理技術(shù)的飛速發(fā)展, 人們對路面不平度的采集、測量和各種試驗方法也在不斷的更新和改進(jìn)[1],這一領(lǐng)域已經(jīng)涌現(xiàn)出了多種測量和試驗分析的新方法。一般按照測量原理的不同可分為直接接觸式測量儀和非接觸式測量儀(響應(yīng)式測量儀)等。國內(nèi)直接接觸式測量儀的典型代表是1979年, 國內(nèi)長春汽車研究所的趙繼海等人發(fā)明的拖車式真實路形儀[1], 通過測量拖車上前后輪與拖臂等部件之間的角度變化來獲得路面的真實路形。

在對一定道路的測量和分析研究中，學(xué)者們發(fā)現(xiàn)路面不平度雖然不能夠用普通數(shù)學(xué)函數(shù)描述，但是其具有隨機(jī)、平穩(wěn)、各態(tài)歷經(jīng)的特征[1],這些特征在統(tǒng)計學(xué)意義上具有完整的理論描述方法，研究發(fā)現(xiàn)路面不平度可以用平穩(wěn)隨機(jī)過程理論來分析描述。一種普遍的做法是選取與道路垂直縱斷面與道路表面的交線作為路面不平度的樣本，取與道路垂直縱斷面相垂直的某一個道路平面作為理想平面，以該理想平面為基準(zhǔn)，路面各點具有高程值，通過樣本的數(shù)字特征方差和功率譜密度函數(shù)可以描述路面的不同特征。選取的理想平面應(yīng)該使樣本的均值為0，這樣樣本的方差可以描述路面平度大小的總體情況；功率譜密度函數(shù)表達(dá)的是路面平度能量在空間頻域的分布，它描述了路面不平度的基本形式和路面總體特征。

在進(jìn)行車輛的平順性等特性仿真分析時，還應(yīng)該考慮車輛前、后輪的相關(guān)性以及左、右輪的相關(guān)性。很多研究嘗試用各種不同的數(shù)學(xué)方法來模擬路面不平度，典型的有ARMA模型方法[2]、諧波疊加法(或稱三角級數(shù)法)[3-4]、積分單位白噪聲法以及濾波器整形白噪聲法[5]等。本研究選用諧波疊加法模擬路面的激勵。

1 諧波疊加法

根據(jù)隨機(jī)路面的統(tǒng)計特征，路面譜為平穩(wěn)、遍歷的均值為零的Gaussian過程，因此可用三角級數(shù)進(jìn)行模擬。諧波疊加法采用以離散譜逼近目標(biāo)隨機(jī)過程的模型，是一種離散化數(shù)值模擬路面的方法，本文將采用正弦波三角級數(shù)對隨機(jī)路面進(jìn)行模擬。

路面不平度通常采用功率譜密度函數(shù)來描述其高程值特征。國際標(biāo)準(zhǔn)協(xié)會文件ISO/TC108/SC2N67制定了路面不平度的功率譜密度表達(dá)式模型和分級方法，我國國家標(biāo)準(zhǔn)也沿用了這種方法。該方法規(guī)定在一段有意義的空間頻率[n1，n2]范圍內(nèi)，根據(jù)路面的不同等級情況給定了路面位移譜密度Gq(n)。Gq(n)可以根據(jù)標(biāo)準(zhǔn)或試驗取不同的具體等級對應(yīng)值。路面不平度功率譜密度擬合表達(dá)式為

(1)

式中：n0為參考空間頻率，一般取為n0=0.1 m-1；Gq(n0)為空間頻率為n0時的路面功率譜密度，也叫路面不平度系數(shù)，與路面等級相對應(yīng)；W為頻率指數(shù)，決定了路面功率譜密度頻率結(jié)構(gòu)，一般分級路面譜的頻率指數(shù)取為W=2；n為空間頻率有效頻帶中的某一空間頻率，其帶寬為[n1，n2]。

n1和n2規(guī)定了有效頻帶的上限和下限，帶寬的確定應(yīng)保證汽車在以平均速度行駛時，路面不平度引起的振動要包括汽車振動的主要固有頻率。

(2)

對于這個積分運算，可以采用離散近似法，將空間頻率范圍[n1，n2]劃分成m個小區(qū)間Δni，當(dāng)劃分區(qū)間個數(shù)m較大時，對于每個小區(qū)間的積分值可用式Gq(nmid,i)*Δni表示，Gq(nmid,i)每個小區(qū)間的中心頻率nmid,i(i=1,2,…,m)處對應(yīng)的功率譜密度值。這樣，可將式(2)改寫成

(3)

(4)

將各個小區(qū)間對應(yīng)的正弦波函數(shù)進(jìn)行疊加，就得到垂直縱剖面時域路面隨機(jī)位移。

(5)

式中：x為時域路面縱向位置；θi為均勻分布在[0，2π]之間的隨機(jī)數(shù)。

國標(biāo)GB7031—87等級路面標(biāo)準(zhǔn)[7-9]根據(jù)路面功率譜密度將路面分為8級，規(guī)定了各級路面不平度系數(shù)Gq(n0)的幾何平均值和各級路面均方根值σq，如表1所示。

表1 GB7031路面等級劃分

根據(jù)式(5)以及表1提供的參數(shù)，取空間頻率范圍為[0.011，2.83]m-1，將空間頻率范圍分成200份，即令m=200，通過編制Matlab程序，就可以得到隨機(jī)路面各激勵點高程值的數(shù)值和分布。圖1，2是B級路面、道路長度為400 m的路面不平度及其功率譜密度曲線圖。

由圖2可知，擬合的隨機(jī)路面滿足路面功率譜密度函數(shù)的要求。計算B級路面不平度均方根值7.385×10-3m與標(biāo)準(zhǔn)路面不平度均方根值的相對誤差為2.95%，在可接受的誤差范圍內(nèi)，故可認(rèn)為按照諧波疊加法建立的二維隨機(jī)路面滿足標(biāo)準(zhǔn)路面譜要求，隨機(jī)路面可以滿足單車轍路面激勵計算的要求。

2 三維隨機(jī)路面的Matlab程序?qū)崿F(xiàn)

式(5)給出了隨機(jī)路面沿道路長度方向高程值的分布，并通過該方法建立了滿足標(biāo)準(zhǔn)路面譜要求的隨機(jī)路面，但是，這是一個二維路面分布，它沒有表達(dá)出橫向剖面的高程值分布情況。要建立三維的隨機(jī)路面，應(yīng)該對表達(dá)式(5)進(jìn)行擴(kuò)展，將其路面高程值的分布擴(kuò)展到橫向剖面。

對空間頻率n的進(jìn)行考究可知，空間頻率表示單位距離內(nèi)包含的單一波長的個數(shù)。參考文獻(xiàn)[10-12]2種擴(kuò)展方法，結(jié)合實際仿真效果，得出本文三維空間路面擴(kuò)展算法為

(6)

式中：x為隨機(jī)路面縱向位置；y為隨機(jī)路面橫向位置；θi(y)為路面上任意行y屬于[0，2π]之間的隨機(jī)數(shù)。

這種擴(kuò)展算法基于下面的考慮：隨機(jī)路面的任意一個縱向剖面都應(yīng)該滿足路面譜功率譜密度函數(shù)要求，并且任意2個縱向剖面之間不具有相關(guān)性，而對橫向剖面沒有這樣的要求；隨機(jī)路面應(yīng)具有隨機(jī)性，不能出現(xiàn)確定性函數(shù)關(guān)系形式的路面。

在進(jìn)行整車平順性仿真分析時，關(guān)注的是整車行進(jìn)方向的路面不平度分布情況以及橫向車輪之間的路面不平度分布是否存在相關(guān)關(guān)系，故需要在計算單一縱向剖面的時候，將θi限定為一組m個相同的隨機(jī)數(shù)，這樣能保證這一組縱向剖面不平度數(shù)據(jù)滿足功率譜密度要求；而為了使不同縱向剖面的路面不平度分布不具有相關(guān)關(guān)系，需要將橫向參數(shù)y加入到公式中并且各組隨機(jī)數(shù)取值不能一樣，這樣在公式中加入的隨機(jī)數(shù)用θi(y)來表達(dá)比較合適，表示路面上隨意行y屬于[0，2π]之間的隨機(jī)數(shù)。

根據(jù)擴(kuò)展算法式(6)，利用Matlab編制程序，計算得到了一定條件下的三維空間路面不平度的空間分布。如圖3給出了8 m×400 m B級三維隨機(jī)路面不平度分布情況。

截取不同縱向剖面數(shù)據(jù)繪制路面不平度曲線。圖4，5為y=0.7 m處縱向路面不平度及其功率譜密度曲線，圖6，7為y=6.5 m縱向路面不平度及其功率譜密度曲線。

計算2條縱向路面不平度均方根值，分別為7.374 3×10-3m，7.356 4×10-3m，與標(biāo)準(zhǔn)B級路面均方根值7.61×10-3m的相對誤差為-3%和-3.3%，這個相對誤差在可接受的范圍內(nèi)，擬合路面滿足均方根值要求；由圖5，7可知，路面不平度與標(biāo)準(zhǔn)路面功率譜密度函數(shù)擬合率較高，縱向路面不平度滿足功率譜密度函數(shù)要求。

3 ADAMS三維隨機(jī)路面模型構(gòu)建

進(jìn)行整車級的ADAMS動力學(xué)仿真時，通常要用到路面文件，ADAMS提供的路面比較簡單，很難滿足車輛性能仿真的要求。因此，建立三維路面文件的通用模型，在各種工況下適配車輛仿真的條件，滿足各種仿真試驗的要求，顯得尤為重要。

ADAMS中自身配備了一些路面文件，它們以.rdf文件的形式存在，他們可以滿足簡單仿真試驗的要求。另外，他們也給出了ADAMS中路面文件編制的格式和要求，從而為進(jìn)行通用性路面模型的建立提供了條件。采用三角網(wǎng)絡(luò)法[13]構(gòu)建的三維路面文件，能夠通用于ADAMS/view和ADAMS/car中，滿足仿真的要求。

3.1 路面模型的一般組成

ADAMS 3D路面文件一般包含6個模塊，分別是表頭(MDI HEADER)、單元(units)、定義方式(definition)、偏移量(offset)、節(jié)點(nodes)、單元(elements)[14-15]。其中，表頭一般包含路面文件的一些類型聲明；如果沒有特別要求，偏移量模塊可以省略。建立隨機(jī)路面模型，需要重點關(guān)注的是路面文件中節(jié)點和單元模塊。節(jié)點模塊中每一個節(jié)點都是一個四維向量矩陣，分別包含節(jié)點序號和節(jié)點的x，y，z方向上的坐標(biāo)。單元模塊的每一個單元都是三角形單元，它是由一個五維向量矩陣，由3個節(jié)點序號和節(jié)點動、靜摩擦系數(shù)組成；一般動、靜摩擦系數(shù)也可以合并為一個，統(tǒng)稱為摩擦系數(shù)，這樣單元也可以是一個四維向量矩陣。圖8為節(jié)點和單元模塊的編寫格式。圖9為節(jié)點和單元的組成。只要確定節(jié)點(nodes)和單元(elements)，按照上述規(guī)則輸入到路面文件中，就可以得到通用的路面文件。

3.2 通用路面模型的編制

為此，分別建立一個四維向量節(jié)點矩陣NodeNum_N×4和五維向量單元矩陣ElementNum_E×5。對于節(jié)點矩陣NodeNum_N×4，由于節(jié)點的排序是按照“先縱后橫，由小到大”的原則，因此，可以考慮按照列向量分別編寫填充節(jié)點矩陣NodeNum_N×4。

第1列為節(jié)點序號，節(jié)點序號標(biāo)識節(jié)點的排序和個數(shù)，直接為由小到大的自然數(shù)編號，由于已經(jīng)確定路面有Num_N個節(jié)點，故可以將(1,Num_N)的順序自然數(shù)直接寫入NodeNum_N×4的第1列；

第2列縱向坐標(biāo)列陣為路面節(jié)點的x方向坐標(biāo)值，把縱向某一行節(jié)點的x方向坐標(biāo)值依次排列并循環(huán)Num_y次填入NodeNum_N×4的第2列；

第3列橫向坐標(biāo)列陣為路面節(jié)點的y方向坐標(biāo)值，把縱向某一行節(jié)點的y方向坐標(biāo)值依次排列并循環(huán)Num_y次填入NodeNum_N×4的第3列；

第4列豎直坐標(biāo)列陣為路面節(jié)點的z方向坐標(biāo)值，把存放高程值的ZNum_x×Num_y矩陣?yán)背闪邢蛄刻钊爰纯桑?/p>

對于單元矩陣ElementNum_E×5，同樣考慮上述方法，按照列向量分別編寫填充。

得到路面的節(jié)點矩陣和單元矩陣之后，將其分別輸入到ADAMS中的.rdf替換相應(yīng)模塊數(shù)據(jù)，就形成了通用ADAMS隨機(jī)路面模型文件。

4 結(jié)束語

本文根據(jù)諧波疊加法生成了二維隨機(jī)路面譜模型，模型滿足標(biāo)準(zhǔn)路面譜功率譜密度曲線要求，將二維路面擴(kuò)展為三維隨機(jī)路面模型；研究了ADAMS中隨機(jī)路面文件編制規(guī)則，提出了編制三維隨機(jī)路面模型文件的方法，編制生成了ADAMS三維隨機(jī)路面模型文件；仿真過程表明，依據(jù)該方法得到了三維路面模型文件可用于ADAMS仿真研究計算中。

參考文獻(xiàn)：

[1] 段虎明,石峰,謝飛,等.路面不平度研究綜述[J].振動與沖擊,2009,28(9):95-101.

DUAN Hu-ming,SHI Feng,XIE Fei,et al.Road Roughness Study Summarize[J].Journal of Vibration and Shock,2009,28(9):95-101.

[2] 唐光武,賀學(xué)鋒,顏永福.路面不平度的數(shù)學(xué)模型及計算機(jī)模擬研究[J].中國公路學(xué)報,2000,13(1):115-117.

TANG Guang-wu,HE Xue-feng,YAN Yong-fu.Mathematical Model of Road Roughness in Time Domain and Its Simulation[J].China Journal of Highway and Transport, 2000,13(1):115-117.

[3] AU F T K,CHENG Y S,CHEUNG Y K.Effects of Random Road Surface Roughness and Long-Termd Eflection of Prestressed Concrete Girder and Cable-Stayed Bridges on Impact due to Moving Vehicles[J].Computers and Structures,2001(79):853-872.

[4] SCHIEHLEN W,HU B.Spectral Simulation and Shock Absorber Identification[J].International Journal of Non-Linear Mechanics,2003(38):161-171.

[5] 周云山,于秀敏.汽車電控系統(tǒng)理論與設(shè)計[M].北京:北京理工大學(xué)出版社,1999.

ZHOU Yun-shan,YU Xiu-min.The Theory and Design of Automobile Electronic Control System[M].Beijing:Beijing Institute of Technology Press,1999.

[6] 徐延海.隨機(jī)路面譜的計算機(jī)模擬[J].農(nóng)業(yè)機(jī)械學(xué)報,2007,38(1):33-36.

XU Yan-hai.Computer Simulation on Stochastic Road Irregularities[J].Transactions of The Chinese Society of Agricultural Machinery,2007,38(1):33-36.

[7] GB/T7031-1986.車輛振動輸入與路面不平度表示方法[S].北京:中國標(biāo)準(zhǔn)出版社,1987.

GB/T7031-1986.Vechicle Vibration Imput and Road Rougness Express Method[S].Beijing:China Standard Publishing Company,1987.

[8] 余志生.汽車?yán)碚揫M].4版.北京：機(jī)械工業(yè)出版社,2007.

YU ZHI-sheng.Automobile Theory[M].4th ed.Beijing:Machinery Industry Press,2007.

[9] 徐昭鑫.隨機(jī)振動[M].北京:國防工業(yè)出版社,1990.

XU Shao-xin.Random Vibration[M].Beijing:National Defense Industry Press,1990.

[10] 梁新成,張軍,徐瑞.三角級數(shù)法的路面重構(gòu)仿真研究[J].汽車工程學(xué)報,2011,1(5):442-447.

LIANG Xin-cheng,ZHANG Jun,XU Rui.Study of Road Reconstruction with Trigonometric Series Method[J].Chinese Journal of Automotive Engineering,2011,1(5):442-447.

[11] 吳參,王維銳,陳穎,等.三維路面譜的仿真與驗證[J].浙江大學(xué)學(xué)報:工學(xué)版,2009,43(10):1935-1938.

WU Can,WANG Wei-rui,CHEN Ying,et al.Simulation and Validation of Three Dimensional Road Surface Spectrum[J].Journal of Zhejiang University:Engineering Science ed,2009,43(10):1935-1938.

[12] 張永林,胡志剛,陳立平.時空相關(guān)車輛道路的高效數(shù)值仿真[J] .農(nóng)業(yè)機(jī)械報,2005,36(9):13-15.

ZHANG Yong-lin,HU Zhi-gang,CHEN Li-ping.A New Approach to Numerical Simulation of Vehicle Road Process with Temporal Spatial Correlation[J].Transactions of the Chinese Society of Agricultural Machinery,2005,36(9):13-15.

[13] 趙治.ADAMS/Car軟件中隨機(jī)路面建立方法[J].客車技術(shù)與研究,2012(1):11-13.

ZHAO Zhi.Method on Building Random Road in ADAMS/Car[J].Bus and Couch Technology Reserch,2012(1):11-13.

[14] 陳軍.MSC.ADAMS技術(shù)與工程分析實例[M].1版.北京：中國水利水電出版社,2008.

CHEN Jun.Examples of Technology and Engineering Analysis by MSC.ADAMS[M].1st ed.Beijing:China Water Conservancy and Hydropower Press,2008.

[15] 鄭建榮.ADAMS:虛擬樣機(jī)技術(shù)入門與提高[M].北京：機(jī)械工業(yè)出版社,2002.

ZHENG Jian-rong.ADAMS:Introduction and Improvement of Virtual Prototyping Technology[M].Beijing:Machinery Industry Press,2002.