薛明華，陳圣仟 ，楊世文

(1.濟(jì)南市勘察測繪研究院，濟(jì)南 250101；2.貴州省建筑設(shè)計(jì)研究院有限責(zé)任公司，貴陽 550000)

砂性土(無黏聚力)是由大量土顆粒構(gòu)成的多體相互作用體系[1]。球形度(form or sphericity)、圓滑度/棱角度(roundness or angularity)、粗糙度(roughness)從3個(gè)不同尺度描述了顆粒的形態(tài)特征[2-3]，被認(rèn)為是顆粒形狀量化的基本指標(biāo)[4]，用于評價(jià)顆粒的形狀。另外，還衍生出諸如凹凸度(non-convexity)、長寬比(aspect ratio)等[5]另外一些顆粒形狀評價(jià)指標(biāo)。顆粒形狀會對砂性土的宏觀性質(zhì)產(chǎn)生影響，如相對密實(shí)度、滲透特性、壓縮特性、剪切強(qiáng)度及變形等[2-10]。而宏觀上展現(xiàn)的力學(xué)特性本質(zhì)上是細(xì)觀尺度(顆粒層面)顆粒間相互作用的累積[4,11]。對于顆粒形狀與砂性土宏觀力學(xué)特性之間的關(guān)系，研究人員開展過大量的室內(nèi)試驗(yàn)或數(shù)值模擬研究工作。試驗(yàn)研究認(rèn)為，幾乎所有砂性土的宏觀力學(xué)特性都會受到顆粒形狀的影響[2-4,11]。比如，顆粒的形狀會影響砂性土的極限孔隙比(相對密實(shí)度)，多棱角、表面粗糙和球形度低的顆粒堆積體具有較大的極限孔隙比區(qū)間長度，相反，較光滑球形顆粒堆積體的極限孔隙比區(qū)間長度較小[2,12-14]。在三軸壓縮試驗(yàn)、直剪試驗(yàn)中，相同密實(shí)度條件下，多棱角的砂性土比表面圓滑的砂性土具有較大的強(qiáng)度，且在峰值應(yīng)力前具有較大變形模量[9,13]。當(dāng)試樣體積變形恒定時(shí)，表面圓滑的砂較表面多棱角的砂穩(wěn)定性差，強(qiáng)度低[15]。臨界狀態(tài)下砂性土的內(nèi)摩擦角、臨界狀態(tài)線的斜率、峰值內(nèi)摩擦角以及峰值應(yīng)力前砂性土的變形模量均隨砂粒棱角度的增加而增加[9,16]。

多數(shù)室內(nèi)或數(shù)值試驗(yàn)的研究闡明了不同顆粒形狀會影響砂性土的強(qiáng)度和變形特性，但顆粒形狀對砂性土強(qiáng)度和變形特性的影響機(jī)理則鮮有研究。室內(nèi)試驗(yàn)具有明顯的局限性，如顆粒的粗糙度、棱角度和球形度對砂性土宏觀力學(xué)性質(zhì)的影響很難做到實(shí)際意義上的分離[4,9,17]。制備由規(guī)則形狀的顆粒構(gòu)成的試樣無疑脫離了天然砂性土中顆粒隨機(jī)不規(guī)則形狀的特性，且放大了顆粒形狀的影響。在這種情況下，諸多學(xué)者對顆粒形狀與砂性土力學(xué)特性的關(guān)系持不同意見，除了顆粒的球形度、圓滑度和粗糙度與砂性土力學(xué)特性的基本認(rèn)識以外，有人認(rèn)為顆粒形狀對砂性土剪切強(qiáng)度的影響甚微[18]，有人指出顆粒形狀是否對砂性土剪切強(qiáng)度產(chǎn)生影響尚不清楚[2]。

在此基礎(chǔ)上，現(xiàn)通過室內(nèi)直剪試驗(yàn)及其離散元模擬試驗(yàn)探索顆粒形狀與砂性土力學(xué)特性，特別是強(qiáng)度與變形特性。研究對象是較為常見的砂性土，其顆粒為隨機(jī)不規(guī)則形狀，且在相同的試驗(yàn)條件下設(shè)計(jì)基于圓形顆粒的對比試驗(yàn)。

1 直剪試驗(yàn)

直剪試驗(yàn)的試樣為天然硅質(zhì)粗砂，其級配曲線如圖1所示，不均勻系數(shù)為1.44。試樣平均粒徑為0.9 mm，密度為2.5 g/cm3，最大孔隙比1.076，最小孔隙比0.67，L/d50≈44>40且L/dmax≈18>15(L為試樣的最小尺寸，d50為平均顆粒粒徑，dmax為最大顆粒粒徑)。根據(jù)文獻(xiàn)[19-20]，試樣的尺寸效應(yīng)很小。

剪切盒總尺寸為60 mm×60 mm×43 mm(長×寬×高)，包括上下兩盒。剪切過程采用應(yīng)變控制式，剪切速率可在0.001～2 mm/min范圍內(nèi)自由調(diào)節(jié)，上盒固定，由電機(jī)推動下盒剪切，垂直壓力通過杠桿懸掛載荷砝碼施加。試驗(yàn)時(shí)，垂直和水平位移分別通過兩個(gè)千分位移表讀取，剪切力通過水平應(yīng)力環(huán)采集，精度達(dá)±2.2 N。

制樣時(shí)按照50%的相對密實(shí)度(中密)換算試樣的質(zhì)量，分層擊實(shí)，最后平整試樣的上表面，以較慢的剪速(約0.05 mm/min)在3種法向應(yīng)力(22、50和100 kPa)下剪切。對應(yīng)每種法向應(yīng)力進(jìn)行3個(gè)試樣剪切，當(dāng)2次試驗(yàn)結(jié)果差異很大時(shí)增加剪切試樣進(jìn)行剪切試驗(yàn)，每種法向應(yīng)力下試驗(yàn)結(jié)果取3次試驗(yàn)的算術(shù)平均值。

2 模擬試驗(yàn)

2.1 粒徑不同、細(xì)觀結(jié)構(gòu)相同試樣生成方法

采用二維顆粒流離散元分析程序(particle flow code，PFC)模擬直剪試驗(yàn)。模擬直剪試驗(yàn)的理論基礎(chǔ)參見文獻(xiàn)[20]。模擬直剪試驗(yàn)的對象可以是具有不同長寬比的橢圓形顆粒[21]、邊數(shù)可控制的多邊形顆粒[22]，還可以是通過顆粒自動填充(computerized particle analysis,CPA)形成的任意形狀的顆粒[23]。由于天然顆粒形狀的復(fù)雜性，要實(shí)現(xiàn)模擬試樣中顆粒的形狀與實(shí)際的完全相同，現(xiàn)有的技術(shù)是做不到的。然而，即使在顆粒形狀上存在差異，砂性土在直剪試驗(yàn)中表現(xiàn)的力學(xué)特性還是可以被模擬試驗(yàn)準(zhǔn)確的再現(xiàn)，以至于追求顆粒形狀的完美匹配顯得沒有必要[7]。

劉清秉等[24]、Ye[25]通過比較樣本容量為25、50和100的顆粒的形狀統(tǒng)計(jì)指標(biāo)，指出砂性土的某些顆粒形狀量化指標(biāo)受樣本容量影響較小。據(jù)此，從直剪試驗(yàn)的試樣中任意選取了25個(gè)砂粒(圖2)，基于顆粒的數(shù)字圖像處理技術(shù)[5]，獲取顆粒的輪廓后，采用CPA填充技術(shù)在PFC中生成目標(biāo)形狀顆粒(圖3)。隨后按照直剪試驗(yàn)試樣的尺寸和級配曲線(圖1)生成不規(guī)則形狀的試樣。為進(jìn)行對比性分析，采用等體積替換法[25]，將不規(guī)則的顆粒等體積的替換為圓形顆粒。為方便表述，稱圓形顆粒試樣為試樣一，不規(guī)則形狀顆粒試樣為試樣二。由圖1可看出替換后試樣的級配幾乎沒變。

圖2 電鏡下的砂粒Fig.2 Sand under the electron microscope

圖3 CPA填充技術(shù)生成的顆粒Fig.3 Particles generated by CPA filling technology

顆粒在沒有摩擦作用時(shí)可以堆積成致密的試樣，而在較大的摩擦作用時(shí)可以堆積成非常松散的試樣[26]。在二維模擬試驗(yàn)中必須確定試樣的相對密實(shí)度，因?yàn)轭w粒材料的強(qiáng)度和變形性質(zhì)受砂性土的相對密實(shí)度的影響，而顆粒的形狀會影響砂性土的極限孔隙比(相對密實(shí)度)[2]。通過足夠多數(shù)量的顆粒即可確定砂性土的極限孔隙率[20]，試樣一最小孔隙率為0.098，最大孔隙率為0.188，試樣二最小孔隙率為0.076，最大孔隙率為0.142。試樣一的極限孔隙比區(qū)間長度比試樣二的大且試樣二最大孔隙比比試樣一的小，上述特征符合具有較低的球形度、較高的棱角度和粗糙的表面的砂礫土相比于圓形砂礫土應(yīng)有較大的極限孔隙比差值的結(jié)論[13-14]?；诖?，為了與室內(nèi)試驗(yàn)對應(yīng)，采用50%的相對密度制成試樣以備剪切。

由于試樣為硅質(zhì)砂，且應(yīng)力水平相對較小(不超過400 kPa)，因此不考慮顆粒破碎對試驗(yàn)結(jié)果造成的影響。且認(rèn)為粒間接觸剛度不隨接觸力的增加而增加，可以采用線性接觸模型代替赫茲接觸模型。考慮到采用圓形顆粒試樣在模擬直剪試驗(yàn)時(shí)顆粒的轉(zhuǎn)動應(yīng)該受到限制，因此，粒間接觸模型采用抗旋轉(zhuǎn)的線性接觸模型[6]。通過顆粒和接觸處設(shè)置阻尼以保持剪切過程在準(zhǔn)靜態(tài)的情況下進(jìn)行。模擬直剪試驗(yàn)采用恒定的剪切速度進(jìn)行剪切(0.05 mm/s)，至剪位移達(dá)到10 mm(剪應(yīng)變約為16%)時(shí)停止。

通過反向擬合室內(nèi)試驗(yàn)中得到的應(yīng)力-應(yīng)變曲線校對模擬試驗(yàn)中的細(xì)觀參數(shù)。考慮到細(xì)觀參數(shù)與砂性土的極限孔隙率的相互影響，反復(fù)調(diào)試了細(xì)觀參數(shù)，最終細(xì)觀參數(shù)見表1和表2。砂性土的剪切強(qiáng)度或變形會因顆?？臻g分布的不同而在一定范圍內(nèi)波動[20]。為此，模擬直剪試驗(yàn)采用3種不同隨機(jī)數(shù)構(gòu)建模型，分別進(jìn)行平行試驗(yàn)，取其平均值作為模擬直剪試驗(yàn)的結(jié)果并計(jì)算其標(biāo)準(zhǔn)差。

表1 接觸的細(xì)觀參數(shù)Table 1 Meso-parameters of contact

表2 邊界墻和顆粒的細(xì)觀參數(shù)Table 2 Mesoscopic parameters of boundary walls and particles

3 試驗(yàn)結(jié)果及討論

3.1 標(biāo)定結(jié)果

直剪試驗(yàn)及其模擬試驗(yàn)的剪應(yīng)力-剪位移曲線、法向位移-剪位移曲線見圖4和圖5(帶誤差線的均為模擬直剪試驗(yàn)的結(jié)果?；诮７椒ê图?xì)觀參數(shù)，模擬直剪試驗(yàn)很好地再現(xiàn)了直剪試驗(yàn)的剪切過程。剪切至8～10 mm時(shí)，試樣的體積基本上處于恒定的狀態(tài)，即臨界狀態(tài)[17]。通過直剪試驗(yàn)和模擬直剪試驗(yàn)獲得的臨界摩擦角分別為33.9°和32.3°。基于剪切后的試樣所做的6組休止角試驗(yàn)，獲得試樣休止角為31.8°±0.6°。直剪試驗(yàn)和模擬直剪試驗(yàn)獲得的臨界摩擦角與試樣休止角的平均值分別相差2.1°和0.5°。Herle等[27]曾指出干砂臨界摩擦角φcv與其休止角非常接近，相差約±0.5°，并提出了用干砂休止角估算其臨界摩擦角的方法。Guo等[16]統(tǒng)計(jì)的磨圓度(roundness)較好的渥太華標(biāo)準(zhǔn)砂(Ottawa standard sand)休止角為31°。因此，直剪試驗(yàn)及其模擬直剪試驗(yàn)均具有很好的可信度。

圖4 直剪試驗(yàn)及其模擬直剪試驗(yàn)的剪應(yīng)力-剪位移曲線Fig.4 The shear stress-shear displacement curve of the direct shear test and its simulated direct shear test

圖5 直剪試驗(yàn)及其模擬直剪試驗(yàn)的法向位移-剪位移曲線Fig.5 Normal displacement-shear displacement curve of direct shear test and its simulated direct shear test

3.2 顆粒形狀對試樣宏觀力學(xué)性質(zhì)的影響

試樣二與試樣一相比，除了顆粒形狀上的差別外，試樣的相對密實(shí)度、顆粒級配、試驗(yàn)條件以及顆粒的數(shù)量均是相同的。為了研究法向應(yīng)力對研究結(jié)果的影響，除了采用22、50和100 kPa法向應(yīng)力外，還增加了法向應(yīng)力為200、300和400 kPa的模擬直剪試驗(yàn)。兩種試樣均在極限孔隙率下(密實(shí)狀態(tài)和松散狀態(tài))的進(jìn)行剪切試驗(yàn)。試樣密實(shí)狀態(tài)下的峰值強(qiáng)度取剪應(yīng)力能達(dá)到的最大值，臨界值強(qiáng)度取8～10 mm段的剪應(yīng)力平均值并計(jì)算其標(biāo)準(zhǔn)差，試樣松散狀態(tài)下的強(qiáng)度取剪應(yīng)力能達(dá)到的最大值。

按照Mohr-Coulomb準(zhǔn)則進(jìn)行直線擬合，得到不同法向應(yīng)力下密實(shí)狀態(tài)的兩種試樣的峰值強(qiáng)度和臨界值強(qiáng)度(圖6和圖7)。結(jié)果表明，當(dāng)法向應(yīng)力小于200 kPa時(shí)，試樣一峰值強(qiáng)度大于試樣二；當(dāng)法向應(yīng)力大于200 kPa時(shí)，試樣一峰值強(qiáng)度小于試樣二。隨著法向應(yīng)力增大，試樣二表現(xiàn)出較高的剪切強(qiáng)度。試樣二臨界值強(qiáng)度以及臨界摩擦角均略比試樣一的大，但差別甚小(試樣一和試樣二臨界摩擦角分別為32.6°和33°)。有文獻(xiàn)報(bào)道稱不規(guī)則的顆粒形狀會增加砂性土的剪切強(qiáng)度，例如，Kandasami等[4]在300 kPa圍壓下分別對由石英砂和玻璃球組成的砂性土做了空心圓柱體抗扭試驗(yàn)，得出石英砂不規(guī)則的形狀增加了試樣的剪切強(qiáng)度。然而,Ng[18]在100 kPa圍壓下基于4種具有不同長細(xì)比的橢球體試樣進(jìn)行三軸壓縮試驗(yàn)，其結(jié)論表明顆粒形狀對砂性土的峰值和臨界摩擦角影響很小。Cho等[2]表示顆粒形狀對砂性土力學(xué)特性的影響機(jī)理仍不清楚。在顆粒形狀影響砂性土極限孔隙比(相對密實(shí)度)的情況下，諸多文獻(xiàn)忽略了對試樣相對密實(shí)度的校核[2]，這也是文獻(xiàn)[10]中沒有注意到的。他們的結(jié)論中沒有提及顆粒形狀對剪切強(qiáng)度大小產(chǎn)生明顯影響的條件。模擬直剪試驗(yàn)結(jié)果表明，對于密實(shí)狀態(tài)的砂性土，只有在較高的應(yīng)力水平下顆粒不規(guī)則形狀才能明顯增加砂性土的剪切峰值強(qiáng)度，在較低的法向應(yīng)力時(shí)其剪切峰值強(qiáng)度反而比規(guī)則圓形砂性土的低。在文中整個(gè)法向應(yīng)力范圍內(nèi)，試樣二臨界值強(qiáng)度以及臨界摩擦角均略比試樣一的大，但差別很小，可以認(rèn)為顆粒不規(guī)則形狀沒有明顯增加砂性土的臨界值強(qiáng)度。

圖6 在密實(shí)狀態(tài)下的峰值強(qiáng)度Fig.6 Peak intensity in the dense state

圖7 在密實(shí)狀態(tài)下的臨界值強(qiáng)度Fig.7 Critical value strength in the dense state

圖8為松散狀態(tài)下的兩種試樣在不同法向應(yīng)力時(shí)的剪切強(qiáng)度?？梢钥闯觯嚇佣膬?nèi)摩擦角明顯比試樣一大，然在較低法向應(yīng)力時(shí)，極值強(qiáng)度上的差別幾乎可以忽略，隨著法向應(yīng)力的增加，兩種試樣極值強(qiáng)度差成比例增加。

圖8 在松散狀態(tài)下的極值強(qiáng)度Fig.8 The extreme strength in loose state

4 顆粒形狀對砂性土力學(xué)性質(zhì)的影響機(jī)理

4.1 直剪試驗(yàn)中試樣應(yīng)力的變化

自開始剪切至出現(xiàn)峰值強(qiáng)度的過程中，顆粒間接觸力方向由開始的隨機(jī)分布[圖9(a)]轉(zhuǎn)變?yōu)槎ㄏ蚍植糩圖9(b)]，即形成柱狀力鏈(column)。當(dāng)達(dá)到峰值強(qiáng)度時(shí)，剪切面上柱狀力鏈Fi方向近似相等，令其與法向夾角為θi，則法向應(yīng)力N、剪應(yīng)力T與Fi之間的關(guān)系式為

(1)

式(1)中：n為剪切面[圖9(b)中虛線]上柱狀力鏈的個(gè)數(shù)。剪切面上n個(gè)Fi矢量和可以用一個(gè)綜合的矢量F表示，其與法向夾角為θ(反算或模擬直剪試驗(yàn)中直接獲得)，則法向應(yīng)力N、剪應(yīng)力T與F的關(guān)系式為

(2)

剪應(yīng)力T等于0時(shí)，θ等于0，意味著沒剪切前，剪切面上綜合矢量F方向是豎直的。剪應(yīng)力逐漸增大到峰值過程中，θ亦逐漸增大到極值。一方面，顆粒緩慢移動，在剪切盒對角線方向上顆粒間的接觸越來越多，到峰值處最明顯；另一方面，對角線方向上顆粒間的接觸力逐漸增大，而其他方向上的接觸力逐漸向剪切盒對角線方向上偏轉(zhuǎn)、減小或消失 [圖9(a)、圖9(b)的插圖]。剪應(yīng)力出現(xiàn)峰值后，對角線方向上接觸力的增大使其法向分量超過加壓板的壓力，從而使剪切面以上的顆粒連同加壓板反向運(yùn)動，同時(shí)剪切面附近形成的柱狀力鏈繞其中心順時(shí)針轉(zhuǎn)動[22]，θ又開始減小，根據(jù)式(1)、式(2)可知，剪應(yīng)力T亦開始減小，剪切帶區(qū)域空隙增大，試樣體積膨脹。當(dāng)θi過小時(shí)，柱狀力鏈會發(fā)生潰曲，新的柱狀力鏈又會生成，剪應(yīng)力逐漸進(jìn)入到臨界值狀態(tài)[圖9(c)]。

圖9 密實(shí)試樣一剪切過程中不同階段顆粒間接觸力狀態(tài)(法向應(yīng)力300 kPa)Fig.9 The state of the contact force between particles at different stages during the shearing process of the compact specimen (normal stress 300 kPa)

根據(jù)式(2)中剪應(yīng)力T、法向應(yīng)力N和綜合矢量F與N夾角θ的關(guān)系，提出一種新的本構(gòu)關(guān)系如圖10所示，剪應(yīng)力T的大小可根據(jù)參數(shù)N和θ做出預(yù)測。根據(jù)Rowe剪脹理論，等粒徑的圓形顆粒極密排列情況下，綜合矢量F與法向夾角可達(dá)60°，因此0°≤θ≤60°。兩種試樣表明，剪切過程中θ的變化受顆粒形狀的影響，剪應(yīng)力達(dá)到峰值和臨界值時(shí)的θ分別為θp和θc。

圖10 不同法向應(yīng)力下剪應(yīng)力T和夾角θ的關(guān)系Fig.10 The relationship between the angle θ and the shear stress T under different normal stresses

兩種試樣的Tp-θp(Tp為峰值剪應(yīng)力)曲線均較好呈現(xiàn)雙曲線形態(tài)，低法向應(yīng)力時(shí)，試樣一θp偏大而試樣二偏小，高應(yīng)力時(shí)相反，θp均落于40°～50°范圍，甚至試樣二的θp區(qū)間更小。而對于Tc-θc(Tc為臨界剪應(yīng)力)曲線而言，試樣二依然呈現(xiàn)出較好的雙曲線形態(tài)，而試樣一在低應(yīng)力時(shí)離散性較大，試樣二總體上具有較大的θc，且在高應(yīng)力情況下更明顯，θc傾向于30°～40°范圍。正如文獻(xiàn)[19]中所述，密實(shí)試樣的θp、θc應(yīng)介于30°～60°。

4.2 剪脹作用的影響

Rowe[28]認(rèn)為顆粒間的摩擦作用和剪脹作用是剛性顆粒材料剪切強(qiáng)度的主要來源，提出應(yīng)力剪脹關(guān)系式為

(3)

式(3)中：σ11、σ22為試樣所受的最大和最小主應(yīng)力，Pa；dεv為體積應(yīng)變增量；dε1為剪切方向上的應(yīng)變增量；φu為顆粒表面的摩擦因數(shù)。式(3)說明，對于同一個(gè)試樣，其應(yīng)力比R=σ11/σ22會與剪脹因子D=1-dεv/dε1同步變化，當(dāng)應(yīng)力比達(dá)到峰值時(shí)，剪脹因子會出現(xiàn)最大值。

直剪試驗(yàn)中，有關(guān)剪脹角ψ的表達(dá)式為

(4)

式(4)中：dεx、dεy分別為試樣在剪切方向和法向上的變形增量；tanψ實(shí)質(zhì)上是法向位移-剪位移曲線的斜率。模擬直剪試驗(yàn)結(jié)果表明，剪應(yīng)力峰值所對應(yīng)的tanψ總是最大的，即最大剪脹角ψmax，不同法向應(yīng)力下兩種試樣的ψmax見圖11。Bolton[29]在砂土的三軸壓縮試驗(yàn)中得出最大剪脹角ψmax與峰值摩擦角φp的關(guān)系式為

圖11 不同法向應(yīng)力下兩種試樣的ψmaxFig.11 Two specimens under different normal stresses ψmax

φp=0.8ψmax+φc

(5)

式(5)中：φc為砂土的臨界摩擦角。

圖12為兩種試樣的θp、θc與ψmax的關(guān)系曲線，虛線為文獻(xiàn)[16]基于Ottawa標(biāo)準(zhǔn)砂做的三軸壓縮試驗(yàn)中的結(jié)果。雖然直剪試驗(yàn)與三軸試驗(yàn)關(guān)于參數(shù)的定義有所區(qū)別，但是它們的物理意義一樣，模擬直剪試驗(yàn)結(jié)果表明，θp、θc與ψmax同樣服從很好的線性關(guān)系，與文獻(xiàn)[16]的三軸試驗(yàn)結(jié)果非常吻合。

當(dāng)法向應(yīng)力較小時(shí)(22 kPa和50 kPa)，在剪應(yīng)力達(dá)到峰值時(shí)，顆粒體尚未出現(xiàn)明顯的剪切帶，體應(yīng)變不明顯，剪應(yīng)變很小。根據(jù)Barton等[30]提出的不平整的巖體節(jié)理面抗剪強(qiáng)度理論，在剪應(yīng)變較小以及法向應(yīng)力較小時(shí)，剪應(yīng)力T、法向應(yīng)力N的關(guān)系式為

T=Ntan(θc+ψ)

(6)

由圖12可知，在低法向應(yīng)力時(shí)，試樣一的峰值剪脹角比試樣二的大，而由圖7可知，兩種試樣的θc幾乎無差別，因而試樣一在低法向應(yīng)力下峰值強(qiáng)度較大主要原因是試樣一的剪脹角比試樣二的大些。

文獻(xiàn)[17]公式中，φp和φcv分別為三軸壓縮試驗(yàn)中試樣的峰值及臨界摩擦角，φmax為最大剪脹角圖12 兩種試樣θp、θc與ψmax的關(guān)系Fig.12 The relationship between two samples θp，θc and ψmax

圖13為圓形顆粒體和橢圓形顆粒體剪脹變形圖解，圓形顆粒與橢圓形顆粒的面積相同，假設(shè)顆粒間均未發(fā)生相對滑動，則柱狀顆粒體會發(fā)生轉(zhuǎn)動，當(dāng)發(fā)生相同剪切位移時(shí)(dεx=1.79)，圓形顆粒體的位移增量(dεy=1.03)要比橢圓形顆粒體的(dεy=0.74)大些，按照剪脹角的定義[式(4)]不難看出圓形試樣剪切時(shí)具有較大的剪脹角的原因。當(dāng)法向應(yīng)力較高(300 kPa和400 kPa)時(shí)，剪應(yīng)力達(dá)到峰值時(shí)發(fā)生了較大的剪切位移，兩種試樣的峰值剪脹角ψmax也趨于相同，此時(shí)式(6)不再適用。文獻(xiàn)[17]認(rèn)為多棱角的顆粒之間可形成鎖定(interparticle locking)結(jié)構(gòu)，增加了顆粒體的剪脹性，因而增加了顆粒體的峰值摩擦角φp。

圖13 橢圓形與圓形顆粒體變形圖解Fig.13 Illustration of deformation of ellipse and round particles

4.3 顆粒間摩擦作用的影響

模擬直剪試驗(yàn)表明，受摩擦作用影響的顆粒間的滑動和滾動及其相互轉(zhuǎn)化機(jī)制一定程度上控制著顆粒體的宏觀強(qiáng)度，隨著顆粒表面摩擦系數(shù)的增高，顆粒體的宏觀強(qiáng)度呈非比例增高[31]。當(dāng)抗旋轉(zhuǎn)阻力矩系數(shù)小于顆粒表面摩擦因數(shù)時(shí)，顆粒體強(qiáng)度會因其增大而增大，當(dāng)抗旋轉(zhuǎn)阻力矩系數(shù)超過顆粒表面摩擦因數(shù)時(shí)，顆粒體強(qiáng)度的增加不再明顯[32]。為此，對兩種試樣在100 kPa法向應(yīng)力下，顆粒表面摩擦因數(shù)在0～1變化時(shí)，密實(shí)狀態(tài)顆粒體的強(qiáng)度(峰值強(qiáng)度、臨界值強(qiáng)度)進(jìn)行了探討，在設(shè)定一個(gè)顆粒表面摩擦因數(shù)時(shí)，同時(shí)給予一個(gè)相同的抗旋轉(zhuǎn)組力矩系數(shù)，試驗(yàn)結(jié)果見圖14。

在沒有摩擦的情況下，兩種試樣的強(qiáng)度值極低，說明摩擦作用是無黏性粗顆粒體(相對于黏土或粉土而言的砂土或粗粒土等)強(qiáng)度的主要來源之一。當(dāng)摩擦因數(shù)為零時(shí)，試樣一中圓形顆粒間接觸力必然通過圓心，這樣在沒有切向力的情況下，圓形顆粒不受彎矩，因而不會發(fā)生轉(zhuǎn)動；相反，試樣二顆粒的形狀是不規(guī)則的，接觸力偏離形心，顆粒受到較大的偏心距，所有顆粒均發(fā)生較大的轉(zhuǎn)動，導(dǎo)致試樣一最終的體積膨脹量(1.77%)遠(yuǎn)比試樣二(0.38%)的大，其強(qiáng)度上的差異可能與此有關(guān)。另外，無摩擦作用時(shí)，兩種試樣內(nèi)顆粒間接觸力分布均勻，均無明顯的應(yīng)力局部化現(xiàn)象(即無強(qiáng)弱力鏈)。

Skinner[33]曾指出臨界狀態(tài)摩擦角幾乎不受顆粒間摩擦作用強(qiáng)弱的影響，然而Abriak[34]通過一系列的雙軸試驗(yàn)反駁了Skinner的結(jié)論，得出當(dāng)顆粒間摩擦角在0°～30°變化時(shí)，相應(yīng)的臨界狀態(tài)摩擦角從20°增加到26°。模擬直剪試驗(yàn)表明(圖14)，兩種試樣的峰值強(qiáng)度和臨界值強(qiáng)度均隨著顆粒表面摩擦因數(shù)和抗旋轉(zhuǎn)阻力矩系數(shù)的升高而非比例升高，這與文獻(xiàn)[1]中結(jié)論類似。在較低的摩擦作用下，試樣二的峰值強(qiáng)度比試樣一的低。當(dāng)顆粒表面摩擦因數(shù)較大時(shí)，試樣二的峰值強(qiáng)度更高，由此可以得出，顆粒間較大的摩擦作用，可使得試樣二峰值強(qiáng)度較試樣一高。當(dāng)顆粒間摩擦作用較弱時(shí)，臨界狀態(tài)剪切強(qiáng)度受其影響明顯，較強(qiáng)時(shí)影響不明顯。但是，摩擦作用的強(qiáng)弱不會改變不同形狀顆粒體臨界摩擦角似乎相等這一特征。

圖14 100 kPa法向應(yīng)力下剪切強(qiáng)度與顆粒表面摩擦因數(shù)和抗旋轉(zhuǎn)阻力矩系數(shù)比值的關(guān)系Fig.14 The relationship between shear strength and the ratio of particle surface friction coefficient and moment coefficient of resistance to rotation at 100 kPa normal stress

5 結(jié)論

采用PFC中Clump模塊構(gòu)建不規(guī)則形狀的顆粒并生成仿真砂樣，顆粒的形狀取自天然砂粒，一定程度上反映了砂性土中砂粒的形狀特征。通過室內(nèi)直剪試驗(yàn)，標(biāo)定了模擬直剪試驗(yàn)所需的細(xì)觀參數(shù)，嚴(yán)格地確定了不同試樣的極限孔隙率。設(shè)定不同的法向應(yīng)力條件，開展了一系列的模擬直剪試驗(yàn)。通過設(shè)置基于圓形顆粒的對比試驗(yàn)，揭示了天然砂性土中顆粒形狀對其宏細(xì)觀力學(xué)性質(zhì)的影響。定義了剪切面上顆粒間接觸力合力與法向應(yīng)力的夾角參數(shù)θ，數(shù)值上其正切函數(shù)值等于剪應(yīng)力T和法向應(yīng)力N的比值?；讦葮?gòu)造了新的本構(gòu)關(guān)系式，揭示了剪應(yīng)力變化的原因?；诩裘浗抢碚摚瑥念w粒層面揭示了法向應(yīng)力較低時(shí)不規(guī)則形狀顆粒體剪切強(qiáng)度較圓形試樣低的主要原因。得到如下結(jié)論。

(1)密實(shí)狀態(tài)下不規(guī)則形狀的砂粒組成的試樣較圓形顆粒試樣具有更大的峰值摩擦角和幾乎相等的臨界摩擦角；松散狀態(tài)下具有更大的內(nèi)摩擦角。

(2)在較低的法向應(yīng)力條件下密實(shí)不規(guī)則形狀顆粒體的剪切峰值強(qiáng)度比圓形顆粒體低，而在較高的法向應(yīng)力條件下比圓形顆粒體高。兩者臨界值強(qiáng)度差別甚小。

(3)法向應(yīng)力較低時(shí)不規(guī)則形狀顆粒體剪切強(qiáng)度較圓形試樣低的主要原因是圓形試樣具有較大剪脹角。