基于級配設計的瀝青混合料動態模量預測模型研究*

劉 彬 羅 蓉 陳 輝 馮光樂

(武漢理工大學交通學院1) 武漢 430063) (湖北省公路工程技術研究中心2) 武漢 430063) (湖北省交通廳工程質量監督局3) 武漢 430014)

0 引 言

動態模量是瀝青混凝土材料設計的關鍵參數之一，可運用其預測瀝青混凝土材料在荷載作用下的力學響應.但是，瀝青混凝土是典型的黏彈性材料，其模量并非固定值，而是隨著加載溫度和頻率的變化而變化.目前，獲得瀝青混凝土動態模量的方法主要有試驗法和經驗公式法.試驗法測試瀝青混凝土的動態模量可以采用重復荷載間拉伸試驗，但是試驗法所需要用到的儀器設備較為昂貴，并且試驗需要在不同溫度和頻率下進行，耗時較長；經驗公式法是在大量試驗數據基礎上，通過軟件擬合，得到含有材料參數的方程對混合料的動態模量進行預測，運用較多的經驗公式有Witczak模型、Hirsch模型等[1].經驗公式的建立所用到的數據是在一定的試驗條件下得到的，若實際的條件與公式建立的條件不一致，那么預測得到的結果與實際的可能相差很大.并且從本質上看，試驗法和經驗法都是建立在瀝青混合料的宏觀力學層面，而瀝青混合料的宏觀力學性質與其材料組成和內部結構是密切相關的，因此，有必要從微觀層面建立合適的模型，使之能夠在設計之初準確預測瀝青混合料力學性質，尤其是預測混合料的動態模量，以達到預測瀝青路面路用性能的目的.

瀝青混合料是典型的多相復合材料，近些年來，不少學者從復合材料微觀力學層面對瀝青混合料的動態模量進行預測.Li等[2]將瀝青混合料視為瀝青基體與集料夾雜相的兩相復合材料，形成一個雙層的嵌入模型，提出了考慮級配與顆粒大小的模型對瀝青混合料的動態模量進行預測；Li等[3]提出兩步法預測混合料的動態模量；朱興一等[4]基于復合材料微觀力學建立了瀝青混合料多相的兩層嵌入式微觀力學模型，通過等效夾雜理論，將各檔集料與空隙逐一均勻化投入基體中，預測混合料的動態模量；郭乃勝等[5]基于復合材料有效模量理論，建立瀝青混合料動態模量預測的三相微觀力學模型；Luo等[6]認為瀝青混合料有瀝青、集料以及空隙三相組分組成，在Hashin等[7]提出的復合材料理論上推導了瀝青混合料自洽微觀力學模型，對混合料動態模量進行預測.

本文基于文獻[6]提出的瀝青混合料自洽微觀力學模型，結合Shu等[8]提出的瀝青混合料動態模量預測模型，建立的考慮集料尺寸、級配，以及空隙的微觀力學模型，對瀝青混合料動態模量進行預測，并通過與室內動態模量試驗結果對比，分析模型預測結果與實際結果的差別.

1 瀝青混合料力學分析模型

1.1 瀝青混合料自洽微觀力學模型

瀝青混合料是一個復雜的多相體系，Luo等基于Hashin-Shtrikman邊界，推導得出了瀝青混合料自洽微觀力學模型，模型結構為

(1)

(2)

式中：K1,K2,K3，K*分別為集料、瀝青、空氣及瀝青混合料的體積模量；G1,G2,G3，G*分別為集料、瀝青、空氣及瀝青混合料的剪切模量；c1,c2,c3分別為混合料中集料、瀝青及空氣的體積分數.

在上述模型中，通過集料、瀝青、空氣三相材料的材料參數，即可求解出混合料材料參數；而通過混合料及任意兩相材料的材料參數，也可求解出第三相材料的參數.基于該思想，Luo等首先通過混合料、瀝青，以及空氣的材料參數計算出集料的參數，再將計算得到的集料參數以及試驗得到的瀝青和空氣參數代入模型，計算出混合料的材料參數，與實際測試得到的混合料參數相比較.通過這種自洽的方法，驗證了該模型是行之有效的.

利用上述自洽的方法，通過混合料、瀝青，及空氣的參數計算得出的集料的體積模量和剪切模量是一個隨著加載時間和加載頻率變化的量，這與普遍認為的集料是彈性材料是不相符的.而在實際預測過程中，要通過三相材料的材料參數預測混合料的參數，只能通過假設集料的模量進行計算.并且，瀝青混合料設計級配對其力學性能影響較大，而該模型中忽略了混合料級配這一關鍵因素，因此，需要進一步改進上述模型，建立考慮混合料級配的微觀力學模型，利用其對混合料動態模量進行預測.

1.2 基于級配的微觀力學模型

瀝青混合料中，集料作為最主要的部分起到骨架和填充的作用，瀝青則作為膠結料使得集料顆粒與顆粒之間相互粘結.為了考慮集料級配對混合料力學性質的影響，首先將混合料視為由無數個裹覆了一定厚度的瀝青膜的集料顆粒及自由分布的空隙組成，并且從混合料中取出一個微單元，在該單元中空隙的分布是隨機分散在顆粒四周，見圖1.

圖1 瀝青混合料內部顆粒單元

集料形狀較為復雜，為了簡化計算過程，將集料視為球形顆粒，并且瀝青混合料中集料尺寸差別較大，即使是同一篩孔上的集料大小也有較大的差別，為了減小這種差別，采取平均半徑的方法計算每一個篩孔上的集料半徑a，即

(3)

式中：ai為第i號篩網上集料顆粒的平均半徑；di，di+1分別是第i號篩網及i+1號篩網的尺寸大小.

假設裹覆在集料顆粒表面的瀝青膜厚度與集料顆粒大小無關，及每一顆集料所裹覆的瀝青膜厚度是相同的，則瀝青膜厚度為[9]

(4)

(5)

(6)

(7)

(8)

瀝青混合料的模量E可以視為由無數個顆粒單元組成的，因此，瀝青混合料的有效彈性模量E可以采用積分的形式為

(9)

式中：amin，amax為混合料中最小和最大集料粒徑；f(a)為集料級配函數.

由于式(9)計算積分較為復雜，可將其簡化為

(10)

2 結果分析

2.1 原材料參數確定

為了驗證上述模型，選用SBS改性瀝青和70#基質瀝青2種瀝青，以及石灰巖和輝綠巖2種集料設計了4種混合料，其中瀝青和集料的各項性質均滿足規范要求.利用旋轉壓實儀分別成型了AC-20石灰巖試件、AC-13輝綠巖試件，初始試件尺寸為直徑150 mm、高度170 mm的圓柱體，經過鉆芯切割后為直徑100 mm、高150 mm的圓柱體試件，其試驗級配見表1.

表1 試驗級配

其中，石灰巖和輝綠巖2種集料的彈性模量E1均假設為40 000 MPa，并且不考慮加載過程中集料泊松比的變化，認為其泊松比為常數，μ1=0.25.通過式(11)～(12)確定集料的剪切模量G1和體積模量K1.通過動態剪切流變試驗(DSR)確定SBS改性瀝青和70#基質瀝青的動態剪切模量G2，并利用CAM模型[10]繪制瀝青剪切模量主曲線，見圖2.通過式(13)計算瀝青的體積模量K2.一般而言，瀝青的泊松比取為0.5，但是當瀝青泊松比取為0.5時，通過式(13)計算得出的瀝青體積模量為無窮大，因此，為了避免這種情況，本文瀝青的泊松比取為μ2=0.49.空氣的剪切模量和體積模量均取為0.

圖2 瀝青剪切模量主曲線

(11)

(12)

(13)

2.2 試驗結果分析

采用單軸壓縮動態試驗測試瀝青混合料的動態模量，試驗采用MTS試驗系統對試件施加半正弦荷載，施加荷載模式采用應力控制，同時控制試件應變小于10-4試件處于無損階段.根據美國AASHTO TP62-07試驗規程，選取試驗溫度分別為-10，4.4，21.1，37.8和54.4 ℃，試驗頻率分別為0.1,0.5,1.0 ,5.0,10和25 Hz.按照上述方法測試得到各溫度各頻率下瀝青混合料的動態模量，根據時溫等效原理，將各溫度頻率下的動態模量進行平移，繪制動態模量曲線.

將以上得到的顆粒單元中的各個參數代入式(1)～(2)，利用excel中的規劃求解功能，將其中一個方程作為目標值，另一個方程作為約束，單元顆粒的剪切模量和體積模量作為變化值，求解出單元顆粒的剪切模量和體積模量，利用式(8)即可求解得到單元顆粒的彈性模量，通過該方法求得各溫度頻率下混合料的動態模量，并繪制曲線.將預測得到的瀝青混合料動態模量曲線與實際測試結果進行對比，見圖3.通過判定系數判斷預測動態模量與實測數據的擬合程度，本文選取的四種混合料類型預測結果與實測結果的判定系數R2見表2.

圖3 瀝青混合料動態模量預測結果與試驗結果

集料類型瀝青類型判定系數R2石灰巖70#基質瀝青0.994 0SBS改性瀝青0.993 9輝綠巖70#基質瀝青0.991 1SBS改性瀝青0.989 9

由表2中的判定系數可知，本文提出的基于集料級配設計的微觀力學模型預測出的瀝青混合料的動態模量與室內試驗測試得到的動態模量較為接近，4種混合料類型的預測結果判定系數都達到了0.98以上，說明該模型能夠較為準確的預測瀝青混合料動態模量.

基質瀝青混合料的預測效果要優于改性瀝青混合料的預測效果，這可能由于預測模型沒有考慮改性瀝青中的改性劑在混合料中起到的作用.而由圖3可知，不論是石灰巖還是輝綠巖，使用改性瀝青均能有效提高混合料的動態模量，結合圖2可知，瀝青作為膠結料在混合料中起到較大的作用，提高瀝青彈性模量，能夠在一定程度上提高混合料的模量.

由圖3可知，4種混合料預測得到的曲線均要低于實測得到的模量曲線.出現這種情況的原因可能有：①本模型在考慮集料形狀、瀝青膜厚度、孔隙分布形式的假設與實際不符導致的；②本文假設的集料的模量、泊松比與實際情況有差別；③考慮瀝青體積分數的時候，沒有考慮集料吸收的瀝青體積，且瀝青泊松比假設存在一定的差異，研究表明 瀝青材料在較高溫度揮著較大的應力水平下或是瀝青產生裂縫時，其泊松比隨著加載時間的延長而增長[11].并且，在頻率較高時，預測與實測的差異更加的明顯.這一現象可以解釋為本文所用的預測模型很大程度強化了復合材料體系的整體性，而弱化了每一相組分在不同溫度和頻率下所表現出的性質.瀝青混合料是典型的黏彈性材料，在低溫或高頻時材料性質主要由瀝青性質起主導作用，而在高溫或低頻時則由集料性質起主導作用.因此，本預測模型弱化了瀝青材料在高頻時對混合料模量的主導作用，因此出現隨著頻率增加，預測結果與實測結果差別增大的現象.

在本次試驗中，集料的模量和泊松比均是假設的，其值為參考文獻中所提供的數值.為了研究集料模量和泊松比大小對混合料模量的影響，本文首先保持集料泊松比μ1=0.25不變，集料的模量取為30,40,50,60,70,80及90 MPa；然后保持集料模量為40 MPa，泊松比分別選為0.22,0.25,0.28,0.31和0.34，并在21.1 ℃,10 Hz的條件下通過預測模型預測瀝青混合料的動態模量，預測結果見圖4.由圖4可知，隨著集料模量的增加，混合料模量變化并不明顯；而隨著集料泊松比的增加，預測得到的混合料模量有所提升，但是總體影響并不明顯，因此，可以認為在本文所建立的預測模型中，假設的集料模量和泊松比是有效的，對于混合料模量預測影響較小.

圖4 混合料模量隨不同因素的變化

3 結 論

1) 根據Luo等提出的自洽微觀力學模型的基礎上，基于混合料中顆粒單元復合模型，提出了考慮集料尺寸以及級配組成設計的瀝青混合料動態模量預測方法，對各溫度頻率下瀝青混合料的動態模量進行預測，并與室內測試結果進行對比.結果表明，雖然其預測值要小于實測值，尤其是在低溫高頻條件下，二者差異較為明顯，但是該方法能夠有效的預測混合料的動態模量，預測結果判定系數達到0.98以上.

2) 在文章所提出的模型中，集料彈性模量以及泊松比對于混合料的模量影響較小，本文所選取的E1=40 000 MPa，μ1=0.25是合理的.

3) 瀝青的彈性模量的增強，能夠提高混合料的動態模量，對于集料表面裹覆的瀝青膜厚度以及被集料吸收的瀝青的體積，對于混合料模量預測值有較大影響,因此，進一步將更為準確的瀝青膜厚度和瀝青真實的體積分數代入微觀力學模型中進行計算式下一步工作的重點.

參考文獻

[1] 胡旭東,張起森,范勇軍.HMA動態模量Witczak和Hirsch預測模型[J].中外公路,2006,26(6):204-207.

[2] LI G, LI Y, METCALF J B, et al. Elastic modulus prediction of asphalt concrete[J]. Journal of Materials in Civil Engineering, 1999,11(3):236-241.

[3] LI Y, METCALF J B. Two-step approach to prediction of asphalt concrete modulus from two-phase micromechanical models[J]. Journal of Materials in Civil Engineering, 2005,17(4):407-415.

[4] 朱興一,黃志義,陳偉球.基于復合材料細觀力學模型的瀝青混凝土彈性模量預測[J].中國公路學報,2010,23(3):29-34.

[5] 郭乃勝,趙穎華.基于細觀力學的瀝青混合料動態模量預測[J].工程力學,2012,29(10):13-19.

[6] LUO R, LYTTON R L. Self-consistent micromechanics models of an asphalt mixture[J]. Journal of Materials in Civil Engineering, 2011,23(1):49-55.

[7] HASHIN Z, SHTRIKMAN S. A variatinal approach to the elastic behavior of multiphase materials[J]. Journal of the Mechanics & Physics of Solids, 1963,11(2):127-140.

[8] SHU X, HUANG B S. Dynamic modulus prediction of HMA mixtures based on the viscoelastic micromechanical model [J]. Journal of Material in Civil Engineering, 2008,20(8):530-538.

[9] 黃曉勇.瀝青混合料動態模量預測方法研究[J].公路交通技術,2015(4):32-36.

[10] YUSOFF N I M, SHAW M T, AIREY G D. Modelling the linear viscoelastic rheological properties of bituminous binders[J]. Construction & Building Materials, 2012,25(5):2171-2189.

[11] 張裕卿,黃曉明.基于微觀力學的瀝青混合料黏彈性預測[J].吉林大學學報(工學版),2010,40(1):52-57.