秦愛芳,呂康立

傳統(tǒng)的土力學研究基本上都是以飽和土為主,而在實際工程中,大多數(shù)的土處于非飽和狀態(tài),在干旱、半干旱地區(qū)尤為如此.對于工程中的非飽和土路堤和土壩等問題,若仍然用飽和土的理論來指導施工的話,將會影響工程建造的質(zhì)量以及施工的進度.

對于非飽和土,典型的固結(jié)理論分別由Blight[1],Scott[2],Barden[3]和Fredlund等[4-5]提出,其中Fredlund固結(jié)理論比較完善且具有權(quán)威性.

對于飽和黏彈性地基固結(jié)問題,王奎華等[6]、蔡袁強等[7-8]、謝康和等[9]分別采用Merchant和Kelvin模型進行了大量的深入研究.對于非飽和黏彈性土地基在變荷載作用下的固結(jié)問題,國內(nèi)外相關(guān)研究還較少.Shan等[10]對恒載和加荷隨時間指數(shù)變化時多層非飽和土彈性地基的固結(jié)進行了研究.Zhou等[11]研究了多種初始條件和邊界條件下,變荷載作用時的非飽和土彈性地基一維固結(jié)問題.Liem等[12]將非飽和土彈性地基的固結(jié)理論拓展到了二維空間.

對于非飽和黏彈性土地基一維固結(jié)問題,秦愛芳等[13-15]依據(jù)Fredlund非飽和土一維固結(jié)理論,并對其進行了一些假設,得到了瞬時加荷下黏彈性土地基一維固結(jié)的半解析解.本工作在上述研究的基礎上,對實際工程中常用的3種循環(huán)荷載(正弦荷載、三角形荷載和梯形荷載)作用下的黏彈性地基,采用Merchant黏彈性模型求得了一維固結(jié)的半解析解,并對其固結(jié)特性進行了分析.研究結(jié)果對工程中的實際問題具有指導作用.

1 半解析解的推導

1.1 計算模型和本構(gòu)方程

本工作中的基本假定與Fredlund的非飽和土一維固結(jié)理論[5]相同,計算模型如圖1所示.非飽和土層的頂面為透水透氣面,底面為不透水不透氣面.土層厚度為H,施加的循環(huán)荷載為σ(t),坐標原點設在地表,深度方向(z軸)以豎直向下為正.取底面積為1、高度為dz的單元體V0=1×1×dz為研究對象.

黏彈性模型采用Merchant模型(三單元模型),該模型由一個Kelvin體(由一個彈性體和一個黏性體并聯(lián)而成)和一個彈性體串聯(lián)組成,如圖2所示.Merchant模型的本構(gòu)方程為

式中,σ為應力,η為黏性體的黏滯系數(shù),ε為應變,E0和E1為彈性體的彈性模量.對式(1)作關(guān)于時間t的Laplace變換,得到

式中,s為Laplace變量.

圖1 計算模型Fig.1 Calculation model

圖 2 Merchant模型Fig.2 Merchant model

1.2 任意加荷時Laplace變換下的解

本工作依據(jù)Fredlund一維固結(jié)理論的本構(gòu)方程[5],針對Merchant黏彈性地基模型,采用李氏比擬法,引入Laplace變換下的柔度系數(shù)V(s)代替線彈性模型中的常數(shù)1/E,再通過線彈性情況下的液相和氣相的控制方程、Darcy定律及Fick定律,采用文獻[13]中的方法,得到Laplace變換下的4個微分方程.

式中:

其中

將式(3),(4),(5),(6)寫為矩陣形式的偏微分方程[13],即

式中,~X(z,s)=[~ua(z,s),~uw(z,s),~Ja(z,s),~vw(z,s)]T,

根據(jù)Cayley-Hamilton理論[16],式(7)的解的一般形式為

式中,T(z,s)=exp[zA(s)],S(z,s)=0T(z ? ?)B(s)d?.由式(8)可得,土層頂面的狀態(tài)向量與任意深度處狀態(tài)向量的傳遞關(guān)系為

式中,T11～T44及S1～S4的含義同文獻[13].

1.3 邊界條件和初始條件

因為土層上表面為透氣透水面,底面為不透氣不透水面,所以邊界條件為

初始條件為

將邊界條件和初始條件代入式(9),可解得Laplace變換下任意深度處的超孔隙氣壓力~ua(z,s)、超孔隙水壓力~uw(z,s)和土層壓縮量~w(z,s),方法同文獻[13].

2 循環(huán)荷載及其Laplace變換下的解

2.1 正弦荷載

式中:σ0為外荷載的峰值,取σ0=100 kPa;w為圓頻率,取w=10?6rad/s;t為荷載加載時間.

圖3 正弦荷載Fig.3 Sinusoidal load

對式(15)作關(guān)于t的Laplace變換,可得

2.2 三角形荷載

三角形荷載如圖4所示,表達式為

圖4 三角形荷載Fig.4 Triangle load

對式(17)作關(guān)于時間t的Laplace變換,可得

式中,σ0=100 kPa,T=106s.

2.3 梯形荷載

梯形荷載如圖5所示,表達式為

對式(19)作關(guān)于時間t的Laplace變換,可得

式中,σ0=100 kPa,T=106s.

將上述3種循環(huán)荷載情況下的~σ(s),σ(0)分別代入式(12),(13)和(14),可得Laplace變換下的解,然后再進行Laplace逆變換,即可求得3種循環(huán)荷載作用下任意深度的超孔隙氣壓力、超孔隙水壓力及土層沉降量.

由于本工作中的~ua(z,s),~uw(z,s)和~w(z,s)非常復雜,通過Laplace逆變換無法直接獲得超孔隙氣壓力、超孔隙水壓力及土層沉降的數(shù)學表達式,因而需要通過數(shù)值法求解逆變換.本工作運用較常用且精度較高的Durbin和Crump方法[17],通過編制程序?qū)ua(z,s),~uw(z,s)和 ~w(z,s)進行Laplace逆變換.

3 算例分析

下面分析3種循環(huán)荷載作用下非飽和黏彈性地基固結(jié)過程中的超孔隙水壓力、超孔隙氣壓力以及固結(jié)度隨時間的變化規(guī)律.固結(jié)度為各時刻的土體沉降與土體最終沉降的比值,即Us=~w(z,t)/~w(z,∞).當加荷到一定時間,土體固結(jié)趨于穩(wěn)定,在循環(huán)荷載作用下所能達到的最大沉降趨近于一固定值,即認為此最大沉降為土體固結(jié)的最終沉降~w(z,∞).

3.1 正弦荷載

圖6為ka/kw=100時正弦荷載作用下超孔隙氣壓力、超孔隙水壓力和固結(jié)度隨時間的變化規(guī)律.由圖可以看出,超孔隙氣壓力、超孔隙水壓力的消散和土體的固結(jié)均隨荷載呈周期性變化.超孔隙氣壓力在第一個周期內(nèi)就消散到趨近于0,然后在0附近上下振動;超孔隙水壓力在第一個周期內(nèi)消散到0.2,然后在0.2附近上下振動,形成一段“平臺期”,最后又消散到0,并在0附近上下振動;固結(jié)度在第一個周期前半段一直增大,隨后上下振動(整體呈現(xiàn)逐漸固結(jié)趨勢),最后在Us=1附近上下振動.

圖6正弦荷載作用下非飽和黏彈性地基的固結(jié)特性隨時間的變化(ka/kw=100)Fig.6 Change in consolidation characteristics of viscoelastic unsaturated soils for sinusoidal loading(ka/kw=100)

圖7 為正弦荷載作用下ka/kw不同時超孔隙氣壓力、超孔隙水壓力及固結(jié)度隨時間的變化規(guī)律.由圖可以看出,ka/kw的大小對超孔隙氣壓力的消散影響非常明顯.當ka/kw較大時,超孔隙氣壓力在第一個周期就消散到趨近于0,之后在0附近上下振動;ka/kw越大,超孔隙水壓力達到的峰值越小,超孔隙水壓力消散越快;當ka/kw較小時,超孔隙氣壓力和超孔隙水壓力的消散曲線的切線斜率幾乎相同;ka/kw越大,固結(jié)度在周期內(nèi)達到的峰值越大,固結(jié)速度越快.

圖7 正弦荷載作用下ka/kw不同時非飽和黏彈性地基的固結(jié)特性隨時間的變化Fig.7 Change in consolidation characteristics of viscoelastic unsaturated soils for sinusoidal loading at diあerent values of ka/kw

3.2 三角形荷載

圖8三角形荷載作用下ka/kw不同時非飽和黏彈性地基的固結(jié)特性隨時間的變化Fig.8 Change in consolidation characteristics of viscoelastic unsaturated soils for triangle loading at diあerent values of ka/kw

圖8 為三角形荷載作用下ka/kw不同時超孔隙氣壓力、超孔隙水壓力及土層沉降隨時間的變化.由圖可以看出,三角形荷載作用下超孔隙氣壓力、超孔隙水壓力和固結(jié)度的變化規(guī)律與正弦荷載作用下的規(guī)律基本相同.ka/kw較大時,超孔隙氣壓力很快消散到0,然后在0附近上下振蕩;ka/kw越大,超孔隙氣、水壓力消散越快,超孔隙氣、水壓力的峰值越小.超孔隙水壓力雖不類似超孔隙氣壓力很快消散至0附近,但整體呈消散趨勢(中心線向下傾斜).另外,由圖8(c)可以看出,ka/kw越大,固結(jié)速度越快.

3.3 梯形荷載

圖9為梯形荷載作用下ka/kw不同時超孔隙氣壓力、超孔隙水壓力及固結(jié)度隨時間的變化規(guī)律.由圖可以看出,梯形荷載作用下的變化規(guī)律與正弦荷載、三角形荷載作用下的變化規(guī)律類似.另外,在每個周期的中間段,均有一個帶凹陷的“平臺期”,而在“平臺期”的兩端各有一個凸起,這是由于土體黏滯性的存在,土中孔隙壓力的消散尚未適應荷載的突變引起的.另外,ka/kw越大,超孔隙氣壓力和超孔隙水壓力消散越快,土體固結(jié)速度越快.

圖9 梯形荷載作用下ka/kw不同時非飽和黏彈性地基的固結(jié)特性隨時間的變化Fig.9 Change in consolidation characteristics of viscoelastic unsaturated soils for trapezoidal loading at diあerent values of ka/kw

4 結(jié)束語

本工作基于已有的非飽和黏彈性地基一維固結(jié)理論,應用Merchant黏彈性模型,運用Laplace變換和Cayley-Hamilton方法,求得了正弦荷載、三角形荷載和梯形荷載3種情況下的非飽和黏彈性地基一維固結(jié)時超孔隙氣壓力、超孔隙水壓力以及土層沉降的解;然后,采用Crump及Durbin方法實現(xiàn)Laplace逆變換,獲得了半解析解;最后,通過算例分析了3種循環(huán)荷載作用下非飽和黏彈性地基的一維固結(jié)特性.研究結(jié)果表明:①循環(huán)荷載作用下超孔隙氣壓力、超孔隙水壓力和固結(jié)度的變化與荷載類似,呈周期性變化;②循環(huán)荷載作用下非飽和黏彈性地基固結(jié)時的超孔隙氣壓力和超孔隙水壓力的變化規(guī)律基本相同,當ka/kw越大時,超孔隙氣壓力越快消散到0,超孔隙氣壓力和超孔隙水壓力消散越快;③循環(huán)荷載作用下非飽和土黏彈性地基的固結(jié)過程中,ka/kw越大,固結(jié)速度越快;④循環(huán)荷載作用下超孔隙氣壓力、超孔隙水壓力和固結(jié)度呈上下振動趨勢,但從整體上看,變化曲線的中心線向下傾斜,超孔隙壓力越來越小(一直在消散),而固結(jié)度越來越大(土體一直在固結(jié)).

[1]BLIGHT G E.Strength and consolidation characteristics of compacted soils[D].London:University of London,1961.

[2]SCOTT R F.Principles of soil mechanics[M].Massachusetts:Addison Wesley Publishing Company,1963.

[3]BARDEN L.Consolidation of compacted and unsaturated clays[J].Geotechnique,1965,15(3):267-286.

[4]FREDLUND D G,HASAN J U.One-dimensional consolidation theory:unsaturated soils[J].Canadian Geotechnical Journal,1979,17(3):521-531.

[5]FREDLUND D G,RAHARDjO H.Soil mechanics for unsaturated soils[M].New York:John Wiley and Sons,1993.

[6]王奎華,謝康和,曾國熙.雙面半透水邊界的一維黏彈性固結(jié)理論[J].巖土工程學報,1998,20(2):34-36.

[7]蔡袁強,徐長節(jié),袁海明.任意荷載下成層黏彈性地基的一維固結(jié)[J].應用數(shù)學和力學,2001,22(3):307-313.

[8]蔡袁強,梁旭,鄭灶鋒,等.半透水邊界的黏彈性土層在循環(huán)荷載下的一維固結(jié)[J].土木工程學報,2003,36(8):86-90.

[9]謝康和,周瑾,董亞欽.循環(huán)荷載作用下地基一維非線性固結(jié)解析解[J].巖土力學與工程學報,2006,25(1):21-26

[10]SHAN Z D,LING D S,DING H J,et al.Analytical solution for the 1D consolidation of unsaturated multi-layered soil[J].Computers and Geotechnics,2014,57:17-23.

[11]ZHOU W H,ZHAO L S.One-dimensional consolidation of unsaturated soil subjected to timedependent loading with various initial and boundary conditions[J].International Journal of Geomechanics,2014,14(2):291-301.

[12]LIEM H,BEHZAD F.Analytical solution for the two-dimensional plane strain consolidation of an unsaturated soil stratum subjected to time-dependent loading[J].Computers and Geotechnics,2015,67:1-16.

[13]QIN A F,CHEN G J,TAN Y W,et al.Analytical solution to one dimensional consolidation in unsaturated soils[J].Applied Mathematics and Mechanics(English Edition),2008,29(10):1329-1340.

[14]QIN A F,SUN D A,Tan Y W.Analytical solution to one-dimensional consolidation in unsaturated soils under loading changed exponentially with time[J].Computers and Geotechnics,2010,37(1):233-238.

[15]秦愛芳,羅坤,孫德安.非飽和土黏彈性地基一維固結(jié)特性分析[J].上海大學學報(自然科學版),2010,16(2):203-209.

[16]KORN G A,KORN T M.Mathematical handbook for scientists and engineerings[M].New York:McGraw-Hill Inc,1968.

[17]DURBIN F.Numerical inversion of the Laplace transforms:an eきcient improvement to Dubner and Abate’s method[J].The Computer Journal,1974,17(4):371-376.