斜入射平面SH波在层状半空间中沉积谷地周围的三维散射

采用直接刚度法计算自由场动力响应,以层状半空间中移动斜线均布荷载动力格林函数模拟散射波场,采用间接边界元方法求解了层状半空间中沉积谷地对斜入射平面SH波的三维散射问题.由于文中采用的层状场地三维动力刚度矩阵是精确的,且用于模拟散射波场的均布移动斜线荷载可以直接施加在沉积交界面处而不存在奇异性,所以本文方法具有很高精度.文中以均匀半空间和基岩上单一土层中沉积谷地对入射SH波的散射为例进行了数值计算,研究表明,沉积谷地对地震波的三维散射与二维散射之间存在本质差别;层状半空间中沉积谷地与均匀半空间中沉积谷地附近地表位移存在显著差异.     

含断层破碎带场地对平面SH波的放大作用

借助Wolf的土层和半空间精确动力刚度矩阵和斜线荷载格林函数,采用间接边界元法在频域内求解了含断层破碎带场地对SH波的放大作用,着重分析了破碎带较窄断层两侧围岩动力响应的基本规律及场地动力特性对散射的影响.结果表明,断层破碎带对SH波有显著的影响,即使破碎带很窄,也可对SH波产生很大的放大作用;层状场地动力特性对放大作用有显著的影响.     

弹性层状半空间中无限长洞室对斜入射 平面SH波的三维散射(Ⅱ) 数值结果与分析

以基岩上单一土层场地为例, 计算分析了在斜入射平面SH波作用下弹性层状半空间中无限长洞室附近的地表位移. 研究表明, 层状半空间中地下洞室对波的散射与均匀半空间情况存在显著差别. 层状场地由于考虑了场地自身的动力特性, 使得洞室附近地表位移幅值的空间变化更为复杂, 基岩与土层刚度比、 土层厚度对散射效应均有着重要影响. 随着基岩与土层刚度比的增大, 地表位移幅值整体上逐渐增大; 随着土层厚度的增大, 土层对地表位移幅值的影响逐渐减小. 在频域解答的基础上, 给出了层状半空间中洞室对斜入射SH波散射的时域解答, 并以Ricker波为例进行了数值计算.      

弹性层状半空间中无限长洞室对斜入射平面SH波的三维散射(Ⅱ)——数值结果与分析

以基岩上单-土层场地为例,计算分析了在斜入射平面SH波作用下弹性层状半空间中无限长洞室附近的地表位移.研究表明,层状半空间中地下洞室对波的散射与均匀半空间情况存在显著差别.层状场地由于考虑了场地自身的动力特性,使得洞室附近地表位移幅值的空间变化更为复杂,基岩与土层刚度比、土层厚度对散射效应均有着重要影响.随着基岩与土层刚度比的增大,地表位移幅值整体上逐渐增大;随着土层厚度的增大,土层对地表位移幅值的影响逐渐减小.在频域解答的基础上,给出了层状半空间中洞室对斜入射SH波散射的时域解答,并以Ricker波为例进行了数值计算.     

弹性层状半空间中凸起地形对入射平面SH波的放大作用

对Wolf理论进行拓展,使之可解决凸起地形对波的散射问题,进而利用间接边界元法,求解了弹性层状半空间中凸起地形对入射平面SH波的放大作用问题。本文模型的显著特点之一是考虑了层状半空间的动力特性以及层状半空间和凸起地形的阻尼;特点之二是计算精度高。文中以基岩上单一土层中半圆凸起地形对入射平面SH波的放大作用为例进行了数值计算分析。研究表明,基岩上单一土层中凸起地形对入射平面SH波放大作用和均匀半空间中凸起地形有着本质的差别;土层动力特性不仅影响凸起地形地表位移的幅值,还会影响地表位移的频谱;阻尼会显著降低凸起地形对高频波的放大作用。     

横观各向同性层状半空间中凹陷地形对平面SH波的散射

沉积风化作用下土体表现出明显的成层特性和横观各向同性(TI)特性。将更为符合实际的TI介质模型引入到局部地形对地震波的散射问题研究有着重要的理论意义和工程实用价值。本文在TI土层平面外动力刚度矩阵的基础上,推导了斜线均布荷载动力格林函数,进而以该格林函数为基本解建立间接边界元方法研究了TI层状半空间中凹陷地形对平面SH波的散射问题。文中验证了方法的正确性,并分别在频域和时域内进行了数值计算分析。研究表明,TI介质中凹陷地形对SH波的散射与各向同性介质情况存在较大差异,凹陷附近地表位移幅值依赖于SH波的入射角度、入射频率和竖向与水平剪切模量比值;另外土体竖向与水平剪切模量比值对基岩上单一土层场地动力特性有着显著影响,进而使得TI层状场地凹陷附近位移幅值与各向同性情况差异变得更为复杂。     

瑞雷波斜入射下层状半空间中沉积谷地周围的三维散射研究

针对层状半空间中沉积谷地对斜入射瑞雷波的三维散射问题, 采用直接刚度法计算自由场波场, 以层状半空间中移动斜线均布荷载动力格林影响函数求解三维散射波场, 建立了求解该问题的间接边界元方法. 通过与已有结果的比较, 验证了该方法的正确性, 并以均匀半空间以及弹性基岩上单一土层场地中沉积谷地为例进行了计算分析. 研究结果表明: 层状半空间与均匀半空间中沉积谷地对瑞雷波的散射存在显著差别; 层状半空间中瑞雷波的振动模态对沉积附近位移幅值有着重要影响; 土层刚度和厚度等参数也对沉积附近位移幅值大小及空间分布有着显著的影响.      

平面SV波在层状半空间中沉积谷地周围的散射

采用间接边界元方法求解了入射平面SV波在层状半空间中沉积谷地周围的散射问题.问题的解答包含自由场和散射场两部分.自由场可由直接刚度法求得,散射场由层状半空间中斜线荷载动力格林函数来模拟.文中以入射平面SV波在基岩上单一土层中沉积谷地周围散射为例研究了土层和沉积谷地周围的影响.结果表明,由于考虑了土层的动力特性,平面SV...     

弹性层状半空间中沉积谷地对入射平面SH波的放大作用

对Wolf理论进行拓展,使之可解决沉积谷地对波的散射问题,进而利用间接边界元法,求解了弹性层状半空间中沉积谷地对入射平面SH波的放大作用问题。通过自由场反应分析,求得假想边界上各点的位移和各单元的应力反应。施加虚拟均布荷载,求得位移和应力的格林函数。根据应力和位移连续边界条件确定虚拟分布荷载,将自由场反应和虚拟分布荷载产生的反应叠加起来,即得到问题的解答。最后,以基岩上单一土层中沉积谷地对入射平面SH波的放大作用为例进行了数值计算分析,结果表明,层状半空间中沉积谷地和均匀半空间中沉积谷地对波的放大作用存在显著差别。文中分析了造成差别的原因,并讨论了覆盖层厚度和刚度对放大作用的影响。     

层状半空间中洞室对平面SH波的放大作用

利用间接边界元法,求解了弹性层状半空间中洞室对入射平面SH波的放大作用问题,并以基岩上单一土层中洞室对入射平面SH波的放大作用为例进行了数值计算分析。本文模型的特点之一是考虑了层状场地的动力特性,因而更接近于实际工程;特点之二是计算精度非常高。研究表明,层状半空间中洞室对波的放大作用与均匀半空间中情况有着本质的差别;层状半空间中洞室附近地表动力响应由土层动力特性和洞室对波的散射二者共同决定。土层动力特性不仅影响洞室附近地表位移的幅值,还会影响地表位移的频谱。在土层的前几阶共振频率附近,随着基岩与土层剪切波速比的增大,土层的影响随之增大,而随着土层厚度的增加,土层的影响随之减小,并逐渐趋于均匀半空间情况。     

Diffraction of SH-waves by topographic features in a layered transversely isotropic half-space

The scattering of plane SH-waves by topographic features in a layered transversely isotropic (TI) half-space is investigated by using an indirect boundary element method (IBEM). Firstly, the anti-plane dynamic stiffness matrix of the layered TI half-space is established and the free fields are solved by using the direct stiffness method. Then, Green’s functions are derived for uniformly distributed loads acting on an inclined line in a layered TI half-space and the scattered fields are constructed with the deduced Green’s functions. Finally, the free fields are added to the scattered ones to obtain the global dynamic responses. The method is verified by comparing results with the published isotropic ones. Both the steady-state and transient dynamic responses are evaluated and discussed. Numerical results in the frequency domain show that surface motions for the TI media can be significantly different from those for the isotropic case, which are strongly dependent on the anisotropy property, incident angle and incident frequency. Results in the time domain show that the material anisotropy has important effects on the maximum duration and maximum amplitudes of the time histories.     

3D seismic response of a 2D hill-valley staggered topography modeled by a 2.5D multi-domain IBEM

A two and a half dimensional(2.5 D) multidomain indirect boundary element method(IBEM) is developed to study the wave scattering of obliquely incident P-, SVand SH-waves by a hill-valley staggered topography in a multi-layered half-space. The IBEM algorithm includes using2.5 D full-space and half-space Green's functions to construct scattered fields in decomposed closed and opened half-space regions, respectively, and using the dynamic stiffness method to solve the free fields. All regions are finally integrated by introducing the compatibility conditions to obtain the total wave fields. The proposed 2.5 D IBEM has the flexibility in dealing with complex boundaries by directly applying the fictitious loads on the regions' boundaries, with a less storage requirement compared to the full 3 D models. Besides, by combining the specific advantages of the two kinds of Green's functions, the method is well suitable for handling coupled topographies with high accuracy. The method is validated by comparison with published results for a single valley as well as a single hill topography. The effects of height-to-width ratio of hill and layering on dynamic responses are further parametrically investigated by numerical implementations in frequency domain. Results show that the interaction between valley and hills can lead to a more significant amplification within the valley region, and dynamic responses are deeply influenced by the height-towidth of hill, simultaneously depending on incident angle and frequency. Besides, the site effects become more complex when the stratification feature is taken into account.     

2.5D scattering of incident plane SV waves by a canyon in layered half-space

This paper presents a 2.5D scattering of incident plane SV waves by a canyon in a layered half-space by using the indirect boundary element method (IBEM). A free field response analysis is performed to provide the displacements and stresses on the boundary of the canyon where fictitious uniform moving loads are applied to calculate the Green’s functions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The free field displacements are added to the fic...     

Scattering and diffraction of plane SH-waves by periodically distributed canyons

A new method is presented to study the scattering and diffraction of plane SH-waves by periodically distributed canyons in a layered half-space. This method uses the indirect boundary element method combined with Green’s functions of uniformly distributed loads acting on periodically distributed inclined lines. The periodicity feature of the canyons is exploited to limit the discretization effort to a single canyon, which avoids errors induced by the truncation of the infinite boundary, and the computational complexity and the demand on memory can be significantly reduced. Furthermore, the total wave fields are decomposed into the free field and scattered field in the process of calculation, which means that the method has definite physical meaning. The implementation of the method is described in detail and its accuracy is verified. Parametric studies are performed in the frequency domain by taking periodically distributed canyons of semi-circular and semi-elliptic cross-sections as examples. Numerical results show that the dynamic responses of periodically distributed canyons can be quite different from those for a single canyon and significant dynamic interactions exist between the canyons.     

层状半空间中周期分布凸起地形对平面SH波的散射

提出了一种新的以层状半空间中周期分布斜线荷载动力格林函数为基本解的间接边界元方法,研究了周期分布凸起地形对平面SH波的散射问题.方法将散射波场分解为凸起内部散射波场和凸起外部散射波场.凸起内部散射波场通过在凸起闭合边界上施加虚拟斜线荷载产生的动力响应来模拟,而凸起外部散射波场则通过在凸起与半空间交界面上施加虚拟周期分布斜线荷载产生的动力响应来模拟.周期分布斜线荷载动力格林函数的引入,使得本文方法仅需针对一个凸起进行边界单元的离散和求解,便可完成问题的求解,避免了通过截断无限边界求解而引入的误差,方法具有较高精度的同时显著降低了求解自由度.文中通过与已有结果的比较,验证了方法的正确性,并以均匀半空间和基岩上单一土层中周期分布凸起为例进行了数值计算分析.研究表明,凸起间距对凸起地形间的动力相互作用有着显著的影响,同时层状半空间中周期分布凸起地形对SH波的散射与均匀半空间情况也有着显著的差别.     

The scattering of plane P_1 wave by a 3-D sedimentary basin in a fluid-saturated poroelastic half-space using IBEM

The indirect boundary element method(IBEM) is applied to investigate the scattering of elastic waves around a 3-D sedimentary basin filled with fluidsaturated poroelastic medium. Based on this method, the free field and scattered field can be solved according to the boundary conditions. And the numerical accuracy has been verified. The effects of parameters on elastic wave scattering are studied, such as boundary condition, incident frequency,incident angle and porosity of medium. Numerical results illustrate that the amplification effect of surface displacement near poroelastic sedimentary basin is notable. In addition, for the case of large porosity the drainage condition has a significant impact on the response amplitude. Due to the fluid exchange at the interface under the drained condition, the displacement amplitude can be much larger than that under the undrained condition in present study. The study can provide a theoretical basis for the anti-seismic design of engineering structures located in sedimentary basin.     

In-plane soil–structure interaction in layered,fluid-saturated,poroelastic half-space I: Structural response

Linear in-plane soil–structure interaction in two dimensions (2D) is studied in fluid-saturated, poroelastic, layered half-space using the Indirect Boundary Element Method (IBEM). The structure is a shear wall supported by a rigid embedded foundation. Exact stiffness matrices for the soil layer and half-space, and Green׳s functions of uniformly distributed loads and pore pressure on an inclined line are derived. Results of the system response in the frequency domain are presented for the special case of single soil layer over bedrock, semi-circular foundation and zero seepage force. The effects of water saturation, soil porosity, depth of soil layer, rigidity contrast between layer and bedrock are investigated in the frequency domain for incident plane P- and SV waves. The results suggest that water saturation may cause increase of the system frequency by more than 10%.     

Earthquake site effect modeling in the Granada basin using a 3-D indirect boundary element method

In this work, we develop the indirect boundary element method (IBEM) to simulate the seismic site response in a realistic, large-scale 3-D sedimentary basin. Most previous applications of boundary element method have used full-space Green’s functions for wave propagation between element points. We use half-space Green’s functions, which include the seismic wavefield interactions at the free surface and require only the boundary elements of the basin interface. In this way, the size of the matrix equation for solution in the IBEM can be reduced to approximately a quarter of that using full-space Green’s functions. The site response modeling of the Granada basin in southern Spain using the IBEM shows that the basin-induced scattering waves were identified as propagating back and forth inside the basin. The scattered waves also generate surface waves that are weakly propagated outside of the basin. The wave propagation inside and outside of the basin shows different patterns. We observe that the scattered wave is locally amplified, and its propagation direction deviates from that of the incident waves propagation direction. Therefore, the computed seismic response in the basin could provide us with good estimates of the seismic motion.