王 磊，王渭明

(山東科技大學(xué)a.土木建筑學(xué)院;b.山東省土木工程防災(zāi)減災(zāi)重點(diǎn)實(shí)驗(yàn)室，山東 青島 266590)

1 研究背景

弱膠結(jié)軟巖廣泛分布于西部煤系地層，由于其特定的沉積環(huán)境，具有明顯的層狀結(jié)構(gòu)［1－3］。從巖體力學(xué)角度分析，這類巖層表現(xiàn)為明顯的各向異性和塑性變形特征。在這類地層中進(jìn)行巷道開挖支護(hù)，圍巖應(yīng)力場分布規(guī)律、變形特征及破壞模式同一般巖體有明顯差異。

如何在巖體工程研究中考慮材料的各向異性對(duì)應(yīng)力場、位移場以及破壞模式的影響，國內(nèi)外學(xué)者進(jìn)行了大量的研究。Yoshinaka［4－5］研究了4種不同軟巖的變形特征，分析了軟巖的強(qiáng)度參數(shù)對(duì)應(yīng)力分布和變形特征的影響;Nasseri等［6－8］通過室內(nèi)試驗(yàn)研究了變質(zhì)巖的單軸壓縮與三軸壓縮變形特征，并提出了考慮各向異性的本構(gòu)模型;李鏡培［9］以Maxwell黏性模型推導(dǎo)了軟巖的各向異性黏彈塑性本構(gòu)關(guān)系，并分析了模型各參數(shù)對(duì)軟巖變形分布規(guī)律的影響;Adachi［10］以室內(nèi)試驗(yàn)為基礎(chǔ)分析軟巖的應(yīng)力應(yīng)變特征，建立了考慮軟巖應(yīng)變硬化和應(yīng)變軟化的彈塑性本構(gòu)關(guān)系;Mikenina等［11］從微觀角度分析了材料各向異性彈塑性變形特征，并進(jìn)行了各向異性彈塑性本構(gòu)關(guān)系二次開發(fā)，分析了軟巖隧道的變形分布特征;李曉紅等［12］根據(jù)層狀巖體中隧道施工變形規(guī)律，以橫觀各向異性彈塑性模型分析了深埋隧道施工變形特征;王渭明等［13－14］通過將各向異性材料分層計(jì)算，分析了凍結(jié)軟巖的黏塑性變形。

上述研究成果從不同角度研究了材料各向異性對(duì)其變形特征的影響，本文嘗試以空間滑動(dòng)面理論為基礎(chǔ)，通過室內(nèi)試驗(yàn)確定模型參數(shù)，推導(dǎo)基于統(tǒng)一強(qiáng)度理論的各向異性彈塑性本構(gòu)模型，并通過FLAC3D二次開發(fā)，分析軟巖巷道變形分布及其支護(hù)控制技術(shù)。

2 空間滑動(dòng)面理論

空間滑動(dòng)面理論(SMP理論)考慮了中間主應(yīng)力的影響(如圖1所示)，能較好地反映巖土材料的實(shí)際變形狀態(tài)［15］。張連衛(wèi)［16］、羅汀［17］、鄒博［18］、連鎮(zhèn)營［19］等分別基于空間滑動(dòng)面理論研究了材料的破壞模式。

圖1 π平面上SMP準(zhǔn)則Fig.1 SMP criterion in π plane

考慮內(nèi)黏聚力的空間滑動(dòng)面準(zhǔn)則可以表示為

式中:σ1，σ2，σ3為第一、第二、第三主應(yīng)力;μ 為材料摩擦系數(shù)，μ=tanφ，φ為材料的摩擦角，c為材料黏聚力。

3 基于SMP的橫觀各向異性模型

3.1 橫觀各向異性彈性本構(gòu)關(guān)系

在如圖2所示的橫觀各向異性體中，材料的彈性本構(gòu)關(guān)系可以表示為

式中aij為柔度系數(shù)，可以由工程彈性常數(shù)表示:

式中:E1，μ1為各向同性面內(nèi)的彈性模量和泊松比;E2，μ2為垂直于各向同性面內(nèi)的彈性模量和泊松比;G2為垂直于各向同性面內(nèi)的剪切模量。式(3)中各參數(shù)均可由室內(nèi)試驗(yàn)獲取。

圖2 橫觀各向異性材料Fig.2 Transversely anisotropic materials

3.2 基于SMP的各向異性模型推導(dǎo)

由式(1)可構(gòu)造屈服函數(shù)如下:

由式(4)，根據(jù)FLAC3D中增量算法，非關(guān)聯(lián)流動(dòng)法則條件下的塑性勢函數(shù)可以表示為

材料的應(yīng)變增量為

各向異性條件下，由廣義胡克定律應(yīng)力增量可以表示為

式中矩陣Si表示彈性條件下應(yīng)力應(yīng)變關(guān)系。

根據(jù)塑性力學(xué)理論，塑性應(yīng)變增量可以表示為

式中λ為屈服塑性因子。

由式(6)至式(8)可得在塑性條件下應(yīng)力增量為

新的應(yīng)力狀態(tài)仍然在塑性屈服面上，即滿足

將式(9)代入式(10)可求得塑性因子λ。

根據(jù)FLAC3D理論中的彈性應(yīng)力假設(shè)和彈性應(yīng)力猜想，任一時(shí)步的應(yīng)力可以表示為:

由式(10)至式(12)可得新的應(yīng)力狀態(tài)為

由式(13)在各向異性條件下根據(jù)前一時(shí)步的應(yīng)力狀態(tài)求得當(dāng)前時(shí)步的應(yīng)力狀態(tài)。

以上為材料受壓屈服推導(dǎo)，在受拉條件下，材料的屈服函數(shù)為

式中σt為材料的抗拉強(qiáng)度。

3.3 FLAC3D二次開發(fā)流程

根據(jù)上述推導(dǎo)關(guān)系，在橫觀各向異性條件下FLAC3D二次開發(fā)基本流程為如圖3所示。

圖3 FLAC3D二次開發(fā)流程Fig.3 Process of secondary development in FLAC3D

根據(jù)圖3，各向異性模型二次開發(fā)主要包括函數(shù)Ini()和ANSMP()的編譯，可通過C++2005編程調(diào)試并形成dll文件，直接配置到FLAC3D中即可。

4 模型分析

4.1 橫觀各向異性參數(shù)室內(nèi)測試

目前測定巖體各向異性參數(shù)的常用方法是不同傾角下單軸壓縮試驗(yàn)，根據(jù)文獻(xiàn)Talesnick［20］和Tien Yong Ming［21－22］等研究結(jié)果，可取 β =0°，β =45°，β=90°三組傾角巖體試件進(jìn)行室內(nèi)單軸壓縮試驗(yàn)確定巖體的橫觀各向異性參數(shù)，如圖4所示。

圖4 取樣與試件Fig.4 The rock samples

根據(jù)測試結(jié)果，巖體試件的各向異性參數(shù)分別為

4.2 模型驗(yàn)證分析

結(jié)合內(nèi)蒙某礦井底車場大巷，建立三維數(shù)值模型如圖5所示。

圖5 三維數(shù)值模型Fig.5 Three-dimensional numerical model

深度范圍內(nèi)巖層分2層，采用本文基于SMP的橫觀各向異性模型，各向異性參數(shù)根據(jù)試驗(yàn)結(jié)果和現(xiàn)場變形測試進(jìn)行了修正，各層力學(xué)參數(shù)取值如表1所示。

表1 巖層1、巖層2力學(xué)參數(shù)Table 1 Mechanical parameters of rock layer 1 and rock layer 2

巷道開挖模擬塑性區(qū)分布如圖6(a)所示。從圖6(a)中可以看出，由于各向異性的影響，巷道開挖后塑性區(qū)分布呈現(xiàn)非對(duì)稱性，而在各向同性條件下，應(yīng)力分布于塑性區(qū)分布是對(duì)稱的。

當(dāng)材料在各個(gè)方向的彈性模量和泊松比分別相同時(shí)，模型退化為各向同性彈塑性本構(gòu)關(guān)系，此時(shí)塑性區(qū)分布應(yīng)對(duì)稱。為進(jìn)一步驗(yàn)證模型的可靠性，設(shè)計(jì)數(shù)值模擬試驗(yàn)，材料參數(shù)取值如表2所示。

表2 力學(xué)參數(shù)Table 2 Mechanical parameters

模擬計(jì)算塑性區(qū)分布如圖6(b)所示。

圖6 塑性區(qū)分布Fig.6 Distribution of plastic zone of anisotropic and isotropic materials in the simulation

從圖6(b)中可以看出，將各個(gè)方向材料參數(shù)取值相同后，塑性區(qū)分布大致對(duì)稱，符合各向同性彈塑性本構(gòu)關(guān)系下塑性區(qū)分布規(guī)律，因此模型是可靠的。

4.3 巷道施工支護(hù)分析

由于巖體分布的各向異性，使得巷道開挖后最大應(yīng)力集中區(qū)發(fā)生變化，而傳統(tǒng)分析模型塑性區(qū)對(duì)稱出現(xiàn)在拱肩部，對(duì)兩肩同時(shí)加強(qiáng)支護(hù)雖然也可以有效控制施工變形，但是造成極大的材料浪費(fèi)。

經(jīng)模擬計(jì)算，原設(shè)計(jì)錨桿受力如圖7(a)所示。

由圖7(a)可以看出，由于地層各向異性分布特征，使得錨桿受力呈現(xiàn)明顯偏壓狀態(tài)，需要修改錨桿支護(hù)參數(shù)。

根據(jù)塑性區(qū)分布狀態(tài)，將右側(cè)肩部錨桿長度改為3.5 m，其他參數(shù)不變，經(jīng)模擬計(jì)算錨桿受力如圖7(b)所示。

圖7 錨桿受力分布Fig.7 Distribution of anchor force in the original design and after the optimization

由圖7可知，經(jīng)塑性區(qū)錨桿長度修正后，圍巖塑性區(qū)明顯減小，且錨桿受力趨向于均勻。

經(jīng)該礦現(xiàn)場應(yīng)用，優(yōu)化支護(hù)方案后有效控制了巷道表面變形，噴層也沒有出現(xiàn)開裂現(xiàn)象。

5 結(jié)論

(1)在各向異性條件下，巷道開挖圍巖破壞區(qū)分布呈現(xiàn)非對(duì)稱性，而應(yīng)用各向同性模型模擬計(jì)算時(shí)并不能反映巷道施工非對(duì)稱破壞的現(xiàn)象。

(2)根據(jù)將各向異性參數(shù)在各主向上取相同值，模型退化為各向同性彈塑性模型，圍巖塑性區(qū)分布呈現(xiàn)對(duì)稱分布的特征，設(shè)計(jì)模擬方案，驗(yàn)證了模型的可靠性。

(3)結(jié)合實(shí)際巷道監(jiān)測進(jìn)行了橫觀各向異性條件下錨桿支護(hù)分析，根據(jù)模擬結(jié)果，錨桿受力呈現(xiàn)明顯不對(duì)稱。

(4)經(jīng)修改錨桿支護(hù)參數(shù)，將破壞區(qū)錨桿加長，根據(jù)模擬結(jié)果，圍巖塑性區(qū)分布明顯減小，且錨桿受力趨向均勻。

(5)根據(jù)巷道施工現(xiàn)場應(yīng)用，基于各向異性的關(guān)鍵區(qū)錨桿支護(hù)對(duì)控制圍巖變形是有效的。

［1］周應(yīng)麟，邱喜華.層狀巖層圍巖隧道穩(wěn)定性的探討［J］.地下空間與工程學(xué)報(bào)，2006，2(2):345－348.(ZHOU Ying-lin，QIU Xi-hua.Study on Stability of Stratified Rock Tunnel［J］.Chinese Journal of Underground Space and Engineering，2006，2(2):345－348.(in Chinese))

［2］孫 鈞.海底隧道工程設(shè)計(jì)施工若干關(guān)鍵技術(shù)的商榷［J］.巖石力學(xué)與工程學(xué)報(bào)，2006，25(8):1514－1522.(SUN Jun.Discussion on Some Key Technical Issues for Design and Construction of Undersea Tunnels［J］.Chinese Journal of Rock Mechanics and Engineering，2006，25(8):1514 －1522.(in Chinese))

［3］張坤勇，殷宗澤，梅國雄.土體各向異性研究進(jìn)展［J］.巖土力學(xué)，2004，25(9):1503－1509.(ZHANG Kun-yong，YIN Zong-ze，MEI Guo-xiong.Development of Soil’s Anisotropy Study ［J］.Rock and Soil Mechanics，2004，25(9):1503 －1509.(in Chinese))

［4］YOSHINAKA R，TRAN T V，OSADA M.Non-linear Stress and Strain-Dependent Behaviour of Soft Rocks under Cyclic Triaxial Conditions［J］.International Journal of Rock Mechanics and Mining Sciences，1998，35(7):941－955.

［5］YOSHINAKA R，TRAN T V，OSADA M.Deformation Behaviour of Soft Rocks During Consolidated Undrained Cyclic Triaxial Testing［J］.International Journal of Rock Mechanics and Mining Sciences，1998，33(6):557 －572.

［6］NASSERI M H B，RAO K S，RAMAMURTHY T.Anisotropic Strength and Deformational Behavior of Himalayan Schists［J］.International Journal of Rock Mechanics and Mining Sciences，2003，40(1):3 －23.

［7］TALESNICK M L，BLOCH-FRIEDMAN E A.Compatibility of Different Methodologies for the Determination of Elastic Parameters of Intact Anisotropic Rocks［J］.International Journal of Rock Mechanics and Mining Sciences，1999，36:919－940.

［8］KRIVOI A F，POPOV G Y.Features of the Stress Field near Tunnel Inclusions in an Inhomogeneous Anisotropic Space［J］.International Applied Mechanics，2008，44(6):626－634.

［9］李鏡培，饒平平，張述濤.各向異性黏彈性隧道開挖位移時(shí)變解析解［J］.哈爾濱工業(yè)大學(xué)學(xué)報(bào)，2011，43(4):94 －98.(LI Jing-pei，RAO Ping-ping，ZHANG Shutao.Time-Varying Displacement Analytical Solution of Anisotropic Visco-Elastic Tunnel Excavation［J］.Journal of Harbin Institute of Technology，2011，43(4):94 －98.(in Chinese))

［10］ADACHI T，OKA F.An Elasto-Plastic Constitutive Model for Soft Rock with Strain Softening［J］.International Journal for Numerical and Analytical Methods in Geomechanics，1995，19(4):233 －247.

［11］LAVRIKOV S V，MIKENINA O A，REVUZHENKO A F.A Non-Archimedean Number System to Characterize the Structurally Inhomogeneous Rock Behavior Nearby a Tunnel［J］.Journal of Rock Mechanics and Geotechnical Engineering，2011，3(2):153－160.

［12］李曉紅，夏彬偉，李丹，等.深埋隧道層狀圍巖變形特征分析［J］.巖土力學(xué)，2010，31(4):1163－1167.(LI Xiao-hong，XIA Bin-wei，LI Dan，et al.Deformation Characteristics Analysis of Layered Rockmass in Deep Buried Tunnel［J］.Rock and Soil Mechanics，2010，31(4):1163 －1167.(in Chinese))

［13］王渭明，路林海.正交各向異性復(fù)合井壁應(yīng)力變形分析與應(yīng)用［J］.力學(xué)與實(shí)踐，2009，31(1):52－56.(WANG Wei-ming，LU Lin-hai.Analysis of Stresses and Deformations of a Compound and Orthotropic Shaft with Applications［J］.Mechanics in Engineering，2009，31(1):52 －56.(in Chinese))

［14］王渭明，王 磊，代春泉.基于強(qiáng)度分層計(jì)算的弱膠結(jié)軟巖凍結(jié)壁變形分析［J］.巖石力學(xué)與工程學(xué)報(bào)，2011，30(增 2):4110 － 4116.(WANG Wei-ming，WANG Lei，DAI Chun-quan.Frozen Wall Deformation Analysis in Weakly Cemented Soft Rock Based on Layered Calculation of Strength［J］.Chinese Journal of Rock Mechanics and Engineering，2011，30(Sup.2):4110 －4116.(in Chinese))

［15］MTSUOKA H，HOSHIKAWA T，UENO K.A General Failure Criterion and Stress-Strain Relation for Granular Materials to Metals［J］.Soils and Foundations，1990，30(2):119－127.

［16］張連衛(wèi)，張建民，張嘎.基于SMP的粒狀材料各向異性強(qiáng)度準(zhǔn)則［J］.巖土工程學(xué)報(bào)，2008，30(8):1107－1111.(ZHANG Lian-wei，ZHANG Jian-min，ZHANG Ga.SMP-Based Anisotropic Strength Criteria of Granular Materials［J］.Chinese Journal of Geotechnical Engineering，2008，30(8):1107 －1111.(in Chinese))

［17］羅 汀，李 萌，孔玉俠，等.基于SMP的巖土各向異性強(qiáng)度準(zhǔn)則［J］.巖土力學(xué)，2009，30(增2):127－131.(LUO Ting，LI Meng，KONG Yu-xia，et al.Failure Criterion Based on SMP for Anisotropic Geomaterials［J］.Rock and Soil Mechanics，2009，30(Sup.2):127 －131.(in Chinese))

［18］鄒 博，姚仰平，路德春.變換應(yīng)力三維化方法在清華模型中的應(yīng)用［J］.巖石力學(xué)與工程學(xué)報(bào)，2005，24(23):4303 －4307.(ZOU Bo，YAO Yang-ping，LU Dechun.Qinghua Model Revised by SMP Criterion［J］.Chinese Journal of Rock Mechanics and Engineering，2005，24(23):4303 －4307.(in Chinese))

［19］連鎮(zhèn)營，韓國城，姚仰平.基于SMP準(zhǔn)則的改進(jìn)劍橋模型及其在基坑工程中應(yīng)用［J］.大連理工大學(xué)學(xué)報(bào)，2002，42(1):93 － 97.(LIAN Zhen-ying，HAN Guocheng，YAO Yang-ping.Revised Cam-clay Model Based on SMP Criterion and Its Application to Foundation Engineering［J］.Journal of Dalian University of Technology，2002，42(1):93 －97.(in Chinese))

［20］HOULSBY G T.A General Failure Criterion for Frictional and Cohesive Materials［J］.Soils and Foundations，1986，26(2):97－101.

［21］TIEN Y M，KUO M C，JUANG C H.An Experimental Investigation of the Failure Mechanism of Simulated Transversely Isotropic Rocks［J］.International Journal of Rock Mechanics and Mining Sciences，2006，43(8):1163－1181.

［22］PIETRUSZCZAKA S，LYDZBA D，SHAO J F.Modelling of Inherent Anisotropy in Sedimentary Rocks［J］.International Journal of Solids and Structures，2002，39(3):637－648.