陳 鋮劉小清羅正東寧志強(qiáng)

（①湘潭大學(xué)土木工程與力學(xué)學(xué)院 湘潭 411105）

（②巖土力學(xué)與工程國家重點(diǎn)實(shí)驗(yàn)室，中國科學(xué)院武漢巖土力學(xué)研究所 武漢 430071）

顆粒級配對粗粒土強(qiáng)度與變形特性影響的研究*

陳 鋮①②劉小清①羅正東①寧志強(qiáng)①

（①湘潭大學(xué)土木工程與力學(xué)學(xué)院 湘潭 411105）

（②巖土力學(xué)與工程國家重點(diǎn)實(shí)驗(yàn)室，中國科學(xué)院武漢巖土力學(xué)研究所 武漢 430071）

采用三維顆粒流軟件PFC3D的內(nèi)置FISH語言，進(jìn)行二次開發(fā)，模擬粗粒土級配，建立一個研究高鐵填方路基粗粒土變形特性的顆粒流模型。通過粗粒土三軸試驗(yàn)，標(biāo)定了表征粗粒土細(xì)觀力學(xué)性質(zhì)的模型參數(shù)，并驗(yàn)證了模型的有效性。考慮粗粒土曲率系數(shù)和不均勻系數(shù)對顆粒級配的影響，模擬并研究了5種不同級配下粗粒土在3種不同圍壓下的變形特性。試驗(yàn)結(jié)果表明，曲率系數(shù)對粗顆粒土變形特性影響較大，當(dāng)其有效粒徑與中值粒徑及限制粒徑的差距過大時，粗粒土中顆粒填充及致密性將變差，變形增大，強(qiáng)度降低。以較大曲率系數(shù)與較大不均勻系數(shù)組合的粗粒土受力性能較好，但此種級配的粗粒土壓縮性大。級配良好且不均勻系數(shù)較大，其應(yīng)力鏈分布越均勻，顆粒的擠壓效果越顯著，力的傳遞和分配也越均勻。

粗粒土 曲率系數(shù) 不均勻系數(shù) 應(yīng)力鏈

0 引 言

粗粒土泛指沖洪積砂卵石、采挖堆石、礫質(zhì)土以及各種軟巖分化料。目前，在高速鐵路填方路基中，填料一般采用粗顆粒填料。根據(jù)郭慶國（1998）在文獻(xiàn)中的結(jié)論，粗顆粒是指粒徑5mm以上的顆粒。幾十年來，專家學(xué)者們對粗粒土的強(qiáng)度、應(yīng)力－應(yīng)變關(guān)系、復(fù)雜應(yīng)力狀態(tài)下的力學(xué)表現(xiàn)做了重點(diǎn)研究，但基本上都習(xí)慣于借助傳統(tǒng)的室內(nèi)試驗(yàn)來研究粗粒土的強(qiáng)度及變形特性。Ahad et al.（2009）研究了粗粒土最大粒徑對抗剪強(qiáng)度的影響；秦紅玉等（2004）通過三軸試驗(yàn)研究了粗粒料中不同泥巖含量對堆石壩粗粒料強(qiáng)度的影響；張嘎等（2004）對粗粒土與結(jié)構(gòu)接觸面在循環(huán)荷載作用下的力學(xué)特性進(jìn)行了系統(tǒng)的試驗(yàn)研究，討論了循環(huán)荷載作用下接觸面受力變形的規(guī)律和機(jī)理；程展林等（2007）通過CT三軸試驗(yàn)，從細(xì)觀角度研究了粗粒土顆粒之間的相互作用以及顆粒的相對位置對其強(qiáng)度及變形性質(zhì)的影響；姜景山等（2009）通過大型三軸壓縮試驗(yàn)研究了密度、圍壓對粗粒土力學(xué)性質(zhì)的影響，并分析了圍壓對應(yīng)力－應(yīng)變曲線的影響；馮大闊等（2009）通過實(shí)驗(yàn)研究，分析了粗粒土與結(jié)構(gòu)的接觸面間的相互作用現(xiàn)象，總結(jié)了在常剛度法向邊界條件下兩者相互作用的基本規(guī)律；周健等（2013）利用空心圓柱儀對礫石土進(jìn)行三軸壓縮、扭剪、三軸－扭剪聯(lián)合3種單調(diào)加載路徑試驗(yàn)以及應(yīng)力主軸旋轉(zhuǎn)路徑試驗(yàn)，研究了應(yīng)力路徑對礫石土力學(xué)性質(zhì)的影響；鄧國棟等（2014）采用大型動三軸儀對高速鐵路路基粗粒土填料進(jìn)行實(shí)驗(yàn)研究，分析應(yīng)力水平、固結(jié)圍壓、振動頻率以及振動次數(shù)對動彈性模量和阻尼比的影響。同時，一些學(xué)者通過改變粗粒土中某顆粒級配段粗粒的含量來研究粗粒含量對粗粒土力學(xué)及變形性質(zhì)的影響。王瑞甫（2014）利用單軸壓縮蠕變試驗(yàn)，探討粗粒土中細(xì)顆粒含量與蠕變參數(shù)之間的關(guān)系。石熊等（2014）通過向素土中摻加不同比例的碎石，對改良后的高鐵路基粗粒土進(jìn)行大型三軸剪切試驗(yàn)，結(jié)果表明級配改良后的粗粒土壓實(shí)特性顯著提高，各主要檢測指標(biāo)也有較大幅度的提高。進(jìn)入21世紀(jì)，離散元方法被國內(nèi)學(xué)者引入來研究巖土材料的力學(xué)性質(zhì)。耿麗等（2011）通過室內(nèi)試驗(yàn)，基于三維離散元顆粒流理論，從細(xì)觀角度出發(fā)，得到按級配生成的粗粒土三軸試驗(yàn)的顆粒流數(shù)值模型，通過對不同圍壓下顆粒流模型與室內(nèi)三軸試驗(yàn)的力學(xué)關(guān)系曲線的對比，研究了影響粗粒土強(qiáng)度的顆粒細(xì)觀參數(shù)；徐肖峰等（2013）通過大型直剪試驗(yàn)和離散元模擬方法研究了粗粒含量對粗粒土宏觀及細(xì)觀性質(zhì)的影響機(jī)理；董輝等（2015）采用顆粒離散單元方法，對影響堆積碎石土宏觀變形特征的細(xì)觀參數(shù)進(jìn)行虛擬試驗(yàn)敏感度分析。粗顆粒填料的選擇優(yōu)化是高速鐵路填方路基工程建設(shè)中的關(guān)鍵技術(shù)之一，而顆粒級配優(yōu)化又是粗顆粒填料選擇優(yōu)化中的重要一項(xiàng)。進(jìn)入21世紀(jì)，顆粒流數(shù)值方法被引入我國，越來越多的學(xué)者用這一方法來研究巖土的力學(xué)性質(zhì)，特別是在對砂土、黏性土及堆石體的研究方面取得了不少的成果，但是利用顆粒流數(shù)值方法對高鐵路基粗粒土力學(xué)性質(zhì)的研究并不多見，尤其是利用離散元法研究顆粒級配對高鐵路基粗粒土受力及變形性質(zhì)的研究尚屬少見。

本文選取武廣高鐵湖南境內(nèi)某地區(qū)的鐵路路基粗粒土進(jìn)行了室內(nèi)大型三軸試驗(yàn)，利用離散元方法，借鑒前人建立顆粒流模型的經(jīng)驗(yàn)，建立一個三軸試驗(yàn)?zāi)Ｐ?，該顆粒流模型能夠較好地模擬粗粒土室內(nèi)三軸試驗(yàn)時的力學(xué)表現(xiàn)。在此基礎(chǔ)上，重新編程，開發(fā)出一種級配生成程序，利用該程序能使生成的粗粒土顆粒流模型的顆粒級配與室內(nèi)試驗(yàn)的顆粒級配一致，并且能夠根據(jù)需要生成任意指定級配的粗粒土模型試樣。參考室內(nèi)三軸試驗(yàn)結(jié)果，標(biāo)定了粗粒土顆粒流模型的細(xì)觀模型參數(shù)，驗(yàn)證了模型的有效性。采用標(biāo)定的細(xì)觀參數(shù)，研究顆粒級配對高鐵路基粗粒土強(qiáng)度與變形性質(zhì)的影響，為高鐵路基粗粒土級配選擇提供一定依據(jù)，也為研究粗顆粒填料在動荷載作用下的力學(xué)表現(xiàn)奠定了一定基礎(chǔ)。本文假定不考慮顆粒破碎以及顆粒形狀，用球形顆粒模型代表粗粒土顆粒。

1 試驗(yàn)方案與結(jié)果分析

1.1 試驗(yàn)方案

本文中大型三軸實(shí)驗(yàn)的原料是粗粒土，選自武漢—廣州高速鐵路湖南段某鐵路填方路基的粗粒土填料。圓柱形三軸試驗(yàn)試樣高度H＝645mm，直徑D＝254mm，粗粒土顆粒粒徑范圍為5～38mm。該實(shí)驗(yàn)為排水實(shí)驗(yàn)，徑向變形通過綁設(shè)在距試樣兩端1／3及1／2高度處的3個測量儀器測得，徑向變形取3個量測值的平均值，體積變形由軸向變形和徑向變形算得。試樣直徑和最大顆粒粒徑比為6.68，試樣尺寸與平均顆粒半徑之比為60.09。根據(jù)劉海濤等（2009）的結(jié)論：當(dāng)試樣尺寸與平均顆粒半徑之比大于40時，試樣尺寸對粗粒土抗剪強(qiáng)度沒有影響。根據(jù)中華人民共和國國家發(fā)展與改革委員會（2007）規(guī)定，試樣直徑與試樣中最大顆粒粒徑之比不應(yīng)小于5，試樣尺寸的選取消除了尺寸效應(yīng)對試驗(yàn)結(jié)果的影響。試樣級配曲線（圖1）。

圖1 顆粒級配曲線圖Fig.1 Grain size distribution curve

本文采用偏應(yīng)力q和靜水應(yīng)力不變量p來描述材料的應(yīng)力-應(yīng)變關(guān)系。應(yīng)力、應(yīng)變變量如下所示：

1.2 試驗(yàn)結(jié)果分析

室內(nèi)大型三軸實(shí)驗(yàn)結(jié)果（圖2），當(dāng)圍壓分別為10.3kPa、41.3kPa、68.9kPa時，試件的應(yīng)力比－q／p與偏應(yīng)變κ、體應(yīng)變εvol與偏應(yīng)變κ的關(guān)系曲線具有相同的走勢。如圖2a所示，隨著偏應(yīng)變的增加，應(yīng)力比值增大，最終趨于穩(wěn)定，偏應(yīng)變在0％～1％的范圍內(nèi)，應(yīng)力比增加很快，偏應(yīng)變在1％之后，應(yīng)力比的增幅幾近于零。不同圍壓下的應(yīng)力比－偏應(yīng)變曲線的形狀走勢相同，隨著圍壓的增大，應(yīng)力比－偏應(yīng)變曲線的穩(wěn)定峰值降低。從圖2b可以看出，圍壓對試樣體應(yīng)變的變化幅度有影響：圍壓越大，體應(yīng)變的變化幅度越小。在偏應(yīng)變κ增加的過程中，粗粒土變形由剪縮變形轉(zhuǎn)為剪脹變形，隨著圍壓的增大，粗粒土剪縮變形段所對應(yīng)的偏應(yīng)變κ范圍變大，即體積變化變緩。

圖2 室內(nèi)試驗(yàn)的應(yīng)力比－偏應(yīng)變、體應(yīng)變－偏應(yīng)變關(guān)系曲線Fig.2 Stress-strain curves in indoor testa.應(yīng)力比－偏應(yīng)變曲線；b.體應(yīng)變－偏應(yīng)變曲線

2 顆粒流數(shù)值模擬研究

2.1 數(shù)值模型生成過程

編寫程序，生成符合試驗(yàn)土樣顆粒級配的數(shù)值模型。首先，將級配曲線分段，得到分段點(diǎn)的粒徑及對應(yīng)的質(zhì)量累計百分?jǐn)?shù)，兩個分段點(diǎn)對應(yīng)的質(zhì)量累計百分?jǐn)?shù)差即為粒徑區(qū)間段土樣占整個土樣的質(zhì)量百分?jǐn)?shù)。然后，分段生成?；就亮５拿芏认嗤|(zhì)量含量即是體積含量，以體積量為控制參數(shù)，以類似半徑膨脹法的方法隨機(jī)生成每個區(qū)段所對應(yīng)的土體顆粒，這些顆粒的位置可以重疊。最后，基于動力學(xué)基本原理，運(yùn)行程序，顆粒相互排斥充斥整個空間形成土樣模型。級配生成法的基本過程如圖3所示，最終生成的是符合試驗(yàn)土樣顆粒級配的顆粒流數(shù)值模型，其中相同級配段的顆粒用同種顏色表示，此模型尺寸與室內(nèi)試驗(yàn)一致，高645mm，直徑254mm。

耿麗（2011）分析了顆粒流細(xì)觀參數(shù)對粗粒土強(qiáng)度的影響，標(biāo)定了粗粒土三軸試驗(yàn)?zāi)Ｐ偷募?xì)觀參數(shù)，本文參考了其標(biāo)定方法及徐金明等（2010）標(biāo)定細(xì)觀參數(shù)的方法。首先調(diào)節(jié)摩擦系數(shù)，通過與試驗(yàn)曲線比較，確定數(shù)值模型的強(qiáng)度峰值；然后通過調(diào)整法向切向剛度比，通過數(shù)值擬合的方法，進(jìn)一步標(biāo)定與室內(nèi)三軸試驗(yàn)結(jié)果相對應(yīng)的細(xì)觀參數(shù)；最終選取一組最能成功擬合室內(nèi)試驗(yàn)結(jié)果（圖4a）的細(xì)觀參數(shù)（表1）。顆粒流數(shù)值模擬的結(jié)果與室內(nèi)三軸試驗(yàn)結(jié)果相比，前者應(yīng)力比最大值之差為后者應(yīng)力比最大值的0.5％，穩(wěn)定后兩者應(yīng)力比值基本相等。同時，對比兩者的體應(yīng)變和偏應(yīng)變的關(guān)系曲線（圖4b），兩條曲線走勢一致，基本吻合，數(shù)值模型的有效性和細(xì)觀參數(shù)的可靠性得到驗(yàn)證。

表1 顆粒流數(shù)值模擬參數(shù)設(shè)置表Table 1 Numerical simulation of particle flow parameters

表2 數(shù)值模型級配信息表Table 2 Information of grain size distribution of numericalmodels

圖4 不同圍壓下顆粒流數(shù)值模型應(yīng)力比－偏應(yīng)變、體應(yīng)變－偏應(yīng)變曲線Fig.4 Stress-strain curves under different confining pressures in numericalmodelsa.應(yīng)力比－偏應(yīng)變曲線；b.體應(yīng)變－偏應(yīng)變曲線

圖3 級配生成法生成模型過程Fig.3 Process of generatingmodel with distribution generation method

2.2 不同顆粒級配下，粗粒土強(qiáng)度特性的顆粒流模擬分析

采用2.1節(jié)的模型生成方法以及標(biāo)定的細(xì)觀參數(shù)，生成其他4種不同級配的粗粒土顆粒流數(shù)值模型，其級配信息見表2，生成的顆粒流模型（圖5）。其中，樣本1代表的是2.1節(jié)所描述的室內(nèi)三軸試驗(yàn)的顆粒流模型，每種顏色代表一段級配，此5個樣本的模型尺寸和顆粒粒徑范圍與室內(nèi)試驗(yàn)一致。樣本生成后，選取3種不同的圍壓（10.3kPa、41.3kPa和68.9kPa）進(jìn)行常規(guī)三軸壓縮試驗(yàn)的顆粒流模擬。計算結(jié)果如圖6、圖7、圖8所示。

圖5 不同級配的數(shù)值模型Fig.5 Numericalmodelswith different grain size distribution

圖6 圍壓為10.3kPa，不同顆粒級配下材料的體應(yīng)變－偏應(yīng)變、應(yīng)力比－偏應(yīng)變關(guān)系曲線Fig.6 Stress-strain curves with different grain composition under confining pressuresσc＝10.3kPa

圖7 圍壓為41.3kPa，不同顆粒級配下材料的體應(yīng)變－偏應(yīng)變、應(yīng)力比－偏應(yīng)變關(guān)系曲線Fig.7 Stress-strain curves with different grain composition under confining pressuresσc＝41.3kPa

2.3 不同顆粒級配下，粗粒土變形特性的顆粒流模擬分析

如圖9所示，圍壓等于10.3kPa時，樣本1兩端的顆粒密集分布，位移帶最長，說明樣本1在兩端加載板附近有較多小球產(chǎn)生了較大的運(yùn)動，同樣的加載條件下，樣本1最易被壓縮，故而加載初期樣本1的體縮現(xiàn)象最為明顯（圖6a）。樣本2的位移場分布最為稀疏且位移矢量模量較其他樣本小，說明在加載過程中樣本2的顆粒移動較小，比較能保持最初的顆粒組構(gòu)，表現(xiàn)出較大強(qiáng)度，故而出現(xiàn)了如圖6b所示的“強(qiáng)度趕超”現(xiàn)象。并且，粗粒土中顆粒位移場的分布規(guī)律遵從其不均勻系數(shù)的大小排列順序，即樣本2的位移場最為稀疏，隨著不均勻系數(shù)的增大，樣本中部位置顆粒的運(yùn)動越活躍，位移場分布越均勻。粗粒土的顆粒模型中其顆粒位移場分布越均勻，說明顆粒的移動范圍越大；位移場疏密程度越均勻，說明發(fā)生相互錯動的顆粒范圍越廣泛，顆粒原先的結(jié)構(gòu)越容易被改變，體變也越大。隨著圍壓的加大，顆粒的位移場分布變得越均勻，在位移場中部，發(fā)生徑向位移的范圍變小，顆粒的軸向運(yùn)動明顯。并且由于圍壓的加大，顆粒移動變得困難，體變速率明顯變小。從圖9可以看出，樣本1體現(xiàn)出較好的受力性能，樣本5體現(xiàn)出較大的壓縮性，表明較大曲率系數(shù)（1、3之間，靠近3）與較大不均勻系數(shù)組合后粗粒土的強(qiáng)度越高，但是其壓縮性也相應(yīng)地增大。

圖8 圍壓為68.9kPa，不同顆粒級配下材料的體應(yīng)變－偏應(yīng)變、應(yīng)力比－偏應(yīng)變關(guān)系曲線Fig.8 Stress-strain curveswith different grain composition under confining pressuresσc＝68.9kPa

圖9 不同圍壓下，不同顆粒流數(shù)值模型中顆粒的位移場變化圖Fig.9 Displacement field distribution of various numerical models under different confining pressures

從表2可以看出，樣本3、4、5的有效粒徑和中值粒徑及限制粒徑的差距較樣本1、2大得多。有效粒徑與中值粒徑及限制粒徑差距很大，表明土體中的細(xì)小顆粒所占比例較大。同一圍壓下，樣本3、4、5的位移場較樣本1、2均勻，圍壓加大，位移場分布越均勻，樣本兩端位移集中區(qū)域變小，說明顆粒相互錯動明顯，粗粒土容易發(fā)生體縮變形，其強(qiáng)度增高。細(xì)小顆粒占的比例過大，粒度分布連續(xù)性差，將會導(dǎo)致空隙空間分布‘不均勻化’，粗粒土中容易出現(xiàn)較大空隙，將不利于空隙的有效填充，不利于致密土顆粒骨架的構(gòu)建，因此就出現(xiàn)了如樣本4這樣級配良好的粗粒土的強(qiáng)度和變形特性反而不如其他級配不良的樣本。根據(jù)上述試驗(yàn)結(jié)果，在1和3之間，隨著曲率系數(shù)的減小，可以近似得出粗粒土變形性能增強(qiáng)的結(jié)論。由此可知，要想得到較致密的粗粒土填料，有效粒徑與中值粒徑及限制粒徑的差距不宜過大。

2.4 不同顆粒級配下，粗粒土顆粒間力鏈變化規(guī)律的顆粒流模擬分析

圖10 不同圍壓下，不同顆粒流數(shù)值模型中顆粒間應(yīng)力鏈的分布圖Fig.10 Stress chain distribution of various numericalmodels under different confining pressures

圖10 所示的是顆粒流模擬試驗(yàn)中顆粒間的應(yīng)力分布圖。其中，紅色部分代表拉應(yīng)力，黑色部分代表壓應(yīng)力，應(yīng)力分布在幾何上呈現(xiàn)鏈狀分布，我們稱之為力鏈。顆粒間應(yīng)力的性質(zhì)、強(qiáng)弱以及其分布狀況在力鏈土中都能得到直觀的體現(xiàn)。圍壓增加，模型試樣的徑向束縛加強(qiáng)，顆粒的徑向運(yùn)動的約束也相應(yīng)的加強(qiáng)，在加載條件不變的情況下，顆粒的擠壓程度將隨圍壓的加大而加強(qiáng)，這一過程表現(xiàn)到力鏈圖上就是圍壓加大，拉應(yīng)力鏈減少，壓應(yīng)力鏈增加。在同一圍壓下，隨著不均勻系數(shù)的增大，拉力鏈減少，壓力鏈增多，且級配越接近良好級配條件，其力鏈分布越均勻，這說明級配越好，顆粒的擠壓效果越顯著，力的傳遞和分配也越均勻。

3 結(jié) 論

（1）采用顆粒流方法模擬分析了粗粒土的強(qiáng)度與變形特性，建立了符合材料真實(shí)級配曲線的顆粒流數(shù)值模型，通過與室內(nèi)三軸試驗(yàn)結(jié)果比較，驗(yàn)證了本文所建顆粒流數(shù)值模型的正確性與有效性。該方法能夠靈活改變試驗(yàn)條件，方便而有效地開展不同試驗(yàn)條件下的數(shù)值模擬。

（2）選取5種不同的顆粒級配，分別建立其相應(yīng)的顆粒流模型，模擬分析了顆粒級配對粗粒土力學(xué)性質(zhì)的影響：曲率系數(shù)對粗粒土的強(qiáng)度與變形特性影響較大，當(dāng)其有效粒徑與中值粒徑及限制粒徑的差距過大時，粗粒土中顆粒填充及組合致密性較差，強(qiáng)度與變形性能下降。

（3）通過顆粒流程序，分析了顆粒的運(yùn)動變化規(guī)律和顆粒間應(yīng)力鏈的變化規(guī)律，較好地描述了粗粒土的細(xì)觀力學(xué)性質(zhì)。隨著圍壓的增加，顆粒位移場的分布變得均勻，在位移場中部，顆粒發(fā)生徑向位移的范圍變小，顆粒的軸向運(yùn)動明顯，體變速率降低。顆粒間壓應(yīng)力鏈逐漸增加，向更加密實(shí)的擠壓狀態(tài)轉(zhuǎn)變。較大曲率系數(shù)（1、3之間，靠近3）與較大不均勻系數(shù)組合后粗粒土的強(qiáng)度與變形性能較好，但此種級配的粗粒土壓縮性增大。不均勻系數(shù)增大，拉力鏈減少，壓力鏈增多，且級配越接近良好級配條件，其力鏈分布越均勻，顆粒的擠壓效果越顯著，力的傳遞和分配也越均勻。

（4）本文的研究，沒有考慮粗粒土顆粒形狀的影響，而事實(shí)上粗粒土顆粒具有隨機(jī)的形狀，顆粒間除了存在摩擦力外還存在較大的“咬合”力。此外，在粗粒土顆粒接觸處，由于顆粒形狀不同，容易產(chǎn)生應(yīng)力集中，會出現(xiàn)某點(diǎn)的接觸力極大，顆粒受壓產(chǎn)生破碎，顆粒破碎后，粗粒土中應(yīng)力將產(chǎn)生重分布，勢必也會影響粗粒土的強(qiáng)度與變形性質(zhì)。所以，綜合考慮顆粒破碎、顆粒形狀以及顆粒級配的共同作用將是今后研究粗粒土細(xì)觀力學(xué)性質(zhì)的一個主要方向。

Ahad Bagherzadeh-khalkhali，Mirghasemi A A.2009.Numerical and experimental direct shear tests for coarse-grained soils［J］. Particuology，7（1）：83～91.

Cheng Z L，Ding H S，Wu L P.2007.Experimental study on mechanical behavior of granular material［J］.Chinese Journal of Geotechnical Engineering，29（8）：1151～1158.

Deng G D，Zhang JS，Wang Q Y，et al.2014.Experimental research on dynamic parameters of high-speed railway coarse-grained padding［J］.Journal of Railway Science and Engineering，2（11）：76～83.

Dong H，Ma Y Y，F(xiàn)u H L，et al.2015.Particles discrete element simulation accumulation of gravel soil deformation parameter sensitivity analysis［J］.Chinese Journal of Computational Mechanics，32（2）：192～199.

Feng D K，Zhang G，Zhang JM，et al.2009.Experimental study on 3D cyclic behaviors of soil-structure interface under constant normal stiffness condition［J］.Chinese JournalofGeotechnical Engineering，31（10）：1571～1577.

Geng L，Huang Z Q，Miao Y.2011.Meso-mechanics simulation triaxial test of coarse-grained soil［J］.Journal of Civil Engineering and Management，28（4）：24～29.

Guo Q G.1998.The engineering properties and application of coarse grained soil［M］.Henan：The Yellow River Water Conservancy Press：6～9.

Jiang JS，Liu H L，Cheng Z L，et al.2009.Influences of density and confining pressure on mechanical properties for coarse grained soils［J］.Journal of Yangtze River Scientific Research Institute，26（8）：46～50.

Liu H T，Cheng X H.2009.Discrete element analysis for size effects of coarse-grained soils［J］.Rock and Soil Mechanics，30（S）：287～292.

Li JZ，Xu JM，Huang D Y.2015.Particle flow simulation of deformation and failure mechanism of granite based on actual distributions of meso-compositions［J］.Journal of Engineering Geology，23（S）：84～90.

National Development and Reform Commission.2007.Code for coarsegrained soil tests for hydropower and water conservancy engineering（DL／T5356—2006）［S］.Beijing：China Electric Power Press.

Qin H Y，Liu H L，Gao Y F，et al.2004.Research on strength and deformation behavior of coarse aggregates based on large-scale triaxial tests［J］.Rock and Soil Mechanics，25（10）：1575～1580.

Shi X，Zhang JS，Liu B，et al.2014.Test of high-speed railway course grained filler of improved particle size distribution［J］.Journal ofCentral South University（Science and Technology），45（11）：3964～3969.

Wang L F.2014.Experimental research on influence of content of fine particles on creep properties of coarse-grained soil［J］.Technology of Highway and Transport，1（5）：1～4.

Xu JM，Xie Z L，Jia H T.2010.Simulation ofmesomechanical properties of limestone using particle flow code［J］.Rock and SoilMechanics，31（S2）：390～395.

Xu X F，Wei H Z，Meng Q S，et al.2013.Dem simulation on effect of coarse gravel content to direct shear strength and deformation characteristics of coarse-grained soil［J］.Journal of Engineering Geology，21（2）：311～316.

Zhang G，Zhang J M.2004.Experimental study on cyclic behavior of interface between soil and structure［J］.Chinese Journal of Geotechnical Engineering，26（2）：254～258.

Zhou J，Wang Z H，Zhang J，et al.2013.Macro-meso research on mechanical behavior of a gravelly soil under various stress paths［J］. Chinese Journal of Rock Mechanics and Engineering，32（8）：1721～1728.

程展林，丁紅順，吳良平.2007.粗粒土試驗(yàn)研究［J］.巖土工程學(xué)報，29（8）：1151～1158.

鄧國棟，張家生，王啟云，等.2014.高速鐵路粗粒土填料動力參數(shù)試驗(yàn)研究［J］.鐵道科學(xué)與工程學(xué)報，2（11）：76～83.

董輝，馬一躍，傅鶴林，等.2015.顆粒離散元模擬堆積碎石土變形的參數(shù)靈敏度分析［J］.計算力學(xué)學(xué)報，32（2）：192～199.

馮大闊，張嘎，張建民，等.2009.常剛度條件下粗粒土與結(jié)構(gòu)接觸面三維力學(xué)特性的試驗(yàn)研究［J］.巖土工程學(xué)報，31（10）：1571～1577.

耿麗，黃志強(qiáng)，苗雨.2011.粗粒土三軸試驗(yàn)的細(xì)觀模擬［J］.土木工程與管理學(xué)報，28（4）：24～29.

郭慶國.1998.粗粒土的工程特性及應(yīng)用［M］.河南：黃河水利出版社：6～9.

姜景山，劉漢龍，程展林，等.2009.密度和圍壓對粗粒土力學(xué)性質(zhì)的影響［J］.長江科學(xué)院院報，26（8）：46～50.

劉海濤，程曉輝.2009.粗粒土尺寸效應(yīng)的離散元分析［J］.巖土力學(xué)，30（增）：287～292.

李進(jìn)昭，徐金明，黃大勇.2015.考慮細(xì)觀組分實(shí)際分布的花崗巖變形破壞過程顆粒流模擬［J］.工程地質(zhì)學(xué)報，23（增）：84～90.

秦紅玉，劉漢龍，高玉峰，等.2004.粗粒料強(qiáng)度和變形的大型三軸試驗(yàn)研究［J］.巖土力學(xué)，25（10）：1575～1580.

石熊，張家生，劉蓓，等.2014.高速鐵路粗粒土填料級配改良試驗(yàn)［J］.中南大學(xué)學(xué)報（自然科學(xué)版），45（11）：3964～3969.

王瑞甫.2014.細(xì)粒含量對粗粒土蠕變特性影響的試驗(yàn)研究［J］.公路交通技術(shù)，1（5）：1～3.

徐金明，謝芝蕾，賈海濤.2010.石灰?guī)r細(xì)觀力學(xué)特性的顆粒流模擬［J］.巖土力學(xué)，31（S2）：390～395.

徐肖峰，魏厚振，孟慶山，等.2013.粗粒含量對礫類土直剪過程中強(qiáng)度與變形特性影響的離散元模擬研究［J］.工程地質(zhì)學(xué)報，21（2）：311～316.

張嘎，張建民.2004.循環(huán)荷載作用下粗粒土與結(jié)構(gòu)接觸面變形特性的試驗(yàn)研究［J］.巖土工程學(xué)報，26（2）：254～258.

中華人民共和國國家發(fā)展與改革委員會.2007.電水利工程粗粒土試驗(yàn)規(guī)程（DL／T5356—2006）［S］.北京：中國電力出版社.

周健，王子寒，張姣，等.2013.不同應(yīng)力路徑下礫石土力學(xué)特性的宏細(xì)觀研究［J］.巖石力學(xué)與工程學(xué)報，32（8）：1721～1728.

STUDY OF STRENGTH AND DEFORMATION CHARACTERISTIC OF UNBOUND GRANUAL MATERIAL W ITH DIFFERENT GRAIN COMPOSITION

CHEN Cheng①②LIU Xiaoqing①LUO Zhengdong①NING Zhiqiang①
（①College of Civil Engineering and Mechanics，Xiangtan University，Xiangtan 411105）
（②State Key Laboratory of Geomechanics and Geotechnical Engineering，Institute of Rock and Soil Mechanics，Chinese Academy of Sciences，Wuhan 430071）

To study the strength and deformation characteristic of unbound granular materials of high-speed railway，a granular flow model is established with the internal FISH language of three-dimensional particle flow code PFC3D.Meanwhile，themicro parameters of unbound granularmaterials are calibrated with the conventional triaxial test results.Considering the influence of the curvature coefficient and the uniformity coefficient on the grain size distribution，five kinds of different grain composition of unbound granular materials are simulated under three different confining pressures.The results show that：There is a predominant factor in the curvature coefficient of unbound granularmaterials.The differences between the effective size，median diameter and constrained size can't be allowed too large，because the strength behaviors of the unbound granular materials with large curvaturecoefficient and uneven coefficient performances are higher than the others，and the deformation behaviors become more and more large.The non-uniformity coefficient is another important factor on the strength and deformation behavior of unbound granular materials.With the larger non-uniformity coefficient and well-graded granular materials won't changing the force chain distribution under the different stress condition，and the extrusion becomes more obvious and the transfer and distribution of force aremore uniform.

Unbound granularmaterial，Curvature coefficient，Non-uniformity coefficient，Stress chain

TU43

：A

10.13544／j.cnki.jeg.2016.06.019

2015－09－06；

2015－12－03.

國家自然科學(xué)基金項(xiàng)目（41372314），湘潭大學(xué)第七屆大學(xué)生研究性學(xué)習(xí)和創(chuàng)新實(shí)驗(yàn)計劃項(xiàng)目資助.

陳鋮（1977－），男，博士，副教授，主要從事巖土材料的本構(gòu)關(guān)系、微細(xì)觀力學(xué)特性與路基填料動力特性等方面的研究.Email：

chencheng401＠sina.com

劉小清（1989－），男，碩士生，主要從事巖土微細(xì)觀力學(xué)特性與路基填料動力特性方面的研究.Email：1295025071＠qq.com