考慮加載速率效應(yīng)的黏土損傷特性研究

摘"要：對黏土的強度與變形特性的研究是深基坑及礦井建設(shè)的施工設(shè)計基礎(chǔ)，揭示不同固結(jié)壓力和加載速率對應(yīng)的黏土強度與損傷特性有助于更好地服務(wù)于地下結(jié)構(gòu)的施工設(shè)計。對不同固結(jié)壓力的黏土試樣進行不同加載速率的三軸試驗，分析黏土的應(yīng)力-應(yīng)變曲線及兩種因素（固結(jié)壓力、加載速率）對黏土強度的影響規(guī)律；對破壞后試樣CT掃描，獲取的試樣內(nèi)部裂隙的損傷分布揭示了三軸壓縮試樣呈現(xiàn)腰鼓狀的原因；在試驗的基礎(chǔ)上，假設(shè)黏土的強度服從Weibull分布，引入加載速率對強度參數(shù)的影響規(guī)律，得到了考慮加載速率效應(yīng)的黏土損傷本構(gòu)模型。

關(guān)鍵詞：黏土；土工試驗；加載速率；損傷本構(gòu)；CT

Investigation of clay damage characteristics with effect of loading rates considered

ZHU Hangqi^1，2， WANG Bo³， WANG Tao^1，2

（1.Beijing Urban Construction Survey and Design Institute Co.， Ltd.， Beijing 10008， China；

2.Beijing Key Laboratory of Deep Foundation Pit Geotechnical Engineering for Urban Rail Transit， Beijing 100084， China；

3.Suzhou Transportation Culture Tourism Investment Group Co.， Ltd.， Suzhou 234000， Anhui， China）

Abstract： The study of the strength and deformation characteristics of clay is the foundation of construction design for deep foundation pits and mines. Revealing the strength and damage characteristics of clay corresponding to different consolidation pressures and loading rates can help better serve the construction design of underground structures. This research conducts triaxial tests on clay samples with different consolidation pressures and loading rates to analyze the stress-strain curve of clay and the influence of two factors （consolidation pressure and loading rate） on clay strength. The CT scan of the damaged sample reveals the distribution of internal cracks and the reason for the bulging appearance of the triaxial compression sample. Based on the experiment and assuming that the strength of clay follows a Weibull distribution， the influence of loading rate on strength parameters is introduced， and a clay damage constitutive model considering the effect of loading rate is obtained.

Keywords： clay; geotechnical testing; loading rate; damage constitutive; CT

黏土作為一種常見的巖土材料，廣泛分布于表土層，其力學(xué)特性受到各種外部因素的影響，并與加載條件密切相關(guān)。完全彈性體材料的應(yīng)力應(yīng)變關(guān)系符合廣義胡克定律，因此其應(yīng)力應(yīng)變關(guān)系只取決于初始和最終的應(yīng)力狀態(tài)，與加載過程無關(guān)。但對于具有內(nèi)摩擦特性的彈塑性材料土體而言，其力學(xué)特性不僅取決于初始和結(jié)束的應(yīng)力狀態(tài)，還受到加載過程中各種因素的影響。多年來，諸多學(xué)者試驗研究了黏土的含水率、加載方式、溫度以及不同圍壓環(huán)境對黏土強度與變形特征的影響，提出了反映相關(guān)力學(xué)特性的若干數(shù)學(xué)模型（王海波等，2022；張衛(wèi)兵等，2021；張芳枝等，2009）。加載速率是影響?zhàn)ね亮W(xué)性質(zhì)的重要因素之一，目前的研究多集中于地表土的加載速率效應(yīng)（加載速率對土體力學(xué)性質(zhì)的影響）研究（De Cola et al.，2018；Omidvar et al.，2012；黨進謙等，2009），對固結(jié)黏土力學(xué)性質(zhì)的研究較少。

鄧亞虹等（2016）研究了晚更新世黃土受不同加載速率荷載作用產(chǎn)生的變形與強度特性，發(fā)現(xiàn)加載速率會顯著影響黃土的應(yīng)力-應(yīng)變特性，同時對黃土抗剪強度參數(shù)的影響規(guī)律不盡相同；李建中等（2008）通過黏土的單軸加載試驗，發(fā)現(xiàn)黏土的加載速率效應(yīng)隨含水量的降低有所降低。楊愛武等（2020）基于不同卸荷路徑和卸荷速率的三軸卸荷試驗，探究了卸荷路徑和卸荷速率這2種因素對吹填土孔壓、應(yīng)力應(yīng)變關(guān)系，強度及破壞特征的影響，并根據(jù)試驗數(shù)據(jù)構(gòu)建了考慮卸荷速率影響的卸荷強度破壞公式。朱寶龍等（2015）利用CT-多功能土工三軸儀，對處于固結(jié)排水試驗中的合肥地區(qū)非飽和黏土進行了CT掃描，發(fā)現(xiàn)土體在經(jīng)過短暫的壓密后，原有孔隙消失，產(chǎn)生的裂紋寬度隨時間擴展較快。從以上學(xué)者的研究成果可以看出，土在三軸試驗中表現(xiàn)出的力學(xué)特性與加載速率及固結(jié)應(yīng)力密切相關(guān)。

本文在分析不同因素影響?zhàn)ね翉姸忍匦缘幕A(chǔ)上，對淮南地區(qū)黏土試樣進行了不同固結(jié)應(yīng)力和加載速率的三軸壓縮試驗和破壞試樣的CT掃描，并根據(jù)試驗結(jié)果建立了黏土的損傷本構(gòu)模型，對施工過程中土體結(jié)構(gòu)失穩(wěn)現(xiàn)象的預(yù)測和防治有一定的參考意義。

1 "固結(jié)黏土的力學(xué)特性

1.1 "加載速率和圍壓的影響

試驗用土取自淮南地區(qū)的勘探鉆孔，其顆粒級配如圖1所示，其中黏土的級配依據(jù)GB/T 50123—2019。

為便于試驗條件的控制和重復(fù)，試驗采用重塑土制備試樣。參照GB/T50123—2019《土工試驗方法標準》關(guān)于試樣制備的要求，將取樣所得的原狀散土烘干、打散、篩選，分多次往配置好的干土中加水拌和均勻，稱取一定質(zhì)量的土放入標準土樣制樣模具，將樣品加工成直徑39.1 mm、高度80 mm的圓柱狀試樣，試樣的尺寸誤差控制在0.15 mm以內(nèi)。試樣最大干密度為2.135 g·cm^-3，飽和度為81.6%。

考慮到試驗參數(shù)的精確度，主要選用南京土壤儀器廠的TSZ型三軸儀進行試驗。由于重塑土的高壓固結(jié)時間較長，本文開展了3種圍壓（0.5 MPa、2 MPa和4"MPa）、5種應(yīng)變加載速率（0.001～0.05 mm·s^-1）下的三軸抗壓強度實驗。為了降低試樣差異對試驗結(jié)果的影響，每組試驗重復(fù)3次，剔除離散性較大的試驗結(jié)果后取平均值。試樣采用k₀固結(jié)方式，根據(jù)前人研究經(jīng)驗，固結(jié)系數(shù)取0.34（周國慶等，2016；陳國慶等，2014），加載應(yīng)力路徑如圖2所示。

首先分別以1 × 10^-3 mm·s^-1、5 × 10^-3 mm·s^-1、1 × 10^-2"mm·s^-1、5 × 10^-2 mm·s^-1、1 × 10^-1 mm·s^-1的加載速率對土樣加壓至設(shè)計壓力；然后穩(wěn)壓固結(jié)（軸壓、圍壓維持不變）至軸向應(yīng)變速率小于0.01 mm·h^-1；最后，保持圍壓穩(wěn)定不變，以初始加壓階段（固結(jié)階段前的加壓段）相同的加載速率施加軸向壓力至軸向應(yīng)變達到20%。

圖3為不同圍壓和加載速率試驗獲得的黏土試樣應(yīng)力-應(yīng)變曲線。從圖中可以看出，對應(yīng)不同圍壓，在彈性階段應(yīng)變加載速率對黏土的應(yīng)力應(yīng)變曲線影響較小，隨著軸向應(yīng)變的增加差異逐漸體現(xiàn)出來，在應(yīng)變屈服點處尤為明顯（應(yīng)變5%附近）。低圍壓（0.5 MPa）對應(yīng)的黏土應(yīng)力應(yīng)變曲線呈現(xiàn)出了峰值偏應(yīng)力（圖3 a），峰值后偏應(yīng)力呈下降趨勢，隨著圍壓的增加，曲線形式變?yōu)槔硐霃椝苄?；圍壓相同時，隨著加載速率的增加，黏土的屈服強度逐漸增大，此時黏土的應(yīng)變加載速率效應(yīng)逐漸顯現(xiàn)，隨著圍壓的增大，不同加載速率對應(yīng)的應(yīng)力應(yīng)變曲線差異逐漸減小，說明圍壓的增加會減弱黏土的應(yīng)變速率效應(yīng)。

當應(yīng)力應(yīng)變曲線存在峰值時，取峰值強度作為黏土的三軸壓縮強度；當應(yīng)力應(yīng)變曲線呈理想彈塑性時，取軸向應(yīng)變15%時的偏應(yīng)力作為黏土的抗壓強度。從圖4中可以看出黏土的抗壓強度隨著加載速率的增大而增大；同一加載速率對應(yīng)的黏土的抗壓強度隨著圍壓的增大而增大。在計算泊松比時，試樣的徑向變形取平均值，剪切時試樣面積A_t的校正，按照下式計算：

（1）

式中："是剪切過程中在t時刻測得的總體變；是加載桿侵入壓力室引起的液壓油體變，通過計算加載桿侵入體積來計算；"是相同壓力條件下無試樣時卸荷過程中液壓油的體變；以上數(shù)據(jù)均通過對三軸儀的標定獲得。

不同圍壓和加載速率對應(yīng)的黏土強度（T=σ₁-σ₃）、彈性模量（E）和泊松比（ν）見表1。

1.2 "黏土的CT掃描實驗

計算機層析成像（computerized tomography，CT）技術(shù)是一種無損傷探測物體內(nèi)部結(jié)構(gòu)的技術(shù)，它可以動態(tài)、定量以及無損傷地探測物體在受力過程中內(nèi)部結(jié)構(gòu)的變化。20世紀90年代以來，越來越多的科研人員將CT掃描技術(shù)運用到巖土材料研究中，并取得了一系列的成果。本文采用工業(yè)CT對處于2 MPa圍壓、加載速率為0.01 mm·s^-1的黏土試樣破壞前后的中心區(qū)域進行了掃描，掃描參數(shù)見表2。

通過土樣CT掃描實驗獲得的原始實驗數(shù)據(jù)是從CT機里讀取的以TIFF格式存儲的文件，圖像中的像素點是用十六位的信息表示的，實際上存儲的是掃描層面內(nèi)各點的灰度值，其值從0到65 536之間，這是一種特殊的圖像文件格式。本文采用MATLAB程序讀取TIFF格式的圖像信息和具有實用價值的特定信息，將圖像導(dǎo)入MATLAB程序中，根據(jù)提取的特定信息，確定黏土、礦物夾雜和裂隙3種成分對應(yīng)的灰度閾值，按圖像中3種成分閾值區(qū)域各自所占的像素數(shù)量，即可得到掃描層面中各成分所占的百分比。

剔除偽影嚴重的土樣外側(cè)部分，選取試樣中間區(qū)域重構(gòu)三維圖像，得到試樣和裂隙的三維重構(gòu)模型（圖5）。

對比原圖發(fā)現(xiàn)，經(jīng)過處理后的圖像，較好地保留了圖像的細節(jié)信息。其中綠色的區(qū)域為裂隙。從圖5中可以看出，CT掃描可以很好的重建試樣的結(jié)構(gòu)信息，裂隙可以經(jīng)過閾值分割技術(shù)完整的提取出來，這也有助于利用CT掃描重構(gòu)圖像對試樣破壞前后的細觀結(jié)構(gòu)進行定量分析。從圖5 c可以看出，試樣在破壞后，試樣上半部分的三維裂隙發(fā)育較試樣下半部分更為顯著，中間區(qū)域裂隙發(fā)育擴展集中，由此導(dǎo)致試樣中心區(qū)域密度變小，體積變大，這也是三軸壓縮試樣在破壞后呈現(xiàn)腰鼓狀的原因。

2 "考慮應(yīng)變速率影響的損傷本構(gòu)模型

2.1 "損傷變量的定義

黏土在高壓三軸加載條件下的破壞主要為塑性破壞，因此，黏土的強度準則應(yīng)為土體內(nèi)部出現(xiàn)塑性破壞時的應(yīng)力狀態(tài)?？梢詫ね猎谌S試驗中的破壞看作一個土樣內(nèi)部損傷連續(xù)發(fā)展的過程，土樣內(nèi)部孔隙、微裂隙在荷載作用下不斷發(fā)展、貫通，進而形成宏觀裂隙直至破壞（Shang et al.，2015；Sun et al.，2011；趙曉東等，2009；商翔宇等，2009）。圍壓和加載速率均會對黏土試樣的細觀結(jié)構(gòu)產(chǎn)生影響，如何合理地反映受圍壓和加載速率雙重影響的試樣損傷規(guī)律，關(guān)鍵在于損傷變量的選取。有關(guān)巖土材料損傷的研究表明，用彈性模量表示的損傷變量適用于彈塑性材料（張超等，2019；高娟等，2017；朱振南等，2019），因此，本文也采用彈性模量來定義圍壓和加載速率引起的黏土損傷變量。

目前，按損傷力學(xué)原理建立巖土材料的本構(gòu)模型主要有兩種途徑，一是從能量守恒的角度出發(fā)，以材料受損前后能量等價為基礎(chǔ)建立損傷模型；二是按變形前后應(yīng)變等價的原則，即按照Lemaitre（1985）的理論建立損傷模型。變形前后應(yīng)變等價的原理如式（2）所示：

（2）

式中，、分別表示未損傷及加載損傷后黏土的有效應(yīng)力，為基準圍壓下的黏土彈模，E_d為受損后試樣的彈性模量。用彈性模量來定義圍壓所致?lián)p傷，并結(jié)合虎克定律可得：

（3）

式中：為刻畫圍壓增加導(dǎo)致?lián)p傷（體現(xiàn)為初始彈性模量的減小）變量；現(xiàn)將圍壓增加導(dǎo)致的損傷定義為初始損傷，在圍壓施加的基礎(chǔ)上繼續(xù)按一定的加載速率"ε 施加軸壓，土樣在軸向應(yīng)力σ_ε作用下發(fā)生軸向變形ε_ε。再次利用推廣后的應(yīng)變等價原理得到第二損傷狀態(tài)的參數(shù)：

（4）

式中，、分別為黏土第二損傷狀態(tài)的有效應(yīng)力和等效彈性模量，為第二損傷變量。將式（3）代入式（4），得到：

（5）

式中，為黏土速率和圍壓耦合對應(yīng)的總損傷變量，為試樣損傷時所承受荷載。由式（5）可以看出，圍壓和加載速率共同作用使黏土總損傷表現(xiàn)出明顯的非線性特征，式中的反映了圍巖和加載速率的耦合作用。

本文借鑒曹文貴等（2012）和姚兆明等（2020）在巖土材料損傷本構(gòu)模型中所采用的處理方法，假定黏土由無數(shù)形狀和強度各異的微結(jié)構(gòu)單元（微元）組成，在試樣受到荷載作用時，強度較低的微元首先損傷與破壞，并服從某種統(tǒng)計規(guī)律，據(jù)此，黏土的損傷破壞過程可以采用連續(xù)介質(zhì)力學(xué)中的損傷理論進行描述。設(shè)黏土的損傷變量D為外界荷載作用下破壞的微元數(shù)量N_d與土體未受損的總微元數(shù)量N之比，則有：

（6）

式中，D為黏土考慮圍壓和加載速率時的損傷變量。D=0時，黏土處于無損狀態(tài)；D=1時，黏土處于完全損傷狀態(tài)；0lt;Dlt;1時，黏土處于不同程度的損傷狀態(tài)。

在荷載的作用下，黏土內(nèi)部土微元的破壞是隨機的，與土微元的強度T有關(guān)。假定土微元的破壞概率隨微元強度T的概率密度函數(shù)為P（T），當外界荷載達到某一應(yīng)力水平T時，黏土的微元破壞數(shù)量為：

（7）

假設(shè)土微元的破壞概率密度函數(shù)P（T）服從代入Weibull分布，即：

（8）

式中，m及F₀為Weibull分布參數(shù)，反映了分布曲線的的形態(tài)。將式（8）代入式（7）積分求出N_d，進一步代入式（6），由此得到了三軸壓縮條件下基于Weibull分布的荷載損傷變量表達式：

（9）

2.2 "建立基于Weibull分布的損傷本構(gòu)模型

土微元強度隨機分布變量的形式由土的破壞機理和強度準則決定，不同的強度準則，T有不同的表達形式。常用的經(jīng)典強度理論包括Mohr-Cloulomb屈服準則，Drucker-prager準則，Duncan準則等，其中Mohr-Cloulomb屈服準則是以實驗為基礎(chǔ)建立的強度準則，能很好的反映土的強度對荷載的敏感性，本文通過試驗測得的黏土強度可以采用基于加載速率修正的Mohr-Cloulomb屈服準則進行描述。

利用主應(yīng)力表示的Mohr-Cloulomb屈服準則可以表示為：

（10）

式中，φ為黏土的內(nèi)摩擦角，c為黏土的黏聚力。用式（10）黏土三軸壓縮試驗結(jié)果進行擬合，可以得到黏土的黏聚力和內(nèi)摩擦角隨加載速率的變化，見表3。

將表3中的強度參數(shù)（c、φ）相對于加載速率進行擬合，擬合結(jié)果如圖6所示。

由此可得到考慮加載速率的黏土抗壓強度預(yù)測方程。結(jié)合式（5）、式（9）和式（10）可得基于Weibull分布的黏土加載速率損傷演化方程為：

（11）

根據(jù)應(yīng)變等價原理及廣義的胡克定律，聯(lián)立式（5）可以得到黏土三軸壓縮條件下的統(tǒng)計損傷本構(gòu)方程：

（12）

3 "結(jié)論

1）加載速率對黏土的應(yīng)力應(yīng)變曲線具有顯著的影響。在彈性階段加載速率對黏土的應(yīng)力應(yīng)變曲線影響較??；隨著軸向應(yīng)變的增加差異逐漸體現(xiàn)出來，在應(yīng)變屈服點處尤為明顯（應(yīng)變5%附近）。低圍壓（0.5 MPa）對應(yīng)的黏土應(yīng)力應(yīng)變曲線存在峰值偏應(yīng)力，大于0.5"MPa的圍壓對應(yīng)的應(yīng)力應(yīng)變曲線變?yōu)槔硐霃椝苄?。隨圍壓的增加，不同加載速率對應(yīng)的應(yīng)力應(yīng)變曲線差異逐漸減小，因此圍壓的增加會減弱黏土的應(yīng)變速率效應(yīng)。

2）CT掃描可以有效揭示出黏土破壞后的裂隙發(fā)育情況。試樣在破壞后，上部裂隙發(fā)育比下部更密，中間區(qū)域裂隙發(fā)育擴展集中，由此導(dǎo)致試樣中心區(qū)域密度小，體積變大，這也是三軸壓縮試樣在破壞后呈現(xiàn)腰鼓狀的原因。

3）通過試驗得到的黏土強度可以采用基于加載速率修正的Mohr-Cloulomb屈服準則進行描述。通過考慮加載速率對黏土彈性模量的影響，建立了一種基于Weibull分布的黏土損傷本構(gòu)模型。

參考文獻

曹文貴，趙衡，李翔，張玲，2012.基于殘余強度變形階段特征的巖石變形全過程統(tǒng)計損傷模擬方法［J］.土木工程學(xué)報，45（6）：139-145.

陳國慶，王衍森，冷陽光，張春虎，2014.長時高壓K0固結(jié)黏土的卸載力學(xué)特性試驗研究［J］. 巖石力學(xué)與工程學(xué)報，33（S1）：2 996-3 002.

黨進謙，蔣倉蘭，吉中亮，2009.剪切速率對結(jié)構(gòu)性黃土力學(xué)性狀的影響［J］.地下空間與工程學(xué)報，5（3）：459-462.

鄧亞虹，魏寶華，王晗，李喜安，宋焱勛，2016.原狀黃土變形與強度特性的加載速率效應(yīng)［J］.中國公路學(xué)報，29（7）：22-29.

高娟，賴遠明，2017.凍結(jié)鹽漬土三軸剪切試驗過程中的損傷及壓融分析［J］.巖土工程學(xué)報，40（4）：707-715.

李建中，徐力生，2008.低含水量粘土的加載速率效應(yīng)與蠕變變形［J］.煤田地質(zhì)與勘探（1）：41-44.

商翔宇，2009. 不同應(yīng)力水平深部粘土力學(xué)特性研究［D］.徐州：中國礦業(yè)大學(xué).

王海波，呂偉華，武荘，朱文波，2022.不同溫度應(yīng)力路徑下飽和黏土剪切特性［J/OL］.巖土力學(xué)，43（3）：679-687.

楊愛武，楊少坤，張振東，2020.基于不同卸荷速率與路徑影響下吹填土力學(xué)特性研究［J］.巖土力學(xué)，41（9）：2 891-2 900.

姚兆明，余文，麻世壟，郭夢圓，2020.凍結(jié)黏土單軸抗壓試驗及損傷本構(gòu)模型分析［J］.力學(xué)與實踐，42（2）：196-201.

張超，楊期君，曹文貴，2019.考慮峰值后區(qū)應(yīng)力跌落速率的脆巖損傷本構(gòu)模型研究［J］.巖土力學(xué)，40（8）：3 099-3 106.

張芳枝，陳曉平，2009.非飽和黏土變形和強度特性試驗研究［J］.巖石力學(xué)與工程學(xué)報，28（S2）：3 808-3 814.

張衛(wèi)兵，莫劍鵬，容賢沂，唐明亮，2021.基于三軸壓縮試驗的黏土抗剪強度影響因素研究［J］.西部交通科技（8）：60-62.

趙曉東，周國慶，王博，2009.深部重塑土高壓應(yīng)力路徑試驗研究［J］.中國礦業(yè)大學(xué)學(xué)報，38（4）：471-475.

周國慶，況聯(lián)飛，馬金榮，商翔宇，趙曉東，2016.深土力學(xué)特性研究現(xiàn)狀及發(fā)展展望［J］.中國礦業(yè)大學(xué)學(xué)報，45（2）：195-204.

朱寶龍，巫錫勇，李曉寧，魏久平，2015.合肥地區(qū)重塑黏性土細觀結(jié)構(gòu)演化三軸CT試驗［J］.西南交通大學(xué)學(xué)報，50（1）：144-149.

朱振南，蔣國盛，田紅，吳文兵，梁榮柱，竇斌，2019.基于Normal分布的巖石統(tǒng)計熱損傷本構(gòu)模型研究［J］.中南大學(xué)學(xué)報（自然科學(xué)版），50（6）：1 411-1 418.

De Cola F，"Pellegrino A，"Gloessner C，"Penumadu D，"Petrinic N，"2018. Effect of Particle Morphology，"Compaction，"and Confinement on the High Strain Rate Behavior of Sand［J］. Experimental Mechanics，"58（2）：223-242.

Lemaitre J，"1985. A Continuous Damage Mechanics Model for Ductile Fracture ［J］. Transactions of the Asme Journal of Engineering Journal of Engineering Materials and Technology，"107：83-89.

Omidvar M，"Iskander M，"BLESS S，"2012. Stress-strain behavior of sand at high strain rates［J］. International Journal of Impact Engineering，"49：192-213.

SHANG Xiangyu，"ZHOU Guoqing，"LU Yong，"2015. Stress-dependent undrained shear behavior of remolded deep clay in East China ［J］. Journal of Zhejiang University-Science A（Applied Physics amp; Engineering），"16（3）：171-181.

SUN Ruhua，"LI Wenping，"WANG Weili，"2011. Research on deep clay high pressure triaxial unloading and structural constitutive model ［J］. International symposium on geomechanics and geotechnics："from micro to macro. Geomechanics and geotechnics："from micro to macro（1）：381-394.