張 帥，程曉輝，王天麟

(清華大學(xué)土木工程系，北京 100084)

作為典型的顆粒材料，砂土具有顯著的壓硬性，其剛度表現(xiàn)受到實(shí)際應(yīng)力水平的影響。土體剛度的這種特性在工程應(yīng)用中需要得到充分的重視，尤其是在土體應(yīng)力有較大變化時(shí)的應(yīng)用場(chǎng)景中，例如基坑的開(kāi)挖/回填[1-2]、土石壩的蓄水/泄水[3]、地震波沿深度方向的傳遞[4-5]和地下結(jié)構(gòu)抗震[6]等?；娱_(kāi)挖相關(guān)的研究表明，如果不考慮土體剛度隨深度(即應(yīng)力水平)的變化，將大幅度低估擋土墻的側(cè)向變形[1]。膨脹土地基中的樁-土共同作用研究發(fā)現(xiàn)，土體剛度對(duì)于樁身軸力和樁頭位移都有很大影響[7]。

室內(nèi)實(shí)驗(yàn)(波速法、共振柱、扭剪和三軸等)研究表明，土體在極小應(yīng)變(小于10-5)條件下處于準(zhǔn)彈性狀態(tài)，其模量為常數(shù)[8-10]?？捎洿藭r(shí)的剪切模量為極小應(yīng)變剪切模量G0或最大剪切模量Gmax。土體的極小應(yīng)變模量是深入研究小應(yīng)變模量或者動(dòng)模量的前提，具有極為重要的科研及應(yīng)用價(jià)值。

土體的極小應(yīng)變模量受到許多因素的影響，主要有固結(jié)應(yīng)力水平、密實(shí)度、顆粒形狀、粒徑級(jí)配、微觀結(jié)構(gòu)等。對(duì)于同一種砂土而言，主要是密實(shí)度和固結(jié)應(yīng)力水平的影響，最常用的回歸公式是Hardin-Richart公式[11]：

式中，F(xiàn)(e)是孔隙率e的函數(shù)，主要反映了堆積密度對(duì)于土體剛度的影響，常取為(eg-e)2/(1+e)。A、m和eg是土體參數(shù)，綜合反映了顆粒形狀、粒徑級(jí)配、微觀結(jié)構(gòu)等因素的影響。(p′)m主要反映了土體剛度隨應(yīng)力水平而非線性增長(zhǎng)的特征。理想彈性球的經(jīng)典赫茲接觸給出G0～(p′)13的關(guān)系，但是天然土體顆粒之間的接觸與理想彈性球體相互間的接觸有一定區(qū)別，土力學(xué)實(shí)驗(yàn)得出的m值大多分布在1/3~3/5的區(qū)間內(nèi)[12-13]。

實(shí)驗(yàn)室實(shí)驗(yàn)中的相關(guān)研究大多是基于等向固結(jié)狀態(tài)，而實(shí)際地基大都處于非等向固結(jié)狀態(tài)。因此，基于實(shí)驗(yàn)室實(shí)驗(yàn)的Hardin-Richart公式直接應(yīng)用于實(shí)際工程中，會(huì)帶來(lái)一定的誤差。許多實(shí)驗(yàn)研究表明，非等向固結(jié)下的土體剛度并不符合Hardin-Richart公式的預(yù)測(cè)，從而需要增加一定的修正系數(shù)。

Hardin-Richart等經(jīng)驗(yàn)公式屬于Cauchy彈性模型，在循環(huán)荷載等復(fù)雜應(yīng)力路徑下可能表現(xiàn)出能量非保守特性，即彈性能量是與應(yīng)力路徑相關(guān)的，這在理論上不正確，同時(shí)也會(huì)引起不正確的殘余變形[14]。而彈性力學(xué)中的超彈性模型也稱(chēng)為Green彈性模型，其應(yīng)力-應(yīng)變關(guān)系始終從一個(gè)自由能函數(shù)的微分形式給出，完全滿(mǎn)足熱力學(xué)基本定律的要求，在不同材料的本構(gòu)關(guān)系研究中得到廣泛應(yīng)用[15-16]。另外，超彈性模型在處理正剪耦合時(shí)比較方便，例如剪脹性、應(yīng)力引起的各向異性、主應(yīng)力旋轉(zhuǎn)等[17-18]。因此研究非等向固結(jié)，尤其是高度非等向狀態(tài)下的土體剛度，采用超彈性模型更為合適。Houlsby等研究了能反映土體剛度隨平均應(yīng)力變化的超彈性模型[16]，但非等向固結(jié)狀態(tài)的影響并沒(méi)有得到討論。

本文首先對(duì)以往學(xué)者完成的非等向固結(jié)狀態(tài)下的砂土極小應(yīng)變剛度實(shí)驗(yàn)進(jìn)行研究分析，得出有較多實(shí)驗(yàn)數(shù)據(jù)支持的經(jīng)驗(yàn)規(guī)律；然后針對(duì)3種能夠反映土體剛度隨應(yīng)力水平增大規(guī)律的超彈性模型進(jìn)行研究，發(fā)現(xiàn)其中僅有HE1模型對(duì)砂土非等向固結(jié)狀態(tài)下的剛度規(guī)律趨勢(shì)預(yù)測(cè)正確；最后，對(duì)于HE1模型進(jìn)行了參數(shù)影響分析。

1 非等向固結(jié)實(shí)驗(yàn)研究分析

關(guān)于非等向固結(jié)狀態(tài)的研究，一般是考慮常規(guī)三軸固結(jié)狀態(tài)，即記此時(shí)的固結(jié)應(yīng)力比非等向固結(jié)狀態(tài)下的土體極小應(yīng)變剪切模量，可以由下列形式的修正Hardin-Richart公式給出：

其中：R為修正系數(shù)，取值由k13決定。實(shí)驗(yàn)研究發(fā)現(xiàn)，土體剪切剛度主要受相應(yīng)平面內(nèi)的兩個(gè)主應(yīng)力的影響[19]，因此Yu和Richart[20]建議采用來(lái)代替p′對(duì)剛度做歸一化：

并基于Ottawa砂和Toyoura砂實(shí)驗(yàn)數(shù)據(jù)給出了經(jīng)驗(yàn)公式：

在使用p′進(jìn)行歸一化時(shí)，Yu-Richart公式(式(3)和式(4))換算為式(2)后給出的修正系數(shù)R隨著固結(jié)應(yīng)力比k13的增大先是略微增大，而后較大幅度下降?？偟膩?lái)說(shuō)，當(dāng)固結(jié)應(yīng)力比小于某閾值時(shí)，修正系數(shù)R基本為1；當(dāng)固結(jié)應(yīng)力比相對(duì)較大時(shí)，修正系數(shù)R迅速降低，降低幅值可超過(guò)20%。這與Kuribayashi等關(guān)于Toyoura砂[21](圖中三角符號(hào))、Tatsuoka等關(guān)于Toyoura砂[22](圖中未列出)以及Bellotti等關(guān)于Ticino砂[23](圖中方塊符號(hào))的實(shí)驗(yàn)研究結(jié)果相一致。另外，DEM數(shù)值模擬[24]、極小應(yīng)變條件下的彈性模量的實(shí)驗(yàn)研究[25]和循環(huán)荷載條件下的剪切模量的實(shí)驗(yàn)研究[26]同樣得出了類(lèi)似的結(jié)論。

表1 前人實(shí)驗(yàn)總結(jié)Table 1 Summary of previous experiments

圖1 不同固結(jié)應(yīng)力比下的土體極小應(yīng)變模量修正系數(shù)Fig.1 Comparison of the relationships between shear modules ratio and stress ratio by various investigations

雖然也有實(shí)驗(yàn)研究認(rèn)為修正系數(shù)R隨著k13的增大而增大[27-29]，如Payan等[27]關(guān)于Blue砂的實(shí)驗(yàn) (圖中圓形符號(hào))，但這些研究的k13值范圍普遍較小，一般小于2.5。鑒于Yu-Richart的研究系統(tǒng)性更強(qiáng)，且得到多數(shù)相關(guān)研究的重復(fù)驗(yàn)證，因此我們主要采用其研究結(jié)果來(lái)作為參考?？紤]到Y(jié)u-Richart建議了新的歸一化方式，因此，本文中兩種歸一化方式的結(jié)果都予以呈現(xiàn)。

2 超彈性模型

2.1 超彈性模型簡(jiǎn)介

本文所有應(yīng)力均指有效應(yīng)力，所有應(yīng)變均指彈性應(yīng)變(極小應(yīng)變條件下，彈性應(yīng)變等于總應(yīng)變)。σij為Cauchy有效應(yīng)力張量，εij為彈性應(yīng)變張量；δij為Kronecker delta記號(hào)，上述2階張量均遵循啞標(biāo)的愛(ài)因斯坦求和約定。應(yīng)力不變量有平均應(yīng)力及剪應(yīng)力其中為偏應(yīng)力張量。同樣，應(yīng)變不變量有體積應(yīng)變及剪應(yīng)變其中為偏應(yīng)變張量。在常規(guī)三軸條件下，一般采用如下形式的簡(jiǎn)寫(xiě)；此時(shí)的應(yīng)力和應(yīng)變不變量有

超彈性模型是指材料具有Helmholtz自由能F(又可稱(chēng)為彈性勢(shì)能)或Gibbs自由能E(又可稱(chēng)為應(yīng)變余能)，其應(yīng)力-應(yīng)變關(guān)系可以從相應(yīng)的自由能函數(shù)的微分中推導(dǎo)得出。這兩種形式的自由能可以通過(guò)勒讓德變換來(lái)互相推導(dǎo)，但是對(duì)于某些函數(shù)形式的自由能，可能并不存在簡(jiǎn)單函數(shù)形式的另一種對(duì)應(yīng)的自由能函數(shù)形式。

對(duì)于Helmholtz自由能形式的超彈性模型有：

進(jìn)而得出其應(yīng)力-應(yīng)變關(guān)系為：

對(duì)于Gibbs自由能形式的超彈性模型有：

進(jìn)而得出其應(yīng)力-應(yīng)變關(guān)系為：

考慮Hardin-Richart公式給出的剛度隨平均應(yīng)力的冪律關(guān)系，其對(duì)應(yīng)的超彈性模型大致有3種[16]，如表2所示。其中應(yīng)力和應(yīng)變不變量的定義及其與三軸狀態(tài)下的不變量的關(guān)系參見(jiàn)2.1節(jié)；參數(shù)n反映了土體剛度隨應(yīng)力水平的增長(zhǎng)，與Hardin-Richart公式中的m有關(guān)，具體關(guān)系2.1.1節(jié)、2.1.2節(jié)、2.1.3節(jié)將推導(dǎo)給出；ξ或?與泊松比有關(guān)；B或C影響著土的剛度模量的大小。

表2 3種能反映Hardin-Richart公式的超彈性模型Table 2 Three types of hyperelastic models compatible to the Hardin-Richart formula

2.1.1 HE1模型

首先得出HE1超彈性模型的應(yīng)力不變量為：

進(jìn)而可以得出其切向剛度矩陣，其中：

如果考慮等向固結(jié)狀態(tài)下：

將土體剪切模量中的應(yīng)變項(xiàng)消去，得出：

與式(2)對(duì)比，可以得出HE1模型給出的Hardin-Richart公式修正系數(shù)為：

和式(14)共同繪圖可以得出其隨固結(jié)應(yīng)力比變化的規(guī)律。

2.1.2 HE2模型

對(duì)于HE2超彈性模型，其兩種自由能形式給出的結(jié)果是等價(jià)的，在此僅給出Helmholtz自由能形式的結(jié)果。首先得出其應(yīng)力不變量為：

進(jìn)而得出其切向剛度矩陣，其中：

在等向固結(jié)狀態(tài)下同樣可以得出與HE1類(lèi)似的符合經(jīng)典Hardin-Richart公式的結(jié)果，在此不再重復(fù)。

2.1.3 HE3模型

首先得出HE3超彈性模型的應(yīng)變不變量為

進(jìn)而得出其切向柔度矩陣，求逆可得出其切向剛度矩陣，其中

式中：-n相當(dāng)于Hardin-Richart公式中的m，因此n的取值一般為-1/3~3/5，即小于0?？紤]到體應(yīng)變不能為負(fù)，則由式(20)可知即可以看出HE3模型對(duì)于土體能達(dá)到的最大剪切應(yīng)力比也有一個(gè)限值。

2.2 模擬對(duì)比分析

式(16)、式(19)、式(21)中給出了修正Hardin-Richart公式中的修正系數(shù)R，將其隨固結(jié)應(yīng)力比k13的變化繪制出來(lái)如圖2所示。圖中同時(shí)給出了經(jīng)典Hardin-Richart公式以及Yu-Richart經(jīng)驗(yàn)公式的結(jié)果。

HE1模型給出的修正系數(shù)R隨著固結(jié)應(yīng)力比k13的增大而衰減，與Yu-Richart經(jīng)驗(yàn)公式給出的趨勢(shì)較為一致；而HE2和HE3模型給出的R都隨著固結(jié)應(yīng)力比k13的增大而持續(xù)增大，與實(shí)驗(yàn)規(guī)律不一致。需要說(shuō)明的是，圖示曲線具體數(shù)值仍然受到模型參數(shù)的影響，但并不改變這種明顯不同的趨勢(shì)。從模型本身角度而言，這些不同的趨勢(shì)主要受到兩個(gè)因素的影響：一個(gè)是彈性剪脹或剪縮特性，另一個(gè)是剪切剛度的影響因素。從式(10)可以看出，HE1模型具有彈性剪脹特性，即保持平均應(yīng)力不變條件下剪切，能使得彈性體積應(yīng)變減??；而式(11)表明其剪切剛度僅受彈性體積應(yīng)變影響，因而非等向固結(jié)這種剪切能夠使得剪切剛度變小。式(17)顯示HE2模型同樣具有彈性剪脹特性，但式(18)顯示其剪切剛度同時(shí)受到彈性體積應(yīng)變和彈性剪切應(yīng)變的影響，非等向固結(jié)狀態(tài)下這兩者共同作用使得剪切剛度增大。而式(20)顯示HE3模型則具有彈性剪縮特性，不能反映非等向固結(jié)狀態(tài)下土體剛度的衰減規(guī)律。

圖2 不同固結(jié)應(yīng)力比下的土體極小應(yīng)變模量修正系數(shù)Fig.2 The relationships between shear modules ratio and stress ratio

考慮到HE2和HE3預(yù)測(cè)的趨勢(shì)與實(shí)驗(yàn)結(jié)果有較大差異，因此我們僅對(duì)HE1進(jìn)行定量的參數(shù)分析。在做此分析時(shí)，主要對(duì)比一下HE1的預(yù)測(cè)與Yu-Richart公式的結(jié)果，因此在繪圖時(shí)采用Yu-Richart的建議：縱坐標(biāo)剛度的歸一化時(shí)采用即為式(3)中的R*；橫坐標(biāo)固結(jié)應(yīng)力比的歸一化時(shí)采用式(4)中的kn。HE1模型對(duì)土體能達(dá)到的最大固結(jié)應(yīng)力比有一個(gè)限值，我們采用這個(gè)限值對(duì)其固結(jié)應(yīng)力比進(jìn)行歸一化，從而得到HE1模型中相應(yīng)的kn。

HE1模型的3個(gè)參數(shù)中：B只影響土體剛度的大小，不影響該相對(duì)規(guī)律；n的具體數(shù)值可通過(guò)實(shí)驗(yàn)數(shù)據(jù)很精確地標(biāo)定，且對(duì)于表1所總結(jié)的砂土都大致為1；ξ的具體數(shù)值難以通過(guò)常規(guī)土力學(xué)實(shí)驗(yàn)來(lái)精確標(biāo)定，所以有必要進(jìn)行參數(shù)分析。將Yu-Richart公式與HE1模型的結(jié)果對(duì)比，如圖3所示。HE1模型的結(jié)果受到其參數(shù)ξ的影響：ξ越大時(shí)，修正系數(shù)R*整體越小。HE1的結(jié)果與Yu-Richart公式總體比較吻合。在較低應(yīng)力比時(shí)，HE1的結(jié)果要略小于Yu-Richart公式；在較高應(yīng)力比時(shí)，HE1的結(jié)果要略大于(ξ=2或ξ=3)或略小于(ξ=4)Yu-Richart公式。

圖3 HE1模型不同固結(jié)應(yīng)力比下的模量修正系數(shù)Fig.3 The relationships between shear modules ratio and stress ratio by HE1 model

3 結(jié)論

(1) 3種類(lèi)型的超彈性模型能夠反映等向固結(jié)狀態(tài)下土體剛度隨應(yīng)力水平增大而增大的冪律關(guān)系，但僅有HE1模型對(duì)土體非等向固結(jié)狀態(tài)下的剛度規(guī)律趨勢(shì)預(yù)測(cè)正確。這主要是因?yàn)槠渚哂邪◤椥约裘浽趦?nèi)的特殊的正剪耦合特性。

(2) HE1模型對(duì)土體非等向固結(jié)狀態(tài)下的剛度規(guī)律預(yù)測(cè)結(jié)果受到其參數(shù)ξ取值的影響，但總體而言與實(shí)驗(yàn)結(jié)果總結(jié)的規(guī)律一致。

(3) 下一步工作，考慮具有非線性剛度地基與地下結(jié)構(gòu)的彈性動(dòng)力相互作用，以便反映在邊界條件下的地基非等向固結(jié)極小應(yīng)變剛度對(duì)地下結(jié)構(gòu)抗震性能的影響。