邊坡上部基坑開挖數值模擬分析

魏 達 劉思奇 呂美群 申楊凡

(東北林業大學土木工程學院，黑龍江 哈爾濱 150040)

邊坡上部基坑開挖數值模擬分析

魏 達 劉思奇 呂美群 申楊凡

(東北林業大學土木工程學院，黑龍江 哈爾濱 150040)

為了解邊坡上部基坑開挖對于邊坡穩定性的影響，利用有限元軟件ansys對邊坡上部不同深度的基坑開挖進行了數值模擬，通過計算得到了邊坡在開挖前后水平方向的變位以及邊坡塑性應變的差異，所得結果可為邊坡支擋防護等相關工程的設計提供參考。

邊坡，穩定性，基坑開挖，有限元計算軟件

0 引言

在巖土與水利建設中，邊坡工程是十分重要的，邊坡的穩定性會嚴重影響工程的施工和運營安全。邊坡上部進行一定深度的基坑開挖既破壞了原有邊坡的穩定性能，同時也是一種上部結構的卸載過程，因此邊坡上部的基坑開挖對于邊坡的穩定性有一定的影響。

本文通過有限元模型對于邊坡進行不同深度的基坑開挖數值模擬，可以得到不同開挖深度的基坑對于邊坡x方向的變位以及塑性開展區。

1 計算模型

本文采用大型通用有限元計算軟件ansys建立邊坡平面模型，由于邊坡可以認為是縱向較長的實體，因此將幾何模型簡化為平面問題進行簡化是合理的。模型單元采用材料庫中常見的Plane82單元，此平面單元由于具有8個節點，具有較大的精準度，可以很好的適應不規則網格。定義邊坡截面中邊坡的高度為90 m，邊坡坡率為1.5。

本次邊坡模型選用D-P(Drucker-Prager)模型作為土體的本構模型，ansys中的D-P模型是使用Von Mises屈服作為相應判斷條件。在D-P模型中，其等效應力σe以及Von Mises屈服法則的計算公式如式(1)～式(3)所示[1]。

σe=3βσm+[{s}T[M]{s}/2]1/2

(1)

F=3βσm+[{s}T[M]{s}/2]1/2-σy=0

(2)

(3)

式中：σm——靜水壓力；

β——材料參數；

{s}——偏差應力矩陣；

C——粘聚力；

φ——土體的剪脹角。

本例中，土體采用的粘聚力為123.13 kPa，內摩擦角為14.96°，剪脹角設置為0。

建立相關幾何模型以后對模型進行網格劃分并施加約束，對有限元模型的下部土體進行完全約束，并對邊坡左右兩邊進行橫向膨脹約束，對模型進行自重荷載加載，通過殺死一定單元的方式模擬邊坡上部的基坑開挖進行計算。其有限元模型示意如圖1所示。

2 計算結果分析

對邊坡上部進行不同深度的基坑開挖，其中所開挖的基坑距離斜坡為20 m，基坑寬度為20 m，所開挖的邊坡基坑深度分別為5 m，10 m，15 m，20 m。經過計算后進入通用后處理可以查看得到邊坡的x方向以及y方向的位移變形云圖、塑性應變云圖等。如果邊坡的x變形較大，塑性區持續開展形成貫通區域時，就容易發生邊坡塌陷等災害，因此通過邊坡的x變形以及邊坡塑性變形可以綜合判斷邊坡穩定性。

通過邊坡上部不同深度的基坑開挖進行求解得到邊坡x方向的位移云圖以及塑性云圖，與未開挖的邊坡模型計算進行對比，其中x方向隨著基坑開挖深度的不斷增大邊坡上部位移最大值變化如圖2所示，未開挖時的邊坡x方向位移云圖與上部基坑開挖20 m時的x方向位移云圖如圖3所示。

通過圖2以及圖3可以看出，基坑開挖前后，邊坡在-x方向的位移值要遠大于在x方向的位移值，且-x方向位移最大值以及x方向位移最大值的出現位置均相同，-x方向位移最大值均出現在坡腳處，而x方向位移最大值均出現在坡頂處。隨著邊坡開挖深度的不斷增大，-x方向的位移值呈現四次拋物線的形式下滑，這是由于基坑開挖對于邊坡而言是一個卸載的過程，此時對于邊坡穩定性是有利的。隨著邊坡開挖深度的不斷增大，x方向的位移值呈現四次拋物線的形式上升，這是由于出現x方向的位移最大值在坡頂處，而在坡頂處進行基坑開挖則會破壞局部結構的穩定性，使得基坑附近的土體更容易塌陷。

通過有限元計算所得到的塑性區開展的塑性應變最大值如圖4所示，基坑開挖前后的邊坡塑性應變云圖如圖5所示。

通過圖4以及圖5可以看出，無論基坑開挖前后，邊坡的塑性區開展都是從坡腳處開始，其塑性應變最大值都在坡腳處。隨著開挖深度的不斷增大，其塑性應變最大值出現了先降再升再降的雙波谷式變化，這是由于最開始的基坑開挖對于邊坡土體是一個卸載作用，使得坡腳的塑性應變減小，從圖5b)中可以看出，當基坑深度達到20 m時，此時的基坑附近也出現了較大的塑性區，如基坑深度不斷增大，其塑性區開展容易形成貫穿區域，不利于邊坡的穩定性。

3 結語

通過有限元計算軟件ansys對不同深度的邊坡上部基坑進行開挖的數值模擬，可以得到以下幾點結論：

1)不同深度的基坑開挖對邊坡的穩定性有一定的影響。

基坑開挖前后，邊坡x方向最大值出現的位置是一定的，均出現在坡腳處。但是隨著邊坡開挖深度的不斷增大，-x方向的位移值呈現四次拋物線的形式下滑，而x方向的位移值呈現四次拋物線的形式上升，因此在邊坡上部進行開挖時，盡量采用相應的支擋結構。

2)隨著邊坡開挖深度的不斷增大，邊坡的塑性應變最大值出現了先降再升再降的雙波谷式變化，當基坑深度達到20 m時，此時的基坑附近也出現了較大的塑性區，因此較大的基坑開挖不利于邊坡的穩定性。

[1] 李 婕.基于ANSYS加筋路堤沉降的有限元分析[J].山西建筑,2008,34(26):262-263.

[2] 劉世強，丁佳攀，何鈺龍，等.基坑開挖對邊坡穩定性的影響分析[J].建材發展導向,2014(5):28.

[3] 王恒星,孫 碩,劉 潔,等.凍融作用對路堤邊坡穩定性的影響[J].北方交通,2014(12):47-49，54.

[4] 閻 波，方 云.基于ANSYS的邊坡開挖模擬[J].山西建筑,2008,34(1):93-94.

[5] 包 磊.基于ansys的邊坡穩定性強度折減有限元分析[J].土工基礎，2012，26(8)：48-51.

Numerical simulation analysis of excavation in the upper part of slope

Wei Da Liu Siqi Lv Meiqun Shen Yangfan

(SchoolofCivilEngi.,NortheastForeastryUniv.,Harbin150040,China)

In order to solution on the upper part of slope excavation influence on slope stability, using finite element software ansys, the upper slope of different depth of foundation pit excavation numerical simulation and calculation of the slope in the horizontal direction before and after the excavation of the displacement and slope plastic strain differences. The results can provide a reference for the design of slope retaining and protection.

slope， stability， foundation pit excavation， finite element software

2015-05-23

魏 達(1992- )，男，在讀本科生； 劉思奇(1993- )，男，在讀本科生； 呂美群(1993- )，女，在讀本科生； 申楊凡(1993- )，男，在讀本科生

1009-6825(2015)22-0098-02

TU463

A