Influences of excavation of foundation pit of underground pipe gallery on mechanical deformation characteristics of CFG piles
-
摘要: 在先施工工程桩再进行开挖的基坑工程中，坑底土体在基坑开挖卸荷的过程中产生的回弹变形会导致基坑底部土体产生相对于桩体向上的位移，此时桩体可能因较大的拉力断裂，从而大大减弱地基承载力。以贵州某地下综合管廊项目为依托，利用有限元计算软件MIDAS GTS NX模拟了基坑开挖的施工过程，研究了桩长和桩体弹模对CFG桩受力及变形的影响规律。结果表明：桩长对CFG桩受力及变形有着较为明显的影响，而桩长超过一定范围后影响显著降低；增加桩体弹模会让桩体承受更大的拉应力，但是能使桩体处于更加安全的状态。Abstract: In the foundation pit in which the engineering piles are constructed first and then excavated, the rebound deformation of the soil at the bottom of the pit in the process of excavation and unloading will lead to the upward displacement of the soil at the bottom of the pit relative to the piles. At this time, the piles may be broken due to large tensile force, which will greatly reduce the bearing capacity of the foundation. Based on an underground comprehensive pipe gallery project in Guizhou, the construction process of excavation of foundation pit is simulated by using the finite element software MIDAS GTS NX, and the influence laws of length and elastic modulus of the piles on the stress and deformation of the CFG piles are studied. The results show that the length of the piles has obvious influences on the stress and deformation of the CFG piles, and the influences decrease significantly when the pile length exceeds a certain range. Increasing their elastic modulus will make the piles bear greater tensile stress, but it can make them be in a safer state.
-
0. 引言
黄土在中国西北地区覆盖广泛,随着西部大开发的进展和“一路一带”倡议的落实,黄土地区的重大工程不断出现,与此同时也带来了一系列的工程问题。黄土作为一种典型的特殊土,具有复杂的力学特性,因此建立恰当的本构模型预测黄土的强度变形特性是解决工程问题的基础和关键；而且本构模型的研究也是土力学研究的热点和难点问题之一。
国内外学者对黄土的本构模型进行了大量研究,建立黄土的弹塑性本构模型[1]、损伤本构模型[2]、基于Desai的扰动状态本构模型[3],建立相应的黄土扰动状态本构模型[4]、破损力学模型[5]。而其中沈珠江所提出的破损力学模型[5]把土的复杂结构简化为结构性较强的结构块和结构性较弱的软弱带来描述,荷载由结构块和软弱带二者承担,从而将土视为一种二元介质材料,进而建立一种新的二元介质模型,进而描述土变形过程中实际的破损现象。国内外学者已将该模型成功的应用于黄土[6-9],这些模型的建立为更进一步的研究提供了思路。
统一强度理论[10]不仅可以反映材料的三维强度特性,覆盖现有强度准则的上限和下限,并且可以线性逼近很多的强度准则。在本构模型中考虑统一强度理论,无疑将使得模型有着更为广泛的适用性。因此本文将分析黄土强度演化过程,在研究黄土受载作用下强度的演变规律的基础上,结合二元介质模型,引入统一强度理论,建立黄土的统一二元介质本构模型。
1. 统一强度理论的不变量表达式
统一强度理论在岩土工程领域应用广泛,该强度理论从双剪单元体出发,考虑双剪应力单元体上的所有剪应力及其对应的主应力对材料破坏不同的贡献[10],其统一强度理论的应力不变量表达式为
(1) 当
(2a) 当
(2b) 式中,b为中间主应力系数；c, φ分别为黏聚力和内摩擦角。统一强度理论不再是一个单独的强度准则,它不仅能够解释材料的破坏机理,而且当中间主应力系数b变化时,能够代表一系列的强度准则。
2. 黄土二元介质模型
2.1 黄土二元介质模型理论
黄土在变形过程中,具有明显的非线性性质,在变形过程中同时伴随着弹性、塑性及脆性等力学行为特性。为描述黄土的这些力学行为特性,根据各自不同的力学性质,按照二元介质模型理论[6],可以假定采用弹簧、滑片以及滑片等基本元件描述弹性、塑性及脆性等力学性质,采用不同的力学元件组合,表示不同的介质行为,从而解释材料的不同力学行为特性。
黄土的强度可以认为由摩擦强度和黏聚强度组成,研究表明黏聚力和摩擦力在变形中不会同步发生作用,只有在黏聚作用破坏后,摩擦强度才逐渐发挥作用；摩擦强度和黏聚强度所发挥作用的强弱,决定了黄土的应变硬化、应变软化等力学行为特性。借助二元介质模型,可以分别采用摩擦元和胶结元分别表示黄土强度的构成。在荷载作用下,胶结元率先发挥作用,随着变形的发展,其荷载的承担逐渐向摩擦元过渡。
因此可以借助二元介质模型进一步将黄土抽象成由弹脆性元(胶结元的结构体)和弹塑性元(摩擦元的破碎体)组成的二元介质型材料。在加载变形过程中,一个微单元体积内弹脆性元逐渐破裂并向弹塑性元转化。弹脆性元向弹塑性元转变的过程可用破损参数B描述,B在黄土变形过程中,可以表示为二元介质荷载的分担比,类似于损伤力学中的损伤变量,也可以理解为胶结元的破损程度,因此其取值范围为0≤B≤1,假定材料未受到损伤,则B=0；若B=1,则表示材料内部损伤演化到极限状态。
结构性岩土材料在变形过程中,弹脆性元向弹塑性元转化,在一个代表性微单元体积内应用非均匀介质的均匀化理论,可以得到单元内应力情况。按照等效应变的假设,二元介质模型的应力关系可以表示为
(3) 式中,上标α表示黄土实际变形或应力；下标1和2分别表示弹脆性元和弹塑性元。其增量形式可以表示为
(4) 由式(4)可知单元总应力增量
2.2 弹脆性元部分本构关系
黄土在变形过程中,黏聚力首先发挥作用,其力学行为采用弹性理论进行描述。因此,认为弹脆性元在单元内未破损之前,只有弹性变形,其力学特性按弹性元件受力而发生形变,即符合Hooke定律,则弹脆性元的弹性应力应变增量关系为
(5) 式中,
2.3 弹塑性元部分本构关系
当弹脆性元破损之后,其荷载逐渐由弹塑性元进行承担。假定弹塑性元的变形和一般岩土材料的弹塑性变形一样,包括弹性变形和塑性变形：
(6) 弹性变形可根据胡克定律求得,塑性变形可以按照传统岩土塑性理论进行求解。考虑到统一强度理论在岩土中有较为广泛的应用,且可以反映中间主应力的影响,因此在建立弹塑性元力学模型时,采用统一强度理论确定其本构模型。
根据统一强度理论,弹塑性元的屈服准则可以表示为
(7) 式中,ε为硬化参数。
根据硬化准则可知,材料的屈服面不是固定不变的,材料在发生塑性变形后,其后继屈服面常发生改变,即发生硬化现象。但在本文中主要对二元介质模型进行研究分析,硬化参数的影响较为次要,为了突出二元介质模型的优越性,假定硬化参数为零,即不考虑黄土屈服的后续硬化行为,在变形过程中保持屈服面的大小不发生改变。
岩土材料是摩擦性材料,采用非相关联流动法则考虑剪胀性和体积变化的影响,引入剪胀角ψ2代替摩擦角φ2,得到塑性势函数：
(8) 式中,ψ2=ζφ2, (0≤ζ≤1),若ζ=0,表示无剪胀现象,若ζ=1,表示发生最大剪胀。一般情况下ψ2<φ2,即剪胀角小于摩擦角。
根据弹塑性理论,可以确定弹塑性元的本构模型为
(9) 当不考虑虑硬化行为时,式(9)中硬化函数H=0。
弹塑性元屈服面在偏平面上的极限线如图1所示,从图中可以看出,弹塑性元模型存在奇异点,其屈服函数和塑性势函数在偏平面上不是光滑曲线,因此在进行数值计算时,必须对该模型的奇异性进行处理。
常见的处理方法,采用矢量平均值来代替该奇异点的值,即其值为角点在左右不同边趋近所得值的均值,如当θσ=arctan(sinφ2/
此外,屈服函数和塑性势函数在θσ=-30°, 30°处由于中间主应力影响系数b的取值可能会发生奇异性。当0≤b<1时,将θσ=-30°, 30°代入相应的条件式中即可。当b=1时,屈服函数和塑性势函数在θσ=-30°, 30°的值趋于无穷大,为了消除这种奇异性的影响,采用对流动矢量取极限的数学方法。
2.4 破损参数的确定
黄土在加载过程中破坏是一个连续变化的过程,用微观物理量描述破损过程更能反映黄土的破坏。本文采用黄土微元体分析黄土的渐进破损过程,假定弹脆性元在为破损之前其应力应变关系服从广义胡克定律,其弹性模量为初始切线模。弹脆性元破损之后变为弹塑性元,假定弹脆性元的破坏分布是受应力、应变和破坏概率影响的,是一随机、连续的变化量,其弹破坏服从Weibull分布统计规律,则可以得到破损参数为
(10) 式中,A为反映结构破损速率的参数,m为试验参数。
2.5 黄土二元介质本构模型的建立
由式(3)根据非均匀介质的均匀化理论和破损参数的定义,可得到
(11) 式中,
因此,黄土的二元介质本构模型为
(12) 2.6 黄土二元介质模型参数的确定
根据二元介质损伤本构模型的建立过程可以确定在该本构模型中有3种类型的模型参数：弹脆性元及弹塑性元两种单元的模型参数、以及损伤函数参数中的相关参数。
弹脆性元的弹性参数可通过试样进行试验得到,弹脆性元弹性模量E1、泊松比μ1可通过常规三轴试验测得。弹塑性元的参数弹性模量E2、泊松比μ2和内摩擦角φ2可通过重塑土的常规三轴试验测得,当缺少重塑土的三轴试验时,弹塑性元的弹性参数还可以通过弹脆性元的相应系数适当的折减得到,或者与弹脆性元保持一致。中间主应力系数b可根据验证的具体试验条件来确定。剪胀角ψ2可以通过相应的常规三轴试验结果确定。
对于损伤参数B的确定,损伤参数可通过试验得到的应力应变关系关系曲线和破损参数的定义求得,对其公式进行变换,曲线变换为直线求解,推导求得,将式(10)移项并取两次自然对数变换,可以得到
(13) 3. 本构模型的验证
为验证本文所建立的二元介质模型的可行性,采用饱和黄土的常规三轴CU试验进行验证。本次试验所用黄土土样采自西安市地铁四号线工程D4KC-9标段行政中心车站。根据工程地质勘察资料可知,所取的黄土试样为Q4黄土,对其分别进行了等向压缩和固结不排水三轴试验,其三轴压缩试验结果如图2所示。
在常规三轴压缩试验中,中间主应力系数b=0。对于所建立的二元介质模型,其弹塑性元部分本构不考虑硬化作用。模型弹性模量等参数采用三轴试验的初始切线模量作为弹脆性元的弹性模量,在这里不考虑弹塑性元和弹脆性元初始模量之间的差别,两种参考介质弹性模型认为相同。
破损函数B的参数A和m和围压相关,随围压的不同而改变,所以需要根据不同围压下常规三轴CU试验所得数据和确定破损函数参数的方法求得各个围压下损伤参数的值,如表1所示。
表 1 损伤参数的数值Table 1. Damage parameters of binary-medium model参数 围压/kPa 100 200 300 A 2.64 2.52 2.07 m 0.413 0.418 0.437 图2为二元介质模型和试验结果比较图。从图中可以看到二元介质模型理论结果和试验结果基本趋势一致,所建立的本构模型较好地预测了黄土的应力应变关系。
为更进一步比较分析二元介质模型在预测黄土强度变形特性的优势,结合黄土的强度变形特性和统一屈服面的特点,确定统一弹塑性本构模型的硬化函数,建立一个基于统一强度理论的弹塑性本构模型,其硬化准则表示为
(14) 该模型预测结果如图2中虚线所示,从图2中比较可以看出,考虑硬化的统一弹塑性本构模型在预测黄土的应力应变关系时误差较大,尤其是在变形的初始阶段,误差更大。
相比统一硬化模型模拟曲线,二元介质模型能更好地预测黄土的应力应变关系。无论是在变形的初始阶段还是黄土接近破坏的阶段,二元介质模型理论结果和试验结果都能较好地吻合,说明在二元介质模型比统一硬化模型具有更好的适用性。
分析原因,土体抗剪强度主要由黏聚力和摩擦力组成,在变形的初始阶段,土体相对比较完整,黏聚力对抵抗变形其主要作用,摩擦力发挥的作用较小。由于在统一硬化模型主要考虑了摩擦的硬化效应,未考虑黏聚力的影响,且初始摩擦力较小,计算理论强度较低,因此在初始阶段统一硬化的理论结果与试验曲线差别较大。而随着变形的增加,黏聚力逐渐丧失,摩擦力不断增大,此时摩擦力对变形的抵抗起到主要作用,因此在变形的后期阶段,统一硬化模型预测结果误差减小。而在二元介质模型中,用弹脆性元模拟黄土的黏聚作用,用弹塑性元模拟黄土的摩擦作用,变形过程用破损参数B来描述和量化弹脆性元向弹塑性元转化,来清晰描述土体破坏过程,因此二元介质模型能更好地预测黄土的应力应变关系。
4. 结论
(1)依据统一强度理论的应力不变量形式,确定了统一屈服准则形式。分析了黄土强度构成及其发挥特点,结合二元介质模型本构模型,确定了黄土二元参考介质模型及其破损参数,采用等效应变原理,确定参考介质应力和应变关系,建立了黄土的二元介质本构模型,并分析了模型参数的确定方法。
(2)采用黄土的三轴试验对所建立的二元介质模型进行了分析和验证,结果表明,所建立的模型能够较好的预测黄土的应力应变关系；为比较分析,进异同一步确定了黄土的统一硬化弹塑性本构模型,并采用试验进行分析,研究了二元介质模型和统一硬化弹塑性本构模型在预测土的强度变形特性的异同,结果表明二元介质模型考虑了黄土的黏聚作用和摩擦作用,能够更真实地模拟土的强度发挥过程,可以提供更好的预测结果。
-
![]()
图 4 基坑底部(土体)回弹量变化曲线
Figure 4. Variation curves of rebound amount at bottom of foundation pit (soil)
![]()
图 7 基坑横截面最大拉应力变化曲线
Figure 7. Variation curves of maximum tensile stress of cross section of foundation pit
![]()
图 8 基坑底部(土体)回弹量变化曲线
Figure 8. Variation curve of rebound amount at the bottom of foundation pit (soil)
![]()
图 10 基坑横截面最大拉应力变化曲线
Figure 10. Variation curves of maximum tensile stress of cross section of foundation pit
表 1 土层参数取值表
Table 1 Values of parameters of soil layers
土层 重度/(kN·m-3) Eref50/MPa Erefoed/MPa Erefur/MPa 泊松比ν 黏聚力c/kPa 摩擦角φ/(°) 杂填土 19.0 30 30 90 0.32 5 28 黏土(可塑) 18.0 43 43 129 0.34 25 18 黏土(软塑) 17.3 22 22 66 0.34 20 15 强风化泥岩 24.0 120 120 360 0.37 30 20 注：表中Eref50为三轴实验割线刚度；Erefoed为主压密加载试验的切线刚度；Erefur为卸载弹性模量。 
下载: 导出CSV
表 2 计算步骤设置一览表
Table 2 List of calculation steps
施工阶段 施工阶段描述 第1步 初始地应力平衡 第2步 施工基坑围护结构与CFG桩 第3步 施工冠梁 第4步 施工内支撑（位移清零） 第5~10步 开挖至-6 m（每步开挖1 m）
下载: 导出CSV
-
[1] 郭德强. 软土地区基坑回弹机理及工程桩受力变形研究[D]. 天津: 天津大学, 2017. GUO De-qiang. The study of pit resilience mechanism and behavior of pile in soft region[D]. Tianjin: Tianjin University, 2017. (in Chinese)
[2] IWASAKI Y, WATANABE H, FUKUDA M, et al. Construction control for underpinning piles and their behaviour[J]. Géotechnique, 1994, 44(4): 681–689. doi: 10.1680/geot.1994.44.4.681
[3] 朱火根, 孙加平. 上海地区深基坑开挖坑底土体回弹对工程桩的影响[J]. 岩土工程界, 2005(3): 43–46. https://www.cnki.com.cn/Article/CJFDTOTAL-YSJS20050300I.htm ZHU Huo-gen, SUN Jia-ping. Influence of soil rebound at the bottom of deep foundation pit in Shanghai area on engineering piles[J]. Geological Exploration for Non-Ferrous Metals, 2005(3): 43–46. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSJS20050300I.htm
[4] 郑刚, 刁钰, 吴宏伟. 超深开挖对单桩的竖向荷载传递及沉降的影响机理有限元分析[J]. 岩土工程学报, 2009, 31(6): 837–845. doi: 10.3321/j.issn:1000-4548.2009.06.004 ZHENG Gang, DIAO Yu, ZHENG Gang, DIAO Yu, NG C W W. Finite element analysis on mechanism of effect of extra-deep excavation on vertical load transfer and settlement of a single pile[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(6): 837–845. (in Chinese) doi: 10.3321/j.issn:1000-4548.2009.06.004
[5] 胡琦, 凌道盛, 陈云敏, 等. 深基坑开挖对坑内基桩受力特性的影响分析[J]. 岩土力学, 2008, 29(7): 1965–1970. doi: 10.3969/j.issn.1000-7598.2008.07.043 HU Qi, LING Dao-sheng, CHEN Yun-min, et al. Study of loading characters of pile foundation due to unloading of deep foundation pit excavation[J]. Rock and Soil Mechanics, 2008, 29(7): 1965–1970. (in Chinese) doi: 10.3969/j.issn.1000-7598.2008.07.043
[6] 张翔, 刘松玉, 吴恺, 等. 工程桩对基坑回弹变形影响的数值模拟分析[J]. 岩土工程学报, 2021, 43(增刊2): 11–14. doi: 10.11779/CJGE2021S2003 ZHANG Xiang, LIU Song-yu, WU Kai, et al. Numerical analysis of influences of engineering piles on rebound deformation of foundation pit[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(S2): 11–14. (in Chinese) doi: 10.11779/CJGE2021S2003
[7] 罗耀武. 大面积深开挖对抗拔桩承载性状的影响研究[D]. 杭州: 浙江大学, 2010. LUO Yao-wu. The influence of large and deep excavation on the bearing behavior of uplift piles[D]. Hangzhou: Zhejiang University, 2010. (in Chinese)
[8] 郑刚, 张立明, 刁钰. 开挖条件下坑底工程桩工作性状及沉降计算分析[J]. 岩土力学, 2011, 32(10): 3089–3096. doi: 10.3969/j.issn.1000-7598.2011.10.034 ZHENG Gang, ZHANG Li-ming, DIAO Yu. Analysis of working performance of piles beneath excavation bottom and settlement calculation[J]. Rock and Soil Mechanics, 2011, 32(10): 3089–3096. (in Chinese) doi: 10.3969/j.issn.1000-7598.2011.10.034
[9] 杜一鸣, 郑刚, 张立明. 考虑深基坑开挖效应的超长灌注桩桩身压缩综合系数的理论分析[J]. 岩土力学, 2014, 35(7): 2019–2028. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201407036.htm DU Yi-ming, ZHENG Gang, ZHANG Li-ming. Theoretical analysis of comprehensive coefficient of pile compression for super-long bored pile considering effect of deep excavation[J]. Rock and Soil Mechanics, 2014, 35(7): 2019–2028. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201407036.htm
[10] 苟尧泊, 俞峰, 杨予. 基坑开挖引起既有桩基侧摩阻力中性点位置分析[J]. 岩土力学, 2015, 36(9): 2681–2687. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201509032.htm GOU Yao-bo, YU Feng, YANG Yu. Locating neutral point of excavation-induced skin friction on existing piles[J]. Rock and Soil Mechanics, 2015, 36(9): 2681–2687. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201509032.htm
-
期刊类型引用(6)
1. 王番，郅彬，刘恩龙，王小婵，邓博团，李金华，张辉. 考虑局部胶结破损热力学行为的结构性黄土二元介质本构模型. 岩土力学. 2025(01): 97-109 . 
百度学术
2. 武颖利，任红磊，郭万里. 胶结砂砾石料剪切二元强度准则与本构模型研究. 岩土工程学报. 2025(03): 516-524 . 本站查看
3. 万旭升，颜梦宇，路建国，晏忠瑞，刘凤云. 冻结硫酸钠粉质黏土强度与本构模型研究. 冰川冻土. 2023(01): 201-211 . 百度学术
4. 张革，刘恩龙. 基于CT动态扫描的冻土细观二元介质本构模型. 岩土工程学报. 2023(09): 1888-1896 . 本站查看
5. 张树明，蒋关鲁，叶雄威，蔡俊峰，袁胜洋，罗斌. 基于破损参数简化的二元介质冻结粉细砂土本构模型. 岩土力学. 2022(07): 1854-1864 . 百度学术
6. 侯乐乐，翁效林，崔艺铖，韩巍巍，陈禹勋. K_0固结结构性黄土的各向异性本构模型. 岩石力学与工程学报. 2022(10): 2124-2134 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 152
- HTML全文浏览量: 35
- PDF下载量: 30
- 被引次数: 10
