王 博 呂 果 李 江

（建設(shè)綜合勘察研究設(shè)計院有限公司，北京 100007）

0 引言

砂土宏觀上具有離散性和非均勻性，微觀上具有多樣性和易變性[1]，在路基工程、山區(qū)邊坡防護等工程領(lǐng)域中較為常見，保證其穩(wěn)定性尤為重要，因此分析砂土的抗剪能力極為關(guān)鍵。

目前，關(guān)于砂土剪切特性研究主要從實際試驗和模擬試驗兩方面開展，對于實際試驗而言，Been 等[2]通過采用狀態(tài)參數(shù)來描述Kogyuk 砂的剪脹性。郝冬雪等[3]進行砂土三軸剪切試驗，發(fā)現(xiàn)高壓下砂土應力-應變曲線后期無應變軟化現(xiàn)象。對于模擬試驗而言，劉方成等[4]在砂土中摻入橡膠，進行了直剪試驗離散元模擬，從力鏈和顆粒位移等分析了細觀機制的變化；李 航等[5]從顆粒形狀和孔隙率出發(fā)，模擬了標準砂剪切；楊忠平等[6]、張振平等[7]利用顆粒離散元軟件進行了土石混合體直剪試驗模擬，均得出了混合體抗剪強度隨含石量增加而增加的結(jié)論。

盡管眾多學者在砂土剪切特性方面取得了一定研究成果，但所采用砂土均具有單一性，沒有考慮砂土粒徑對其力學特性的影響。河砂作為一種在工程領(lǐng)域常見的土體，通常粒徑變化范圍較大，研究其不同粒徑下的剪切特性對于砂土穩(wěn)定性研究尤為關(guān)鍵。為此，利用應變控制式直剪儀進行3 種不同平均粒徑砂（分別對應細砂、中砂和粗砂）的直剪試驗，分析了剪切應力-剪切位移關(guān)系和界面抗剪強度指標等宏觀力學變化，建立了砂土抗剪強度預測模型。此外，開展離散元PFC2D 進行砂土直剪試驗，對不同平均粒徑下的剪切帶厚度變化規(guī)律進行分析，針對通過改變砂土粒徑提升砂土抗剪能力進行研究。

1 試驗設(shè)備、材料、方案

1.1 試驗設(shè)備及材料

試驗設(shè)備為ZJ 型直剪儀（見圖1），上下剪切盒通過插銷構(gòu)成整體，試驗準備前先放置透水石，再放置砂土，最后放置透水石。砂土試樣直徑6.18 cm、高度4 cm。試驗所用砂土為某工地河砂，將砂土烘干后分為3 個粒組，分別為0.075～0.25 mm（細砂）、0.25～0.5 mm（中砂）和0.5～2 mm（粗砂），如圖2 所示。砂土的土粒比重Gs為2.64，三種砂土取平均粒徑dⅠ、dⅡ和dⅢ分別為0.15 mm、0.36 mm 和1.22 mm，初始密度ρ分別為1.37 g/cm3、1.40 g/cm3和1.45 g/cm3，剪切試樣直徑為61.8 mm、高為20 mm。

圖1 應變控制式直剪儀示意圖

圖2 試驗用砂

1.2 試驗方案

試驗剪切速率為0.8 mm/min，考慮因素為3 種不同的平均粒徑，4 種法向應力（50 kPa、100 kPa、150 kPa 和200 kPa），共12 組試驗。

2 宏觀試驗結(jié)果分析

2.1 剪切應力-剪切位移特征曲線

剪切應力-剪切位移特征曲線見圖3，通過對圖3 分析發(fā)現(xiàn)，剪切應力隨剪切位移先快速增大后緩慢增大，到達抗剪強度后基本穩(wěn)定。匯總出各情況下抗剪強度（見表1），由表1 可見，在同一平均粒徑下，法向應力越大，抗剪強度越大，但增大幅度基本一致。在同一法向應力下，砂的抗剪強度隨砂土平均粒徑的增加而增加。

表1 各情況下抗剪強度表

圖3 剪切應力-剪切位移特征曲線

2.2 抗剪強度的預測

以法向應力σ和平均粒徑d作為自變量，以抗剪強度τf作為因變量，進行多元線性回歸分析，可以進行抗剪強度的預測，公式如下：

式中：β0、β1、β2為回歸系數(shù)。將2.1 節(jié)中數(shù)據(jù)代入式（1），回歸后得到多元線性回歸關(guān)系為：

將試驗抗剪強度與預測抗剪強度繪制到圖4中，發(fā)現(xiàn)預測與試驗走勢相同，數(shù)據(jù)基本一致，說明該多元線性模型回歸效果良好，利用法向應力σ和平均粒徑d來測算抗剪強度τf較為準確，可以應用于實際工程。

圖4 模型與試驗數(shù)據(jù)對比

2.3 抗剪強度指標分析

在巖土領(lǐng)域，有多種土體破壞準則，例如：Tresca破壞準則、Drucke-Mises 破壞準則、Mohr-Coulomb 破壞準則和雙剪破壞準則等[8]。對于干砂而言，黏聚力c為0，因此適合采用Coulomb 破壞準則對抗剪強度τf與法向應力σ間的關(guān)系進行表達，具體表達式如下：

式中：φ為內(nèi)摩擦角。

圖5 為抗剪強度-法向應力關(guān)系，可以看出：相關(guān)系數(shù)R2均在0.999 以上，Coulomb 破壞準則表達準確。同樣發(fā)現(xiàn)抗剪強度隨法向應力和平均粒徑的增加而增大。通過式（1）得出平均粒徑dⅠ、dⅡ和dⅢ對應的內(nèi)摩擦角分別為27.88°、29.55°和32.46°，可見內(nèi)摩擦角隨平均粒徑的增大而增大。

圖5 抗剪強度-法向應力關(guān)系

3 剪切帶厚度規(guī)律分析

砂土在剪切結(jié)束后產(chǎn)生剪切帶，剪切帶厚度是巖土界關(guān)注的重點，而常規(guī)的室內(nèi)試驗無法觀測到剪切帶，也無法確定其厚度，只能借助離散元等微觀軟件進行判斷。離散元法作為描述顆粒土力學特性的數(shù)值分析手段被廣泛應用于巖土工程領(lǐng)域[9-10]，常用的離散元模擬軟件為PFC2D 或3D 模擬軟件，本文采用PFC2D 進行不同平均粒徑下剪切帶厚度分析。

3.1 模型建立及試驗方案

3.1.1 模型假設(shè)與尺寸建立

模型假設(shè)顆粒均為剛性體，且為圓形；顆粒與顆粒間通過點進行連接接觸，且受力可以發(fā)生平移、轉(zhuǎn)動。設(shè)置①-⑥號廣義墻，其內(nèi)部放置土體顆粒，如圖6 所示。對于①-③廣義墻模擬上剪切盒、④-⑥廣義墻模擬下剪切盒，為了與室內(nèi)常規(guī)直剪儀盛放土體的尺寸相同，剪切盒整體尺寸為6.18 cm（長）×4 cm（高）。設(shè)定④-⑥固定不動，為了消除慣性力對剪切試驗的影響，通過控制應變率使①-③墻向右運動緩慢均勻運動。為與室內(nèi)試驗結(jié)合，平均粒徑分為dⅠ、dⅡ和dⅢ三種，即為0.15 mm、0.36 mm 和1.22 mm，但采用改顆粒尺寸計算時間長、效率慢，因此采用粒徑放大法[11]，將平均粒徑dⅠ、dⅡ和dⅢ尺寸擴大3倍，即為0.45 mm、1.08 mm 和3.66 mm。

圖6 土體直剪試驗PFC2D 離散元模型

3.1.2 模型參數(shù)設(shè)置

本模擬試驗中接觸問題均采用線性剛度模型，公式為：

式中：kn、ks分別為法向、切向接觸剛度，k[A]、k[B]表示任意兩個接觸實體單元的顆粒剛度，對于顆粒間和顆粒-墻體間的法向接觸剛度，取值分別為1×108N/m和4×108N/m，同樣的，切向接觸剛度取值分別為1×108N/m 和4×108N/m。

3.1.3 模擬方案

為與室內(nèi)試驗對比，平均粒徑分為dⅠ、dⅡ和dⅢ三種，法向應力分50 kPa、100 kPa、150 kPa 和200 kPa 四種，剪切至位移4 mm 結(jié)束試驗，剪切速率基本與室內(nèi)試驗一致。由于室內(nèi)宏觀試驗與模擬試驗土體摩擦系數(shù)有較大差別，因此結(jié)合室內(nèi)摩擦系數(shù)進行參數(shù)調(diào)整，即為調(diào)參過程。最終確定3 種平均粒徑下的土體摩擦系數(shù)分別為0.22、0.25 和0.30。

3.2 細觀結(jié)果分析

3.2.1 模型合理性驗證

由于篇幅有限，因此僅展示部分試驗與模型的剪切應力-剪切位移關(guān)系（見圖7）。通過對圖7 分析發(fā)現(xiàn)，通過PFC2D 模擬得出的數(shù)據(jù)與試驗數(shù)據(jù)差距不大，兩者曲線走勢基本一致，這充分說明模型建立合理，可以用于剪切帶微觀分析。

圖7 剪切應力-剪切位移試驗與模型對比

3.2.2 不同平均粒徑下剪切帶規(guī)律分析

剪切后期，剪切應力處于穩(wěn)定狀態(tài)，此時可進行剪切帶厚度的確定。為了便于剪切帶厚度分析，因此對顆粒進行染色處理，處理為10 條綠色和10 條紅色，以平均粒徑Ⅰ進行說明，見圖8（a），剪切試驗初始狀態(tài)，10 條綠色和10 條紅色均為垂直分布，當剪切完畢后，在剪切界面處顆粒發(fā)生錯動，以錯動的垂直距離定義為剪切帶厚度L（見圖8（b））。

圖8 剪切帶厚度觀條設(shè)置示意圖

圖9 為不同平均粒徑下剪切帶厚度變化，發(fā)現(xiàn)砂土平均粒徑越大，剪切帶厚度越大，其中平均粒徑dⅠ到dⅡ的增加幅度小于dⅡ到dⅢ。法向應力越大，剪切帶厚度同樣越大，但增加幅度越來越小。此外，發(fā)現(xiàn)剪切帶厚度L與法向應力σ存在一元二次函數(shù)關(guān)系：

圖9 不同平均粒徑下剪切帶厚度變化

式中：a、b和c為擬合參數(shù)，匯總函數(shù)見表2，發(fā)現(xiàn)相關(guān)系數(shù)R2均在0.99 以上，擬合關(guān)系良好。

4 結(jié)論

利用應變控制式直剪儀，進行了3 種不同平均粒徑砂的室內(nèi)直剪試驗和離散元模擬試驗，得出以下結(jié)論：

（1）砂的抗剪強度隨砂土平均粒徑的增加而增加。以法向應力σ和平均粒徑d為自變量，建立了砂土抗剪強度τf多元線性預測模型，且模型預測準確，對工程有一定實際意義。此外，內(nèi)摩擦角隨平均粒徑的增大而增大。

（2）剪切帶厚度隨平均粒徑的增大而增加，其中平均粒徑dⅠ到dⅡ的增加幅度小于dⅡ到dⅢ。剪切帶厚度與法向應力之間符合一元二次函數(shù)關(guān)系。

綜上兩點結(jié)論，在實際工程施工中，合理采用粒徑較大的砂土有利于提高抗剪能力，增加工程穩(wěn)定性。