崔 溦 魏 杰 李國(guó)棟

(1天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室， 天津 300350)(2天津大學(xué)中國(guó)地震局地震工程綜合模擬與城鄉(xiāng)抗震韌性重點(diǎn)實(shí)驗(yàn)室， 天津 300350)(3中國(guó)電建市政建設(shè)集團(tuán)有限公司， 天津 300384)

混凝土具有原材料豐富、抗水性好、抗壓強(qiáng)度高、成本低等優(yōu)點(diǎn)，被廣泛應(yīng)用于土木、水利、橋梁等建筑工程中，是當(dāng)代應(yīng)用最廣泛的人工建筑材料之一[1].在細(xì)觀層次上，混凝土被視為由砂漿、粗骨料及二者之間的界面過(guò)渡區(qū)(ITZ)組成的三相復(fù)合材料.粗骨料的體積占混凝土試件總體積的 40%～50%，對(duì)混凝土的彈性模量、密度等起著決定性的作用[2].試驗(yàn)證明，粒徑越大、形態(tài)細(xì)長(zhǎng)和扁平的粗骨料所占比例越大，則水分越容易在粗骨料表面聚集，從而造成 ITZ 的強(qiáng)度減弱，即粗骨料的形狀和強(qiáng)度會(huì)間接影響混凝土試件的力學(xué)特性[3].現(xiàn)有混凝土力學(xué)特性研究中，為了簡(jiǎn)化模型，將粗骨料視為剛性材料進(jìn)行數(shù)值模擬[4]，然而對(duì)混凝土試件進(jìn)行單軸壓縮試驗(yàn)發(fā)現(xiàn)，存在粗骨料破碎和主裂縫穿過(guò)粗骨料延伸擴(kuò)展的現(xiàn)象.因此，研究粗骨料破碎對(duì)混凝土力學(xué)性能的影響具有重要意義.

Wriggers等[5]使用有限元法對(duì)混凝土進(jìn)行建模,基于蒙特卡洛方法隨機(jī)生成粒徑、尺寸、位置不同的球形粗骨料，通過(guò)塑性損傷本構(gòu)關(guān)系預(yù)測(cè)裂縫的擴(kuò)展，研究表明裂縫首先發(fā)生在骨料周圍，然后在砂漿內(nèi)傳播.周正峰等[6]采用黏聚力單元模型模擬混凝土中裂縫的萌生和擴(kuò)展，結(jié)果表明，混凝土強(qiáng)度隨ITZ強(qiáng)度的增加而增加.近年來(lái)，鑒于離散元方法在模擬裂縫的便捷性，學(xué)者們使用離散元方法(DEM)對(duì)混凝土進(jìn)行中尺度建模[7-8].Nitka等[9]采用不同力學(xué)性能的球形黏結(jié)顆粒模擬砂漿和骨料，而裂縫則由黏結(jié)鍵破壞表示.目前，關(guān)于顆粒破碎的研究大多集中在堆石料、鈣質(zhì)砂、鐵路道砟等方面[10]，混凝土研究則涉及較少，粗骨料破碎對(duì)混凝土力學(xué)特性的影響有待進(jìn)一步研究.關(guān)于顆粒破碎的模擬方法主要有連續(xù)-離散耦合法(Combined FDEM)和離散元法(DEM)[11]，且基于離散元的顆粒破碎法具有計(jì)算效率高、能考慮顆粒的復(fù)雜形狀等優(yōu)點(diǎn)[12].

本文基于離散元軟件 EDEM，建立了粗骨料-砂漿-ITZ的三相混凝土離散元模型.采用黏結(jié)顆粒模型(BPM)建立粗骨料破碎模型，以模擬粗骨料和砂漿之間的力學(xué)特性，從細(xì)觀層次探究粗骨料破碎對(duì)混凝土變形和破壞特性的影響.

1 數(shù)值模型

本文將混凝土視為砂漿、粗骨料和ITZ組成的三相復(fù)合材料，將ITZ建模為砂漿和骨料顆粒之間的弱界面，其實(shí)現(xiàn)方式為均勻削弱骨料黏結(jié)力或刪除部分黏結(jié)鍵[13].考慮骨料與骨料、骨料與砂漿、砂漿與砂漿之間的黏結(jié)特性，采用三維掃描技術(shù)對(duì)粗骨料進(jìn)行掃描，并重建其三維幾何數(shù)字模型.基于黏結(jié)顆粒模型(BPM)，模擬粗骨料與砂漿之間的力學(xué)行為，建立混凝土三相離散元模型.

1.1 BPM模型

BPM模型能實(shí)現(xiàn)顆粒之間的黏結(jié)與顆粒破碎模擬.在BPM模型中，黏結(jié)鍵將小顆粒單元黏結(jié)成為大塊的物料模型[14].在外力作用下小顆粒單元間的黏結(jié)鍵破碎，從而模擬物料破碎、產(chǎn)生裂縫的效果，適用于模擬巖石破碎、混凝土產(chǎn)生裂縫等脆性材料破壞的過(guò)程.

BPM模型黏結(jié)計(jì)算示意圖見(jiàn)圖1.圖中，R1和R′1分別為顆粒1的物理半徑和黏結(jié)半徑；R2和R′2分別為顆粒2的物理半徑和黏結(jié)半徑.由于黏結(jié)半徑等于或略大于顆粒半徑，因此生成黏結(jié)鍵的兩顆?？赡軟](méi)有發(fā)生真正接觸.在生成黏結(jié)鍵的瞬間，顆粒的合力和合力矩會(huì)重置為0，計(jì)算每一時(shí)間步長(zhǎng)下合力和合力矩的變化值，計(jì)算公式為

(1)

式中，ΔFn、ΔFt分別為法向和切向力變化值；ΔMn、ΔMt分別為法向和切向力矩變化值；vn、vt分別為法向和切向速度；ωn、ωt分別為法向和切向角速度；Kn、Kt分別為法向和切向剛度；Δt為計(jì)算時(shí)間步長(zhǎng)；S為黏結(jié)區(qū)域面積；J為黏結(jié)區(qū)域橫截面的極慣性矩.對(duì)于球形顆粒單元，有

(2)

(3)

圖1 BPM模型黏結(jié)計(jì)算示意圖

BPM模型中的黏結(jié)區(qū)域?yàn)閳A盤(pán)形狀，其外圍的最大拉應(yīng)力σmax和最大剪應(yīng)力τmax分別為

(4)

(5)

式中，F(xiàn)n、Ft分別為法向力和切向力；Mn和Mt分別為法向力矩和切向力矩.

當(dāng)σmax大于抗拉強(qiáng)度或τmax大于抗剪強(qiáng)度時(shí)，顆粒間的黏結(jié)鍵斷裂.黏結(jié)鍵斷裂后，對(duì)應(yīng)的力、力矩、剛度等不再作用于顆粒單元.

1.2 黏結(jié)顆粒BPM模型的參數(shù)標(biāo)定

ITZ是指骨料和砂漿之間寬度為15～30 μm的低強(qiáng)度、高孔隙率的薄層.為簡(jiǎn)化處理，參考文獻(xiàn)[9,13]的方法，將ITZ建模為骨料和水泥砂漿顆粒之間的薄弱接觸，此時(shí)ITZ不存在物理寬度.骨料與骨料、骨料與砂漿以及砂漿與砂漿之間的接觸均采用BPM模型模擬.

選用強(qiáng)度為C30的混凝土，尺寸為 150 mm×150 mm×150 mm.利用PLS-500T型微機(jī)控制電液伺服加載試驗(yàn)機(jī)對(duì)混凝土試件進(jìn)行單軸壓縮試驗(yàn)，壓板的加載速率為1 mm/min.將加載板位移除以試件初始長(zhǎng)度可得軸向應(yīng)變，最終得到混凝土單軸壓縮的應(yīng)力-應(yīng)變曲線.在EDEM軟件中建立混凝土的BPM模型(見(jiàn)圖2).砂漿顆粒的半徑為3.5 mm，隨機(jī)生成的粗骨料顆粒粒徑范圍為10～31.5 mm.將粗骨料三維模型導(dǎo)入EDEM軟件后，軟件將粗骨料顆粒默認(rèn)為剛體，在外荷載的作用下無(wú)法變形和破壞，因此將粗骨料作為不可破碎的剛體.

圖2 混凝土單軸壓縮離散元模型

BPM模型中包含如下3類參數(shù)：① 材料的本征參數(shù)，如密度、剪切模量、泊松比等；② 接觸參數(shù)，如碰撞恢復(fù)系數(shù)、靜摩擦系數(shù)、滾動(dòng)摩擦系數(shù)等；③ 黏結(jié)參數(shù)，如單位面積的法向剛度、單位面積的切向剛度、臨界法向力、臨界切向力、黏結(jié)半徑等.C30混凝土中的材料為砂漿和粗骨料；模擬混凝土的單軸壓縮試驗(yàn)時(shí)，壓板和底板為鋼材.因此，需要確定砂漿、粗骨料、鋼材的本征參數(shù)及三者之間的基本接觸參數(shù).此外，還需對(duì)砂漿顆粒與砂漿顆粒、砂漿顆粒與粗骨料顆粒間的黏結(jié)參數(shù)進(jìn)行參數(shù)標(biāo)定.

本征參數(shù)是材料自身的特性參數(shù)，參考文獻(xiàn)[7，15]，砂漿的泊松比為0.2，密度為2 000 kg/m3，剪切模量為0.56 GPa，彈性模量為1.3 GPa；粗骨料的泊松比為0.3，密度為2 630 kg/m3，剪切模量為0.72 GPa，彈性模量為5.3 GPa；鋼材的泊松比為0.3，密度為7 800 kg/m3，剪切模量為0.92 GPa，彈性模量為6.5 GPa.

基本接觸參數(shù)可由室內(nèi)試驗(yàn)或者現(xiàn)場(chǎng)試驗(yàn)測(cè)得.文獻(xiàn)[16]利用室內(nèi)試驗(yàn)對(duì)混凝土的基本接觸參數(shù)進(jìn)行了研究，本文選取的接觸參數(shù)見(jiàn)表1.

表1 接觸參數(shù)

黏結(jié)參數(shù)與顆粒粒徑、粗骨料形狀、混凝土配合比等有關(guān)，且對(duì)混凝土的宏觀力學(xué)性能存在較大影響，需要對(duì)其進(jìn)行標(biāo)定.參數(shù)標(biāo)定時(shí)，對(duì)C30混凝土試件進(jìn)行單軸壓縮試驗(yàn)，并在EDEM軟件中模擬，不斷更改黏結(jié)參數(shù)，直至數(shù)值模擬得到的應(yīng)力-應(yīng)變曲線與試驗(yàn)結(jié)果相符.

標(biāo)定過(guò)程中初始單位面積的法向剛度kn和單位面積的切向剛度kt分別為

(6)

(7)

式中，E、ν分別為材料的彈性模量和泊松比.

在kn和kt的基礎(chǔ)上，反復(fù)修改黏結(jié)參數(shù)，直到符合混凝土宏觀力學(xué)特性，結(jié)果見(jiàn)表2.

表2 混凝土的黏結(jié)參數(shù)

圖3給出了室內(nèi)試驗(yàn)和數(shù)值模擬得到的應(yīng)力-應(yīng)變曲線對(duì)比圖.由圖可知，試驗(yàn)結(jié)果和模擬結(jié)果的曲線變化趨勢(shì)相似.在加載初期，室內(nèi)試驗(yàn)曲線和數(shù)值模擬曲線呈現(xiàn)非線性變化.隨著加載的進(jìn)行，混凝土進(jìn)入線彈性階段，曲線均近似線性變化，直至達(dá)到比例極限.隨后，混凝土開(kāi)始軟化，應(yīng)力-應(yīng)變曲線開(kāi)始彎曲，且隨著加載的進(jìn)行，曲線的彎曲程度不斷增加.當(dāng)模擬值和試驗(yàn)值分別達(dá)到峰值應(yīng)力fc和f′c時(shí)，曲線開(kāi)始迅速下降.當(dāng)應(yīng)力值達(dá)到收斂值時(shí)，混凝土完全破壞.

表3給出了室內(nèi)試驗(yàn)和數(shù)值模擬的混凝土強(qiáng)度指標(biāo)對(duì)比.由表可知，抗壓強(qiáng)度、抗拉強(qiáng)度、抗折強(qiáng)度的誤差都在10%以內(nèi)，說(shuō)明數(shù)值模擬的準(zhǔn)確性較高.圖4為混凝土最終破壞形態(tài)的對(duì)比圖.由圖可知，在外荷載的作用下，試件底部和頂部生成細(xì)小裂縫；隨著荷載的不斷增大，裂縫延伸擴(kuò)展，最終造成混凝土試件的破壞.數(shù)值模擬得到的最終破壞形態(tài)與實(shí)際單軸壓縮試驗(yàn)結(jié)果相似性較好.

圖3 混凝土單軸壓縮試驗(yàn)的應(yīng)力-應(yīng)變曲線對(duì)比圖

表3 室內(nèi)試驗(yàn)和數(shù)值模擬的強(qiáng)度指標(biāo)對(duì)比 MPa

(a) 室內(nèi)試驗(yàn)

(b) 數(shù)值模擬

1.3 顆粒破碎BPM模型的參數(shù)標(biāo)定

離散元模擬顆粒破碎的方法包括基于BPM模型的顆粒破碎方法和碎片置換法(FRM)[11]. 本文采用基于BPM模型的顆粒破碎方法模擬粗骨料破碎，其原理是通過(guò)黏結(jié)鍵將小顆粒黏結(jié)成的團(tuán)聚體替代原始顆粒(即母顆粒).當(dāng)黏結(jié)鍵承受的剪切力或者法向力超過(guò)極限值時(shí)，黏結(jié)鍵斷裂，顆粒也會(huì)隨之破碎.

選擇4種代表性形狀的碎石作為粗骨料，編號(hào)分別為Ⅰ～Ⅳ，其產(chǎn)地為天津，母巖為花崗巖，粒徑為10.0～31.5 mm.為反映混凝土粗骨料的三維結(jié)構(gòu)特征，采用三維掃描技術(shù)對(duì)粗骨料進(jìn)行掃描并重建其三維幾何數(shù)字模型.為使仿真模擬效果更接近于實(shí)際情況，粗骨料用非球形顆粒表示，在EDEM軟件中構(gòu)建相應(yīng)的BPM模型，結(jié)果見(jiàn)圖5.

(a) 顆粒Ⅰ

(b) 顆粒Ⅱ

(c) 顆粒Ⅲ

(d) 顆粒Ⅳ

計(jì)算并對(duì)比真實(shí)顆粒和離散元模型顆粒形狀參數(shù)，包括球度s、凸度Cx、主尺度(長(zhǎng)軸長(zhǎng)度a、中軸長(zhǎng)度b和短軸長(zhǎng)度c)等.球度S用于衡量顆粒與球體的接近程度，定義為等同球體的表面積與顆粒表面積的比值；凸度反映了三維顆粒的緊湊程度.球度和凸度的表達(dá)式分別為

(8)

(9)

式中，V為顆粒體積；As為顆粒表面積；VCH為凸包的體積.

表4給出了顆粒Ⅰ的形狀參數(shù)對(duì)比.由表可知，三維模型對(duì)粗骨料的還原性較高，離散元模型能反映出粗骨料的形狀特性.

表4 顆粒Ⅰ的形狀參數(shù)對(duì)比

利用WDW-300型電子式萬(wàn)能材料試驗(yàn)機(jī)對(duì)顆粒Ⅰ進(jìn)行單軸壓縮試驗(yàn)，加載速度設(shè)置為2 mm/min.粗骨料參數(shù)標(biāo)定結(jié)果見(jiàn)表5.

基于標(biāo)定參數(shù)得到的力-位移曲線見(jiàn)圖6.由圖可知，力-位移曲線中包含多個(gè)峰值點(diǎn).隨著加載的進(jìn)行，在達(dá)到每個(gè)峰值點(diǎn)后力均會(huì)急劇下降，表示顆粒發(fā)生一次破碎.對(duì)比數(shù)值模擬和室內(nèi)試驗(yàn)中第1次發(fā)生破碎時(shí)的力和位移，其誤差分別為3.1%和7.6%，故可認(rèn)為表5中的黏結(jié)參數(shù)能夠反映出粗骨料真實(shí)的力學(xué)性質(zhì).

表5 粗骨料的黏結(jié)參數(shù)

圖6 顆粒Ⅰ在單軸壓縮試驗(yàn)下的力-位移曲線

圖7對(duì)比了顆粒Ⅰ在室內(nèi)試驗(yàn)和數(shù)值模擬中的最終破壞模式.由圖可知，無(wú)論是室內(nèi)試驗(yàn)還是數(shù)值模擬中，顆粒Ⅰ均發(fā)生了剪切破壞.裂縫從粗骨料中部萌生并向四周擴(kuò)展，直至貫穿于整個(gè)粗骨料，從而造成粗骨料的破碎，說(shuō)明數(shù)值模擬得到的破壞模式與試驗(yàn)結(jié)果吻合較好.

(a) 室內(nèi)試驗(yàn)

(b) 數(shù)值模擬

2 骨料破碎對(duì)混凝土力學(xué)特性的影響

為使數(shù)值模型更接近于實(shí)際試驗(yàn)，將部分粗骨料設(shè)置為可破碎狀態(tài)，探究可破碎粗骨料的質(zhì)量分?jǐn)?shù)w對(duì)混凝土試件的應(yīng)力-應(yīng)變曲線、抗壓強(qiáng)度及最終破壞模式的影響.受計(jì)算效率限制，且實(shí)際混凝土受壓過(guò)程中只有部分粗骨料發(fā)生破壞，故將w設(shè)置為0%～40%.為減少可破碎骨料隨機(jī)分布對(duì)模擬結(jié)果的影響，每個(gè)工況下進(jìn)行3組試驗(yàn)，將結(jié)果取平均值.w=10%時(shí)的混凝土模型見(jiàn)圖8.

圖8 混凝土模型(w=10%)

2.1 應(yīng)力-應(yīng)變曲線

圖9為含不同可破碎粗骨料質(zhì)量分?jǐn)?shù)的混凝土的應(yīng)力-應(yīng)變曲線.由圖可知，當(dāng)混凝土試件中的w=0%時(shí)，混凝土試件的峰值應(yīng)力fc= 33.17 MPa，對(duì)應(yīng)的峰值應(yīng)變?chǔ)與=1.98×10-3.不同可破碎粗骨料質(zhì)量分?jǐn)?shù)下的混凝土應(yīng)力-應(yīng)變曲線形狀相似，變化趨勢(shì)基本相同.粗骨料的破碎會(huì)影響混凝土的峰值應(yīng)力和應(yīng)變.當(dāng)0%5%時(shí)峰值應(yīng)力和應(yīng)變隨w的增加逐漸減小.

圖9 不同可破碎粗骨料質(zhì)量分?jǐn)?shù)混凝土應(yīng)力-應(yīng)變曲線

綜上可知，粗骨料破碎對(duì)混凝土的變形特性存在顯著影響.粗骨料的破碎會(huì)增強(qiáng)混凝土的硬化特性和延性，改變峰值應(yīng)力和應(yīng)變，且峰值應(yīng)力和峰值應(yīng)變整體上隨可破碎粗骨料質(zhì)量分?jǐn)?shù)的增加而減少.

2.2 抗壓強(qiáng)度

混凝土的抗壓強(qiáng)度受粗骨料粒徑、水灰比、水泥級(jí)配等多種因素影響，是混凝土結(jié)構(gòu)最重要的力學(xué)特性之一[17].圖10為可破碎粗骨料質(zhì)量分?jǐn)?shù)對(duì)抗壓強(qiáng)度的影響曲線.由圖可知，隨著混凝土中可破碎粗骨料質(zhì)量分?jǐn)?shù)的增大，抗壓強(qiáng)度先增大后減小，最后趨于穩(wěn)定.

圖10 粗骨料破碎對(duì)抗壓強(qiáng)度的影響

圖10中的曲線可分為3個(gè)階段：① 當(dāng)w<3%時(shí)，混凝土的抗壓強(qiáng)度隨w的增大而略微增加.w=0%時(shí)混凝土的抗壓強(qiáng)度為33.17 MPa；w=3%時(shí)，混凝土的抗壓強(qiáng)度達(dá)到最大值35.23 MPa.產(chǎn)生這一現(xiàn)象的原因在于，采用離散元方法模擬混凝土?xí)r，離散元顆粒間會(huì)存在一定的孔隙，w較小時(shí)粗骨料破碎產(chǎn)生的顆粒會(huì)對(duì)孔隙進(jìn)行填充，從而增加了混凝土試件的密實(shí)性，使混凝土的抗壓強(qiáng)度增大.②混凝土抗壓強(qiáng)度隨w的增大而減小，w=5%時(shí)的抗壓強(qiáng)度為33.24 MPa，與w=0%時(shí)基本一致.隨著w的增加，混凝土抗壓強(qiáng)度迅速下降，可破碎粗骨料對(duì)裂縫的阻礙作用減弱，導(dǎo)致混凝土更易破壞.當(dāng)5%

2.4 破壞模式

單軸壓縮試驗(yàn)過(guò)程中，混凝土內(nèi)部局部區(qū)域出現(xiàn)應(yīng)力集中并萌生裂縫，在外荷載的作用下裂縫會(huì)繼續(xù)擴(kuò)展，其擴(kuò)展方向與外荷載方向基本平行.當(dāng)裂縫擴(kuò)展到一定程度后，裂縫之間相互貫穿形成主裂縫.主裂縫決定了混凝土試件的破壞形式.圖11給出了考慮和不考慮粗骨料破碎時(shí)混凝土試件的破壞模式.由圖可知，不考慮粗骨料破碎時(shí)，裂縫沿著粗骨料顆粒與砂漿顆粒的交界面延伸擴(kuò)展，粗骨料沒(méi)有破碎.對(duì)于可破碎粗骨料質(zhì)量分?jǐn)?shù)為15%的混凝土試件，部分粗骨料發(fā)生破碎且存在主裂縫穿過(guò)粗骨料的現(xiàn)象，位于主裂縫位置上的可破碎粗骨料破碎最為嚴(yán)重，其他位置損傷程度較小.

圖12給出了混凝土破壞前后的黏結(jié)鍵.試件破壞后，斷裂黏結(jié)鍵個(gè)數(shù)為24 108.斷裂黏結(jié)鍵顯示了混凝土中裂縫的延伸趨勢(shì)和易破碎位置，主裂縫貫穿整個(gè)試件，這與圖12中的破壞模式吻合較好.

(a) w=0%

(b) w=15%

(a) 破壞前

(b) 破壞后

3 結(jié)論

1) 采用三維掃描技術(shù)構(gòu)建了粗骨料三維幾何模型.通過(guò)BPM建立了粗骨料破碎模型，以模擬粗骨料和砂漿之間的力學(xué)特性，進(jìn)而實(shí)現(xiàn)了考慮粗骨料破碎的混凝土粗骨料-砂漿-ITZ的三相離散元模型.

2) 含不同可破碎粗骨料質(zhì)量分?jǐn)?shù)的混凝土的應(yīng)力-應(yīng)變曲線形狀相似，發(fā)展趨勢(shì)相同.峰值應(yīng)力和應(yīng)變整體上隨可破碎粗骨料質(zhì)量分?jǐn)?shù)的增加而減少.

3) 粗骨料的破碎會(huì)影響混凝土的抗壓強(qiáng)度.當(dāng)w<3%時(shí)，少量粗骨料破碎產(chǎn)生的顆粒會(huì)填充于離散元模型中顆粒間的孔隙中，增加密實(shí)性，混凝土抗壓強(qiáng)度出現(xiàn)短暫的上升；隨著w的增加，混凝土更易產(chǎn)生裂縫，抗壓強(qiáng)度迅速下降.當(dāng)w達(dá)到臨界值32%后，抗壓強(qiáng)度幾乎不變.

4) 從細(xì)觀角度解釋了裂縫不僅沿著ITZ拓展，還存在著貫穿粗骨料的現(xiàn)象.