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

(1.西安科技大學(xué) 安全科學(xué)與工程學(xué)院，陜西 西安 710054；2.西安科技大學(xué) 能源學(xué)院，陜西 西安 710054；3.西安科技大學(xué) 理學(xué)院，陜西 西安710054)

0 引 言

煤巖體中微裂紋結(jié)構(gòu)和變形局部化行為會造成其失穩(wěn)破壞，從而導(dǎo)致嚴(yán)重的礦井災(zāi)害事故[1]。目前針對堅硬煤巖變形局部化演化已取得大量成果。而軟煤作為一種特殊煤巖，其力學(xué)特性和變形特征與堅硬煤巖有較大區(qū)別[2]，因此從變形局部化角度研究軟煤失穩(wěn)破壞過程具有重要意義。

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

文中將石膏與水以質(zhì)量比為7∶3混合制成含裂隙軟煤試樣，進行單軸壓縮試驗并采用數(shù)字散斑相關(guān)測量方法(digital speckle correlation measurement method，DSCM)獲取試樣加載過程中表面應(yīng)變與位移信息。分析不同傾角試樣表面變形局部化演化特征及對應(yīng)裂紋起裂、擴展?fàn)顟B(tài)，分階段探究不同傾角試樣局部化帶位移和能量演化規(guī)律。

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

1.1 試樣制備

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

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

表1 軟煤單軸壓縮結(jié)果Table 1 Soft coal sample specification parameters in uniaxial compression test

表2 試樣尺寸及試驗參數(shù)Table 2 Sample size and test parameters

1.2 試驗系統(tǒng)

試驗加載系統(tǒng)采用DNS200電子萬能試驗機，加載方式為位移控制，加載速率為0.5 mm/min，數(shù)據(jù)采集頻率為1 Hz。為降低端部效應(yīng)影響，在試樣與上下壓頭之間涂抹凡士林。

試驗采用DSCM系統(tǒng)捕捉全場變形圖像信息[15]。該測量系統(tǒng)包括光源、CCD攝像機、圖像采集系統(tǒng)和VIC-3D軟件。試驗時將測量物體放置在加載設(shè)備中，由圖像采集系統(tǒng)控制采樣頻率和分辨率，由CCD攝像機采集圖像，通過VIC-3D軟件計算試樣加載過程中全場應(yīng)變、位移信息。為防止頻閃對數(shù)據(jù)采集的影響，用兩架OSRAM?55460光學(xué)儀器燈對物體表面進行照射[16]。試驗中圖像采集速率為1.0幀/s，圖像分辨率為1 600 pixel×1 200 pixel，物面分辨率為0.1 316 mm/pixel。系統(tǒng)構(gòu)成如圖1所示。試驗要求加載系統(tǒng)與圖像采集系統(tǒng)時間同步，開始加載的同時采用VIC-3D軟件的連續(xù)采集功能采集試樣表面變形圖像，直至試樣發(fā)生破壞。

圖1 試驗系統(tǒng)構(gòu)成Fig.1 Test system composition diagram

1.3 試驗結(jié)果

通過試驗得到D0與D60組試樣應(yīng)力-應(yīng)變曲線，如圖2(a)和(b)所示。同組裂隙傾角試樣均一性較好。

D0組試樣在0.007<ε<0.009內(nèi)達(dá)到峰值強度，其峰值強度范圍在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內(nèi)達(dá)到峰值強度，其峰值強度范圍在7.60～8.70 MPa之間，平均值8.10 MPa，最大相差1.08 MPa，彈性模量在742.00～786.63 MPa之間，平均值766.75 MPa，最大相差44.63 MPa。

圖2 試樣應(yīng)力-應(yīng)變曲線Fig.2 Stress-strain curves of samples

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

2 試驗結(jié)果分析

2.1 變形局部化區(qū)域劃分

采用主應(yīng)變差值法[18]來確定局部化區(qū)域，即劃分區(qū)的主應(yīng)變差值較周圍區(qū)域呈現(xiàn)明顯局部性集中。如圖5(a)和(b)所示，根據(jù)試樣最終破壞模式及破壞前變形場，將變形局部化帶所在位置(也是最終宏觀破壞出現(xiàn)位置)按出現(xiàn)先后順序用Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ標(biāo)注。

為對比分析不同裂隙傾角軟煤試樣變形場演化特征，將D0-2與D60-3試樣加載全程應(yīng)力按照公式(1)進行歸一化處理。同時，以峰值應(yīng)力對應(yīng)的應(yīng)變εc為基準(zhǔn)，采用(2)式對應(yīng)變進行處理。圖3為應(yīng)力-應(yīng)變水平曲線。

(1)

式中σstd為歸一化的應(yīng)力水平(0<σstd<1)；σ為瞬時應(yīng)力，MPa；σc為峰值應(yīng)力,MPa。

(2)

式中εstd為應(yīng)變水平；ε為瞬時應(yīng)變；εc為峰值應(yīng)力對應(yīng)的應(yīng)變。

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

圖3 試樣應(yīng)力-應(yīng)變水平曲線Fig.3 Stress level-strain curves of samples

2.2 裂隙軟煤試樣變形場演化分析

為分析不同裂隙傾角軟煤試樣加載全程變形場演化特征，采用DSCM系統(tǒng)計算得到加載全程應(yīng)變場云圖。圖6為試樣峰值后(標(biāo)試點5)應(yīng)變云圖，其中拉應(yīng)變?yōu)檎?，壓?yīng)變?yōu)樨?fù)。不同顏色代表不同應(yīng)變區(qū)域，紅色代表局部高拉應(yīng)變區(qū)，紫色代表局部高壓應(yīng)變區(qū)[19]。選取圖6應(yīng)變云圖中50 mm×80 mm區(qū)域進行放大，取圖3同一應(yīng)力水平下標(biāo)識點1～5時刻所對應(yīng)云圖為參考，從細(xì)觀層面分析試樣應(yīng)變演化過程。

D0-3試樣應(yīng)變演化如圖4(a)～(d)所示。在0.2σc應(yīng)力水平下，如圖4(a)所示，試樣處于壓密階段，應(yīng)變云圖顏色基本一致，說明試樣表面變形均勻。在0.6σc應(yīng)力水平下，如圖4(b)所示，試樣處于彈性階段，試樣應(yīng)變值較上一應(yīng)力水平有微小增大，預(yù)制裂隙中央位置開始出現(xiàn)變形局部化現(xiàn)象。根據(jù)朱泉企等[20]的大量實驗，該應(yīng)力水平前試樣表面沒有裂紋萌生。繼續(xù)加載至0.79σc應(yīng)力水平下，如圖4(c)所示，此時試樣變形集中區(qū)域范圍和變形量值都有所增加，局部化帶Ⅰ基本穩(wěn)定成形。當(dāng)應(yīng)力水平達(dá)到0.92σc，如圖4(d)所示，變形量最大值達(dá)到0.003 46，變形進一步集中。隨著加載繼續(xù)進行，局部化帶Ⅰ向上端延伸，試樣最終發(fā)生整體破壞。D60-2試樣應(yīng)變演化如圖5(e)～(h)所示。試樣前3個應(yīng)力水平下，如圖5(a)、(b)、(c)所示。與D0-3試樣區(qū)別在于其變形局部化并不明顯。這是由于D60-2試樣結(jié)構(gòu)主要受裂隙面控制，加載前期結(jié)構(gòu)較為穩(wěn)定，未有明顯應(yīng)力集中。當(dāng)應(yīng)力水平達(dá)到0.92σc時，如圖5(d)所示，在預(yù)制裂隙端部出現(xiàn)2個明顯的變形集中區(qū)域，并有交匯趨勢，應(yīng)變量值達(dá)到0.024 4。隨著加載的進行，分別向右上角和左側(cè)擴展貫通，最終導(dǎo)致試樣破壞。

圖4 D0-3試樣應(yīng)變場演化過程Fig.4 Evolution process of strain field of 0°specimen

圖5 D60-2試樣應(yīng)變場演化過程Fig.5 Evolution process of strain field of 60°specimen

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

2.3 裂隙軟煤局部化帶位移演化分析

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

(3)

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

圖6 不同傾角試樣最終破裂形態(tài)與峰值后 變形場云圖對照Fig.6 Comparison of final fracture morphology of samples with different dip angles and the cloud map of deformation field after peak value

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

(4)

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

圖7 測點布置示意Fig.7 Schematic diagram of measuring points arrangement

由圖可以看出，不同裂隙傾角試樣局部化帶位移速率演化過程相似，均存在緩慢演化和加速演化階段。加載初期，局部化帶張開與錯動均不明顯。加載至應(yīng)力峰值附近時，D0-3沿局部化帶Ⅰ發(fā)生張開，Ⅱ、Ⅲ發(fā)生錯動，D60-2試樣沿局部化帶Ⅰ、Ⅱ發(fā)生張開，Ⅲ發(fā)生錯動。應(yīng)力峰值后，由于局部化帶之間擴展、交匯、貫通導(dǎo)致D0-3與D60-2試樣張開與錯動速率加速演化。對比圖8(a)和(b)發(fā)現(xiàn)D0-3試樣錯動速率量值高于錯動速率值一個量級，對比圖8(c)和(d)發(fā)現(xiàn)D60-2試樣張開與錯動速率量級相當(dāng)。這表明D0-3試樣應(yīng)力峰值前以張拉破壞為主，峰值后以張拉剪切混合破壞為主，D60-2試樣整體以張拉剪切混合破壞為主。

由2.2節(jié)應(yīng)變場演化可知加載至應(yīng)力峰值前后，試樣發(fā)生宏觀裂紋擴展，宏觀裂紋擴展導(dǎo)致裂隙試樣剪切錯動速率演化明顯，在此對試樣錯動速率進一步分析。

如圖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-4mm/s；加載時間為59 s(峰后)時，局部化帶Ⅳ與Ⅴ迅速增加至39.9×10-4mm/s，36.7×10-4mm/s，此時試樣發(fā)生失穩(wěn)破壞。D0-3試樣局部化帶錯動速率在峰值前后迅速增加，主要原因是試樣應(yīng)變場在應(yīng)力峰值前后劇烈演化，各局部化帶迅速貫通交匯導(dǎo)致。

2018年年初，16歲華裔天才少女朱易奪得全美花滑錦標(biāo)賽新人組女單冠軍，超出第2名高達(dá)35分，成為令西方矚目的新星。然而，她卻毅然放棄在美國的發(fā)展機會，9月底加入中國國家隊，成為偶像——“冰上蝴蝶”陳露的弟子，力爭沖刺2022年北京冬奧會，為祖國奪金。而在這背后，是“科學(xué)大神”父親的夢想，一直在激勵著她！

圖8(d)為D60-2試樣錯動速率演化曲線。從圖可以得到加載時間為203 s(0.92σc)時，局部化帶Ⅱ錯動速率值為2.7×10-4mm/s；加載時間為225 s時，局部化帶Ⅰ錯動速率達(dá)到最大值4.7×10-4mm/s，加載時間為233(0.98σc)時，局部化帶Ⅳ發(fā)生剪切錯動導(dǎo)致錯動速率迅速增大至13.5×10-4mm/s，加載時間為261 s與271 s(峰值后)時，局部化帶Ⅲ和Ⅴ相繼達(dá)到最大值3.9×10-4mm/s，9.5×10-4mm/s，此時試樣發(fā)生失穩(wěn)破壞。D60-2試樣局部化帶錯動速率在峰后迅速增加，主要原因是D60-2試樣應(yīng)變場在應(yīng)力峰值后劇烈演化，各局部化帶迅速貫通交匯導(dǎo)致。

對比(b)圖與(d)圖，發(fā)現(xiàn)D60-2試樣錯動速率量值明顯低于D0-3試樣。其原因在于，在相同加載條件下，同種材料的裂隙擴展力相同，但裂隙傾角增大，擴展阻力隨之增大[22]，導(dǎo)致其擴展程度降低，表現(xiàn)為D60-2試樣錯動速率量值上低于D0-3試樣一個量級。

2.4 裂隙軟煤試樣變形能量演化分析

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

(5)

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

圖8 位移速率演化曲線Fig.8 Evolution curves of displacement rate

圖9(a)和(b)分別為D0-3與D60-2試樣局部化帶變形能密度演化曲線?？梢钥闯鲈嚇友莼?jīng)歷3各階段，分別能量累積的A階段、能量累積與釋放交替的A-R階段、能量釋放的R階段。在應(yīng)力加載峰值前，由于預(yù)制裂隙試樣處于壓密階段，變形局部化帶以能量積累為主。應(yīng)力峰值前后，變形局部化帶既有能量積累又有釋放。加載至后期能量表現(xiàn)為釋放。由于兩試樣能量演化規(guī)律具有相似性，故在此僅以D0-3試樣演化過程做詳細(xì)分析，其余不予贅述。

D0-3試樣加載至峰值前階段，局部化帶Ⅰ與Ⅱ幾乎同時達(dá)到能量峰值后開始釋放能量，這是由于局部化帶Ⅰ與Ⅱ首先發(fā)生微小變形而產(chǎn)生影響。局部化帶Ⅲ、Ⅳ、Ⅴ還在進行能量積累。在加載應(yīng)力峰值后階段，局部化帶Ⅳ與Ⅴ快速增加，并達(dá)到能量峰值2.02×106J/m3，2.23×106J/m3，同時，局部化帶Ⅰ與Ⅱ能量下降到9.6×105J/m3，6.4×105J/m3，局部化帶Ⅲ能量積累達(dá)到峰值3.25×105J/m3，然后迅速釋放能量，進一步擴展。其中局部化帶Ⅲ能量積累速率最慢，累積能量最高，主要原因是局部化帶Ⅱ向下擴展，以及局部化帶Ⅳ向上擴展均與其連接貫通而對其產(chǎn)生影響，這圖4云圖結(jié)果一致。從圖中可以看出，各局部化帶錯動速率突增時刻與能量密度累積與釋放交替時刻相對應(yīng)，表明裂隙面錯動狀態(tài)變化是引起變形能密度出現(xiàn)波動的主要原因。

圖9 試樣變性能密度Fig.9 Deformation energy density of sample

對比圖9(a)和(b)可看出，D60-2試樣A-R階段相對于D0-3試樣有一定程度后移。產(chǎn)生這種現(xiàn)象的原因是裂隙擴展能量釋放率(即為裂隙擴展力)大于消耗率(即為裂隙擴展阻力)能量，表現(xiàn)為裂隙將擴展[22]。而D60-2試樣擴展阻力較大，相比而言D0-3試樣裂紋更易擴展從而引起能量較早釋放。這與唐禮忠數(shù)值計算出隨裂隙傾角增大，試樣發(fā)生能量釋放相對后移的結(jié)果一致[27]。

3 結(jié) 論

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

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

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