第4章北翼上向進(jìn)路護(hù)壁承載層穩(wěn)定性數(shù)值模擬

白象山鐵礦床主礦體賦存在閃長巖與砂頁巖接觸帶部位。其形態(tài)主要受礦區(qū)

背斜構(gòu)造控制，橫向呈平緩拱形。礦體走向長度1780m,傾向以4B線為界以南

向南傾，傾角一般25°左右；4B線以北向北傾，傾角一般5°?25°,總體以13°,

左右呈波狀向北傾伏。礦體賦存標(biāo)高-200?-600m。平均厚度34.41m,最大厚度

可達(dá)121.72m,水平厚度最大達(dá)300m。

水平厚度最厚處在礦體的西翼Y=81800-Y=82200（-400m標(biāo)高上下），水平

最大厚度達(dá)300?400m。其他部分厚度變化在10?120多米的范圍內(nèi)。礦體上部

收縮連成一片，拱頂部分也相對(duì)較厚，垂高最大達(dá)96m（2#勘探線）。

礦區(qū)地層主要有三疊系.匕統(tǒng)黃馬青組，侏羅系中下統(tǒng)象山群，白堊系下統(tǒng)上

火山巖組，以及第四系沖坡積層。礦體的直接頂?shù)装鍑鷰r主要為閃長巖、角頁巖，

其次為砂頁巖。礦體頂?shù)装褰V圍巖均不同程度地存在一個(gè)不穩(wěn)定帶，其厚度一

般為10?20m。

礦區(qū)內(nèi)斷裂構(gòu)造較為發(fā)育，勘探共發(fā)現(xiàn)10條斷層。可分為成礦前斷層及成

礦后斷層。其中地下開采影響較大的有3條斷層，分別為Fi、F2和F3。

目前，針對(duì)礦山所采取的機(jī)械化盤區(qū)點(diǎn)柱式上向分層充填采礦法，在不鋪設(shè)

人工假底的情況下對(duì)頂?shù)字M(jìn)行回采，對(duì)護(hù)壁厚度及回收方式對(duì)圍巖進(jìn)行穩(wěn)定性

進(jìn)行數(shù)值模擬。本研究將利用FLAC3D分析軟件，模擬不同回采率、護(hù)壁厚度

條件，承載層及相鄰充填體的應(yīng)力、應(yīng)變、塑性變形等情況，進(jìn)而評(píng)價(jià)目前采場(chǎng)

的穩(wěn)定性，并確定最合理的結(jié)構(gòu)參數(shù)。

4.1FLAC3D軟件簡介

FLAC3D采用的“顯式拉格朗日”算法和“混合-離散分區(qū)”技術(shù)，能夠非常

準(zhǔn)確的模擬材料的塑性破壞和流動(dòng)，由于無須形成剛度矩陣，因此，基于較小內(nèi)

存空間就能夠求解大范圍的三維巖土工程問題。由于采用了自動(dòng)慣量和自動(dòng)阻尼

系數(shù)，克服了顯式公式存在的小時(shí)間步長的限制以及阻尼問題。所以，F(xiàn)LAC3D

是一個(gè)求解三維巖土問題的最理想工具之一。FLAC3D有良好的前處理功能，計(jì)

算時(shí)程序自動(dòng)將模型剖分成六面體單元，每個(gè)單元都可以有自己的材料模型，材

料可以在外力及應(yīng)力場(chǎng)的作用下發(fā)生屈服流動(dòng)，節(jié)點(diǎn)的位置也會(huì)隨著模型的屈服

流動(dòng)而發(fā)生改變（大變形時(shí)），因此,FLAC3D在模擬大變形問題上有其獨(dú)到之處。

FLAC3D包含了三維網(wǎng)格自動(dòng)生成器，通過匹配、連接，由網(wǎng)格生成器生成所需

要的網(wǎng)格。生成器還能夠自動(dòng)產(chǎn)生交叉結(jié)構(gòu)區(qū)域網(wǎng)格（如相交巷道）。另外，為方

便用戶劃分網(wǎng)格，F(xiàn)LAC3D提供了12種原始模型（primitivemesh）,使得用戶能

根據(jù)工程實(shí)際更為快捷方便的建立理想的計(jì)算模型。FLAC3D的后處理功能也非

常強(qiáng)大，用戶可以根據(jù)需要通過相應(yīng)的操作命令打印或者繪制出自己所需的數(shù)據(jù)

或者圖形。FLAC3D近具有記錄追蹤功能（history）,使得用戶可以很方便的得到

所需參數(shù)在計(jì)算過程中的發(fā)展曲線或者得到多個(gè)參數(shù)之間的相互關(guān)系曲線。針對(duì)

不同的材料，F(xiàn)LAC3D軟件提供多達(dá)10種材料模型，能更真實(shí)地模擬實(shí)際材料

的力學(xué)行為。同時(shí)，F(xiàn)LAC3D中包含靜力、動(dòng)力、蠕變、滲流、溫度等5種計(jì)算

模式，并且還可以進(jìn)行多模式的耦合分析,FLAC3D可以模擬多種巖十工程地質(zhì)

不連續(xù)面，包括斷層、節(jié)理等以及常見的多種支護(hù)形式，例如梁（beam）、錨桿

（索）（cable）、樁（pile）殼（shell）等。FLAC3D還包含一個(gè)強(qiáng)有力的程序語言?F1SH,

能夠使用戶定義新的變量和函數(shù)。為特定需要的用戶提供了一個(gè)專門工具。

4.2計(jì)算模型構(gòu)建

4.2.1基本假定

鑒于井下開采技術(shù)條件復(fù)雜多變，為便于建模和分析計(jì)算，需作出如下假設(shè):

（1）開挖礦體的水平厚度和傾角為固定值；

（2）礦巖體假設(shè)為理想彈塑性體，在屈服點(diǎn)以后，隨著塑性流動(dòng)，材料強(qiáng)

度和體積無改變；

（3）礦體和圍巖為局部均質(zhì)、各向同性的材料，塑性流動(dòng)不改變材料各向

同性；

（4）考慮到巖石的脆性，分析中涉及到的所有物理量均與時(shí)間無關(guān)；

（5）不考慮應(yīng)變硬化（或軟化）；

（6）模型側(cè)面限制水平位移，底面限制垂直位移；

（7）考慮到有限元程序的局限性，假設(shè)場(chǎng)地內(nèi)無構(gòu)造活動(dòng)的影響，原巖地

應(yīng)力為大地靜力場(chǎng)型，各巖層之間為整合接觸，巖層內(nèi)部為連續(xù)介質(zhì)，模型中不

考慮地下水活動(dòng)的影響，也不考慮巖層和礦體中的結(jié)構(gòu)面、裂隙和軟弱層的存在

與影響；

(8)計(jì)算選定的載荷將不隨單元方向變化而改變，始終保持它們的最初的

方向，表面載荷作用在變形單元表面法向，且可被用來模擬“跟隨”力，大小就是

上覆巖層或充填體層的重量。

4.2.2材料力學(xué)參數(shù)確定

根據(jù)礦山實(shí)際情況以及地質(zhì)報(bào)告，對(duì)巖石及巖體進(jìn)行力學(xué)參數(shù)分析。并對(duì)充

填體進(jìn)行室內(nèi)試驗(yàn)得到彈性模量、抗壓強(qiáng)度、泊松比、容重、粘結(jié)力、內(nèi)摩擦角

等礦巖體力學(xué)參數(shù)，其抗拉強(qiáng)度的獲得參考業(yè)內(nèi)經(jīng)驗(yàn)m約是其單軸抗壓強(qiáng)度而

的15%?19%,從而得出計(jì)算礦巖體力學(xué)參數(shù)如表4-1。

表4-1相關(guān)礦巖基本力學(xué)參數(shù)

彈性模抗壓強(qiáng)度許用應(yīng)力內(nèi)摩

泊松容重丫粘結(jié)力

類別量EjalGt擦角

比P(kN-m-3)(MPa)

(GPa)(MPa)(MPa)a(°)

上盤56403.40.2224.899.349.8

礦體59606.430.234.999.5448.5

下盤32151.20.2222.343.2649.6

普通

0.211.730.210.1518.130.47也7

充填體

人工假

0.240.390.2218.051.2348.

底

4.2.3模型建立

頂?shù)字厥罩饕扇∩舷蜻M(jìn)路充填法為主的采礦方法組合方案，即：礦巖穩(wěn)

固性差、涌水量大的地段采用普通進(jìn)路充填法(進(jìn)路規(guī)格3?4mx3?4m)；穩(wěn)固性

稍好、涌水量中等地段，采用預(yù)控頂上向進(jìn)路充填采礦法(進(jìn)路寬度6m、進(jìn)路

高度6-7m)；穩(wěn)固性較好、涌水量較小地段則采用小分段空?qǐng)鏊煤蟪涮畈色E法

FLAC3D5.00

>52012ItascaConsultingGroup.Inc

Step276507

2017/10/2214:21:06

ZGroup

GroupSlot:Any

Ictt

1Defaiit

1hb

wl

Iw2

Z

YX

c)75%回收率模型正視圖

FLAC3D5.00

€2012ItascaConsultingGroi4).Inc

d)75%回收率回收過程模型正視圖

圖4-1數(shù)值模型圖

4.3模擬結(jié)果及分析

本次數(shù)值模擬分析主要研究進(jìn)路采場(chǎng)護(hù)壁及相鄰充填體的應(yīng)力分布規(guī)律、位

移分布規(guī)律和最佳的承載層厚度及進(jìn)路斷面尺寸。數(shù)值模擬主要結(jié)果數(shù)據(jù)列于表

4-2,應(yīng)力、變形情況見圖4-2至37,由結(jié)果數(shù)據(jù)及圖表可知：

表4-2數(shù)值模擬結(jié)果主要數(shù)據(jù)表

結(jié)果數(shù)據(jù)

承載層相鄰礦柱

模型

Z方向位最大拉應(yīng)最大壓應(yīng)最大拉應(yīng)最大壓應(yīng)

移/mm力/MPa力/MPa力/MPa力/MPa

50%-0.971.040.110.831.98

50%過程-1.201.240.180.63.24

75%-5.625.891.661.753.74

75%過程-5.015.961.731.483.49

上表為模擬結(jié)果數(shù)據(jù)，下圖為拉應(yīng)力結(jié)果圖:

FLAC3D5.00

?2012ItascaConsultingGroupme.

Step220988

2017/1(^229:47:00

ContourofXX-Stress

CaicuiaMWVoUneencAmm

1.6096E*06

1.5000E^06

1.2500E*06

1OOOOE*06

7.5000E*05

5.0000E*05

25000E>05

O.OOOOE?OO

-2.5000E*05

-5.0000E*05

-7.5000E*05

-1.0000E*06

-1.2500E*06

-1.5000E406

-1.7500E*06

-2.0000E*06

-2.2500E*06

-2.5000E*06

-2.7500E*06

-29749E*06

圖4-250%回收率拉應(yīng)力正視圖

FLAC3D5.00

必12ItascaConsultingGroupme

Step220988

2017/10/229:47:13

C()n(ourofXX-Stress

FIneon

CakutetAdbyVoUnecnc卸ng

1I1.6102E*06

1.5000E*06

12500E>06

1.0000E*06

/.ouwt*ue>

50000E>05

25000E*05

O.OOOOE*OO

-2.5000E*05

-5.0000E*05

-7.5000E*05

-1.0000E+06

?1.2500E+06

-1.5000E*06

-17500E*06

-20000E*06

-2.2500E*06

-2.5000E*06

-27500E*06

-2.9737E*06

圖4-350%回收率拉應(yīng)力切片圖」

FLAC3D5.00

?2012ItascaConsultingGroupme.

Step220988

2017/10/229:47:44

ContourofXX-Stress

PUne:on

CalcuMdbyVaUrwfinc

9.1852E>05

90000E*05

80000E*05

7.0000E*05

6.0000E*05

50000E*05

4.0000E*05

3.0000E*05

20000E*05

1.0000E*05

O.OOOOE*OO

-1.0000E?05

?2,0000E25

-3.0000E*05

-4.0000E*05

-5.0000E*05

60000E25

-7.0000E*05

30000E25

-9.0000E*05

-9.1018E*05

圖4-450%回收率拉應(yīng)力切片圖_2

FLAC3D5.00

?2012ItascaConsultingGroupme.

Step208627

2017/1(^2213:50:48

ContourofXX-Stross

CaicuiaMWVoiumeencAmm

I1.7877E*06

1.7500E>06

1.5000E*06

12500E>06

1.0000E*06

7.5OOOE*O5

50000E-05

2.5000E*05

O.OCOOE*00

-2.5000E*05

-5.0000E*05

-7.5000E*05

-1.0000E+06

-1.2500E*06

-1.5000E*06

-1.7500E*06

-2.0000E*06

-2.2500E*06

-2.5000E*06

?27243E+06

圖4-550%回收率回采過程拉應(yīng)力切片正視圖

FLAC3D5.00

€^012ItascaConsult)“GroupInc

Step208627

2017/10/2213:49:46

ContourofXX-S(ress

Planeon

CakMcdbyVaUncfiricAvar9^ng

1.7877E*06

L7500E.06

1.5000E*06

1.2500E*06

1OOOOE*Q6

7.5000E-05

5.0000E*05

25000E*05

O.OOOOE*OO

-2.5000E*05

-5.0000E-05

-7.5000E*05

-1.0000E*06

-1.2500E406

-1.5000E*06

-1.7500E*06

-2.0000E406

-2.2500E*06

-2.5000E*06

-27239E*06

圖4-650%回收率回采過程拉應(yīng)力切片圖」

FLAC3D5.00

?2012ItascaConsultingGroupme.

Step208627

2017/1(W213:50:29

ContourofXX-Stress

PUne:on

CateuWWdbyVaUrwfinc

I1.56XE-06

1.4000E*06

12000E*06

1.0000E*06

8.0000E*05

60000E*05

?40000E.05

'2.0000E*05

O0COOE*0O

.-2.0000E*05

I-4.0000E*05

■COOOOE25

|77175E25

圖4-750%回收率回采過程拉應(yīng)力切片圖_2

FLAC3D5.00

€2012ItascaConsultingGroup.Inc

Step276507

2017/10/2213:42:12

ContourofXX-Str?K?

CalcuiaMWVokmtncAvera0ng

6.0488E*06

6.0000E*06

5.0000E*06

40000E*06

3.0000E*06

2.0000E*06

1OOOOE*06

O.OOOOE?OO

-1.0000E*06

.2,0000E*06

-3.0000E*06

-4.0000E-K)6

-5.0000E*06

-60000E*06

-7.0000E*06

-8.0000E+06

-90000E*06

-1.0000E*07Z

-1.1000E*07

-11859E*O7

LX

圖4-875%回收率拉應(yīng)力正視圖

FLAC3D5.00

物＞12ItascaConsultingGroup.Inc

Step276507

2017/10/2213:42:37__________

ContourofXX-Stress

Planeon

CafcuWwlbyMoSemcAveQgng

I6.0472E—06

16.0000E*06

150000E*06

■4.0000E*06

13.0000E*06

20000E+06

1.0000E*06

O.OOOOE?OO

?-1.0000E*06

-2.0000E*06

|-3.0000E*06

-4.0000E*06

-50000E*06

60000E+06

i-7.0000E+06

-80000E*06

-9.0000E*06z

-1.0000E+07

-1.1000E*07

-1.1859E*07Ju

圖4-975%回收率拉應(yīng)力切片圖」

FLAC3D5.00

€2012ItascaConsultingGroup.Inc

Step276507

2017/10/2213:43:17

ContourofXX-Stress

Planeon

ColcukaMVoUrwmcAvoraging

3.1125E+06

3.0000E*06

25000E*06

2.0000E*06

1.5000E*06

1OOOOE*Q6

I5.0000E*05

0.0000E*00

-50000E*05

-1.0000E*06

-1.5000E*06

-2.0000E*06

-2.5000E*06

-3.0000E*06

-3.5000E+06

-40000E*06

45000E26Y

-4.7004E+06

圖4-1075%回收率拉應(yīng)力切片圖_2

FLAC3D5.00

?2012ItascaConsultingGroup.Inc.

Step259900

2017/10/2213:5805

ContourofXX-Stress

CafcutatedbyVoUnrtncAveraoinQ

I60983E*06

6.0000E—06

5.0000E*06

40000E*06

3.0000E*06

2.0000E+06

10O00E+06

O.OOOOE?OO

-1.0000E+06

-2.0000E*06

1-3.0000E*06

■4.0000E-K)6

-5.0000E*06

60000E+06

-7.0000E*06

-8.0000E+06

-O.OOOOE^Oe

-1.0000E*07

-1.1000E*07

-1J463E*O7

圖4-1175%回收率回采過程拉應(yīng)力正視圖

FLAC3D5.00

€C012ItascaConsultingGroup.Inc

Step259900

2017/10/2213:58:37

ContourofXX-Stress

Planeon

CakukatodbyVoUnolricAveraging

I6.0977E*06

I60000E*06

I5.0000E*06

■4.0000E*06

30000E*06

2.0000E—06

1.0000E*06

OOOOOE?OO

-1.0000E*06

20000E26

1-3,0000E*06

-4.0000E*06

-5.0000E*06

60000E+06

-7.0000E*06

-8.0000E+06

-9.0000E406

-1.0000E*07

-1.1000E*07

-1.1462E*07L

圖4-1275%回收率回來過程拉應(yīng)力切片圖」

FLAC3D5.00

?2012ItascaConsultingGroup.Inc

Step259900

2017/10/2213:59:02

ContourofXX-Stress

Planeon

ColcukaMbyVoUrwmcAvAraging

I6.0662E*06

I6.0000E*06

■50000E*06

4.0000E*06

3.0000E+06

20000E+06

1.0000E*06

O.OOOOE?OO

-1.0000E*06

-2.0000E*06

-3.0000E*06

-4.0000E?06

-5.0000E*06

60000E26

-6.1274E+O6

Y

圖4-1375%回收率回采過程拉應(yīng)力切片圖_2

根據(jù)上圖關(guān)于拉應(yīng)力的結(jié)果圖，可看出，50%回收率時(shí)其拉應(yīng)力要小于75%

回收率的拉應(yīng)力，假底層的存在有效的降低了護(hù)壁層拉應(yīng)力的大小。根據(jù)75%回

收率的回收順序，先“隔一采一”也即是50%回收率的過程，充填后再進(jìn)行二步回

收，再進(jìn)行二步充填。由于二步充填中的“相鄰礦柱''變成一步回采過程中的充

填體，這使得75%回收率時(shí)總體拉應(yīng)力大于50%回收率。綜合以上結(jié)果圖可知

5%回收率時(shí)總體拉應(yīng)力大于5()%回收率，跨中拉應(yīng)力總體大于相鄰其他區(qū)域拉

應(yīng)力，回采過程中由于空區(qū)的存在使得拉應(yīng)力略有增加。

FLAC3D5.00

?2012ItascaConsultingGroupme.

Step220988

2017/KV229:40:59

ContourofZZ-Stress

CalcutetedbyVoune<nc但電由

2.3447E*05

O.OOOOE?OO

-2.5000E*05

-5.0000E*05

-7,5000E*05

-1.0000E*06

I-12500E*06

-1.5000E*06

-1.7500E*06

-2.0000E*06

-2.2500E406

I呢

?2.411+06

圖4-1450%回收率壓應(yīng)力正視圖

FLAC3D5.00

?2012ItascaConsultingGroupme.

Step220988

2017/10/229:46:05

ContourofZZ-Stress

PUne:on

CafcuWWdbyVaUrwfinc

2.3442E>05

O.QOOOE*00

-2.5000E*05

I-5.0000E25

-7.5000E*05

-10000E*06

-1.2500E*06

1-1.5000E*06

-1.7500E*06

-2.0000E406

-2.2500E*06

-2.412OE?O6

圖4-I550%回收率壓應(yīng)力切片圖」

FLAC3D5.00

?2012ItascaConsultingGroupme.

Step220988

2017/1(W29:46:36

ContourofZZ-Stress

Plane:on

CefcuMtodbyVoUrwfinc

1.4972E>05

0.00O0E*0O

?2,0000E25

I-4.0000E*05

60000E*05

30000E25

I-1.0000E<06

■-1.2000E*06

1-1.4000E*06

.-1,6000E406

I-1.8000E*06

I-2.0000E*06

|Q1786E26

圖4-1650%回收率壓應(yīng)力切片圖_2

FLAC3D5.00

?2012ItascaConsultingGroupme.

Step208627

2017/10/2213:51:46

ContourofZZ-Stress

CafcutetedbyVoUne<ncAve*—"

2.4163E*O5

O.OOOOE?OO

-2.5000E*05

-5.0000E*05

-7.5000E*05

-1.0000E*06

?12500E"

-1.5000E<06

-1.7500E*06

-2.0000E*06

-2.2500E*06

-2.5000E*06

-2.7500E*06

-3.0000E*06

-3.2500E*06

-3.4295E*06

圖4-1750%回收率回采過程壓應(yīng)力正視圖

FLAC3D5.00

?2012ItascaConsultingGroupme.

Step208627

2017/1(W213:52那

ContourofZZ-Stress

PUne:on

CafcuWWdbyVaUrwfinc

i2.4162E>05

O.QOOOE*00

-2.5000E*05

-5.0000E25

?7.5000E25

..10000E*06

-1.2500E*06

-1.5000E*06

-1.7500E*06

-2.0000E406

-2.2500E*06

-2.5000E?06

-27500E*06

-3.0000E*06

-3.2500E*06

-3.4293E*06

圖4-1850%回收率回采過程壓應(yīng)力切片圖」

FLAC3D5.00

?2012ItascaConsultingGroupme.

Step208627

2017/10/2213:54:07

ContourofZZ-Stress

Planeon

ColcuMMbyVoUrwfinc

I-1.2753E*O6

-1.4000E*06

-16000E*06

-1.8000E*06

-2.0000E+06

-22000E*06

-2.4000E*06

-2.6000E*06

-2.8000E*06

I---333...023009001007EEE***000666

圖4-1950%回收率回采過程壓應(yīng)力切片圖_2

FLAC3D5.00

€2012ItascaConsultingGroup.Inc

Step276507

2017/10/2213:4447

ContourofZZ-Stress

CafcutatedOyVoUneincAver紿ng

6.4517E*O4

O.OOOOE?OO

-2.5000E+05

-50000E*05

-7.5000E*05

-1.0000E*06

.12$00E*O6

-1.5000E*06

-1.7500E*06

.2,0000E*06

-2.2500E*06

-2.5000E-KJ6

.2.7500E*06

.3,0000E*06

-3.2500E-K)6

-3.5000E+06

-37435E*06

圖4-2075%回收率壓應(yīng)力正視圖

FLAC3D5.00

€2012ItascaConsultingGroup.Inc.

Step276507

2017/10/2213:4422

ContourorZZ-Stress

Planeon

CafcMRbyVoUWicAMA附OQ

6.4370EQ

O.OOOOE*(X)

25000EM5

-5.0000E*05

I-7.5000E+05

-1.0000E*06

-1.2500E*06

-1.5000E*06

-1.7500E*06

.2.0000E*06

-2.2500E*06

-2.5000E?06

-2.7500E406

-3.0000E*06

-3.2500E+06

-3.5000E*06

?3.7436E26

圖4-2I75%回收率壓應(yīng)力切片圖」

FLAC3D5.00

€2012ItascaConsultingGroup.Inc

Step276507

2017/10/2213:43:43

ContourorZZ-Stress

Planeon

ColcukaMbyVoljnwlncAveracjng

3.6975E*04

O.OOOOE*(X)

?25000EM5

-5.0000E*05

-7.5000E+05

-1.0000E*06

-1.2500E*06

-1.5000E*06

-1.7500E*06

-2.0000E*06

-2.2500E*06

-25000E*06

.27500E*06

-3.0000E*06

-3.2500E+06

-3.5000E*06

-3.5711E*06

圖4-2275%回收率壓應(yīng)力切片圖_2

FLAC3D5.00

€2012ItascaConsultingGroup.Inc

Step259900

2017/10/2214:00:27

ContourofZZ-Stress

Cafcutated?VoUneincAver紿ng

34192E*O4

O.OOOOE?OO

-2.5000E+05

-50000E*05

-7.5000E*05

-1.0000E*06

.12$00E*O6

-1.5000E*06

-1.7500E*06

.2,0000E*06

-2.2500E*06

-2.5000E-KJ6

.2.7500E*06

.3,0000E*06

-3.2500E-K)6

-3.5000E+06

-37216E*O6

圖4-2375%回收率回采過程壓應(yīng)力正視圖

FLAC3D5.00

€2012ItascaConsultingGroup.Inc.

Step259900

2017/10/2214:00:04

ContourorZZ-Stress

Planeon

ColcukaMt)yVoUrwlncAvoracjng

3.4103EQ

O.OOOOE*(X)

25000EM5

-5.0000E*05

-7.5000E+05

-1.0000E*06

-1.2500E*06

-1.5000E*06

I-1.7500E*06

.2.0000E*06

-2.2500E*06

-2.5000E?06

-2.7500E406

-3.0000E*06

-3.2500E+06

-3.5000E*06

-3.7219E*O6L

圖4-2475%回收率回采過程壓應(yīng)力切片圖」

FLAC3D5.00

€2012ItascaConsultingGroup.Inc

Step259900

2017/10/2213:5940

ContourofZZ-Stress

Planeon

CofcgtedbyVoljnwlncAveracjng

-3.1889E*05

-4.0000E+05

60000EM5

-8.0000E*05

-1.0000E*06

-12000E*06

-1.4000E*06

-1,eOOOE*06

-18000E*06

-2.0000E*06

-2.2000E*06

-2.4000E*06

-26000E*06

1-2.6684E*O6

圖4-2575%回收率回采過程壓應(yīng)力切片圖_2

根據(jù)壓應(yīng)力圖可知，承載層壓應(yīng)力小于相鄰f柱壓應(yīng)力，并且75%回收率時(shí)

總體壓應(yīng)力大于50%回收率的壓應(yīng)力?；夭蛇^程中的壓應(yīng)力由于采空區(qū)的存在,

使得其壓應(yīng)力略大于充填后的壓應(yīng)力。壓應(yīng)力總體大小均未達(dá)到巖石破壞條件,

根據(jù)模擬結(jié)果，兩種回收率都不會(huì)引起由壓應(yīng)力主導(dǎo)而發(fā)生破壞。

FLAC3D5.00

€2012ItascaConsultingGroup.Inc.

Step220988

2017/10/229:38:26

ConOto.OuOrOOOEf?ZOO-Displacenient

2.208707J0940EU-90Q23

2.2000E-02

24.0000E-032

11.8000E-02

1.6000E.02

61..04000E-023

1.2000E-02

1.0000E.02

I8.0000E-03

圖4-2650%回收率位移正視圖

FLAC3D5.00

>32012ItascaConsultingGroup.Inc

Step220988

2017/10/229:39:03___________

ContourOfZ-Displacement

Planeon

I2.2755E-02

I2.2000E-02

I2.0000E-02

/1.8000E-02

1.6000E-02

I1.4000E92

?12000E-02

[1.0000E-02

18.0000E-03

6OOOOE-03

.4.0000E-03

I2.0000E-03

IOOOOOE-OO

1-2.5795E-O4

LE

圖4-2750%回收率位移切片圖」

FLAC3D5.00

€2012ItascaConsultingGroup.Inc

Step220988

2017/10/229:40:18

ContourOfZ-Displacement

Plane:on

2.1292E-03

2.1000E-03

2.0000E-03

1.9000E-03

1.8000E-03

1.7000E03

16000E-03

1.5000E-03

1.4000E-03

13000E.03

1.2000E-03

1.1000E93

1OOOOE-03

9.0000E-04

8.0000EQ

7.0000EQ

62265E.04

圖4-2850%回收率位移切片圖_2

FLAC3D5.00

€2012ItascaConsultingGroup.Inc.

Step208627

2017/10/2213:55:12

ContourOfZ-Displacenient

22975E92

2.2000E-02

_2.0000E-02

Ia1.8000E-02

1.6000E.02

?1.4000E-02

J1.2000E-02

11.0000E.02

8.0000E-03

.6.0000E-03

I40000E-03

■2.0000E-03

|O.OOOOE?OO

圖4-2950%回收率回采過程位移正視圖

FLAC3D5.00

>32012ItascaConsultingGroup.Inc

Step208627

2017/10/2213:56:54__________

ContourOfZ-Displacement

Planeon

I2.2975E-02

I2.2000E-02

B2.0000E-02

1.8000E-02

1.6000E-02

■1.4000E92

I12000E-02

11.0000E-02

8.0000E-03

.60000E-03

I4.0000E-03

I2.0000E-03

IOOOOOE-OO

z

一

XB

圖4-3050%回收率回采過程位移切片圖」

FLAC3D5.00

€2012ItascaConsultingGrovp.Inc.

Step208627

2017/10/2213:5637

ContourOfZ-Displacement

Plane:on

1.6812E-03

1.6500E-03

1.6000E93

15500E.03

1.5000E-03

1.4500E-03

1.4000E-03

1.3500E-03

1.3000E-03

12500E-03

1.2000E-03

1.1500E-03

1.1000E.03

1.0500E-03

1.0000E-03

9.5000EQ

9.0000E-04

8.5000E-04

8.0000EQ

75000E-04

7.0572E-04

圖4-3150%回收率回采過程位移切片圖_2

FLAC3D5.00

€2012ItascaConsulti