劉占宏

(陜西國防工業(yè)職業(yè)技術學院 建筑與熱能分院,陜西 西安 710302)

?

基于有限元原理對建筑基坑邊坡支護穩(wěn)定性的計算研究

劉占宏

(陜西國防工業(yè)職業(yè)技術學院 建筑與熱能分院,陜西 西安 710302)

選取某一實際護坡支護工程進行有限元建模，并采用MATLAB8.0軟件以及ANSYS9.0軟件找到滑裂面，計算出滑裂面上切應力的數值以及最大切應力的作用位置，再通過有限元法計算出邊坡支護穩(wěn)定性安全系數，并與傳統(tǒng)方法計算得出的安全系數進行對比。結果顯示：有限元法計算得出的安全系數較大，且?guī)追N方法的安全系數變化趨勢是相同的。

有限元; 安全系數; 滑裂面; 穩(wěn)定性

隨著國內工程建設的不斷增多，許多建設部門都會涉及到大量的邊坡問題，邊坡一旦發(fā)生失穩(wěn)破壞，將會發(fā)生嚴重危害，因此，把邊坡失穩(wěn)災害降到最低一直是工程設計人員所不懈追求的目標。鑒于此，本研究擬通過對某實地的有限元建模以探尋建筑工地中人造邊坡的穩(wěn)定性問題，旨在發(fā)現其中規(guī)律，并將其推廣到更廣泛的邊坡穩(wěn)定性研究領域之中。

1 有限元建模原理

有限元方法基本求解步驟：根據實際問題確定求解域的幾何區(qū)域；將求解域劃分為有限個網格單元的離散域(網格單元越多，計算結果越精確)；推導有限單元的列式，建立單元試函數，進而形成單元矩陣；采取迭代法或者直接法等對聯(lián)立方程組進行求解，獲得單元節(jié)點處狀態(tài)變量近似值。單元的各種性質是根據單元位移來建立的，假設單元位移分量由坐標線性函數表示，單元應變?yōu)椋?/p>

{ε}=[εxεyεxy]T

(1)

根據廣義的胡克定律，針對線彈性提來說，單元應變分量{ε}與單元應力分量[σ]存在線性關系：

[σ]=[D]{ε}=[D][B]{δ}e

(2)

在公式(2)中，[D]為彈性矩陣，僅與材料性質有關，在平面應力狀況下為：

(3)

{F}e=?[B]T{δ}tdxdy=?[B]T[D][B]tdxdy{δ}e=[k]{δ}e

(4)

其中，[K]為單元剛度矩陣，結合公式(4)可知

[K]=?[B]T[D][B]tdxdy

(5)

單元分析完畢之后，即可結合結構力學方法組合所有單元特征，推導出整個單元體系的性質，這一過程被稱為總體分析。結合公式(5)，將結點i的平衡方程，集合結構各個結點對應的平衡方程，得到整體結構的平衡方程為：

[K]{δ}={R}

(6)

采用有限元法進行邊坡穩(wěn)定的分析與傳統(tǒng)極限平衡法相比，主要具備以下優(yōu)勢：能夠充分考慮到土體本構的非線性關系作為邊坡穩(wěn)定性分析的基礎；可適應復雜多變的邊坡條件，例如存在結構物作用以及支護、幾何形狀復雜的邊坡。

2 基于有限元法的安全系數計算

結合寶雞市一在建民用住宅樓護坡支護工程實際情況，根據經驗可知邊坡發(fā)生失穩(wěn)破壞的主要方式為在塑性區(qū)形成一個滑裂帶(見圖1)。

圖1 邊坡失穩(wěn)滑裂帶分布Fig.1 Distribution of slope failure zone

(7)

結合摩爾庫倫準則可知，抗剪強度τf與正應力σ間存在如下關系：

τf=σtanφ+c

(8)

因此整個滑動面安全系數為：

(9)

結合有限元原理可知，邊坡破壞是一個漸進的過程，破壞區(qū)域從一點開始面，逐步拓展到全局，土體應力時刻變化，直到失穩(wěn)破壞，在邊坡沿某一滑動面達到極限平衡條件發(fā)生破壞時，對滑動面上每一點均有：

τ=τf/F

(10)

圖2 微元體受力圖Fig.2 Infinitesimal body diagram

根據力的平衡條件可知

Ti+Ri=0

(11)

即有：

(12)

其中Δli為微元曲線長度，那么結合微積分基本思想可知，針對整個曲面則有：

(13)

進一步可得：

F=∫lτfdl/∫lτdl

(14)

在公式(14)中F即為安全系數。

針對本研究實際問題，假設如圖3所示邊坡土體ABC為脫離體。

圖3 滑動體受力示意圖Fig.3 Schematic diagram of force of sliding body

在圖3中，邊坡應力分布為σij(x,y)，滑裂面正應力為σ，剪應力為τ，脫離體受到外力σ與τ合力平衡，則在滑裂面上可由摩爾—庫倫準則得到如下公式：

(15)

(16)

3 算例分析

3.1 計算流程設計

選取寶雞市某一施工基坑，支護方式為護坡支護，具體力學參數為：土容重γ=19.60kN/m3、坡角θ=45°、坡高h=20m、內摩擦角φ=20°、彈性模量E=2.8×107Pa、粘聚力c=22.5kPa、泊松比μ=0.3，坡體底邊固定約束，左右邊界水平約束，見圖4。

圖4 約束條件示意圖Fig.4 Sketch map of constraint condition

結合有限元原理，研究按照較密的劃分原則對圖4進行網格劃分(見圖5)，結點總數3 925個，單元總數為2 431個。

圖5 網格劃分Fig.5 Grid division

借助MATLAB8.0軟件以及ANSYS9.0軟件對其進行計算，計算流程見圖6。

圖6 滑裂面與安全系數計算流程圖Fig.6 Flow chart of calculation of slip surface and safety factor

3.2 計算結果與結構分析

3.2.1 滑裂面位置確定以及滑裂面切應力分析

在有限元計算過程中，自重等級分30個等級勻速加載，結合ANSYS9.0軟件計算塑性區(qū)塑性應變等值曲線(見圖7)，其中紅線位置為邊坡最危險滑動面，并得到最危險滑動面的擬合曲線(圖8(A))以及各節(jié)點對應的切應力數據(圖8(B))。

圖7 塑性應變等值線與最危險滑動面Fig.7 Plastic strain contour and the most dangerous slip surface

圖8 滑裂面曲線擬合以及滑裂面切應力擬合Fig.8 Curve fitting of slip surface and fitting of the shear stress on the slip surface

由圖8(A)可知滑裂面近似于圓弧，與滑裂面為圓弧假設具有一致性，結合圖8(B)的擬合曲線可算出滑裂面上的每一處所對應的切應力，其最大切應力大概在結點4、5、6之間，因此在實際工程護坡設計中，要保障最大切應力處符合設計要求，以免發(fā)生護坡失穩(wěn)滑裂破壞。

3.2.2 安全系數的計算及對比分析

影響邊坡穩(wěn)定的參數有許多，分別選取2個參數為變量(坡角θ、粘聚力c)，根據圖6中的程序流程圖，進行ANSYS程序設計，計算了每個參數變化情況下的安全系數值，并與傳統(tǒng)的瑞典法、畢肖普法、簡布法進行對比分析，結果見表1和表2。

表1 坡角θ為變量(H=20m、φ=20°、c=22.6kPa)

表2 粘聚力c為變量(H=20m、φ=20°、θ=45°)

由表1和表2可知：幾種方法得出的安全系數變化趨勢都是一致的，說明了該研究計算方法科學且有效。此外，采取有限元法尋找滑裂面得到的安全系數比傳統(tǒng)的瑞典法、畢肖普法、簡布法計算所得安全系數值要大，究其緣由，可能是傳統(tǒng)的瑞典法、畢肖普法以及簡布法加入了過多假設，鋼化了坡體結構，造成計算結果偏于保守，而采取有限元法尋找滑裂面由于更加接近實際工程情況，從而使得計算結果更加精確。

4 結論

1)滑裂面切應力近似于拋物線(見圖8)，最大切應力產生在滑裂面中間位置(節(jié)點4、5、6)，在實際工程中，若護坡底部產生細微隆起，實際產生滑動位置為護坡中部位置，要注意對護坡中部加強錨桿鎖固，加大約束力度。

2)有限元法研究安全系數偏大，這是由于有限元法更加貼近實際，傳統(tǒng)方法鋼化了坡體結構。因此，在實際護坡支護過程中，可根據有限元計算結果適當簡化支護方案。

3)目前有限元原理在建筑基坑邊坡支護穩(wěn)定性的計算方面還不夠成熟，本研究將有限元法融入到邊坡支護的計算中，可為日后的支護方案設計提供了新的計算思路。

[1] 趙杰.邊坡穩(wěn)定有限元分析方法中若干應用問題研究[D/OL].大連理工,2006:28[2016-05-20].http://www.cnki.net/KCMS/detail/detail.aspx?.

[2] 常崇義.有限元輪軌滾動接觸理論及其應用研究[D/OL].中國鐵道科學研究院,2010:10-14[2016-5-16].http://www.cnki.net/KCMS/detail/detail.aspx?.

[3] 鄧岳保.豎井地基固結解析理論與有限元分析[D/OL].浙江大學,2013:37-39[2016-5-12].http://www.cnki.net/KCMS/detail/detail.aspx?.

[4] 趙建軍.公路邊坡穩(wěn)定性快速評價方法及應用研究[D/OL].成都理工大學,2007:12-13[2016-4-26].http://www.cnki.net/KCMS/detail/detail.aspx?.

[5] 倪衛(wèi)達.基于巖土體動態(tài)劣化的邊坡時變穩(wěn)定性研究[D/OL].中國地質大學,2014:46-49[2016-5-18].http://www.cnki.net/KCMS/detail/detail.aspx?.

[6] 孫志彬.邊坡穩(wěn)定性上限分析方法及參數反演研究[D/OL].中南大學,2013:36-38[2016-5-25].http://www.cnki.net/KCMS/detail/detail.aspx?.

[7] 吳振君,王水林,湯華,等.一種新的邊坡穩(wěn)定性因素敏感性分析方法—可靠度分析方[J].巖石力學與工程學報,2010,29(10):2050-2055.

[8] 黃潤秋.巖石高邊坡發(fā)育的動力過程及其穩(wěn)定性控制[J].巖石力學與工程學報,2008,27(8):1525-1544.

[9] 鄧東平,李亮.基于滑動面搜索新方法對地震作用下邊坡穩(wěn)定性擬靜力分析[J].巖石力學與工程學報,2012,31(1):86-98.

[10]HISHAM T,EID K H,RABIE D W.Drained residual shear strength at effective normal stresses relevant to soil slope stability analyses[J].Engineering Geology,2016,38(11):106-112.

[11]LARISA B L.Domain decomposition finite element/finite difference method for the conductivity reconstruction in a hyperbolic equation[J].Communications in Nonlinear Science and Numerical Simulation,2016,5(12):83-86.

[12]DANIEL J E,ELENA B F,JAIME C G,et al.Decision aid system founded on nonlinear valuation,dispersion-based weighting and correlative aggregation for wire rope selection in slope stability cable nets[J].Expert Systems With Applications,2016,76(5):121-126.

[13]VAFAEIAN B,LE L H,TRAN N H T,et al.Micro-scale finite element modeling of ultrasound propagation in aluminum trabecular bone-mimicking phantoms:A comparison between numerical simulation and experimental results[J].Ultrasonics,2016,4(7):35-37.

[14]HONG F F,GUO H,HOU J,et al.A stabilized mixed finite element method for steady and unsteady reaction-diffusion equations[J].Computer Methods in Applied Mechanics and Engineering,2016,48(5):52-58.

Stability calculation of slope support in building foundation pit based on finite element theory

LIU Zhanhong

(Shaanxi institute of Technology Architecture and Thermal Energy Branch, Xi'an, Shaanxi 710302,China)

Based on a practical slope supporting engineering for finite element modeling, we applied MATLAB8.0 software and ANSYS9.0 software to found slip surface, calculated slip surface shear force value and maximum shear force position. Then, calculated supporting the stability of the slope safety coefficient by the finite element method, and compared with the traditional method. The results showed that the safety factor calculated by the finite element method was larger, and the change trend of the safety factor of several methods was the same.

finite element; safety factor; slip surface; stability

1004—5570(2016)06-0093-05

2016-06-04

劉占宏(1977-),男,講師,研究方向:土木工程技術,E-mail:137477424@qq.com.

TU470

A