Stresses and deformations of support structures of deep foundation pit based on centrifuge and numerical simulation
-
摘要: 依据实际基坑开挖和支护结构工程特性进行模型概化,得出适合的离心模型试验方案,采用离心模型试验分析超深基坑开挖过程支护结构工程的受力和变形分布规律特征。同时建立三维有限元模型对圆形基坑开挖过程中支护结构的变形进行数值模拟,并将试验和数值计算结果进行对比分析。结果表明:支护结构水平位移模式为两头小,中间大的“涨肚型变形”;并且随着基坑深度的加深,地下连续墙水平位移最大值点逐渐下移;地表沉降呈凹槽形沉降形式,随着开挖的进行,沉降槽底部向着远离基坑的方向发展;地连墙背后土压力变化值呈非线性,开挖初期,土压力沿深度变化很小,但随着开挖的进行,土压力变化量逐渐增大;数值分析与离心模型试验的开挖支护结构变形特征结果较相近。研究成果可为基坑开挖设计切实有效的支护结构提供科学依据。Abstract: Based on the actual characteristics of excavation and support structures of foundation pit, a suitable centrifugal model test scheme is obtained, and the centrifugal model test is used to analyze the distribution characteristics of stress and deformation of the support structures of the ultra-deep foundation pit during the excavation process. Simultaneously a three-dimensional finite element model is established to simulate the deformation of the support structures during the excavation process of a circular foundation pit, and the test and numerical results are compared and analyzed. The results show that the horizontal displacement mode of the support structures is the "bulge shaped deformation" with small ends and large middle. As the depth of the foundation pit increases, the point with the maximum horizontal displacement of the underground diaphragm wall gradually moves downwards. The surface settlement is in the form of grooves, and as excavation continues, the bottom of the settlement groove develops towards the direction away from the foundation pit. The variation of soil pressure behind the diaphragm wall is non-linear. At the early stage of excavation, the variation of soil pressure along depth is small, but as excavation continues, the variation of soil pressure gradually increases. As the excavation continues, the bending moment value gradually increases, reaching the maximum positive bending moment near the depth of 17 m, about 2200 kN·m/m. The deformation characteristics of the support structures obtained from the numerical analysis and centrifugal model tests are relatively similar. The research results contribute to providing scientific basis for proposing practical and effective protection measures of the support structures for excavation of foundation pits.
-
Keywords:
- deep foundation pit /
- excavation /
- support structure /
- deformation /
- centrifuge /
- numerical simulation
-
0. 引言
福州长乐国际机场位于长乐—漳州诏安断裂带,地震形势十分复杂[1]。机场场地在地震作用下的稳定性是影响机场建设和使用安全的重要因素之一[2]。该机场场地地下水位较高,土质多为砂土,饱和砂土对于动荷载具有特殊的敏感性,在动荷载作用下易发生强度降低甚至液化引起结构失稳。因此,分析土的动力特性参数(动剪切模量比、阻尼比)和抗液化能力是分析重要及重大工程场地地震效应必不可少的基础[3]。
自19世纪起,许多学者对于土体动剪切模量与阻尼比的影响因素进行了探讨,主要聚焦于围压、固结比、颗粒形状、细颗粒含量及颗粒破碎等方面[4]。土体动剪切模量比和阻尼比随剪应变的变化曲线的不确定性有重要影响[5]。在工程实践中,重塑土试样测得性质与原位土性质存在一定出入,但已有研究表明[6],对于动剪切模量随剪应变关系进行归一化处理,可消除颗粒级配与相对密度对G/Gmax - γ曲线的影响。孔宪京等[7]对不同级配堆石料展开大型动三轴试验,对比切模量比及阻尼比随剪应变变化范围,证实动三轴室内试验测定的准确性。蔡袁强等[8]分析了荷载频率在不同应力路径下对路基粗粒填料长期动力特性的影响。李瑞山等[9]分析了荷载频率等对土体剪切模量比、阻尼比发展规律的影响。
值得注意的是,除去砂土动力特性参数分析,砂土液化对场地稳定性的影响也不容忽视。Seed等[10]基于饱和砂土液化势的概念,结合室内动三轴试验与现场勘测数据,对土体液化进行评估。石兆吉等[11]收集唐山地震土层液化数据,建立了饱和砂土、粉土的临界剪切波速判别方法。陈文化等[12]基于中国《水利水电工程地质勘察规范》中标准贯入试验,建立了考虑结构特性的临界标准贯入锤击数计算公式,为场地设计和稳定性判别提供新思路。
本文利用室内动三轴试验,分析福州饱和砂土动剪切模量比和阻尼比变化特性及其与围压之间的关系,测定场地土体动强度;并基于剪切波速法及Seed简化法结合现场剪切波速试验,判别场地液化情况,为机场场地的抗震设计和防液化措施提供指导。
1. 室内动三轴试验
1.1 试验方法
本研究采用英国GDS公司生产的振动三轴试验系统,轴向振动频率为0.01~10 Hz;最大轴向荷载20 kN;动态位移100 mm,该仪器采用自动数据采集系统,能够精准捕捉各个应变范围内动力参数变化。
试验材料分别取自10个样点的不同深度,共计36组试样。试验步骤主要有制样、试样饱和、加载。试样直径为39.1 mm,高为80 mm。采用反压饱和方法,施加300~400 kPa反压使试样达到较高的饱和度,本试验中所有试样在饱和结束后B系数均大于0.95。试验采用固结不排水试验,试样的固结应力比Kc为1.0,围压大小根据取样深度确定并在试验过程中保持恒定。通过设定轴向应力加载,加载波形为正弦波,荷载振动频率采用1 Hz。土体在轴向受大主应力循环变化,完成一个周期循环后,土体应力应变曲线将形成一个滞回圈,根据土工试验规范计算土体的动力特性参数。
动弹性模量及阻尼比试验采用振动频率为1 Hz,加载波形为正弦波,采用逐级施加动应力幅的方法,每个试样采用5~6级动应力,后一级的动应力幅一般为前一级的2倍左右,每级的振动次数不大于10次。
动强度试验时,采用动荷载振动频率采用1 Hz,加载波形为正弦波,对每个试样试验控制动荷载的大小,使试样破坏周数分别为10周次、20~30周次、100周次左右。当孔隙水压力与围压相等试验停止。
根据土工试验规范[13],本文试样的土体动力特性参数计算公式如下:
动弹性模量Ed:
Ed=σdεd×100。 (1) 式中:σd为轴向动应力(kPa),εd为轴向动应变(%)。
阻尼比λ:
λ=Az4π As。 (2) 式中:Az为图 1中滞回圈ABCDA面积(cm2);As为三角形OAB面积(cm2)。
动剪切模量、动剪应变:
Gd=Ed2(1+ν), (3) γd=εd(1+ν)。 (4) 式中:ν为泊松比。
1.2 动剪切模量比及阻尼比试验结果分析
图 2分别为不同取样点处部分试样在不同围压条件下砂土G/Gmax - γ,λ - γ关系图。为更直观简洁的描述动剪切模量随剪应变的衰减关系,将动剪切模量G与最大动剪切模量Gmax归一化处理,可以看到,当初始相对密度一定时,围压越大,土体受到的侧向约束力越大,其抵抗剪切变形的能力也越强,试样的动剪切模量比也越大。随着剪应变的增加,砂土逐渐慢慢受到破坏,土体抵抗剪切变形的能力减弱,同时试样的强度也变小,动剪切模量比随动剪应变的增加而减小。阻尼比随着剪切应变γ的增加呈现三段式发展。当剪切应变小于0.0001时,孔隙变化较小,土体阻尼比变化不明显。剪应变从0.0001到0.001时,其阻尼比随动剪切应变迅速增长,当剪应变达到0.001时,阻尼比曲线逐渐趋于平缓,此时阻尼比接近最大值。在剪应变增大的过程中,试样孔隙逐渐变小,此时克服的颗粒间的摩擦力需要耗散更多能量。
1.3 动强度试验结果分析
图 3为不同取样点位、不同深度处试样动强度试验结果。AZK20、AZK32、SZK38和SZK32取样点位试样,在不同取样深度即不同围压下,同一取样点位土体液化应力比均随着振次增加而减小,但不同取样深度下的土体液化应力差别较大,可能是由于同一取样点位不同深度下试样级配、孔隙比、颗粒堆叠排列方式不同,从而引起液化应力比的变化。可以看出,围压并非影响土体液化应力比的单一因素,土体场地稳定性研究可进一步探讨级配、孔隙比、颗粒排列形式的影响。
2. 场地稳定性判别方法研究
2.1 剪切波速判别法
中国自唐山地震后就开始注意到用地震剪切波速预测液化势问题,宏观震害资料和现场剪切波速Vs测定值表明,是否液化与Vs密切相关;并且经验表明,剪切波速法预测液化势的结果与用标准贯入试验预测结果基本一致。石兆吉等[11]提出了以有效应力为零作为液化判别标准,结合剪切波速预测砂土液化势公式如下:
Vscr=198√αmaxg(ds−0.0133d2s)。 (5) 式中:当烈度为7,8,9度时,αmax分别取0.1,0.2,0.4;ds为土层深度。
基于上述理论,取地震烈度为7度,求得场地不同深度处饱和砂土层的临界剪切波速,将土体实测剪切波速曲线与临界剪切波速曲线绘制于图 4。对比发现在埋深为8 m时,SZK38、SZK32、ASZK40点位实测波速曲线与临界剪切波速曲线相交,表明在埋深8 m内场地土层不发生液化,在8~10 m范围内,砂土在地震作用下会发生液化。
2.2 地震抗液化剪应力法
抗液化剪应力法是Seed等[10]提出一个确定土层中任一点的剪应力的近似方法。通过假设深度为h的单元土柱为刚体,并已知地面最大地震加速度,则可计算出土柱底面上的地震剪应力,由于土体为非刚体,求出剪应力应大于实际最大剪应力,需对其进行修正。由于地震时产生的动应力亦随时间而不规则地变化。其中最大地震剪应力只作用于一瞬间,以它代表地震期间的地震剪应力是不合理的。为了便于将其与室内以等幅值动应力试验求得的抗液化剪应力作比较,将地震期间这种不规则的地震剪应力,简化成一种等效的、有一定循环作用次数的均匀剪应力。通常简称为等效循环均匀剪应力。根据动力试验分析,可以近似地取最大动剪应力的65%作为等效循环均匀剪应力。
(1) 确定等效地震剪应力τe,根据场地设计地震参数计算,可计算出不同取样深度土体在地震作用下等效循环剪应力。
τe=0.65Kd∑γΔhgαmax。 (6) 式中:Δh为土层厚度;Kd为深度修正系数;γ为土体天然重度;g为重力加速度;αmax为最大地面加速度。
(2) 确定抗液化剪应力τd,天然场地下土层抗液化剪应力。
τd=Cr⋅σ′(σd2σ0)N。 (7) 式中:τd为抗液化剪应力;Cr一般为0.54~0.58;σ′与N有关;σ′为有效上覆压力;(σd/2σ0)N为抗液化剪应力比,根据福州长乐市抗震设防烈度确定等效循环振次N。
(3) 场地液化判别。将计算得到的τe和τd进行比较。若τe>τd则可能发生液化,计算结果如表 1,表中AZK20、ASZK32等代表不同取样点位,L1、L2等代表不同取样深度。
表 1 抗液化剪力计算结果Table 1. Calculated results of liquefaction shear resistance编号 干重度/ (kN·m-3) 取样深度/m Kd Cr σ′ (σd2σ0)N 等效地震剪应力τe 抗液化
剪应力τd是否液化 AZK20-L1 1.46 5.60~6.00 0.945 0.56 83.95 0.222 3.867 10.421 否 AZK20-L2 1.46 7.50~8.00 0.920 0.56 102.54 0.168 5.411 9.675 否 AZK20-L3 1.46 9.50~10.0 0.880 0.56 114.85 0.192 7.287 12.356 否 ASZK32-L1 1.46 3.50~4.00 0.970 0.56 54.75 0.232 2.589 7.114 否 ASZK32-L2 1.46 6.50~7.00 0.932 0.56 98.55 0.222 4.478 12.234 否 ASZK32-L3 1.46 7.50~8.00 0.920 0.56 105.65 0.201 5.411 11.909 否 ASZK32-L4 1.46 8.50~9.00 0.900 0.56 110.25 0.188 6.373 11.588 否 SZK38-L1 1.46 3.50~4.0 0.970 0.56 54.75 0.164 2.589 5.029 否 SZK38-L2 1.46 5.50~6.00 0.945 0.56 83.95 0.154 3.867 7.226 否 SZK38-L3 1.49 8.50~9.00 0.900 0.56 110.25 0.124 6.376 7.664 否 SZK32-L1 1.40 3.50~4.00 0.970 0.56 52.50 0.210 2.483 6.174 否 SZK32-L2 1.40 5.50~6.00 0.945 0.56 80.50 0.222 3.709 9.993 否 SZK32-L3 1.40 8.50~9.00 0.900 0.56 106.05 0.195 6.296 11.581 否 ASZK40-L1 1.44 3.50~4.00 0.970 0.56 54.00 0.209 2.554 6.327 否 ASZK40-L2 1.44 6.5~7.00 0.932 0.56 97.20 0.299 4.416 16.260 否 ASZK40-L3 1.44 7.50~8.00 0.920 0.56 104.10 0.352 5.342 20.544 否 ASZK40-L4 1.44 8.50~9.00 0.900 0.56 108.50 0.388 6.296 23.588 否 通过对比剪切波速法和抗液化剪应力法的判别标准可知,在取样深度小于8 m,其预测结果具有良好的一致性,对于深层土体仍需进一步判别。
3. 结论
本文通过对福州长乐国际机场为例,对现场10个取样点不同深度开展了36组动三轴试验,并对其中5个点位开展了动强度试验。并采用室内试验和现场测试两种方法探究场地砂土的抗液化能力,主要得到以下3点结论。
(1) 不同取样点试样在围压50 kPa,动剪切模量比约为0.23,阻尼比约为0.17,当围压增大为100 kPa和150 kPa时,动剪切模量比较为分散,土体动力特性随取样深度变化而发生。
(2) 不同取样深度下围压并非土体动强度单一影响因素,试样级配、相对密实度、颗粒堆叠方式等对土体动力特性的影响不容忽略。相同级配下,围压越大,颗粒间接触和排列更为紧密,其抗液化能力越大。
(3) 基于剪切波速法和抗液化剪应力法对场地进行判别,对于埋深小于8 m,判别结果具有良好的一致性;对于较深土层,两种判别结果存在一定出入,适用于较深土层液化判别方法亟待研究。
-
[1] 刘俊岩. 深基坑工程[M]. 北京: 中国建筑工业出版社, 2001. LIU Junyan. Deep Foundation Pit Engineering[M]. Beijing: China Architecture & Building Press, 2001. (in Chinese)
[2] 王莹, 夏才初, 陈孝湘, 等. 基于离心试验下软黏土地区窄基坑稳定性分析[J]. 地下空间与工程学报, 2022, 18(6): 1853-1862, 1872. WANG Ying, XIA Caichu, CHEN Xiaoxiang, et al. Stability analysis of narrow foundation pit in saturated soft clay area based on centrifugal test[J]. Chinese Journal of Underground Space and Engineering, 2022, 18(6): 1853-1862, 1872. (in Chinese)
[3] 马险峰, 曹明洋. 上海软黏土深基坑在机开挖离心模型试验研究[J]. 土木工程学报, 2023, 56(8): 131-139. MA Xianfeng, CAO Mingyang. Centrifuge model tests on excavation in Shanghai soft soil layers using in-flight excavation tools[J]. China Civil Engineering Journal, 2023, 56(8): 131-139. (in Chinese)
[4] 陈仁朋, 刘书伦, 孟凡衍, 等. 软黏土地层基坑开挖对旁侧隧道影响离心模型试验研究[J]. 岩土工程学报, 2020, 42(6): 1132-1138. doi: 10.11779/CJGE202006018 CHEN Renpeng, LIU Shulun, MENG Fanyan, et al. Centrifuge modeling of excavation effects on a nearby tunnel in soft clay[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(6): 1132-1138. (in Chinese) doi: 10.11779/CJGE202006018
[5] 郭海庆, 陶善之, 张泉. 基于离心机模型试验的超深基坑受力变形研究[J]. 地下空间与工程学报, 2020, 16(1): 177-186. GUO Haiqing, TAO Shanzhi, ZHANG Quan. Research on stress and deformation of super deep foundation pits based on the centrifugal model test[J]. Chinese Journal of Underground Space and Engineering, 2020, 16(1): 177-186. (in Chinese)
[6] 冯志, 沈正龙, 孟庆亮, 等. 泵站深基坑支护结构的离心模型试验研究[J]. 水资源与水工程学报, 2019, 30(4): 183-188. FENG Zhi, SHEN Zhenglong, MENG Qingliang, et al. Centrifugal model test of supporting structure of pump station deep foundation pit[J]. Journal of Water Resources and Water Engineering, 2019, 30(4): 183-188. (in Chinese)
[7] 王国辉, 陈文化, 聂庆科, 等. 深厚淤泥质土中基坑开挖对基桩影响的离心模型试验研究[J]. 岩土力学, 2020, 41(2): 399-407. WANG Guohui, CHEN Wenhua, NIE Qingke, et al. Impacts of pit excavation on foundation piles in deep silty soil by centrifugal model tests[J]. Rock and Soil Mechanics, 2020, 41(2): 399-407. (in Chinese)
[8] 鲁泰山, 何欢, 刘松玉, 等. 基于SCPTU测试的软土基坑开挖数值模拟研究[J]. 土木工程学报, 2023, 56(增刊2): 141-148. LU Taishan, HE Huan, LIU Songyu, et al. Study on numerical simulation of soft soil foundation pit excavation based on SCPTU test[J]. China Civil Engineering Journal, 2023, 56(S2): 141-148. (in Chinese)
[9] 张鹏, 邓智平, 王磊, 等. 基于Midas GTS的某近海深基坑开挖三维有限元数值模拟分析[J]. 甘肃科学学报, 2024, 36(1): 125-129. ZHANG Peng, DENG Zhiping, WANG Lei, et al. 3D finite element numerical simulation analysis of deep foundation pit excavation based on Midas GTS[J]. Journal of Gansu Sciences, 2024, 36(1): 125-129. (in Chinese)
[10] 包承纲, 饶锡保. 土工离心模型的试验原理[J]. 长江科学院院报, 1998, 15(2): 1-3, 7. BAO Chenggang, RAO Xibao. Test principle of geotechnical centrifugal model[J]. Journal of Changjiang River Scientific Research Institute, 1998, 15(2): 1-3, 7. (in Chinese)