羅 勇

(中南大學(xué) 資源與安全工程學(xué)院，湖南 長沙 410083)

0 引言

基于M-C準(zhǔn)則的Kastner方程，假定圍巖為理想彈塑性材料，被廣泛應(yīng)用于巷道彈塑性分析[1-2]。近年來，對(duì)巷道彈塑性分析開展了許多相關(guān)研究，得到了一些有益的成果。例如，蔣斌松等[3]利用非關(guān)聯(lián)流動(dòng)法則和M-C準(zhǔn)則對(duì)巷道變形進(jìn)行分析；孟龍等[4]基于M-C準(zhǔn)則，得到了圍巖應(yīng)力解；有研究表明[5-6]中間主應(yīng)力可使巖體強(qiáng)度增大30%左右；侯公羽等[7]利用D-P準(zhǔn)則推導(dǎo)了理想彈塑性巷道的解析解；胡小榮等[8]、曾開華等[9]考慮中間主應(yīng)力對(duì)變形的影響；劉萬榮等[10]、劉志欽等[11]分析了中間主應(yīng)力巷道應(yīng)力的影響解析式。上述這些研究雖分析了中間主應(yīng)力對(duì)彈塑性解的影響，但均沒有考慮巖體的應(yīng)變軟化行為。為此，基于非關(guān)聯(lián)流動(dòng)法則和D-P準(zhǔn)則，考慮中間主應(yīng)力、圍巖軟化和膨脹，推導(dǎo)了巷道彈性解，并進(jìn)行應(yīng)用分析。

1 力學(xué)模型與屈服準(zhǔn)則

1.1 巷道力學(xué)模型

設(shè)在原巖應(yīng)力為p0的地層中開挖一半徑為R0的圓形巷道，表面作用有支護(hù)阻力pi，圍巖依次形成彈性區(qū)、塑性區(qū)和破裂區(qū)，如圖1。

圖1 巷道力學(xué)模型Fig.1 Mechanical model of roadway

圖中Rb，Rp分別為破裂區(qū)和塑性區(qū)半徑。設(shè)水平無限長巷道深度大于20倍硐徑；原巖應(yīng)力各向等壓，不計(jì)影響范圍內(nèi)巖體重力，圍巖視為連續(xù)、均質(zhì)和各向同性材料。

1.2 應(yīng)變軟化

為了描述巖體的軟化行為，設(shè)軟化階段φ和c與塑性應(yīng)變?chǔ)舙滿足二次函數(shù)關(guān)系[12]，如圖2所示。

圖2 c和φ軟化模型Fig.2 Softening model of c and φ

設(shè)二次函數(shù)關(guān)系式為：

c=F(εp-ε0)2+K(εp-ε0)+D

(1)

式中：F，K，D為常數(shù)。在峰值強(qiáng)度(ε0，c0)和殘余強(qiáng)度(εb，cb)處，由圖2可得：

c0=D，cb=F(εb-ε0)2+K(εb-ε0)+D

(2)

在c減至cb后，塑性應(yīng)變不再使其變化，即：

(dc/dεp)εp=εb=2F(εb-ε0)+K=0

(3)

聯(lián)立式(1)～(3)可得：

c=Sc(εp-ε0)2/(εb-ε0)-2Sc(εp-ε0)+c0

(4)

同理：

φ=Sφ(εp-ε0)2/(εb-ε0)-2Sφ(εp-ε0)+φ0

(5)

式中：Sc=(c0-cb)/(εb-ε0)，為c的軟化模量；Sφ=(φ0-φb)/(εb-ε0)，為φ的軟化模量；ε0，φ0，c0分別是峰值處的應(yīng)變、內(nèi)摩擦角和黏聚力，單位分別為無量綱、(°)和MPa。

1.3 破裂膨脹

圍巖破壞后具有破裂膨脹行為[13-14]，可利用非關(guān)聯(lián)流動(dòng)法則確定切向和徑向塑性應(yīng)變的關(guān)系[9]：

塑性區(qū)：

(6)

破裂區(qū)：

(7)

式中：β1和β2是與巖體破裂膨脹有關(guān)的系數(shù)。

1.4 Drucker-Prager屈服準(zhǔn)則

D-P準(zhǔn)則表達(dá)式為：

(8)

n=(σ2-σ3)/(σ1-σ3)

(9)

σθ-σr=(1-M)σr-N

(10)

考慮應(yīng)變軟化，則在塑性區(qū)有：

σθ-σr=(1-Mp)σr-Np

(11)

在破裂區(qū)有：

σθ-σr=(1-Mb)σr-Nb

(12)

式中：σr和εr，σθ和εθ，σz分別為徑向應(yīng)力和應(yīng)變、切向應(yīng)力和應(yīng)變以及巷道軸向應(yīng)力，應(yīng)力單位為MPa。

2 巷道圍巖應(yīng)力與位移

2.1 基本方程

平衡微分方程：

dσr/dr+(σr-σθ)/r=0

(13)

幾何方程：

εr=du/dr，εθ=u/r

(14)

物理方程：

(15)

式中：u，v，E分別為位移、泊松比和彈性模量，單位分別為mm、無量綱和MPa。

2.2 彈性區(qū)應(yīng)力與位移

對(duì)于彈性區(qū)，當(dāng)r→時(shí)，有σr=σθ=p0；r=Rp時(shí)，σr和σθ滿足式(10)和σθ+σr=2p0，且有求解得彈性區(qū)應(yīng)力：

(16)

將式(16)代入式(15)求解，得彈性區(qū)應(yīng)變：

(17)

式(17)第2式代入式(14)，得彈性區(qū)位移：

(18)

2.3 塑性區(qū)應(yīng)力和位移

在塑性區(qū)內(nèi)，σr和σθ滿足式(11)和式(13)，把式(11)代入式(13)，積分解得：

σr=CrMp-1+Np/(1-Mp)

(19)

(20)

在塑性區(qū)，圍巖總的應(yīng)變由2部分組成，即塑性應(yīng)變和彈性應(yīng)變，由此得到：

(21)

在塑性區(qū)，巖體表現(xiàn)出破裂膨脹特性，聯(lián)立式(6)、(14)和(21)得：

(23)

式(23)代入式(14)，可得到塑性區(qū)內(nèi)應(yīng)變：

(24)

2.4 破裂區(qū)應(yīng)力與位移

(25)

由式(21)，同理可以得到破裂區(qū)總應(yīng)變，即：

(26)

(27)

2.5 塑性區(qū)和破裂區(qū)半徑

在塑性軟化區(qū)內(nèi)，由式(24)可以得：

(28)

由于應(yīng)變軟化，塑性區(qū)cp，φp與εp分別滿足式(4)和(5)，將式(28)代入得：

(29)

式中：H=(1+ν)[p0(M-1)+N]/[E(1+M)]。

當(dāng)r=Rb時(shí)，c=cb，φ=φb，代入式(29)，可以得到塑性區(qū)與破裂區(qū)的半徑比值，即：

(30)

(31)

聯(lián)立式(30)和(31)得破裂區(qū)半徑。如圍巖無破裂膨脹和應(yīng)變軟化特性，此時(shí)cb=c0，φb=φ0，β1=β2=1代入式(27)和(31)，則得到只考慮中間主應(yīng)力下的一般解答，與文獻(xiàn)[7]所得解一致。

3 算例分析

假設(shè)開挖一半徑R0為3 m的圓形巷道，p0為15 MPa，巷道表面均勻作用支護(hù)阻力pi，大小為1.2 MPa。圍巖基本參數(shù)為：E=2 GPa；μ=0.25；c0=5.3 MPa；cb=2.1 MPa；φb=16.3°；φ0=32°；Sφ=1 963°；Sc=400 MPa；β1=1.42；β2=1.4。

3.1 中間主應(yīng)力對(duì)Rp，Rb及u的影響

圖3為中間主應(yīng)力對(duì)Rp和Rb的影響，圖4為中間主應(yīng)力對(duì)巷道變形u的影響。由圖3和圖4得，n逐漸增大時(shí)，Rp，Rb及u均呈現(xiàn)先減后增的變化規(guī)律。當(dāng)0≤n≤0.5，Rp，Rb及u隨n增大減小明顯；當(dāng)0.5≤n≤0.7，Rp，Rb及u隨n增大呈先減后增的規(guī)律，但變化不明顯，且當(dāng)n=0.7左右有最小值；當(dāng)0.7≤n≤1時(shí)，Rp，Rb及u隨n增大增加越明顯，但n=1比n=0時(shí)值小。因此，中間主應(yīng)力對(duì)Rp，Rb及u的影響具有區(qū)間性。由圖3、圖4可知考慮中間主應(yīng)力時(shí)，最大可使Rp大約減少20%，Rb減少12%，u減少22%。由此可見，中間主應(yīng)力對(duì)巷道的變形和破壞影響較大。

圖3 中間主應(yīng)力對(duì)Rp和Rb的影響Fig.3 Influence of intermediate principal stress on Rp and Rb

圖4 中間主應(yīng)力對(duì)巷道變形u的影響Fig.4 Influence of intermediate principal stress on u

3.2 應(yīng)變軟化對(duì)Rp，Rb及u的影響

圖5為c軟化對(duì)Rp和Rb的影響，圖6為內(nèi)聚力軟化對(duì)位移的影響。由圖5和圖6可知，Sc對(duì)Rp不產(chǎn)生影響；Sc恒定不變時(shí)，隨cb的增大，Rp，Rb和u均逐漸減小，且有cb較小時(shí)，對(duì)Rp，Rb及u的影響相對(duì)更顯著，當(dāng)cb增大到某一值后繼續(xù)增加，則對(duì)Rp，Rb及u的影響不明顯；cb恒定不變時(shí)，Sc的值越大，其對(duì)u的影響較小，而對(duì)Rb的影響越大，但隨Sc增大，Rb受其影響程度降低。圖7為φ軟化對(duì)Rp和Rb的影響，圖8為φ軟化對(duì)u的影響。分析圖7和圖8得，Sφ對(duì)Rp沒有影響，對(duì)u影響??；Sφ恒定不變時(shí)，Rp，Rb和u隨φb增大均呈現(xiàn)降低趨勢(shì)，且有φb越小，對(duì)Rp，Rb和u的影響程度越小，而φb越大，Rp和Rb兩者越接近；當(dāng)φb不變，Sφ越大，對(duì)Rb影響越大。

圖5 c軟化對(duì)Rp和Rb的影響Fig.5 Influence of cohesive force-softening on Rp and Rb

圖6 內(nèi)聚力軟化對(duì)位移的影響Fig.6 Influence of cohesive force-softening on displacement

圖7 φ軟化對(duì)Rp和Rb的影響Fig.7 Influence of φ on plastic and fractured zone radius

圖8 φ軟化對(duì)u的影響Fig.8 Influence of φ on displacement

3.3 破裂膨脹對(duì)Rp，Rb和u的影響

由式(27)和(31)可知，β1對(duì)Rp沒有影響，對(duì)Rb和u有影響；β2對(duì)Rp和Rb都沒有影響，但對(duì)u有影響。膨脹系數(shù)β1和β2對(duì)Rb及u的影響如圖9和圖10所示。

圖9 破裂膨脹對(duì)Rb的影響Fig.9 Influence of expansion on fractured zone radius

圖10 破裂膨脹對(duì)位移的影響Fig.10 Influence of expansion on displacement

分析圖9和圖10可以得出，β2恒定時(shí)，隨β1增大，Rb和u呈現(xiàn)逐漸增大趨勢(shì)，每當(dāng)β1增大0.5，破裂區(qū)半徑大約增大1%～2%，u增大7%～8%。β1恒定時(shí)，隨β2增大，變形u逐漸增大，且有β1越大，β2對(duì)變形u的影響越顯著。因此，膨脹對(duì)Rb影響不大，但對(duì)變形影響顯著。

3.4 支護(hù)阻力對(duì)Rp和Rb及位移u的影響

pi對(duì)Rp和Rb的影響見圖11，pi對(duì)u的影響見圖12。分析圖11可得，pi恒定時(shí)，隨p0增大，Rp和Rb均增大；p0不變時(shí)，隨pi增大，Rp和Rb均降低，當(dāng)pi增加0.5 MPa，Rp和Rb大約減少5%左右。分析圖12可知，巷道變形u隨pi增加逐漸減小，且p0越大，pi對(duì)巷道變形影響就變得越顯著。

圖11 pi對(duì)Rp和Rb的影響Fig.11 Influence of pi on plastic and fracture zone radius

圖12 pi對(duì)u的影響Fig.12 Influence of pi on displacement

4 工程應(yīng)用

潘一礦東區(qū)的-842 m東翼大巷位于二疊系上石盒子組含煤系地層中。該區(qū)域最大主應(yīng)力為34.67 MPa，巷道揭露巖性主要為砂質(zhì)泥巖、細(xì)砂巖、11-3煤及煤線等，圍巖單軸抗壓強(qiáng)度大都在20～40 MPa之間，屬于典型的高應(yīng)力軟巖巷道。巷道開挖后，圍巖強(qiáng)度不斷劣化，表現(xiàn)出明顯的應(yīng)變軟化行為。

文獻(xiàn)[15]在該巷道中選取2段進(jìn)行現(xiàn)場(chǎng)實(shí)驗(yàn)。2段均采用“錨桿+錨索+鐵絲網(wǎng)+混凝土”支護(hù)，但一段采用中空錨注錨桿，對(duì)圍巖注漿，減小圍巖軟化程度，另一段采用普通端錨錨桿，不進(jìn)行注漿，其他參數(shù)均相同。觀測(cè)結(jié)果發(fā)現(xiàn)，注漿段圍巖變形和破裂范圍明顯小于未注漿段。因此，降低圍巖軟化程度，能有效減小軟巖巷道的破裂范圍和變形。

5 結(jié)論

1)中間主應(yīng)力對(duì)巷道變形和破壞的影響具有區(qū)間性。塑性區(qū)、破裂區(qū)半徑及位移隨中間主應(yīng)力增加先減少后增大，且有|n-0.7|越大，影響程度越大。

2)巷道變形、塑性區(qū)和破裂區(qū)大小均隨軟化程度增大而增大；軟化模量對(duì)巷道變形影響不明顯，但對(duì)破裂區(qū)和塑性區(qū)大小受其影響比較顯著。

3)塑性區(qū)膨脹行為不影響塑性區(qū)半徑，但影響巷道變形和破裂區(qū)范圍；破裂區(qū)膨脹行為對(duì)破裂區(qū)和塑性區(qū)范圍不產(chǎn)生影響，但會(huì)影響巷道的變形。

4)支護(hù)阻力對(duì)圍巖破壞和變形影響明顯。當(dāng)支護(hù)提供的阻力每提高0.5 MPa，破裂區(qū)和塑性區(qū)的寬度縮小5%左右，且巷道所處地應(yīng)力越高，支護(hù)對(duì)巷道變形產(chǎn)生的作用效果就越顯著。