單軸壓縮下橫觀各向同性巖石破裂過程聲發(fā)射特性的離散元模擬

劉京鑠，范金星，李　娟，劉　斌

（湖南水利水電職業(yè)技術(shù)學(xué)院 水利水電工程系，湖南 長沙 410131）

單軸壓縮下橫觀各向同性巖石破裂過程聲發(fā)射特性的離散元模擬

劉京鑠，范金星，李娟，劉斌

（湖南水利水電職業(yè)技術(shù)學(xué)院 水利水電工程系，湖南 長沙 410131）

采用離散單元法研究了橫觀各向同性巖石在不同層面傾角條件下的單軸壓縮破壞過程及聲發(fā)射特性。結(jié)果表明：由于層面傾角的變化，導(dǎo)致橫觀各向同性巖石破裂過程具有不同的聲發(fā)射特性。聲發(fā)射特性與應(yīng)力存在一定的耦合關(guān)系，且聲發(fā)射空間響應(yīng)集中在2種巖石的交界面上。互層巖體由于層間力學(xué)屬性不同，極有可能在交界面處引起應(yīng)力集中，使得在交界面附近裂紋最先萌生，裂紋進一步向交界面兩側(cè)巖體中擴展，從而引起宏觀裂紋的出現(xiàn)，最終引起巖體的破壞。隨著層面傾角的增大，巖石的單軸抗壓強度和彈性模量呈先減小后增大，層面傾角90°時的強度甚至超過了0°時的強度。

橫觀各向同性巖石；軟硬互層；層面傾角；單軸壓縮；聲發(fā)射；離散單元法

0　引言

層狀巖體是廣泛存在的地質(zhì)材料，其各向異性性狀對巖體的力學(xué)行為有重要的影響。巖石彈塑性變形只是各種微結(jié)構(gòu)之間相互作用的宏觀體積響應(yīng)。變形的主要機理是微裂紋成核和增長及其聚集而生成細(xì)觀裂紋和裂紋的擴展與傳播［1］。但是，由于巖體中包含的微觀顆粒數(shù)目龐大，表現(xiàn)形態(tài)各異，要從理論上精確地描述和分析各個損傷實體，幾乎是不可能的。隨著計算機技術(shù)的發(fā)展，數(shù)值模擬手段越來越多地應(yīng)用到巖石力學(xué)研究中。

巖體內(nèi)部微缺陷和微裂紋等的成核、擴展并快速釋放能量，是巖體損傷、失效與破壞的主要機制。聲發(fā)射技術(shù)可以實時監(jiān)測巖石的微觀損傷和演化過程，是一種行之有效的研究手段［2-8］。離散元通過顆粒間的接觸和粘結(jié)來儲存巖體彈性應(yīng)變能。在應(yīng)力作用下，當(dāng)顆粒間平行黏結(jié)的應(yīng)力大于其黏結(jié)強度時，黏結(jié)發(fā)生破壞斷裂，并釋放應(yīng)變能。所以，采用離散元模擬聲發(fā)射行為是可行的［9］。

本文借助離散元軟件，用2種不同的巖石材料來組成不同層面傾角的橫觀各向同性的巖石試件，通過單軸加載數(shù)值模擬試驗，模擬橫觀各向同性巖石破裂的整個過程，結(jié)合聲發(fā)射監(jiān)測技術(shù)，研究單軸壓縮條件下橫觀各向同性巖石試件的力學(xué)特性。

1　離散元數(shù)值分析

離散元細(xì)觀參數(shù)的標(biāo)定是：以巖石常規(guī)單軸壓縮應(yīng)力-應(yīng)變數(shù)值模擬曲線和室內(nèi)試驗曲線吻合為細(xì)觀參數(shù)標(biāo)定的準(zhǔn)則；基本過程是，不斷調(diào)整顆粒法向、切向接觸剛度，同時控制其峰值強度在一定的范圍內(nèi)，使顆粒流的細(xì)觀數(shù)值模擬得到的結(jié)果逼近材料的宏觀力學(xué)響應(yīng)。細(xì)觀單元力學(xué)參數(shù)間雖相互獨立，而又同時影響著試件的某種宏觀力學(xué)響應(yīng)。與材料強度有關(guān)的細(xì)觀參數(shù)有接觸剛度、黏結(jié)強度和摩擦因數(shù)，同時變化這些參數(shù)來得到宏觀強度變化規(guī)律是非常困難的，目前尚無有效的解決方法。因此，選擇顆粒細(xì)觀結(jié)構(gòu)、力學(xué)參數(shù)與物理模型的過程是非常復(fù)雜而又困難的，有一個大量試錯的過程。

表1　介質(zhì)顆粒參數(shù)Table 1　Parameters of medium particles

表2　介質(zhì)黏結(jié)參數(shù)Table 2　Bonding parameters of mediums

圖1　不同傾角的數(shù)值模擬試件示意圖Fig.1　A diagram of numerical simulation specimens with different orientation angles

2　結(jié)果與分析

為了探討橫觀各向同性巖石的力學(xué)性質(zhì)，利用不同層面傾角的橫觀各向同性巖石試件進行單軸壓縮試驗?zāi)M。圖2為不同層面傾角試件的應(yīng)力-應(yīng)變和累計AE數(shù)曲線。由圖可知，不同層面傾角試件的應(yīng)力-應(yīng)變和累計AE數(shù)曲線的形態(tài)基本一致。應(yīng)力-應(yīng)變曲線可分為5個階段。

1）初始壓密階段：應(yīng)力-應(yīng)變關(guān)系曲線微向上彎曲。此階段幾乎沒有裂紋產(chǎn)生，累計AE數(shù)維持在較低水平。

2）彈性階段：應(yīng)力-應(yīng)變曲線近似于直線，隨著荷載的增加，巖石內(nèi)部開始產(chǎn)生微裂紋。此階段為裂紋的穩(wěn)定擴展段，對應(yīng)累計AE數(shù)穩(wěn)定增加。

3）應(yīng)變硬化階段：此階段為裂紋的加速擴展階段，應(yīng)力-應(yīng)變曲線呈顯著非線性，隨著法向荷載的繼續(xù)增加，黏結(jié)斷裂后相應(yīng)的受力構(gòu)架會發(fā)生應(yīng)力調(diào)整，導(dǎo)致臨近單元應(yīng)力增加，臨近單元應(yīng)力水平達(dá)到黏結(jié)強度時隨之破壞，應(yīng)力傳遞路徑不斷變化，受力構(gòu)架不斷削弱。同一時間破壞的黏結(jié)數(shù)急劇增加，表現(xiàn)為累計AE數(shù)急劇增加。這說明內(nèi)部破壞急劇產(chǎn)生，裂紋加速發(fā)展，逐漸出現(xiàn)一定的貫通趨向性。

4）應(yīng)變軟化階段：巖體承載力達(dá)到峰值強度后，其內(nèi)部結(jié)構(gòu)遭到破壞，但試件基本保持整體狀。到本階段，應(yīng)力迅速衰減，出現(xiàn)明顯的應(yīng)力降，裂紋快速發(fā)展，內(nèi)部損傷不斷累積，最終形成宏觀斷裂面，使試件發(fā)生失穩(wěn)破壞。累計AE數(shù)呈近線性增長趨勢。

5）殘余強度階段：試件變形主要表現(xiàn)為沿宏觀斷裂面的塊體滑移，試件承載力隨變形增大迅速下降，并降到0。此時，累計AE數(shù)增加速率降低并最終趨于穩(wěn)定。

圖2　不同層面傾角試件的應(yīng)力-應(yīng)變和累計AE數(shù)曲線Fig.2　Diagrams of stress-strain and accumulative AE curves of specimens with different orientation angles

從圖3還可以看出， 軸向應(yīng)變從0.2%增大到0.6%，累計AE數(shù)逐漸增加。而軸向應(yīng)變?yōu)?.6%時的累計AE數(shù)與軸向應(yīng)變?yōu)?.8%時的累計AE數(shù)幾乎相同。這說明軸向應(yīng)變?yōu)?.6%～0.8%已處于殘余強度階段，累計AE數(shù)增加速率降低并最終趨于穩(wěn)定。層面傾角從0°增大到75°時，不同軸向應(yīng)變時的累計AE數(shù)的差值減小。層面傾角90°時，差值增大。

圖3　不同軸向應(yīng)變時的累計AE數(shù)-層面傾角曲線Fig.3　A curve diagram of accumulative AE event-orientation angles at different axial strains

目前，基于聲發(fā)射特性的研究往往是從聲發(fā)射-應(yīng)變的角度進行分析，然而應(yīng)力是導(dǎo)致材料變形破壞的主動因素。而應(yīng)變、聲發(fā)射是應(yīng)力作用下的被動響應(yīng)。所以，聲發(fā)射與應(yīng)力之間的耦合關(guān)系對于揭示聲發(fā)射的潛在力學(xué)本質(zhì)更有幫助。圖4為加載過程中累計AE數(shù)-應(yīng)力曲線。由圖可知，不同層面傾角試件的聲發(fā)射特性與應(yīng)力關(guān)系都呈現(xiàn)出一定的規(guī)律，二者存在一定的耦合關(guān)系。

圖4　不同層面傾角試件的累計AE數(shù)-偏應(yīng)力曲線Fig.4　A curve diagram of accumulative AE event-deviator stresses of rock samples with different orientation angles

巖石的穩(wěn)定性與裂紋的分布關(guān)系密切。巖石工程的失穩(wěn)破壞往往是由于原有裂紋的演化、擴展和貫通造成的，微裂紋成核時瞬間釋放的彈性應(yīng)變能是聲發(fā)射最主要的組成部分。聲發(fā)射空間響應(yīng)對應(yīng)著微裂紋的空間分布。不同層面傾角條件下巖石試件破壞時的裂紋分布如圖 5所示。由圖可知，試件的聲發(fā)射空間響應(yīng)集中在材料交界面上。交界面上的微裂紋開始貫通成較大的宏觀斷裂面，裂紋首先出現(xiàn)在交界面處。由于橫觀各向同性巖石的非連續(xù)性與非均質(zhì)性，在2種材料的交界面處形成復(fù)雜的應(yīng)力-應(yīng)變關(guān)系，極有可能在交界面處引起應(yīng)力集中，導(dǎo)致裂紋最有可能產(chǎn)生于交界面處。且層面傾角的存在使得交界面處存在剪切應(yīng)力。所以，交界面就成為裂紋優(yōu)先萌生和擴展的位置［10-11］。

圖5　不同層面傾角條件下試件破壞時的聲發(fā)射空間響應(yīng)Fig.5　Spatial responses of acoustic emissions at the failure for specimens with different orientation angles

圖6為不同傾角試件的單軸抗壓強度和彈性模量關(guān)系曲線。由圖可知，隨著層面傾角的增大，巖石的單軸抗壓強度先減小后增大。層面傾角從0°增大到30°，單軸抗壓強度降幅較小。=15°和30°時，單軸抗壓強度從14.67 MPa分別減小到13.40 MPa和12.74 MPa，降幅分別為8.67%和13.16%。之后單軸抗壓強度加速減小，=60°時，巖石試件的單軸抗壓強度減小到最小值4.53 MPa。隨后，單軸抗壓強度隨著層面傾角的增大而增大。當(dāng)=75°時，單軸抗壓強度增大到5.06 MPa，并在=90°時達(dá)到最大值，甚至高于=0°時的強度，這也說明層狀巖體具有強烈的各向異性。相比于單軸抗壓強度，彈性模量的變化范圍較小，尤其是層面傾角在0°～45°的范圍內(nèi)，彈性模量隨著層面傾角的增大而單調(diào)緩慢減小。說明層面傾角對彈性模量的影響較小。雖然單軸抗壓強度是在=60°時最小，但彈性模量卻是在=45°時最小。=75°時單軸抗壓強度較小，但彈性模量卻是最大。這說明了橫觀各向同性巖石，在彈性階段所表現(xiàn)出的彈性模量的大小并不能代表其最終峰值強度的大小。

由以上各傾角試件的單軸壓縮模擬試驗獲得的應(yīng)力-應(yīng)變關(guān)系發(fā)現(xiàn)，低傾角軸向應(yīng)變值有時會比較大，弱相巖石B的彈性模量較低、強度也較?。ㄝ^易發(fā)生破壞），因此低傾角試件的力學(xué)行為主要由弱相巖石B主導(dǎo)。中傾角試件（為45°～75°）破壞形式主要為沿交界面產(chǎn)生滑動破壞，試件在受到軸向應(yīng)力時，交界面剪力強度承受不住剪應(yīng)力的持續(xù)增大導(dǎo)致交界面破壞。因此，中傾角試件的力學(xué)行為主要由交界面性質(zhì)所主導(dǎo)。高傾角（=90°）試件由巖石材料A和B一同承受載荷，強度較高的巖石材料A會分擔(dān)更多的應(yīng)力，故高傾角試件的力學(xué)行為由巖石材料A所主導(dǎo)的比例較多。

圖6　單軸抗壓強度和彈性模量與層面傾角的關(guān)系曲線Fig.6　A contrast diagram of UCS andelastic modulus with different orientation angles

3　結(jié)論

1）由于層面傾角的變化，導(dǎo)致橫觀各向同性巖石破裂過程具有不同的聲發(fā)射特性，且聲發(fā)射空間響應(yīng)集中在材料交界面上。

2）層狀巖體由于層間力學(xué)屬性不同，極有可能在交界面處引起應(yīng)力集中，使得在交界面附近裂紋最先萌生，裂紋進一步向交界面兩側(cè)巖體中擴展，從而引起宏觀裂紋出現(xiàn)，最終導(dǎo)致巖體破壞。

3）應(yīng)力是導(dǎo)致聲發(fā)射事件產(chǎn)生的主動因素，聲發(fā)射特性與應(yīng)力存在一定的耦合關(guān)系。

4）隨著層面傾角的增大，巖石的單軸抗壓強度和彈性模量呈先減小后增大的變化趨勢。

［1］劉卡丁，張玉軍.層狀巖體剪切破壞面方向的影響因素［J］.巖石力學(xué)與工程學(xué)報，2002，21（3）：335-339.LIU Kading，ZHANG Yujun.Influence Factors on Shear Failure Orientation of Layered Rocks［J］.Chinese Journal of Rock Mechanics and Engineering，2002，21（3）：335-339.

［2］ BAUD P，KLEIN E，WONG T F.Compaction Localization in Porous Sandstones: Spatial Evolution of Damage and Acoustic Emission Activity［J］.Journal of Structural Geology，2004，26（4）：603-624.

［3］INGRAHAM M D，ISSEN K A，HOLCOMB D J.Use of Acoustic Emissions to Investigate Localization in High-Porosity Sandstone Subjected to True Triaxial Stresses［J］.Acta Geotechnica，2013，8（6）：645-663.

［4］CHITRALA Y，MORENO C，SONDERGELD C，et al.An Experimental Investigation into Hydraulic Fracture Propagation Under Different Applied Stresses in Tight Sands Using Acoustic Emissions［J］.Journal of Petroleum Science and Engineering，2013，108：151-161.

［5］ISHIDA T.Acoustic Emission Monitoring of Hydraulic Fracturing in Laboratory and Field［J］.Construction and Building Materials，2001，15（5/6）：283-295.

［6］DONG L J，LI X B，ZHOU Z L，et al.Three-Dimensional Analytical Solution of Acoustic Emission Source Location for Cuboid Monitoring Network Without Pre-Measured Wave Velocity［J］.Transactions of Nonferrous Metals Society of China，2015，25（1）：293-302.

［7］PRADHAN S，STROISZ A M，F(xiàn)J R E，et al.Stress-Induced Fracturing of Reservoir Rocks: Acoustic Monitoring and CT Image Analysis［J］.Rock Mechanics and Rock Engineering，2015，48（6）：2529-2540.

［8］穆康，俞縉，李宏，等.水-力耦合條件下砂巖聲發(fā)射和能量耗散的顆粒流模擬［J］.巖土力學(xué)， 2015， 36（5）：1496-1504.MU Kang，YU Jin，LI Hong，et al.Acoustic Emission of Sandstone with Hydro-Mechanical Coupling and PFC-Based Modelling of Energy Dissipation［J］.Rock and Soil Mechanics，2015，36（5）：1496-1504.

［9］王云飛，黃正均，崔芳.煤巖破壞過程的細(xì)觀力學(xué)損傷演化機制［J］.煤炭學(xué)報，2014，39（12）：2390-2396.WANG Yunfei，HUANG Zhengjun，CUI Fang.Damage Evolution Mechanism in the Failure Process of Coal Rock Based on Mesomechanics［J］.Journal of China Coal Society，2014，39（12）：2390-2396.

［10］LI Y P，LIU W，YANG C H，et al.Experimental Investigation of Mechanical Behavior of Bedded Rock Salt Containing Inclined Interlayer［J］.International Journal of Rock Mechanics and Mining Sciences，2014，69：39-49.

［11］LI Y P，CHEN L Z，WANG Y H.Experimental Research on Pre-Cracked Marble Under Compression［J］.International Journal of Solids and Structures，2005，42（9/10）：2505-2516.

（責(zé)任編輯：鄧光輝）

DEM Simulation of Acoustic Emission Characteristics in the Failure Process of Transversely Isotropic Rocks Under Uniaxial Compression

LIU Jingshuo，F(xiàn)ANG Jinxing， LI Juan，LIU Bin（Department of Hydraulic Engineering，Hunan Polytechnic of Water Resources and Electric Power，Changsha 410131，China）

A research has been conducted， with the discrete element method introduced， on the failure process under uniaxial compression and the acoustic emission characteristics of transversely isotropic rocks with different orientation angles.The experimental results show that the difference in orientation angle of transversely isotropic rocks results in various acoustic emission characteristics.A coupling relation has been found between acoustic emission characteristics and the corresponding stresses， with the spatial response of acoustic emissions occurring mainly at the interface of two kinds of rocks.Differences between mechanical properties of transversely isotropic rocks result in a stress concentration at the interface of inter-stratified rock mass where a crack initiation may first occur， which expands into the nearby area and then causes macroscopic cracks to form， finally resulting in the total failure of the entire specimen.With the increase of orientation angles， the uniaxial compression strength and modulus of elasticity decrease initially and then increase.The compression strength reaches its peak when=90°， where it exceeds the corresponding compression strength when=0° .

transversely isotropy rock；inter-layered rock；orientation angle；uniaxial compression；acoustic emission （AE）； discrete element method

TU452

A

1673-9833（2016）04-0005-05

10.3969/j.issn.1673-9833.2016.04.002

2016-06-20

國家自然科學(xué)基金資助項目（51174228），湖南省教育廳科學(xué)研究基金資助項目（14C0746，15C0895）

劉京鑠（1979-），男，湖南邵東人，湖南水利水電職業(yè)技術(shù)學(xué)院講師，中南大學(xué)博士生，主要研究方向為巖土工程和水工結(jié)構(gòu)，E-mail：ljs7327@163.com