A coupled Legendre-Laguerre polynomial method with analytical integration for the Rayleigh waves in a quasicrystal layered half-space with an imperfect interface

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.

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.

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 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.

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

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.

Green's function of multi-layered poroelastic half-space for models of ground vibration due to railway traffic

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.

Dynamic soil-tunnel interaction in layered half-space for incident plane SH waves

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.

Dynamic soil-tunnel interaction in layered half-space for incident P-and SV-waves

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.

Non-axisymmetrical vibration of elastic circular plate on layered transversely isotropic saturated ground

The non-axisymmetrical vibration of elastic circular plate resting on a layered transversely isotropic saturated ground was studied. First, the 3-d dynamic equations in cylindrical coordinate for transversely isotropic saturated soils were transformed into a group of governing differential equations with 1-order by the technique of Fourier expanding with respect to azimuth, and the state equation is established by Hankel integral transform method, furthermore the transfer matrixes within layered media are derived based on the solutions of the state equation. Secondly, by the transfer matrixes, the general solutions of dynamic response for layered transversely isotropic saturated ground excited by an arbitrary harmonic force were established under the boundary conditions, drainage conditions on the surface of ground as well as the contact conditions. Thirdly, the problem was led to a pair of dual integral equations describing the mixed boundaryvalue problem which can be reduced to the Fredholm integral equations of the second kind solved by numerical procedure easily. At the end of this paper, a numerical result concerning vertical and radical displacements both the surface of saturated ground and plate is evaluated.

Surface Motion of Alluvial Valley in Layered Half-Space for Incident Plane P-Waves

The indirect boundary element method (IBEM) is used to study the surface motion of an alluvial valley in layered half-space for incident plane P-waves based on Wolf’s theory. Firstly, the free field response can be solved by the direct stiffness method, and the scattering wave response is calculated by Green’s functions of distributed loads acting on inclined lines in a layered half-space. The method is verified by comparing its results with literature and numerical analyses are performed by taking the amplification of incident plane P-waves by an alluvial valley in one soil layer resting on bedrock as an example. The results show that there exist distinct differences between the wave amplification by an alluvial valley embedded in layered half-space and that in homogeneous half-space and there is interaction between the valley and the soil layer. The amplitudes are relatively large when incident frequencies are close to the soil layer’s resonant frequencies.

Time domain method for calculating free field motion of a layered half-space subjected to obliquely incident body waves

In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.

Layer-element analysis of multilayered saturated soils subject to axisymmetric vertical time-harmonic excitation

The analytical layer-elements for a single poroelastic soil layer and the underlying half-space are established using an algebraic manipulation and Hankel trans- form. According to the boundary conditions and adjacent continuity conditions of general stresses and displacements, a global matrix equation in the transform domain for multi- layered saturated soil media is assembled and solved. Solutions in the frequency domain can be further obtained with an inverse Hankel transform. Numerical examples are used to examine accuracy of the present method and demonstrate effects of soil parameters and load conditions on dynamic responses of the multilayered poroelastic saturated soils.

Dynamic Stress and Pore Pressure Around Circular Cavity in Saturated Porous Elastic Half-space Under Harmonic Plane Dilatational Waves

Dynamic stress concentration and pore pressure concentration around an infinitely long cylindrical cavity of circular cross-section subjected to harmonic plane dilatational waves in fluid-saturated porous elastic half-space were obtained by a complex function method based on potential function and multi-polar coordinate. The steady state Biot’s dynamic field equations of porous elastic solid with a viscous liquid were uncoupled into Helmholtz equations via given potential functions. A circular cavity with large radius is used to replace the straight boundary of the saturated porous elastic half-space. The stresses and pore pressures were obtained by using complex functions in multi-polar coordinates with certain boundary conditions of the solid matrix and the fluid matrix. The approximate solutions were compared to existing numerical solutions. Then the variations of the coefficients of dynamic stress concentration and the pore pressures concentration on boundaries of the cavity were discussed with different parameter conditions. The results of the given numerical example indicate that the method used is useful and efficient to the scattering and dynamic stress concentration of plane dilatational waves in saturated porous elastic half-space.

Solution of a Green's function for a saturated porous medium in a half-space subjected to a torsional force

Based on the solutions of the Green＇s function for a saturated porous medium obtained by the authors, and using transformation of axisymmetric coordinates, Sommerfeld integrals and superposition of the influence field on a free surface, the authors have obtained displacement solutions of a saturated porous medium subjected to a torsional force in a half-space. The relationship curves of the displacement solutions and various parameters （permeability, frequency, etc.） under action of a unit of torque are also given in this paper. The results are consistent with previous Reissner＇s solutions, where a two-phase medium decays to a single-phase medium. The solution is useful in solving relevant dynamic problems of a two- phase saturated medium in engineering.

3-D Scattering of Obliquely Incident Plane P Waves by Alluvial Valley Embedded in Layered Half-Space

The indirect boundary element method(IBEM) was established to solve the problem of 3-D seismic responses of 2-D topographies,by calculating the free-field responses with the direct-stiffness method and simulating the scattering wave fields with the dynamic Green's functions of moving distributed loads.The proposed method yields accurate results,because the 3-D dynamic stiffness matrixes used are exact and the fictitious moving distributed loads can be acted directly on the interface between the alluvial valley and the layered half-space without singularity.The comparison with the published methods verifies the validity of the proposed method.And the numerical analyses are performed to give some beneficial conclusions.The study shows that 3-D scattering by an alluvial valley is essentially different from the 2-D case,and that the presence of soil layer affects not only the amplitude value of surface displacements but also the distribution of surface displacements.

Three-dimensional Scattering of Obliquely Incident Plane SH Waves by an Alluvial Valley in a Layered Half-space

The indirect boundary element method （IBEM） is used to study three-dimensional scattering of obliquely incident plane SH waves by an alluvial valley embedded in a layered half-space. The free-field response of the layered half-space is calculated by the direct stiffness method, and dynamic Green＇s functions of moving distributed loads acting on inclined lines in a layered half-space are calculated to simulate the scattering wave field. The presented method yields very accurate results since the three-dimensional dynamic stiffness matrix is exact and the moving distributed loads can act directly on the valley boundary without singularity. Numerical results and analyses are performed for amplification of obliquely incident plane SH waves around an alluvial valley in a uniform half-space and in single layer over half-space. The results show that the three-dimensional responses are distinctly different from the two-dimensional responses, and the displacement amplitudes around alluvial valleys in a uniform haft-space are obviously different from those in a layered half-space.

Study on Electrical Potential by Buried Source Electrode within the Horizontally Layered Half-space Model

The paper first studies the analytical expression of electrical potential by a point current source in a uniform half-space medium,and then focuses on the distribution of electrical potential in a horizontally layered half-space model by a point current source within the surface layer or the bottom layer.Finally,the electrical potential by a source electrode in any layer of a layered half-space model is presented.

Effects of Solid Matrix and Porosity of Porous Medium on Heat Transfer of Marangoni Boundary Layer Flow Saturated with Power-Law Nanofluids

The effect of the solid matrix and porosity of the porous medium are first introduced to the study of power-law nanofluids, and the Marangoni boundary layer flow with heat generation is investigated. Two cases of solid matrix of porous medium including glass balls and aluminum foam are considered. The governing partial differential equations are simplified by dimensionless variables and similarity transformations, and are solved numerically by using a shooting method with the fourth-fifth-order Runge-Kutta integration technique. It is indicated that the increase of the porosity leads to the enhancement of heat transfer in the surface of the Marangoni boundary layer flow.

3-D DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED SOILS

A study on dynamic response of transversely isotropic saturated poroelastic media under a circular non-axisymmetrical harmonic source has been presented by Huang Yi et al. using the technique of Fourier expansion and Hankel transform. However, the method may not always be valid. The work is extended to the general case being in the rectangular coordinate. The purpose is to study the 3-d dynamic response of transversely isotropic saturated soils under a general source distributing in arbitrary rectangular zoon on the medium surface. Based on Biot＇s theory for fluid- saturated porous media, the 3-d wave motion equations in rectangular coordinate for transversely isotropic saturated poroelastic media were transformed into the two uncoupling governing differential equations of 6-order and 2-order respectively by means of the displacement functions. Then, using the technique of double Fourier transform, the governing differential equations were easily solved. Integral solutions of soil skeleton displacements and pore pressure as well as the total stresses for poroelastic media were obtained. Furthermore, a systematic study on half-space problem in saturated soils was performed. Integral solutions for surface displacements under the general harmonic source distributing on arbitrary surface zone, considering both case of drained surface and undrained surface, were presented.

Complementary Sliding Mode Speed Control with Saturation Function Boundary Layer Optimization

A simple method is proposed to optimize the thickness parameter of the boundary layer of the saturation function in the Complementary Sliding Mode Control (CSMC). A pair of complementary sliding surfaces are constructed. And the Taylor series is used to estimate the steady state error of CSMC system to optimize the parameter value for the boundary layer of the saturation function without artificial settings. This proposed CSMC strategy is applied to the speed regulation in permanent magnet synchronous motor with lump uncertainties. The experimental results show that the proposed CSMC strategy can obtain an optimal value for the boundary layer parameter effectively suppressing the chattering and keep an excellent system performance.

Modeling of Surface Waves in a Fluid Saturated Poro-Elastic Medium under Initial Stress Using Time-Space Domain Higher Order Finite Difference Method

In this present context, mathematical modeling of the propagation of surface waves in a fluid saturated poro-elastic medium under the influence of initial stress has been considered using time dependent higher order finite difference method (FDM). We have proved that the accuracy of this finite-difference scheme is 2M when we use 2nd order time domain finite-difference and 2M-th order space domain finite-difference. It also has been shown that the dispersion curves of Love waves are less dispersed for higher order FDM than of lower order FDM. The effect of initial stress, porosity and anisotropy of the layer in the propagation of Love waves has been studied here. The numerical results have been shown graphically. As a particular case, the phase velocity in a non porous elastic solid layer derived in this paper is in perfect agreement with that of Liu et al. (2009).