張民生，申志聰，王振豪，孔令明

(1. 北京航空航天大學(xué)交通科學(xué)與工程學(xué)院，北京 100191；2. 中國(guó)海洋大學(xué)環(huán)境科學(xué)與工程學(xué)院，山東 青島 266100；3. 北京建筑大學(xué)土木與交通工程學(xué)院，北京 100044)

砂土是一種結(jié)構(gòu)簡(jiǎn)單但在高壓下易顆粒破碎的散粒組合體，其基本力學(xué)特性主要有壓硬性、摩擦性及剪脹性。由于顆粒間結(jié)構(gòu)松散，極易擾動(dòng)，制備土樣孔隙比變化較大，導(dǎo)致砂土有許多壓縮線，這是砂土區(qū)別于粘土的特征之一。砂土顆粒間接觸面積極小，壓縮過(guò)程中應(yīng)力大于顆粒強(qiáng)度，從而產(chǎn)生顆粒破碎現(xiàn)象，這一力學(xué)特征也是砂土區(qū)別于粘土的顯著標(biāo)志。

在寬應(yīng)力變化范圍內(nèi)，砂土的等向壓縮具有以下特征：低壓條件下，砂土壓縮性較差，幾乎不可壓縮，體積變化主要以顆粒間的彈性變形為主，在V-lgp空間內(nèi)壓縮曲線變化緩慢；隨著壓力增加，砂土的壓縮性顯著增大，在一定應(yīng)力區(qū)間壓縮曲線斜率變化明顯，該區(qū)間被認(rèn)為是土顆粒破損或破碎的標(biāo)志；在高應(yīng)力區(qū)，不同密實(shí)度的試樣壓縮曲線趨向歸并為同一條曲線[1-5]。

現(xiàn)有模型在描述砂土壓縮時(shí)主要采用單對(duì)數(shù)方法、雙對(duì)數(shù)方法、冪次方法及非線性函數(shù)方法。單對(duì)數(shù)方法假設(shè)孔隙比與壓力的對(duì)數(shù)呈線性關(guān)系，這種方法對(duì)于粘土有很好的適用性，劍橋模型及修正劍橋模型[6]均采用該方法，雖然一些學(xué)者也采用此方法描述顆粒材料[7-8]，但該方法在描述低應(yīng)力時(shí)存在一定的局限性。而且，單對(duì)數(shù)方法中孔隙比與固結(jié)(壓縮)壓力呈對(duì)數(shù)關(guān)系，當(dāng)壓力足夠大時(shí)，孔隙比出現(xiàn)負(fù)值，這不符合試驗(yàn)規(guī)律。為克服單對(duì)數(shù)方法無(wú)法真實(shí)反映物理試驗(yàn)現(xiàn)象的問(wèn)題，一些學(xué)者[2,9-10]提出了雙對(duì)數(shù)法描述粘土固結(jié)特性，而Sheng等[2]利用該方法很好地描述了砂土低應(yīng)力下的非線性壓縮特性及臨界狀態(tài)線。冪次方法起初基于砂土彈性模量與壓縮壓力或垂直壓力的冪次方呈正比[11]。由于砂土彈性模量的冪次表達(dá)方法在工程中受到廣泛認(rèn)可，且真實(shí)砂土臨界狀態(tài)線也符合冪次規(guī)律[5]，因此建立一個(gè)基于冪次變化規(guī)律的壓縮曲線將對(duì)本構(gòu)模型的建立提供極大方便。同時(shí)冪次函數(shù)用于描述低應(yīng)力下砂土壓縮特性[12]也具有一定的適用性，而且基于冪次函數(shù)描述砂土等向壓縮特性已取得了一定的效果[1,3-4,13-16]。

本文在分析砂土的等向壓縮特征基礎(chǔ)上，依托臨界狀態(tài)線建立砂土等向壓縮參考線，并通過(guò)塑性修正方法獲取其它壓縮線的斜率。建立的模型能夠同時(shí)考慮塑性變形與彈性模量。利用已有砂土試驗(yàn)數(shù)據(jù)對(duì)模型進(jìn)行驗(yàn)證。

1 砂土壓縮模型

1.1 高壓下砂土的極限孔隙比

常壓下制備砂土存在最小孔隙比，隨著壓力增加最小孔隙比也會(huì)逐漸減小，極高壓力下砂土顆粒發(fā)生破碎形成多個(gè)細(xì)小顆粒，導(dǎo)致顆粒間接觸空間減小。當(dāng)顆粒破碎到一定程度時(shí)不再發(fā)生破碎，接觸空間也不再減小，從而形成了高壓下的極限孔隙比。由于顆粒力學(xué)特性不同，不同的砂土破碎后形成的粒徑大小及形狀各不相同，因此高壓下的極限孔隙比沒(méi)有統(tǒng)一值。通過(guò)分析已有關(guān)于高壓下不同砂土壓縮試驗(yàn)[17]及剪切試驗(yàn)[18-19](試驗(yàn)壓力大于50 MPa，見圖1)，可以看出試驗(yàn)停止時(shí)，孔隙比主要在0.05～0.2區(qū)間內(nèi)分布，即使在1 000 MPa下砂土的極限孔隙比仍然不為0，因此描述砂土壓縮規(guī)律時(shí)應(yīng)參考極限孔隙比影響。

圖1 極限孔隙比隨壓力分布Fig.1 Limited void ratio vary with pressure in V-lgp

1.2 砂土參考?jí)嚎s線

飽和粘土的正常固結(jié)線是描述其在不同應(yīng)力歷史下固結(jié)特性的參考線。由于粘土顆粒的物理化學(xué)作用，正常固結(jié)線相對(duì)較為穩(wěn)定且唯一，即土的應(yīng)力狀態(tài)與土的孔隙比之間存在著唯一性關(guān)系，而且這種關(guān)系可以通過(guò)粘土的等向固結(jié)試驗(yàn)獲取。

而砂土顆粒構(gòu)成則較為簡(jiǎn)單，容易擾動(dòng)，導(dǎo)致低應(yīng)力下制備的試樣孔隙比變化較大，進(jìn)而使得砂土的壓縮線不唯一，也就是土的應(yīng)力狀態(tài)與砂土孔隙比之間唯一性關(guān)系不成立，從而增加了砂土壓縮特性數(shù)學(xué)描述的困難。因此，尋找合適的等向壓縮線作為參考線是有效描述砂土等向壓縮特性關(guān)鍵，本文將該曲線稱為參考?jí)嚎s線。通過(guò)砂土的等向壓縮試驗(yàn)可以發(fā)現(xiàn)不同孔隙比的砂土等向壓縮線在高壓下最終均會(huì)匯聚于同一條線[1]。參考?jí)嚎s線作為其中的一條等向壓縮線，最終也會(huì)與其它等向壓縮線在高壓下相匯聚。同時(shí)，Javanmardi等[20]與Sheng等[2]認(rèn)為在雙對(duì)數(shù)空間內(nèi)同樣存在一條參考狀態(tài)線,該線與臨界狀態(tài)線(CSL)在壓力為1 kPa時(shí)相交，并且在壓力無(wú)限大時(shí)與臨界狀態(tài)線近似平行。姚仰平等[21]認(rèn)為在壓力為0時(shí)，粘土的固結(jié)狀態(tài)與臨界狀態(tài)應(yīng)該為同一狀態(tài)，這意味著兩者在壓力為0時(shí)相交?；谇叭说难芯炕A(chǔ)[2,20-21]，本文建議砂土參考?jí)嚎s線特征：在壓力為0時(shí)與臨界狀態(tài)線相交，在高壓時(shí)與其它壓縮線匯聚。

砂土等向壓縮線及臨界狀態(tài)線(CSL)在V-lgp空間內(nèi)的非線性變化特征促使學(xué)者尋找更好的數(shù)學(xué)函數(shù)來(lái)描述其壓縮規(guī)律。為考慮極限孔隙比對(duì)等向壓縮影響，本文采用V體積(V=1+e)對(duì)數(shù)與壓力的冪次關(guān)系來(lái)描述砂土臨界狀態(tài)線和參考?jí)嚎s線。以Toyoura[22]砂為例，采用以下公式擬合砂土CSL線為：

(1)

式中:p為壓縮壓力；Vc為p對(duì)應(yīng)的背景狀態(tài)體積；Vmin為高壓下的極限體積；Vc0為臨界狀態(tài)線上p為0時(shí)的體積；λc為曲線斜率；ξ為材料常數(shù),對(duì)于Toyoura砂取值2/3。從式(1)可以看出CSL線在ln(V-Vmin)-(p/pa)ξ空間內(nèi)為直線段且斜率為λc。從圖2中可以看出利用式1擬合的Toyoura砂臨界狀態(tài)線數(shù)據(jù)離散性較低，擬合效果較為理想。

根據(jù)假設(shè)，砂土的參考?jí)嚎s線與臨界狀態(tài)線在0點(diǎn)相交，且數(shù)學(xué)表達(dá)方式與臨界狀態(tài)線相似，因此參考?jí)嚎s線的表達(dá)式可表示為：

(2)

式中:Vr為參數(shù)壓縮線上p對(duì)應(yīng)的體積；Vr0為參考?jí)嚎s線上壓力為0時(shí)的體積；且Vr0=Vc0；λr為曲線斜率。如圖2，取等體積下臨界狀態(tài)線和壓縮線上A點(diǎn)和B點(diǎn)對(duì)應(yīng)的壓力pA與pB，分別代入式(1)與式(2)，并聯(lián)立可得

(3)

圖2 V-lgp空間內(nèi)臨界狀態(tài)線及參考?jí)嚎s線變化規(guī)律[22]Fig.2 The critical state line and reference compression curve in V-lgp space[22]

由式3可見，參考?jí)嚎s線與CSL線在任意體積下所對(duì)應(yīng)的壓力比值恒定，這與粘土的性質(zhì)是完全一致，說(shuō)明用對(duì)數(shù)-冪次函數(shù)描述砂土壓縮線具有合理性。圖3為參考?jí)嚎s線在V-lgp空間內(nèi)分布規(guī)律。從圖中可以看出：低壓時(shí)曲線斜率接近于0，體積隨壓力變化極??；隨著壓力增加曲線斜率急劇變化后并在一定壓力范圍內(nèi)保持相對(duì)穩(wěn)定；在較大高壓階段曲線又逐漸變?yōu)榫徠?，此時(shí)已趨近臨界體積。

圖3 V-lgp平面內(nèi)參考?jí)嚎s線變化規(guī)律

根據(jù)參考?jí)嚎s線與臨界狀態(tài)線之間的關(guān)系，在臨界狀態(tài)線已知條件下，可確定臨界狀態(tài)在壓力為0時(shí)的體積Vc0。制備體積為Vc0的砂土進(jìn)行等向壓縮即可得到參考?jí)嚎s線。此外，根據(jù)參考?jí)嚎s線的定義，可根據(jù)不同壓縮線在高壓下匯聚于一條線這一基本試驗(yàn)規(guī)律而進(jìn)行擬合得到。圖4為Cambria砂[23]的壓縮試驗(yàn)及三軸剪切試驗(yàn)結(jié)果。在假定Vmin為0.17，ξ為2/3前提下，利用公式1對(duì)臨界狀態(tài)線進(jìn)行擬合得到Vc0，然后利用公式2調(diào)整參考?jí)嚎s線斜率λr使與其在高壓下與壓縮曲線匯聚，從而得到λr的確定值。砂土在高壓下會(huì)發(fā)生顆粒破碎，圖中Cambria砂的臨界狀態(tài)和等向壓縮試驗(yàn)均包含了砂土的顆粒破碎，本文建立的參考?jí)嚎s線函數(shù)基于等向壓縮試驗(yàn)和三軸剪切試驗(yàn)，因此曲線函數(shù)已將顆粒破碎及壓縮作為一個(gè)整體來(lái)考慮。從對(duì)比結(jié)果可以看出參考?jí)嚎s線在描述包含砂土顆粒破碎時(shí)仍然具有適用性。

圖4 Cambria砂的臨界狀態(tài)線及參考?jí)嚎s線[23]

1.3 其它壓縮線

對(duì)于其它壓縮線的變化規(guī)律，本文采用增量形式表達(dá)，其思想是通過(guò)狀態(tài)參量來(lái)修正參考?jí)嚎s線。對(duì)式(2)進(jìn)行微分可得：

(4)

假設(shè)砂土的彈性變形同樣符合式(2)變化規(guī)律，則參考?jí)嚎s線上當(dāng)前體積對(duì)應(yīng)的彈性增量可表示為：

(5)

式中κ為ln(V-Vmin)-(p/pa)ξ空間內(nèi)的砂土彈性變形系數(shù)。式5表明砂土的彈性變形不僅與壓力有關(guān)而且與孔隙比有關(guān)。Pestana等[1]及Chang等[24]采用過(guò)類似的方式來(lái)描述砂土的彈性變形。本文進(jìn)一步考慮了極限體積的影響。砂土塑性變形可表示為：

(6)

由式6可以看出參考?jí)嚎s線上塑性變形的增量同樣與壓力和體積有關(guān)。其它體積下砂土的塑性變形可利用塑性修正函數(shù)來(lái)修正，即：

(7)

式中：kn為恒大于0的材料參數(shù)；ψ為e-lgp空間內(nèi)當(dāng)前體積到參考?jí)嚎s線的豎直距離(如圖3所示)，即：

ψ=V-Vr。

(8)

將式7代入式6可得

(9)

最終得到不同體積下的壓縮增量表達(dá)式為：

(10)

砂土體積應(yīng)變?cè)隽縿t可表示為：

(11)

式中采用塑性修正方法求取其它體積下的壓縮增量具有明顯優(yōu)勢(shì)。式10中，exp(knψ)>0，λr-κ>0，則exp(knψ)(λr-κ)>0，這就保證了修正后的壓縮曲線斜率恒大于κ，從而避免了數(shù)值計(jì)算過(guò)程中壓縮變形呈現(xiàn)負(fù)增長(zhǎng)的情形，保證數(shù)值計(jì)算穩(wěn)定性。當(dāng)前孔隙比位于參考?jí)嚎s線上方時(shí)ψ>0，由于kn為恒大于0則exp(knψ)>1，(λr-κ)exp(knψ)+κ>λr，說(shuō)明壓縮線向參考?jí)嚎s線趨近。當(dāng)前孔隙比位于參考?jí)嚎s線上方時(shí)ψ<0，exp(knψ)<1，(λr-κ)exp(knψ)+κ<λr，則等向壓縮線相對(duì)于參考?jí)嚎s線較為平緩。

1.4 彈性模量

假設(shè)砂土彈性變形在ln(V-Vmin)-(p/pa)ξ空間內(nèi)只與壓縮壓力有關(guān)且符合線性關(guān)系，則彈性變形增量可表示為：

(12)

則體積彈性模量K及剪切模量G則可表示為：

(13)

(14)

雖然假設(shè)砂土彈性變形在ln(V-Vmin)-(p/pa)ξ空間內(nèi)只與壓縮壓力有關(guān)，但式13已體現(xiàn)出體積對(duì)體積彈性模量K的影響。從式13可以看出體積彈性模量K不僅隨壓縮壓力升高而增加，也隨著體積的減小而增加，很好地反映了砂土的壓硬性及彈性模量的密實(shí)度相關(guān)性。

1.5 模型參數(shù)

建立的砂土壓縮模型包括6個(gè)參數(shù)，Vmin、Vr0、λr、κ、ξ及kn。其中ξ為材料常數(shù)，當(dāng)材料為粘土?xí)r取值為1[25]，對(duì)于各向同性壓縮的顆粒材料，ξ趨近于2/3[1,24]，有些學(xué)者將ξ取值于0.7[5]或1/2[3,26]，本文取2/3。Vmin可取200 MPa壓力對(duì)應(yīng)的體積，對(duì)于Toyoura砂本文取值為1.2(見表1)。

表1 模型參數(shù)列表

κ為回彈曲線斜率，可在ln(V-Vmin)-(p/pa)ξ空間內(nèi)利用砂土回彈曲線及最密壓縮線擬合得到(如圖5)。通過(guò)對(duì)Toyoura砂[22]壓縮試驗(yàn)回彈曲線及最密壓縮線的擬合得到曲線斜率分別為0.005、0.004 6及0.005 9(見圖5)，本文取0.005 3。

Vr0與λr為參考?jí)嚎s線所對(duì)應(yīng)的初始體積和斜率，結(jié)合Toyoura砂臨界狀態(tài)線等向壓縮試驗(yàn)，確定λr為0.013，Vr為1.944。

參數(shù)kn的數(shù)值可由姚仰平等[27]的硬化參量及Dafalias等[28]提出的狀態(tài)參量確定。kn=4(nf-nd)，由文獻(xiàn)nf為3.5，nd為-1.7，得到kn為20.8。

圖5 κ標(biāo)定結(jié)果

2 模型驗(yàn)證

圖6為Toyoura砂土的壓縮試驗(yàn)結(jié)果與模型預(yù)測(cè)對(duì)比，試驗(yàn)數(shù)據(jù)來(lái)源于Verdugo等[22]的飽和砂土的等向壓縮試驗(yàn)，初始孔隙比分別為0.6、0.81、0.9及1.02對(duì)應(yīng)的初始體積分別為1.6、1.81、1.9及2.02。從圖6中可以看出，當(dāng)前體積位于參考?jí)嚎s線上時(shí)，壓縮曲線斜率大于參考?jí)嚎s線并向參考?jí)嚎s線趨近。當(dāng)前體積位于參考?jí)嚎s下之下時(shí)，曲線相對(duì)于參考?jí)嚎s線較為平緩，當(dāng)前體積與參考?jí)嚎s線的豎直距離越大則曲線平緩程度越高。模型預(yù)測(cè)的不同體積所對(duì)應(yīng)的壓縮線最終都匯聚于參考?jí)嚎s線。從圖中可以看出。預(yù)測(cè)體積變化趨勢(shì)與試驗(yàn)結(jié)果也較為吻合，說(shuō)明模型能夠很好地描述不同體積和壓力下砂土的壓縮特性。

圖6 Toyoura砂壓縮試驗(yàn)結(jié)果與預(yù)測(cè)對(duì)比[22]

圖7為Miura 等[29]開展的Toyoura砂的三軸壓縮試驗(yàn)。試驗(yàn)砂土比重為2.64，最大顆粒直徑為0.25 mm，不均勻系數(shù)為1.5，最大孔隙比為0.92，最小孔隙比為0.58，試驗(yàn)試樣初始孔隙比與體積分別為0.83和1.83。整個(gè)試驗(yàn)進(jìn)行了3次加載和2次卸載過(guò)程：首先由300 kPa的起始?jí)毫虞d至10 MPa(體積約為1.728)，然后卸載至200 kPa(體積約為1.77)；繼續(xù)加載至17 MPa(體積約為1.666)，再次卸載至488 kPa(體積約為1.71)；最后再次加載至20 MPa。本文預(yù)測(cè)了3次加載過(guò)程，從圖中可以看出由于初始體積位于參考?jí)嚎s線之下，加載時(shí)壓縮線斜率相對(duì)較為平緩，當(dāng)前體積距參考?jí)嚎s線越大，則曲線越平緩。圖7顯示模型預(yù)測(cè)結(jié)果與砂土壓縮試驗(yàn)結(jié)果較為一致，同樣很好地反映了砂土的壓縮特性。

圖7 Toyoura砂加卸載試驗(yàn)結(jié)果與模型預(yù)測(cè)對(duì)比[29]

現(xiàn)有彈性模量試驗(yàn)手段除壓縮試驗(yàn)外，還有彎曲元試驗(yàn)、共振柱試驗(yàn)及現(xiàn)場(chǎng)剪切波速試驗(yàn)，這些方法均屬于小變形范疇。而本文所采用的砂土模量是基于中等變形范疇的壓縮試驗(yàn)規(guī)律得到的。由于砂土的彈性模量與其應(yīng)變水平相關(guān)，應(yīng)變?cè)叫椥阅Ａ吭酱?，并且彈性模量隨著應(yīng)變水平的增加呈現(xiàn)規(guī)律變化。為驗(yàn)證體積對(duì)彈性模量變化趨勢(shì)影響，分別采用應(yīng)變因子FRC(γ)=3和FBE(γ)=3.85對(duì)共振柱試驗(yàn)(RC)與彎曲元(BE)試驗(yàn)數(shù)據(jù)進(jìn)行歸一化處理。圖8為預(yù)測(cè)Toyoura砂的剪切模量與彎曲元及共振柱試驗(yàn)數(shù)據(jù)[30]對(duì)比。試驗(yàn)數(shù)據(jù)采用圍壓為100、200和500 kPa時(shí)Toyoura砂的共振柱和彎曲元試驗(yàn)結(jié)果，通過(guò)試驗(yàn)數(shù)據(jù)可以看出在相同圍壓下剪切模量隨著體積的減小呈非線性增加，而預(yù)測(cè)的剪切模量變化趨勢(shì)與實(shí)測(cè)模量變化趨勢(shì)較為一致，說(shuō)明本文提出的壓縮模型同樣能夠較為可靠地預(yù)測(cè)砂土彈性變形趨勢(shì)。

圖8 Toyoura砂土剪切模量與體積變化關(guān)系[30]

通過(guò)分析模型特點(diǎn)以及模型預(yù)測(cè)與試驗(yàn)對(duì)比結(jié)果，基于體積的對(duì)數(shù)與壓力的冪次關(guān)系構(gòu)建的砂土壓縮模型，不僅表現(xiàn)了砂土在低應(yīng)力段的小變形特征及高壓下的近似線性關(guān)系，而且能夠表現(xiàn)出極限體積附近的低壓縮性特征，能夠很好地描述砂土全壓力下的壓縮特性。與現(xiàn)有單對(duì)數(shù)函數(shù)描述的砂土壓縮曲線相比，模型在低壓下的變形特征及極限體積附近的變形特征描均具有明顯優(yōu)點(diǎn)。同時(shí)，在該關(guān)系函數(shù)框架下建立的彈性模量表達(dá)式能夠很好地描述體積彈性模量隨體積變化趨勢(shì)。

相對(duì)來(lái)說(shuō)，利用CSL線來(lái)確定參考?jí)嚎s線進(jìn)而描述砂土等向壓縮特性的工作稍顯復(fù)雜，且砂土CSL線的建立需要一定數(shù)量的高質(zhì)量三軸剪切試驗(yàn)數(shù)據(jù)。但從參考?jí)嚎s線的建立過(guò)程可以看出，本文已經(jīng)建立CSL線與參考?jí)嚎s線的相關(guān)關(guān)系，這是建立砂土二維本構(gòu)模型最為重要的工作之一，因此建立的CSL線與參考?jí)嚎s線的相關(guān)關(guān)系不僅適用于砂土等向壓縮描述，而且能夠?yàn)樯巴炼S模型建立提供理論依據(jù)。從長(zhǎng)遠(yuǎn)角度，這種復(fù)雜的工作還是值得去探索的。

3 結(jié)論

砂土等向壓縮線不唯一且在e-lgp空間內(nèi)非線性變化及為砂土的數(shù)學(xué)描述帶來(lái)很大的困難，為解決這一問(wèn)題，本文基于臨界狀態(tài)提出了砂土參考?jí)嚎s線的數(shù)學(xué)函數(shù)，主要研究成果可以概括為以下幾個(gè)方面：

(1) 通過(guò)分析砂土等向壓縮特性及力學(xué)特性，建立了與臨界狀態(tài)線相關(guān)的砂土的參考?jí)嚎s線，并采用ln(V-Vmin)-(p/pa)ξ空間內(nèi)線性函數(shù)進(jìn)行描述。

(2) 利用基于當(dāng)前體積與參考?jí)嚎s線的狀態(tài)參量建立修正函數(shù)以獲取其它體積下的壓縮曲線斜率；修正函數(shù)僅對(duì)砂土塑性變形進(jìn)行修正以保證數(shù)學(xué)計(jì)算穩(wěn)定性。

(3) 構(gòu)建的壓縮模型共6個(gè)參數(shù)且均能夠通過(guò)物理試驗(yàn)或簡(jiǎn)單數(shù)學(xué)擬合獲得。

(4) 基于Toyoura砂土的等向壓縮試驗(yàn)、循環(huán)壓縮試驗(yàn)與模型預(yù)測(cè)結(jié)果，提出的砂土壓縮模型能夠很好地反映砂土的壓縮特性；基于模型預(yù)測(cè)的剪切模量與彎曲元和共振柱試驗(yàn)結(jié)果對(duì)比，模型同樣能夠較為可靠地預(yù)測(cè)砂土彈性變形規(guī)律。