劉 學 王來貴 李宏艷

(1.煤炭科學研究總院礦山安全技術研究分院；2.煤炭資源高效開采與潔凈利用國家重點實驗室；3.遼寧工程技術大學力學與工程學院)

不同均質度下巷道壓剪損傷算法與模擬研究*

劉 學1，2王來貴3李宏艷1，2

(1.煤炭科學研究總院礦山安全技術研究分院；2.煤炭資源高效開采與潔凈利用國家重點實驗室；3.遼寧工程技術大學力學與工程學院)

基于韋伯分布理論，依據非均質巖石剪切破壞準則，采用彈性全量法分級加載，通過Fortran語言編寫出非均質巖石剪切破壞有限元程序,模擬不同均質度下巷道頂板單軸受壓作用下的破壞過程。研究結果為：應力場和變形場明顯受到巖體非均質性的影響，其損傷區(qū)域分布表現出無序、隨機的特點；均質度是影響巷道圍巖宏觀強度的重要參數，覆壓作用下的最終破壞形式以剪切破壞為主。該研究成果為預測巷道結構整體穩(wěn)定性提供了方法和依據。

巷道結構 非均質 剪切破壞 有限元

巖石是一種天然材料，是由礦物晶粒集合體和膠結材料的非晶體顆粒組成，在其形成的漫長地質年代中遭受到復雜的地質作用，使其力學性質在空間上表現出極大的非均勻性。巖石作為一種具有內部結構的非均質材料，在外力作用下，內部微缺陷的成核、擴展及相互之間的作用,決定了巖石的宏觀變形、破壞特性。盡管人們在研究巖石破裂方面取得了諸多有益的成果[1-3]，對均質性研究已經成熟，但對非均質性巖體破壞的研究尚處于起步階段[4-7]，很少涉及到巖石自身非均質性引起的壓剪損傷破壞。而在工程破壞中是否考慮巖體非均質性，將會得到不同的結果。因此，上覆巖層的非均質性作為引起井巷頂板壓剪破壞原因之一，是不容忽視的重要因素，其破壞前兆和破裂后承載能力變化是采礦工作者重點關心的問題。

1 Weibull分布

Weibull率先提出了用統(tǒng)計數學描述材料非均勻性的方法。假定組成材料單元的力學性質滿足Weibull分布(或其它隨機分布形式)，對尺度效應、強度理論的研究起到了重要作用[8-12]。

韋伯分布的概率密度函數為：

(1)

式中，x為材料介質力學性質參數(如巖石強度、彈性模量等)；m為分布函數的形狀參數，決定分布密度曲線的基本形狀；α為分布函數的尺度參數，為介質力學性質的平均值。

隨著m值的增大，巖石介質的力學性質將集中在一個狹窄的范圍內，這表明巖石介質的性質較均勻；而當m值減小時，則巖石介質的力學性質分布范圍變寬，表明巖石介質的性質趨于非均質。如圖1所示。

圖1 Weibull概率密度函數(雙參數)

2 均質度與巖體損傷關系研究

2.1 荷載、材料參數和邊界條件

巷道結構如圖2所示。彈性模量型按Webull隨機分布(如圖3)，最大值為1.53E10，最小值為3.28E9，泊松比為0.3，容重為3 kN/m3。模型底邊假設無限遠邊界，施加為全約束，頂部施加均布荷載，從171 MPa施加載荷，每時步增量為1 MPa。

2.2 不同均質度下巷道頂板壓剪損傷破壞模擬

采用Weibull統(tǒng)計分布函數描述巷道模型的細觀非均質性，利用FEPG模擬不同均質度模型在單軸受壓作用下的損傷破壞過程。由以上分析可知，Weibull分布中的均質度系數m是描述力學參數的均勻程度。

圖2 巷道結構力學模型

圖3 隨機彈性模量

選取m=0.5,m=1,m=2,m=5,分析不同均質度對巷道模型壓剪損傷的變化影響。見圖4～圖7所示。

圖4 m=0.5模型損傷演化

圖5 m=1模型損傷演化

圖6 m=2模型損傷演化

圖7 m=5模型損傷演化

由不同m值巷道模型損傷演化圖可以看出，均質度很小時，損傷區(qū)域隨著載荷步的施加會逐漸擴張聯(lián)通，直至整個巷道結構大面積出現破壞，巷道模型最容易出現損傷。當m=0.5時，模型在加載第一步就出現了損傷域，說明均質度越小，巖體介質的力學性質分布范圍變寬，表明巖體介質的性質趨于非均質，損傷區(qū)域越大；當均質度很大時，巖體介質的力學性質將集中在一個狹窄的范圍內，性質趨于均質，當m=5時，巷道模型出現了很狹小的損傷區(qū)域，模型壓剪破壞較慢。

為了從細觀角度分析不同m值對巷道巖體損傷度的影響，選取了具有代表性的監(jiān)測點1、2進行跟蹤監(jiān)測(如圖8)，并進行曲線擬合對比分析(如圖9，圖10)。

由圖9、圖10可以看出，隨著m值增大，損傷出現的越慢，說明巖體結構整體性質越均質，宏觀強度越大；而當m值越小時，巖體結構均勻性越差，力學性質分布不一，造成巖體宏觀強度越低，在受壓荷載作用下，越容易出現損傷，巷道抵抗破壞的能力越低。

圖8 監(jiān)測點布置

圖9 不同m值下監(jiān)測點1損傷演化

圖10 不同m值下監(jiān)測點2損傷演化

3 結 論

(1)通過模擬受壓荷載作用下的非均質巷道破壞狀態(tài)，得出應力場和變形場明顯受到巖體非均質性的影響，其損傷區(qū)域分布表現出無序、隨機的特點。

(2)充分體現了覆壓作用下的最終破壞形式以剪切破壞為主，其損傷破壞起源于模型內部。由于模型內部無應力集中，即無應力釋放現象，為后續(xù)巷道圍巖壓剪破壞的計算奠定理論基礎。

(3)出現了非對稱損傷區(qū)域，且損傷區(qū)域為不光滑絮狀，是由于彈性模量不均勻所致；隨著加載步的遞增，巷道損傷區(qū)域在前一步損傷破壞的基礎上繼續(xù)沿45°延展，頂板上部呈現漏斗狀損傷區(qū)，與實際井巷開采形成的沉陷區(qū)域吻合。

(4)均質度是影響巖體宏觀強度的重要參數，當均質度越大時，巖體結構宏觀強度越大，力學性質越均勻，抗壓剪能力越強；當均質度越小時，巖體結構宏觀強度越小，力學性質越不均勻，抗壓剪能力越差。因此，監(jiān)測巷道圍巖的均質度變化，可為預測預報巷道結構整體穩(wěn)定性提供方法和依據。

[1]Muneohori,JunochiroKubo.Analysisofprobabilityandrangeofaveragestressineachphaseofheterogeneitymaterials[J].JournaloftheMechanicsandPhysicsofSolids,1998,46(3):537-556.

[2]TangCA.Numericalsmulationofprogressiverockfailureandassociatedseismicity[J].InternationalJournalofRockMechanicsandMiningSciences,1997,34(2):249-261.

[3] 唐春安.巖石破裂過程中的災變[M].北京:煤炭工業(yè)出版社，1993.

[4] 陳永強.非均勻材料有效力學性能和破壞過程的數值模擬[D].北京:清華大學，2001.

[5] 梁正召.三維條件下的巖石破裂過程分析及其數值試驗方法研究[D].沈陽:東北大學，2005.

[6] 曹文貴，趙明華，劉成學.基于Weibull分布的巖石損傷軟化模型及其修正方法研究[J].巖石力學與工程學報，2004,23(19):3223-3231.

[7] 唐春安，劉紅元，秦四清，等.非均勻性對巖石介質中裂紋擴展模式的影響[J].地球物理學報，2000，43(1)：116-121.

[8] 王學濱.材料缺陷對巖樣變形局部化影響的數值模擬[J].巖土力學，2006,27(8):1241-1247.

[9] 王學濱，潘一山.加載速度對隨機缺陷巖樣破壞過程的影響[J].巖石力學與工程學報,2007(zl):3493-3497.

[10] 王來貴,趙 娜,初 影,等. 不同面積載荷作用下的巖石試件破裂數值模擬[J].沈陽建筑大學學報：自然科學版,2007(6):44-47.

[11] 李宏艷. 煤巖物理力學性質與沖擊傾向性關系[J].煤礦開采,2011(3):43-46.

[12] 謝和平，彭瑞東，周宏偉，等.基于斷裂力學與損傷力學的巖石強度理論研究進展[J].地下空間與工程學報,2004(6):1076-1081.

Simulation Research of the Compression-shear Damage Algorithm of Roadway with Different Degree of Heterogeneity

Liu Xue1,2Wang Laigui3Li Hongyan1,2

(1.Mine Safety Technology Branch, China Coal Research Institute;2.State Key Laboratory of Coal Mining and Clean Utilization;3.School of Mechanics and Engineering,Liaoning Technical University)

Based on the Wei-bull distribution theory and the shear failure criterion of the heterogeneous rock stage with loading elastic total method under step loading, the finite element program of heterogeneous shear failure of rock is wrote with FORTRAN language so as to simulate the characteristics of roadway roof failure process under unixial compression homogeneity. The simulation results show that the stress and deformation fields are influenced significantly by the rock mass heterogeneity;characteristics of the damaged areas distribution is disordered and random; homogeneity degree is the important parameter that affect the macroscopic strength of roadway surrounding rock,it is fully embodies that the damage forms is given priority to shear failure in the end.The above research results can provide the basis and method for predicting the overall structural stability of roadway.

Roadway structure, Heterogeneity, Shear failure, Finite element

*國家重點基礎研究發(fā)展計劃(973)資助項目(編號:2012CB724208)；國家自然科學基金項目(編號:51174112)。

2015-04-24)

劉 學(1987—)，女，助理工程師，碩士，100013 北京市朝陽區(qū)和平里青年溝東。