孫　超　史迪菲　原利明

(吉林建筑大學(xué)測(cè)繪與勘查工程學(xué)院，長(zhǎng)春　130118)

擋土墻主動(dòng)土壓力理論與有限元理論分析比較

孫超史迪菲原利明

(吉林建筑大學(xué)測(cè)繪與勘查工程學(xué)院，長(zhǎng)春130118)

摘要：本文根據(jù)朗金主動(dòng)土壓力理論，通過研究彈性半空間體內(nèi)的應(yīng)力狀態(tài)，根據(jù)土的極限平衡條件得出的土壓力計(jì)算方法.利用大型通用軟件ABAQUS計(jì)算擋土墻主動(dòng)土壓力，并與朗金土壓力理論下的主動(dòng)土壓力進(jìn)行比較，說明不同模型選取對(duì)試驗(yàn)結(jié)果的影響.

關(guān)鍵詞：彈塑性；擋土墻；朗金主動(dòng)土壓力理論

擋土墻作為防止土體坍塌的構(gòu)筑物廣泛應(yīng)用于房屋建筑、水利、鐵路以及公路和橋梁工程中.由于擋土墻的土壓力計(jì)算十分復(fù)雜，目前計(jì)算多采用古典的朗金(Rankine，1857年)和庫倫(Coulomb，1773年)土壓力理論[1].隨著當(dāng)前有限元類軟件發(fā)展日新月異，大型通用軟件ABAQUS基于彈塑性理論對(duì)主動(dòng)土壓力理論的研究迅速發(fā)展，許多理論模擬可以通過計(jì)算得到，從而達(dá)到節(jié)省人力、物力的目的.因此，利用ABAQUS土壓力擁有廣泛地發(fā)展前景.

1模型及相關(guān)參數(shù)

1.1　Mohr-Coulomb模型與界面接觸模型[2]

土的Mohr-Coulomb破壞條件為：

(1)

式中：p與q分別為平均應(yīng)力和廣義偏差應(yīng)力；θσ為應(yīng)力Lode角；c與φ分別為土的粘聚力和內(nèi)摩擦角.

將總應(yīng)變速率分解為彈性應(yīng)變速率{εe}和塑性速率應(yīng)變{εp}兩部分：

(2)

其中，彈性應(yīng)變速率與塑性應(yīng)變速率分別由廣義Hook定律和塑性流動(dòng)法則確定，即

(3)

(4)

綜合考慮上列各式，采用塑性一致性條件，得到下列彈塑性本構(gòu)方程：

(5)

式中：[Dep]為彈塑性系數(shù)矩陣，在塑性加載條件下

(6)

而在卸載與中性變載條件下，[Dep]=[De].

擋土墻與土體在剛度上存在較大差異，其交界面會(huì)產(chǎn)生相對(duì)滑動(dòng).因此，在有限元數(shù)值分析中，擋土墻和土體之間設(shè)置接觸類型為表面與表面接觸，接觸屬性為切向行為，摩擦公式定義為罰函數(shù)算法，摩擦系數(shù)μ=0.001.

1.2　朗金土壓力模型

朗金土壓力理論通過研究彈性半空間體內(nèi)的應(yīng)力狀態(tài)，根據(jù)土的極限平衡條件而得出的土壓力計(jì)算方法，適用于擋土墻墻背豎直光滑，填土面水平.根據(jù)朗金土壓力理論，當(dāng)土體中某點(diǎn)處于極限平衡狀態(tài)時(shí)，大小主應(yīng)力σ1和σ3應(yīng)滿足以下關(guān)系式：

粘性土：

(7)

(8)

設(shè)墻背豎直光滑，填土面水平，當(dāng)擋土墻偏離土體時(shí)，墻背土體中離地表任意深度z處豎向應(yīng)力σz為大主應(yīng)力σ1，σx為小主應(yīng)力σ3，故可得朗金主動(dòng)土壓力強(qiáng)度σa為：

粘性土：

(9)

1.3　有限元計(jì)算模型

有限元方法的基礎(chǔ)是變分原理和加權(quán)余量法[3]，其基本求解思想是把計(jì)算域劃分為有限個(gè)互不重疊的單元，在每個(gè)單元內(nèi)，選擇一些合適的節(jié)點(diǎn)作為求解函數(shù)的插值點(diǎn)，將微分方程中的變量改寫成由各變量或其導(dǎo)數(shù)的節(jié)點(diǎn)值與所選用的插值函數(shù)組成的線性表達(dá)式，借助于變分原理或加權(quán)余量法，將微分方程離散求解.采用不同的權(quán)函數(shù)和插值函數(shù)形式，便構(gòu)成不同的有限元方法.

對(duì)于擋土墻和土體，采用8節(jié)點(diǎn)平面應(yīng)變單元C3D8R;擋土墻墻高7m，頂部寬0.6m，底部寬1m，墻背豎直，光滑，墻后填土面水平，并作用有均布荷載q=200kPa，擋土墻示意圖見圖1，各土層物理力學(xué)指標(biāo)見表1.

圖1 擋土墻示意圖

表1 各土層物理力學(xué)指標(biāo)

2朗金主動(dòng)土壓力理論，ABAQUS模擬彈性、彈塑性主動(dòng)土壓力強(qiáng)度結(jié)果及分析

2.1　朗金主動(dòng)土壓力強(qiáng)度計(jì)算結(jié)果

通過計(jì)算，雜填土頂面處的土壓力強(qiáng)度為81.20kPa，雜填土底面處的土壓力強(qiáng)度為107.66kPa，粉質(zhì)粘土頂面處的土壓力強(qiáng)度為91.57kPa，粉質(zhì)粘土底面處的土壓力強(qiáng)度為119.96kPa.

2.2　ABAQUS彈性模型、彈塑性模型[4]主動(dòng)土壓力強(qiáng)度計(jì)算結(jié)果

通過模型模擬結(jié)果，彈性理論下雜填土頂面處的土壓力強(qiáng)度為17.24kPa，粉質(zhì)粘土底面處的土壓力強(qiáng)度為148.79kPa，見圖2.彈塑性理論下雜填土頂面處的土壓力強(qiáng)度為17.84kPa，粉質(zhì)粘土底面處的土壓力強(qiáng)度為164.35kPa，見圖3.

圖2 彈性主動(dòng)土壓力應(yīng)力云圖　　　　　圖3 彈塑性主動(dòng)土壓力應(yīng)力云圖　　　　　圖4 壓應(yīng)力隨深度變化曲線

2.3　朗金主動(dòng)土壓力理論、ABAQUS彈性模型、彈塑性模型主動(dòng)土壓力強(qiáng)度分析比較

通過應(yīng)力云圖可以得出結(jié)論，土在彈塑性狀態(tài)下的壓應(yīng)力與土在彈性狀態(tài)下的壓應(yīng)力相比要小，擋土墻在主動(dòng)土壓力作用下的應(yīng)力云圖與理論基本一致，壓應(yīng)力最大值發(fā)生在擋土墻外表面與土接觸的部位.通過受力分析可知，這是由于主動(dòng)土壓力作用在擋土墻下表面上部三分點(diǎn)處，對(duì)擋土墻底產(chǎn)生力臂受壓所致.

通過圖4比較各應(yīng)力曲線可以看到，土在彈塑性與彈性兩種不同狀態(tài)下的對(duì)擋土墻壓應(yīng)力變化曲線形狀基本一致.在埋深較小時(shí)，與朗金理論不同的是地表壓應(yīng)力值，隨著深度的加大，三條曲線開始靠近，在土層變化處開始出現(xiàn)差別，土在彈性狀態(tài)下的壓應(yīng)力增量逐漸加大，彈塑性狀態(tài)下的壓應(yīng)力增量開始減小，但基本上呈線性增長(zhǎng)，與朗金主動(dòng)土壓力理論一致.由此可知，對(duì)土體的研究應(yīng)采用彈塑性模型，這樣所得的結(jié)果偏保守，比較安全.

3結(jié)論

通過建立彈性模型和彈塑性模型，在相同的條件下模擬，可以得出如下結(jié)論：

(1) 朗金理論計(jì)算得到的主動(dòng)土壓力偏保守，彈塑性理論和彈性理論計(jì)算得到的主動(dòng)土壓力有所增加；

(2) 彈性理論模型和彈塑性理論模型的壓應(yīng)力基本一致；

(3) 彈性模型下的主動(dòng)土壓力應(yīng)力值大于彈塑性模型下的主動(dòng)土壓力應(yīng)力值；

(4) 彈塑性理論的計(jì)算結(jié)果更接近于工程實(shí)際.

參考文獻(xiàn)

[1] 趙明華,俞曉.土力學(xué)與基礎(chǔ)工程[M].武漢:武漢理工大學(xué)出版社,2009:124-129.

[2] 費(fèi)康,張建偉.ABAQUS在巖土工程中的應(yīng)用[M].北京:中國水利水電出版社,2010:67-70.

[3] 崔春義,李順群,奕茂田,趙穎華.基于ABAQUS的連續(xù)墻支護(hù)深基坑開挖彈塑性數(shù)值分析[J].巖土工程學(xué)報(bào),2008(S1):86-90.

[4] 石亦平,周玉蓉.ABAQUS有限元分析實(shí)例詳解[M].北京:機(jī)械工業(yè)出版社,2006:86-98.

Retaining Wall Active Earth Pressure Theory and the Theory of Finite Element Analysis

SUN Chao，SHI Di-fei，YUAN Li-ming

(SchoolofGeometricsandProspectingEngineering,JilinJianzhuUniversity,Changchun,China130118)

Abstract：Rankine active earth pressure theory through the study of stress state,the elastic half space body root soil conditions of limit equilibrium of soil pressure calculation method.The common software ABAQUS is used to calculate active earth pressure of retaining,with Rankine earth pressure theory of active earth pressure under comparison,to illustrates the influence of different selection model on the test results.

Keywords:elastic-plastic；retaining wall；rankine active earth pressure theory

中圖分類號(hào)：TU 413

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：2095-8919(2015)06-0019-03

作者簡(jiǎn)介：孫超(1978~)，男，黑龍江省東寧縣人，副教授，博士.

收稿日期：2015-01-14.