羅　群　劉維森

(1.長(zhǎng)沙市中等城鄉(xiāng)建設(shè)職業(yè)技術(shù)學(xué)校,湖南 長(zhǎng)沙　410000；　2.中國(guó)人壽財(cái)產(chǎn)保險(xiǎn)長(zhǎng)沙市中心支公司，湖南 長(zhǎng)沙　410000)

1　概述

水對(duì)邊坡穩(wěn)定性具有重要的影響。一方面，孔隙水壓力使土體的有效應(yīng)力降低，影響了土體的抗剪強(qiáng)度，從而降低了邊坡的穩(wěn)定性。另一方面，孔隙水的重力產(chǎn)生使邊坡土體下滑的作用效果，這部分的力對(duì)邊坡穩(wěn)定性的影響更為明顯。通常，孔隙水壓力采用滲流場(chǎng)進(jìn)行模擬。蔣中明等[1]采用FLAC3D軟件分析了邊坡在非飽和條件下的穩(wěn)定性計(jì)算方法。對(duì)于礦山邊坡，于振等[2]采用FLAC3D軟件計(jì)算了應(yīng)力和滲流耦合作用下的穩(wěn)定性。王偉興等[3]結(jié)合Geo-Studio和FLAC2D分析了邊坡的流固耦合效應(yīng)。基于離散單元法，關(guān)千軍[4]采用UDEC軟件研究了滲流對(duì)巖質(zhì)裂隙邊坡的穩(wěn)定性影響。

邊坡穩(wěn)定性計(jì)算的商業(yè)軟件均嵌入了分析滲流場(chǎng)作用的模塊。例如，F(xiàn)LAC3D采用configure fluid命令進(jìn)入滲流求解模式[5]；Geo-Studio中的Seep/W模塊模擬滲流場(chǎng)。在實(shí)際的使用中，不同邊界條件的處理對(duì)滲流的計(jì)算效果有較為明顯的差別。本文基于FLAC3D有限差分軟件，計(jì)算分析了不同滲流邊界條件對(duì)計(jì)算結(jié)果的影響，以期對(duì)邊坡的設(shè)計(jì)和施工提供一些參考。

2　FLAC3D滲流耦合的計(jì)算方法

流固耦合作用的模擬是數(shù)值模擬方法中比較復(fù)雜的問(wèn)題。FLAC3D在滲流計(jì)算方面具有一些優(yōu)勢(shì)，可以在設(shè)定流體介質(zhì)和土體介質(zhì)滲流參數(shù)的情況下，模擬得到合理的滲流場(chǎng)。同時(shí)，借助于耦合模型，能夠?qū)⒖讐鹤兓鸬膽?yīng)力場(chǎng)影響進(jìn)行分析計(jì)算。此外，滲流邊界條件對(duì)流場(chǎng)的影響也較為顯著。

2.1　滲流分析的參數(shù)設(shè)置

在FLAC3D滲流模式下，需要設(shè)置的參數(shù)包括：滲透系數(shù)、流體密度。除此之外，對(duì)于顆?？梢詨嚎s的土體，需要設(shè)置比奧系數(shù)和比奧模量；對(duì)于不可壓縮的土體，需要設(shè)置流體的體積模量和孔隙率。如軟土地基堆載等情況，需將土體視為可壓縮土體。邊坡穩(wěn)態(tài)分析中，通常將土體簡(jiǎn)化為不可壓縮土體，主要反映出孔隙水的重力作用對(duì)邊坡穩(wěn)定性的影響。此外，穩(wěn)態(tài)分析中，采用較低的流體模量可以加快滲流計(jì)算的收斂。

2.2　流體邊界

通常采用的流體邊界有兩種：不透水邊界和透水邊界。兩者最大的區(qū)別在于不透水邊界上的孔壓值可以自由變化，而透水邊界上的孔壓值是固定的。FLAC3D默認(rèn)邊界為不透水邊界。而在實(shí)際穩(wěn)定滲流情況下，邊坡的邊界上的孔壓值經(jīng)常是是固定不變的，這就需要將邊界設(shè)置為透水邊界。在FLAC3D中，透水邊界可以采用fix或apply命令進(jìn)行設(shè)置，命令如下：

{fix pp孔壓值(可省略)range節(jié)點(diǎn)范圍}；或{apply pp孔壓值range節(jié)點(diǎn)范圍}。

fix命令需要和initial pp命令聯(lián)合使用。在apply命令中可以采用grad關(guān)鍵詞來(lái)表示孔壓的變化情況，或者采用hist來(lái)定義隨時(shí)間變化的孔壓值。需要注意的是，在坡面(斜面)為透水邊界的情況下，節(jié)點(diǎn)范圍指定在該坡面上的命令如下：

{apply pp 0 range plane normal-1 0 1 origin 0 0 0 z 0 10}。

該命令中，坡面平面采用normal和origin這兩個(gè)向量參數(shù)進(jìn)行指定，范圍限制在z向0 m～10 m范圍內(nèi)的平面。

3　安全系數(shù)計(jì)算結(jié)果對(duì)比

在分析邊坡穩(wěn)定性分析中，可以采用固定孔壓的方法進(jìn)行有效應(yīng)力的分析，或者采用數(shù)值模擬所得到的穩(wěn)定流場(chǎng)進(jìn)行計(jì)算分析。相比較滲流耦合的方法，固定孔壓方法能夠節(jié)省大量的計(jì)算時(shí)間。本文分別采用固定孔壓方式和穩(wěn)定流場(chǎng)兩種模式進(jìn)行分析邊坡穩(wěn)定性。

建立如圖1所示的模型。邊坡高度為10 m，坡角為45°。h為右側(cè)邊界水位距坡趾的距離。分別選取三種計(jì)算條件，即h=-0.25 m，0 m，5 m。在h=5 m的情況下，即坡體外水壓對(duì)坡面具有一定的作用效果。這部分力對(duì)邊坡產(chǎn)生反壓作用。本文采用applynstress命令，將水壓力視為法向力作用于坡面，從而模擬這部分的作用效果。

表1　計(jì)算參數(shù)取值表

固定孔壓情況下的計(jì)算過(guò)程為：

1)生成初始應(yīng)力場(chǎng)，位移場(chǎng)清零。2)采用water table命令生成水面，水面以下的孔壓由程序自動(dòng)計(jì)算完成。該計(jì)算條件下，不需要打開(kāi)config fluid模式。3)采用強(qiáng)度折減法和二分算法，計(jì)算邊坡的安全系數(shù)，程序自動(dòng)計(jì)算得到邊坡在臨界狀態(tài)下的滑移面。穩(wěn)定滲流場(chǎng)情況下的計(jì)算過(guò)程為：1)生成初始應(yīng)力場(chǎng)，位移場(chǎng)清零。2)打開(kāi)config fluid滲流模式。采用initial pp命令指定邊界孔壓分布情況，通過(guò)fix或apply命令將邊界設(shè)置為透水條件，即固定邊界處的孔壓。3)關(guān)閉應(yīng)力計(jì)算模式，計(jì)算穩(wěn)定滲流場(chǎng)。命令為：set fluid on mech off。4)關(guān)閉滲流場(chǎng)計(jì)算，打開(kāi)應(yīng)力場(chǎng)模式，同時(shí)將流體模量設(shè)置為0，避免應(yīng)力變化再次引起孔壓的變化。

兩種模式下，均采用強(qiáng)度折減法和二分算法，計(jì)算得到邊坡的安全系數(shù)。在不同水頭情況下，邊坡的安全系數(shù)計(jì)算結(jié)果如表2所示。從表2可以看出，采用滲流耦合方法計(jì)算得到的安全系數(shù)略小于固定孔壓模式。其原因是，滲流模式的孔壓值略大于固定孔壓模式。此外，隨著邊坡外水面的升高，邊坡的安全系數(shù)在增大。特別是h=5 m的情況，邊坡外水的反壓作用比較明顯，安全系數(shù)有明顯的提高。圖2，圖3為流場(chǎng)的等值線云圖。圖4，圖5分別是h=0和h=5 m情況下的剪切應(yīng)變?cè)隽吭茍D，圖中存在明顯的剪切帶，即為邊坡的滑移面。從圖4和圖5中可以看出，隨著坡體外水位的增高，滑移面具有明顯的變化。h=5 m時(shí)，邊坡的滑移范圍擴(kuò)大，滑移面有可能越過(guò)坡趾，形成垮塌。

表2　安全系數(shù)計(jì)算結(jié)果對(duì)比

4　結(jié)語(yǔ)

基于FLAC3D分析軟件，模擬了邊坡的滲流場(chǎng)。對(duì)于滲流場(chǎng)邊界的設(shè)置方法以及需要注意的問(wèn)題進(jìn)行了說(shuō)明。分別計(jì)算并討論了不同流場(chǎng)邊界條件對(duì)滲流場(chǎng)的影響。結(jié)合強(qiáng)度折減法，計(jì)算了固定孔壓模式和穩(wěn)定滲流模式兩種情況下邊坡的安全系數(shù)。計(jì)算結(jié)果表明滲流耦合方法計(jì)算得到的安全系數(shù)略小于固定孔壓模式。此外，坡體外水壓對(duì)坡面的反壓效果比較明顯，安全系數(shù)增大，且潛在的滑移面變化較大。

參考文獻(xiàn):