Dynamic response of utility tunnel during the passage of Rayleigh waves

The modeling methodologies and calculation of dynamic response of underground structure under Rayleigh waves is investigated in this paper. First the free field responses under Rayleigh waves are analyzed and the numerical results agree well with the theoretical results. Then, the approximate Rayleigh waves are put forward based on the preliminary re-search, and Rayleigh wave field is obtained through fast Fourier transform technique. Taking a utility tunnel as an example, its dynamic responses under Rayleigh waves is calculated by ABAQUS. The results demonstrate that bending deformation is the main component of structural deformation and the deformation at the top of the structure is about twice as much as that at bottom of the structure. The effect of soil-structure interface and the buried depth of underground structure are also investi-gated via parameter analysis. For the shallow buried underground structures, Rayleigh waves can be the key factor to control the responses and damage of the structure.     

A plane strain soil-structure interaction model

A plane strain model for dynamic soil-structure interaction problems under harmonic state is presented. The boundary element method is used to study the response of a homogeneous isotropic linear elastic soil. The far field displacement at the free surface is approximated by an outgoing Rayleigh wave. The finite element method is used to describe the response of the building, of the foundation and possibly of a finite part of the inhomogeneous non-linear soil. Two coupling procedures are described. The model is applied to a problem previously studied in the antiplane case. Incident P, SV and Rayleigh waves are considered. The results show an amplification and an attenuation of the structure motion with frequency when incident Rayleigh waves and P, SV body waves are respectively considered.     

Flexible support of a pile embedded in unsaturated soil under Rayleigh waves

Under the action of Rayleigh waves, pile head is easy to rotate with a concrete pile cap, and pure fixed-head condition is rarely achieved, which is a common phenomenon for it usually occurs on the precast piles with insufficient anchorage. In addition, the propagation characteristics of Rayleigh wave have been changed significantly due to the existence of capillary pressure and the coupling between phases in unsaturated soil, which significantly affects the pile-soil interaction. In order to study the above problems, a coupled vibration model of unsaturated soil–pile system subjected to Rayleigh waves is established on the basis that the pile cap is equivalent to a rigid mass block. Meanwhile, the soil constitution is simplified to linear-elastic and small deformations are assumed to occur during the vibration phase of soil–pile system. Then, the horizontal dynamic response of a homogeneous free-field unsaturated soil caused by propagating Rayleigh waves is obtained by using operator decomposition theory and variable separation method. The dynamic equilibrium equation of a pile is established by using the dynamic Winkler model and the Timoshenko beam theory, and the analytical solutions of the horizontal displacement, rotation angle, bending moment and shear force of pile body are derived according to the boundary conditions of flexible constraint of pile top. Based on the present solutions, the rationality of the proposed model is verified by comparing with the previous research results. Through parametric study, the influence of rotational stiffness and yield bending moment of pile top on the horizontal dynamic characteristics of Rayleigh waves induced pile is investigated in detailed. The analysis results can be utilized for the seismic design of pile foundation under Rayleigh waves.     

Seismic response of pile groups under oblique-shear and Rayleigh waves

A simple analytical solution is presented to calculate the pile-soil-pile interaction and eventually the dynamic response of pile groups when excited by the passage of Rayleigh waves and obliquely incident SH-waves. A dynamic Winkler model, with realistic frequency-dependent ‘springs’ and ‘dashpots’ in conjunction with physically motivated approximations is utilized to compute the wave field radiating from an oscillating pile and the effect of this field on an adjacent pile. The coupled rocking motion of the pile group resulting from Rayleigh waves and the torsional motion of the pile group resulting from SH-waves is accurately predicted by a simple mathematical expression. The results of the presented method can be obtained with ‘hand calculations’ and are in excellent agreement with results from ‘rigorous’ solutions based on integral equation formulations. It is found that the group response is primarily affected from the phase difference of the input seismic motion at the location of each pile (wave-passage effect). Pile-soil-pile interaction has insignificant effect and can be neglected.     

地下综合管廊土-结构接触面参数对地震动力响应的影响数值分析

为研究土-结构接触面参数对地下综合管廊地震动力响应特征的影响,建立动力有限元数值模型,模型边界采用激励侧固定边界、远离激励侧黏性边界、其余侧自由场边界的优化组合动力边界,土体本构采用HSS模型,接触面采用改进Goodman单元,动力荷载考虑三种情况(Rayleigh波的作用、底部激励了美国加利福尼亚Upland地震波以及前两者的共同作用),分别研究不同地震动输入、接触面折减系数的改变对综合管廊内力及加速度的影响。研究结果表明:在相同的折减系数条件下,与静力作用相比,动力作用下的结构内力明显增大,综合管廊设计时应考虑地震荷载作用下内力增大的情况;随着界面折减系数的增加,正弯矩极值减小,负弯矩极值增大,加速度峰值增大;在相同接触面折减系数条件下,底部地震波输入产生的结构内力极值显著高于仅有Rayleigh波输入的情况;考虑Rayleigh波和地震波共同作用条件下,引起的管廊结构内力极值与仅考虑底部地震波输入时的结构内力极值差异不大。研究成果可供地下综合管廊结构地震响应精细化数值模拟及抗震设计参考。     

A comparative analysis of seismic response of shallow buried underground structure under incident P,SV and Rayleigh waves

In this study, A time-domain seismic response analysis method and a calculation model of the underground structure that can realize the input of seismic P, SV and Rayleigh waves are established, based on the viscoelastic artificial boundary elements and the boundary substructure method for seismic wave input. After verifying the calculation accuracy, a comparative study on seismic response of a shallow-buried, double-deck, double-span subway station structure under incident P, SV and Rayleigh waves is conducted. The research results show that there are certain differences in the cross-sectional internal force distribution characteristics of underground structures under different types of seismic waves. The research results show that there are certain differences in the internal force distribution characteristics of underground structures under different types of seismic waves. At the bottom of the side wall, the top and bottom of the center pillar of the underground structure, the section bending moments of the underground structure under the incidences of SV wave and Rayleigh wave are relatively close, and are significantly larger than the calculation result under the incidence of P wave. At the center of the side wall and the top floor of the structure, the peak value of the cross-sectional internal force under the incident Rayleigh wave is larger than the calculation result under SV wave. In addition, the floor of the underground structure under Rayleigh waves vibrates in both the horizontal and vertical directions, and the magnification effect in the vertical direction is more significant. Considering that the current seismic research of underground structures mainly considers the effect of body waves such as the shear waves, sufficient attention should be paid to the incidence of Rayleigh waves in the future seismic design of shallow underground structures.     

地震槽波的数学-物理模拟初探

针对地震槽波在低速层的传播特性,开展了煤层内地震槽波勘探的数值模拟和物理模拟研究的初探工作.在数值模拟研究方面,采用交错网格有限差分法对煤层中的地震槽波进行三分量全波场模拟.基于波场快照和人工合成地震记录研究了不同模型中的波场特征和各种波型的传播规律.在物理模拟方面,通过选用不同配比的环氧树脂和硅橡胶类材料构建地震槽波物理模型,利用透射法和反射法观测系统获得了清晰的地震槽波记录以研究槽波的地震学特征.研究表明,在煤层内槽波的地震波场中,Love型槽波的能量小于Rayleigh型槽波的SV分量,大于Rayleigh型槽波的SH分量.相对于Love型槽波和Rayleigh型槽波的SH分量,Rayleigh型槽波的SV分量在围岩中的泄露能量较强.在煤层界面附近的围岩中,地震波仍以槽波形式传播,随着距离的增加能量逐渐衰减.随着煤层变薄,煤层槽波主频向高频方向移动,频散现象增强,传播速度增大.     

Forward modeling of Rayleigh surface waves for analytical characterization of dominant dispersion trends

Forward modeling is of critical importance for inversion analysis of surface wave methods to obtain shear-wave velocity (V_S) profiles of soil sites. The dynamic stiffness matrix (DSM) method can provide forward modeling of Rayleigh surface waves to simulate complex wave propagation in layered soil sites. However, contamination from body waves and interference of multiple Rayleigh wave modes can reduce the accuracy of theoretical dispersion curves, especially at irregular soil sites with embedded low-velocity or high-velocity layers. An analytical method is developed herein to combine the techniques of the multichannel analysis of surface waves method with the DSM method to improve the accuracy of the theoretical dispersion analysis for soil sites. The proposed method implements multichannel analysis of the analytical displacement responses to capture dominant dispersion trends. Comparison of the results obtained with the new method against those from the transfer matrix method and the literature indicates that the new method can (1) effectively minimize the effects of contamination caused by body waves and interference from several Rayleigh wave modes, and (2) generate accurate dominant dispersion trends for soil sites with various stiffness profiles, especially for the high-frequency dispersion characteristics of the profiles with embedded low-velocity layers.     

多道瞬态瑞雷波技术在公路采空塌陷区探测中应用

地下采空塌陷区岩体结构具有明显的结构不连续、松散、断裂、空洞等异常结构特征,本文在分析公路采空塌陷区地球物理场特征的基础上,详细讨论了应用多道瞬态瑞雷波技术探测采空塌陷区的可行性以及相关技术,深入分析基于瑞雷波波场特征、频谱结构、能量衰减、频散特征、相速度以及多阶模等特征探测和评价地层结构的均匀性、地下空洞及异常结构体等地下障碍物的存在及其空间分布的解释标志.滚动式多道瑞雷波探测技术可以有效地探测地下采空区、软弱结构体等异常结构体的赋存状态及空间分布.     

On the interaction of rayleigh surface waves with structures

A two-dimensional soil-structure interaction analysis is carried out for transient Rayleigh surface waves that are incident on a structure. The structure is modelled by a three-degree of freedom rigid basemat to which is attached a flexible superstructure, modelled by a single mass-spring system. The structural responses to a given Rayleigh wave train are compared with those that would have been obtained if the free-field acceleration-time history had been applied as a normally incident body wave. The results clearly exhibit the [frequency filtering] effects of the rigid basemat on the incident Rayleigh waves. It is shown that, if seismic excitation of a structure is, in fact, due to Rayleigh surface waves, then an analysis assuming normally incident body waves can considerably over-estimate structural response, both at basemat level for horizontal and vertical motions and for vertical oscillations of the superstructure. However, in the examples considered here, relatively large rocking effects were induced by the Rayleigh waves, thus giving maximum horizontal accelerations in the superstructure that were of comparable magnitude for Rayleigh and normally incident body waves.     

处于海洋环境的桩柱在瑞利波作用下的响应

本文以一海洋深水环境桩柱为例,考虑到海洋地基为两相饱和土介质,基于Biot两相饱和多孔介质理论,通过数值算例分析了饱和介质中瑞利波加速度、位移分布规律并与单相弹性介质地基条件下进行了对比,结果表明,瑞利波主要作用于地表附近且振动幅度随深度衰减很快,饱和地基下的瑞利波影响深度要大于单相弹性地基,在瑞利波地震作用下,考虑两相饱和介质下的桩柱位移响应大于单相弹性介质地基。并且由于桩柱处海洋深水环境,本文对环境水体对桩柱振动响应的影响进行了研究,结果表明,流固耦合效应可明显降低桩柱运动响应,在海洋深水环境地震响应分析时,有必要考虑饱和地基以及流固耦合效应对桩柱运动响应的影响。     

Magneto-elastic waves and disturbances in initially stressed conducting media

Summary The object of the present paper is to investigate magneto-elastic waves and disturbances in initially stressed conducting media. Firstly, the theory of magneto-elastic surface waves in a conducting medium under an initial uniaxial tension has been deduced and then it has been employed to investigate the particular cases of surface waves such as Rayleigh, Love and Stoneley waves. Secondly, propagation of waves in an elastic layer has been considered using the fundamental equations of motion for magneto-elastic waves in conductors under an initial uniaxial tension or compression. This is followed by the case of plane Lamb's problem in a magneto-elastic semi-space under the same initial tension. The final results obtained in the above cases are in agreement with the corresponding classical problems when the initial tension is zero and the magnetic field is absent.     

小应变硬化模型(HSS)在Rayleigh波作用下场地响应分析中的应用

为研究冲击荷载或地震作用下产生的,以Rayleigh波为主的面波对浅层地表土体动力响应特征以及数值模拟中土层阻尼的设置方法,以厦门地区浅层的素填土及粉质黏土为研究对象,采用有限元动力分析,土体本构采用小应变硬化模型(HSS),利用模型本身的滞回环特性,输入变化的小应变参数,考察HSS模型的小应变参数对场地动力响应的影响,并与土体采用摩尔-库伦模型结合Rayleigh阻尼("MC+Rayleigh阻尼")的计算结果进行对比。研究表明:当采用带有滞回环的HSS模型时,波速随初始剪切模量Gref0的增大而增大,但振幅减小,残余变形量也有所减小;小应变参数γ0.7对波的影响较小;HSS模型能够给出残余变形量,而"MC+Rayleigh阻尼"由于本构模型为理想弹塑性模型,在卸载重加载条件下表现为纯弹性行为,无法反映出卸载重加载过程中塑性应变的积累及其累积阻尼效应;但HSS模型还不能够全面反映循环加载作用下塑性体积应变的累积,因此在考虑滞回阻尼的基础上,仍然建议借助Rayleigh阻尼来更加全面地模拟土体的实际阻尼特性。     

Time-domain solution for transient dynamic response of a large-diameter thin-walled pipe pile

The propagation of stress waves in a large-diameter pipe pile for low strain dynamic testing cannot be explained properly by traditional 1D wave theories. A new computational model is established to obtain a wave equation that can describe the dynamic response of a large-diameter thin-walled pipe pile to a transient point load during a low strain integrity test. An analytical solution in the time domain is deduced using the separation of variables and variation of constant methods. The validity of this new solution is verifi ed by an existing analytical solution under free boundary conditions. The results of this time domain solution are also compared with the results of a frequency domain solution and fi eld test data. The comparisons indicate that the new solution agrees well with the results of previous solutions. Parametric studies using the new solution with reference to a case study are also carried out. The results show that the mode number affects the accuracy of the dynamic response. A mode number greater than 10 is required to enable the calculated dynamic responses to be independent of the mode number. The dynamic response is also greatly affected by soil properties. The larger the side resistance, the smaller the displacement response and the smaller the refl ected velocity wave crest. The displacement increases as the stress waves propagate along the pile when the pile shaft is free. The incident waves of displacement and velocity responses of the pile are not the same among different points in the circumferential direction on the pile top. However, the arrival time and peak value of the pile tip refl ected waves are almost the same among different points on the pile top.     

基于萤火虫算法的雷瑞波非线性反演(英文)

雷瑞波具有强振幅、低频和低速的特点,在反射地震勘探中通常是需要被压制的强噪声。本文研究如何利用雷瑞波获取近地表地层的横波速度和地下结构,选取萤火虫优化算法进行面波的反演,萤火虫优化算法是一种新的粒子群算法理论,具有稳定、快捷、全局搜索等特点。针对萤火虫优化算法优缺点进行了讨论和改进,通过对理论模型和野外数据的测试应用,将提取的瑞利面波频散曲线反演得到横波速度信息。结果表明萤火虫优化算法能实现面波非线性反演,并具有分辨率高、抗干扰能力强等优点和实际使用前景。     

煤巷小构造Rayleigh型槽波超前探测数值模拟

对煤巷小构造地震波场进行了数值模拟研究,分析了层状煤层中地震波的传播特征.研究表明:(1)在煤巷迎头前方煤层内以纵波震源激发的Rayleigh型槽波相对于体波能量较强,波列较长,波速较低.(2)沿煤层传播的Rayleigh型槽波在小构造面上产生Rayleigh型槽波反射波,反射Rayleigh型槽波垂直分量相对于水平分量能量较强.沿煤层反向传播的反射Rayleigh型槽波在煤巷迎头面上转换为沿煤巷底板传播的Rayleigh面波.沿煤巷底板可以接收到能量较强的反射Rayleigh型槽波产生的Rayleigh面波,其可以作为超前探测小构造面的特征波.在地震记录上反射Rayleigh型槽波产生的Rayleigh面波波至最迟,在时间域与其他波列时间间隔较大,其垂直分量能量相对于水平分量较强,在地震记录上容易识别.(3)在相同的地质条件下应用反射地震超前探测方法,标志煤巷迎头前方存在小构造面的反射地震波能量较弱,受煤巷顶、底板界面和采煤迎头面的强反射波干扰,在地震记录中难以识别.     

Rayleigh波场的数值模拟及其应用

本文总结了Rayleigh波场数值模拟的4种思路和各自特点,根据Lamb问题的理论分析成果提出了基于地表激振的Rayleigh波场数值模拟技术,论述了地表集中震源作用下引起地表波动场的特点及主要影响因素。在此基础上,利用Plaxis 2D有限元软件实现了Rayleigh波场的数值模拟,并结合算例验证了所提方法的可行性和结果的合理性。最后,针对多层建筑结构,研究了不同地震动输入模式下结构动力反应的特点。结果表明,Rayleigh波作用下结构的动力反应特性明显区别于在底部输入剪切波时的结果,不同地震动输入模式对结构的振动形态和破坏模式有着明显的影响。     

Spatial behaviour of Rayleigh waves in layered half‐spaces under active surface sources

In the free state, Rayleigh waves are assumed to travel in the form of planar wavefronts. Under such an assumption, the propagation behaviour of the modes of Rayleigh waves in layered half‐spaces is only frequency dependent. The frequency behaviour, which is often termed as dispersion, is determined by the shear wave velocity profile of layered soils within the depth related to wavelength (or frequency). According to this characteristic, the shear wave velocity profile can be back‐analysed from the dispersion. The technique is widely used in the surface wave testing. However, the wavefronts of Rayleigh waves activated by the surface sources are non‐planar. The geometric discrepancy could result in Rayleigh waves manifesting distance‐dependent behaviour, which is referred to as spatial behaviour in this paper. Conventional analysis ignoring this spatial behaviour could introduce unexpected errors. In order to take the effects of sources on the propagation behaviour into account, a new mathematical model is established for Rayleigh waves in layered elastic media under vertical disc‐like surface sources using the thin‐layer method. The spatial behaviour of the activated modes and the apparent phase velocity, which is the propagation velocity of Rayleigh waves superposed by the multiple modes, are then analysed. Aspects of the spatial behaviour investigated in this paper include the equilibrium path, the particle orbit, and the geometric attenuation of the activated Rayleigh waves. The results presented in this paper can provide some guidelines for developing new inverse mathematical models and algorithms.     

Soil–pile interaction during the passage of rayleigh waves: An analytical solution

A simple analytical solution is presented to calculate the single-pile response when excited by the passage of Rayleigh seismic waves. Closed-form expressions for the horizontal and vertical displacement distributions are presented for piles with finite or infinite length. The analytical results for both free-head and fixed-head piles are obtained through a dynamic Winkler model, with realistic frequency-dependent ‘springs’ and ‘dashpots’. The results of the presented method are in excellent agreement with results of a rigorous solution. It is shown that in vertical motion, the differences between pile and soil displacements are far more significant than in horizontal motion, and therefore, further work is needed to investigate the importance of pile-soil-pile interaction (group effects), because of the vertical component of Rayleigh seismic waves.     

Rocking vibrations of a rigid circular foundation on poroelastic half-space to elastic waves

A study is carefully conducted for the rocking response of a rigid circular foundation resting on a poroelastic half-space when subjected to seismic waves under the framework of Biot’s theory. The free-field waves, rigid-body scattering field waves and radiation scattering field waves are introduced to consider the complex behavior of the soil owing to the scattering phenomena caused by the existence of the foundation. The contact surface between the soil and the foundation is supposed to be perfectly bonded and fully permeable. Combining with the divided wave fields, two sets of dual integral equations elaborating the mixed boundary-value conditions are established, and then reduced to Fredholm integral equations. Therefore, with a semi-analytical method, the expressions of the rocking displacements are obtained. The numerical results of the rocking vibration of the foundation for incident P, SV and Rayleigh waves are presented. The influences of certain parameters, such as the permeability of the soil, the incident angle, Poisson’s ratio and the mass of the foundation, on the rocking vibration of the foundation are explored and studied. Different reactions are found when the foundation is excited by different waves.