郭彪，韓穎，龔曉南，盧萌盟

(1. 重慶市市政設(shè)計研究院，重慶，400020；2. 浙江大學(xué) 軟弱土與環(huán)境土工教育部重點實驗室，浙江 杭州，310058；3. 機械工業(yè)第三設(shè)計研究院，重慶，400039；4. 中國礦業(yè)大學(xué) 深部巖土力學(xué)與地下工程國家重點實驗室，江蘇 徐州，221008)

砂井聯(lián)合堆載預(yù)壓是提高軟土地基承載力、減小地基沉降的有效方法，對常荷載下砂井地基的固結(jié)問題，國內(nèi)外眾多學(xué)者進行了大量的研究，取得了豐富的成果[1-16]。但在實際工程中，地基荷載通常是隨時間變化的，對于變荷載下的砂井地基固結(jié)問題，也有不少研究者進行了研究，如：Olson[16]假設(shè)砂井為理想井，給出了荷載分級線性施加情況下的解析解；Tang等[17-18]進一步考慮了涂抹和井阻作用；Zhu等[19]用數(shù)值方法對砂井地基進行了研究，并考慮了大變形、非線性等因素的影響；耿雪玉等[20]運用 Laplace變換求得了復(fù)雜荷載作用下未打穿砂井地基的固結(jié)解。上述解答要么考慮的因素不夠全面，要么過于復(fù)雜，很難應(yīng)用于工程實踐。另外，大量的研究[12-14]表明：涂抹區(qū)土體由于受到的施工擾動程度不同，其滲透性是連續(xù)變化的，離砂井越近，擾動程度越高，土的滲透性越小。為此，本文作者考慮涂抹效應(yīng)、井阻作用、地基土橫豎向組合滲流等因素，并且涂抹區(qū)的水平滲透系數(shù)采用更符合實際的線性連續(xù)變化模式，將上部荷載通過傅里葉變換轉(zhuǎn)化為三角級數(shù)的形式[21]，得到更便于工程應(yīng)用的分析變荷載下砂井地基固結(jié)的解析解，并編制了應(yīng)用程序。

1 固結(jié)方程及求解

圖1所示為隨時間任意變化荷載下砂井地基示意圖。在推導(dǎo)過程中，進行以下假定：

① 等應(yīng)變條件成立；

② 土中水的滲流服從Darcy定律；

③ 荷載在地基中引起的附加應(yīng)力沿深度均勻分布；

④ 在任一深度 z處從土體流入樁體的水量等于樁體中向上水流的增量。

根據(jù)等應(yīng)變條件，有：

式中：εv為地基豎向應(yīng)變；Es為土體壓縮模量；q(t)為上部荷載；u為土體任一深度處的平均超靜孔壓，

圖1 隨時間任意變化荷載下砂井地基簡化模型Fig.1 Model for consolidation of soft soil with vertical drains and general time-dependent loading

rw為砂井半徑；re為砂井影響區(qū)半徑；u為土體任一點的超靜孔壓。

根據(jù)謝康和等[9]的研究，砂井地基的固結(jié)方程為：

式中：kv為土體豎向滲透系數(shù)；kw為砂井滲透系數(shù)；kmf(r)為土體水平向滲透系數(shù)；km為最大水平滲透系數(shù)；f(r)為描述滲透系數(shù)隨r變化的函數(shù)。樁周流量連續(xù)方程為：

式中：uw為砂井任一深度處的超靜孔壓。邊界條件為：

任何加載都需要一個過程，即使是瞬時加載也是如此，只不過時間很短。因此，初始條件可寫為：

⑤ t=0時，u=0。

[10，17]，可以得到固結(jié)控制方程為：

式中：

假設(shè)地基土體涂抹區(qū)水平滲透系數(shù)隨r增大而線性增大，未擾動區(qū)保持不變，不妨稱其為模式 1，如圖2(a)所示。參考文獻[10]，可以得到：

其中：

圖2 影響區(qū)土體水平滲透系數(shù)Fig.2 Horizontal permeability coefficient in influence zone

s=rs/rw；rs為涂抹區(qū)半徑；δ為最小滲透系數(shù)和最大滲透系數(shù)之比，0＜δ=kn/km＜1。

為了便于比較分析，本文也給出涂抹區(qū)水平滲透系數(shù)不變情況下(模式2)的解，如圖2(b)所示。此時，F(xiàn)a的表達式為：

方程(5)是非齊次微分方程，其邊界條件式③和④以及初始條件式⑤是齊次的，滿足沖量定理條件，可設(shè)：

由式(6)，υ和υ滿足關(guān)系式：

其中，wυ方程(5)對應(yīng)齊次微分方程的解：

邊界條件和初始條件為：

采用分離變量法求解，設(shè)

代入式(18)可得：

即

利用邊界條件⑥和⑦以及式(17)和(19)，可得：

常微分方程(21)的特征方程為：

其解可表示為：

式(21)的通解可表示為：

由邊界條件⑧和⑨可得：

則 Z ( z) = asin(λ z)。由z=H和Z′(z)=0可得：

將k=λi代入式(23)，可得：

方程(22)的通解可表示為：

由于其初始時刻為t=τ+0，可設(shè) T (τ+0)= Δ ，Δ 為1個常數(shù)，則：

方程(18)的解可寫成：

將式(29)代入式(17)可得：

根據(jù)三角函數(shù)的正交性和初始條件⑧可得：

將式(29)和(30)分別代入式(15)和(16)可得固結(jié)方程的解為：

在 q(t)比較復(fù)雜的情況下，上面的解難以應(yīng)用，為方便工程應(yīng)用，需要將其簡化。對荷載 q(t)，一般來說是滿足狄利克雷(Dirichlet)條件的，可將其展開成傅里葉級數(shù)：

式中：ωn= 2 n π/ T；T為q(t)的周期。在q0的作用下，dq0=dt=0，因此，對地基孔壓沒有影響。

在變載qn(t)作用下，將式(39)代入式(32)和(33)并積分可以得到：

運用疊加原理，地基在荷載 q(t)作用下，砂井和地基土的平均超靜孔壓可寫成：

求得了u( z, t)后，可以得到地基任意深處的固結(jié)度為：

式中：qu為最大荷載。

進一步可以得到地基的總平均估計度為：

式中：

2 an和bn的求解

系數(shù)an和bn的求得是上述求解過程中的關(guān)鍵。在分級線性施加荷載情況下(如圖 3所示)，對其中任意一級加載，q(t)可寫為：

將式(47)代入式(36)和(37)可以得到：

其中：

圖3 分級線性施加荷載示意圖Fig.3 Loading process consisting of a combination of linear and constant functions with time

在荷載隨時間曲線變化的情況下，可以用直線段來模擬，如圖4所示。當直線段足夠多時，就可以達到足夠工程應(yīng)用的精度要求。這樣，可以同樣采用式(48)～(51)求得 an和 bn。

圖4 荷載隨時間曲線變化示意圖Fig.4 Relationship between loading and time

3 程序的編制

利用計算機編程，上面的求解過程很容易實現(xiàn)。步驟如下。

(1) 輸入地基土參數(shù)和砂井參數(shù)，包括：砂井半徑rw，影響區(qū)半徑re，涂抹區(qū)半徑rs，土層厚度H，砂井滲透系數(shù)kw，土體豎向滲透系數(shù)kv，未擾動土水平滲透系數(shù)km，擾動程序參數(shù)δ，土體壓縮模量Es等。

(2) 輸入上部荷載參數(shù)，若q(t)由如圖3所示的一系列直線段組合而成，則輸入各拐點處的時間ti和qi；若荷載為如圖4所示的隨時間曲線變化，則根據(jù)需要取一系列t，再由荷載函數(shù)q(t)計算對應(yīng)的q，時間間隔越小，精度便越高；輸入q(t)的周期T，若q(t)不是周期函數(shù)，則可將最大分析時間取作T。

(3) 由式(26)求得 βm。

(4) 利用式(48)，(49)和(50)求得一組系數(shù) an和 bn。

(5) 利用式(40)和(41)求得在荷載qn(t)作用下砂井和地基土的超靜孔壓

(6) 遞增 n，重復(fù)步驟(4)和(5)，將其結(jié)果疊加求得在荷載 q(t)作用下砂井和地基土的超靜孔壓 uw(z,t)和u( z, t)。

(7) 利用式(45)求得地基的平均固結(jié)度。

4 解的驗證

設(shè)荷載單級線性施加，如圖5所示。圖6所示為本文解與唐曉武解[17]的比較結(jié)果。由于唐曉武解假設(shè)涂抹區(qū)水平滲透系數(shù)不變，為了使兩者條件一致，本例解中采用式(14)表達的Fa，即模式2。

圖5 單級線性施加荷載Fig.5 Single ramp loading

圖6 本文解與唐曉武解的比較Fig.6 Comparison of the present solution with that by Tang X W

由圖6可見：2孔壓曲線基本上重合，由此可知本文解是可靠的。計算參數(shù)為：wr=0.2 m，re=4 m，rs=2 m，H=20 m，kw=1×10-3cm/s，km=kv=1×10-6cm/s，δ=0.5，Es=20 MPa，達到最大荷載時間t1=30 d，最終荷載qu=20 kN。

5 幾種荷載形式下的固結(jié)性狀

5.1 多級線性施加荷載

荷載分級線性施加示意圖如圖7所示。其中，荷載分兩級線性施加，t1=30 d，t2=60 d，t3=100 d，q1=20 kN，qu=50 kN。其他計算參數(shù)與上例的相同。2種滲透系數(shù)模式下超靜孔壓曲線及固結(jié)曲線如圖8和圖9所示。從圖8和圖9可以看出：考慮滲透系數(shù)呈線性變化時固結(jié)要快。

圖7 荷載分級線性施加示意圖Fig.7 Multi-ramp loading

圖8 荷載兩級線性施加時地基的超靜孔壓曲線Fig.8 Excess pore-water pressure of sand drains foundation under double-ramp loading

5.2 三角形循環(huán)荷載

三角形循環(huán)荷載如圖 10所示，其中 T=60 d，qu=50 kPa，其他計算參數(shù)與前面的相同。地基超靜孔壓變化曲線如圖11所示。從圖11可以看出：超靜孔壓呈震蕩減小趨勢，震蕩幅度并不隨時間減小，但其震蕩中線逐漸趨于一定值；在循環(huán)荷載作用下，當上部荷載減小即卸荷時，地基土回彈體積膨脹，從而產(chǎn)生負超靜孔壓。滲透系數(shù)為模式1時的孔壓曲線在模式 2時孔壓曲線的下方，即孔壓為正時，模式 1的孔壓小于模式2的孔壓，但模式1的負孔壓要比模式2的大。

圖9 荷載兩級線性施加時地基的固結(jié)度曲線Fig.9 Consolidation rate of sand drains foundation under double-ramp loading

圖10 三角形循環(huán)荷載Fig.10 Triangular cyclic loading

圖 12所示為三角形循環(huán)荷載作用下地基有效應(yīng)力曲線。從圖12可以看出：有效應(yīng)力呈震蕩上升趨勢，同樣其震蕩幅度不隨時間減小，但其震蕩中線逐漸趨于一定值。模式1的有效應(yīng)力要比模式2的大。

5.3 梯形循環(huán)荷載

地基在如圖 13所示的梯形循環(huán)荷載作用下，T=180 d，t1=30 d，t2=90 d，t3=120 d，qu=50 kPa，地基超靜孔壓曲線和附加有效應(yīng)力曲線如圖14和圖15所示，可以看出其與三角形循環(huán)荷載作用時有相似的規(guī)律。

圖11 三角形循環(huán)荷載作用下地基超靜孔壓曲線Fig.11 Excess pore-water pressure of sand drains foundation under triangular cyclic loading

圖12 三角形循環(huán)荷載作用下地基有效應(yīng)力曲線Fig.12 Effective stress of sand drains foundation under triangular cyclic loading

圖13 梯形循環(huán)荷載Fig.13 Trapezium cyclic loading

5.4 拋物線形循環(huán)荷載

如圖16所示為二次拋物線形循環(huán)荷載，拋物線方程為：

圖14 梯形循環(huán)荷載下地基超靜孔壓曲線Fig.14 Excess pore-water pressure of sand drains foundation under trapezium cyclic loading

圖15 梯形循環(huán)荷載下地基有效應(yīng)力曲線Fig.15 Effective stress of sand drains foundation under trapezium cyclic loading

圖16 拋物線形循環(huán)荷載Fig.16 Parabolic cyclic loading

如t1=60 d，T=180 d，qu=50 kPa，在循環(huán)荷載作用下，考慮滲透系數(shù)呈線性變化。地基的超靜孔壓曲線和附加有效應(yīng)力曲線如圖17和圖18所示。其中，X為逼近拋物線的直線段段數(shù)。

從圖17和圖18可以看出：它們均具有很好的收斂性，當X為2時4時有較大差別，X為4時與10時差別不大，而X為10時與為100時的孔壓曲線和有效應(yīng)力曲線基本上重合，這說明采用直線段模擬荷載隨時間曲線變化是很有效的。

圖17 拋物線形循環(huán)荷載下地基超靜孔壓曲線Fig.17 Excess pore-water pressure of sand drains foundation under parabolic cyclic loading

圖18 拋物線形循環(huán)荷載下地基有效應(yīng)力曲線Fig.18 Effective stress of sand drains foundation under parabolic cyclic loading

6 結(jié)論

(1) 在荷載分級施加時，考慮涂抹區(qū)滲透系數(shù)線性變化時地基固結(jié)速度比假設(shè)其不變時的快。

(2) 在循環(huán)荷載作用下，砂井地基超靜孔壓呈震蕩減小趨勢，在卸荷時可能會產(chǎn)生負壓，震蕩幅度并不隨時間減小，但其震蕩中線逐漸趨于一定值；考慮涂抹區(qū)滲透系數(shù)線性變化時的正孔壓比假設(shè)其不變時的小，而負超靜孔壓更大。

(3) 在循環(huán)荷載作用下，地基有效應(yīng)力呈震蕩上升趨勢；同樣，其震蕩幅度不隨時間減小，震蕩中線逐漸趨于一定值?？紤]滲透系數(shù)線性變化時有效應(yīng)力比假設(shè)其不變時的大。

(4) 當荷載隨時間曲線變化時，采用多段直線段逼近荷載曲線是求解此問題的有效方法。

參考文獻：

[1] Barron R A. Consolidation of fine grained soil by drains wells[J].Transactions of ASCE, 1948, 113: 718-742.

[2] Richart, Jr F E. A review of the theories for sand drains[J]. Proc ASCE, 1957, 83(S3): 1-38.

[3] Yoshikuni H, Nakanodo H. Consolidation of soils by vertical drain wells with finite permeability[J]. Soils and Foundations,1974, 14(2): 35-46.

[4] Hansbo S, Jamiolkowski M, Kok L. Consolidation by vertical drains[J]. Geotechnique, 1981, 31(1): 45-46.

[5] 曾國熙，王鐵儒，顧堯章. 砂井地基的若干問題[J]. 巖土工程學(xué)報，1981, 3(3): 74-81.ZENG Guo-xi, WANG Tie-ru, GU Yao-zhang. Several problems of sand drains[J]. Chinese Journal of Geotechnical Engineering,1981, 3(3): 74-81.

[6] Hart E G, Kindner R L, Boyer W C. Analysis for partially penetrating sand drains[J]. Journal of Soil Mechanics and Foundation Division, ASCE, 1958, 84(S4): 1-15.

[7] 陳根媛. 多層地基的一維固結(jié)計算方法與砂井地基計算的改進建議[J ]. 水利水運科學(xué)研究, 1984, 5(2) :18-29.CHEN Gen-yuan. One-dimensional consolidation of layered systems and a suggestion for improving calculation of sand-drained subsoil[J]. Journal of Nanjing Hydraulic Research Institute, 1984, 5(2): 18-29.

[8] 謝康和. 砂井地基: 固結(jié)理論、數(shù)值分析與優(yōu)化設(shè)計[D]. 杭州: 浙江大學(xué)建筑工程學(xué)院, 1987: 1-62.XIE Kang-he. Sand drained ground: Analytical & numerical solutions of consolidation and optimal design[D]. Hangzhou:Zhejiang University. College of Civil Engineering and Architecture, 1987: 1-62.

[9] 謝康和, 周開茂. 未打穿豎向排水井地基固結(jié)理論[J]. 巖土工程學(xué)報, 2006, 28(6): 679-684.XIE Kang-he, ZHOU Kai-mao. Consolidation theory for soft soil with partially penetrated vertical drains[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(6): 679-684.

[10] 盧萌盟, 謝康和, 張玉國, 等. 考慮土體水平滲透系數(shù)變化的復(fù)合地基固結(jié)解[J]. 浙江大學(xué)學(xué)報: 工學(xué)版, 2008, 41(11):1996-2001.LU Meng-meng, XIE Kang-he, ZHANG Yu-guo, et al.Analytical solution for consolidation of composite foundation accounting for variation of soil horizontal permeability coefficient[J]. Journal of Zhejiang University: Engineering Science, 2008, 41(11): 1996-2001.

[11] Tang X W, Onitsuka K. Consolidation of double- layered ground with vertical drains[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2001, 25: 449-465.

[12] Bergado D T, Asakami H, Alfaro M C, et al. Smear effects of vertical drains on soft Bangkok clay[J]. Journal of Geotechnical Engineering, 1991, 117(10): 1509-1530.

[13] Indraratna B, Redana I W. Babaoratory determination of smear zone due to vertical drain installation[J]. Journal of Geotechnical and Geoenvironmental Engineering, 1998, 124(2): 180-184.

[14] Sharma J S, Xiao D. Characterization of a smear zone around vertical drains by large-scale laboratory tests[J]. Canadian Geotechnical Journal, 2000, 37(6): 1265-1271.

[15] XIE Kang-he, LU Meng-meng, LIU Gan-bin. Equal strain consolidation for stone columns reinforced foundation[J].International Journal for Numerical and Analytical Methods in Geomechanics, 2009, 33(15): 1721-1735.

[16] Olson R E. Consolidation under time dependent loading[J].Journal of the Geotechnical Engineering Division, 1977,103(GT1): 55-60.

[17] Tang X W, Onitsuka K. Consolidation by vertical drains under time-dependent loading. International[J]. Journal for Numerical and Analytical Methods in Geomechanics, 2000, 24(9):739-751.

[18] Leo C J. Equal strain consolidation by vertical drains[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2004,130(3): 316-327.

[19] Zhu G F, Yin J H. Finite element analysis of consolidation of soils with vertical drain[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2000, 24(4):337-366.

[20] 耿雪玉, 蔡袁強, 徐長節(jié). 復(fù)雜荷載作用下考慮下臥層三維滲流的未打穿豎井地基固結(jié)分析[J]. 巖土力學(xué), 2008, 29(6):1521-1529.GENG Xue-yu, CAI Yuan-qiang, XU Chang-jie. Soil of consolidation with partially penetrated vertical drains under complicated loadings considering three-dimensional seepage in soil beneath vertical drain[J]. Rock and Soil Mechanics, 2008,29(6): 1521-1529.

[21] Conte E, Troncone A. Radial consolidation with vertical drains and general time-dependent loading[J]. Canadian Geotechnical Journal, 2009, 46(1): 25-36.