天津濱海軟黏土蠕變特性及其經(jīng)驗?zāi)Ｐ脱芯?/h1>
2014-05-17 00:57:50王元戰(zhàn)董焱赫
水道港口 2014年3期
王元戰(zhàn)，董焱赫
（天津大學(xué) 建筑工程學(xué)院，天津市港口與海洋工程重點實驗室，天津 300072）
天津濱海軟黏土蠕變特性及其經(jīng)驗?zāi)Ｐ脱芯?/p>
王元戰(zhàn)，董焱赫
（天津大學(xué) 建筑工程學(xué)院，天津市港口與海洋工程重點實驗室，天津 300072）
在室內(nèi)三軸蠕變試驗的基礎(chǔ)上，分析了天津濱海地區(qū)重塑飽和軟黏土的非線性蠕變特性。蠕變試驗采用分級加載的方式進(jìn)行，利用“陳氏法”對實測數(shù)據(jù)進(jìn)行處理，得到飽和軟黏土在不同應(yīng)力狀態(tài)下的蠕變試驗曲線。在此基礎(chǔ)上，依據(jù)黏土的應(yīng)力-應(yīng)變、應(yīng)變-時間關(guān)系，分別用Singh?Mitchell和Mesri經(jīng)驗?zāi)Ｐ蛿M合試驗數(shù)據(jù)?？偨Y(jié)擬合結(jié)果，結(jié)合濱海軟黏土的蠕變特性，提出一種分段擬合的Mesri模型，在蠕變前期，模型的應(yīng)力-應(yīng)變關(guān)系采用雙曲線函數(shù)，應(yīng)變-時間關(guān)系采用冪函數(shù)；在蠕變后期，模型的應(yīng)力-應(yīng)變關(guān)系采用雙曲線函數(shù)，應(yīng)變-時間關(guān)系采用雙曲線函數(shù)，以兩條擬合曲線的交點作為分段擬合的分界點。對比分析表明，該方法能更好地描述具有衰減穩(wěn)定特性的蠕變曲線。選取天津濱海地區(qū)原狀土蠕變試驗數(shù)據(jù)對該模型進(jìn)行驗證，擬合結(jié)果較好，由此建立了適用于具有衰減穩(wěn)定蠕變特性的濱海軟黏土非線性蠕變模型。
軟黏土；蠕變；衰減穩(wěn)定；經(jīng)驗?zāi)Ｐ?；天津濱海
Biography：WANG Yuan?zhan（1958-），male，professor.
在中國經(jīng)濟(jì)發(fā)展活躍的長江口、珠江口、天津濱海等海岸地區(qū)，廣泛分布著軟黏土地基，這些地區(qū)經(jīng)濟(jì)發(fā)展迅速，在這些軟土地基上興建了大量的建筑物，以天津港為例，港口規(guī)模不斷擴(kuò)大，吞吐量已經(jīng)躍居世界第四，大量的陸域設(shè)施和圍海造陸工程正在建設(shè)中。但由于該種土類流變效應(yīng)較為明顯，靜載條件下變形緩慢增加，即發(fā)生蠕變效應(yīng)，使得這些建在工程病害層上的建筑物存在著不同程度的沉降，使之有失穩(wěn)、滑塌的危險。如近年天津北疆高樁碼頭多個轉(zhuǎn)角處由于軟黏土岸坡土體的蠕變發(fā)生了基樁連同樁帽與上部結(jié)構(gòu)之間出現(xiàn)較大相對變位的現(xiàn)象，導(dǎo)致橫梁的擱置長度不足，有的只剩下了幾厘米，使得碼頭結(jié)構(gòu)處于危險狀態(tài)之中［1］。開展濱海軟黏土的蠕變試驗研究，建立其合適的蠕變模型，是亟待開展的工作，對研究濱海地區(qū)港口碼頭和防波堤結(jié)構(gòu)的長期穩(wěn)定性有著重要意義。關(guān)于長期荷載下軟黏土的蠕變特性的問題，學(xué)者們展開了大量的研究。
在國外，Zhen［2］等在Perzyna超限應(yīng)力理論的基礎(chǔ)上，建立了一個各向異性的黏彈塑性軟黏土模型，并針對法國St.Herblain地區(qū)的軟土進(jìn)行了等應(yīng)變率蠕變試驗，對所建立的模型進(jìn)行了驗證；Leoni［3-5］等和Ver?meer等考慮土體各向異性，建立軟土蠕變模型，并對軟土的蠕變行為進(jìn)行數(shù)值模擬研究。Kelln［6］等嘗試采用粘彈力學(xué)原理揭示軟粘土的蠕變行為。在國內(nèi)，李軍世［7］等采用Singh?Mitchell應(yīng)變速率-時間關(guān)系方程描述土體的蠕變特性，結(jié)合上海淤泥質(zhì)黏土室內(nèi)試驗，給出了上海淤泥質(zhì)黏土的蠕變形式；王常明［8］等的三軸蠕變試驗研究表明，濱海相沉積的軟土具有非線性蠕變的特性，雙曲線模型更適合其應(yīng)力-應(yīng)變關(guān)系，修改了Singh?Mitchell模型的應(yīng)力-應(yīng)變關(guān)系，建議了濱海軟土應(yīng)力-應(yīng)變-時間關(guān)系模型。王元戰(zhàn)［9］等基于屈服理論，研究了天津濱海淤泥質(zhì)軟黏土的蠕變變形規(guī)律，建立了反映濱海軟黏土的非線性流變本構(gòu)模型；閆澍旺［10］等，采用大直徑原狀試樣進(jìn)行了三軸不固結(jié)不排水蠕變試驗，研究了天津濱海新區(qū)吹填黏土的蠕變變形的規(guī)律性，在此基礎(chǔ)上建立天津濱海新區(qū)軟黏土的無屈服面蠕變模型；王元戰(zhàn)［11］等采用Merchant蠕變模型三軸蠕變試驗的結(jié)果曲線進(jìn)行擬合，建立了適用于天津濱海地區(qū)的軟黏土本構(gòu)模型。
雖然目前有很多關(guān)于土體蠕變模型的研究成果，也不乏對天津濱海軟黏土蠕變特性的研究，但是這些成果不足以完全詮釋天津濱海地區(qū)軟黏土的蠕變特性，原因有2個：一是蠕變模型自身復(fù)雜性和局限性。二是以往研究得出的模型大多只能描述自身的試驗數(shù)據(jù)，沒有用其他試驗數(shù)據(jù)加以驗證，缺乏綜合比較與分析。
本文以天津港病害層重塑淤泥質(zhì)黏土為例進(jìn)行蠕變試驗研究。依據(jù)試驗數(shù)據(jù)，結(jié)合Singh?Mitchell和Mesri經(jīng)驗蠕變模型，提出一種分段擬合的Mesri模型，在蠕變前期，模型的應(yīng)力-應(yīng)變關(guān)系采用雙曲線函數(shù)，應(yīng)變-時間關(guān)系采用冪函數(shù)；在蠕變后期，模型的應(yīng)力-應(yīng)變關(guān)系采用雙曲線函數(shù)，應(yīng)變-時間關(guān)系采用雙曲線函數(shù)，以兩條擬合曲線的交點作為分段擬合的分界點。對比分析表明，該模型能更好地描述具有衰減穩(wěn)定特性的蠕變曲線。選取以往天津濱海地區(qū)兩種典型軟黏土蠕變試驗數(shù)據(jù)對該模型進(jìn)行驗證，擬合結(jié)果較好，由此建立了適用于具有衰減穩(wěn)定蠕變特性的濱海軟黏土非線性蠕變模型。
1 三軸蠕變試驗
1.1 試驗設(shè)備及試樣
設(shè)備采用的是無級加載應(yīng)力控制式三軸蠕變試驗儀，該三軸蠕變試驗儀由普通應(yīng)變式三軸儀通過改裝軸向加載設(shè)備而成，儀器主要由軸向加載設(shè)備、圍壓控制測量系統(tǒng)、反壓控制測量系統(tǒng)、體積變化和孔隙水壓力測量系統(tǒng)等組成，將軸向加載設(shè)備由原來的應(yīng)變控制式改造成應(yīng)力控制式加載。這樣改造試驗儀器簡單易行，且既可以利用原有應(yīng)變式三軸儀器設(shè)備優(yōu)良穩(wěn)定的圍壓控制測量系統(tǒng)、反壓控制測量系統(tǒng)、體積變化和孔隙水壓力測量系統(tǒng)，又可以通過受力平衡計算，利用干砂和砝碼施加各級軸向荷載，獲得恒定且準(zhǔn)確的軸向偏應(yīng)力。
黏土試樣為天津港地區(qū)5～10 m土層飽和重塑軟黏土，土樣呈流塑狀態(tài)。依據(jù)土工試驗規(guī)程（SL237-1999）制作重塑試樣。試樣直徑39.1 mm，高為80 mm。土樣的基本指標(biāo)如表1所示。
1.2 試驗方法及加載方案
在土的室內(nèi)流變試驗中，對于蠕變試驗有兩種不同的加載方式，即分別加載法和分級加載法。所謂分級加載，就是在同一試樣上逐級加上不同的應(yīng)力，即在一級應(yīng)力水平下蠕變經(jīng)歷給定的時間或達(dá)到穩(wěn)定后，將應(yīng)力水平提高到下一級的水平，直到所需的應(yīng)力水平［12］。為了克服分別加載的局限，減小采用不同土樣進(jìn)行試驗對試驗結(jié)果帶來的離散性，本試驗加載方式采用分級加載法。
由于軟黏土具有滲透性差的特點，試驗采用不排水蠕變試驗，圍壓分別取50 kPa，100 kPa，150 kPa。在蠕變試驗進(jìn)行前，應(yīng)選取同一批次重塑土樣進(jìn)行圍壓50 kPa、100 kPa和150 kPa下的固結(jié)不排水剪切試驗，以試樣應(yīng)變達(dá)到15%為破壞標(biāo)準(zhǔn)，確定破壞偏應(yīng)力qf分別為70.2 kPa，120.5 kPa和178.8 kPa。根據(jù)qf確定每級加載偏應(yīng)力，取加載級數(shù)為5級，加載方案如表2所示，再通過受力平衡計算各級加砂質(zhì)量。
表1 土樣基本特征表Tab.1 Basic physical parameters of soil samples
蠕變試驗步驟為：第1步，將制好的試樣放入真空室中，將室內(nèi)抽氣呈負(fù)壓狀態(tài)，10 h以后打開排水閥，使土樣進(jìn)水飽和。第2步，飽和完成后，將試樣安放在三軸壓力室進(jìn)行固結(jié)，由于軟黏土排水性差，故固結(jié)時間設(shè)為24 h。第3步，固結(jié)基本完成以后，設(shè)置量測豎向變形的百分表初始讀數(shù)，關(guān)閉排水閥，再施加第一級偏應(yīng)力ΔF1，同時通過百分表量測軸向變形。在偏應(yīng)力作用下蠕變穩(wěn)定或經(jīng)歷給定的時間（本次試驗每級偏應(yīng)力施加時間確定為48 h）后，再增加配重，使偏壓力達(dá)到ΔF2的水平，以此直到五級荷載施加完畢。
表2 土樣分級加載方案Tab.2 Schemes of step loading kPa
1.3 試驗結(jié)果分析
限于篇幅，以下以圍壓100 kPa下的蠕變試驗為例，進(jìn)行分析。試驗得到的曲線為分級加載下的應(yīng)變-時間曲線，該流變曲線的可用性較小，需要分級進(jìn)行整理，通常的整理方法有線性法和陳氏法。線性法是將土體視為線性流變體，此時土體蠕變響應(yīng)滿足線性疊加原理，蠕變曲線可以簡單地由每一級荷載的蠕變響應(yīng)疊加而成。
若要更準(zhǔn)確的研究土的蠕變，應(yīng)將土體視為非線性流變體。由陳宗基教授提出并由他的學(xué)生所發(fā)展的“陳氏法”通過采用適當(dāng)?shù)膶嶒灱夹g(shù)與方法，用作圖法建立真實變形過程的疊加關(guān)系，在考慮土體為非線性流變體或線性流變體時均可適用，在國內(nèi)外巖土流變學(xué)研究中得到越來越廣泛的應(yīng)用。陳氏法可以使蠕變試驗通過采用適當(dāng)?shù)募虞d程序從單一試樣得到更多的試驗資料，這對試驗周期長，并且試驗價格昂貴的巖土蠕變性質(zhì)研究有重大的實際意義。本試驗對蠕變試驗結(jié)果采用的整理方法就是陳氏法。
用“陳氏法”處理得到不同應(yīng)力狀態(tài)下的軸向應(yīng)變與加載時間的蠕變關(guān)系曲線如圖1所示，由圖1可知：（1）在不排水條件下，當(dāng)主應(yīng)力差達(dá)到破壞偏應(yīng)力前，變形呈現(xiàn)出明顯的衰減穩(wěn)態(tài)蠕變過程。在各級應(yīng)力下，應(yīng)變隨著時間的增長逐漸趨于穩(wěn)定。當(dāng)主應(yīng)力差較低時，應(yīng)變呈現(xiàn)出明顯的衰減蠕變趨勢，蠕變現(xiàn)象不明顯；當(dāng)應(yīng)力水平較高時，試樣蠕變曲線包括衰減階段和穩(wěn)態(tài)蠕變階段，試樣穩(wěn)態(tài)蠕變階段的斜率隨主應(yīng)力差的增加而增加；（2）在不排水條件下，試樣的應(yīng)變隨著偏應(yīng)力水平的提高而越來越大，且增加幅度越來越大。在試樣破壞之前，變形曲線并未出現(xiàn)加速蠕變階段。由蠕變曲線得到的應(yīng)力-應(yīng)變等時曲線如圖2所示：不同時刻，除初始時刻外，各個時刻的應(yīng)力-應(yīng)變曲線為彼此相似的曲線，可以用同一類函數(shù)描述；在同一時刻下，應(yīng)力-應(yīng)變曲線具有明顯屈服特性，曲線呈非線性，這說明，需要用非線性的蠕變模型來描述濱海軟黏土蠕變特性。
圖1 應(yīng)變-時間關(guān)系圖Fig.1 ε?t curves of creep test
圖2 應(yīng)力-應(yīng)變等時曲線Fig.2 σ?ε isochronal curves
2 蠕變經(jīng)驗?zāi)Ｐ?/h2>
蠕變的經(jīng)驗?zāi)Ｐ停菑牧鲃永碚摶蚶匣碚摮霭l(fā)，直接給出土蠕變方程的應(yīng)力-應(yīng)變-時間函數(shù)形式。對每種不同的材料，甚至不同的條件，有多樣的經(jīng)驗?zāi)Ｊ?，和元件模型、屈服面模型、?nèi)時模型等相比，它比較直觀，參數(shù)簡單且容易獲取，亦為工程設(shè)計人員所樂意采用。目前常用的具有代表性的包括Singh?Mitchell蠕變速率和蠕變函數(shù)，常速率軸向應(yīng)變條件下的雙曲線計算模型以及Mesri應(yīng)力-應(yīng)變-時間蠕變函數(shù)等。
2.1 Singh?Mitchell模型
1968年Singh等提出的著名的Singh?Mitchell經(jīng)驗蠕變模型［13-14］是在總結(jié)了單級常應(yīng)力加載，排水與不排水三軸壓縮試驗數(shù)據(jù)的基礎(chǔ)上，采用指數(shù)函數(shù)表示應(yīng)力-應(yīng)變關(guān)系、采用冪函數(shù)表示應(yīng)變-時間關(guān)系的經(jīng)典的經(jīng)驗蠕變模型，該模型能夠很好地反映土的蠕變特性。Singh?Mitchell蠕變模型寫成應(yīng)變率-應(yīng)力水平-時間的關(guān)系式，可表示為
式中：ε˙為任一時刻t的軸向應(yīng)變速率；t為試樣受荷時間；Dr=（σ1-σ3）/（σ1-σ3）f為偏應(yīng)力水平，（σ1-σ3）f為破壞偏應(yīng)力，可由常規(guī)三軸排水壓縮實驗得到；Ar是在單位參考時刻tr，偏應(yīng)力水平為零時的應(yīng)變速率；m由ln˙-lnt擬合直線的斜率確定；α為應(yīng)變速率對數(shù)與剪應(yīng)力關(guān)系圖中線性段的斜率。
對式（1）進(jìn)行積分，在不考慮初始應(yīng)變且m不為1時。經(jīng)推導(dǎo)可得
式中：B=Artr/（1-m）；β=α；λ=1-m。該模型需要確定的3個參數(shù)是B，β，λ。
2.2 Mesri模型
針對Singh?Mitchell模型存在的局限性，Mesri.G［15-16］于1981年提出了新的經(jīng)驗?zāi)Ｐ?。該模型?yīng)力-應(yīng)變關(guān)系采用雙曲線函數(shù)，應(yīng)變-時間關(guān)系仍采用冪函數(shù)。經(jīng)典的Mesri模型公式如下
式中有3個參數(shù)：參考應(yīng)變2/（Eu/Su），擬合比Rf和蠕變參數(shù)λ。
用Singh?Mitchell模型和Mesri模型分別擬合試驗曲線，擬合結(jié)果如圖3和圖4所示。
2.3 分段擬合的Mesri模型
從擬合結(jié)果可以看出，Mesri模型和Singh?Mitchell模型都存在后期擬合不準(zhǔn)確的情況，原因在于兩個模型的應(yīng)變-時間關(guān)系都采用了冪函數(shù)的形式，而冪函數(shù)是非衰減函數(shù)，當(dāng)時間增大時，曲線斜率較大，對應(yīng)的應(yīng)變值遠(yuǎn)大于實測值，對于尚處于衰減穩(wěn)定期的試驗數(shù)據(jù)是不適用的。
2.3.1變形函數(shù)（應(yīng)力-應(yīng)變關(guān)系）的確定
當(dāng)取某一時刻t1時，公式（2）和（3）分別轉(zhuǎn)化為t1時刻的應(yīng)力-應(yīng)變方程。以100 kPa圍壓下的實測數(shù)據(jù)為例，取t1=1 d，分別用兩種函數(shù)去擬合試驗數(shù)據(jù)，通過應(yīng)變對數(shù)-偏應(yīng)力水平曲線（圖5）可以看出，當(dāng)偏應(yīng)力水平處于中間段時，實測數(shù)據(jù)呈線性關(guān)系，而在偏應(yīng)力水平趨近于1或趨近于0時，實測數(shù)據(jù)則呈現(xiàn)出明顯的非線性；指數(shù)函數(shù)關(guān)系只能描述偏應(yīng)力水平處于中間段（20%～80%偏應(yīng)力水平）的應(yīng)力-應(yīng)變關(guān)系，而雙曲函數(shù)能夠準(zhǔn)確地描述全部應(yīng)力水平下的狀態(tài)。故應(yīng)力-應(yīng)變關(guān)系應(yīng)采用雙曲線函數(shù)描述。
圖3 試驗與Sing?Mitchell模型計算結(jié)果比較Fig.3 Comparison between test data and calculated results using S?M model
圖4 試驗數(shù)據(jù)與Mesri模型計算結(jié)果比較Fig.4 Comparison between test data and calculated results using Mesri model
圖5 應(yīng)力-應(yīng)變等時曲線Fig.5 σ?ε isochronal curves
2.3.2 蠕變函數(shù)的確定
結(jié)合上文研究，為了修正Mesri模型的誤差，由C.C.維亞洛夫提出的通用的時間函數(shù)［17］入手，本文建議了一種分段擬合的經(jīng)驗?zāi)Ｐ?，該模型前期仍采用Mesri函數(shù)的形式；在擬合后期數(shù)據(jù)時將Mesri函數(shù)應(yīng)變-時間關(guān)系修正為雙曲線函數(shù)，該函數(shù)的形式為
式中：T為蠕變參數(shù)，2/（Eu/Su）和Rf為t=∞時應(yīng)力-應(yīng)變方程的參考應(yīng)變和擬合比，故該修正模型共有3個參數(shù)需要確定。
2.3.3 Mesri模型參數(shù)的確定
令t1=1 440 min（1 d），此時的應(yīng)力-應(yīng)變函數(shù)表示t=1 d時的應(yīng)力-應(yīng)變關(guān)系，即
以ε/Dr為縱坐標(biāo)，ε為橫坐標(biāo)，將公式（5）轉(zhuǎn)換為線性方程的形式，對試驗數(shù)據(jù)進(jìn)行線性擬合。直線的截距即為t1=1 440 min（1 d）時刻下的參考應(yīng)變2/（Eu/Su），斜率即為擬合比Rf。擬合結(jié)果如圖6所示。
應(yīng)變-時間雙對數(shù)曲線如圖7所示。圖7中，各個偏應(yīng)力下的曲線呈線性分布，對各個曲線進(jìn)行線性擬合，各條直線斜率的平均值即為λ。
2.3.4 修正的Mesri模型參數(shù)的確定
由公式（4），當(dāng)t=∞時，公式變?yōu)?/p>
將公式（6）代入公式（4）得
設(shè)Y=t/ε，a=1/ε∞，b=T/ε∞，將式（7）變換為Y=ax+b的線性形式，所得t/ε–t關(guān)系如圖8所示。將圖中數(shù)據(jù)進(jìn)行線性擬合，分別得到5級荷載下T的值，并可以求出各級偏應(yīng)力下的ε∞。
以ε∞/Dr為縱坐標(biāo)，ε∞為橫坐標(biāo)，將公式（6）轉(zhuǎn)換為線性方程的形式，對t=∞時刻的應(yīng)力-應(yīng)變關(guān)系進(jìn)行線性擬合。直線的截距即為t=∞時刻下的參考應(yīng)變2/（Eu/Su），斜率即為擬合比Rf。擬合結(jié)果如圖9所示。
2.3.5 分段擬合的蠕變模型
通過以上步驟分別得到Mesri模型和修正Mesri模型的擬合曲線，由于Mesri模型能夠較好地描述蠕變的衰減階段，而修正的Mesri模型能夠較好地描述低蠕變速率的穩(wěn)態(tài)蠕變階段。分別用兩函數(shù)擬合試驗數(shù)據(jù)，用修正的Mesri函數(shù)修正經(jīng)典Mesri函數(shù)的擬合曲線。經(jīng)過分析，以兩條曲線的交點作為分段擬合的分界點，交點對應(yīng)時刻為tj，當(dāng)0＜t＜tj時，黏土的蠕變特性用Mesri模型描述，當(dāng)t≥tj時，黏土的蠕變特性用修正的Mesri模型描述。最終的擬合參數(shù)如表3所示。用分段擬合的Mesri模型對試驗曲線進(jìn)行擬合，所得結(jié)果如圖10～圖12所示。比較圖11與圖3和圖4，分段擬合的Mesri模型能夠更好地描述試驗數(shù)據(jù)的蠕變特性，尤其是對蠕變后期的預(yù)測要優(yōu)于經(jīng)典的Mesri模型和Singh?Mitchell模型。
圖6 ε/Dr?ε關(guān)系曲線Fig.6 ε/Drwith ε under confining pressure of 100 kPa
圖7 lnε-lnt關(guān)系曲線Fig.7 Relationships between lnε and lnt
圖8 100 kPa圍壓下蠕變試驗數(shù)據(jù)換算曲線Fig.8 Conversion of creep curves under 100 kPa confining pressure
圖9 ε∞/Drε∞關(guān)系曲線Fig.9 ε∞/Drwith ε∞ under confining pressure of 100 kPa
3 模型在濱海原狀土試驗數(shù)據(jù)的驗證
天津濱海地區(qū)的淤泥質(zhì)、粉質(zhì)黏土層多呈灰色，有機(jī)質(zhì)含量豐富，具有高含水率、高孔隙比、低強(qiáng)度、高壓縮性、低滲透性等特點。
為了驗證分段擬合的Mesri模型，選取天津港北疆碼頭淤泥質(zhì)黏土及粉質(zhì)黏土的原狀土蠕變試驗數(shù)據(jù)，試驗數(shù)據(jù)來自文獻(xiàn)［14］，試樣是現(xiàn)場鉆孔取得的原狀土，土性參數(shù)如表4。
表3 分段擬合的各個參數(shù)值Tab.3 Values of parameters for piecewise fitting
依據(jù)土層深度，分別選取100 kPa圍壓下淤泥質(zhì)黏土和200 kPa圍壓下粉質(zhì)黏土的蠕變試驗數(shù)據(jù)。加載等級分為4級，對于淤泥質(zhì)黏土，各級偏壓分別為25 kPa、50 kPa、75 kPa和100 kPa，固結(jié)不排水（CU）試驗的破壞偏壓力為134.4 kPa；對于粉質(zhì)黏土，各級偏壓力分別為50 kPa、100 kPa、150 kPa和200 kPa，破壞偏壓力為244.9 kPa。
用分段擬合的Mesri模型去擬合試驗曲線，擬合結(jié)果如圖13、圖14所示。各項參數(shù)見表5。
由圖13、圖14，比較擬合曲線與實測曲線，發(fā)現(xiàn)對于淤泥質(zhì)黏土，分段擬合的Mesri模型能夠準(zhǔn)確地描述其蠕變特性，擬合精度較高；對于粉質(zhì)黏土，除較高應(yīng)力狀態(tài)下擬合結(jié)果存在一定誤差外，其余擬合曲線基本與實測曲線一致，且最大誤差能夠控制在10%以內(nèi)。蠕變函數(shù)能夠很好的描述天津濱海淤泥質(zhì)黏土和粉質(zhì)黏土的蠕變特性。
表4 土樣基本特征表Tab.4 Basic physical parameters of soil samples
圖10 圍壓50 kPa下蠕變擬合曲線Fig.10 Fitting creep curves under 50 kPa confining pressure
表5 分段擬合模型的各個參數(shù)值Tab.5 Values of parameters for modified Mesri model
圖11 圍壓100 kPa下蠕變擬合曲線Fig.11 Fitting creep curves under 100 kPa confining pressure
圖12 圍壓150 kPa下蠕變擬合曲線Fig.12 Fitting creep curves under 150 kPa confining pressure
圖13 圍壓100 kPa下淤泥粘土蠕變擬合曲線Fig.13 Fitting creep curves of muddy clay under 100 kPa confining pressure
圖14 圍壓200 kPa下粉質(zhì)粘土蠕變擬合曲線Fig.14 Fitting creep curves of silt clay under 200 kPa confining pressure
4 結(jié)論
通過室內(nèi)三軸蠕變試驗，得到了天津港地區(qū)重塑飽和軟黏土在不同圍壓狀態(tài)下的應(yīng)變-時間關(guān)系及應(yīng)力-應(yīng)變等時關(guān)系，建立了分段擬合的Mesri經(jīng)驗蠕變模型，通過對試驗結(jié)果和模型的分析得到以下結(jié)論。
（1）在不同圍壓，以及各級偏應(yīng)力作用下，蠕變試驗曲線具有明顯的衰減穩(wěn)定特性，沒有加速蠕變階段；應(yīng)力-應(yīng)變等時曲線具有明顯的非線性特征，除初始時刻外，其余曲線為彼此相似曲線，可以用同一類非線性應(yīng)力-應(yīng)變函數(shù)描述。
（2）在總結(jié)試驗結(jié)果和現(xiàn)有經(jīng)驗蠕變模型的基礎(chǔ)上，建立了修正的Mesri模型，該模型將應(yīng)變-時間關(guān)系修正為雙曲線關(guān)系，用于對Mesri模型的擬合曲線加以修正。在此基礎(chǔ)上，建議了一種分段擬合的Mesri模型，用Mesri模型描述蠕變的衰減階段，用修正的Mesri模型描述蠕變的穩(wěn)態(tài)階段，以兩條擬合曲線的交點作為分段擬合的分界點。該模型更符合文中試驗數(shù)據(jù)，能夠更好地反映和預(yù)測具有衰減穩(wěn)定特征的濱海軟黏土的蠕變特性。
（3）選取天津濱海地區(qū)典型淤泥質(zhì)黏土和粉質(zhì)黏土的原狀土蠕變試驗數(shù)據(jù)對分段擬合的Mesri模型進(jìn)行驗證，擬合結(jié)果較好，由此建立了適用于具有衰減穩(wěn)定蠕變特性的濱海軟黏土非線性蠕變模型。
［1］李越松,趙沖久,趙利平,等.天津港岸坡土體蠕變對高樁碼頭的影響［J］.中國港灣建設(shè),2009（4）：5-8.
LI Y S,ZHAO C J,ZHAO L P,et al.Influence of bank soil creep deformation on high-piled wharf in Tianjin port［J］.China Har?bor Engineering,2009（4）：5–8.
［2］YIN Zhen?Yu,CHANG Ching S,MINNA Karstunen,et al.An anisotropic elastic?viscoplastic model for soft clays［J］.Interna?tional Journal of Solids and Structures,ASCE,2010,47（5）：665-677.
［3］LEONI M,VERMEER P A,KARSTUNEN M.Validation of Anisotropic Creep Model for Soft Soils［C］//KARSTUNEN M,LEO?NI M.Geotechnics of Soft Soils：Focus on Ground Improvement.Glasgow：Taylor&Francis,2008：165-171.
［4］LEONI M,KARSTUNEN M,VERMEER P A.Anisotropic creep model for soft soils［J］.Geotechnique,2008,58（3）：215-226.
［5］VERMEER P A,LEONI M,et al.Modeling and Numerical Simulation of Creep in Soft Soils［C］//CHAN D H,LAW K T.Soft Soil Engineering.Vancouver：Taylor&Francis,2007：57-71.
［6］ KELLN C,SHARMA J,HUGHES D,et al.An Improved Elastic?viscoplastic Soil Model［J］.Canadian Geotechnical Journal,2008,45（10）：1 356-1 376.
［7］ 李軍世,孫鈞.上海淤泥質(zhì)黏土的Mesri蠕變模型［J］.土木工程學(xué)報,2001,34（6）：74-79.
LI J S,SUN J.Mesri′s creep model for Shanghai silt?clay［J］.China Civil Engineering Journal,2001,34（6）：74-79.
［8］王常明,王清,張淑華.濱海軟土蠕變特性及蠕變模型［J］.巖石力學(xué)與工程學(xué)報,2004,23（2）：227-230.
WANG C M,WANG Q,ZHANG S H.Creep characteristics and creep model of marine soft soils［J］.Chinese Journal of Rock Me?chanics and Engineering,2004,23（2）：227－230.
［9］王元戰(zhàn)，王婷婷，王軍.濱海軟土非線性蠕變模型及其工程應(yīng)用研究［J］.巖土力學(xué),2009,30（9）：2 679-2 685.
WANG Y Z,WANG T T,WANG J.A nonlinear rheological model of soft clay and its application to Tianjin littoral area［J］.Rock and Soil Mechanics,2009,30（9）：2 679-2 685.
［10］閆澍旺,劉克瑾,李偉,等.天津濱海新區(qū)軟黏土的蠕變特性及無屈服面模型探究［J］.巖土力學(xué),2010,31（5）：1 431-1 436.
YAN S W,LIU K J,LI W,et al.Study of creep properties of soft clay in Tianjin Binhai New Area and no?yield?surface constitu?tive model［J］.Rock and Soil Mechanics,2010,31（5）：1 431-1 436.
［11］王元戰(zhàn),黃東旭,肖忠,等.天津濱海兩種典型軟黏土蠕變特性試驗研究［J］.巖土工程學(xué)報,2012,34（2）：379-384.
WANG Y Z,HUANG D X,XIAO Z.Experimental research on creep properties of two typical soft clays in coastal region of Tian?jin［J］.Chinese Journal of Geotechnical Engineering,2012,34（2）：379-384.
［12］孫鈞.巖土材料流變及其工程應(yīng)用［M］.北京：中國建筑工業(yè)出版社,1999.
［13］SINGH A,MITCHELL J K.General stress?strain?time function for clay［J］.Journal of the clay mechanics and foundation divi?sion,ASCE,1968，94（SM1）：21-46.
［14］田雙珠，李越松，及春寧.天津港典型軟黏土蠕變試驗研究［J］.水道港口，2009,30（6）：440-443.
TIAN S Z，LI Y S，JI C N.Creep test of soft clay for Tianjin Port［J］.Journal of Waterway and Harbor,2009,30（6）：440-443.
［15］MESRI G，F(xiàn)EBRES?CORDERO E,SHIELDS D R,et al.Shear stress?strain?time behavior of clays［J］.Geotechnique,1981,31（4）：537-552.
［16］王琛,劉浩吾,許強(qiáng).三峽泄灘滑坡滑動帶土的改Mesri蠕變模型［J］.西南交通大學(xué)學(xué)報,2004,39（1）：15-19.
WANG C,LIU H W,XU Q.Modified Mesri′s creep model for soils in sliding zone of Xietan landslide in Three Gorges［J］.Jour?nal of Southwest Jiaotong University,2004,39（1）：15-19.
［17］C C維亞諾夫.土力學(xué)的流變學(xué)原理［M］.北京：科學(xué)出版社，1987：152-156.
Experimental research on creep properties of soft clays in coastal region of Tianjin and its empirical creep model
WANG Yuan?zhan,DONG Yan?he
（School of Civil Engineering,Tianjin Key Laboratory of Port and Ocean Engineering,Tianjin University,Tianjin300072,China）
Based on the indoor triaxial creep tests,the nonlinear creep properties of the remoulded saturated soft clay in coastal area of Tianjin were analyzed.The creep tests were carried out by the step loading method,and the creep curves of the soft clay under different stress states were got by processing the measured data by Mr.Chen′s method.In addition,Singh?Mitchell model and Mesri model were respectively used to fit the test curves which are the stress?strain isochronal curves and strain?time curves.Combining with the fitting results and the creep proper?ties of the coastal clay,a piecewise fitting Mesri model was put forward.In the early stage of the creep,the hyperbol?ic function was taken to describe the stress?strain relationship,and the power function was used to describe the strain?time relationship.In the late stage of the creep,two different kinds of hyperbolic functions were respectively used to matching the stress?strain relationship and the strain?time relationship.The intersection of the fitted curves was taken as the cut?off point.Compared with Singh?Mitchell and Mesri model,this model is more suitable to reflect the attenuation and stable state of the creep characteristics.The test data of two kinds of typical undisturbed clay in Tianjin which are muddy clay and silt clay were taken to verify this model.The fitting results were good.So a non?linear creep model which is suitable for the soft clay in the coastal area of Tianjin is established.
soft clay;creep;attenuation and stable state property;empirical creep model;coastal region of Tianjin
TV 443；O 242.1
A
1005-8443（2014）03-0209-08
2013-09-22；
2013-10-18
國家自然科學(xué)基金項目（51279128）
王元戰(zhàn)（1958-），男，教授，博士生導(dǎo)師，主要從事港口海岸及近海工程結(jié)構(gòu)設(shè)計理論和方法、土與結(jié)構(gòu)相互作用等方面的研究。
猜你喜歡
一類連續(xù)和不連續(xù)分段線性系統(tǒng)的周期解研究數(shù)學(xué)物理學(xué)報(2021年4期)2021-08-30 08:28:02 濱海白首烏華人時刊(2020年21期)2021-01-14 01:33:28 不一般的黏土插畫少年漫畫(藝術(shù)創(chuàng)想)(2020年6期)2020-10-23 09:20:52 濱海頂層公寓現(xiàn)代裝飾(2019年12期)2020-01-19 02:59:16 岳濱海 藏石欣賞寶藏(2018年12期)2019-01-29 01:50:56 分段計算時間小學(xué)生學(xué)習(xí)指導(dǎo)(低年級)(2018年11期)2018-12-03 05:05:00 黏土多肉植物少年漫畫(藝術(shù)創(chuàng)想)(2018年5期)2018-09-11 01:57:32 報紙“黏土”等好孩子畫報(2018年1期)2018-04-14 02:04:16 一起讀吧娃娃樂園·繪本(2016年8期)2016-12-26 21:49:13 3米2分段大力士“大”在哪兒?太空探索(2016年9期)2016-07-12 10:00:04
王元戰(zhàn)，董焱赫
（天津大學(xué) 建筑工程學(xué)院，天津市港口與海洋工程重點實驗室，天津 300072）
天津濱海軟黏土蠕變特性及其經(jīng)驗?zāi)Ｐ脱芯?/p>
王元戰(zhàn)，董焱赫
（天津大學(xué) 建筑工程學(xué)院，天津市港口與海洋工程重點實驗室，天津 300072）
在室內(nèi)三軸蠕變試驗的基礎(chǔ)上，分析了天津濱海地區(qū)重塑飽和軟黏土的非線性蠕變特性。蠕變試驗采用分級加載的方式進(jìn)行，利用“陳氏法”對實測數(shù)據(jù)進(jìn)行處理，得到飽和軟黏土在不同應(yīng)力狀態(tài)下的蠕變試驗曲線。在此基礎(chǔ)上，依據(jù)黏土的應(yīng)力-應(yīng)變、應(yīng)變-時間關(guān)系，分別用Singh?Mitchell和Mesri經(jīng)驗?zāi)Ｐ蛿M合試驗數(shù)據(jù)?？偨Y(jié)擬合結(jié)果，結(jié)合濱海軟黏土的蠕變特性，提出一種分段擬合的Mesri模型，在蠕變前期，模型的應(yīng)力-應(yīng)變關(guān)系采用雙曲線函數(shù)，應(yīng)變-時間關(guān)系采用冪函數(shù)；在蠕變后期，模型的應(yīng)力-應(yīng)變關(guān)系采用雙曲線函數(shù)，應(yīng)變-時間關(guān)系采用雙曲線函數(shù)，以兩條擬合曲線的交點作為分段擬合的分界點。對比分析表明，該方法能更好地描述具有衰減穩(wěn)定特性的蠕變曲線。選取天津濱海地區(qū)原狀土蠕變試驗數(shù)據(jù)對該模型進(jìn)行驗證，擬合結(jié)果較好，由此建立了適用于具有衰減穩(wěn)定蠕變特性的濱海軟黏土非線性蠕變模型。
軟黏土；蠕變；衰減穩(wěn)定；經(jīng)驗?zāi)Ｐ?；天津濱海
Biography：WANG Yuan?zhan（1958-），male，professor.
在中國經(jīng)濟(jì)發(fā)展活躍的長江口、珠江口、天津濱海等海岸地區(qū)，廣泛分布著軟黏土地基，這些地區(qū)經(jīng)濟(jì)發(fā)展迅速，在這些軟土地基上興建了大量的建筑物，以天津港為例，港口規(guī)模不斷擴(kuò)大，吞吐量已經(jīng)躍居世界第四，大量的陸域設(shè)施和圍海造陸工程正在建設(shè)中。但由于該種土類流變效應(yīng)較為明顯，靜載條件下變形緩慢增加，即發(fā)生蠕變效應(yīng)，使得這些建在工程病害層上的建筑物存在著不同程度的沉降，使之有失穩(wěn)、滑塌的危險。如近年天津北疆高樁碼頭多個轉(zhuǎn)角處由于軟黏土岸坡土體的蠕變發(fā)生了基樁連同樁帽與上部結(jié)構(gòu)之間出現(xiàn)較大相對變位的現(xiàn)象，導(dǎo)致橫梁的擱置長度不足，有的只剩下了幾厘米，使得碼頭結(jié)構(gòu)處于危險狀態(tài)之中［1］。開展濱海軟黏土的蠕變試驗研究，建立其合適的蠕變模型，是亟待開展的工作，對研究濱海地區(qū)港口碼頭和防波堤結(jié)構(gòu)的長期穩(wěn)定性有著重要意義。關(guān)于長期荷載下軟黏土的蠕變特性的問題，學(xué)者們展開了大量的研究。
在國外，Zhen［2］等在Perzyna超限應(yīng)力理論的基礎(chǔ)上，建立了一個各向異性的黏彈塑性軟黏土模型，并針對法國St.Herblain地區(qū)的軟土進(jìn)行了等應(yīng)變率蠕變試驗，對所建立的模型進(jìn)行了驗證；Leoni［3-5］等和Ver?meer等考慮土體各向異性，建立軟土蠕變模型，并對軟土的蠕變行為進(jìn)行數(shù)值模擬研究。Kelln［6］等嘗試采用粘彈力學(xué)原理揭示軟粘土的蠕變行為。在國內(nèi)，李軍世［7］等采用Singh?Mitchell應(yīng)變速率-時間關(guān)系方程描述土體的蠕變特性，結(jié)合上海淤泥質(zhì)黏土室內(nèi)試驗，給出了上海淤泥質(zhì)黏土的蠕變形式；王常明［8］等的三軸蠕變試驗研究表明，濱海相沉積的軟土具有非線性蠕變的特性，雙曲線模型更適合其應(yīng)力-應(yīng)變關(guān)系，修改了Singh?Mitchell模型的應(yīng)力-應(yīng)變關(guān)系，建議了濱海軟土應(yīng)力-應(yīng)變-時間關(guān)系模型。王元戰(zhàn)［9］等基于屈服理論，研究了天津濱海淤泥質(zhì)軟黏土的蠕變變形規(guī)律，建立了反映濱海軟黏土的非線性流變本構(gòu)模型；閆澍旺［10］等，采用大直徑原狀試樣進(jìn)行了三軸不固結(jié)不排水蠕變試驗，研究了天津濱海新區(qū)吹填黏土的蠕變變形的規(guī)律性，在此基礎(chǔ)上建立天津濱海新區(qū)軟黏土的無屈服面蠕變模型；王元戰(zhàn)［11］等采用Merchant蠕變模型三軸蠕變試驗的結(jié)果曲線進(jìn)行擬合，建立了適用于天津濱海地區(qū)的軟黏土本構(gòu)模型。
雖然目前有很多關(guān)于土體蠕變模型的研究成果，也不乏對天津濱海軟黏土蠕變特性的研究，但是這些成果不足以完全詮釋天津濱海地區(qū)軟黏土的蠕變特性，原因有2個：一是蠕變模型自身復(fù)雜性和局限性。二是以往研究得出的模型大多只能描述自身的試驗數(shù)據(jù)，沒有用其他試驗數(shù)據(jù)加以驗證，缺乏綜合比較與分析。
本文以天津港病害層重塑淤泥質(zhì)黏土為例進(jìn)行蠕變試驗研究。依據(jù)試驗數(shù)據(jù)，結(jié)合Singh?Mitchell和Mesri經(jīng)驗蠕變模型，提出一種分段擬合的Mesri模型，在蠕變前期，模型的應(yīng)力-應(yīng)變關(guān)系采用雙曲線函數(shù)，應(yīng)變-時間關(guān)系采用冪函數(shù)；在蠕變后期，模型的應(yīng)力-應(yīng)變關(guān)系采用雙曲線函數(shù)，應(yīng)變-時間關(guān)系采用雙曲線函數(shù)，以兩條擬合曲線的交點作為分段擬合的分界點。對比分析表明，該模型能更好地描述具有衰減穩(wěn)定特性的蠕變曲線。選取以往天津濱海地區(qū)兩種典型軟黏土蠕變試驗數(shù)據(jù)對該模型進(jìn)行驗證，擬合結(jié)果較好，由此建立了適用于具有衰減穩(wěn)定蠕變特性的濱海軟黏土非線性蠕變模型。
1 三軸蠕變試驗
1.1 試驗設(shè)備及試樣
設(shè)備采用的是無級加載應(yīng)力控制式三軸蠕變試驗儀，該三軸蠕變試驗儀由普通應(yīng)變式三軸儀通過改裝軸向加載設(shè)備而成，儀器主要由軸向加載設(shè)備、圍壓控制測量系統(tǒng)、反壓控制測量系統(tǒng)、體積變化和孔隙水壓力測量系統(tǒng)等組成，將軸向加載設(shè)備由原來的應(yīng)變控制式改造成應(yīng)力控制式加載。這樣改造試驗儀器簡單易行，且既可以利用原有應(yīng)變式三軸儀器設(shè)備優(yōu)良穩(wěn)定的圍壓控制測量系統(tǒng)、反壓控制測量系統(tǒng)、體積變化和孔隙水壓力測量系統(tǒng)，又可以通過受力平衡計算，利用干砂和砝碼施加各級軸向荷載，獲得恒定且準(zhǔn)確的軸向偏應(yīng)力。
黏土試樣為天津港地區(qū)5～10 m土層飽和重塑軟黏土，土樣呈流塑狀態(tài)。依據(jù)土工試驗規(guī)程（SL237-1999）制作重塑試樣。試樣直徑39.1 mm，高為80 mm。土樣的基本指標(biāo)如表1所示。
1.2 試驗方法及加載方案
在土的室內(nèi)流變試驗中，對于蠕變試驗有兩種不同的加載方式，即分別加載法和分級加載法。所謂分級加載，就是在同一試樣上逐級加上不同的應(yīng)力，即在一級應(yīng)力水平下蠕變經(jīng)歷給定的時間或達(dá)到穩(wěn)定后，將應(yīng)力水平提高到下一級的水平，直到所需的應(yīng)力水平［12］。為了克服分別加載的局限，減小采用不同土樣進(jìn)行試驗對試驗結(jié)果帶來的離散性，本試驗加載方式采用分級加載法。
由于軟黏土具有滲透性差的特點，試驗采用不排水蠕變試驗，圍壓分別取50 kPa，100 kPa，150 kPa。在蠕變試驗進(jìn)行前，應(yīng)選取同一批次重塑土樣進(jìn)行圍壓50 kPa、100 kPa和150 kPa下的固結(jié)不排水剪切試驗，以試樣應(yīng)變達(dá)到15%為破壞標(biāo)準(zhǔn)，確定破壞偏應(yīng)力qf分別為70.2 kPa，120.5 kPa和178.8 kPa。根據(jù)qf確定每級加載偏應(yīng)力，取加載級數(shù)為5級，加載方案如表2所示，再通過受力平衡計算各級加砂質(zhì)量。
表1 土樣基本特征表Tab.1 Basic physical parameters of soil samples
蠕變試驗步驟為：第1步，將制好的試樣放入真空室中，將室內(nèi)抽氣呈負(fù)壓狀態(tài)，10 h以后打開排水閥，使土樣進(jìn)水飽和。第2步，飽和完成后，將試樣安放在三軸壓力室進(jìn)行固結(jié)，由于軟黏土排水性差，故固結(jié)時間設(shè)為24 h。第3步，固結(jié)基本完成以后，設(shè)置量測豎向變形的百分表初始讀數(shù)，關(guān)閉排水閥，再施加第一級偏應(yīng)力ΔF1，同時通過百分表量測軸向變形。在偏應(yīng)力作用下蠕變穩(wěn)定或經(jīng)歷給定的時間（本次試驗每級偏應(yīng)力施加時間確定為48 h）后，再增加配重，使偏壓力達(dá)到ΔF2的水平，以此直到五級荷載施加完畢。
表2 土樣分級加載方案Tab.2 Schemes of step loading kPa
1.3 試驗結(jié)果分析
限于篇幅，以下以圍壓100 kPa下的蠕變試驗為例，進(jìn)行分析。試驗得到的曲線為分級加載下的應(yīng)變-時間曲線，該流變曲線的可用性較小，需要分級進(jìn)行整理，通常的整理方法有線性法和陳氏法。線性法是將土體視為線性流變體，此時土體蠕變響應(yīng)滿足線性疊加原理，蠕變曲線可以簡單地由每一級荷載的蠕變響應(yīng)疊加而成。
若要更準(zhǔn)確的研究土的蠕變，應(yīng)將土體視為非線性流變體。由陳宗基教授提出并由他的學(xué)生所發(fā)展的“陳氏法”通過采用適當(dāng)?shù)膶嶒灱夹g(shù)與方法，用作圖法建立真實變形過程的疊加關(guān)系，在考慮土體為非線性流變體或線性流變體時均可適用，在國內(nèi)外巖土流變學(xué)研究中得到越來越廣泛的應(yīng)用。陳氏法可以使蠕變試驗通過采用適當(dāng)?shù)募虞d程序從單一試樣得到更多的試驗資料，這對試驗周期長，并且試驗價格昂貴的巖土蠕變性質(zhì)研究有重大的實際意義。本試驗對蠕變試驗結(jié)果采用的整理方法就是陳氏法。
用“陳氏法”處理得到不同應(yīng)力狀態(tài)下的軸向應(yīng)變與加載時間的蠕變關(guān)系曲線如圖1所示，由圖1可知：（1）在不排水條件下，當(dāng)主應(yīng)力差達(dá)到破壞偏應(yīng)力前，變形呈現(xiàn)出明顯的衰減穩(wěn)態(tài)蠕變過程。在各級應(yīng)力下，應(yīng)變隨著時間的增長逐漸趨于穩(wěn)定。當(dāng)主應(yīng)力差較低時，應(yīng)變呈現(xiàn)出明顯的衰減蠕變趨勢，蠕變現(xiàn)象不明顯；當(dāng)應(yīng)力水平較高時，試樣蠕變曲線包括衰減階段和穩(wěn)態(tài)蠕變階段，試樣穩(wěn)態(tài)蠕變階段的斜率隨主應(yīng)力差的增加而增加；（2）在不排水條件下，試樣的應(yīng)變隨著偏應(yīng)力水平的提高而越來越大，且增加幅度越來越大。在試樣破壞之前，變形曲線并未出現(xiàn)加速蠕變階段。由蠕變曲線得到的應(yīng)力-應(yīng)變等時曲線如圖2所示：不同時刻，除初始時刻外，各個時刻的應(yīng)力-應(yīng)變曲線為彼此相似的曲線，可以用同一類函數(shù)描述；在同一時刻下，應(yīng)力-應(yīng)變曲線具有明顯屈服特性，曲線呈非線性，這說明，需要用非線性的蠕變模型來描述濱海軟黏土蠕變特性。
圖1 應(yīng)變-時間關(guān)系圖Fig.1 ε?t curves of creep test
圖2 應(yīng)力-應(yīng)變等時曲線Fig.2 σ?ε isochronal curves
2 蠕變經(jīng)驗?zāi)Ｐ?/h2>
蠕變的經(jīng)驗?zāi)Ｐ停菑牧鲃永碚摶蚶匣碚摮霭l(fā)，直接給出土蠕變方程的應(yīng)力-應(yīng)變-時間函數(shù)形式。對每種不同的材料，甚至不同的條件，有多樣的經(jīng)驗?zāi)Ｊ?，和元件模型、屈服面模型、?nèi)時模型等相比，它比較直觀，參數(shù)簡單且容易獲取，亦為工程設(shè)計人員所樂意采用。目前常用的具有代表性的包括Singh?Mitchell蠕變速率和蠕變函數(shù)，常速率軸向應(yīng)變條件下的雙曲線計算模型以及Mesri應(yīng)力-應(yīng)變-時間蠕變函數(shù)等。
2.1 Singh?Mitchell模型
1968年Singh等提出的著名的Singh?Mitchell經(jīng)驗蠕變模型［13-14］是在總結(jié)了單級常應(yīng)力加載，排水與不排水三軸壓縮試驗數(shù)據(jù)的基礎(chǔ)上，采用指數(shù)函數(shù)表示應(yīng)力-應(yīng)變關(guān)系、采用冪函數(shù)表示應(yīng)變-時間關(guān)系的經(jīng)典的經(jīng)驗蠕變模型，該模型能夠很好地反映土的蠕變特性。Singh?Mitchell蠕變模型寫成應(yīng)變率-應(yīng)力水平-時間的關(guān)系式，可表示為
式中：ε˙為任一時刻t的軸向應(yīng)變速率；t為試樣受荷時間；Dr=（σ1-σ3）/（σ1-σ3）f為偏應(yīng)力水平，（σ1-σ3）f為破壞偏應(yīng)力，可由常規(guī)三軸排水壓縮實驗得到；Ar是在單位參考時刻tr，偏應(yīng)力水平為零時的應(yīng)變速率；m由ln˙-lnt擬合直線的斜率確定；α為應(yīng)變速率對數(shù)與剪應(yīng)力關(guān)系圖中線性段的斜率。
對式（1）進(jìn)行積分，在不考慮初始應(yīng)變且m不為1時。經(jīng)推導(dǎo)可得
式中：B=Artr/（1-m）；β=α；λ=1-m。該模型需要確定的3個參數(shù)是B，β，λ。
2.2 Mesri模型
針對Singh?Mitchell模型存在的局限性，Mesri.G［15-16］于1981年提出了新的經(jīng)驗?zāi)Ｐ?。該模型?yīng)力-應(yīng)變關(guān)系采用雙曲線函數(shù)，應(yīng)變-時間關(guān)系仍采用冪函數(shù)。經(jīng)典的Mesri模型公式如下
式中有3個參數(shù)：參考應(yīng)變2/（Eu/Su），擬合比Rf和蠕變參數(shù)λ。
用Singh?Mitchell模型和Mesri模型分別擬合試驗曲線，擬合結(jié)果如圖3和圖4所示。
2.3 分段擬合的Mesri模型
從擬合結(jié)果可以看出，Mesri模型和Singh?Mitchell模型都存在后期擬合不準(zhǔn)確的情況，原因在于兩個模型的應(yīng)變-時間關(guān)系都采用了冪函數(shù)的形式，而冪函數(shù)是非衰減函數(shù)，當(dāng)時間增大時，曲線斜率較大，對應(yīng)的應(yīng)變值遠(yuǎn)大于實測值，對于尚處于衰減穩(wěn)定期的試驗數(shù)據(jù)是不適用的。
2.3.1變形函數(shù)（應(yīng)力-應(yīng)變關(guān)系）的確定
當(dāng)取某一時刻t1時，公式（2）和（3）分別轉(zhuǎn)化為t1時刻的應(yīng)力-應(yīng)變方程。以100 kPa圍壓下的實測數(shù)據(jù)為例，取t1=1 d，分別用兩種函數(shù)去擬合試驗數(shù)據(jù)，通過應(yīng)變對數(shù)-偏應(yīng)力水平曲線（圖5）可以看出，當(dāng)偏應(yīng)力水平處于中間段時，實測數(shù)據(jù)呈線性關(guān)系，而在偏應(yīng)力水平趨近于1或趨近于0時，實測數(shù)據(jù)則呈現(xiàn)出明顯的非線性；指數(shù)函數(shù)關(guān)系只能描述偏應(yīng)力水平處于中間段（20%～80%偏應(yīng)力水平）的應(yīng)力-應(yīng)變關(guān)系，而雙曲函數(shù)能夠準(zhǔn)確地描述全部應(yīng)力水平下的狀態(tài)。故應(yīng)力-應(yīng)變關(guān)系應(yīng)采用雙曲線函數(shù)描述。
圖3 試驗與Sing?Mitchell模型計算結(jié)果比較Fig.3 Comparison between test data and calculated results using S?M model
圖4 試驗數(shù)據(jù)與Mesri模型計算結(jié)果比較Fig.4 Comparison between test data and calculated results using Mesri model
圖5 應(yīng)力-應(yīng)變等時曲線Fig.5 σ?ε isochronal curves
2.3.2 蠕變函數(shù)的確定
結(jié)合上文研究，為了修正Mesri模型的誤差，由C.C.維亞洛夫提出的通用的時間函數(shù)［17］入手，本文建議了一種分段擬合的經(jīng)驗?zāi)Ｐ?，該模型前期仍采用Mesri函數(shù)的形式；在擬合后期數(shù)據(jù)時將Mesri函數(shù)應(yīng)變-時間關(guān)系修正為雙曲線函數(shù)，該函數(shù)的形式為
式中：T為蠕變參數(shù)，2/（Eu/Su）和Rf為t=∞時應(yīng)力-應(yīng)變方程的參考應(yīng)變和擬合比，故該修正模型共有3個參數(shù)需要確定。
2.3.3 Mesri模型參數(shù)的確定
令t1=1 440 min（1 d），此時的應(yīng)力-應(yīng)變函數(shù)表示t=1 d時的應(yīng)力-應(yīng)變關(guān)系，即
以ε/Dr為縱坐標(biāo)，ε為橫坐標(biāo)，將公式（5）轉(zhuǎn)換為線性方程的形式，對試驗數(shù)據(jù)進(jìn)行線性擬合。直線的截距即為t1=1 440 min（1 d）時刻下的參考應(yīng)變2/（Eu/Su），斜率即為擬合比Rf。擬合結(jié)果如圖6所示。
應(yīng)變-時間雙對數(shù)曲線如圖7所示。圖7中，各個偏應(yīng)力下的曲線呈線性分布，對各個曲線進(jìn)行線性擬合，各條直線斜率的平均值即為λ。
2.3.4 修正的Mesri模型參數(shù)的確定
由公式（4），當(dāng)t=∞時，公式變?yōu)?/p>
將公式（6）代入公式（4）得
設(shè)Y=t/ε，a=1/ε∞，b=T/ε∞，將式（7）變換為Y=ax+b的線性形式，所得t/ε–t關(guān)系如圖8所示。將圖中數(shù)據(jù)進(jìn)行線性擬合，分別得到5級荷載下T的值，并可以求出各級偏應(yīng)力下的ε∞。
以ε∞/Dr為縱坐標(biāo)，ε∞為橫坐標(biāo)，將公式（6）轉(zhuǎn)換為線性方程的形式，對t=∞時刻的應(yīng)力-應(yīng)變關(guān)系進(jìn)行線性擬合。直線的截距即為t=∞時刻下的參考應(yīng)變2/（Eu/Su），斜率即為擬合比Rf。擬合結(jié)果如圖9所示。
2.3.5 分段擬合的蠕變模型
通過以上步驟分別得到Mesri模型和修正Mesri模型的擬合曲線，由于Mesri模型能夠較好地描述蠕變的衰減階段，而修正的Mesri模型能夠較好地描述低蠕變速率的穩(wěn)態(tài)蠕變階段。分別用兩函數(shù)擬合試驗數(shù)據(jù)，用修正的Mesri函數(shù)修正經(jīng)典Mesri函數(shù)的擬合曲線。經(jīng)過分析，以兩條曲線的交點作為分段擬合的分界點，交點對應(yīng)時刻為tj，當(dāng)0＜t＜tj時，黏土的蠕變特性用Mesri模型描述，當(dāng)t≥tj時，黏土的蠕變特性用修正的Mesri模型描述。最終的擬合參數(shù)如表3所示。用分段擬合的Mesri模型對試驗曲線進(jìn)行擬合，所得結(jié)果如圖10～圖12所示。比較圖11與圖3和圖4，分段擬合的Mesri模型能夠更好地描述試驗數(shù)據(jù)的蠕變特性，尤其是對蠕變后期的預(yù)測要優(yōu)于經(jīng)典的Mesri模型和Singh?Mitchell模型。
圖6 ε/Dr?ε關(guān)系曲線Fig.6 ε/Drwith ε under confining pressure of 100 kPa
圖7 lnε-lnt關(guān)系曲線Fig.7 Relationships between lnε and lnt
圖8 100 kPa圍壓下蠕變試驗數(shù)據(jù)換算曲線Fig.8 Conversion of creep curves under 100 kPa confining pressure
圖9 ε∞/Drε∞關(guān)系曲線Fig.9 ε∞/Drwith ε∞ under confining pressure of 100 kPa
3 模型在濱海原狀土試驗數(shù)據(jù)的驗證
天津濱海地區(qū)的淤泥質(zhì)、粉質(zhì)黏土層多呈灰色，有機(jī)質(zhì)含量豐富，具有高含水率、高孔隙比、低強(qiáng)度、高壓縮性、低滲透性等特點。
為了驗證分段擬合的Mesri模型，選取天津港北疆碼頭淤泥質(zhì)黏土及粉質(zhì)黏土的原狀土蠕變試驗數(shù)據(jù)，試驗數(shù)據(jù)來自文獻(xiàn)［14］，試樣是現(xiàn)場鉆孔取得的原狀土，土性參數(shù)如表4。
表3 分段擬合的各個參數(shù)值Tab.3 Values of parameters for piecewise fitting
依據(jù)土層深度，分別選取100 kPa圍壓下淤泥質(zhì)黏土和200 kPa圍壓下粉質(zhì)黏土的蠕變試驗數(shù)據(jù)。加載等級分為4級，對于淤泥質(zhì)黏土，各級偏壓分別為25 kPa、50 kPa、75 kPa和100 kPa，固結(jié)不排水（CU）試驗的破壞偏壓力為134.4 kPa；對于粉質(zhì)黏土，各級偏壓力分別為50 kPa、100 kPa、150 kPa和200 kPa，破壞偏壓力為244.9 kPa。
用分段擬合的Mesri模型去擬合試驗曲線，擬合結(jié)果如圖13、圖14所示。各項參數(shù)見表5。
由圖13、圖14，比較擬合曲線與實測曲線，發(fā)現(xiàn)對于淤泥質(zhì)黏土，分段擬合的Mesri模型能夠準(zhǔn)確地描述其蠕變特性，擬合精度較高；對于粉質(zhì)黏土，除較高應(yīng)力狀態(tài)下擬合結(jié)果存在一定誤差外，其余擬合曲線基本與實測曲線一致，且最大誤差能夠控制在10%以內(nèi)。蠕變函數(shù)能夠很好的描述天津濱海淤泥質(zhì)黏土和粉質(zhì)黏土的蠕變特性。
表4 土樣基本特征表Tab.4 Basic physical parameters of soil samples
圖10 圍壓50 kPa下蠕變擬合曲線Fig.10 Fitting creep curves under 50 kPa confining pressure
表5 分段擬合模型的各個參數(shù)值Tab.5 Values of parameters for modified Mesri model
圖11 圍壓100 kPa下蠕變擬合曲線Fig.11 Fitting creep curves under 100 kPa confining pressure
圖12 圍壓150 kPa下蠕變擬合曲線Fig.12 Fitting creep curves under 150 kPa confining pressure
圖13 圍壓100 kPa下淤泥粘土蠕變擬合曲線Fig.13 Fitting creep curves of muddy clay under 100 kPa confining pressure
圖14 圍壓200 kPa下粉質(zhì)粘土蠕變擬合曲線Fig.14 Fitting creep curves of silt clay under 200 kPa confining pressure
4 結(jié)論
通過室內(nèi)三軸蠕變試驗，得到了天津港地區(qū)重塑飽和軟黏土在不同圍壓狀態(tài)下的應(yīng)變-時間關(guān)系及應(yīng)力-應(yīng)變等時關(guān)系，建立了分段擬合的Mesri經(jīng)驗蠕變模型，通過對試驗結(jié)果和模型的分析得到以下結(jié)論。
（1）在不同圍壓，以及各級偏應(yīng)力作用下，蠕變試驗曲線具有明顯的衰減穩(wěn)定特性，沒有加速蠕變階段；應(yīng)力-應(yīng)變等時曲線具有明顯的非線性特征，除初始時刻外，其余曲線為彼此相似曲線，可以用同一類非線性應(yīng)力-應(yīng)變函數(shù)描述。
（2）在總結(jié)試驗結(jié)果和現(xiàn)有經(jīng)驗蠕變模型的基礎(chǔ)上，建立了修正的Mesri模型，該模型將應(yīng)變-時間關(guān)系修正為雙曲線關(guān)系，用于對Mesri模型的擬合曲線加以修正。在此基礎(chǔ)上，建議了一種分段擬合的Mesri模型，用Mesri模型描述蠕變的衰減階段，用修正的Mesri模型描述蠕變的穩(wěn)態(tài)階段，以兩條擬合曲線的交點作為分段擬合的分界點。該模型更符合文中試驗數(shù)據(jù)，能夠更好地反映和預(yù)測具有衰減穩(wěn)定特征的濱海軟黏土的蠕變特性。
（3）選取天津濱海地區(qū)典型淤泥質(zhì)黏土和粉質(zhì)黏土的原狀土蠕變試驗數(shù)據(jù)對分段擬合的Mesri模型進(jìn)行驗證，擬合結(jié)果較好，由此建立了適用于具有衰減穩(wěn)定蠕變特性的濱海軟黏土非線性蠕變模型。
［1］李越松,趙沖久,趙利平,等.天津港岸坡土體蠕變對高樁碼頭的影響［J］.中國港灣建設(shè),2009（4）：5-8.
LI Y S,ZHAO C J,ZHAO L P,et al.Influence of bank soil creep deformation on high-piled wharf in Tianjin port［J］.China Har?bor Engineering,2009（4）：5–8.
［2］YIN Zhen?Yu,CHANG Ching S,MINNA Karstunen,et al.An anisotropic elastic?viscoplastic model for soft clays［J］.Interna?tional Journal of Solids and Structures,ASCE,2010,47（5）：665-677.
［3］LEONI M,VERMEER P A,KARSTUNEN M.Validation of Anisotropic Creep Model for Soft Soils［C］//KARSTUNEN M,LEO?NI M.Geotechnics of Soft Soils：Focus on Ground Improvement.Glasgow：Taylor&Francis,2008：165-171.
［4］LEONI M,KARSTUNEN M,VERMEER P A.Anisotropic creep model for soft soils［J］.Geotechnique,2008,58（3）：215-226.
［5］VERMEER P A,LEONI M,et al.Modeling and Numerical Simulation of Creep in Soft Soils［C］//CHAN D H,LAW K T.Soft Soil Engineering.Vancouver：Taylor&Francis,2007：57-71.
［6］ KELLN C,SHARMA J,HUGHES D,et al.An Improved Elastic?viscoplastic Soil Model［J］.Canadian Geotechnical Journal,2008,45（10）：1 356-1 376.
［7］ 李軍世,孫鈞.上海淤泥質(zhì)黏土的Mesri蠕變模型［J］.土木工程學(xué)報,2001,34（6）：74-79.
LI J S,SUN J.Mesri′s creep model for Shanghai silt?clay［J］.China Civil Engineering Journal,2001,34（6）：74-79.
［8］王常明,王清,張淑華.濱海軟土蠕變特性及蠕變模型［J］.巖石力學(xué)與工程學(xué)報,2004,23（2）：227-230.
WANG C M,WANG Q,ZHANG S H.Creep characteristics and creep model of marine soft soils［J］.Chinese Journal of Rock Me?chanics and Engineering,2004,23（2）：227－230.
［9］王元戰(zhàn)，王婷婷，王軍.濱海軟土非線性蠕變模型及其工程應(yīng)用研究［J］.巖土力學(xué),2009,30（9）：2 679-2 685.
WANG Y Z,WANG T T,WANG J.A nonlinear rheological model of soft clay and its application to Tianjin littoral area［J］.Rock and Soil Mechanics,2009,30（9）：2 679-2 685.
［10］閆澍旺,劉克瑾,李偉,等.天津濱海新區(qū)軟黏土的蠕變特性及無屈服面模型探究［J］.巖土力學(xué),2010,31（5）：1 431-1 436.
YAN S W,LIU K J,LI W,et al.Study of creep properties of soft clay in Tianjin Binhai New Area and no?yield?surface constitu?tive model［J］.Rock and Soil Mechanics,2010,31（5）：1 431-1 436.
［11］王元戰(zhàn),黃東旭,肖忠,等.天津濱海兩種典型軟黏土蠕變特性試驗研究［J］.巖土工程學(xué)報,2012,34（2）：379-384.
WANG Y Z,HUANG D X,XIAO Z.Experimental research on creep properties of two typical soft clays in coastal region of Tian?jin［J］.Chinese Journal of Geotechnical Engineering,2012,34（2）：379-384.
［12］孫鈞.巖土材料流變及其工程應(yīng)用［M］.北京：中國建筑工業(yè)出版社,1999.
［13］SINGH A,MITCHELL J K.General stress?strain?time function for clay［J］.Journal of the clay mechanics and foundation divi?sion,ASCE,1968，94（SM1）：21-46.
［14］田雙珠，李越松，及春寧.天津港典型軟黏土蠕變試驗研究［J］.水道港口，2009,30（6）：440-443.
TIAN S Z，LI Y S，JI C N.Creep test of soft clay for Tianjin Port［J］.Journal of Waterway and Harbor,2009,30（6）：440-443.
［15］MESRI G，F(xiàn)EBRES?CORDERO E,SHIELDS D R,et al.Shear stress?strain?time behavior of clays［J］.Geotechnique,1981,31（4）：537-552.
［16］王琛,劉浩吾,許強(qiáng).三峽泄灘滑坡滑動帶土的改Mesri蠕變模型［J］.西南交通大學(xué)學(xué)報,2004,39（1）：15-19.
WANG C,LIU H W,XU Q.Modified Mesri′s creep model for soils in sliding zone of Xietan landslide in Three Gorges［J］.Jour?nal of Southwest Jiaotong University,2004,39（1）：15-19.
［17］C C維亞諾夫.土力學(xué)的流變學(xué)原理［M］.北京：科學(xué)出版社，1987：152-156.
Experimental research on creep properties of soft clays in coastal region of Tianjin and its empirical creep model
WANG Yuan?zhan,DONG Yan?he
（School of Civil Engineering,Tianjin Key Laboratory of Port and Ocean Engineering,Tianjin University,Tianjin300072,China）
Based on the indoor triaxial creep tests,the nonlinear creep properties of the remoulded saturated soft clay in coastal area of Tianjin were analyzed.The creep tests were carried out by the step loading method,and the creep curves of the soft clay under different stress states were got by processing the measured data by Mr.Chen′s method.In addition,Singh?Mitchell model and Mesri model were respectively used to fit the test curves which are the stress?strain isochronal curves and strain?time curves.Combining with the fitting results and the creep proper?ties of the coastal clay,a piecewise fitting Mesri model was put forward.In the early stage of the creep,the hyperbol?ic function was taken to describe the stress?strain relationship,and the power function was used to describe the strain?time relationship.In the late stage of the creep,two different kinds of hyperbolic functions were respectively used to matching the stress?strain relationship and the strain?time relationship.The intersection of the fitted curves was taken as the cut?off point.Compared with Singh?Mitchell and Mesri model,this model is more suitable to reflect the attenuation and stable state of the creep characteristics.The test data of two kinds of typical undisturbed clay in Tianjin which are muddy clay and silt clay were taken to verify this model.The fitting results were good.So a non?linear creep model which is suitable for the soft clay in the coastal area of Tianjin is established.
soft clay;creep;attenuation and stable state property;empirical creep model;coastal region of Tianjin
TV 443；O 242.1
A
1005-8443（2014）03-0209-08
2013-09-22；
2013-10-18
國家自然科學(xué)基金項目（51279128）
王元戰(zhàn)（1958-），男，教授，博士生導(dǎo)師，主要從事港口海岸及近海工程結(jié)構(gòu)設(shè)計理論和方法、土與結(jié)構(gòu)相互作用等方面的研究。
猜你喜歡
數(shù)學(xué)物理學(xué)報(2021年4期)2021-08-30 08:28:02
華人時刊(2020年21期)2021-01-14 01:33:28
少年漫畫(藝術(shù)創(chuàng)想)(2020年6期)2020-10-23 09:20:52
現(xiàn)代裝飾(2019年12期)2020-01-19 02:59:16
寶藏(2018年12期)2019-01-29 01:50:56
小學(xué)生學(xué)習(xí)指導(dǎo)(低年級)(2018年11期)2018-12-03 05:05:00
少年漫畫(藝術(shù)創(chuàng)想)(2018年5期)2018-09-11 01:57:32
好孩子畫報(2018年1期)2018-04-14 02:04:16
娃娃樂園·繪本(2016年8期)2016-12-26 21:49:13
太空探索(2016年9期)2016-07-12 10:00:04
