, ,

(1.海口市市政工程設(shè)計(jì)研究院,海南 ?？?570100；2.石家莊鐵道大學(xué) 土木工程學(xué)院,河北 石家莊 050043)

0 引言

公路工程主要是人-車(chē)-路-環(huán)境所構(gòu)成的非常復(fù)雜的系統(tǒng)工程，行車(chē)安全性及乘車(chē)舒適性等是由人-車(chē)-路-環(huán)境的相互作用所決定的。公路線(xiàn)形一經(jīng)確定就不能隨便更改，而線(xiàn)形的優(yōu)劣是公路設(shè)計(jì)及其作用的主要評(píng)價(jià)方法。在公路的運(yùn)營(yíng)期間，由于設(shè)計(jì)的缺陷和疏漏，往往造成公路幾何特性與車(chē)輛行駛特性、駕駛員特性三者之間的不匹配，甚至出現(xiàn)事故多發(fā)路段(交通黑點(diǎn))。據(jù)統(tǒng)計(jì)，2012年9月30日至10月7日長(zhǎng)假期間，全國(guó)共發(fā)生道路交通事故68 422起，同比下降24.1%。涉及人員傷亡的道路交通事故2 164起，造成794人死亡、2 473人受傷，直接財(cái)產(chǎn)損失1 325萬(wàn)元。由可見(jiàn)此，對(duì)于公路線(xiàn)形的檢驗(yàn)和評(píng)價(jià)應(yīng)該首先著手于設(shè)計(jì)階段，在設(shè)計(jì)完成以后，利用線(xiàn)形檢驗(yàn)和評(píng)價(jià)方法，發(fā)現(xiàn)它的不合理之處進(jìn)行及時(shí)改正，使得公路線(xiàn)形的投入和產(chǎn)出達(dá)到協(xié)調(diào)與統(tǒng)一。

鑒于汽車(chē)工程研發(fā)的車(chē)輛動(dòng)力學(xué)仿真技術(shù)，對(duì)于機(jī)械工業(yè)中產(chǎn)品的“設(shè)計(jì)—樣品—試驗(yàn)—修改設(shè)計(jì)”設(shè)計(jì)理念逐漸地受到重視，進(jìn)行類(lèi)似于工業(yè)產(chǎn)品那樣的虛擬行駛試驗(yàn)，來(lái)評(píng)價(jià)道路線(xiàn)形組合設(shè)計(jì)的質(zhì)量，而其發(fā)展則是一個(gè)緩慢的過(guò)程。Alan Nicholson[1]于1998年研究了公路平豎組合線(xiàn)形的一致性與汽車(chē)運(yùn)行速度之間的關(guān)系；2007年李磊等[2]利用軟件ADAMS/Car建立了汽車(chē)運(yùn)動(dòng)學(xué)模型和動(dòng)力學(xué)模型，并進(jìn)行動(dòng)力學(xué)仿真分析，能較全面地評(píng)價(jià)和預(yù)測(cè)動(dòng)力學(xué)性能；2010年朱亮亮等[3]采用ADAMS軟件建立實(shí)體仿真模型，將理論計(jì)算和仿真計(jì)算的內(nèi)容分析整理，為以后的公路設(shè)計(jì)與研究提供借鑒；2011年游克思等[4]利用MATLAB/Simulink建立了3自由度整車(chē)動(dòng)力學(xué)模型，通過(guò)數(shù)值仿真模擬，分析了汽車(chē)在不同超高、縱坡坡度以及合成縱坡下的動(dòng)力響應(yīng)，在此基礎(chǔ)上進(jìn)一步分析了在不同車(chē)速以及前輪轉(zhuǎn)角輸入下，道路幾何參數(shù)對(duì)車(chē)輛操縱穩(wěn)定性的影響。

1 UM軟件簡(jiǎn)介

Universal Mechanism軟件由俄羅斯布良斯克國(guó)立理工大學(xué)(Bryansk State Technical University)的著名數(shù)學(xué)家、計(jì)算力學(xué)專(zhuān)家Dmitry Pogorelov教授帶領(lǐng)的計(jì)算力學(xué)實(shí)驗(yàn)室開(kāi)發(fā)，簡(jiǎn)稱(chēng)UM。它是一款非常優(yōu)秀的通用多體系統(tǒng)動(dòng)力學(xué)仿真分析軟件，并具有輪式車(chē)輛、履帶車(chē)輛、機(jī)車(chē)車(chē)輛、列車(chē)、輪軌磨耗和車(chē)橋耦合等諸多專(zhuān)用分析模塊，已廣泛應(yīng)用于相關(guān)行業(yè)。

UM軟件基于現(xiàn)代計(jì)算機(jī)輔助工程技術(shù)，從建模到仿真和后處理融合了許多優(yōu)秀的算法和程序，從各個(gè)方面增強(qiáng)軟件的實(shí)用性。UM使用簡(jiǎn)單方便。用戶(hù)圖形界面友好，設(shè)計(jì)非常人性化，將多體系統(tǒng)的數(shù)學(xué)模型都隱藏在程序內(nèi)核，便于用戶(hù)操作。同時(shí)，UM也是一個(gè)開(kāi)放的系統(tǒng)，其良好的開(kāi)放性和兼容性，便于與其他軟件(如ANSYS、ADAMS/CAR等)聯(lián)合使用，具有一套自己的程序風(fēng)格和標(biāo)準(zhǔn)語(yǔ)言，可以根據(jù)相關(guān)語(yǔ)法就能做進(jìn)一步的開(kāi)發(fā)。

2 車(chē)輛-道路動(dòng)力學(xué)模型

2.1 汽車(chē)動(dòng)力學(xué)模型

UM中包含有一個(gè)特殊模塊UM Automotive用于車(chē)輛動(dòng)態(tài)分析。這個(gè)模塊包括集成到程序核心的附加工具，還有特有的懸架和傳動(dòng)裝置庫(kù)(獨(dú)立傳遞)。UM包括以下主要組件：道路幾何形狀生成與可視化；道路不平順文件和功率譜密度文件庫(kù)；道路微觀外形(不平順)生成與可視化工具；精確的輪胎力模型(輪胎與道路接觸力)；驅(qū)動(dòng)模型；特有的動(dòng)態(tài)實(shí)驗(yàn)裝置。UM Automotive模塊提供了兩種駕駛員模型：MacAdam’s model和Second order preview model。輪胎參數(shù)在UM Simulation里設(shè)置，利用輪胎向?qū)Чぞ呖梢钥焖賱?chuàng)建所需模型，用于表達(dá)模型參數(shù)存儲(chǔ)文件為(*.tr)。UM Automotive模塊中提供了兩種車(chē)型：小客車(chē)與貨車(chē)。本文主要研究的是公路平縱組合線(xiàn)形對(duì)行車(chē)動(dòng)力影響，通過(guò)比較確定哪種線(xiàn)形是好的，從而能夠滿(mǎn)足人們?nèi)找嬖鲩L(zhǎng)的舒適性和安全性要求。因此，直接選用UM中的車(chē)輛模型進(jìn)行車(chē)輛的動(dòng)態(tài)仿真，并在必要時(shí)對(duì)模型進(jìn)行了參數(shù)修改，而把重點(diǎn)放在了道路的平、縱、橫組合設(shè)計(jì)過(guò)程中。車(chē)輛動(dòng)力學(xué)模型主要有客車(chē)與貨車(chē)，在仿真分析時(shí)主要采用vaz21_09模型，模型具體參數(shù)見(jiàn)表1。

表1 汽車(chē)vaz21_09模型主要參數(shù)參數(shù)數(shù)值參數(shù)數(shù)值車(chē)體質(zhì)量/t1.269車(chē)輪搖頭轉(zhuǎn)動(dòng)慣量/(kg·m2)2.5×10-1車(chē)體側(cè)滾轉(zhuǎn)動(dòng)慣量/(kg·m2)4.0×102車(chē)輪橫向阻尼/(kN·s·m-1)6.0×101車(chē)體點(diǎn)頭轉(zhuǎn)動(dòng)慣量/(kg·m2)1.2×103車(chē)輪縱向阻尼/(kN·s·m-1)1.0×101車(chē)體搖頭轉(zhuǎn)動(dòng)慣量/(kg·m2)1.2×103車(chē)輪垂向阻尼/(kN·s·m-1)3.0×101車(chē)體材料密度/(kg·m-3)7.8×103輪胎內(nèi)寬/m1.5×10-1輪胎內(nèi)直徑/m5.6×10-1前輪輪距/m1.39輪胎截面高度/m1.2×10-1后輪輪距/m1.37輪胎截面寬度/m1.742×10-1軸距/m2.468車(chē)輪質(zhì)量(除去轉(zhuǎn)臂)/kg2.0×101車(chē)輪橫向剛度/(kN·m-1)1.0×102車(chē)輪側(cè)滾轉(zhuǎn)動(dòng)慣量/(kg·m2)2.5×10-1車(chē)輪縱向剛度/(kN·m-1)1.0×102車(chē)輪點(diǎn)頭轉(zhuǎn)動(dòng)慣量/(kg·m2)5.0×10-1車(chē)輪垂向剛度/(kN·m-1)6.25×102

2.2 道路模型

UM Automotive模塊中的路面工具包括兩部分：Macrogeometry(用于宏觀線(xiàn)路設(shè)計(jì))和Irregularities(生成路面不平度)。順著車(chē)輛的運(yùn)行方向，即為道路縱斷面的路面不平度，道路外形主要分為三種：宏觀(長(zhǎng)波)、微觀(短波)和粗糙度(粗糙紋理)。本文道路模型的建立主要是運(yùn)用緯地軟件(道路三維動(dòng)態(tài)可視化幾何設(shè)計(jì)軟件)，并從運(yùn)動(dòng)學(xué)角度出發(fā)，結(jié)合影響平豎曲線(xiàn)組合的主要因素，考慮汽車(chē)速度變化的影響，建立相關(guān)的3D道路線(xiàn)形設(shè)計(jì)文件，通過(guò)離散平、縱曲線(xiàn)，并將離散點(diǎn)坐標(biāo)轉(zhuǎn)化為UM軟件道路模型。

3 仿真計(jì)算

3.1 道路平縱橫斷面設(shè)計(jì)資料及參數(shù)

本汽車(chē)動(dòng)力仿真試驗(yàn)采用的是?？谑醒蛏铰肪W(wǎng)一期工程項(xiàng)目(主干道北段)。道路起點(diǎn)為?？谑欣@城高速公路博片村路口處，本設(shè)計(jì)起點(diǎn)坐標(biāo)為：X=204 627.916，Y=197 947.960(坐標(biāo)系統(tǒng)為海南?？讵?dú)立坐標(biāo)系統(tǒng))，途徑羊山地區(qū)本年度國(guó)際高爾夫球賽事項(xiàng)目區(qū)，終點(diǎn)與規(guī)劃路在永興鎮(zhèn)段內(nèi)銜接，本次設(shè)計(jì)終點(diǎn)坐標(biāo)為(X=199 895.948，Y=192 877.071)， 全長(zhǎng)為9 405.025 m。在本設(shè)計(jì)中，道路紅線(xiàn)寬度為30 m，按雙向6車(chē)道設(shè)置。道路等級(jí)為城市主干道Ⅰ級(jí)，設(shè)計(jì)車(chē)速為60 km/h，最大縱坡為3.253%，最小縱坡為0.5%。為了在UM中建模的方便，將設(shè)計(jì)路段進(jìn)行離散化，并選取其中特殊的一段(此段同時(shí)含有平曲線(xiàn)、緩和曲線(xiàn)、豎曲線(xiàn)以及超高的設(shè)置)進(jìn)行汽車(chē)動(dòng)力學(xué)仿真計(jì)算，具體用樁號(hào)表示為：K1+834.797～K4+534.797。將設(shè)計(jì)路段的平面設(shè)計(jì)資料利用緯地道路CAD系統(tǒng)軟件進(jìn)行重新設(shè)計(jì)，輸出線(xiàn)形的逐樁坐標(biāo)表，并將其用二維直角坐標(biāo)系繪制出來(lái)，具體如圖1、圖2所示。

圖1 仿真路段平曲線(xiàn) 圖2 仿真路段豎曲線(xiàn)

由于本設(shè)計(jì)資料中的超高過(guò)渡是按三次拋物線(xiàn)漸變的，通過(guò)設(shè)計(jì)橫斷面資料可得汽車(chē)左右輪的超高，并對(duì)離散的點(diǎn)進(jìn)行三次拋物線(xiàn)擬合，其得到的數(shù)據(jù)經(jīng)過(guò)軟件自動(dòng)修正后如圖3所示。然后通過(guò)點(diǎn)生成譜的方式，并且令諧波數(shù)分別等于10 000隨機(jī)生成路段的左右輪不平順，將超高與不平順疊加導(dǎo)入U(xiǎn)M中，經(jīng)過(guò)UM自動(dòng)修正后，生成仿真路段超高不平順如圖4所示。

圖3 仿真路段超高修正圖 圖4 仿真路段超高不平順

對(duì)于汽車(chē)仿真時(shí)的最大行車(chē)速度，取設(shè)計(jì)速度為60 km/h(16.7 m/s)。為了使仿真結(jié)果接近實(shí)際情況，通過(guò)統(tǒng)計(jì)身邊開(kāi)車(chē)駕駛?cè)嗽谵D(zhuǎn)彎時(shí)的駕駛速度，合理確定該路段的行車(chē)速度，從而符合實(shí)際駕駛操作。另外，為了減小軟件對(duì)超高修正的影響，起始路段采用勻加速模型，而曲線(xiàn)路段的車(chē)速適當(dāng)降低為36 km/h(10 m/s)。根據(jù)車(chē)速變化計(jì)算出各段行駛時(shí)間得到時(shí)間與速度關(guān)系如圖5所示。

圖5 汽車(chē)行駛速度圖

3.2 動(dòng)態(tài)仿真分析

由于實(shí)際工程中道路同側(cè)內(nèi)采用的是單向橫坡形式，因此仿真時(shí)將道路模型簡(jiǎn)化為單車(chē)道形式。通過(guò)對(duì)上述數(shù)據(jù)及參數(shù)的設(shè)置后，汽車(chē)動(dòng)態(tài)仿真時(shí)采用MacAdam駕駛員模型，參數(shù)Preview time取為1 s，Reaction time delay取為0.1 s，Neuromuscular lag取為0.1 s，Number of control steps取為4。首先需要得出橫向合成加速度的數(shù)據(jù)，然后利用2D Curve editor求導(dǎo)解出橫向加速度變化率的數(shù)據(jù)，橫向合成加速度值為其空間向量的模，即X、Y、Z軸加速度投影的模的平方和的二次方根。而其指向可以通過(guò)Animation Window實(shí)時(shí)監(jiān)控，非常方便。

汽車(chē)某時(shí)刻的行駛狀態(tài)，經(jīng)過(guò)一系列的汽車(chē)動(dòng)力學(xué)仿真，首先得到路段的各加速度在坐標(biāo)軸上的投影變化圖以及x軸、y軸和z軸加速度合成圖，即模m方向的加速度圖，如圖6所示。

圖6 路段仿真分析結(jié)果數(shù)據(jù)圖(忽略不平順)

通過(guò)分析圖6可以得出以下結(jié)論：平縱橫合成加速度值除了初始階段(主要是由于UM軟件對(duì)超高的自動(dòng)修正，此段無(wú)實(shí)際意義可以忽略)，其余段均小于0.6 m/s2，采用文獻(xiàn)[5]的橫向力加速度標(biāo)準(zhǔn)來(lái)評(píng)價(jià)舒適性，發(fā)現(xiàn)道路設(shè)計(jì)成果的動(dòng)態(tài)性能比較好。路段內(nèi)，平縱橫合成加速度出現(xiàn)拉升的位置基本上都是汽車(chē)加減速與勻速突變的位置。由于車(chē)輛運(yùn)行時(shí)與x、y軸成一定的角度，因此，從圖6中可以發(fā)現(xiàn)在直線(xiàn)路段上，x軸向加速度和y軸向的加速度要受到路拱坡度的影響，產(chǎn)生的加速度分別向各軸進(jìn)行投影，在曲線(xiàn)段內(nèi)，這種影響尤為明顯。此外，將上述圖中兩個(gè)路段的平縱橫合成加速度對(duì)時(shí)間求導(dǎo)，得到了合成加速度變化率圖形，具體如圖7所示。

圖7 仿真路段合加速度變化率(忽略不平順)

從圖7中可以看出，在加減速路段與勻速路段交接的位置加速度變化率很大，影響了乘車(chē)的舒適性。圖7中雖然將超高修正段的數(shù)據(jù)設(shè)為0，但是由于車(chē)輛行駛過(guò)程中，相鄰的曲線(xiàn)段之間存在耦合的特點(diǎn)，受修正段的影響，其相鄰段內(nèi)的加速度變化率也是比較大的。而在路段尾端由圓曲線(xiàn)進(jìn)入緩和曲線(xiàn)、由緩和曲線(xiàn)進(jìn)入直線(xiàn)時(shí)，加速度變化率雖然比較小，但是波動(dòng)還是比較密集的。另外，由于路段一在中間兩段緩和曲線(xiàn)之間沒(méi)有加入直線(xiàn)段，從圖7中可以發(fā)現(xiàn)拐點(diǎn)之后的加速度變化率波動(dòng)很大，直接影響了乘車(chē)的舒適性。

以上是在沒(méi)有加入不平順的情況下分析得到的結(jié)果，而實(shí)際情況，道路是存在不平順的，因此，在動(dòng)力仿真時(shí)加入一定的不平順更符合實(shí)際情況，具體仿真結(jié)果如圖8所示。

圖8 仿真路段合成加速度變化(考慮不平順)

從圖8中可以看出，在超高加入不平順之后，各軸向加速度都出現(xiàn)了不同程度的波動(dòng)，尤其是豎向加速度波動(dòng)最為明顯，而縱橫軸向加速度雖然有一定的波動(dòng)，但是基本上都是在未加不平順的加速度曲線(xiàn)上下一定范圍內(nèi)波動(dòng)，其中在曲線(xiàn)路段，不平順使得各加速度的波動(dòng)顯著增大。所以，實(shí)際公路運(yùn)營(yíng)期間，當(dāng)曲線(xiàn)路段出現(xiàn)問(wèn)題時(shí)，要及時(shí)進(jìn)行修整，否則會(huì)造成嚴(yán)重的不良后果。此外，汽車(chē)行駛時(shí)還可以對(duì)其輪胎的受力與滑動(dòng)進(jìn)行監(jiān)測(cè)，對(duì)公路線(xiàn)形存在不合理的位置進(jìn)行修改，減小車(chē)輛的損耗，進(jìn)而提高汽車(chē)行駛的安全性和穩(wěn)定性。

4 結(jié)語(yǔ)

將公路線(xiàn)形設(shè)計(jì)軟件與汽車(chē)動(dòng)力仿真分析軟件進(jìn)行結(jié)合，應(yīng)用緯地道路CAD系統(tǒng)軟件將前人進(jìn)行研究應(yīng)用的公路資料，重新設(shè)計(jì)，并得到平縱橫線(xiàn)形數(shù)據(jù)資料，將其編輯成UM軟件中仿真時(shí)所需要的文件，完成道路模型的建立。從2010年中國(guó)市政工程中南設(shè)計(jì)研究院設(shè)計(jì)的?？谑醒蛏铰肪W(wǎng)一期工程項(xiàng)目(主干道北段)設(shè)計(jì)資料中，選取了一段線(xiàn)形特殊的道路，進(jìn)行仿真分析，用平縱橫合成加速度及其變化率為評(píng)價(jià)指標(biāo)進(jìn)行了相關(guān)的分析。在仿真分析時(shí)考慮了路面的不平順，利用UM軟件隨機(jī)生成了不平順，并與超高相疊加得到左右輪不平順文件。在汽車(chē)動(dòng)力學(xué)仿真時(shí)，加入不平順來(lái)研究組合線(xiàn)形的舒適性問(wèn)題。由于本人時(shí)間和能力有限，動(dòng)力仿真分析時(shí)采用的是UM中現(xiàn)有的車(chē)輛模型，研究時(shí)會(huì)產(chǎn)生一定的誤差，所以，對(duì)于車(chē)輛模型方面，有待于車(chē)輛工程等專(zhuān)業(yè)的研究人員建立我國(guó)現(xiàn)實(shí)情況下的車(chē)輛模型。在道路模型方面，仿真時(shí)采用的是剛性路面，從而忽視了實(shí)際中路面對(duì)車(chē)輛的減震作用。如果車(chē)輛通過(guò)道路時(shí)，路面能夠產(chǎn)生相應(yīng)的變形，此時(shí)的結(jié)果將更加符合實(shí)際。不過(guò)對(duì)于剛性路面來(lái)說(shuō)，如果進(jìn)行水泥路面的仿真分析還是比較符合實(shí)際情況的。

參 考 文 獻(xiàn)

[1]ALAN NICHOLSON. Super elevation, side friction and roadway consistency[J].Journal of Transportation Engineering, 1998, 124(5): 411-418.

[2]李磊, 任勇生, 孫愛(ài)芹. ADAMS/Car在汽車(chē)動(dòng)力學(xué)仿真分析中的應(yīng)用[J]. 現(xiàn)代制造技術(shù)與裝備, 2007(1): 25-27.

[3]朱亮亮, 李向國(guó), 李木松. 高速公路平豎曲線(xiàn)組合行車(chē)舒適性研究[J].石家莊鐵道大學(xué)學(xué)報(bào):自然科學(xué)版,2010, 23(3): 68-72.

[4]游克思, 武加恒, 孫璐,等.道路幾何參數(shù)對(duì)車(chē)輛操縱穩(wěn)定性影響[J]. 公路交通科技, 2011, 28(10): 109-117.

[5]伍劍奇.高速公路行車(chē)舒適性評(píng)價(jià)研究[J].華東公路, 2009(6): 61-62.