沈 凱，肖曉春，徐 軍

( 1. 紹興文理學(xué)院 土木工程學(xué)院，浙江 紹興 312000；2. 紹興文理學(xué)院 浙江省巖石力學(xué)與地質(zhì)災(zāi)害重點(diǎn)實(shí)驗(yàn)室，浙江 紹興 312000；3. 遼寧工程技術(shù)大學(xué) 力學(xué)與工程學(xué)院，遼寧 阜新 123000 )

隨著煤炭開采深度的不斷增加，出現(xiàn)許多不同于淺部開采的問題。深部開采條件下，在賦存環(huán)境、地質(zhì)構(gòu)造以及開采擾動(dòng)所致的應(yīng)力場(chǎng)遷移和傳遞等因素共同影響下，極易發(fā)生超低摩擦型沖擊地壓，而軟弱夾層是頂板巷道圍巖的關(guān)鍵控制層位，因此，研究軟弱夾層對(duì)超低摩擦效應(yīng)的影響規(guī)律具有重要的實(shí)踐意義。

對(duì)于深部巖體超低摩擦現(xiàn)象，國(guó)內(nèi)外諸多專家學(xué)者做了大量的研究工作。國(guó)外最早由俄羅斯科學(xué)家KURLENYA M V[1-2]等提出深部巖體超低摩擦理論，并開展了塊系巖體受沖擊作用時(shí)的實(shí)驗(yàn)室模擬試驗(yàn)，驗(yàn)證了超低摩擦現(xiàn)象的存在；國(guó)內(nèi)錢七虎院士[3-6]等在研究擺型波產(chǎn)生機(jī)理時(shí)發(fā)現(xiàn)了巖體介質(zhì)中的超低摩擦現(xiàn)象，對(duì)深部巖體工程響應(yīng)提出了一系列新的特征科學(xué)現(xiàn)象；對(duì)于超低摩擦試驗(yàn)研究，吳昊[7]、王洪亮[8-9]、王德榮[10]、許瓊萍[11]等進(jìn)行了深入的研究。在理論研究方面許多學(xué)者提出過不同的理論，吳昊[12]等基于KURLENYA M V的模型試驗(yàn)，得出工作塊體水平位移的解析表達(dá)式，驗(yàn)證了超低摩擦現(xiàn)象的存在；潘一山[13-15]、王凱興[16-19]等基于擺型波理論對(duì)沖擊荷載作用下塊系巖體介質(zhì)動(dòng)力模型進(jìn)行求解，得到各巖塊的加速度響應(yīng)，通過比較得出了塊體中心頻率和頻帶之間的規(guī)律，通過分析動(dòng)力作用下巖塊間的相對(duì)位移最大值，得到了塊系巖體超低摩擦啟動(dòng)發(fā)生的條件；李利萍[20-21]等從塑性動(dòng)力學(xué)角度解釋了巖體超低摩擦效應(yīng)，提出超低摩擦型沖擊地壓這一新概念，并利用FLAC3D模擬得到塊系巖體同時(shí)受垂直沖擊載荷和水平靜力時(shí)的動(dòng)力響應(yīng)。在軟弱夾層研究方面，唐禮忠[22]等對(duì)于動(dòng)力擾動(dòng)下含軟弱夾層圍巖巷道力學(xué)響應(yīng)特性進(jìn)行了研究，發(fā)現(xiàn)有無軟弱夾層對(duì)于巷道的拱頂和底板應(yīng)力波幅值、巷道圍巖變形等都有較大的影響；RUTTER E H[23]等通過試驗(yàn)研究軟弱夾層對(duì)斷層泥摩擦效應(yīng)的影響，發(fā)現(xiàn)軟弱夾層對(duì)斷層超低摩擦滑動(dòng)起重要作用。

上述對(duì)于超低摩擦效應(yīng)主要從試驗(yàn)和理論方面進(jìn)行研究，其中理論方面研究主要對(duì)建立的一維剛性塊系塊體模型進(jìn)行力學(xué)響應(yīng)分析，通過求解運(yùn)動(dòng)方程并對(duì)法向力進(jìn)行分析。但一些深部巖體開采中發(fā)生超低摩擦現(xiàn)象時(shí)存在軟弱夾層，現(xiàn)有研究中關(guān)于弱層對(duì)煤巖體在靜、動(dòng)載荷作用下的動(dòng)力響應(yīng)分析鮮有研究。因此，采用數(shù)值模擬方法，以水平加速度、水平位移、殘余位移等作為表征超低摩擦效應(yīng)的特征量，分析不同加載方式、不同沖擊載荷強(qiáng)度以及弱層對(duì)超低摩擦效應(yīng)影響的規(guī)律，以期為超低摩擦型沖擊地壓提供理論參考。

1 含弱層的組合煤巖數(shù)值模擬

1.1 數(shù)值模型的建立

根據(jù)文獻(xiàn)[7]對(duì)超低摩擦效應(yīng)的試驗(yàn)研究，筆者模擬含弱層煤巖體的超低摩擦效應(yīng)。數(shù)值模型參考遼寧紅陽三礦的超低摩擦型沖擊地壓，模型示意如圖1所示。

圖1 超低摩擦型沖擊地壓示意Fig. 1 Schematic diagram of rock-burst with ultra-low friction

為了研究弱層對(duì)于超低摩擦效應(yīng)的影響，應(yīng)用有限差分軟件FLAC3D對(duì)含弱層煤巖體的超低摩擦效應(yīng)進(jìn)行數(shù)值模擬。模型尺寸為砂巖300 mm( 長(zhǎng) )×120 mm( 寬 )×90 mm( 高 )；弱層300 mm( 長(zhǎng) )×120 mm( 寬 )×10 mm( 高 )；煤塊300 mm( 長(zhǎng) )×120 mm( 寬 )×50 mm( 高 )。塊體之間設(shè)置接觸面，建立三維含弱層組合煤巖的超低摩擦數(shù)值模擬模型，模型的網(wǎng)格劃分與監(jiān)測(cè)點(diǎn)布置如圖2所示。

圖2 數(shù)值模擬模型Fig. 2 Numerical simulation model diagram

監(jiān)測(cè)點(diǎn)p1，p2分別在弱層的上下面，p3，p4，p5在煤塊的中線上，間距50 mm，其中p3在煤塊中心處。

1.2 模型參數(shù)選擇

應(yīng)用FLAC3D中巖石類材料模型，由于FLAC3D中自帶的非線性本構(gòu)模型主要用于準(zhǔn)靜態(tài)和具有單調(diào)性的動(dòng)態(tài)響應(yīng)情況，因此動(dòng)態(tài)響應(yīng)分析選擇摩爾庫侖本構(gòu)模型。在FLAC3D動(dòng)力計(jì)算時(shí)，模型邊界對(duì)入射波產(chǎn)生反射而影響計(jì)算的真實(shí)性，為了減少波的反射選擇靜態(tài)邊界條件，這與理論模型邊界條件一致。煤巖力學(xué)材料參數(shù)見表1[24]。

表1 含弱層組合煤巖力學(xué)參數(shù)[24]Table 1 Mechanical parameters of coal rock with weak beds [24]

2 數(shù)值模擬結(jié)果與分析

2.1 可行性驗(yàn)證

目前對(duì)超低摩擦試驗(yàn)研究主要是對(duì)塊系塊體施加靜、動(dòng)載荷研究其動(dòng)力響應(yīng)，為了驗(yàn)證本文模擬的可靠性，根據(jù)文獻(xiàn)[7]的塊體波動(dòng)特性試驗(yàn)研究數(shù)據(jù)，對(duì)含弱層煤巖體施加垂直動(dòng)載荷與水平靜載荷，在靜、動(dòng)載荷共同作用下，獲得塊體的加速度時(shí)程曲線。

圖3為本文數(shù)值模擬結(jié)果與文獻(xiàn)[7]試驗(yàn)結(jié)果的對(duì)比。由圖3可知，本文數(shù)值模擬結(jié)果與文獻(xiàn)[7]的試驗(yàn)結(jié)果并不完全相同，但都具有類似的變化趨勢(shì)，均在經(jīng)歷了一段時(shí)間的震蕩后趨于穩(wěn)定。因此，試驗(yàn)與數(shù)值模擬結(jié)果對(duì)比驗(yàn)證了本文數(shù)值模擬的可靠性。由于超低摩擦試驗(yàn)與數(shù)值模擬存在的局限性，以及所取模型尺寸的差異性，因此在數(shù)值上存在差異，但整體變化趨勢(shì)相同，驗(yàn)證了本文數(shù)值模擬的合理性和可行性。

圖3 數(shù)值模擬與文獻(xiàn)[7]試驗(yàn)結(jié)果對(duì)比Fig. 3 Comparison between numerical simulation and Reference [7] experiments

2.2 沖擊載荷對(duì)超低摩擦效應(yīng)的影響

在深部開采過程中會(huì)受到人工爆破、周期來壓、礦震等擾動(dòng)，在遇到強(qiáng)沖擊載荷時(shí)極易發(fā)生超低摩擦現(xiàn)象，繼而誘發(fā)沖擊地壓等一系列災(zāi)害，因此研究強(qiáng)沖擊載荷作用下含弱層煤巖體的超低摩擦效應(yīng)具有重要意義。為了研究強(qiáng)沖擊載荷作用下的動(dòng)力響應(yīng)，先對(duì)含弱層煤巖體進(jìn)行垂直沖擊載荷的數(shù)值模擬研究。對(duì)含弱層煤巖體模型施加500，1 000，1 500 N垂直沖擊載荷vp，不同垂直沖擊載荷作用下監(jiān)測(cè)點(diǎn)p1，p2，p3的水平加速度時(shí)程曲線如圖4所示。

圖4 不同垂直沖擊載荷作用下監(jiān)測(cè)點(diǎn)水平加速度時(shí)程曲線Fig. 4 Time history curves of horizontal acceleration of monitoring point under different vertical impact loads

由圖4可知，在不同垂直沖擊載荷作用下頂板和煤層的水平加速度變化是多次劇烈的波動(dòng)，隨后逐漸衰減呈正弦波，且煤塊存在受迫振動(dòng)和自由振動(dòng)2個(gè)階段。由監(jiān)測(cè)點(diǎn)的加速度變化可知，水平加速度峰值沿頂板、弱層和煤層依次減弱，動(dòng)載擾動(dòng)后巷道圍巖應(yīng)力、位移均發(fā)生了顯著變化，模型內(nèi)的垂直應(yīng)力由頂板向底部轉(zhuǎn)移，導(dǎo)致加速度峰值逐漸減??；加速度經(jīng)過弱層后減弱效果更加顯著，說明弱層對(duì)于應(yīng)力波的傳播有一定的衰減作用。水平加速度幅值隨著沖擊載荷強(qiáng)度的增強(qiáng)而增大，在強(qiáng)沖擊載荷作用下煤巖體內(nèi)儲(chǔ)存的彈性能釋放越多，塊體在水平方向越容易發(fā)生滑移，超低摩擦效應(yīng)越顯著。為了進(jìn)一步研究弱層對(duì)超低摩擦效應(yīng)的影響，圖5為在相同沖擊載荷作用下有無弱層時(shí)不同位置處測(cè)點(diǎn)的加速度時(shí)程曲線。

由圖5可知，無弱層煤巖體在較短時(shí)間內(nèi)完成振蕩過程，隨后達(dá)到穩(wěn)定；有弱層煤巖體則在較長(zhǎng)時(shí)間內(nèi)完成上下振蕩，隨后趨于穩(wěn)定，這是由于弱層具有一定的流變特性，在應(yīng)力擾動(dòng)下產(chǎn)生塑性變形以及蠕變，使應(yīng)力波的傳播速度有所下降，并使應(yīng)力波傳播時(shí)間延長(zhǎng)，需要經(jīng)歷較長(zhǎng)時(shí)間才能達(dá)到穩(wěn)定的平衡。含弱層煤巖體中的煤塊水平加速度幅值明顯低于無弱層煤巖體，應(yīng)力波經(jīng)過弱層后頂板與煤層加速度峰值的差值變小，說明弱層的存在對(duì)于應(yīng)力波有一定的衰減作用。

圖5 相同沖擊載荷作用下有無弱層時(shí)監(jiān)測(cè)點(diǎn)加速度時(shí)程曲線Fig. 5 Time history curves of acceleration at monitoring points with or without weakly layers under the same impact loads

2.3 靜載對(duì)超低摩擦效應(yīng)的影響

在對(duì)含弱層煤巖體進(jìn)行靜動(dòng)載響應(yīng)模擬前，需要對(duì)作用在煤塊上的最大靜摩擦力進(jìn)行測(cè)定。對(duì)數(shù)值模型依次施加水平靜力(F)100～1 400 N，圖6為在不同水平靜力作用下煤塊水平位移時(shí)程曲線。

圖6 不同水平靜力作用下水平位移時(shí)程曲線Fig. 6 Time history curves of horizontal displacement under different horizontal static forces

由圖6可知，僅水平靜力作用時(shí)，水平位移變化會(huì)迅速達(dá)到峰值然后趨于穩(wěn)定，最初煤塊產(chǎn)生的水平位移較小。當(dāng)水平靜力較小時(shí)煤塊的水平位移趨于穩(wěn)定的時(shí)間基本一致；當(dāng)水平靜力達(dá)到1 200 N時(shí)，煤塊水平位移達(dá)到穩(wěn)定的時(shí)間提前了近1 ms，且產(chǎn)生了較大的殘余位移。

圖7為煤塊水平殘余位移時(shí)程曲線。由圖7可知，水平靜力為100～300 N時(shí)，煤塊的水平殘余位移幾乎沒有變化，水平靜力為100 N時(shí)的水平殘余位移僅為0.594 μm；水平靜力為400～1 200 N時(shí)，水平殘余位移的增長(zhǎng)速度加快，水平靜力為1 200 N時(shí)的水平殘余位移達(dá)10.70 μm，說明煤塊與相鄰塊體間產(chǎn)生了相對(duì)滑動(dòng)，但塊體間的摩擦力仍為靜摩擦力；水平靜力達(dá)1 200 N后水平殘余位移呈指數(shù)增長(zhǎng)，水平靜力由1 200 N增加至1 300 N時(shí)，水平殘余位移增長(zhǎng)了42%，監(jiān)測(cè)點(diǎn)的水平殘余位移發(fā)生了突變，煤塊發(fā)生了整體滑移，此時(shí)煤塊與弱層之間發(fā)生的摩擦力為滑動(dòng)摩擦力，因此，可知最大水平靜摩擦力為1 200～1 300 N。

圖7 水平靜力作用下煤塊水平殘余位移時(shí)程曲線Fig. 7 Time history curve of coal residual displacement under horizontal static force

2.4 靜、動(dòng)載荷對(duì)超低摩擦效應(yīng)的影響

為了研究含弱層煤巖體在靜、動(dòng)載荷作用下的超低摩擦效應(yīng)，先對(duì)煤塊施加水平靜力F，使之達(dá)到平衡后在頂板巖塊中心位置施加半正弦波形式的垂直沖擊載荷pv，其作用形式為F(t) =pvsinωt，其中，t為作用時(shí)間；ω為沖擊載荷頻率。為了研究水平靜力對(duì)煤塊超低摩擦效應(yīng)的影響，對(duì)數(shù)值模型施加相同垂直沖擊載荷與不同的水平靜力，監(jiān)測(cè)點(diǎn)p3的水平位移變化曲線及位移云圖如圖8所示。

圖8 靜、動(dòng)載荷作用下煤塊水平位移曲線及位移云圖Fig. 8 Horizontal displacement curve and cloud image of coal block under dynamic and static combination action

由圖8可知，靜、動(dòng)載荷作用下煤塊水平方向的位移首先經(jīng)歷了“上下”周期振動(dòng)，隨后振動(dòng)幅值減小逐漸趨于新的平衡并產(chǎn)生殘余位移；當(dāng)水平靜力較小時(shí)，煤塊主要受垂直沖擊載荷的作用，此時(shí)煤塊的殘余位移較小(F=100 N 時(shí)殘余位移僅0.594 μm )，此時(shí)煤塊水平位移響應(yīng)曲線如圖8( a )所示，最大水平位移僅0.64 μm，為了解整個(gè)煤巖塊體水平位移變化情況，得到計(jì)算結(jié)束時(shí)煤巖塊體的水平位移云圖為圖8( b )，位移變化從煤塊中部向四周擴(kuò)散分布，大致呈帶狀對(duì)稱分布，在煤塊與弱層接觸的部位有較大水平位移，此處最易發(fā)生超低摩擦現(xiàn)象，此現(xiàn)象與文獻(xiàn)[7]做的超低摩擦試驗(yàn)所得結(jié)論一致，這是因?yàn)樵谒届o力較小時(shí)煤巖體內(nèi)應(yīng)力分布較為均勻；水平靜力增大至特征摩擦力Fx0時(shí)，水平位移曲線未衰減至0，而是出現(xiàn)了殘余位移，殘余位移隨著水平靜力的增大而增大；當(dāng)F=300 N時(shí)，水平殘余位移為1.84 μm；當(dāng)F=500 N時(shí)，水平殘余位移為3.10 μm，水平位移最大值為3.19 μm，煤巖塊體水平位移曲線和計(jì)算結(jié)束時(shí)的水平位移云圖如圖8( e )～( f )所示，此時(shí)水平位移變化不是由中間向兩邊擴(kuò)散，而是整體呈帶狀分布，在底板巖塊中部水平位移出現(xiàn)負(fù)值，這是由于水平靜力增大后開采空間附近移近量增大，對(duì)于煤層的夾持作用更加明顯，沖擊破壞呈整體性；外部擾動(dòng)的增大使弱層蓄積大量的應(yīng)變能，并使弱層的流變特性改變，應(yīng)力沿頂板向深部擴(kuò)展，弱層進(jìn)一步發(fā)生剪切破壞，塊體間的摩擦減弱效應(yīng)更加明顯；當(dāng)水平靜力增大至超低摩擦臨界力cF時(shí)(F=1 000 N )，煤塊水平位移逐漸增加且運(yùn)動(dòng)速率加快，水平位移動(dòng)力響應(yīng)曲線如圖8( g )所示，由于煤塊受到上下塊體周期性的壓緊與脫離，在脫離狀態(tài)時(shí)塊體間接觸最弱，此時(shí)最易發(fā)生滑移，在壓緊狀態(tài)時(shí)水平位移變化緩慢，出現(xiàn)了階梯式上升變化，此條件下整個(gè)煤巖塊體水平位移云圖如圖8( h )所示，水平位移呈帶狀分布，從左至右水平位移逐漸增大，煤塊的最大水平位移出現(xiàn)在與弱層交界處，最大值為6.44 μm。

綜上分析可知，在垂直沖擊載荷和水平靜力共同作用下，水平靜力具有2個(gè)特征值：一是特征摩擦力Fx0，即在垂直沖擊載荷一定時(shí)，煤塊開始產(chǎn)生殘余位移時(shí)的水平靜力；二是超低摩擦臨界力cF，即在一定垂直沖擊載荷作用下，煤塊開始發(fā)生滑動(dòng)的水平靜力。因此，可以把煤塊的運(yùn)動(dòng)分為3個(gè)階段：

( 1 ) 當(dāng)0 ≤F＜Fx0時(shí)，塊體間為靜摩擦力，煤塊此時(shí)發(fā)生水平移動(dòng)。

( 2 ) 當(dāng)Fx0≤F≤Fc時(shí)，由于塊體間摩擦作用效應(yīng)減弱，煤塊開始有相對(duì)運(yùn)動(dòng)趨勢(shì)，并且緩慢移動(dòng)，最后趨于穩(wěn)定。

( 3 ) 當(dāng)F>Fc時(shí)，此時(shí)塊體間的摩擦效應(yīng)大大減弱，煤塊滑移速度加快，發(fā)生滑動(dòng)。

為了分析不同加載方式對(duì)于超低摩擦效應(yīng)的影響，將水平靜載荷作用與靜、動(dòng)載荷同時(shí)加載時(shí)煤塊的水平位移曲線進(jìn)行對(duì)比分析，圖9為煤塊在不同加載方式下水平殘余位移變化情況。

圖9 不同加載方式下水平殘余位移變化情況Fig. 9 Variation of horizontal residual displacement under different loading modes

不同加載方式下巷道圍巖的應(yīng)力重分布以及應(yīng)力傳遞路徑都會(huì)發(fā)生改變，產(chǎn)生沖擊載荷后應(yīng)力向煤層深部轉(zhuǎn)移顯著，使得深部煤層大面積失穩(wěn)而發(fā)生沖擊地壓。由圖9可知，水平位移變化整體呈對(duì)數(shù)增長(zhǎng)，靜、動(dòng)載荷同時(shí)作用時(shí)的水平殘余位移顯著大于水平靜力作用，當(dāng)水平靜力從100 N依次增大至1 000 N時(shí)，水平殘余位移差值分別為0.264，0.735，1.230，0.790，1.020，1.230，1.470，1.780，2.680 μm，差值隨著水平靜力的增大也逐漸增大，由此可知，水平靜力是超低摩擦現(xiàn)象發(fā)生的充分條件，而垂直沖擊載荷充當(dāng)了“催化劑”的角色，其通常會(huì)改變塊體間接觸行為，使塊體間摩擦減弱效應(yīng)增強(qiáng)，加速煤塊的滑移，水平靜力越大超低摩擦效應(yīng)越明顯。

3 超低摩擦啟動(dòng)條件

在動(dòng)載應(yīng)力波作用下，當(dāng)垂直應(yīng)力波對(duì)煤巖體起到拉伸作用時(shí)，頂?shù)装迮c煤層間的摩擦作用最弱，此時(shí)產(chǎn)生超低摩擦現(xiàn)象，在沖擊載荷引起的水平推力作用下煤體呈大位移滑移突變而誘發(fā)沖擊。因此，外界動(dòng)載荷對(duì)于超低摩擦效應(yīng)影響較大，筆者通過煤塊的水平位移變化這個(gè)特征量來研究超低摩擦的啟動(dòng)條件，根據(jù)文獻(xiàn)[21]對(duì)于塊系塊體的超低摩擦效應(yīng)研究可知，影響超低摩擦效應(yīng)的因素主要有水平靜力、沖擊載荷和沖擊載荷頻率。筆者在改變沖擊載荷大小而保持其他因素不變的情況下研究其動(dòng)力響應(yīng)，通過監(jiān)測(cè)煤塊水平位移變化得出超低摩擦效應(yīng)的啟動(dòng)條件。

以水平靜力F=500 N為例，對(duì)煤巖體弱層數(shù)值模型施加p(t) =pvsin 2 000πt形式的垂直沖擊載荷，記錄監(jiān)測(cè)點(diǎn)水平位移變化。不同沖擊載荷下監(jiān)測(cè)點(diǎn)水平位移變化如圖10所示。

由圖10可知，不同沖擊載荷作用下煤塊產(chǎn)生多次劇烈波動(dòng)，隨后驟變衰減為正弦波，時(shí)間節(jié)點(diǎn)分界非常顯著，煤塊經(jīng)歷了受迫振動(dòng)和自由振動(dòng)2個(gè)階段，受迫振動(dòng)的時(shí)間隨著沖擊載荷的增大而延長(zhǎng)，在深部煤巖體內(nèi)最初儲(chǔ)存了大量的應(yīng)變能，在沖擊載荷作用下大量應(yīng)變能釋放形成破壞沖擊，沖擊載荷越大，水平位移在受迫振動(dòng)階段極值點(diǎn)越多，振動(dòng)頻率越快，受迫振動(dòng)時(shí)間延長(zhǎng)，這與李利萍[21]等的研究成果相符。由圖10可知，水平位移的峰值在沖擊載荷較小時(shí)變化不大，這也印證了上述沖擊載荷相當(dāng)于“催化劑”的結(jié)論，水平靜力是發(fā)生超低摩擦的主要因素，只有當(dāng)沖擊載荷增大至較大時(shí)會(huì)引起較大的水平位移，如圖10( d )所示當(dāng)沖擊載荷增大至800 N時(shí)煤塊發(fā)生較大水平位移，煤塊水平位移呈階梯式上升且將開始滑移，此時(shí)發(fā)生了滑動(dòng)摩擦。

圖10 不同沖擊載荷下煤塊水平位移時(shí)程曲線Fig. 10 Time history curves of coal horizontal displacement under different impact amplitudes

為了研究煤巖體弱層超低摩擦效應(yīng)的啟動(dòng)條件，圖11給出了不同沖擊載荷下的水平殘余位移關(guān)系曲線。

圖11 不同沖擊載荷下水平殘余位移Fig. 11 Horizontal residual displacements under different impact amplitudes

由圖11可知，水平殘余位移隨著沖擊載荷的增強(qiáng)不斷增大，并呈對(duì)數(shù)增長(zhǎng)，說明沖擊能量增大使得深部巖體中儲(chǔ)存的大量應(yīng)變能釋放，更多的應(yīng)變能轉(zhuǎn)化為動(dòng)能，產(chǎn)生較大的水平位移，沖擊載荷也會(huì)使塊體間更易相互脫離，此時(shí)超低摩擦效應(yīng)越強(qiáng)，這與實(shí)際工程中外界擾動(dòng)越大越容易發(fā)生摩擦滑移現(xiàn)象一致。當(dāng)沖擊載荷在200 N內(nèi)時(shí)殘余位移增加較明顯；當(dāng)沖擊載荷大于200 N后殘余位移近似線性增加，且增加較為緩慢。沖擊載荷從50 N逐級(jí)增加至1 000 N時(shí)，殘余位移分別增加3.4%，12.0%，3.0%，1.8%，2.4%，1.7%，1.8%，1.3%，1.2%，0.5%。沖擊載荷由100 N增加至200 N時(shí)，殘余位移增長(zhǎng)量遠(yuǎn)遠(yuǎn)大于其他增加量，在這個(gè)變化過程中煤塊剛發(fā)生摩擦滑移，此時(shí)即為超低摩擦的啟動(dòng)時(shí)刻，隨著沖擊載荷的不斷增大，煤塊與相鄰塊體的脫離越來越顯著，開始慢慢滑動(dòng)最后發(fā)生顯著的大位移整體傾出，即超低摩擦型沖擊地壓。

4 結(jié) 論

( 1 ) 在垂直沖擊載荷作用下監(jiān)測(cè)點(diǎn)加速度呈準(zhǔn)正弦周期變化，煤塊的變形強(qiáng)度隨沖擊載荷強(qiáng)度的增強(qiáng)而增大，且煤塊發(fā)生微小的水平位移；通過對(duì)比有無弱層時(shí)數(shù)值模型的動(dòng)力響應(yīng)發(fā)現(xiàn)，弱層的存在改變了煤巖體的變形特性，由于弱層的流變特性，在應(yīng)力擾動(dòng)下發(fā)生塑性變形以及蠕變使得應(yīng)力波的傳播速度有所下降，且弱層對(duì)應(yīng)力波具有一定的消減作用。

( 2 ) 在水平靜力作用下煤塊水平位移峰值隨著水平靜力的增加不斷增大，在水平靜力達(dá)到一定值時(shí)煤塊提前發(fā)生滑移，通過分析水平殘余位移隨水平靜力變化的情況，得出煤塊最大靜摩擦力的特征區(qū)間為1 200～1 300 N，此時(shí)煤塊將要發(fā)生滑動(dòng)。

( 3 ) 在靜、動(dòng)載荷同時(shí)作用下煤塊經(jīng)歷了受迫振動(dòng)和自由振動(dòng)2個(gè)階段，經(jīng)過一段時(shí)間后趨于穩(wěn)定，煤塊運(yùn)動(dòng)狀態(tài)大致分3個(gè)階段；通過對(duì)比有無弱層時(shí)數(shù)值模型動(dòng)力響應(yīng)，得出水平方向力是超低摩擦發(fā)生的主要原因，而垂直沖擊載荷只是改變塊體間摩擦減弱效應(yīng)，充當(dāng)超低摩擦中的“催化劑”，應(yīng)力波傳遞過程中主要向垂直方向傳遞，弱層的黏彈性質(zhì)會(huì)改變煤巖體整體剛度，導(dǎo)致塊體間摩擦效應(yīng)減弱，煤塊加速發(fā)生滑移突變，且弱層對(duì)應(yīng)力波具有一定的消減作用。

( 4 ) 煤塊的水平殘余位移隨著沖擊載荷的增大而增大，呈對(duì)數(shù)變化，增長(zhǎng)率呈先增大后減小的規(guī)律，沖擊載荷在100～200 N變化時(shí)，煤塊發(fā)生滑移突變，此特征區(qū)間即為超低摩擦的啟動(dòng)條件。

參考文獻(xiàn)( References )：

[1] KURLENYA M V，OPARIN V N，VOSTRIKOV V I. Effect of anomalously low friction in block media[J]. Journal of Applied Mechanics & Technical Physics，1999，40( 6 )：1116-1120.

[2] KURLENYA M V，OPARIN V N，BALMASHNOVA E G. On dynamic behavior of ''self-stressed'' block media，Part Ⅱ：Comparison of theoretical and experimental data[J]. Journal of Mining Science，2001，37( 5 )：455-461.

[3] 錢七虎. 非線性巖石力學(xué)的新進(jìn)展——深部巖體力學(xué)的若干關(guān)鍵問題[A]. 中國(guó)巖石力學(xué)與工程學(xué)會(huì)，第八次全國(guó)巖石力學(xué)與工程學(xué)術(shù)大會(huì)論文集[C]. 成都：中國(guó)巖石力學(xué)與工程學(xué)會(huì)，2004.

[4] 錢七虎. 深部地下空間開發(fā)中的關(guān)鍵科學(xué)問題[A]. 第230次香山科學(xué)會(huì)議：深部地下空間開發(fā)中的基礎(chǔ)研究關(guān)鍵技術(shù)問題[C].北京，2004.

[5] 錢七虎. 深部巖體工程響應(yīng)的特征科學(xué)現(xiàn)象及“深部”的界定[J]. 東華理工學(xué)院學(xué)報(bào)，2004，27( 1 )：1-5.QIAN Qihu. The characteristic scientific phenomena of engineering response to deep rock mass and the implication of ''deepness''[J].Journal of East China Institute of Technology，2004，27( 1 )：1-5.

[6] 王明洋，戚承志，錢七虎. 深部巖體塊系介質(zhì)變形與運(yùn)動(dòng)特性研究[J]. 巖石力學(xué)與工程學(xué)報(bào)，2005，24( 16 )：2825-2830.WANG Mingyang，QI Chengzhi，QIAN Qihu. Study on deformation and motion characteristics of blocks in deep rock mass[J]. Chinese Journal of Rock Mechanics and Engineering，2005，24( 16 )：2825-2830.

[7] 吳昊，方秦，張亞棟，等. 一維塊系地質(zhì)塊體波動(dòng)特性的試驗(yàn)和理論研究[J]. 巖土工程學(xué)報(bào)，2010，32( 4 )：600-611.WU Hao，F(xiàn)ANG Qin，ZHANG Yadong，et al. Propagation properties of stress waves in one-dimensional geo-block medium[J]. Chinese Journal of Geotechnical Engineering，2010，32( 4 )：600-611.

[8] 王洪亮，葛濤，王德榮，等. 塊系巖體動(dòng)力特性理論與實(shí)驗(yàn)對(duì)比分析[J]. 巖石力學(xué)與工程學(xué)報(bào)，2007，26( 5 )：951-958.WANG Hongliang，GE Tao，WANG Derong，et al. Comparison of theoretical and experimental analyses of dynamic characteristics of block mass[J]. Chinese Journal of Rock Mechanics and Engineering，2007，26( 5 )：951-958.

[9] 王洪亮，周澤平，王明洋. 物體界面動(dòng)力性能的研究[J]. 防護(hù)工程，2006，28( 2 )：30-34.WANG Hongliang，ZHOU Zeping，WANG Mingyang. Study on dynamic properties of object interface[J]. Protective Engineering，2006，28( 2 )：30-34.

[10] 王德榮，陸渝生，馮淑芳. 深部巖體動(dòng)態(tài)特性多功能試驗(yàn)系統(tǒng)的研制[J]. 巖石力學(xué)與工程學(xué)報(bào)，2008，27( 3 )：601-606.WANG Derong，LU Yusheng，F(xiàn)ENG Shufang. Development of multipurpose test system for dynamic behaviors of deep rock mass[J].Chinese Journal of Rock Mechanics and Engineering，2008，27( 3 )：601-606.

[11] 許瓊萍，陸渝生，王德榮. 深部巖體塊系摩擦減弱效應(yīng)試驗(yàn)[J].解放軍理工大學(xué)學(xué)報(bào)( 自然科學(xué)版 )，2009，10( 3 )：285-289.XU Qiongping，LU Yusheng，WANG Derong. Experimental study on friction weakened effect of deep block rock mass[J]. Journal of PLA University of Science and Technology( Natural Science Edition )，2009，10( 3 )：285-289.

[12] 吳昊，方秦，王洪亮. 深部塊系巖體超低摩擦現(xiàn)象的機(jī)理分析[J].巖土工程學(xué)報(bào)，2008，30( 5 )：769-775.WU Hao，F(xiàn)ANG Qin，WANG Hongliang. Mechanism of anomalously low friction phenomenon in deep block rock mass[J]. Chinese Journal of Geotechnical Engineering，2008，30( 5 )：769-775.

[13] 潘一山，王凱興. 巖塊尺度對(duì)擺型波傳播影響研究[J]. 巖石力學(xué)與工程學(xué)報(bào)，2012，31( S2 )：3459-3465.PAN Yishan，WANG Kaixing. Study effect of block-lock scale on pendulum-type wave propagation[J]. Chinese Journal of Rock Mechanics and Engineering，2012，31( S2 )：3459-3465.

[14] 潘一山，王凱興，肖永惠. 基于擺型波理論的防沖支護(hù)設(shè)計(jì)[J].巖石力學(xué)與工程學(xué)報(bào)，2013，32( 8 )：1537-1543.PAN Yishan，WANG Kaixing，XIAO Yonghui. Design of anti-scour support based on theory of pendulum-type wave[J]. Chinese Journal of Rock Mechanics and Engineering，2013，32( 8 )：1537-1543.

[15] 潘一山，王凱興. 巖體間超低摩擦發(fā)生機(jī)理的擺型波理論[J]. 地震地質(zhì)，2014，36( 3 )：833-844.PAN Yishan，WANG Kaixing. Pendulum-type waves theory on the mechanism of anomalously low friction mechanism between rock masses[J]. Seismology and Geology，2014，36( 3 )：833-844.

[16] 王凱興，潘一山. 擺型波傳播過程塊系巖體頻域響應(yīng)反演巖塊間黏彈性性質(zhì)[J]. 煤炭學(xué)報(bào)，2013，38( S1 )：19-24.WANG Kaixing，PAN Yishan. Frequency domain response of blockrock mass inversion partings viscoelastic property on pendulum-type wave propagation[J]. Journal of China coal society，2013，38( S1 )：19-24.

[17] 王凱興，孟村影，楊月，等. 塊系覆巖中擺型波傳播對(duì)巷道支護(hù)動(dòng)力響應(yīng)影響[J]. 煤炭學(xué)報(bào)，2014，39( 2 )：347-352.WANG Kaixing，MENG Cunying，YANG Yue，et al. Dynacmic response of roadway support on pendulum type waves propagation in overburden block rock mass[J]. Journal of China coal society，2014，39( 2 )：347-352.

[18] 王凱興，潘一山，曾祥華，等. 塊系巖體間黏彈性性質(zhì)對(duì)擺型波傳播的影響[J]. 巖土力學(xué)，2013，34( S2 )：174-179.WANG Kaixing，PAN Yishan，ZENG Xianghua，et al. Effect of viscoelasticity in block-rock mass partings to the propagation of pendulum waves[J]. Rock and Soil Mechanics，2013，34( S2 )：174-179.

[19] 王凱興，潘一山，竇林名. 擺型波傳播過程塊系巖體能量傳遞規(guī)律研究[J]. 巖土工程學(xué)報(bào)，2016，38( 12 )：2309-2314.WANG Kaixing，PAN Yishan，DOU Linming. Energy transfer in blockrock mass during propagation of pendulum-type waves[J]. Chinese Journal of Geotechnical Engineering，2016，38( 12 )：2309-2314.

[20] 李利萍，潘一山，章夢(mèng)濤. 基于簡(jiǎn)支梁模型的巖體超低摩擦效應(yīng)理論分析[J]. 巖石力學(xué)與工程學(xué)報(bào)，2009，28( S1 )：2715-2720.LI Liping，PAN Yishan，ZHANG Mengtao. Theoretical analysis of effect of anomalously low friction on rock mass based on simply supported beam model[J]. Chinese Journal of Rock Mechanics and Engineering，2009，28( S1 )：2715-2720.

[21] 李利萍，王淋，邰英樓. 動(dòng)靜組合加載作用下塊系巖體動(dòng)力響應(yīng)分析[J]. 中國(guó)安全生產(chǎn)科學(xué)技術(shù)，2016，12( 11 )：5-10.LI Liping，WANG Lin，TAI Yinglou. Dynamic response analysis on block rock mass under combined static-dynamic load[J]. Journal of Safety Science and Technology，2016，12( 11 )：5-10.

[22] 唐禮忠，高龍華，王春，等. 動(dòng)力擾動(dòng)下含軟弱夾層巷道圍巖力學(xué)響應(yīng)特性[J]. 科技導(dǎo)報(bào)，2014，32( 21 )：56-61.TANG Lizhong，GAO Longhua，WANG Chun，et al. Mechanical response features of roadway surrounding rock with weak interlayer under dynamic disturbance[J]. Science and Technology Review，2014，32( 21 )：56-61.

[23] RUTTER E H，HACKSTONA J，YEATMAN E. Reduction of friction on geological faults by weak-phase smearing[J]. Journal of Structural Geology，2013，51：52-60.

[24] 何富連. 巨高巷道幫不穩(wěn)定三角塊滑移及控制[J]. 巖石力學(xué)與工程學(xué)報(bào)，2015，31( 4 )：3238-3247.HE Fulian. Slip and control of coal side unstable triangle block in super-high roadway[J]. Chinese Journal of Rock Mechanics and Engineering，2015，31( 4 )：3238-3247.