李飛,2，孫俊2，張飛*，何淼，許龔成

(1.鹽城工學(xué)院土木工程學(xué)院， 江蘇鹽城224001; 2.蘇州科技大學(xué)土木工程學(xué)院， 江蘇蘇洲215009)

0 引言

隨著城市建設(shè)的不斷發(fā)展，高層、超高層建筑以及地下工程層出不窮，伴隨而來的深基坑工程也越來越多。同時(shí)，城市土地資源日趨緊張，新建基坑不得不緊鄰已有建(構(gòu))筑物或緊鄰河流湖泊，此時(shí)基坑一側(cè)為有限土體，若仍按經(jīng)典土壓力理論來計(jì)算此時(shí)作用在支護(hù)上的土壓力，可能偏于浪費(fèi)[1]。

目前國(guó)內(nèi)外已有一些學(xué)者針對(duì)有限土體土壓力展開研究，采用的方法主要為極限平衡法和薄層單元法[2-4]。胡明義[5]分析了擋土墻后有限土體破裂面在不受限、部分受限以及完全受限三種情況下土壓力的分布。高印立[6-7]基于塑性上限理論求得有限土體主被動(dòng)土壓力的上限解，并與朗肯土壓力的計(jì)算值比較。王文杰等[8]利用微分方程，較為系統(tǒng)地考慮了墻背黏聚力，阻力等因素，建立平衡方程，并通過簡(jiǎn)單算例進(jìn)行驗(yàn)證。李峰等[9]考慮了滑動(dòng)楔形體與墻背和土體之間的粘著力，通過滑裂體的極限平衡狀態(tài)建立方程并求解，其結(jié)果與朗肯土壓力理論值相比較，認(rèn)為有限土體成立的條件為寬深比在0.5～0.75之間。王洪亮等[10-12]推導(dǎo)的計(jì)算方法考慮了墻土界面摩擦，同時(shí)還考慮了粘土的黏聚力過大和有限土體寬度過小的特殊情況，推導(dǎo)出對(duì)應(yīng)的計(jì)算方法。王閆超等[13]依據(jù)薄層單元法，對(duì)隔離體進(jìn)行受力分析，推導(dǎo)出主動(dòng)土壓力合力及作用位置，并分析了破裂角θ與有限土體寬度的關(guān)系。劉鑫[14]對(duì)考慮放坡情況下的有限土體進(jìn)行了數(shù)值模擬，模擬了當(dāng)擋墻垂直和擋墻角度為75°時(shí)土壓力的分布。胡敏云等[15]針對(duì)深基坑排樁支護(hù)有限土壓力展開研究，指出排樁嵌入深度與土壓力的關(guān)系。

目前大多數(shù)研究均是針對(duì)擬開挖基坑緊鄰已有建筑的情況[16-17]，對(duì)基坑臨河情況下的土壓力研究甚少。我國(guó)東部大部分城市，市內(nèi)河流縱橫，開挖基坑往往不得不緊鄰河流，一些建筑為了景觀效應(yīng)，也會(huì)依河、湖而建。工程上大多采用半經(jīng)驗(yàn)半理論的方式來計(jì)算此時(shí)基坑支護(hù)的土壓力，難以進(jìn)行驗(yàn)證與預(yù)測(cè)基坑開挖形狀。本文擬針對(duì)臨河情況下，對(duì)黏性土基坑與河流之間形成的有限土體進(jìn)行受力分析，建立方程，得出其主動(dòng)土壓力的計(jì)算方法，為實(shí)際工程中臨河基坑支護(hù)土壓力的計(jì)算提供參考。

1 計(jì)算模型

1.1 計(jì)算假定與考慮因素

有限土體是指當(dāng)擋墻與河流岸邊或已有建筑基礎(chǔ)的距離小于一定值時(shí)，土體發(fā)生滑裂的破裂面被河流或者已有建筑基礎(chǔ)所截?cái)啵荒芡耆l(fā)展，如圖1所示。本文截取有限寬度土體，對(duì)其滑裂土體進(jìn)行受力分析，此有限土體的滑動(dòng)模式已被各學(xué)者廣泛采用。為了簡(jiǎn)化計(jì)算，本文僅考慮滑裂土體與下部不動(dòng)土體之間的摩擦，不考慮墻土界面摩擦，即不考慮擋土墻上的切向力；不考慮河床的影響，即認(rèn)為河岸是垂直于地表的；針對(duì)坑底土體為滲透性低的黏土的情況，將河流及地下水位對(duì)有限土體的作用力視為靜止水壓力，忽略滲流與動(dòng)水壓力的影響，采用水土合算法。本文計(jì)算基本假設(shè)條件如下：

①擋土墻背面光滑，與土體之間沒有摩擦；

②有限土體為單層各向同性的c，φ均勻的粘性土；

③河岸與地平面垂直，即不考慮河床的影響；

④河流足夠深，破裂面不會(huì)被河底以下所截?cái)唷?/p>

1.2 受力分析與方程的建立

受力簡(jiǎn)圖如圖2所示，有限土體處于極限平衡狀態(tài)時(shí)的土壓力即為主動(dòng)土壓力，此時(shí)即將發(fā)生滑裂的土體為楔形體ABCD，基坑開挖深度為H，有限土體寬度b，滑裂面與水平方向的夾角為θ，土的黏聚力為c，內(nèi)摩擦角為φ，河面與地面的距離為h。下部不動(dòng)土體對(duì)滑裂體的作用力為R。

滑動(dòng)土體每延度重力：

(1)

CD面上的黏聚力：

(2)

圖1 形成臨河有限土體的基坑
Fig.1 Pit near the river

圖2 滑裂體受力簡(jiǎn)圖
Fig.2 Stress diagram of slipped soiltake shape the limited soil

河流中水的作用范圍為：

CE=H-h-btanθ。

(3)

河流中水對(duì)有限土體的水平向作用力：

(4)

對(duì)水平向和豎直向建立平衡方程：

(5)

將W，C，F(xiàn)代入上式：

(6)

通過上述公式求Ea，一種方法是Ea對(duì)θ求導(dǎo)，求出Ea的極值，即為有限土體的主動(dòng)土壓力，另一種是試算法，將不同的θ(0～90°)代入式中，取得的最大的Ea值即為主動(dòng)土壓力。本文利用MATLAB軟件將Ea對(duì)θ求導(dǎo)：

(7)

當(dāng)dEa/dθ=0時(shí)，此時(shí)的θ即為要求的破裂角，由上式可知，θ是一個(gè)與重度γ、開挖深度H、黏聚力c、內(nèi)摩擦角φ、河面離地面的深度h、有限土體寬度b有關(guān)的復(fù)雜函數(shù)，并不是朗肯土壓力理論中認(rèn)為的(45+φ/2)，很難求出θ顯式。通過室內(nèi)土工試驗(yàn)及現(xiàn)場(chǎng)測(cè)量，能得到相關(guān)參數(shù)，將其代入dEa/dθ=0，利用MATLAB、EXCEL等軟件便能求出對(duì)應(yīng)的θ。

但上式是針對(duì)有限土體的破裂面是被河流截?cái)嗟那闆r，如圖2所示，此時(shí)0

1.3 特殊情況分析

以上計(jì)算是針對(duì)滑裂土體為楔形體且破裂面是被河流水下部分所截?cái)啵?dāng)有限土體寬度b足夠小或開挖深度H偏大時(shí)，可能會(huì)出現(xiàn)滑裂土體截面為三角形或者破裂面是被河面上部與地面之間的土體所截?cái)嗟膬煞N特殊情況，如圖3所示：

(a) 破裂面在地面上

(b) 破裂面在水面以上地面以下

圖3(a)所示情形，滑裂土體截面為三角形，適用條件為tanθ>H/b，此時(shí)：

對(duì)水平向和豎直向建立平衡方程：

(8)

將W，C代入式(8)，得：

(9)

利用MATLAB將Ea對(duì)θ求導(dǎo)，得：

(10)

同樣求得當(dāng)dEa/dθ=0時(shí)，此時(shí)的θ為極限平衡狀態(tài)下的破裂角，代入式(9)，即可求出主動(dòng)土壓力。

圖3(b)所示情形，土體的滑裂面并不是被河流水面以下所截?cái)?，而是被河面之上，地面以下這部分所截?cái)?，適用條件為(H-h)/b

每沿度重力：

(11)

CD面上的黏聚力：

(12)

對(duì)水平向和豎直向建立平衡方程：

(13)

將W、C代入式(13)：

(14)

利用MATLAB將Ea對(duì)θ求導(dǎo)，得：

(15)

1.4 計(jì)算方法小結(jié)

如上所述，當(dāng)破裂角θ不同時(shí)，采用的有限土體主動(dòng)土壓力的計(jì)算方法也是不同的，具體如下：

根據(jù)式(7)、(10)、(15)對(duì)滑動(dòng)土體破裂角進(jìn)行試算，然后根據(jù)tanθ在上式中所在的范圍，選擇對(duì)應(yīng)的計(jì)算公式，求得其主動(dòng)土壓力。

2 算例分析

假設(shè)擬開挖基坑臨河距離為6 m，基坑開挖深度為20 m，土層為均質(zhì)的粉質(zhì)粘土，其重度γ=18 kN/m3,黏聚力c=10 kPa，內(nèi)摩擦角φ=20°。河流深5 m，河面離地面距離為2 m。

例如要計(jì)算基坑開挖到10 m深的時(shí)候的有限土體主動(dòng)土壓力，先將各參數(shù)代入式(7)，利用MATLAB軟件求得當(dāng)dEa/dθ=0時(shí)，θ為53.86°，tanθ=1.37?？雌涫欠駶M足破裂面被河面以下所截?cái)嗟臈l件：0

分別按上述過程計(jì)算在各開挖深度時(shí)的主動(dòng)土壓力合力，各深度土壓力的分布可以用主動(dòng)土壓力合力的差分形式算出。當(dāng)開挖深度過低時(shí)，b/H過大，不構(gòu)成有限土體，文獻(xiàn)[5]認(rèn)為當(dāng)b/H>1，破裂角θ<45°，有限土體計(jì)算不適用，這也與本文計(jì)算方法吻合。表1列出了b/H<1時(shí)本文方法與朗肯土壓力理論在各深度的土壓力數(shù)值：

表1 各深度土壓力Tab.1 Earth pressureat different depths

將本文的計(jì)算方法與朗肯土壓力公式得出的數(shù)據(jù)作比較，如圖4所示：

從圖4中可以看出有限土體土壓力呈非線性分布，總體小于朗肯土壓力，隨著深度增大，有限土體土壓力和朗肯土壓力的差別越大，但開挖到一定深度前，存在有限土體土壓力大于朗肯土壓力的情況。當(dāng)H=10，即b/H=0.6時(shí)，有限土體土壓力與朗肯土壓力接近相等。有限土體土壓力的增值隨著深度的增大而慢慢變小，當(dāng)開挖至基坑底部時(shí)，其土壓力約為朗肯土壓力的60 %。

分析水位高度h對(duì)有限土體土壓力的影響，改變上述算例中水位高度h，分別取1 m，2 m，3 m時(shí)，將得到有限土體的主動(dòng)土壓力作比較，結(jié)果如圖5所示，當(dāng)開挖深度小于擬開挖深度的一半時(shí)，水位高低對(duì)有限土體的土壓力影響很小，當(dāng)開挖深度大于擬開挖深度的一半時(shí)，隨著水位高度h的增大，有限土體土壓力逐漸變大。

圖4 本文方法土壓力與朗肯土壓力對(duì)比圖
Fig.4 Comparison diagram of the earth pressure ofthis article’s method and Rankine’s theory

圖5 不同水位對(duì)有限土體土壓力的影響
Fig.5 Influence for earth pressure oflimited soils on different water levels

3 數(shù)值模擬分析

3.1 計(jì)算模型及參數(shù)

圖6 基坑計(jì)算斷面圖
Fig.6 Cross section of the foundation pit

圖7 網(wǎng)格劃分示意圖
Fig.7 Meshing diagram

3.2 計(jì)算工況模擬

模擬基坑采取邊開挖邊降水的方式，開挖與支護(hù)總共分為7步進(jìn)行，如下表2所示，考慮到土體的初始面并非是完整的水平面，需采用重力加載的方式使其完成先期固結(jié)，河流水壓力作用的河床也需使其邊界條件固定。

表2 基坑開挖數(shù)值模擬計(jì)算工況Tab.2 Numerical simulation calculation of foundation pit excavation

3.3 計(jì)算結(jié)果分析

本文模擬了臨河距離為6 m，8 m，10 m三種情況下，按上述工序進(jìn)行基坑的開挖與支護(hù)，臨河部分有限土體作用在擋土墻上的主動(dòng)土壓力分布情況，計(jì)算結(jié)果如圖8所示，有限土體的土壓力隨著開挖深度的增加而增大，增幅逐漸減小，在與朗肯土壓力理論中土壓力隨深度線性增大存在著較大差異。隨著臨河距離(有限土體寬度)的增大，有限土體土壓力也逐漸增大，并趨近于朗肯土壓力值，符合實(shí)際情況。

為了驗(yàn)證本文所述計(jì)算方法的正確性，在一定水位條件下，將基坑開挖形成的臨河有限土體的土壓力理論計(jì)算結(jié)果與數(shù)值計(jì)算結(jié)果進(jìn)行對(duì)比分析。由圖9可以得到，理論解與模擬解呈非線性關(guān)系，與朗肯土壓力結(jié)果相比差異較大。理論解與數(shù)值解總體是小于朗肯土壓力結(jié)果，但在開挖深度為7～10 m，即0.6

圖8 不同臨河距離對(duì)土壓力的影響
Fig.8 Influence for earth pressure ondifferent distance to the river

圖9 計(jì)算結(jié)果對(duì)比
Fig.9 Comparison of calculation results

4 結(jié)論

基于目前對(duì)基坑支護(hù)有限土體的研究均針對(duì)緊鄰建筑物的情況，本文基于極限平衡理論，對(duì)開挖基坑臨近河流的情況進(jìn)行了有限土壓力計(jì)算分析。建立了臨河條件下有限寬度土體土壓力的計(jì)算模型，推導(dǎo)出有限土體主動(dòng)土壓力的計(jì)算公式，并探究其應(yīng)用條件與范圍。通過算例分析與有限元軟件數(shù)值模擬，與朗肯土壓力理論對(duì)比，得出以下結(jié)論：

①基于極限平衡理論，推導(dǎo)得出臨河狀態(tài)下有限土體的主動(dòng)土壓力計(jì)算公式，公式推導(dǎo)較為嚴(yán)密，具有一定的準(zhǔn)確性。

②有限土體的土壓力呈非線性分布，隨著基坑開挖，b/H逐漸減小，有限土體土壓力增量逐漸減小，開挖至坑底時(shí)，有限土體主動(dòng)土壓力比朗肯土壓力的計(jì)算值小得多。

③有限土體主動(dòng)土壓力總體小于朗肯土壓力計(jì)算值，但當(dāng)0.6

此外，以下情況也值得研究：①考慮河床的影響，即河岸并非垂直于地平面，滑動(dòng)土體并非楔形體，而是五邊形的情況。②河流深度小于擋土墻的高度的情況，即滑動(dòng)土體的破裂面被河流底部以下位置所截?cái)?，此類情況也應(yīng)著重考慮。