摘要：在應(yīng)力跌落模型的基礎(chǔ)上引進(jìn)軟化閾值，建立了彈塑脆性模型，模擬巖土材料的脆性軟化性質(zhì)?；贛ohr-Coulomb準(zhǔn)則，考慮了巖土材料屈服后的塑性軟化和體積膨脹，推導(dǎo)了圓形巷道圍巖的軟化區(qū)半徑、塑性區(qū)半徑、洞周位移及圍巖內(nèi)應(yīng)力表達(dá)式，最后通過算例分析了剪脹、軟化程度和彈模劣化對(duì)破裂區(qū)范圍的影響，為巷道的穩(wěn)定分析以及支護(hù)設(shè)計(jì)提供理論依據(jù)。

關(guān)鍵詞：巷道；彈塑性；剪脹；軟化；解析解

中圖分類號(hào)：TU451.2

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1674-4764（2013）03-0007-05

Elastic-Plastic Analysis of Surrounding Rocks Considering Material

Dilatancy and Softening

Wei Jianjun1，2

（1.State Key Laboratory for Geomechanics and Deep Underground Engineering， China University of Mining Technology，

Xuzhou 221008， Jiangsu，P.R.China；2. RD Centre of Jiangsu Construction Safety and Disaster Mitigation，

Jiangsu Institute of Architectural Technology，Xuzhou 221116， Jiangsu，P.R.China）

Abstract：

Based on stress dropping model， an elastoplastic-brittle model was established by using a softening threshold value to model the features of the geotechnical materials. According to the Mohr-Coulomb criterion， the radius of the rupture and plastic zone， the displacement and stress expressions of surrounding rock were obtained with the consideration of the strain-softening and dilatancy of geotechnical materials. Furthermore， the model was verified by case study and the influence of strain-softening， dilatancy and elastic modulus deterioration on calculation result were presented.

Key words：

circular opening；elastoplastic；dilation；strain-softening；analytical solution

隨著礦山開采的大規(guī)模進(jìn)行和采深的增加，造成的災(zāi)害日益增多，對(duì)深部資源的安全高效開采造成巨大威脅。深部巖體處于高地應(yīng)力、高溫、高滲透壓的惡劣環(huán)境中，使得巖體的組織結(jié)構(gòu)、力學(xué)性態(tài)和工程響應(yīng)均發(fā)生根本性變化，圍巖表現(xiàn)出特殊的非線性力學(xué)行為[1-3]。深部巷道開挖后二次應(yīng)力場(chǎng)形成引起的高地應(yīng)力集中導(dǎo)致圍巖壓剪應(yīng)力超過其強(qiáng)度進(jìn)而破壞處于剪脹狀態(tài)；同時(shí)，開挖后的近表圍巖內(nèi)空隙水壓力降低引起的圍巖有效應(yīng)力增大加劇了圍巖的破裂；高地溫帶來的附加應(yīng)力對(duì)圍巖破裂擴(kuò)展帶了不可忽視的影響。巖石進(jìn)入峰后階段由于應(yīng)力跌落方式的不確定性，難以用經(jīng)典理論來描述[4-5]。Hoek等[6]于1997年提出了估算巖石峰后力學(xué)行為的方法，給出了巖石峰后力學(xué)模型。隨后，學(xué)者對(duì)圍巖的彈塑性分析進(jìn)行了大量研究。其他國(guó)家的研究集中于基于Hoek-Brown強(qiáng)度準(zhǔn)則采用彈脆塑性模型進(jìn)行分析[7-10]，中國(guó)的研究則考慮巖石的非線性軟化和破壞過程，但在解析計(jì)算方面為了求解方便往往設(shè)置不合理的限制，如在圍巖破壞過程中體應(yīng)變?cè)O(shè)置為零[11-12]、彈性模量設(shè)為常數(shù)[13]等，這些方法均不符合巷道受擾變形的機(jī)制。

圓形硐室的軸對(duì)稱彈塑性分析是個(gè)經(jīng)典問題，有許多學(xué)者進(jìn)行過研究。文獻(xiàn)[8]對(duì)Mohr-Coulomb介質(zhì)中的應(yīng)變軟化過程進(jìn)行了研究，其中考慮了黏結(jié)力和剪脹角的下降，內(nèi)摩擦角在軟化變形中保持不變；文獻(xiàn)[7]將圓形硐室圍巖內(nèi)形成的塑性區(qū)劃分為一系列的同心圓環(huán)，采用差分法求解得到軟化區(qū)與殘余強(qiáng)度區(qū)尺寸以及相應(yīng)的位移分布；文獻(xiàn)[14-16]考慮了統(tǒng)一強(qiáng)度理論，軟化段采用雙折線模型，得到了硐室圍巖彈塑性分析的統(tǒng)一解。

魏建軍：考慮剪脹和軟化的巷道圍巖彈塑性分析

本文根據(jù)巖土體材料在出現(xiàn)應(yīng)力跌落以前或多或少要經(jīng)過一段塑性變形的過程，考慮引入應(yīng)力跌落條件閾值來判斷是否出現(xiàn)應(yīng)力跌落，構(gòu)建一個(gè)彈塑脆性模型，考慮圍巖的應(yīng)變軟化和剪脹特性，采用Mohr-Coulomb屈服準(zhǔn)則對(duì)圓形硐室圍巖進(jìn)行彈塑性分析。

1彈塑脆性模型及問題描述

前人所采用的應(yīng)力跌落模型的應(yīng)力應(yīng)變曲線由兩段組成，即屈服前的線彈性階段和屈服后的殘余強(qiáng)度階段，即認(rèn)為材料在達(dá)到彈性極限之前為線彈性，一旦達(dá)到彈性極限，便會(huì)屈服，并且強(qiáng)度降為殘余強(qiáng)度。其實(shí)，真正彈脆性材料是不存在的，材料在出現(xiàn)應(yīng)力跌落以前，或多或少要經(jīng)過一段塑性變形的過程。因此，本文采用彈塑脆性模型，如圖1所示。巷道圍巖應(yīng)力應(yīng)變曲線可分為處于彈性區(qū)、塑性發(fā)展區(qū)和殘余變形區(qū)3種狀態(tài)。

模型以微元體中與最大主應(yīng)力和最小主應(yīng)力都成45°的面上的最大剪應(yīng)變來描述巖石不可逆的塑性變形，當(dāng)最大剪應(yīng)變達(dá)到某一值時(shí)，材料即出現(xiàn)脆性軟化，該最大剪應(yīng)變的臨界值即為軟化閾值，其值在下文用γf表示。

本文研究巖土材料中平面應(yīng)變條件下的圓孔問題。如圖2所示，在彈塑脆性材料的無限空間中作用有靜水壓力p0，其中有一無限長(zhǎng)圓形硐室，開挖半徑為a；支護(hù)荷載為pi；隨著圍巖卸荷發(fā)展，硐室周圍巖土體進(jìn)入塑性，貼近洞口部分區(qū)域剪應(yīng)變達(dá)到軟化閾值，材料出現(xiàn)軟化形成環(huán)形軟化區(qū)，這樣巷道圍巖自內(nèi)向外由軟化區(qū)、塑性區(qū)和彈性區(qū)3部分組成，為書寫方便，其對(duì)應(yīng)的物理量分別附下標(biāo)s、p、e以示區(qū)別。

2基本假定與方程

做如下基本假定：

1）圍巖材料均勻、連續(xù)、各向同性。

2）材料屈服服從Mohr-Coulomb準(zhǔn)則。

3）圍巖材料具有明顯的剪脹特性，塑性及軟化區(qū)體積應(yīng)變不等于0。采用線性剪脹模型，位移以指向洞內(nèi)為正，設(shè)剪脹參數(shù)為βi，則塑性應(yīng)變?chǔ)舙1與εp3的關(guān)系為：

4算例與討論

采用文獻(xiàn)[17]的計(jì)算實(shí)例。貴州錦屏二級(jí)水電站大水溝廠房支1洞位于東雅礱江右岸，地質(zhì)巖性主要為大理巖，截面尺寸為3.26 m×3.30 m（寬×高），近似為圓形硐室。經(jīng)過地應(yīng)力實(shí)測(cè)表明探洞斷面上豎直主應(yīng)力約為22.9 MPa，水平主應(yīng)力約為19.8 MPa。

結(jié)合現(xiàn)場(chǎng)測(cè)試與設(shè)計(jì)院試驗(yàn)結(jié)果確定本研究段的巖體力學(xué)計(jì)算參數(shù)為：彈性模量Ee=Ep=16 GPa，Es=2 GPa；泊松比υ=0.3；在彈性區(qū)和塑性區(qū)粘聚力cp=2.5 MPa，摩擦角φp=35°；在塑性軟化區(qū)粘聚力cs=1.4 MPa，摩擦角φ=23°。取剪脹角ψp=18°，ψs=13°，巖石的應(yīng)變軟化閾值分別取0.001、0.004、0.01。地應(yīng)力近似為兩向等壓，取p0=21 MPa；支護(hù)荷載pi=0。通過計(jì)算得到其對(duì)應(yīng)的軟化區(qū)及塑性區(qū)范圍見表1。圍巖的應(yīng)力與應(yīng)變分布如圖3所示。

從圖3可以看出，軟化閾值對(duì)分析結(jié)果影響顯著，隨著軟化閾值減小，軟化區(qū)和塑性區(qū)范圍增大，在相同松動(dòng)圈半徑位置，應(yīng)力隨之減小，應(yīng)變則增大。

圖4為軟化閾值取0.004時(shí)，圍巖彈模劣化系數(shù)k分別取0.125、0.300、0.600三種情況下的應(yīng)變曲線。從結(jié)果可以看出，近巷圍巖彈模劣化是導(dǎo)致深部巷道圍巖變形較大的因素之一，不考慮圍巖彈模劣化計(jì)算結(jié)果偏危險(xiǎn)。

5結(jié)論

1）在理想彈脆性模型的基礎(chǔ)上，加上一個(gè)軟化閾值，建立了彈塑脆性模型，研究巖土開挖中的硐室圍巖彈塑性問題，基于M-C強(qiáng)度理論和非關(guān)聯(lián)流動(dòng)法則，得到了圓形硐室圍巖應(yīng)力與變形解析表達(dá)式。

2）軟化閾值變化對(duì)圍巖應(yīng)力場(chǎng)影響顯著，隨著軟化閾值的增大，塑性區(qū)圍巖應(yīng)力增大，變形減小，反之亦然。

3）分析了巖石彈模劣化系數(shù)對(duì)圍巖變形的影響。隨著劣化程度的加大，破裂區(qū)范圍變化不大，圍巖變形顯著，是地下工程近巷圍巖大變形的重要原因。

參考文獻(xiàn)：

[1]周宏偉，謝和平，左建平.深部高地應(yīng)力下巖石力學(xué)行為研究進(jìn)展[J].力學(xué)進(jìn)展，2005，35（1）：91-99.

Zhou H W，Xie H P，Zuo J P.Developments in researches on mechanical behaviors of rocks under the condition of high ground pressure in the depths [J].Advances in Mechanics， 2005，35（1）： 91-99.

[2]王德榮，甄樹新，汪新紅，等.深部坑道圍巖壓力與變形分析[J].巖土力學(xué)，2007，28（3）：70-77. Wang D R，Zhen S X，Wang X H，et al. Analysis of pressure and deformation of wall rock near deep tunnel [J].Rock and Soil Mechanics，2007，28（3）：70-77.

[3]劉泉聲，盧興利.煤礦深部巷道破裂圍巖非線性大變形及支護(hù)對(duì)策研究[J].巖土力學(xué)，2010，31（10）：3273-3280.

Liu Q S，Lu X L.Research on nonlinear large deformation and support measures for broken surrounding rocks of deep coal mine roadway [J].Rock and Soil Mechanics，2010，31（10）：3273-3280.

[4]張黎明，王在泉，李華峰，等.粉砂巖峰后破壞區(qū)應(yīng)力脆性跌落的試驗(yàn)和本構(gòu)方程研究[J].實(shí)驗(yàn)力學(xué)，2008，23（3）：234-240. Zhang L M，Wang Z Q，Li H F，et al. Theoretical and experimental study on siltstone brittle stress drop in post failure region [J].Journal of Experimental Mechanics，2008，23（3）：234-240.

[5]王水林，王威，吳振君.巖土材料峰值后區(qū)強(qiáng)度參數(shù)演化與應(yīng)力-應(yīng)變曲線關(guān)系研究[J].巖石力學(xué)與工程學(xué)報(bào)，2010，29（8）：1524-1529. Wang S L，Wang W，Wu Z J.Study of relationship between evolution of post-peak strength parameters and stress-strain curves of geomaterlals [J].Chinese Journal of Rock Mechanics and Engineering，2010，29（8）： 1524-1529.

[6]Hoek E，Brown T E.Practical estimates of rock mass strength [J].International Journal of Rock Mechanics and Mining Sciences，1997，34（8）：1165-1186.

[7]Sharan K S.Eastic-brittle-plastic analysis of circular openings in Hoek-Brown media [J].International Journal of Rock Mechanics and Mining Sciences，2003， 40（6）： 817-824.

[8]Lee K Y，Pietruszczak S.A new numerical procedure for elasto-plastic analysis of a circular opening excavated in a strain-softening rock mass [J].Tunnelling and Underground Space Technology，2008，23（5）：588-599.

[9]Carranza T C，F(xiàn)airhurst C.The elasto-plastic response of underground excavations in rock masses that satisfy the Hoek-Brown failure criterion [J].International Journal of Rock Mechanics and Mining Sciences，1999，36（6）：777-809.

[10]Chandra S，Nilsen B， Ming L. Predicting excavation methods and rock support： a case study from the Himalayan region of India [J].Bulletin of Engineering Geology and the Environment，2009，69（2）：257-266.

[11]潘岳，王志強(qiáng).基于應(yīng)變非線性軟化的圓形硐室圍巖彈塑性分析[J].巖石力學(xué)與工程學(xué)報(bào)，2005，24（6）：915-921.

Pan Y， Wang Z Q.Elasto-plastic analysis on surrounding rock of circular chamber based on strain nonlinear softening [J]. Chinese Journal of Rock Mechanics and Engineering，2005，24（6）：915-921.

[12]王亮，趙均海，李小偉.巖質(zhì)圓形隧洞圍巖應(yīng)力場(chǎng)彈塑性新解[J].工程地質(zhì)學(xué)報(bào)，2007，15（3）：424-428. Wang L，Zhao J H，Li X W.New elasto-plastic solution of rock circular tunnel surrounding rockmass stress field [J].Journal of Engineering Geology，2007，15（3）：424-428.

[13]蔣斌松，張強(qiáng)，賀永年，等.深部圓形巷道破裂圍巖的彈塑性分析[J].巖石力學(xué)與工程學(xué)報(bào)，2007，26（5）：983-986. Jiang B S， Zhang Q， He Y N， et al. Elastoplastic analysis of cracked surrounding rocks in deep circular openings[J]. Chinese Journal of Rock Mechanics and Engineering， 2007，26（5）： 983-986.

[14]范文，俞茂宏，陳立偉，等.考慮剪脹及軟化的洞室圍巖彈塑性分析的統(tǒng)一解[J].巖石力學(xué)與工程學(xué)報(bào)，2004，23（19）：3213-3221. Fan W，Yu M H，Chen L W，et al.Unified elastoplastic solution for surrounding rocks of openings with consideration of material dilatancy and softening [J].Chinese Journal of Rock Mechanics and Engineering， 2004，23（19）：3213-3221.

[15]李同錄，陳立偉，俞茂宏，等.考慮材料軟化的洞室圍巖彈塑性分析的統(tǒng)一解[J].長(zhǎng)安大學(xué)學(xué)報(bào)：自然科學(xué)版，2004，24（3）：49-52. Li T L，Chen L W， Yu M H， et al. Unified solution of elastic-plastic surrounding rocks of cavity considering material softening [J].Journal of Changan University： Natural Science Edition， 2004， 24（3）：49-52.

[16]張常光，胡云世，趙均海，等.深埋圓形水工隧洞彈塑性應(yīng)力和位移統(tǒng)一解[J].巖土工程學(xué)報(bào)，2010，32（11）：1738-1746. Zhang C G，Hu Y S，Zhao J H，et al.Elastic-plastic unified solutions for stresses and displacements of a deep buried circular hydraulic tunnel [J].Chinese Journal of Geotechnical Engineering， 2010，32（11）：1738-1746.

[17]江權(quán).高地應(yīng)力下硬巖彈脆塑性劣化本構(gòu)模型與大型地下洞室群圍巖穩(wěn)定性分析[D].武漢：中國(guó)科學(xué)院武漢巖土力學(xué)研究所，2007.