基于雙對數壓縮模型的真空預壓非線性固結解

韓文君 劉松玉 章定文 劉維正

(1東南大學交通學院, 南京 210096)(2常州市軌道交通發展有限公司, 常州 213022)(3中南大學土木工程學院, 長沙 410075)

真空預壓法是一種優良的地基加固方法．工程設計中,常采用砂井固結理論研究真空預壓固結過程[1-2]．Schiffman等[3]給出了砂井非線性固結解答;Hansbo等[4]采用水力梯度和滲透速度的指數關系來反應土體固結過程中滲透系數變化過程;Indraratna等[5]引入孔隙比e與有效應力σ′、孔隙比e與水平滲透系數kh的關系式,得到了正常固結土和超固結土的非線性固結解答;周琦等[6]采用e-σ′及e-kh關系式,研究了砂井地基正常固結土固結性狀;董志良[7]建立了真空預壓加固地基固結方程;Indraratna等[8]假定砂井內超靜孔壓(負壓)沿豎向呈線性衰減分布,求得整個地基的平均超靜孔壓;彭劼等[9-11]提出了考慮真空荷載實際邊界條件的固結解答．但是,已有的解答都未能同時考慮真空荷載實際邊界和土體固結參數變化的特性．

本文在文獻[5-6,10]的基礎上,考慮真空荷載沿深度衰減的實際邊界條件,引入能對土體壓縮曲線進行線性化的雙對數坐標lg(1+e)-lgσ′及lg(1+e)-lgkh,建立了真空預壓砂井地基非線性固結近似解答．

1 真空預壓砂井固結計算模型

1.1 基本假定及求解條件

完全打穿砂井地基固結的計算簡圖見圖1．假定土體滿足等應變條件,土中水的滲流服從達西定律;任意深度z處,從土體中流入砂井的水量等于砂井中向上水流量的增量．涂抹區與未擾動區的土體除滲透系數不同外,其他性質相同．真空荷載為瞬時荷載,固結過程中為常量．圖1中,dw,rw分別為砂井直徑和半徑;ds,rs分別為涂抹區直徑和半徑;de,re分別為砂井影響區直徑和半徑;r為距砂井中心距離;H為土層厚度;p0為膜下真空荷載的絕對值;kl為真空荷載衰減系數．

圖1 砂井地基徑向固結軸對稱問題的計算簡圖

根據文獻[12],砂井地基僅考慮徑向固結時的基本控制方程為

(1)

(2)

(3)

井周的流量連續方程為

(4)

式中,uw為排水板內任一深度的超靜孔壓;kw為排水板滲透系數．

式(1)～(4)的邊界條件如下:① 當r=re時,?ur/?r=0;② 當r=rs時,ks?us/?r=kh?ur/?r;③ 當r=rs時,ur=us;④ 當r=rw時,us=uw;⑤ 當z=0時,uw=-p0;⑥ 當z=H時,?uw/?z=p0(1-kl)/H．

1.2 固結方程的推導

整個地基的平均超靜孔隙水壓力值為

(5)

將式(1)和(2)兩邊均對r積分,結合邊界條件,將積分結果代入式(5)可得[8]

(6)

整個土層平均固結度為

(7)

式中,Ru為無量綱的中間變量．

將式(6)代入式(7),整理后可得[5]

(8)

式中,σ為土體總應力．

(9)

式中,Th為時間因子;ch0為土體的固結系數;e0為土體的初始孔隙比;mv0為土體的初始體積壓縮系數;kh0為土體的初始水平滲透系數．

引入反應土體固結非線性的雙對數模型,以便在同一坐標系下建立固結方程,即[13-14]

(10)

(11)

對于正常固結土,將式(10)對σ′求導,可得[6]

(12)

將式(12)聯合式(10)、(11),代入式(9),可得

(13)

(14)

式(14)是含有變量Ru的非線性偏微分方程,該方程沒有通解．由于0≤Ru≤1,則

(15)

式中，Wav為W的極限平均值.

(16)

古箏合奏曲在歷史中有很多，但大多是和其他的一些民樂合奏。現在古箏合奏的形式多樣，有和民樂合奏的《枉凝眉》還有和鋼琴合奏的《臨安遺恨》等等，在演奏形式上取得了重大突破。由于其豐富的音色和演奏技巧，所以有了很多不同的演奏形式。在獨奏和合奏方面不僅繼承了傳統還取得了新的成就，彈唱因為難度希望在將來有進一步的發展。

(17)

(18)

(19)

(20)

按沉降變形計算地基的平均固結度Us為

(21)

根據文獻[6,13-14]可知,t時刻土體沉降St及最終沉降S∞的計算公式如下:

正常固結土

(22)

(23)

(24)

(25)

(26)

(27)

2 實例驗證

表1 不同方法計算得到的最終沉降量 mm

由表1可知,利用不同方法計算的最終沉降量基本相近,且本文解答的計算結果與實測值最為接近,與Indraratna解[5]、Tuan解[10]及Hansbo解[12]的差別在5%以內;Tuan解[10]相對于Hansbo解[12]的沉降量偏小．

利用不同方法計算得到的超靜孔隙水壓力以及任意時刻的沉降量分別見圖2和圖3．由圖2可知,固結初期基于不同解答所得的超靜孔壓與實測數據都有不同程度的偏差,整體來說都是偏小的．隨著時間的增長,考慮固結過程土體參數變化的解答與超靜孔隙水壓力實測曲線吻合,基于Indraratna解[5]、周琦解[6]所得的結果略大于實測值．由圖3可知,不考慮非線性時,基于Tuan解[10]及Hansbo解[12]計算得到的沉降明顯小于實測值及其他計算值,尤其是在固結初期,這種差別較為顯著．周琦解[6]由于僅考慮了正常固結狀態時的沉降,因此其沉降計算值明顯小于實測值．基于Indraratna解[5]及本文解答所得的結果與實測值較為吻合．由此可見,基于本文解答的計算結果更加接近實測結果．

圖2 超靜孔隙水壓力

圖3 不同時刻的沉降量分布

3 固結性狀影響因素分析

圖5 /對固結度的影響

4 結語

本文考慮真空荷載沿深度衰減分布的實際邊界條件,引入能對土體壓縮曲線線性化的雙對數坐標lg(1+e)-lgσ′,建立了真空預壓砂井地基非線性固結近似解答,計算得到的沉降、超靜孔壓與實測值很接近．通過參數計算分析可知,對于超固結土而言,根據變形計算的固結度隨著壓縮指數與滲透指數比值的增大而減小,隨著真空度沿深度衰減系數的減小而減小;當壓縮指數與滲透指數比值大于1時,正常固結狀態中基于孔壓計算的固結度隨真空荷載衰減系數的減小而增大．隨著壓縮指數與回彈指數比值的增大,地基超靜孔壓消散速率逐漸增大,超固結狀態下的沉降速率逐漸減小,正常固結狀態下的沉降速率逐漸增大,土體由超固結狀態進入正常固結狀態的時間逐漸增加．

)

[1] 龔曉南. 地基處理手冊[M]. 3版.北京:中國建筑工業出版社,2008.

[2] 趙維炳. 控制工后沉降處理深厚軟土地基[M]. 北京:人民交通出版社,2006.

[3] Schiffman R L, Fungaroli A A. Consolidation due to tangential loads[C]//Proceedingsofthe6thInternationalConferenceonSoilMechanicsandFoundationEngineering. Montreal, Canada, 1965: 188-192.

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

[5] Indraratna B, Rujikiatkamjorn C, Sathananthan L. Radial consolidation of clay using compressibility indices and varying horizontal permeability[J].CanadianGeotechnicalJournal, 2005,42(1): 330-341.

[6] 周琦,劉漢龍,陳志波. 考慮固結參數變化時砂井地基的非線性徑向固結[J].巖土力學,2007,28(S1):855-858.

Zhou Qi, Liu Hanlong, Chen Zhibo. Non-linear radial consolidation of sand-drained ground with varying volume compressibility and horizontal permeability[J].RockandSoilMechanics, 2007,28(S1): 855-858. (in Chinese)

[7] 董志良. 堆載及真空預壓砂井地基固結解析理論[J]. 水運工程,1992(9):1-7.

Dong Zhiliang. Consolidation theory on heaped load and vacuum preloading of sand drain foundation[J].PortandWaterEngineering, 1992 (9): 1-7. (in Chinese)

[8] Indraratna B, Rujikiatkamjorn C, Sathananthan L. Analytical and numerical solutions for a single vertical drain including the effects of vacuum preloading[J].CanadianGeotechnicalJournal, 2005,42(4):994-1014.

[9] 彭劼,董江平,宋恩潤,等. 考慮加載過程的真空預壓軸對稱解析解[J]. 巖土力學,2010,31(S1):79-85.

Peng Jie, Dong Jiangping, Song Enrun,et al. Axisymmetric analytical solution of vacuum preload considering the loading process[J].RockandSoilMechanics, 2010,31(S1):79-85. (in Chinese)

[10] Tuan A T, Toshiyuki M. Equivalent plane strain modeling of vertical drains in soft ground under embankment combined with vacuum preloading[J].ComputersandGeotechnics, 2008,35(5): 655-672.

[11] 張儀萍,嚴露,俞亞南,等. 真空預壓加固軟土地基變形與固結計算研究[J]. 巖土力學,2011,32(S1):150-154.

Zhang Yiping, Yan Lu, Yu Yanan,et al. Solutions for settlement and consolidation of soft ground with vacuum preloading[J].RockandSoilMechanics, 2011,32(S1):150-154. (in Chinese)

[12] Hansbo S. Consolidation of fine-grained soils by prefabricated drains[C]//Proceedingsofthe10thInternationalConferenceonSoilMechanicsandFoundationEngineering. Stockholm, Sweden, 1981: 677-682.

[13] Butterfield R. A natural compression law for soils [J].Géotechnique, 1979,29(4): 469-480.

[14] 劉維正. 飽和粘性土原位力學特性分析與應用研究[D]. 南京:東南大學交通學院,2011.

[15] 韓文君. 劈裂真空法加固軟土地基理論與工程應用研究 [D]. 南京:東南大學交通學院,2012.