易軍艷，SHEN Shihui，MUHUNTHAN Balasingam，馮德成

(1.哈爾濱工業(yè)大學(xué)交通科學(xué)與工程學(xué)院，150090 哈爾濱;2.華盛頓州立大學(xué)土木與環(huán)境工程學(xué)院，99164 美國華盛頓州)

多孔瀝青混合料粘彈塑性損傷模型

易軍艷1，SHEN Shihui2，MUHUNTHAN Balasingam2，馮德成1

(1.哈爾濱工業(yè)大學(xué)交通科學(xué)與工程學(xué)院，150090 哈爾濱;2.華盛頓州立大學(xué)土木與環(huán)境工程學(xué)院，99164 美國華盛頓州)

為合理描述多孔瀝青混合料在中低溫度外界荷載作用下的力學(xué)特性，基于增量型本構(gòu)方程，采用Weibull損傷函數(shù)、廣義Maxwell粘彈模型與D-P塑性模型，構(gòu)建了粘彈塑性損傷模型.以此模型為分析手段，對不同溫度和加載速率下的單軸壓縮應(yīng)力-應(yīng)變曲線進行擬合，并分析溫度與加載速率對模型參數(shù)的影響規(guī)律.分析結(jié)果表明:多孔瀝青混合料粘彈參數(shù)隨著溫度的降低逐步退化成彈性參數(shù)，塑性模型中的體積模量和剪切模量也隨溫度呈現(xiàn)出明顯的粘彈特性，塑性應(yīng)變產(chǎn)生時對應(yīng)的應(yīng)變值與損傷應(yīng)變闕值基本保持一致，溫度及加載速率對于混合料的損傷擴展也有顯著影響.構(gòu)建的理論模型可以有效表征多孔瀝青混合料在常溫和低溫下受荷時的力學(xué)損傷特性.

道路工程;多孔瀝青混合料;粘彈性;塑性;損傷

由于抗滑、降噪、雨天行駛安全等優(yōu)點，多孔瀝青混合料愈來愈得到國內(nèi)外道路工作者的重視［1-3］.但是，由于空隙較大，容易堵塞，且由于結(jié)構(gòu)組成與常規(guī)密級配瀝青混合料不同，強度僅靠內(nèi)摩阻力與相對更弱的粘聚力提供，耐久性問題一直是制約多孔瀝青混合料進一步應(yīng)用的瓶頸，尤其在季節(jié)性冰凍地區(qū)使用規(guī)模更?。?-8］.多孔瀝青混合料是一種典型的粘彈性材料，在常溫下受荷后瀝青或膠漿可以通過松弛作用消散掉部分應(yīng)力.但當(dāng)溫度降低時，瀝青類材料逐步開始顯示出彈脆性，一旦破壞就會導(dǎo)致嚴重的結(jié)構(gòu)松散［9-10］.

鑒于多孔瀝青混合料的以上特點，在季凍區(qū)應(yīng)用時首先需要明確其力學(xué)特性與損傷機制.目前對材料的損傷機理研究可通過宏微觀分析試驗手段得到，如圖像處理、工業(yè) CT掃描技術(shù)等［11-12］.此外，理論模型也可用于損傷特性分析［13-14］.但是，以往的模型分析大多試驗溫度為中高溫，對于多孔瀝青混合料在較低溫度下的損傷特點，國內(nèi)外研究仍較少涉及.

由此，本文基于粘彈塑性損傷模型，對多孔瀝青混合料在季節(jié)性冰凍地區(qū)最不利季節(jié)的溫度作用下力學(xué)行為進行分析，以期更好地理解其破壞特性，為進一步的應(yīng)用提供參考.

1 粘彈塑性損傷本構(gòu)方程

在構(gòu)建本構(gòu)方程之前，作如下假設(shè):1)假設(shè)多孔瀝青混合料初始損傷為0;2)多孔瀝青混合料的損傷主軸與應(yīng)力主軸和應(yīng)變主軸重合.

1.1 損傷模型

瀝青混合料在成型后內(nèi)部結(jié)構(gòu)即存在不同程度的缺陷與裂紋，盡管這些缺陷與裂紋在混合料內(nèi)部是離散分布的，但在連續(xù)損傷力學(xué)中，材料的這些缺陷可使用連續(xù)的內(nèi)部損傷場變量表示.在巖土及路面領(lǐng)域的研究中，Weibull分布函數(shù)被一些學(xué)者用以描述材料內(nèi)部的損傷過程，其表達式［13-14］為

其中:m，n為模型參數(shù)，γ為損傷時應(yīng)變門檻值.

假設(shè)損傷因子D變化率滿足Weibull分布，即dD/dε=φ(ε)，則損傷因子可寫為

在本文中，引入損傷力學(xué)中的有效應(yīng)力概念，即

1.2 粘彈性模型

國內(nèi)外一般使用簡單的機械模型組合來描述瀝青混合料的粘彈特性，如 Burgers、Kelvin、Maxwell和廣義Maxwell模型等.本文對多孔瀝青混合料進行單軸壓縮試驗，采用的是位移控制模式，即保持應(yīng)變速率不變，探討應(yīng)力隨時間與應(yīng)變的變化情況.已有研究表明，廣義Maxwell模型適用于描述恒應(yīng)變率加載下的材料粘彈性為.為簡化模型參數(shù)，選擇2個Maxwell模型與1個彈簧并聯(lián)，模型如圖1所示.

圖1 廣義Maxwell模型

該模型的連續(xù)方程可使用積分形式表示［15］為

式中:t1、t2為松弛時間，t1=λ1/E1，t2=λ2/E2;˙ε為應(yīng)變率;G 為松弛模量，MPa;E0、E1、E2為彈性模量分量;λ1、λ2為粘度分量.

積分型本構(gòu)方程當(dāng)進行數(shù)據(jù)擬合時，需對其離散化，即構(gòu)建增量型本構(gòu)方程，因此改寫式(4)為

其中:σ2為由λ1、E1組成的 Maxwell模型應(yīng)力，σ3為由λ2、E2組成的Maxwell模型應(yīng)力，σ為模型整體應(yīng)力.

1.3 塑性模型

在選擇塑性模型時，現(xiàn)有常用的屈服函數(shù)包括 Tresca、Von Mises、Mohr-Coulomb 與 Drucker-Prager屈服準(zhǔn)則等.其中Tresca屈服函數(shù)即為最大剪應(yīng)力理論，其函數(shù)曲線存在尖角，不利于數(shù)值計算，這也是Mohr-Coulomb屈服函數(shù)的問題.因此最終屈服函數(shù)選擇應(yīng)用廣泛的Drucker-Prager模型，該模型屈服函數(shù)在π平面的投影為光滑連續(xù)的圓形，且考慮了靜水壓力的影響.其屈服函數(shù)形式為［16］

其中:I1為第一主應(yīng)力張量不變量，J2為第二偏應(yīng)力張量不變量，k為各向同性硬化參數(shù)為累計塑性應(yīng)變，a為模型參數(shù).

在屈服面上，df=0，即上式可寫為

在選取塑性流動規(guī)則時，已有研究表明瀝青混合料的力學(xué)行為符合非關(guān)聯(lián)流動法則［17-19］，即塑性勢函數(shù)g并不等于屈服函數(shù)f.為了得到非關(guān)聯(lián)塑性模型，Drucker-Prager塑性勢函數(shù)采用與屈服函數(shù)類似的結(jié)構(gòu)形式，僅將其中的參數(shù)a變?yōu)閰?shù) b，有

因此，流動法則可表示為

式中:H為硬化模量，dλ為正比例標(biāo)量，δij為克羅內(nèi)克符號，Cijkl為彈性常數(shù)張量.

如前所述，為了更好地驗證模型，引入有效應(yīng)力與有效應(yīng)變的概念.對于本文的驗證試驗，加載模式為單軸壓縮.因此有效應(yīng)力ˉσe可表示為

而有效塑性應(yīng)變可用增量形式表達為

其中:G為剪切模量，MPa;K為體積彈性模量，MPa.

2 試 驗

2.1 試驗材料

試驗用集料為玄武巖，其粗集料技術(shù)指標(biāo)為:針片狀含量10.1%，小于0.075 mm的顆粒的質(zhì)量分數(shù)為 0.3%，壓碎值 10.0%，粘附性等級為4級.

瀝青為橡膠瀝青，其25℃針入度為62.3(0.1 mm)，軟化點 θR＆B為 63.6 ℃，RTFOT 后殘留針入度比為61.1%，180 ℃粘度為2.31 Pa·s，25℃彈性恢復(fù)86.3.

試驗多孔瀝青混合料參考公路行業(yè)規(guī)范推薦級配范圍，經(jīng)過對比體積指標(biāo)與性能要求確定級配［20］.其關(guān)鍵篩孔通過率如下:16 mm篩孔通過率為95%，4.75 mm篩孔通過率為 20%，2.36 mm篩孔通過率為15%，0.075 mm篩孔通過率為4%.瀝青用量為5.2%.

2.2 試驗方法

試驗用多孔瀝青混合料試樣高度與直徑均為100 mm，采用旋轉(zhuǎn)壓實儀直接成型，空隙率控制在18%±1%.加載采用單軸壓縮模式，試驗溫度分別為-10、0、10、20 ℃，在每種溫度下施加不同應(yīng)變率荷載，直至試樣完全破壞.具體加載速率如下:-10 ℃下為 2、1、0.5 mm/min;0 ℃下為 5、2、1 mm/min;10 ℃下為12.5、5、2 mm/min;20 ℃下為 50、5、2 mm/min.

試驗過程中采集荷載與位移，并轉(zhuǎn)化為應(yīng)力與應(yīng)變.每種試驗條件下制備4或5個平行試樣，在加載之前，對試樣施加小應(yīng)力預(yù)壓.試驗系統(tǒng)為MTS810型電液伺服試驗設(shè)備.在模型中，一共有12個模型參數(shù)需要確定，即損傷模型中的m、n、γ，粘彈性模型中的 E0、E1、E2、λ1、λ2，塑性模型中的K、G、a、b.由于粘彈性模型中松弛時間 t1=λ1/E1，因此粘彈性模型中的5個參數(shù)選用 E0、E1、E2、t1、t2.參數(shù)的確定過程包括首先根據(jù)經(jīng)驗對其擬定初始值，然后采用最小平方法進行非線性擬合，擬合工具采用Excel中規(guī)劃求解功能，設(shè)置迭代次數(shù)與允許誤差，控制模型擬合與試驗結(jié)果誤差平方和最小，以確定模型參數(shù).

3 試驗結(jié)果與模型分析

3.1 試驗結(jié)果與模型擬合

圖2為不同溫度下2 mm/min(0.000 3/s)加載速率時的試驗結(jié)果與模型擬合.可以看到對于不同溫度時的應(yīng)力應(yīng)變曲線，模型可很好地進行模擬.

圖2 不同溫度0.000 3/s加載速率下壓縮試驗結(jié)果與模型擬合

圖3、4分別為20℃與-10℃時的試驗結(jié)果與模型擬合.由于篇幅所限，本文并沒有給出10℃與0℃下的試驗數(shù)據(jù).但從4種溫度下的試驗數(shù)據(jù)與模型擬合結(jié)果看，模型在0℃以上與試驗數(shù)據(jù)擬合較好，當(dāng)溫度下降到-10℃時，由于瀝青混合料將顯示出更多的彈性，其力學(xué)性能對加載速率的依賴性減弱.因此，圖4中，無論是試驗結(jié)果還是模型擬合，多孔瀝青混合料的加載應(yīng)力應(yīng)變曲線都受加載速率影響較小.

圖3 20℃不同加載速率下單軸壓縮強度試驗結(jié)果與模型擬合

圖4 －10℃不同加載速率下單軸壓縮強度試驗結(jié)果與模型擬合

3.2 模型有效性分析

由圖2～4可知當(dāng)溫度較高或加載速率不大時，模型擬合情況較好，但是模型在描述多孔瀝青混合料在低溫下力學(xué)行為的有效性方面還需進一步分析.模型參數(shù)的取值和變化是影響模型精度的重要因素.參數(shù)的變化意味著混合料某些力學(xué)特性的改變，對于判斷材料屬性的變化規(guī)律有較大的幫助.

在前文的粘彈性模型中，假設(shè)在固定溫度下的粘彈屬性不變，加載速率的影響可以通過模型反映.最終得到溫度與各項參數(shù)之間的關(guān)系如圖5所示.可以看到隨著溫度的降低，E0、E1、E2、t1與t2都有不同程度的上升.但其中E1和E2當(dāng)溫度降低到0℃以下時變化趨勢減弱，表明粘彈性模型中粘彈構(gòu)件作用降低，材料開始顯示出更多的彈性屬性，也因此彈性模量盡管有所波動，但維持在同一水平.t1與t2表征了材料松弛應(yīng)力的能力，時間越短，表明材料具有更好松弛特性，也顯示出更明顯的粘彈特征.因此總的來看，材料參數(shù)的變化表明:溫度的下降使多孔瀝青混合料由粘彈性材料向彈性材料轉(zhuǎn)化，材料應(yīng)力松弛能力減弱.

圖5 粘彈參數(shù)隨溫度變化

在塑性模型中，經(jīng)典彈塑性理論假設(shè)剪切模量G和體積彈性模量K為常量，不隨溫度或加載速率的變化而改變.但是對于多孔瀝青混合料這種典型的粘彈性材料，其模量也應(yīng)具有溫度或頻率敏感性.圖6、7分別為G與K在溫度與加載速率變化時的取值.盡管試驗結(jié)果為離散的數(shù)據(jù)點，但通過數(shù)據(jù)回歸可以繪出其變化規(guī)律.

圖6 剪切模量G隨溫度與加載速率變化

圖7 體積彈性模量K隨溫度與加載速率變化

粘彈性材料的一個重要特征就是時溫等效性，即溫度變化導(dǎo)致的性能改變也可以通過改變加載時間(或加載頻率)得到.同時，粘彈性材料的力學(xué)特性(動態(tài)模量等)隨溫度或頻率的變化呈“S”形狀，即在固定溫度下存在一個敏感頻率范圍，在這個范圍內(nèi)材料屬性變化明顯，當(dāng)頻率過低或過高時，材料屬性隨頻率變化的幅度減小.從圖6、7可以看到，多孔瀝青混合料的剪切模量G和體積彈性模量K隨溫度或加載速率呈明顯的S型變化，證明了模型中G和K取值的合理性.

損傷模型中的γ為損傷時應(yīng)變門檻值，即當(dāng)應(yīng)變大于γ時，認為損傷開始發(fā)展.對于受荷變形的線粘彈性材料，一般認為只要時間足夠長其力學(xué)特性就可以恢復(fù).在塑性模型中，每一應(yīng)變增量后，通過屈服函數(shù)判斷是否出現(xiàn)塑性變形，并以此作為塑性應(yīng)變產(chǎn)生.本文認為塑性應(yīng)變的出現(xiàn)即意味著損傷的開始，而在初始線粘彈性階段，損傷等于0.為了驗證這個假設(shè)，對塑性應(yīng)變產(chǎn)生時對應(yīng)的總應(yīng)變與損傷時應(yīng)變門檻值γ進行分析，結(jié)果，見圖8.可以看到兩者之間有著很好的線性關(guān)系，圖中數(shù)據(jù)點包括不同溫度與不同加載速率下的結(jié)果，但是總的趨勢依然滿足所作的假設(shè).

圖8 塑性應(yīng)變產(chǎn)生γ時對應(yīng)的總應(yīng)變與損傷時應(yīng)變門檻值比較

對于典型的粘彈性材料，溫度和加載速率的不同，將使材料呈現(xiàn)出不同的損傷行為.低溫和高加載速率下材料顯示出彈脆性，而高溫和低加載速率下材料更傾向于粘塑性.對上述不同溫度及加載速率下多孔瀝青混合料的損傷特性進行分析，結(jié)果如圖9、10所示.為更好地比較溫度的影響，圖中僅給出20℃與-10℃下的曲線.

圖9 20℃3種加載速率下混合料損傷曲線

圖10 －10℃3種加載速率下混合料損傷曲線

可以看到，當(dāng)溫度為20℃時，0.000 33/s(即2 mm/min)的加載速率下多孔瀝青混合料顯示出更多的粘彈性，在這種情況下?lián)p傷從一開始加載就出現(xiàn)并緩慢增長，直到混合料不再能抵抗外荷載.而當(dāng)加載速率增大時，混合料逐漸顯示出明顯的彈脆性，即在一定彈性范圍內(nèi)變形可以恢復(fù)，損傷為0，但是一旦超過這個彈性范圍，損傷迅速發(fā)展，混合料很快出現(xiàn)脆性開裂破壞.

對于-10℃時的損傷曲線，由于溫度較低，瀝青混合料更趨向于顯示出彈性，此時混合料特性受加載速率的影響相對要小很多，因此，在不同加載速率下，損傷曲線基本相同.

對于彈脆性材料，直觀上看，材料越硬(即溫度越低，加載速率越大)，相對應(yīng)的損傷起始時間就越大，并且損傷擴展的速率也越快.圖11為4種溫度時0.000 33/s(即2 mm/min)加載速率下混合料的損傷曲線，即單純考慮溫度對損傷特性的影響.由于本文假設(shè)塑性變形的出現(xiàn)即為損傷的初始，因此溫度越低，出現(xiàn)塑性變形的時間將越長，導(dǎo)致初始損傷發(fā)生時的應(yīng)變也越大.此外，溫度的高低直接決定了損傷的擴展速率.低溫下的彈脆特性將帶來較快的損傷速率，而溫度較高時，材料的損失擴展行為更為緩慢.

圖11 4種溫度固定加載速率下混合料損傷曲線

4 結(jié) 論

1)不同溫度與加載速率下多孔瀝青混合料的應(yīng)力-應(yīng)變曲線顯示出明顯的粘彈效應(yīng).在高溫與低應(yīng)變率下，多孔瀝青混合料呈現(xiàn)出粘彈塑性破壞;而當(dāng)溫度較低或應(yīng)變率高時，其表現(xiàn)為彈脆性破壞.

2)構(gòu)建的粘彈塑性損傷模型可以有效反映多孔瀝青混合料在較低溫度下的力學(xué)行為.

3)由理論模型參數(shù)隨溫度及加載速率的變化規(guī)律可知，隨著溫度的降低，多孔瀝青混合料模量逐步增加，但是到一定值后升高幅度降低，同時材料應(yīng)力松弛能力減弱.

4)隨著溫度及加載速率的變化，多孔瀝青混合料剪切模量和體積彈性模量也呈現(xiàn)出“S”型變化.損傷曲線的變化規(guī)律也證明了構(gòu)建理論模型的合理性.

［1］ALVAREZ A E，MARTIN A E，ESTAKHRI C K，et al.Synthesis of current practice on the design，construction，and maintenance of porous friction course［R］.Texas:Texas Transportation Institute and The Texas A＆M University System，2006.

［2］曹東偉，劉清泉，唐國奇.排水瀝青路面［M］.北京:人民交通出版社，2010.

［3］MUNDEN H J，PIERCE L M，MAHONEY J P，et al.Low noise pavement evaluation in Washington State［C］//TRB Annual Meeting.Washington:Transportation Research Board，2008:1-14.

［4］POULIKAKOS L，TAKAHASHI S，PARTL M.A comparison of Swiss and Japanese porous asphalt through various mechanical tests［C］//3rd Swiss Transport Research Conference.Monte Verità/Ascona:Swiss Transport Research，2004:1-12.

［5］Danish Road Institute.Noise reducing pavements-state of the art in Denmark［R］.Denmark:Road Directorate，Ministry of Transport-Denmark，2005.

［6］YILDIRIM Y，DOSSEY T，F(xiàn)ULTS K，et al.Cold weather performance of new generation open graded friction courses［R］.Texas:Center for Transportation Research，The University of Texas at Austin，2007.

［7］諸永寧.排水性瀝青路面排水性能研究與排水設(shè)施的設(shè)計［D］.南京:東南大學(xué)，2004.

［8］SCHAUS L，TIGHE S，UZAROWSKI L.Porous asphalt pavement designs:Canadian climate use［C］//TRB Annual Meeting.Washington:Transportation Research Board，2008:1-21.

［9］MOHAN S.Winter damage of porous asphal-case study using a meso-mechanics based tool for lifetime optimization of PA［D］.Netherland:Delft University of Technology，2010.

［10］LEBENS M.Porous asphalt pavement performance in cold regions［R］.Minnesota:Minnesota Department of Transportation Research Services，2009.

［11］WANG L B，F(xiàn)ROST J D，VOYIADJIS G Z，et al.Quantification of damage parameters using X-ray tomography images［J］.Mechanics of Materials，2003，35(8):777-790.

［12］肖源杰，倪富健，蒯海東，等.基于圖像的粗集料形態(tài)對瀝青面層抗剪性能的影響［J］.鄭州大學(xué)學(xué)報:工學(xué)版，2006，27(4):44-48.

［13］鄭健龍，呂松濤，田小革.基于蠕變試驗的瀝青粘彈性損傷特性［J］.工程力學(xué)，2008，25(2):193-196.

［14］邵臘庚，周曉青，李宇峙，等.基于直接拉伸試驗的瀝青混合料粘彈性損傷特性研究［J］.土木工程學(xué)報，2005，38(4):125-128.

［15］易軍艷.基于界面行為的多孔瀝青混合料凍融損傷特性研究［D］.哈爾濱:哈爾濱工業(yè)大學(xué)，2012.

［16］CHEN W F，HAN D J.Plasticity for structural engineers［M］.FL Fort Lauderdale:J Ross Publishing Inc，2007.

［17］DARABI M K，ABU A R K，MASAD E A，et al.A thermo-viscoelastic-viscoplastic-viscodamage constitutive model for asphaltic materials［J］.International Journal of Solids and Structures，2011，48(1):191-207.

［18］FLOREA D.Nonassociated elastic/viscoplastic model for bituminous concrete［J］.International Journal of Engineering Science，1994，32(1):87-93.

［19］CHEHAB G.Characterization of asphalt concrete in tension using a viscoelastoplastic model［D］.North Carolina:North Carolina State University，2002.

［20］JTG F40—2004.公路瀝青路面施工技術(shù)規(guī)范［S］.北京:中華人民共和國交通部，2004.

Viscoelastic-plastic damage model of porous asphalt mixtures

YI Junyan1，SHEN Shihui2，MUHUNTHAN Balasingam2，F(xiàn)ENG Decheng1
(1.School of Transportation Science and Engineering，Harbin Institute of Technology，150090 Harbin，China;2.Department of Civil and Environmental Engineering，Washington State University，99164 WA，USA)

To characterize the mechanical behaviors of porous asphalt mixtures under loading at normal and low temperatures，a viscoelastic-plastic damage model，which includes Weibull damage function，Generalized Maxwell and Drucker-Prager model，was presented based on incremental constitutive equation.Experimental data from uniaxial compressive strength tests conducted at different strain rate and different temperature were used to validate the efficacy of the model，and the change rules of model parameters in different temperatures and loading rates were investigated.The test and analysis results show that parameters in viscoelastic model begin to behave like the elastic materials.The volume modulus and shear modulus also show obvious viscoelastic properties along with the changes of temperature and load rate.Additionally the strain plastic strain initializes keep almost same to the critical damage strain，which accommodate to the presented hypothesis.The damage curves of porous asphalt mixtures can also reveal the effects of different temperature and load rate.This model can successfully characterize the damage and mechanical behaviors of porous asphalt mixtures at normal and low temperatures.

road engineering;porous asphalt mixture;viscoelasticity;plasticity;damage

U416.2

A

0367-6234(2014)02-0066-06

2013-02-25.

交通運輸建設(shè)科技項目(2011318 801670).

易軍艷(1983—)，男，博士，講師;

馮德成(1968—)，男，教授，博士生導(dǎo)師.

易軍艷，yijunyan@hit.edu.cn.

(編輯 魏希柱)