秦擁軍，田盼盼，于 江，羅 玲

(新疆大學 建筑工程學院，新疆 烏魯木齊 830047)

再生混凝土細觀層次的數(shù)字圖像模型分析

秦擁軍，田盼盼，于 江，羅 玲

(新疆大學 建筑工程學院，新疆 烏魯木齊 830047)

針對再生混凝土已有細觀骨料模型算法復雜、建模理想化和處理技術(shù)手段傳統(tǒng)等問題，基于真實骨料圖像，結(jié)合數(shù)字圖像處理和邊界矢量化技術(shù)，通過ANSYS軟件的參數(shù)化語言編程，建立了再生混凝土數(shù)字圖像細觀模型。應用數(shù)字圖像細觀模型模擬了單軸受壓立方體試塊的應變云圖，分析了宏觀力學響應和再生粗骨料對再生混凝土抗壓強度的影響。研究結(jié)果表明：大多數(shù)再生混凝土裂紋從大骨料邊界產(chǎn)生、發(fā)展和形成。再生粗骨料摻量(質(zhì)量分數(shù))在50%內(nèi)的再生混凝土，不同骨料模型模擬的抗壓強度值不同。本文提出的數(shù)字圖像細觀模型模擬結(jié)果均符合試驗現(xiàn)象和力學規(guī)律，證明了建模方法的合理性和有效性。

細觀層次；再生混凝土；數(shù)字圖像模型；粗骨料；力學性能

0 引言

再生混凝土是發(fā)展生態(tài)綠色混凝土、實現(xiàn)建筑資源可持續(xù)發(fā)展的重要途徑之一[1-4]。目前的研究從宏觀層次上很難客觀合理地解釋再生混凝土的裂紋擴展規(guī)律、破壞過程及靜、動態(tài)破壞機理，且難以描述混凝土材料的軟化行為。因此，研究在不同應力路徑下再生混凝土破壞過程中的損傷微裂紋孕育、發(fā)展和貫通過程[5]，明確細觀層次機理是一個不可缺少的關(guān)鍵環(huán)節(jié)[6]。

目前，國內(nèi)外學者對混凝土細觀層次的力學模型和數(shù)值模擬進行了研究，取得了一定成果[7-10]，但仍缺乏成熟的細觀模型[7]，原有的格構(gòu)模型[8]、M-H模型[9]和隨機力學特性模型[10]等均存在一些不足，如幾何建模過程復雜、建模理想化和處理技術(shù)手段傳統(tǒng)等。為確保再生混凝土結(jié)構(gòu)的安全，需要研究細觀力學模型，為再生混凝土裂縫擴展方向及其變化規(guī)律的研究提供重要的理論基礎(chǔ)[11]。

數(shù)字圖像處理是指將圖像信號利用計算機進行處理轉(zhuǎn)換成數(shù)字信號的過程[12]，具有信息量大、傳輸速度快、作用距離遠等優(yōu)點，飛速發(fā)展的計算機技術(shù)使得圖像處理技術(shù)成為現(xiàn)代社會快速獲取信息的重要手段之一。圖像分析技術(shù)和理論也逐漸應用于再生混凝土材料的研究中。

本文針對再生混凝土細觀骨料模型算法復雜和處理技術(shù)手段傳統(tǒng)等問題，基于真實骨料圖像和邊界矢量化技術(shù)提出了再生混凝土數(shù)字圖像細觀模型，通過試驗驗證了該模型，進一步研究了再生粗骨料摻量和形狀對再生混凝土的影響。

1 再生混凝土數(shù)字圖像模型的建立

數(shù)字圖像模型是指基于真實骨料圖像，利用數(shù)字圖像處理技術(shù)和邊界矢量化技術(shù)建立的再生混凝土細觀模型。再生混凝土數(shù)字圖像模型生成流程如圖1所示。

1.1 幾何模型的生成

根據(jù)試驗方案，設(shè)計了再生骨料摻量(質(zhì)量分數(shù)，下同)分別為0%、30%和50%的R0、R3和R53組試件，每組3個。由于篇幅有限，本文主要選取R3組第2個試件(R3-2)進行分析，其尺寸為150 mm×150 mm,再生骨料摻量為30%。

應用MATLAB軟件做圖像處理，應用AutoCAD軟件在關(guān)鍵點拾取骨料邊界，形成不規(guī)則多邊形的骨料。邊界矢量化借助Auto CAD軟件中“面域”功能，將每個骨料矢量處理成整體。

圖1 再生混凝土數(shù)字圖像模型生成流程圖

由于Walaraven公式[13]建立了骨料級配和骨料面積的關(guān)系，使得細觀骨料模型能夠在二維平面內(nèi)進行分析[14]，本文模型創(chuàng)建基于二維平面模型開展。

將試件切割處理后，得到如圖2所示的R3-2試件切割圖片，導入ANSYS軟件中，得到如圖3所示的R3-2試件骨料幾何模型。生成幾何模型是進行數(shù)字圖像模型模擬的基本條件。

圖2 R3-2試件切割圖片圖3 R3-2試件骨料幾何模型圖

1.2 材料屬性及網(wǎng)格劃分

將模型導入ANSYS軟件，用參數(shù)化設(shè)計語言(APDL)建立試塊邊界面，通過布爾減運算從邊界面“減去”已有骨料，得到混凝土砂漿面。模型有3種材料：砂漿、原生骨料和再生骨料。自由網(wǎng)格劃分采用網(wǎng)格尺寸為0.002 m 的PLANE82三角形單元。

由于再生混凝土建筑物拆除過程中將不同強度等級、不同種類和不同來源的混凝土混在一起，很難保證再生混凝土成分的單一性，導致再生混凝土性能測試結(jié)果有差異[15]，因此很難統(tǒng)一模型中的材料屬性。本文中，數(shù)字圖像模型材料屬性參數(shù)一方面來自于本課題組試驗數(shù)據(jù)[16-18]；另一方面借鑒大量其他學者的研究成果[14,19-20]。3種材料的基本力學參數(shù)如表1所示。

表1 3種材料的基本力學參數(shù)

2 再生混凝土單軸受壓數(shù)值模擬

2.1 應變云圖

建立骨料幾何模型后，對不同材料賦相應的材料屬性，模擬R3-2試件再生混凝土單軸受壓過程，得到 R3-2試件應變云圖，如圖4所示。圖4中，MIN和MAX分別表示計算結(jié)果的應變最小值和最大值，負值代表壓應變，正值代表拉應變。

圖4 R3-2試件應變云圖

從圖4中可以看出：第1荷載步(見圖4a)，應力較小，再生混凝土變形主要來自骨料和水泥砂漿的較小彈性變形，應變值較小。第17荷載步(見圖4b)，應變逐漸增大，應變區(qū)域逐漸增加、擴展、相連，主要出現(xiàn)在骨料與骨料距離較近的砂漿間，但也有少部分應變區(qū)域是朝著骨料與骨料距離較大的砂漿間發(fā)展。第34荷載步(見圖4c)，砂漿的應變逐漸增大，增大程度比第17荷載步明顯，骨料應變比砂漿應變小。第50荷載步(見圖4d)，砂漿的應變繼續(xù)增加，骨料應變也在增加，骨料彈性模量遠大于砂漿，根據(jù)胡克定律，在彈性范圍內(nèi)，應變與彈性模量成反比，應變貫穿再生混凝土內(nèi)部，最終使得再生混凝土破裂直至完全破壞。

2.2 宏觀響應及裂縫形態(tài)

根據(jù)滑動裂紋模型理論，在受壓狀態(tài)下微裂紋的進一步擴展需要壓應力的持續(xù)增長來維持[21]。再生混凝土宏觀裂紋的開展，來自試件應力持續(xù)增長的內(nèi)部損傷發(fā)展，是壓應力持續(xù)增長的宏觀響應，圖5為3個試件破壞圖和裂紋圖。圖5中：R0-2為R0組第2個試件，R3-2為R3組第2個試件，R5-3為R5組第3個試件。由圖5可看出兩個規(guī)律：

(1)大多數(shù)裂紋從大骨料邊界產(chǎn)生發(fā)展，最終大多數(shù)裂紋也是沿著大骨料邊界形成。因為大骨料與砂漿接觸的面積小于同等體積小骨料與砂漿接觸的面積，大骨料相對小骨料握裹力小，所以裂縫首先從大骨料邊界產(chǎn)生。

(2)骨料-骨料間距較小的地方產(chǎn)生發(fā)展的裂紋較多。因為骨料-骨料間距較小的地方，受力面積小，應力增長快，加之兩骨料擠壓效果明顯，所以應變增大較快，較多裂紋從骨料-骨料間距較小的地方產(chǎn)生發(fā)展。

3 再生粗骨料對再生混凝土抗壓強度的影響

3.1 再生粗骨料摻量對再生混凝土抗壓強度的影響

再生粗骨料摻量分別為0%、30%和50%時，對R0、R3 和R5這3組試件分別進行抗壓強度的研究。再生粗骨料摻量對再生混凝土抗壓強度的影響如圖6所示。由圖6可以看出：再生粗骨料摻量越高，再生混凝土抗壓強度越低。產(chǎn)生上述現(xiàn)象的原因是：再生混凝土是將廢棄混凝土經(jīng)過破碎、篩分、清洗，并按一定級配和比例配制的混凝土[22]，再生粗骨料外包裹舊砂漿，形成了“石子-舊砂漿”的“薄弱環(huán)節(jié)”，使得粗骨料性能略弱于原生骨料性能，再生粗骨料替代原生骨料摻入混凝土中，反映到宏觀性能上，表現(xiàn)為混凝土抗壓性能降低。這符合試驗現(xiàn)象和力學規(guī)律。

圖5 3個試件破壞圖和裂紋圖

3.2 再生粗骨料形狀對再生混凝土抗壓強度的影響

再生粗骨料摻量分別為0%、30%和50%時，對圓形、多邊形和數(shù)字圖像這3種細觀骨料模型進行研究，可得在相同材料參數(shù)及軸心受壓條件下，骨料形狀不同，混凝土抗壓強度值也有所不同。再生混凝土試件抗壓強度從大到小依次是：圓形骨料模型>數(shù)字圖像骨料模型>多邊形骨料模型，如圖7所示。數(shù)字圖像骨料模型是基于真實骨料切片，骨料形狀介于理想的圓形和多邊形之間。圓形骨料的邊緣光滑連續(xù)，應力分布均勻，而多邊形骨料尖角及突出部分容易產(chǎn)生應力集中現(xiàn)象[23]。當有力作用于試件時，多邊形骨料尖角及突出部分首先發(fā)生破壞，其符合力學試驗現(xiàn)象和規(guī)律。

圖6 再生粗骨料摻量對再生混凝土抗壓強度的影響圖7 不同骨料模型對再生混凝土抗壓強度的影響

4 結(jié)論

(1)本文建立的數(shù)字圖像細觀模型以真實骨料切片為基礎(chǔ)，通過有限元軟件和試驗進行驗證和分析，驗證了此建模方法在再生混凝土單軸壓縮試驗范圍內(nèi)是可借鑒和使用的。

(2)大多數(shù)裂紋從大骨料邊界產(chǎn)生發(fā)展，最終大多數(shù)裂紋也是沿著大骨料邊界形成。

(3)再生粗骨料摻量50%內(nèi)的再生混凝土，不同骨料形狀對應的細觀骨料模型模擬的再生混凝土抗壓強度也不同，抗壓強度從大到小依次是：圓形骨料模型>數(shù)字圖像骨料模型>多邊形骨料模型。

[1] 楊醫(yī)博,鄭子麟,郭文瑛,等.全再生細骨料的制備及其對混凝土性能影響的試驗研究[J].功能材料,2016,47(4):157-163.

[2] 常彥錚.再生混凝土抗壓試驗的數(shù)值模擬與等效抗壓強度折減系數(shù)研究[D].西安:西安建筑科技大學，2011.

[3] 汪振雙,譚曉倩.鋼纖維再生粗集料混凝土的力學性能和抗凍性研究[J].硅酸鹽通報,2016,34(4):1184-1187.

[4] 王磊,李泉龍,陳杏.再生粗骨料混凝土的變形性能及影響因素分析[J].建筑結(jié)構(gòu),2016,46(2):30-34.

[5] 王欽亭,王云飛.砂巖單軸壓縮與劈裂破壞過程的損傷研究[J].河南科技大學學報(自然科學版),2014,35(4):54-58.

[6] 金瀏.細觀混凝土分析模型與方法研究[D].北京:北京工業(yè)大學,2014.

[7] 張曉龍.荷載損傷混凝土損傷指標選取及其對耐久性的影響[D].杭州:浙江大學,2015.

[8] VAN M J.Fracture processes of concrete-assessment of material parameters for fracture model[M].Florida:CRC Press,1997.

[9]MOHAMED A R,HANSEN W.Micromechanical modeling of concrete response under static loading-part1:model development and validation[J].ACI materials journal,1999,96(2):196-203.

[10] 王玉琳,王向東,韓金啟,等.基于損傷理論的混凝土Ⅰ-Ⅱ復合型裂縫擴展角[J].河南科技大學學報(自然科學版),2015,36(6):51-60.

[11] 唐春安,朱萬成.混凝土損傷與斷裂：數(shù)值試驗[M].北京:科學出版社,2003.

[12] 趙宇杰,姚娟文,趙彥敏.基于MATLAB的數(shù)字圖像色彩空間轉(zhuǎn)換與應用[J].電子技術(shù)與軟件工程,2015(5):113.

[13] WALARAVEN J C,REINHARDT H W.Theory and experiments on the mechanical behavior of cracks in plain and reinforced concrete subjected to shear loading[J].HERON,1991,26(1A):26-35.

[14] 嵇紹華.纖維混凝土細觀結(jié)構(gòu)及破裂過程數(shù)值模擬研究[D].西安:長安大學,2014.

[15] 常彥錚.應用ANSYS建立再生混凝土隨機骨料模型的方法研究[J].價值工程,2010(35):158-159.

[16] 于江,梁芮,秦擁軍.C30 再生粗骨料混凝土的性能試驗研究[J].建筑科學,2015,31(1):59-63.

[17] 李蕾.新疆寒冷地區(qū)礦物摻合料再生混凝土配合比設(shè)計及其性能研究[D].烏魯木齊:新疆大學,2014.

[18] 梁芮.烏魯木齊市廢棄混凝土再生粗骨料性能及應用研究[D].烏魯木齊:新疆大學,2014.

[19] 叢卓紅,閆紅光,鄭南翔.水泥穩(wěn)定碎石用硬質(zhì)砂巖關(guān)鍵性指標[J].長安大學學報(自然科學版),2011,31(5):12-17.

[20] 曹文平.高強高性能混凝土細觀數(shù)值演化分析[D].西安:西安建筑科技大學,2012.

[21] 蒲成志.巖體斷裂與蠕變損傷破壞機理研究[D].長沙:中南大學,2014.

[22] 陳宗平,應武擋,陳宇良,等.短齡期再生混凝土三軸受壓力學性能及其本構(gòu)關(guān)系[J].建筑材料學報,2015,18(6):935-940.

[23] 宋來忠,沈濤,余波.混凝土二維參數(shù)化骨料模型的創(chuàng)建方法[J].工程力學,2013,30(10):5-13.

國家自然科學基金項目(51668061,51668060)

秦擁軍(1970-),男,江蘇太倉人,副教授,碩士,碩士生導師,主要從事再生混凝土材料方面的研究.

2016-09-05

1672-6871(2017)02-0054-05

10.15926/j.cnki.issn1672-6871.2017.02.010

TU528

A