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

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

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

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

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

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

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

0　引言

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

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

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

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

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

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

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

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

2　結果與分析

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

3　結論

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

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

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

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

［1］劉卡丁，張玉軍.層狀巖體剪切破壞面方向的影響因素［J］.巖石力學與工程學報，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］.巖土力學， 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］王云飛，黃正均，崔芳.煤巖破壞過程的細觀力學損傷演化機制［J］.煤炭學報，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.

（責任編輯：鄧光輝）

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

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

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