彭安平，李亮，張智

離散元表征的級配碎石加州承載比細(xì)觀力學(xué)行為研究

彭安平1, 2，李亮1，張智1

(1.中南大學(xué) 土木工程學(xué)院，湖南 長沙 410075；2. 湖南路橋建設(shè)集團(tuán)有限責(zé)任公司，湖南 長沙 410075)

基于離散元法構(gòu)建符合真實(shí)級配的石灰石碎石二維細(xì)觀模型，確定模型的細(xì)觀參數(shù)與加載速率，研究級配碎石的加州承載比(CBR)試驗(yàn)。通過對試驗(yàn)和模擬中貫入力?貫入量曲線以及CBR值進(jìn)行對比驗(yàn)證，驗(yàn)證模型的合理性與正確性。分析加載后的主應(yīng)力方向、接觸力鏈和位移矢量等細(xì)觀機(jī)制變化規(guī)律。分析加載速率、摩擦因數(shù)與剛度比等參數(shù)對石灰石碎石CBR值的影響。研究結(jié)果表明：隨著加載速率的降低，CBR數(shù)值基本呈緩慢降低的趨勢；集料的縱向變化明顯于橫向變化，離散元模型能夠較好地反映級配碎石的細(xì)觀作用機(jī)制；顆粒間摩擦因數(shù)的增加能夠明顯提高級配碎石的抗壓能力，碎石體系的抗變形能力隨著法向剛度大小增大而增大。

級配碎石；加州承載比；PFC2D；離散元模型；細(xì)觀結(jié)構(gòu)

級配碎石廣泛用于道路工程建設(shè)當(dāng)中，很多國家都采用加州承載比(CBR)作為衡量路基路面材料抵抗集中荷載的強(qiáng)度設(shè)計(jì)指標(biāo)。然而，由于級配碎石存在顯著的非線性力學(xué)性能和非連續(xù)介質(zhì)結(jié)構(gòu)，受到當(dāng)前力學(xué)理論的限制，無法從理論上詳細(xì)分析級配碎石的力學(xué)參數(shù)、實(shí)驗(yàn)條件與抗壓力變形性能之間的關(guān)系。另外，傳統(tǒng)的室內(nèi)試驗(yàn)方法要求通過重復(fù)大量的測試數(shù)據(jù)，才能夠分析出集料性能與加州承載比之間的規(guī)律來優(yōu)化級配碎石的集料選擇與設(shè)計(jì)，從而使得級配碎石的設(shè)計(jì)周期長、資源浪費(fèi)大[1]。在級配構(gòu)成與CBR方面，長安大學(xué)馬骉 等[2]進(jìn)行了相當(dāng)詳細(xì)的室內(nèi)試驗(yàn)，包括對提高級配碎石的抗變形能力、動(dòng)態(tài)回彈模量、抗剪切性能等方面進(jìn)行了研究；Nunes等[3]利用加州承載比試驗(yàn)研究了用于路面施工的二次碎石材料，得到礦石和石片碎料的加州承載比分別為5%和35%；Mousa等[4]對混合不同含量碎石的再生瀝青路面的回彈模量和加州承載比進(jìn)行了研究；龔璐[5]通過對回彈模量、加州承載比、塑性變形等參數(shù)的影響因素進(jìn)行分析并確定了級配碎石的優(yōu)先級配設(shè)計(jì)等。在級配碎石力學(xué)特性的數(shù)值模擬方面，王豐勝[6]基于KENLAYER的計(jì)算模型對碎石級配的宏觀力學(xué)性能進(jìn)行了一系列研究探索；任皎龍[7]使用有限元數(shù)值模擬技術(shù)，對級配碎石的加州承載比的影響因素進(jìn)行了較為全面的分析。Bigi等[8]研究了不同級配組成的級配碎石在不同含水量狀態(tài)下的動(dòng)三軸試驗(yàn)，同時(shí)對級配碎石在凍融過程中的力學(xué)特性進(jìn)行了分析；Janoo等[9]研究了級配碎石中粗集料的三維形狀、棱角性和表面紋理等因素對其性能的影響，規(guī)定了對集料等級進(jìn)行評價(jià)的相關(guān)指標(biāo)，并分析了不同尺寸式樣的級配組成、集料規(guī)格、含水率和空隙率等參數(shù)對級配碎石性能的影響。鑒于級配碎石復(fù)雜而又相互影響的力學(xué)參數(shù)和加州承載比室內(nèi)實(shí)驗(yàn)無法避免的缺點(diǎn)。以PFC2D為代表的離散元數(shù)值軟件廣泛應(yīng)用于道路材料的力學(xué)特性研究[10?15]。因此，本文研究基于離散元法，對級配碎石CBR試驗(yàn)進(jìn)行模擬和影響因素進(jìn)行研究分析，探討研究級配碎石微細(xì)觀結(jié)構(gòu)的手段。

1 室內(nèi)試驗(yàn)

本文采用CBR值評價(jià)?；鶎硬牧铣休d力，其表征了材料的豎向剛度和抗剪切能力。選用石灰?guī)r質(zhì)集料用作級配碎石其級配如圖1所示，其原因如下：1) 石灰?guī)r在軋制過程中易軋制成所需的形狀，且較易滿足級配要求。2)石灰?guī)r中所含的碳酸鹽在級配碎石基層施工過程中能夠與水發(fā)生水化反應(yīng)，可以在集料中發(fā)揮膠結(jié)作用，從而能夠改善級配碎石力學(xué)特性。通過重型擊實(shí)試驗(yàn)以獲得其最大干密度為2.346 g/cm3及最佳含水率為4.85。最后，在此基礎(chǔ)上進(jìn)級配碎石的CBR試驗(yàn)，試件成型和室內(nèi)試驗(yàn)詳細(xì)描述可以參考文獻(xiàn)[16]。

圖1 石灰?guī)r碎石級配圖

根據(jù)CBR定義可知，其貫入規(guī)定量(2.5 mm和5 mm)時(shí)所需荷載與標(biāo)準(zhǔn)碎石貫入相同量時(shí)所需要荷載的比值，以百分率表示其所得到的CBR值。CBR計(jì)算公式如下：

式中：CBR-加州承載比；-貫入量為2.5 mm或5 mm時(shí)的荷載；P-貫入標(biāo)準(zhǔn)碎石2.5 mm時(shí)的標(biāo)準(zhǔn)荷載(7 000 kPa)或貫入標(biāo)準(zhǔn)碎石5 mm時(shí)的標(biāo)準(zhǔn)荷載(10 500 kPa)。

2 級配碎石CBR數(shù)值模擬方法

2.1 幾何模型的建立

利用顆粒流軟件PFC2D所構(gòu)建的級配碎石加州承載比試驗(yàn)?zāi)Ｐ腿鐖D2所示，首先按照實(shí)際試驗(yàn)的預(yù)定尺寸創(chuàng)建矩形邊界區(qū)域。在其邊界范圍內(nèi)按照級配碎石的設(shè)計(jì)級配創(chuàng)建各檔的集料顆粒。另外，考慮到集料的不規(guī)則形狀和模型計(jì)算速率的影響，將粒徑大于2.36 mm的粗集料利用clump模擬其不規(guī)則的形狀特征，其中clump形狀采用多種不同形狀的集料二維圖形，利用PFC2D內(nèi)置命令將集料圖形導(dǎo)入模型中并生成單一不規(guī)則集料以體現(xiàn)集料的不規(guī)則形態(tài)特性，同時(shí)將代表級配碎石圓盤顆粒隨機(jī)替換成不規(guī)則的clump，以體現(xiàn)集料的隨機(jī)分布以及不規(guī)則形狀特性。而粒徑小于2.36 mm的細(xì)集料統(tǒng)一采用粒徑為2 mm的顆粒模擬其力學(xué)特性。

模型尺寸為150 mm×120 mm，由3 637個(gè)圓盤顆粒和898個(gè)clump組成。由于clump中的pebble顆粒相互之間無接觸作用，故模型中實(shí)際存在4 535個(gè)計(jì)算單元。為模擬CBR試驗(yàn)的加載邊界條件，在頂部生成尺寸大小為50 mm×50 mm的模型如圖2所示，然后在壓頭兩側(cè)生成50 mm寬的clump模擬頂部荷載板作用。在生成期間，不可避免存在大量的顆粒重疊現(xiàn)象，造成顆粒之間的接觸不平衡力過大影響模型計(jì)算的收斂。為減小顆粒間的不平衡力的影響，模型通過不斷循環(huán)作用使顆粒之間的重疊量達(dá)到最小值。

圖2 CBR試驗(yàn)數(shù)值模型

圖3 測量圓分布圖

模型加載后顆粒間會(huì)出現(xiàn)應(yīng)力重分布的現(xiàn)象，應(yīng)力重分布必然會(huì)導(dǎo)致顆粒間應(yīng)力主方向發(fā)生偏轉(zhuǎn)。因此為研究顆粒間應(yīng)力主方向規(guī)律，通過measure命令在模型中設(shè)置一定數(shù)量的半徑為5 mm的測量圓(如圖3)以監(jiān)測各處應(yīng)力大小，以此測量模型中各部分的應(yīng)力大小狀況。

2.2 細(xì)觀力學(xué)參數(shù)的選取

PFC模擬需要通過細(xì)觀尺度顆粒之間的相互作用來反映宏觀尺度級配碎石整體的力學(xué)特性，而模型中利用顆粒間接觸模型反映顆粒間的力學(xué)特性如圖4所示。CBR數(shù)值模型中石灰石碎石可根據(jù)粒徑大小分為粗集料(＞2.36 mm)和細(xì)集料(≤2.36 mm)2種類型。因此，模型中主要存在3種接觸：細(xì)集料內(nèi)部接觸、粗細(xì)集料之間的接觸以及粗集料之間的接觸。

圖4 顆粒間接觸界面

石灰?guī)r碎石中含有碳酸鹽，其在級配碎石施工過程中能夠與水發(fā)生水化反應(yīng)，致使細(xì)集料在模型中能夠起到膠結(jié)作用。因此，利用linear(圖5)模擬粗集料之間的力學(xué)特性，而利用linearcbond(圖6)接觸模型模擬細(xì)集料之間的包括膠結(jié)等力學(xué)作用。

圖5 線性模型

(a) 未黏結(jié)；(b) 黏結(jié)

黏結(jié)時(shí)線性接觸黏結(jié)模型線性部分的力?位移法則：如果接觸受到的法向力超過其抗拉強(qiáng)度時(shí)，法向黏結(jié)破壞且接觸的法向和切向力都會(huì)被重新置零。如果切向力超過其具有的抗剪強(qiáng)度時(shí)，切向黏結(jié)破壞并且摩擦因數(shù)和法向力決定的摩阻力會(huì)限制顆粒的滑動(dòng)。相反，如果接觸所受到的法向和切向力未超過其強(qiáng)度限制，即接觸處于黏結(jié)狀態(tài)時(shí)其力學(xué)特性表現(xiàn)為線性接觸模型。

本文根據(jù)文獻(xiàn)[16]中的級配碎石的宏細(xì)觀參數(shù)。經(jīng)過多次的試算和調(diào)整，選取的模型中接觸細(xì)觀參數(shù)如表1所示。級配碎石的CBR試驗(yàn)屬于準(zhǔn)靜態(tài)問題，一般采用較大的阻尼(β，β)，比如本文采用法向與切向阻尼比均為0.7，以提高計(jì)算的收斂速度，使系統(tǒng)迅速達(dá)到穩(wěn)定狀態(tài)。

2.3 加載條件的模擬

級配碎石CBR試驗(yàn)的加載速率為1~1.25 mm/min，如果模型無法按照真實(shí)的速率計(jì)算。因此，為提高計(jì)算效率，選擇最佳加載速率以保證數(shù)值模擬的準(zhǔn)確性，通過后續(xù)加載速率的分析，選擇0.015 m/s作為模型的加載速率。在加載期間，壓頭以0.015 m/s的恒定速率向集料施加力的作用，直到貫入深度為3.5 mm時(shí)停止加載。另外，為模擬加載試件頂面的6 kg荷載板，給荷載板添加力，使其與重力之等效于6 kg荷載板的作用。

表1 模型中材料細(xì)觀力學(xué)參數(shù)

2.4 可靠性驗(yàn)證

為研究碎石的豎向剛度以及剪切能力，進(jìn)行了級配碎石的CBR室內(nèi)試驗(yàn)和數(shù)值模擬試驗(yàn)。同時(shí)，分別記錄了2類試驗(yàn)各自的貫入力隨貫入量的變化曲線以及貫入量為2.5 mm時(shí)的CBR值。圖7表示貫入深度為3.5 mm時(shí)的模型狀態(tài)。由圖7可知，模擬壓頭下方有限區(qū)域(如圖位移值大的集料)內(nèi)發(fā)生移動(dòng)，并且壓頭的擾動(dòng)作用明顯。

圖7 模型破壞狀態(tài)

圖8表示試驗(yàn)與模擬得到的曲線趨勢對比圖。由圖可知，模擬與試驗(yàn)曲線的發(fā)展規(guī)律基本吻合，其中模擬出現(xiàn)幾次波動(dòng)，主要由于初期作用在顆粒上的力與荷載板所受的力反復(fù)平衡所致。另外，由表2中CBR值相差不大且保持在4.75%。由于是二維模擬，故差值沒有三維模擬小一些，但是二維模擬其變化規(guī)律很有意義。表明通過建立正確的離散元模型和輸入合理的參數(shù)，使用離散法可以準(zhǔn)確的模擬級配碎石的CBR試驗(yàn)。

圖8 CBR試驗(yàn)與數(shù)值模擬曲線對比

表2 CBR室內(nèi)實(shí)驗(yàn)結(jié)果與數(shù)值模擬結(jié)果的對比

3 結(jié)果與分析

3.1 加載速率

由于受到計(jì)算效率的影響，模擬加載速率難以與實(shí)際試驗(yàn)相一致。利用不同的加載速率將影響模擬結(jié)果對試驗(yàn)對比精度從而影響模擬的準(zhǔn)確性。為了選擇最佳加載速率，故研究不同速率對CBR值的影響，其中選取的加載速度為：0.5，0.1，0.05，0.015和0.01 m/s。圖9表示不同加載速率的貫入力與貫入量曲線，圖中隨著加載速率的減小曲線的變化趨勢逐漸相近，尤其是速率0.015 m/s與0.01 m/s變化趨勢較為相近。說明隨著加載速率的減小，貫入力受速度影響所導(dǎo)致的變化逐漸趨于穩(wěn)定。

圖9 不同加載速率的貫入力與貫入量曲線

由圖10可知，級配碎石的CBR值隨著速率的減小逐漸表現(xiàn)出緩慢降低的趨勢。當(dāng)速率大于0.015 m/s時(shí)，貫入量為2.5 mm時(shí)的CBR值隨著速率的減小逐漸降低，說明當(dāng)速率大于0.015 m/s時(shí)，速率明顯會(huì)對CBR值產(chǎn)生影響。其次，當(dāng)速率小于等于0.015 m/s時(shí)，隨著速率的減小拉應(yīng)力峰值大小變化基本趨于平緩，說明當(dāng)速率小于等于0.015 m/s時(shí)，速率的變化量對應(yīng)力峰值影響甚微，甚至可以忽略不計(jì)。因此，基于模型計(jì)算效率和準(zhǔn)確性等方面考慮，在后續(xù)研究瀝青混合料的細(xì)觀參數(shù)對力學(xué)特性的影響中，模型采用0.015 m/s作為最佳加載速率。

圖10 速度對CBR的影響

3.2 模型細(xì)觀機(jī)制分析

當(dāng)壓頭貫入模型時(shí)，顆粒間會(huì)發(fā)生嵌擠、摩擦作用，集料內(nèi)部會(huì)發(fā)生應(yīng)力重分布的現(xiàn)象，應(yīng)力重分布必然會(huì)導(dǎo)致應(yīng)力主方向發(fā)生偏轉(zhuǎn)。因此，基于應(yīng)力主方向的偏轉(zhuǎn)情況、量化力鏈分布以及粗骨料量化速度矢量分布情況，對貫入后模型中的受力形態(tài)以及演變規(guī)律進(jìn)行分析。應(yīng)力主方向0的計(jì)算表達(dá)式為[17]：

式中：，和分別表示水平、切向和豎向應(yīng)力。測量圓能夠監(jiān)測出各個(gè)測量圓范圍內(nèi)水平應(yīng)力、切向應(yīng)力和豎向應(yīng)力大小，然后根據(jù)式(2)計(jì)算出應(yīng)力主方向。為了能方便形象的表示每個(gè)測量圓內(nèi)應(yīng)力主方向分布情況，模型中采用離散元裂隙網(wǎng)絡(luò)表示應(yīng)力主方向的分布狀況。

圖11表示接觸力以及應(yīng)力主方向分布規(guī)律，圖中接觸線條粗細(xì)表示作用力的大小即線條越粗作用力越大。其中圖11(a)表示初始狀態(tài)的接觸力及主應(yīng)力方向分布規(guī)律，此時(shí)壓頭尚未作用于試樣表面，壓頭下部區(qū)域主應(yīng)力方向無規(guī)律分布。當(dāng)壓頭貫入量為2.5 mm時(shí)，壓頭底部區(qū)域會(huì)出現(xiàn)明顯的應(yīng)力集中現(xiàn)象且逐漸向下延伸，而模型兩側(cè)作用力相對較小，接觸力的縱向變化比橫向變化相對明顯。壓頭底部(圖中矩形虛線內(nèi))的主應(yīng)力方向在貫入量為0時(shí)，主應(yīng)力方向分布表現(xiàn)不規(guī)則，而貫入量為2.5 mm時(shí)，由于壓頭貫入的作用集料已經(jīng)形成主要受力骨架，主應(yīng)力方向表現(xiàn)為是豎直規(guī)則排列，壓頭兩側(cè)反而表現(xiàn)為不規(guī)則分布狀態(tài)。是由于壓頭下部區(qū)域主要承壓，而兩側(cè)荷載板的存在顆粒能夠向兩側(cè)向上稍稍移動(dòng)。因此，主應(yīng)力方向能夠在加載中較好的反映模型內(nèi)部的受力方向分布 規(guī)律。

(a) 0.0 mm；(b) 3.5 mm

圖12表示級配碎石在不同貫入量時(shí)主應(yīng)力方向以及粗集料速度矢量分布規(guī)律，圖中箭頭表示顆粒速度矢量大小，箭頭方向?yàn)轭w粒運(yùn)動(dòng)方向，箭頭大小代表顆粒速度大小，箭頭長度越大也就位移越大。當(dāng)壓頭貫入量為0.5 mm時(shí)，速度矢量僅存在壓頭下方及兩側(cè)有限區(qū)域內(nèi)，其它區(qū)域均為點(diǎn)狀分布，壓頭下部的主應(yīng)力方向僅在頂部一小部分區(qū)域規(guī)則豎直排放，表明此時(shí)的集料僅在表明發(fā)生輕微擾動(dòng)；當(dāng)貫入量為1.5 mm時(shí)，壓頭底部區(qū)域的集料速度矢量自壓頭下呈擴(kuò)散型向下方擴(kuò)展并且主應(yīng)力方向豎直向上分布，表面顆?？焖龠w移的同時(shí)也存在主要承載，另外，壓頭兩側(cè)區(qū)域的速度矢量減小，由于荷載板壓力與顆粒受力達(dá)到平衡致使速度放緩。當(dāng)貫入量為2.5 mm時(shí)，速度矢量與主應(yīng)力方向與貫入量為1.5時(shí)類似，致使速度大小在減緩，這是由于集料在貫入力作用下逐漸密實(shí)，同時(shí)模型兩側(cè)速度矢量式中處于點(diǎn)狀分布，集料擾動(dòng)不明顯，表明速度矢量的縱向變化比橫向變化明顯。

(a) 0.005 mm；(b) 0.015 mm；(c) 0.002 5 mm

3.3 摩擦因數(shù)的影響

為研究級配碎石之間摩擦因數(shù)對CBR值的力學(xué)影響，進(jìn)行了一系列不同摩擦因數(shù)(0.12，0.22，0.32，0.42，0.62)的數(shù)值模擬試驗(yàn)。不同摩擦因數(shù)下級配碎石貫入應(yīng)力與貫入量曲線如圖13所示。圖中貫入力曲線隨著摩擦因數(shù)的增大其上升趨勢由平緩逐漸變?yōu)槎盖汀Ｓ蓤D14可知，摩擦因數(shù)在0.12~0.62范圍內(nèi)，CBR值隨著摩擦因數(shù)的增大而逐漸增大。這是由于摩擦因數(shù)越大，集料表面越粗糙，集料間越難以發(fā)生相對移動(dòng)，故集料體系更加穩(wěn)定，且在模型集料體系越穩(wěn)固，級配碎石抵抗荷載的能力就越強(qiáng)，因此級配碎石的CBR值相對也就越大。同時(shí)也間接證明了集料的表面紋理粗糙度對提高級配碎石CBR具有積極意義。

圖13 不同摩擦因數(shù)的貫入力與貫入量曲線

圖14 CBR隨摩擦因數(shù)的變化趨勢

3.4 剛度比的影響

為研究顆粒間法向與切向剛度比對級配碎石細(xì)觀力學(xué)特性的影響，進(jìn)行一系列不同剛度比(0.41，2.41，3.41，4.41，5.41和6.41)的CBR模擬。圖15表示不同剛度比條件下貫入力隨貫入量變化曲線。與摩擦因數(shù)相比，剛度比的變化對貫入力曲線形態(tài)的變化影響相對較小。圖16中CBR值隨著剛度比的增加而逐漸下降，且兩者之間存在一定的線性關(guān)系。是由于在集料彈性模量不變的情況下，CBR隨著剛度比的減小而減小，表明碎石的法向剛度在減小，集料在荷載作用下的抗變形能力降低，導(dǎo)致CBR值降低。但是剛度比代表著碎石本身的力學(xué)特性，并不影響顆粒的表面接觸狀態(tài)，進(jìn)而不會(huì)對顆粒間的相對移動(dòng)有所影響，故貫入力曲線形態(tài)及CBR值大小相差不大。

圖15 不同剛度比的貫入量與貫入力曲線

圖16 CBR隨剛度比的變化趨勢

4 結(jié)論

1) 級配碎石CBR值隨著速率的減小逐漸呈緩慢降低的趨勢，當(dāng)速率小于0.015 m/s時(shí)，速率對間接拉伸強(qiáng)度影響很小。故0.015 m/s作為最佳速率能夠節(jié)省計(jì)算時(shí)間。

2) 隨著貫入量的增大，主應(yīng)力方向以豎直向上的方式規(guī)則排列在壓頭作用區(qū)域，細(xì)觀接觸力在壓頭處應(yīng)力集中并向底部擴(kuò)展，速度矢量自壓頭下呈擴(kuò)散型向下方擴(kuò)展。表明縱向變化明顯于橫向，模型能夠較好地反映級配碎石的細(xì)觀作用機(jī)制。

3) 貫入力隨著摩擦因數(shù)的增大而增大，其趨勢由平緩變陡峭；CBR值隨著摩擦因數(shù)的增大逐漸增大且兩者存在一定的線性關(guān)系。間接證明集料的表面紋理粗糙度對提高級配碎石CBR具有積極意義。

4) 隨著剛度比的增大，級配碎石CBR值逐漸下降且幅度不大，但貫入力曲線形態(tài)的變化不大。表明碎石體系的抗變形能力主要與法向剛度大小有關(guān)。

由于計(jì)算效率的限制，難以細(xì)致的模擬2 mm以下的集料和礦粉等微小顆粒，從而無法研究兩者的最佳用量，將在以后通過高性能計(jì)算機(jī)加以解決，這對推動(dòng)級配碎石科學(xué)技術(shù)的發(fā)展具有重要意義。

[1] 謝遠(yuǎn)勇, 李改. 級配碎石動(dòng)態(tài)CBR試驗(yàn)研究[J]. 公路工程, 2014, 39(2): 5?8. XIE Yuanyong, LI Gai. Experimental research on dynamic CBR of aggregate crushed stone[J]. Highway Engineering, 2014, 39(2): 5?8.

[2] 馬骉, 王秉綱. 基于抗變形能力的級配碎石組成設(shè)計(jì)方法[J]. 長安大學(xué)學(xué)報(bào)(自然科學(xué)版), 2007, 27(5): 1?5. MA Biao, WANG Binggang. Graded and broken stone composing design method based on deformation resistance[J]. Journal of Chang’an University (Natural Science Edition), 2007, 27(5): 1?5.

[3] Nunes M C M, Bridges M G, Dawson A R. Assessment of secondary materials for pavement construction: Technical and environmental aspects[J]. Waste Management, 1996, 16(1/3): 87?96.

[4] Mousa E, Azam A, Elshabrawy M, et al. Laboratory characterization of reclaimed asphalt pavement for road construction in Egypt[J]. Canadian Journal of Civil Engineering, 2017, 44(2): 134?137.

[5] 龔璐. 級配碎石基層級配設(shè)計(jì)及應(yīng)用研究[D]. 長沙: 長沙理工大學(xué), 2008. GONG Lu. The research on grading design and application of unbound stone aggregates base[D]. Changsha: Changsha University of Science & Technology, 2008.

[6] 王豐勝. 影響級配碎石模量的結(jié)構(gòu)因素敏感性分析[J]. 合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2009, 32(9): 1414? 1417. WANG Fengsheng. Sensitivity analysis of the factors influencing the modulus of graded broken stone[J]. Journal of Hefei University of Technology (Natural Science Edition), 2009, 32(9): 1414?1417.

[7] 任皎龍. 級配碎石CBR數(shù)值試驗(yàn)方法及其應(yīng)用[D]. 西安: 長安大學(xué), 2011. REN Jiaolong. CBR simulation test method of graded broken stone and its appliance[D]. Xi’an: Chang’an University, 2011.

[8] Bigi S R, Berg R L. Material testing and initial pavement design modeling[J]. Material Testing & Initial Pavement Design Modeling, 1996, 23(3): 123?128.

[9] Janoo V C, Bayer J J. The effect of aggregate angularity on base course performance[J]. Effect of Aggregate Angularity on Base Course Performance, 2001, 33(1): 23?29.

[10] 王端宜, 趙熙. 瀝青混合料單軸壓縮試驗(yàn)的離散元仿真[J]. 華南理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2009, 37(7): 37?41. WANG Duanyi, ZHAO Xi. Simulation of uniaxial compression test for asphalt mixture with discrete element method[J]. Journal of South China University of Technology (Natural Science Edition), 2009, 37(7): 37? 41.

[11] Khattak M J, Khattab A, Rizvi H R, et al. Imaged-based discrete element modeling of hot mix asphalt mixtures[J]. Materials & Structures, 2015, 48(8): 2417?2430.

[12] 蔣應(yīng)軍, 李思超, 王天林. 級配碎石動(dòng)三軸試驗(yàn)的數(shù)值模擬方法[J]. 東南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 43(3): 604?609. JIANG Yingjun, LI Sichao, WANG Tianlin. Numerical method of dynamic triaxial test on graded crushed rock[J]. Journal of Southeast University (Natural Science Edition), 2013, 43(3): 604?609.

[13] Zeghal M. Discrete element method investigation of the resilient behavior of granular materials[J]. Journal of Transportation Engineering, 2004, 130(4): 503?509.

[14] XIAO Y, Tutumluer E. Gradation and packing characteristics affecting stability of granular materials: Aggregate imaging-based discrete element modeling approach[J]. International Journal of Geomechanics, 2016, 17(3): 0401?0464.

[15] 蔣應(yīng)軍, 任皎龍, 李頔, 等. 礦料強(qiáng)嵌擠骨架密實(shí)級配的PFC2D數(shù)值試驗(yàn)研究[J]. 華南理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 40(2): 92?98. JIANG Yingjun, REN Jiaolong, LI Di, et al. Investigation into close-grained gradation of aggregates skeleton under strong interlocked force via PFC2D numerical test[J]. Journal of South China University of Technology (Natural Science Edition), 2012, 40(2): 92?98.

[16] ZHANG J U, LI J, YAO Y S. Geometric anisotropy modeling and shear behavior evaluation of graded crushed rocks[J]. Construction and Building Materials, 2018, 183(3): 346?355.

[17] 賴漢江. 樁承式加筋路堤土拱形成—演化機(jī)理及承載特性宏細(xì)觀分析[D]. 武漢: 華中科技大學(xué), 2016. LAI Hanjian. Analysis of formation-evolution mechanism and bearing characteristics of soil arching within a piled embankment[D]. Wuhan: Huazhong University of Science & Technology, 2016.

Study on the meso-mechanical behavior of California bearing ratio of graded gravel with discrete element method

PENG Anping1, 2, LI Liang1, ZHANG Zhi1

(1. School of Civil Engineering, Central South University, Changsha 410075, China;2. Hunan Road & Bridge Construction Group Co., Ltd, Changsha 410075, China)

A two-dimensional discrete element model of a graded gravel material with real graduation is developed based on discrete element software PFC2D to simulate and analyze the California bearing ratio (CBR) of a graded gravel material. The results of the force-penetration value curve and CBR value from the discrete element simulation were compared with actual test results and analyze the error and its source of the numerical model, which validated reasonability of the model. Then, the model was applied to study the micro-mechanism variation rule of the main stress direction, contact force chain and displacement vector of the CBR test of a graded gravel material, e.g., it also study the effects of various loading speeds, different friction coefficients and stiffness ratio of particles on the CBR of the material. The results indicate that the CBR value decreases slowly with the decrease of loading rate. The longitudinal variation of aggregate is more obvious than the transverse variation, and the discrete element model can better reflect the micro-mechanism of graded aggregates In addition, the increase of friction coefficient between particles cannot obviously improve the compressive ability of the graded gravel and the deformation resistance of gravel system increases with the increase of normal stiffness.

graded gravel; California bearing ratio; PFC2D; discrete element model; micro-structure

U414.103

A

1672 ? 7029(2019)10? 2467 ? 08

10.19713/j.cnki.43?1423/u.2019.10.012

2019?01?16

湖南省交通廳交通項(xiàng)目(2016-22)；國家自然科學(xué)基金資助項(xiàng)目(51308554)

彭安平(1974?)，男，湖南邵陽人，高級工程師，從事路基路面工程方面的科學(xué)研究；E?mail：663357762@qq.com

(編輯 涂鵬)