Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate ...Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate for actual situations.Although soil exhibits material nonlinearity under earthquake excitation,existing research on the TI medium is limited to soil linearity and neglects the nonlinear response of TI sites.A 2D equivalent linear model for a layered TI half-space subjected to seismic waves is derived in the transformed wave number domain using the exact dynamic stiffness matrix of the TI medium.This study introduces a method for determining the effective shear strain of TI sites under oblique wave incidence,and further describes a systematic study on the effects of TI parameters and soil nonlinearity on site responses.Numerical results indicate that seismic responses of the TI medium significantly differ from those of isotropic sites and that the responses are highly dependent on TI parameters,particularly in nonlinear cases,while also being sensitive to incident angle and excitation intensity.Moreover,the differences in peak acceleration and waveform for various TI materials may also be amplified due to the strong nonlinearity.The study provides valuable insights for improving the accuracy of seismic response analysis in engineering applications.展开更多
The solution of static deformation of two homogeneous, isotropic, and perfectly elastic half-spaces in welded contact due to a long dip-slip fault with uniform slip is well known. The aim of the present article is to ...The solution of static deformation of two homogeneous, isotropic, and perfectly elastic half-spaces in welded contact due to a long dip-slip fault with uniform slip is well known. The aim of the present article is to obtain the static deformation of two homogeneous, isotropic, and perfectly elastic half-spaces in welded contact due to nonuniform slip along a long dip-slip fault. The fault is vertical and lies entirely in the lower half-space, extending up to the interface. Four slip profiles are considered： lin- ear, parabolic, elliptic, and cubic. Closed-form expressions for the elastic residual field have been ob- tained for the different slip profiles. The displacement field due to four nonuniform slip profiles is com- pared with the displacement field due to uniform slip. For comparison, we have assumed the source potency for different slip profiles to be the same, which is achieved by taking the fault slip at the inter- face constant and varying the fault width. It is found that the displacement field varies significantly at the fault width rather than at the interface. Moreover, uniform slip along the fault makes the edges singular. This singularity at the lower edge is not present in the case of nonuniform slip along a long dip-slip fault.展开更多
In this paper,a transversely isotropic piezoelectric half-space with the isotropy axis parallel to the z-axis is considered under rotation on a rigid circular disk bonded to the surface of the piezoelectric medium.Thi...In this paper,a transversely isotropic piezoelectric half-space with the isotropy axis parallel to the z-axis is considered under rotation on a rigid circular disk bonded to the surface of the piezoelectric medium.This is a type of Reissner-Sagoci mixed boundary value problem.By utilizing the Hankel transform,the mixed boundary value problem is simplified into solving a pair of dual integral equations.Full-field analytical expressions for displacement,stresses,and electric displacement inside the half-space are obtained.The shear stresses and electric displacement on the surface are found to be singular at the edge of the rigid circular disk,and the stress intensity factors and electric displacement intensity factor are defined.Numerical results show that material properties and geometric size have significant effects on displacement,shear stresses,and electric displacement.展开更多
In this article, under quite weaker conditions, theorems of Helly’s type for(not necessarily finite) family of partially closed half-spaces are presented in both analytic and geometric forms. Examples are also provid...In this article, under quite weaker conditions, theorems of Helly’s type for(not necessarily finite) family of partially closed half-spaces are presented in both analytic and geometric forms. Examples are also provided to show that these conditions cannot be omitted in general.展开更多
The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when th...The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.展开更多
A solution for the lithospheric deformation of a uniform,elastically isotropic layer(EIL)of uniform thickness welded with a uniform,elastically orthotropic half-space(EOHS)due to a vertical tensile line fault(VTLF)wit...A solution for the lithospheric deformation of a uniform,elastically isotropic layer(EIL)of uniform thickness welded with a uniform,elastically orthotropic half-space(EOHS)due to a vertical tensile line fault(VTLF)with an opening in the horizontal direction located in the isotropic layer is derived in the integral form by employing Airy’s stress function approach for the plane strain case.The linear combination of exponential terms appearing in the denominator of the integral expressions of the deformation field of the EIL is expressed as a finite sum of exponential terms(FSET)by applying the method of least squares to analytically compute the deformation field.The displacement field is discussed in detail and computed numerically by considering the EOHS as olivine or barytes material or considering half-space to be isotropic.展开更多
The wave propagation behavior in an elastic wedge-shaped medium with an arbitrary shaped cylindrical canyon at its vertex has been studied.Numerical computation of the wave displacement field is carried out on and nea...The wave propagation behavior in an elastic wedge-shaped medium with an arbitrary shaped cylindrical canyon at its vertex has been studied.Numerical computation of the wave displacement field is carried out on and near the canyon surfaces using weighted-residuals(moment method).The wave displacement fields are computed by the residual method for the cases of elliptic,circular,rounded-rectangular and flat-elliptic canyons,The analysis demonstrates that the resulting surface displacement depends,as in similar previous analyses,on several factors including,but not limited,to the angle of the wedge,the geometry of the vertex,the frequencies of the incident waves,the angles of incidence,and the material properties of the media.The analysis provides intriguing results that help to explain geophysical observations regarding the amplification of seismic energy as a function of site conditions.展开更多
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 d...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 fi.mctions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The free field displacements are added to the fictitious uniform moving loads induced displacements and the total response is obtained. Numerical calculations are performed for a canyon with homogenous and in one layer over bedrock. The effects of the thickness and stiffness of the layer on the amplification are studied and discussed.展开更多
This paper presents 2.5D scattering of incident plane SH waves by a canyon in layered half-space by the indirect boundary element method (IBEM). The free field response is carried out to give the displacements and str...This paper presents 2.5D scattering of incident plane SH waves by a canyon in layered half-space by the indirect boundary element method (IBEM). The free field response is carried out to give the displacements and stresses on the line which forms boundary of the canyon. The fictitious uniform moving loads are applied to the same line to calculate the Green's functions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The displacements due to the free field and from the fictitious uniform moving loads have to be added to obtain the whole motion. The numerical results are carried out for the cases of a canyon in homogenous and in one layer over bedrock. The results show that the 2.5D wave scattering problem is essentially different from the 2D case, and there exist distinct differences between the wave amplification by a canyon in layered half-space and that in homogeneous half-space. The reasons for the distinct difference are explored, and the effects of the thickness and stiffness of the layer on the amplification are discussed.展开更多
This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffract...This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated （either drained or undrained） cases. There are large phase shift between the dry case and the saturated cases, as well as slightly longer resultant wavelengths for the undrained case than those for the drained case and longer resultant wavelengths for the drained case than those for the dry case. For small porosity the surface displacement amplitudes for the saturated cases are almost identical to those for the dry case; while for large porosity, the effect of drainage condition becomes significant, and the surface displacement amplitudes for the undrained case are larger than those for the drained case. As the incident frequency increases, the effect of porosity becomes significant, and more significant for the undrained case than that for the drained case. As the porosity increases, the pore pressures increase significantly but their oscillations become smoother. As the incident frequency increases, the pore pressures become more complicated.展开更多
This study proposes a Green＇s function, an essential representation of water-saturated ground under moving excitation, to simulate ground borne vibration from trains. First, general solutions to the governing equatio...This study proposes a Green＇s function, an essential representation of water-saturated ground under moving excitation, to simulate ground borne vibration from trains. First, general solutions to the governing equations of poroelastic medium are derived by means of integral transform. Secondly, the transmission and reflection matrix approach is used to formulate the relationship between displacement and stress of the stratified ground, which results in the matrix of the Green＇s function. Then the Green＇s function is combined into a train-track-ground model, and is verified by typical examples and a field test. Additional simulations show that the computed ground vibration attenuates faster in the immediate vicinity of the track than in the surrounding area. The wavelength of wheel-rail unevenness has a notable effect on computed displacement and pore pressure. The variation of vibration intensity with the depth of ground is significantly influenced by the layering of the strata soil. When the train speed is equal to the velocity of the Rayleigh wave, the Mach cone appears in the simulated wave field. The proposed Green＇s function is an appropriate representation for a layered ground with shallow ground water table, and will be helpful to understand the dynamic responses of the ground to complicated moving excitation.展开更多
This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green＇s functions of compressive and shear wave sources are deriv...This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green＇s functions of compressive and shear wave sources are derived based on Biot＇s theory. The scattered waves are constructed using fictitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verification of the accuracy is performed by： （1） checking the satisfaction extent of the boundary conditions, （2） comparing the degenerated solutions of a single-phased case with well- known solutions, and （3） examining the numerical stability of the solutions. The nature of diffraction of plane SV waves around a cavity in a poroelastic half-space is investigated by numerical examples.展开更多
The dynamic soil-tunnel interaction is studied by indirect boundary element method （IBEM）, using the model of a rigid tunnel in layered half-space, which is simplified to a single soil layer on elastic bedrock, subj...The dynamic soil-tunnel interaction is studied by indirect boundary element method （IBEM）, using the model of a rigid tunnel in layered half-space, which is simplified to a single soil layer on elastic bedrock, subjected to incident plane SH waves. The accuracy of the results is verified through comparison with the analytical solution. It is shown that soil-tunnel interaction in layered half-space is larger than that in homogeneous half-space and this interaction mechanism is essentially different from that of soil-foundation-superstructure interaction.展开更多
This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green＇s functions of compressional a...This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green＇s functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot＇s theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.展开更多
In the past two decades numerous studies were made to develop and improve the theory and practical computation techniques of synthesizing theoretical seismograms for the model of multi-layered half-space. Today, synth...In the past two decades numerous studies were made to develop and improve the theory and practical computation techniques of synthesizing theoretical seismograms for the model of multi-layered half-space. Today, synthesizing theoretical seismograms in multi-layered half-space is an important tool for understanding the structure of the Earth’s interior as well as the seismic source process from well-recorded seismograms data. As part of a review of the state-of-the-art, in this article I shall present a systematic and self-contained theory of elastic waves in multi-layered half-space media based on the developments in the past two decades.展开更多
The dynamic stiffness method combined with the Fourier transform is utilized to derive the in-plane Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic（TI）half-s...The dynamic stiffness method combined with the Fourier transform is utilized to derive the in-plane Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic（TI）half-space.The loaded layer is fixed to obtain solutions restricted in it and the corresponding reactions forces,which are then applied to the total system with the opposite sign.By adding solutions restricted in the loaded layer to solutions from the reaction forces,the global solutions in the wavenumber domain are obtained,and the dynamic Green’s functions in the space domain are recovered by the inverse Fourier transform.The presented formulations can be reduced to the isotropic case developed by Wolf（1985）,and are further verified by comparisons with existing solutions in a uniform isotropic as well as a layered TI halfspace subjected to horizontally distributed loads which are special cases of the more general problem addressed.The deduced Green’s functions,in conjunction with boundary element methods,will lead to significant advances in the investigation of a variety of wave scattering,wave radiation and soil-structure interaction problems in a layered TI site.Selected numerical results are given to investigate the influence of material anisotropy,frequency of excitation,inclination angle and layered on the responses of displacement and stress,and some conclusions are drawn.展开更多
Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green＇s fimctions of line sourc...Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green＇s fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.展开更多
The dynamic soil-tunnel interaction is studied by the model of a rigid tunnel embedded in layered half-space, which is simplified as a single soil layer on elastic bedrock to the excitation of P- and SV-waves. The ind...The dynamic soil-tunnel interaction is studied by the model of a rigid tunnel embedded in layered half-space, which is simplified as a single soil layer on elastic bedrock to the excitation of P- and SV-waves. The indirect boundary element method is used, combined with the Green＇ s function of distributed loads acting on inclined lines. It is shown that the dynamic characteristics of soil-tunnel interaction in layered half-space are different much from that in homoge- neous half-space, and that the mechanism of soil-tunnel interaction is also different much from that of soil-founda- tion-superstructure interaction. For oblique incidence, the tunnel response for in-plane incident SV-waves is com- pletely different from that for incident SH-waves, while the tunnel response for vertically incident SV-wave is very similar to that of vertically incident SH-wave.展开更多
By introducing the equivalent stiffness of an elastic half-space interacting with a Timoshenko beam, the displacement solution of the beam resting on an elastic half-space subjected to a moving load is presented. Base...By introducing the equivalent stiffness of an elastic half-space interacting with a Timoshenko beam, the displacement solution of the beam resting on an elastic half-space subjected to a moving load is presented. Based on the relative relation of wave velocities of the half-space and the beam, four cases with the combination of different parameters of the half-space and the beam, the system of soft beam and hard half-space, the system of sub-soft beam and hard half- space, the system of sub-hard beam and soft half-space, and the system of hard beam and soft half-space are considered. The critical velocities of the moving load are studied using dispersion curves. It is found that critical velocities of the moving load on the Timoshenko beam depend on the relative relation of wave velocities of the half-space and the beam. The Rayleigh wave velocity in the half-space is always a critical velocity and the response of the system will be infinite when the load velocity reaches it. For the system of soft beam and hard half-space, wave velocities of the beam are also critical velocities. Besides the shear wave velocity of the beam, there is an additional minimum critical velocity for the system of sub-soft beam and hard half-space. While for systems of （sub-） hard beams and soft half-space, wave velocities of the beam are no longer critical ones. Comparison with the Euler-Bernoulli beam shows that the critical velocities and response of the two types of beams are much different for the system of （sub-） soft beam and hard half-space but are similar to each other for the system of （sub-） hard beam and soft half space. The largest displacement of the beam is almost at the location of the load and the displacement along the beam is almost symmetrical if the load velocity is smaller than the minimum critical velocity （the shear wave velocity of the beam for the system of soft beam and hard half-space）. The largest displacement of the beam shifts behind the load and the asymmetry of the displacement along the beam increases with the increase of the load velocity due to the damping and wave racliation. The displacement of the beam at the front of the load is very small if the load velocity is larger than the largest wave velocity of the beam and the half space. The results of the present study provide attractive theoretical and practical references for the analysis of ground vibration induced by the high-speed train.展开更多
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...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.展开更多
基金National Natural Science Foundation of China under Grant No.U2139208。
文摘Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate for actual situations.Although soil exhibits material nonlinearity under earthquake excitation,existing research on the TI medium is limited to soil linearity and neglects the nonlinear response of TI sites.A 2D equivalent linear model for a layered TI half-space subjected to seismic waves is derived in the transformed wave number domain using the exact dynamic stiffness matrix of the TI medium.This study introduces a method for determining the effective shear strain of TI sites under oblique wave incidence,and further describes a systematic study on the effects of TI parameters and soil nonlinearity on site responses.Numerical results indicate that seismic responses of the TI medium significantly differ from those of isotropic sites and that the responses are highly dependent on TI parameters,particularly in nonlinear cases,while also being sensitive to incident angle and excitation intensity.Moreover,the differences in peak acceleration and waveform for various TI materials may also be amplified due to the strong nonlinearity.The study provides valuable insights for improving the accuracy of seismic response analysis in engineering applications.
基金supported by the Major Research Project and Junior Research Fellowship from the University Grant Commission,New Delhi
文摘The solution of static deformation of two homogeneous, isotropic, and perfectly elastic half-spaces in welded contact due to a long dip-slip fault with uniform slip is well known. The aim of the present article is to obtain the static deformation of two homogeneous, isotropic, and perfectly elastic half-spaces in welded contact due to nonuniform slip along a long dip-slip fault. The fault is vertical and lies entirely in the lower half-space, extending up to the interface. Four slip profiles are considered： lin- ear, parabolic, elliptic, and cubic. Closed-form expressions for the elastic residual field have been ob- tained for the different slip profiles. The displacement field due to four nonuniform slip profiles is com- pared with the displacement field due to uniform slip. For comparison, we have assumed the source potency for different slip profiles to be the same, which is achieved by taking the fault slip at the inter- face constant and varying the fault width. It is found that the displacement field varies significantly at the fault width rather than at the interface. Moreover, uniform slip along the fault makes the edges singular. This singularity at the lower edge is not present in the case of nonuniform slip along a long dip-slip fault.
基金supported by the National Natural Science Foundation of China(11872203)the Creative Research Groups(51921003).
文摘In this paper,a transversely isotropic piezoelectric half-space with the isotropy axis parallel to the z-axis is considered under rotation on a rigid circular disk bonded to the surface of the piezoelectric medium.This is a type of Reissner-Sagoci mixed boundary value problem.By utilizing the Hankel transform,the mixed boundary value problem is simplified into solving a pair of dual integral equations.Full-field analytical expressions for displacement,stresses,and electric displacement inside the half-space are obtained.The shear stresses and electric displacement on the surface are found to be singular at the edge of the rigid circular disk,and the stress intensity factors and electric displacement intensity factor are defined.Numerical results show that material properties and geometric size have significant effects on displacement,shear stresses,and electric displacement.
文摘In this article, under quite weaker conditions, theorems of Helly’s type for(not necessarily finite) family of partially closed half-spaces are presented in both analytic and geometric forms. Examples are also provided to show that these conditions cannot be omitted in general.
基金Project supported by the National Natural Science Foundation of China(No.12102131)the Natural Science Foundation of Henan Province of China(No.242300420248)the International Science and Technology Cooperation Project of Henan Province of China(No.242102521010)。
文摘The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.
文摘A solution for the lithospheric deformation of a uniform,elastically isotropic layer(EIL)of uniform thickness welded with a uniform,elastically orthotropic half-space(EOHS)due to a vertical tensile line fault(VTLF)with an opening in the horizontal direction located in the isotropic layer is derived in the integral form by employing Airy’s stress function approach for the plane strain case.The linear combination of exponential terms appearing in the denominator of the integral expressions of the deformation field of the EIL is expressed as a finite sum of exponential terms(FSET)by applying the method of least squares to analytically compute the deformation field.The displacement field is discussed in detail and computed numerically by considering the EOHS as olivine or barytes material or considering half-space to be isotropic.
文摘The wave propagation behavior in an elastic wedge-shaped medium with an arbitrary shaped cylindrical canyon at its vertex has been studied.Numerical computation of the wave displacement field is carried out on and near the canyon surfaces using weighted-residuals(moment method).The wave displacement fields are computed by the residual method for the cases of elliptic,circular,rounded-rectangular and flat-elliptic canyons,The analysis demonstrates that the resulting surface displacement depends,as in similar previous analyses,on several factors including,but not limited,to the angle of the wedge,the geometry of the vertex,the frequencies of the incident waves,the angles of incidence,and the material properties of the media.The analysis provides intriguing results that help to explain geophysical observations regarding the amplification of seismic energy as a function of site conditions.
基金National Natural Science Foundation of China Under Grant No.50908156 and 50978183Tianjin Natural Science Foundation Under Grant No. 07JCZDJC10100
文摘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 fi.mctions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The free field displacements are added to the fictitious uniform moving loads induced displacements and the total response is obtained. Numerical calculations are performed for a canyon with homogenous and in one layer over bedrock. The effects of the thickness and stiffness of the layer on the amplification are studied and discussed.
基金supported by National Natural Science Foundation of China (Nos. 50908156 and 50978183)
文摘This paper presents 2.5D scattering of incident plane SH waves by a canyon in layered half-space by the indirect boundary element method (IBEM). The free field response is carried out to give the displacements and stresses on the line which forms boundary of the canyon. The fictitious uniform moving loads are applied to the same line to calculate the Green's functions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The displacements due to the free field and from the fictitious uniform moving loads have to be added to obtain the whole motion. The numerical results are carried out for the cases of a canyon in homogenous and in one layer over bedrock. The results show that the 2.5D wave scattering problem is essentially different from the 2D case, and there exist distinct differences between the wave amplification by a canyon in layered half-space and that in homogeneous half-space. The reasons for the distinct difference are explored, and the effects of the thickness and stiffness of the layer on the amplification are discussed.
基金support from the Program for New Century Excellent Talents in University (NCET-05-0248)the Key Program for Applied Basic Research of Tianjin Municipality (07JCZDJC10100)
文摘This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated （either drained or undrained） cases. There are large phase shift between the dry case and the saturated cases, as well as slightly longer resultant wavelengths for the undrained case than those for the drained case and longer resultant wavelengths for the drained case than those for the dry case. For small porosity the surface displacement amplitudes for the saturated cases are almost identical to those for the dry case; while for large porosity, the effect of drainage condition becomes significant, and the surface displacement amplitudes for the undrained case are larger than those for the drained case. As the incident frequency increases, the effect of porosity becomes significant, and more significant for the undrained case than that for the drained case. As the porosity increases, the pore pressures increase significantly but their oscillations become smoother. As the incident frequency increases, the pore pressures become more complicated.
基金National Natural Science Foundation of China Key Project,under Grant No.50538030Postdoctoral Science Foundation of China under Grant No.2013M531084Natural Science Foundation of Heilongjiang Province of China under Grant No.E201221
文摘This study proposes a Green＇s function, an essential representation of water-saturated ground under moving excitation, to simulate ground borne vibration from trains. First, general solutions to the governing equations of poroelastic medium are derived by means of integral transform. Secondly, the transmission and reflection matrix approach is used to formulate the relationship between displacement and stress of the stratified ground, which results in the matrix of the Green＇s function. Then the Green＇s function is combined into a train-track-ground model, and is verified by typical examples and a field test. Additional simulations show that the computed ground vibration attenuates faster in the immediate vicinity of the track than in the surrounding area. The wavelength of wheel-rail unevenness has a notable effect on computed displacement and pore pressure. The variation of vibration intensity with the depth of ground is significantly influenced by the layering of the strata soil. When the train speed is equal to the velocity of the Rayleigh wave, the Mach cone appears in the simulated wave field. The proposed Green＇s function is an appropriate representation for a layered ground with shallow ground water table, and will be helpful to understand the dynamic responses of the ground to complicated moving excitation.
基金Program for New Century Excellent Talents in University Under Grant No. NCET-05-0248the Key Program for Applied Basic Research of Tianjin Municipality Under Grant No. 07JCZDJC10100
文摘This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green＇s functions of compressive and shear wave sources are derived based on Biot＇s theory. The scattered waves are constructed using fictitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verification of the accuracy is performed by： （1） checking the satisfaction extent of the boundary conditions, （2） comparing the degenerated solutions of a single-phased case with well- known solutions, and （3） examining the numerical stability of the solutions. The nature of diffraction of plane SV waves around a cavity in a poroelastic half-space is investigated by numerical examples.
基金National Natural Science Foundation of China under Grant 51378384Key Project of Natural Science Foundation of Tianjin Municipality under Grant 12JCZDJC29000
文摘The dynamic soil-tunnel interaction is studied by indirect boundary element method （IBEM）, using the model of a rigid tunnel in layered half-space, which is simplified to a single soil layer on elastic bedrock, subjected to incident plane SH waves. The accuracy of the results is verified through comparison with the analytical solution. It is shown that soil-tunnel interaction in layered half-space is larger than that in homogeneous half-space and this interaction mechanism is essentially different from that of soil-foundation-superstructure interaction.
基金support from the Program for New Century Excellent Talents in University (NCET-05-0248)the Key Program for Applied Basic Research of Tianjin Municipality (07JCZDJC10100)
文摘This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green＇s functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot＇s theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.
文摘In the past two decades numerous studies were made to develop and improve the theory and practical computation techniques of synthesizing theoretical seismograms for the model of multi-layered half-space. Today, synthesizing theoretical seismograms in multi-layered half-space is an important tool for understanding the structure of the Earth’s interior as well as the seismic source process from well-recorded seismograms data. As part of a review of the state-of-the-art, in this article I shall present a systematic and self-contained theory of elastic waves in multi-layered half-space media based on the developments in the past two decades.
基金National Natural Science Foundation of China under grant No.51578373 and 51578372the Natural Science Foundation of Tianjin Municipality under Grant No.16JCYBJC21600
文摘The dynamic stiffness method combined with the Fourier transform is utilized to derive the in-plane Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic（TI）half-space.The loaded layer is fixed to obtain solutions restricted in it and the corresponding reactions forces,which are then applied to the total system with the opposite sign.By adding solutions restricted in the loaded layer to solutions from the reaction forces,the global solutions in the wavenumber domain are obtained,and the dynamic Green’s functions in the space domain are recovered by the inverse Fourier transform.The presented formulations can be reduced to the isotropic case developed by Wolf（1985）,and are further verified by comparisons with existing solutions in a uniform isotropic as well as a layered TI halfspace subjected to horizontally distributed loads which are special cases of the more general problem addressed.The deduced Green’s functions,in conjunction with boundary element methods,will lead to significant advances in the investigation of a variety of wave scattering,wave radiation and soil-structure interaction problems in a layered TI site.Selected numerical results are given to investigate the influence of material anisotropy,frequency of excitation,inclination angle and layered on the responses of displacement and stress,and some conclusions are drawn.
基金supported by National Natural Science Foundation of China (50978183)
文摘Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green＇s fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.
基金supported by the National Natural Science Foundation of China(No.51378384)the Key Project of Natural Science Foundation of Tianjin Municipality(No. 12JCZDJC29000)
文摘The dynamic soil-tunnel interaction is studied by the model of a rigid tunnel embedded in layered half-space, which is simplified as a single soil layer on elastic bedrock to the excitation of P- and SV-waves. The indirect boundary element method is used, combined with the Green＇ s function of distributed loads acting on inclined lines. It is shown that the dynamic characteristics of soil-tunnel interaction in layered half-space are different much from that in homoge- neous half-space, and that the mechanism of soil-tunnel interaction is also different much from that of soil-founda- tion-superstructure interaction. For oblique incidence, the tunnel response for in-plane incident SV-waves is com- pletely different from that for incident SH-waves, while the tunnel response for vertically incident SV-wave is very similar to that of vertically incident SH-wave.
基金Project supported by the National Natural Science Foundation of China (No.50538010) the Doctoral Education of the State Education Ministry of China (No.20040335083) Encouragement Fund for Young Teachers in University of Ministry of Education.
文摘By introducing the equivalent stiffness of an elastic half-space interacting with a Timoshenko beam, the displacement solution of the beam resting on an elastic half-space subjected to a moving load is presented. Based on the relative relation of wave velocities of the half-space and the beam, four cases with the combination of different parameters of the half-space and the beam, the system of soft beam and hard half-space, the system of sub-soft beam and hard half- space, the system of sub-hard beam and soft half-space, and the system of hard beam and soft half-space are considered. The critical velocities of the moving load are studied using dispersion curves. It is found that critical velocities of the moving load on the Timoshenko beam depend on the relative relation of wave velocities of the half-space and the beam. The Rayleigh wave velocity in the half-space is always a critical velocity and the response of the system will be infinite when the load velocity reaches it. For the system of soft beam and hard half-space, wave velocities of the beam are also critical velocities. Besides the shear wave velocity of the beam, there is an additional minimum critical velocity for the system of sub-soft beam and hard half-space. While for systems of （sub-） hard beams and soft half-space, wave velocities of the beam are no longer critical ones. Comparison with the Euler-Bernoulli beam shows that the critical velocities and response of the two types of beams are much different for the system of （sub-） soft beam and hard half-space but are similar to each other for the system of （sub-） hard beam and soft half space. The largest displacement of the beam is almost at the location of the load and the displacement along the beam is almost symmetrical if the load velocity is smaller than the minimum critical velocity （the shear wave velocity of the beam for the system of soft beam and hard half-space）. The largest displacement of the beam shifts behind the load and the asymmetry of the displacement along the beam increases with the increase of the load velocity due to the damping and wave racliation. The displacement of the beam at the front of the load is very small if the load velocity is larger than the largest wave velocity of the beam and the half space. The results of the present study provide attractive theoretical and practical references for the analysis of ground vibration induced by the high-speed train.
基金supported by the National Natural Science Foundation of China (50908156,50978183)Tianjin Research Program of Application Foundation and Advanced Technology (12JCQNJC04700)
文摘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.
