任意荷載下連續(xù)排水邊界分數(shù)階黏彈性地基一維固結(jié)模型

王玨　童立紅　金立　徐長節(jié)

摘 要：以基于Caputo分數(shù)階導數(shù)的彈壺元件修正Kelvin模型來描述飽和土體一維固結(jié)的力學行為，并引入連續(xù)排水邊界條件，通過Laplace變換，聯(lián)立求解得到任意荷載下連續(xù)排水邊界分數(shù)階黏彈性地基有效應力及固結(jié)沉降的解析解。采用Laplace逆變換，獲得了其時域內(nèi)的理論解，并分析了梯形循環(huán)荷載及施工荷載作用下相關(guān)參數(shù)對固結(jié)沉降的影響。研究結(jié)果表明：循環(huán)荷載作用下，黏土地基的沉降變化呈振蕩增長，且振蕩幅值隨著邊界透水性的增大而增大;分數(shù)階次α增大，使固結(jié)前期沉降速率減慢，而在固結(jié)后期，α值對沉降的影響正好相反;循環(huán)荷載下沉降變化曲線的振蕩幅值隨著分數(shù)階次α的增大而減小。此外，一維固結(jié)沉降的發(fā)展還與土體力學參數(shù)及荷載參數(shù)相關(guān)，彈性模量E越大，最終沉降量越小;黏彈性體的延遲時間F越大，固結(jié)沉降變化越慢。

關(guān)鍵詞：任意荷載;連續(xù)排水邊界;分數(shù)階導數(shù);黏彈性;一維固結(jié)

中圖分類號：TU431 文獻標志碼：A 文章編號：2096-6717（2020）01-0056-08

Abstract：The Kelvin constitutive model is modified by the spring-pot element based on the Caputo fractional derivative to describe the mechanical behavior of one-dimensional consolidation of saturated soil. After introducing the continuous drainage boundary condition， the analytical solutions of the effective stress and the settlement under time-dependent loading are derived by performing Laplace transformation. The Laplace inverse transformation is used to obtain the theoretical solutions in time domain， and the influences of relevant parameters on the settlement under trapezoidal cyclic loading and construction loading are studied. The results show that the settlement of viscoelastic soil under cyclic loading increases in an oscillating manner， and the amplitude of the oscillation increases with the boundary permeability. A higher value of the fractional order α slows the development of settlement in the early stage of consolidation. However， in the later stage of consolidation， the effect of α on settlement is reversed. The oscillation amplitude of the settlement under cyclic loading decreases with increase of α. Furthermore， detailed analysis indicates that the development of one-dimensional consolidation settlement is also related to mechanical properties of soil and loading parameters. The larger the elastic modulus E is， the smaller the final settlement， and the greater the delay time of viscoelastic is， the slower the settlement occurs.

Keywords：time-dependent loading; continuous drainage boundary; fractional order derivative; viscoelastic; one-dimensional consolidation

在Terzaghi固結(jié)理論中，土體被處理為線彈性模型，而流變特性是軟土的一種重要的工程特性[1]。因此，考慮軟黏土的流變特性，將土體視為黏彈性介質(zhì)通常更符合實際工程[2]。Taylor等[3]首先引入Kelvin模型來描述土骨架的黏彈性變形;Tan[4]基于Maxwell模型對受側(cè)限土體的固結(jié)和滯流進行了研究。此后，金問魯?shù)萚5]提出了固結(jié)方程的一個近似解法，并給出了各種條件下簡單問題的解答;趙維炳[6]基于廣義Voigt模型，推導了飽和土體一維固結(jié)問題的普遍理論解答;Xie等[7-8]引入Merchant模型及四元件模型到固結(jié)理論中，分析了軟黏土的固結(jié)特性;蔡袁強等[9]求解了任意荷載下成層粘彈性地基一維變形問題。然而，上述經(jīng)典流變模型不能很好地與實驗數(shù)據(jù)相吻合[10]，主要是由于整數(shù)階微分算子的性質(zhì)決定了經(jīng)典流變模型的核函數(shù)通常是指數(shù)函數(shù)的組合，欲精確描述實驗數(shù)據(jù)，常常不得不取消高階的微分項或者以降低本構(gòu)模型的應用范圍為代價[11]。

Gement[12]首先提出了黏彈性材料的分數(shù)階導數(shù)本構(gòu)模型，而后一些學者將其引入到固結(jié)理論中，并指出分數(shù)階導數(shù)流變模型可以有效克服經(jīng)典模型的缺點。Koeller[13]用基于分數(shù)導數(shù)的彈壺元件替換牛頓黏壺，研究分析了多種模型的流變特性;孫海忠等[14]采用含分數(shù)導數(shù)的Kelvin模型對珠江三角洲南沙地區(qū)典型軟土的流變試驗數(shù)據(jù)進行擬合，得到很好的結(jié)果;Yin等[15]對分數(shù)階軟土蠕變過程中的力學性能進行了系統(tǒng)的研究;汪磊等[16]基于分數(shù)階導數(shù)理論引入Kelvin-Voigt模型，獲得了任意荷載情況下一維固結(jié)問題的半解析解;劉忠玉等[17]求得了恒載下基于分數(shù)階Kelvin模型飽和軟黏土一維固結(jié)理論解，并通過對比一維流變固結(jié)試驗曲線及整數(shù)階模型理論曲線，指出基于分數(shù)階Kelvin模型模擬的孔壓消散曲線更接近試驗曲線。

另一方面，實際工程中土體的邊界往往是處于透水與不透水之間的一種中間狀態(tài)[18]。蔡袁強等[19]、汪磊等[20]研究了半透水邊界條件下一維固結(jié)問題。但是半透水邊界計算相對復雜，且不能嚴格滿足初始條件，限制了土體固結(jié)方程解的適用性[21]。基于此，梅國雄等[18]提出了一個從不透水到透水的雙面不對稱連續(xù)排水邊界。目前，關(guān)于變荷載、連續(xù)排水邊界及分數(shù)階導數(shù)黏彈性模型耦合的一維固結(jié)理論分析很少見諸于文獻。筆者針對Caputo分數(shù)階導數(shù)的彈壺元件修正Kelvin模型黏彈性地基，引入連續(xù)排水邊界條件，推導了任意荷載下連續(xù)排水邊界分數(shù)階黏彈性地基一維固結(jié)方程的半解析解，并分析了相關(guān)參數(shù)對軟黏土固結(jié)沉降特性的影響。

4 結(jié) 論

基于Caputo分數(shù)階導數(shù)的彈壺元件修正Kelvin模型，引入連續(xù)排水邊界條件，利用Laplace變換求得考慮連續(xù)排水邊界條件時分數(shù)階導數(shù)黏彈性地基在任意隨時間變化的荷載下有效應力及沉降的解析解，運用Laplace逆變換得到其時域內(nèi)的數(shù)值解。通過系統(tǒng)的算例分析，可以得到如下結(jié)論：

1）循環(huán)荷載作用下，黏土地基的沉降變化呈振蕩增長，但滯后于荷載的變化，且振蕩幅值隨著邊界透水性的增大而增大。

2）分數(shù)階次α增大，使固結(jié)前期沉降發(fā)展速率減慢，但在固結(jié)后期，α值對沉降的影響正好相反，最終固結(jié)沉降達到穩(wěn)定的時間隨著α的增大而縮短。另外，隨著分數(shù)階次α的增大，循環(huán)荷載下沉降變化曲線的振蕩幅值明顯減小。

3）分數(shù)階黏彈性地基一維固結(jié)沉降的發(fā)展還與土體力學參數(shù)及荷載參數(shù)相關(guān)。彈性模量E越大，最終沉降量越小，固結(jié)沉降達到穩(wěn)定的時間越短，且循環(huán)荷載下固結(jié)沉降的振蕩幅值越小;黏彈性體的延遲時間F越大，固結(jié)沉降變化速率越慢。

參考文獻：

[1] 潘曉東， 蔡袁強， 徐長節(jié)， 等. 循環(huán)荷載下黏彈性飽和土層的一維固結(jié)[J]. 巖土力學， 2006， 27（2）： 272-276.

PAN X D， CAI Y Q， XU C J， et al. One-dimensional consolidation of viscoelastic soil layer under cyclic loading [J]. Rock and Soil Mechanics， 2006， 27（2）： 272-276.（in Chinese）

[2] 解益， 李培超， 汪磊， 等. 分數(shù)階導數(shù)黏彈性飽和土體一維固結(jié)半解析解[J]. 巖土力學， 2017， 38（11）： 3240-3246.

XIE Y， LI P C， WANG L， et al. Semi-analytical solution for one-dimensional consolidation of viscoelastic saturated soil with fractional order derivative [J]. Rock and Soil Mechanics， 2017， 38（11）： 3240-3246.（in Chinese）

[3] TAYLOR D W， MERCHANT W. A theory of clay consolidation accounting for secondary compression [J]. Journal of Mathematics and Physics， 1940， 19（1-4）： 167-185.

[4] TAN T K. Secondary time effects and consolidation of clays [J]. Science China Mathematics， 1958， 1（11）： 1060-1075.

[5] 金問魯， 門福錄. 粘土固結(jié)與次時間效應單維問題的近似解[J]. 水利學報， 1964（2）： 69-74.

JIN W L， MENG F L. An approximate solution to the unidimensional problem of clay consolidation and sub-time effects [J]. Journal of Hydraulic Engineering， 1964（2）： 69-74.（in Chinese）

[6] 趙維炳. 廣義Voigt模型模擬的飽水土體一維固結(jié)理論及其應用[J]. 巖土工程學報， 1989， 11（5）： 78-85.

ZHAO W B. Gen-dimensional consolidation theory of saturated soil with generalized Voigt model simulation and its application [J]. Journal of Geotechnical Engineering， 1989， 11 （5）： 78-85.（in Chinese）

[7] XIE K H， LIU X W. A study on one dimensional consolidation of soils exhibiting rheological characteristics[C]// Proceedings of International Symposium on Compression and Consolidation of Clayey Soils， Rotterdam， 1995： 388.

[8] XIE K H， XIE X Y， LI X B. Analytical theory for one-dimensional consolidation of clayey soils exhibiting rheological characteristics under time-dependent loading[J]. International Journal for Numerical and Analytical Methods in Geomechanics， 2008， 32（14）： 1833-1855.

[9] 蔡袁強， 徐長節(jié)， 袁海明. 任意荷載下成層粘彈性地基的一維固結(jié)[J]. 應用數(shù)學和力學， 2001， 22（3）： 307-313.CAI Y Q， XU C J， YUAN H M. One-dimensional consolidation of layered viscoelastic foundations under arbitrary loading [J]. Applied Mathematics and Mechanics， 2001， 22（3）： 307-313.（in Chinese）

[10] 張為民. 一種采用分數(shù)階導數(shù)的新流變模型理論[J]. 湘潭大學（自然科學學報）， 2001， 23（1）： 30-36.

ZHANG W M. A new rheologicai model theory with fractional order derivative [J]. Natural Science Journal of Xiangtan University， 2001， 23（1）： 30-36.（in Chinese）

[11] 劉林超， 張衛(wèi). 具有分數(shù)Kelvin模型的粘彈性巖體中水平圓形硐室的變形特性[J]. 巖土力學， 2005， 26（2）： 121-123.

LIU L C， ZHANG W. Deformation properties of horizontal round adits in viscoelastic rocks with fractional Kelvin model [J]. Rock and Soil Mechanics， 2005， 26（2）： 121-123.（in Chinese）

[12] GEMANT A. A method of analyzing experimental results obtained from elasto-viscous bodies [J]. Physics， 1936， 7（8）： 311-317.

[13] KOELLER R C. Applications of fractional calculus to the theory of viscoelasticity [J]. Journal of Applied Mechanics， 1984， 51（2）： 299-307.

[14] 孫海忠， 張衛(wèi). 一種分析軟土黏彈性的分數(shù)導數(shù)開爾文模型[J]. 巖土力學， 2007， 28（9）： 1983-1986.SUN H Z， ZHANG W. Analysis of soft soil withviscoelastic fractional derivative Kelvin model [J]. Rock and Soil Mechanics， 2007， 28（9）： 1983-1986.（in Chinese）

[15] YIN D S， LI Y Q， WU H， et al. Fractional description of mechanical property evolution of soft soilsduring creep [J]. Water Science and Engineering， 2013， 6（4）： 446-455.

[16] 汪磊， 孫德安， 解益， 等. 任意荷載下分數(shù)階導數(shù)黏彈性飽和土體一維固結(jié)[J]. 巖土工程學報， 2017， 39（10）： 1823-1831.

WANG L， SUN D A， XIE Y， et al. One-dimensional consolidation of fractional order derivative viscoelastic saturated soils under arbitrary loading [J]. Journal of Geotechnical Engineering， 2017， 39（10）： 1823-1831.（in Chinese）

[17] 劉忠玉， 楊強. 基于分數(shù)階Kelvin模型的飽和黏土一維流變固結(jié)分析[J]. 巖土力學， 2017， 38（12）： 3680-3687， 3697.

LIU Z Y， YANG Q. One-dimensional rheological consolidation analysis of saturated clay using fractional order Kelvins model [J]. Rock and Soil Mechanics， 2017， 38（12）： 3680-3687， 3697.（in Chinese）

[18] 梅國雄， 夏君， 梅嶺. 基于不對稱連續(xù)排水邊界的太沙基一維固結(jié)方程及其解答[J]. 巖土工程學報， 2011， 33（1）： 28-31.

MEI G X， XIA J， MEI L. Terzaghis one-dimensional consolidation equation and its solution based on asymmetric continuous drainage boundary [J]. Chinese Journal of Geotechnical Engineering， 2011， 33（1）： 28-31.（in Chinese）

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

CAI Y Q， LIANG X， ZHENG Z F， et al. One-dimensional consolidation of viscoelastic soil layer with semi-permeable boundary under cyclic loadings [J]. Chinese Journal of Civil Engineering， 2003， 36（8）： 86-90.（in Chinese）

[20] 汪磊， 李林忠， 徐永福， 等. 半透水邊界下分數(shù)階黏彈性飽和土一維固結(jié)特性分析[J]. 巖土力學， 2018， 39（11）： 4142-4148.

WANG L， LI L Z， XU Y F， et al. Analysis of one‐dimensional consolidation of fractional viscoelastic saturated soils with semi-permeable boundary [J]. Rock and Soil Mechanics， 2018， 39（11）： 4142-4148.（in Chinese）

[21] 宋劍鵬， 謝晨， 朱洵， 等. 廣義透水邊界下壩基非線性一維固結(jié)分析[J]. 人民黃河， 2018， 40（9）： 100-103， 107.

SONG J P， XIE C， ZHU X， et al. Generalized drainage boundary applied in dam foundation with one-dimensional non-linear consolidatio [J]. Yellow River， 2018， 40（9）： 100-103， 107.（in Chinese）

[22] DURBIN F. Numerical inversion of Laplace transforms： an efficient improvement to Dubner and Abates method [J]. The Computer Journal， 1974， 17（4）： 371-376.

（編輯 胡玲）