王章瓊，晏鄂川，劉毅學(xué)，王忠杰

（1.中國(guó)地質(zhì)大學(xué) 工程學(xué)院，武漢 430074；2.湖北省十房高速公路建設(shè)指揮部，湖北 十堰 442000）

1 引 言

鄂西北地區(qū)為國(guó)內(nèi)典型的變質(zhì)巖發(fā)育區(qū)，以武當(dāng)群片巖為主，區(qū)域地質(zhì)構(gòu)造較為復(fù)雜，片理面極其發(fā)育。十房高速公路通省隧道以深埋長(zhǎng)大成為該地區(qū)片巖公路隧道的控制性工程。隧道建設(shè)期間遇到大變形、縱向裂縫、鋼拱架被剪斷等問題，嚴(yán)重阻礙了隧道的正常建設(shè)。因此，研究片巖強(qiáng)度各向異性特征及其機(jī)制十分必要。

片巖屬典型的層狀巖體，其變形有著明顯的各向異性，這對(duì)工程巖體的穩(wěn)定性會(huì)產(chǎn)生特殊的影響[1－4]。王宏圖等[5]對(duì)層狀復(fù)合煤巖在三軸不等壓應(yīng)力狀態(tài)下的變形特性和強(qiáng)度特性進(jìn)行了較為系統(tǒng)的試驗(yàn)研究；劉卡丁等[6]提出了層狀巖體抗剪強(qiáng)度參數(shù)的經(jīng)驗(yàn)表達(dá)式，并通過解析法分析了層面與外荷載相對(duì)方位對(duì)剪切破壞面方位的影響；張玉軍等[7－8]針對(duì)層狀巖體抗剪強(qiáng)度具有方向性這一特點(diǎn)，提出了層狀巖體在不同方向上抗剪強(qiáng)度的表達(dá)式；高春玉等[9]通過室內(nèi)單軸和三軸試驗(yàn)，研究了某隧道圍巖中細(xì)微層理對(duì)巖石變形特性、強(qiáng)度特性及其參數(shù)的影響，結(jié)果表明，砂板巖力學(xué)特性的各向異性特征顯著；劉勝利等[10]通過室內(nèi)單軸壓縮和間接拉伸試驗(yàn)，研究了武當(dāng)群綠泥石片巖各向異性的力學(xué)特性。國(guó)外學(xué)者也對(duì)層狀巖體的變形、強(qiáng)度特征及其機(jī)制開展了大量研究[11－14]。

綜上所述，目前對(duì)片巖各向異性特征研究較多，但對(duì)各向異性機(jī)制研究較少，且已有研究多通過試樣在室內(nèi)試驗(yàn)后的變形破壞特征或數(shù)值模擬的方法進(jìn)行分析。上述研究?jī)H從整體、宏觀的角度研究層狀巖體力學(xué)特性和分析各向異性機(jī)制，從微觀角度進(jìn)行的研究較少。

本文通過對(duì)十房高速通省隧道片巖不同片理面傾角的室內(nèi)單軸壓縮試驗(yàn)，對(duì)比分析片巖變形參數(shù)的各向異性特征；在此基礎(chǔ)上，通過試樣宏觀破壞現(xiàn)象、試驗(yàn)前后片巖微觀結(jié)構(gòu)電鏡掃描以及數(shù)值模擬切片功能查看試樣內(nèi)部變形情況等3個(gè)層面和角度的分析，探討片巖變形參數(shù)各向異性的機(jī)制。

2 片巖變形參數(shù)特性試驗(yàn)研究

2.1 試樣制備

按不同片理面傾角加工試樣，加工后的試樣為標(biāo)準(zhǔn)圓柱樣，巖樣尺寸為50 mm× 100 mm?？紤]到片巖特殊的片理特性控制著結(jié)構(gòu)面的形態(tài)，為了能更好地反映隧道實(shí)際開挖過程中巖體的變形破壞情況，試樣按照?qǐng)D1所示片理面與試驗(yàn)時(shí)加載方向的相對(duì)關(guān)系進(jìn)行加工。圖中試樣a、b、c片理面傾角分別為 0°、45°、90°。

2.2 試驗(yàn)條件與方法

單軸壓縮試驗(yàn)采用INSTRON-1346型電液伺服巖石試驗(yàn)系統(tǒng)來完成。試驗(yàn)時(shí)按1 MPa/s的加荷速率加載圍壓，軸向加載速率為0.01 mm/s，直至試樣完全破壞。

2.3 應(yīng)力-應(yīng)變曲線

由圖2可知，片理面傾角為90°時(shí)，單軸極限抗壓強(qiáng)度和變形模量均較大；片理面傾角為0°時(shí)，單軸極限抗壓強(qiáng)度和變形模量略小于片理面傾角為90°時(shí)；片理面傾角為45°時(shí)，單軸極限抗壓強(qiáng)度和變形模量均較小。不同片理面傾角的片巖，其單軸極限抗壓強(qiáng)度和變形模量表現(xiàn)出明顯的各向異性特征。

圖2 試樣全應(yīng)力-應(yīng)變曲線Fig.2 Complete stress-strain curves of specimens

同時(shí)，不同片理面傾角的試樣，其應(yīng)力-應(yīng)變關(guān)系曲線特征也存在明顯的差異性。在孔隙裂隙壓密階段、彈性變形至微破裂穩(wěn)定發(fā)展階段和非穩(wěn)定破裂發(fā)展階段，曲線趨勢(shì)基本相同，但片理面傾角為0°和90°的試樣應(yīng)變量明顯大于片理面傾角為45°時(shí)，主要原因是，片理面傾角為45°時(shí)，單軸抗壓強(qiáng)度主要由片理面控制，試樣在較小的軸向力作用下即發(fā)生破壞。在峰值后階段，曲線的差異性更為明顯。片理面傾角為0°時(shí)，曲線出現(xiàn)多級(jí)震蕩，但整體趨勢(shì)較陡，表明其達(dá)到抗壓強(qiáng)度峰值后仍具有一定的承載力；片理面傾角為 45°時(shí)，峰值后應(yīng)力急劇降低或出現(xiàn)小幅度震蕩，曲線較陡；片理面傾角為90°時(shí)，峰值后曲線出現(xiàn)“臺(tái)階”，直至應(yīng)變?cè)隽窟_(dá)到0.001左右才開始跌落，其破壞表現(xiàn)出明顯的延性特征。不同片理面傾角的試樣，應(yīng)力-應(yīng)變曲線峰值后變化趨勢(shì)不同，這主要與各自的破壞機(jī)制有關(guān)，將在隨后的機(jī)制分析中作詳細(xì)論述。

3 片巖變形參數(shù)各向異性機(jī)制分析

3.1 變形參數(shù)分析

由單軸抗壓試驗(yàn)得到的彈性模量和泊松比等參數(shù)見表 1，其中，E1、E2、E3分別表示片理面傾角為 0°、45°、90°時(shí)試樣的彈性模量，v1、v2、v3分別表示片理面傾角為0°、45°、90°時(shí)試樣的泊松比。由表可知，片巖變形參數(shù)具有明顯的各向異性特征。按照彈性模量E由大到小的順序，試樣的片理面傾角依次為：90°、0°、45°，其比值E1: E2:E3=1.20:1.00:1.75；按照泊松比v由大到小的順序，試樣的片理面傾角依次為：45°、0°、90°，其比值v1:v2:v3= 0.91:1.00:0.73。

表1 片巖各向異性強(qiáng)度參數(shù)Table 1 Anisotropic strength parameters of schist

3.2 試樣變形破壞模式分析

不同片理面傾角的片巖單軸壓縮變形破壞型式如圖3所示。由試驗(yàn)結(jié)果可知，受片理面傾角影響，片巖表現(xiàn)出不同的破壞形式。

（1）片理面傾角為 0°時(shí)，試樣以剪斷破壞為主，同時(shí)在破壞面附近出現(xiàn)一些不規(guī)則的、沿軸向發(fā)展的剪斷面，試樣的破壞是片理面間巖石剪斷及微裂隙發(fā)展共同作用的結(jié)果。

（2）片理面傾角為45°時(shí)，試樣多沿片理面發(fā)生純剪切破壞，同時(shí)還存在少量次級(jí)剪切破壞面，表明試樣的變形破壞主要受片理面控制。

（3）片理面傾角為90°時(shí)，主要表現(xiàn)為沿軸向的劈裂破壞，局部出現(xiàn)錯(cuò)落擠壓，該變形破壞形態(tài)主要是片理面間巖石發(fā)生的“壓桿”失穩(wěn)破壞所導(dǎo)致。

上述不同片理面傾角試樣表現(xiàn)出不同變形破壞形式的主要原因是片理面法線方向與軸向力方向相對(duì)關(guān)系的不同。當(dāng)片理面與軸向力正交或平行時(shí)，片巖的變形主要受片理面間巖石控制；當(dāng)片理面與軸向力斜交時(shí)，片巖的變形主要受片理面控制。

圖3 試樣試驗(yàn)后破壞形態(tài)Fig.3 Typical failure mode of samples after test

3.3 片巖微觀結(jié)構(gòu)分析

為了從微觀角度分析片理面對(duì)片巖變形特性的影響及試驗(yàn)前后試樣內(nèi)部結(jié)構(gòu)的變化情況，采用電鏡掃描方法得到了不同片理面傾角試樣受力前后的內(nèi)部微觀結(jié)構(gòu)圖，見圖4。

由圖4可知，試驗(yàn)前片理面間巖石礦物呈定向排列，表面較平整、光滑，片理面張開度較??；試驗(yàn)后片理面張開度增大，新鮮破裂面上可見大量巖石碎屑及脫落的黏土礦物。片理面傾角為0°時(shí)，受力后由于片理面垂直方向受壓，片理面間巖石發(fā)生彎曲折斷破壞，并沿垂直片理面方向形成貫通裂隙，裂隙附近可見由于巖石斷裂產(chǎn)生的大量片狀碎屑。片理面傾角為45°時(shí)，受力后片理面間層狀巖石有明顯的錯(cuò)動(dòng)現(xiàn)象，且錯(cuò)動(dòng)方向基本與片理面平行，巖石表面可見云母受擠壓、摩擦后形成的粉末附著于片狀巖石上。片理面傾角為90°時(shí)，受力后片理面間層狀巖石產(chǎn)生了沿軸力方向分布的張拉裂隙，裂隙張開度較大，巖石表面光滑，無明顯擦痕。

圖4 試樣試驗(yàn)前后微觀形態(tài)對(duì)比Fig.4 Micromorphology of samples before and after tests

4 數(shù)值模擬驗(yàn)證

數(shù)值模擬方法已廣泛應(yīng)用到巖土體變形及力學(xué)特性的研究中[15]，并取得了較理想的效果。黃書嶺[4]、左宇軍[16]、王安明[17]、周斌[18]等采用數(shù)值試驗(yàn)方法，研究了層狀巖體的各向異性特性。為揭示巖試樣在壓縮試驗(yàn)過程中內(nèi)部應(yīng)力、應(yīng)變變化情況，本文采用FLAC3D進(jìn)行片巖單軸抗壓數(shù)值試驗(yàn)。

4.1 數(shù)值模擬模型及參數(shù)

根據(jù)室內(nèi)三軸試驗(yàn)試樣特點(diǎn)，建立了如圖5所示的片理面傾角分別為 0°、45°、90°的數(shù)值模擬模型，模型高為4 m、直徑為2 m，用等厚材料（10 cm）分別對(duì)巖石和片理面進(jìn)行等效處理。材料物理力學(xué)參數(shù)見表 2。試驗(yàn)中本構(gòu)模型采用應(yīng)變軟化模型，具體參數(shù)取值見表 3。試驗(yàn)采用位移加載的方式，速度為2 ′ 1 0-3mm/步。

圖5 數(shù)值試驗(yàn)巖樣Fig.5 Samples of numerical experiment

表2 計(jì)算參數(shù)Table 2 Calculation parameters

表3 應(yīng)變軟化參數(shù)Table 3 Strain softening parameters

4.2 數(shù)值模擬結(jié)果分析

為了形象描述模型內(nèi)部的信息，F(xiàn)LAC3D后處理中設(shè)置了切片操作功能[19]，即在模型中設(shè)置剖面，輸出剖面上的響應(yīng)結(jié)果?，F(xiàn)利用該功能對(duì)達(dá)到峰值強(qiáng)度且軸向力趨于穩(wěn)定后的各試樣進(jìn)行切片，剖面上軸向位移云圖如圖6所示。

由圖可知：

（1）片理面傾角為0°時(shí)，在平行于x軸的剖面上，位移云圖在沿軸向呈層狀分布的基礎(chǔ)上出現(xiàn)錯(cuò)動(dòng)（巖層壓斷），錯(cuò)動(dòng)帶在軸向上連續(xù)并形成兩個(gè)大的貫通面，貫通面與試樣軸向呈一定夾角，且沿軸向?qū)ΨQ。在平行于y軸的剖面上，位移云圖也出現(xiàn)一定的錯(cuò)動(dòng)，但規(guī)模相對(duì)較小，且未形成連續(xù)的貫通面。這表明片理面傾角為0°時(shí)，試樣的破壞形式為巖石折斷面貫通導(dǎo)致的剪切破壞。

（2）片理面傾角為45°時(shí)，在平行于x軸的剖面上，試樣內(nèi)部沿片理面發(fā)生明顯錯(cuò)動(dòng)。在平行于y軸的剖面上，僅上、下邊界附近局部位移較大，可能是受邊界效應(yīng)影響的結(jié)果。這表明，片理面傾角為45°時(shí)試樣主要發(fā)生沿片理面的剪切破壞，其破壞形態(tài)主要受片理面控制。

（3）片理面傾角為90°時(shí)，在平行于x軸的剖面上，位移沿軸向大體呈層狀分布，但沿徑向有一定幅度的錯(cuò)動(dòng)現(xiàn)象，且錯(cuò)動(dòng)規(guī)律在軸向上具有較好的一致性。在平行于y軸的剖面上（平行于片理面），位移的分布情況與平行于x軸的剖面上大體相同，但軸向位移云圖錯(cuò)動(dòng)幅度明顯偏小。這表明片理面傾角為90°時(shí)，試驗(yàn)后試樣主要發(fā)生平行片理面方向的劈裂破壞，其原因是，片理面強(qiáng)度遠(yuǎn)低于片理面間巖石，易沿片理面發(fā)生劈裂破壞。

對(duì)比室內(nèi)試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果可知，兩種試驗(yàn)條件下的試樣變形破壞模式基本相同，但數(shù)值試驗(yàn)結(jié)果中試樣的上、下端部變形情況與室內(nèi)試驗(yàn)結(jié)果存在一定差異，主要原因可能是數(shù)值試驗(yàn)時(shí)的邊界效應(yīng)所致。

5 各向異性機(jī)制分析

圖6 剖面上軸向位移云圖(單位: m)Fig.6 Contour of axial displacements in sections(unit: m)

通過上述片巖力學(xué)性質(zhì)及變形破壞模式的單軸壓縮試驗(yàn)及試驗(yàn)前、后片巖內(nèi)部微觀結(jié)構(gòu)特征的電鏡掃描以及單軸壓縮試驗(yàn)過程中試樣內(nèi)部變形特征數(shù)值模擬等3個(gè)層面的綜合分析，可初步歸納總結(jié)出片巖變形參數(shù)各向異性特征及機(jī)制，具體分析如下：

（1）片理面傾角為0°時(shí)，主要發(fā)生片理面間巖體的壓斷破壞，在宏觀上表現(xiàn)為剪切破壞，其剛度主要由束狀構(gòu)造的造巖礦物提供，片理面間巖體剛度得到了較好的發(fā)揮，因此，片巖彈性模量較大、泊松比較小。

（2）片理面傾角為45°時(shí)，片巖的變形主要由軟弱結(jié)構(gòu)面——片理面控制，片巖易發(fā)生沿片理面方向的剪切破壞，而片理面的抗變形能力遠(yuǎn)低于片理面間巖石。因此，片巖巖彈性模量最小、泊松比最大。

（3）片理面傾角為90°時(shí)，其力學(xué)效應(yīng)類似平行于片理面間巖石的壓桿失穩(wěn)，破壞形式主要表現(xiàn)為張拉劈裂，片理面間巖石的抗壓強(qiáng)度得到了充分發(fā)揮。因此，片巖巖彈性模量最大、泊松比最小。

6 結(jié) 論

（1）不同片理面傾角的片巖試樣，單軸壓縮試驗(yàn)應(yīng)力-應(yīng)變關(guān)系曲線特征存在明顯的差異。片理面傾角為45°時(shí)，峰值后曲線表現(xiàn)出明顯的脆性特征；片理面傾角為90°時(shí)，次之；片理面傾角為0°時(shí)，表現(xiàn)出一定的延性特征。

（2）片巖變形參數(shù)具有明顯的各向異性特征。按照彈性模量E由大到小的順序，試樣的片理面傾角依次為：90°、0°、45°，其比值 E1:E2:E3=1.20:1.00:1.75；按照泊松比 v由大到小的順序，試樣的片理面傾角依次為：45°、0°、90°，其比值 v1:v2:v3=0.91:1.00:0.73。

（3）綜合分析結(jié)果表明：片理面傾角為0°時(shí)，主要發(fā)生片理面間巖體的壓斷破壞，在宏觀上表現(xiàn)為剪切破壞；片理面傾角為45°時(shí)，控制片巖強(qiáng)度的主要因素是軟弱結(jié)構(gòu)面——片理面，片巖易發(fā)生沿片理面方向的剪切破壞；片理面傾角為90°時(shí)，其力學(xué)效應(yīng)類似平行于片理面間巖石的壓桿失穩(wěn)，破壞形式主要表現(xiàn)為張拉劈裂。

（4）片巖變形特性及破壞模式表現(xiàn)出明顯各向異性的根本原因是，片理面空間方位與軸向力方向不同的組合，使得片理面對(duì)片巖強(qiáng)度的貢獻(xiàn)差異明顯：片理面與軸向力斜交時(shí)，對(duì)片巖強(qiáng)度影響最大，平行時(shí)（90°）次之，垂直時(shí)（0°）最小。

[1]SERRANO A,OLALLA C.Ultimate bearing capacity of an anisotropic discontinuous rock mass,Part I: Basic modes of failure[J].International Journal of Rock Mechanics &Mining Sciences,1998,35(3): 301－324.

[2]黃書嶺,徐勁松,丁秀麗,等.考慮結(jié)構(gòu)面特性的層狀巖體復(fù)合材料模型與應(yīng)用研究[J].巖石力學(xué)與工程學(xué)報(bào),2010,29(4): 743－756.HUANG Shu-ling,XU Jin-song,DING Xiu-li,et al.Study of layered rock mss composite model based on characteristics of structural plane and its application[J].Chinese Journal of Rock Mechanics and Engineering,2010,29(4): 743－756.

[3]何忠明,彭振斌,曹平.層狀巖體單軸壓縮室內(nèi)試驗(yàn)分析與數(shù)值模擬[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2010,41(5): 1906－1912.HE Zhong-ming,PENG Zhen-bin,CAO Ping.Test and numerical simulation for stratified rock mass under uniaxial compression[J].Journal of Central South University(Science and Technology),2010,41(5): 1906－1912.

[4]黃書嶺,丁秀麗,鄔愛清,等.層狀巖體多節(jié)理本構(gòu)模型與試驗(yàn)驗(yàn)證[J].巖石力學(xué)與工程學(xué)報(bào),2012,31(8):1627－1635.HUANG Shu-ling,D1NG Xiu-li,WU Ai-qing,et al.Study of multi-joint constitutive model of layered rock mass and experimental verification[J].Chinese Journal of Rock Mechanics and Engineering,2012,31(8):1627－1635.

[5]王宏圖,鮮學(xué)福,賀建民.層狀復(fù)合煤巖的三軸力學(xué)特性研究[J].礦山壓力與頂板管理,1999,(1): 81－83,88.WANG Hong-tu,XIAN Xue-fu,HE Jian-min.Study of the triaxial mechanical properties of a stratified coal-rock union[J].Journal of Mine Pressure and Strata Control,1999,(1): 81－83,88.

[6]劉卡丁,張玉軍.層狀巖體剪切破壞面方向的影響因素[J].巖石力學(xué)與工程學(xué)報(bào),2002,21(3): 335－339.LIU Ka-ding,ZHANG Yu-jun.Influence factors on shear failure orientation of layered rocks[J].Chinese Journal of Rock Mechanics and Engineering,2002,21(3): 335－339.

[7]張玉軍,劉誼平.層狀巖體抗剪強(qiáng)度的方向性及剪切破壞面的確定[J].巖土力學(xué),2001,22(3): 254－257.ZHANG Yu-jun,LIU Yi-ping.Anisotropy of shear strength of layered rocks and deformation of shear failure plane[J].Rock and Soil Mechanics,2001,22(3): 254－257.

[8]張玉軍,劉誼平.層狀巖體的三維彈塑性有限元分析[J].巖石力學(xué)與工程學(xué)報(bào),2002,21(11): 1615－1619.ZHANG Yu-jun,LIU Yi-ping.3D elastoplastic FEM analysis on layered rockmass[J].Chinese Journal of Rock Mechanics and Engineering,2002,21(11): 1615－1619.

[9]高春玉,徐進(jìn),李忠洪,等.雪峰山隧道砂板巖各向異性力學(xué)特性的試驗(yàn)研究[J].巖土力學(xué),2011,32(5):1360－1364.GAO Chun-yu,XU Jin,LI Zhong-hong,et al.Experimental study of anisotropically mechanical characteristics of sandy slate in Xuefeng mountain tunnel[J].Rock and Soil Mechanics,2011,32(5): 1360－1364.

[10]劉勝利,陳善雄,余飛,等.綠泥石片巖各向異性特性研究[J].巖土力學(xué),2012,33(12): 3616－3623.LIU Sheng-li,CHEN Shan-xiong,YU Fei,et al.Anisotropic properties study of chlorite schist[J].Rock and Soil Mechanics,2012,33(12): 3616－3623.

[11]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.

[12]SITHARAM T G,SRIDEVI J,SHIMIZU N.Practical equivalent continuum characterization of jointed rock masses[J].International Journal of Rock Mechanics and Mining Sciences,2001,38(3): 437－448.

[13]NOVA R.The failure of transversely isotropic rocks in triaxial compression[J].International Journal of Rock Mechanics and Mining Sciences,1980,17(6): 325－332.

[14]JAEGER J C.Shear failure of anisotropic rocks[J].Geology Magazine,1960,97(1): 65－72.

[15]唐輝明,晏鄂川,胡新麗.工程地質(zhì)數(shù)值模擬的理論與方法[M].武漢: 中國(guó)地質(zhì)大學(xué)出版社,2001.

[16]左宇軍,唐春安,朱萬成,等.動(dòng)載荷作用下巖石破壞過程的數(shù)值試驗(yàn)研究[J].巖土力學(xué),2008,29(4): 887－892.ZUO Yu-jun,TANG Chun-an,ZHU Wan-cheng,et al.Numerical test on failure process of rock subjected to dynamic loading[J].Rock and Soil Mechanics,2008,29(4): 887－892.

[17]王安明,楊春和,黃誠(chéng),等.層狀鹽巖力學(xué)和變形特性數(shù)值試驗(yàn)研究[J].巖土力學(xué),2009,30(7): 2173－2178.WANG An-ming,YANG Chun-he,HUANG Cheng,et al.Numerical experiment study of deformation and mechanical properties of layered salt rock[J].Rock and Soil Mechanics,2009,30(7): 2173－2178.

[18]周斌,張可能,柳群義.考慮應(yīng)變軟化特征的層狀巖體三軸壓縮數(shù)值試驗(yàn)分析[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,40(3): 828－832.ZHOU Bin,ZHANG Ke-neng,LIU Qun-yi.Numerical test analysis for stratified rock mass under triaxial compression considering strain soften characteristic[J].Journal of Central South University(Science and Technology),2009,40(3): 828－832.

[19]陳育民,徐鼎平.FLAC/FLAC3D基礎(chǔ)與工程實(shí)例[M].北京: 中國(guó)水利水電出版社,2009.