李榮華

（楊凌職業(yè)技術(shù)學院陜西咸陽712100）

改進Mohr-Coulomb算法在湖泊地區(qū)基坑開挖中的應用

李榮華

（楊凌職業(yè)技術(shù)學院陜西咸陽712100）

針對湖泊地區(qū)的復雜富水軟土地基坑開挖施工時經(jīng)常發(fā)生失穩(wěn)坍塌問題。利用最大拉應力準則分析巖土介質(zhì)的抗張強度以此改進Mohr-Coulomb算法策略，并運用該算法策略對湖泊地區(qū)基坑開挖鉆孔壁穩(wěn)定性進行分析，結(jié)合FLAC3D數(shù)值仿真模擬軟件，重點分析了不同地下水深、穩(wěn)定液比重、開挖直徑條件下的基坑開挖的變形規(guī)律進行了一系列的數(shù)值模擬研究與評價，結(jié)果表明：基坑開挖鉆孔深度達到10m至15m之間時，土層孔壁的橫向變形將顯著增大。該改進算法策略為湖泊地區(qū)基坑開挖的設(shè)計與施工提供參考。

基坑開挖；Mohr-Coulomb算法；塌陷問題；穩(wěn)定性

在湖泊地區(qū)廣泛分布著河相沉積的軟弱黏性土層［1］。含水量大、強度低、壓縮性高、透水性差、埋藏深厚是該類土的特點，在荷載作用下會產(chǎn)生相當大的沉降和沉降差，其承載力和穩(wěn)定性也往往不能滿足工程要求［2］?；娱_挖以承載力高、無擠土、施工時無振動、噪音小和宜于在城市建筑物密集地區(qū)使用等優(yōu)點，被廣泛應用于湖泊地區(qū)現(xiàn)代基礎(chǔ)工程建設(shè)中［3-5］。然而，基坑開挖失穩(wěn)問題是施工過程中普遍存在并一直備受困擾的重大工程問題，在施工過程中稍有不慎，極易發(fā)生塌陷等工程質(zhì)量事故［6］?；娱_挖穩(wěn)定性問題的出現(xiàn)在很大程度上影響施工的進度、增加施工成本、降低樁基的承載力，嚴重制約著開挖的應用與發(fā)展。

Mohr-Coulomb強度理論能較好地描述巖土材料的強度特性和破壞行為，然而Mohr-Coulomb屈服面在主應力空間是一個六角錐面，角點處的導數(shù)不連續(xù)，給計算帶來了不便，導致Mohr-Coulomb算法不便應用到有限元數(shù)值計算中［7］。本研究結(jié)合拉伸型和剪切型Mohr-Coulomb屈服準則［8］，運用最大拉應力準則改進Mohr-Coulomb算法，并以FLAC3D

1　改進的Mohr-Coulomb算法

以應力不變量形式表示的剪切型Mohr-Coulomb屈服準則可以表示為［9］：

式中：c、φ分別為粘聚力和內(nèi)摩擦角；σm=I1/3為平均應力；為等效應力；I1、J2分別為應力第1不變量和應力偏量第2不變量；θ為Lode角；本研究規(guī)定拉應力為正，壓應力為負。其中，

以應力不變量形式表示的拉伸型Mohr-Coulomb屈服準則可以表示為：

式中：ft為巖土介質(zhì)的抗拉強度。拉伸型Mohr-Coulomb屈服準則是一個等邊三角形，在主應力空間屈服面由3個分別垂直于主應力軸的平面組成［10］；剪切型Mohr-Coulomb屈服準則是一個不等角的六邊形，在主應力空間為一個棱錐面，中心軸線與等傾線重合［11］。

采用雙曲線方程對拉伸型Mohr-Coulomb屈服準則和剪切型Mohr-Coulomb屈服準則進行擬合。通過調(diào)整參數(shù)m的大小來反映巖土介質(zhì)的抗張強度的大小，取值范圍為0≤m≤1。當m=0時，復合屈服準則退化為剪切型Mohr-Coulomb屈服準則；當m=1時，巖土介質(zhì)的抗張強度為0。同時，參數(shù)m可以修正屈服面上的尖頂，使尖角變的光滑，避免了數(shù)值計算的發(fā)散和收斂的緩慢。

由于Mohr-Coulomb屈服面存在6個棱角這些奇異點，使數(shù)值計算變繁和收斂緩慢［12］。為了盡量逼近和接近屈服面，使得改進后的屈服面盡量接近Mohr-Coulomb屈服面，并且在棱角處得到光滑連續(xù)的處理。則改進的Mohr-Coulomb屈服準則的表達式為

2　數(shù)值計算分析模型的建立

在建立開挖模型中，模擬實際開挖，不考慮施工機械或人為因素對基坑壁穩(wěn)定性的影響，地面無超載，土體模型均采用改進的Mohr-Coulomb模型，使用非對稱求解器［13］進行求解，根據(jù)最大拉應力準則，調(diào)整參數(shù)m=0.5以平衡巖土介質(zhì)之間的抗張強度。假設(shè)各層土體土質(zhì)為均質(zhì)各向同性的彈塑性體，目的是為了反映出一些具有普遍性的規(guī)律。

2.1數(shù)值模擬方案和基本條件的設(shè)定

為了深入分析基坑開挖穩(wěn)定性在實際開挖施工不同工況條件下的變化規(guī)律，采用數(shù)值軟件FLAC3D［14］，在不同地下水埋深、穩(wěn)定液比重、開挖直徑、開挖間距的施工工況下模擬開挖，針對工程施工區(qū)域的地質(zhì)背景，具體設(shè)計了四種不同工況下開挖模擬方案，如表1所示。

表1　模擬方案設(shè)定

在對基坑開挖工程區(qū)域進行基坑數(shù)值模擬分析時，計算模型均采用三維計算模型，采用模型長寬均為15m，大于十倍的開挖直徑以此來消除邊界約束對基坑壁穩(wěn)定性的影響，深度55m（大于開挖深度加十倍的開挖半徑半徑），為了更好地反映開挖的實際受力狀態(tài)和變形規(guī)律，充分考慮邊界條件并將其邊界條件約束如下［15］：1）單約束開挖模型的左右邊界，取u1=0，u2≠0，u3≠0（u1為X軸方向的位移，u2為Y軸方向的位移，u3為z軸方向的位移）；2）單約束開挖模型底邊邊界，取u1≠0，u2≠0，u3=0；3）單約束開挖模型前后邊界，取u1≠0，u2=0，u3≠0；4）初始基坑水壓力場的設(shè)定，根據(jù)工程實際施工地下水位，設(shè)定地下水位埋深1.5m處初始孔隙水壓力為0；在模擬降水施工時設(shè)定降水一側(cè)同降水深度的基坑水壓力場為0；5）模型的上邊界不予約束為自由邊界。

2.2模型計算參數(shù)的設(shè)定

文中選取具有代表性的參數(shù)，地層主要由第四系人工填土層（Q4ml）、第四系全新統(tǒng)的風-水堆積層（Q4eol+m）、沼澤相沉積層（Q4h）、海陸相交互沉積層（Q4mc）、第四系上更新統(tǒng)的海陸相交互沉積層（Q3mc）、沖、洪積層（Q3al+pl）、殘積層（Q3el）以及燕山期花崗巖（γ52-3）構(gòu)成［16］。由于該場區(qū)地層的特殊性和項目裝置的特殊要求，場區(qū)重要裝置及管廊全部采用樁基基礎(chǔ)。但由于該工程區(qū)地基特殊的松散地層地質(zhì)條件與地下水動力等因素，在基坑開挖鉆孔灌注樁施工中，常常出現(xiàn)比較大的鉆孔孔壁變形現(xiàn)象，有時也會演變成較為嚴重塌孔問題。

因此，蠟油加氫污水提升裝置的灌注樁基坑開挖施工區(qū)域進行系統(tǒng)的模擬分析。該聯(lián)合區(qū)域內(nèi)，地層條件復雜多變，砂層厚度較厚，而且在施工中部分基坑已經(jīng)出現(xiàn)塌孔現(xiàn)象。根據(jù)巖土體地層結(jié)構(gòu)分布特點，并結(jié)合勘察報告內(nèi)的場區(qū)主要土層分布圖和各土層等高線圖，將周圍的巖土體作為主要的研究對象，并把復雜的巖土體性質(zhì)、材料參數(shù)等單一化和理想化，得到鉆孔灌注樁基坑壁穩(wěn)定性分析的地層分布簡圖，如圖1所示。

圖1　數(shù)值模擬底層分布簡圖

假設(shè)各地層的物理力學參數(shù)詳見表2。

2.3計算單元網(wǎng)格劃分

在研究模型中，巖土體均采用八節(jié)點六面體等參實體單元，為了更好地模擬研究鉆孔灌注開挖基坑壁的變形情況及穩(wěn)定性，將鉆孔附近的單元進行局部加密，圖4為鉆孔三維模型單元網(wǎng)格劃分圖。

表2　地層參數(shù)

圖2　單元劃分的三維模擬

3　基坑開挖穩(wěn)定性的數(shù)值模擬分析

由于該研究區(qū)域的灌注樁基坑開挖工程施工現(xiàn)場臨近湖泊，其地下水埋深較淺，地下水對孔壁土體的靜水壓力作用較大，為了分析地下水埋深對鉆孔開挖基坑壁穩(wěn)定性的影響規(guī)律，文中基于研究區(qū)域灌注樁基坑開挖工程施工場地的地質(zhì)背景，針對不同地下水埋深條件下進行模擬鉆孔并分析其基坑壁的穩(wěn)定性。模型建立以后，首先加載重力并設(shè)置3種地下水埋深條件分別為0.5m、1.5m、2.0m。然后分別在3個模型中開挖直徑是800 mm的鉆孔基坑開挖，鉆孔基坑開挖的設(shè)計深度為45m，采用比重為1.15的護壁穩(wěn)定液，分別得出地下水埋深為0.5m、1.5m、2.0m的工況下基坑開挖鉆孔壁橫向位移變形云圖，如圖3至圖5所示。

圖3　水深0.5m時孔壁橫向位移

圖4　水深1.5m時孔壁橫向位移

圖5　水深2m時孔壁橫向位移

由圖3至圖5可見，不同地下水埋深件下基坑開挖鉆孔孔壁的橫向位移云圖可以發(fā)現(xiàn)，在該研究區(qū)域地層條件下模擬施工鉆孔孔壁呈現(xiàn)非線性的變形，在砂土層及其夾層中基坑壁的橫向變形顯著大于其他土層，且在其他施工條件一定的情況下，隨著地下水埋深的增加，基坑壁變形總體上有減小的趨勢。根據(jù)改進的Mohr-Coulomb屈服準則，土體模型在拉應力的作用下，基坑側(cè)面對開挖鉆孔產(chǎn)生橫向壓力，得到不同穩(wěn)定液比重條件下基坑開挖鉆孔壁橫向位移值，如表3所示。

表3　不同地下水埋深條件下基坑開挖鉆孔橫向位移

由表3可見，在基坑開挖鉆孔深度10m至15m之間的土層孔壁的橫向變形顯著增大，而該土層區(qū)間內(nèi)主要成分是砂土，由此可見砂土層對于鉆孔施工的擾動極其敏感。砂土具有獨特的物理力學特征，其粘聚力的理論值為零，實際工程地質(zhì)條件下砂土成份中含有一定量的黏土顆粒，對砂土顆粒起到一定的膠結(jié)作用，但仍然不足以提高其土體結(jié)構(gòu)的整體穩(wěn)定性。因此，開挖鉆孔施工引起的重分布應力作用于砂土層時其土體結(jié)構(gòu)容易發(fā)生失穩(wěn)，在數(shù)值仿真模型的基坑開挖鉆孔中表現(xiàn)為孔壁向臨空面的橫向位移變形；而在實際的灌注樁基坑開挖施工鉆孔中表現(xiàn)為該土層孔壁土體的脫落、滑塌現(xiàn)象，甚至引起整個開挖的塌陷。由此可見，砂土層的穩(wěn)定對于灌注樁基坑開挖施工孔壁的穩(wěn)定性起著舉足輕重的作用，在含砂土層中灌注樁基坑開挖鉆孔施工應予以足夠的重視，采取相應防治措施，有針對性的防治灌注樁基坑鉆孔施工基坑壁失穩(wěn)問題。

隨著地下水埋深的增加，基坑開挖孔壁位移變形總體上呈現(xiàn)減小的趨勢，在地下水位埋深分別為0.5m、1.5m、2m的條件下，數(shù)值模擬孔壁的橫向位移最大值分別為15.37 mm、14.19 mm和13.67 mm，可見地下水埋深越淺，對灌注樁基坑開挖施工孔壁的穩(wěn)定性越不利。地下水埋深即反映地下水位的高低，根據(jù)鉆孔灌注樁孔壁穩(wěn)定性力學評價模型可知，地下水對灌注樁基坑開挖施工孔壁產(chǎn)生靜壓力作用。根據(jù)液體壓強的基本原理，地下水位越高，其對與基坑開挖孔壁土體產(chǎn)生的壓強越大，即地下水對基坑壁引起較大的致塌力。因此地下水埋深越淺，基坑壁土體越不穩(wěn)定。因此，在含水豐富的地基土中進行灌注樁基坑開挖鉆孔施工前，對施工區(qū)域采取一定的降水措施降低地下水位，即增加地下水埋深，對于開挖施工基坑壁失穩(wěn)問題的防治具有積極的作用。

4　結(jié)論

本研究運用最大拉應力準則改進Mohr-Coulomb算法，以FLAC3D數(shù)值仿真模擬軟件為平臺，針對湖泊地區(qū)地質(zhì)分層和地下水埋深的特點，采用灌注樁基坑開挖鉆孔進行施工。利用FLAC3D數(shù)值分析軟件對巖石、土層和支護結(jié)構(gòu)等建立高級三維模型，運用改進的Mohr-Coulomb算法計算出不同地下水埋深條件下基坑開挖鉆孔橫向位移，并重點分析了湖泊地區(qū)的復雜富水軟土地基坑開挖過程中，不同地下水深、穩(wěn)定液比重、開挖直徑條件下的基坑開挖的變形規(guī)律進行了一系列的數(shù)值模擬，為湖泊地區(qū)基坑開挖的設(shè)計與施工提供參考。

［1］張燦，周愛鋒，張曉楠，等.湖泊沉積記錄的古洪水事件識別及與氣候關(guān)系［J］.地理科學進展，2015，34（7）:898-908.

［2］彭霜霜，王志銳.東非高原沼澤勘察方法研究與成果分析［J］.中外公路，2015（1）:18-23.

［3］張陳蓉，俞劍，黃茂松.基坑開挖對鄰近地下管線影響的變形控制標準［J］.巖土力學，2012，33（7）:2027-2034.

［4］陳國良，劉修國，尚建嘎，等.三維地質(zhì)結(jié)構(gòu)模型的切割分析技術(shù)及方法［J］.計算機工程，2007，33（20）:184-186.

［5］奚雪峰，班建民，付保川.基坑支護工程的三維模擬仿真系統(tǒng)研究［J］.計算機工程與應用，2011，47（12）:213-216.

［6］李文科，張研，張寶民，等.川中震旦系—二疊系古巖溶塌陷體成因、特征及意義［J］.石油勘探與開發(fā)，2014，41（5）:513-522.

［7］李翠華，姜清輝，周創(chuàng)兵，Mohr-Coulomb準則角點問題的主應力空間互補算法［J］.工程力學，2014（4）:134-140.

［8］WANG Fei，HUANG Xing-chun.Analysis of ThreedimensionalCrackPropagationbyUsingDisplacement Discontinuity Method［J］.Journal of Donghua University：English Edition，2010，27（6）:835-840.

［9］賈善坡，陳衛(wèi)忠，楊建平，等.基于修正Mohr-Coulomb準則的彈塑性本構(gòu)模型及其數(shù)值實施［J］.巖土力學，2010，31（7）: 2051-2058.

［10］范雙雙，楊燦軍，彭時林，等.水下滑翔機關(guān)鍵承壓系統(tǒng)設(shè)計與試驗研究［J］.浙江大學學報（工學版），2014（4）:633-640.

［11］李慧，楊雪，趙廣生.基于ABAQUS柔性模型的ADAMS模型研究［J］.科技創(chuàng)新與應用，2012（11）:20-21.

［12］孫偉，張琦，孫豐，等.工程機械虛擬訓練動態(tài)地形仿真［J］.系統(tǒng)仿真學報，2009（6）:1776-1780.

［13］趙永華，遲學斌，王武.HPSEPS軟件包及其千核應用［J］.華中科技大學學報（自然科學版），2011，39（1）:59-62.

［14］姜慶紅，范慧鶴，李紅英，等.FLAC3D中INTERFACE建模新方法在采空區(qū)模擬中的應用［J］.計算機與現(xiàn)代化，2008（10）:113-115.

［15］田樹昆，曹蘭柱.基于CAD與ANSYS的FLAC3D邊坡模擬分析［J］.微計算機信息，2010，26（33）:259-260.

［16］徐雪戰(zhàn)，孟祥瑞，何葉榮，等.基于三維可視化與虛擬仿真技術(shù)的綜采工作面生產(chǎn)仿真研究［J］.中國安全生產(chǎn)科學技術(shù)，2014（1）:26-32.

Study on the improved Mohr-Coulomb algorithm for excavation of foundation pit in the lake area

LI Rong-hua
（Yangling Vocational&Technical College，Xianyang 712100，China）

In view of the lake district rich water complex in construction of soft soil foundation pit excavation often buckling collapse problem.Using maximum tensile stress criterion analysis the tensile strength of rock and soil medium so as to improve the Mohr-Coulomb algorithm strategy，and the strategies of using the algorithm to the lake district of borehole wall stability of foundation pit excavation is analyzed，combining with FLAC3D numerical simulation software，with different groundwater depth，stable liquid proportion，the excavation diameter under the condition of the deformation law of foundation pit excavation research and evaluation，conducted a series of numerical simulation results show that the excavation hole depth up to 10m to 15m，between the lateral deformation of the soil hole wall will increase significantly.The lake district of the strategy of the improved algorithm provides reference for design and construction of foundation pit excavation.

excavation of foundation pit；Mohr-Coulomb algorithm；collapse problem；stability

TN473.2

A

1674－6236（2016）17-0011-04

2016-03-08稿件編號：201603098

國家自然科學基金項目（41301442）；楊凌職業(yè)技術(shù)學院科學研究基金項目（A2013051）

李榮華（1986—），女，山東鄆城人，碩士，講師。研究方向:道路與鐵道工程。