楊宏印，來家玉，肖尊群,，王 鑫

(1.武漢工程大學土木工程與建筑學院，湖北 武漢 430074；2.武漢工程大學興發(fā)礦業(yè)學院，湖北 武漢 430074)

0 引 言

石英砂巖的碎屑物質(zhì)中90%以上為石英(可包括燧石和硅質(zhì)巖)碎屑，可有少數(shù)長石、巖屑等，重礦物很少，常為穩(wěn)定的、磨圓的重礦物。微風化石英砂巖是典型硬巖代表，宜巴高速公路穿過的主要區(qū)域的微風化石英砂巖巖樣青色、質(zhì)地堅硬，節(jié)理裂隙不發(fā)育，這樣的巖體形成的巖質(zhì)邊坡穩(wěn)定性好，邊坡巖體通常處于自然穩(wěn)定狀態(tài)。為了探究這一類巖石，目前常用的手段是通過室內(nèi)試驗獲得宏觀力學特性[1-3]。不同圍壓下巖樣的雙軸壓縮試驗可以從一定程度上反映不同深度處邊坡巖體強度特性，并獲得巖體邊坡穩(wěn)定評價重要參數(shù)黏聚力c和內(nèi)摩擦角Φ[4-6]，因此，雙軸試驗結(jié)果對于評價巖質(zhì)邊坡穩(wěn)定具有十分重要的意義。除此之外，該類巖質(zhì)邊坡在降雨條件下是否還能保持穩(wěn)定狀態(tài)，還需要進行飽和試樣的雙軸試驗進一步研究。宏觀物理參數(shù)只能反應(yīng)巖體某個特定的狀態(tài)，而邊坡失穩(wěn)是一個逐步演化的過程，要得到這個細觀演化過程，目前還沒有很好的室內(nèi)試驗可以做到。基于顆粒流的離散單元法可以實現(xiàn)對試樣細觀演化過程的模擬。顆粒流理論的PFC2D離散元程序是研究巖樣壓縮試驗微觀組構(gòu)演化規(guī)律的重要手段[7-9]。

顆粒間平行黏結(jié)接觸模型能傳遞力和力矩，被廣泛用于模擬巖石內(nèi)部顆粒相互作用。HAZZARD等[10]采用接觸黏結(jié)模型分析了花崗巖的破壞特征；余華中等[11]對由lump組成的大理巖進行雙軸試驗研究；周杰等[12]采用平行黏結(jié)模型研究砂巖三軸破壞演化宏細觀機理；李守巨等[13]對非均質(zhì)巖石-中晶大理巖單軸壓縮試驗破壞過程細觀模擬。

本文以微風化石英砂巖巖樣在飽和、天然兩種狀態(tài)下的雙軸壓縮試驗宏觀力學響應(yīng)數(shù)據(jù)為基礎(chǔ)，建立雙軸壓縮試驗PFC2D數(shù)值模型，重點研究不同浸水條件、不同圍壓下試樣法向、切向接觸力、配位數(shù)等主要微觀組構(gòu)參數(shù)演化規(guī)律的影響，探尋巴東組石英砂巖雙軸壓縮條件下破壞的細觀機理，進而獲得石英砂巖巖質(zhì)邊坡失穩(wěn)細觀演化機制。

1 數(shù)值試驗設(shè)計

1.1 數(shù)值試樣的制備

進行顆粒流模擬試驗時，首先生成試樣模型墻體將顆粒包圍，然后通過設(shè)置孔隙率，形成給定孔隙率的試樣，將試樣所有顆粒速度清零，得到初始試樣。初始試樣需經(jīng)伺服達到指定圍壓后才能進行加載，伺服圍壓與室內(nèi)試驗圍壓完全一致，分為3 MPa、6 MPa、9 MPa、12 MPa四個圍壓。試樣形成后，對試樣顆粒設(shè)置平行黏結(jié)接觸，對上下加載板施加速度實現(xiàn)試樣的加載過程，并記錄整個試驗過程中墻體的位移、不平衡力，通過后處理得到試件的宏觀變形曲線。顆粒流模型尺寸與室內(nèi)試驗試件尺寸一致，試樣寬50 mm，高100 mm，顆粒半徑取0.05～0.07 mm，初始孔隙率為0.15，密度2.64×103kg/m3。

1.2 細觀參數(shù)標定

不同圍壓下的石英砂巖試樣宏觀物理力學參數(shù)見表1。需要指出的是，泊松比v采用對應(yīng)試樣的單軸抗壓強度試驗結(jié)果。根據(jù)宏觀力學參數(shù)標定得到試樣平行黏結(jié)細觀參數(shù)見表2。

1.3 破壞模式控制

圖1 飽和石英砂巖試樣摩爾應(yīng)力圓對比Fig.1 Comparison of molar stress circle of saturated quartz sandstone sample

圖2 天然石英砂巖試樣摩爾應(yīng)力圓對比Fig.2 Comparison of molar stress circle of natural quartz sandstone sample

表1 石英砂巖室內(nèi)雙軸試驗結(jié)果Table 1 Quartz sandstone indoor biaxial test results

表2 石英砂巖顆粒接觸細觀參數(shù)標定結(jié)果Table 2 Calibration results of quartz sandstone particle contact mesoscopic parameters

圖3 飽和石英砂巖試樣破壞對比Fig.3 Destruction comparison of saturated quartz sandstone samples

圖4 天然石英砂巖試樣破壞對比Fig.4 Destruction comparison of natural quartz sandstone samples

圖5 飽和石英砂巖應(yīng)力應(yīng)變曲線對比Fig.5 Comparison of stress-strain curves of saturated quartz sandstone

圖6 天然石英砂巖應(yīng)力應(yīng)變曲線對比Fig.6 Comparison of stress-strain curves of natural quartz sandstone

1.4 試驗結(jié)果對比

將校正后的細觀物理力學性質(zhì)參數(shù)用于PFC2D顆粒流模型中，進行雙軸壓縮試驗?zāi)M。圖5和圖6分別為飽和試樣和天然試樣PFC2D模型與室內(nèi)雙軸壓縮試驗應(yīng)力-應(yīng)變對比曲線。由圖5和圖6可以看出，不同圍壓下的巖樣對應(yīng)的峰值應(yīng)力差異較大，峰值應(yīng)力與圍壓呈正比，圍壓越大，峰值應(yīng)力越大。不同含水條件下對應(yīng)的峰值應(yīng)力也有一定的差異，相同圍壓下天然試樣的峰值應(yīng)力較飽和試樣高。無論是飽和試樣還是天然試樣，雙軸壓縮試驗的初期，都存在非線性變形階段。為了模擬該階段，在加載初期設(shè)置線性模量和平行黏結(jié)模量隨軸向應(yīng)變變化函數(shù)，通過Fish語言程序?qū)崿F(xiàn)。室內(nèi)試驗與數(shù)值模擬在達到峰值前，應(yīng)力-應(yīng)變曲線的擬合度很好。

2 石英砂巖的微觀組構(gòu)分析

2.1 不同圍壓下顆粒間接觸力對比

以天然石英砂巖為例，繪制不同圍壓下巖樣達到破壞時法向、切向接觸力分布玫瑰花圖如圖7所示。在圖7中，徑向極軸表示某個方向統(tǒng)計范圍內(nèi)所有接觸的平均接觸力，N；環(huán)向為統(tǒng)計范圍的方位角，(°)。

由圖7(a)可以看出，不同統(tǒng)計范圍，圍壓對法向接觸力的影響不同。對于試樣的上下兩端，法向接觸力由大到小排列順序依次為12 MPa、9 MPa、6 MPa、3 MPa，這與不同圍壓下試樣的峰值應(yīng)力順序一致，圍壓越大，試樣破壞時峰值應(yīng)力越大，試樣上下兩端的法向接觸力越大，表明圍壓對試樣破壞時上下兩端的法向接觸力影響較小。對于試樣的左右兩側(cè)，法向接觸力大小呈現(xiàn)無規(guī)律性，表明圍壓對試樣破壞時左右兩側(cè)的法向接觸力影響較大。

由圖7(b)可以看出，圍壓不能改變顆粒間切向接觸力的分布范圍，只能改變接觸力的大小。圍壓為12 MPa時，切向接觸力最大，其次是9 MPa、6 MPa和3 MPa，這與不同圍壓下試樣的峰值應(yīng)力順序一致，圍壓越大，試樣破壞時峰值應(yīng)力越大，試樣的切向接觸力越大，表明圍壓對試樣破壞時切向接觸力影響較小。

2.2 不同狀態(tài)試樣接觸力對比分析

以3 MPa圍壓為例，繪制試樣達到破壞時天然和飽和石英砂巖法向、切向接觸力統(tǒng)計玫瑰花圖如圖8所示。相同圍壓下，天然狀態(tài)下石英砂巖顆粒接觸力比飽和狀態(tài)下大，宏觀表現(xiàn)為天然石英砂巖抗壓強度比飽和石英砂巖高。試樣浸水以后，宏觀表現(xiàn)上，試樣雙軸抗壓強度發(fā)生損失，從細觀角度看，可以理解為浸水條件降低了組成巖樣顆粒之間的平行黏結(jié)接觸參數(shù)。浸水條件僅能改變試樣顆粒接觸力大小，不能改變顆粒接觸力分布范圍。

圖7 天然石英砂巖接觸力分布Fig.7 Contact force distribution of natural quartz sandstone

圖8 石英砂巖不同狀態(tài)下接觸力分布對比Fig.8 Comparison of contact force distribution in different states of quartz sandstone

2.3 配位數(shù)演化

配位數(shù)為單個顆粒平均接觸數(shù)，配位數(shù)的演化關(guān)系直接反映壓縮過程中顆粒之間相對移動程度。圖9為天然石英砂巖配位數(shù)隨時間、裂隙數(shù)的演化關(guān)系。由圖9可知，壓縮前期，配位數(shù)隨時間呈水平曲線，隨著壓縮的進行，試樣內(nèi)部產(chǎn)生裂隙，配位數(shù)緩慢減小，隨著加載的繼續(xù)，裂隙數(shù)量繼續(xù)增多，同時原有的微裂隙進一步擴展，配位數(shù)急劇減小，最終達到相對穩(wěn)定的狀態(tài)。不同圍壓下配位數(shù)進入下降階段的時間點不同，沒有一定的規(guī)律，下降階段的曲線斜率基本相同，說明圍壓對配位數(shù)進入下降階段的時間點影響很大，對配位數(shù)下降速率影響很小。不同圍壓下配位數(shù)的初始值是不一樣的，圍壓越大，試樣的起始配位數(shù)越高。

2.4 孔隙度演化

圖9 天然石英砂巖配位數(shù)演化曲線Fig.9 Coordination evolution curve of natural quartz sandstone

圖10 飽和石英砂巖孔隙度演化曲線Fig.10 Porosity evolution curve of saturated quartz sandstone

圖11 天然石英砂巖孔隙度演化曲線Fig.11 Porosity evolution curve of natural quartz sandstone

孔隙度是研究試樣內(nèi)部相對位移及裂隙發(fā)育的一個重要參數(shù)。圖10和圖11分別為不同圍壓下飽和、天然石英砂巖孔隙度與加載時間、裂隙數(shù)演化關(guān)系。不同圍壓下孔隙度的初始值是不一樣的，圍壓越大，試樣的起始孔隙度越小，試樣更密實。隨著加載的進行，孔隙被壓縮，孔隙度減小，隨著進一步加載，試樣內(nèi)部產(chǎn)生裂隙，孔隙度出現(xiàn)不同程度的升高。對于飽和試樣，當圍壓為6 MPa、9 MPa時，孔隙度曲線在上升時表現(xiàn)迅速，試樣趨向于脆性破壞；而圍壓為3 MPa、12 MPa時，孔隙度曲線表現(xiàn)平穩(wěn)，試樣更趨向于塑性破壞。但天然試樣在不同圍壓下孔隙度曲線均上升迅速，試樣趨向于脆性破壞。說明浸水條件對試樣孔隙度隨時間的演化影響較大，含水情況可以改變試樣的破壞模式。

3 結(jié) 論

1) 圍壓和浸水條件只能改變數(shù)值模型統(tǒng)計角度范圍內(nèi)顆粒接觸力的大小，不能改變接觸力的分布形態(tài)。無論是法向接觸力還是切向接觸力，在試樣發(fā)生破壞時，試樣內(nèi)各個統(tǒng)計方向都存在一定的接觸力。

2) 圍壓對試樣破壞時上下兩端的法向接觸力影響較小，對左右兩側(cè)的法向接觸力影響較大，不同圍壓下試樣達到破壞時上下兩端的法向接觸力大小與抗壓成強度正比，而左右兩側(cè)的法向接觸力呈現(xiàn)出無規(guī)律性。圍壓對試樣破壞時切向接觸力影響較小。不同圍壓下試樣達到破壞時的切向接觸力大小與抗壓強度成正比。

3) 圍壓對起始配位數(shù)影響較大，圍壓越高，起始配位數(shù)越大。對配位數(shù)進入下降階段的時間點影響很大，對配位數(shù)下降速率影響很小。

4) 圍壓越大，試樣起始孔隙度越小，隨著加載的進行，試樣孔隙度下降，下降速率基本一致，含水情況可以改變試樣的破壞模式。