摘要
为保证软土地层中地铁线路建造及运营的稳定性,可通过注浆加固提高结构周边土体强度,而根据加固程度不同,土体介质的非连续性将逐渐转变,对应动力传递特性随之变化,因此行车状态下加固区及周边地层的动力响应分布模式有待进一步探索. 基于解析理论建立含加固区的地铁线路-地层动力耦合模型,将加固土及天然土考虑为具有不同物理力学属性的饱和多孔介质,利用Biot理论分别推导运动方程,地铁结构视为均匀厚度的无限长圆柱壳体,采用无扭矩条件下薄壁圆柱壳理论进行描述. 根据动力系统各交界面间的位移及应力连续条件,将以上各部位振动方程构建形成耦合动力方程组,计算了地铁周边含加固区地层的振动情况,并进一步对比了加固区刚度及黏滞阻尼对地层动力响应分布及传递的影响. 结果表明:稳态作用下,地铁周边地层切向加速度最大值在结构斜下方两侧30°处呈蝉翼状分布,径向加速度最大值出现在荷载作用方向土体内部;由于输入波形的简谐特性,加速度在地层中呈波动式衰减,但在加固区范围内切向加速度出现明显的放大现象,振动进入天然土后,受地层阻尼影响,动力响应显著衰减;提高加固区刚度有利于降低对切向加速度的传导,而径向加速度变化不大,黏滞阻尼减小时,所传递动力响应有所提升.
随着居民对市内出行需求的逐步提高,地铁建设正在快速发展. 为便于隧道开挖,城市地铁线路大多建设在较松散地层,对于土质较差地区,尤其是临近重要结构物或易产生较大扰动变形的富水软土地层,可在开挖前通过注浆方式在线路周围形成一定范围的加固区,以保证施工稳定性. 进入运营阶段后,列车荷载引起的环境振动问题逐渐受到关注,加固区的存在必将影响地铁线路周边近场地层的振动效应,有必要进一步开展含加固区地层动力响应传递特性的研究.
已有针对地铁线路周边地层动力响应的研究中,马蒙
由以上研究成果可知,数值方法在地铁线路环境振动研究中发挥了重要作用,随着模型的逐步完善,可以从多方面反映埋置列车荷载引起的地层动力响应. 但各模型尚未考虑加固作用对线路近场动力特性的影响,为了探索影响机理,需借助解析方法加以深入分析. 胡安峰
上述各类解析方法同样未涉及对地层局部加固线路周边土体动力响应的求解. 而对于地铁预加固处理的描述,大多集中于加固机理探
1 计算模型
考虑到地铁区间线路为沿轴线延伸的长条状结构,列车运行的动力响应在各断面较为相似,因此选取代表性断面建立含加固区的地铁线路-地层动力耦合模型,如

图1 含加固区动力计算模型
Fig. 1 Dynamic calculation model with reinforcement area
模型中将隧道结构简化为埋置于均质饱和土体中的无限长空心圆柱壳体,隧道内径为R1,外径为R2,隧道外部存在厚度为h1的加固层. 将列车荷载模拟为作用在隧道底拱处的径向集中简谐荷载,考虑隧道衬砌、注浆层及外部饱和土体三者的动力协调关系展开求解. 模型假定如下:
1)忽略隧道及周围土体沿线路轴向的变形及应力,将荷载作用下模型的动力响应看作平面应变问题.
2)隧道埋深可不考虑上部地表界限.
3)模型各层间保持位移和应力协调条件.
2 控制方程及求解
2.1 饱和介质波动方程
由Biot理
(1) |
(2) |
式中:λ和μ为介质拉梅常数;α和M为Biot常数,即土颗粒和孔隙流体的压缩性参数;ρ=ns ρf+(1-ns)ρs,ns为土体孔隙率,ρs和ρf分别为土骨架和流体密度;m=ρf/ns;b为与流体动力黏滞系数ηd和土体动力渗透系数k相关的土体渗透参数;ui和wi分别为土骨架位移及孔隙水相对土骨架位移,两者利用标量势函数Φ和矢量势函数Ψ可表示为u│w=▽ϕ+▽×ψ,并在极坐标下展开为法向及切向位移分量,表达如下:
(3) |
简谐荷载(Fn
(4) |
(5) |
式中:为Laplace算子,在极坐标中.
从
(6) |
(7) |
式中:k1、k2、k3分别表示饱和土中快纵波、慢纵波和横波的复波数,即,,其中, ,,.
另外,由
(8) |
式中:;.
饱和多孔介质应力应变关系为:
(9) |
式中:δij为Kronecker delta符号;e=ui,i为土骨架体积应变;εij=(ui,j+uj,i)/2为土骨架应变;孔压pf=-αMe+Mζ,其中ζ=-wi,i.
通过对势函数的谐振表达,将
取加固层(R2≤r≤R2+h1)中n阶模态下土体位移、应力响应可表示为:
(10) |
(11) |
r=∞时,土体介质远离激振点,考虑无限远处位移、应力及孔压分量均衰减为0,即A1n=B1n=En=0. 因此加固层外部r>R2+h1处饱和土体在n阶模态下的位移、应力响应可表示为:
(12) |
(13) |
式中:[U3×3]、[T3×3]为[U3×6]、[T3×6]中去除A1n、B1n和En对应列的剩余矩阵.
2.2 衬砌壳体运动方程
基于无扭矩条件下薄壁圆柱壳理论,n阶模态下衬砌的运动方程可表示为:
(14) |
(15) |
式中:Ur、Uθ分别为隧道衬砌径向及环向位移;;,,,El、ρl、vl分别为衬砌的弹性模量、密度和泊松比;h为衬砌厚度;q、s分别为隧道衬砌所受径向与切向净应力,即衬砌内表面的径向和切向外荷载Q、S与衬砌所受加固层的径向和环向土压力σrr、σrθ的差值.
2.3 边界条件
在r=R2处,根据隧道衬砌与周围土体交界面上的位移连续条件,有:
(16) |
若隧道底拱θ处荷载为径向点荷载Q,则:
(17) |
式中:F荷载幅值;δ(·)为狄拉克函数.
不考虑隧道边界透水,则:
(18) |
在r =R2+h1处,根据加固层与饱和土体接触面位移、应力连续条件有:
(19) |
联合式(10)~
3 动力传递特性分析
以常规软土地层地铁隧道为例,计算参数取值如
参数 | 数值 |
---|---|
拉梅常数μ/Pa |
7.7×1 |
拉梅常数λ/Pa |
1.9×1 |
土颗粒压缩性参数α | 1.0 |
孔隙流体压缩性参数M |
5.0×1 |
流体密度ρf/(kg· |
1.0×1 |
土骨架密度ρs/(kg· |
2.0×1 |
饱和土体孔隙率ns | 0.4 |
土体渗透参数b/(N·s· |
1.0×1 |
衬砌弹性模量El/MPa |
3.0×1 |
衬砌密度ρl/(kg· |
2.5×1 |
衬砌泊松比νl | 0.2 |
隧道外径R2/m | 3.0 |
衬砌厚度h/cm | 25.0 |
加固区厚度h1/m | 3.0 |
作为对比,首先针对未加固及不同加固层厚度下的地层动力响应分布规律进行分析,以说明加固区对动力传递的影响. 围绕3 m加固厚度,计算单位竖向荷载作用下,土体未加固及加固范围为2、3、4 m时,地铁线路周边地层的竖向加速度响应如

(a) 线路侧向地层

(b) 线路竖向地层
图2 不同加固范围下地层振动加速度分布
Fig.2 Distribution of stratum vibration acceleration under different reinforcement ranges
为验证模型的可靠性,利用ABAQUS软件依照上述参数建立外形尺寸为100 m×80 m×60 m的数值模型,并对隧道周围3 m范围内土体进行加固,如

图3 含加固区动力有限元模型
Fig. 3 Dynamic finite element model with reinforcement area

图4 计算结果对比
Fig. 4 Comparison of calculation results
由
3.1 加固区刚度
隧道外围加固区可以看作土体与浆液的混合体,弹性模量与注入浆液量正相关. 考虑土体不同加固程度,弹性模量可在40~120 MPa间取值,不同加固区刚度对应地层动力响应如

(a) r=6 m处切向加速度

(b) r=6 m处径向加速度

(c) r=10 m切向加速度

(d) r=10 m处径向加速度
图5 加固区刚度对地层动力响应的影响
Fig.5 Influence of reinforcement area stiffness on stratum dynamic response
从

(a) θ=150°方向切向加速度

(b) θ=180°方向径向加速度
图6 加固区刚度影响下加速度幅值衰减曲线
Fig.6 Acceleration amplitude attenuation curves under the influence of reinforcement area stiffness
由
3.2 加固区黏滞阻尼
土体黏滞阻尼比β受到加固材料选取及加固程度的共同作用而产生变化,为实现单因素分析目的,保证其他参数不变的情况下,考虑β取0.1、0.05、0.025,分别计算各工况隧道周围地层动力响应,如

(a) r=6 m处切向加速度

(b) r=6 m处径向加速度

(c) r=10 m处切向加速度

(d) r=10 m处径向加速度
图7 加固区黏滞阻尼对地层动力响应的影响
Fig.7 Influence of viscous damping of reinforcement area on stratum dynamic response

(a) θ=150°处切向加速度

(b) θ=180°处径向加速度
图8 加固区黏滞阻尼影响下加速度幅值衰减曲线
Fig.8 Acceleration amplitude attenuation curves under the influence of viscous damping of reinforcement area
4 结 论
本文基于Biot饱和多孔介质理论建立土体运动方程,模拟加固及天然状态下地层动力传导特性,同时采用圆柱壳体来模拟地铁隧道结构,进而以多层介质交界面位移及应力协调关系建立含加固区的隧道-地层动力耦合解析模型. 利用上述模型探索了稳态作用下地铁周边地层加固参数对动力系统振动传递的影响,所得结论如下:
1) 竖向简谐荷载作用下,地铁周边地层切向加速度呈蝉翼状分布,最大值出现在结构中心竖轴两侧各30°斜下方土体中. 径向加速度最大值沿荷载作用方向出现在地铁下方地层. 在远离地铁结构过程中,地层振动加速度呈波动式衰减.
2) 地层加固对结构周边土体切向加速度有显著放大作用,最大值出现在加固区边缘位置. 径向加速度在加固区内部迅速下降,未体现放大趋势.
3) 加固区刚度越大,对切向动力响应的传导性能越弱,而径向加速度受加固程度的影响较小. 随着加固区阻尼减小,对应动力阻滞作用降低,地层响应幅值随之增大.
参考文献
马蒙,刘维宁,刘卫丰.列车引起环境振动预测方法与不确定性研究进展[J].交通运输工程学报,2020,20(3):1-16. [百度学术]
MA M, LIU W N,LIU W F. Research progresses of prediction method and uncertainty of train-induced environmental vibration[J]. Journal of Traffic and Transportation Engineering, 2020, 20(3):1-16.(in Chinese) [百度学术]
HOU C Y,CUI Z D,YUAN L.Accumulated deformation and microstructure of deep silty clay subjected to two freezing-thawing cycles under cyclic loading[J].Arabian Journal of Geosciences,2020, 13(12): 452. [百度学术]
ZHONG Z L,ZOU H,HU X X,et al.Experimental study on stiffness softening of soil-rock mixture backfill under metro train cyclic load[J].Advances in Materials Science and Engineering,2021, 2021(1): 3024490. [百度学术]
QU S, YANG J J, ZHU S Y, et al.Experimental study on ground vibration induced by double-line subway trains and road traffic[J]. Transportation Geotechnics, 2021,29:100564. [百度学术]
BASHIR S,CHOWDHARY A R,AKHTAR N.Prediction and attenuation of ground vibrations generated by moving trains[J]. Current Science, 2022, 124(2): 202-209. [百度学术]
WU L,ZHANG X D,KUANG F M,et al.Dynamic characteristics of adjacent tunnel and surroundings under the vibration of trains in the first tunnel with different buried depths[J]. Arabian Journal for Science and Engineering, 2021, 46(5): 5105-5119. [百度学术]
赵江涛, 牛晓凯, 苏洁, 等. 地铁致凸形地貌地表振动局部放大现象实测及机理分析[J]. 振动与冲击, 2022,41(4):287-293. [百度学术]
ZHAO J T,NIU X K,SU J,et al. In-situ measurement and mechanism analysis for local amplification phenomena of ground vibration in the protruding topography induced by subway[J].Journal of Vibration and Shock,2022,41(4):287-293.(in Chinese) [百度学术]
马龙祥, 刘维宁, 刘卫丰, 等. 地铁列车振动环境影响预测的薄片有限元-无限元耦合模型[J]. 岩石力学与工程学报,2016,35(10): 2131-2141. [百度学术]
MA L X,LIU W N,LIU W F,et al.Sliced finite element-infinite element coupling model for predicting environmental vibration induced by metro train[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(10): 2131-2141.(in Chinese) [百度学术]
雷华阳, 张磊, 许英刚, 等. 快速地铁列车荷载下软土地基沉降数值模拟研究[J]. 岩土工程学报, 2019, 41(增刊1): 45-48. [百度学术]
LEI H Y,ZHANG L,XU Y G,et al.Study on numerical simulation of soft soil foundation settlement under the load of express subway train[J]. Chinese Journal of Geotechnical Engineering,2019, 41(Sup.1): 45-48.(in Chinese) [百度学术]
王力东 ,朱志辉, 韩艳, 等. 地铁车致隧道与土体振动的高效时-频混合预测方法[J]. 振动工程学报, 2022, 35(2):359-368. [百度学术]
WANG L D,ZHU Z H,HAN Y,et al. Efficient time-frequency hybrid method for predicting the subway train induced vibrations of tunnel and ground[J]. Journal of Vibration Engineering, 2022,35(2): 359-368.(in Chinese) [百度学术]
GU Z F,WEI H L,LIU Z K,et al. Dynamic response mechanism of silt ground under vibration load[J]. Sustainability, 2022, 14(16): 10335. [百度学术]
胡安峰, 孙波, 谢康和. 下卧基岩饱和地基在移动荷载作用下的动力响应[J]. 振动与冲击, 2012, 31(4): 151-156. [百度学术]
HU A F, SUN B, XIE K H. Dynamic response of saturated subgrade with rock substratum subjected to moving loads[J].Journal of Vibration and Shock, 2012, 31(4): 151-156.(in Chinese) [百度学术]
ZHOU S H, HE C,GUO P J,et al. Modeling of vehicle-track-tunnel-soil system considering the dynamic interaction between twin tunnels in a poroelastic half-space[J]. International Journal of Geomechanics, 2020, 20(1): 04019144. [百度学术]
袁宗浩, 戴凯鑫, 潘晓东, 等. 黏弹性边界对圆形衬砌隧道动力响应的影响[J]. 浙江工业大学学报, 2020, 48(5): 532-540. [百度学术]
YUAN Z H,DAI K X,PAN X D,et al. Effect of viscoelastic boundary on dynamic response of circular tunnel[J]. Journal of Zhejiang University of Technology, 2020, 48(5): 532-540.(in Chinese) [百度学术]
YUAN Z H, CAI Y Q,SUN H L,et al. The influence of a neighboring tunnel on the critical velocity of a three-dimensional tunnel-soil system[J]. International Journal of Solids and Structures, 2021, 212: 23-45. [百度学术]
黄强, 冯青松, 黄宏伟,等.移动简谐荷载下浮置板轨道-隧道-地基模型二维振动响应研究[J]. 应用基础与工程科学学报, 2023, 31(3): 663-674. [百度学术]
HUANG Q, FENG Q S, HUANG H W, et al. 2D dynamic response of floating slab track-tunnel-ground model under moving harmonic load[J]. Journal of Basic Science and Engineering,2023, 31(3): 663-674.(in Chinese) [百度学术]
何超, 贾缘平, 周顺华. 隧道-地基-流体耦合系统动力响应的三维解析计算方法[J]. 力学学报, 2023, 55(6): 1329-1341. [百度学术]
HE C, JIA Y P, ZHOU S H. Three-dimensional analytical method for calculating vibrations of a coupled tunnel-soil-fluid system[J]. Chinese Journal of Theoretical and Applied Mechanics,2023,55(6): 1329-1341.(in Chinese) [百度学术]
DING W T,HUANG X H,DAI Z Y,et al.Analysis of the collapse mechanism and stabilization optimization of the composite stratum at the boundary between prereinforced and unreinforced areas near a shield launching area[J]. International Journal of Geomechanics, 2023, 23(5): 05023002. [百度学术]
王鑫, 罗富荣, 张飞, 等. 深孔注浆条件下地铁施工地层位移预测研究[J]. 地下空间与工程学报,2020, 16(5): 1521-1528. [百度学术]
WANG X,LUO F R,ZHANG F,et al. Study on stratum displacement prediction method induced by metro construction under deep hole grouting condition[J]. Chinese Journal of Underground Space and Engineering, 2020, 16(5):1521-1528.(in Chinese) [百度学术]
SHI C H,SUN X H, LIU S L,et al. Analysis of seepage characteristics of a foundation pit with horizontal waterproof curtain in highly permeable strata[J]. Water, 2021,13(9):1303. [百度学术]
侯乐乐,翁效林,黄文鹏,等.湿陷性黄土地铁隧道基底注浆加固处治试验[J].长安大学学报(自然科学版),2022,42(2):91-102. [百度学术]
HOU L L,WENG X L,HUANG W P,et al.Experimental on grouting reinforcement and treatment of metro tunnel basement in collapsible loess[J].Journal of Chang’an University (Natural Science Edition),2022,42(2):91-102.(in Chinese) [百度学术]
BIOT M A.Theory of propagation of elastic waves in a fluid-saturated porous solid. I. low-frequency range[J]. Journal of the Acoustical Society of America,1956,28(2):168-178. [百度学术]
BIOT M A. Theory of propagation of elastic waves in a fluid-saturated porous solid.Ⅱ.higher frequency range[J]. The Journal of the Acoustical Society of America,1956,28(2):179-191. [百度学术]