車輛多輪隨機動載作用下柔性瀝青路面的應(yīng)變分析

陳 洋， 戴宗宏， 陳煥明， 劉大維

(青島大學(xué) 機電工程學(xué)院，山東 青島 266071)

?

車輛多輪隨機動載作用下柔性瀝青路面的應(yīng)變分析

陳 洋， 戴宗宏， 陳煥明， 劉大維

(青島大學(xué) 機電工程學(xué)院，山東 青島 266071)

為進行車輛多輪隨機動載作用下柔性瀝青路面的應(yīng)變分析，采用有限元分析軟件ABAQUS建立考慮柔性瀝青路面黏彈性的三維有限元模型，仿真計算車輛多輪隨機動載作用下柔性瀝青路面三個方向的應(yīng)變響應(yīng)，分析車輛多軸作用下路面應(yīng)變時程變化特性，得到路面各點的最大應(yīng)變。結(jié)果表明， 在路面某一位置，三種瀝青材料組成的瀝青各層頂面和下面層底面產(chǎn)生的三個方向的應(yīng)變大小及變化規(guī)律均不相同；瀝青上面層頂面和下面層底面無論是路面長度方向還是寬度方向，各點的三個方向的最大應(yīng)變均不相同，沿路面長度方向，中、后軸車輪荷載產(chǎn)生的應(yīng)變大于前軸車輪荷載產(chǎn)生的應(yīng)變，且橫向應(yīng)變大于縱向和垂向應(yīng)變。研究結(jié)果可為柔性瀝青路面結(jié)構(gòu)設(shè)計與路面壽命分析及預(yù)測提供參考依據(jù)。

隨機動載；柔性瀝青路面；黏彈性；應(yīng)變；有限元

隨著我國高等級公路建設(shè)的不斷發(fā)展，交通量迅速增長，車輛大型化、車輛超載、車輛渠化行駛等問題日趨嚴重。瀝青路面常出現(xiàn)龜裂、橫向裂縫、縱向裂縫、開裂、坑槽、車轍等破壞現(xiàn)象，嚴重降低了道路使用性能。在設(shè)計路面結(jié)構(gòu)、分析路面性能和進行路面疲勞壽命分析時，必須了解車輛荷載作用下路面的應(yīng)變變化，最大應(yīng)變值的產(chǎn)生位置等，才能確保路面的使用壽命。多年來，線彈性層狀理論作為經(jīng)典的分析方法被大量學(xué)者應(yīng)用于路面分析。但對柔性瀝青路面而言，瀝青混合料的黏彈性特性與車輛輪軸數(shù)、輪胎動荷載、作用時間等有關(guān)，線彈性分析模型計算車輛荷載作用下路面結(jié)構(gòu)中應(yīng)力和應(yīng)變，已無法精確地進行路面力學(xué)行為的定性和定量的描述[1-2]。為了更加準確的模擬、預(yù)測車輛引起的路面動態(tài)響應(yīng)，許多學(xué)者考慮了柔性瀝青路面結(jié)構(gòu)的黏彈性和非線性[3-9]。目前在研究車輛動載作用下路面動態(tài)響應(yīng)特性方面，大多對車輛荷載進行了簡化[10-15]，難以全面分析車輛多輪荷載對柔性瀝青路面性能的影響?；诖?，本文建立多軸車輛多輪隨機動載作用下考慮柔性瀝青路面黏彈性的三維有限元模型，分析多軸車輛多輪隨機動載作用下柔性瀝青路面應(yīng)變的動態(tài)變化特性，為柔性瀝青路面結(jié)構(gòu)設(shè)計與路面壽命分析及預(yù)測提供參考依據(jù)。

1 車輛動載作用下路面結(jié)構(gòu)動力學(xué)有限元方程

根據(jù)彈性動力學(xué)理論，車輛多輪動載作用下路面多層結(jié)構(gòu)動力學(xué)有限元方程為：

(1)

式中：M為路面系統(tǒng)質(zhì)量矩陣；C為路面系統(tǒng)阻尼矩陣；K為路面系統(tǒng)剛度矩陣；u為節(jié)點位移向量；F(t)為車輛多輪動載荷矩陣。

路面系統(tǒng)阻尼矩陣采用瑞利阻尼假設(shè)求解：

{C}=α{M}+β{K}

(2)

式中：α和β是與路面結(jié)構(gòu)材料、固有頻率和阻尼比有關(guān)的比例常數(shù)[16]。

2 柔性瀝青路面結(jié)構(gòu)有限元模型

柔性瀝青路面結(jié)構(gòu)如圖1所示，由瀝青上面層(改性瀝青SMA13)、中面層(改性瀝青Sup20)、下面層(普通瀝青Sup25)、基層(級配碎石)和土基組成。

圖1 瀝青路面結(jié)構(gòu)Fig.1 The structure of the asphalt pavement

在建立柔性瀝青路面三維有限元模型時，選取路面長為100 m，寬為18 m，高為5 m，采用ABAQUS有限元軟件中的C3D8R六面體單元對路面結(jié)構(gòu)進行網(wǎng)格劃分?？紤]到計算效率和計算精度，在寬度方向?qū)④囕v車輪接觸區(qū)域網(wǎng)格細化，網(wǎng)格尺寸為0.015 m～0.085 m，車輛車輪非接觸區(qū)域網(wǎng)格尺寸為0.6 m；在長度方向(車輛前進方向)網(wǎng)格尺寸均為0.2 m；在深度方向，瀝青上面層、中面層、下面層分別劃分為2層、3層、4層，基層和土基劃分為4層和15層。所建立的路面結(jié)構(gòu)各層之間的接觸面為完全連續(xù)，整個路面結(jié)構(gòu)單元數(shù)目為588 000個，節(jié)點數(shù)目為624 747個，瀝青路面結(jié)構(gòu)有限元模型如圖2所示。

柔性瀝青路面結(jié)構(gòu)模型各層參數(shù)如表1所示[17]。

圖2 瀝青路面結(jié)構(gòu)有限元模型Fig.2 The finite element model of the asphalt pavement

層類型層厚δ/m彈性模量E/MPa密度ρ/(kg·m-3)改性瀝青SMA130.04黏彈性2450改性瀝青Sup200.06黏彈性2450普通瀝青Sup250.08黏彈性2450級配碎石0.42502000土基4.42601800

對于瀝青層的黏彈性屬性，可根據(jù)瀝青混合料不同溫度和荷載頻率下復(fù)數(shù)模量實驗結(jié)果，通過基本黏彈性關(guān)系轉(zhuǎn)化確定其松弛模量，在ABAQUS中用Prony級數(shù)系列來表征。各瀝青混合料Prony系列的系數(shù)(25℃)如表2所示[17]。

表2 松弛模量Prony系列系數(shù)/MPa

3 車輛多輪動載荷及加載

隨機動載是車輛在不平路面上行駛時，因路面不平度而引起車輛振動產(chǎn)生的。隨機動載用與路面相互作用的車輛輪胎動載荷來描述，可利用質(zhì)量-彈簧-阻尼單元組成的車輛振動模型，或利用多體動力學(xué)方法建立整車行駛動力學(xué)模型，仿真計算輪胎動載荷。為較為準確地獲得多軸車輛多輪輪胎動載荷，本文采用基于路面有理函數(shù)功率譜密度的諧波疊加法建立雙輪轍空間域B級路面不平度模型，并采用多體動力學(xué)仿真軟件SIMPACK建立剛?cè)狁詈系恼囆旭倓恿W(xué)模型(圖3)，計算車輛各軸輪胎動載荷[18-19]。

車輛為某重型自卸汽車，滿載質(zhì)量27 000 kg，前軸載荷70 kN，中、后軸載荷各為100 kN，前軸和中軸距為3.865 m，中軸和后軸距為1.36 m，前軸輪距為2.02 m，中、后軸輪距為1.84 m，輪胎壓力為0.95 MPa。

圖3 整車行駛動力學(xué)模型Fig.3 Dynamic multi-body model of heavy vehicle

圖4為車輛滿載，車速為60 km/h時，前軸、中軸、后軸兩側(cè)輪胎的法向作用力隨行駛距離的變化曲線。

圖4 車輛各軸輪胎法向力Fig.4 The normal force of the vehicle wheels

為了模擬車輛行駛過程中車輛各軸車輪動載的作用效果，在ABAQUS軟件中用子程序vdload施加移動載荷。vdload子程序可以用來定義一組關(guān)于位置時間和速度的函數(shù)，能很好地模擬垂直壓力在路面上的移動加載過程。當車速為60 km/h時，求解的輪胎法向作用力頻率為170.5 Hz，時間間隔為0.005 8 s，也就是每一次取得的載荷在這個時間間隔內(nèi)按不變處理，在下一個時間間隔內(nèi)取另一組數(shù)值。由于vdload子程序?qū)虞d區(qū)域施加的是面載荷，因此需對各個車輪法向載荷進行處理，得到相應(yīng)面載荷數(shù)值，再將其加載到瀝青路面有限元模型上。對于車輛各軸車輪如此眾多的載荷數(shù)據(jù)直接編寫子程序十分困難，本文利用Python腳本語言編寫程序，將獲得的車輪載荷數(shù)據(jù)填充到vdload子程序。

試驗表明，重型車輛接地印跡近似為矩形[20]，本文所用車輛前、后輪胎與路面的接觸印跡分別為0.158 m×0.229 m和0.194 m×0.272 m的矩形，接地面積分別為0.036 1 m2和0.052 6 m2[21]。

4 仿真結(jié)果及分析

圖5為車輛左側(cè)車輪隨機動載作用下路面某一位置的瀝青層三個方向應(yīng)變時程變化曲線。

由圖5可知，三軸車輛通過路面時，當車輪接近與離開該位置時，三個方向的應(yīng)變均發(fā)生了較大變化，三種瀝青材料組成的瀝青各層頂面和下面層底面產(chǎn)生的三個方向的應(yīng)變大小及變化規(guī)律均不相同，三個層頂面的三個方向應(yīng)變變化規(guī)律基本相似，垂向應(yīng)變?yōu)槔瓚?yīng)變，縱向應(yīng)變出現(xiàn)拉壓應(yīng)變變化，橫向應(yīng)變?yōu)閴簯?yīng)變，且橫向壓應(yīng)變最大；瀝青下面層底面(Sup25)應(yīng)變變化與各層頂面的應(yīng)變變化不同，垂向應(yīng)變?yōu)閴簯?yīng)變，縱向應(yīng)變出現(xiàn)壓拉應(yīng)變變化、橫向應(yīng)變?yōu)槔瓚?yīng)變，且橫向拉應(yīng)變最大。由此可見，瀝青路面的破壞一方面是由于瀝青上面層受到較大的壓應(yīng)變而產(chǎn)生的自上而下的破壞；另一方面，由于瀝青下面層底面受拉應(yīng)變作用，將會向上拱起，使得瀝青下面層底面與粒料基層發(fā)生分離，產(chǎn)生一種自下而上的破壞。

從垂向應(yīng)變的時程變化來看(見圖5(a))，瀝青頂面各層在車輪荷載作用過程中均呈現(xiàn)拉應(yīng)變狀態(tài)，且隨著車輛前輪接近該位置，瀝青各層頂面的垂向拉應(yīng)變逐漸變大；當車輪作用在該位置時，垂向拉應(yīng)變變小；在車輛后輪離開后，需要較長時間才能恢復(fù)到初始狀態(tài)；而瀝青下面層底面(Sup25)始終處于壓應(yīng)變狀態(tài)，且變化趨勢與瀝青各層頂面正好相反。

從縱向應(yīng)變的時程變化來看(見圖5(b))，當車輛前輪和中輪接近該位置時，瀝青各層頂面縱向應(yīng)變?yōu)槔瓚?yīng)變；當車輪作用在該位置時，縱向應(yīng)變變?yōu)閴簯?yīng)變；當車輛前輪離開該位置時，縱向應(yīng)變變?yōu)槔瓚?yīng)變；當中輪離開時，由于中、后軸軸距較小，縱向壓應(yīng)變雖有所減小，但隨著后輪荷載的再次作用，縱向壓應(yīng)變再次增大；當后輪離開該位置時，縱向應(yīng)變又變?yōu)槔瓚?yīng)變；而瀝青下面層底面(Sup25)的縱向應(yīng)變變化趨勢與瀝青各層頂面正好相反。上述現(xiàn)象也說明柔性瀝青路面材料具有典型的黏彈性特性。

圖5 瀝青層應(yīng)變時程變化曲線Fig.5 The strain time history curves of the asphalt layers

從縱向應(yīng)變的時程變化可以看出，在車輛駛過的整個過程中，瀝青層經(jīng)歷了拉-壓-拉應(yīng)變的交變變化，在實際瀝青路面結(jié)構(gòu)中，車輛隨機荷載的反復(fù)作用會使路面出現(xiàn)反復(fù)的拉壓交變作用，而這種長期的交變作用必然會導(dǎo)致瀝青路面疲勞破壞的產(chǎn)生。

從橫向應(yīng)變的時程變化來看(見圖5(c))，瀝青各層頂面在車輪荷載作用的過程中均呈現(xiàn)壓應(yīng)變狀態(tài)，而瀝青下面層底面(Sup25)始終處于拉應(yīng)變狀態(tài)。

圖6為車輛前軸和中、后軸左側(cè)車輪通過路面時，沿路面長度方向瀝青上面層頂面(SMA13)和下面層底面(Sup25)最大應(yīng)變變化曲線。

由圖6可看出，在路面長度方向的不同位置，路面三個方向最大應(yīng)變均不相同，中、后軸車輪荷載產(chǎn)生的應(yīng)變(因中、后軸輪距相同，中、后軸車輪荷載產(chǎn)生的最大應(yīng)變?yōu)槎咧械淖畲笾?大于前軸車輪荷載產(chǎn)生的應(yīng)變，且橫向應(yīng)變大于縱向和垂向應(yīng)變。

圖6 路面長度方向瀝青上面層頂面和下面層底面最大應(yīng)變變化曲線Fig.6 The maximum strain curves of asphalt surface and bottom layer along the length direction

圖7為車輛隨機動載作用下沿路面長度和寬度方向瀝青上面層頂面(SMA13)和下面層底面(Sup25)最大應(yīng)變變化曲線。

由圖7可以看出，車輛通過路面時，在瀝青上面層頂面和瀝青下面層底面無論是路面長度方向還是寬度方向，各點的三個方向的最大應(yīng)變均不相同。

5 結(jié) 論

(1) 建立了考慮黏彈性的柔性瀝青路面三維有限元模型，仿真計算了重型車輛多輪隨機動載作用下柔性瀝青路面的動態(tài)應(yīng)變響應(yīng)。

(2) 三軸車輛通過路面時，當車輪接近與離開該位置時，三種瀝青材料組成的瀝青各層頂面和下面層底面產(chǎn)生的三個方向的應(yīng)變大小及變化規(guī)律均不相同，三個層頂面的三個方向應(yīng)變變化規(guī)律基本相似，垂向應(yīng)變?yōu)槔瓚?yīng)變，縱向應(yīng)變出現(xiàn)拉壓應(yīng)變變化，橫向應(yīng)變?yōu)閴簯?yīng)變，且橫向壓應(yīng)變最大；瀝青下面層底面應(yīng)變變化與各層頂面的應(yīng)變變化不同，垂向應(yīng)變?yōu)閴簯?yīng)變，縱向應(yīng)變出現(xiàn)壓拉應(yīng)變變化、橫向應(yīng)變?yōu)槔瓚?yīng)變，且橫向拉應(yīng)變最大。

(3) 在路面長度方向的不同位置，路面三個方向的最大應(yīng)變均不相同，中、后軸車輪荷載產(chǎn)生的應(yīng)變大于前軸車輪荷載產(chǎn)生的應(yīng)變，且橫向應(yīng)變大于縱向和垂向應(yīng)變。

(4) 車輛通過路面時，瀝青上面層頂面和瀝青下面層底面無論是路面長度方向還是寬度方向，各點的三個方向的最大應(yīng)變均不相同。

(5) 研究結(jié)果可為柔性瀝青路面結(jié)構(gòu)設(shè)計與路面壽命分析及預(yù)測提供參考依據(jù)。

圖7 瀝青上面層頂面和下面層底面最大應(yīng)變變化曲線Fig.7 The maximum strain curves of asphalt surface and bottom layer

[1] AL-QADI I L, ELSEIFI M A. Viscoelastic modeling and field validation of flexible pavements[J]. Journal of Engineering Mechanics, 2006, 132(2): 172-178.

[2] WANG H, Al-QADI I L. Combined effect of moving wheel loading and three-dimensional contact stresses on perpetual pavement responses[J]. Transportation Research Record Journal of the Transportation Research Board, 2009, 2095(1):53-61.

[3] WANG Hao, AL-QADI I L．Evaluation of surface-related pavement damage due to tire braking[J]．Road Materials & Pavement Design，2011，11(1)：101-121．

[4] ALI B，SHAHROUR M，SHAHROUR I．Elasto-viscoplastic finite element analysis of the long-term behavior of flexible pavements[J]．Road Materials & Pavement Design，2011，9(3)：463-479．

[5] DE ARAUJ O, JUNIOR P C，SOARES J B，et al．Dynamic viscoelastic analysis of asphalt pavements using a finite element formulation[J]．Road Materials & Pavement Design，2011，11(2)：409-433．

[6] KHAVASSEFAT P. Dynamic response of flexible pavements at vehicle-road interaction[J]. Road Materials & Pavement Design, 2015, 16(2):256-276.

[7] LIN J H．Variations in dynamic vehicle load on road pavement[J]．International Journal of Pavement Engineering，2014，15(6)：558-563．

[8] KHAVASSEFAT P，JELAGIN D，BIRGISSON B．The non-stationary response of flexible pavements to moving loads[J]．International Journal of Pavement Engineering，2015：1-13．

[9] SARKAR A．Numerical comparison of flexible pavement dynamic response under different axles[J]．International Journal of Pavement Engineering，2015：1-11．

[10] BESKOU N D，THEODORAKOPOULOS D D．Dynamic effects of moving loads on road pavements：A review[J]．Soil Dynamics and Earthquake Engineering，2011，31(4)：547-567．

[11] ZHANG Y，ZHANG X．Dynamic response analysis of pavement and subgrade of highmay[C]//Multimedia Technology (ICMT)，2011 International Conference on．IEEE，2011：1034-1036．

[12] 曹衛(wèi)鋒，呂彭民．黏彈性瀝青路面車輛動力響應(yīng)仿真分析研究[J]．計算機仿真，2013，30(4)：173-177．

CAO Weifeng，Lü Pengmin．Research on simulation analysis of dynamic response in viscoelastic asphalt pavement under moving vehicle loads[J]．Computer Simulation，2013，30(4)：173-177．

[13] 董忠紅，呂彭民．移動荷載下黏彈性層狀瀝青路面動力響應(yīng)模型[J]．工程力學(xué)，2011(12)：153-159．

DONG Zhonghong，Lü Pengmin．A model to study the dynamic response of visco-elastic layered system under moving load[J]．Engineering Mechanics，2011(12)：153-159．

[14] 張洪亮，胡長順，許偉清．移動荷載作用下柔性路面的動力響應(yīng)[J]．長安大學(xué)學(xué)報(自然科學(xué)版)，2005，25(5)：6-10．

ZHANG Hongliang，HU Changshun，XU Weiqing．Dynamic response of flexible pavement under moving load[J]．Journal of Chang’an University(Natural Science Edition)，2005，25 (5)：6-10．[15] 陳恩利，劉永強，趙進寶．移動荷載作用下路面動力響應(yīng)試驗研究[J]．振動與沖擊，2014，33(16)：62-67．

CHEN Enli，LIU Yongqiang，ZHAO Jinbao．Experiments on dynamic response of pavement under moving load[J]．Journal of Vibration and Shock，2014，33(16)：62-67．

[16] 廖公云，黃曉明．ABAQUS有限元軟件在道路工程中的應(yīng)用[M]．南京：東南大學(xué)，2008．

[17] 趙延慶，鐘陽．瀝青路面動態(tài)黏彈性響應(yīng)分析[J]．振動與沖擊，2009，28(9)：159-162．

ZHAO Yanqing，ZHONG Yang．Dynamic viscoelastic analysis of asphalt pavements[J]．Journal of Vibration and Shock，2009，28(9)：159-162．

[18] 劉祥銀，陳洋，高攀，等．雙輪轍激勵下多軸重型車輛動載特性仿真分析[J]．振動與沖擊，2015，34(13)：48-52．

LIU Xiangyin，CHEN Yang，GAO Pan，et al．Simulation analysis for multiaxial heavy vehicles’ dynamic load characteristics under bilateral tracks’ road excitation[J]． Journal of Vibration and Shock，2015，34(13)：48-52．

[19] LIU Dawei, JIANG Rongchao, YANG Yuedong，et al．Simulation study of heavy vehicle road-friendliness under bilateral tracks’ excitation[J]．Advanced Materials Research，2013, 680：422-428．

[20] WEISSMAN S. Influence of tire-pavement contact stress distribution on development of distress mechanisms in pavements[J]. Transportation Research Record Journal of the Transportation Research Board, 1999, 1655(1):161-167.

[21] 陳靜. 車輛與路面相互作用的基礎(chǔ)研究[D]. 長春：吉林大學(xué), 2002.

Strain analysis of a flexible asphalt pavement under multi-wheel random dynamic loads of vehicles

CHEN Yang, DAI Zonghong, CHEN Huanming, LIU Dawei

(College of Mechanical & Electronic Engineering, Qingdao University, Qingdao 266071, China)

In order to analyze the strain of a flexible asphalt pavement under vehicle random dynamic loads, the finite element model of the flexible asphalt pavement considering its viscoelasticity was established with the finite element analysis software ABAQUS. The flexible asphalt pavement’s strain responses in three directions under multi-wheel random dynamic loads of vehicles were simulated. The strains’ time histories features were analyzed under random dynamic loads of a multiaxial vehicle and the maximum strain of each point on the pavement was obtained. The results showed that the strains of each asphalt layer’s top surface and bottom surface for 3 layers of different asphalt materials are different in three directions and their variations are also different; on the top surface of the upper; layer and the bottom surface of the lowest layer in the longitudinal direction or the width direction, the maximum strains in 3 directions of each point are different; the strains caused by the middle and rear wheel axle loads are larger than those caused by the front wheel axle loads in the longitudinal direction and the transverse strains are bigger than the vertical and longitudinal strains. The results provided a reference for the design of flexible asphalt pavement structures and the analysis of pavement life.

random dynamic loads; flexible asphalt pavement; viscoelastic; strain; FEM

國家自然科學(xué)基金資助項目(51475248)

2016-02-14 修改稿收到日期：2016-03-01

陳洋 女，碩士生，1991年10月生

劉大維 男，博士，教授，1957年3月

U416.2

A

10.13465/j.cnki.jvs.2016.19.003