單軸應力下的孔洞–裂隙擴展規(guī)律數值模擬

郭修成

單軸應力下的孔洞–裂隙擴展規(guī)律數值模擬

郭修成

(長城鉆探工程有限公司工程技術研究院，遼寧 盤錦 124010)

為研究不同孔洞–裂隙(簡稱“孔–隙”)賦存條件下的裂紋擴展規(guī)律，利用RFPA軟件，對不同裂紋傾角及不同非均質系數下的巖體破壞進行數值模擬分析，獲得其裂紋擴展過程、聲發(fā)射規(guī)律、應力–應變曲線，同時與原試驗結果進行對比驗證。結果表明：完整試樣裂紋沿著剪切方向產生，含孔–隙試樣裂紋沿裂隙尖端及孔口側邊產生；翼裂紋貫穿試件的同時，在預制裂紋尖端或孔口側邊產生水平方向的次生裂紋，并產生分叉，非均質系數影響次生裂紋走向；壓載前期試樣以拉破壞為主，壓載后期以拉–剪組合破壞為主，次生裂紋的產生與剪切破壞有關；聲發(fā)射累計能量與聲發(fā)射累計數前期緩慢增大，后期迅速增大，預制裂紋傾角越小，非均質系數越大，聲發(fā)射累計能量越大；不同裂紋傾角及不同非均質系數試件的應力–應變曲線均經歷3個階段：彈性變形階段、非線性變形階段及殘余變形階段，孔–隙的存在降低了試樣的峰值強度，影響試件的脆性度。研究結果為進一步認識孔–隙相互作用規(guī)律提供了參考。

單軸應力；孔–隙作用；聲發(fā)射；裂紋擴展；應力–應變曲線；數值模擬

孔洞、裂隙廣泛存在于巖體內部[1-3]，如深埋隧洞在孔洞周圍存在大量的裂隙，在復雜應力下裂隙與裂隙、裂隙與孔洞及孔洞與孔洞之間會擴展貫通[4-6]，當達到極限狀態(tài)時，巖體結構將發(fā)生災變，因此，對巖體缺陷的力學規(guī)律研究將有益于認識及防治巖體災變。

國內外對含缺陷(包括裂隙及孔洞)巖體的力學規(guī)律研究較多，陳雪峰等[6]開展了不同厚度節(jié)理爆破模型試驗，得到節(jié)理模型的裂紋擴展規(guī)律；王飛等[7]對平行節(jié)理相互作用下巖體的力學性能進行了室內實驗；岳中文等[8]對含切槽炮孔的脆性巖體在爆破作用下的動態(tài)裂紋擴展進行了研究；李竟艷等[9]基于Abaqus軟件對動態(tài)拉伸載荷下巖石材料泛形裂紋擴展進行了數值模擬研究；李宏等[10]基于離散元顆粒流程序PFC對含非貫通型表面裂隙組的試件進行了分析，得到含多組裂隙相互貫通擴展的規(guī)律；周喻等[11]對含不同方位角的雙圓孔類巖石試件進行了單軸壓縮試驗，得到雙孔的相互作用規(guī)律。以往研究均僅針對裂隙與裂隙之間或孔與孔之間的裂紋擴展規(guī)律，對裂隙與孔之間裂紋的擴展規(guī)律研究較少。

本文利用RFPA軟件，采用數值模擬方法，對不同非均質系數、不同方位角下的裂紋擴展過程、應力–應變曲線及聲發(fā)射規(guī)律進行分析，將模擬結果與原試驗結果進行對比，驗證數值模擬結果的合理性，研究成果為認識孔–隙相互作用規(guī)律提供參考。

1 數值模擬計算理論

1.1 單軸壓縮下脆性材料的破壞準則

裂隙在受到雙軸應力作用下，Griffith理論強度[12-18](裂隙發(fā)展理論)可以寫成：

式中：1、3分別為最大、最小主應力，MPa；t為抗拉強度，MPa。

若巖體受到單向壓縮，則3=0，1=c，c表示單軸抗壓強度，MPa，可得：

將Griffith理論寫成Moh-Coulomb形式，可以將式(2)用正應力和剪應力表示，即：

式中：、分別為裂紋上的正應力與切應力，RFPA軟件正是根據這一理論進行模擬計算。

1.2 巖石介質非均勻性的統(tǒng)計理論

非均質系數可描述巖體材料的不均勻程度。對于同一種巖石材料來說，礦物晶體、膠結物晶體及各種微缺陷等具有各自的排列方式，相互之間結合強度有差異，導致其物理性質不能用同一特征值描述。1939年，Weibull率先提出用統(tǒng)計數學描述材料非均勻性的方法，Weibull統(tǒng)計分布函數表達式[11]可表示成：

式中：為巖石介質的基元體力學性質參數(強度、彈性模量等)；0為基元體力學性質的平均值；為分布函數的非均質系數，反映巖石介質的均質性；()為巖石基元體力學性質的統(tǒng)計分布密度。

本文對于巖石不均勻程度的描述采用本節(jié)非均質系數定義。

1.3 RFPA系統(tǒng)的基元狀態(tài)定義及計算流程

1.3.1 RFPA基元狀態(tài)定義

RFPA中的基元系統(tǒng)指在模型中當前功能的實體介質，它的本構關系由應力應變關系及破壞準則式(3)所描述。脆性介質最大特點之一是介質的變形達到某個閾值時，介質的力學性質將發(fā)生突變，即介質由于損傷而導致其承載力急劇下降，甚至完全喪失承載力，基元將發(fā)生相變。

RFPA中將壓、剪破壞引起的基元破壞稱為Ⅰ類相變，將拉應力引起的基元破壞稱為Ⅱ類相變。當基元介質發(fā)生相變后，RFPA引入空氣單元的概念，由于空氣單元的基元彈性模量極低，可近似認為實體介質的力學行為已不再存在，這樣，在不改變模型數學結構的前提下，使模型在總體特性上可以反映因基元破裂而引起的物理特性改變。值得注意的是，對于已形成的裂紋面(即空氣單元)，當裂隙兩面的介質在壓應力作用下出現壓密或壓實的現象時，將空氣單元改為接觸單元，使其剛度增大，起到傳遞應力的作用，這就是所謂的接觸單元特性。

1.3.2 RFPA計算流程

a. 建模及網格 可以根據用戶需要選擇基元(分為實體、支護、孔洞)，網格采用自動剖分技術。

b. 應力應變計算 根據輸入的邊界條件，進行有限元計算，輸出應力、應變及位移信息。

c. 相變分析 根據式(1)—式(3)判斷單元是否發(fā)生相變，進行基元的轉化處理。

2 裂紋擴展數值計算模型

2.1 模型及參數

為研究不同裂紋方位角及不同非均質系數下含孔洞裂隙試件在單軸壓力下的裂紋擴展規(guī)律，根據文獻[19]的試驗結果進行有限元建模計算。數值模擬采用二維試樣模型，寬60 mm，長120 mm，采用平面應力模型，圓孔布置于試件中心，圓孔半徑為5 mm，孔洞側邊裂隙長度為12 mm，裂隙傾角定義為裂隙方向與水平方向的夾角，計算模型如圖1a所示；不同材料非均質系數=1.1、3、5、10、50的彈性模量分布如圖1b—圖1f所示。試樣采用單軸壓縮應力加載，加載速率設置為0.004 mm/步；材料的基礎參數見表1。

圖1 試樣設計及不同非均值系數下的彈性模量分布

表1 材料基礎參數

2.2 模擬計算工況

為研究不同裂隙傾角及不同非均質系數下的孔–隙相互作用規(guī)律及對材料物理力學性質的影響，設置5個計算工況：①工況1，完整對比組，目的是與含孔洞及裂隙試件的對比研究；②工況2，不同裂隙傾角組，目的是對比不同裂隙傾角下試件物理力學特性差異；③工況3，不同非均質系數組，目的是對比不同非均質系數對試件物理力學特性的影響；④工況4，含單孔試樣組，目的是與多孔試樣進行對比分析；⑤工況5，含單裂隙試樣組，目的是與含孔–隙試樣進行對比分析，具體計算工況詳見表2。

3 模擬計算結果分析

3.1 裂紋擴展過程

3.1.1 工況1、工況2及工況4

完整試樣的裂紋擴展計算過程與含孔洞及裂隙試樣的裂紋擴展過程如圖2所示。由圖2可知，孔–隙的存在及裂紋傾角對試件破壞形態(tài)有較大影響，主要規(guī)律體現如下。

表2 不同裂隙傾角及非均質系數的計算工況

對于完整試樣來說，試件呈現典型的剪切破壞模式，即剪切破壞面沿試件約60°方向產生(圖2a)。對單含孔洞試樣而言，裂紋首先在孔洞上下出現，隨后出現剪切裂紋，使試樣呈現典型的剪切破壞模式(圖2b)。預制裂紋傾角較小時(=0°,15°,30°)，裂紋首先沿著孔口上下端產生，并持續(xù)擴展；當裂隙傾角為30°時，在預制裂紋下端產生翼裂紋，最終裂紋沿著最大主應力方向發(fā)展并貫穿試件(圖2c—圖2e)；值得注意的是，裂紋擴展過程中產生一定的分叉，同時，在裂紋擴展后期預制裂紋尖端產生次生裂紋，次生裂紋沿水平方向發(fā)展，并產生更大的裂紋分叉。當預制裂紋傾角較大時(=45°,60°)，翼裂紋首先沿預制裂紋尖端產生，并沿最大主應力發(fā)展，未在孔口上下頂端產生；最終翼裂紋貫穿試件的同時在預制裂紋尖端產生水平方向的次生裂紋，并產生更大的分叉(圖2f,圖2g)。當預制裂紋傾角很大時(=75°,90°)，翼裂紋產生于預制裂紋尖端，但次生裂紋產生于孔口左右側頂端(圖2h—圖2i)。

圖2 工況1、工況2與工況4裂紋擴展過程

3.1.2 工況3

不同非均質系數下的試件裂紋擴展過程如圖3所示。由圖3可知，不同非均質系數下的孔–隙裂紋擴展規(guī)律也有一定差異。

不同非均質系數下，翼裂紋首先在預制裂紋尖端處產生，并逐漸沿著最大主應力方向發(fā)展，非均質系數越大，裂紋分叉程度就越弱。不同非均質系數下的次生裂紋形態(tài)也有不同，當非均質系數較小時(=1.1)，次生裂紋形態(tài)為反翼裂紋，即次生裂紋擴展方向與翼裂紋擴展方向相反(圖3a)；而非均質系數較大時(=3,5,10,50)，次生裂紋擴展則主要沿著水平方向進行擴展(圖3b—圖3e)。

3.1.3 工況5

含單裂隙不同傾角下的試件裂紋擴展過程如圖4所示。由圖4可知，不同裂隙傾角下的裂紋擴展規(guī)律與含孔–隙試樣的裂紋擴展規(guī)律大有不同，裂隙的分叉效應更加明顯，當裂隙的傾角較小時(=0°,15°)，裂紋從預制裂隙的中部產生(圖4a)；當預制裂隙的傾角較大時(大于30°)，裂紋將從預制裂隙的尖端產生，隨后朝著最大主應力方向發(fā)展，最終使得試樣破壞。

圖3 不同非均質系數下的裂紋擴展過程

圖4 不同裂隙傾角下的裂紋擴展過程

3.2 聲發(fā)射規(guī)律分析

3.2.1 工況1、工況2及工況4

完整試件與含不同傾角裂隙及孔洞試件的聲發(fā)射過程及聲發(fā)射曲線如圖5—圖6所示。在RFPA中，白色圓圈代表發(fā)生剪切破壞，紅色圓圈代表發(fā)生拉伸破壞。

由圖5可知，聲發(fā)射規(guī)律與3.1節(jié)裂紋擴展規(guī)律呈現高度的一致性，主要得出以下規(guī)律。

對于完整試樣來說，聲發(fā)射事件首先無規(guī)律分布于試件的不同部位，這是由試件的非均質性導致；在壓載后期，沿剪切方向產生拉剪復合破壞，聲發(fā)射事件也集中于此區(qū)域(圖5a)。

對于單含孔洞試樣而言，聲發(fā)射事件首先在孔口上下集中，為拉破壞；隨后沿剪切方向產生，為拉剪復合破壞(圖5b)。

對于含孔–隙試件來說，聲發(fā)射事件首先在孔口兩側或預制裂紋尖端產生拉破壞，在壓載一段時間后產生剪切破壞，且剪切破壞主要沿著次生裂紋的擴展方向，可見次生裂紋主要是由剪切破壞產生(圖5c—圖5g)；對于傾角較大的情況而言(圖5h—圖5i)，由于預制裂隙的傾角較大，剪切裂紋則從孔洞周圍產生。

聲發(fā)射能量與聲發(fā)射事件隨計算步呈現前期緩慢增大，后期迅速增大的規(guī)律，且預制裂紋的傾角越小，聲發(fā)射的累計能量越大(圖6)。

3.2.2 工況3

不同非均質系數下的聲發(fā)射規(guī)律如圖7—圖8所示。由圖可見，不同非均質系數的聲發(fā)射規(guī)律和累計能量曲線與前述規(guī)律較為類似，值得注意的是，非均質系數越大，聲發(fā)射的累計能量也越大。

3.3 應力應變分析

不同試樣的應力–應變曲線如圖9所示。

由圖9可知，不同裂紋傾角及不同非均質系數試件的應力–應變曲線均經歷3個階段：

a.彈性變形階段 此階段應力隨應變基本上呈線性分布；

b.非線性變形階段 此階段應力–應變曲線偏離線性變形，主要是由于隨著巖體內部隨機分布的相變基元數量的增加和局部裂紋的形成，巖石試件的彈性模量逐漸降低；

圖5 工況1、工況2及工況4聲發(fā)射規(guī)律

圖6 聲發(fā)射累計能量與聲發(fā)射累計數(m=10)

圖7 不同非均質系數下的聲發(fā)射規(guī)律

圖8 聲發(fā)射累計能量與聲發(fā)射累計數(45°)

c.殘余變形階段 此階段應力水平較低，但是變形較大，試件形成宏觀剪切面。

不同非均質系數對試樣的應力–應變曲線影響較小，但對不同孔–隙形式影響較大：

①孔–隙的存在極大地降低了試樣的峰值強度，且預制裂隙傾角越小，試樣的強度降幅越大；

②孔–隙的存在同時降低了試樣整體的脆性度，完整試樣的應力–應變曲線呈典型的脆性破壞模式，即試樣到達峰值荷載時應力呈現“斷崖式”跌落(圖9a)，含孔–隙試樣則具有一定的峰后特征，且預制裂隙傾角越小，延性破壞特征越明顯。

4 與現有研究成果對比

為驗證本文數值模擬的合理性，將模擬結果與本文計算模型的原試驗結果[19]進行對比。文獻[19]選用典型的砂巖作為試驗材料，試樣尺寸為120 mm× 60 mm×30 mm，孔洞直徑=10 mm，孔洞兩側裂隙長度2=2=12 mm，裂隙傾角為45°，采用單軸壓縮加載方式，加載速率為4×10–6m/s，加載至試樣破壞，最終破壞形態(tài)如圖10a所示；圖10b為本文的數值模擬結果。由圖10可知，數值模擬與原試驗結果呈較好的一致性，即產生了翼裂紋與次生裂紋；同時，裂紋產生位置也有較好的對應性，從而驗證了本文數值模擬的合理性。

圖9 不同試樣的應力–應變曲線

5 結論

a. 剪切裂紋產生于完整試樣，裂紋尖端及孔口拉裂紋出現在含孔–隙試樣，裂紋傾角與非均質系數對裂紋擴展影響巨大。

b.拉伸破壞發(fā)生在受載前期，而拉–剪破壞則出現在后期，聲發(fā)射數及能量前期增速較慢，后期迅速加快，預制裂紋傾角與聲發(fā)射累計能量呈負相關，而非均質系數則相反。

c. 受載下的應力–應變曲線均經歷3個階段：彈性變形階段、非線性變形階段及殘余變形階段，孔–隙的存在降低了試樣的峰值強度，影響試件的脆性度。

d. 將數值模擬結果與原試驗結果進行對比，裂紋擴展發(fā)生在預制裂紋尖端與孔口周圍，試驗結果與數值模擬結果一致，驗證了數值模擬結果的可靠性。

圖10 數值模擬與試驗結果的對比

e. 預制裂紋傾角小于30°時，裂紋從預制裂隙的中部產生，大于30°時裂紋從預制裂紋尖端產生；非均質系數越大，裂紋擴展分叉越少。

請聽作者語音介紹創(chuàng)新技術成果等信息，歡迎與作者進行交流

[1] 馮彥軍，康紅普. 受壓脆性巖石Ⅰ–Ⅱ型復合裂紋水力壓裂研究[J]. 煤炭學報，2013，38(2)：226–232.FENG Yanjun，KANG Hongpu. The initiation ofⅠ-Ⅱmixed mode crack subjected to hydraulic pressure in brittle rock under compression[J]. Journal of China Coal Society，2013，38(2)：226–232.

[2] 嚴鵬，陳拓，盧文波，等. 巖爆動力學機理及其控制研究進展[J]. 武漢大學學報(工學版)，2018，51(1)：1–14. YAN Peng，CHEN Tuo，LU Wenbo，et al. A review of dynamic mechanism and controlling of rock burst[J]. Journal of Wuhan University(Engineering Edition)，2018，51(1)：1–14.

[3] 郭建強，劉新榮. 強度準則與巖爆判據統(tǒng)一的研究[J]. 巖石力學與工程學報，2018，37(增刊1)：3340–3352. GUO Jianqiang，LIU Xinrong. Study on the uniformity between strength criterion and rock burst criterion[J]. Journal of Rock Mechanics and Engineering，2018，37(S1)：3340–3352.

[4] 郭品坤. 煤與瓦斯突出層裂發(fā)展機制研究[D]. 徐州：中國礦業(yè)大學，2014. GUO Pinkun. Research on laminar spallation mechanism of coal and gas outburst propagation[D]. Xuzhou：China University of Mining and Technology，2014.

[5] 劉瑞峰，朱哲明，劉邦，等. 爆炸載荷下砂巖動態(tài)斷裂特性的試驗研究[J]. 巖石力學與工程學報，2019，38(3)：445–454.LIU Ruifeng，ZHU Zheming，LIU Bang，et al. Experimental study on dynamic fracture characteristics of sandstones under blasting[J]. Journal of Rock Mechanics and Engineering，2019，38(3)：445–454.

[6] 陳雪峰，趙孝學，汪海波，等. 節(jié)理充填巖體爆炸應力波傳播規(guī)律模型試驗與應用研究[J]. 中國安全生產科學技術，2018，14(12)：130–134. CHEN Xuefeng，ZHAO Xiaoxue，WANG Haibo，et al. Model tests and application research on propagation laws of blasting stress wave in jointed and filled rock mass[J]. Journal of Safety Science and Technology，2018，14(12)：130–134.

[7] 王飛，曹平，曹日紅，等. 平行節(jié)理相互作用對節(jié)理巖體的力學行為影響[J]. 中南大學學報(自然科學版)，2018，49(10)：2498–2507. WANG Fei，CAO Ping，CAO Rihong，et al. Influence of parallel joint interaction on mechanical behavior of jointed rock mass[J]. Journal of Central South University(Science and Technology)，2018，49(10)：2498–2507.

[8] 岳中文，邱鵬，楊仁樹，等. 切槽炮孔偏心裝藥爆源近區(qū)裂紋動態(tài)力學特征實驗研究[J]. 巖石力學與工程學報，2017，36(增刊2)：3792–3798.YUE Zhongwen，QIU Peng，YANG Renshu，et al. Experimental study on crack dynamic mechanical characteristic near blasting source of grooved borehole with eccentric decouple charge[J]. Journal of Rock Mechanics and Engineering，2017，36(S2)：3792–3798.

[9] 李竟艷，高文學，宋肖龍. 巖石材料動載下泛形裂紋擴展數值模擬[J]. 振動與沖擊，2018，37(22)：88–91. LI Jingyan，GAO Wenxue，SONG Xiaolong. Numerical simulation on the extension of a ubiquitiformal crack of rock materials under dynamic loading[J]. Journal of Vibration and shock，2018，37(22)：88–91.

[10] 李宏，張麗萍，鄧清海，等. 非貫通裂隙巖體的單軸壓縮顆粒流數值模擬[J]. 水力發(fā)電，2018，44(11)：26–30.LI Hong，ZHANG Liping，DENG Qinghai，et al. Numerical simulation of uniaxial compression particle flow in non-penetrating fractured rock mass[J]. Water Power，2018，44(11)：26–30.

[11] 周喻，劉冰，王莉，等. 單軸壓縮條件下含雙圓孔類巖石試樣力學特性的細觀研究[J]. 巖石力學與工程學報，2017，36(11)：2662–2671. ZHOU Yu，LIU Bing，WANG Li，et al. Mesoscopic mechanical properties of rock-like material containing two circular holes under uniaxial compression[J]. Journal of Rock Mechanics and Engineering，2017，36(11)：2662–2671.

[12] 章廣成，胡靜. 非貫通單節(jié)理巖體裂紋擴展方向研究[J]. 煤田地質與勘探，2011，39(4)：43–48. ZHANG Guangcheng，HU Jing. Study on crack propagation direction of intermittent single jointed rock mass[J]. Coal Geology & Exploration，2011，39(4)：43–48.

[13] 張力民，張慧，劉紅巖. 單軸壓縮荷載下非貫通閉合節(jié)理巖體損傷本構模型[J]. 煤田地質與勘探，2016，44(1)：79–84. ZHANG Limin，ZHANG Hui，LIU Hongyan. A damage constitutive model for rock mass with non-persistently closed joints under uniaxial compression load[J]. Coal Geology & Exploration，2016，44(1)：79–84.

[14] 姚改煥，宋戰(zhàn)平，余賢斌. 石灰?guī)r聲發(fā)射特性的試驗研究[J].煤田地質與勘探，2006，34(6)：44–46. YAO Gaihuan，SONG Zhanping，YU Xianbin. Experimental study on acoustic emission characteristics of limestone[J]. Coal Geology & Exploration，2006，34(6)：44–46.

[15] 駱中山，曹函，孫平賀，等. 水力壓裂脈沖頻率對煤巖預制孔密封效果影響[J]. 煤田地質與勘探，2018，46(2)：190–195. LUO Zhongshan，CAO Han，SUN Pinghe，et al. Influence of hydraulic fracturing pulse frequency on precasting sealing effect of coal rock[J]. Coal Geology & Exploration，2018，46(2)：190–195.

[16] 張帆，馬耕，劉曉，等. 煤巖水力壓裂起裂壓力和裂縫擴展機制實驗研究[J]. 煤田地質與勘探，2017，45(6)：84–89. ZHANG Fan，MA Geng，LIU Xiao，et al. Experimental study on initiation pressure and mechanism of fracture propagation of hydraulic fracturing in coal and rock mass[J]. Coal Geology & Exploration，2017，45(6)：84–89.

[17] 李地元，成騰蛟，周韜，等. 沖擊載荷作用下含孔洞大理巖動態(tài)力學破壞特性試驗研究[J]. 巖石力學與工程學報，2015，34(2)：249–260. LI Diyuan，CHENG Tengjiao，ZHOU Tao，et al. Experimental study of the dynamic strength and fracturing characteristics of marble specimens with a single hole under impact loading[J].Journal of Rock Mechanics and Engineering，2015，34(2)：249–260.

[18] 王蒙，朱哲明，謝軍. 巖石Ⅰ–Ⅱ復合型裂紋動態(tài)擴展SHPB實驗及數值模擬研究[J]. 巖石力學與工程學報，2015，34(12)：2474–2485. WANG Meng，ZHU Zheming，XIE Jun. Experimental and numerical studies of the mixed-modeⅠandⅡcrack propagation under dynamic loading using SHPB[J]. Journal of Rock Mechanics and Engineering，2015，34(12)：2474–2485.

[19] 程龍，楊圣奇，劉相如. 含缺陷砂巖裂紋擴展特征試驗與模擬研究[J]. 采礦與安全工程學報，2012，29(5)：719–724. CHENG Long，YANG Shengqi，LIU Xiangru. Experimental and numerical investigation on crack expansion of sandstone containing flaws[J]. Journal of Mining & Safety Engineering，2012，29(5)：719–724.

Numerical simulation of propagation law of pore-fracture under uniaxial stress

GUO Xiucheng

(Engineering and Technology Research Institute of Great Wall Drilling Engineering Co. Ltd., Panjin 124010, China)

In order to study the law of crack propagation under different pore-fracture conditions, the crack propagation process, acoustic emission law and stress-strain curve under different crack inclination angle and different heterogeneous coefficient were numerically simulated by using RFPA software. At the same time, the results are compared with the experimental results. The results show that the cracks of the intact specimen are produced along the shear direction, while the cracks of the specimens with pores occur along the side of the crack tip and the orifice. The secondary cracks in horizontal direction and bifurcation at the prefabricated crack tip or the side of the orifice at the same time are produced through the specimens with wing cracks, and the heterogeneity coefficient affects the strike of the secondary cracks. Pre-ballast specimen with tensile failure in the late stage of the main and ballast stage, the tensile-shear combination is mainly destroyed, and the generation of secondary cracks is related to the shear failure. The cumulative energy of the acoustic emission and the acoustic emission increased slowly in the early stage; and rapidly in the later stage, the smaller the inclination angle of the pre-cast crack, the larger the heterogeneous coefficient, the greater the accumulated energy of the acoustic emission, and the stress time curves of the test pieces at different crack inclination angles and different mean values are all subjected to three stages. The existence of the porosity reduces the peak strength of the samples and affects the brittleness of the specimens. The results of the study provide some reference for further understanding the law of the pore-gap interaction.

uniaxial stress; pore-fracture action; acoustic emission; crack propagation; stress-strain curve; numerical simulation

TD31

A

10.3969/j.issn.1001-1986.2020.02.027

1001-1986(2020)02-0179-08

2019-05-20；

2019-11-18

長城鉆探工程有限公司科研項目(GWDC2017B-01)

Science and Research Project of Engineering and Technology Research Institute of Great Wall Drilling Engineering Co. Ltd.(GWDC2017B-01)

郭修成，1986年生，男，四川廣安人，工程師，從事鉆井工程等相關研究工作. E-mail：guoxc.gwdc@cnpc.com.cn

郭修成. 單軸應力下的孔洞–裂隙擴展規(guī)律數值模擬[J]. 煤田地質與勘探，2020，48(2)：179–186.

GUO Xiucheng. Numerical simulation of propagation law of pore-fracture under uniaxial stress[J]. Coal Geology & Exploration，2020，48(2)：179–186.

(責任編輯 周建軍)