黃 睿，夏唐代

(浙江大學軟弱土與環(huán)境土工教育部重點實驗室，310058 杭州)

無粘性土的地震主動土壓力水平層法改進

黃 睿，夏唐代

(浙江大學軟弱土與環(huán)境土工教育部重點實驗室，310058 杭州)

為了研究主應力偏轉(zhuǎn)對土壓力的影響，基于擬靜力法，考慮平移模式下?lián)跬翂笸馏w的主應力偏轉(zhuǎn)和土層間的水平剪力，改進水平層分析法，根據(jù)應力狀態(tài)和靜力平衡條件建立基本方程組，通過迭代搜索程序進行求解，計算得出滑裂面形狀、地震主動土壓力分布、土壓力系數(shù)和合力作用點，并對其進行參數(shù)討論.結(jié)果表明:對于無粘性土，搜索得到的滑裂面曲線形狀近似直線，故簡化計算時，采用平面滑裂面的假設比較合理;填土內(nèi)摩擦角和水平地震作用對土壓力分布的影響較大，墻土摩擦角和豎向地震作用影響相對較小.與傳統(tǒng)擬靜力法對比，兩種方法所得的土壓力分布基本吻合，但是，改進水平層法所得的曲線呈現(xiàn)非線性變化趨勢.

地震;土壓力;水平層法;主應力偏轉(zhuǎn)

擋土墻的土壓力作為土力學的經(jīng)典問題，一直受到很多工程和科研人員的關(guān)注，其計算的關(guān)鍵問題在于確定作用在墻背上的側(cè)土壓力性質(zhì)、大小、方向和作用點，經(jīng)典的庫倫理論和朗肯理論已經(jīng)可以基本滿足工程需要，但是，地震作用時，擋土墻的土壓力計算卻是一個比較復雜的問題，相關(guān)的理論研究還不完善.目前，國內(nèi)外普遍采用的計算方法是基于庫倫理論演化出的Mononobe-Okabe擬靜力法，即用慣性力考慮地震作用，以靜力學的方法來近似處理動力學的問題，該方法思路清晰，計算簡便，被看作研究地震作用下土壓力的基本理論.文獻［1－3］提出采用水平層分析法計算土壓力和側(cè)壓力系數(shù)，并基于擬靜力法，對地震下的土壓力分布進行探討.文獻［4－5］在此基礎上利用水平層法對傾斜擋土墻背等復雜工況的地震土壓力進行了深入研究.文獻［6－8］根據(jù)水平土層的靜力和力矩平衡建立基本方程組，對擋土墻和邊坡的地震土壓力問題進行了求解.

然而，已有的地震土壓力研究很少考慮土拱效應引起的主應力偏轉(zhuǎn)影響，而當土體與結(jié)構(gòu)之間存在摩擦作用時，必定有土拱效應，進而會改變土中主應力的方向，影響土壓力的分布.文獻［9－11］通過預先假設土拱曲線的形狀，對土拱效應下平面滑裂面的主動、被動土壓力分布進行了一系列理論性分析.文獻［12－13］針對上述方法中土層單元無法滿足靜力平衡的不足，將土體的應力分析和水平層的靜力平衡條件相結(jié)合，提出考慮主應力偏轉(zhuǎn)的土壓力計算改進方法，并與經(jīng)典土壓力和試驗結(jié)果進行對比.

本文基于前人研究成果，以平移模式下的剛性擋土墻后的無粘性土為研究對象，考慮土拱效應對主應力偏轉(zhuǎn)的影響，事先并不假定滑裂面的形狀，而通過土體的應力分析和水平單元層的靜力平衡條件建立基本方程，得到每個水平單元層關(guān)于各層的主應力和滑裂面傾角的隱式方程組，由于在土拱的實際形成過程中，擋土墻和基底土體對滑動土體有約束作用，水平土層之間的剪應力必定客觀存在，故基本方程中加入水平剪力部分.通過編程，對多層水平單元的土壓力進行迭代搜索計算，求出滑裂面的形狀，土壓力的分布，土壓力系數(shù)和合力作用點，并討論了土體內(nèi)摩擦角、墻土外摩擦角、地震加速度影響系數(shù)等參數(shù)對計算結(jié)果的影響，再通過與傳統(tǒng)擬靜力法的計算結(jié)果對比，驗證本文方法的合理性.

1 理論模型

1.1 計算假定

本文土壓力理論模型見圖1，計算采用以下基本假定［12－13］:

1)擋土墻為剛性，墻背豎直，填土面水平;

2)墻后填土為各向同性且均質(zhì)的無粘性土;

3)當填土處于極限平衡狀態(tài)時，墻背處和滑裂面土體達到應力極限狀態(tài);

4)滑動土體有兩條滑裂面，一條沿著墻背，另一條位于墻后填土內(nèi)，滑裂面形狀由水平土層的靜力平衡方程確定;

5)每個水平土層單元上的應力沿土層長度的變化不大，可由每層對應的墻背處D和滑裂面處E的應力平均值來表示;

6)地震作用按照擬靜力法思路，以慣性力的形式施加在水平土體單元層上.

圖1 擋土墻地震主動土壓力計算模型

1.2 應力分析

在平移模式下，擋土墻的位移達到一定程度時，墻背面和滑裂面的土體達到極限應力狀態(tài)，由于摩擦力的存在，主應力的方向發(fā)生偏轉(zhuǎn)，φ為填土內(nèi)摩擦角，δ為墻土間外摩擦角，θ為該水平層的滑裂面傾角，墻背處D和滑裂面處E的應力摩爾圓見圖2.

圖2 墻背和滑裂面處土體的應力狀態(tài)

墻背面D點的應力狀態(tài):

1)水平應力

式中σ1為大主應力，Ka為朗肯主動土壓力系數(shù)，

ψD為墻背面處主應力偏轉(zhuǎn)角，

2)剪應力

3)豎向應力

4)墻背面反應力的水平和豎直分量

墻后填土中滑裂面E點的應力狀態(tài):

1)水平應力

式中ψE為填土中滑裂面處的主應力偏轉(zhuǎn)角，由文獻［13］分析可得

2)剪應力

3)豎向應力

4)滑裂面外法線方向的正應力

5)平行于滑裂面的切應力

6)滑裂面處反力的水平分量

7)滑裂面處反力的豎直分量

1.3 水平層受力分析

從滑動土體中，任取一水平單元薄層DE，見圖1，其上表面長Li，下表面長 Li+1，上表面深度zi，厚度為Δz，作用于水平土層上的力包括:上下表面的豎向應力pi、pi+1和剪力Vi、Vi+1;墻面處反力合力為Qi，其水平和豎向分量Qxi、Qzi;每層的墻面處反力應力的水平和豎向分量qxi、qzi;滑裂面處反力合力為Ri，其水平和豎向分量Rxi、Rzi;每層的滑裂面處反力應力的水平和豎向分量rxi、rzi;水平土層自重Wi，水平和豎向地震慣性力khWi、kvWi，其中，kh、kv分別為水平和豎向地震加速度影響系數(shù).

根據(jù)基本假設5)，水平單元土層上的豎向應力和剪應力都用滑動土層兩端土體的應力平均值表示，則上下表面的剪力表示為

平均豎向應力

墻背面反力的水平分量

墻背面反力的豎向分量

滑裂面反力合力的水平分量，用應力沿著水平層厚度的平均值表示

同理，滑裂面反力合力的豎向分量

水平單元土層的自重

其中γ為土重度，

1.4 基本方程和計算方法

水平單元土層DE在各力作用下靜力平衡.水平方向上有

豎直方向上有

地震主動土壓力水平合力可表示為

土壓力合力為

合力作用點可表示為

地震主動土壓力系數(shù)為

當滑動土體被劃分為n個水平單元層時，式(23)、(24)就組成了2n個基本方程，通過整理后，方程組未知量為每個水平單元土層的滑裂面傾角θ0～θn和大主應力σ10～σ1n，共有2n+2個未知量，需要另外補充兩個條件方能求解，對于填土的邊界面(z=0)處，設墻后填土有均布超載p，必須滿足第一個水平單元層的上表面豎向平均應力與超載相等，即p0=p;另外，在水平單元層較薄時，相鄰土層的滑裂面傾角變化很小，可近似取θ0= θ1.

由于基本方程組對θ是隱式方程，不能直接求解，本文采用迭代搜索算法編程計算，首先假定填土頂面的水平單元土層長L0和滑裂面傾角θ0(可采用庫倫土壓力理論所得的滑裂面角為初值)，并由邊界面超載p表示出第一水平層的σ10，將σ10代入式(24)表示出σ11，然后把θ0和σ11代入式(23)驗證是否滿足水平方向平衡條件，若不滿足則調(diào)整θ0，重復上述步驟直到滿足平衡.與第一層類似，每個水平土層采用此算法搜索滿足平衡條件的θi，并表示出Li，當逐層計算到滑動土體的底部時，若計算出的Ln≈0，即滿足預定的計算精度，則輸出計算結(jié)果，否則，設定步長調(diào)整L0，Ln＞0就減小L0，Ln＜0就增大L0，步長的設置需要反復地試算.重新搜索計算，直到滿足精度要求為止.當?shù)囁氵^程收斂時，計算結(jié)果一旦滿足要求，程序就結(jié)束計算，所得的結(jié)果在精度范圍內(nèi)是唯一的.最后，將計算結(jié)果按照式(25)～(28)可求出主動土壓力合力、合力作用點和土壓力系數(shù)，計算流程見圖3.

圖3 主動土壓力的計算流程

2 參數(shù)分析

現(xiàn)按照所述的計算思路，對參數(shù)進行分析，統(tǒng)一選取土重γ=18 kN/m3，墻高H=10 m，填土面超載p=0，討論土的內(nèi)摩擦角φ、墻土外摩擦角δ、水平和豎向地震加速度影響系數(shù)kh和kv，對滑裂面、土壓力分布、土壓力系數(shù)及合力作用點的影響.

2.1 滑裂面形狀

本文方法并沒有預先假設滑裂面的形狀，而是根據(jù)對每一水平土層的靜力平衡，搜索該層對應的傾角值θi，由此得到每個土層的長度Li，從而可繪出滑裂面形狀，不同參數(shù)所得的滑裂面形狀見圖 4，工況情況分別為:(a)φ =25°、30°、35°、40°，δ=0.5φ，kh=0.2，kv=0.5kh;(b)φ =40°，δ =0、1/3φ、1/2φ、2/3φ，kh=0.2，kv=0.5kh;(c)φ =40°，δ=0.5φ，kh=0、0.1、0.15、0.2、0.25，kv=0.5kh;(d)φ =40°，δ=0.5φ，kh=0.2，kv=0、1/3kh、1/2kh、2/3kh、kh.

由圖4可知，通過搜索繪出的滑裂面形狀曲線近似為直線狀，對于一般工程計算，地震作用下，假設擋土墻后非粘性土的滑裂面為平面是比較合理的.參數(shù)討論表明，土的內(nèi)摩擦角φ和水平地震加速度系數(shù)kh對滑裂面傾斜角度的影響非常明顯，而墻土摩擦角δ和豎向加速度系數(shù)kv的影響相對較小.滑裂面相對水平面的總體傾斜角度隨著地震加速度的增大而減小，隨著內(nèi)外摩擦角的增大而增大.土拱引起的主應力偏轉(zhuǎn)主要由參數(shù)δ體現(xiàn)，計算表明，δ對滑裂面形狀的影響主要體現(xiàn)在墻高的中上部，在墻高中部以下，滑裂面形狀差別不大.值得注意的是，基于本文假設所得的滑裂面形狀，在加速度作用較大或摩擦角較小時，曲線出現(xiàn)了略微反拱，其原因主要是根據(jù)靜力平衡條件搜索滑裂面傾角時，擋土墻的中上部土層大主應力較小，由大主應力表示的水平剪力也較小，水平向剪力和地震力對土層的影響不如擋土墻下部土層的大，靜力平衡時的滑裂面角就出現(xiàn)上面小下面大的情況.

圖4 滑裂面形狀曲線的參數(shù)分析

2.2 土壓力分布

按流程圖，選取與滑裂面分析時類似的工況，得到的土壓力分布見圖5.由于不同墻土摩擦角δ時，土壓力分布曲線的區(qū)分度較小，故圖5只給出不同φ、kh、kv時的土壓力分布曲線.結(jié)果同樣表明，土的內(nèi)摩擦角和水平地震加速度影響系數(shù)對土壓力分布的影響較為顯著，豎向加速度系數(shù)影響較小.土壓力隨著水平和豎向地震加速度影響系數(shù)的增大而增大，隨著內(nèi)摩擦角的增大而減小.計算深度在0.8H以內(nèi)時，土壓力的分布近似為線性變化，在0.8H到H范圍內(nèi)，與非地震作用時的土壓力先增大后減小變化不同，地震作用下，土壓力出現(xiàn)了增大趨勢陡增的非線性變化，并且水平地震作用越大，非線性化的陡增趨勢越明顯，可見，地震會引起土壓力的重分布.與文獻［10］相比，本文針對平動模式的擋土墻，采用不同的滑裂面假定，并且利用靜力平衡條件考慮土拱影響，再加之地震作用，故土壓力的分布形式存在較大不同.

圖5 土壓力分布曲線的參數(shù)分析

2.3 土壓力系數(shù)和合力作用點

按照式(27)、(28)，求得不同摩擦角和水平地震加速度系數(shù)下，地震主動土壓力系數(shù)Kae和合力作用點hae的結(jié)果見表1.相同地震加速度時，Kae隨內(nèi)摩擦角φ的增大和水平地震作用的減小而顯著減小，隨墻土摩擦角δ的增大而略微減小;與此不同，合力作用點hae對內(nèi)摩擦角的變化并不敏感，內(nèi)摩擦角相同時，δ較大，合力作用點略微向上升高;摩擦力情況相同時，作用點隨地震作用的增大而下降.

表1 土壓力系數(shù)和合力作用點

3 理論對比

采用本文方法和傳統(tǒng)的Mononobe－Okabe擬靜力法計算地震主動土壓力的分布，結(jié)果對比見圖 6，計算參數(shù)為 φ =40°，δ=0.5φ，kh=0.2，kv=0.5kh.

圖6 土壓力分布曲線對比

由圖6可看出，本文方法所得結(jié)果和擬靜力法整體上基本吻合，但由于傳統(tǒng)擬靜力法是由庫倫土壓力理論直接演化而來，土壓力分布曲線是由合力計算公式沿深度求導所得，而合力是關(guān)于深度的二次函數(shù)形式，故曲線為線性變化.本文則是將擬靜力的思路結(jié)合水平層法計算土壓力，同時考慮了土拱效應造成的主應力偏轉(zhuǎn)，所得結(jié)果在計算深度0.80H以內(nèi)與傳統(tǒng)擬靜力法差別不大，而在0.8H到H段，土壓力值非線性增大趨勢變強，整個土壓力分布為曲線形式.

4 結(jié)論

1)基于擬靜力法的思路，考慮土拱效應引起的主應力偏轉(zhuǎn)和土層間的剪應力，對平移模式下的擋土墻地震主動土壓力的水平層法進行了改進.通過土體的應力分析和水平層的靜力平衡條件建立基本方程組，采用編程的方法進行搜索計算，得到極限平衡狀態(tài)下的滑裂面形狀，進而求出土壓力分布曲線，土壓力系數(shù)和合力作用點.

2)對于非粘性土的情況，按照本文方法逐層搜索得到的水平層滑裂面傾角θi，沿深度變化的曲線近乎為直線，故對于工程中的簡化計算，假設平面滑裂面是比較合理的.

3)填土內(nèi)摩擦角和水平地震作用對土壓力分布有顯著影響，而墻土間摩擦角和豎向地震作用的影響相對較小;土壓力系數(shù)Kae隨著內(nèi)摩擦角增大和水平地震作用的減小而顯著減小;隨墻土摩擦角的增大而略微減小;合力作用點高度隨內(nèi)摩擦角變化不大，隨墻土摩擦角的增大和地震作用的減小而略微升高.

4)與傳統(tǒng)擬靜力法計算的地震主動土壓力分布相比，本文計算所得結(jié)果基本相符，但是土壓力曲線呈非線性變化.

［1］WANG Yuanzhan.Distribution of earth pressure on a retaining wall［J］.Geotechnique，2000，50(1):83－88.

［2］王元戰(zhàn)，李新國，陳楠楠.擋土墻主動土壓力分布與側(cè)壓力系數(shù)［J］.巖土力學，2005，26(7):1019－1022.

［3］王立強，王元戰(zhàn)，遲麗華.擋土墻地震土壓力及其分布［J］.中國港灣建設，2007(5):1－5.

［4］楊劍，高玉峰，程永鋒，等.地震條件下傾斜擋土墻被動土壓力研究［J］.巖土工程學報，2009，31(9):1391－1397.

［5］林宇亮，楊果林，趙煉恒，等.地震動土壓力水平層分析法［J］.巖石力學與工程學報，2010，29(12):2581－2591.

［6］NOURI H，JONES C.Development of horizontal slice method for seismic stability analysis of reinforced slopes and walls［J］.Geotextile and Geomembranes，2006，24(3):175－187.

［7］AZAD A，YASROBI S，PAK A.Seismic active pressure distribution history behind rigid retaining walls［J］.Soil Dynamics and Earthquake Engineering，2008，28(5):365－375.

［8］NOURI H，F(xiàn)AKHER A，JONES C.Evaluating the effects of the magnitude and amplification of pseudostatic acceleration on reinforced soil slopes and walls using the limit equilibrium horizontal slices method［J］.Geotextile and Geomembranes，2008，26(3):263－278.

［9］PAIK K H.Estimation of active earth pressure against rigid retaining walls considering arching effects［J］.Geotechnique，2003，53(7):643－653.

［10］應宏偉，蔣波，謝康和.考慮土拱效應的擋土墻主動土壓力分布［J］.巖土工程學報，2007，29(5):717－722.

［11］吳明，彭建兵，徐平，等.考慮土拱效應的擋墻后土壓力研究［J］.工程力學，2011，28(11):89－95.

［12］章瑞文.擋土墻主動土壓力理論研究［D］.杭州:浙江大學，2007.

［13］章瑞文，徐日慶.土拱效應原理求解擋土墻土壓力方法的改進 ［J］.巖土力學，2008，29(4):1057－1060.

Improved horizontal slice method for seismic active earth pressure of cohesionless soil

HUANG Rui，XIA Tangdai
(MOE Key Laboratory of Soft Soils and Geoenvironmental Engineering，Zhejiang University，310058 Hangzhou，China)

Based on pseudo-static theory，improved horizontal slice method is derived for retaining walls in translational mode to investigate the influence of principle stress deflection on the earth pressure.The method considers the principle stress deflection of soil behind walls and the horizontal shear force between slices.The basic equations are established according to stress state and static equilibrium condition，and solved by using iterative search algorithm to calculate the slip surface shape，seismic active earth pressure distribution，earth pressure coefficient and resultant force point，then parametric discussion is taken.Numerical results indicate that the slip surface shape for cohesionless soil is approximately linear，so it’s quite reasonable to use the hypothesis of planar slip surface in simplified calculation.In addition，soil friction angle and horizontal seismic action significantly impact earth pressure distribution，instead，the influence of wall friction angle and vertical seismic action are relatively smaller.Compared with the traditional pseudo-static method，earth pressure distributions by two different methods basically agree with each other，but the curve from the improved horizontal slice method trends towards nonlinear increase.

earthquake;earth pressure;horizontal slice method;principle stress deflection

TU432

A

0367－6234(2014)02－0085－07

2013－06－24.

國家自然科學基金與高速鐵路基礎研究聯(lián)合基金資助項目(U1234204).

黃 睿(1987—)，男，博士研究生;

夏唐代(1965—)，男，教授，博士生導師.

黃 睿，zjuhr@hotmail.com.

(編輯 趙麗瑩)