固體充填開(kāi)采頂板多跨超靜定結(jié)構(gòu)分析與充填支架控制作用研究

劉建功，趙家巍

(1.河北工程大學(xué),河北 邯鄲 056038; 2.中國(guó)礦業(yè)大學(xué)(北京) 力學(xué)與建筑工程學(xué)院，北京 100083)

在我國(guó)當(dāng)前生態(tài)文明建設(shè)的大背景下，煤炭工業(yè)向綠色發(fā)展轉(zhuǎn)型的步伐逐漸加大[1]，全國(guó)主要產(chǎn)煤地區(qū)陸續(xù)出臺(tái)了有關(guān)煤礦矸石地面排放的限制政策，為保證礦區(qū)的可持續(xù)發(fā)展，同時(shí)為解決我國(guó)巨量的“三下”壓煤開(kāi)采與地表變形控制問(wèn)題[2-3]，矸石固體充填開(kāi)采的需求日益增加。

當(dāng)前充填開(kāi)采工藝和裝備取得了長(zhǎng)足的發(fā)展[4-7]，但關(guān)于巖層的運(yùn)動(dòng)與控制理論還尚需完善。關(guān)于充填開(kāi)采巖層運(yùn)動(dòng)問(wèn)題，文獻(xiàn)[8-12]應(yīng)用彈性地基梁理論分析了充填體對(duì)巖層變形的控制作用;文獻(xiàn)[13]通過(guò)FLAC3D內(nèi)置Fish語(yǔ)言編制了相應(yīng)于矸石壓實(shí)實(shí)驗(yàn)測(cè)試結(jié)果的充填體非線性壓實(shí)程序，用數(shù)值模型模擬了矸石充填采空區(qū)后上覆巖層移動(dòng)規(guī)律;文獻(xiàn)[14]通過(guò)相似物理模擬試驗(yàn)?zāi)M了不同充填率情況下上覆巖層應(yīng)力變化和巖層移動(dòng)特征;文獻(xiàn)[15]通過(guò)分析充填開(kāi)采支架與圍巖關(guān)系和上覆巖層移動(dòng)特征，確立了頂板載荷估算方法，并在給定直接頂懸距的條件下探討了圍巖與支架的關(guān)系;文獻(xiàn)[16]根據(jù)充填采煤液壓支架結(jié)構(gòu)原理及控頂作用，建立了支架頂梁的力學(xué)模型，分析得出了頂梁受力情況及立柱受力之間的關(guān)系，并運(yùn)用Pro/E軟件對(duì)支架進(jìn)行三維建模和運(yùn)動(dòng)學(xué)仿真分析，得出了各主要部件在工作過(guò)程中的運(yùn)動(dòng)特征曲線。上述研究對(duì)推動(dòng)充填開(kāi)采巖層運(yùn)動(dòng)與控制理論的發(fā)展起到了積極的作用，但對(duì)頂板、支架、充填體三者相互作用關(guān)系還有待更深層次的探討，尤其是缺少在頂板載荷、支架支撐力、充填體壓實(shí)穩(wěn)定長(zhǎng)度等諸多參數(shù)未知條件下有效分析三者之間互相作用、互相影響量化關(guān)系的方法。

充填開(kāi)采上覆巖層應(yīng)力場(chǎng)分布特征對(duì)充填支架工作狀態(tài)和充填體力學(xué)行為有著較大的影響，而充填支架和充填體又反作用于頂板應(yīng)力分布，充填支架承受的載荷不僅與原巖地質(zhì)與開(kāi)采條件有關(guān)，還與充填體力學(xué)性能密切相關(guān)，因此，充填開(kāi)采巖層應(yīng)力場(chǎng)的重新分布過(guò)程是頂板、支架、充填體三者相互作用、相互影響的結(jié)果，若單方面給定某一變量的數(shù)值則失去了各因素之間的互動(dòng)性。文獻(xiàn)[17-18]在大量的工程觀測(cè)基礎(chǔ)上，提出了“連續(xù)曲形梁”結(jié)構(gòu)模型(圖1)，該模型闡述了連續(xù)曲形梁的形成過(guò)程及其幾何與力學(xué)特征，由于結(jié)構(gòu)的連續(xù)而保持了力的傳遞，從而減輕了對(duì)原巖應(yīng)力場(chǎng)的擾動(dòng)，從本質(zhì)上揭示了充填開(kāi)采巖層運(yùn)動(dòng)形式。本文進(jìn)一步對(duì)頂板、支架、充填體三者相互作用機(jī)理進(jìn)行深入研究，為充填開(kāi)采工程中充填支架和充填體關(guān)鍵力學(xué)參數(shù)的選取與確定提供一個(gè)合理的分析和研究方法。

圖1 充填開(kāi)采上覆巖層“連續(xù)曲形梁”結(jié)構(gòu)形態(tài)Fig.1 Structure of “continuous curved beam” of overlying strata in filling mining

1 頂板-支架-充填體相互作用力學(xué)分析

1.1 力學(xué)模型及其假設(shè)條件

對(duì)于垮落法開(kāi)采的巖層運(yùn)動(dòng)，長(zhǎng)期的實(shí)踐觀測(cè)和研究表明，隨工作面的推進(jìn)，采空區(qū)上方堅(jiān)硬巖層斷裂成排列整齊的巖塊，巖塊間受水平推力作用而形成鉸接關(guān)系，“砌體梁”結(jié)構(gòu)模型由此而提出[19]，為我國(guó)煤礦開(kāi)采礦山壓力與巖層控制奠定了理論基礎(chǔ)(圖2)。

圖2 “砌體梁”結(jié)構(gòu)受力模型[19]Fig.2 “Voussoir Beam” structure model[19]

圖2中q1～q8為頂板載荷，r1為煤壁支承反力，r2～r8為已冒落矸石支承反力。假設(shè)圖2模型中第一巖塊為懸露巖塊，則此結(jié)構(gòu)的自由度ω[19]為

ω=3n-2(n-1)-(n+2)=0(1)

其中，n為巖塊的塊數(shù);n-1為鉸鏈數(shù);n+2為鏈桿數(shù)。由此可知，此結(jié)構(gòu)為靜定結(jié)構(gòu)。

本文在“砌體梁”結(jié)構(gòu)基礎(chǔ)上，根據(jù)充填開(kāi)采時(shí)堅(jiān)硬巖層僅發(fā)生微小變形、結(jié)構(gòu)依然連續(xù)而不再斷裂這一現(xiàn)象，建立“頂板-支架-充填體”相互作用力學(xué)模型(圖3)，從結(jié)構(gòu)力學(xué)上看，其自由度ω<0，為多跨超靜定連續(xù)結(jié)構(gòu)。在該力學(xué)模型中，煤壁和充填體作為主要承載體，由于充填開(kāi)采對(duì)巖層變形有要求，不允許頂板出現(xiàn)較大的充前下沉，因此充填支架應(yīng)對(duì)頂板起限定變形作用。

圖3 充填開(kāi)采“頂板-支架-充填體”相互作用力學(xué)模型Fig.3 “Roof-support-filling body” interaction mechanics model of filling mining

“頂板-支架-充填體”相互作用力學(xué)模型基于如下假設(shè)條件而提出:

(1)以基本頂為分析對(duì)象，承受上覆巖層非均布載荷，將直接頂簡(jiǎn)化為均布載荷。

(2)基本頂伸入煤壁端為固定支承，充填體支撐端因存在豎向位移和水平約束而假設(shè)為定向支承，在充填體壓實(shí)穩(wěn)定前還承受轉(zhuǎn)動(dòng)約束。

(3)將工作面前方煤壁和后方充填體假設(shè)為彈性地基支承，兩者彈性地基系數(shù)不同。

(4)為便于結(jié)構(gòu)力學(xué)分析和計(jì)算，將分布彈性地基支承離散為若干跨等效彈性支承。

1.2 多跨超靜定連續(xù)結(jié)構(gòu)力學(xué)分析

將兩端約束以未知力和彎矩表示，支座支承以反力Ri代替，得出“頂板-支架-充填體”相互作用力學(xué)模型的計(jì)算簡(jiǎn)圖，如圖4所示。

圖4 力學(xué)模型計(jì)算簡(jiǎn)圖Fig.4 Calculation sketch of above mechanical model

圖4中R1,R2為煤壁等效彈性支承反力；R3，R4為充填支架支承反力;R5～Rn為充填體等效彈性支承反力;lm，lz，lc分別為煤壁支承區(qū)、支架控頂區(qū)以及充填體控頂區(qū)長(zhǎng)度;MO，F(xiàn)OX，F(xiàn)OY分別為煤壁端頂板彎矩、水平力、垂直支承力，MC，F(xiàn)CX分別為充填體端頂板彎矩、水平力。

煤壁和充填體等效彈性剛度按式(2)計(jì)算:

(2)

式中，Kmi和Kci分別為煤壁和充填體的等效彈性剛度;kmi和kci分別為煤壁和充填體的彈性地基系數(shù);lmi和lci分別為煤壁和充填體的等效彈性支承間距;b為力學(xué)模型的單位厚度。

取計(jì)算簡(jiǎn)圖4中的某一跨度為li作為隔離體，其受力分析如圖5所示。

圖5 隔離單跨受力Fig.5 Stress analysis of an isolated single span structure

根據(jù)支座i兩側(cè)剪力Fsi和彎矩Mi平衡條件，可得

Fsi+Ri=q(x)+Fsi+1-Ri+1(3)

Fi-Fi+1=0(4)

式中l(wèi)ix為至第i跨左端的距離。

煤壁和充填體等效彈性支承反力Ri與其壓縮量Di的關(guān)系為

(6)

當(dāng)i=3,4時(shí)，Ri為充填支架在限定變形狀態(tài)Di下的恒定支撐作用力，D4為給出值。

對(duì)于該超靜定連續(xù)結(jié)構(gòu)問(wèn)題的求解，因超靜定次數(shù)為n+2，共n+5個(gè)未知約束，為得到問(wèn)題的解還需引入變形連續(xù)(變形協(xié)調(diào))條件。將連續(xù)結(jié)構(gòu)沿第i個(gè)支承處的截面切開(kāi)，去除3個(gè)內(nèi)部多余約束，并代以三對(duì)多余未知力X1,X2,X3，每一對(duì)未知力分別作用在截面切口兩側(cè)，截面兩側(cè)沿3個(gè)未知力方向的位移為相對(duì)水平位移Δ1、相對(duì)豎向位移Δ2和相對(duì)轉(zhuǎn)角Δ3，則第i個(gè)截面切口處的形變分量為

(r=1,2,3;s=1,2,3)(7)

式中，X1(i)=Fi，X2(i)=Fsi，X3(i)=Mi;δrs(i)為單位力作用時(shí)在截面切口i處產(chǎn)生的位移;ΔrRi為荷載Ri單獨(dú)作用時(shí)產(chǎn)生的3個(gè)未知力方向的位移。

式(7)中所有單位力產(chǎn)生的位移系數(shù)和自由項(xiàng)均可由變形位能式(8)求出:

(8)

為保證截面兩側(cè)是連續(xù)的而不發(fā)生任何相對(duì)位移，則要求基本結(jié)構(gòu)在多余未知力X1，X2，X3及外載荷共同作用下的位移形態(tài)應(yīng)與原結(jié)構(gòu)的完全相同，于是具體的位移條件為Δr(i)=0，將式(7)展開(kāi)得到第i個(gè)截面切口形變分量的矩陣形式

(9)

式中，ΔHi,ΔVi,ΔTi分別為垂直力、水平力和彎矩左截面處引起的形變分量。

將式(9)中的矩陣寫(xiě)成如下形式

式中，ΔHiRi，ΔViRi,ΔTiRi分別為支座反力Ri單獨(dú)作用時(shí)產(chǎn)生的3個(gè)未知力方向的形變分量。

則第1～n個(gè)截面切口形變分量的矩陣表達(dá)式為

這樣就由式(10)得到3n-2(n-1)個(gè)變形協(xié)調(diào)方程，與前面3個(gè)平衡方程構(gòu)成n+5個(gè)求解方程，對(duì)應(yīng)n+5個(gè)未知約束，所求問(wèn)題全部得解，即當(dāng)確定頂板載荷q(x)后可以解出方程中所有支座的反力、梁體內(nèi)力、梁端彎矩等未知量。

對(duì)于q(x)的計(jì)算，文獻(xiàn)[19]在進(jìn)行巖層受力分析時(shí)將頂板載荷視為線性分布載荷，文獻(xiàn)[20]根據(jù)組合梁理論推導(dǎo)得出了頂板載荷的式(11)，計(jì)算得出的為均布載荷，常用來(lái)進(jìn)行巖梁極限跨距的求解。

(11)

其中，E1，E2，…，En為各層巖層的彈性模量;n為巖層數(shù);h1，h2，…，hn為各層巖層的厚度;γ1，γ2，…，γn為各層巖層的容重。當(dāng)計(jì)算到(qn+1)1<(qn)1時(shí)，則以(qn)1作為作用于第1層巖層的單位面積上的載荷。

采用上述2種方法估算頂板載荷對(duì)某些問(wèn)題的求解均具有一定的適用性，但在處理充填支架的受力分析問(wèn)題時(shí)與實(shí)際差距較大。用力學(xué)模型求解采動(dòng)支承應(yīng)力比較困難，而數(shù)值模型則不能較好地反映充填支架的受力狀態(tài)，為了更真實(shí)地反映頂板載荷的分布，更合理地分析頂板、支架和充填體的相互作用關(guān)系，本文嘗試將力學(xué)模型與數(shù)值模型相結(jié)合，將數(shù)值模擬得出的覆巖應(yīng)力曲線作為力學(xué)模型的頂板載荷q(x)，為該類(lèi)問(wèn)題的分析提供一個(gè)新的解決思路。

2 覆巖應(yīng)力分布數(shù)值模擬分析

2.1 數(shù)值模型建立

建立有限元平面應(yīng)變數(shù)值模型(圖6)，模型尺寸長(zhǎng)×高為500 m×300 m，上邊界施加5 MPa地層載荷，下邊界固定，兩側(cè)邊界設(shè)置水平約束，模型自上而下概化為較軟巖層、關(guān)鍵層、軟弱巖層、基本頂、直接頂、煤層、底板等7個(gè)巖層，煤層走向開(kāi)挖長(zhǎng)度200 m，充填體地基系數(shù)15 MN/m3，充填支架限定變形量150 mm，巖石本構(gòu)模型采用M-C模型，煤巖體物理與力學(xué)參數(shù)見(jiàn)表1。

圖6 開(kāi)挖與充填數(shù)值模擬計(jì)算模型Fig.6 Numerical model of filling and mining

巖層名稱(chēng)巖性厚度/m彈性模量/GPa泊松比黏聚力/MPa內(nèi)摩擦角/(°)抗拉強(qiáng)度/MPa抗剪切強(qiáng)度/MPa較軟巖層砂質(zhì)頁(yè)巖603.50.304.5328.215.0關(guān)鍵層砂巖6020.00.257.52715.023.0軟弱巖層泥巖602.50.353.2342.45.5基本頂細(xì)砂巖129.50.286.2299.816.0直接頂粉砂巖55.00.265.5267.513.0煤層煤43.00.332.3283.66.7底板粉砂巖1006.00.324.6308.414.0

2.2 充填開(kāi)采應(yīng)力場(chǎng)分布

工作面在距左側(cè)邊界350 m處開(kāi)切眼，自右向左開(kāi)采、充填，推進(jìn)至30，100，200 m時(shí)的Von Mises應(yīng)力場(chǎng)如圖7所示，通過(guò)圖7可以直觀地判斷出應(yīng)力擾動(dòng)程度和范圍以及充填采動(dòng)后擾動(dòng)應(yīng)力場(chǎng)的演化過(guò)程。

圖8 垮落法開(kāi)采覆巖Von Mises應(yīng)力場(chǎng)分布云圖Fig.8 Von Mises stress field distribution cloud map with mining method

從圖8可以看出，由于充填體的支撐作用，隨工作面向前推進(jìn)，工作面和開(kāi)切眼造成的應(yīng)力擾動(dòng)區(qū)被限制在很小的范圍內(nèi)，工作面后方充填區(qū)域的覆巖應(yīng)力場(chǎng)很快恢復(fù)平衡，對(duì)比垮落法開(kāi)采(圖7)，充填開(kāi)采顯著減弱了對(duì)覆巖應(yīng)力場(chǎng)的擾動(dòng)。將直接頂、基本頂下邊界和關(guān)鍵層中部的垂向應(yīng)力值提取出來(lái)，得到如圖9所示的充填開(kāi)采覆巖垂直應(yīng)力場(chǎng)的分布曲線，基本頂下邊界垂向應(yīng)力即為直接頂承受的載荷。

圖9 充填開(kāi)采覆巖垂直應(yīng)力場(chǎng)分布曲線Fig.9 Vertical stress field distribution of overlying strata

工作面推進(jìn)長(zhǎng)度200 m，改變充填體地基系數(shù)，其值依次為5，10，20，25 MN/m3時(shí)的基本頂下邊界垂直應(yīng)力場(chǎng)分布曲線如圖10所示。

圖10 不同充填體地基系數(shù)頂板應(yīng)力分布Fig.10 Stress distribution on the roof with different foundation coefficient of filling bodies

3 頂板-支架-充填體力學(xué)模型算例分析

由數(shù)值模型計(jì)算得出的直接頂載荷與直接頂自重之和即為所求q(x)，將煤壁前方峰值點(diǎn)至充填體穩(wěn)定壓實(shí)區(qū)這一段的應(yīng)力曲線提取出來(lái)，增加直接頂自重引起的均布載荷0.125 MPa，對(duì)數(shù)值模擬得到的計(jì)算應(yīng)力曲線進(jìn)行擬合，并將數(shù)值模型坐標(biāo)(190～210 m)與力學(xué)模型坐標(biāo)(0～60 m)對(duì)應(yīng)，得到圖11所示的應(yīng)力分布擬合曲線和式(12)所示的數(shù)學(xué)擬合公式。

(12)

將直接頂覆巖載荷q(x)代入前述力學(xué)模型中進(jìn)行計(jì)算，煤壁彈性地基系數(shù)取km=600 MN/m3，充填體彈性地基系數(shù)kc=15 MN/m3，頂板限定變形量150 mm，分布彈性地基離散跨間距5 m，求解力學(xué)模型中支座反力、梁體內(nèi)力、梁端彎矩等各項(xiàng)未知量，具體計(jì)算過(guò)程借助Matlab計(jì)算程序?qū)崿F(xiàn)，需指明的是計(jì)算得出的支座反力為等效集中力，最后還應(yīng)根據(jù)跨間距進(jìn)行線性分布載荷還原，以反映真實(shí)的受力情況。

圖11 頂板應(yīng)力擬合曲線Fig.11 Fitting curve of roof stress

采用同樣的計(jì)算過(guò)程和計(jì)算方法，依次計(jì)算出頂板限定變形量150 mm時(shí)分別對(duì)應(yīng)于充填體彈性地基系數(shù)為5，10，20，25 MN/m3以及充填體彈性地基系數(shù)15 MN/m3時(shí)對(duì)應(yīng)于頂板限定變形量分別為50，100，200，250 mm時(shí)的支座反力，得到不同q(x)條件下的“頂板-支架-充填體”三者相互作用力學(xué)關(guān)系(圖12)，因最大值與最小值的差值較大，圖12(a)采用對(duì)數(shù)坐標(biāo)顯示。

圖12 “頂板-支架-充填體”相互作用關(guān)系Fig.12 Interaction relationship of “roof-support-filling body”

通過(guò)改變充填體彈性地基系數(shù)(頂板限定變形量固定為150 mm)，計(jì)算結(jié)果表明(圖12(a))，當(dāng)充填體彈性地基系數(shù)由5 MN/m3增加至25 MN/m3時(shí)，充填支架前頂梁支護(hù)強(qiáng)度由1.27 MPa降低至0.54 MPa，后頂梁支護(hù)強(qiáng)度由1.34 MPa降低至0.61 MPa;當(dāng)充填體承載能力增強(qiáng)后，充填支架需要提供的支撐力減小，頂板同時(shí)轉(zhuǎn)移部分載荷到前方煤壁，意味著頂板傳遞力的作用有所增強(qiáng)。

通過(guò)改變頂板的限定變形量(充填體彈性地基系數(shù)固定為15 MN/m3)，計(jì)算結(jié)果表明(圖12(b))，當(dāng)限定變形量由250 mm減小至50 mm時(shí)，充填支架前頂梁需提供的支撐力由0.45 MPa提高至1.77 MPa，后頂梁需提供的支撐力由0.55 MPa提高至1.96 MPa;隨限定變形量減小，充填支架支撐力提高，而前方煤壁與后方充填體承受的載荷值有所回落。

因充填開(kāi)采對(duì)充填前頂板變形量有要求，不允許頂板產(chǎn)生較大的充前下沉，但根據(jù)圖12(b)的分析，當(dāng)對(duì)頂板充前下沉控制要求越高，充填支架需提供的支撐力更大，甚至超過(guò)當(dāng)前支架強(qiáng)度設(shè)計(jì)能力和裝備制造水平，因此應(yīng)將頂板的充前下沉控制在合理的范圍之內(nèi)，當(dāng)限定頂板變形100 mm時(shí)，改變充填體彈性地基系數(shù)，得出基本頂?shù)膬?nèi)力(剪力和彎矩)分布曲線，如圖13所示。

圖13 不同充填體彈性地基系數(shù)基本頂內(nèi)力分布Fig.13 Internal pressure distribution of basic roof with dif-ferent foundation coefficient of filling bodies

由式(13)計(jì)算基本頂破斷的臨界條件:

(13)

式中，F(xiàn)s為梁的截面剪力;Rs為梁的抗剪強(qiáng)度;A為截面面積;M為梁的彎矩;Rt為梁的抗拉強(qiáng)度;h為梁的厚度。

根據(jù)表1中的基本頂力學(xué)參數(shù)計(jì)算得出基本頂?shù)呐R界剪力條件為Fsmax=192 MN,臨界彎矩條件為Mmax=235.2 MN·m,則在上述條件下基本頂保持連續(xù)的參數(shù)為充填體彈性地基系數(shù)不小于15 MN/m3，根據(jù)該條件下充填支架與頂板相互作用關(guān)系(圖14)，充填支架前頂梁支護(hù)強(qiáng)度不小于0.85 MPa，后頂梁支護(hù)強(qiáng)度不小于0.92 MPa。

從上述分析來(lái)看，由于頂板彎曲下沉?xí)r后頂梁位移量大于前頂梁，因此后頂梁承受著更多的載荷，而這一結(jié)果與目前某些煤礦充填支架實(shí)際設(shè)計(jì)有所差異，尤其對(duì)于5 m以上大采高充填支架其后頂梁設(shè)計(jì)支護(hù)強(qiáng)度略低于前頂梁(例如邢東煤礦ZC5160/30/50型支架前頂梁支護(hù)強(qiáng)度0.80 MPa，后頂梁強(qiáng)度0.68 MPa)，后頂梁支護(hù)強(qiáng)度不足會(huì)造成更大的頂板充前下沉，同時(shí)造成前頂梁支護(hù)強(qiáng)度較大的富裕系數(shù)。因此根據(jù)本文的分析結(jié)果，充填支架后頂梁的設(shè)計(jì)支護(hù)強(qiáng)度應(yīng)大于前頂梁。充填支架型式如圖15所示。

以控制地表沉降和保護(hù)地表附著物為目標(biāo)的充填開(kāi)采，必須限制關(guān)鍵層變形，當(dāng)基本頂形成連續(xù)曲形梁以后，上覆巖層呈整體彎曲下沉變形狀態(tài)，關(guān)鍵層得到有效保護(hù)，從而避免地表大規(guī)模塌陷和建筑物的破壞，保護(hù)地下含水層和地面水系，實(shí)現(xiàn)煤炭生態(tài)保護(hù)性開(kāi)采。

4 結(jié) 論

(1)根據(jù)充填開(kāi)采工程實(shí)踐中巖層的實(shí)際變形形態(tài)，建立了基于多跨超靜定連續(xù)結(jié)構(gòu)的“頂板-支架-充填體”相互作用力學(xué)模型，將分布彈性地基支承離散為若干跨支承，從而根據(jù)支承基礎(chǔ)的反力求得充填支架和充填體對(duì)頂板的支撐力。

(2)數(shù)值模型的計(jì)算結(jié)果表明，由于充填體的支撐作用，隨工作面向前推進(jìn)，工作面和開(kāi)切眼造成的應(yīng)力擾動(dòng)區(qū)被限制在很小的范圍內(nèi)，工作面后方充填區(qū)域的覆巖應(yīng)力場(chǎng)很快恢復(fù)平衡，對(duì)比垮落法開(kāi)采，充填開(kāi)采顯著減弱了對(duì)覆巖應(yīng)力場(chǎng)的擾動(dòng)。將覆巖垂直應(yīng)力場(chǎng)中的基本頂下邊界垂向應(yīng)力曲線提取出來(lái)并進(jìn)行數(shù)學(xué)公式擬合，從而得出了力學(xué)模型計(jì)算所需的頂板載荷q(x)。

(3)通過(guò)算例分析得出了“頂板-支架-充填體”3者之間的量化關(guān)系，當(dāng)充填體承載能力增強(qiáng)后，充填支架需要提供的支撐力減小，頂板同時(shí)轉(zhuǎn)移部分載荷到前方煤壁，說(shuō)明頂板傳遞力的作用有所增強(qiáng);當(dāng)充填體彈性地基系數(shù)固定時(shí)，充填支架支撐力隨頂板限定變形量減小而提高，同時(shí)前方煤壁與后方充填體承受的載荷值有所回落。

(4)在文中給出的開(kāi)采與地質(zhì)條件下，當(dāng)限定頂板變形量為100 mm時(shí)，基本頂保持連續(xù)的充填參數(shù)為充填體彈性地基系數(shù)不小于15 MPa，充填支架前頂梁支護(hù)強(qiáng)度不小于0.85 MPa，后頂梁支護(hù)強(qiáng)度不小于0.92 MPa。