孫志林，卓斌，廖中平

（1. 長沙理工大學 道路災變防治及交通安全教育部工程研究中心，湖南 長沙，410114；2. 湖南省高速公路管理局，湖南 長沙，410016）

瀝青路面非線性疲勞損傷特性及應力狀態(tài)演變規(guī)律

孫志林1，卓斌2，廖中平1

（1. 長沙理工大學 道路災變防治及交通安全教育部工程研究中心，湖南 長沙，410114；2. 湖南省高速公路管理局，湖南 長沙，410016）

運用通用有限元軟件ABAQUS及其二次開發(fā)平臺，建立考慮路面材料非線性疲勞損傷的瀝青路面結構數(shù)值計算模型，分析瀝青路面結構在車輛荷載反復作用下路面結構損傷的空間分布、演化規(guī)律以及路面結構內(nèi)部水平正應力的空間分布與演化規(guī)律。研究結果表明：路面結構損傷主要分布在雙輪中心線下靠近基層與底基層層底的區(qū)域，隨著荷載作用次數(shù)增加，基層層底與底基層層底損傷度均增加，且增加幅度逐漸增大；雙輪中心線下靠近基層與底基層層底區(qū)域，考慮損傷的路面結構相比無損路面結構，層底水平拉應力均有所減小，且隨荷載作用次數(shù)增加，水平拉應力逐漸減小，減小的幅度逐漸增大。研究結果可用于路面維修養(yǎng)護中路面破壞區(qū)域及程度的判斷，以及路面設計研究中設計指標的確定。

道路工程；路面結構；損傷力學；疲勞損傷模型；有限元

瀝青路面疲勞破壞機理是長期研究的熱點，有關學者利用斷裂力學、損傷力學、現(xiàn)象學等方法針對瀝青路面疲勞破壞進行了大量研究[1?4]。目前，針對瀝青路面結構疲勞損傷過程的研究特別是基于損傷力學－有限元全耦合方法對瀝青路面結構在交通荷載作用下疲勞損傷過程的相關研究較少[5?6]。為了更準確地分析瀝青路面結構疲勞損傷過程，了解瀝青路面結構在交通荷載下力學狀態(tài)變化規(guī)律，本文作者將運用損傷力學－有限元全耦合方法，基于ABAQUS[7]計算平臺及相應疲勞損傷用戶材料子程序UMAT，分析瀝青路面結構在車輛荷載作用下路面結構損傷度及應力狀態(tài)的演化規(guī)律[8?9]，以便為瀝青路面結構設計、養(yǎng)護維修等提供依據(jù)。

1 疲勞損傷有限元方法

在疲勞損傷力學理論中，一個廣泛應用的非線性損傷演化模型即為Chaboche模型，其損傷演化方程為[10]

式中：D為損傷度，在本研究中取0≤D≤0.5；N為荷載作用次數(shù)（疲勞壽命）；a*，p和q為材料的疲勞損傷特性參數(shù)；σ為拉伸應力。本文采用此模型作為材料疲勞損傷分析模型。

有限元計算模型選擇平面應變模型，此時，有限元方法的本構方程為

式中：σx和σz分別為x和z方向的正應力；τxz為剪應力；E為彈性模量；μ為泊松比；D為損傷度；εx和εz分別為x和z方向的正應變；γxz為剪應變。

為了準確地進行損傷力學有限元疲勞損傷累積分析，采用全耦合方法，即每隔一定的應力循環(huán)次數(shù)對單元剛度矩陣重新計算，以反映疲勞損傷累積效應對單元剛度矩陣的影響。在有限元軟件ABAQUS提供的用戶子程序接口上，采用FORTRAN77編寫用戶子程序UMAT，以反映疲勞損傷對單元剛度的影響，即耦合疲勞損傷的材料模型進行相應的編譯連接，用于進行疲勞損傷分析。所用計算方法中，用戶子程序UMAT中的單元剛度矩陣在每個增量步中都進行調(diào)用，每次調(diào)用之后按照疲勞損傷演化規(guī)律及時更新剛度矩陣中的疲勞損傷度D，因此，該方法屬于全耦合解法，解法準確度要高于全解耦和半耦合解法，具有明確的物理意義。

2 結構計算模型

2.1 路面結構

采用我國目前常用的半剛性路面結構形式。路面結構各組成部分的材料如下：面層材料為瀝青混凝土，基層材料為水泥穩(wěn)定碎石，底基層材料為二灰土。各層材料的參數(shù)見表1。

表1 路面結構計算參數(shù)Table 1 Calculation parameters of pavement structure

有限元計算模型為平面應變模型，其長×寬為10 m×10 m。假定模型兩側與底部完全約束，路表為自由面，沒有約束?？拷s束邊界的位置各種力學響應很小，可忽略不計，故假定模型兩側與底部完全約束合理。

2.2 材料疲勞損傷參數(shù)

張行在有關拉伸疲勞試驗中，針對Chaboche提出的模型，得出：

式中：KT為應力集中系數(shù)，在本文中取1；Ncr為裂縫形成疲勞壽命。

依據(jù)相關疲勞方程和試驗數(shù)據(jù)[11?12]，可以得到p和c的值，并取q=0，由式（3）可計算出a*。整理計算結果如表2所示。

2.3 荷載條件

本文中，選取車輛荷載寬度W=156 mm，雙輪中心距B=312 mm。選取荷載集度P=0.16 MPa，計算可得施加荷載為標準軸載，以保證各層結構的應力水平和路表彎沉在通常范圍內(nèi)，不影響路面結構力學響應的規(guī)律性分析。

表2 Chaboche疲勞損傷模型參數(shù)Table 2 Chaboche fatigue damage model parameters

3 損傷場分析

3.1 損傷度分布云圖

圖1所示為上述路面結構在車輛荷載作用600萬次后，損傷度在路面結構中的分布云圖。從圖1可以看出面層沒有損傷，基層與底基層均發(fā)生疲勞損傷，在雙輪下靠近層底的位置損傷最嚴重。

圖1 損傷度分布云圖Fig. 1 Cloud picture of damage degree distribution

3.2 損傷度空間分布規(guī)律

利用ABAQUS有限元計算軟件，計算可得瀝青路面結構在車輛荷載作用600萬次后，損傷度空間分布規(guī)律如圖2和圖3所示，其中，圖2所示為雙輪中心線下垂直方向基層與底基層損傷度分布規(guī)律，圖3所示為基層與底基層層底沿水平方向損傷度分布規(guī)律。

從圖2可以看出：在雙輪中心線下，無論是基層還是底基層，越靠近層底，損傷度越大，從層底往上損傷迅速減小；對于基層，距層底略5 cm處，損傷度已接近于0；對于底基層，損傷自下往上逐步遞減，直到接近距層底約20 cm的位置損傷度接近0?？梢姷谆鶎訐p傷范圍比基層的要廣，但損傷程度比基層的低。

圖3所示為基層層底與底基層層底沿水平方向損傷度分布規(guī)律。從圖3可以看出：損傷最大的位置出現(xiàn)在雙輪中心線下，往兩側損傷度減小，說明該位置是整個結構最先出現(xiàn)破壞的地方；當損傷達到一定程度時，便會在此位置逐漸形成裂紋直至開展形成宏觀裂縫；對于基層，距中心線25 cm的位置損傷度已經(jīng)接近于0；對于底基層相對要緩慢點，不過在接近距中心線40 cm的位置損傷度也已經(jīng)接近于0。該結果與損傷度云圖反映的結果也一致。這種損傷區(qū)域集中的原因主要是非線形損傷模型中考慮了當前損傷度對損傷速率的影響。

圖2 雙輪中心線下基層與底基層損傷度Fig. 2 Damage degree of base and sub-base along middle line between double wheels

圖3 基層與底基層層底損傷度Fig. 3 Damage degree at bottom of base and sub-base

3.3 結構層層底損傷度隨荷載作用次數(shù)的變化

層底損傷度隨荷載作用次數(shù)的變化規(guī)律如圖4所示。從圖4可以看出：基層層底的損傷度隨著荷載作用次數(shù)的增加而增加；當荷載作用次數(shù)從0增加到50萬次時，基層層底損傷度從0增加到0.018；當荷載作用次數(shù)從650萬次增加到700萬次時，損傷度增加了0.071；隨著荷載作用次數(shù)的增加損傷度增加的幅度越來越快。底基層損傷度在初始階段其增長速度比基層的小，但后期增長速度比基層的大。

圖4 層底損傷度隨荷載作用次數(shù)的變化規(guī)律Fig. 4 Change law of damage degree at layer bottom with load cycles

選用Chaboche非線性疲勞損傷演化模型進行分析。該模型中損傷度演化率與當前損傷度有關，當前損傷度越大，損傷演化率越大；同時，隨著損傷度的增大，分析點的水平拉應力會減小，有減小損傷演化率的趨勢；當損傷度導致的損傷演化率增大值大于水平拉應力減小導致的損傷演化率減小值時，損傷度便會快速增大。

4 應力場分析

4.1 路面結構水平正應力空間分布規(guī)律

瀝青路面結構在車輛荷載作用600萬次后，基層和底基層雙輪中心線下水平正應力和層底水平正應力分布規(guī)律分別如圖5和圖6所示。

從圖5可以看出：雙輪中心線下基層的水平正應力在整個厚度內(nèi)不再呈線性分布；在靠近層底的一側，由于受拉應力損傷的影響，與無損結構相比有所減小；但當距層底距離大于5 cm時，由于水平拉應力較小或者只有水平壓應力，結構層的損傷很小或者沒有，因此，水平拉應力保持線性變化。底基層內(nèi)水平正應力在垂直方向的分布與基層的類似，在靠近層底的位置，由于受拉應力損傷的影響有所減小，當距底基層層底的距離大于7 cm以后則呈線性變化。但在底基層內(nèi)水平正應力始終為拉應力。

從圖6可以看出：基層層底水平正應力沿水平方向的分布呈現(xiàn)出下凹的應力分布圖形。從雙輪中心線的位置往兩側水平拉應力逐步增大，之后增大幅度有所減小，直到距中心線20 cm處水平拉應力不再增加，再往兩側水平拉應力開始減小，且減小的幅度很快增大；底基層水平正應力沿水平方向的分布與基層的類似，但底基層在中心線及其附近水平正應力減小幅度較大。

圖5 雙輪中心線下基層與底基層水平正應力Fig. 5 Horizontal normal stress of base and sub-base along middle line between double wheels

圖6 基層與底基層層底水平正應力Fig. 6 Horizontal normal stress at bottom of base and sub-base

4.2 結構層層底水平正應力隨荷載作用次數(shù)的變化

瀝青路面結構雙輪中心線下層底水平正應力隨荷載作用次數(shù)的變化規(guī)律如圖7所示。

從圖7可以看出：隨著荷載作用次數(shù)的增加，基層層底與底基層層底的水平正應力均減小，且減小的幅度逐步增大。這是因為隨著荷載次數(shù)的增加，損傷減小的越來越快，損傷處模量的減小也會變快，由于應力的重分布作用，相應的水平正應力的減小也會逐漸加快。

圖7 層底水平正應力隨荷載作用次數(shù)的變化規(guī)律Fig. 7 Change law of horizontal normal stress at layer bottom with load cycles

5 結論

（1） 路面結構的損傷主要分布在雙輪中心線下靠近層底的區(qū)域；隨著車輛荷載作用次數(shù)增加，基層層底與底基層層底的損傷度增加，且增長的幅度逐漸增大。依據(jù)結構損傷度的分析，可以實時了解結構在不同時期、不同位置損傷度的分布和損傷程度，從而確定對路面結構進行養(yǎng)護維修的時間和位置以及所采用的措施。

（2） 考慮損傷的路面結構相比不考慮損傷的路面結構，在雙輪中心線下靠近層底的區(qū)域內(nèi)，水平拉應力均有所減?。辉陔p輪中心線下基層與底基層層底的位置，隨著荷載作用次數(shù)增加，水平拉應力逐漸減小，且減小的幅度逐漸增大。路面結構層底應力水平并不是固定不變的，因此，有必要確定更合理的路面結構設計控制指標進行路面結構設計。

[1] Carpenter S H, Ghuzlan K, Shen S. Fatigue endurance limit for highway and airport pavements[J]. Transportation Research Record, 2003, 18（32）： 131?138.

[2] Castro M, Sanchez J A. Estimation of asphalt concrete fatigue curves： A damage theory approach[J]. Construction and Building Materials, 2008, 22（6）： 1232?1238.

[3] SI Zhiming, Little D N, Lytton R L. Characterization of microdamage and healing of asphalt concrete mixtures[J]. Journal of Materials in Civil Engineering, 2202, 14（6）： 461?470.

[4] 葛折圣, 黃曉明. 運用損傷力學理論預測瀝青混合料的疲勞性能[J]. 交通運輸工程學報, 2003, 3（1）： 40?42.

GE Zhesheng, HUANG Xiaoming. Prediction of asphalt mixtures fatigue properties using damage mechanics[J]. Journal of Traffic and Transportation Engineering, 2003, 3（1）： 40?42.

[5] SHU Xiang, HUANG Baoshan, Dragon V. Laboratory evaluation of fatigue characteristics of recycled asphalt mixture[J]. Construction and Building Materials, 2008, 22（7）： 1323?1330.

[6] 關宏信, 鄭健龍, 張起森. 瀝青混合料黏彈性疲勞損傷模型研究[J]. 力學與實踐, 2007, 29（2）： 50?53.

GUAN Hongxin, ZHENG Jianlong, ZHANG Qisen. Viscoelastic fatigue damage model of bituminous mixture[J]. Mechanics in Engineering, 2007, 29（2）： 50?53.

[7] 石亦平, 周玉蓉. ABAQUS有限元分析實例詳解[M].北京： 機械工業(yè)出版社, 2006： 23?36.

SHI Yiping, ZHOU Yurong. Detailed description of ABAQUS finite element analysis example[M]. Beijing： China Machine Press, 2006： 23?36.

[8] Tschegg E K, Jamek M, Lugmayr R. Fatigue crack growth in asphalt and asphalt-interfaces[J]. Engineering Fracture Mechanics, 2011,78（6）： 1598?1602.

[9] LI Qiang, Lee H J, Kim T W. A simple fatigue performance model of asphalt mixtures based on fracture energy[J]. Construction and Building Materials, 2012, 27（1）： 605?611.

[10] 鄭健龍, 呂松濤. 瀝青混合料非線性疲勞損傷模型[J].中國公路學報, 2009, 22（5）： 21?28.

ZHENG Jianlong, Lü Songtao. Nonlinear fatigue damage model for asphalt mixtures[J]. China Journal of Highway and Transport, 2009, 22（5）： 21?28.

[11] 孫榮山, 汪水銀. 級配變化對水泥穩(wěn)定碎石材料疲勞性能影響的研究[J]. 公路交通科技, 2007, 24（6）： 57?61.

SUN Rongshan, WANG Shuiyin. Research on fatigue performance of the different gradation of cement stabilized crushed stone[J]. Journal of Highway and Transportation Research and Development, 2007, 24（6）： 57?61.

[12] 謝軍, 郭忠印. 瀝青混合料疲勞響應模型試驗研究[J]. 公路交通科技, 2007, 24（5）： 21?25.

XIE Jun, GUO Zhongyin. Researching on fatigue model of asphalt mixtures[J]. Journal of Highway and Transportation Research and Development, 2007, 24（5）： 21?25.

（編輯 陳燦華）

Nonlinear fatigue damage characteristics and stress evolution law of asphalt pavement

SUN Zhilin1, ZHUO Bin2, LIAO Zhongping1

（1. Engineering Research Center of Catastrophic Prophylaxis and Treatment of Road & Traffic Safety, Changsha University of Science &Technology, Changsha 410114, China;

2. Hunan Expressway Administration Bureau, Changsha 410016, China）

Using general finite element software ABAQUS and the secondary development platform, an asphalt pavement numerical model which considers the pavement material nonlinear fatigue damage was established to analyze the asphalt pavement structure damage’s spatial distribution and evolution, the spatial distribution and evolution of the horizontal stress within the pavement structure. The results show that the pavement structural damage mainly distributes in a specific damage zone located at the bottom of the base and sub-base and along the middle line between the double wheels. With the increase of the load cycles, the damage degree at the bottom of the base and the sub-base increases, and the rate of which gradually increases. The horizontal tensile stress in the specific damage zone of a pavement structure considering damage is less than that of the pavement structure without considering damage. Simultaneously, this horizontal tensile stress gradually decreases, and the rate increases with the increase of load cycles. The results can be used to judge failure region and degree of the pavement structure on the asphalt pavement maintenance, and to determine the design index of the asphalt pavement design research.

road engineering; pavement structure; damage mechanics; fatigue damage model; finite element

U416.01

A

1672?7207（2014）02?0576?05

2013?05?25；

2013?09?10

國家自然科學青年基金資助項目（51008038，51208066）

孫志林（1979?），男，湖南慈利人，博士，講師，從事路面工程研究；電話：15084967758；E-mail：sunzhilin1979@qq.com