王崔林,許 強,劉文德
(1.成都理工大學(xué),四川成都610059;2.地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護國家重點實驗室,四川成都610059)
材料在受力時發(fā)生的變形或裂紋的擴展會以彈性波的形式釋放應(yīng)變能的現(xiàn)象稱為聲發(fā)射。巖石試件在不同的加載階段都有聲發(fā)射產(chǎn)生,且在不同的階段有不同的聲發(fā)射特征[1]。聲發(fā)射技術(shù)現(xiàn)已廣泛被用來研究巖石、混凝土等脆性材料,它能連續(xù)、實時地監(jiān)測材料內(nèi)部微裂紋的萌生、擴展和貫通,進而研究材料漸進失穩(wěn)破壞的全過程[2],為監(jiān)測巖體動力災(zāi)害提供依據(jù)。
分形這一名詞是Benoit Mandelbrot在20世紀70年代為了表征復(fù)雜圖形和復(fù)雜過程首次將拉丁文Fractus轉(zhuǎn)化后引入自然科學(xué)領(lǐng)域的,其在之后的幾年的研究中將分形理論及其應(yīng)用推向一個全新的階段[3- 4];Himura Naoto、Sato Kazushi等采用彈性模型模擬了非均勻材料的破壞,通過與巖石聲發(fā)射試驗對比分析發(fā)現(xiàn),巖石在破壞過程中聲發(fā)射b值與分形維值D變化保持一致[5];Biancolinia等研究了在疲勞狀態(tài)下裂縫形成和擴展時所產(chǎn)生的聲發(fā)射事件,同時進行盒維數(shù)的計算,發(fā)現(xiàn)分維值可以很好地描述巖石加載過程中裂紋的動態(tài)演化情況[6];趙凱華、朱照宣、黃畇等在80年代末將分形理論及其應(yīng)用引入到國內(nèi),在國內(nèi)掀起了分形教學(xué)和研究的高潮[7]。謝和平將巖石力學(xué)和分形幾何相結(jié)合,應(yīng)用于巖石破裂全過程的分析中,建立了巖石分形理論[8]。至此,國內(nèi)巖石分形理論的研究已取得了一定的成果。趙奎等分析了單軸加載條件下巖石聲發(fā)射試驗Kaiser點及其相鄰點聲發(fā)射能量的關(guān)聯(lián)分維數(shù)的特征,得到了 Kaiser 點特征頻帶能量百分比大于相鄰點的重要結(jié)論[9];尹賢剛、李庶林、唐海燕等研究了巖石聲發(fā)射平靜期分形維數(shù)的變化,并對平靜期聲發(fā)射參數(shù)進行了量化[10];梁忠雨等分析了大理巖和紅砂巖在單軸壓縮條件下聲發(fā)射參數(shù)的分形特征,得出了巖石聲發(fā)射過程參數(shù)的分形特征具有一定的尺度范圍[11];柴肇云、康天合等通過對泥巖進行變角剪切壓模試驗,研究泥巖壓剪破壞裂隙演化規(guī)律和破裂塊體分布的分形特征,得到了隨著剪切角的增大,破裂塊體分布的分形維數(shù)值呈對數(shù)關(guān)系遞減的結(jié)論[12];劉京紅等運用分形理論建立了聲發(fā)射參數(shù)分形維數(shù)計算模型,對沖擊傾向性煤巖單軸壓縮聲發(fā)射試驗信號分形特征進行了研究,可有效預(yù)報煤巖破裂所引起的巖爆等煤巖動力災(zāi)害現(xiàn)象[13]。
本文對灰?guī)r進行了單軸壓縮及劈裂加載試驗,研究了灰?guī)r在這2種加載模式下漸進破壞全過程中聲發(fā)射事件的變化,并用分形理論分析了巖石聲發(fā)射參數(shù)序列分形特征,對比分析了灰?guī)r在2種不同加載條件下全過程不同應(yīng)力比的聲發(fā)射事件數(shù)的分形維數(shù)特征,為巖石聲發(fā)射動力災(zāi)害監(jiān)測提供依據(jù)。
取青川大光包滑坡的灰?guī)r,按照試驗內(nèi)容和國家試驗規(guī)范的要求加工成標準圓柱形試件。單軸壓縮試驗試件尺寸大小100 mm×50 mm(高度×直徑),試件端面的平整度在0.02 mm以內(nèi),共3個;劈裂試驗試件尺寸大小50 mm×50 mm(高度×直徑),試件端面的不平整度不超過0.1 mm,共3個。對試樣編號后,量測試樣的高度、直徑、質(zhì)量等物理參數(shù),統(tǒng)計結(jié)果見表1。
本次試驗采用MTS815 Flex Test GT巖石力學(xué)試驗設(shè)備對灰?guī)r進行單軸壓縮和巴西劈裂試驗,聲發(fā)射測試和采集由PCI-2 AE實時三維定位監(jiān)測系統(tǒng)完成。聲發(fā)射傳感器頻率為200 kHz,最大信號可達100 dB,動態(tài)范圍大于85 dB。本試驗中,設(shè)置門檻電壓為0.5 V和放大器增益為45 dB時,可達到最佳效果,用八通道聲發(fā)射儀記錄聲發(fā)射時間序列參數(shù)和原始波形數(shù)據(jù)。聲發(fā)射試驗檢測系統(tǒng)見圖1。
圖1 聲發(fā)射監(jiān)測系統(tǒng)
圖2 聲發(fā)射事件數(shù)-應(yīng)力與時間關(guān)系
本文室內(nèi)試驗對灰?guī)r試樣分別進行了單軸壓縮和巴西劈裂聲發(fā)射特性參數(shù)試驗,測試巖體受力變形過程中的聲發(fā)射特性,選取典型試樣D-3和P-2進行分析。通過記錄的數(shù)據(jù)繪出應(yīng)力-應(yīng)變曲線,并與聲發(fā)射事件數(shù)擬合,分析巖石破壞不同階段所對應(yīng)的聲發(fā)射特性,研究巖石在破壞過程中的聲發(fā)射活動特性。圖2是灰?guī)r單軸壓縮試驗和巴西劈裂試驗過程中的聲發(fā)射事件數(shù)-應(yīng)力與時間關(guān)系。圖3是聲發(fā)射能量-累積能量-應(yīng)力與時間關(guān)系。
圖3 聲發(fā)射能量-累積能量-應(yīng)力與時間關(guān)系
由圖2、3可知,相比于劈裂加載,單軸壓縮加載過程中灰?guī)r的聲發(fā)射活動在整個加載過程中都較為強烈。加載初期的聲發(fā)射活動是由于單軸壓縮或劈裂加載開始時巖石和加載板之間的摩擦而產(chǎn)生的,此后聲發(fā)射活動趨于穩(wěn)定,聲發(fā)射事件數(shù)和能量基本在很小的范圍變化。隨著荷載的進一步增大,單軸壓縮過程中聲發(fā)射事件數(shù)和能量都出現(xiàn)了3次較大的波動,反映了試樣內(nèi)部3次較大的裂紋擴展。隨著巖石進入破壞階段,單軸壓縮過程中巖石載荷能力下降,聲發(fā)射活動也隨之減少,但聲發(fā)射事件數(shù)和能量依然維持在很高的水平,這反映了灰?guī)r在壓縮破壞之后,內(nèi)部應(yīng)力還沒有馬上達到平衡,故而聲發(fā)射活動還未結(jié)束;但在劈裂加載過程中,達到峰值強度過后,巖石喪失載荷能力,聲發(fā)射活動、聲發(fā)射事件數(shù)和能量同時消失,這反映了灰?guī)r在劈裂破壞時,巖樣內(nèi)部主破裂面貫通,巖石立即失去了強度。
劈裂加載條件下聲發(fā)射的活動性相比于單軸壓縮條件下弱很多,且它們各變形破壞階段的聲發(fā)射特性也不同。因此,巖石聲發(fā)射的特性可被用來定性分析巖石的破壞模式,探測巖石的微觀損傷演化,預(yù)測判斷現(xiàn)場工程巖體的宏觀斷裂失穩(wěn)過程。
用分形理論研究巖石在整個加載過程中的破裂問題已經(jīng)取得一定的成果,分形維數(shù)在整個加載過程中的大小,可以反映巖石破裂的尺度。用分形理論探究灰?guī)r在單軸壓縮和劈裂加載條件下分形維數(shù)的變化,以此來反演在2種不同加載條件下的巖石內(nèi)部的破裂,結(jié)合宏觀巖石破壞形態(tài),能更加清楚地了解巖石的破壞過程。
分形維數(shù)是分形理論的基本量,其中關(guān)聯(lián)維數(shù)是最常用的分形維數(shù)之一。1983年,Grassberger和Procaccia根據(jù)嵌入理論和重構(gòu)相空間思想,提出了從時間序列直接計算關(guān)聯(lián)維數(shù)D的G-P算法[14]。該算法將聲發(fā)射基本參數(shù)序列作為研究對象,任何1個聲發(fā)射基本參數(shù)序列對應(yīng)1個容量為n的序列集,即
X={x1,x2,…,xn}
(1)
式中,X為容量為n的序列集。
該序列集合可以構(gòu)成1個m維的相空間(m X1={x1,x2,…,xm} (2) 然后右移1個數(shù)據(jù)再取m個數(shù)構(gòu)造第2個向量,依次類推,構(gòu)成N=n-m+1個向量,它們的關(guān)聯(lián)函數(shù)為C(r),即 (3) 圖5 不同加載應(yīng)力比的關(guān)聯(lián)維數(shù)D 圖4 相空間維數(shù)m與關(guān)聯(lián)維數(shù)D的關(guān)系 對每1個給定的尺度r都可得到1個C(r)。在雙對數(shù)坐標系中可得到一系列點{lnC[r(k)],lnr(k)},對這些點進行數(shù)據(jù)擬合,若回歸結(jié)果為直線,則表明聲發(fā)射序列在給定的尺度r范圍內(nèi)具有分形特征,直線的斜率就是聲發(fā)射參數(shù)序列的關(guān)聯(lián)維數(shù)D,即 (4) 本文在MATLAB中實現(xiàn)了該算法的程序編寫,能對巖石損傷破壞過程的聲發(fā)射時間序列數(shù)據(jù)進行歸一化處理、相空間重構(gòu)計算和關(guān)聯(lián)維數(shù)的計算。 重構(gòu)相空間維數(shù)m在計算關(guān)聯(lián)維數(shù)時的取值不同,得到的關(guān)聯(lián)維數(shù)D也不同,因此在計算巖石不同加載階段聲發(fā)射基本參數(shù)序列的關(guān)聯(lián)維數(shù)時,應(yīng)采用相同的相空間維數(shù)。用MATLAB計算得到的相空間維數(shù)m和關(guān)聯(lián)維數(shù)D的關(guān)系見圖4。從圖4可知,m取值為2~6時,關(guān)聯(lián)維數(shù)大致呈線性增加,意味著關(guān)聯(lián)維數(shù)梯度接近于穩(wěn)定;m>6時,關(guān)聯(lián)維數(shù)D隨相空間維數(shù)m的增加基本不變。本文取相空間維數(shù)m=4。 圖6 不同加載應(yīng)力比的關(guān)聯(lián)維數(shù)D 選取典型單軸灰?guī)r試樣D-3的聲發(fā)射計數(shù)時間序列,計算在加載應(yīng)力比為(0~0.2)σc、(0.2~0.4]σc、(0.4~0.6)σc、(0.6~0.8)σc、(0.8~0.9)σc和(0.9~1.0)σc(σc為峰值荷載)時的關(guān)聯(lián)維數(shù)D,結(jié)果見圖5。 從圖5可知,在單軸壓縮加載的整個過程中,關(guān)聯(lián)維數(shù)D整體上呈下降趨勢。在(0~0.2)σc應(yīng)力區(qū)間內(nèi),由于剛性試驗機和試樣表面發(fā)生橫向摩擦和初始裂紋被壓密,而導(dǎo)致有大量聲發(fā)射事件產(chǎn)生,所以此時關(guān)聯(lián)維數(shù)D不是最大的,而后進入彈性階段,沒有新的裂紋產(chǎn)生,巖石內(nèi)部沒有破裂,幾乎沒有聲發(fā)射事件產(chǎn)生,聲發(fā)射關(guān)聯(lián)維數(shù)D達到最大值,而后開始出現(xiàn)下降走勢,在(0.4~0.6)σc應(yīng)力區(qū)間內(nèi)下降得最快,這意味著聲發(fā)射大事件所占的比例增加,巖石試樣內(nèi)較大尺度的微破裂增多,巖石內(nèi)部裂紋擴展從無序逐漸向有序發(fā)展。隨著載荷的繼續(xù)增加,試樣破裂聲發(fā)射關(guān)聯(lián)維數(shù)D持續(xù)下降,巖石內(nèi)部損傷程度增加,在破壞前均下降到最小值,表明此時巖石內(nèi)部裂紋逐漸沿著主破裂面連接貫通,并最終形成宏觀破壞面導(dǎo)致巖石失穩(wěn)破壞。 選取典型劈裂灰?guī)r試樣P-2的聲發(fā)射計數(shù)時間序列,計算巖石劈裂全過程中加載應(yīng)力比為(0~0.2)σt、(0.2~0.4)σt、(0.4~0.6)σt、(0.6~0.8)σt、(0.8~0.9)σt和(0.9~1.0)σt(σt為劈裂破壞荷載)時的關(guān)聯(lián)維數(shù)D,結(jié)果見圖6。 從圖6可知,在劈裂加載的整個過程中,關(guān)聯(lián)維數(shù)D整體上呈下降趨勢。在劈裂加載初期(0~0.2)σt應(yīng)力區(qū)間內(nèi),由于鋼絲墊條與巖石、鋼絲墊條與剛性試驗機之間存在摩擦,聲發(fā)射活動較強,所以此時的關(guān)聯(lián)維數(shù)D不是最大的,在(0~0.2)σt~(0.2~0.4)σt應(yīng)力區(qū)間內(nèi),隨著載荷的增加,關(guān)聯(lián)維數(shù)D逐漸增大,這是因為在較小應(yīng)力水平時,巖石試樣內(nèi)部的破裂是以小尺度的微破裂為主,巖石試樣初始裂紋分布較為均勻,聲發(fā)射小事件所占的比例較多且在不斷增加,巖石處于相對穩(wěn)定階段。當載荷增加到一定程度時,聲發(fā)射關(guān)聯(lián)維數(shù)D達到最大值,而后開始出現(xiàn)下降走勢,這意味著聲發(fā)射活動開始變強,巖石試樣內(nèi)較大尺度的微破裂逐漸增多,且?guī)r石內(nèi)裂紋擴展逐漸從無序向有序發(fā)展。在破壞前((0.9~1.0)σt)關(guān)聯(lián)維數(shù)D下降到最小值,表明巖石內(nèi)部微破裂在劈裂作用下逐漸連接貫通,并形成最終宏觀破壞面導(dǎo)致巖石失穩(wěn)破壞。 不同加載模式下分形維數(shù)D變化見圖7。從圖7可知,劈裂加載條件下聲發(fā)射分形維數(shù)相較于單軸加載條件下普遍大一些,表明了單軸加載條件下聲發(fā)射活動性強于劈裂加載,反映了劈裂加載條件下巖石的破裂尺度較小,裂紋的擴展分布較均勻;2種加載模式下聲發(fā)射分形維數(shù)D整體上均呈現(xiàn)下降趨勢,且在巖石試樣破壞前下降到最低值,這反映了巖石破裂的分維規(guī)律性,表明了巖石失穩(wěn)破壞是一個降維有序的過程,可將其作為巖體失穩(wěn)破壞的前兆,從而有效預(yù)測判斷現(xiàn)場工程巖體的宏觀斷裂失穩(wěn)過程。 圖7 不同加載模式下分形維數(shù)D變化 本文通過對灰?guī)r在單軸壓縮和劈裂加載過程中聲發(fā)射分形特性試驗的研究,得出以下結(jié)論: (1)單軸壓縮條件下聲發(fā)射的活動性相比于劈裂加載條件下強很多,且兩者各變形破壞階段的聲發(fā)射特性也不同。因此,巖石聲發(fā)射的特性可被用來定性分析巖石的破壞模式、分析巖石的破壞演化的各個階段、預(yù)測判斷現(xiàn)場工程巖體的宏觀斷裂失穩(wěn)過程。 (2)劈裂加載條件下的聲發(fā)射分形維數(shù)普遍大于單軸加載,表明了單軸加載條件下聲發(fā)射活動性強于劈裂加載,反映了劈裂加載條件下巖石的破裂尺度較小,裂紋的擴展分布較均勻。 (3)不同加載條件下的破裂過程中聲發(fā)射關(guān)聯(lián)維數(shù)D整體上均呈現(xiàn)下降趨勢,在巖石試樣破壞前下降到最低值,這表明巖石失穩(wěn)破壞是一個降維有序的過程,可將其作為巖體失穩(wěn)破壞的前兆,從而有效預(yù)測判斷現(xiàn)場工程巖體的宏觀斷裂失穩(wěn)過程。3.2 相空間維數(shù)的確定
3.3 單軸壓縮條件下分形維數(shù)演化特征
3.4 劈裂條件下分形維數(shù)演化特征
4 結(jié) 語