張?zhí)燔? 王喜娜 景晨 張磊 潘紅宇 張碩

摘 要：為探究含裂隙軟煤的局部化演化特征，將石膏與水以質量比7：3混合制成類軟煤含裂隙試樣，進行單軸壓縮試驗并采用數(shù)字散斑相關測量方法采集加載全程試樣表面應變、位移信息，分析其變形位移和變形能密度演化特征。結果表明：不同傾角試樣局部化帶開始出現(xiàn)位置和擴展方向不同，傾角為0°時局部化帶沿預制裂隙中央向試樣上端和下端部擴展，60°時局部化帶沿預制裂隙端部向試樣對角方向擴展，隨傾角增大，局部化帶穩(wěn)定時所對應的應力水平接近1。試樣破壞由局部化帶的張開和錯動共同作用導致，局部化帶張開與錯動速率存在緩慢演化和加速演化階段，應力峰值前后，局部化帶之間擴展交匯引起速率加速演化。局部化帶的能量演化經歷累積（A）、累積與釋放交替（A-R）、釋放（R）3個階段。在加載初始階段，總體表現(xiàn)出A過程，加載至應力峰值前后，能量發(fā)生A-R波動式演化過程，加載至應力峰值后期，表現(xiàn)為R過程，A-R階段隨裂隙傾角增大發(fā)生相對后移現(xiàn)象。關鍵詞：數(shù)字散斑;局部化帶;錯動速率;變形能密度中圖分類號：TD 315

文獻標志碼：A

文章編號：1672-9315（2021）02-0203-10

DOI：10.13800/j.cnki.xakjdxxb.2021.0203開放科學（資源服務）標識碼（OSID）：

Localized deformation characteristicsof soft coal

with fractures at different dip angles

ZHANG Tianjun1，WANG Xina1，JING Chen1，ZHANG Lei2，PAN Hongyu1，ZHANG Shuo3

（1.College of Safety Science and Engineering，Xian University of Science and Technology，Xian 710054，China;

2.College of Energy Science and Engineering，Xian University of Science and Technology，Xian 710054，China;

3.College ofSciences，Xian University of Science and Technology，Xian 710054，China）

Abstract：In order to explore the localized evolution characteristics of soft coal with fissures，uniaxial compression tests on soft coal samples with single fissures at different inclination angles were carried out.The digital speckle correlation measurement method was used to collect the surface strain and displacement information of the sample during the entire loading process，with its localized displacement and evolution characteristics of deformation energy density analyzed.The results show that：The localization zone of the sample with different dip angles began to appear different positions and expansion directions.When the dip Angle was 0，the localization zone expanded step by step along the center of the prefabricated crack to the upper and lower ends of the sample.When the dip Angle was 60，the localization zone expanded along the end of the prefabricated crack to the diagonal direction of the sample.The failure of the specimen is caused by the joint action of the opening and dislocation of the localized zone，and there is a slow evolution of the localized zone opening and dislocation rate.And the accelerated evolution stage，before and after the stress peak，the expansion and intersection between the localized zones cause the rate to accelerate.The energy evolution of the localized zone of the soft coal samples with different fracture angles undergoes accumulation（A），alternating accumulation and release（A-R），release（R）Three stages.In the initial stage of loading，the overall A process is shown.Before and after loading to the peak stress，the energy undergoes an A-R fluctuating evolution process.When the load reaches the late stress peak，it appears as an R process.The A-R phase occurs relative to the increase of the crack inclination angle showing backward movement phenomenon.Key words：digital speckle;localized zone;dislocation rate;deformation energy density

0 引 言煤巖體中微裂紋結構和變形局部化行為會造成其失穩(wěn)破壞，從而導致嚴重的礦井災害事故[1]。目前針對堅硬煤巖變形局部化演化已取得大量成果。而軟煤作為一種特殊煤巖，其力學特性和變形特征與堅硬煤巖有較大區(qū)別[2]，因此從變形局部化角度研究軟煤失穩(wěn)破壞過程具有重要意義。目前國內外學者在煤巖的變形局部化方面做了大量研究。SHIOYA等采用立體攝影測量技術對單軸壓縮下煤巖體的變形局部化進行研究[3]。趙豫生等應用激光全息干涉法觀測單軸壓縮條件下巖石局部化變形現(xiàn)象[4]。為實現(xiàn)應變的全局觀測，孫慶平等綜合運用高密度光柵、形貌以及激光Raman微探針技術對陶瓷材料進行單軸以及三點彎曲下塑性變形局部化行為研究[5]。為量化描述材料局部化演化特征，潘一山、馬少鵬等采用白光散斑相關方法對巖土材料變形局部化進行研究，測定了巖石變形局部化開始時刻、寬度、傾角以及演化規(guī)律[6-7]。上述成果主要是針對完整試樣在受荷變形過程中的局部化演化規(guī)律進行分析，且研究材料集中在巖體和硬煤領域。而軟煤在煤體強度、變形特征和破裂特征等方面與中硬煤有顯著差異[8]。目前針對軟煤變形局部化領域研究成果較少，其中毛靈濤等基于數(shù)字散斑相關方法得到煤樣變形局部化區(qū)域與煤樣最終破壞區(qū)域一致的結果[9]。呂玉凱等進一步研究發(fā)現(xiàn)單軸壓縮下煤樣的變形局部化帶先后經歷了相對滑動主導、相對拉伸主導以及相對拉伸和滑動共同主導3個階段，但未深入分析其主導因素[10]。王學濱等探究了單軸壓縮下煤樣最大剪切應變變異系數(shù)與變形局部化帶的關系，得到局部化帶變清晰的過程中會引起變異系數(shù)的突增[11]。張?zhí)燔姷葟牟煌式嵌忍骄苛怂绍浢簶釉趩屋S壓縮下的表面變形局部化演化特征，認為裂隙擴展速率隨含水率增高而降低[12]。這些研究成果對于了解軟煤試樣變形局部化具有積極意義。但由于實驗條件限制，大多數(shù)研究局限于定性分析，且研究材料未考慮含裂隙軟煤。因此對含裂隙軟煤試樣的局部化演化規(guī)律和破壞特征方面有待進一步研究。文中將石膏與水以質量比為7∶3混合制成含裂隙軟煤試樣，進行單軸壓縮試驗并采用數(shù)字散斑相關測量方法（digital speckle correlation measurement method，DSCM）獲取試樣加載過程中表面應變與位移信息。分析不同傾角試樣表面變形局部化演化特征及對應裂紋起裂、擴展狀態(tài)，分階段探究不同傾角試樣局部化帶位移和能量演化規(guī)律。

1 裂隙軟煤試樣單軸壓縮試驗

1.1 試樣制備軟煤取樣困難且不易加工成型，可采用相似材料對其研究[13]。石膏材料與煤體的力學性質相似，且價格低廉。表1中試樣1～3及S1～S3分別是山西和新疆某礦松軟煤樣，文中所制作的試樣各參數(shù)均在松軟煤的參數(shù)范圍內，可進行松軟煤體的相關研究。

首先，將石膏與水以7∶3質量比混合配置成試樣漿液。接著，將所得漿液澆注于70 mm×70 mm×140 mm方形試樣盒中。采用預埋抽條法在試樣中心部位預制單條裂隙。大量實驗表明0°與60°單裂隙對試樣影響較大，因此預制φ=0°及φ=60°這2個典型裂隙傾角，裂隙寬度d=0.2 mm，長度2L=10 mm，深度c=35 mm。最后使用黑白啞光噴漆對干燥試樣含裂隙面進行處理。所制試樣情況見表2。

1.2 試驗系統(tǒng)試驗加載系統(tǒng)采用DNS200電子萬能試驗機，加載方式為位移控制，加載速率為0.5 mm/min，數(shù)據(jù)采集頻率為1 Hz。為降低端部效應影響，在試樣與上下壓頭之間涂抹凡士林。試驗采用DSCM系統(tǒng)捕捉全場變形圖像信息[15]。該測量系統(tǒng)包括光源、CCD攝像機、圖像采集系統(tǒng)和VIC-3D軟件。試驗時將測量物體放置在加載設備中，由圖像采集系統(tǒng)控制采樣頻率和分辨率，由CCD攝像機采集圖像，通過VIC-3D軟件計算試樣加載過程中全場應變、位移信息。為防止頻閃對數(shù)據(jù)采集的影響，用兩架OSRAM55460光學儀器燈對物體表面進行照射[16]。試驗中圖像采集速率為1.0幀/s，圖像分辨率為1 600 pixel×1 200 pixel，物面分辨率為0.1 316 mm/pixel。系統(tǒng)構成如圖1所示。試驗要求加載系統(tǒng)與圖像采集系統(tǒng)時間同步，開始加載的同時采用VIC-3D軟件的連續(xù)采集功能采集試樣表面變形圖像，直至試樣發(fā)生破壞。

1.3 試驗結果通過試驗得到D0與D60組試樣應力-應變曲線，如圖2（a）和（b）所示。同組裂隙傾角試樣均一性較好。D0組試樣在0.007<ε<0.009內達到峰值強度，其峰值強度范圍在3.09～3.21 MPa之間，平均值3.16 MPa，最大相差0.12 MPa，彈性模量在230.00～271.10 MPa之間，平均值247.87 MPa，最大相差41.1 MPa。D60組試樣在0.022<ε<0.028內達到峰值強度，其峰值強度范圍在7.60～8.70 MPa之間，平均值8.10 MPa，最大相差1.08 MPa，彈性模量在742.00～786.63 MPa之間，平均值766.75 MPa，最大相差44.63 MPa。

不同裂隙傾角試樣應力-應變曲線存在差異，D60組試樣彈性模量和峰值強度均大于D0組試樣，60°試樣峰值強度是0°試樣的2.3～2.8倍，彈性模量是0°試樣的3.1～3.7倍。這與汪杰等試驗得到的結果吻合度較高[17]。由于預制裂隙的存在會降低試樣的強度，預制傾角為0°時，加載方向與預制裂隙面垂直，試樣強度受預制裂隙的影響很小;預制傾角為60°時，試樣強度主要受裂隙面控制，預制裂隙的作用較明顯。當受到加載后，裂隙首先產生閉合，由于裂隙傾角越大，裂隙面之間產生摩擦阻力越大，提高了試樣強度，因此D60組試樣強度高于D0組試樣。D60組試樣彈性模量高于D0組試樣是由于0°預制裂隙的壓縮變形在試樣總變形中所占比重較大而引起的。

2 試驗結果分析

2.1 變形局部化區(qū)域劃分采用主應變差值法[18]來確定局部化區(qū)域，即劃分區(qū)的主應變差值較周圍區(qū)域呈現(xiàn)明顯局部性集中。如圖5（a）和（b）所示，根據(jù)試樣最終破壞模式及破壞前變形場，將變形局部化帶所在位置（也是最終宏觀破壞出現(xiàn)位置）按出現(xiàn)先后順序用Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ標注。為對比分析不同裂隙傾角軟煤試樣變形場演化特征，將D0-2與D60-3試樣加載全程應力按照公式（1）進行歸一化處理。同時，以峰值應力對應的應變εc為基準，采用（2）式對應變進行處理。圖3為應力-應變水平曲線。

σstd=σσc

（1）

式中 σstd為歸一化的應力水平（0<σstd<1）;σ為瞬時應力，MPa;σc為峰值應力，MPa。

εstd=εεc

（2）

式中 εstd為應變水平;ε為瞬時應變;εc為峰值應力對應的應變。取同一應力水平下標識點所對應時刻應變云圖進行變形局部化分析。如圖3所示，以應變值出現(xiàn)突變?yōu)橐罁?jù)選取標識點1～4，分別對應應力水平0.2σc、0.6σc、0.79σc、0.92σc，標識點5位于應力峰值后破壞階段。

2.2 裂隙軟煤試樣變形場演化分析為分析不同裂隙傾角軟煤試樣加載全程變形場演化特征，采用DSCM系統(tǒng)計算得到加載全程應變場云圖。圖6為試樣峰值后（標試點5）應變云圖，其中拉應變?yōu)檎?，壓應變?yōu)樨?。不同顏色代表不同應變區(qū)域，紅色代表局部高拉應變區(qū)，紫色代表局部高壓應變區(qū)[19]。選取圖6應變云圖中50 mm×80 mm區(qū)域進行放大，取圖3同一應力水平下標識點1～5時刻所對應云圖為參考，從細觀層面分析試樣應變演化過程。D0-3試樣應變演化如圖4（a）～（d）所示。在0.2σc應力水平下，如圖4（a）所示，試樣處于壓密階段，應變云圖顏色基本一致，說明試樣表面變形均勻。在0.6σc應力水平下，如圖4（b）所示，試樣處于彈性階段，試樣應變值較上一應力水平有微小增大，預制裂隙中央位置開始出現(xiàn)變形局部化現(xiàn)象。根據(jù)朱泉企等[20]的大量實驗，該應力水平前試樣表面沒有裂紋萌生。繼續(xù)加載至0.79σc應力水平下，如圖4（c）所示，此時試樣變形集中區(qū)域范圍和變形量值都有所增加，局部化帶Ⅰ基本穩(wěn)定成形。當應力水平達到0.92σc，如圖4（d）所示，變形量最大值達到0.003 46，變形進一步集中。隨著加載繼續(xù)進行，局部化帶Ⅰ向上端延伸，試樣最終發(fā)生整體破壞。D60-2試樣應變演化如圖5（e）～（h）所示。試樣前3個應力水平下，如圖5（a）、（b）、（c）所示。與D0-3試樣區(qū)別在于其變形局部化并不明顯。這是由于D60-2試樣結構主要受裂隙面控制，加載前期結構較為穩(wěn)定，未有明顯應力集中。當應力水平達到0.92σc時，如圖5（d）所示，在預制裂隙端部出現(xiàn)2個明顯的變形集中區(qū)域，并有交匯趨勢，應變量值達到0.024 4。隨著加載的進行，分別向右上角和左側擴展貫通，最終導致試樣破壞。

試樣最終破壞形態(tài)和采用數(shù)字散斑相關方法分析得到試樣峰后變形場云圖對比結果如圖6所示。從圖6可以看出計算得出的變形集中區(qū)域與試樣最終破壞區(qū)域相吻合。D0-3試樣最終破壞形態(tài)如圖6（a）所示。變形集中區(qū)域由預制裂隙中央開始向上下兩端孕育演化成局部化帶，最后形成一條貫穿預制裂隙的宏觀裂紋，并與沿左側以及預制裂隙端部的微裂紋貫通，導致試樣最終破壞。D60-2試樣最終破壞形態(tài)如圖6（b）所示。與D0-3試樣不同之處在于其變形局部化帶由預制裂隙兩端分別向試樣對角方向孕育演化，最終與左側形成的微裂紋貫通導致試樣破壞。

2.3 裂隙軟煤局部化帶位移演化分析為量化不同裂隙傾角試樣變形局部化特征，通過局部化帶位移演化速率進行分析。將位移速率分為張開位移速率和錯動位移速率，其計算過程為公式（3）。

Vo=Cot，

Vd=Cdt

（3）

式中 Co為裂紋張開位移，mm;Cd為裂紋錯動位移，mm;Vo為裂紋張開速率，mm/s;Vd為裂紋錯速率，mm/s;t為時間，s。

而上式中的局部化帶張開位移與錯動位移，采用如圖7所示方法計算。在變形局部化帶兩側分別對稱的選取5組測點（以測點1、測點2為中心對稱區(qū)域），位移分解為x方向（水平）位移u和y方向（豎直）位移v。相對張開與錯動位移計算公式為[21]

式中 u′為兩側點x方向位移差，mm;v′為兩測點y方向位移差，mm;θ為兩側點連線與x軸夾角。按照上述公式得到加載全程D0-3和D60-2試樣各條局部化帶張開與錯動位移速率隨時間演化曲線，如圖7（a）～（d）所示。

由圖可以看出，不同裂隙傾角試樣局部化帶位移速率演化過程相似，均存在緩慢演化和加速演化階段。加載初期，局部化帶張開與錯動均不明顯。加載至應力峰值附近時，D0-3沿局部化帶Ⅰ發(fā)生張開，Ⅱ、Ⅲ發(fā)生錯動，D60-2試樣沿局部化帶Ⅰ、Ⅱ發(fā)生張開，Ⅲ發(fā)生錯動。應力峰值后，由于局部化帶之間擴展、交匯、貫通導致D0-3與D60-2試樣張開與錯動速率加速演化。對比圖8（a）和（b）發(fā)現(xiàn)D0-3試樣錯動速率量值高于錯動速率值一個量級，對比圖8（c）和（d）發(fā)現(xiàn)D60-2試樣張開與錯動速率量級相當。這表明D0-3試樣應力峰值前以張拉破壞為主，峰值后以張拉剪切混合破壞為主，D60-2試樣整體以張拉剪切混合破壞為主。由2.2節(jié)應變場演化可知加載至應力峰值前后，試樣發(fā)生宏觀裂紋擴展，宏觀裂紋擴展導致裂隙試樣剪切錯動速率演化明顯，在此對試樣錯動速率進一步分析。如圖8（b）所示，分析D0-3試樣錯動速率演化曲線。從圖可以得到加載時間為52 s（0.92σc）時，局部化帶Ⅰ與Ⅱ錯動速率值為9.17×10-4mm/s，19.4×10-4mm/s。加載時間為55 s（σc）時，局部化帶Ⅲ發(fā)生剪切錯動引起錯動速率迅速增大至25.6×10-4 mm/s;加載時間為59 s（峰后）時，局部化帶Ⅳ與Ⅴ迅速增加至39.9×10-4 mm/s，36.7×10-4 mm/s，此時試樣發(fā)生失穩(wěn)破壞。D0-3試樣局部化帶錯動速率在峰值前后迅速增加，主要原因是試樣應變場在應力峰值前后劇烈演化，各局部化帶迅速貫通交匯導致。

圖8（d）為D60-2試樣錯動速率演化曲線。從圖可以得到加載時間為203 s（0.92σc）時，局部化帶Ⅱ錯動速率值為2.7×10-4mm/s;加載時間為225 s時，局部化帶Ⅰ錯動速率達到最大值4.7×10-4 mm/s，加載時間為233（0.98σc）時，局部化帶Ⅳ發(fā)生剪切錯動導致錯動速率迅速增大至13.5×10-4 mm/s，加載時間為261 s與271 s（峰值后）時，局部化帶Ⅲ和Ⅴ相繼達到最大值3.9×10-4 mm/s，9.5×10-4 mm/s，此時試樣發(fā)生失穩(wěn)破壞。D60-2試樣局部化帶錯動速率在峰后迅速增加，主要原因是D60-2試樣應變場在應力峰值后劇烈演化，各局部化帶迅速貫通交匯導致。對比（b）圖與（d）圖，發(fā)現(xiàn)D60-2試樣錯動速率量值明顯低于D0-3試樣。其原因在于，在相同加載條件下，同種材料的裂隙擴展力相同，但裂隙傾角增大，擴展阻力隨之增大[22]，導致其擴展程度降低，表現(xiàn)為D60-2試樣錯動速率量值上低于D0-3試樣一個量級。

2.4 裂隙軟煤試樣變形能量演化分析學者們認為是局部化帶能量驅動煤巖體失穩(wěn)破壞[23-25]。為研究軟煤試樣破裂釋放能量的條件，采用變形能密度研究試樣破壞變形過程能量演化特征，對比不同裂隙傾角對能量演化的影響。由于工程中大多情況下是圍巖參與釋放彈性能，因此能量分析區(qū)域選為圖5（a）和（b）所示變形局部化帶外區(qū)域，用分析區(qū)域各個點應變分量平均

值表示此區(qū)域應變分量值，計算得到軟煤試樣變形局部化帶區(qū)域變形能密度量值。其計算公式如下[22]

U=E2（ε21+ε22-2με1ε2）

（5）

式中 E為試件彈性模量，MPa;μ為試樣泊松比;ε1為第一主應變;ε2為第二主應變;U為變形能密度，J/m3。

圖9（a）和（b）分別為D0-3與D60-2試樣局部化帶變形能密度演化曲線。可以看出試樣演化經歷3各階段，分別能量累積的A階段、能量累積與釋放交替的A-R階段、能量釋放的R階段。在應力加載峰值前，由于預制裂隙試樣處于壓密階段，變形局部化帶以能量積累為主。應力峰值前后，變形局部化帶既有能量積累又有釋放。加載至后期能量表現(xiàn)為釋放。由于兩試樣能量演化規(guī)律具有相似性，故在此僅以D0-3試樣演化過程做詳細分析，其余不予贅述。D0-3試樣加載至峰值前階段，局部化帶Ⅰ與Ⅱ幾乎同時達到能量峰值后開始釋放能量，這是由于局部化帶Ⅰ與Ⅱ首先發(fā)生微小變形而產生影響。局部化帶Ⅲ、Ⅳ、Ⅴ還在進行能量積累。在加載應力峰值后階段，局部化帶Ⅳ與Ⅴ快速增加，并達到能量峰值2.02×106 J/m3，2.23×106 J/m3，同時，局部化帶Ⅰ與Ⅱ能量下降到9.6×105 J/m3，64×105 J/m3，局部化帶Ⅲ能量積累達到峰值325×105 J/m3，然后迅速釋放能量，進一步擴展。其中局部化帶Ⅲ能量積累速率最慢，累積能量最高，主要原因是局部化帶Ⅱ向下擴展，以及局部化帶Ⅳ向上擴展均與其連接貫通而對其產生影響，這圖4云圖結果一致。從圖中可以看出，各局部化帶錯動速率突增時刻與能量密度累積與釋放交替時刻相對應，表明裂隙面錯動狀態(tài)變化是引起變形能密度出現(xiàn)波動的主要原因。

結合試樣變形場、位移場演化特征，含裂隙軟煤試樣能量演化主要體現(xiàn)出下列2個特征：一是含裂隙軟煤試樣能量演化經歷A、A-R、R這3個階段。二是在A-R階段，變形能開始發(fā)生累積和釋放的波動式變化，局部化帶錯動速率出現(xiàn)最值時刻與能量演化交替時刻相對應，這與試樣破壞和局部化帶錯動速率密切相關的結論是一致的[26]。對比圖9（a）和（b）可看出，D60-2試樣A-R階段相對于D0-3試樣有一定程度后移。產生這種現(xiàn)象的原因是裂隙擴展能量釋放率（即為裂隙擴展力）大于消耗率（即為裂隙擴展阻力）能量，表現(xiàn)為裂隙將擴展[22]。而D60-2試樣擴展阻力較大，相比而言D0-3試樣裂紋更易擴展從而引起能量較早釋放。這與唐禮忠數(shù)值計算出隨裂隙傾角增大，試樣發(fā)生能量釋放相對后移的結果一致[27]。

3 結 論

1）含不同傾角試樣局部化帶出現(xiàn)位置和擴展方向不同。傾角為0°時局部化帶沿預制裂隙中央向試樣上端和下端部擴展，為60°時局部化帶沿預制裂隙端部向試樣對角方向擴展。傾角為0°時試樣在0.79σc應力水平下開始出現(xiàn)穩(wěn)定的變形局部化帶，傾角為60°時在0.92σc應力水平下開始出現(xiàn)穩(wěn)定的應變局部化帶。

2）試樣破壞是由局部化帶的張開和錯動共同作用導致的，其張開與錯動速率存在緩慢演化和加速演化階段。應力峰值前后，變形局部化帶之間擴展、交匯、貫通引起局部化帶張開與錯動速率加速演化，該階段傾角為0°試樣錯動速率值高于60°試樣一個量級。

3）含不同裂隙傾角試樣能量演化經歷累積（A）、累積與釋放交替（A-R）、釋放（R）3個階段。在應力加載初始階段，總體表現(xiàn)出A過程，加載至應力峰值前后能量發(fā)生A-R的波動式演化過程，加載至應力峰值后期，表現(xiàn)為R過程。A-R階段隨裂隙傾角增大產生相對后移現(xiàn)象。

參考文獻（References）：

[1] 呂玉凱，姜耀東.不同沖擊傾向性煤樣變形場演化特征[J].煤炭學報，2014，39（11）：2172-2176.LYU Yukai，JIANG Yaodong.Deformation field evolution characteristics of coal samples with different shock tendencies[J].Journal of China Coal Society，2014，39（11）：2172-2176.

[2]伍永平，高喜才.不同加載模式軟硬煤巖側向變形特征的對比實驗研究[J].煤炭學報，2010，35（S1）：44-48.WU Yongping，GAO Xicai.Comparative experimental study on lateral deformation characteristics of soft and hard coal with different loading modes[J].Journal of China Coal Society，2010，35（S1）：44-48.

[3]SHIOYA T.Analysis of crack propagation in brittle material[J].JSME International Journal.Ser.1，Solid Mechanics，Strength of Materials，1994，37（1）：1-6.

[4]趙豫生，陳顒，王仁.巖石裂紋附近的變形局部化與破壞——激光全息干涉法的應用[J].地球物理學報，1988，31（3）：294-305.ZHAO Yusheng，CHEN Yong，WANG Ren.Deformation localization and failure near rock cracks：Application of laser holographic interferometry[J].Journal of Geophysics，1988，31（3）：294-305.

[5]孫慶平，趙智軍，卿新林，等.ZrO2相變多晶體塑性變形局部化行為的宏觀—細觀實驗研究[J].中國科學（A輯·數(shù)學·物理學·天文學·技術科學），1994，33（4）：383-388.SUN Qingping，ZHAO Zhijun，QING Xinlin，et al.Macroscopic and microscopic experimental study on the localization behavior of ZrO2 phase transition polycrystal plastic deformation[J].Science of China（Series A，Mathematics，Physics，Astronomy and Technical Science），1994，33（4）：383-388.

[6]馬少鵬，王來貴，趙永紅.巖石圓孔結構破壞過程變形場演化的實驗研究[J].巖土力學，2006，27（7）：1082-1086.MA Shaopeng，WANG Laigui，ZHAO Yonghong.Experimental study on deformation field evolution during failure process of rock circular pore structure[J].Rock and Soil Mechanics，2006，27（7）：1082-1086.

[7]潘一山，楊小彬.巖石變形破壞局部化的白光數(shù)字散斑相關方法研究[J].實驗力學，2001，27（2）：220-225.PAN Yishan，YANG Xiaobin.Study on the correlation method of white light digital speckle for rock deformation and failure localization[J].Experimental Mechanics，2001，27（2）：220-225.

[8]高保彬，李回貴，李林，等.同組軟硬煤煤樣聲發(fā)射及分形特征研究[J].巖石力學與工程學報，2014，33（S2）：3498-3504.

GAO Baobin，LI Huigui，LI Lin，et al.Research on acoustic emission and fractal characteristics of soft and hard coal samples in the same group[J].Journal of Rock Mechanics and Engineering，2014，33（S2）：3498-3504.

[9]毛靈濤，連秀云，郝耐，等.基于數(shù)字體散斑法煤樣內部三維應變場的測量[J].煤炭學報，2015，40（1）：65-72.MAO Lingtao，LIAN Xiuyun，HAO Nai，et al.Measurement of the internal threedimensional strain field of coal samples based on digital volume speckle method[J].Journal of China Coal Society，2015，40（1）：65-72.

[10]呂玉凱，蔣聰，成果，等.不同沖擊傾向煤樣表面溫度場與變形場演化特征[J].煤炭學報，2014，39（2）：273-279.LYU Yukai，JIANG Cong，CHENG Guo，et al.Evolution characteristics of surface temperature field and deformation field of coal samples with different shock tendency[J].Journal of China Coal Society，2014，39（2）：273-279.

[11]王學濱，侯文騰，潘一山，等.基于數(shù)字圖像相關方法的單軸壓縮煤樣應變局部化過程試驗[J].煤炭學報，2018，43（4）：984-992.WANG Xuebin，HOU Wenteng，PAN Yishan，et al.Experiment on the strain localization process of uniaxial compression coal sample based on digital image correlation method[J].Journal of China Coal Society，2018，43（4）：984-992.

[12]張?zhí)燔姡o翔，張磊，等.瓦斯抽采鉆孔孔周裂隙演化及等效裂紋寬度試驗研究[J].巖石力學與工程學報，2019，38（S2）：3625-3633.ZHANG Tianjun，JI Xiang，ZHANG Lei，et al.Experimental study on crack evolution and equivalent crack width around gas drainage boreholes[J].Chinese Journal of Rock Mechanics and Engineering，2019，38（S2）：3625-3633.

[13]張?zhí)燔姡瑥埨?，李樹剛，?瓦斯抽采鉆孔孔周裂紋擴展規(guī)律[J].遼寧工程技術大學學報（自然科學版），2018，37（3）：499-507.ZHANG Tianjun，ZHANG Lei，LI Shugang，et al.Crack propagation law around borehole in gas drainage[J].Journal of Liaoning Technical University（Natural Science Edition），2018，37（3）：499-507.

[14]李勝利，靳鐘銘，魏錦平.軟煤大煤樣的壓裂特性研究[J].太原理工大學學報，1999，30（6）：567-570.LI Shengli，JIN Zhongming，WEI Jinping.Study of crack characteristic of big soft coal samples[J].Journal of Taiyuan University of Technology，1999，30（6）：567-570.

[15]伍永平，高喜才.不同加載模式軟硬煤巖側向變形特征的對比實驗研究[J].煤炭學報，2010，35（S1）：44-48.WU Yongping，GAO Xicai.Experimental comparative study on lateral deformation characteristics of coal sample in different loading path[J].Journal of China Coal Society，2010，35（S1）：44-48.

[15]GAO G，YAO W，XIA K，et al.Investigation of the rate dependence of fracture propagation in rocks using digital image correlation（DIC）method[J].Engineering Fracture Mechanics，2015，138：146-155.

[16]NIRMAL J J S L.Tampering detection using hybrid local and global features in wavelettransformed space with digital images[J].Soft Computing：A Fusion of Foundations，Methodologies and Applications，2020，24（1-3）：5427-5443.

[17]汪杰，李楊，宋衛(wèi)東，等.不同傾角節(jié)理巖體損傷演化特征分析[J].哈爾濱工業(yè)大學學報，2019，51（8）：143-150.WANG Jie，LI Yang，SONG Weidong，et al.Analysis of damage evolution characteristics of jointed rock masses with different dip angles[J].Journal of Harbin Institute of Technology，2019，51（8）：143-150.

[18]鄭捷，姚孝新，陳顒.巖石變形局部化的試驗研究[J].地球物理學報，1982，26（6）：554-563.ZHENG Jie，YAO Xiaoxin，CHEN Yong.Experimental study on rock deformation localization[J].Chinese Journal of Geophysics，1982，26（6）：554-563.

[19]潘紅宇，葛迪，張?zhí)燔?，?應變率對巖石裂隙擴展規(guī)律的影響[J].煤炭學報，2018，43（3）：675-683.PAN Hongyu，Ge Di，ZHANG Tianjun，et al.The influence of strain rate on the law of rock fracture propagation[J].Journal of China Coal Society，2018，43（3）：675-683.

[20]朱泉企，李地元，李夕兵.含預制橢圓形孔洞大理巖變形破壞力學特性試驗研究[J].巖石力學與工程學報，2019，38（S1）：2724-2733.

ZHU Quanqi，LI Diyuan，LI Xibing.Experimental study on the deformation and failure mechanical characteristics of marble with prefabricated elliptical holes[J].Chinese Journal of Rock Mechanics and Engineering，2019，38（S1）：2724-2733.

[21]宋義敏，鄧琳琳，呂祥鋒，等.巖石摩擦滑動變形演化及聲發(fā)射特征研究[J].巖土力學，2019，40（8）：2899-2913.SONG Yimin，DENG Linlin，LYU Xiangfeng，et al.Research on the evolution of rock friction and sliding deformation and acoustic emission characteristics[J].Rock and Soil Mechanics，2019，40（8）：2899-2913.

[22]宋義敏，姜耀東，馬少鵬，等.巖石變形破壞全過程的變形場和能量演化研究[J].巖土力學，2012，33（5）：1352-1365.SONG Yimin，JIANG Yaodong，MA Shaopeng，et al.Research on deformation field and energy evolution during the whole process of rock deformation and failure[J].Rock and soil mechanics，2012，33（5）：1352-1365.

[23]

LI X Y，ZENG J，DING Y J，et al.Research on lateral scale effect and constitutive model of rock damage energy evolution[J].Geotechnical and Geological Engineering，2018，36（4）：2415-2424.

[24]MIAO S J，LAI X P，ZHAO X G，et al.Simulation experiment of AEbased localization damage and deformation characteristic on covering rock in minedout area[J].International Journal of Minerals Metallurgy and Materials，2009，16（3）：255-260.

[25]王學濱，潘一山，馬瑾.剪切帶內部應變（率）分析及基于能量準則的失穩(wěn)判據(jù)[J].工程力學，2003，20（2）：111-115.WANG Xuebin，PAN Yishan，MA Jin.Analysis of strain（rate）in shear zone and instability criterion based on energy criterion[J].Engineering Mechanics，2003，20（2）：111-115.

[26]宋義敏，邢同振，趙同彬，等.巖石單軸壓縮變形場演化的聲發(fā)射特征研究[J].巖石力學與工程學報，2017，36（3）：534-542.SONG Yimin，XING Tongzhen，ZHAO Tongbin，et al.Acoustic emission characteristics of rock uniaxial compression deformation field evolution[J].Chinese Journal of Rock Mechanics and Engineering，2017，36（3）：534-542.

[27]唐禮忠，宋徉霖.含非共面重疊型微裂隙類巖石試樣單軸受壓宏細觀力學特性顆粒流模擬[J].巖石力學與工程學報，2019，38（11）：2161-2171.

TANG Lizhong，SONG Yulin.Particle flow simulation of macroand mesomechanical characteristics of rock specimens containing noncoplanar overlapping microfractures under uniaxial compression[J].Chinese Journal of Rock Mechanics and Engineering，2019，38（11）：2161-2171.