Shaking table tests on seismic liquefaction characteristics of soil surrounding tunnels in saturated loess stratum
-
摘要: 以兰合铁路某在建隧道为工程背景,对其黄土段开展了1∶20大比例尺饱和黄土围岩地震液化的系列振动台试验,揭示了土-结构相互作用对隧道衬砌周围土体液化势的影响规律。研究结果表明:土-结构相互作用显著提高了隧道衬砌周围土体的液化势,当遭遇地震荷载时,围岩会先于同一水平地层土体而产生液化;激励地震加速度的较小时(= 1.0 m/s2),同一水平高度靠近衬砌结构土体的动剪应力峰值、动孔隙水压力峰值、和加速度峰值都要比远离衬砌结构的大,具有明显的放大效应;随着激励地震加速度的增大至接近围岩液化阶段时,动孔隙水压力的放大效应明显,而剪应力和加速度的放大效应减小或者消失。同时上述结果中围岩土体加速度的变化特征,揭示了前期相关数值分析研究存在的不足。Abstract: Based on the previous numerical and theoretical analyses, a series of shaking table tests with scale of 1∶20 for seismic liquefaction analysis of saturated surrounding loess are carried out to simulate the loess section of a tunnel under construction on the Lanzhou-Hezuo Railway in China. The results show that the soil–structure interaction significantly increases the liquefaction potential of the soil near the tunnel linings, and the surrounding soil will produce liquefaction before the soil of the same level when subjected to seismic loads. When the peak acceleration of seismic loads is minor (i. e., = 1.0 m/s2), the peaks of dynamic shear stress , dynamic pore water pressure , and acceleration near the lining structure are larger than those far away from the lining at the same level height, and have obvious amplification effects. With the increase of to near the liquefaction stage of surrounding soil, the amplifying effects of dynamic pore water pressure are obvious, while the amplifying effects of shear stress and acceleration decreases or disappear. In particular, the variation characteristics of the acceleration of the surrounding soil reveal the deficiency of the previous numerical simulations.
-
Keywords:
- loess tunnel /
- liquefaction potential /
- shaking table test /
- dynamic shear stress
-
0. 引言
随着城市化进程不断加速,城市地下空间的开发越来越受到人们的关注。大规模的新建地铁盾构隧道在掘进过程中,不可避免地出现穿越既有建筑桩基的工况[1]。隧道开挖引起的地层应力释放,造成既有桩基周围地层的应力场发生改变,从而导致桩基产生附加变形和附加内力,并造成桩基承载力发生变化,而过度的承载力损失将危及上部结构安全[2-3]。因此,准确评价隧道近接施工对既有摩擦桩基的影响是当前城市地下空间开发所面临的较为紧迫的问题。
当前,针对隧道施工引起的邻近桩基响应问题,学者们基于现场实测、模型试验[4-5]、数值模拟和理论解析等方法已取得了较为丰硕的成果。如李雪等[6]针对饱和砂土地区盾构隧道超近接高铁桥墩摩擦桩的工程问题,对既有高铁桥墩桩基变形进行了现场测试;Loganathan等[7]和孙庆等[8]采用离心机试验研究了隧道开挖对邻近桩基的瞬态和长期影响,指出隧道开挖影响下邻近桩基最大横向变形和附加弯矩均发生在隧道轴线附近。
隧道施工引起的邻近桩基响应问题,实质为隧道−土体−桩基之间相互作用。学者们已构建了十分丰富的解析计算模型。如ZHANG等[9]基于Pasternak地基模型及桩侧土体影响,研究了被动桩在盾构掘进下水平向变形响应;程康等[10]等进一步考虑桩基的横向剪切性能,将桩简化成Pasternak地基上的Timoshenko梁,进一步研究了盾构掘进引起的桩基水平变形;冯国辉等[11-13]分别将既有桩基视为Winkler地基和Pasternak地基中的Euler-Bernoulli梁,分析了盾构掘进下邻近桩基的横向附加响应。
本文在既有研究的基础上,结合隧道开挖前后地层应力场的变化,建立盾构近接施工下邻近摩擦桩基承载力解析计算模型,分析了隧道施工对既有摩擦桩基承载特性的影响机制,并进一步基于三维有限元分析论证了隧道洞周地层加固对既有桩基变形控制的有效性。
1. 盾构近接施工扰动下摩擦桩基承载特性分析
1.1 基本思路
盾构施工对邻近摩擦桩基承载特性的影响,实质在于隧道施工会引起周围土体应力场重分布,而导致摩擦桩侧法向正应力和切向剪应力发生变化,继而影响桩侧摩阻力。因此对于盾构施工扰动下邻近摩擦桩基承载力分析的基本思路:
(1)隧道施工前,地层内任一点的应力状态由围岩自重生成竖向应力σv0、水平向应力σh0,将该值作为桩侧地层应力状态初始值,可求得初始桩侧阻力Qs0。
(2)计算隧道施工后地层内任一点的竖向应力σv、水平应力σh和剪应力τxy,根据此时的应力状态计算隧道施工后的桩侧阻力Qs。
(3)结合步骤(1),(2),求得摩擦桩基承载力变化值ΔQs=Qs0−Qs。
1.2 解析模型的建立及基本假定
为获得盾构近接施工扰动下临近摩擦桩基的承载特性,建立如图 1所示的隧道-地层-桩基相互作用模型。图 1中x,y为弹性半空间任意一点的直角坐标,r,θ为相应的极坐标;h为隧道轴线埋深;σr,σθ,τrθ分别为任意一点土体径向应力、环向应力和剪应力;l为桩基长度;x0为桩基轴线与隧道轴线的水平间距;r0为隧道开挖半径。
在理论计算时作如下基本假定:
(1)为避免地层应力场计算时的奇异,隧道覆径比大于1.5。
(2)忽略土体与衬砌之间的摩擦力,即接触面处的剪应力满足:
τ|r=r0=0。 (1) (3)采用间隙参数w表征盾构开挖引起的等效地层损失,计算方法为
w=Vlπr022πr0=12Vlr0。 (2) 式中:Vl为地层损失率。
1.3 初始桩侧阻力Qs0
任在原岩应力状态下,桩侧任意一点的应力状态为
σv0=γ(h−y),σh0=kσv0,τxy0=0。} (3) 式中:γ为土体有效重度;k为侧压力系数。
根据摩擦桩基承载理论,此时桩基承载力为
Qs0=π d∫lfs0(y)dy=π d∫l(τxy0+ksσh0)dy。 (4) 式中:d为桩基直径;fs0为初始状态下的桩侧摩阻力分布函数;ks为桩侧界面摩擦系数,ks=tanδφ;δφ为界面摩擦角。
1.4 开挖扰动后桩侧阻力Qs
以干砂地层为例,BOBET基于应力函数法给出了此类地层中隧道开挖后任意一点的应力,即
σr=a0r2−12γh(1+k)+[14γr(k+3)−2c′1r−3+12c1r−1]sinθ+[12γh(1−k)−6a′2r−4−2b′2r−2]⋅cos2θ+[−14γr(1−k)−12c′3r−5−6d′3r−3]sin3θ, (5a) σθ=−a0r2−12γh(1+k)+[14γr(3k+1)+2c′1r−3−12c1r−1]sinθ+[−12γh(1−k)+6a′2r−4+2b′2r−2]⋅cos2θ+[14γr(1−k)+12c′3r−5+6d′3r−3]sin3θ, (5b) τrθ=[14γr(1−k)+2c′1r−3]cosθ−[12γh(1−k)+6a′2r−4+2b′2r−2]sin2θ−[14γr(1−k)−12c′3r−5−6d′3r−3]cos3θ。 (5c) 式中:a0,a′2,b′2,c1,c′1,c′3,d′1,d′3为常数,可根据问题的边界条件确定,其他符号意义同前。
进一步,可通过坐标变换将地层中一点应力状态(式(5))由极坐标系转换到直角坐标系,即
σx=σr+σθ2+σr−σθ2cos2θ−τrθsin2θ, (6a) τxy=σr−σθ2sin2θ+τrθcos2θ。 (6b) 将式(5)代入式(6),结合摩擦桩基承载理论,得到隧道开挖扰动后的桩基承载力,即
Qs=πd∫lfs(y)dy=πd∫l(τxy+ksσx)dy。 (7) 式中:fs为隧道开挖后桩侧摩阻力分布函数;σ′x,τ′xy分别表示直角坐标系下任意一点的有效水平应力和剪应力。
1.5 隧道开挖对摩擦桩基承载力的影响分析
以某地铁盾构隧道连续侧穿老旧建筑群施工为工程背景,基于上述隧道开挖前后摩擦桩基承载力理论分析,进行隧道开挖对邻近摩擦桩基承载特性的影响分析。选取典型工况:既有桩基及开挖隧道位于砂土地层中,地层重度约为18 k N/m3,侧压力系数约为0.32,根据地勘资料,计算时地层弹性模量和泊松比分别取30 MPa,0.3;盾构隧道开挖半径为3.14 m,轴线埋深为15 m,隧道衬砌采用C50钢筋混凝土管片,管片厚度为0.3 m,弹性模量和泊松比分别为34500 MPa,0.2;既有建筑物桩基直径为0.5 m。为反映隧道施工对摩擦桩基承载能力的影响,定义摩擦桩基承载力影响因子RQ:
RQ=QsQs0。 (8) RQ>1表示承载力提高,RQ<1表示承载力损失,当RQ=1时表示承载力无变化。
如图 2所示为桩隧横向间距为1.5r0,地层损失率为0.8%时,摩擦桩基承载力影响因子随桩长的变化曲线。由图可见:承载力影响因子变化曲线随桩长的增大而呈漏斗状分布,当l/h<1.0时,承载力影响因子随桩长的增大而迅速减小,而当l/h>1.0时,承载力影响因子随桩长的增大而又逐渐增大,并逐渐趋近于1.0;桩隧相对埋深越趋近1.0,承载力影响因子变化梯度越大,承载力损失越多;当桩基端部与隧道轴线相平时,即l/h=1.0时,其承载力损失最多,在本算例中此时承载力影响因子RQ=0.17,桩基承载力损失82.9%。
如图 3所示为当桩隧横向间距为3.0 r0,桩隧相对埋深l/h = 1.0时,摩擦桩基承载力影响因子随地层损失率的变化曲线。由图 3可见:摩擦桩基承载力损失因子与地层损失率呈线性相关关系,随着地层损失率的增大,摩擦桩基承载力影响因子随之线性减小;在本算例中,当地层损失率仅为0.4%时,承载力影响因子为0.98,而当地层损失率增大至1.8%时,摩擦桩基承载力影响因子减小至0.15,变化84.7%。
2. 盾构侧穿摩擦桩基建筑物控制技术
盾构隧道施工对周围岩土体造成扰动,应力平衡状态改变造成围岩应力释放,导致临近桩基产生沉降变形。综合前述计算分析,从地层加固措施入手减小盾构施工对邻近桩基的影响,通过注浆加固地层的方式减小盾构施工扰动在地层中的传播。
2.1 计算模型及计算参数
所建立的盾构-围岩-建筑物的三维数值模型如图 4所示,模型尺寸为长60 m (沿隧道开挖方向),宽120 m (垂直隧道开挖方向),高50 m (沿地层深度方向)。模型侧面约束水平位移,底部为固定边界,上表面为自由边界,共划分56972个单元,41673个节点。数值模型中盾构隧道管片外径6.0 m,内径5.4 m,厚度300 mm,宽度1.2 m,采用C50混凝土,考虑衬砌接头对管片刚度的削弱效应,管片刚度需折减15%。
模型中地层采用摩尔库仑模型,盾壳、隧道衬砌、同步注浆区(即后文所述的等代层)、袖阀管注浆加固区、建筑物梁板柱、桩基础均采用线弹性模型。地层、等代层、袖阀管注浆加固区采用实体单元模拟,盾壳、隧道衬砌、建筑楼板采用板单元模拟,建筑框架梁、立柱、建筑桩基采用梁单元模拟。参数取值见表 1。
表 1 结构参数取值Table 1. Values of structural parameters名称 密度
γ/(kN·m-3)弹性模量
E/GPa泊松比 盾壳 78.5 210 0.30 隧道衬砌 25.0 29.3 0.20 建筑框架
(梁、板、柱)25.0 20.4 0.20 桩基础 25.0 24.0 0.20 等代层 20.0 0.06 0.25 袖阀管注浆加固区 22.0 0.25 0.25 2.2 盾构隧道施工过程模拟
(1)同步注浆
在盾构机施工过程中,在盾尾管片脱出后,由于盾构机刀盘的超挖间隙、装配管片所需的操作间隙、盾构机壳的厚度3个因素将会在盾构管片和土体间产生盾尾间隙。为了减小施工扰动,控制地层变形,在盾构掘进过程中将会采取同步注浆对盾尾间隙进行填充。但是实际工程中,由于地层向盾尾间隙道的移动程度、注浆体注浆后的分布情况、隧道围岩壁面受扰动的范围和程度等情况相对复杂,因此,在建模过程中通常将注浆体简化为均质的、弹性的、等厚的等代层处理,等代层的厚度δ可为
δ=μΔ。 (9) 式中:Δ为盾尾间隙的理论计算值;μ为折减系数,取值范围为0.7~2.0,对于极软土层取上限值,对于硬土层取下限值。
根据盾构穿越地层条件及施工扰动范围及程度取等代层厚度为100 mm,弹性模量为60 MPa,泊松比为0.25,采用线弹性体模拟,并根据实际施工情况取注浆压力为300 kPa。
(2)土仓压力
对于土压平衡盾构机而言,掌子面的稳定一般依靠于密封土仓的土仓压力平衡掌子面前方的水土压力,并且考虑到盾构刀盘切削土体使土体剥落而导致土仓压力不稳定的情况,为了防止掌子面坍塌,根据施工经验,一般会预留10~20 kPa的预备压力,所以土仓压力一般按下式确定:
P0=Pc+Pw+p。 (10) 式中:P0为土仓压力;Pc为土压力;Pw为水压力;P为预备压力。其中Pc=∑K0γizi,Pc=γwh,K0为土体静止侧压力系数,γi为第i层土的重度,zi为第i层土的厚度,γw为水的重度,h为地下水位都盾构刀盘的高度。
通过计算取土仓压力为200 kPa均布施加于开挖面土体上。
(3)盾壳与地层摩擦力
盾构掘进过程,盾构机和地层发生相对运动而产生摩擦力,会使得盾构机周围土体由于剪切作用而向前运动,沿伸到地表则表现为地表隆起或者沉降。盾壳和地层间的摩擦力计算公式为
f=μs∑γizi。 (11) 式中:μs为盾构机与地层的摩擦系数,通过计算得到的摩擦力为153 kPa。
3. 计算结果分析
3.1 地层位移
在整个开挖过程中,注浆加固区附近的土体变形明显变小,当双线隧道开挖完成后,对于未进行注浆加固的工况来说,沉降发生区域在左线隧道和右线隧道间呈长条带状分布,而注浆加固后,沉降发生区域沿开挖方向呈现“两头大,中间小”分布,凸向隧道左线。选取隧道轴线纵向30 m处的横断面作为分析断面,对加固地层前后的地表沉降进行分析,计算结果如图 5所示。
由图 5可知,地表沉降曲线呈“U”字形,符合“Peck”曲线分布,并与实测数据有较好的拟合度。当右线隧道开挖完成时,地表沉降最大值出现在右线隧道正上方,当左线隧道开挖完成后,地层受二次开挖扰动,沉降值变大,沉降槽宽度加宽,最低点往左移动。对于距起始开挖面30 m处的横断面而言,处于加固区,受右线隧道开挖影响减小,沉降最大值出现在左线隧道正上方附近,注浆加固前后沉降明显变小,最大值分别为14.8,10.7 mm,降低27.8%。
3.2 桩基变形
选取距离起始开挖面最近的一排3根桩进行分析,距离隧道由近到远分别为桩基P1,P2,P3。由图 6可知,双线隧道开挖完后,由于土体的卸荷作用,桩基发生挠曲变形,桩基顶部向隧道方向移动,桩基底部远离隧道移动。离隧道越远,桩基水平位移越小,由于上部结构的约束作用,各单桩桩顶位移基本相同。在注浆加固后,桩基水平位移明显减小,桩顶最大水平位移由2.6 mm减少至1.4 mm,降低46%。
如图 7所示为加固地层和未加固地层工况下左右线贯通后桩基竖向位移沿桩身埋深分布曲线。由图 7可见:桩身沉降沿桩身基本相同,离隧道越远,沉降越小,桩基P1,P2沉降受隧道开挖影响较大,在注浆加固后,桩基沉降由5.7 mm减少至2.6 mm,降低54%。
4. 结论
针对盾构近接施工对邻近摩擦桩基的影响问题,分别建立了隧道开挖扰动下摩擦桩基承载力解析计算模型和三维数值分析模型,探明了隧道开挖对既有摩擦桩基承载力的影响机制,并论证了地层加固对桩基变形控制的有效性,得到3点结论。
(1)基于摩擦桩基承载力理论和隧道开挖前后地层应力场的分析,建立了隧道开挖影响下摩擦桩基承载力计算方法,并定义摩擦桩基承载力影响因子分析了不同工况下隧道施工对桩基承载力的影响。
(2)隧道施工对不同长度的摩擦桩基影响不同,当桩端埋深与隧道轴线埋深一致时,桩基承载力损失最为严重,而当桩长远大于隧道埋深时其承载力损失较小;隧道施工时的地层损失率越大,桩基承载损失越严重,二者呈正相关。
(3)建立了盾构侧穿浅层摩擦桩基建筑物施工的三维有限元模型,分析了盾构侧穿桩基施工时的桩基沉降变形规律。计算结果表明在采取洞周地层加固措施后,隧道围岩注浆加固后土体变形明显变小,可有效减少临近桩基水平位移和桩身的沉降。
-
表 1 隧道现场黄土物理特性
Table 1 Physical properties of in-site loess
孔隙比e 干密度ρd/(g·cm-3) 含水率w/% 黏粒d ≤ 0.005 粉粒0.005<d≤0.075 (mm) 砂粒0.075<d ≤ 0.25 (mm) 0.86 1.45 17.2% 16.8% 67.7% 14.5% -
[1] WANG Lanmin. Loess Dynamics[M]. Beijing: Seismological Press, 2003. (in Chinese) [2] 王兰民. 黄土地层大规模地震液化滑移的机理与风险评估[J]. 岩土工程学报, 2020, 42(1): 1-19. doi: 10.11779/CJGE202001001 WANG Lanmin. Mechanism and risk evaluation of sliding flow triggered by liquefaction of loess deposit during earthquakes[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(1): 1-19. (in Chinese) doi: 10.11779/CJGE202001001
[3] 李德武, 马为功. 二次衬砌施作时机的弹粘塑性有限元分析[J]. 现代隧道技术, 2012, 49(4): 6-9. https://www.cnki.com.cn/Article/CJFDTOTAL-XDSD201204004.htm LI Dewu, MA Weigong. Elasto-viscoplastic finite element analysis of the installation time of a secondary lining[J]. Modern Tunnelling Technology, 2012, 49(4): 6-9. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-XDSD201204004.htm
[4] GUAN Baoshu. Key Points of Tunnel Engineering Design[M]. Beijing: China Communications Press, 2003. (in Chinese) [5] 许新桩. 陕北黄土地区铁路隧道基底病害机理分析及治理措施[J]. 公路交通科技(应用技术版), 2015, 11(9): 17-19. https://www.cnki.com.cn/Article/CJFDTOTAL-GLJJ201509008.htm XU Xinzhuang. Mechanism analysis and treatment measures of railway tunnel foundation disease in loess area of northern Shaanxi[J]. Journal of Highway and Transportation Research and Development, 2015, 11(9): 17-19. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GLJJ201509008.htm
[6] 来弘鹏, 谭智鹏, 孙玉坤, 等. 富水黄土隧道施工过程围岩水分迁移规律研究[J]. 中国公路学报, 2023, 36(1): 150-161. doi: 10.3969/j.issn.1001-7372.2023.01.013 LAI Hongpeng, TAN Zhipeng, SUN Yukun, et al. Study on law of water migration in surrounding rock during construction of water-rich loess tunnel[J]. China Journal of Highway and Transport, 2023, 36(1): 150-161. (in Chinese) doi: 10.3969/j.issn.1001-7372.2023.01.013
[7] MA W G, WANG L M, WANG P, et al. The variation characteristics of dynamic shear stress reduction coefficient rd in loess tunnel stratum[J]. Soil Dynamics and Earthquake Engineering, 2023, 174: 108192. doi: 10.1016/j.soildyn.2023.108192
[8] SEED H B, IDRISS I M. Simplified procedure for evaluating soil liquefaction potential[J]. Journal of the Soil Mechanics and Foundations Division, 1971, 97(9): 1249-1273. doi: 10.1061/JSFEAQ.0001662
[9] ZHANG Kexu. Geotechnical Earthquake Engineering and Engineering Vibration[M]. Beijing: Science Press, 2016. (in Chinese) [10] DING X M, ZHANG Y L, WU Q, et al. Shaking table tests on the seismic responses of underground structures in coral sand[J]. Tunnelling and Underground Space Technology, 2021, 109: 103775.