陳 余，李傳勛

(江蘇大學(xué)土木工程與力學(xué)學(xué)院，江蘇，鎮(zhèn)江 212013)

經(jīng)典太沙基固結(jié)計(jì)算時(shí)，常將排水邊界看作理想化的完全透水邊界，事實(shí)上，當(dāng)砂墊層透水性不良時(shí)，土層的邊界條件大多數(shù)介于透水與不透水之間。同時(shí)，太沙基一維固結(jié)理論中也存在初始條件與邊界條件相矛盾的邏輯問題，太沙基固結(jié)理論中的排水邊界在初始的瞬間超靜孔壓有突變，忽略了固結(jié)過程本身是一個(gè)連續(xù)的過程?；诖耍穱鄣萚1]提出了介于透水和不透水之間的連續(xù)排水邊界條件，對(duì)太沙基一維固結(jié)理論進(jìn)行了有益的補(bǔ)充。

與完全透水邊界、完全不透水邊界及第三類排水邊界相比較，連續(xù)排水邊界具有適定性，適用于大多數(shù)軟黏土地基的邊界排水情況，能嚴(yán)格滿足固結(jié)模型的初始條件，且其在一定條件下能夠退化為完全透水邊界，因此進(jìn)行基于連續(xù)排水邊界下土體固結(jié)理論的研究具有一定的理論意義和實(shí)際意義。張國英等[2]建立了連續(xù)排水邊界下一維均質(zhì)地基固結(jié)模型，采用有限傅里葉變換方法求得超靜孔壓解答，分析了不同排水邊界界面參數(shù)下固結(jié)性狀的異同。蔡烽等[3]基于ABAQUS有限元軟件開發(fā)了連續(xù)排水邊界條件下固結(jié)計(jì)算的子程序，并將計(jì)算結(jié)果與連續(xù)排水邊界下固結(jié)解析解對(duì)比驗(yàn)證子程序的可靠性，從而在技術(shù)上提供了連續(xù)排水邊界應(yīng)用于實(shí)際有限元計(jì)算的途徑。在此基礎(chǔ)上，馮健雪等[4]給出了單級(jí)加載下考慮連續(xù)排水邊界的土層一維固結(jié)解析解，其在求解過程中連續(xù)排水邊界是基于瞬時(shí)加載建立的。Feng等[5]建立了同時(shí)考慮自重與外荷載引起的初始孔壓非均布和連續(xù)排水邊界的地基一維固結(jié)模型，并給出了該模型的解析解。需要指出的是，以上基于連續(xù)排水邊界的固結(jié)理論均沒有考慮土體的非線性固結(jié)特性。

事實(shí)上土體的非線性壓縮和滲透特性對(duì)土體固結(jié)性狀的影響不可忽視。基于此宗夢(mèng)繁等[6]在初始有效應(yīng)力和固結(jié)系數(shù)保持不變的基本假定下，得到了連續(xù)排水邊界下土體一維非線性固結(jié)模型解析解。宗夢(mèng)繁等[7]進(jìn)一步考慮土體的流變固結(jié)特性，得到瞬時(shí)和單級(jí)加載下考慮連續(xù)排水邊界的土體一維流變固結(jié)解析解。田乙等[8]在連續(xù)排水邊界下采用分?jǐn)?shù)階Kelvin模型來描述飽和土體的流變特性，得到了分?jǐn)?shù)階導(dǎo)數(shù)黏彈性飽和土體一維固結(jié)半解析解。李稱等[9]在已有的簡化降解固結(jié)模型基礎(chǔ)上引入連續(xù)排水邊界條件，建立了不同排水邊界下飽和城市固廢一維降解固結(jié)模型及其解析解。江文豪等[10]采用大應(yīng)變幾何假定，研究了考慮連續(xù)排水邊界下飽和軟土的一維大變形固結(jié)問題。Zhang等[11]分析了連續(xù)排水邊界條件下豎井地基的超孔隙水壓力和平均固結(jié)度，探討了頂面排水能力、涂抹和井阻效應(yīng)對(duì)固結(jié)性狀的影響。Tian等[12]在多級(jí)荷載和連續(xù)排水邊界條件下，利用超孔隙水壓力疊加公式，得到了超孔隙水壓力和固結(jié)度的解析解，得到了任意時(shí)變荷載作用下更一般的連續(xù)排水邊界。由此可見，盡管目前對(duì)連續(xù)排水邊界下土體的線性和非線性固結(jié)問題均有一定研究，但前述的對(duì)連續(xù)排水邊界下固結(jié)理論的研究均假定土中滲流遵從達(dá)西定律，而實(shí)際工程中軟黏土中的滲流在低水力坡降下常常會(huì)出現(xiàn)偏離達(dá)西定律的現(xiàn)象。試驗(yàn)和實(shí)測(cè)均表明軟黏土中滲流在低水力坡降下可能會(huì)出現(xiàn)偏離達(dá)西定律的現(xiàn)象[13 ? 17]。其中Mitchell和Younger[14]的實(shí)驗(yàn)數(shù)據(jù)顯示，在低水力坡降下的飽和黏土中滲流的實(shí)測(cè)流速很小甚至沒有滲流發(fā)生，這是在水力坡降較小的情況下，由于低滲透性黏性土中有粘滯阻力的存在所導(dǎo)致的。為進(jìn)一步簡化飽和軟黏土中水的這一滲流現(xiàn)象，Miller和Low[18]提出了考慮初始水力坡降的滲流模型，即當(dāng)水力坡降i小于起始水力坡降i0時(shí)，土中無滲流發(fā)生；只有當(dāng)水力坡降i大于i0時(shí)，土中才會(huì)發(fā)生滲流。需說明的是，即使土中滲流在低水力坡降下沒有完全停止，但考慮起始水力坡降的滲流模型可看作是低水力坡降下流速極其微小的一種簡化模型，故有必要對(duì)考慮起始水力坡降的固結(jié)理論開展研究[19 ? 27]。目前考慮起始水力坡降的固結(jié)理論均認(rèn)為土層的排水邊界為完全透水或完全不透水，但實(shí)際中土層邊界的透水性并不是完全理想化的，其排水邊界孔壓隨時(shí)間指數(shù)衰減的連續(xù)排水邊界更符合實(shí)際工況，尤其對(duì)存在起始水力坡降的黏性土的排水邊界更適用。

本文同時(shí)考慮軟黏土中滲流存在起始水力坡降以及土層表面的連續(xù)排水邊界特性，對(duì)軟黏土地基一維固結(jié)理論展開研究，并著重分析邊界系數(shù)b和起始水力坡降i0對(duì)固結(jié)性狀的影響。本文提供的固結(jié)理論模型適用于大多數(shù)軟黏土地基的排水固結(jié)計(jì)算，尤其是黏粒含量較多導(dǎo)致的土中滲流在低水力坡降下存在起始值的工況，進(jìn)而為該類工程地基的固結(jié)變形計(jì)算提供重要的理論支持。

1 邊界條件的分析及固結(jié)模型的建立

如圖1所示，無限均布荷載q0施加于均質(zhì)黏土層，黏土層厚度為H，土層頂面是連續(xù)排水邊界，其界面排水系數(shù)為b，土層底面為完全不透水邊界，土體壓縮模量為Es。

圖1 軟土地基一維固結(jié)模型Fig.1 One-dimensional consolidation model of soft ground

本文在建立固結(jié)模型時(shí)采用了Terzaghi固結(jié)理論中除軟黏土中水滲流服從達(dá)西定律以外的其他基本假定。

土中滲流考慮存在的起始水力坡降，其值為i0，以重力方向向下為正，則土中滲流規(guī)律[18]可以表示為：

式中：i為水力坡降；kv為黏土層的滲透系數(shù)；v為黏土層中水的流速；i0為起始水力坡降，其值可以由黏性土的滲透試驗(yàn)來確定。滲透系數(shù)的室內(nèi)測(cè)定可采用常水頭法獲得，通過測(cè)試不同水頭下的滲透系數(shù)，繪制出滲流速度和水力坡降的關(guān)系曲線，曲線與水力梯度坐標(biāo)軸的交點(diǎn)可確定起始水力坡降。土層頂面連續(xù)排水邊界處的孔壓應(yīng)滿足[1]：

式中：u為超靜孔隙水壓力；t為時(shí)間；b為連續(xù)排水邊界界面系數(shù)，反映土體排水能力，界面參數(shù)是通過孔隙水壓力傳感器測(cè)量出的土體不同位置和時(shí)間的超靜孔隙水壓力再反演獲得[28]。b值越大邊界排水性越好，當(dāng)b→∞，連續(xù)排水邊界退化為完全透水邊界。土中滲流存在的起始坡降必然引起固結(jié)過程中的邊界移動(dòng)，記t時(shí)刻移動(dòng)邊界距離透水面的距離為h(t)，此時(shí)滲流鋒面z=h(t)處超靜孔隙水壓力u及水力坡降i應(yīng)該滿足的邊界條件為[19]：

式中：z為深度； γw為水的重度。

根據(jù)側(cè)限條件下孔隙比變化與豎向有效應(yīng)力變化間的關(guān)系，可建立滲流前鋒以上土體的一維固結(jié)控制方程：

式中，cv為固結(jié)系數(shù)，cv=kvEs/γw?？刂品匠虒?duì)應(yīng)的連續(xù)排水邊界為式(2)，起始水力坡降引起的移動(dòng)邊界表達(dá)式為式(3)和式(4)。模型的初始條件為：

2 固結(jié)模型的解析解答

2.1 滲流移動(dòng)邊界未到土層底面的解析解

至此，該問題對(duì)應(yīng)的固結(jié)模型已建立完備。

為方便求解做如下變量代換，令：

將式(7)代入式(3)～式(6)中，則固結(jié)控制微分方程、邊界條件及初始條件轉(zhuǎn)化為：

將w(z,t) 和q0e?bt按關(guān)于z的固有函數(shù)進(jìn)行傅里葉正弦級(jí)數(shù)展開，即：

進(jìn)而t時(shí)刻z≤h(t)的超靜孔壓表達(dá)式為：

當(dāng)h(t)

當(dāng)t趨近于無窮(固結(jié)已完成)時(shí)，如果滲流鋒面仍不能至土層底面，此時(shí)土中殘留孔壓沿深度分布為：

根據(jù)式(3)表達(dá)的移動(dòng)邊界處超靜孔壓需滿足的條件，結(jié)合式(14)超靜孔壓在z=h(t)處的表達(dá)式，可得到滲流移動(dòng)邊界隨時(shí)間變化規(guī)律為：

應(yīng)用無量綱變量，式(15)對(duì)應(yīng)的無量綱表達(dá)式為：

本研究初步提示，早期足量肝素用于STEMI患者，有助于改善冠脈血流、改善心功能，且并不增加出血事件。但，本研究為單中心研究，樣本數(shù)量有限，還需要以后更大樣本量的研究。

按t時(shí)刻孔壓定義的平均固結(jié)度Upt為：

將超靜孔壓表達(dá)式代入，平均固結(jié)度可進(jìn)一步表達(dá)為：

根據(jù)有效應(yīng)力原理，可將有效應(yīng)力表達(dá)為：

進(jìn)而可得到t時(shí)刻土體沉降量St為：

式中： σ′為有效應(yīng)力；Es為土體壓縮模量，在無側(cè)向變形條件下，豎向應(yīng)力和應(yīng)變的比值，反映土體對(duì)豎向壓縮變形的抵抗能力。令固結(jié)時(shí)間向無窮遠(yuǎn)處發(fā)展，根據(jù)殘留孔壓分布可得到移動(dòng)邊界能到達(dá)土層底面與移動(dòng)邊界始終不能到達(dá)土層底面的最終沉降量，其表達(dá)式為：

根據(jù)按變形定義的計(jì)算公式Ust=St/St→∞，最終得到按變形定義的平均固結(jié)度Ust。

2.2 滲流移動(dòng)邊界到達(dá)土層底面的解析解

當(dāng)滲流移動(dòng)邊界到達(dá)土層底面時(shí)，可在式(14)中令h(t)=H獲得t時(shí)刻超靜孔壓表達(dá)式為：

在滲流鋒面能發(fā)展至底面的條件下，當(dāng)t趨近于無窮(固結(jié)完成)時(shí)，此時(shí)殘留孔壓沿深度分布為：

根據(jù)平均固結(jié)度的定義：

得滲流移動(dòng)邊界已到達(dá)土層底部的平均固結(jié)度Upt：

根據(jù)不考慮側(cè)向變形條件下的壓縮量計(jì)算公式算得土層沉降量St為：

進(jìn)而得到按變形定義的平均固結(jié)度Ust為：

3 解析解的退化驗(yàn)證

3.1 基于連續(xù)邊界條件的固結(jié)模型解析解的退化

前文提到b是排水邊界的界面系數(shù)，其反映地基土體邊界的排水性能，其值越大，邊界的透水性越好；相反，其值越小，邊界的透水性越差。如果令b趨近于無窮大，那么連續(xù)排水邊界條件下的固結(jié)模型就退化為完全透水邊界下考慮起始比降的一維固結(jié)模型。滲流移動(dòng)邊界未到達(dá)土層底面的超靜孔壓計(jì)算式退化為：

當(dāng)0

當(dāng)h(t)

式(29)和式(30)正是王坤等[21]初始孔壓均布條件下考慮起始比降的一維固結(jié)解解析解。滲流移動(dòng)邊界X與時(shí)間因子Tv的關(guān)系式退化為：

滲流移動(dòng)邊界未到達(dá)土層底面時(shí)按孔壓定義的平均固結(jié)度和沉降量的計(jì)算式退化為：

滲流移動(dòng)邊界到達(dá)土層底面時(shí)按孔壓定義平均固結(jié)度和沉降量的計(jì)算式分別退化為：

3.2 不考慮起始水力坡降時(shí)解析解的退化

如果忽略考慮起始水力坡降的存在，考慮起始水力坡降的滲流定律就退化為達(dá)西定律，此時(shí)如果令i0=0(R=0)，超靜孔壓的解析表達(dá)式退化為：

式(36)恰好是張國英等[2]連續(xù)排水邊界條件下均質(zhì)地基單面排水一維固結(jié)解析解。按孔壓定義的平均固結(jié)度和沉降量的計(jì)算式分別退化為：

以上退化說明，現(xiàn)有的僅考慮起始水力坡降而忽略復(fù)雜排水邊界的固結(jié)解析解及達(dá)西定律下考慮連續(xù)排水邊界的固結(jié)解析解均是本文解析解的特例。

3.3 本文解析解與已有文獻(xiàn)的對(duì)比

由3.1小節(jié)知，當(dāng)B取值足夠大可近似看作完全透水邊界，再根據(jù)文獻(xiàn)[27]，將完全透水邊界下的超靜孔壓隨時(shí)間變化的關(guān)系曲線與B=100 000時(shí)本文模型的超靜孔壓與時(shí)間的關(guān)系曲線進(jìn)行對(duì)比，發(fā)現(xiàn)兩者曲線幾乎重合，如圖2所示，再次驗(yàn)證了本文固結(jié)解析解的正確性與可靠性。

圖2 已有文獻(xiàn)解與本文解析解的對(duì)比(B=100 000)Fig.2 Comparison between existing solution in literature and analytical solution (B=100 000)

4 固結(jié)性狀分析

影響本文固結(jié)模型性狀的主要因素是土中滲流存在的起始水力坡降和邊界透水性能兩方面。因此，本文將重點(diǎn)分析無量綱變量R及反映界面參數(shù)b的無量綱變量B對(duì)固結(jié)性狀的影響。

如圖3所示，某工程地基土體為均質(zhì)飽和軟土，其上堆載了5 m高的路基填土，路基填土對(duì)地基土層頂面瞬時(shí)施加了100 kPa的恒荷載，整個(gè)地基土層處于正常固結(jié)狀態(tài)(土層自重作用下的固結(jié)變形已完成)，地基土層底部為完全不透水層。土層固結(jié)計(jì)算的相關(guān)參數(shù)見表1。

圖3 黏性土層的邊界條件Fig.3 Boundary conditions of cohesive soil

表1 計(jì)算分析所用參數(shù)Table 1 Parameters adopted in the following analysis

4.1 基于R對(duì)固結(jié)性狀的影響

對(duì)同一地基土層而言，R值反映了外荷載和起始水力坡降的影響。故研究不同R值下地基土體的固結(jié)性狀也就反映了荷載和起始水力坡降對(duì)固結(jié)性狀的影響。

起始水力坡降i0的存在致使土中超靜孔壓不能完全消散，且當(dāng)i0取值大于q0/(γwH)(R>1)時(shí)，移動(dòng)邊界始終不能達(dá)到土層底面，在施加荷載后瞬間從排水面至滲流移動(dòng)邊界處的土層發(fā)生滲流。當(dāng)i0取值小于q0/(γwH) (R<1)時(shí)，移動(dòng)邊界以下至不排水面的土層不發(fā)生滲流，但隨固結(jié)時(shí)間的推移，滲流移動(dòng)邊界將最終發(fā)展至土層底面，整個(gè)土層均會(huì)發(fā)生滲流，即滲流鋒面最終能發(fā)展至土層底面。圖4分別給出當(dāng)R>1和R<1情況下的移動(dòng)邊界隨時(shí)間的發(fā)展過程。從圖中可發(fā)現(xiàn)，滲流移動(dòng)邊界向土層底部發(fā)展的快慢與R值的大小密切相關(guān)，R值越大，移動(dòng)邊界下移速度越慢，R值越小，滲流移動(dòng)邊界下移速度則越快。當(dāng)R值大于1時(shí)，滲流前鋒則不能到達(dá)土層底面，其最終停留位置為1/R。在連續(xù)排水邊界下，R值對(duì)固結(jié)性狀的影響與完全透水邊界下相比并未發(fā)生明顯改變，其固結(jié)性狀與文獻(xiàn)[24]描述的固結(jié)性狀基本一致。

圖4 R對(duì)移動(dòng)邊界的影響Fig.4 Influence of R on moving boundary

圖5為Tv=0.3時(shí)超靜孔隙水壓力隨深度變化曲線。可發(fā)現(xiàn)同一深度處達(dá)西定律下(R=0)土中超靜孔壓殘留值最小，超靜孔壓殘留值會(huì)隨R值增大而增大。這說明R值影響超靜孔壓的消散速率，其消散速率隨R值增大而逐漸減慢。同時(shí)，達(dá)西定律(R=0)下隨著固結(jié)時(shí)間的推移，超靜孔壓能夠在固結(jié)完成時(shí)消散完成，即超靜孔壓的最終殘留值為0。圖6為不同R值下z/H=0.5處超靜孔壓隨時(shí)間消散過程，可發(fā)現(xiàn)達(dá)西定律下超靜孔壓消散最快，且其能夠最終消散至0。R值一旦大于0，則超靜孔壓就不能最終消散至0。且R值越大，超靜孔壓消散速率越慢，同時(shí)刻下超靜孔壓殘留值越大。

圖5 R對(duì)z/H -u/q0曲線的影響(Tv=0.3)Fig.5 Influence of R on z/H -u/q0(Tv=0.3)

圖6 R對(duì)u/q0-Tv曲線的影響(z/H=0.5)Fig.6 Influence of R on u/q0-Tv(z/H=0.5)

以上超靜孔壓消散規(guī)律直接反應(yīng)至土層固結(jié)度表現(xiàn)為，當(dāng)R=0(達(dá)西定律)時(shí)，由于最終超靜孔壓能夠完全消散，此時(shí)平均固結(jié)度能到達(dá)100%。一旦R值不等于0，鑒于超靜孔壓最終不能完全消散，則土層平均固結(jié)度最終值就不能到達(dá)100%。同一時(shí)刻下平均固結(jié)度隨著R值的增大而減小。且R值越大，孔壓消散完成的時(shí)間會(huì)越短，因此固結(jié)度達(dá)到穩(wěn)定值的時(shí)間也越短，如圖7所示。

圖7 R對(duì)固結(jié)度Up的影響Fig.7 Influence of R on average consolidation degree Up

4.2 基于B對(duì)固結(jié)性狀的影響

界面系數(shù)b的大小反映了排水面邊界的透水性情況，通過對(duì)b值大小的控制，可以得到不同透水性邊界條件的固結(jié)模型。事實(shí)證明，排水邊界透水性的強(qiáng)弱對(duì)固結(jié)性狀有著顯著影響。

隨著界面排水能力的降低(B值變小)，土中超靜孔隙水壓力消散減慢，滲流移動(dòng)邊界下移速率亦會(huì)減慢，如圖8所示。B值越小，移動(dòng)邊界到達(dá)底部的時(shí)間越長；B值越大，移動(dòng)邊界到達(dá)土層底部的時(shí)間越少，當(dāng)B=100 000時(shí)，連續(xù)排水邊界可以近似看作完全透水邊界，此時(shí)移動(dòng)邊界隨時(shí)間的下移速度最快。

圖8 B對(duì)滲流移動(dòng)邊界的影響Fig.8 Influence of B on moving boundary

圖9為不同界面參數(shù)下Tv=0.8時(shí)超靜孔壓隨深度分布曲線圖。通過計(jì)算發(fā)現(xiàn)滲流移動(dòng)邊界在Tv=0.8時(shí)已移動(dòng)至土層底面。從分布曲線中可發(fā)現(xiàn)，同一深度處超靜孔壓值隨B值的減小而逐漸增大。說明該時(shí)刻土中超靜孔壓在完全透水邊界下消散最快，殘留的超靜孔壓值最小。隨著B值的減小，超靜孔壓消散速率逐漸減慢，土中殘留的超靜孔壓逐漸增大。

圖9 B對(duì)z/H-u/q0曲線的影響(Tv=0.8)Fig.9 Influence of B on z/H-u/q0(Tv=0.8)

同樣的，固結(jié)性狀在某深度處超靜孔壓隨時(shí)間消散曲線中也有所展現(xiàn)。圖10為在z/H=0.1處的超靜孔壓隨時(shí)間變化的曲線，該組曲線進(jìn)一步描述了無量綱變量B對(duì)超靜孔壓消散的影響。B值越大，某深度處的超靜孔壓消散速率越快，同時(shí)也越接近于完全透水邊界下的u/q0-Tv曲線。這樣的固結(jié)性狀在固結(jié)度隨時(shí)間發(fā)展曲線中進(jìn)一步展現(xiàn)，如圖11所示，B值越大，相同時(shí)間下地基固結(jié)度越高，且完全透水邊界下土層的平均固結(jié)度最大。

圖10 B對(duì)u/q0-Tv曲線的影響(z/H=0.1)Fig.10 Influence of B on u/q0(z/H=0.1)

圖11 B對(duì)固結(jié)度Up的影響Fig.11 Influence of B on average consolidation degree Up

5 結(jié)論

考慮軟黏土地基中存在的不同透水性邊界及土中滲流存在的起始水力坡降，在太沙基一維固結(jié)理論的基礎(chǔ)上重新建立并推導(dǎo)了一維固結(jié)控制方程，得到其解析解，結(jié)論如下：

(1)給出了基于連續(xù)排水邊界和起始水力坡降的軟土一維固結(jié)解析解。該解答為同時(shí)考慮起始水力坡降和不同透水性邊界的實(shí)際軟土固結(jié)計(jì)算提供了可供參考的計(jì)算方法。

(2)本文解可退化為考慮起始水力坡降的一維線性固結(jié)解析解，也可退化為達(dá)西定律下考慮連續(xù)排水邊界的軟土一維固結(jié)解析解。這兩種情況下的解析解均是本文解析解的特例。

(3)連續(xù)排水邊界下起始水力坡降對(duì)軟黏土固結(jié)性狀影響與完全透水邊界下相比，并未發(fā)生明顯改變。

(4)連續(xù)排水邊界對(duì)起始水力坡降所引起的移動(dòng)邊界下移速度影響較大，透水性越差，移動(dòng)邊界下移速度越慢。但總體而言，考慮起始水力坡降后連續(xù)排水邊界排水能力對(duì)超靜孔壓消散的影響與達(dá)西定律下相比沒有發(fā)生明顯改變。