塊石強度及尺寸對土石混合體壓縮特性的影響分析

武小菲　張樂　江松

【摘 要】文章利用二維離散元程序建立了細觀結構均勻的土石混合體簡化顆粒流模型，模擬雙軸壓縮試驗分析塊石強度及尺寸對土石混合體壓縮過程和壓壞特性的影響規律。結果表明：塊石尺寸足以在土石混合體中形成骨架結構時，塊石才能充分發揮作用。土石混合體含塊石骨架時，其壓縮破壞過程可分為彈性變形階段、塊石骨架承力階段和整體破壞失穩階段。土石混合體中，塊石形成骨架結構時，隨塊石強度增大，彈性變形階段有所延長，塊石骨架結構承力作用顯著提升，整體失穩破壞的陡變性明顯減弱，殘余強度有效提高;塊石未形成骨架結構時，塊石強度對混合體強度影響不大，在結構破壞后對整體失穩有一定的延緩。大粒徑塊石會導致土石混合體在密實狀態下壓縮破壞后體積增加，塊石強度越高，尺寸越大，混合體壓縮破壞后體積增量越大。

【關鍵詞】土石混合體; 雙軸壓縮; 壓縮過程; 壓壞特性; PFC2D

【中圖分類號】TU411.5【文獻標志碼】A

1 土石混合體

土石混合體作為特殊的工程地質材料，其材料成分復雜、顆粒離散性大，且被廣泛應用于道路、堤壩等工程[1-2]。土石混合體的壓實過程，常伴隨著內部結構破壞與重組，而填料壓實的目的在于保證土石混合體具有足夠的強度與穩定性[3]。因此，有必要研究土石混合體的壓縮破壞性質。

當前，關于土石混合體的壓縮破壞有諸多研究。董云[4]通過室內試驗研究了含石量等對土石混合體強度的影響，發現含石量小于30 %時對強度影響較小，70 %左右強度達最大值。孫遜等[5]通過現場試驗發現，土石混合體壓實過程中，粗顆粒破碎成細顆粒后填充顆粒間孔隙，使粗細顆粒間彼此咬合，從而導致壓實度增大。劉勇等[6]發現塊石含量等因素對壓實效果的影響權重為：大塊石含量>小塊石含量>最大粒徑。柴賀軍等[7]認為土石混合體中顆粒最大粒徑對混合體的強度及變形特性均存在一定影響。曹文貴等[8]通過室內大型三軸試驗發現，巖性是土石混合體強度的主要影響因素之一。楊冰等[9]利用PFC3D研究了含石量對土石混合體骨架效應的影響，提出可進一步研究塊石形狀的作用及混合體壓碎過程。蔣明鏡等[10-12]發現，多種土體基于離散元的平面雙軸壓縮試驗模擬與室內三軸試驗結果在宏觀力學響應上近似。張學朋等[13]借助PFC2D中的Bonded particle model模型建立雙軸壓縮模型，比對室內試驗校正參數后，再現了花崗巖壓縮試驗全過程。

顯然，塊石強度和塊石尺寸對土石混合體的強度特性有較大影響，顆粒流模擬能有效地獲取混合體內部破壞特征。本文利用PFC2D建立簡化土石混合體二維顆粒流模型，考慮塊石尺寸和塊石強度2個因素各3種水平，模擬了9種工況下土石混合體的壓縮破壞過程及壓壞特性，以期獲取這2個因素對土石混合體壓縮特性的影響規律。

2 顆粒流模擬

2.1 模型建立

使顆粒整齊排列，模型尺寸取2.5 m×2.5 m，孔隙率僅8.52 %，以便于簡化分析，模型如圖1。模型邊界設置為4面足夠長的墻體，法向剛度設定為5×108 N/m，切向剛度和摩擦系數設置為0。小粒徑無黏性砂礫由粒徑分別為2.071 cm和5 cm的2種圓盤顆粒模擬，緊密接觸排列，模型中顆粒共4 901顆。大粒徑塊石通過設定黏結強度的長方體顆粒集合模擬，含不同尺寸大粒徑塊石的土石混合體模型見圖1。

雙軸壓縮過程中，通過設定橫向墻體的速度進行加載，豎向墻體施加隨圍壓（橫向偏應力）變化的速度進行圍壓伺服控制。計算過程中記錄模型的橫向偏應力-應變、豎向偏應力-應變和模型體積應變來分析土石混合體的壓縮破壞性質。無大粒徑塊石的細粒土顆粒模型中，設置圍壓100 kPa，顆粒間摩擦系數0.5，進行雙軸壓縮模擬。空白組模型的偏應力與豎向應變的關系曲線見圖2，模型壓縮破壞后見圖3。顯然，圍壓基本保持在100 kPa左右，僅在模型臨界破壞狀態時產生輕微波動。模型加載及伺服控制有效，混合體破壞形式表現為模型邊緣的細觀結構破壞和中部裂隙貫通。

2.2 模型參數

PFC模擬時，需設置顆粒單元、墻體單元的接觸剛度，顆粒間的摩擦系數和黏結強度等，從而設定接觸模型，與宏觀物質的力學性質建立聯系。本文分析時顆粒之間設置線性接觸剛度模型和接觸滑動模型，大粒徑塊石內部除上述2模型外，增設點黏結模型。顆粒法向及切向接觸剛度均設置為1×108 N/m。摩擦系數和黏結強度通過試算擬定。

調整顆粒模型中摩擦系數，在相同加載速度、不同圍壓下進行試算，豎向偏應力最大值見圖4。摩擦系數小于0.7時，呈現摩擦系數越大，豎向偏應力最大值越小的規律。摩擦系數等于0.9時，豎向偏應力最大值突增，且圍壓200 kPa時增幅更為顯著。圍壓100 kPa、摩擦系數為0.1和0.9時，空白組模型的豎向偏應力-應變曲線見圖5。結合圖2分析可知，摩擦系數越大，豎向偏應力最大值對應的豎向應變越小，即破壞前彈性模量越大;且隨著摩擦系數增加，豎向偏應力最大值后應力逐級降低的現象越明顯，即破壞后殘余變形增量越大。

這是因為，模型中孔隙率小，顆粒排列整齊且能夠相互咬合，內摩擦力大，細觀結構穩定。摩擦系數不大于0.7時，摩擦系數越大，顆粒模型彈性模量越大，相同加載速度時呈現出摩擦系數越大，豎向偏應力最大值越小的規律。因此，模擬時忽略細顆粒間摩擦系數的差異，統一取粒間摩擦系數為0.5，以分析多種工況下土石混合體的壓縮破壞性質。

假定土石混合體中大粒徑塊石內部顆粒的黏結強度在法向與切向相等，以便于簡化分析。圍壓100 kPa和200 kPa作用下，分別設置顆粒黏結強度為50 kPa、100 kPa和200 kPa，進行雙軸壓縮模擬，模型豎向偏應力最大值見圖6。顯然，2種圍壓下模型豎向偏應力最大值接近，且隨黏結強度的增大，豎向偏應力最大值近乎呈線性增長。100 kPa圍壓下，3種顆粒黏結強度對應的塊石強度水平為6.94 MPa、11.35 MPa和19.69 MPa。

圍壓100 kPa時，黏結強度分別取50 MPa、100 MPa和200 MPa，空白組模型的豎向偏應力-應變曲線見圖7。顯然，黏結強度越大，空白組模型的豎向偏應力峰值及峰值豎向應變越大。對比圖5可知，不同黏結強度下豎向偏應力可能存在2次峰值，且峰值前的彈性模量相近。隨黏結強度的增大，偏應力最大值較第1次峰值增幅減小，黏結強度為200 kPa時僅有1次峰值。據此，因此，文中后續取豎向偏應力最大值來分析多工況下土石混合體的壓縮破壞性質。

2.3 模擬方案

考慮塊石強度與塊石尺寸2種因素各3種水平，研究土石混合體的壓縮破壞性質，模擬工況見表1。塊石強度根據圖6，取100 kPa下3個顆粒黏結強度對應水平取值。通過不同的塊石尺寸反映粒徑和含量的差異，3個水平下模型見圖8，B1、B2和B3水平所對應的塊石含量為78 %、50 %和28 %。

3 模擬結果分析

3.1 豎向偏應力最大值

2種圍壓下，9種工況在雙軸壓縮下豎向偏應力峰值見圖9。可見，僅塊石尺寸B1（50 cm×75 cm）水平下的豎向偏應力峰值隨塊石強度的增大有明顯增長。結合圖8可知，土石混合體中未能形成塊石骨架結構時，3個塊石強度水平下土石混合體的強度幾乎不變。大粒徑塊石相互接觸形成骨架結構后，塊石才能充分發揮作用，隨塊石強度增大，土石混合體強度近乎呈線性提高。此外，100 kPa和200 kPa圍壓條件下混合體的豎向偏應力最大值規律一致，為簡化分析，下文壓縮破壞過程及體變特征分析僅取圍壓為100 kPa工況。

3.2 壓縮破壞過程

以工況A2B1為例，分析含塊石骨架結構土石混合體壓縮破壞過程，豎向偏應力-豎向應變曲線見圖10，壓縮破壞過程中不同階段的破壞特征見圖11。可知，含塊石骨架結構的土石混合體壓縮破壞過程中，首先出現小粒徑砂礫的剪切破壞，相應的豎向偏應力有輕微降低。然后，主要由大粒徑塊石形成的骨架結構負壓，荷載達一定值后，細顆粒發生細觀結構的整體破壞，大粒徑塊石發生剪切破壞。豎向偏應力達最大值后，模型整體失穩，荷載持續作用下大粒徑塊石逐步發生移位、變形和破壞。據此，土石混合體含塊石骨架結構時，其壓縮破壞過程可分為彈性變形階段、塊石骨架承力階段和整體破壞失穩階段。

不同塊石強度下含塊石骨架土石混合體的豎向偏應力-豎向應變曲線見圖12，整體破壞失穩階段見圖13。對比3種工況（A1B1、A2B1、A3B1）易知，土石混合體在壓縮過程中，隨塊石強度增大，彈性變形階段有所延長，塊石骨架結構承力作用顯著提升，整體破壞失穩階段陡變性明顯減弱，殘余強度有效提高。不同塊石強度下土石混合體整體破壞失穩特征差異較大，強度較低時，混合體內部形成傾斜剪切破壞面;隨塊石強度增長，混合體破壞呈現小粒徑的細觀破碎和大粒徑塊石的變形、移位等。

以A3B1、A3B2、A3B3為例，分析3種工況下相同塊石強度、不同塊石尺寸的土石混合體的壓縮破壞過程。豎向偏應力-豎向應變見圖14，整體破壞失穩階段見圖15。可知，大粒徑塊石在土石混合體內部不足以形成骨架結構時，豎向偏應力-應變曲線中僅出現1次峰值。塊石排列相對較密時，土石混合體的破壞主要為砂礫的細觀結構破壞和大粒徑塊石的邊角輕微破損，塊石的存在對土石混合體強度無明顯增強。塊石排列相對松散時，土石混合體的破壞主要表現為塊石間砂礫的剪切破裂，塊石的存在雖不能增強混合體的峰值強度，但殘余強度有一定的小幅回升，說明該類工況下大粒徑塊石可延緩土石混合結構的整體失穩。

3.3 體積應變

取土石混合體體積應變減小為正，建立各工況下體積應變與豎向應變的關系曲線見圖16。顯然，大粒徑塊石的強度越小，塊石尺寸越小，相應工況下的體積應變越小。各工況下模型整體破壞失穩后，曲線多出現小幅回升現象，即混合體的體積增加，與豎向偏應力-應變曲線表征規律一致。因此，大粒徑塊石會導致密實狀態的土石混合體在破壞后體積增加，且塊石強度越高，尺寸越大，體積應變增值越大。該規律與現場試驗結論一致，從側面論證了本文所得規律的合理性。

4 結束語

本文利用離散元程序PFC2D進行雙軸壓縮模擬，研究了塊石強度及尺寸對土石混合體壓縮過程與壓壞特性的影響規律，所得主要結論如下：

（1）大粒徑塊石相互接觸形成骨架結構后，塊石才能充分發揮作用，隨塊石強度增大，土石混合體強度近乎呈線性提高。

（2）含塊石骨架結構的土石混合體壓壞過程可分為3個階段：彈性變形階段、塊石骨架承力階段及整體破壞失穩階段。

（3）當塊石形成骨架時，土石混合體在壓縮過程中，隨塊石強度增大，彈性變形階段有所延長，塊石骨架結構承力作用顯著提升，整體破壞失穩階段陡變性明顯減弱，殘余強度有效提高;大粒徑塊石在混合體內部未形成骨架時，土石混合體的破壞形式與塊石排列的緊密程度相關，塊石的存在對土石混合體強度提升較小，但對混合體的整體失穩破壞有一定的延緩作用。

（4）大粒徑塊石會導致密實狀態下的土石混合體壓縮破壞后體積增加，且塊石強度越高，尺寸越大，混合體模型壓縮破壞后的體積增量越大。

參考文獻

[1]胡瑞林，李曉，王宇，等.土石混合體工程地質力學特性及其結構效應研究[J].工程地質學報，2020，28（2）：255-281.

[2]XU W J，HU R L，TAN R J.Some geomechanical properties of soil-rock mixtures in Tiger-leaping Gorge Area，China[J]. Geotechnique，2007，57（3）：255-264.

[3]江松.大粒徑無黏性土石混填路基壓實特性研究[D].成都：西南交通大學，2018.

[4]董云.土石混合料強度特性的試驗研究[J].巖土力學，2007，28（6）：1269-1274.

[5]孫遜，李志清，周應新，等.土石混填路基現場壓實試驗及其固結機理研究[J].工程地質學報，2018，26（2）：467-474.

[6]劉勇，劉暉，蒲文明，等.土石混合填料壓實特性的顆粒流模擬研究[J].路基工程，2020，208（1）：44-48.

[7]柴賀軍，閻宗嶺，賈學明.土石混填路基修筑技術[M].北京：人民交通出版社，2009.

[8]曹文貴，黃文健，王江營，等.土石混填體變形力學特性大型三軸試驗研究[J].湖南大學學報：自然科學版，2016，43（3）：142-148.

[9]楊冰，楊軍，常在，等.土石混合體壓縮性的三維顆粒力學研究[J].巖土力學，2010，31（5）：1645-1650.

[10]蔣明鏡，付昌，劉靜德，等.不同沉積方向各向異性結構性砂土離散元力學特性分析[J].巖土工程學報，2016，38（1）：138-146.

[11]蔣明鏡，李秀梅.雙軸壓縮試驗中砂土剪切帶形成的離散元模擬分析[J].山東大學學報：工學版，2010，40（2）：52-58.

[12]蔣明鏡，李濤，胡海軍.結構性黃土雙軸壓縮試驗的離散元數值仿真分析[J].巖土工程學報，2013，35（S2）：241-246.

[13]張學朋，王剛，蔣宇靜，等.基于顆粒離散元模型的花崗巖壓縮試驗模擬研究[J].巖土力學，2014，（z1）：99-105.

[定稿日期]2021-01-22

[基金項目]四川省重點研發項目（項目編號：2018JY0076）;四川省科技計劃重點項目（項目編號：2012YFS0321）

[作者簡介]武小菲（1974～），女，碩士，講師，主要研究方向為巖土工程及市政工程。