黃雄立，韋 慧

(1.長沙理工大學交通運輸工程學院，湖南 長沙 410004;2.湖南理大交通科技發(fā)展有限公司，湖南 長沙 410015)

針對軟土地基上修建建筑物不能滿足工程要求的穩(wěn)定與變形問題，盡管目前有很多種地基處理方法，但排水固結(jié)法仍是實際工程中常用的有效方法。然而，在豎向排水體施工過程中，施工導管擠土作用會對地基土產(chǎn)生較大的擾動。砂井施工擾動很早就有學者認識到，并指出完全排土式砂井地基沉降穩(wěn)定時間甚至要比不處理的地基要長［1］。工程實踐結(jié)果同樣表明:豎井打設(shè)使土體的結(jié)構(gòu)發(fā)生擾動破壞，導致土體強度降低，地基承載力降低［2－4］;增加了建筑物的沉降［5－6］;影響豎井的排水效果。但上述更多的是集中在工程實錄階段，針對擾動機理的研究較少。因此，有必要對豎井施工擾動進行分析和評價。目前對豎井施工擾動效應的分析，廣泛采用Vesic［7］圓柱形孔擴張彈塑性解理論，分析中大多將土體視為理想彈塑性材料，服從Tresca屈服準則或Mohr－Coulomb 屈服準則［8－10］。能夠較好地反映軟黏土應力應變關(guān)系的本構(gòu)模型應為基于臨界土力學概念的修正劍橋模型。Cao等在修正劍橋模型基礎(chǔ)上給出了不排水擴張過程中土體應力的精確解析解，并分析了旁壓試驗、靜力觸探等土工原位測試問題［11－12］。但該精確解需要通過數(shù)值積分實現(xiàn)，實際應用不方便。本文根據(jù)文獻［12］精確解的一種簡化形式，對豎井施工過程中土體應力狀態(tài)進行分析，并在此基礎(chǔ)上評價了豎井施工擾動效應，建立了一個修正的擾動函數(shù)。

1 豎井周圍土體應力變化

1.1 計算模型及基本方程

豎井的打設(shè)一般在短時間內(nèi)就可完成，將此過程視為不排水條件下平面應變軸對稱的圓柱形孔擴張問題，計算模型參見圖1。假設(shè)飽和軟黏土中圓孔的初始應力場、孔徑為σ0和a0。當均勻分布的內(nèi)壓力從σ0增大到σa時，孔徑從a0增大到a。隨著內(nèi)壓力值的增加，圍繞著圓孔的圓柱形區(qū)域?qū)⒂蓮椥誀顟B(tài)進入塑性狀態(tài)，塑性區(qū)隨σa的增加而不斷擴大，rp為塑性區(qū)半徑。分析井周應力狀態(tài)時，彈性階段本構(gòu)方程為廣義虎克定律;塑性區(qū)的屈服準則采用修正劍橋模型進行描述，表達式為:

圖1 豎井的圓孔擴張分析模型Fig.1 Analytical model of cylindrical expansion for vertical drain

式中:q和p'分別為偏應力和有效平均應力;p'y為等向固結(jié)屈服應力;M為臨界狀態(tài)線的斜率。

由平面應變軸對稱問題的平衡微分方程、幾何方程和上述的屈服條件，圓柱形孔擴張問題塑性區(qū)的有關(guān)解答如下［12］。

有效應力q、p'與徑向距離r的關(guān)系為:

偏應力比q與平均有效應力p'的關(guān)系:

式中:p'0為初始平均有效應力，p'0=σ0。

因此，聯(lián)立式(2)和(3)，可得到q與r的之間關(guān)系，再根據(jù)應力平衡方程、邊界條件和有效應力原理，可到孔壓u與 r的關(guān)系［12］:

然而，式(4)的積分不能顯式地表達，需通過數(shù)值積分計算，在實際應用不夠方便。本文根據(jù)塑性區(qū)的偏應力與極限偏應力相差不大的事實［12］，假定塑性區(qū)的偏應力q等于極限偏應力qu，可得與上述相關(guān)的近似解答。

塑性區(qū)半徑與孔徑的關(guān)系為:

式中:σr和σθ分別為徑向應力和環(huán)向應力。

1.2 孔周應力分布規(guī)律

修正劍橋模型相對于Tresca屈服準則和Mohr－Coulomb屈服準則可以考慮土體屈服應力(應力歷史)的影響。根據(jù)等向平均有效應力p'和豎向有效應力σ'v的關(guān)系，可得到等向超固結(jié)比R與傳統(tǒng)的豎向超固結(jié)比ROC的相互關(guān)系為［11］:

假設(shè) G/p(0=30，M=1.2，Λ =0.75，取 ROC分別為1.0，1.5 和 3.0，根據(jù)公式(6)，(7)和(8)，分析了ROC對圓柱形孔周應力分布的影響。選取豎井施工對井周應力狀態(tài)影響的特征量——超靜孔壓u和徑向應力σr，其歸一化的應力與徑向距離的關(guān)系參見圖2。由圖可見:ROC越大，在孔周產(chǎn)生的u和σr越大，但衰減的速度越快，因此塑性區(qū)半徑rp和影響范圍越小;不同ROC的u和σr均隨徑向距離增加而呈對數(shù)衰減。

另外，對于天然沉積軟黏土，ROC越大，表明屈服應力越大，土結(jié)構(gòu)性越強，靈敏度就越大。因此，從另一方面表明靈敏度同樣對孔周應力的分布有影響。關(guān)于靈敏度對擠土產(chǎn)生的超靜孔壓的影響，研究表明:沉樁過程中產(chǎn)生超靜孔隙水壓力與土的靈敏度有關(guān)，靈敏度越高的土中產(chǎn)生的孔壓也越高［13］。

2 施工擾動效應評價

2.1 擾動評價及修正擾動函數(shù)的建立

圖2 超固結(jié)比(ROC)對圓柱形孔周應力分布的影響Fig.2 Effect of ROCon stress around the cylindrical cavity

豎井周圍土體不排水強度因受到不同大小的施工擾動而不同程度的降低。文獻［10］根據(jù)孔周應力在徑向上呈對數(shù)衰減規(guī)律，假定擾動度在徑向上亦為對數(shù)衰減。并認為在塑性區(qū)內(nèi)，小孔壁即rm處擾動最大即完全擾動，此時飽和軟黏土的不排水強度為重塑土的不排水強度;在彈塑性邊界上，即rp處擾動度為0，不排水強度為原狀土的強度Su;井周土體強度Sd與擾動度FD的關(guān)系為:

根據(jù)上式，結(jié)合靈敏度定義與上述假定，文獻［10］提出一種基于飽和軟黏土靈敏度St的擾動度函數(shù)，表達式為:

由前面的分析可知，擾動引起的超靜孔壓u在徑向上為對數(shù)衰減，且受靈敏度的影響。因此，本文的擾動函數(shù)也采用基于靈敏度St的對數(shù)形式。但文獻［10］定義的擾動函取值范圍不是0～1，而是0～1－1/St，即在完全擾動的孔壁rm處擾動度不是其實際定義的100%，而是隨靈敏度St變化而變化的，物理意義不夠明確;另外，僅把塑性區(qū)當作豎井施工擾動區(qū)，只考慮塑性區(qū)的擾動，沒有考慮彈性區(qū)的擾動，但實際中也有部分彈性區(qū)存在擾動［5，14，15］。鑒于此，對于靈敏性土，為了更加準確地描述豎井施工產(chǎn)生的附加擾動度，本文對文獻［10］提出的擾動度函數(shù)進行修正，令其在擾動區(qū)的取值從1變化到0，表達式為:

式中:FD(Field disturbance)表示豎井施工造成的原位擾動;rp為塑性區(qū)半徑，可由式(5)確定;rd為擾動區(qū)半徑，根據(jù)此處擾動度FD=0的定義，得其表達式為:

從式(12)可看出，擾動區(qū)半徑rd隨St的增大而減小:這與圖2顯示的規(guī)律是一致的。進而對擾動區(qū)的不排水強度重新定義:

根據(jù)擾動度FD是徑向距離r對數(shù)函數(shù)的假定，修正的擾動函數(shù)進一步簡化為:

圖3 施工擾動度FD隨徑向距離的變化Fig.3 Variation of FDwith radical distance

2.2 擾動函數(shù)的驗證

豎井施工產(chǎn)生的附加擾動會降低井周土的水平向滲透系數(shù)。傳統(tǒng)固結(jié)理論認為井周附近由于施工擾動存在一個涂抹區(qū)，并假定涂抹區(qū)內(nèi)滲透系數(shù)是不變的，且涂抹區(qū)半徑是2～4倍豎井等效半徑。但最近有研究者根據(jù)室內(nèi)試驗和現(xiàn)場試驗結(jié)果表明:施工擾動區(qū)不僅僅是涂抹區(qū)，且滲透系數(shù)在擾動區(qū)內(nèi)不是不變的，而隨擾動度的大小變化［16－18］。擾動區(qū)的滲透系數(shù)與擾動度FD的關(guān)系定義為:

式中:kh為擾動區(qū)水平滲透系數(shù);khr為完全重塑土水平滲透系數(shù);khu為未擾動土水平滲透系數(shù)。

將式(14(a))代入式(15)，經(jīng)化簡后，得到歸一化后的滲透系數(shù)與徑向距離的關(guān)系為:

因此，根據(jù)式(16)計算的滲透系數(shù)結(jié)果，與本文收集到的一些文獻中的試驗結(jié)果進行比較，可以驗證提出的擾動函數(shù)合理性。但相關(guān)文獻中沒有提供計算系數(shù)A和B所需的完整參數(shù)，本文只能對文獻中的試驗數(shù)據(jù)進行擬合，定性地檢驗擾動函數(shù)合理性。擬合結(jié)果見表1和圖4。從表1和圖4可以看出:不管是室內(nèi)試驗還是現(xiàn)場試驗的結(jié)果，對數(shù)擬合系數(shù)均較高，與式 (16)所表示的趨勢是一致的，表明本文所定義的擾動函數(shù)是合理、有效的。

圖5所示為根據(jù)豎井(砂井和PVD)地基現(xiàn)場實測沉降反演得到的水平向固結(jié)系數(shù)Ch(e)經(jīng)室內(nèi)測得的原狀樣豎直向固結(jié)系數(shù)Cv歸一化后，與井徑比(dm為打設(shè)豎井套管按周長等效的直徑)的關(guān)系。從圖5可見:歸一化的固結(jié)系數(shù)對數(shù)值與徑向距離呈線性關(guān)系，即說明砂井地基固結(jié)系數(shù)隨距井中心的徑向距離增大而對數(shù)式增大。盡管圖5所示為根據(jù)實測沉降反演的結(jié)果，分析的是豎井地基平均擾動度與徑向距離的關(guān)系，但固結(jié)系數(shù)是隨擾動度而變化的，且兩者存在對數(shù)關(guān)系。因此，這也從另一方面說明了本文所給出的在徑向上的對數(shù)擾動函數(shù)存在合理性。

表1 試驗數(shù)據(jù)擬合結(jié)果Table 1 Fitted results of test data

圖4 歸一化的滲透系數(shù)與徑向距離的關(guān)系Fig.4 Relationship between normalized permeability coefficient and radical distance

圖5 歸一化固結(jié)系數(shù)與修正的井徑比關(guān)系［19］Fig.5 Relationship between Normalized consolidation coefficient and modified drain spacing ratio for vertical drains

3 結(jié)論

(1)采用修正劍橋模型作為屈服條件的圓柱孔擴張理論，得出豎井施工對井周應力狀態(tài)影響的特征量即超靜孔壓u和徑向應力σr在徑向上呈對數(shù)衰減，且受地基土靈敏度的影響。

(2)在準確分析豎井地基擾動區(qū)范圍和界定擾動函數(shù)取值范圍的基礎(chǔ)上，提出了修正的基于靈敏度考慮的對數(shù)擾動函數(shù)來評價豎井施工擾動影響。

(3)根據(jù)文獻中豎井打設(shè)的室內(nèi)模型試驗與實際工程現(xiàn)場試驗測得的固結(jié)指標隨徑向距離的變化結(jié)果，驗證了所提出的擾動函數(shù)合理性，故式(16)為準確分析豎井地基的固結(jié)奠定了很好的基礎(chǔ)。

［1］Casagrande L，Poulos S.On the effectiveness of sand drains［J］.Canadian Geotechnical Journal，1969，6(3):286－326.

［2］Akagi T.Effect of mandrel driven sand drains on strength［C］//Proceeding of the Tenth International Conference on Soil Mechanics and Foundation Engineering.1981.

［3］施建勇，趙維炳，艾英缽，等.砂井施工對軟黏土擾動的研究［J］.河海大學學報，1997 ，25(2):30－33.SHI Jian-yong，ZHAO Wei-bing， AI Ying-bo，et al.Effect of sand drainage construction on soft clay disturbance［J］.Journal of Hohai University:Natural Sciences，1997，25(2):30 －33.

［4］鄒新軍，楊 眉，趙明華.基于室內(nèi)模型試驗的砂井復合地基作用機理［J］.鐵道科學與工程學報，2009，6(3):61－64.ZOU Xin-jun，YANG Mei，ZHAO Ming-hua.Behavior of the sand well composite foundation by indoormodel test［J］.Journal of Railway Science and Engineering，2009，6(3):61－64.

［5］劉吉福，楊春林.珠江三角洲地區(qū)豎向排水體施工擾動初探［J］.巖石力學與工程學報，2003，22(1):142－147.LIU Ji-fu，YANG Chun-lin.Preliminary study onconstruction disturbance of vertical drainage bodies in the pearl river delta［J］.Chinese Journal of Rock Mechanics and Engineering，2003，22(1):142 －147.

［6］王立忠，李玲玲.結(jié)構(gòu)性土體的施工擾動及其對沉降的影響［J］.巖土工程學報，2007，29(5):697 －704.WANG Li-zhong，LI Ling-ling.Field disturbance of structured clay and its effect on settlements of soil foundation［J］.Chinese Journal of Geotechnical Engineering，2007，29(5):697 －704.

［7］Vesic A S.Expansion of cavities in an infinite soil mass［J］.Journal of Soil Mechanics and Foundation ，ASCE，1972 ，98(3):265－291.

［8］徐永福，傅德明.結(jié)構(gòu)性軟土中打樁引起的超孔隙水壓力［J］.巖土力學，2000，21(1):53－55.XU Yong-fu，F(xiàn)U De-ming.Excess pore pressure induced in piling in saturated structural soft soils［J］.Rock and Soil Mechanics，2000，21(1):53 －55.

［9］張儀萍，王 偉.豎井施工對軟土地基的擾動效應研究［J］.浙江大學學報:工學版，2006，40(7):1239－1243.ZHANG Yi-ping， WANG Wei. Study on disturbance effect on soft soil ground during vertical drain construction［J］.Journal of Zhejiang University:Engineering Science，2006，40(7):1239 －1243.

［10］徐書平，劉祖德，司馬軍.考慮擾動的飽和軟黏土球形孔擴張彈塑性解［J］.巖石力學與工程學報，2006，25(1):3212－3217.XU Shu-ping，LIU Zu-de，SI Ma-jun.Elato－plastic solution of spherical cavity expansion considering disturbed saturated soft clay［J］.Chinese Journal of Rock Mechanics and Engineering，2006，25(1):3212 －3217.

［11］Cao L F，Chang M F，Teh C I.Cavity expansion in modified Cam Clay and its application to the interpretation of piezocone tests［R］.Singapore:Geotechnical Research ReportNTU/GT96－03，Nanyang Technological University，1996.

［12］Cao L F，Teh C I，Chang M F.Undrained cavity expan-sion in modified cam clay.I:Theoretical analysis［J］.Geotechnique，2001，51(4):323 －334.

［13］Poulos H G，Davis E H.Pile foundation analysis and design［M］.London:John Wiley and Sons，1980.

［14］王 軍，高玉峰，高紅珍.結(jié)構(gòu)性軟土地基施工擾動定量分析［J］.巖土力學，2005，26(5):789－794.WANG Jun，GAO Yu-feng，GAO Hong-zhen.Quantitative analysis of foundation treatment disturbance of structured soft clay ground［J］.Rock and Soil Mechanics，2005，26(5):789 －794.

［15］雷金山，楊秀竹，夏力農(nóng)，等.軟土地層花管壓密注漿的有限元模擬［J］.鐵道科學與工程學報，2009，6(3):28－30.LEI Jin-shan，YANG Xiu-zhu，XIA Li-nong，et al.Finite element smiulation on compaction grouting with bored pipe in soft soil［J］.Journal of Railway Science and Engineering，2009，6(3):28 －30.

［16］Onoue A，Ting N H，Germaine J T，et al.Permeability of disturbed zone around vertical drains［C］//Proceedings of the Congress of the Geotechnical Engineering Division，American Society of Civil Engineers.New York，1991:879－890.

［17］Madhav M R，Park Y M，Miura N.Modeling and study of smear zones around band shaped drains［J］.Soils and Foundations，1993，33(4):135 －147.

［18］Sathananthan I，Indraratna B.Laboratory evaluation of smear zone and correlation between permeability and moisture content［J］.Journal of Geotechnical and Geoenvironmental Engineering，2006，132(7):942 －945.

［19］Saye S R.Assessment of soil disturbance by the installation of displacement sand drains and prefabricated vertical drains［J］.Geotechnical Special Publication，2003:372－393.