何富連，何文瑞，陳冬冬，謝生榮，李 輝，何常西

(中國礦業(yè)大學(xué)(北京) 能源與礦業(yè)學(xué)院，北京 100083)

工作面上覆基本頂?shù)钠茢辔恢?、破斷順序和破斷形態(tài)對(duì)指導(dǎo)工作面安全回采以及鄰側(cè)區(qū)段巷道的煤柱寬度合理確定等意義顯著[1-2]。

要研究基本頂?shù)钠茢嘁?guī)律，需要明確兩類問題:① 采用梁結(jié)構(gòu)模型，還是板結(jié)構(gòu)模型;② 模型的邊界條件是固支、簡支、彈性基礎(chǔ)邊界還是彈-塑性基礎(chǔ)邊界，這2個(gè)方面選擇對(duì)科學(xué)合理的制定措施至關(guān)重要。

多年來，礦業(yè)科技工作者分別建立了基本頂在固支邊界、簡支邊界、文克爾彈性基礎(chǔ)邊界(錢鳴高院士等根據(jù)煤體抗剪切能力弱等特點(diǎn)提出煤體近似符合文克爾彈性地基假設(shè)[1])等條件下的梁結(jié)構(gòu)模型，據(jù)此分析了基本頂破斷機(jī)理、破斷形態(tài)、破斷時(shí)的反彈壓縮擾動(dòng)，采場初次來壓的結(jié)構(gòu)特征及來壓步距等問題[1,3-7];建立了基本頂在固支、簡支及自由邊界條件下的板結(jié)構(gòu)模型，分析了基本頂?shù)钠茢辔恢?、破斷順序及“O-X”型破斷形態(tài)特征等[8-15];初次破斷前的基本頂四周處于上覆較軟巖層和下伏直接頂與煤層的夾支狀態(tài)，由于煤層的剛度小于甚至遠(yuǎn)小于基本頂?shù)膭偠龋敲聪嗤芰l件下煤層的壓縮變形量大于甚至是遠(yuǎn)大于基本頂?shù)膲嚎s變形量，可見下伏煤層限制基本頂下沉的能力很弱，特別是煤層厚度較大且較軟時(shí)，此時(shí)遠(yuǎn)無法滿足固支邊界條件假設(shè)，為了彌補(bǔ)此缺陷，構(gòu)建并系統(tǒng)研究了文克爾彈性基礎(chǔ)邊界基本頂板結(jié)構(gòu)首采工作面初次破斷、周期破斷、一側(cè)采空與兩側(cè)采空時(shí)的破斷規(guī)律及影響因素等[16-20]，并得到了固支邊界模型得不到的重要結(jié)論，這些研究從不同角度，基于不同的模型假設(shè)得到的研究結(jié)論對(duì)采礦工程實(shí)踐起到了重要的指導(dǎo)作用。

但是，煤層開采后，基本頂懸頂區(qū)域周邊一定深度范圍內(nèi)的煤體必然進(jìn)入塑性狀態(tài)，即全區(qū)域煤體并不能全部滿足彈性基礎(chǔ)邊界條件假設(shè)?？梢?，全面考慮采空區(qū)周邊煤體的彈-塑性變形，才能更符合基本頂初次破斷特征，得到的結(jié)論才能更有效的指導(dǎo)實(shí)踐。筆者通過構(gòu)建彈-塑性基礎(chǔ)邊界條件下的基本頂板結(jié)構(gòu)初次破斷模型，基于有限差分方法研究煤體的塑化程度和塑化范圍、基本頂厚度、彈性模量及未塑化煤體的彈性基礎(chǔ)系數(shù)、懸頂跨度對(duì)基本頂板結(jié)構(gòu)破斷位置、破斷順序及整體形態(tài)特征的影響，對(duì)于認(rèn)識(shí)基本頂初次破斷過程及其結(jié)構(gòu)形態(tài)具有重要價(jià)值。

1 基本頂邊界條件類型

煤層開采后，懸頂區(qū)域的基本頂處于上覆與下伏巖層的夾支狀態(tài)。若基本頂上下覆巖層的剛度為無窮大，可構(gòu)建固支邊界模型，如圖1(a)所示;實(shí)際上，基本頂上下覆巖層的剛度小于甚至遠(yuǎn)小于基本頂?shù)膭偠?，特別是煤層厚度較大且較軟時(shí)是無法滿足固支邊界條件要求的，所以如圖1(b)所示，可假設(shè)基本頂受到上下覆彈性巖層的夾支，即可構(gòu)建彈性基礎(chǔ)邊界模型;對(duì)于實(shí)際的采礦工程問題，支撐基本頂?shù)拿后w淺部必然處于塑性狀態(tài)，而圍巖表層甚至處于破碎狀態(tài)，深部煤體處于彈性狀態(tài)，所以如圖1(c)所示，可構(gòu)建彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)模型，研究所得結(jié)論顯然更符合實(shí)際。

圖1 基本頂邊界條件對(duì)比Fig.1 Comparison ofmain roof boundary conditions

2 彈-塑性邊界基本頂板結(jié)構(gòu)力學(xué)模型

2.1 力學(xué)模型建立

根據(jù)彈性薄板力學(xué)假設(shè)[21]

(1)

式中，h為板厚度,m;l為板短邊長度,m。

一般條件下采場上覆基本頂均滿足上述要求，所以滿足彈性薄板假設(shè)[21-22]。

開采區(qū)域周邊的基本頂四周受到上下覆巖層的夾支，尤其是下伏煤層的基礎(chǔ)系數(shù)是限制基本頂在煤體支撐區(qū)變形與破斷的關(guān)鍵。通常條件下淺部煤體處于塑性狀態(tài)，深部煤體處于彈性狀態(tài)，即支撐基本頂?shù)臑橹苓厪?塑性煤體，而煤體塑化程度和范圍必然對(duì)基本頂?shù)钠茢嗵卣饔杏绊憽?/p>

據(jù)此建立如圖2所示的考慮煤體彈塑性變形的彈-塑性基礎(chǔ)邊界基本頂薄板初次破斷力學(xué)模型，其中，ABCD區(qū)域?yàn)橐呀?jīng)開采而基本頂懸頂未斷區(qū)域、A1B1C1D1區(qū)域之內(nèi)與ABCD區(qū)域之外為塑化煤體支撐區(qū)、A1B1C1D1區(qū)域之外而A2B2C2D2區(qū)域之內(nèi)為彈性煤體支撐區(qū)，邊A1B1,B1C1,C1D1及A1D1為彈-塑性煤體分界線，工作面的長度AB為2a、開采懸頂區(qū)的跨度AD為2b、外邊界A2B2長度為2a2、外邊界A2D2長度為2b2，長邊塑性區(qū)的寬度為bs(即AB與A1B1或CD與C1D1之間的距離)、短邊塑性區(qū)的寬度為as(即AD與A1D1,B1C1與BC之間的距離)，長邊彈-塑性分界邊A1B1長度為2as+2a、短邊彈-塑性分界邊A1D1長度為2bs+2b;開采區(qū)域上覆基本頂承擔(dān)載荷為q，其值為基本頂自重與上覆軟巖層載荷之和。

圖2 彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)力學(xué)模型Fig.2 Mechanical model of main roof structure with elastic-plastic foundation boundary

S懸頂區(qū)的基本頂撓度微分方程為

(2)

式中，ω(x,y)為S區(qū)域基本頂撓度函數(shù);D為基本頂剛度。

S1塑性區(qū)的基本頂撓度微分方程為

(3)

式中，ωs(x,y)為S1區(qū)基本頂撓度函數(shù);ks為塑性區(qū)煤體基礎(chǔ)系數(shù)(支撐系數(shù))。

S2彈性區(qū)的基本頂撓度微分方程為

(4)

式中，ωt(x,y)為S2區(qū)基本頂撓度函數(shù);kt為彈性區(qū)煤體基礎(chǔ)系數(shù)(支撐系數(shù))。

(5)

式中，ν為泊松比;E為基本頂彈性模量，GPa;h為基本頂厚度，m。

上述參數(shù)中，彈性煤體區(qū)的基礎(chǔ)系數(shù)kt為煤體的彈性模量與煤體厚度的比值[7,20];煤體塑性區(qū)的寬度由文獻(xiàn)[1]式(8-5)計(jì)算得到。ABCD邊界線即淺部煤體基礎(chǔ)系數(shù)(支撐系數(shù))設(shè)為k0，k0代表了塑性區(qū)煤體的塑化程度，由于開采區(qū)淺部的煤體必然塑化甚至破碎，所以淺部煤體的基礎(chǔ)系數(shù)k0值較小甚至為0(代表淺部破碎)，而開采區(qū)周邊深部煤體的塑化程度越來越小，即基礎(chǔ)系數(shù)越來越大，直到彈性煤體區(qū)時(shí)煤體基礎(chǔ)系數(shù)達(dá)最大值kt，設(shè)S1塑化區(qū)煤體基礎(chǔ)系數(shù)為ks并滿足關(guān)系式k0

2.2 邊界條件

2.2.1內(nèi)邊的邊界條件

邊DC,BC,AD及AB為基本頂塑性基礎(chǔ)起始邊界(也是已開采區(qū)與煤體區(qū)分界邊)，該邊界上各點(diǎn)撓度同時(shí)滿足懸頂區(qū)S的微分方程(2)與塑化區(qū)S1的微分方程(3)，且各邊上的撓度、截面法向線轉(zhuǎn)角、彎矩與剪力分別連續(xù)，即滿足:

(6)

2.2.2彈-塑性分區(qū)的邊界條件

邊D1C1,B1C1,A1D1及A1B1為基本頂下伏煤體的彈-塑性分界邊，邊上各點(diǎn)撓度同時(shí)滿足塑化區(qū)S1微分方程(3)與彈性區(qū)S2微分方程(4)，且各邊上的撓度、截面法向線轉(zhuǎn)角、彎矩與剪力分別連續(xù)，即滿足:

(7)

2.2.3模型外邊界的邊界條件

圖2中，對(duì)于首采面的基本頂板結(jié)構(gòu)初次破斷來說，開采區(qū)域的外邊界A2B2C2D2距離開采懸頂區(qū)的距離要求不受或者基本不受開采區(qū)ABCD的擾動(dòng)影響，由于外邊界C2D2,B2C2,A2B2及A2D2不受開采擾動(dòng)的影響，那么這些邊不僅滿足撓度為0，而且滿足截面法向線轉(zhuǎn)角為0，即滿足式:

(8)

3 彈-塑性基礎(chǔ)邊界板結(jié)構(gòu)模型計(jì)算法

要研究彈-塑性基礎(chǔ)邊界條件下的基本板結(jié)構(gòu)破斷規(guī)律，就需要得到偏微分方程(2)～(4)在邊界條件(6)～(8)條件下的解，但是求出精確解極為困難，即便是彈性基礎(chǔ)邊界條件下的求解難度也很大[1,17-20]。眾所周知，獲得精確解并不是采礦工程問題所追求的，能求出近似解滿足分析并解決采礦工程實(shí)際需求即可，而有限差分法就是符合該要求的有效方法[17-23]。

3.1 差分節(jié)點(diǎn)編號(hào)

為了便于采用有限差分法具體求解偏微分方程(2)～(4)在邊界條件(6)～(8)條件下的解，需要對(duì)差分節(jié)點(diǎn)進(jìn)行編號(hào)，這樣便于計(jì)算處理，圖3為13節(jié)點(diǎn)編號(hào)圖，中間特征節(jié)點(diǎn)為P點(diǎn)，節(jié)點(diǎn)間距為d。

圖3 差分法節(jié)點(diǎn)編號(hào)Fig.3 Node number of difference method

3.2 撓度方程的差分方程

基于差分節(jié)點(diǎn)編號(hào)圖3可得式(2),(3)及(4)在特征節(jié)點(diǎn)P的差分方程分別為式(9),(10)及(11)。

20ωβ，δ-8(ωβ+1，δ+ωβ-1，δ+ωβ，δ+1+ωβ，δ-1)+

2(ωβ+1，δ+1+ωβ+1，δ-1+ωβ-1，δ+1+ωβ-1，δ-1)+

(9)

ωβ，δ-1)+2(ωβ+1，δ+1+ωβ+1，δ-1+ωβ-1，δ+1+

ωβ-1，δ-1)+ωβ+2，δ+ωβ-2，δ+ωβ，δ+2+ωβ，δ-2=0

(10)

ωβ，δ-1)+2(ωβ+1，δ+1+ωβ+1，δ-1+ωβ-1，δ+1+

ωβ-1，δ-1)+ωβ+2，δ+ωβ-2，δ+ωβ，δ+2+ωβ，δ-2=0

(11)

3.3 外邊界條件方程的差分方程

一般條件下，對(duì)于實(shí)際開采邊界范圍以及采礦工程所要求的精度來說，外邊界的范圍不需要無窮遠(yuǎn)，只要是矩形區(qū)ABCD長邊長度的3～5倍以上[1,17-20]即可近似滿足式(8)要求的條件，且有限邊界對(duì)于具體求解來說是有利的。所以可得外邊界條件式(8)基于特征節(jié)點(diǎn)P的差分方程為式(12)。

(12)

3.4 彈-塑性基礎(chǔ)板結(jié)構(gòu)模型求解及破斷指標(biāo)

根據(jù)基本頂板結(jié)構(gòu)在開采懸頂區(qū)S、塑性區(qū)S1及彈性區(qū)S2的撓度偏微分方程的差分方程(9)～(11)可知，3個(gè)區(qū)域的撓度偏微分方程的差分方程均含有13個(gè)差分節(jié)點(diǎn)且撓度未知，而節(jié)點(diǎn)之間相互關(guān)聯(lián)并通過邊界條件來約束范圍，由此可建立并求解基于邊界條件的基本頂板結(jié)構(gòu)撓度未知的13節(jié)點(diǎn)差分方程組，方程組的解即為各個(gè)節(jié)點(diǎn)的撓度解。顯然，求解差分方程比直接求解偏微分方程要簡單，為了滿足計(jì)算精度要求，需要構(gòu)建足夠數(shù)量的差分方程(設(shè)節(jié)點(diǎn)間距d=0.2)，那么雖然計(jì)算較易，但是方程數(shù)量多，此時(shí)可以采用軟件Matlab[23-24]輔助計(jì)算，具體可采用sparse構(gòu)建系數(shù)為稀疏矩陣聯(lián)合代數(shù)方程組，然后采用函數(shù)gmres函數(shù)求解方程組，從而得到各個(gè)撓度未知節(jié)點(diǎn)的撓度解。

各個(gè)區(qū)域的所有撓度未知節(jié)點(diǎn)的撓度求解出來后，可得分量彎矩值，代入式(13)即可得到各個(gè)節(jié)點(diǎn)的主彎矩值，通過各個(gè)區(qū)域的主彎矩極值與彎矩極限Ms對(duì)比[1,17]，可分析彈-塑性基礎(chǔ)邊界條件下的基本頂板結(jié)構(gòu)是否發(fā)生破斷，破斷時(shí)的位置和破斷順序及影響因素等。

(13)

4 彈-塑性基礎(chǔ)基本頂板結(jié)構(gòu)模型分析

4.1 彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)的主彎矩形態(tài)特征

為系統(tǒng)研究彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)的破斷規(guī)律，選取一特征參數(shù)采用上述計(jì)算方法計(jì)算并著重分析基本頂?shù)闹鲝澗貥O值大小及位置特征。為了研究問題方便，設(shè)塑性區(qū)的煤體基礎(chǔ)系數(shù)ks由淺部基礎(chǔ)系數(shù)k0到彈性煤體區(qū)基礎(chǔ)系數(shù)kt呈正相關(guān)線性增長;取工作面長度132 m(即圖2中的AB長度)、推進(jìn)距離44 m(即圖2中的AD長度或跨度);基本頂?shù)膹椥阅Ａ縀、厚度h、泊松比ν、載荷q分別為32 GPa,6.5 m,0.22,0.32 MPa;煤體淺部基礎(chǔ)系數(shù)k0=0、煤體彈性區(qū)基礎(chǔ)系數(shù)kt為1.5 GN/m3，周邊煤體塑性區(qū)寬度(設(shè)為b0)為4 m。

通過計(jì)算得到各節(jié)點(diǎn)最大與最小主彎矩值并繪制出基本頂全區(qū)域主彎矩形態(tài)特征圖。圖4(a)為彈-塑性基礎(chǔ)邊界基本頂最小主彎矩M3的分布形態(tài)特征云圖，圖4(b)為彈-塑性基礎(chǔ)邊界基本頂最大主彎矩M1的分布形態(tài)特征云圖。

根據(jù)圖4基本頂各個(gè)區(qū)域的最大與最小主彎矩M1與M3的形態(tài)特征和數(shù)值分布規(guī)律可知:

(1)開采區(qū)域長邊與短邊深入彈-塑性煤體區(qū)域基本頂板結(jié)構(gòu)的最大與最小主彎矩M1與M3為負(fù)值，所以深入彈-塑性煤體區(qū)域基本頂板結(jié)構(gòu)的上表面受拉應(yīng)力而下表面受壓應(yīng)力，由巖石抗拉強(qiáng)度遠(yuǎn)小于抗壓強(qiáng)度可知，基本頂深入彈-塑性煤體區(qū)的上表面先于下表面破斷。

(2)開采區(qū)域長邊與短邊區(qū)域的絕對(duì)值最大主彎矩為最小主彎矩M3的相反數(shù)，且位置在深入彈塑性煤體區(qū)域而不是固支邊界模型得到的沿著煤壁。

(3)開采懸頂區(qū)域中部基本頂?shù)淖畲笈c最小主彎矩M1與M3均為正值，所以該區(qū)域基本頂?shù)纳媳砻媸軌簯?yīng)力而下表面受拉應(yīng)力，且下表面先于上表面破斷。

(4)基本頂中部區(qū)域的最大主彎矩在開采區(qū)域的中點(diǎn)(坐標(biāo)為圖2中O點(diǎn)坐標(biāo)(0,0)位置)，且數(shù)值為對(duì)應(yīng)點(diǎn)主彎矩M1。

不同參數(shù)條件下彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)全區(qū)域的最大與最小主彎矩的分布形態(tài)與圖4相似，且開采區(qū)域長邊與短邊主彎矩極值的位置也處于彈-塑性煤體區(qū)，開采區(qū)域中部的最大主彎矩在開采區(qū)的中點(diǎn)，所以要得到彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)的破斷規(guī)律可以通過研究各個(gè)特征區(qū)域的主彎矩極值的大小和位置及隨影響因素的變化特征來分析。

圖4 彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)主彎矩云圖Fig.4 Cloud chart of main bending moment of main roof plate structure with elastic-plastic foundation boundary

根據(jù)上述分析結(jié)果，建立如圖5所示的彈-塑性基礎(chǔ)邊界基本頂主彎矩極值位置圖，其中，中部區(qū)域的主彎矩極值為中點(diǎn)(0,0)位置處的最大主彎矩M1，并設(shè)Msz=M1|(0,0);長邊區(qū)域絕對(duì)值最大主彎矩為主彎矩M3極值的負(fù)數(shù)，其深入煤體距離為Lsc，即距離AB邊或CD邊的長度為Lsc，設(shè)該極值為Msc，那么滿足Msc=|M3|(0，b+Lsc)|與Msc=|M3|(0，-b-Lsc)|;短邊區(qū)域絕對(duì)值最大主彎矩為主彎矩M3極值的負(fù)數(shù)，該極值深入煤體距離為Lsd，即距離AD邊或BC邊的長度為Lsd，設(shè)該主彎矩極值的絕對(duì)值為Msd，那么滿足Msd=|M3|(a+Lsd，0)|，Msd=|M3|(-a-Lsd，0)|?？梢?，彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)的主彎矩極值主要在開采區(qū)中部、長邊彈-塑性煤體區(qū)及短邊彈-塑性煤體區(qū)，分析這幾個(gè)特征區(qū)域的主彎矩Msz,Msc與Msd的大小及所在位置Lsc與Lsd即可分析彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)的初次破斷順序、破斷位置及整體破斷規(guī)律。

圖5 彈-塑性基礎(chǔ)邊界基本頂主彎矩極值位置Fig.5 Location of the extreme value of the main roof main bending moment with the boundary of elastic-plastic foundation

彈-塑性基礎(chǔ)邊界條件下，基本頂板結(jié)構(gòu)的初次破斷規(guī)律主要受基本頂?shù)暮穸萮、彈性模量E、跨度L(L越大代表基本頂強(qiáng)度越大，可懸頂距離越大)、實(shí)體煤區(qū)彈性煤體基礎(chǔ)系數(shù)kt、煤體塑化范圍as,bs及淺部煤體基礎(chǔ)系數(shù)k0影響，本文采用控制變量法著重研究這些因素對(duì)彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)初次破斷特征的影響，最后研究這些因素的權(quán)重關(guān)系。

由于采用控制變量法進(jìn)行研究，以下所采用的計(jì)算基礎(chǔ)參數(shù)未經(jīng)說明更改的均為:工作面長度及跨度分別為132,44 m;基本頂?shù)膹椥阅Ａ縀、厚度h、泊松比ν及載荷q分別為32 GPa,6.5 m,0.22及0.32 MPa;淺部煤體的基礎(chǔ)系數(shù)k0及未塑化彈性煤體基礎(chǔ)系數(shù)kt分別為0及1.5 GN/m3，周邊煤體塑性區(qū)寬度as及bs均為4 m。

4.2 破斷特征的彈性煤體基礎(chǔ)系數(shù)kt效應(yīng)

如圖6所示，未塑化的彈性煤體基礎(chǔ)系數(shù)kt不僅可以改變基本頂?shù)钠茢辔恢们铱筛淖兓卷數(shù)钠茢囗樞颉?/p>

圖6 主彎矩及位置隨kt變化規(guī)律Fig.6 Variation rule of main bending moment and position with kt

(1)基本頂破斷位置方面:隨kt值增大，長邊主彎矩Msc深入煤體距離Lsc由大于煤體塑性區(qū)寬度bs過渡到逐步小于bs，這說明kt較小時(shí)長邊破斷線在未塑化的彈性煤體區(qū)上覆，kt較大時(shí)長邊破斷線在塑化煤體區(qū)上覆;kt值越小，長邊破斷線深入煤體距離Lsc越大，即煤體越軟，約束基本頂整體下沉變形的能力越弱，破斷線越超前煤壁，此時(shí)與固支邊界模型的差距也越大;短邊破斷位置Lsd與長邊Lsc的變化規(guī)律相同。

(2)基本頂破斷順序及形態(tài)方面:隨kt值增大，基本頂中部主彎矩Msz減小，而長邊與短邊的主彎矩Msc與Msd逐漸增大，但增長幅度小于降低幅度;①kt較小時(shí)，Msz>Msc>Msd，主彎矩達(dá)到彎矩極限Ms時(shí)，基本頂破斷順序?yàn)?中部下表面→長邊深入彈性煤體區(qū)上表面→短邊深入彈性煤體區(qū)上表面，最終在彈性煤體區(qū)上方形成“O”型斷裂圈，整體破斷形態(tài)為“O-X”型，如圖7(d),(e)所示;②kt較大時(shí)，Msc>Msz>Msd，基本頂破斷順序?yàn)?長邊深入塑化煤體區(qū)上表面→中部下表面→短邊深入塑化煤體區(qū)上表面，最終在塑化煤體區(qū)上覆形成“O”型斷裂圈，整體破斷形態(tài)為“O-X”型，如圖7(a)～(c)所示;③kt=kts時(shí)，基本頂“O”型斷裂圈在煤體彈-塑性分界線區(qū)域，如圖7(f),(g)所示;④kt=ktj時(shí)，Msz=Msc=Mj，即存在基本頂中部與長邊同時(shí)破斷的情況。

圖7 彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)破斷特征Fig.7 Fracture characteristic diagram of main roof plate structure with elastic-plastic foundation boundary

而固支邊模型得到基本頂斷裂線沿著煤壁與kt值無關(guān)，可見固支邊界模型所得結(jié)論與實(shí)際差距較大，尤其是破斷位置及破斷順序方面。

4.3 破斷特征的基本頂厚度效應(yīng)

如圖8所示，基本頂?shù)暮穸萮不僅可以改變基本頂?shù)钠茢辔恢们铱筛淖兓卷數(shù)钠茢囗樞颉?/p>

(1)基本頂破斷位置方面:隨h值增大，長邊主彎矩Msc深入煤體距離Lsc由小于煤體塑性區(qū)寬度bs過渡到大于bs，這說明h較大時(shí)長邊破斷線在未塑化的彈性煤體區(qū)上覆，h較小時(shí)長邊破斷線在塑化煤體區(qū)上覆;h值越大，長邊破斷線深入煤體距離Lsc越大，即基本頂越厚，需要更大范圍內(nèi)的煤體約束基本頂才能發(fā)生破斷，此時(shí)與固支邊界模型的差距也越大;短邊破斷位置Lsd與長邊Lsc的變化規(guī)律相同。

(2)基本頂破斷順序及形態(tài)方面:隨h值減小，基本頂中部主彎矩Msz減小，而長邊與短邊的主彎矩Msc與Msd逐漸增大;①h較小時(shí)，Msc>Msz>Msd，基本頂破斷順序?yàn)?長邊深入塑化煤體區(qū)上表面→中部下表面→短邊深入塑化煤體區(qū)上表面，最終在塑化煤體區(qū)上覆形成“O”型斷裂圈，整體為“O-X”型破斷形態(tài)，如圖7(a)～(c)所示;②h較大時(shí)，Msz>Msc>Msd，基本頂破斷順序?yàn)?中部下表面→長邊深入彈性煤體區(qū)上表面→短邊深入彈性煤體區(qū)上表面，最終在彈性煤體區(qū)上方形成“O”型斷裂圈，整體破斷形態(tài)為“O-X”型，如圖7(d),(e)所示;③h=hs時(shí)，基本頂“O”型斷裂圈在煤體彈-塑性分界線區(qū)域，如圖7(f),(g)所示;④h=hj時(shí)，Msz=Msc=Mj，即存在基本頂中部與長邊同時(shí)破斷的情況?；卷攺椥阅Ａ縀對(duì)破斷規(guī)律的影響與h相似。

圖8 主彎矩及位置隨基本頂厚度變化規(guī)律Fig.8 Variation rule of main bending moment and position with thickness of main roof

而固支邊模型得到基本頂斷裂線沿著煤壁，與h及E值無關(guān)，可見固支邊界模型所得結(jié)論與實(shí)際差距較大，尤其是破斷位置及破斷順序方面。

4.4 破斷特征的跨度效應(yīng)

圖2中，開采懸頂跨度越大，即代表基本頂?shù)膹?qiáng)度越大，那么研究跨度L對(duì)彈-塑性基礎(chǔ)邊界基本頂破斷規(guī)律的影響就是間接研究基本頂強(qiáng)度對(duì)破斷規(guī)律的影響。由圖9可得:

(1)基本頂破斷位置方面:隨L值增大，長邊主彎矩Msc深入煤體距離Lsc由大于煤體塑性區(qū)寬度bs過渡到逐步小于bs，這說明L較小時(shí)長邊破斷線在未塑化的彈性煤體區(qū)上覆，L較大時(shí)長邊破斷線在塑化煤體區(qū)上覆;L值越小，長邊破斷線深入煤體距離Lsc越大，此時(shí)與固支邊界模型的差距也越大;短邊破斷位置Lsd與長邊Lsc的變化規(guī)律相同。

(2)基本頂破斷順序方面:①L較大時(shí)，Msc>Msz>Msd，基本頂破斷順序?yàn)?長邊深入塑化煤體區(qū)上表面→中部下表面→短邊深入塑化煤體區(qū)上表面，最終在塑化煤體區(qū)上方形成“O”型斷裂圈，整體破斷形態(tài)為“O-X”型，如圖7(a)～(c)所示;②L較小時(shí)，Msz>Msc>Msd，基本頂破斷順序?yàn)?中部下表面→長邊深入彈性煤體區(qū)上表面→短邊深入彈性煤體區(qū)上表面，最終在彈性煤體區(qū)上方形成“O”型斷裂圈，整體破斷形態(tài)為“O-X”型，如圖7(d),(e)所示;③L=Ls時(shí)，基本頂“O”型斷裂圈在煤體彈-塑性分界線區(qū)域，如圖7(f),(g)所示;④L=Lj時(shí)，Msz=Msc=Mj，即存在基本頂中部與長邊同時(shí)破斷的情況。

而固支邊模型得到基本頂斷裂線沿著煤壁，與L值無關(guān)，可見固支邊界模型所得結(jié)論與實(shí)際差距較大，尤其是破斷位置及破斷順序方面。

4.5 破斷特征的煤體塑化范圍b0效應(yīng)

如圖10所示，煤體的塑化范圍b0不僅可以顯著改變基本頂?shù)钠茢辔恢们铱梢愿淖兓卷數(shù)钠茢囗樞颉?/p>

圖9 主彎矩及位置隨跨度L變化規(guī)律Fig.9 Variation rule of main bending moment and position with the span L

圖10 主彎矩及位置隨b0變化規(guī)律Fig.10 Variation rule of main bending moment and position with b0

(1)基本頂破斷位置方面:隨煤體塑化范圍b0值增大，長邊主彎矩Msc深入煤體距離增大;當(dāng)b0=b0s時(shí)，Lsc=L0sc=b0s，即煤體塑化范圍與基本頂長邊斷裂線深入煤體距離相等，可見，b0較小時(shí)長邊破斷線在彈性煤體區(qū)上覆，b0較大時(shí)長邊破斷線在塑化煤體區(qū)上覆;b0越大，長邊破斷線深入煤體距離Lsc越大，即煤體塑化范圍越大，約束基本頂下沉變形的能力越弱，破斷線越深入煤壁，此時(shí)與固支邊界模型的差距也越大;短邊破斷位置Lsd與長邊Lsc的變化規(guī)律相同。

(2)基本頂破斷順序方面:煤體塑化范圍b0增大時(shí)，各個(gè)主彎矩均增大，這相當(dāng)于煤體支撐基本頂?shù)哪芰p弱(或說懸頂程度增大)，所以各個(gè)主彎矩均增大;① 塑化范圍b0較小時(shí)，滿足Msc>Msz>Msd，基本頂破斷順序?yàn)?長邊深入彈化煤體區(qū)上表面→中部下表面→短邊深入彈性煤體區(qū)上表面，最終在彈性煤體區(qū)上方形成“O”型斷裂圈，整體破斷形態(tài)為“O-X”型，如圖7(d),(e)所示;② 塑化范圍b0較大時(shí)，滿足Msz>Msc>Msd，基本頂破斷順序?yàn)?中部下表面→長邊深入塑化煤體區(qū)上表面→短邊深入塑化煤體區(qū)上表面，最終在塑化煤體區(qū)上覆形成“O”型斷裂圈，整體破斷形態(tài)為“O-X”型，如圖7(a)～(c)所示;③b0=b0s=L0sc時(shí)，“O”型斷裂圈在煤體彈-塑性分界線區(qū)域，如圖7(f),(g)所示;④b0=b0j時(shí)，Msz=Msc=Mj，即存在基本頂中部與長邊同時(shí)破斷的情況。

而固支邊模型與彈性基礎(chǔ)邊界模型均無法研究煤體塑化范圍對(duì)基本頂板結(jié)構(gòu)破斷位置及規(guī)律的影響，可見傳統(tǒng)模型所得結(jié)論有局限性。

4.6 破斷特征的煤體塑化程度效應(yīng)

由圖11可知煤體的塑化程度(采用淺部塑化煤體基礎(chǔ)系數(shù)k0來表征)不僅可以顯著改變基本頂?shù)钠茢辔恢们铱梢愿淖兓卷數(shù)钠茢囗樞?

(1)基本頂破斷位置方面:隨煤體塑化程度增大(即k0減小)，長邊主彎矩Msc深入煤體距離Lsc由小于煤體塑化寬度bs過渡到逐步大于bs，這說明k0較小時(shí)長邊破斷線在未塑化的彈性煤體區(qū)上覆;k0較大時(shí)長邊破斷線在塑化煤體區(qū)上覆;k0值越小，長邊破斷線深入煤體距離Lsc越大，即煤體塑化程度越大，約束基本頂下沉變形的能力越弱，破斷線越深入煤體，此時(shí)與固支邊界模型的差距也越大;短邊破斷位置Lsd與長邊Lsc的變化規(guī)律相同。

(2)基本頂破斷順序方面:塑化程度增大時(shí)(即k0減小)，各個(gè)主彎矩均增大，且長邊Msc增長程度大于中部Msz的增長幅度，這相當(dāng)于淺部煤體支撐基本頂?shù)哪芰p弱，基本頂懸頂程度相對(duì)增大，所以各個(gè)主彎矩均增大;① 塑化程度較大時(shí)，滿足Msc>Msz>Msd，基本頂破斷順序?yàn)?長邊深入彈性煤體區(qū)上表面→中部下表面→短邊深入彈性煤體區(qū)上表面，最終在彈性煤體區(qū)上覆形成“O”型斷裂圈，整體破斷形態(tài)為“O-X”型，如圖7(d),(e)所示;② 塑化程度較小時(shí)，滿足Msz>Msc>Msd，基本頂破斷順序?yàn)?中部下表面→長邊深入塑化煤體區(qū)上表面→短邊深入塑化煤體區(qū)上表面，最終在塑化煤體區(qū)上覆形成“O”型斷裂圈，整體為“O-X”型破斷形態(tài)，如圖7(a)～(c)所示;③k=k0s時(shí)，基本頂“O”型斷裂圈在煤體彈-塑性分界線區(qū)域，如圖7(f),(g)所示;④k0=k0j時(shí)，Msz=Msc=Mj，即存在中部與長邊同時(shí)破斷的情況。

圖11 主彎矩及位置隨k0變化規(guī)律Fig.11 Variation rule of main bending moment and position with k0

而固支邊模型與彈性基礎(chǔ)邊界模型均無法研究煤體塑化程度對(duì)基本頂板結(jié)構(gòu)破斷位置及規(guī)律的影響，可見傳統(tǒng)模型所得結(jié)論有局限性。

4.7 破斷特征的影響因素之間存在的關(guān)系(比值不變規(guī)律)

上文采用控制變量法分析了基本頂?shù)暮穸萮、彈性模量E、跨度L、實(shí)體煤區(qū)未塑化的彈性煤體基礎(chǔ)系數(shù)kt、煤體塑化范圍及煤體塑化程度k0對(duì)彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)破斷特征及規(guī)律的影響。下面分析影響因素之間存在的關(guān)系。

如圖12所示，同時(shí)改變kt與h的大小，且滿足kt與k0取任意比值不變(圖中兩者比值為10:1)，kt與h3也取任意比值不變，可以得到彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)各個(gè)主彎矩大小及位置均不改變，破斷順序和位置也不變。

圖12 主彎矩及位置的影響因素權(quán)重關(guān)系Fig.12 Weight relation of main bending moment and position

可見，雖然彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)破斷的影響因素復(fù)雜，但是k0,kt與h3之間滿足上述“比值不變”條件時(shí)，具有破斷規(guī)律不變的特征。

4.8 與彈性基礎(chǔ)邊界模型破斷規(guī)律對(duì)比分析

彈性基礎(chǔ)邊界模型與彈-塑性基礎(chǔ)邊界模型都考慮了煤體可變形的實(shí)際特征，但是考慮煤體塑化程度和塑化范圍的模型得到傳統(tǒng)模型得不到的諸多結(jié)論，更符合實(shí)際且對(duì)實(shí)踐的指導(dǎo)意義顯著，下面為2類模型的對(duì)比說明。

(1)破斷位置方面。2類模型都可以得到基本頂深入煤體斷裂，但是彈性基礎(chǔ)邊界模型只能得到斷裂線在彈性煤體區(qū)，而本文模型得到基本頂?shù)臄嗔丫€位置主要有3種：① 彈性煤體區(qū);② 塑性煤體區(qū);③ 煤體的彈塑性分界線。本文模型更加全面且符合實(shí)際，對(duì)工程實(shí)踐的指導(dǎo)意義更大。

(2)破斷順序方面?；卷?shù)某醮纹茢辔恢枚加?類:① 長邊超前煤壁區(qū)的基本頂上表面;② 基本頂懸頂區(qū)中部的下表面，但煤體的塑化范圍和塑化程度也顯著影響基本頂?shù)钠茢辔恢煤推茢囗樞?，而彈性基礎(chǔ)邊界模型無法得到實(shí)際工程中煤體必然塑化引起基本頂破斷順序改變的結(jié)論。

(3)破斷形態(tài)方面。基本頂巖板整體破斷形態(tài)都是“O-X”型，而本文模型得到基本頂O型圈位置有3類:① 在彈性煤體區(qū)上覆;② 在塑性煤體區(qū)上覆;③ 在煤體的彈塑性分界線上覆，而彈性基礎(chǔ)邊界模型只能得到O型圈在煤體的彈性區(qū)上覆，顯然本文模型得到的結(jié)論更加全面且符合實(shí)際。

5 工程算例分析

研究清楚符合實(shí)際圍巖關(guān)系的基本頂板結(jié)構(gòu)的破斷順序、位置及整體形態(tài)特征，對(duì)于預(yù)防頂板大面積切頂來壓事故及合理確定工作面?zhèn)确较锏赖奈恢玫戎陵P(guān)重要。

基于彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)模型，下面選取2個(gè)方面，即從工作面推進(jìn)方向基本頂斷裂位置和工作面短邊區(qū)域基本頂斷裂位置角度說明本文力學(xué)模型的重要工程意義。

現(xiàn)給出一具體算例進(jìn)行說明，圖2中，一工作面的跨度為42 m、工作面長度為200 m，基本頂厚度為10 m、彈性模量為35 GPa，彈性煤體基礎(chǔ)系數(shù)為1.2 GN/m3，周邊塑性區(qū)煤體寬度為5 m，淺部煤體基礎(chǔ)系數(shù)為0(即淺部煤體基本無承載能力)，塑性區(qū)煤體由淺部到深部的基礎(chǔ)系數(shù)按線性增長，基于本文的計(jì)算方法和破斷準(zhǔn)則，得到如下基本結(jié)論，并與傳統(tǒng)模型進(jìn)行對(duì)比分析。

5.1 工作面推進(jìn)方向基本頂板結(jié)構(gòu)破斷位置的工程意義

本文的力學(xué)模型得到，基本頂超前煤壁的斷裂位置Lsc為6.6 m(與實(shí)際值接近)，斷裂線處在煤體的彈性區(qū)(圖13(b))，而傳統(tǒng)的固支邊界模型得到基本頂?shù)臄嗔丫€沿著長邊煤壁(圖13(a))，顯然在基本斷裂位置方面?zhèn)鹘y(tǒng)模型不能反應(yīng)工程實(shí)際。

圖13 工作面長邊區(qū)基本頂斷裂位置對(duì)比Fig.13 Comparison of main roof fracture positions at the long side area of working face

工作面推進(jìn)方向上，基本頂是超前煤壁破斷的，那么破斷時(shí)，前方有彈-塑性煤體的支撐，斷裂時(shí)工作面并不會(huì)顯著來壓，當(dāng)支架推進(jìn)到斷裂線區(qū)域時(shí)，壓力顯著增大，這說明，工作面來壓是個(gè)過程且破斷與顯著來壓之間有時(shí)間差，那么基本頂剛斷裂時(shí)即預(yù)警，提前采取措施加強(qiáng)支架控制，防止工作面大面積切頂事故的發(fā)生;而傳統(tǒng)模型得到，基本頂?shù)臄嗔丫€沿著煤壁，顯然不能從斷裂位置角度得到有助于實(shí)踐的結(jié)論，且不能有效指導(dǎo)提前采取措施預(yù)警工作面大面積切頂事故的發(fā)生。

同時(shí)由上文得到的煤體塑化程度和范圍越大時(shí)，開采區(qū)域周邊基本頂?shù)闹鲝澗卦酱螅抑鲝澗厣钊朊后w的距離越大的結(jié)論明確表明了在實(shí)踐過程中，采取對(duì)煤體的弱化措施不僅可以加速基本頂?shù)钠茢?減小來壓步距)，而且破斷線更加深入煤體，基本頂剛破斷時(shí)可以減弱對(duì)支架的沖擊，這對(duì)工作面頂板的安全控制有重要指導(dǎo)意義。

5.2 工作面?zhèn)认?短邊區(qū))基本頂板結(jié)構(gòu)破斷位置的工程意義

本文的力學(xué)模型得到，工作面?zhèn)认蚧卷斏钊朊后w的斷裂位置Lsd為6.9 m(與實(shí)際值接近)，處在煤體的彈性區(qū)(圖14(b))，而傳統(tǒng)的固支邊界模型得到基本頂?shù)臄嗔丫€沿著短邊煤壁(圖14(a))，這說明傳統(tǒng)模型在基本斷裂位置方面不能有效反應(yīng)工程實(shí)際。

基本頂側(cè)向斷裂線的位置對(duì)沿空巷道位置選擇意義重大，直接決定了沿空巷道頂板結(jié)構(gòu)穩(wěn)定性程度。傳統(tǒng)的固支邊界模型無法得到短邊區(qū)域基本頂?shù)臄嗔盐恢迷谏钊朊后w區(qū)域;彈性基礎(chǔ)邊界板結(jié)構(gòu)模型雖然可以得到基本頂?shù)臄嗔盐恢迷谏钊朊后w處，但是卻得不到基本頂斷裂線在煤體的彈性區(qū)還是塑性區(qū)或者彈-塑性分界面。

圖14 工作面短邊區(qū)基本頂斷裂位置對(duì)比Fig.14 Comparison of main roof fracture positions in short side of working face

沿空巷道位置選擇不僅與基本頂在側(cè)方的斷裂線位置有關(guān)，還與基本頂斷裂線在煤體的彈性區(qū)還是塑性區(qū)有關(guān)，如圖14(b)所示，若在采空區(qū)鄰側(cè)掘進(jìn)下一區(qū)段回采巷道，煤柱寬度Lm與基本頂在側(cè)方的斷裂位置Lsd的大小關(guān)系(有3種:Lm>Lsd,Lm

可見，本文構(gòu)建的彈-塑性基礎(chǔ)邊界模型可以有效彌補(bǔ)傳統(tǒng)固支邊界與彈性基礎(chǔ)邊界板結(jié)構(gòu)模型的缺陷，不僅可以有效得到基本頂深入煤體的斷裂位置，而且可以得到斷裂位置在彈性煤體區(qū)還是塑性煤體區(qū)，這對(duì)于巷道位置選擇，煤柱寬度留設(shè)等均有重要指導(dǎo)意義。

6 結(jié) 論

(1)彈-塑性基礎(chǔ)邊界基本頂板結(jié)構(gòu)的破斷位置、順序及形態(tài)，均隨煤體的塑化范圍和塑化程度、彈性煤體的基礎(chǔ)系數(shù)、懸頂跨度、基本頂?shù)暮穸群蛷椥阅Ａ慷淖儭?/p>

(2)煤體塑化程度和塑化范圍增大，基本頂各個(gè)區(qū)域的主彎矩絕對(duì)值均增大，即基本頂更易發(fā)生破斷，且基本頂在懸頂區(qū)周邊破斷位置深入煤體距離增大，這對(duì)于工程實(shí)踐中控制基本頂斷裂線深入煤體位置和來壓步距具有指導(dǎo)意義。

(3)煤體塑化范圍增大，中部主彎矩增長程度大于長邊的，即塑化范圍大時(shí)，基本頂中部區(qū)域更易先破斷;煤體塑化程度增大，基本頂長邊主彎矩增長程度大于中部的，即煤體塑化程度大時(shí)，基本頂長邊更易先破斷。

(4)煤體塑化范圍及淺部塑化煤體基礎(chǔ)系數(shù)較大時(shí)，基本頂?shù)钠茢辔恢眉绊樞驗(yàn)?開采懸頂區(qū)中部→長邊深入塑化煤體區(qū)→短邊深入塑化煤體區(qū)，在塑化煤體區(qū)上覆形成“O”型斷裂圈;反之，基本頂?shù)钠茢辔恢眉绊樞驗(yàn)?長邊深入彈性煤體區(qū)→開采懸頂區(qū)中部→短邊深入彈性煤體區(qū)，在彈性煤體區(qū)上覆形成“O”型斷裂圈，最終均形成“O-X”型破斷形態(tài)。

(5)淺部塑化煤體基礎(chǔ)系數(shù)k0與彈性煤體基礎(chǔ)系數(shù)kt為任意比值不變且kt與h3為任意比值不變時(shí)，基本頂?shù)闹鲝澗卮笮〖拔恢镁蛔?，基本頂?shù)钠茢嘁?guī)律不變。

本文的力學(xué)模型，彌補(bǔ)了傳統(tǒng)的固支邊界板結(jié)構(gòu)模型和彈性基礎(chǔ)邊界板結(jié)構(gòu)模型的缺陷。