陳 芳 ，王新剛 ，吳永風(fēng)

（1.江西省水利科學(xué)研究院，江西 南昌 330029；2.江西省水工安全工程技術(shù)研究中心，江西 南昌 330029；3.中國地質(zhì)大學(xué) （武漢）工程學(xué)院，湖北 武漢 430074）

0 引 言

碎石土在自然界分布十分廣泛[1]，其邊坡的穩(wěn)定性受到工程界高度重視。本文室內(nèi)試驗所取土邊坡碎石土含量多，下滑力大，邊坡體內(nèi)碎石土所受剪應(yīng)力大，準(zhǔn)確確定碎石土本構(gòu)模型對建立邊坡數(shù)值模型、分析邊坡穩(wěn)定性有著重要的意義。

本構(gòu)模型是土力學(xué)的核心問題[2]。目前描述土的應(yīng)力應(yīng)變關(guān)系的數(shù)學(xué)模型有許多種，歸納起來有兩大類:一類是彈性模型，較典型是Duncan-Chang模型和K-G模型[3]；另一類是彈塑性模型，較典型的有Cambridge模型、清華模型、沈珠江的雙屈服面模型等。碎石土的力學(xué)性質(zhì)和本構(gòu)模型是巖土工程研究的一個重要課題，建立適用于碎石土的本構(gòu)模型一直是工程界、學(xué)術(shù)界努力的方向，Duncan提出的E-B模型[3]能很好地反映碎石土在破壞前的應(yīng)力應(yīng)變關(guān)系，但不能考慮碎石土的剪縮性，計算所得的水平位移偏大；沈珠江雙屈服面模型由于其模型簡單[4-6]、參數(shù)易于確定的特點，在工程中得到了較為廣泛的應(yīng)用，但在高圍壓或較大剪應(yīng)變條件下仍存在模擬土的剪脹性偏大的不足。筆者根據(jù)文獻(xiàn)[2]提出的結(jié)合兩種模型優(yōu)點的改進(jìn)模型，利用VC++編程在FLAC3D軟件中實現(xiàn)了改進(jìn)的鄧肯-張E-B模型的開發(fā)。

1 改進(jìn)模型在FLAC3D軟件中的開發(fā)

本文以 FLAC中的摩爾-庫倫本構(gòu)模型為框架，進(jìn)行改進(jìn)鄧肯-張本構(gòu)模型的二次開發(fā)，改進(jìn)的鄧肯-張本構(gòu)模型能較好地描述碎石土的非線性、反映土的剪脹性，且不受土體剪應(yīng)變較大的影響[2]。FLAC3D軟件所提供的udm.zip文件中包含摩爾-庫倫模型的頭文件 （.h）和C++源文件（.cpp）。改進(jìn)的鄧肯-張模型主程序開發(fā)流程見圖1。

圖1 改進(jìn)鄧肯-張模型主程序開發(fā)流程

2 工程應(yīng)用

2.1 模型

模型邊坡海拔4426～5341 m，相對高差 910 m，開采后人工邊坡加自然斜坡高1065 m，屬于高陡邊坡。北東部剖面線走向約北偏東30°，近乎垂直臺階坡面。為便于分析該邊坡開采后的穩(wěn)定性問題，采用ANSYS軟件前處理，建立邊坡北東部邊坡剖面計算模型 （見圖2）。共劃分單元88758個、節(jié)點78909個，在FLAC3D軟件中進(jìn)行計算，其中碎石土層采用本文開發(fā)的改進(jìn)鄧肯-張模型，其他巖土體采用摩爾-庫倫模型。

圖2 北東部剖面計算模型

2.2 模型參數(shù)

該邊坡的穩(wěn)定性分析與數(shù)值模擬計算所必需的物理力學(xué)參數(shù)由室內(nèi)土力學(xué)試驗獲取，具體數(shù)值見表1。

表1 數(shù)值計算模型參數(shù)

2.3 碎石土邊坡開挖后分析

北東部邊坡開挖位移分析見圖3。從圖3可知，采用3種不同模型的邊坡最大位移均出現(xiàn)在邊坡中部附近，高程約在4980～5000 m之間，且在此高程區(qū)間，邊坡開挖不穩(wěn)定，需要進(jìn)行護(hù)坡加固。采用摩爾-庫倫模型邊坡最大位移約為4.17 m，方向指向坡體外，計算結(jié)果偏安全，不適合土石混合材料。采用鄧肯-張模型邊坡最大位移約為5.41 m，由于該模型不能考慮碎石土的剪縮性，計算得出的水平位移偏大，計算結(jié)果偏危險。采用本文開發(fā)的改進(jìn)鄧肯-張模型邊坡最大位移約為4.69 m，由于該模型既考慮碎石土的體縮特性，又考慮了碎石土的剪脹性，適合于邊坡高剪應(yīng)力情況，避免邊坡在剪應(yīng)變較大時，剪脹體變偏大，計算結(jié)果較合理。

圖3 北東部邊坡開挖位移分析（單位:m）

為進(jìn)一步對比3種本構(gòu)模型剪應(yīng)力、水平位移的計算結(jié)果，在該邊坡計算模型碎石土區(qū)域及其附近關(guān)鍵位置設(shè)置14個監(jiān)測點，監(jiān)測點位置見圖4。

由監(jiān)測點剪應(yīng)力 （見圖5）可以看出:①采用摩爾-庫倫模型計算結(jié)果在1～7號監(jiān)測點剪應(yīng)力較大，7號監(jiān)測點以后剪應(yīng)力較??；②采用鄧肯-張模型在1～7號監(jiān)測點剪應(yīng)力較小，7號監(jiān)測點以后剪應(yīng)力較大；③采用本文改進(jìn)的鄧肯-張模型在各監(jiān)測點剪應(yīng)力均較為合理；④3種不同的模型計算得出的剪應(yīng)力在8號監(jiān)測點處達(dá)到最大，約為600～650 kPa。

圖4 1～14號監(jiān)測點位置示意

由監(jiān)測點水平位移 （見圖6）可以看出:①采用摩爾-庫倫模型計算結(jié)果在1～4號監(jiān)測點位移較大，在7號監(jiān)測點之后位移約1 m，整體位移較采用其他2種模型小，計算結(jié)果偏于安全；②采用鄧肯-張模型計算結(jié)果整體偏大，在7號監(jiān)測點之后位移約1.3 m；③采用本文改進(jìn)的鄧肯-張模型在7號監(jiān)測點之后位移約1.2 m，計算結(jié)果比較合理，既考慮了土的剪脹性，又適宜剪應(yīng)力較大的情況。

圖5 監(jiān)測點剪應(yīng)力

圖6 監(jiān)測點水平位移

3 結(jié) 論

（1）在Visual Studio環(huán)境中以FLAC中的摩爾-庫倫本構(gòu)模型為框架，進(jìn)行改進(jìn)鄧肯-張本構(gòu)模型的二次開發(fā)流程，并將VC++語言開發(fā)的本構(gòu)模型編譯成duncan-changshen.dll文件 （動態(tài)鏈接庫），在調(diào)用時加載上去。改進(jìn)的鄧肯-張本構(gòu)模型能較好地描述碎石土的非線性，反映土的剪脹性，且不受土體剪應(yīng)變較大的影響。

（2）以某碎石土邊坡為例，分別調(diào)用摩爾-庫倫模型、鄧肯-張模型、二次開發(fā)的改進(jìn)鄧肯-張模型，數(shù)值模型的計算結(jié)果進(jìn)行對比分析表明，改進(jìn)的鄧肯-張模型計算結(jié)果較合理，具有一定的應(yīng)用前景。

［1］王光進(jìn)，楊春和，張 超，等.粗粒含量對散體巖土顆粒破碎及強(qiáng)度特性試驗研究［J］.巖土力學(xué)，2009，30（12）:3649-3654.

［2］羅剛，張建民.鄧肯-張模型和沈珠江雙曲面模型的改進(jìn)［J］.巖土力學(xué)，2004，25（6）:887-890.

［3］DUNCAN J M，CHANG C Y.Nonlinear analysis of stress and strain in soils［J］.Journal of the Soil Mechanics and Foundations Division,ASCE,1970,96（5）:1629-1653.

［4］沈珠江.考慮剪脹性的土和石料的非線性應(yīng)力應(yīng)變模式［J］.水利水運科學(xué)研究，1986（4）:1-14.

［5］沈珠江，王劍平.橫山水庫土石壩有效應(yīng)力應(yīng)變分析［J］.水利學(xué)報，1990（4）:59-65.

［6］沈珠江.結(jié)構(gòu)性粘土的彈塑性損傷力學(xué)模型［J］.巖土工程學(xué)報，1993，15（3）:21-28.