江祥林 錢振東 宋 鑫
(1長安大學(xué)公路學(xué)院,西安 610100)(2江西省橋梁檢(監(jiān))測及加固重點(diǎn)實(shí)驗(yàn)室,南昌 330038)(3東南大學(xué)教育部智能運(yùn)輸系統(tǒng)工程研究中心,南京 210096)
基于離散單元法的環(huán)氧瀝青混凝土虛擬斷裂試驗(yàn)研究
江祥林1,2錢振東3宋 鑫3
(1長安大學(xué)公路學(xué)院,西安 610100)
(2江西省橋梁檢(監(jiān))測及加固重點(diǎn)實(shí)驗(yàn)室,南昌 330038)
(3東南大學(xué)教育部智能運(yùn)輸系統(tǒng)工程研究中心,南京 210096)
摘 要:為明晰環(huán)氧瀝青混凝土的斷裂力學(xué)行為特征,基于數(shù)字圖像處理技術(shù),采用離散元方法建立了切口小梁二維模型,開展了環(huán)氧瀝青混凝土的虛擬彎曲斷裂試驗(yàn).從細(xì)觀角度分析了環(huán)氧瀝青混凝土斷裂過程中的力學(xué)響應(yīng),探討了其斷裂機(jī)理及裂紋擴(kuò)展路徑,并與室內(nèi)試驗(yàn)結(jié)果進(jìn)行了對比.結(jié)果表明,采用離散元方法可較好地反映環(huán)氧瀝青混凝土的斷裂力學(xué)特性,數(shù)值模擬過程中虛擬試件的力學(xué)響應(yīng)與理論結(jié)果相符.虛擬試驗(yàn)獲取的材料強(qiáng)度參數(shù)與室內(nèi)試驗(yàn)的測試結(jié)果接近,誤差僅為0.25%,但前者所得的勁度模量因破壞撓度偏大,較后者低16.56%.在裂紋擴(kuò)展過程中,破壞易產(chǎn)生于黏結(jié)較為薄弱的集料-砂漿界面.虛擬試驗(yàn)與室內(nèi)試驗(yàn)的對比分析結(jié)果驗(yàn)證了離散元模型以及參數(shù)取值的正確性.
關(guān)鍵詞:環(huán)氧瀝青混凝土;離散單元法;斷裂模型;虛擬試驗(yàn)
裂縫是鋼橋面用環(huán)氧瀝青混凝土(EAC)的主要病害形式.目前,針對EAC鋪裝裂縫行為的研究多偏重于宏觀尺度下鋪面材料的力學(xué)及使用性能研究;但作為一種準(zhǔn)脆性、非均質(zhì)復(fù)合材料,EAC的破壞行為特性很大程度上取決于細(xì)觀層次的結(jié)構(gòu)特性.因此,如何研究細(xì)觀尺度下EAC的裂縫失效行為、提出針對性防治技術(shù),是研究人員當(dāng)前關(guān)注的熱點(diǎn)問題.
近年來,離散元法在處理非連續(xù)介質(zhì)問題上得到了較大的發(fā)展和應(yīng)用.Abbas[1]采用PFC2D軟件模擬了瀝青混合料的簡單性能試驗(yàn)和劈裂試驗(yàn),研究了瀝青混合料高溫條件下的黏彈性響應(yīng)和低溫條件下的開裂原理.Kim等[2-3]首次建立了內(nèi)嵌雙線性內(nèi)聚力接觸本構(gòu)的瀝青混合料圓盤緊湊拉伸(DC(T))離散元模型,研究了加載速率[4]和尺寸效應(yīng)[5]等對混合料斷裂行為的影響.楊軍等[6]為了更好地評價(jià)和預(yù)估瀝青混合料高溫穩(wěn)定性能,對瀝青混合料的三軸剪切試驗(yàn)進(jìn)行了離散元模擬.張德育等[7]利用瀝青混合料離散元虛擬單軸壓縮蠕變試驗(yàn),從材料不連續(xù)的角度研究了瀝青混合料的永久變形能力.
本文采用離散元軟件PFC2D,基于圖像采集處理技術(shù),建立了低溫條件下的EAC二維離散元斷裂模型,從細(xì)觀尺度上探求了EAC的斷裂行為特征.研究結(jié)果有助于指導(dǎo)材料的抗裂設(shè)計(jì)與養(yǎng)護(hù)修復(fù)方案的制定.
低溫條件下的環(huán)氧瀝青混合料可視為線彈性材料.在充分考慮各單元間的接觸關(guān)系和PFC2D軟件提供的接觸模型特點(diǎn)后,選擇如表1所示的接觸模型來描述混合料內(nèi)部各相材料單元間的接觸關(guān)系.
表1 低溫條件下環(huán)氧瀝青混合料離散元模型內(nèi)單元接觸類型
接觸剛度與材料宏觀模量間具有如下的轉(zhuǎn)換關(guān)系[8]:
式中,kn,ks分別為單元間的法向、切向接觸剛度;Ec為單元間的接觸楊氏模量;ν為泊松比;t為單元厚度.
鋼橋面鋪裝采用玄武巖集料,取集料楊氏模量為 56.8 GPa,泊松比為 0.2[7-11].環(huán)氧瀝青砂漿的楊氏模量由抗壓回彈模量試驗(yàn)確定為7.2 GPa,泊松比為0.25.細(xì)觀參數(shù)可由式(1)和(2)計(jì)算確定.玄武巖集料的黏結(jié)強(qiáng)度取為24 MPa[9-10];環(huán)氧瀝青砂漿的黏結(jié)強(qiáng)度由劈裂試驗(yàn)獲得,環(huán)境溫度為-10℃的條件下黏結(jié)強(qiáng)度為12.13 MPa.目前,對于界面黏結(jié)強(qiáng)度,尚缺乏成熟的測試方法,本文采用文獻(xiàn)[10]中反復(fù)試算與真實(shí)物理試驗(yàn)結(jié)果對比的方法來校核估定.細(xì)觀參數(shù)可根據(jù)接觸黏結(jié)模型參數(shù)與宏觀材料強(qiáng)度間的關(guān)系式計(jì)算得到,即
式中,Sn,Ss分別為承受的抗拉力和抗剪力;σc,τc分別為材料的極限抗拉強(qiáng)度和抗剪強(qiáng)度;ˉR為相互接觸的2個(gè)顆粒半徑的平均值.
在確定平行黏結(jié)模型的參數(shù)時(shí),需確定法向剛度ˉkn、切向剛度ˉkn、法向抗拉強(qiáng)度ˉσc、抗剪強(qiáng)度ˉτc和平行黏結(jié)半徑ˉR五個(gè)參數(shù).考慮到小梁斷裂過程中主要承受彎拉(壓)應(yīng)力,而且黏結(jié)接觸模型中小梁在受到彎曲作用時(shí)才會同時(shí)傳遞力和彎矩,因此平行黏結(jié)模型中環(huán)氧瀝青砂漿的法向強(qiáng)度和剛度由小梁彎曲試驗(yàn)得到,環(huán)境溫度為-10℃時(shí)分別取為38.54和5 407 MPa.在平行黏結(jié)模型中,強(qiáng)度與材料宏觀彎拉強(qiáng)度呈對應(yīng)關(guān)系,無需轉(zhuǎn)化,其他模型參數(shù)可由下式得到[8]:
式中,c為平行黏結(jié)的楊氏模量;n為平行黏結(jié)的法向與切向剛度比,此處取為 1.3[5].
EAC小梁經(jīng)切割處理后,采用高精度數(shù)碼設(shè)備采集其側(cè)面圖像(見圖1(a)).小梁長250 mm,高35 mm.將采集到的截面圖轉(zhuǎn)化為灰度圖,進(jìn)行二值化處理.二值化后的邊緣圖像由像素1和0組成,分別表示集料和砂漿.將粗細(xì)集料的劃分粒徑設(shè)定為2.36mm,采用 Image-Pro Plus 6.0軟件對細(xì)集料進(jìn)行濾除處理,將其劃歸環(huán)氧瀝青砂漿部分,處理之后的二維數(shù)字試件如圖1(b)所示.
圖1 EAC小梁截面原始圖與灰度圖
采用離散元軟件PFC2D,建立EAC小梁離散元模型.經(jīng)過多次試算,并綜合國內(nèi)外研究經(jīng)驗(yàn)[11-12],確定顆粒單元半徑為 0.5 mm,厚度為30 mm,共劃分為8 750個(gè)單元,并將確定的細(xì)觀參數(shù)賦予各單元接觸模型.由于EAC的空隙率小(一般為1.5% ~3.0%),在進(jìn)行圖像處理時(shí)識別空隙困難,因此從環(huán)氧瀝青砂漿離散元模型內(nèi)隨機(jī)刪除2.0%的單元作為空隙,以提高模擬的真實(shí)性和準(zhǔn)確度.
鋼橋面瀝青混凝土鋪裝層裂縫類病害以Ⅰ型張開開裂為主.考慮到切口小梁彎曲斷裂過程簡單,受力明確,故在小梁跨中下方刪除部分單元,形成高7 mm、寬2 mm的裂縫.而后,在梁底距跨中左右各100 mm位置處施加豎向約束.跨中上方施加向下的強(qiáng)制位移,加載寬度為10 mm,速率為1 mm/min,得到如圖2所示的離散元模型.離散元模型集料的分布和形狀特征與原試件保持一致.
圖2 環(huán)氧瀝青混凝土切口小梁離散元模型(單位:mm)
采用離散元小梁模型,對EAC切口小梁彎曲試驗(yàn)進(jìn)行仿真.根據(jù)模型計(jì)算結(jié)果,分別提取時(shí)間與荷載、加載點(diǎn)位移(即跨中撓度)間的對應(yīng)數(shù)據(jù)關(guān)系.將荷載-跨中撓度曲線與室內(nèi)試驗(yàn)結(jié)果進(jìn)行對比,結(jié)果見圖3.
圖3 荷載-跨中撓度曲線
由圖3可知,虛擬試驗(yàn)結(jié)果曲線與室內(nèi)試驗(yàn)結(jié)果曲線存在一定差異.低溫條件下忽略了環(huán)氧瀝青砂漿的黏彈性,視其為線彈性材料,故在加載階段,數(shù)值模擬曲線為直線.加載至約1.44 kN時(shí),曲線發(fā)生突變,荷載出現(xiàn)最大值.結(jié)合圖2可知,預(yù)制裂縫上方存在粗集料,由于集料的抗拉強(qiáng)度大于環(huán)氧瀝青砂漿,因此出現(xiàn)了斷裂增韌現(xiàn)象.室內(nèi)試驗(yàn)中,在初始加載階段,跨中撓度隨荷載作用呈非線性增長,而后隨著撓度的增大,荷載線性增加,材料表現(xiàn)為彈性.當(dāng)小梁達(dá)到破壞極限狀態(tài)后,圖3中的2條荷載-跨中撓度曲線均驟然回落.室內(nèi)試驗(yàn)中,試件發(fā)生結(jié)構(gòu)性破壞后,傳感器自動判斷出試件不再具有承載能力,曲線自然平滑趨向于橫軸;而在虛擬試驗(yàn)中,斷裂路徑上單元顆粒間的接觸模型需完全破壞,荷載才會降為0,因此,尚有顆粒單元承擔(dān)部分荷載作用,曲線出現(xiàn)了一定的起伏.
分別讀取虛擬試驗(yàn)與室內(nèi)試驗(yàn)中荷載達(dá)到峰值的時(shí)間T、破壞荷載PB、破壞時(shí)跨中撓度d等數(shù)據(jù),計(jì)算出彎拉強(qiáng)度RB、勁度模量SB等彎曲破壞評價(jià)指標(biāo)以及應(yīng)力強(qiáng)度因子KIC、斷裂能GF等斷裂參數(shù),結(jié)果見表2.
表2 虛擬試驗(yàn)結(jié)果與室內(nèi)試驗(yàn)結(jié)果的對比
由表2可知,在虛擬試驗(yàn)中,經(jīng)歷33.01 s的位移加載后,小梁所受荷載達(dá)到最大值1 617.10 N,跨中撓度為0.55 mm;在室內(nèi)試驗(yàn)中,到達(dá)極值點(diǎn)的荷載時(shí)間為26.78 s,相對應(yīng)的荷載和跨中撓度分別為1 621.15 N和0.46 mm.實(shí)際試驗(yàn)中,到達(dá)峰值荷載的時(shí)間取決于跨中三維斷裂區(qū)的抗彎能力,而不僅是二維區(qū)域.圖4所示的破壞斷面中存在其他粗集料,因此在加載階段,相比虛擬試驗(yàn),室內(nèi)試驗(yàn)中的荷載增長速率較大,小梁的破壞時(shí)間與撓度較小.通過對比可知,虛擬試驗(yàn)獲取的材料強(qiáng)度參數(shù)與室內(nèi)試驗(yàn)的測試結(jié)果接近,誤差僅為0.25%,由此驗(yàn)證了材料參數(shù)取值的合理性.由于應(yīng)力強(qiáng)度因子與臨界荷載成正比例關(guān)系[13],故兩計(jì)算結(jié)果差異較小.虛擬試驗(yàn)中的破壞模量偏大,因此計(jì)算所得的勁度模量較試驗(yàn)結(jié)果低16.56%.對比圖2中撓度曲線的區(qū)域面積,由于虛擬試驗(yàn)所得曲線回落緩慢,故其斷裂能相比于室內(nèi)試驗(yàn)結(jié)果至少提高了近73%.
圖4 EAC切口小梁斷面圖
在模擬過程中選取典型的時(shí)間節(jié)點(diǎn),分析模型內(nèi)部結(jié)構(gòu)的力學(xué)響應(yīng).圖5為裂縫萌生、裂縫啟裂、裂縫擴(kuò)展、接近完全斷裂時(shí)節(jié)點(diǎn)的小梁應(yīng)力分布和單元位移圖.
圖5 斷裂模擬過程的力學(xué)響應(yīng)分析
由圖5可知,斷裂過程中,裂縫尖端始終受拉應(yīng)力的作用.當(dāng)荷載尚未超過最大值時(shí),小梁荷載作用下,上半?yún)^(qū)域受壓應(yīng)力,下半?yún)^(qū)域受拉應(yīng)力,且在有效跨徑范圍內(nèi)形成主應(yīng)力軌跡線(見圖5(a)和(b)).對比不同階段的小梁應(yīng)力分布圖可知,隨著裂紋的發(fā)展,拉應(yīng)力分布范圍逐漸向上移動;當(dāng)荷載超過峰值后,小梁內(nèi)部拉、壓應(yīng)力范圍和數(shù)值均呈減小趨勢,切口裂縫不斷上升,裂紋嘴張開位移明顯增大.綜上可知,在斷裂模擬過程中,小梁內(nèi)部的應(yīng)力分布變化與相關(guān)力學(xué)理論相符,由此表明,建立的離散元模型和選擇的材料參數(shù)是正確有效的.
為了深入認(rèn)識切口小梁的斷裂路徑和機(jī)理,截取3.1節(jié)中典型時(shí)間節(jié)點(diǎn)相對應(yīng)的跨中斷裂過程區(qū)裂紋發(fā)展情況(見圖6).
圖6 細(xì)觀尺度下的裂紋發(fā)展
圖6(a)為荷載增加到1 538.60 N時(shí)的裂紋情況.由于集料的剛度較大,裂紋首先出現(xiàn)于集料內(nèi)部單元之間.當(dāng)荷載達(dá)到峰值時(shí),預(yù)制裂縫上端的粗集料已接近于“穿心”破壞,此刻預(yù)制裂縫右上方集料與砂漿之間也發(fā)生了黏結(jié)破壞,但并沒有與粗集料裂紋相連(見圖6(b)).當(dāng)跨中撓度接近于0.778 mm時(shí),裂紋出現(xiàn)較大擴(kuò)展,貫穿粗集料并近似沿著集料-砂漿界面上升(見圖6(c)).當(dāng)試件接近于完全斷裂時(shí),小梁跨中上部區(qū)域裂縫演化迅速,多個(gè)砂漿單元、界面單元之間均發(fā)生了破壞(見圖6(d)).縱觀小梁斷裂模擬的整個(gè)過程,當(dāng)同時(shí)存在砂漿內(nèi)部單元與界面單元時(shí),裂紋多傾向于出現(xiàn)在界面單元之間,這與界面黏結(jié)強(qiáng)度較小有關(guān).
圖7為室內(nèi)試驗(yàn)中觀測到的小梁裂紋.將其與虛擬試驗(yàn)得到的裂紋進(jìn)行對比,發(fā)現(xiàn)這2條裂紋的擴(kuò)展路徑相似度較高,均經(jīng)過預(yù)制裂紋上方的粗集料,且穿透的位置幾乎一致,在裂紋繼續(xù)發(fā)展演化的過程中,也均繞開細(xì)集料,沿其邊緣向上發(fā)展.2條裂紋在起裂位置以及后期發(fā)展路徑上有所不同,這可能是由于二維模型的局限性所致.由此再次證明了本文對環(huán)氧瀝青混合料切口小梁細(xì)觀斷裂行為的數(shù)值模擬是有效的,可較為準(zhǔn)確地反映低溫條件下切口小梁在集中荷載作用下的斷裂過程.
圖7 室內(nèi)試驗(yàn)中小梁裂縫
1)基于采集的環(huán)氧瀝青混凝土小梁截面圖像,應(yīng)用圖像處理技術(shù),并編制顆粒流命令,生成了環(huán)氧瀝青混凝土切口小梁二維數(shù)字試件.
2)在環(huán)境溫度為-10℃的條件下,環(huán)氧瀝青混凝土呈脆性破壞.虛擬實(shí)驗(yàn)中,裂紋首先出現(xiàn)在粗集料中,而后貫穿預(yù)制裂縫向上發(fā)展;在裂紋擴(kuò)展過程中,破壞易產(chǎn)生于集料-砂漿界面.
3)室內(nèi)試驗(yàn)結(jié)果及理論分析結(jié)果均表明,離散元模擬結(jié)果合理,與室內(nèi)試驗(yàn)結(jié)果相關(guān)性較高、規(guī)律相近,由此驗(yàn)證了離散元模型以及參數(shù)取值的正確性.
4)利用PFC2D軟件模擬了環(huán)氧瀝青混合料切口小梁的斷裂過程.然而,實(shí)際的瀝青混合料屬于三維空間體結(jié)構(gòu),建立的二維模型仍無法突破這一限制,因此,仿真結(jié)果與室內(nèi)試驗(yàn)結(jié)果仍存在一定的差異.
[1]Abbas A R.Simulation of the micromechanical behavior of asphalt mixtures using the discrete element method[D].Washington DC:Washington State University,2004.
[2]Kim H,Buttlar W G.Micromechanical fracture modeling of asphalt mixture using the discrete element method[C]//Advances in Pavement Engineering(GSP130).Austin,TX,USA,2005:209-223.
[3]Kim H,Wagoner M P,Buttlar W G.Simulation of fracture behavior in asphalt concrete using a heterogeneous cohesive zone discrete element model[J].Journal of Materials in Civil Engineering,2008,20(8):552-563.
[4]Kim H,Wagoner M P,Buttlar W G.Rate-dependent fracture modeling of asphalt concrete using the discrete element method[J].Canadian Journal of Civil Engineering,2009,36(2):320-330.
[5]Kim H,Wagoner M P,Buttlar W G.Numerical fracture analysis on the specimen size dependency of asphalt concrete using a cohesive softening model[J].Construction and Building Materials,2009,23(5):2112-2120.
[6]楊軍,張旭,朱浩然.瀝青混合料三軸剪切試驗(yàn)的離散元模擬研究[J].建筑材料學(xué)報(bào),2012,15(1):64-68.
Yang Jun,Zhang Xu,Zhu Haoran.Discrete element simulation on tri-axial shear test of asphalt mixtures[J].Journal of Building Materials,2012,15(1):64-68.(in Chinese)
[7]張德育,黃曉明,高英.瀝青混合料三維離散元虛擬單軸蠕變試驗(yàn)[J].華南理工大學(xué)學(xué)報(bào):自然科學(xué)版,2010,40(7):15-20.
Zhang Deyu,Huang Xiaoming,Gao Ying.Three-dimension virtual uniaxial creep test of asphalt mixture by using discrete element method[J].Journal of South China University of Technology:Natural Science Edition,2010,40(7):15-20.(in Chinese)
[8]Itasca Consulting Group.PFC2D version3.1[M].Minneapolis,MN, USA:ItascaConsultingGroup Inc.,2004:101-112.
[9]Dai Q L.Micromechanical modeling of constitutive and damage behavior ofheterogeneous asphalt materials[D].Kingston,RI,USA:University of Rhode Island,2004.
[10]You Z P,Adhikari S,Dai Q L.Three-dimensional discrete element models for asphalt mixtures[J].Journal of Engineering Mechanics,2008,134(12):1053-1062.
[11]陳俊.基于離散元方法的瀝青混合料虛擬疲勞試驗(yàn)研究[D].南京:東南大學(xué)交通學(xué)院,2010.
[12]Kim H,Buttlar W G.Discrete fracture modeling of asphalt concrete[J].International Journal of Solids and Structures,2009,46(13):2593-2604.
[13]程靳,趙數(shù)山.斷裂力學(xué)[M].北京:科學(xué)出版社,2006:56-79.
Experimental study on virtual fracture of epoxy asphalt concrete based on discrete element method
Jiang Xianglin1,2Qian Zhendong3Song Xin3
(1School of Highway,Chang'an University,Xi'an 610100,China)
(2Jiangxi Key Laboratory of Bridge Test and Reinforcement,Nanchang 330038,China)
(3Intelligent Transportation System Institute,Southeast University,Nanjing 210096,China)
Abstract:In order to clarify the fracture mechanical behaviors of epoxy asphalt concrete(EAC),a single-notched beam two-dimensional model was established by the discrete element method based on digital image processing technology.The virtual bending fracture tests of EAC were performed to study the mechanical response in the fracture process of EAC under meso-scale.The fracture mechanism and cracks propagation paths were analyzed and compared with the laboratory test results.The results show that the discrete element method can be used to exhibit the fracture characteristics of EAC excellently.The mechanical responses of the virtual tests obtained by numerical simulation are coincident with the theoretical results.The material strength parameters acquired by the virtual test are close to those of the indoor tests,with an error of 0.25%.However,the stiffness modulus obtained by the virtual test is 16.56%lower than the test results due to the large deflection.The damage tends to occur in the interface with inferior adhesive performance between the aggregate and asphalt mastic during the cracks propagation.The comparison between the virtual test results and the experimental results verify the availability of the discrete element model and the relevant material parameters.
Key words:epoxy asphalt concrete;discrete element method;fracture model;virtual test
中圖分類號:U416.2
A
1001-0505(2014)01-0173-05
doi:10.3969/j.issn.1001 -0505.2014.01.031
收稿日期:2013-07-12.
江祥林(1975—),男,博士生;錢振東(聯(lián)系人),女,博士,教授,博士生導(dǎo)師,qianzd@seu.edu.cn.
基金項(xiàng)目:國家自然科學(xué)基金資助項(xiàng)目(51178114)、江西省交通運(yùn)輸廳科技資助項(xiàng)目(2010C00005).
江祥林,錢振東,宋鑫.基于離散單元法的環(huán)氧瀝青混凝土虛擬斷裂試驗(yàn)研究[J].東南大學(xué)學(xué)報(bào):自然科學(xué)版,2014,44(1):173-177.[doi:10.3969/j.issn.1001 -0505.2014.01.031]