景宏君 張延青 顧行文 凌華

摘要：為研究土石混合填料在不同土石混合比情況下的力學(xué)響應(yīng)，以國(guó)道316旬陽至安康二級(jí)公路改建工程中的填石高路堤工程為依托，利用大型三軸剪切試驗(yàn)儀，在4種不同圍壓條件下，分別對(duì)土石混合質(zhì)量比為3∶7和3∶2的試樣進(jìn)行固結(jié)不排水剪切試驗(yàn)。發(fā)現(xiàn)當(dāng)土石混合比為3∶7時(shí)，在低圍壓（200 kPa）條件下，試樣應(yīng)力-應(yīng)變曲線為應(yīng)變硬化型，在高圍壓（400，600和800 kPa）條件下，試樣應(yīng)力-應(yīng)變曲線在達(dá)到峰值應(yīng)力后有所下降，表現(xiàn)為弱應(yīng)變軟化。在4種圍壓條件下，試樣體應(yīng)變均表現(xiàn)為隨軸向應(yīng)變先增大后減小，即試樣先剪縮后剪脹，且表現(xiàn)為隨著圍壓的增大，剪縮量不斷增大，然后達(dá)到峰值，隨后體應(yīng)變減小，表現(xiàn)出剪脹趨勢(shì)，且隨著圍壓的增大，剪脹趨勢(shì)減弱。當(dāng)土石混合比為3∶2時(shí)，在4種圍壓條件下，試樣應(yīng)力-應(yīng)變曲線均表現(xiàn)為應(yīng)變硬化，無明顯的峰值強(qiáng)度。試樣體應(yīng)變均隨軸向應(yīng)變?cè)龃蠖龃?，即試樣一直表現(xiàn)為剪縮，且圍壓越大，剪縮趨勢(shì)越明顯，剪縮量越大。體變速率 （dεv/dε1）能較好地反映試樣變形特性和應(yīng)力-應(yīng)變特性。結(jié)果表明：土石混合比的變化對(duì)土石混合料剪脹性影響很大，對(duì)應(yīng)力-應(yīng)變特性的影響主要體現(xiàn)在加載后期。

關(guān)鍵詞：土石混合料;大型三軸剪切試驗(yàn);應(yīng)力-應(yīng)變關(guān)系;剪脹性;體變速率

中圖分類號(hào)：TU 411文獻(xiàn)標(biāo)志碼：A

DOI：10.13800/j.cnki.xakjdxxb.2019.0213文章編號(hào)：1672-9315（2019）02-0270-06

0引言

在山區(qū)修筑道路，由于缺乏優(yōu)良的路基填料，常選用隧道爆破或路塹開挖產(chǎn)生的棄渣。若單純使用棄渣作為填料，則路基難以壓實(shí)且穩(wěn)定性差，因此往往利用棄渣與細(xì)粒土混合形成的土石混合體作為路基填料。土石混合體通常是指由具有一定工程尺度、強(qiáng)度較高的塊石、細(xì)粒土及孔隙構(gòu)成且具有一定含石量的極端不均勻松散巖土介質(zhì)系統(tǒng)[1]。國(guó)內(nèi)眾多學(xué)者對(duì)土石混合體工程力學(xué)特性進(jìn)行研究，并取得了一定成果。曹文貴等考慮含水量、含石量、巖性及土性等因素的影響，進(jìn)行了一系列土石混填體大型三軸試驗(yàn)，研究土石混填體的變形力學(xué)特性[2]。丁自偉等利用改進(jìn)的原位直剪試驗(yàn)裝置，開展礦區(qū)典型巖土體力學(xué)特性研究[3]。寧金成等采用大型三軸試驗(yàn)研究了土石混合體的強(qiáng)度特性[4]。劉麗萍等為研究擊實(shí)功、擊實(shí)方法及等效處理法對(duì)土石混合料擊實(shí)結(jié)果的影響，開展了一系列擊實(shí)試驗(yàn)[5]。趙川等采用直剪試驗(yàn)對(duì)碎石土的抗剪強(qiáng)度與含石量、含水量等的關(guān)系進(jìn)行了研究[6-8]。孫文君等采用大三軸試驗(yàn)研究了黃土與碎石混合體的力學(xué)性質(zhì)[9]。王江營(yíng)等開展一系列大型三軸剪切試驗(yàn)，研究了不同應(yīng)力路徑下土石混填體的變形力學(xué)特性[10]。金磊等為揭示塊石含量和塊石形狀對(duì)膠結(jié)土石混合體力學(xué)特性的影響，開展了一系列室內(nèi)大三軸試驗(yàn)研究[11]。閆云明等為研究粗粒含量對(duì)土石混合料抗剪強(qiáng)度的影響，對(duì)不同粗粒含量的土石混合料試樣進(jìn)行飽和固結(jié)排水三軸壓縮試驗(yàn)[12]。舒志樂等采用大型三軸試驗(yàn)，研究了不同粒度分維不同圍壓下土石混合體的應(yīng)力-應(yīng)變特性、強(qiáng)度參數(shù)特征及峰值強(qiáng)度特征[13]。楊忠平等、詹永祥等基于大型直剪試驗(yàn)和離散元程序，研究了含石量變化對(duì)土石混合體抗剪強(qiáng)度的影響[14-15]。夏加國(guó)等進(jìn)行大型三軸固結(jié)不排水剪切試驗(yàn)，研究了土石混合體在含超徑顆粒情況下的力學(xué)響應(yīng)[16]。王強(qiáng)等開展室內(nèi)大型擊實(shí)試驗(yàn)與數(shù)值模擬，研究不同塊石質(zhì)量分?jǐn)?shù)及不同塊石形狀的軟巖土石混合體的擊實(shí)特性及其細(xì)觀機(jī)理[17]。國(guó)外學(xué)者也對(duì)土石混合料的變形和強(qiáng)度特性做了較多研究，Zhang等開展三軸試驗(yàn)，研究了具有相同含石量但不同粒度分布的3種土石混合體的力學(xué)行為[18]。Shan等基于數(shù)字圖像，提出了一種提出了一種土石混合體微觀結(jié)構(gòu)PFC2D模型的自動(dòng)生成方法[19]。Li等、Xu等進(jìn)行了一系列原位水平推剪試驗(yàn)，研究了土石混合體的強(qiáng)度特性[20]。

目前，通過試驗(yàn)和數(shù)值模擬手段對(duì)土石混合料的強(qiáng)度和變形等力學(xué)性狀進(jìn)行了許多研究。但由于土石混合料顆粒粒徑變化較大，變形力學(xué)特性較為復(fù)雜，在用作路堤填料時(shí)，易出現(xiàn)施工困難、穩(wěn)定性差等問題。因此，深入研究土石混合體變形力學(xué)特性，具有重要的理論與工程實(shí)際意義。文中以國(guó)道316旬陽至安康二級(jí)公路改建工程中的填石高路堤工程為依托，進(jìn)行大型三軸固結(jié)不排水剪切試驗(yàn)，研究土石混合填料力學(xué)特性，為本項(xiàng)目工程及其他類似工程的穩(wěn)定性分析提供參考。

1試驗(yàn)

1.1試驗(yàn)設(shè)備和試驗(yàn)材料

本次三軸剪切試驗(yàn)是在南京水利科學(xué)研究院的大型靜力三軸儀（圖1）上進(jìn)行的，試樣尺寸為300 mm×700 mm，最大限制粒徑為60 mm.

試驗(yàn)材料取自國(guó)道316旬陽至安康二級(jí)公路改建工程中的填石高路堤工程，為黃土和絹云石英片巖的混合體。試驗(yàn)級(jí)配曲線如圖2所示。

1.2試驗(yàn)過程

根據(jù)試驗(yàn)要求，按0～5，5～10，10～20，20～40，40～60 mm粒徑檔分5等份制備試樣。將透水板置于試樣底部，橡皮膜綁扎在底座上，之后安裝成型筒，將橡皮膜外翻在成型筒上并在成型筒外抽氣，使橡皮膜貼緊成型筒內(nèi)壁。放入第1層試樣并均勻鋪平，根據(jù)試樣要求的干容重采用振動(dòng)器進(jìn)行振實(shí)，振動(dòng)器底板靜壓為14 kPa，振動(dòng)頻率為40 Hz.試樣裝好后，加上透水板和試樣帽，扎緊橡皮膜，去掉成型筒（去掉成膜桶后的三軸試樣如圖3所示），安裝圍壓室，打開圍壓室排氣孔，向內(nèi)注水，當(dāng)有水從排氣孔流出時(shí)，關(guān)閉排氣孔，采用滴水飽和法使試樣飽和。按要求對(duì)飽和試樣施加圍壓，進(jìn)行固結(jié)。固結(jié)完成后，進(jìn)行試驗(yàn)。

重復(fù)上述過程分別進(jìn)行圍壓為200，400，600和800 kPa狀態(tài)下的試驗(yàn)研究。在試驗(yàn)過程中，計(jì)算機(jī)自動(dòng)采集試樣的軸向荷載、軸向變形和排水量，并同步繪制應(yīng)力-應(yīng)變曲線，直到試樣破壞或達(dá)到軸向應(yīng)變的15%.當(dāng)應(yīng)力-應(yīng)變曲線有峰值時(shí)，將峰值點(diǎn)作為破壞點(diǎn)，其對(duì)應(yīng)的主應(yīng)力差（σ1-σ3）作為該試樣的破壞強(qiáng)度，反之則取軸向應(yīng)變的15%所對(duì)應(yīng)的點(diǎn)為破壞點(diǎn)，對(duì)應(yīng)的主應(yīng)力差（σ1-σ3）為試樣的破壞強(qiáng)度。剪切完成后的試樣如圖4所示。

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

2.1應(yīng)力-應(yīng)變關(guān)系分析

圖5為不同混合比試樣在不同圍條件下的大型三軸固結(jié)不排水剪切試驗(yàn)的偏應(yīng)力-軸向應(yīng)變關(guān)系曲線。

從圖5可以看出：當(dāng)土石混合比一定且軸向應(yīng)變相同時(shí)，試樣的偏應(yīng)力隨圍壓的增大而增大。土石混合比為3∶7的試樣，當(dāng)圍壓為200 kPa時(shí)，其偏應(yīng)力-軸向應(yīng)變曲線為應(yīng)變硬化型，當(dāng)圍壓為400，600和800 kPa時(shí)，其偏應(yīng)力-軸向應(yīng)變曲線表現(xiàn)為弱應(yīng)變軟化，且軸向應(yīng)變未達(dá)到15%時(shí)，試樣已破壞。土石混合比為3∶2的試樣，在4種圍壓下，其偏應(yīng)力-軸向應(yīng)變曲線均表現(xiàn)為應(yīng)變硬化。

土石混合比為3∶7的試樣，其含石量較高，隨著圍壓的增大，石料被壓碎，導(dǎo)致偏應(yīng)力-軸向應(yīng)變曲線在達(dá)到峰值應(yīng)力后有少許下降。土石混合比為3∶2的試樣，其土料含量較高，顆粒不易被壓碎。

2.2體應(yīng)變分析

對(duì)于三軸排水剪切試驗(yàn)而言，試樣的體應(yīng)變是通過量測(cè)試樣的排水量而得到的，而在不排水剪切條件下，通常認(rèn)為試樣體積不發(fā)生變化。但對(duì)于含石量較高的試樣，由于在剪切過程中塊石間存在滑動(dòng)擠壓、破碎和填充空隙等現(xiàn)象，試樣同樣存在剪脹和剪縮現(xiàn)象。

由圖6可以看出：對(duì)于土石混合比為3∶7的情況，試樣在4種圍壓下，體應(yīng)變均表現(xiàn)為隨軸向應(yīng)變先增大后減小，說明試樣先剪縮后剪脹，且圍壓越大，剪縮量越大，然后達(dá)到峰值，隨后表現(xiàn)出剪脹趨勢(shì)，且圍壓越大剪脹量越小，隨著軸向荷載的不斷增大，試樣發(fā)生破壞。對(duì)于土石混合比為3∶2情況，試樣在4種圍壓條件下，體應(yīng)變均隨軸向應(yīng)變?cè)龃蠖龃螅瑹o明顯的峰值強(qiáng)度，表明試樣一直表現(xiàn)為剪縮，且圍壓越大，剪縮量越大。

為了更直觀的表現(xiàn)試樣的變形性狀，引入體變速率 （dεv/dε1）。圖7為體變速率隨軸向應(yīng)變關(guān)系曲線。體變速率為正表示剪縮，體變速率為負(fù)表示剪脹。

由圖7可看出：對(duì)于土石混合比為3∶7的試樣，在4種圍壓下，其體變速率均隨著軸向應(yīng)變的增大從某一正值先增大后減小為負(fù)值，隨后又有略微增大但仍為負(fù)值，說明試樣是先剪縮后剪脹。其曲線的最小值點(diǎn)基本對(duì)應(yīng)于應(yīng)變軟化型曲線的峰值強(qiáng)度點(diǎn)。當(dāng)圍壓一定時(shí)，隨著軸向應(yīng)變的增大，剪縮趨勢(shì)逐漸減弱，剪縮量減小為0時(shí)，開始出現(xiàn)剪脹現(xiàn)象，隨后剪脹趨勢(shì)不斷加強(qiáng)。當(dāng)圍壓較低（200 kPa）時(shí)，在軸向應(yīng)變較?。s5%）時(shí)開始出現(xiàn)剪脹現(xiàn)象，隨著圍壓的增大，出現(xiàn)剪脹現(xiàn)象時(shí)的軸向應(yīng)變隨之增大。當(dāng)軸向應(yīng)變一定時(shí)，隨著圍壓的增大，剪縮量不斷增大，剪脹量不斷減小。對(duì)于土石混合比為3∶2的試樣，加載初期，其體變速率與混合比為3∶7的試樣具有相同的趨勢(shì)，即體變速率隨著軸向應(yīng)變的增大從某一正值先增大后減小，但最終并未減小為負(fù)值，說明試樣一直表現(xiàn)為剪縮。應(yīng)力-應(yīng)變曲線為硬化型。當(dāng)圍壓一定時(shí)，隨著軸向應(yīng)變的增大，剪縮量不斷減小。當(dāng)軸向應(yīng)變一定時(shí)，隨著圍壓的增大，剪縮量不斷增大。可以明顯看出，隨著軸向應(yīng)變的增大，土石混合比為3∶7的試樣比土石混合比為3∶2的試樣較早地發(fā)生破壞。

當(dāng)土石混合比為3∶2時(shí)，試樣土料含量較高，塊石之間相互接觸較少，加載初期，在圍壓作用下，細(xì)顆粒迅速調(diào)整位置，試樣內(nèi)部孔隙被壓密，表現(xiàn)為剪縮。圍壓增大，試樣被壓密程度增加，剪縮量增大。隨著試樣不斷被壓密，剪縮量增大量逐漸減小，表現(xiàn)為剪縮趨勢(shì)減弱。

當(dāng)土石混合比為3∶7時(shí)，試樣石料含量較高，塊石間互相接觸的機(jī)會(huì)增多，在圍壓作用下，塊石快速調(diào)整位置，彼此間相互支撐形成不易破壞的骨架結(jié)構(gòu)。當(dāng)圍壓較低時(shí)，塊石不易被壓碎，骨架結(jié)構(gòu)不會(huì)破壞，由于骨架結(jié)構(gòu)的支撐作用，試樣內(nèi)部的孔隙不易被填充，因此，在加載初期細(xì)顆粒只填充部分孔隙，試樣的剪縮量較小。隨著軸向荷載的增大，細(xì)顆粒間出現(xiàn)滾動(dòng)爬坡現(xiàn)象，試樣開始表現(xiàn)出剪脹趨勢(shì)，體應(yīng)變減小。隨著圍壓的增大，剪縮量不斷增大，塊石之間出現(xiàn)互相擠壓現(xiàn)象，部分塊石擠入其他塊石之間，骨架結(jié)構(gòu)破壞，試樣表現(xiàn)出剪脹趨勢(shì)。之后，圍壓再增大，塊石間互相接觸的部位出現(xiàn)應(yīng)力集中，部分塊石被壓碎為細(xì)顆粒，而細(xì)顆粒迅速調(diào)整位置填充孔隙，致使試樣密實(shí)，表現(xiàn)為圍壓增大，剪脹趨勢(shì)減弱。

3結(jié)論

1）土石混合料的應(yīng)力-應(yīng)變關(guān)系與土和石所占比例及外荷載大小有很大關(guān)系。土石混合料中石料含量較高時(shí)，低圍壓下，其應(yīng)力-應(yīng)變曲線表現(xiàn)為應(yīng)變硬化，高圍壓下，其應(yīng)力-應(yīng)變曲線表現(xiàn)為弱應(yīng)變軟化。土石混合料中土料含量較高時(shí)，其應(yīng)力-應(yīng)變曲線均表現(xiàn)為應(yīng)變硬化;

2）土石混合料中土料和石料比例的變化對(duì)其剪脹性影響很大，當(dāng)石料含量較高時(shí)，試樣出現(xiàn)明顯的剪脹現(xiàn)象，且剪脹量較大。相比土料含量較高的試樣，石料含量較高的試樣在軸向荷載較小時(shí)發(fā)生破壞;

3）體變速率能較好地反映體變變化和應(yīng)力-應(yīng)變關(guān)系。當(dāng)體變速率從正值先增大后減小為負(fù)值，隨后又有所增大但仍為負(fù)值，則體變?yōu)橄燃艨s后剪脹，應(yīng)力-應(yīng)變曲線為軟化型，體變速率為0時(shí)，開始出現(xiàn)剪脹，體變速率達(dá)到最小值時(shí)，試樣剪脹性最大，對(duì)應(yīng)于應(yīng)力-應(yīng)變曲線的峰值強(qiáng)度。當(dāng)體變速率從正值先增大后減小但仍為正值時(shí)，則體變表現(xiàn)為剪縮，應(yīng)力-應(yīng)變曲線為硬化型。

參考文獻(xiàn)（References）：

[1]徐文杰，胡瑞林.土石混合體概念、分類及意義[J].水文地質(zhì)工程地質(zhì)，2009，36（4）：50-56.

XU Wenjie，HU Ruilin.Conception，classification and significations of soilrock mixture[J].Hydrogeology & Engineering Geology，2009，36（4）：50-56.

[2]曹文貴，黃文健，王江營(yíng)，等.土石混填體變形力學(xué)特性大型三軸試驗(yàn)研究[J].湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2016，43（3）：142-148.

CAO Wengui，HUANG Wenjian，WANG Jiangying，et al.Largescale triaxial test study on deformation and intensity characteristics of soilrock aggregate mixture[J].Journal of Hunan University （Natural Sciences），2016，43（3）：142-148.

[3]丁自偉，錢坤.基于原位直剪試驗(yàn)的巖土體力學(xué)特性研究[J].西安科技大學(xué)學(xué)報(bào)，2017，37（1）：32-37.

DING Ziwei，QIAN Kun.Mechanic characteristics of rock and soil mass based on insitu direct shear test[J].Journal of Xi’an University of Science and Technology，2017，37（1）：32-37.

[4]寧金成，孫久民.土石混合體的力學(xué)性能影響因素研究[J].中外公路，2012，32（2）：207-210.

NING Jincheng，SUI Jiumin.Study on influence factors of mechanical properties of soilrock mixture[J].Journal of China & Foreign Highway，2012，32（2）：207-210.

[5]劉麗萍，折學(xué)森.土石混合料壓實(shí)特性試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào)，2006，25（1）：206-210.

LIU Liping，SHE Xuesen.Study on compaction property of earthrock mixture[J].Chinese Journal of Rock Mechanics and Engineering，2006，25（1）：206-210.

[6]趙川，石晉旭，唐紅梅.三峽庫區(qū)土石比對(duì)土體強(qiáng)度參數(shù)影響規(guī)律的試驗(yàn)研究[J].公路，2006（11）：32-35.

ZHAO Chuan，SHI Jinxu，TANG Hongmei.Test and research on influence law of soilrock ratio in Three Gorges Reservoir Region on soil strength parameters[J].Highway，2006（11）：32-35.

[7]董云.土石混合料強(qiáng)度特性的試驗(yàn)研究[J].巖土力學(xué)，2007，28（6）：1269-1274.

DONG Yun.Experimental study on intensity character of rocksoil aggregate mixture[J].Rock and Soil Mechanics，2007，28（6）：1269-1274.

[8]田恒蛟，薛亞東.路用土石混合料的剪切特性試驗(yàn)研究[J].地下空間與工程學(xué)報(bào)，2014，10（1）：30-35，50.

TIAN Hengjiao，XUE Yadong.Experimental research on shear behavior of soilrock mixture for road construction[J].Chinese Journal of Underground Space and Engineering，2014，10（1）：30-35，50.

[9]孫文君，宋楊，路維，等.黃土-碎石混合體的力學(xué)性質(zhì)試驗(yàn)研究[J].公路，2016，（5）：178-181.

SUN Wenjun，SONG Yang，LU Wei，et al.Experimental study on mechanical properties of loesscrushed stone mixture[J].Highway，2016，（5）：178-181.

[10]王江營(yíng)，曹文貴，蔣中明，等.不同應(yīng)力路徑下土石混填體變形力學(xué)特性大型三軸試驗(yàn)研究[J].巖土力學(xué)，2016，37（2）：424-430.

WANG Jiangying，CAO Wengui，JIANG Zhongming，et al.Largescale triaxial tests on deformation and mechanical behavior of soilrock aggregate mixture under different stress paths[J].Rock and Soil Mechanics，2016，37（2）：424-430.

[11]金磊，曾亞武，張森.塊石含量及形狀對(duì)膠結(jié)土石混合體力學(xué)性能影響的大型三軸試驗(yàn)[J].巖土力學(xué)，2017，38（1）：141-149.

JIN Lei，ZENG Yawu，ZHANG Sen.Large scale triaxial tests on effects of rock block proportion and shape on mechanical properties of cemented soilrock mixture[J].Rock and Soil Mechanics，2017，38（1）：141-149.

[12]閆云明，王雯璐.粗粒摻量對(duì)土石混合料強(qiáng)度特性的影響[J].科技通報(bào)，2017，33（9）：142-145.

YAN Yunming，WANG Wenlu.Effects of coarse fraction content and grain size on strength properties of soil and rock mixtures[J].Bulletin of Science and Technology，2017，33（9）：142-145.

[13]舒志樂，劉新榮，劉保縣，等.基于分形理論的土石混合體強(qiáng)度特征研究[J].巖石力學(xué)與工程學(xué)報(bào)，2009，28（增1）：2651-2656.

SHU Zhile，LIU Xinrong，LIU Baoxian，et al.Study of strength properties of earthrock aggregate based on fractals[J].Chinese Journal of Rock Mechanics and Engineering，2009，28（S1）：2651-2656.

[14]楊忠平，雷曉丹，王雷，等.含石量對(duì)土石混合體剪切特性影響的顆粒離散元數(shù)值研究[J].工程地質(zhì)學(xué)報(bào)，2017，25（4）：1035-1045.

YANG Zhongping，LEI Xiaodan，WANG Lei，et al.Impact of stone content to shear properties of soilrock mixture using particle flow code simulation[J].Journal of Engineering Geology，2017，25（4）：1035-1045.

[15]詹永祥，姚海林.塊石含量對(duì)土石混合體抗剪強(qiáng)度特性的影響[J].公路，2018（2）：20-25.

ZHAN Yongxiang，YAO Hailin.Effect of stone content on shear strength characteristics of soilrock mixture[J].Highway，2018（2）：20-25.

[16]夏加國(guó)，胡瑞林，祁生文，等.含超徑顆粒土石混合體的大型三軸剪切試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào)，2017，36（8）： 2031-2039.

XIA Jiaguo，HU Ruilin，QI Shengwen，et al.Largescale triaxial shear testing of soil rock mixtures containing oversized particles[J].Chinese Journal of Rock Mechanics and Engineering，2017，36（8）：2031-2039.

[17]王強(qiáng)，胡新麗，徐迎，等.軟巖土石混合體的擊實(shí)特性與細(xì)觀機(jī)理研究[J].浙江大學(xué)學(xué)報(bào)（工學(xué)版），2018，52（12）：2995-2305.

WANG Qiang，HU Xinli，XU Ying，et al.Investigation on compaction characteristics and mesomechanismof soilrock mixture with soft rock blocks[J].Journal of Zhejiang University（Engineering Science），2018，52（12）：2995-2305.

[18]ZHANG Haiyang，XU Wenjie，YU Yuzhen.Triaxial tests of soil–rock mixtures with different rock block distributions[J].Soils and Foundations，2016，56（1）：44-56.

[19]SHAN Pengfei，LAI Xingping.Mesoscopic structure PFC2D model of soil rock mixture based on digital image[J].Journal of Visual Communication and Image Representation，2019，58（1）：407-415.

[20]Li X，Liao Q L，He J M.Insitu tests and a stochastic structural model of rock and soil aggregate in the three gorges reservoir area，China[J].Journal of Rock Mechanics and Mining Science，2005，41（3）：494-495.

[21]Xu W J，Hu R L，Tan R J.Some geomechanical properties of soilrock mixtures in Tigerleaping Gorge Area，China[J].Geotechnique，2007，57（3）：255-264.

[22]GONG Jian，LIU Jun.Analysis on the mechanical behaviors of soilrock mixtures using Discrete Element Method[J].Procedia Engineering，2015，102（12）：1783-1792.

[23]XU Wenjie，WANG Shi，ZHANG Haiyang，et al.Discrete element modelling of a soilrock mixture used in an embankment dam[J].International Journal of Rock Mechanics & Mining Sciences，2016，86（7）：141-156.

[24]WEI Houzhen，XU Wenjie，WEI Changfu，et al.Influence of water content and shear rate on the mechanical behavior of soilrock mixtures[J].Technological Sciences，2018，61（8）：1127-1136.

[25]遲明杰，趙成剛，李小軍.砂土剪脹機(jī)理的研究[J].土木工程學(xué)報(bào)，2009，42（3）：99-104.

CHI Mingjie，ZHAO Chenggang，LI Xiaojun.Stressdilation mechanism of sands[J].China Civil Engineering Journal，2009，42（3）：99-104.