馬 剛 ，常曉林 ，劉嘉英 ，周 偉

（1. 武漢大學(xué) 水資源與水電工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430072；2. 武漢大學(xué) 水工巖石力學(xué)教育部重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430072）

1 引 言

非關(guān)聯(lián)的彈塑性巖土材料尤其是顆粒類材料在塑性極限內(nèi)存在多種失穩(wěn)模式，典型的如應(yīng)變局部化失穩(wěn)。當(dāng)應(yīng)力狀態(tài)處于塑性極限面以內(nèi)時，一定情況下也會發(fā)生分散性失穩(wěn)。術(shù)語“分散性失穩(wěn)”對應(yīng)于應(yīng)變局部化失穩(wěn)，強(qiáng)調(diào)了失穩(wěn)模式中沒有出現(xiàn)應(yīng)變局部化或者剪切帶。

相比于對應(yīng)變局部化的研究，對分散性失穩(wěn)模式的研究比較少。Castro等[1]在對飽和松砂進(jìn)行靜力加載試驗(yàn)時發(fā)現(xiàn)，完全不排水條件下偏應(yīng)力q達(dá)到峰值后，在試樣軸向施加微小的應(yīng)力增量，試樣整體馬上發(fā)生急劇的失穩(wěn)，應(yīng)力狀態(tài)尚處于峰值強(qiáng)度破壞面以內(nèi)，稱之為靜力液化，以區(qū)別于飽和砂土在循環(huán)加載中出現(xiàn)的動力液化。隨著試樣相對密度的增加，靜力液化現(xiàn)象越來越弱并最終消失。研究人員通過室內(nèi)試驗(yàn)，研究了排水條件、圍壓、應(yīng)力水平、相對密度等因素對靜力液化起動機(jī)制的影響[2－7]。Lade等[3]將不同試驗(yàn)條件下靜力液化時的應(yīng)力狀態(tài)在應(yīng)力空間內(nèi)連起來形成失穩(wěn)線，當(dāng)試樣的應(yīng)力水平在失穩(wěn)狀態(tài)線以上時，在不排水條件下通過軸向加載會發(fā)生急劇失穩(wěn)。

通過對大量邊坡失穩(wěn)案例的分析，Olson[8]認(rèn)為由于地下水引起的靜力液化可能是邊坡失穩(wěn)的隱含機(jī)制之一。認(rèn)識到靜力液化與滑動失穩(wěn)的密切聯(lián)系后，大大擴(kuò)展了人們對巖土失穩(wěn)機(jī)制的認(rèn)知深度和研究范圍[9－10]。

隨著研究的深入人們發(fā)現(xiàn)，對密實(shí)的顆粒集合體，通過控制試樣在剪切過程中的體積變形也會出現(xiàn)與靜力液化相似的現(xiàn)象，即分散性失穩(wěn)。由于不排水剪切試驗(yàn)時試樣的體積保持不變，可以認(rèn)為松砂在不排水剪切試驗(yàn)中出現(xiàn)的靜力液化是分散性失穩(wěn)的一個特例。通常在室內(nèi)試驗(yàn)或離散元模擬中沿不同的應(yīng)變路徑剪切顆粒集合體，控制其體積變形，研究其在不同體變情況下可能發(fā)生的分散性失穩(wěn)。Lancelot等[11]、Jrad等[12]采用常規(guī)的試驗(yàn)裝置，研究了顆粒材料在應(yīng)變加載路徑下的力學(xué)特性。Darve等[9]和Nicot等圍繞顆粒材料的分散性失穩(wěn)，綜合采用室內(nèi)試驗(yàn)[13－14]、理論建模[15－16]和離散元模擬[17－18]開展了大量富有建設(shè)性的研究工作。

本文采用連續(xù)-離散耦合分析方法（combined finite-discrete element method，F(xiàn)DEM）模擬密實(shí)的顆粒集合體在不同應(yīng)變加載路徑下的力學(xué)特性，研究其在特定應(yīng)變路徑下可能出現(xiàn)的分散性失穩(wěn)。以往的研究表明，F(xiàn)DEM能體現(xiàn)顆粒層面的摩擦、咬合、破碎和重排列對宏觀力學(xué)特性的影響，再現(xiàn)巖土顆粒材料，如粗粒土、堆石的非線性、壓硬性、剪脹和剪縮特性[19－20]。

2 等比例應(yīng)變加載路徑的數(shù)值試驗(yàn)

2.1 數(shù)值試樣

采用隨機(jī)顆粒生成程序生成形狀隨機(jī)的不規(guī)則多面體顆粒，按照體積等效原則，定義顆粒的等效直徑為

生成的數(shù)值試樣如圖 1(a)所示。試樣形狀為300 mm×300 mm×300 mm的立方體，共包含8 927個不規(guī)則凸多面體顆粒，試樣中顆粒的等效粒徑分布如圖 1(b)所示，采用 2階四面體單元離散為142 277個單元，401 590個節(jié)點(diǎn)。

圖1 數(shù)值試樣及試樣中顆粒的等效粒徑分布Fig.1 Numerical sample and equivalent particle size distribution

2.2 細(xì)觀參數(shù)

數(shù)值試驗(yàn)需要的細(xì)觀參數(shù)較多，其中部分參數(shù)可以通過常規(guī)物理力學(xué)試驗(yàn)直接確定。本研究不針對特定的巖土顆粒材料，只是在一般參數(shù)取值范圍內(nèi)選取了一套參數(shù)進(jìn)行數(shù)值模擬，參數(shù)見表1。

2.3 等比例應(yīng)變加載路徑

為了研究堆石料的分散性失穩(wěn)特性，本文設(shè)計(jì)了一組軸對稱的等比例應(yīng)變加載路徑數(shù)值試驗(yàn)，數(shù)值試樣和計(jì)算參數(shù)分別如圖1、表1所示。采用應(yīng)變增量的形式，將軸對稱的等比例應(yīng)變加載路徑：

表1 數(shù)值試驗(yàn)所用的細(xì)觀參數(shù)Table 1 Set of parameters used in simulations

當(dāng)R = 1時，試樣體積變形為0，對應(yīng)于不排水剪切試驗(yàn)；當(dāng)R＞1時，試樣發(fā)生體積收縮，R = 1時對應(yīng)于單軸壓縮試驗(yàn)；當(dāng)0＜R＜1時，試樣發(fā)生體積膨脹。

在軸對稱的等比例應(yīng)變加載路徑數(shù)值試驗(yàn)中，先對試樣施加各向等壓的應(yīng)力狀態(tài)直至預(yù)定的圍壓值，然后按照上述的等比例應(yīng)變加載路徑分別控制試樣軸向和側(cè)向加載板的位移，以保證試樣的應(yīng)變增量滿足式（1）。分別進(jìn)行圍壓為0.8、1.6、2.4、3.2 MPa下的等比例應(yīng)變加載路徑數(shù)值試驗(yàn)，每個圍壓下控制應(yīng)變路徑的參數(shù) R取 1.333、1.000、0.667、0.500和0.400。

3 分散性失穩(wěn)模式

圖2為不同等比例應(yīng)變路徑的剪應(yīng)力p與平均應(yīng)力q關(guān)系曲線。當(dāng)R≥0.667時，p-q曲線沒有出現(xiàn)峰值，剪應(yīng)力持續(xù)增大至試驗(yàn)結(jié)束；當(dāng) R = 0.5時，p-q曲線存在一個局部峰值，剪應(yīng)力在達(dá)到這個局部峰值后經(jīng)過短暫的降低，此后又繼續(xù)增大；當(dāng)R = 0.4時，剪應(yīng)力q在達(dá)到峰值后迅速減小，表明試樣已喪失了承載能力，發(fā)生了分散性失穩(wěn)。

圖2 子午面上的應(yīng)力路徑Fig.2 Stress paths in the meridian plane

圖3 不同等比例應(yīng)變路徑的剪應(yīng)力、應(yīng)力差及2階功演化曲線Fig.3 Evolutions of deviator stress, modified deviator stress and second-order work under different proportional strain path loading tests

圖4 典型的分散性失穩(wěn)模式Fig.4 Unstable state related to a diffuse mode of failure strictly inside the Mohr–Coulomb plastic limit condition.

根據(jù)本文的數(shù)值試驗(yàn)結(jié)果，將p-q平面內(nèi)的應(yīng)力狀態(tài)分成如圖5所示的3個區(qū)域。以R = 1的等體積加載路徑得到的p-q曲線為界，分為剪脹區(qū)和剪縮區(qū)，剪縮區(qū)內(nèi)的應(yīng)力狀態(tài)是穩(wěn)定的，而剪脹區(qū)內(nèi)以R = 0.4的p-q曲線為界，分為剪脹-穩(wěn)定區(qū)和剪脹-非穩(wěn)定區(qū)。剪脹-穩(wěn)定區(qū)內(nèi)試樣的宏觀力學(xué)特性表現(xiàn)為剪應(yīng)力應(yīng)變硬化、體積剪脹，而剪脹-非穩(wěn)定區(qū)內(nèi)試樣的宏觀力學(xué)特性則為剪應(yīng)力應(yīng)變軟化、體積剪脹。

4 失穩(wěn)線

將不同圍壓下發(fā)生分散性失穩(wěn)時的應(yīng)力狀態(tài)在p-q平面內(nèi)連接起來形成失穩(wěn)線，同時對圖 1(a)的顆粒集合體，采用表1參數(shù)，進(jìn)行等p應(yīng)力路徑的常規(guī)三軸剪切試驗(yàn)，得到不同平均應(yīng)力p情況下的峰值和臨界應(yīng)力狀態(tài)。圖6為p-q平面上的峰值強(qiáng)度線、臨界狀態(tài)線和失穩(wěn)線。在相同平均應(yīng)力p下有峰值強(qiáng)度線>臨界狀態(tài)線>失穩(wěn)線，圖中失穩(wěn)線以上的深色部分即為在峰值強(qiáng)度線以內(nèi)的不穩(wěn)定應(yīng)力狀態(tài)區(qū)域。

圖5 基于2階功判據(jù)的應(yīng)力狀態(tài)分區(qū)Fig.5 Zones of stability and instability under forced compaction and dilation

圖6 子午面上的峰值強(qiáng)度線、臨界狀態(tài)線和失穩(wěn)線Fig.6 Failure line, critical state line and instability line in meridian plane

沿等p應(yīng)力路徑將試樣剪切至不同的應(yīng)力水平P1～P9，再沿R = 0.400的等比例應(yīng)變路徑加載，結(jié)果如圖7所示。圖中，P1～P4位于失穩(wěn)線以下，P5位于失穩(wěn)線上，P6～P8位于不穩(wěn)定應(yīng)力狀態(tài)區(qū)域，P9位于峰值強(qiáng)度線上。從圖中可以看出，當(dāng)應(yīng)力水平位于失穩(wěn)線以下時，在等比例應(yīng)變加載過程中，試樣的剪應(yīng)力逐漸增大至峰值，峰值點(diǎn)的應(yīng)力狀態(tài)基本都落在失穩(wěn)線附近。當(dāng)應(yīng)力水平位于失穩(wěn)線及其以上（如P5～P9）時，經(jīng)過等比例應(yīng)變路徑加載，試樣的剪應(yīng)力徑直減小，應(yīng)力狀態(tài)貼著破壞線向下發(fā)展。

將試樣剪切至不同的應(yīng)力水平，經(jīng)過相同的等比例應(yīng)變加載路徑后，試樣最終都會發(fā)生分散性失穩(wěn)。不同應(yīng)力水平的區(qū)別是：當(dāng)應(yīng)力水平位于失穩(wěn)線以下時，在剛開始加載時試樣尚具有一定的抗剪切能力；應(yīng)力水平位于失穩(wěn)線及其以上時，試樣在等比例應(yīng)力加載路徑下逐漸喪失承載能力。

圖7 剪切至不同應(yīng)力水平后的等比例應(yīng)變加載路徑Fig.7 Simulated response of numerical sample after shearing to different stress levels

5 結(jié) 論

（1）給定的應(yīng)變路徑下平均應(yīng)力-剪應(yīng)力關(guān)系曲線存在一個峰值點(diǎn)，達(dá)到這個峰值點(diǎn)后，如果繼續(xù)施加一個無限小的軸向荷載，試樣整體出現(xiàn)急劇失穩(wěn)。采用應(yīng)變控制加載時，試驗(yàn)可以繼續(xù)進(jìn)行下去，此時平均主應(yīng)力和剪應(yīng)力持續(xù)減小直至0。

（2）根據(jù)Hill的材料失穩(wěn)準(zhǔn)則，對一個體積為V材料系統(tǒng)，其總的2階功為負(fù)時，即使不改變施加在材料邊界上的外荷載，也沒有外部能量的輸入，材料系統(tǒng)也會發(fā)生這種不可逆的整體失穩(wěn)。

（3）將子午面的應(yīng)力狀態(tài)分成3個區(qū)域，以等體積加載路徑得到的平均應(yīng)力-剪應(yīng)力曲線為界，分為剪脹區(qū)和剪縮區(qū)。剪縮區(qū)內(nèi)的應(yīng)力狀態(tài)是穩(wěn)定的，而剪脹區(qū)分為剪脹-穩(wěn)定區(qū)和剪脹-非穩(wěn)定區(qū)。

（4）在相同平均應(yīng)力下有峰值強(qiáng)度線高于臨界狀態(tài)線，臨界狀態(tài)線高于失穩(wěn)線。峰值強(qiáng)度線與失穩(wěn)線之間的區(qū)域即為不穩(wěn)定應(yīng)力狀態(tài)區(qū)域。將試樣剪切至不同的應(yīng)力水平，經(jīng)過相同的等比例應(yīng)變加載路徑后試樣最終都會發(fā)生分散性失穩(wěn)。

[1] CASTRO G, POULOS S J. Factors affecting liquefaction and cyclic mobility[J]. Journal of the Geotechnical Engineering Division, 1977, 103(6): 501－506.

[2] KRAMER S L, SEED H B. Initiation of soil liquefaction under static loading conditions[J]. Journal of Geotechnical Engineering, 1988, 114(4): 412－430.

[3] LADE P V. Static instability and liquefaction of loose fine sandy slopes[J]. Journal of Geotechnical Engineering,1992, 118(1): 51－71.

[4] YAMAMURO J A, LADE P V. Static liquefaction of very loose sands[J]. Canadian Geotechnical Journal, 1997,34(6): 905－917.

[5] UTHAYAKUMAR M, VAID Y P. Static liquefaction of sands under multiaxial loading[J]. Canadian Geotechnical Journal, 1998, 35(2): 273－283.

[6] LEONG W K, CHU J, TEH C I. Liquefaction and instability of a granular fill material[J]. Geotechnical Testing Journal, 2000, 23(2): 178－192.

[7] WANATOWSKI D, CHU J. Static liquefaction of sand in plane strain[J]. Canadian Geotechnical Journal, 2007,44(3): 299－313.

[8] OLSON S M. Liquefaction analysis of level and sloping ground using field case histories and penetration resistance[D]. U.S.A.: University of Illinois, 2001.

[9] DARVE F, LAOUAFA F. Instabilities in granular materials and application to landslides[J]. Mechanics of Cohesive-Frictional Materials, 2000, 5(8): 627－652.

[10] BORJA R I. Condition for liquefaction instability in fluid-saturated granular soils[J]. Acta Geotechnica, 2006,1(4): 211－224.

[11] LANCELOT L, SHAHROUR I, AL M M. Instability and static liquefaction on proportional strain paths for sand at low stresses[J]. Journal of Engineering Mechanics,2004, 130(11): 1365－1372.

[12] JRAD M, SUKUMARAN B, DAOUADJI A.Experimental analyses of the behaviour of saturated granular materials during axisymmetric proportional strain paths[J]. European Journal of Environmental and Civil Engineering, 2012, 16(1): 111－120.

[13] NICOT F, DAOUADJI A, HADDA N, et al. Granular media failure along triaxial proportional strain paths[J].European Journal of Environmental and Civil Engineering, 2013, 17(9): 777－790.

[14] DAOUADJI A, Algali H, DARVE F, et al. Instability in granular materials: experimental evidence of diffuse mode of failure for loose sands[J]. Journal of Engineering Mechanics, 2009, 136(5): 575－588.

[15] NICOT F, DARVE F. A micro-mechanical investigation of bifurcation in granular materials[J]. International Journal of Solids and Structures, 2007, 44(20): 6630－6652.

[16] NICOT F, DARVE F. Diffuse and localized failure modes:Two competing mechanisms[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2011, 35(5): 586－601.

[17] NICOT F, HADDA N, BOURRIER F, et al. Failure mechanisms in granular media: A discrete element analysis[J]. Granular Matter, 2011, 13(3): 255－260.

[18] SIBILLE L, DONZé F V, NICOT F, et al. From bifurcation to failure in a granular material: A DEM analysis[J]. Acta Geotechnica, 2008, 3(1): 15－24.

[19] 馬剛, 周偉, 常曉林, 等. 堆石體三軸剪切試驗(yàn)的三維細(xì)觀數(shù)值模擬[J]. 巖土工程學(xué)報(bào), 2011, 33(5): 746－753.MA Gang, ZHOU Wei, CHANG Xiao-lin, et al. 3D mesoscopic numerical simulation of triaxial shear tests for rockfill[J]. Chinese Journal of Geotechnical Engineering,2011, 33(5): 746－753.

[20] 馬剛, 周偉, 常曉林, 等. 考慮顆粒破碎的堆石體三維隨機(jī)多面體細(xì)觀數(shù)值模擬[J]. 巖石力學(xué)與工程學(xué)報(bào),2011, 30(8): 1671－1682.MA Gang, ZHOU Wei, CHANG Xiao-lin, et al.Mesoscopic numerical simulation of rockfill considering particle breakage by using three-dimensional stochastic polyhedrons[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(8): 1671－1682.