Amplification of in-plane seismic ground motion by group cavities in layered half-space (Ⅱ): with saturated poroelastic soil layers

As the continuation study on amplification of in-plane seismic ground motion by underground group cavities in layered half-space, this study extends to the case of poroelastic half-space with dry poroelastic and saturated poroelastic soil layers. The influence of poroelastic layers on the amplification of seismic ground motion is studied both in frequency domain and time domain using indirect boundary element method （IBEM）. It is shown that for the example of a saturated poroelastic site in Tianjin under the excitation of Taft wave and E1 Centro wave, the amplification of seismic ground motion in poroelastic case is slightly smaller than that in the elastic case, and the amplification of PGA （peak ground acceleration） and its PRS （peak response spectrum）.. can be increased up to 38.8% and 64.6%; the predominant period of response spectra in poroelastic case becomes shorter to some extent compared with that in the elastic case. It is suggested that the effect of underground group cavities in poroelastic half-space on design seismic ground motion should be considered.

Seismic analysis of high-speed railway irregular bridge–track system considering V-shaped canyon effect

To explore the effect of canyon topography on the seismic response of railway irregular bridge-track system that crosses a V-shaped canyon, seismic ground motions of the horizontal site and V-shaped canyon site were simulated through theoretical analysis with 12 earthquake records selected from the Pacific Earthquake Engineering Research Center(PEER) Strong Ground Motion Database matching the site condition of the bridge.Nonlinear seismic response analyses of an existing 11-span irregular simply supported railway bridge-track system were performed under the simulated spatially varying ground motions. The effects of the V-shaped canyon topography on the peak ground acceleration at bridge foundations and seismic responses of the bridge-track system were analyzed. Comparisons between the results of horizontal and V-shaped canyon sites show that the top relative displacement between adjacent piers at the junction of the incident side and the back side of the V-shaped site is almost two times that of the horizontal site, which also determines the seismic response of the fastener. The maximum displacement of the fastener occurs in the V-shaped canyon site and is 1.4 times larger than that in the horizontal site. Neglecting the effect of V-shaped canyon leads to the inappropriate assessment of the maximum seismic response of the irregular high-speed railway bridge-track system. Moreover, engineers should focus on the girder end to the left or right of the two fasteners within the distance of track seismic damage.

3D scattering of obliquely incident plane SV waves by an alluvial valley embedded in a fluid-saturated, poroelastic layered half-space

The indirect boundary element method is used to study the 3D dynamic response of an infinitely long alluvial valley embedded in a saturated layered half-space for obli- quely incident SV waves. A wave-number transform is first applied along the valley＇s axis to reduce a 3D problem to a 2D plane strain problem. The problem is then solved in the section perpendicular to the axis of the valley. Finally, the 3D dynamic responses of the valley are obtained by an inverse wave-number transform. The validity of the method is con- firmed by comparison with relevant results. The differences between the responses around the valley embedded in dry and in saturated poroelastic medium are studied, and the effects of drainage conditions, porosity, soil layer stiffness, and soil layer thickness on the dynamic response are dis- cussed in detail resulting in some conclusions.

Amplification of in-plane seismic ground motion by group cavities in layered half-space (Ⅰ)

Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method （IBEM）, and the effect of cavity interval and spectrum of incident waves on the amplification are studied by numerical examples. It is shown that there may be large interaction between cavities, and group cavities with certain intervals may have significant amplification to seismic ground motion. The amplification of PGA （peak ground acceleration） and its PRS （peak response spectrum） can be increased up to 45.2% and 84.4%, for an example site in Tianjin, under the excitation of Taft wave and E1 Centro wave; and group cavities may also affect the spectra of the seismic ground motion. It is suggested that the effect of underground group cavities on design seismic ground motion should be considered.

3D Diffraction of obliquely incident SH waves by twin infinitely long cylindrical cavities in layered poroelastic half-space

Abstract This paper studies three-dimensional diffraction of obliquely incident plane SH waves by twin infinitely long cylindrical cavities in layered poroelastic half-space using indirect boundary element method. The approach is validated by comparison with the literature, and the effects of cavity interval, incident frequency, and boundary drainage condition on the diffraction are studied through numerical examples. It is shown that, the interaction between two cavities is significant and surface displacement peaks become large when two cavities are close, and the surface displacement may be significantly amplified by twin cavities, and the influence range with large amplification can be as wide as 40 times of the cavity radius. Surface displacements in dry poroelastic case and saturated poroelastic cases with drained and undrained boundaries are evidently different under certain circumstances, and the differences may be much larger than those in the free-field response.

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.5D)multi-domain indirect boundary element method(IBEM)is developed to study the wave scattering of obliquely incident P-,SV-and SH-waves by a hill-valley staggered topography in a multi-layered half-space.The IBEM algorithm includes using 2.5D 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.5D 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 3D 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-to-width of hill,simultaneously depending on incident angle and frequency.Besides,the site effects become more complex when the stratification feature is taken into account.

An end-to-end 3D seismic simulation of underground structures due to point dislocation source by using an FK-FEM hybrid approach

Based on the domain reduction idea and artificial boundary substructure method,this paper proposes an FK-FEM hybrid approach by integrating the advantages of FK and FEM(i.e.,FK can efficiently generate high-frequency three translational motion,while FEM has rich elements types and constitutive models).An advantage of this approach is that it realizes the entire process simulation from point dislocation source to underground structure.Compared with the plane wave field input method,the FK-FEM hybrid approach can reflect the spatial variability of seismic motion and the influence of source and propagation path.This approach can provide an effective solution for seismic analysis of underground structures under scenario of earthquake in regions where strong earthquakes may occur but are not recorded,especially when active faults,crustal,and soil parameters are available.Taking Daikai subway station as an example,the seismic response of the underground structure is simulated after verifying the correctness of the approach and the effects of crustal velocity structure and source parameters on the seismic response of Daikai station are discussed.In this example,the influence of velocity structure on the maximum interlayer displacement angle of underground structure is 96.5%and the change of source parameters can lead to the change of structural failure direction.