Loading [MathJax]/jax/output/SVG/fonts/TeX/Size1/Regular/Main.js
  • 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

埋置移动荷载作用下三维黏弹性半空间的动力响应

李怡君, 胡安峰, 李聪, 马玉喜, 徐浩

李怡君, 胡安峰, 李聪, 马玉喜, 徐浩. 埋置移动荷载作用下三维黏弹性半空间的动力响应[J]. 岩土工程学报, 2024, 46(S2): 145-149. DOI: 10.11779/CJGE2024S20021
引用本文: 李怡君, 胡安峰, 李聪, 马玉喜, 徐浩. 埋置移动荷载作用下三维黏弹性半空间的动力响应[J]. 岩土工程学报, 2024, 46(S2): 145-149. DOI: 10.11779/CJGE2024S20021
LI Yijun, HU Anfeng, LI Cong, MA Yuxi, XU Hao. Dynamic response of a visco-elastic half-space subjected to buried moving loads[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(S2): 145-149. DOI: 10.11779/CJGE2024S20021
Citation: LI Yijun, HU Anfeng, LI Cong, MA Yuxi, XU Hao. Dynamic response of a visco-elastic half-space subjected to buried moving loads[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(S2): 145-149. DOI: 10.11779/CJGE2024S20021

埋置移动荷载作用下三维黏弹性半空间的动力响应  English Version

基金项目: 

国家自然科学基金面上项目 52179110

详细信息
    作者简介:

    李怡君(1992—),女,讲师,主要从事交通岩土工程及土动力学等方面的研究。E-mail: yijunli@whpu.edu.cn

  • 中图分类号: TU435

Dynamic response of a visco-elastic half-space subjected to buried moving loads

  • 摘要: 采用积分变换方法对埋置移动荷载作用下三维黏弹性半空间地基振动的衰减规律进行了理论研究。首先,基于理想弹性介质的动力Navier方程,利用Helmholtz分解定理及傅里叶变换方法,在直角坐标系中求得了变换域内弹性介质动力响应分量的表达式。结合相应的边界条件及连续条件,得到了埋置移动荷载作用下黏弹性半空间地基动力响应的积分形式解。当荷载埋深为零时,埋置移动荷载退化为地表移动荷载,对比发现退化解结果与已有结果吻合较好。最后,采用IFFT法对地基土体动力响应进行数值求解,分析了移动荷载埋深、速度等因素对地基振动的分布衰减规律以及频谱分布的影响。
    Abstract: The attenuation of the vibration of a three-dimensional visco-elastic half-space subjected to the embedded moving loads is investigated theoretically. By employing the Fourier integral transform method and the Helmholtz decomposition, the Navier equations for the ideal elastic medium are solved, and the dynamic components in the transformed domain are derived in the Cartesian coordinate. Combined with the boundary conditions and continuity conditions, the analytical solutions to the dynamic response of the viscoelastic half-space due to the embedded moving loads are obtained. When the depth of the moving loads is zero, the results obtained in this study are in good agreement with the published ones. By employing the IFFT method, some numerical examples are selected to discuss the influences of the depth, velocity and viscous damping of the soil on the propagation of the vibration and the amplitude spectra of the displacement.
  • 随着城市轨道交通的快速发展,列车荷载引起的环境振动问题受到越来越多的关注。不少学者对移动荷载作用下地基的动力响应开展了研究。Eason[1]通过积分变换法求得了速度小于瑞丽波速的移动点荷载作用下均质弹性半空间位移、应力响应的数值解。Hung等[2]研究了亚音速、跨音速及超音速3种不同速度范围内的移动荷载作用下黏弹性半空间的动力响应规律,并分析了移动荷载形状分布及加振频率的影响。De等[3]采用传递反射矩阵法求得了成层弹性半空间在移动荷载的动力格林函数。自Biot[4-5]提出饱和多孔介质的动力控制方程后,多数学者都基于该控制方程对饱和土体动力响应进行研究。Jin等[6]利用梯形求积公式获得了振动荷载作用下半无限多孔饱和固体中的应力和孔隙水压力数值解。胡安峰等[7]利用Fourier积分变换求解了下卧基岩饱和半平面在移动线荷载作用下的动力响应。Xu[8]等研究了饱和成层土体上无限长Euler-Bernoulli梁在移动荷载作用下的动力响应。Cai等[9]研究了移动矩形荷载作用下饱和土体的动力响应。

    上述研究中,移动荷载大多作用在地基表面,而对地铁等地下构筑物来说,移动荷载作用在地表以下一定深度处。Pak[10]采用势函数的方法求解了弹性半空间内部点源荷载作用下的动力响应。Senjunctichai等[11]求得了二维均质多孔弹性半空间在内部荷载激励作用下的动力格林函数。陈胜立等[12]通过Hankel变换变换方法,求解了埋置点源荷载的轴对称Lame问题。杜秦文等[13]采用Hankel积分变换法,得到了埋置点源简谐荷载作用下Gibson土体的动力Green函数。上述文献中的埋置荷载多假定为固定位置动荷载。Metrikine等[14]通过建立一个二维弹性层-梁耦合的简化模型,对隧道中列车运行引起的地表振动规律进行了研究。Yuan等 [15]在Metrikine等[14]的基础上,求解了饱和地基-梁模型在移动荷载作用下的振动解析解。上述两个模型均为二维模型,忽略了垂直于荷载移动方向的地基振动。Forrest等[16]通过建立筒中筒(PIP)模型,在柱坐标系下对移动点荷载作用下圆形衬砌隧道埋置于弹性全空间中的动力响应情况进行了求解,但未考虑地表边界的影响。本文建立了三维黏弹性半空间地基模型,并通过引入势函数及利用三维傅里叶积分变换及逆变换方法,在笛卡尔坐标系下求解了埋置移动点荷载作用下黏弹性半空间的动力响应积分形式解。最后,通过数值算例分析了埋置荷载移动速度,埋置深度等因素对地表振动传播分布及衰减规律的影响。

    均质各向同性弹性体的控制方程以张量形式可表示为(不计体力):

    μui,jj+(λ+μ)uj,ji=ρ¨ui
    (1)

    式中:ui为弹性体的位移,i=x,y,zρ为弹性介质的质量密度;λμ为弹性体的Lame常数;(..)为对时间t的二阶导数。

    根据Helmholtz分解定理,引入标量势函数ϕ,矢量势函数ψj=(ψ1, ψ2,ψ3)。对位移场作如下分解:

    ui=ϕ,i+eijkψk,j
    (2)

    将式(2)代入式(1)可得如下关于势函数的表达式:

    2ϕ1cP22ϕt2=0 
    (3a)
    2Ψj1cS22Ψjt2=0
    (3b)

    式中:cPcS分别为弹性半空间体的压缩波与剪切波波速:cP=λ+2μρcS=μρ2为Laplace算子,2=2x2+2y2+2z2

    假定在直角坐标系中矢量势函数满足如下条件:

    μψj,j=0
    (3c)

    定义对xyt的三重Fourier变换及逆变换分别为

    ˆ¯¯g(k1,k2,z,ω)=+++g(x,y,z,t)eik1xeik2yeiωtdxdydt 
    (4a)
    g(x,y,z,t)=1(2π)3+++ˆ¯¯g(k1,k2,z,ω)eik1xeik2yeiωtdk1dk2dω
    (4b)

    将式(3a)~(3c)进行式(4a)中的Fourier积分变换可得:

    d2ˆ¯¯ϕdz2BP2ˆ¯¯ϕ=0  
    (5a)
    d2^¯¯Ψjdz2BS2^¯¯Ψj=0 
    (5b)
    d^¯¯Ψ3dz=(ik1^¯¯Ψ1+ik2^¯¯Ψ2)
    (5c)

    式中:BP2=k12+k22kP2BS2=k12+k22kS2kP=ωcPkS=ωcS。其中BpBs的实部大于零。

    由弹性介质本构方程及傅里叶变换公式可得波数-频率域内位移及应力响应分量表达式为

    ˆ¯¯ux=ik1ˆ¯¯ϕ+ik2ˆ¯¯ψ3ˆ¯¯ψ2z 
    (6a)
    ˆ¯¯uy=ik2ˆ¯¯ϕik1ˆ¯¯ψ3+ˆ¯¯ψ1z 
    (6b)
    ˆ¯¯uz=ˆ¯¯ϕz+ik1ˆ¯¯ψ2ik2ˆ¯¯ψ1 
    (6c)
    ˆ¯¯σzz=λ(2ˆ¯¯ϕz2k12ˆ¯¯ϕk22ˆ¯¯ϕ)+2μ(ik1ˆ¯¯ψ2zik2ˆ¯¯ψ1z+2ˆ¯¯ϕz2) 
    (6d)
    ˆ¯¯τxz=μ(2ik1ˆ¯¯ϕz+k1k2ˆ¯¯ψ1k12ˆ¯¯ψ22ˆ¯¯ψ2z2+ik2ˆ¯¯ψ3z) 
    (6e)
    ˆ¯¯τzy=μ(2ik2ˆ¯¯ϕz+k22ˆ¯¯ψ1+2ˆ¯¯ψ1z2k1k2ˆ¯¯ψ2ik1ˆ¯¯ψ3z) 
    (6f)

    图 1给出了本文的理论研究模型示意图。沿x轴正方向匀速移动的竖向点荷载作用在黏弹性半空间体表面以下深度z=h处的平面上,移动荷载表达式为F=peiω0tδ(xct)δ(y),其在变换域内荷载表达式为ˆ¯¯F=p2πδ(ω0+ω+ck1)。其中p为荷载幅值,c为荷载移动速度,ω0为荷载自振频率,δ为狄拉克函数。考虑荷载由x轴的负无穷处向正无穷处移动,当t=0时,荷载恰好移动到坐标原点处。

    图  1  埋置移动荷载作用于弹性半空间内部示意图
    Figure  1.  Diagram of elastic half-space with embedded moving loads

    根据波的辐射特性,将弹性半空间体分为上下两个区域:区域1(0zh)与区域2(h+z+),如图 1所示。则由式(3a)、(3b)可解得变换域内势函数可表示为

    0zh时,

    ˆ¯¯Φ=a1eBPz+a2eBPz ˆ¯¯Ψ1=a3eBSz+a4eBSz ˆ¯¯Ψ2=a5eBSz+a6eBSz ˆ¯¯Ψ3=ieBSzBS(a3k1+a5k2)+ieBSzBS(a4k1+a6k2)}
    (7a)

    h+z+时,

    ˆ¯¯Φ=a7eBPz ˆ¯¯Ψ1=a8eBSz ˆ¯¯Ψ2=a9eBSz ˆ¯¯Ψ3=ieBSzBS(a8k1+a9k2)}
    (7b)

    式中:aj为未知常数,j=19

    图 1可看出,半空间表面的边界条件及荷载作用平面处的连续条件可表示为

    z=0时,

    σzz(x,y,0)=0 τxz(x,y,0)=0 τyz(x,y,0)=0}
    (8)

    z=h时,

    σzz(x,y,h+)σzz(x,y,h)=F 
    (9a)
    τxz(x,y,h+)τxz(x,y,h)=0 
    (9b)
    τyz(x,y,h+)τyz(x,y,h)=0 
    (9c)
    uzz(x,y,h+)uzz(x,y,h)=0 
    (9d)
    uxz(x,y,h+)uxz(x,y,h)=0 
    (9e)
    uyz(x,y,h+)uyz(x,y,h)=0
    (9f)

    观察式(7)~(9)可看到,边界条件个数与连续条件个数之和等于未知数个数,故可进行求解。

    对边界条件与连续条件进行傅里叶变换后,将式(6a)~(6d)代入,可得到一组关于未知系数aj的线性方程组:Aijaj=fi,其中矩阵Aijfi的表达式见附录1。求得未知系数aj后,代入位移表达式, 即可得到变换域内的位移响应解ˆ¯¯ux,ˆ¯¯uy,ˆ¯¯uz。由Fourier逆变换即可求得稳态响应的表达式:

    {ux,uy,uz}=1(2π)3++ + δ(ω0+ω+ck1)ˆ¯¯ux,ˆ¯¯uy,ˆ¯¯uz                              eik1xeik2yeiωtdk1dk2dω
    (10)

    为考虑土体的黏滞性,引入如下黏滞阻尼模型:λ* = λ(1 + 2iβ)μ* = μ(1 + 2iβ),其中β为土体介质的黏滞阻尼比。由于式(10)中的动力响应解为无穷积分形式,且其被积函数为一个复杂的振荡函数,故采用IFFT方法进行求解。

    h=0时,本文解可退化为移动荷载作用在黏弹性半空间表面的动力响应解。Hung等[2]研究了黏弹性半空间在不同类型的地表移动荷载作用下的动力响应。土体参数取值密度ρ=2000 kg/m3,泊松比ν=0.25,阻尼比β=0.02,横波波速cs=100 m/s,纵波波速cp=173.2 m/s,表面波波速cR=92 m/s。定义观察点(x,y,z)=(0,0,z0)(其中z0=1 m)处的无量纲化竖向位移V=2πμz0uz/p,无量纲化水平纵向位移W=2πμz0ux/p图 2为荷载移动速度c=50 m/s时,本文位移响应退化解与Hung等[2]结果的对比图,可见两者吻合得较好。

    图  2  本文退化解与Hung等[2]的对比图
    Figure  2.  Comparison between degenerated solutions and results of Hung et al[2]

    定义无量纲位移ux=uxλl0/puz=uzλl0/p,其中l0为特征长度(l0=10 m)。土体Lame常数λ=2×107 N/m2μ=2×107 N/m2,土体密度ρ= 2000 kg/m3,黏滞阻尼比β=0.02,移动荷载幅值p=1 N。下面分析埋置移动点荷载作用下黏弹性半空间地表的振动情况。

    (1)地表竖向位移响应分布

    引入移动坐标系xt=xct图 3给出了不同速度的埋置移动荷载作用下,地表无量纲竖向位移uz在移动坐标系下的曲线分布图,其中y=0,移动荷载埋深分别取h为10,15,20 m。由图 3可看到,地表竖向位移响应关于xt=0基本呈轴对称分布,且荷载埋深越深,振动沿水平方向衰减越快,振动影响范围越小,地表振动幅值越小。此外,对比图 3(a)(b)可发现,荷载速度越大,地表最大竖向位移幅值也越大。

    图  3  地表竖向位移uz
    Figure  3.  Vertical displacements of surface uz

    (2)地表纵向位移响应分布

    在埋置移动荷载作用下,地表不仅会产生竖直方向的振动,也会产生水平方向的动力响应。图 4为当荷载埋深h=10 m时,地表x方向的纵向水平位移响应沿移动坐标xt的分布图。由图 4可看到埋置移动点荷载作用下,地表纵向水平位移响应曲线关于移动坐标轴呈反对称分布,且垂直于荷载移动方向,距离荷载作用位置越远,即y值越大,纵向水平位移越小;对比图 4(a)(b)可看出,随荷载速度增大,纵向水平位移幅值降低。同时,对比图 34可看出,同一埋置荷载作用下,地表土体竖向位移响应幅值远大于纵向水平位移响应幅值,这与实测地面振动规律吻合[17]

    图  4  地表纵向位移ux
    Figure  4.  Longitudinal displacements of surface ux

    (3)地表位移响应的频谱分析

    由式(10)可得原点(0,0,0)处的竖向位移响应表达式为

    uz(0,0,0,t)=1(2π)2+ + ˆ¯¯uz(k1,k2,0,ctω0)                     eik1(xctω0)eik2ydk1dk2
    (11)

    则该点处的竖向振动频谱可表示为

    uz(f)= + uz(0,0,0,t)ei2πftdt       =12π+ˆ¯¯uz(2πf02πfc,k2,0,2πf)dk2
    (12)

    同理可得原点处的水平位移频谱表达式ux(f)

    图 5给出了两种不同速度的常值移动荷载(h=10 m)作用下,地表原点处竖向位移与水平纵向位移响应的频谱图。可看到,当荷载速度较小时,位移频谱主要分布在频率较低的范围内,且竖向位移与水平位移几乎在同一频率达到峰值。当荷载速度变大时,位移频谱曲线的频率分布范围变大,位移幅值也增大,位移峰值所对应的频率也越大。

    图  5  地表位移频谱图
    Figure  5.  Displacement spectra of surface

    采用傅里叶变换及Helmholtz分解的方法对埋置移动荷载作用下弹性半空间地基的动力响应进行理论求解,并通过数值计算,分析了荷载埋深、速度等因素对地表位移响应的影响,得到3点结论。

    (1)移动荷载埋深越深,地表水平位移与竖向位移响应均越小;荷载移动速度越大,地表竖向位移响应幅值越大,纵向水平位移反而越小。

    (2)垂直于荷载移动方向,距离荷载作用位置越远,地表水平位移响应越小。

    (3)埋置荷载移动速度越大,位移响应频谱分布范围越大,位移幅值峰值对应的频率也越大。

  • 图  1   埋置移动荷载作用于弹性半空间内部示意图

    Figure  1.   Diagram of elastic half-space with embedded moving loads

    图  2   本文退化解与Hung等[2]的对比图

    Figure  2.   Comparison between degenerated solutions and results of Hung et al[2]

    图  3   地表竖向位移uz

    Figure  3.   Vertical displacements of surface uz

    图  4   地表纵向位移ux

    Figure  4.   Longitudinal displacements of surface ux

    图  5   地表位移频谱图

    Figure  5.   Displacement spectra of surface

  • [1]

    EASON G. The stresses produced in a semi-infinite solid by a moving surface force[J]. International Journal of Engineering Science, 1965, 2(6): 581-609. doi: 10.1016/0020-7225(65)90038-8

    [2]

    HUNG H H, YANG Y B. Elastic waves in visco-elastic half-space generated by various vehicle loads[J]. Soil Dynamics and Earthquake Engineering, 2001, 21(1): 1-17. doi: 10.1016/S0267-7261(00)00078-6

    [3]

    DE BARROS F C P, LUCO J E. Response of a layered viscoelastic half-space to a moving point load[J]. Wave Motion, 1994, 19(2): 189-210. doi: 10.1016/0165-2125(94)90066-3

    [4]

    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. doi: 10.1121/1.1908241

    [5]

    BIOT M A. Mechanics of deformation and acoustic propagation in porous media[J]. Journal of Applied Physics, 1962, 33(4): 1482-1498. doi: 10.1063/1.1728759

    [6]

    JIN B, YUE Z Q, THAM L G. Stresses and excess pore pressure induced in saturated poroelastic halfspace by moving line load[J]. Soil Dynamics and Earthquake Engineering, 2004, 24(1): 25-33. doi: 10.1016/j.soildyn.2003.09.004

    [7] 胡安峰, 孙波, 谢康和. 下卧基岩饱和地基在移动荷载作用下的动力响应[J]. 振动与冲击, 2012, 31(4): 151-156.

    HU Anfeng, SUN Bo, XIE Kanghe. 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)

    [8]

    XU B, LU J F, WANG J H. Dynamic response of an infinite beam overlying a layered poroelastic half-space to moving loads[J]. Journal of Sound and Vibration, 2007, 306(1/2): 91-110.

    [9]

    CAI Y Q, SUN H L, XU C J. Steady state responses of poroelastic half-space soil medium to a moving rectangular load[J]. International Journal of Solids and Structures, 2007, 44(22/23): 7183-7196.

    [10]

    PAK R Y S. Asymmetric wave propagation in an elastic half-space by a method of potentials[J]. Journal of Applied Mechanics, 1987, 54(1): 121-126. doi: 10.1115/1.3172945

    [11]

    SENJUNTICHAI T, RAJAPAKSE R K N D. Dynamic Green's functions of homogeneous poroelastic half-plane[J]. Journal of Engineering Mechanics, 1994, 120(11): 2381-2404. doi: 10.1061/(ASCE)0733-9399(1994)120:11(2381)

    [12] 陈胜立, 张建民, 陈龙珠. 饱和地基中埋置点源荷载的动力Green函数[J]. 岩土工程学报, 2001, 23(4): 423-426. doi: 10.3321/j.issn:1000-4548.2001.04.008

    CHEN Shengli, ZHANG Jianmin, CHEN Longzhu. Dynamic Green's functions of saturated soils subjected to the internal excitation[J]. Chinese Journal of Geotechnical Engineering, 2001, 23(4): 423-426. (in Chinese) doi: 10.3321/j.issn:1000-4548.2001.04.008

    [13] 杜秦文, 朱向荣, 丁伯阳. 埋置点源简谐荷载作用下Gibson土体动力Green函数[J]. 工程力学, 2007, 24(2): 45-49.

    DU Qinwen, ZHU Xiangrong, DING Boyang. Dynamic Green's function for Gibson soil subjected to internal excitation[J]. Engineering Mechanics, 2007, 24(2): 45-49. (in Chinese)

    [14]

    METRIKINE A V, VROUWENVELDER A C W M. Surface ground vibration due to a moving train in a tunnel: two-dimensional model[J]. Journal of Sound Vibration, 2000, 234(1): 43-66. doi: 10.1006/jsvi.1999.2853

    [15]

    YUAN Z H, BOSTRÖM A, CAI Y Q. Benchmark solution for vibrations from a moving point source in a tunnel embedded in a half-space[J]. Journal of Sound Vibration, 2017, 387: 177-193. doi: 10.1016/j.jsv.2016.10.016

    [16]

    FORREST J A, HUNT H E M. A three-dimensional tunnel model for calculation of train-induced ground vibration[J]. Journal of Sound Vibration, 2006, 294(4/5): 678-705.

    [17]

    ZOU C, WANG Y M, WANG P, et al. Measurement of ground and nearby building vibration and noise induced by trains in a metro depot[J]. The Science of the Total Environment, 2015, 536: 761-773. doi: 10.1016/j.scitotenv.2015.07.123

图(5)
计量
  • 文章访问数:  90
  • HTML全文浏览量:  5
  • PDF下载量:  17
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-06-21
  • 刊出日期:  2024-09-30

目录

/

返回文章
返回