基于非局部Biot理论的平面波作用下饱和地基动力响应
Dynamic Responses of Saturated Soil Foundation Subjected to Plane Wave Based on Nonlocal-Biot Theory
通讯作者: 丁海滨,男,博士生;E-mail:hbding@ecjtu.edu.cn.
责任编辑: 陈晓燕
收稿日期: 2020-05-22
基金资助: |
|
Received: 2020-05-22
作者简介 About authors
王宁(1983-),男,山东省东营市人,讲师,主要从事基础动力学与土动力学及动力检测技术应用等研究
为探究孔隙尺寸效应对饱和土地基动力响应的影响,基于非局部-Biot理论,构建了P波及SV波作用下饱和土地基表面动力响应的计算模型.采用波函数展开法,求解了地表位移及应力的解析解.分析了孔隙尺寸(以非局部参数描述)、入射波频率及入射角对地表位移及应力响应的影响.结果表明:低频下,非局部Biot理论计算结果与经典Biot理论计算结果基本一致;高频时,地表位移及应力随非局部参数变化较为明显,即高频率下,孔隙尺寸效应对地表响应影响不可忽略.入射波频率对地表响应的影响与孔隙尺寸有关,孔隙越大,频率影响越为显著.相同幅值的P波及SV波作用下,SV波引起的地表动力响应大于P波,且SV波在入射角为45° 时发生全反射.研究成果可为半空间饱和土中波动问题的研究提供借鉴.
关键词:
To investigate the influence of pore size effect on dynamic response of saturated soil foundation, a model for predicting the dynamic response of the ground surface of the saturated soil foundation under incident P wave and SV wave is proposed based on the nonlocal-Biot theory. The analytical solution is obtained using the wave function expansion method. The influence of pore size described by nonlocal parameter, input frequency, and the incident angle on the dynamic response of displacement and stress is discussed in detail. The results show that at low frequencies, the calculation results of the nonlocal-Biot theory are basically the same as those of the classical Biot theory. At high frequencies, the surface displacement and stress change significantly with nonlocal parameters, that is, at high frequencies, the effect of pore size on the surface response cannot be ignored. The influence of incident wave frequency on the ground-surface response is related to pore size, that is, the larger the pore size, the more significant the frequency effect. The influence of SV wave on the dynamic response of ground surface is larger than that of P wave. Besides, the total reflection phenomenon is observed at an incident angle of 45° for the incident of SV wave. The results obtained in this work can provide reference for studying the problem of wave propagation in half-space saturated soil foundation.
Keywords:
本文引用格式
王宁, 丁海滨, 童立红, 蒋亚龙.
WANG Ning, DING Haibin, TONG Lihong, JIANG Yalong.
为解决上述问题,Chakraborty等[11]结合经典非局部弹性理论[12]和Biot理论[4]预测了饱和多孔材料中波的传播特性,但其由于考虑了流体的非局部效应,对Biot理论进行过修正,而使其预测结果与实验结果产生一定的偏差.Tong等[13]摒弃了流体部分非局部效应,提出了仅考虑土骨架非局部效应的非局部-Biot理论模型,分析了饱和土中波速的传播及其衰减规律,其结果与Bouzidi等[9]的试验结果吻合得很好,以此说明了所提出模型的合理性.随后,徐长节等[14]、Xu等[15]及Ding等[16]基于非局部-Biot理论研究了饱和土中深埋、浅埋圆形衬砌及浅埋复合式衬砌对入射P波的散射问题.结果表明,高频下孔隙尺寸对衬砌动应力影响显著,而低频情况下,非局部-Biot理论解与经典Biot理论解基本一致.Tong等[17]利用非局部-Biot理论,分析了非局部参数对饱和土中Rayleigh波的传播特性的影响,结果仍表明高频下孔隙尺寸对波传播特性的影响不可忽略.由以上分析可知,高频下,饱和土中孔隙尺寸及孔隙动应力会对波的传播特性产生较大的影响,为此,研究孔隙尺寸在平面波作用下对地表动力响应的影响显得尤为重要.
本文基于非局部-Biot理论,构建了P波及SV波入射下,饱和土地基地表响应计算模型,采用波函数展开法,得出了地表位移及水平应力的解析解,通过算例分析研究了非局部参数、入射波频率及入射角对地表位移及应力响应影响的变化规律.
1 计算模型
假设平面P1(SV)波以入射角为
图1
P1波入射时:
SV波入射时:
式中:i、r分别代表入射波及反射波;
2 非局部-Biot理论控制方程
根据文献[13]可知,饱和孔隙弹性介质的非局部-Biot理论的运动方程为
本构方程:
式中:
运动方程:
式中:ρ为土体密度;ρf 为流体密度;m=ρf/n0,b=ηF(ξ)/k,η为流体黏滞系数,k为流体渗透系数;F(ξ)为高频下黏性修正系数,其表达式为[18]
式中:
3 波场的求解
3.1 饱和土中散射波场的求解
为求解式(4),引入土骨架标量和矢量势函数φs和Ψs,及流体部分标量及矢量势函数φf和Ψf,根据Helmholtz分解定理,饱和土中位移场可表示为
将式(6)带入式(4),并消去φf和Ψf,可得:
式中:kb为剪切波波数,满足
参数关系为
其中:ω为圆角频率.
由式(7)可知,饱和多孔介质中存在两种P波,土骨架部分的P波势函数可表示为
同理,由式(4)和(6)可得流体部分P波和SV波势函数为
式中:
3.2 饱和土中总波场
饱和土中总波场由入射波场和散射波场构成,本文考虑入射波场分别为P波和SV波情况.
入射波为P波时:
入射波为SV波时:
此处需要说明的是,入射波为SV波时,由于P1波波速通常大于SV波波速,当SV波的入射角达到临界角θcr1时,P1波的反射角将大于90°,此时,P1波表现出随深度增加,其幅值呈指数减小的非均匀波.与此同时,在软黏土介质中,P2波波速也有可能大于SV波波速,此时将存在第二个临界角θcr2.由波的反射关系,可得两个临界角为
为便于分析,将式(2)表示为
式中:K1、K2 及K3为反射波待定系数,
4 边界条件及待定系数求解
4.1 应力位移与势函数关系
由式(3)可知,非局部-Biot应力张量与经典Biot理论应力张量满足
由于非局部参数τ为很小量,所以,可将高阶无穷小项
即
结合式(3)、(6)及(14),应力、位移与势函数关系可表示为
式中: ux 、uy 分别为x、y方向土骨架位移;wx 、wy 分别为x、y方向流体相对土骨架位移;σxx 、σyy 、σxy 分别为x、y方向正应力与切应力;Pf为孔隙水压力.
4.2 边界条件
考虑到地表零应力及透水特性,可将其边界条件表示为
将式(1)和(13)分别代入式(15),并结合式(16)可求出P波及SV波入射下势函数的待定系数,待定系数求解方程组如下[3].
P波入射:
SV波入射:
系数矩阵表达式见表1.
表1 系数矩阵表达式
Tab.1
系数 | P波 | SV波 |
---|---|---|
G1i | [-1+τ2(k2- | |
G13 | μ | -2jkμνβ[-1+τ2(k2- |
G2i | -μ | 2jkμναi[-1+τ2(k2- |
G23 | μ | μ(k2-k4τ2+ |
G3i | M(k2- |
通过式(17),并结合表1即可求出反射系数,从而可利用式(15)求解出饱和土中位移及应力表达式.
5 结果与讨论
实际工程中人们所关心的问题是地表的位移及应力变化,地表位移可采用式(15)中第1式计算,所得出的位移为与入射波幅值的比值(本文取入射波幅值为φ0=1),因此后续分析中的位移为无量纲位移,针对P波和SV波入射分别采用ka1和kβ对其位移进行无量纲化.地表应力采用无量纲应力,即
式中:
5.1 结果验证
图2
5.2 算例分析
为充分分析P波及SV波作用下饱和地基地表动力响应情况,本文选取饱和土体计算参数,见表2所示.
表2 饱和土体物理力学参数
Tab.2
参数 | 数值 | 参数 | 数值 |
---|---|---|---|
土颗粒密度ρs/(kg·m-3) | 2650 | 渗透系数k/(m-2) | 1.0×10-8 |
流体密度ρf/(kg·m-3) | 1000 | 黏滞系数η/(Pa·簚s) | 1.0×10-3 |
孔隙比n0 | 0.30 | Biot 参数α | 0.9982 |
泊松比ν | 0.30 | Biot 参数M/MPa | 5859.2 |
土颗粒模量Kg/GPa | 36.0 | 弯曲因子ξ | 1 |
流体体积模量Kf/GPa | 2.0 | 孔隙半径a/μm | 23 |
拉梅常数λ/MPa | 49.7 |
图3、4分别为P波及SV波作用下,地表位移及应力随非局部参数的变化曲线(入射角为30°).由图可知,当入射频率f=10 Hz时, 地表位移及应力随非局部参数的增加几乎没有变化,而随着入射波频率的增加,非局部参数对地表位移及应力具有显著的影响.如当入射波频率达到 1000 Hz,P波入射情况下,地表位移随非局部参数出现增加而增大,应力随非局部参数增加而呈现出减小的趋势.SV波入射情况下,地表水平位移及水平应力均随非局部参数的增加而减小,其竖向位移随非局部参数的增加而增大.其原因是非局部参数的增加,意味着饱和土中孔隙尺寸的增大,入射频率较低(如10 Hz)时,波长远大于孔隙尺寸,此时孔隙尺寸对波传播的影响可以忽略,而随着入射波频率(如 1000 Hz)的增加,饱和土中波长减小,此时随着孔隙尺寸的增加,波长与孔隙尺寸更加接近,导致孔隙尺寸对波传播特性产生较大的影响.此外,从图3和4中还可以看出,随着入射波频率的增加,地表水平位移及竖向位移的起始值逐渐增大,这是由于入射波频率增加引起地表位移增大,而入射波频率对地表应力起始值影响不大.
图3
图3
P波入射下地表位移及应力随非局部参数的变化
Fig.3
Displacement and stress of ground surface versus nonlocal parameter subjected to P wave
图4
图4
SV波入射下地表位移及应力随非局部参数变化
Fig.4
Displacement and stress of ground surface versus nonlocal parameter subjected to SV wave
图5、6分别为不同非局部参数时,P波及SV波入射下,地表位移及应力随入射波频率变化曲线,其中P波及SV波入射角均为30°.由图可知,当入射波频率较低时,本文理论计算结果与经典Biot理论计算结果(τ=0.00 m)基本一致,而随着入射波频率的增加,本文理论计算结果与Biot理论计算结果逐渐显现出差异.由此说明,高频下孔隙尺寸对波传播特性的影响不可忽略.倘若以与经典Biot理论偏差达5%为分界点,由图可以看出,频率分界点(与Biot理论偏差达到5%的频率点)随非局部参数的增大而减小.如SV波入射下,地表水平位移在非局部参数为0.06、0.08及0.10 m时,所对应的临界频率分别为892、672及538 Hz,其原因是非局部参数越大,其对应的饱和土孔隙越大,从而导致其对波传播特性的影响越大.
图5
图5
P波入射下地表位移及应力随入射频率变化
Fig.5
Displacement and stress of ground surface versus input frequency subjected to P wave
图6
图6
SV波入射下地表位移及应力随入射频率变化
Fig.6
Displacement and stress of ground surface versus input frequency subjected to SV wave
图7、8分别为P波及SV波作用下,地表应力及位移随入射波入射角度的变化曲线,入射频率为1000 Hz.由图可知,P波入射下,地表水平及竖向位移随非局部参数的增加而增大,而应力随非局部参数的增大而减小.SV波入射时,地表水平位移及应力在入射角为45° 时减小为0,这是由于SV波在45° 入射时, 在地表发生全反射,此结论与经典弹性解相同[19].在45° 之前水平位移随非局部参数增加呈现出增加的趋势,在此之后随非局部参数的增加而呈现出减小趋势.地表竖向位移随非局部参数增加而增大.此外,图中虚线内出现的拐点随非局部参数增大而向左移动,这是由于非局部参数增加使得临界角度减小所致.对比图7、8可知,对于相同幅值的P波、SV波作用下,SV波所引起的地表响应大于P波.
图7
图7
P波入射下地表位移及应力随入射角变化曲线
Fig.7
Displacement and stress of ground surface versus incident angle subjected to P wave
图8
图8
SV波入射下地表位移及应力随入射角变化
Fig.8
Displacement and stress of ground surface versus incident angle subjected to SV wave
6 结论
基于非局部-Biot理论,利用波函数展开法,求解了P波和SV波作用下,饱和土地基地表响应问题的解析解.将本文解退化为经典Biot理论下无黏性解,并与文献[3]对比,结果一致.通过算例分析研究了非局部参数、入射波频率及入射角对地表动力响应的影响,得出如下结论:
(1) 非局部参数对地表动力响应的影响与入射波频率有关,入射波频率较低(10 Hz)时,本文解与经典Biot理论解基本一致,随入射波频率增加(如 1000 Hz),非局部参数对地表影响的影响较为显著.其原因是频率越高,波长越小,孔隙尺寸对波传播特性的影响越大.
(2) 入射波频率对地表响应的影响与非局部参数有关,倘若以与经典Biot偏差5%为分界点,可知非局部参数越大,频率分界点越小.
(3) P波入射时,地表位移随非局部参数的增大而增大,应力随非局部参数的增大而减小.SV波入射时,地表水平位移及应力在入射角为45° 时,其值为0,这是由于SV波在45° 时发生全反射所致.此外,非局部参数通过减小SV波入射的临界角,而使地表响应的拐点向左移动.相同幅值的P波和SV波作用下,SV波所引起的地表响应大于P波.
本文理论及研究成果可为半空间饱和土中波动问题的研究提供借鉴.
参考文献
平面波入射下深水地基场地动力响应分析
[J].
Dynamic responses of an underwater site subjected to plane P- or SV-wave incidence
[J].
Saturation effects on horizontal and vertical motions in a layered soil-bedrock system due to inclined SV waves
[J].DOI:10.1016/S0267-7261(01)00015-X URL [本文引用: 1]
Wave propagation in a poroelastic half-space saturated with inviscid fluid
[D].
Mechanics of deformation and acoustic propagation in porous media
[J].DOI:10.1063/1.1728759 URL [本文引用: 2]
Macroscopic Lagrangian formulation of poroelasticity with porosity dynamics
[J].DOI:10.1016/j.jmps.2004.05.005 URL [本文引用: 1]
Wave propagation in 3-D poroelastic media including gradient effects
[J].DOI:10.1007/s00419-012-0675-8 URL [本文引用: 1]
Transient dynamic analysis of a fluid-saturated porous gradient elastic column
[J].DOI:10.1007/s00707-011-0539-2 URL [本文引用: 1]
On wave propagation in gradient poroelasticity
[J].DOI:10.1016/j.soildyn.2016.05.010 URL [本文引用: 1]
Measurement of the speed and attenuation of the Biot slow wave using a large ultrasonic transmitter
[J].
Frequency dependencies of phase velocity and attenuation coefficient in a water-saturated sandy sediment from 0.3 to 1.0 MHz
[J].DOI:10.1121/1.2713690 URL [本文引用: 1]
Prediction of negative dispersion by a nonlocal poroelastic theory
[J].DOI:10.1121/1.2816576 URL [本文引用: 1]
On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves
[J].DOI:10.1063/1.332803 URL [本文引用: 1]
On wave propagation characteristics in fluid saturated porous materials by a nonlocal Biot theory
[J].DOI:10.1016/j.jsv.2016.05.042 URL [本文引用: 3]
基于非局部Biot理论下饱和土中深埋圆柱形衬砌对平面弹性波的散射
[J].
Scattering wave generated by the cylindrical lining in saturated soil based on the nonlocal-Biot theory
[J].
Scattering of a plane wave by shallow buried cylindrical lining in a poroelastic half-space
[J].DOI:10.1016/j.apm.2019.01.029 URL [本文引用: 1]
Dynamic responses of shallow buried composite cylindrical lining embedded in saturated soil under incident P wave based on nonlocal-Biot theory
[J].DOI:10.1016/j.soildyn.2019.02.018 URL [本文引用: 1]
Nonlocal scale effect on Rayleigh wave propagation in porous fluid-saturated materials
[J].DOI:10.1016/j.ijmecsci.2018.08.028 URL [本文引用: 1]
Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher frequency range
[J].DOI:10.1121/1.1908241 URL [本文引用: 1]
Surface amplitudes of reflected body waves
[J].DOI:10.1190/1.1438425 URL [本文引用: 1]
/
〈 |
|
〉 |
![](../../../richhtml/1006-2467/richHtml_jats1_1/images/1.png)