1. 引言
双线隧道是区间地铁列车通行常采用的方案,伴随着我国城市地铁建设和运营历程的日益增长,地铁列车荷载作用下双线隧道的稳定问题日益突出 [1] [2]。隧道稳定问题与列车荷载作用下隧道周围地层的动力响应特性直接相关,近年来许多专家学者进行了相关的现场实测、室内试验、数值和理论分析研究 [3] [4]。包小华等 [4] 利用水–土耦合数值模型,分析了苏州地铁一号线隧道周围软土地层的动力学响应,获得了土体有效应力、超孔隙水压以及地层沉降等的变化特征。马学勇等 [5] 分析了南京地铁10号线单圆双线隧道下卧粉细砂层的动力响应,结果表明列车荷载激发的超静孔隙水压约为静水压力的1%。Tang等 [6] 根据现场实测和数值计算结果分析了上海地区地铁列车荷载作用下软土地层沉降与孔隙水压发展规律。不难发现,目前关于双线隧道的分析成果并不多,双线隧道列车荷载产生的相互影响将与单线隧道结果有显著区别 [7] [8],因此,本文通过简化建立了平行双线隧道计算模型,模拟交会相向行驶两辆列车荷载作用,对双线隧道下卧土体动力响应规律进行了总结,并与循环荷载作用下单线隧道动力响应计算结果进行了对比分析。
2. 数值建模与计算工况
2.1. 几何模型
如图1所示,采用平面应变有限元模型,隧道中心位于同一水平,模拟隧道外径为3 m,衬砌厚度为0.35 m,隧道埋深10 m,左右边界与隧道距离30 m,下边界与隧道底距离30 m。双线隧道中心间距分别考虑9 m、10 m和11 m三种工况,即双线隧道净间距分别为3 m、4 m和5 m。采用三角形单元对土体和隧道衬砌进行离散,在隧道周边以及下部的土体进行局部细化和加密。
![](//html.hanspub.org/file/14-2750969x10_hanspub.png)
Figure 1. Sketch of the numerical simulation model
图1. 数值模型示意图
2.2. 材料特性与边界条件
将土层简化为均质土层,采用强化岩土模型对粉质黏土层进行模拟。隧道衬砌结构使用弹性模型进行模拟。材料参数详见表1。
数值模型两侧采用水平约束边界条件(ux = 0),底部边界采用全约束边界条件(ux = uy = 0)。模型的两侧及底部设置标准吸收边界以避免杂波反射。在对应边界上被阻尼器吸收的法向和剪应力分量为:
(1)
式中,ρ是材料的密度。Vp 和Vs分别是压缩波波速和剪切波波速。C1和C2是松弛系数,根据经验取C1 = 1和C2 = 0.25可以使波在边界上得到合理的吸收。
2.3. 循环荷载
数值计算中,地铁列车循环荷载采用简化处理方法,即在隧道底部施加单点单周期荷载,采用正弦波形式:
(2)
式中,A为幅值,f为频率,t为时间。根据文献研究成果 [9],本文计算中取幅值A = 41 kN,频率f = 3.5 Hz,循环荷载作用时间为6 s,循环荷载停止后继续计算2 s,如图2所示。
3. 计算结果与分析
为分析双线隧道下卧土体的动力响应规律,并且方便后面与单线隧道进行比较,分别选择T2隧道中线下、双线孔隧道中线下以及T2隧道侧距隧道底一定距离的水平位置的若干点进行结果分析,具体分析点位置示意如图3所示。由于循环荷载对土体竖向动力响应的影响远大于水平方向,在以下分析中仅对竖向位移、竖向加速度和竖向速度的衰减变化规律进行分析。
![](//html.hanspub.org/file/14-2750969x15_hanspub.png)
Figure 3. Location of result analysis points
图3. 结果分析点位置示意图
3.1. 土体动力响应竖向衰减规律
以隧道T2正下方分析点A、B、C、D为例进行说明,分析点距离T2隧道底距离分别为5 m、10 m、15 m、20 m。图4和图5为分析点竖向加速度和速度时程曲线,可以看出,同一周期内加速度峰值正负大小基本相等,而速度峰值的负向绝对值稍大于正向峰值。土体中竖向加速度和速度峰值均随距循环荷载作用点垂直距离的增加而呈衰减趋势(见图6)。整体而言,土体中动力响应在6.9 s左右停止,即荷载停止施加0.9 s后土体振动基本停止。
![](//html.hanspub.org/file/14-2750969x16_hanspub.png)
Figure 4. Vertical acceleration beneath T2 tunnel
图4. T2隧道下方测点竖向加速度时程图
![](//html.hanspub.org/file/14-2750969x17_hanspub.png)
Figure 5. Vertical velocity beneath T2 tunnel
图5. T2隧道下方测点竖向速度时程图
![](//html.hanspub.org/file/14-2750969x18_hanspub.png)
Figure 6. Attenuation curves of vertical acceleration and velocity vs vertical distance
图6. 竖向加速度和速度随竖向距离衰减曲线
双线隧道中线位置土体由于受到两个隧道循环荷载作用的叠加影响,其动力响应值大于任一隧道正下方垂直测点的动力响应值(如图7所示)。分析点A和E与隧道下方垂直距离均为5 m,分析点A位于隧道T2正下方,而分析点E位于双线隧道中线上。对应分析点A处土体竖向加速度和速度峰值分别为0.62 m/s2和0.032 m/s,而分析点E处土体竖向加速度和速度峰值分别为0.72 m/s2和0.037 m/s。
![](//html.hanspub.org/file/14-2750969x19_hanspub.png)
Figure 7. Attenuation curves of vertical acceleration and velocity vs horizontal distance
图7. 竖向加速度和速度随水平向距离衰减曲线
3.2. 土体动力响应水平向衰减规律
与竖向相似,土体动力响应随距离荷载作用点水平距离的增加也呈衰减趋势。图7为以同一水平上分析点E、A、I、J位置处的竖向加速度和速度衰减曲线。由图中可知,循环荷载作用下随着距离的增加土体动力响应衰减越快,且水平向衰减幅度大于竖向的衰减幅度(对比图6)。
3.3. 隧道净距影响
由以上分析可知,分析点E受双线隧道循环荷载叠加作用最显著。以下以分析点E为例,对隧道净间距L的影响进行分析。图8给出了分析点E在不同隧道净间距条件下的竖向加速度时程曲线。结果显示,随着双线隧道净间距的增加,土体竖向加速度峰值有所降低,但减小有限。对比不同净间距条件下的竖向位移云图可见(图9至图10),土体竖向位移最大值均发生在隧道上方,且位移影响区域随着双线隧道净间距的增大而有所扩大。
3.4. 双线隧道与单线隧道对比
选择T2隧道正下方分析点A、B和C与单线隧道对应位置处的动力响应进行比较 [10]。图11和图12分别为竖向加速度和竖向速度峰值随深度的变化曲线。可以看出,双线隧道中振动荷载的叠加作用影响显著,且双线隧道中土体动力响应随距离的衰减要慢。模拟条件下,双线隧道竖向加速度峰值与单线隧道竖向加速度峰值之比介于1.55~1.87,两者竖向速度峰值之比介于1.68~1.70。
4. 结论
利用数值方法对双线地铁隧道下卧土体的动力响应特性进行了分析,可得到如下结论:
(1) 循环荷载作用过程中双线隧道下卧土体的竖向加速度等动力响应幅值基本稳定,波动不大,荷载停止施加后土体自振时间约为0.9 s (近似为单线隧道下卧土体自振时间的2倍)。
![](//html.hanspub.org/file/14-2750969x20_hanspub.png)
Figure 8. Vertical acceleration versus time data of point E with different net distance
图8. 不同净间距下E测点竖向加速度时程图
![](//html.hanspub.org/file/14-2750969x21_hanspub.png)
Figure 9. Contour of vertical displacement (L = 3 m)
图9. 竖向位移云图(L = 3 m)
![](//html.hanspub.org/file/14-2750969x22_hanspub.png)
Figure 10. Contour of vertical displacement (L = 5 m)
图10. 竖向位移云图(L = 5 m)
![](//html.hanspub.org/file/14-2750969x23_hanspub.png)
Figure 11. Comparison of vertical acceleration of single and double-track tunnel
图11. 单双线隧道竖向加速度对比
![](//html.hanspub.org/file/14-2750969x24_hanspub.png)
Figure 12. Comparison of vertical velocity of single and double-track tunnel
图12. 单双线隧道竖向速度对比
(2) 循环荷载作用下双线隧道下卧土体的动力响应随距离增加而逐渐衰减,且沿深度方向衰减速度小于沿水平向的衰减速度(深度方向衰减速率约为水平向的1/2),表明循环荷载对双线隧道下卧土层叠加效应主要体现在深度方向。
(3) 对比双线与单线隧道下卧土体的动力响应,双线隧道中任一隧道正下方分析点的动力响应值大于单线隧道对应分析点的值,其中竖向加速度峰值之比介于1.55~1.87,竖向速度峰值之比介于1.68~1.70。
基金项目
国家自然科学基金(51874286);江苏省自然科学基金项目(BK20140203)。
参考文献
NOTES
*通讯作者。