Undrained unified solutions to cylindrical cavity expansion in soils and sands based on CSUH model
-
摘要: 目前大多解答难以合理考虑密砂和超固结土的诸多力学特性,因而解答与实际情况存在一定偏差,此外可以统一描述砂土和黏土不排水柱孔扩张的解答很少。基于CSUH本构模型(Unified hardening model for clays and sands)推导了砂土和黏土中不排水柱孔扩张的弹塑性刚度矩阵[Dep]。在此基础上,联合大应变理论推导了不排水扩孔的控制方程,最后通过数值方法得到了方程的半数值半解析解。计算结果表明,预测结果不仅可以合理模拟柱孔扩张过程中松砂的液化和正常固结土的硬化,还可以模拟密砂和超固结土的剪胀、峰值强度、密实状态的衰化和特征状态等特性,因此可以合理模拟砂/黏土柱孔扩张过程中土体应力场和应变场的变化规律。Abstract: The current solutions to cylindrical cavity expansion cannot properly consider the mechanical properties of dense sand and overconsolidated soil. Therefore, there is some discrepency between the solution and the practical situation. In addition, there are few answers that can describe the expansion of undrained expansion in sand and clay uniformly. Based on the CSUH model, the elastoplastic matrix and the governing equations for undrained cylindrical cavity expansion are proposed. Based on the boundary conditions, the exact solutions to the governing equations are obtained numerically. The results show that the predicted results capture the stress−strain relationship more reasonably, including liquefaction, shear dilatancy, strain-hardening and softening, damping of dense state and overconsolidation ratio. Therefore, the proposed approach can reasonably simulate the variation laws of stress field and strain field during the expansion process in dense sand and overconsolidated soil.
-
Keywords:
- cylindrical cavity expansion /
- sand/clay /
- CSUH model /
- undrained solution
-
0. 引言
水下抛填散粒料常被用来填筑围堰基础,例如白鹤滩工程围堰、大藤峡二期深水围堰、碾盘山围堰等。水下抛填散粒料的相对密度是围堰和大坝工程设计中最基本的物理参数之一,直接决定了抛填料的强度参数和变形模量,是进行工程变形和稳定分析的关键参数[1],但是如何准确、合理地确定抛填相对密度仍存在困难[2]。
目前,国内外许多学者研究了典型散粒料的抛填密度问题。程永辉等[3]采用静态抛填离心试验方法研究了风化砂、石渣料和砂砾石的水下抛填密度,指出进占式水下抛填深度越大,由于抛填料运动至深部动能也较大,且受后期抛填料的压实作用也越大,因此形成的密度较大。李青云等[4]通过对三峡工程二期深水高土石围堰运行后的实测分析,一般位置的水下部分(66 m高程以下)堰体密度实测值在1.83~1.98 t/m3之间,与用离心模型试验确定的水下抛填风化砂密度(1.75~1.85 t/m3)相比,略有增长,考虑的围堰4 a运行中的固结效应,二者基本一致,有力证明了采用离心模型试验方法确定的风化砂水下抛填密度的科学性和可行性。Mahmoudi等[5]通过开展一系列重塑饱和砂土试样的不同超固结比不排水三轴试验,结果表明初始相对密度以显著的方式控制不排水抗剪强度,初始相对密度的增加导致膨胀性的放大,从而导致不排水抗剪强度的增加[6]。Wei等[7]采用砂子和两种不同粒径的砂砾石混合,砂砾石的粒径分别为5~10,10~20 mm,两种混合料的最大密度和最小密度随着砂砾石含量的增大均逐渐增大,当砂砾石的含量为75%时达到最大值,砂砾石含量相同时粒径较大的砂砾石得到的密度较大。Chang等[8]对松散状态的宽颗粒级配土,通过不同含砾量的间隙级配砂和级配砂的人工配制试样开展系列排水单剪试验,表明砾石土抗剪强度取决于优势颗粒的充填条件,其中颗粒间孔隙比比整体孔隙比更能代表土体的充填条件。Kalman[9]对比分析不同细颗粒掺入粗颗粒时密度变化,混合料的最大密度和最小密度与细颗粒的粒径以及掺入量密切相关,密度基本随着细颗粒掺入量的增加先增大后减小。Chen等[10]指出颗粒材料的岩土力学性能严格依赖于其初始堆积密度、应力水平和粒度分布,二元颗粒土的性能还受其细小颗粒的质量和数量的影响,峰值抗剪强度随细小颗粒含量的变化呈钟形,与初始填充密度和应力水平无关,与粒径比值有关。因此,颗粒材料的抛填密度决定了其物理力学特性,直接影响了围堰构筑物的变形和稳定,但现有研究无法得到统一的研究结论或计算方法,不同区域的抛填料密度跟颗粒大小和级配直接相关。
以汉江某水电站大坝水下抛填工程为原型条件,开展砂砾石和沙子混合料的进占式水下抛填密度离心试验,将相对密度Dr作为衡量标准,主要研究散粒料颗粒级配等因素对水下抛填密度的影响规律。
1. 单戗立堵水下抛填离心试验设计
1.1 原型条件和概化
汉江某水电站项目为Ⅱ等大(2)型工程,二期主河床土石坝最大坝高为22.58 m,且坝基需要在水下进行抛投填筑,最大水深为12 m。水下抛填料主要采用汉江河床内的散粒料,为砂砾石和砂子混合料。通过现场筛分得到典型的颗粒级配曲线如图 1所示。
现场筛分试验结果表明填料的最大粒径小于100 mm,0.5 mm以下的含量分别为24.6%,32.1%,49.1%,而0.5~2 mm之间的含量分别为1.9%,1.8%,2.1%,计算得到三种抛填料的不均匀系数Cu分别为153,169,53,曲率系数Cc分别为4.49、1.14和0.11。其中,不均匀系数Cu是限制粒径与有效粒径的比值,是反映组成土的颗粒均匀程度的一个指标。曲率系数Cc是反映土的粒径级配累计曲线的斜率是否连续的指标系数。当Cu≥5且Cc=1~3的土为级配良好,现场筛分得到的3种砂砾石抛填料颗粒级配曲线。表明XJ2下游围堰为级配良好,而XJ1上游围堰和XJ3下游围堰料场均为级配不良。
1.2 试验方案
共开展3种抛填料的水下抛填试验,抛填料主包括:XJ-1、XJ-2和XJ-3分别为上游围堰、下游围堰及下游围堰附近的砂砾石。研究内容主要为现场实际12 m水深条件以及12 m水深+14.3 m上覆堆载的水下抛填密度;对比分析不同水深(5~20 m)的水下抛填密度规律。单戗立堵水下抛填离心试验方案如表 1所示。
表 1 单戗立堵水下抛填离心试验方案Table 1. Centrifugal test plans for single embankment vertical blocking underwater dumping and filling编号 抛填料 试验内容 XJ-1 砂砾石,上游围堰级配 设定最大离心加速度40g,模拟原型进占式抛填方式
试验工况:12 m水深条件下水下抛填密度;12 m水深+14.3 m上覆堆载下的水下抛填体密度;不同水深(5~20 m)的水下抛填密度变化规律XJ-2 砂砾石,下游围堰级配 XJ-3 砂砾石,下游围堰堤头附近级配 1.3 模型抛填料
由于原型填料的粒径尺寸较大,在离心模型试验中不能按原型进行试验,需进行缩尺模拟。综合考虑模型箱和模型断面尺寸,以模型箱最小尺寸的1/20即2 cm作为离心模型试验中砂砾石料的最大限制粒径。本次试验采用等量代替法模拟原型抛填料,即用2~20 mm的各级土料等量代替原型料中20 mm以上的填料,2 mm以下用原型料,模型抛填料级配曲线如图 2所示。
对于粗粒土的相对密度试验,最大密度试验采用振动台法,最小密度试验采用固定体积法,试验结果如表 2所示。3种模型抛填料的最大密度和最小密度分别为2.12,1.88 g/cm3;2.10,1.85 g/cm3;1.99,1.73 g/cm3。
表 2 模型抛填料的最大密度和最小密度Table 2. Maximum and minimum densities of model filling materials编号 抛填料 最大密度/(g·cm-3) 最小密度/(g·cm-3) XJ-1 砂砾石,上游围堰级配 2.12 1.88 XJ-2 砂砾石,下游围堰级配 2.10 1.85 XJ-3 砂砾石,下游围堰堤头附近级配 1.99 1.73 1.4 离心场中水下抛填装置
试验在长江科学院CKY-200型现代化多功能土工离心机上完成。离心机的主要性能参数:有效容量为200g·t;最大加速度200g,无级调速,调速精度0.1g;有效半径3.7 m;吊蓝净空尺寸1.2 m×1.0 m×1.5 m;模型箱尺寸(长×宽×高)为100 cm×40 cm×80 cm。
试验研制了单戗立堵水下抛填离心试验装置,能够实现离心机运行(离心加速度不大于100g)过程中模拟水下抛填的单戗立堵进占方式,断面设计图如图 3(a)所示,模型图如图 3(b)所示。当离心加速度逐级增大至设计加速度时,启动电动提升装置,通过钢丝绳和滑轮连接,将模型箱中部隔板缓慢拉出,隔板上部填料随着开口增大依次落入模型箱下部的水中。该装置主要特点是通过隔板与滑槽之间采用滑轮接触,尽可能减少摩阻力;设置溢流孔和集水箱,保证抛填模拟过程中水位保持不变。
2. 水下抛填离心试验结果及分析
通过离心试验得到的三种抛填料水下抛填密度结果如表 3和图 4所示。表 3为单戗立堵水下抛填后得到平均相对密度,图 4为12 m水深以及12 m水深+14.3m厚上覆堆载条件相对密度-深度的关系。本次试验抛填密度采用灌水法测得,首先通过预留排水管将模型中的水排出,依次将5,12,15,20 m四个深度处一定厚度的抛填料取出,采用灌水法测量取出抛填料的体积,计算得到不同深度处平均干密度,进而计算相对密度Dr。相对密度Dr是指无黏性土处于最松状态的孔隙比与天然状态或给定孔隙比之差和最松状态孔隙比与最紧孔隙比之差的比值。
表 3 单戗立堵水下抛填离心试验结果Table 3. Centrifugal test results for single embankment vertical blocking underwater dumping and filling编号 抛填料 水下抛填密度 12 m水深条件 12 m水深+14.3 m上覆堆载 干密度/ (g·cm-3) 相对密度 干密度/(g·cm-3) 相对密度 XJ-1 砂砾石,上游围堰级配 1.951 0.33 2.012 0.58 XJ-2 砂砾石,下游围堰级配 1.936 0.36 2.016 0.68 XJ-3 砂砾石,下游围堰堤头附近级配 1.788 0.27 1.852 0.52 表 3分析表明,水下抛填密度与颗粒级配是否良好密切相关,12 m水深条件下XJ1~XJ3水下抛填相对密度分别为0.33,0.36,0.27,基本都处于松散~中密的过渡阶段。其中,XJ2的水下抛填密度相对大一些,主要原因是XJ1和XJ3的曲率系数Cc分别为4.49和0.11,为级配不良,而XJ2级配良好。当增加14.3 m厚上覆荷载后,XJ1~XJ3水下抛填相对密度均增大,分别为0.58,0.68和0.52,XJ2达到密实状态,上覆荷载作用能够有效增大抛填区域尤其是浅部抛填料的相对密度。
图 4表明一定水深条件下抛填料的相对密度随着深度而逐渐增大,与包承纲[1]、程永辉等[3]得到的规律一致。XJ2相对密度显著均大于XJ1和XJ3,主要是因为XJ2级配良好,粒径尺寸逐级减小,较小粒径能够尽可能填补上一级粒径的孔隙。当抛填料上部增加14.3 m上覆堆载时,水下抛填体的相对密度均增大,上部增幅较大,可达0.4以上,随深度增幅逐渐减小,深度20 m时增幅仅为0.03左右。结果表明上覆荷载作用可有效增大相对密度,尤其对于水下抛填形成相对密度较小的浅部抛填料增幅最为显著。当抛填料的级配良好时,水下抛填联合上部堆载作用可使得抛填体达到密实状态(Dr>0.67)。
3. 散粒料抛填相对密度简化分析
基于二元颗粒“填隙”理论[11],建立了散粒料抛填相对密度简化计算方法,主要考虑抛填料级配、抛填深度、上覆堆载等影响因素。两粒径理想颗粒的孔隙比被定义为
e=e0+Δeh+Δep−kfill。 (1) 则相对密度Dr:
Dr=emax−eemax−emin。 (2) 式中,e0为初始孔隙比,为最大粒径组的初始孔隙比;Δeh和Δep分别为抛填料深度和上覆荷载作用引起的孔隙比变化,可根据固结试验e-p曲线得到;kfill为颗粒级配(不均匀系数Cu和曲率系数Cc)引起的孔隙比变化。
本次离心模型试验填料为砂砾石,测试得到颗粒相对质量密度Gs为2.65;最大粒径组为10~20 mm,相对密度约为0.38;级配填料的最大孔隙比emax和最小孔隙比emin如表 4所示。
表 4 模型抛填料的最大和最小孔隙比Table 4. Maximum and minimum void ratios of model filling materials编号 抛填料 emax emin XJ-1 砂砾石,上游围堰级配 0.41 0.25 XJ-2 砂砾石,下游围堰级配 0.43 0.26 XJ-3 砂砾石,下游围堰堤头附近级配 0.53 0.33 根据填料初始为最松散状态的固结试验e-p曲线绘制孔隙比随压力的关系曲线,得到相对密度-竖向压力的关系,如图 5所示。
参照Rosin-Rammler采用指数函数表示粒度分布的关系式,粒度分布宽度w为
w=lg(DmaxDmin)=lgDen√ln100R(Dmax)Den√ln100R(Dmin)。 (3) 式中:De为特征粒径;R为筛余量;Dmax和Dmin分别为最大和最小粒度。XJ2级配良好,颗粒级配引起的孔隙比变化为0.05;XJ1和XJ3级配不良,其值分别为0.14和0.12。
采用本文提出的散粒料抛填相对密度简化分析方法得到的相对密度与离心试验结果对比如图 6所示。结果表明,计算值与离心试验结果基本一致,最大相对误差小于5%。
4. 结论
(1)通过单戗立堵水下抛填砂砾石离心模型试验,研究了汉江某水电站典型级配条件水下抛填密度分布规律和主要影响因素,得到抛填密度随深度的分布规律,分析了颗粒级配特征(不均匀系数Cu和曲率系数Cc)对密度的影响。
(2)研制了离心场中单戗立堵水下抛填试验装置,该装置使用最大加速度为100g,能够模拟最大抛填水深为30 m,水平抛填范围40 m,实现了抛填过程中水深条件保持不变。
(3)对于汉江典型砂砾石抛填料,颗粒级配条件对抛填相对密度起着控制作用。当抛填料级配良好时,相对密度较大,当抛填料上部施加一定的堆载后,相对密度会得到进一步提高,平均值可达到相对密实状态。研究成果为围堰和土石坝中水下抛填散粒料的设计和施工提供试验依据。
-
表 1 Toyoura砂参数
Table 1 Parameters for Toyoura sand
M N σ′r0 /kPaσ′θ0 /kPaσ′z0 /kPaK0 e0 1.25 1.973 1000 1000 1000 1 0.933 1.25 1.973 1000 1000 1000 1 0.856 1.25 1.973 850 850 1250 0.68 0.735 表 2 不同种类黏土参数
Table 2 Parameters for different types of clay
OCR σ′r0 /kPaσ′θ0 /kPaσ′z0 /kPaK0 V0 cu/kPa 1.0 100 100 160 0.625 2.09 54.3 1.2 100 100 160 0.625 2.06 62.8 3.0 120 120 120 1.000 1.97 115.0 10.0 144 144 72 2.000 1.80 360.2 -
[1] VESIC A S. Expansion of cavities in infinite soil mass[J]. Journal of the Soil Mechanics and Foundations Division, 1972, 98(SM3): 265-269.
[2] RANDOLPH M F, CARTER J P, WROTH C P. Driven piles in clay-the effects of installation and subsequent consolidation[J]. Géotechnique, 1979, 29(4): 361-393. doi: 10.1680/geot.1979.29.4.361
[3] SHEN S L, MIURA N, KOGA H. Interaction mechanism between deep mixing column and surrounding clay during installation[J]. Can Geotech J, 2003, 40(2): 293-307. doi: 10.1139/t02-109
[4] SALGADO R, MITECHELL J K, JAMIOLKOWSKI M B. Cavity expansion and penetration resistance in sand[J]. Journal of Geotechnical and Geoenviromental Engineering, ASCE, 1997, 123(4): 344-354. doi: 10.1061/(ASCE)1090-0241(1997)123:4(344)
[5] CHANG M F, TEH C I, CAO L F. Undrained cavity expansion in modified Cam clay II: application to the interpretation of the piezocone test[J]. Géotechnique, 2001, 51(4): 335-350. doi: 10.1680/geot.2001.51.4.335
[6] CARTER J P, BOOKER J R, YEUNG S K. Cavity expansion in cohesive frictional soils[J]. Geotechnique, 1986, 36(3): 349-358. doi: 10.1680/geot.1986.36.3.349
[7] YU H S, HOULSBY G T. Finite cavity expansion in dilatants soils: loading analysis[J]. Géotechnique, 1991, 41(2): 173-183. doi: 10.1680/geot.1991.41.2.173
[8] MANTARAS F M, SCHNAID F. Cylindrical cavity expansion in dilatant cohesive-frictional materials[J]. Géotechnique, 2002, 52(5): 337-348. doi: 10.1680/geot.2002.52.5.337
[9] CAO L F, TEH CI, CHANG M F. Undrained cavity expansionin modified Cam clay I: theoretical analysis[J]. Géotechnique, 2001, 51(4): 323-334. doi: 10.1680/geot.2001.51.4.323
[10] CAO L F, TEH C I, CHANG M F. Analysis of undrained cavity expansion in elasto-plastic soils with non-linear elasticity[J]. Int J Numer Anal Meth Geomech, 2002, 26(1): 25-52. doi: 10.1002/nag.189
[11] SILVESTRI V, ABOU-SAMRA G. Application of the exact constitutive relationship of modified Cam clay to the undrained expansion of a spherical cavity[J]. Int J Numer Analyt Meth Geomech, 2011, 35(1): 53-66. doi: 10.1002/nag.892
[12] VINCENZO S, GHASSAN A S. Analytical solution for undrained plane strain expansion of a cylindrical cavity in modified cam clay[J]. Geomech Eng, 2012, 4(1): 19-37. doi: 10.12989/gae.2012.4.1.019
[13] CHEN S L, ABOUSLEIMAN Y N. Exact undrained elasto-plastic solution for cylindrical cavity expansion in modified Cam clay soil[J]. Géotechnique, 2012, 62(5): 447-456. doi: 10.1680/geot.11.P.027
[14] CHEN S L, ABOUSLEIMAN Y N. Exact drained solution for cylindrical cavity expansion in modified Cam clay soil[J]. Géotechnique, 2013, 63(6): 510-517. doi: 10.1680/geot.11.P.088
[15] ZHOU H, LIU H L, KONG G Q. Analytical solution of undrained cylindrical cavity expansion in saturated soil under anisotropic initial stress[J]. Comput and Geotech, 2014, 55: 232-239. doi: 10.1016/j.compgeo.2013.09.011
[16] ZHOW H, LIU H L, ZHA Y H, et al. A general semi-analytical solution for consolidation around an expanded cylindrical and spherical cavity in modified Cam Clay[J]. Comput and Geotech, 2017, 91: 71-81. doi: 10.1016/j.compgeo.2017.07.005
[17] 李镜培, 唐剑华, 李林, 等. 饱和黏土中柱孔三维弹塑性扩张机制研究[J]. 岩石力学与工程学报, 2016, 35(2): 378-386. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201602019.htm LI Jing-pei, TANG Jian-hua, LI Lin, et al. Mechanism of three dimensional elastic-plastic expansion of cylindrical cavity in saturated clay[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(2): 378-386. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201602019.htm
[18] LI L, LI J P, SUN D A. Anisotropically elasto-plastic solution to undrained cylindrical cavity expansion in K0-consolidated clay[J]. Comput and Geotech, 2016, 73: 83-90. doi: 10.1016/j.compgeo.2015.11.022
[19] LI J P, GONG W B, LI L, et al. Drained elastoplastic solution for cylindrical cavity expansion in K0-consolidated anisotropic soil[J]. J Eng Mech, 2017, 143(11): 04017133.
[20] CHEN S L, LIU K. Undrained cylindrical cavity expansion in anisotropic critical state soils[J]. Géotechnique, 2019, 69(3): 189-202. doi: 10.1680/jgeot.16.P.335
[21] LIU K, CHEN S L. Analysis of cylindrical cavity expansion in anisotropic critical state soils under drained conditions[J]. Can Geotech J, 2019, 56(5): 675-686. doi: 10.1139/cgj-2018-0025
[22] CHEN H H, LI L, LI J P, et al. Elastoplastic solution to drained expansion of a cylindrical cavity in anisotropic critical-state soils[J]. J Eng Mech, 2020, 146(5): 04020036.
[23] CHEN H H, LI L, LI J P. Stress transform method to undrained and drained expansion of a cylindrical cavity in anisotropic modified Cam-Clay soils[J]. Comput and Geotech, 2019, 106: 128-142. doi: 10.1016/j.compgeo.2018.10.016
[24] SUN D A, SHENG D, SLOAN S W. Elastoplastic modeling of hydraulic and stress-strain behaviour of unsaturated soils[J]. Mechanics of Materials, 2007, 39(3): 212-221. doi: 10.1016/j.mechmat.2006.05.002
[25] YAO Y P, LIU L, LUO T, et al. Unified hardening (UH)model for clays and sands[J]. Comput and Geotech, 2019, 110: 326-343. doi: 10.1016/j.compgeo.2019.02.024
[26] YAO Y P, HOU W, ZHOU A N. UH model: three-dimensional unified hardening model for overconsolidated clays[J]. Géotechnique, 2009, 59(5): 451-469. doi: 10.1680/geot.2007.00029
[27] YAO Y P, GAO Z W, ZHAO J D, et al. Modified UH model: constitutive modeling of overconsolidated clays based on a parabolic Hvorslev envelope[J]. J Geotech Geoenviron Eng, 2012, 138(7): 860-868. doi: 10.1061/(ASCE)GT.1943-5606.0000649
[28] MATSUOKA H, SUN D A. The SMP Concept-Based 3D Constitutive Models for Geomaterials[M]. Taylor and Francis, 2006.
[29] 姚仰平. UH模型系列研究[J]. 岩土工程学报, 2015, 35(2): 193-216. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201502002.htm YAO Yang-ping. Advanced UH models for soils[J]. Chinese Journal of Geotechnical Engineering, 2015, 35(2): 193-216. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201502002.htm
[30] VERDUGO R, ISHIHARA K. The steady state of sandy soils[J]. Soils Found, 1996, 36(2): 81-91. doi: 10.3208/sandf.36.2_81
[31] ISHIHARA K. Liquefaction and flow failure during earthquakes[J]. Géotechnique, 1993, 43(3): 351-415. doi: 10.1680/geot.1993.43.3.351