凍融受荷協(xié)同作用下瀝青混合料損傷模型研究

張 倩，李 澤，溫志廣，侯振華

(1. 西安建筑科技大學(xué)土木工程學(xué)院，陜西 西安 710055；2. 包頭市公路局，內(nèi)蒙古 包頭 014030)

寒區(qū)瀝青混凝土路面常年處于正負(fù)溫交變狀態(tài)，在外荷載和凍融環(huán)境因素綜合作用下，內(nèi)部攜帶的各類初始微缺陷會(huì)不斷萌生、發(fā)展，當(dāng)累積損傷達(dá)到材料失效門檻值時(shí)，裂縫、松散和坑槽等宏觀物理性能逐漸響應(yīng)[1]．近年來，隨著寒區(qū)基礎(chǔ)設(shè)施工程建設(shè)不斷興起，越來越多的瀝青混凝土路面凍融耦合損傷問題亟待解決．因此，建立瀝青混合料凍融受荷協(xié)同損傷模型，探究其損傷劣化機(jī)制，進(jìn)而提出凍融損傷后壽命預(yù)估方法，將對凍區(qū)瀝青路面設(shè)計(jì)、性能評估及養(yǎng)護(hù)時(shí)機(jī)選擇具有重要的理論價(jià)值和現(xiàn)實(shí)意義．

近年來，國內(nèi)外眾多學(xué)者對瀝青混合料凍融破壞研究主要集中在環(huán)境或荷載等單因素對混合料宏觀性能指標(biāo)及凍融損傷特性的影響方面．AMINI等[1-3]研究了瀝青混合料含水量、飽和度、含鹽濃度、溫度等外在因素和級配、油石比、空隙率，對其強(qiáng)度和模量的耐凍性能影響，Si Wei[4]考慮了溫度和凍融循環(huán)對寒區(qū)瀝青混合料壓縮特性的影響，張倩等[5]基于凍融劈裂前后試件的 CT圖像，研究了瀝青混合料微觀結(jié)構(gòu)損傷特性，譚憶秋等[6-7]建立了瀝青混合料凍融損傷本構(gòu)方程，吳志勇[8]獲得了與加載次數(shù)，荷載大小及勁度模量相關(guān)的冪函數(shù)疲勞累積損傷模型．文獻(xiàn)[9-10]通過試驗(yàn)闡釋了瀝青混合料的性能凍融衰減特征及規(guī)律．事實(shí)上，寒區(qū)瀝青路面損傷是多因素全方位、多層次、多角度綜合作用的結(jié)果，不僅受限于材料內(nèi)部不同尺度的宏細(xì)觀影響，而且還有外部環(huán)境和荷載的作用．但現(xiàn)有大部分研究只考慮凍融或荷載等單因素作用下造成的損傷，不能真實(shí)反映瀝青路面實(shí)際工作狀態(tài)．嘗試在凍融循環(huán)條件下研究受荷瀝青混合料多因素的協(xié)同損傷問題，將是寒區(qū)瀝青路面研究的熱點(diǎn)課題．

基于損傷力學(xué)方法，借助宏觀統(tǒng)計(jì)損傷模型，構(gòu)建瀝青混合料凍融受荷協(xié)同作用下?lián)p傷模型，研究凍融次數(shù)與應(yīng)變對瀝青混合料力學(xué)特性及損傷劣化的影響，揭示材料凍融受荷損傷演變規(guī)律，探討不同油石比條件下凍融破壞機(jī)理，并通過凍融彎曲試驗(yàn)結(jié)果驗(yàn)證模型的合理性，以期為瀝青路面抗凍融受荷設(shè)計(jì)和性能評估提供理論支撐．

1 凍融受荷協(xié)同損傷模型

1.1 瀝青混合料凍融受荷協(xié)同破損機(jī)制

凍融與外荷載以不同的力學(xué)機(jī)理作用于瀝青混合料微裂隙等薄弱部位相連通的損傷域，耦合效應(yīng)加劇了空隙中水冰相變、水分遷移，使得荷載與材料損傷之間不斷反饋調(diào)整，導(dǎo)致材料物理力學(xué)性能劣化失效．

考慮到瀝青混合料凍結(jié)和融化的含冰混合料、水、空氣等多相介質(zhì)具有不同的熱物理特性和縮脹率，當(dāng)溫度降低時(shí)，飽水混合料體積收縮，而空隙水凝冰膨脹，易引起跨集料邊界縮脹失調(diào)；由于混合料的非連續(xù)特性，在集料及微空隙間產(chǎn)生顯著的凍脹力不能得到有效傳遞釋放而在粗集料界面附近易發(fā)生應(yīng)力集中，當(dāng)凍脹力超過瀝青-集料界面粘結(jié)力失效閾值時(shí)會(huì)迫使內(nèi)部微裂紋萌生、擴(kuò)展，損傷也相繼形成；當(dāng)溫度上升時(shí)，混合料內(nèi)部凝冰開始融化，同時(shí)伴隨著凍脹力的消散、遷移水對瀝青-集料界面反復(fù)侵蝕．多次凍融循環(huán)后，瀝青混合料累積耗散能逐漸降低，導(dǎo)致其物理力學(xué)性能不可逆衰減，殘余壽命相繼縮短．尤其在消融時(shí)段，外荷載會(huì)對混合料產(chǎn)生負(fù)壓泵吸作用，瀝青砂漿受空隙中正負(fù)動(dòng)水壓力反復(fù)沖刷效應(yīng)加劇，進(jìn)一步加速集料界面之間滑移與錯(cuò)動(dòng)，瀝青砂漿顆粒逐漸損失，集料表面瀝青膜減??；大量微裂隙、孔洞產(chǎn)生、擴(kuò)展、匯合貫通，材料性能迅速退化，最終演化為一定尺度的宏觀裂縫、剝落、松散、坑槽等病害．

1.2 凍融荷載協(xié)同損傷變量

瀝青混合料經(jīng)歷多次凍融后內(nèi)部細(xì)微觀缺陷損傷規(guī)律非常復(fù)雜，表征凍融損傷的彎曲勁度模量等物理力學(xué)指標(biāo)不斷退化．宏觀唯象損傷力學(xué)理論[11]認(rèn)為，瀝青混合料性能指標(biāo)彎曲勁度模量的衰減能夠反映材料的內(nèi)部損傷程度，且在凍融受荷試驗(yàn)過程中便于量測和分析．因此可將瀝青混合料凍融損傷變量表示為

式中：Dn為瀝青混合料n次凍融后的損傷變量；En、E0分別為瀝青混合料n次凍融后和凍融前初始彎曲勁度模量．

由于大多數(shù)瀝青混合料碾壓成型過程中皆攜帶初始微損傷，故將其凍融前初始損傷狀態(tài)定義為基準(zhǔn)損傷狀態(tài)，根據(jù)Lemaitre提出的應(yīng)變等效假設(shè)及有效應(yīng)力概念，可得瀝青混合料損傷本構(gòu)關(guān)系為

式中：Dl為荷載作用下?lián)p傷變量．

根據(jù)張全勝等[12]推廣后的應(yīng)變等價(jià)原理，材料在不同損傷狀態(tài)下的本構(gòu)關(guān)系形式相同，只替換式(2)中的E0即可．因此將瀝青混合料凍融損傷后的狀態(tài)作為第一種損傷狀態(tài)，凍融荷載協(xié)同作用下的損傷狀態(tài)作為第二種損傷狀態(tài)，從而可得損傷本構(gòu)關(guān)系為

將式(1)代入式(3)，得到瀝青混合料凍融荷載協(xié)同損傷本構(gòu)關(guān)系為

其中:

式中：Dt為凍融受荷總損傷變量，Dl Dn為耦合項(xiàng)．

由式(5)可以看出，總損傷變量并不等價(jià)于凍融變量與荷載變量的簡單代數(shù)疊加，二者耦合效應(yīng)加劇了瀝青混合料的損傷程度，呈非線性遞增趨勢．但耦合項(xiàng)前的負(fù)號表示在荷載下集料顆粒的搓揉、滑移削弱了由凍融引起的內(nèi)部損傷，凍脹效應(yīng)得到有效緩解，使總損傷有所弱化．

1.3 凍融荷載協(xié)同損傷演化方程

瀝青混合料是一種內(nèi)含各類隨機(jī)分布細(xì)觀缺陷的非均質(zhì)材料，初始攜帶的各類微空隙、微裂縫可視為連續(xù)分散于材料內(nèi)部的一種細(xì)觀結(jié)構(gòu)損傷場，其損傷程度在多因素綜合作用下也呈概率性分布，假設(shè)瀝青混合料微元體強(qiáng)度服從Weibull分布，其概率密度函數(shù)[13]為

式中，ε為瀝青混合料應(yīng)變值；m、F為表征材料物理力學(xué)性質(zhì)的參數(shù)，分別反映材料對外荷載不同響應(yīng)下的統(tǒng)計(jì)特征，通常與模量E、泊松比υ等材料參數(shù)有關(guān)．

假定瀝青混合料的損傷由局部微元體不均勻破壞引起，將某一荷載下已破壞微元體數(shù)目Nc與總微元體數(shù)目N之比定義為統(tǒng)計(jì)損傷變量Dl．在任意應(yīng)變區(qū)間內(nèi)，產(chǎn)生破壞的微元體數(shù)目為NP(x)dx，當(dāng)加載到某一應(yīng)變水平ε時(shí)，已破壞的微元體數(shù)目可表示為

聯(lián)立式(2)、(8)，根據(jù)應(yīng)力-應(yīng)變曲線上的峰值強(qiáng)度可推導(dǎo)確定損傷統(tǒng)計(jì)本構(gòu)模型參數(shù)m和F分別為

由式(1)、(5)、(8)、(9)和(10)，得出以凍融次數(shù)和應(yīng)變?yōu)閾p傷控制變量的總損傷演化方程為

式(11)表明，當(dāng)僅考慮凍融損傷時(shí)，即ε=0，則Dt=Dn；當(dāng)只考慮受荷損傷時(shí)，即E0=En，則Dt=Dl．

將式(11)代入式(4)得到瀝青混合料凍融受荷損傷本構(gòu)方程為

當(dāng)凍融次數(shù)和受荷應(yīng)變變化時(shí)，由式(1)、 (5)、(8)、(11)可得瀝青混合料總損傷率演化方程為

當(dāng)凍融次數(shù)恒定時(shí)，受荷損傷演化率為

當(dāng)受荷應(yīng)變恒定時(shí)，凍融損傷演化率為

由式(13)、(14)和(15)可知，瀝青混合料總損傷分別沿凍融和荷載兩途徑演化，揭示了凍融循環(huán)次數(shù)和應(yīng)變對其損傷演變相互耦合、相互影響的特性，客觀清晰地再現(xiàn)了瀝青混合料損傷力學(xué)行為及損傷破壞規(guī)律．

2 計(jì)算實(shí)例與模型驗(yàn)證

文獻(xiàn)[14]選取AC-13瀝青混合料，在不同凍融循環(huán)周期下進(jìn)行低溫彎曲試驗(yàn)．其中凍融試驗(yàn)條件為：?20 ℃低溫冷凍箱中凍結(jié) 12 h，25 ℃恒溫水浴中融化12 h，利用低溫冷凍箱與恒溫水浴模擬凍融作用，凍結(jié)時(shí)用塑料袋密封試件，并注水30 mL，使其充分飽水．低溫彎曲試驗(yàn)條件為：試驗(yàn)溫度為?10 ℃，加載速率為50 mm/min．表1為AC-13瀝青混合料級配組成，試驗(yàn)所得相關(guān)數(shù)據(jù)見表2～4，根據(jù)式(9)計(jì)算獲得模型參數(shù)m見表5，圖1是由式(11)計(jì)算得到的 5種油石比瀝青混合料凍融受荷損傷模型演化曲線．

表1 AC-13瀝青混合料級配Tab.1 Gradation of AC-13 asphalt mixture

表2 瀝青混合料凍融彎拉強(qiáng)度Tab.2 The tensile strength of asphalt mixture under freeze-thaw conditions

表3 瀝青混合料凍融破壞彎曲應(yīng)變Tab.3 The ultimate tensile strain of asphalt mixture under freeze-thaw conditions

表4 瀝青混合料凍融破壞勁度模量Tab.4 The limit stiffness modulus of asphalt mixture under freeze-thaw conditions

表5 凍融受荷損傷模型參數(shù)Tab.5 Parameters of damage model under freeze-thaw and loading

圖1 瀝青混合料凍融受荷損傷模型演化曲線Fig.1 Damage model evolution curves of asphalt mixture under freeze-thaw and loading

圖1計(jì)算結(jié)果表明：

不同油石比條件下，瀝青混合料的凍融損傷劣化隨著凍融循環(huán)次數(shù)的增加而加劇，但在經(jīng)歷20次凍融循環(huán)后，損傷變量總體變化幅度相對較小，基本趨于恒定，表明初期凍融循環(huán)對寒區(qū)瀝青混合料的抗彎拉特性影響較大．

在施荷初始階段，當(dāng)凍融次數(shù)相同時(shí)，隨著油石比增加，凍融損傷程度不斷得到弱化，在經(jīng)過20次凍融循環(huán)后，油石比達(dá)到6.5%時(shí)，損傷變量衰減超過50%，說明瀝青混合料凍融損傷程度受油石比影響顯著；然而，油石比的大小又與瀝青混合料的空隙率密切相關(guān)，由于凍融循環(huán)過程試件一直處于飽水狀態(tài)，當(dāng)油石比較低時(shí)，其內(nèi)部較大的空隙完全由冰水相控制，尤其凝冰會(huì)引起內(nèi)部空隙形成巨大的凍脹壓力，導(dǎo)致混合料損傷域不斷擴(kuò)展斷裂．隨著油石比的增大，其空隙率降低，內(nèi)部產(chǎn)生的膨脹力低于瀝青-集料粘結(jié)力失效閾值，瀝青混合料凍融損傷程度減緩并逐漸穩(wěn)定．

凍融循環(huán)次數(shù)不變時(shí)，不同油石比瀝青混合料凍融損傷劣化程度隨應(yīng)變的增大而加劇，這是由于受荷初期表現(xiàn)為損傷弱化階段，內(nèi)部微孔隙、微裂紋受壓閉合，凍融產(chǎn)生的空隙逐漸減少，整體呈“密實(shí)”狀態(tài)；當(dāng)應(yīng)變達(dá)到損傷閾值時(shí)，瀝青混合料開始損傷演化，不斷加速擴(kuò)展，期間伴隨著內(nèi)部微裂縫的萌生、擴(kuò)展、貫通，最終表現(xiàn)為宏觀裂紋的出現(xiàn)，到此損傷變量趨近1，材料產(chǎn)生破壞．且隨著油石比的加大，瀝青混合料壓密段逐漸變長，這與自由瀝青富余有關(guān)，多余的自由瀝青在集料之間發(fā)生相對滑動(dòng)過程中有潤滑作用，外荷載的施加，壓密階段持續(xù)時(shí)間相應(yīng)延長．

圖2 瀝青混合料沿受荷損傷演化率分布曲線Fig.2 Damage evolution rate distribution curves of asphalt mixture along loading

圖 2給出了根據(jù)式(14)計(jì)算得到的瀝青混合料沿受荷損傷演化率分布曲線．由圖可知，不同油石比下的瀝青混合料經(jīng)歷相同凍融次數(shù)時(shí)，隨著應(yīng)變的增加，損傷率呈非線性增大趨勢，達(dá)到峰值后逐漸減小．當(dāng)油石比大于5.5%(最佳油石比)時(shí)，存在應(yīng)變值臨界區(qū)間0.002～0.002 5，若應(yīng)變值小于下限，損傷率幾乎為零；若應(yīng)變值大于上限時(shí)，損傷率呈線性遞增趨勢．在最佳油石比條件下，不同凍融周期的瀝青混合料損傷率較小且演化幅度大致相同，說明按最佳油石比成型的路面可改善其凍融受荷耐久性能．

瀝青混合料沿受荷損傷率演化曲線特征基本一致，在峰值應(yīng)變處存在最大值，但在最佳油石比時(shí)損傷率相對較小，也進(jìn)一步佐證了在此狀態(tài)下瀝青混合料耐久性良好的觀點(diǎn)．相較于油石比為4.5%，損傷率演化曲線峰后斜率在6.5%時(shí)較小，表明油石比較小的瀝青混合料脆性破壞較顯著．因此，為避免寒區(qū)瀝青路面引起脆性失效，可在最佳油石比基礎(chǔ)上適當(dāng)上浮，可有效提高其抗凍融受荷破損的能力．

圖3 瀝青混合料沿凍融損傷演化率分布曲線Fig.3 Damage evolution rate distribution curves of asphalt mixture along freeze-thaw

根據(jù)式(14)和(15)計(jì)算得到瀝青混合料分別沿凍融和荷載兩條途徑損傷演化率曲線，如圖2和圖3所示．當(dāng)凍融次數(shù)一定時(shí)，受荷損傷率呈先增大，達(dá)到某一峰值后逐漸減小的趨勢；當(dāng)受荷應(yīng)變恒定時(shí)，凍融損傷演化率隨著凍融次數(shù)的增加逐漸減小，但在凍融次數(shù)超過 20時(shí)，無論沿凍融損傷還是沿受荷途徑損傷率降低幅度都較大且基本趨于恒定．由圖可知在油石比小于6.5%時(shí)，受荷應(yīng)變越大，凍融損傷率貢獻(xiàn)越小，反之亦然，說明凍融和荷載對瀝青混合料損傷率的貢獻(xiàn)變化是一個(gè)此消彼漲的過程，因此合理評估瀝青混合料的凍融受荷損傷劣化程度須在特定凍融周期和實(shí)際荷載工況下進(jìn)行．

圖4 彎拉強(qiáng)度與凍融次數(shù)關(guān)系Fig.4 Relation curves between tensile strength and freeze-thaw cycles

由于不同油石比下瀝青混合料凍融受荷損傷演化趨勢相同，為驗(yàn)證不同凍融次數(shù)條件下本構(gòu)方程的合理性，分別選取油石比為6.0%和6.5%的試驗(yàn)數(shù)據(jù)按式(12)進(jìn)行計(jì)算，圖 4給出了試驗(yàn)結(jié)果與模型計(jì)算結(jié)果．從圖4可看出，兩種油石比下的模型計(jì)算曲線與試驗(yàn)值均非常接近，吻合較好．

根據(jù)本文凍融受荷模型得到的瀝青混合料損傷演化規(guī)律與文獻(xiàn)[14]中試驗(yàn)現(xiàn)象及實(shí)際分析結(jié)論相一致．說明提出的表征瀝青混合料凍融受荷的損傷本構(gòu)方程和演化模型較為合理可靠．

3 結(jié)論

(1) 考慮凍融與荷載雙因素協(xié)同作用，基于損傷力學(xué)分析方法，建立了表征瀝青混合料凍融受荷損傷模型及演化方程，且驗(yàn)證了其合理性，為從細(xì)觀損傷方向研究宏觀凍融受荷協(xié)同破損機(jī)制提供了新視角．

(2) 凍融與荷載協(xié)同作用可使瀝青混合料總損傷程度加劇，呈非線性增長趨勢，但耦合效應(yīng)使總損傷弱化．二者以不同的損傷機(jī)理作用于瀝青混合料，使得兩種損傷之間不斷反饋調(diào)整，最終導(dǎo)致材料發(fā)生破壞．在油石比恒定條件下，瀝青混合料凍融損傷變量隨凍融次數(shù)的增加而增大，但在 20次之后，其值基本恒定；在凍融損傷變量不變情況下，隨著凍融次數(shù)的增加，應(yīng)變呈降低趨勢；當(dāng)應(yīng)變一定時(shí)，凍融損傷劣化程度隨凍融次數(shù)增加而加?。?/p>

(3) 瀝青混合料沿凍融和荷載兩條途徑損傷演化率曲線呈不同分布特征．沿受荷應(yīng)變損傷時(shí)，損傷率呈先增大后減小趨勢，存在峰值；沿凍融循環(huán)損傷時(shí)，損傷率逐漸減小并最終保持穩(wěn)定，但在兩種途徑下?lián)p傷率在超過 20次凍融循環(huán)后降低幅度都較大且趨于恒定．

References

[1]AMINI B, TEHRANI S S. Simultaneous effects of salted water flow on asphalt concrete pavement deterioration under freeze-thaw cycles [J]. International Journal of Pavement Engineering, 2014, 15(5): 383-391.

[2]SéBASTIEN L, DANIEL P, HERVé D B. Contraction and expansion of partially saturated hot mix asphalt samples exposed to freeze-thaw cycles[J]. Road Materials and Pavement Design, 2015, 16(2): 277-299.

[3]OZGAN E, SERIN S. Investigation of certain engineering characteristics of asphalt concrete exposed to freeze-thaw cycles [J]. Cold Reg. Sci. Technol., 2013（85）: 131-136.

[4]SI Wei, MA Biao, WANG Hainian, et al. Analysis of compressive characteristics of asphalt mixture under freeze-thaw cycles in cold plateau regions[J]. Journal of Highway and Transportation Research and Development,2013, 7(4): 17-22.

[5]張倩, 李創(chuàng)軍. 瀝青混合料凍融劈裂微觀結(jié)構(gòu)損傷特性分析[J]. 公路交通科技, 2010, 27(2): 6-9.ZHANG Qian, LI Chuangjun. Analysis of micro structural damage characteristic of freeze thaw split asphalt mixtures[J]. Journal of Highway and Transportation Research and Development, 2010, 27(2): 6-9.

[6]譚憶秋, 趙立東, 藍(lán)碧武, 等. 瀝青混合料凍融損傷模型及壽命預(yù)估研究[J]. 公路交通科技, 2011,28(6): 1-7.TAN Yiqiu, ZHAO Lidong, LAN Biwu, et al．Research on freeze-thaw damage model and life prediction of asphalt mixture[J].Journal of Highway and Transportation Research and Development, 2011, 28(6): 1-7.

[7]李兆生, 譚憶秋. 瀝青混合料凍融損傷特性研究[J].中國科技論文, 2014, 9(11): 1279-1281.LI Zhaosheng, TAN Yiqiu. Research on asphalt mixture freezing-thawing damage performance[J]. China Science Paper. 2014, 9(11): 1279-1281.

[8]吳志勇. 基于多級等幅荷載下的瀝青混合料損傷累積和瀝青面層疲勞損傷破壞研究[D]. 廣州: 華南理工大學(xué), 2014.WU Zhiyong. Research on cumulative fatigue damage of asphalt mixture and asphalt layer based on multi-level amplitude loading[D].Guangzhou: South China University of Technology, 2014.

[9]司偉, 馬骉, 任君平, 等. 基于可靠度的凍融循環(huán)下瀝青路面性能[J].吉林大學(xué)學(xué)報(bào)(工學(xué)版), 2016, 46(1):126-132.SI Wei, MA Biao, REN Junping, et al. Analysis of asphalt pavement performance under freeze-thaw cycles using reliability method[J].Journal of Jilin University (Engineering and Technology Edition), 2016, 46(1):126-132.

[10]顏可珍, 游凌云, 王曉亮. SMA 混合料凍融循環(huán)性能試驗(yàn)及超聲波法評價(jià)[J].中國公路學(xué)報(bào),2015,28(11):8-14.YAN Kezhen, YOU Linyun, WANG Xiaoliang. Test on characteristics of SMA mixture under freeze-thaw cycle and ultrasonic evaluation method[J]. China Journal of Highway and Transport, 2015, 28(11): 8-14.

[11]李兆霞. 損傷力學(xué)及其應(yīng)用[M]. 北京: 科學(xué)出版社,2002.LI Zhaoxia. Damage mechanics and its application[M].Beijing: Science Press, 2002.

[12]張全勝, 楊更社, 任建喜. 巖石損傷變量及本構(gòu)方程的新探討[J]. 巖石力學(xué)與工程學(xué)報(bào), 2003, 22(1):30-34.ZHANG Quansheng, YANG Gengshe, REN Jianxi. New study of damage variable and constitutive equation of rock[J].Chinese Journal of Rock Mechanics and Engineering, 2003, 22(1): 30-34.

[13]張我華, 孫林柱, 王軍, 等. 隨機(jī)損傷力學(xué)與模糊隨機(jī)有限元[M]. 北京: 科學(xué)出版社, 2011.ZHANG Wohua, SUN Linzhu, WANG Jun, et al. The stochastic damage mechanics and fuzzy-random finite element method[M]. Beijing: Science Press, 2011.

[14]虎見. 凍融循環(huán)作用下的瀝青混合料力學(xué)參數(shù)研究[D].西安: 長安大學(xué), 2011.HU Jian. Study on the mechanical parameters of asphalt mixture under the action of freeze-thaw cycling [D].Xi?an: Chang?an University, 2011.