混凝土單軸壓縮分形損傷演化的CT試驗(yàn)研究

劉書波 劉 婺

(1.沈陽(yáng)建都工程建設(shè)監(jiān)理公司，遼寧沈陽(yáng) 110000;2.遼寧工程技術(shù)大學(xué)建筑工程學(xué)院，遼寧阜新 123000)

地震災(zāi)害荷載作用下，建筑物宏觀方面會(huì)產(chǎn)生裂縫，微觀方面是由于在地震作用下直材料出現(xiàn)突發(fā)性斷裂破壞，此類破壞沒(méi)有明顯預(yù)兆，有較大的危險(xiǎn)性。隨著社會(huì)經(jīng)濟(jì)的進(jìn)步，高層建筑、大跨橋梁及港口建筑結(jié)構(gòu)大量涌現(xiàn)，高強(qiáng)高性能混凝土的使用也日漸廣泛起來(lái)，其雖在強(qiáng)度和性能上較以往的普通混凝土有所提高和改進(jìn)，但制造過(guò)程中的初始缺陷卻同樣難以避免。

眾多實(shí)驗(yàn)顯示，材料的最終宏觀斷裂破壞與其內(nèi)部微裂隙(也即初始損傷)的發(fā)育和聚集有著密切的關(guān)系。近年來(lái)的研究表明，采用先進(jìn)的實(shí)驗(yàn)技術(shù)，引入損傷、非均勻、不連續(xù)的概念研究其內(nèi)部結(jié)構(gòu)，并分析高強(qiáng)高性能混凝土受力破損過(guò)程中的結(jié)構(gòu)變化，對(duì)高強(qiáng)高性能混凝土的配合比設(shè)計(jì)、性能評(píng)價(jià)具有重要意義。CT技術(shù)的出現(xiàn)為研究高強(qiáng)高性能混凝土材料細(xì)觀破壞機(jī)理提供了較以往更先進(jìn)的方法。國(guó)內(nèi)也開(kāi)展了CT技術(shù)在高強(qiáng)高性能混凝土內(nèi)部結(jié)構(gòu)分析和巖石的細(xì)觀結(jié)構(gòu)損傷研究中的應(yīng)用，這些研究為認(rèn)識(shí)高強(qiáng)高性能混凝土材料細(xì)觀損傷的演化提供了有力的支持［1］。借助于室內(nèi)單軸壓縮試驗(yàn)，對(duì)高強(qiáng)高性能混凝土材料細(xì)觀結(jié)構(gòu)損傷演化進(jìn)行了基于CT圖像分析的細(xì)觀力學(xué)研究，重點(diǎn)探討了混凝土材料微觀的損傷演化過(guò)程，對(duì)其力學(xué)性能的損傷演化進(jìn)行了系統(tǒng)的研究。

1 混凝土單軸壓縮的試驗(yàn)設(shè)備

本試驗(yàn)在遼寧工程技術(shù)大學(xué)實(shí)驗(yàn)中心進(jìn)行，其加載方式主要采用軸向壓縮。主要儀器設(shè)備有:TYA-2000型電液式壓力試驗(yàn)機(jī)、西門子64層醫(yī)用螺旋X射線CT斷層診斷儀(Computed Tomography)等。混凝土立方體試塊尺寸為100 mm×100 mm×100 mm，水泥為秦嶺牌525號(hào)硅酸鹽水泥，粗骨料為粒徑10 mm～15 mm的一級(jí)配花崗巖人工碎石，細(xì)骨料為灞河河砂，以寶雞第二熱電廠一級(jí)粉煤灰作為摻合料，標(biāo)準(zhǔn)條件下養(yǎng)護(hù)28 d。

2 混凝土斷裂面的分形特征

1975年，任IBM(International Business Machine)公司研究中心物理部研究員暨哈佛大學(xué)數(shù)學(xué)系教授美籍法國(guó)數(shù)學(xué)家Mandelbrot首先提出了分形(Fractal)的概念［2］，它與混沌理論和耗散理論被稱為20世紀(jì)70年代科學(xué)界三大發(fā)現(xiàn)，是非線性科學(xué)研究領(lǐng)域的重要成果。研究表明，混凝土斷裂要經(jīng)過(guò)斷裂面萌生、斷面擴(kuò)展和動(dòng)態(tài)斷裂的復(fù)雜過(guò)程，該過(guò)程同時(shí)又受應(yīng)力狀態(tài)、材料強(qiáng)度、原始缺陷等眾多因素的影響，從而使混凝土損傷斷裂過(guò)程中出現(xiàn)一系列復(fù)雜的斷面運(yùn)動(dòng)現(xiàn)象，而混凝土本身就是具有多層次自相似的混沌體，混凝土斷裂面也具有分形特征［3］。

3 試驗(yàn)過(guò)程

1)高強(qiáng)高性能混凝土的單軸壓縮試驗(yàn)。該試驗(yàn)采用應(yīng)力控制，加壓過(guò)程中不斷觀察壓力試驗(yàn)機(jī)表盤值，在觀測(cè)試件受壓過(guò)程變化的同時(shí)，記錄下壓力試驗(yàn)機(jī)表盤數(shù)據(jù)，每一次加載對(duì)應(yīng)一個(gè)表盤值，為了能夠盡量使得試件全截面均勻受載以及消除兩受壓端的摩擦約束造成的圍壓作用，控制加載速率為1.2×10－3kN/s，并且在加載時(shí)，試件頂部和底部涂抹潤(rùn)滑劑，直至試件破壞后終止試驗(yàn)。2)CT掃描。選定掃描層厚1.0 mm，層間隔23.8 mm，掃描電壓120 kV，掃描電流300 μA，掃描典型圖像如圖1所示。

圖1 混凝土斷裂面的CT掃描

4 分形特性分析

盡管一個(gè)真實(shí)空間中的分形表面映射到圖像的灰度值表示的方法尚未得到解決，但是利用數(shù)學(xué)函數(shù)方法來(lái)表示圖像的灰度值是可行的。基于分析維數(shù)的數(shù)學(xué)計(jì)算方法可以較為理想的解決裂縫的灰度表面裂縫問(wèn)題。

4.1 分形維數(shù)的數(shù)學(xué)計(jì)算理論

差分盒維數(shù)法(DBC)、盒維數(shù)法(DB)、雙毯(Peleg-blanket)覆蓋法等方法是計(jì)算灰度圖像分形維數(shù)的基本方法，根據(jù)本試驗(yàn)具體的試驗(yàn)數(shù)據(jù)，采用盒維數(shù)法(DB)來(lái)分析混凝土立方體試塊的損傷演化趨勢(shì)。

首先對(duì)分析單元進(jìn)行離散化，保證離散單位能夠基本覆蓋一個(gè)損傷演化區(qū)域，并建立起混凝土斷裂面的裂縫所占的區(qū)域與分得的離散單元的數(shù)學(xué)方程式和網(wǎng)格邊長(zhǎng)r的倒數(shù)，繪制關(guān)系曲線:

其中，r為離散單元的網(wǎng)格邊長(zhǎng);N為裂縫區(qū)域所占的網(wǎng)格數(shù)。

該方程如果呈現(xiàn)線性，則說(shuō)明分形理論可以用來(lái)研究混凝土構(gòu)件表面裂縫的損傷演化發(fā)展。根據(jù)Mandelbrot的盒維數(shù)基本定義，混凝土構(gòu)件表面裂縫的損傷維數(shù)為:

其中，K為常數(shù);D為分形維數(shù)。

4.2 CT圖像分形維數(shù)計(jì)算和分析

圖1中的CT圖像在計(jì)算機(jī)上進(jìn)行不同灰度處理后，采用盒維數(shù)法(DB)，分別取對(duì)數(shù)后擬合得到直線，再根據(jù)該直線的斜率即可求出每幅CT圖像的分形維數(shù)，計(jì)算結(jié)果顯示，高強(qiáng)高性能混凝土斷裂面的分形維數(shù)與普通混凝土的相比，前者要小一些，這是由于:根據(jù)分形理論，分形維數(shù)D與表面形貌的幅值變化劇烈程度有關(guān)，D值大則表面微觀細(xì)節(jié)豐富，D值小則表面相對(duì)平緩，不同于普通混凝土單軸壓縮斷裂時(shí)的沿集料破壞。

通過(guò)加入摻合料來(lái)提高高強(qiáng)混凝土材料的力學(xué)性能，導(dǎo)致試件的脆性性能加大，同時(shí)導(dǎo)致粗細(xì)骨料的破裂減小，對(duì)于試件的破裂面而言，其平整度增大，有效降低了其分形維數(shù)。所以隨著其加載的增大，試件的微觀損傷不斷加劇，進(jìn)而增大了其分形維數(shù)。

通過(guò)分析混凝土試件不同應(yīng)力水平的分形維數(shù)可知:荷載加載初期，分形維數(shù)的變化不大，說(shuō)明此時(shí)混凝土尚未發(fā)生破壞，當(dāng)分形維數(shù)急劇增加的時(shí)候說(shuō)明這一階段高強(qiáng)高性能混凝土材料的損傷演化呈現(xiàn)明顯加速趨勢(shì)，損傷分布更加復(fù)雜化(見(jiàn)圖2)。上述結(jié)論表明，用分形維數(shù)表征高強(qiáng)高性能混凝土斷裂面的破壞程度以及材料內(nèi)部損傷發(fā)展各階段狀態(tài)的可行性［4］。

4.3 損傷演化方程的建立

將高強(qiáng)高性能混凝土內(nèi)部結(jié)構(gòu)圖像的分形維數(shù)變化率定義為一個(gè)如下的損傷變量［5］:

其中，Δw為分形維數(shù)變化的總量;w0，wi分別為初始和第i幅圖像的分形維數(shù)。初始狀態(tài)時(shí)，P=0表示試件無(wú)損傷;加載到一定階段，P=1表示試件宏觀裂紋已經(jīng)形成，試件被破壞。我們可以通過(guò)計(jì)算試驗(yàn)中不同階段CT掃描圖像的分形維數(shù)來(lái)獲取其相應(yīng)階段的損傷變量。用損傷變量和不同加壓階段的應(yīng)力值為坐標(biāo)軸建立二者的關(guān)系，得到如圖3所示的曲線圖。觀察所得損傷變量與不同加壓階段的應(yīng)力值曲線，其呈類似指數(shù)函數(shù)關(guān)系，擬采用指數(shù)函數(shù)進(jìn)行擬合，擬合關(guān)系如下:

其中，a=0.947 7;b=158.730 2;R2=0.910 3;a，b 均為參數(shù);R2為相關(guān)系數(shù)。將參數(shù)代入，所得方程即為損傷演化方程:

圖2 試件應(yīng)力變化與分形維數(shù)的關(guān)系

圖3 試件的損傷演化曲線

5 結(jié)語(yǔ)

1)高強(qiáng)混凝土的斷裂特征與普通混凝土有顯著不同。從斷裂面的形態(tài)看，普通混凝土斷裂時(shí)，其集料一般不斷裂，斷裂面是繞行集料的，從界面區(qū)斷開(kāi)。而高強(qiáng)混凝土則不然，其斷裂面是切斷巖石集料的，使其斷裂面的曲折度減小，真實(shí)斷裂面積變小，斷裂速度加快，脆性比普通混凝土加大。由于分形維數(shù)是用來(lái)表征斷裂面復(fù)雜程度的計(jì)算工具，斷裂面曲折程度越大，分形維數(shù)越大，這一點(diǎn)從本文高強(qiáng)混凝土斷裂面分形維數(shù)的計(jì)算結(jié)果也可以很好地看出。2)在高強(qiáng)高性能混凝土單軸壓縮試驗(yàn)的全過(guò)程中，隨著損傷的不斷積累和演化，由盒維數(shù)法計(jì)算的CT圖像的分形維數(shù)不斷增加，由此可得結(jié)論:分形維數(shù)可作為觀測(cè)研究高強(qiáng)高性能混凝土細(xì)觀損傷演化的定量標(biāo)準(zhǔn)。3)通過(guò)對(duì)高強(qiáng)高性能混凝土損傷演化過(guò)程中不同加壓階段CT圖像分形維數(shù)的分析，定義相應(yīng)的損傷變量，得到了該種材料在單軸受壓時(shí)的損傷演化方程。

［1］張淑娟，賴遠(yuǎn)明，孫志忠，等.CT技術(shù)應(yīng)用于巖石細(xì)觀損傷特性試驗(yàn)研究現(xiàn)狀［J］.甘肅科學(xué)學(xué)報(bào)，2004，16(1):96-100.

［2］張濟(jì)忠.分形［M］.北京:清華大學(xué)出版社，1997.

［3］劉書賢，魏曉剛，王 偉，等.煤礦采動(dòng)與地震耦合作用下建筑物災(zāi)變分析［J］.中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，2013，42(4):526-534.

［4］徐曉鵬，彭瑞東，謝和平，等.基于SEM圖像分維估算的脆性材料細(xì)觀結(jié)構(gòu)演化方法研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2004，23(21):3600-3603.

［5］魏曉剛，劉書賢，麻鳳海，等.動(dòng)力參數(shù)對(duì)核筒懸掛隔震結(jié)構(gòu)減震效果影響［J］.工業(yè)建筑，2013，43(2):34-38，58.