段虎明，石 鋒，馬 穎，張開斌，謝 飛

(1.中國汽車工程研究院博士后科研工作站，重慶 400039;2.中國汽車工程研究院，重慶 400039;3.同濟大學(xué)汽車學(xué)院，上海 201804)

道路路面不平度數(shù)據(jù)是車輛工程試驗中研究外部輸入激勵的主要成分，影響車輛行駛的平順性、乘坐舒適性、操縱穩(wěn)定性、零部件疲勞壽命、運輸效率、油耗等各個方面［1－2］，因此針對道路路面平整度的準(zhǔn)確評價和獲取有代表性的路面特征參數(shù)就顯得非常重要。在汽車工程領(lǐng)域中，常常使用路面不平度曲線的功率譜密度來反映路面的質(zhì)量好壞，然而由于功率譜密度是一條曲線，在分析和評價不同路段之間差異時，使用大量的曲線來分析就顯得很困難，因而在處理這類數(shù)據(jù)時，需要進(jìn)一步的分析路面的特征，提取出若干可以反映路面平整度水平的特征參數(shù)，用于統(tǒng)計分析計算和比較。本文正是在路面不平度的功率譜密度曲線不便于評價道路路面這一問題的基礎(chǔ)上，提出了一種可以通過若干參數(shù)來簡潔評價路面質(zhì)量的方法，為海量路面數(shù)據(jù)統(tǒng)計分析研究提供基礎(chǔ)。

1 道路路面數(shù)據(jù)的測量方法介紹

道路路面的測量從上個世紀(jì)60年代就開始了，通過不斷的改進(jìn)和演變，形成一系列測量設(shè)備和系統(tǒng)，目前國內(nèi)外使用較多的是基于慣性基準(zhǔn)的道路路面數(shù)據(jù)測量系統(tǒng)［1－3］，下面以國家863計劃重點項目中開發(fā)的中國典型道路譜測量系統(tǒng)為例進(jìn)行介紹。

道路譜測量系統(tǒng)主要是通過激光位移傳感器和加速度傳感器等配合測量得到道路路面的不平度曲線信號，系統(tǒng)的構(gòu)成如圖1所示［2－6］。該系統(tǒng)主要有激光斷面儀(內(nèi)置激光位移傳感器、高頻和低頻加速度傳感器等用來測量多路激光信號和加速度信號，二者通過復(fù)雜的數(shù)學(xué)計算即可得到路面的高程曲線)、精密陀螺儀(測量車身姿態(tài)變化信息參數(shù))、GPS測量設(shè)備(實時測量道路的經(jīng)緯海拔高等信息)、高速攝像機(獲取試驗路段的視頻信息)、車速及距離傳感器(測量試驗路段的距離脈沖信息)以及一些輔助的測量設(shè)備(采集控制箱、工控機、UPS不間斷電源、存儲海量數(shù)據(jù)使用的磁盤陣列等)組成。該道路譜測量設(shè)備，不但測量速度快(按照正常的行駛車速測量，不需要交通管制)、而且測量精度較高，測量頻段范圍也寬，較好的彌補了現(xiàn)有測量設(shè)備的不足，得到了較高精度的道路路面的剖面曲線信號。

圖1 測量系統(tǒng)的構(gòu)成Fig.1 Instrument integration of test system vehicle

通過圖1中路面數(shù)據(jù)測量系統(tǒng)的測量和解算，就可以得到初步的路面不平度數(shù)據(jù)，再通過若干的預(yù)處理計算，如信號的毛刺剔除、信號趨勢項的提取、分段試驗數(shù)據(jù)的平滑過渡連接等［7－10］，就可以得到真實有效的道路路面不平度曲線。

2 基于功率譜密度的路面特征參數(shù)提取

在路面不平度數(shù)據(jù)的使用過程中，首先就是要評估路面的平整度水平，充分獲取路面的特征參數(shù)信息，然后才能在大量的路面數(shù)據(jù)庫中提取出需要的道路路段數(shù)據(jù)。常規(guī)的路面特征參數(shù)有:均值(反映路面不平度數(shù)據(jù)的平均水平和集中趨勢)、方差或標(biāo)準(zhǔn)差(反映路面不平度數(shù)據(jù)的發(fā)散程度和離散趨勢)等，除了這些常用的統(tǒng)計參數(shù)，在公路部門行業(yè)常用的路面平整度評價參數(shù)為國際平整度指數(shù)［11］，而在汽車工程領(lǐng)域，常用的路面評價方法是功率譜密度曲線，在國家標(biāo)準(zhǔn)中就有明確的闡述［12］，然而這種功率譜密度的表現(xiàn)形式是一條隨著空間頻率變化的譜密度曲線，而不是若干指標(biāo)數(shù)據(jù)，其局限性在于針對大量路段對比分析其差異時不便于有效進(jìn)行，因此需要通過一種合理的方法從中提煉出若干代表性的參數(shù)，以方便在大量路面數(shù)據(jù)的統(tǒng)計分析和對比研究中使用。

2.1 路面不平度的功率譜分析

由于路面不平度曲線是沿著道路路面前進(jìn)方向變化的路面垂直高度曲線，即它是以等距離間隔為橫軸，以不斷變化的路面垂直高度為縱軸的二維曲線(相當(dāng)于等距離間隔的路面縱向的剖面曲線)。這樣的曲線數(shù)據(jù)，對其計算功率譜密度，得到的是以空間頻率(單位為:1/m)為橫軸，以單位空間頻率范圍內(nèi)的有限均方值為縱軸的譜密度曲線，稱為路面垂直位移功率譜密度，簡稱道路譜或路譜。當(dāng)然路譜也可以使用速度功率譜密度表示。

路譜的計算過程與傳統(tǒng)以等時間間隔的數(shù)據(jù)序列類似，需要設(shè)定很多計算參數(shù)，在GB/T7031－2005中有比較詳細(xì)的說明［12］，比如分析方法一般采用FFT(快速傅里葉變換)方法，窗函數(shù)使用漢寧窗，修正系數(shù)采用8/3［13］，另外需要根據(jù)采集條件設(shè)定分析空間頻率的范圍、采樣頻率、抗混濾波等。例如使用北京北六環(huán)附近省道S321高麗營六元橋段數(shù)據(jù)為例，計算得到如圖2所示的功率譜密度曲線(雙對數(shù)坐標(biāo))。

圖2 功率譜密度曲線Fig.2 Power spectrum density of road profile

從圖2中可以看出雙對數(shù)坐標(biāo)下，由于是定帶寬分析方法來計算功率譜密度，因此在高頻區(qū)域就出現(xiàn)了非常豐富的頻率分量，容易給人產(chǎn)生錯覺［12］。因此在GB/T7031－2005中提出了使用多倍頻程的方法來平滑該曲線，主要分三個區(qū)間:

(1)倍頻程分析:從最低計算頻率到中心頻率為0.031 2 1/m的頻率帶寬;

(2)1/3倍頻程分析:從中心頻率為0.049 6 1/m到0.25 1/m的頻率帶寬;

(3)1/12倍頻程分析:從中心頻率為0.2806 1/m到最高計算頻率。

如表1所示為使用倍頻程平滑計算路譜的中心頻率與截止頻率，nc表示中心頻率，nl和nh分別表示上下截止頻率，1/3倍頻程與1/12倍頻程的中心頻率與截止頻率詳見參考文獻(xiàn)［12］中的表2。

表1 倍頻程帶寬的中心頻率與截止頻率Tab.1 Center and cut－off frequency of octave bandwidth

通過表1所示的多倍頻程區(qū)間，在給定的頻帶內(nèi)，平均功率譜密度按下式計算:

式(1)中Gs(i)為平滑帶寬i上的平滑功率譜密度，Be表示空間頻率分辨率，nH和nL分別為:

式(2)和式(3)中，INT表示取整數(shù)。

通過以上計算即可完成對原始功率譜密度曲線的平滑處理，針對圖2中的曲線使用式(1)進(jìn)行平滑濾波可得到如圖3所示的光滑曲線，該曲線就是工程中常用來評價和分析路面特性的功率譜密度曲線。如圖4所示為在平滑濾波計算時用到的路段描述參數(shù)和倍頻程分析參數(shù)。

2.2 路面分級與評價方法討論

在目前的汽車工程領(lǐng)域，使用路譜來分析和處理路面數(shù)據(jù)時，就是在圖3平滑濾波以后的功率譜密度曲線上進(jìn)行的。該功率譜密度曲線在雙對數(shù)坐標(biāo)下，可使用最小二乘方法擬合出一條直線，用來反映曲線的總體趨勢。擬合公式為:

其中的Gd(n)表示位移功率譜密度，n0表示參考空間頻率(n0=0.1 1/m)，w為擬合功率譜密度的指數(shù)，圖3中光滑曲線的擬合結(jié)果如圖5所示。

在GB/T7031－2005中，道路路面的八級分級圖正是以式(4)為基礎(chǔ)的，具體的分級方法是:設(shè)定w=2，然后以參考文獻(xiàn)［12］中的表C.2中的Gd(n0)的上下限分為A級到H級，共八級，如圖6所示為將功率譜密度曲線繪制在標(biāo)準(zhǔn)八級分級圖中的結(jié)果。

從圖6中可以看出，該路段基本上在A級、B級和C級中，路面質(zhì)量相對較好。但是該分級方法也不能定量描述路面等級情況和質(zhì)量好壞，因此只能籠統(tǒng)的描述和分析路面。

圖3 平滑濾波以后的功率譜密度曲線Fig.3 Power spectrum density after smoothing filter

圖4 平滑濾波時的倍頻程濾波參數(shù)Fig.4 Parameters of octave-band filter

圖5 平滑功率譜密度曲線的直線擬合Fig.5 Fitting a straight line of smoothing PSD

圖6 八級分級圖中的功率譜密度曲線Fig.6 PSD in regionalism map

2.3 道路路面的特征參數(shù)提取

在很多的工程項目中，需要對大量的道路路段進(jìn)行統(tǒng)計分析，綜合所有路面的總體特征，這種需求下，使用上面的常規(guī)功率譜密度分析方法就很難實現(xiàn)，因為每條試驗路面都會生成一條功率譜密度曲線。為了解決該問題，這里提出了一種路面等級比例分析的方法。該方法主要原理是:平滑以后的功率譜密度曲線與標(biāo)準(zhǔn)的八級分級的邊界線都有交點(在圖6中幾個小圓圈的標(biāo)注點)，以這些交點作為分界點，分別統(tǒng)計落在每一個分級區(qū)間的曲線點，統(tǒng)計其在所有數(shù)據(jù)點中所占的比例，這樣我們每條功率譜密度曲線就轉(zhuǎn)化為A級到H級的八個比例參數(shù)，這八個比例參數(shù)的和為100%，圖6中曲線的比例系數(shù)見表2所示。

表2 功率譜密度曲線在路面等級中的比例Tab.2 Percents of a PSD curve in regionalism map

從表2的統(tǒng)計結(jié)果也可看出，該試驗路段主要集中在A、B和C三個等級。通過這樣的處理，路面數(shù)據(jù)平滑以后的功率譜密度曲線就可以用這8個比例參數(shù)替代，用來評估和分析路面的等級高低及質(zhì)量好壞，方便了在后續(xù)的數(shù)據(jù)分析處理中，大量路段試驗數(shù)據(jù)的統(tǒng)計分析和對比研究。

如果在大量路面數(shù)據(jù)處理的統(tǒng)計分析中，需要進(jìn)一步壓縮統(tǒng)計參量，可以考慮使用平滑以后功率譜密度曲線的擬合直線來分析，兩個參數(shù):截距和斜率。截距為擬合直線在0.1 1/m處的截距值(與八級分級的截距位置相同)，其值的大小反映了路面的等級。斜率為擬合直線的斜率，反映道路路面一致性程度，其斜率越接近－2(八級分級的界限斜率)，表明其一致性越好，越偏離－2，表示一致性越差(即對應(yīng)分析路段所跨越和包含的等級越多)。圖5中的擬合直線其斜率和截距見表3所示。

表3 功率譜密度曲線的擬合直線參數(shù)Tab.3 Parameters of fitting a straight line of PSD

由表3可知，由于截距值僅為6.16×1e－6 m3，所以該路段主要集中在A級(A級與B級的分界線截距為32×1e－6 m3)，而斜率值－1.01與－2相差較大，表明該路段的一致性不好，會落在(跨越)多個分級區(qū)間。

通過以上的分析可以看出，一條試驗路段的不平度數(shù)據(jù)，可通過計算其功率譜密度曲線，進(jìn)而得到表征該路段特征的八個比例參數(shù)和兩個擬合直線參數(shù)，這兩種特征參數(shù)的提取方法各有優(yōu)劣:①比例參數(shù)提取方法參數(shù)較多，不利用大量路段的統(tǒng)計分析，但其描述和刻畫路段的特征比較清楚和具體;②擬合直線方法特征參數(shù)比較少，便于統(tǒng)計分析，但是對路段的描述比較籠統(tǒng)(僅使用截距反映路面等級，不準(zhǔn)確，誤差大)?？傮w來看，推薦使用比例參數(shù)方法提取路面的特征參數(shù)，這樣特征參數(shù)的代表性更準(zhǔn)確。

3 實例計算

通過以上的道路路面評價和特征參數(shù)提取方法介紹，我們已從一條試驗路面的采樣數(shù)據(jù)中獲得了表征路面特點的特征參數(shù)，下面隨機選擇幾條不同的試驗路段進(jìn)行實例分析，計算結(jié)果如表4所示。

通過表4的比例特征參數(shù)和擬合直線的特征參數(shù)列表都可以看出，不同等級的路面中，其基本趨勢隨著路面技術(shù)等級(或行政等級)的下降路面的質(zhì)量逐漸變差，這與實際情況比較相符。

從項目開發(fā)中得到的海量路段的試驗數(shù)據(jù)來看，大部分道路都集中在前面四個等級(A級到D級)，很少有路面可以達(dá)到F級以后，這也說明，這些年來隨著國家經(jīng)濟的快速發(fā)展，我國的道路路面整體變好，也說明A級到H級的八級分級方法比較陳舊，應(yīng)適時的修訂國家標(biāo)準(zhǔn)，舍棄后面用不到的道路等級，并細(xì)分前面的道路等級，用于更準(zhǔn)確的評價道路路面。

表4 幾條試驗路面的特征參數(shù)Tab.4 Parameters of some test road profile

4 結(jié)論

基于道路譜課題采集的海量路面試驗數(shù)據(jù)，提出了道路路面評價和特征參數(shù)提取的方法，主要通過對路面不平度曲線的空間域功率譜分析，從平滑的功率譜密度曲線中，提取出了可以表征道路路面特性的若干特征參數(shù)(比例系數(shù)參數(shù)與擬合直線參數(shù))，為大量路面數(shù)據(jù)的后續(xù)處理和分析計算提供基礎(chǔ)。

這兩種道路路面特征參數(shù)的提取方法使用簡單方便，參數(shù)意義明確，提取效果明顯，因此有較好的實用和推廣價值。目前這兩種路面特征參數(shù)提取方法已形成程序模塊，在中國典型道路譜數(shù)據(jù)庫［14－15］中海量路面數(shù)據(jù)處理中得到了很好的應(yīng)用。這些方法不但可用于路面數(shù)據(jù)統(tǒng)計分析，而且也可在其它行業(yè)的工程振動信號的統(tǒng)計分析中推廣應(yīng)用。

［1］段虎明，石 鋒，謝 飛，等.路面不平度研究綜述［J］.振動與沖擊，2009，28(09):95－101.

DUAN Hu-ming，SHI Feng，XIE Fei，et al.A Survey of road roughness study［J］.Journal of Vibration and Shock，2008，28(9):95－101.

［2］段虎明，石 鋒，謝 飛，等.道路譜測量檢測技術(shù)研究綜述［J］.電子測量與儀器學(xué)報，2010，24(1):72－79.

DUAN Hu-ming，SHI Feng，XIE Fei，et al.Review on research of road spectrum measurement［J］.Journal of Electronic Measurement and Instruments，2010，24(1):72－79.

［3］韓 毅，楊殿閣，連小珉.基于路形測量的車載道路譜檢測儀研究與實現(xiàn)［C］.2009中國汽車工程學(xué)會年會論文集，2009:770－774.

［4］段虎明，石 鋒，謝 飛，等.道路譜測量系統(tǒng)的誤差與精度分析［J］.汽車工程，2010，32(9):783－788.

DUAN Hu-ming，SHI Feng，XIE Fei，et al.Error and accuracy analysis on road spectra measurement system［J］.Automotive Engineering，2010，32(9):783－788.

［5］段虎明，石 鋒，趙永杰，等.道路路面的剖面曲線的測量研究與實踐［J］.振動與沖擊，2011，30(3):155－160.

DUAN Hu-ming， SHIFeng， ZHAO Yong-jie， etal.Research and practice of road surface profile measurement［J］.Journal of Vibration and Shock，2011，30(3):155－160.

［6］ Duan Huming，Ma Ying，Shi Feng，et al.Measurement and Application of Road Profile［C］.The 3rd International Symposium on Test Automation＆ Instrumentation，Xiamen，China，2010，05:119－124.

［7］ Duan Huming，Xie Fei，Zhang Kaibin，et al.Signal trend extraction of road surface profile measurement［C］.The 2nd InternationalConference on SignalProcessing System，Dalian，China，2010，07:V2－694–V2－698.

［8］馬 穎，段虎明，謝 飛，等.道路路面測量數(shù)據(jù)中奇異信號識別和修正方法［J］.電子測量與儀器學(xué)報，2010，24(9):878－884.

MA Ying，DUAN Hu-ming，XIE Fei，et al.A method of signal singularity identification and correction in road surface measurement data［J］.Journal of Electronic Measurement and Instruments，2010，24(9):878 －884.

［9］ Duan H M，Ma Y，Xie Fei，et al.Smooth connection method of segment test data in road surface profile measurement［C］. 2010 InternationalConference on Software and Computing Technology， Kunming， China， 2010， 10，18－19.

［10］段虎明，謝 飛，張開斌，等.海量道路路面測量數(shù)據(jù)的若干預(yù)處理方法研究［J］.振動與沖擊，2011，30(8):101－106.

DUAN Hu-ming，XIE Fei，ZHANG Kai-bin，et al.Methods of signal pre-processing with massive road surface measurement data［J］.Journal of Vibration and Shock，2011，30(8):101－106.

［11］ Sayers M W，et al.Guidelines for conducting and calibrating road roughness measurements［R］.World Bank Technical Paper Number 46，1986.

［12］GB/T 7031－2005.機械振動——道路路面譜測量數(shù)據(jù)報告［S］.2005.

［13］柴立群，徐建程，許 喬.加窗后波前功率譜密度的計算值修正［J］.強激光與粒子束，2005，17(12):1835－1838.

CHAI Li-qun，XU Jian-cheng，XU Qiao.Modification of power spectrum density of wavefront after applying window［J］.High Power Laser and Particle Beams，2005，17(12):1835－1838.

［14］謝 飛，段虎明，馬 穎，等.中國典型道路譜數(shù)據(jù)庫的開發(fā)與應(yīng)用［J］.振動與沖擊，2010，29(11):160－164.

XIE Fei，DUAN Hu-ming，MA Ying，et al.Development and application of typical road spectrum database in China［J］.Journal of Vibration and Shock，2010，29(11):148－152.

［15］段虎明，石 鋒，張開斌，等.中國典型汽車道路譜研究的構(gòu)建及應(yīng)用［J］.中國測試，2009，35(4):76－79.

DUAN Hu-ming，SHI Feng，ZHANG Kai-bin，et al.Review on construction and application of Chinese typical vehicle road spectra［J］.China Measurement and Test，2009，35(4):76－79.