糾永志，劉忠玉，樂金朝，孫麗云

（1.同濟(jì)大學(xué) 巖土及地下工程教育部重點(diǎn)實(shí)驗(yàn)室，上海200092；2.同濟(jì)大學(xué) 地下建筑與工程系，上海200092；3.鄭州大學(xué) 土木工程學(xué)院，河南 鄭州450001；4.鄭州大學(xué) 水利與環(huán)境學(xué)院，河南 鄭州450001）

飽和軟黏土非線性1維固結(jié)理論的研究始于上世紀(jì)60年代.Davis和 Raymond等［1－2］分別求得了特定情況下的1維非線性問題的解析解.而對(duì)于一般情況通常難以得到解析解，因此Chen和Zhou等［3－4］分別采用有限差分或有限體積法等對(duì)土體1維非線性固結(jié)理論進(jìn)行了求解.不過這些理論都建立在假定地基中的初始有效應(yīng)力為恒值（即初始有效應(yīng)力沿土體深度是不變的）的基礎(chǔ)上的，但是現(xiàn)實(shí)土體中存在的是沿深度變化的自重應(yīng)力.Gibson等［5］的研究表明，當(dāng)土層較厚時(shí)，固結(jié)土體的自重應(yīng)力將會(huì)對(duì)其固結(jié)過程產(chǎn)生很大的影響，且考慮自重應(yīng)力所得出的土體的固結(jié)速率快于Terzaghi理論值.在計(jì)算地基最終沉降量時(shí)，已有一些方法體現(xiàn)了自重應(yīng)力水平的影響，比如，丁洲祥等［6］在修正割線模量法時(shí)就引入了修正系數(shù)α，并指出該系數(shù)有隨埋深增加而逐漸減小的規(guī)律.所以有必要把土體自重應(yīng)力沿土體深度的實(shí)際分布的因素引入到非線性固結(jié)理論的研究工作中.因此，竇宜等［7］用多種土體進(jìn)行了室內(nèi)離心試驗(yàn)，文獻(xiàn)［8－9］用數(shù)值方法在考慮土體自重應(yīng)力的基礎(chǔ)上對(duì)特定情況進(jìn)行了求解.但是上述對(duì)Terzaghi 1維固結(jié)理論的改進(jìn)都是基于飽和黏土中滲流符合Darcy定律的假定之上的.

在孔隙水滲流方面，Darcy定律因其簡(jiǎn)潔而被廣泛用于多種飽和土，但其對(duì)滲透性較低的致密黏土以及某些軟土的有效性卻一直受到質(zhì)疑.Gibson等［10］指出Darcy定律是有限適用的，土體的滲透與土骨架形狀的改變以及超靜孔隙水壓力減少的梯度有關(guān).另外，文獻(xiàn)［11－13］在滲透試驗(yàn)中都發(fā)現(xiàn)了相關(guān)黏土中的滲流都不同程度地存在著對(duì)Darcy定律偏離的現(xiàn)象，并分別提出了不同形式的非Darcy滲流數(shù)學(xué)擬合式.其中，應(yīng)用較廣的當(dāng)屬 Hansbo［14］公式

式中：q為滲透速度；m為由實(shí)驗(yàn)確定的常數(shù)；k，K分別為指數(shù)形式和線性關(guān)系表達(dá)式中的滲透系數(shù)；i為水力梯度；i1為直線滲流起始水力梯度；i0為直線滲流計(jì)算起始水力梯度.

曲線在i＝i1處是連續(xù)的，因此i0，k滿足下式：

式（2）中，如令m＝1，則i0＝0，k＝K，式（1）可退化為Darcy定律的表達(dá)形式.如忽略式（1）中的曲線段，該式中的i0可視為滲流起始水力梯度［15］.劉忠玉等［16－17］基于忽略式（1）曲線段的直線型非Darcy滲流對(duì)Terzaghi 1維固結(jié)理論進(jìn)行了修正，強(qiáng)調(diào)了起始水力梯度的影響.謝海瀾、劉忠玉、鄂建等［18－21］采用指數(shù)－直線型非 Darcy滲流模型式（1）對(duì)Terzaghi 1維固結(jié)方程進(jìn)行了修正，并分別采用了半解析法、有限體積法和有限差分法對(duì)相關(guān)固結(jié)方程進(jìn)行了求解.這些研究雖然都不同程度地考慮了滲流的非Darcy特性，但沒考慮固結(jié)過程中土體壓縮變形的非線性和土體初始有效應(yīng)力的實(shí)際分布.

本文擬把式（1）所示的 Hansbo［14］非Darcy滲流模型引入到考慮土體自重應(yīng)力的1維非線性固結(jié)理論中，進(jìn)一步探討土體自重應(yīng)力和土體滲流非Darcy特性對(duì)正常固結(jié)飽和軟黏土非線性固結(jié)的影響.

1 問題的描述

本文的分析模型如圖1所示，均質(zhì)正常固結(jié)飽和黏土層的厚度為H，頂面為透水面，底面為不透水面，假定頂面瞬時(shí)加載且該土層在自重應(yīng)力作用下已經(jīng)完成固結(jié).對(duì)于正常固結(jié)飽和黏土其初始有效應(yīng)力可取其自重應(yīng)力σ′cz.采用文獻(xiàn)［18－20］相似的推導(dǎo)過程可得

圖1 分析模型Fig.1 Schematic diagram of analytical model

由e—lgσ′和e—lgK的經(jīng)驗(yàn)關(guān)系式［1－2，8－9］可得

式（3）—（5）中：u＝u（z，t）為t時(shí)刻距黏土層頂面深度z處的孔壓；γw為水的重度；e0和K0分別與某一應(yīng)力狀態(tài)σ′0相對(duì)應(yīng)的孔隙比和滲透系數(shù)；e和K分別為與某時(shí)刻土層任意深度z處的有效應(yīng)力σ′（z，t）相應(yīng)的孔隙比和滲透系數(shù)；cc為壓縮指數(shù)；ck為滲透指數(shù).

在瞬時(shí)均布荷載p0的作用下深度z處總應(yīng)力為

式中：σ′cz0為壓縮土層頂面的初始有效應(yīng)力，其大小反映了土體的埋深；γ′為土的有效重度.

由式（6）可得

通過對(duì)式（7）兩邊關(guān)于z求偏導(dǎo)可得

將式（4），（5）和（8）代入式（3）可得控制方程為

其初始條件和邊界條件分別為

將式（9）量綱一化可得

這里，P≠0意味著考慮土體的自重，P＝0即不考慮自重.同時(shí)，I可改寫為

相應(yīng)的初始條件和邊界條件為

式中：d＝σ′cz0／σ′0，可反映土體的埋深；Q＝p0／σ′0為量綱一化后的土體頂面荷載.

這樣，式（12），（13）—（15）構(gòu)成封閉方程組，通過迭代求解便可得j＋1時(shí)刻任一深度處的有效應(yīng)力.需要注意的是，由于和中含有，因此迭代格式（12）為隱式格式，可從透水邊界開始逐步向不透水邊界移動(dòng)迭代求解.

當(dāng)量綱一化后的有效應(yīng)力F求出后，由式（7）可得量綱一化的超孔隙水壓

2 方程的有限差分求解

令

則方程（10）變?yōu)?/p>

目前方程（11）尚無法求得其解析解，因此本文用有限差分法來對(duì)此固結(jié)方程進(jìn)行求解.即以ΔZ從上到下將土層離散為n個(gè)薄層，同時(shí)以步長ΔT對(duì)時(shí)間進(jìn)行離散，這樣方程（11）離散為

式中：λ＝ΔT／ΔZ2；i為空間節(jié)點(diǎn)數(shù)，i＝1，2，…，n－1；j為空間節(jié)點(diǎn)數(shù)，j＝0，1，2，….當(dāng)1≤i≤n－1時(shí)

相應(yīng)的初始條件和頂面邊界條件離散為

考慮到底面不排水，此時(shí)最底層的排水量應(yīng)等于體積壓縮量，經(jīng)推導(dǎo)可得

這樣可以獲得如下按孔壓和變形定義的兩種平均固結(jié)度公式：

式中：Up，Us分別為量綱一化的孔壓固結(jié)度和變形固結(jié)度；e1，ef和e分別為第i層土初始孔隙比、最終孔隙比和j時(shí)刻的孔隙比.

3 解法驗(yàn)證

當(dāng)c＝1，土體頂面自重應(yīng)力d＝1，且不考慮土體固結(jié)過程中的非Darcy滲流特性，并令土體初始有效應(yīng)力沿深度不變，即P＝0時(shí)，本文模型便退化為Davis和Raymond固結(jié)方程［1］.按前述有限差分法對(duì)符合Davis和Raymond假定的特例計(jì)算了平均固結(jié)度Up（計(jì)算中，取ΔZ＝0.02，ΔT＝10－6，迭代精度為10－6），結(jié)果示于圖2.為便于對(duì)比，圖中也給出了Davis和Raymond解析解結(jié)果［1］.很明顯，本文有限差分法計(jì)算結(jié)果曲線與解析解之間的誤差極小，這證明了本文解法的有效性.

圖2 平均孔壓固結(jié)度Up與時(shí)間因子T關(guān)系曲線（d＝1，c＝1，m＝1，I1＝0）Fig.2 Relationship between Upand T（d＝1，c＝1，m＝1，I1＝0）

4 計(jì)算與分析

4.1 非Darcy滲流參數(shù)對(duì)土體固結(jié)的影響

首先考察在考慮自重應(yīng)力情況下直線滲流起始水力梯度i1的影響.圖3和4分別繪出了當(dāng)c＝0.5和1.2、m＝2.0以及不同I1（即不同i1值）值時(shí)考慮土體自重時(shí)的Up和Us隨時(shí)間的變化曲線.由圖中可以看出，在考慮土體自重的情況下，非Darcy滲流特性和不考慮自重時(shí)［21］相似，即基于非Darcy滲流的孔壓消散和地基沉降速率都明顯慢于基于Darcy滲流的情形（圖中I1＝0的曲線）.

圖3 I1值對(duì)平均孔壓固結(jié)度Up的影響Fig.3 Influence of I1on degree of consolidation Up

為了考察m值的影響，圖5和圖6分別繪出了當(dāng)c＝0.5和c＝1.2時(shí)基于非Darcy滲流（I1＝1.0但不同m值）且考慮土體自重時(shí)的孔壓消散和地基沉降隨時(shí)間變化曲線.很明顯，m值對(duì)孔壓消散速率和地基沉降速率的影響，與I1相似.即m值越大孔壓消散速率和地基沉降速率都越慢，并且隨著時(shí)間的發(fā)展非Darcy滲流參數(shù)的影響越來越顯著.

4.2 考慮非Darcy滲流時(shí)埋深對(duì)固結(jié)度的影響

基于Darcy滲流時(shí)埋深對(duì)固結(jié)度的影響，馬崇武［9］已做過探討，本文只討論非Darcy滲流情況下埋深對(duì)固結(jié)度的影響.

為了分析在非Darcy滲流條件下并考慮土體自重時(shí)土層埋深對(duì)土體的孔壓消散和地基沉降速率的影響，本文分別取c＝0.8和c＝1.5時(shí)計(jì)算了不同埋深時(shí)Up和Us隨時(shí)間的變化曲線，如圖7和圖8所示.這里的結(jié)論與文獻(xiàn)［9］相似，即埋深越大，地基沉降速度就越慢，并且c值越大，土體埋深對(duì)地基沉降速度影響就越大；但土體埋深對(duì)孔壓消散速率的影響規(guī)律與其對(duì)地基沉降的影響不同：當(dāng)c值較小時(shí)（例如c＝0.8）土體的埋深越大，孔壓消散速度越快；當(dāng)c值較大時(shí)（例如c＝1.5），c值對(duì)孔壓消散速度的影響體現(xiàn)出和孔壓較小時(shí)不同的規(guī)律：即在固結(jié)的前期，埋深越大，孔壓消散速度越快，而在固結(jié)的后期，情況正好相反.但和文獻(xiàn)［9］不同的是，由于滲流非Darcy特性的影響，孔壓消散和地基沉降速度都明顯慢于符合Darcy滲流時(shí)的孔壓消散和地基沉降速率.

圖6 m值對(duì)平均變形固結(jié)度Us的影響Fig.6 Influence of mon degree of consolidation Us

4.3 土體自重及非Darcy滲流特性對(duì)土體固結(jié)的影響

為了系統(tǒng)探討土體變形的非線性、土體自重應(yīng)力、土體滲流的非Darcy特性等對(duì)孔壓消散和地基沉降的影響，圖9和圖10分別給出了c＝0.5和c＝1.2時(shí)基于Darcy滲流或非Darcy滲流、考慮土體自重應(yīng)力與不考慮土體自重應(yīng)力時(shí)的孔壓消散曲線和地基沉降曲線.

從圖中可以看出，不論是否考慮非Darcy滲流，考慮土體自重應(yīng)力時(shí)的孔壓消散速率和地基沉降速率都要大于不考慮土體自重時(shí)的情況，這與文獻(xiàn)［7］的結(jié)論相似.另外，不論是否考慮土體自重應(yīng)力，土體滲流的非Darcy特性都延緩了孔壓消散速率和地基沉降速率，這與文獻(xiàn)［18－21］相似.但綜合考慮土體自重應(yīng)力和非Darcy滲流特性的影響時(shí)，這里得出了與上述結(jié)論都有所不同的結(jié)果.由圖10可以看出，在固結(jié)的前期，考慮非Darcy滲流且考慮土體自重（圖中P＝1.0，m＝2.0，I1＝1.0的曲線）的地基沉降速率要快于考慮Darcy滲流且不考慮土體自重（圖中P＝0，Darcy的曲線）的地基沉降速率，但在固結(jié)后期則出現(xiàn)相反的情況，即考慮非Darcy滲流且考慮土體自重的孔壓消散速率逐漸開始慢于考慮Darcy滲流但不考慮土體自重的地基沉降速率.同時(shí)，圖9所示的孔壓消散速率的規(guī)律與之相似.因此不考慮滲流非Darcy特性和土體自重的固結(jié)分析在土體固結(jié)的前期低估了孔壓消散速率和地基沉降速率，而在固結(jié)后期則正好相反.

由圖9可以看出，土體自重對(duì)固結(jié)前期的孔壓消散無明顯的影響，但隨著時(shí)間的發(fā)展考慮土體自重的孔壓消散逐漸快于不考慮土體自重的情形，但到固結(jié)的后期，基于Darcy滲流考慮或不考慮土體自重的孔壓消散曲線逐漸趨于一致，基于非Darcy滲流時(shí)也是如此，因此可以認(rèn)為，土體自重對(duì)孔壓消散的影響主要表現(xiàn)在固結(jié)的中期，而對(duì)前期和后期無明顯影響，而土體滲流的非Darcy特性對(duì)孔壓消散的影響主要表現(xiàn)在中后期，對(duì)前期無明顯影響.

由圖10可以看出，考慮土體自重時(shí)的地基沉降速率一開始就明顯快于不考慮土體自重時(shí)的地基沉降速率，但到固結(jié)的后期二者逐漸趨于一致.因此可以認(rèn)為，土體自重對(duì)地基沉降的影響主要表現(xiàn)在前期和中期，而對(duì)后期的影響較??；土體滲流的非Darcy特性對(duì)地基沉降的影響主要表現(xiàn)在中后期，而對(duì)前期的影響較小.

4.4 和Terzaghi解的對(duì)比

為了和Terzaghi解進(jìn)行對(duì)比，圖9和圖10中以實(shí)線標(biāo)出了Terzaghi解的固結(jié)曲線.Gibson等［5］認(rèn)為考慮自重應(yīng)力所得出的土體的固結(jié)速率快于Terzaghi理論值，但這里通過綜合考慮非達(dá)西滲流、土體自重應(yīng)力和土體非線性變形時(shí)得出了與此不同的結(jié)論.由圖9可以看出，當(dāng)c值較小時(shí)（例如c＝0.5），在考慮土體自重的情況下，基于Darcy滲流（圖中P＝1.0，Darcy的曲線）的孔壓消散速率快于Terzaghi解，這點(diǎn)與文獻(xiàn)［5］相同；但基于非Darcy滲流（圖中P＝1.0，m＝2.0，I1＝1.0的曲線）的解答，則可能出現(xiàn)在固結(jié)的中期孔壓消散速率快于Terzaghi解，而在固結(jié)后期，由于非Darcy特性的影響其孔壓消散速率慢于Terzaghi解的情況.當(dāng)c值較大時(shí)（例如c＝1.2），在考慮土體自重的情況下基于Darcy或非Darcy滲流的孔壓消散速率都要慢于Terzaghi解.

對(duì)比圖9和10可以看出地基沉降速率和孔壓消散速率有著不同的情況.當(dāng)c值較小時(shí)（例如c＝0.5），基于Darcy滲流考慮土體自重和不考慮土體自重時(shí)的地基沉降速率要明顯快于Terzaghi解，基于非Darcy滲流時(shí)由于非Darcy特性的影響只在固結(jié)快要完成的時(shí)候出現(xiàn)了固結(jié)沉降速率慢于Terzaghi解的情況.當(dāng)c值較大時(shí)（例如c＝1.2），考慮土體自重應(yīng)力時(shí)的地基沉降速率和Terzaghi解相對(duì)比則可能出現(xiàn)前期快而后期慢的相反情況.

5 結(jié)論

本文在考慮土體變形非線性的基礎(chǔ)上引入Hansbo滲流模型，并同時(shí)考慮土體自重和埋深的影響，以有效應(yīng)力為求解對(duì)象，推導(dǎo)了飽和黏土的1維固結(jié)方程，并給出了隱式有限差分格式.數(shù)值分析結(jié)果表明了分析飽和黏土非線性固結(jié)時(shí)，考慮滲流非Darcy特性且同時(shí)考慮土體自重和埋深影響的必要性，并得出下列結(jié)論：

（1）無論是基于Darcy還是基于非Darcy滲流，考慮土體自重應(yīng)力時(shí)的孔壓消散速率和地基沉降速率都要大于不考慮土體自重時(shí)的情況.無論是考慮或不考慮土體自重應(yīng)力，土體滲流的非Darcy特性都延緩了孔壓消散速率和地基沉降速率.

（2）土體自重對(duì)孔壓消散的影響主要表現(xiàn)在固結(jié)的中期，而對(duì)前期和后期無明顯影響；對(duì)地基沉降的影響主要表現(xiàn)在前期和中期，而對(duì)后期的影響較小.土體滲流的非Darcy特性對(duì)地基沉降和孔壓消散的影響相似，主要表現(xiàn)為在中后期影響較大而對(duì)前期的影響較小.

（3）當(dāng)c值較小時(shí)（例如c＝0.5），在考慮土體自重的情況下，基于Darcy滲流的孔壓消散速率快于Terzaghi解；但基于非Darcy滲流的解答，則可能出現(xiàn)在固結(jié)的中期孔壓消散速率快于Terzaghi解，而在固結(jié)后期，由于非Darcy特性的影響其孔壓消散速率慢于Terzaghi解的情況.

（4）當(dāng)c值較大時(shí)（例如c＝1.2），在考慮土體自重的情況下，基于Darcy或非Darcy滲流孔壓消散速率都要慢于Terzaghi解.但考慮土體自重應(yīng)力時(shí)的地基沉降速率和Terzaghi解相對(duì)比則可能出現(xiàn)前期快而后期慢的相反情況.

［1］Davis E H，Raymand G P.A non－linear theory of consolidation［J］.Geotechnique，1965，15（2）：161.

［2］Xie K H，Xie X Y，Jiang W.A study on one－dimensional nonlinear of consolidation of bouble－layered soil ［J］.Computers and Geotechnics，2002，29：151.

［3］Chen R P，Zhou W H，Wang H Z，et al.One－dimensional nonlinear consolidation of multi－layered soil by differential quadrature method ［J］.Computers and Geotechnics，2005，32：358.

［4］Abbasi N，Rahimi H，Javadi A A，et al.Finite difference approach for consolidation with variable compressibility and permeability ［J］.Computers and Geotechnics，2007，34（1）：41.

［5］Gibson R E，Schiffman R L，Cargill K W.The theory of onedimensional of saturated clay.Ⅱ： finite non－linear consolidation of think homogeneous layers ［J］.Candian Geotechnical Journal，1981，18：280.

［6］丁洲祥，龔曉南，李又云，等.割線模量法在沉降計(jì)算中存在的問題及改進(jìn)探討［J］.巖土工程學(xué)報(bào)，2005，27（3）：313.DING Zhouxiang，GONG Xiaonan，LI Youyun.Study on some problems of settlement calculation with secant modulus method and its modification ［J］.Chinese Journal of Geotechnical Engineering，2005，27（3）：313.

［7］竇宜，蔡正銀.自重應(yīng)力作用下飽和黏土的固結(jié)變形特性［J］.巖土工程學(xué)報(bào)，1992，14（6）：29.DOU Yi， CAI Zhengyin， SHENG Shuxin.Consolidation Characteristics of Saturated Clays ［J］.Chinese Journal of Geotechnical Engineering，1992，14（6）：29.

［8］李冰河，謝康和，應(yīng)宏偉，等.初始有效應(yīng)力沿深度變化的非線性一維固結(jié)半解析解［J］.土木工程學(xué)報(bào)，1999，32（6）：47.LI Binghe，XIE Kanghe，YING Hongwei，et al.Semianalytical solution of 1－D nonlinear consolidation considering the initial effcetive stress distribution ［J］.China Civil Engineering Journal，1999，32（6）：47.

［9］馬崇武，劉忠玉.考慮飽和黏土埋深影響的一維非線性固結(jié)［J］.巖石力學(xué)與工程學(xué)報(bào)，2007，26（增2）：4372.MA Chong wu， LIU Zhongyu.One－dimensional nonlinear consolidation considering the buried depth of saturated clay layer［J］.Chinese Journal of Rock Mechanics and Engineering，2007，26（S2）：4372.

［10］Gibson R E，England G L，Hussey M J L.The theory of one dimensional consolidation of saturated clays.I：finite nonlinear consolidation of thin homogeneous layers ［J］.Geotechnique，1967，17：261.

［11］王秀艷，劉長禮.對(duì)黏性土孔隙水深流規(guī)律本質(zhì)的新認(rèn)識(shí)［J］.地球?qū)W報(bào)，2003，24（1）：91.WANG Xiuyan，LIU Changli.New understanding of the regularity of water seepage in cohesive soil［J］.Acta Geoscientica Sinica，2003，24（1）：91.

［12］齊添，謝康和，胡安峰，等.蕭山黏土非達(dá)西滲流性狀的試驗(yàn)研究［J］.浙江大學(xué)學(xué)報(bào)：自然科學(xué)版，2007，41（6）：1023.QI Tian，XIE Kanghe，HU Anfeng，et al.Laboratorial study on non－Darcy seepage in Xiaoshan clay［J］.Journal of Zhejiang University：Engineering Science，2007，41（6）：1 023.

［13］Deng Y E，Xie H P，Huang R Q，et al.Law of nonlinear flow in saturated clays and radial consolidation ［J］.Applied Mathematics and Mechanics，2007，28（11）：1427－1436.

［14］Hansbo S.Consolidation of clay with special reference to influence of vertical sand drains［J］.Swedish Geotechnical Institute，1960，18：45.

［15］Law K T，Lee C F.Initial gradient in a dense glacial till［C］／／Proc 10th Int Conf on Soil Mechanics and Foundation Engineering，Stockholm：［s.n.］，1981，441－446.

［16］劉忠玉，劉忠廣，馬崇武.考慮起始水力梯度時(shí)飽和黏土的一維固結(jié)［J］.鄭州大學(xué)學(xué)報(bào)：工學(xué)版，2006，27（3）：21.LIU Zhongyu， LIU Zhongguang， MA Chongwu.Onedimensional consolidation of saturated clays considering intial hydraulic gradient［J］.Journal of Zhengzhou University：Engineering Science，2006，27（3）：21.

［17］劉忠玉，張?zhí)旌?，馬崇武.起始水力梯度對(duì)飽和黏土一維固結(jié)的影響［J］.巖土力學(xué)，2007，29（3）：467.LIU Zhongyu，ZHANG Tianhang，MA Chongwu.Effect of initial hydraulic gradient on one－dimensional consolidation of saturated clays［J］.Rock and Soil Mechanics，2007，29（3）：467.

［18］謝海瀾，武 強(qiáng)，趙增敏，等.考慮非達(dá)西流的弱透水層固結(jié)計(jì)算［J］.巖土力學(xué)，2007，28（5）：1061.XIE Hailan，WU Qiang，ZHAO Zengmin，et al.Consolidation computation of aquitard considering non－Darcy flow［J］.Rock and Soil Mechanics，2007，28（5）：1061.

［19］劉忠玉，孫麗云，樂金朝，等.基于非Darcy滲流的飽和黏土一維固結(jié)理論［J］.巖石力學(xué)與工程學(xué)報(bào)，2009，28（5）：973.LIU Zhongyu，SUN Liyun， YUE Jinchao， et al.Onedimensional consolidation theory of saturated clay based on non－Darcy flow ［J］.Chinese Journal of Rock Mechanics and Engineering，2009，28（5）：973.

［20］鄂建，陳剛，孫愛榮.考慮低速非Darcy滲流的飽和黏性土一維固結(jié)分析［J］.巖土工程學(xué)報(bào)，2009，31（7）：1115.E Jian， CHEN Gang， SUN Airong.One－dimensional consolidation of saturated cohesive soil considering non－Darcy flows［J］.Chinese Journal of Geotechnical Engineering，2009，31（7）：1115.

［21］劉忠玉，糾永志，樂金朝，等.基于非Darcy滲流的飽和黏土一維非線性固結(jié)分析［J］.巖石力學(xué)與工程學(xué)報(bào)，2010，29（11）：2348.LIU Zhongyu，JIU Yongzhi，YUE Jinchao，et al.Onedimensional nonlinear consolidation of saturated clay based on non－Darcy flow［J］.Chinese Journal of Rock Mechanics and Engineering，2010，29（11）：2348.