張 帥，程曉輝，王天麟

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

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

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

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

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

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

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

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

1 非等向固結(jié)實驗研究分析

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

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

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

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

表1 前人實驗總結(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ù)R隨著k13的增大而增大[27-29]，如Payan等[27]關(guān)于Blue砂的實驗 (圖中圓形符號)，但這些研究的k13值范圍普遍較小，一般小于2.5。鑒于Yu-Richart的研究系統(tǒng)性更強，且得到多數(shù)相關(guān)研究的重復(fù)驗證，因此我們主要采用其研究結(jié)果來作為參考?？紤]到Y(jié)u-Richart建議了新的歸一化方式，因此，本文中兩種歸一化方式的結(jié)果都予以呈現(xiàn)。

2 超彈性模型

2.1 超彈性模型簡介

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

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

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

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

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

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

考慮Hardin-Richart公式給出的剛度隨平均應(yīng)力的冪律關(guān)系，其對應(yīng)的超彈性模型大致有3種[16]，如表2所示。其中應(yīng)力和應(yīng)變不變量的定義及其與三軸狀態(tài)下的不變量的關(guān)系參見2.1節(jié)；參數(shù)n反映了土體剛度隨應(yī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)力不變量為：

進而可以得出其切向剛度矩陣，其中：

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

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

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

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

2.1.2 HE2模型

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

進而得出其切向剛度矩陣，其中：

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

2.1.3 HE3模型

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

進而得出其切向柔度矩陣，求逆可得出其切向剛度矩陣，其中

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

2.2 模擬對比分析

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

HE1模型給出的修正系數(shù)R隨著固結(jié)應(yīng)力比k13的增大而衰減，與Yu-Richart經(jīng)驗公式給出的趨勢較為一致；而HE2和HE3模型給出的R都隨著固結(jié)應(yīng)力比k13的增大而持續(xù)增大，與實驗規(guī)律不一致。需要說明的是，圖示曲線具體數(shù)值仍然受到模型參數(shù)的影響，但并不改變這種明顯不同的趨勢。從模型本身角度而言，這些不同的趨勢主要受到兩個因素的影響：一個是彈性剪脹或剪縮特性，另一個是剪切剛度的影響因素。從式(10)可以看出，HE1模型具有彈性剪脹特性，即保持平均應(yīng)力不變條件下剪切，能使得彈性體積應(yīng)變減小；而式(11)表明其剪切剛度僅受彈性體積應(yīng)變影響，因而非等向固結(jié)這種剪切能夠使得剪切剛度變小。式(17)顯示HE2模型同樣具有彈性剪脹特性，但式(18)顯示其剪切剛度同時受到彈性體積應(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ù)測的趨勢與實驗結(jié)果有較大差異，因此我們僅對HE1進行定量的參數(shù)分析。在做此分析時，主要對比一下HE1的預(yù)測與Yu-Richart公式的結(jié)果，因此在繪圖時采用Yu-Richart的建議：縱坐標剛度的歸一化時采用即為式(3)中的R*；橫坐標固結(jié)應(yīng)力比的歸一化時采用式(4)中的kn。HE1模型對土體能達到的最大固結(jié)應(yīng)力比有一個限值，我們采用這個限值對其固結(jié)應(yīng)力比進行歸一化，從而得到HE1模型中相應(yīng)的kn。

HE1模型的3個參數(shù)中：B只影響土體剛度的大小，不影響該相對規(guī)律；n的具體數(shù)值可通過實驗數(shù)據(jù)很精確地標定，且對于表1所總結(jié)的砂土都大致為1；ξ的具體數(shù)值難以通過常規(guī)土力學(xué)實驗來精確標定，所以有必要進行參數(shù)分析。將Yu-Richart公式與HE1模型的結(jié)果對比，如圖3所示。HE1模型的結(jié)果受到其參數(shù)ξ的影響：ξ越大時，修正系數(shù)R*整體越小。HE1的結(jié)果與Yu-Richart公式總體比較吻合。在較低應(yīng)力比時，HE1的結(jié)果要略小于Yu-Richart公式；在較高應(yīng)力比時，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種類型的超彈性模型能夠反映等向固結(jié)狀態(tài)下土體剛度隨應(yīng)力水平增大而增大的冪律關(guān)系，但僅有HE1模型對土體非等向固結(jié)狀態(tài)下的剛度規(guī)律趨勢預(yù)測正確。這主要是因為其具有包括彈性剪脹在內(nèi)的特殊的正剪耦合特性。

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

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