基于級配方程的粗粒料級配演化預(yù)測模型

郭萬里，朱俊高，王青龍，余挺

?

基于級配方程的粗粒料級配演化預(yù)測模型

郭萬里1，朱俊高1，王青龍2，余挺2

(1. 河海大學(xué) 巖土力學(xué)與堤壩工程教育部重點(diǎn)實(shí)驗(yàn)室，江蘇 南京，210098；2. 中國電力建設(shè)集團(tuán)成都勘測設(shè)計研究院有限公司，四川 成都，610072)

為了揭示粗粒料級配隨應(yīng)力應(yīng)變演化的規(guī)律，建立粗粒料“應(yīng)力應(yīng)變→破碎指標(biāo)→級配分布”的數(shù)學(xué)模型(SBG模型)。首先，引入1個適用性良好的級配方程來描述土體的級配，其級配參數(shù)為和；同時，定義1個新的破碎指標(biāo)W，并將W和破碎指標(biāo)g建議為該模型所需的2個破碎指標(biāo)。其次，推導(dǎo)并且驗(yàn)證W和g與級配參數(shù)和的數(shù)學(xué)關(guān)系，實(shí)現(xiàn)“破碎指標(biāo)→級配分布”這一目標(biāo)。根據(jù)粗粒料的三軸試驗(yàn)結(jié)果總結(jié)三軸試樣破壞時W和g與圍壓3之間的關(guān)系、三軸試樣在剪切過程中W和g與平均正應(yīng)力和廣義剪應(yīng)變s之間的關(guān)系，實(shí)現(xiàn)“應(yīng)力應(yīng)變→破碎指標(biāo)”這一目標(biāo)。將以上2個部分聯(lián)合，即可得到“應(yīng)力應(yīng)變→破碎指標(biāo)→級配分布”的數(shù)學(xué)模型(SBG模型)。研究結(jié)果表明：本文提出的模型可以用來描述三軸應(yīng)力狀態(tài)下粗粒料級配隨應(yīng)力應(yīng)變演化的規(guī)律。

粗粒料；顆粒破碎；級配方程；級配演化；三軸試驗(yàn)

粗粒料在加載過程中會發(fā)生顆粒破碎，從而顯著影響土體的強(qiáng)度、變形、滲透性等性質(zhì)[1?3]，因此，對于粗粒料顆粒破碎的研究一直以來都是土力學(xué)的熱 點(diǎn)[4?9]。目前研究多集中于顆粒破碎指標(biāo)。張季如等[3]分析了相對破碎指標(biāo)r與應(yīng)力水平之間的關(guān)系，賈宇峰等[5]通過三軸試驗(yàn)總結(jié)了剪切過程中r隨剪應(yīng)變的變化規(guī)律。不少研究者基于分形理論提出了顆粒破碎之后的分形維數(shù)與圍壓之間的經(jīng)驗(yàn)公式[1?2]。此外，也有研究者嘗試將顆粒破碎規(guī)律融入到本構(gòu)模型之中，但是現(xiàn)有本構(gòu)模型大多是在宏觀層面上反映顆粒破碎對粗粒料應(yīng)力變形的影響[4]?？偟膩碚f，目前研究大多能夠建立破碎指標(biāo)與應(yīng)力或應(yīng)變之間的關(guān)系[5?6]，但是無法進(jìn)一步反映級配的演化。值得注意的是，“破碎指標(biāo)”與“級配分布”并不是呈一一對應(yīng)的關(guān)系：由級配分布可以計算破碎指標(biāo)，但是由破碎指標(biāo)無法推算出級配分布。雖然破碎指標(biāo)在一定程度上能夠反映土體性質(zhì)的變化，但是對于某些土體僅能預(yù)測破碎指標(biāo)是不夠的。例如，在大壩填筑過程中，反濾土石壩中的反濾料的受力是不斷增加的，所引發(fā)的顆粒破碎會導(dǎo)致其級配的變化，而級配又能顯著影響其滲透性。若能夠建立反映粗粒料的級配隨應(yīng)力應(yīng)變演化的數(shù)學(xué)模型，則這對于土石壩等工程的設(shè)計、監(jiān)測與維護(hù)都具有重要意義。為此，本文作者引入對各種級配曲線描述效果良好的級配方程[10?12]，并且定義1個新的顆粒破碎指標(biāo)，推導(dǎo)破碎指標(biāo)與級配參數(shù)之間的關(guān)系，并依此建立1個反映顆粒級配隨應(yīng)力應(yīng)變演化的數(shù)學(xué)模型。最后，以粗粒料的三軸試驗(yàn)結(jié)果為例，驗(yàn)證該模型的合理性。

1 模型的建立

1.1 級配的定量表示方法

本文作者的最終目標(biāo)是建立“應(yīng)力應(yīng)變→級配分布”的數(shù)學(xué)模型，即通過土體的應(yīng)力應(yīng)變狀態(tài)可以預(yù)測當(dāng)前的級配分布。但是，現(xiàn)有的研究大都止步于“應(yīng)力應(yīng)變→破碎指標(biāo)”，即通過應(yīng)力應(yīng)變可以預(yù)測出顆粒破碎的程度，無法進(jìn)一步推測出當(dāng)前的級配分布。

分析其原因，土體的應(yīng)力和應(yīng)變是用定量的數(shù)值表示的，而目前級配分布都只是定性描述，而不是定量的數(shù)學(xué)表示(主要有級配曲線和粒組質(zhì)量分?jǐn)?shù)2種表示方法，雖然各個粒組的質(zhì)量分?jǐn)?shù)都有具體的數(shù)值，但是對于級配分布整體而言，依然是定性描述)。因而“定量的”應(yīng)力應(yīng)變與“定性的”級配之間無法建立數(shù)學(xué)關(guān)系。

因此，本文作者試圖引入之前提出的一個適用性良好的級配方程[10?12]，將級配曲線量化表示為關(guān)于參數(shù)和的函數(shù)。當(dāng)已知級配分布時，一般采用最優(yōu)化擬合[12]來確定參數(shù)和；若已知參數(shù)和，則可以代入式(1)繪出級配曲線。因此，在本文中級配方程與級配分布是等價的。該級配方程為

式中：為任意粒徑；max為最大粒徑，通常為已知量；為粒徑小于的顆粒所占百分?jǐn)?shù)；和為參數(shù)，其中＜1且＞0，以下稱為級配參數(shù)。

1.2 模型的統(tǒng)一表達(dá)式

通過粗粒料在三軸應(yīng)力狀態(tài)下的顆粒破碎試驗(yàn)，初步發(fā)現(xiàn)級配參數(shù)和與圍壓3之間并不存在規(guī)律性關(guān)系，因此，直接建立“應(yīng)力應(yīng)變→級配參數(shù)(級配分布)”的關(guān)系是比較困難的。

引入破碎指標(biāo)來作為應(yīng)力應(yīng)變和級配參數(shù)之間的過渡條件，即建立“應(yīng)力應(yīng)變→破碎指標(biāo)→級配參數(shù)”的關(guān)系。一方面，“應(yīng)力應(yīng)變→破碎指標(biāo)”已有較多的研究成果可供參考[1?2]；另一方面，“破碎指標(biāo)→級配參數(shù)”的關(guān)系可以通過破碎指標(biāo)的定義進(jìn)行數(shù)學(xué)推導(dǎo)。因此，建立這樣的數(shù)學(xué)模型在理論上是切實(shí)可行的。

該數(shù)學(xué)模型的本質(zhì)是通過建立方程組求解當(dāng)前的級配參數(shù)和。由于級配方程有和這2個參數(shù)作為未知數(shù)，因此，需要聯(lián)立2個不同的方程來求解。基于此，引入2個破碎指標(biāo)W和g。

首先，根據(jù)破碎指標(biāo)的定義，只需要進(jìn)行數(shù)學(xué)上的推導(dǎo)即可確定“破碎指標(biāo)→級配參數(shù)”的函數(shù)關(guān)系，設(shè)其統(tǒng)一表達(dá)式為

其次，“應(yīng)力應(yīng)變→破碎指標(biāo)”可以通過試驗(yàn)總結(jié)理論或經(jīng)驗(yàn)公式，設(shè)其統(tǒng)一表達(dá)式為

式中：為應(yīng)力；為應(yīng)變。

最后，聯(lián)立式(2)和式(3)即為本文所建立的“應(yīng)力應(yīng)變→破碎指標(biāo)→級配參數(shù)”的數(shù)學(xué)模型。當(dāng)式(2)和式(3)中的函數(shù)和確定之后，根據(jù)土體所處的應(yīng)力應(yīng)變狀態(tài)即可求出當(dāng)前的級配參數(shù)和，即確定新的級配分布。

1.3 破碎指標(biāo)的選擇

由于式(2)和式(3)都是關(guān)于破碎指標(biāo)的函數(shù)，因此，選擇合適的破碎指標(biāo)W和g成為模型建立的重要環(huán)節(jié)。目前，顆粒破碎指標(biāo)總的來說可以分為3類：第1類是基于特定粒徑或系數(shù)，第2類是基于粒組質(zhì)量分?jǐn)?shù)，第3類是基于破碎勢。

第1類用單個粒徑(比如10和15等)或系數(shù)(比如u或c)的變化來衡量顆粒破碎程度[13]，在計算方面簡單實(shí)用，但描述的級配特征比較單一，難免以偏 概全。

第2類用粒組質(zhì)量分?jǐn)?shù)的變化來衡量顆粒破碎程度，以Marsal[7]定義的g為代表。若用0i和w分別表示試驗(yàn)前、后第個粒組的質(zhì)量分?jǐn)?shù)，則該粒組的質(zhì)量分?jǐn)?shù)變化量為(w?0i)，其中，0i和w示意圖如圖1所示。Marsal定義的破碎指標(biāo)g為所有(w?0i)的正值之和[7]。g能夠反映土顆粒整體破碎的特點(diǎn)，且簡單實(shí)用，因此，將作為模型所需的破碎指標(biāo)之一。

圖1 w0i和wi示意圖

第3類以HARDIN等[8?9]定義的破碎率為代表，是以顆粒破碎前后的級配曲線所圍成的面積作為破碎量，再除以各自所定義的破碎潛能，得到的比值即為破碎率。其中，HARDIN[8]將粒徑=0.074 mm與初始級配曲線圍成的面積作為破碎潛能；EINAV[9]則認(rèn)為顆粒破碎會存在1個極限級配，因此，將=0.074 mm的限制取消，改為將初始級配曲線與極限級配曲線圍成的面積作為破碎潛能，在理論上更為合理。將初始級配曲線、試驗(yàn)后的級配曲線和極限級配曲線與最大粒徑線=max所圍成的面積分別表示為0，1和2，則EINAV[9]定義的破碎指標(biāo)可表示為

事實(shí)上，E在實(shí)際運(yùn)用中并不方便：相比于其他破碎指標(biāo)，確定極限級配曲線的試驗(yàn)屬于額外試驗(yàn)，加大了試驗(yàn)量，而且該試驗(yàn)需要在高應(yīng)力狀態(tài)下進(jìn)行，難度較大。因此，E雖然在理論上為人們所認(rèn)可，但目前被作為破碎指標(biāo)應(yīng)用于實(shí)際研究的不多。由式(4)可知：對于給定的某種土料，初始級配曲線確定之后(即0為定值)，對應(yīng)的極限級配曲線也隨之確定(即2也為定值)。因此，EINAV破碎指標(biāo)E的本質(zhì)是用破碎量(1?0)除以一個同量綱的定值，得到量綱一的百分比即為破碎指標(biāo)。基于此，將式(4)的分母由定值(2?0)替換為定值0，從而定義1個新的破碎指標(biāo)W，如圖2所示。W表達(dá)式為

圖2 BW示意圖

在理論的完整性上W不如E：E的變化區(qū)間是0~1，適合用來對比評價不同土體的顆粒破碎程度；W的理論變化區(qū)間則是0~無窮大。但是，W反映的是級配整體的變化，且簡單實(shí)用，因此，將其作為模型所需的另外1個破碎指標(biāo)。

至此，模型要求的2個破碎指標(biāo)已選定，即g和W，其優(yōu)勢是都能夠反映級配整體的變化，且定義簡單，計算方便。

2 破碎指標(biāo)與級配參數(shù)的關(guān)系推導(dǎo)

2.1 BW與級配參數(shù)的關(guān)系

HARDIN[8]和EINAV[9]定義的破碎指標(biāo)在實(shí)際應(yīng)用時，一般通過梯形分割法求取級配曲線所圍面積。若利用本文的級配方程，則可以直接利用方程積分求取面積。級配方程計算面積示意圖如圖3所示。

由于粒徑的最小值不能為0 mm，故式(1)中的最小值只能趨向于0，最大值為100%。因此，在利用積分求面積時可以先對在區(qū)間[，1]內(nèi)進(jìn)行積分，然后使趨向于0并求極限(見圖3)。

圖3 級配方程計算面積示意圖

級配曲線與=max所圍的面積可表示為

式中d為級配曲線上=時所對應(yīng)的粒徑。

結(jié)合式(6)和式(1)可得：

當(dāng)趨向于0時，將式(7)代入式(6)可得

特別地，當(dāng)=0時，式(8)可進(jìn)一步表示為

因此，若初始級配參數(shù)為0和0，顆粒破碎后的級配參數(shù)為和，根據(jù)式(8)和式(5)可得到本文定義的破碎指標(biāo)W：

當(dāng)趨向于0時對應(yīng)的粒徑為0 mm，而實(shí)際上土顆粒的最小粒徑不可能為0 mm，因此，將設(shè)定為一個較小的值，但不為0，這樣更符合級配分布的實(shí)際情況。

為約定的取值，此處以某砂板巖的三軸試驗(yàn)為例[14]，試驗(yàn)所得各圍壓下的級配參數(shù)如表1所示(其中利用式(10)所求得的面積和破碎指標(biāo)W也列于表1)。當(dāng)分別取為1.0%，0.5%和0.1%時，相同圍壓下的面積和破碎指標(biāo)W基本保持不變。由此可見，當(dāng)較小時，利用式(10)計算得到的破碎指標(biāo)W是穩(wěn)定的，不隨的變化而發(fā)生顯著變化。因此，可以將約定為0.1%。

綜上所述，利用本文的級配方程積分求取破碎指標(biāo)W是可行的。同理，該方法也可以推廣到計算HARDIN[8]和EINAV[9]定義的破碎指標(biāo)。

2.2 Bg與級配參數(shù)的關(guān)系

根據(jù)g的定義，由于粒組質(zhì)量分?jǐn)?shù)的增加量與其余粒組質(zhì)量分?jǐn)?shù)的減少量是相等的，因此，g又可以表示為各粒組質(zhì)量分?jǐn)?shù)變化值之和的1/2，即

0i和w利用級配參數(shù)可以分別表示為

式中：函數(shù)為級配方程(見式(1))；d和d?1分別為某個粒組的最大和最小粒徑。在GB T50123—1999“土工試驗(yàn)方法標(biāo)準(zhǔn)”[15]中，顆粒分析試驗(yàn)分別采用孔徑為60.00，40.00，20.00，10.00，5.00，2.00，1.00，0.50，0.25和0.075 mm的土工篩。

式(11)和式(12)分別確定了g與級配參數(shù)和的函數(shù)關(guān)系。值得注意的是，級配曲線是通過篩分試驗(yàn)確定的，而級配方程是對篩分所得級配曲線進(jìn)行數(shù)學(xué)擬合，兩者之間存在擬合誤差。圖4所示為級配方程與篩分試驗(yàn)所得w示意圖(其中，公式w為利用式(12)計算得到的各個粒組的質(zhì)量分?jǐn)?shù)；篩分w為篩分試驗(yàn)得到的粒組質(zhì)量分?jǐn)?shù)。二者之間存在誤差，因此，利用式(11)所計算的破碎指標(biāo)g與篩分試驗(yàn)確定的g也存在誤差。

表1 面積S，破碎指標(biāo)BW與k的關(guān)系

表2 級配方程計算Bg的誤差分析

圖4 級配方程與篩分試驗(yàn)所得wi示意圖

朱俊高等[10]的研究表明式(1)對于級配曲線的擬合相關(guān)系數(shù)2一般都在0.95以上，由此可知式(12)計算得到的各個粒組質(zhì)量分?jǐn)?shù)與篩分試驗(yàn)相比雖然存在誤差，但是誤差應(yīng)該在合理范圍之內(nèi)。以某堆石料的三軸試樣顆粒破碎數(shù)據(jù)為例，通過篩分試驗(yàn)得到各粒組質(zhì)量分?jǐn)?shù)，通過級配方程擬合得到參數(shù)和。同時利用篩分試驗(yàn)數(shù)據(jù)和式(11)對各圍壓下的g進(jìn)行計算。進(jìn)一步地，以篩分試驗(yàn)計算的g為基準(zhǔn)，算出式(11)計算結(jié)果的相對誤差。級配方程計算g的誤差分析如表2所示。由表2可知：兩者相對誤差都在10%以內(nèi)，屬于合理范圍。因此，可以認(rèn)為利用式(11)計算破碎指標(biāo)g是合理的。

綜上所述，模型中的2個破碎指標(biāo)W和g與級配參數(shù)和之間的函數(shù)關(guān)系式已經(jīng)確定，即模型中的式(2)已確定。

3 破碎指標(biāo)與應(yīng)力應(yīng)變的關(guān)系

3.1 三軸試樣破壞時的顆粒破碎規(guī)律

目前有關(guān)顆粒破碎的研究主要是基于三軸試樣破壞時的試驗(yàn)結(jié)果，因此，下面以普通三軸試驗(yàn)為例，總結(jié)三軸應(yīng)力狀態(tài)下式(3)的函數(shù)關(guān)系1和2，初步驗(yàn)證本模型的合理性。

實(shí)際上，已有較多的經(jīng)驗(yàn)公式可供參考[1, 3, 5]，其中大部分是描述三軸試驗(yàn)破碎時的破碎指標(biāo)與應(yīng)力或應(yīng)變的單因素關(guān)系。對于三軸試樣破壞時的破碎規(guī)律，本文作者建議了破碎指標(biāo)與圍壓3之間的經(jīng)驗(yàn)公式為

式中：1，1，2和2為擬合參數(shù)。

對于式(13)的適用性，現(xiàn)以1組三軸試樣破壞時的數(shù)據(jù)來檢驗(yàn)其擬合效果。三軸試樣破壞時破碎指標(biāo)與圍壓的關(guān)系如圖5所示(其中，3為圍壓；1=9.16，1=0.465；2=5.81，2=0.349)。從圖5可以看出：式(13)對于破碎指標(biāo)與圍壓的關(guān)系描述效果較好。

3.2 三軸試樣剪切過程中的顆粒破碎規(guī)律

式(13)只是針對試樣在不同圍壓下破壞時的這個“點(diǎn)”，而無法推廣到整個加載剪切的“過程”。三軸剪切試驗(yàn)過程中目前也有一些經(jīng)驗(yàn)公式可供參考。比如，賈宇峰等[5]將r表示為廣義剪應(yīng)變s的函數(shù)；劉恩龍等[4]在壓縮試驗(yàn)中將破碎指標(biāo)表示為平均正應(yīng)力的函數(shù)；SALIM等[6]在其彈塑性模型中將破碎指標(biāo)g表示為塑性剪應(yīng)變sp、臨界狀態(tài)應(yīng)力cs(i)和的函數(shù)。本文作者則建議了1個將破碎指標(biāo)表示為廣義剪應(yīng)變s和平均正應(yīng)力的經(jīng)驗(yàn)公式：

圖5 三軸試樣破壞時破碎指標(biāo)與圍壓的關(guān)系

式中：s為粗粒土的固相硬度，一般通過壓縮試驗(yàn)確定[16]。1，1，2和2均為擬合參數(shù)。

對于式(14)的適用性，現(xiàn)以1組三軸剪切過程中的顆粒破碎數(shù)據(jù)[5]來檢驗(yàn)其擬合效果。圖6所示為用式(14)擬合的三軸剪切過程中破碎指標(biāo)與應(yīng)力應(yīng)變的關(guān)系(其中，為破碎率；s=19.1 MPa，1=0.61，1=32，2=0.265，2=39.1)。由圖6可知：式(14)能夠較好地描述破碎指標(biāo)與應(yīng)力應(yīng)變之間的關(guān)系。

土的應(yīng)力狀態(tài)有多種，比如普通三軸應(yīng)力狀態(tài)、0應(yīng)力狀態(tài)以及其他多種復(fù)雜應(yīng)力狀態(tài)等，因此，總結(jié)出適用于一般應(yīng)力狀態(tài)的函數(shù)關(guān)系1和2是本文的最終目的。式(13)和式(14)則都是基于三軸應(yīng)力狀態(tài)，其中式(13)描述的是三軸試樣破壞時粗粒料的破碎規(guī)律，式(14)則描述的是三軸剪切過程中粗粒料的破碎規(guī)律。

圖6 三軸剪切過程中破碎指標(biāo)與應(yīng)力應(yīng)變的關(guān)系

4 模型初步驗(yàn)證

模型的適用性主要取決于2個環(huán)節(jié)：首先是“應(yīng)力應(yīng)變→破碎指標(biāo)”之間經(jīng)驗(yàn)公式的適用性。這一部分是通過顆粒破碎試驗(yàn)結(jié)果來進(jìn)行總結(jié)的。雖然目前尚無適用于一般應(yīng)力狀態(tài)下的經(jīng)驗(yàn)公式，但適用于某一類應(yīng)力狀態(tài)的經(jīng)驗(yàn)公式有較多文獻(xiàn)可供參考[1, 3, 5]。因此，這一部分雖然是模型的重點(diǎn)，但并不是模型的難點(diǎn)。其次是“破碎指標(biāo)→級配參數(shù)”是否可行，這是模型是否成立的關(guān)鍵。

下面主要證明已知W和g時，是否可以求得和。以某粗粒料的三軸試驗(yàn)為例，當(dāng)圍壓分別為 0.6 MPa和1.5 MPa時，試樣篩分實(shí)驗(yàn)破碎指標(biāo)W1=22.2%，g1=10.1%(圍壓為0.6 MPa)；W2=34.3%，g2=15.3% (圍壓為1.5 MPa)。將W和g試驗(yàn)值分別代入式(10)和式(11)聯(lián)立方程求解，得到級配參數(shù)分別如下：1=0.607和1=0.891(圍壓為0.6 MPa)；2=0.638和2=0.794(圍壓為1.5 MPa)。將和代入式(1)，繪制得到的級配曲線預(yù)測值如圖7所示。從圖7可見：模型預(yù)測的級配曲線與篩分試驗(yàn)結(jié)果較吻合，這說明了模型中“破碎指標(biāo)→級配參數(shù)”是可行的。

圖7 模型對級配分布的預(yù)測

5 結(jié)論

1) 定義1個新的破碎指標(biāo)W，連同MARSAL定義的破碎指標(biāo)g一起建議為該模型中所需的2個破碎指標(biāo)，同時推導(dǎo)并證明了W和g與級配參數(shù)和之間的數(shù)學(xué)關(guān)系式，即完成模型中“破碎指標(biāo)→級配參數(shù)”的轉(zhuǎn)換。

2) 對于三軸試樣破壞時顆粒破碎指標(biāo)與圍壓的關(guān)系、三軸試樣剪切過程中破碎指標(biāo)與應(yīng)力和應(yīng)變的關(guān)系，建議合適的經(jīng)驗(yàn)公式，并利用試驗(yàn)數(shù)據(jù)初步證明其適用性，即完成模型中“應(yīng)力應(yīng)變→破碎指標(biāo)”的轉(zhuǎn)換。

3) 以三軸試驗(yàn)為例，已知試驗(yàn)所得的破碎指標(biāo)W和g，成功預(yù)測三軸試樣在不同圍壓下破壞時的級配分布，初步驗(yàn)證了模型的合理性。

[1] 吳愛祥, 王建棟, 彭乃兵. 顆粒級配對粗骨料充填料漿離析的影響[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2016, 47(9): 3201?3207. WU Aixiang, WANG Jiandong, PENG Naibing. Effect of grain composition on coarse aggregate filling slurry segregation[J]. Journal of Central South University(Science and Technology), 2016, 47(9): 3201?7.

[2] 陳镠芬, 高莊平, 朱俊高, 等. 粗粒土級配及顆粒破碎分形特性[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2015(9): 3446?3453. CHEN L, GAO Zhuangping, ZHU Jungao, et al. Gradation of coarse grained soil and fractal geometry character of particle breakage[J]. Journal of Central South University(Science and Technology), 2015(9): 3446?3453.

[3] 張季如, 張弼文, 胡泳, 等. 粒狀巖土材料顆粒破碎演化規(guī)律的模型預(yù)測研究[J]. 巖石力學(xué)與工程學(xué)報, 2016, 35(9): 1898?1905. ZHANG Jiru, ZHANG Biwen, HU Yong, et al. Predicting the particle breakage of granular geomaterials[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(9): 1898?1905.

[4] 劉恩龍, 陳生水, 李國英, 等. 堆石料的臨界狀態(tài)與考慮顆粒破碎的本構(gòu)模型[J]. 巖土力學(xué), 2011, 32(S2): 148?154. LIU Enlong, CHEN Shengshui, LI Guoying, et al. Critical state of rockfill materials and a constitutive model considering grain crushing[J]. Rock & Soil Mechanics, 2011, 32(S2): 148?154.

[5] 賈宇峰, 王丙申, 遲世春. 堆石料剪切過程中的顆粒破碎研究[J]. 巖土工程學(xué)報, 2015, 37(9): 1692?1697. JIA Yufeng, WANG Bingshen, CHI Shichun. Particle breakage of rockfill during triaxial tests[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(9): 1692?1697.

[6] Salim, W., Indraratna, B. A new elastoplastic constitutive model for coarse granular aggregates incorporating particle breakage[J]. Canadian Geotechnical Journal, 2004, 41(4), 657?671.

[7] MARSAL R J. Large-scale testing of rockfill materials[J]. Journal of the Soil Mechanics and Foundation Division, 1967, 93(2): 27?44.

[8] HARDIN B O. Crushing of soil particles[J]. Journal of Geotechnical Engineering, 1985, 111(10): 1177?1192.

[9] EINAV I. Breakage mechanics-part I: theory[J]. Journal of the Mechanics and Physics of Solids, 2007, 55(6): 1274?1297.

[10] 朱俊高, 郭萬里, 王元龍, 等. 連續(xù)級配土的級配方程及其適用性研究[J]. 巖土工程學(xué)報, 2015, 37(10): 1931?1936. ZHU Jungao, GUO Wanli, WANG Yuanlong, et al. The research on equation of soil gradation curve and its applicability[J]. Chinese Journal of Geotechnical, 2015, 37(10): 1931?1936.

[11] 郭萬里, 朱俊高, 溫彥鋒. 對粗粒料4種級配縮尺方法的解釋[J]. 巖土工程學(xué)報, 2016, 38(8): 1473?1480.GUO Wanli, ZHU Jungao, WEN Yanfeng. Unified description for four grading scale methods for coarse aggregate[J]. Chinese Journal of Geotechnical, 2016, 38(8): 1473?1480.

[12] ZHU Jungao, GUO Wanli, WEN Yanfeng, et al. A new gradation equation and its applicability for particle size distributions of various soils[J]. International Journal of Geomechanics, 2018, 18(2): 1?8.

[13] 柏樹田, 崔亦昊. 堆石的力學(xué)性質(zhì)[J]. 水力發(fā)電學(xué)報, 1997, 16(3): 21?30. BAI Shutian, CUI Yihao. The mechanical properties of rockfill[J]. Journal of Hydroelectric Engineering, 1997, 16(3): 21?30.

[14] 傅華, 凌華, 蔡正銀. 粗顆粒土顆粒破碎影響因素試驗(yàn)研究[J]. 河海大學(xué)學(xué)報(自然科學(xué)版), 2009, 37(1): 75?79. FU Hua, LIN Hua, CAI Zhengyin. Influencing factors for particle breakage of coarse grained soil[J]. Journal of Hohai University(Natural Sciences), 2009, 37(1): 75?79.

[15] GB 50123—1999, 土工試驗(yàn)方法標(biāo)準(zhǔn)[S]. GB 50123—1999, Standard for geotechnical test methods[S].

[16] 陳生水, 傅中志, 韓華強(qiáng),等. 一個考慮顆粒破碎的堆石料彈塑性本構(gòu)模型[J]. 巖土工程學(xué)報, 2011, 33(10): 1489?1495. CHEN Shengshui, FU Zhongzhi, HAN H Q, et al. An elastoplastic model for rockfill materials considering particle breakage[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(10): 1489?1495.

(編輯 伍錦花)

Mathematical model based on the gradation equation for predicting gradation evolution of coarse-grained soils

GUO Wanli1, ZHU Jungao1, WANG Qinglong2, YU Ting2

(1. Key Laboratory of Geomechanics and Embankment Engineering of Ministry of Education, Hohai University, Nanjing 210098, China; 2. Chengdu Survey and Design Institute Limited, Power Construction Corporation of China, Chengdu 610072, China)

In order to study the evolution law of gradation for coarse-grained soils, a mathematical model describing the relationship between stress and strain state, breakage indicators and gradation(SBG model) was established. Firstly, a gradation equation in which the two gradation parameters were namedandwas proposed to describe the grain size distribution(GSD) curve. Then, a new grain breakage indicatorWwas defined. Together with another breakage indicatorg, they were both suggested as the two required breakage indicators for the proposed model. On one hand, the relationship betweenWandgand the gradation parametersandwere derived and verified. On the other hand, an empirical equation ofWandgand confining pressure at the failure of the specimens, and an empirical equation ofWandgand the mean normal stress and general shear strain during shearing of the specimens were proposed. The SBG model was thereby established by combining the above two parts. The results show that the proposed SBG model can be used to predict change of breakage indicators and GSD of coarse-grained soils with stress and strain states in triaxial test.

coarse-grained soil; particle crushing; gradation equation; gradation evolution; triaxial test

10.11817/j.issn.1672?7207.2018.08.030

TU43

A

1672?7207(2018)08?2076?07

2017?08?01；

2017?09?21

國家重點(diǎn)研發(fā)計劃項目(2017YFC0404801)；國家自然科學(xué)基金資助項目(51479052)(Project(2017YFC0404801) supported by the National Key Research and Development Program of China; Project(51479052) supported by the National Natural Science Foundation of China

朱俊高，教授，博士生導(dǎo)師，從事土體基本性質(zhì)及本構(gòu)關(guān)系、土石壩工程研究；E-mail: zhujungao@hhu.edu.cn