貴州某邊坡穩定性分析及數值模擬研究

陳信宏

（貴州大學資源與環境工程學院，貴州 貴陽 550025）

貴州由于獨特的自然地理條件和環境，具有獨特的地形地貌特征及降雨特點，雨水入滲使得斜坡土體的飽和度增加，降低土體抗剪強度，從而坡體發生滑動；另一方面，降雨使地下水位上升，含水量增加、水壓力增大，從而降低坡體穩定性。邊坡穩定性分析方法有定性和定量分析，Dong Tang 等（2020）分別考慮了滑移面處于滯水帶中和滯水帶外，通過列舉算例進行驗證，表明所提出的計算方法在圓弧滑面和復合滑面邊坡安全系數計算的合理性。隨著科學技術發展，借助于計算軟件的數值模擬方式不斷推陳出新，相比于物理模擬試驗，選用適宜軟件進行數值模擬具有速率高、對時間和人力物力的要求稍低，也從一定程度上解決了物理試驗缺乏的條件等特點。劉海軍（2012）通過UDEC、FLAC3D 對皖南山區特大型滑坡形成過程進行驗證，展示滑坡經歷彎曲變形、彎曲-傾倒、傾倒-折斷和滑移-剪斷失穩破壞的動態演化過程，發現傾倒變形體坡頂后緣主要以拉張破壞為主，坡腳沿折斷面往坡體內部主要以剪切破壞為主。王浩宇、于永強（2020）以三峽古滑坡為例，利用GEO-slope 和FLAC-3D 軟件研究滑坡可能產生深層和淺層滑動的多種工況，并進一步分析斜坡穩定性。高相波等（2019 年）采用FLAC3D 對浙江神仙居景區神龍瀑邊坡進行定性分析和定量計算，發現神龍瀑邊坡在受拉狀態下危巖體存在整體位移的特征，提出景區危巖體是神龍瀑邊坡潛在破壞區的預測。

本文依托貴陽市某基坑邊坡，就邊坡開挖卸荷狀態，先由理論公式定量分析邊坡穩定性，然后通過數值模擬對卸荷條件下邊坡巖體強度參數的變化進行研究，結果發現坡體第一階梯坡腳處出現較大位移，滑動面主要位于坡體內部。

1 邊坡地質概況

1.1 地形地貌及地質構造

研究區邊坡處貴陽市，屬亞熱帶濕潤溫和型氣候，高原丘陵地形。區內構造較為簡單，無斷層。巖石地層產狀變化不大，下伏基巖為三疊系大冶組（T1d）灰白色薄白云巖，巖層產狀122°∠20°。

1.2 地層巖性

區內巖土構成自上而下依次為：研究區上覆土層為第四系素填土（Q4ml）、第四系殘積紅粘土層（Q4el+dl），下伏基巖為三疊系大冶組（T1d）石灰巖。

2 邊坡穩定性計算

2.1 計算剖面選取

邊坡開挖后基底標高1275.3～1279，形成高6.2～11m 的多段基坑邊坡，周邊條件較復雜，地下水位高，選取其中較具代表性的F—G 段進行計算，工程剖面圖如圖1 所示，邊坡總長約69m，坡高6.4～11.1m。上部覆蓋紅粘土層厚0.5～1m，下部為石灰巖，為巖質順向坡。坡頂受林地保護限制，無放坡空間。邊坡下部屬于臨時坡（1179～1184.85m 部分），上部屬于永久邊坡（1184.85m 以上），基坑邊坡周圍有建筑物存在，破壞后果嚴重，據《建筑邊坡工程技術規范》（GB 50330—2013），邊坡安全等級為一級，重要性系數γ0=1.1。

圖1 工程剖面圖

圖2 網格劃分及按土層分組圖示

2.2 穩定性系數計算

采用極限平衡法計算邊坡抗滑穩定性，所用公式為《建筑邊坡工程技術規程》（GB 50330—2013）附錄A公式（公式略）。依據相關規范以及工程地質手冊查閱取得，所取參數如下表1 所示：

表1 邊坡巖土層參數建議值表

計算得F—G 總下滑力3773.410 KN，總抗滑力3962.136 KN，穩定系數Fs=1.05。由于基坑邊坡開挖影響范圍較大，造成穩定性下降，存在安全隱患，故分析其穩定性尤為重要。

3 采用FLCA 3D 數值模擬分析

3.1 模型的建立及網格劃分

利用FLAC 3D 內置Extrusion 板塊，依據剖面圖建立的模型簡化并建模并劃分網格和單元劃分，巖土體參數見表1。

對模型底面、兩側及前后側進行位移約束和應力邊界條件設置，重力加速度為10，初始應力狀態以默認參數計算，較短時間內就收斂，說明自然狀態下該邊坡并未滑動，仍然處于穩定狀態，經過初始地應力的唯一和速度清零后，只保留應力，本構模型設置為摩爾庫倫模型進行計算。

3.2 計算結果及分析

使用FLAC 3D 進行模擬，研究邊坡發生滑動時變形特征的過程中，采用強度折減法將模型材料的強度參數進行連續折減，直到邊坡失穩發生滑動，邊坡失穩時邊坡整體位移變化情況如圖3 所示。

圖3 總位移云圖

由圖 3 總位移云圖可見，邊坡總位移主要體現在坡中上部，整個滑體在自重作用下向下發生蠕動變形，在水和上部荷載作用下，坡體蠕變加劇，若在坡腳開挖，導致坡腳產生新的應力集中區，發生剪出導致邊坡破壞。填土層與粘土層接觸的第一階梯坡腳位置發生變形且位移較大；滑坡應為牽引式破壞，滑動發生于坡體之內。此外，粘土層前緣在滑坡變形影響范圍內，滑坡發生滑動時，前端可能發生變形。

坡體水平方向的位移如圖4，趨勢與總位移一致，填土層與粘土層交接處水平位移最大，水平向位移的變化規律和分布規律為：在橫向上的位移主要在坡中位置，其次在坡腳位置，破壞由后緣先發生，破壞主要方向朝臨空面方向。大小均具有向坡體前緣呈遞增趨勢。

圖4 水平方向（X 方向）上的位移云圖

圖5 為豎直方向位移云圖，從該位移云圖可以看出在豎直方向上發生的位移主要集中于坡坡體的填土層，呈現逐步向下方的粘土層及坡腳擴散的現象，斜坡往后再發展可能朝著巖土層接觸面附近發生圓弧滑動，并最后會由坡腳剪出。

圖5 豎直方向（Z 方向）上的位移云圖

圖6 和圖7 為最大剪應變率和豎直方向上應變速率云圖。結合前文總位移云圖不難看出，當坡體產生較大位移時，滑動帶位置主要呈現為剪應變區，即為所推測的斜坡發生圓弧滑動的區域。隨著滑坡體淺層滑帶發展延伸至坡體后緣，為坡體主要塑性區，當潛在滑動帶發生貫通后，斜坡失穩。

圖6 最大剪應變率云圖

圖7 豎直方向（Z 方向）應變率云圖

綜上所述，在經開挖和削方減載后，該斜坡可能發生失穩破壞，且存在圓弧滑動面，在變形不斷發展下應力不斷變化，坡體發生滑動時總位移與水平位移均為填土層與粘土層交接部位前緣較大，向后緣和前端坡腳部位逐漸遞減，而豎向位移則主要集中于坡體上部陡緩交接部位及中前部粘土層坡體部位，滑動帶逐步發展至貫通并最后造成坡體失穩。

4 結論

結合相關參考資料及規范標準，通過定性分析與定量分析相結合，主要用公式計算和數值模擬，進行貴陽市某基坑邊坡穩定分析，采用FLAC 3D 進行不發生降雨情況下邊坡變形過程模擬研究，可得到以下認識：

（1）通過公式計算邊坡安全系數，F-G 段邊坡典型剖面安全系數為1.05，邊坡處于不穩定狀態。

（2）針對不穩定斜坡，采用 FLAC3D 模擬并分析邊坡實際的位移特征：邊坡水平方向上的位移主要在坡中位置，其次位移明顯處在坡腳位置，破壞由后緣先發生，產生破壞，相應位移變化在圖中顯示，破壞主要方向朝臨空面方向。

（3）潛在滑動面近似于圓弧形，滑面經兩個圖層中部位置及巖土分界面，剪出口滑體與基巖接觸面。土層分界面處發生了主要的變形特征，且位移較大，該滑坡為牽引式破壞，滑動面位于坡體之內。