馬立杰 韓佩鋼 韓子昂

摘要 為了探究瀝青路面結構層厚度對車轍變形的影響規(guī)律，利用ABAQUS有限元軟件對不同結構層厚度情況下以及溫度變化等外界因素進行綜合模擬分析，根據正交試驗原理，系統(tǒng)分析了不同結構層厚度情況下的車轍的變形結果，并對設計年限內的每公里造價對比分析，得出了最優(yōu)的結構層厚度的組合（SMA-13采用60 mm、AC-20采用60 mm、AC-25采用80 mm、水泥穩(wěn)定碎石300 mm），從而預防車轍的發(fā)生降低維修成本。

關 鍵 詞 瀝青路面；車轍；正交分析；ABAQUS軟件；結構層厚度

中圖分類號 U416.217? ? ?文獻標志碼 A

Rutting simulation analysis based on orthogonal combination of ABAQUS asphalt pavement structure layer thickness

MA Lijie， HAN Peigang， HAN Ziang

（College of Civil and Architectural Engineering， North China University of Science and Technology， Tangshan Hebei 063210， China）

Abstract Simulates the external factors such as different thickness of structural layer and the change of temperature with ABAQUS finite element software in order to explore the effect law of asphalt pavement structural layer thickness on rut deformation. According to the principle of orthogonal test， the deformation results of rut under different structural layer thickness are systematically analyzed， and the cost per kilometer in the design year is compared and analyzed， finally the optimal combination of structural layer thickness is obtained （SMA-13 adopts 60mm， AC-20 adopts 60mm， AC-25 adopts 80mm， and cement stabilized macadam 300mm）， which will prevent the car from rutting and reduce the maintenance cost.

Key words asphalt pavement; rut; orthogonal analysis; ABAQUS software; structure layer thickness

0 引言

隨著經濟的快速發(fā)展，以連接城市各主要分區(qū)，交通功能為主的主干路客貨交通運輸增長迅速。交通量和軸載的增加對道路質量要求也提出了更大挑戰(zhàn)。而目前道路破壞中最明顯的屬車轍破壞。車轍是由于瀝青路面受到溫度的影響且在車輛的長期反復碾壓作用下，出現蠕變現象，隨著荷載的重復作用，使路面面層、基層、底基層等結構層發(fā)生豎向及側向的位移，并且使路面發(fā)生渠化，形成高低不平的凹槽，降低道路使用壽命，并且影響到車輛行駛的舒適度。

國內外學者通過有限元軟件在道路工程設計與驗算中投入了大量的研究，而有限元軟件具有強大的建模功能，可以有效地分析出路面的車轍情況。張?zhí)m峰[1]通過建立路面結構模型，耦合了隨時間溫度的變化，對車轍進行了模擬分析。黃菲[2]通過有限元軟件建模并結合在建項目對路面的車轍變形量進行了合理的預估。邱自萍[3]對高模量瀝青混凝土進行了有限元分析，證明了高模量瀝青混凝土可有效地降低路面車轍破壞。何金龍[4]以長沙市道路為模型采用ABAQUS進行分析，得出在氣溫超過25 ℃時，開始影響瀝青混合料的抗車轍性能。Wichtmann等[5]建立了HCA模型，并用該模型預測道路的永久變形。Ali等[6]采用有限元模進行了交通條件對道路損害分析。以上學者均考慮了溫度對車轍的影響，但未對結構層厚度變化對車轍的影響進行模擬分析，所以文中針對瀝青路面結構設計時的結構層厚度進行組合分析，結合溫度變化等外界因素，根據正交試驗原理利用有限元軟件分析不同結構層厚度組合的車轍發(fā)展規(guī)律[7]。從車轍的變形情況和經濟合理的角度考慮，選擇最佳的結構層厚度組合。

1 瀝青路面結構層層厚確定

本研究方法采用正交因素控制變量法，研究上面層（SMA-13）、中面層（AC-20）、下面層（AC-25）、基層（水泥穩(wěn)定碎石基層）的層厚對于車轍深度的影響。分別用A、B、C、D代表SMA-13、AC-20、AC-25、水泥穩(wěn)定碎石基層，并根據交通等級和《公路瀝青路面設計規(guī)范》對不同粒徑瀝青混合料厚度的規(guī)定，無機結合料類基層（無機結合料穩(wěn)定作為底基層時）在重交通時厚度范圍在200～450 mm。結合現有的工程經驗和成本的合理性，SMA的適宜厚度為35～60 mm、AC-20的適宜厚度為60～100 mm、AC-25的適宜厚度為80～120 mm。將其分為3個水平，進行正交試驗進行分析如表1。

2 瀝青路面結構層有限元模型與參數

2.1 瀝青路面結構層模型

由于本次模擬主要針對城市道路的的車轍模擬，而在城市道路中，車道寬度3.5 m居多，所以本模型選用模型為3.5 m寬，3 m深，石灰土厚度0.2 m，具體劃分如圖1所示。

2.2 溫度場參數

瀝青路面的溫度受到氣溫、太陽輻射、熱傳遞及熱交換過程等因素的影響，所以在模擬過程中加入溫度場參數，可更真實地反應在連續(xù)變化溫度的影響下，路面各結構層的溫度變化情況[8-9]。

溫度場結合唐山市高溫季節(jié)24 h內的溫度變化作為代表性氣溫，本模型以唐山市7月24日的溫度變化為代表，最高溫度38 ℃、最低溫度22 ℃、風速3.5 m/s。將溫度參數輸入ABAQUS子程序（DFLUX、FILM）進行溫度場的模擬分析。唐山市7月24日的溫度變化曲線如圖2所示。

瀝青路面溫度場分析熱特性參數表參考廖公云《ABAQUS有限元軟件在道路工程中的應用》[10]一書中的參數如表2所示。

2.3 材料參數

SMA-13、AC-20、AC-25、水泥穩(wěn)定碎石、石灰土、土基等材料的熱特性參數、溫度變化下彈性參數和蠕變參數等均參考廖公云《ABAQUS有限元軟件在道路工程中的應用》一書中的參數如表3所示。

2.4 荷載參數

城市中的主干路設計時速為40～60 km/h，由于存在交叉路口、商業(yè)密集區(qū)、紅綠燈等情況所以選取40 km/h作為模擬參數，本模型采用以靜代動的方式，采用標準雙輪組單軸100 kN作為軸載，輪壓0.7 MPa。雙輪中心間距31.4 cm，輪間距12.8 cm，輪寬18.6 cm。根據《城鎮(zhèn)道路設計規(guī)范》的要求重交通下累計當量軸次應為1 200～2 500萬次/車道，所以本次模擬取其中值1 500萬次當量軸載次數作為模擬參考，對瀝青混合料進行數值模擬分析[11]。用荷載累積時間公式計算速度為40 km/h、累積軸載次數1 500萬次作用下的累積作用時間為259 216 s。把荷載累積作用時間疊加在一天24 h之內，并根據不同時間段的交通量情況將荷載累積作用時間進行分配如表4所示。

3 正交模擬數值結果分析

通過ABAQUS模擬計算了正交因素組合的9種方案，并耦合了一天內溫度變化等的外界影響。通過ABAQUS數值分析計算得到的瀝青路面的變形云圖，總體的變形效果大致相同，但車轍深度上均有所不同，以方案1的各結構層的變形圖為例，如圖3所示。

其余所有的方案的正交因素分析表及車轍模擬結果如表6所示。對9種組合方案的車轍深度進行極差分析，以此確定最優(yōu)的結構層厚度組合。9種組合方案隨時刻的變化的車轍深度如圖4所示。

根據圖4可以看出，9種組合方案總體發(fā)展過程都呈現出：在10：00之前，車轍深度發(fā)展緩慢。在10：00～18：00之間車轍深度呈線性增長趨勢，這段時間內氣溫超過30 ℃，且客貨運量也隨之增大，瀝青路面由于蠕變產生橫向和豎向的推移變形，且短時間內的車輪對瀝青路面的反復荷載作用路面造成了路面的車轍。雖然18：00之后累積作用時間相比10：00之前要高，但氣溫已經開始下降，車轍深度的增長緩慢趨于平穩(wěn)，可以判斷溫度對瀝青路面車轍的影響是一個重要的因素[12]。根據圖5所示，通過不同方案的組合對比，可以發(fā)現改變各個結構層厚度的組合對車轍有明顯的影響[13]。

考慮到城市道路瀝青路面工程造價，通過廣聯達計價軟件查得SMA-13、AC-20、AC-25、水泥穩(wěn)定碎石的工程指導價如表5所示，并根據不同方案組合計算每公里造價，記錄在表6中。

從表7可以看出通過調整瀝青路面結構層的厚度對路面的車轍深度有一定影響且各結構層厚度的對車轍的影響次序為：A、B、D、C，即A（SMA-13）的厚度對車轍的影響最大，其次是B（AC-20）、D（水泥穩(wěn)定碎石）、C（AC-25）。各層位中A對車轍的影響最大，其次是B對車轍的影響較大，C和D的影響均較小，但D的影響程度又高于C。因為在溫度變化時A的抗壓回彈模量明顯優(yōu)于B和C，具有較好的溫度穩(wěn)定性抵抗車轍的變形。而基層是結構中的承重層，主要承受車輛荷載的豎向力，并將上層應力擴散到土基中。且D屬于半剛性基層具有較高的剛度、擴散應力強，在整體的瀝青路面結構中起到一定抵抗豎向變形的能力。

通過車轍深度的極差分析得到的最佳組合為A3、B1、C1、D2，即SMA-13為60 mm，AC-20為60 mm，AC-25為80 mm，水泥穩(wěn)定碎石為300 mm這種組合的抗車轍性能最好。但考慮到工程造價，所以針對不同方案組合的每公里造價進行了極差分析如表8所示。

通過每公里造價的極差分析可以看出，A在整體的城市道路建設造價中起決定性的作用，且高于B和C合起來的造價。所以在考慮選取B、C、D最佳厚度情況下，在A的合理范圍內調整厚度，結合車轍深度和工程造價的對比進行分析。如表9所示。

結合表9中的數據分析，SMA-13的厚度從40 mm增加到50 mm厚度時，車轍深度同比減小3.42%，但成本將同比增加12%。SMA-13的厚度40 mm增加到60 mm厚度時，車轍深度同比減小12.75%，但成本將同比增加24.08%。從50 mm增加到60 mm時車轍深度顯著降低，同比減小了9.65%。所以SMA-13選擇60 mm對預防車轍效果更為明顯。按照城市道路設計規(guī)范主干路瀝青路面設計使用年限為15 a。當SMA-13采用60 mm所產生的車轍深度49.97 mm，在設計年限內，年均產生3.33 mm的車轍，同比采用40 mm所產生的車轍深度57.27 mm，可增加2.19 a的使用年限，即相比可達到17.19 a使用年限。將成本在其設計年限內進行均攤，采用40 mm時年均每公里造價7.425萬元，而采用60 mm時年均每公里造價8.038 9萬元，相比采用40 mm時，每年每公里僅增加0.613 9萬元的成本。雖然增加了成本，但卻延緩了車轍發(fā)展速度，增加了瀝青路面的全壽命周期。

在后期的養(yǎng)護方面，路面養(yǎng)護費用包括大中修、日常養(yǎng)護費等。根據有關學者的調查[14]，養(yǎng)護部門對每段路每年每平米投入9.24～11.13元的日常維修費用，以3.5 m寬路面每公里的年維修費用將高達到3.234～3.895 5萬元，相比每年每公里增加0.613 9萬元的建設成本來說維修費用是一筆很大的開銷。而且延緩車轍發(fā)展的速度，可延后大修養(yǎng)護的時機，在全壽命周期內相應減少大修養(yǎng)護的次數，這將節(jié)省一大筆的大修養(yǎng)護費用。另外從社會角度，減少維修次數的同時可以相應減少由于維修期間引起的交通延誤費、交通繞行費用以及減少由維修期間堵車而造成的交通事故。

所以綜合分析，SMA-13選取60 mm雖然增加了一定的建設成本，但預防車轍產生的效果明顯，增加了瀝青路面的全壽命周期以及再維修養(yǎng)護費用，在社會經濟角度仍具有一定的優(yōu)勢。

4 最佳組合的的抗車轍性能驗證

選取SMA-13為60 mm，AC-20為60 mm，AC-25為80 mm，水泥穩(wěn)定碎石為300 mm。其他參數均不變，在導入ABAQUS模型進行驗算可得到不同結構層隨時間變化的豎向位移和不同深度的最大豎向位移。如圖6、圖7所示。

根據圖6和圖7可以看出最佳組合的最大車轍深度為49.97 mm，相比9種方案中最差組合的車轍深度明顯下降了17.02%，同時也明顯優(yōu)于其他方案組合。通過圖6和圖7可以看出，上面層的豎向位移變形量最大，瀝青路面車轍變形主要發(fā)生在上面層，其次是中面層和下面層，基層、底基層、地基的豎向位移變均較小。

5 結論

文中主要以城市道路中的瀝青路面結構，耦合了唐山市本地的溫度變化條件等因素，運用ABAQUS軟件模擬對城市中瀝青路面的車轍進行了正交分析，得出如下結論。

1）模擬分析了道路在一天不同時段內的變化規(guī)律，在溫度較低的時段，車轍的增長較為緩慢，而在溫度較高和車輛行車荷載密集的時間段，車轍的深度呈線性增長的總體趨勢。得出溫度是車轍產生的重要因素。

2）通過正交模擬分析得出在各結構層材料不變的情況下，通過改變各結構的厚度也可以有效的抵抗的車轍的產生。

3）從車轍減小程度與造價和長期維修成本的綜合分析，確定最佳結構層厚度組合的上面層60 mm、中面層60 mm、下面層80 mm、水泥穩(wěn)定碎石基層300 mm。

4）瀝青路面車轍變形主要發(fā)生在上面層，其次是中面層和下面層，基層、底基層、地基的豎向位移變均較小。

參考文獻：

[1]? ? 張?zhí)m峰. 連續(xù)變溫瀝青路面車轍變形數值模擬[J]. 公路交通科技，2018，35（2）：15-24.

[2]? ? 黃菲. 瀝青路面永久變形數值模擬及車轍預估[D]. 南京：東南大學，2006.

[3]? ? 邱自萍. 高模量瀝青混凝土路面抗車轍行為數值模擬[D]. 西安：長安大學，2009.

[4]? ? 何金龍. 溫度場下城市瀝青路面車轍成因力學機理分析[D]. 長沙：中南大學，2014.

[5]? ? WICHTMANN T，ROND?N H A，NIEMUNIS A，et al. Prediction of permanent deformations in pavements using a high-cycle accumulation model[J]. Journal of Geotechnical and Geoenvironmental Engineering，2010，136（5）：728-740.

[6]? ? ALI B，SADEK M，SHAHROUR I. Finite-element model for urban pavement rutting：analysis of pavement rehabilitation methods[J]. Journal of Transportation Engineering，2009，135（4）：235-239.

[7]? ? 郭芳. 基于時間硬化蠕變模型的組合式基層瀝青路面結構車轍分析[J]. 公路工程，2015，40（6）：214-217，222.

[8]? ? 李輝，黃曉明，張久鵬，等. 基于連續(xù)變溫的瀝青路面車轍模擬分析[J]. 東南大學學報（自然科學版），2007，37（5）：915-920.

[9]? ? 張巧慧，崔凱. 混合瀝青路面結構溫度耦合演變與仿真模擬研究[J]. 計算機應用與軟件，2020，37（7）：320-323.

[10]? 廖公云，黃曉明. 有限元軟件在道路工程中的應用[M]. 南京：東南大學出版社，2008：200-225.

[11]? 李清波. 基于動態(tài)模量的瀝青路面車轍數值模擬研究[D]. 哈爾濱：東北林業(yè)大學，2019.

[12]? 劉俊卿，劉紅，李倩. 變溫條件下考慮車輛-路面相互作用的車轍分析[J]. 應用數學和力學，2017，38（2）：170-180.

[13]? 高虎，楊新華，唐鵬. 瀝青混合料車轍變形的離散元數值模擬[J]. 固體力學學報，2018，39（3）：277-283.

[14]? 張羽. 基于全壽命周期費用理論的瀝青路面大修時機研究[J]. 公路交通科技（應用技術版），2019，15（12）：197-199.