基于FLAC3D的煤巖體單軸壓縮數(shù)值模擬研究

褚佳琪

（安徽理工大學(xué) 礦業(yè)工程學(xué)院，安徽 淮南 232001）

0 引言

近年來，隨著我國煤炭資源開采深度和開采強(qiáng)度的增加，發(fā)生沖擊地壓的礦井?dāng)?shù)量快速增多，事故頻度和強(qiáng)度也明顯提高。由于我國礦山工程地質(zhì)災(zāi)害問題日漸突出，社會的生產(chǎn)活動以及生態(tài)環(huán)境均受到嚴(yán)重制約，而礦山工程地質(zhì)災(zāi)害受到人為活動條件、氣候植被條件、地形地貌條件、地質(zhì)條件及歷史災(zāi)害發(fā)生的規(guī)模、頻次、密度等許多因素的影響，近年來發(fā)生了不同規(guī)模大小的地質(zhì)災(zāi)害，其中煤巖動力災(zāi)害問題尤為突出。煤巖動力災(zāi)害的主要表現(xiàn)為地震、火山噴發(fā)、山崩滑坡、煤與瓦斯突出、沖擊地壓等，其中沖擊地壓是一種與地應(yīng)力有關(guān)的地質(zhì)災(zāi)害。據(jù)悉，礦山進(jìn)入1 000 m 以下進(jìn)行深部開采時，由于高地應(yīng)力、高地溫、高滲透壓以及強(qiáng)烈開采擾動的影響，煤巖體往往發(fā)生沖擊地壓。

沖擊地壓的發(fā)生，與煤巖體的力學(xué)特性、組合形式以及受力環(huán)境等各種因素密切相關(guān)，眾多學(xué)者針對煤巖組合體的特性進(jìn)行了研究探討，付斌等利用數(shù)值模擬軟件，對不同圍壓和不同組合傾角條件下的煤巖組合體進(jìn)行了力學(xué)特性和聲發(fā)射特征研究；胡順銀等對煤巖模型破壞時的應(yīng)力及變形情況進(jìn)行了模擬分析；張凱文等模擬分析了真三軸加載過程中巖石的損傷與破裂演化過程及特征；蔣中明等對不同圍壓下對巖石進(jìn)行了不同變形參數(shù)及強(qiáng)度參數(shù)的數(shù)值模擬研究；周天白等采用數(shù)值模擬反演了煤巖材料的參數(shù)，對煤巖材料變形破壞的主控參數(shù)進(jìn)行了研究；賈善坡等研究了層狀巖體單軸壓縮情況下的變形破壞特征及其影響其破損的規(guī)律；周元超等針對不同高度比的煤巖組合體，模擬研究了不同組合方式下煤巖組合體的力學(xué)特性和聲發(fā)射特征。

本文基于FLAC3D 的數(shù)值模擬研究，采用摩爾-庫倫力學(xué)模型，建立煤巖體力學(xué)數(shù)值計算模型，進(jìn)行單軸加載數(shù)值模擬，對比分析不同加載時步下（step200、step400、step600）的煤巖體應(yīng)力及位移演化規(guī)律，同時對煤巖組合體的應(yīng)力—應(yīng)變曲線進(jìn)行分析，從數(shù)值模擬的角度對煤巖單軸加載下的力學(xué)特性進(jìn)行認(rèn)識。

1 煤巖力學(xué)機(jī)理

礦山在進(jìn)行深部開采時，由于高地應(yīng)力的存在，使得煤巖體在一定的環(huán)境條件下往往發(fā)生沖擊地壓。煤巖體的變形是衡量深部礦山穩(wěn)定性的重要指標(biāo)，故探討煤巖體的力學(xué)機(jī)理有著重要的意義，若煤巖體的變形破壞較大，對深部礦山的危害也會相對較嚴(yán)重。在深部礦山中，煤巖體內(nèi)任一點的應(yīng)力狀態(tài)可由正應(yīng)力σ、σ、σ和剪應(yīng)力τ、τ、τ、τ、τ、τ組成的應(yīng)力分量表示，且可以用應(yīng)力張量σ表示：

式中，，=，，，且σ可用力學(xué)模型表示，如圖1所示。

圖1 煤巖體力學(xué)模型示意圖

由圖1 可知，煤巖體在深部礦山中處于三向受力的狀態(tài)，對于煤巖體所受到來自、、方向的正應(yīng)力與剪應(yīng)力，符合廣義的胡克定律，即煤巖體在受力之后，煤巖體的應(yīng)力與應(yīng)變之間成線性相關(guān)關(guān)系：同時，煤巖體在深部礦山中的應(yīng)力及應(yīng)變也可用張量分量表示：

式中，、為拉梅常數(shù)；煤巖體的力學(xué)本構(gòu)關(guān)系可通過拉梅常數(shù)表示：

2 模擬過程及分析

2.1 模型建立

本次數(shù)值模擬采用由美國ITASCA 公司開發(fā)的仿真計算軟件—FLAC3D 軟件，該軟件采用三維連續(xù)體的快速拉格朗日分析，同時其包含的11 種材料力學(xué)模型可以更準(zhǔn)確地對要研究的對象進(jìn)行數(shù)值模擬及力學(xué)分析。本次采用的力學(xué)數(shù)值模型是基于在進(jìn)行試驗時的實際試樣，運用FLAC3D 數(shù)值模擬軟件的導(dǎo)入命令模式，將編寫的模型命令導(dǎo)入FLAC3D。基于煤巖體的受力變形特征，建立一個高度為100 mm（Y 方向），直徑為50 mm 的標(biāo)準(zhǔn)巖-煤-巖組合體，其中煤體受到巖體的夾持。在FLAC3D 中，首先利用generate zone 的命令將該煤巖體模型總共劃分了1 000 個單元，并生成標(biāo)準(zhǔn)圓柱形的模型，同時保證各個網(wǎng)格節(jié)點的連接，使得在后續(xù)的模擬過程中保證其受力均勻。該模型與試驗使用的標(biāo)準(zhǔn)試樣一致，能夠滿足模型的模擬。在完成了初步的試驗?zāi)Ｐ鸵院螅疚膯屋S模擬試驗采用摩爾-庫倫塑性模型，該模型可應(yīng)用于煤巖體的實際承載能力和失效荷載的計算，以及其他以煤巖體破壞為關(guān)鍵因素的數(shù)值計算，其中g(shù)roup1 及group3 為巖體，group2 為煤體。初步模型如圖2所示。

圖2 初步模型圖

對模型底部施加約束條件=0，其余方向為自由邊界，對該模型進(jìn)行單軸壓縮數(shù)值模擬，并利用FISH 語言自動記錄應(yīng)力—應(yīng)變的數(shù)值大小，煤巖體中各物理力學(xué)參數(shù)如表1所示。

表1 煤巖物理力學(xué)參數(shù)

在完成煤巖體數(shù)值模型建立和參數(shù)賦值以后進(jìn)行模擬計算，數(shù)值模擬加載方式采用沿著軸線方向施加0.5 MPa/s的軸向應(yīng)力，在其余邊界上不施加任何應(yīng)力，即在軸向上σ=，而在環(huán)向及其他方向上==0；這樣的加載方式有利于保證與實際的實驗加載方式一致，以確保實驗的準(zhǔn)確性，同時確保數(shù)值模擬計算能夠順利進(jìn)行。

2.2 煤巖組合體應(yīng)力及位移分析

以200 時步遞增的順序?qū)γ簬r體進(jìn)行單軸壓縮試驗，并以穿過煤巖體試件軸線的平面對模擬結(jié)果進(jìn)行切片，分別記錄step200、step400 以及step600 時的應(yīng)力與位移演化云圖，具體模擬過程如圖3、圖4所示，其中圖3 為應(yīng)力分布切片，圖4為位移分布切片?？傮w來看，煤巖體試件應(yīng)力及應(yīng)變變化呈現(xiàn)以圓柱試件中心軸為對稱軸的高度對稱變化，煤巖體試件在Y軸方向上的位移上部呈下沉式形變，底部呈現(xiàn)上升式形變，且頂?shù)酌嫘巫兞孔畲?，向煤巖體試件中部位置依次衰減，隨著加載時間的增加，高位移區(qū)逐漸向試件中心部分延展。

圖3 Y 方向應(yīng)力演化云圖

圖4 Y 方向位移演化云圖

煤巖體試件應(yīng)力變化情況大致可分為Ⅰ、Ⅱ、Ⅲ三個階段：

（1）當(dāng)加載運行至200 時步時，此時為煤巖組合體單軸加載的初期，即變化的第Ⅰ階段。受上下部荷載影響，巖體試件在Y 軸方向的應(yīng)力方向為Y 軸負(fù)方向，且?guī)r體中心上部中心部位應(yīng)力最大，煤體內(nèi)部由于受巖體應(yīng)力的作用，呈現(xiàn)與巖體應(yīng)力方向相反的應(yīng)力，并在煤體中心區(qū)域形成較大的應(yīng)力集中區(qū)，此時煤體內(nèi)部的應(yīng)力略大于巖體內(nèi)部的應(yīng)力，而煤體內(nèi)部暫無位移變化，巖體沿著Y 方向開始向煤體進(jìn)行位移，且上下偏移量的數(shù)值相同。

（2）當(dāng)加載運行至400 時步時，此時為煤巖體組合單軸加載的中期，即變化的第Ⅱ階段。隨著加載時間的增加，巖體試件的高應(yīng)力區(qū)范圍逐漸增加，且影響范圍逐漸向煤巖組合體的中心處延展。煤體中心處應(yīng)力轉(zhuǎn)變?yōu)榕c巖體相同的Y 軸負(fù)方向，且煤體兩側(cè)的應(yīng)力達(dá)到20.47 MPa，煤體中心處應(yīng)力最大，最大可達(dá)到62.3 MPa。

（3）當(dāng)加載運行至600 時步時，此時為煤巖體組合體單軸壓縮的后期，即變化的第Ⅲ階段，此時煤體上下部巖體的應(yīng)力集中逐漸消失，加載對巖體產(chǎn)生的影響能夠較快速度的進(jìn)行傳遞，煤體兩側(cè)的應(yīng)力達(dá)到20.32 MPa，相比較于第Ⅱ階段，煤體兩側(cè)的應(yīng)力逐漸向煤體中間靠攏，煤巖組合體整體中間部分應(yīng)力最大，此時最大可達(dá)到72.9 MPa，巖體的位移變化也逐漸向內(nèi)部擴(kuò)展。

煤巖組合體在三個階段都具有明顯的變化特征，即Ⅰ階段煤體與巖體試件具有應(yīng)力方向相反的應(yīng)力集中區(qū)域，Ⅱ階段煤體應(yīng)力與Ⅰ階段媒體應(yīng)力相反，Ⅲ階段巖體應(yīng)力集中區(qū)域消失，同時，應(yīng)力呈現(xiàn)出由巖體傳遞向煤體，由煤體兩側(cè)傳遞向煤體中心的趨勢，通過圖4 同樣可以看出，隨著加載時步的增加，煤巖體Y 方向的位移由巖體逐漸向煤體轉(zhuǎn)移，大小呈現(xiàn)遞增的趨勢。

圖4 中國獨立動畫理論研究突顯詞視圖（左）與控制面板參數(shù)設(shè)置（右）

2.3 煤巖組合體應(yīng)力與應(yīng)變變化關(guān)系分析

強(qiáng)度作為衡量煤巖體破壞狀態(tài)的重要力學(xué)性質(zhì)之一，指煤巖體受到外荷載作用時抵抗破壞時的最大能力，為研究模擬過程中煤巖體所承受的最大破壞能力以及煤巖體加載過程中應(yīng)力與應(yīng)變之間的關(guān)系，以模擬過程中煤巖體的應(yīng)變?yōu)樽宰兞浚簬r體的應(yīng)力作為因變量，得出煤巖體單軸加載過程中的應(yīng)力—應(yīng)變曲線如圖5所示。

由圖5 可以看出煤巖體應(yīng)力隨應(yīng)變大致呈現(xiàn)先快速增加后逐漸減緩的趨勢。隨著加載時間的增加，固定荷載對煤巖體施加軸向應(yīng)力使得煤體產(chǎn)生彈性形變，煤巖體內(nèi)部應(yīng)力隨著應(yīng)變快速增加，但隨著應(yīng)變的逐漸增大，應(yīng)力的增長速度逐漸減緩，在應(yīng)變量達(dá)到5.0×10之前，整個煤巖體處于一個彈性變形階段，當(dāng)應(yīng)變達(dá)到1.0×10時，煤體應(yīng)力隨應(yīng)變增加具有先小幅下降后小幅上升的趨勢，此時煤巖體處于屈服階段，當(dāng)應(yīng)變增加至1.5×10時，應(yīng)力基本維持在一個較高的水平，其數(shù)值大致為32.5 MPa，即該煤巖體所承受的最大破壞能力約為32.5 MPa，若超過該值，則煤巖體發(fā)生延性破壞。

圖5 應(yīng)力—應(yīng)變曲線圖

3 結(jié)論

通過FlAC3D 的數(shù)值模擬結(jié)果可以看出，隨著加載時間的增加，煤巖組合體內(nèi)部應(yīng)力及位移變化具有較強(qiáng)的規(guī)律性：

（1）煤巖組合體試件在Y 軸方向上的位移上部呈下沉式形變，底部呈現(xiàn)上升式形變，且頂?shù)酌嫘巫兞孔畲?，向試件中部位置依次衰減。

（2）煤巖組合體在加載過程中應(yīng)力變化呈現(xiàn)特征較為明顯的Ⅰ、Ⅱ、Ⅲ三個階段的變化，且Ⅰ階段煤體與巖體試件具有應(yīng)力方向相反的應(yīng)力集中區(qū)域，Ⅱ階段煤體應(yīng)力與Ⅰ階段媒體應(yīng)力相反，Ⅲ階段巖體應(yīng)力集中區(qū)域消失，同時應(yīng)力呈現(xiàn)出由巖體傳遞向煤體，由煤體兩側(cè)傳遞向煤體中心的趨勢。

（3）煤巖組合在加載的過程中，煤巖體應(yīng)力隨應(yīng)變大致呈現(xiàn)先快速增加后逐漸減緩的趨勢，且隨著加載時間的增加，煤巖體經(jīng)歷了彈性變形、屈服階段、延性破壞三個變形階段。