An IBEM solution to the scattering of plane SH-waves by a lined tunnel in elastic wedge space

The indirect boundary element method （IBEM） is developed to solve the scattering of plane SH-waves by a lined tunnel in elastic wedge space. According to the theory of single-layer potential, the scattered-wave field can be constructed by applying virtual uniform loads on the surface of lined tunnel and the nearby wedge surface. The densities of virtual loads can be solved by establishing equations through the continuity conditions on the interface and zero-traction conditions on free surfaces. The total wave field is obtained by the superposition of free field and scattered-wave field in elastic wedge space. Numerical results indicate that the IBEM can solve the diffraction of elastic wave in elastic wedge space accurately and effi- ciently. The wave motion feature strongly depends on the wedge angle, the angle of incidence, incident frequency, the location of lined tunnel, and material parameters. The waves interference and amplification effect around the tunnel in wedge space is more significant, causing the dynamic stress concentration factor on rigid tunnel and the displacement amplitude of flexible tunnel up to 50.0 and 17.0, respectively, more than double that of the case of half-space. Hence, considerable attention should be paid to seismic resistant or anti-explosion design of the tunnel built on a slope or hillside.

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

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

Antiplane response of two scalene triangular hills and a semi-cylindrical canyon by incident SH-waves

Antiplane response of two scalene triangular hills and a semi-cylindrical canyon by SH-waves is studied using wave function expansion and complex function method. Firstly, the analytical model is divided into three parts, and the displacement solutions of wave fields are constructed based on boundary conditions in the three regions. Three domains are then conjoined to satisfy the ＂conjunction＂ condition at shared boundary. In addition, combined with the zero-stress condition of semi-cylindrical canyon, a series of infinite algebraic equations for the problem are derived. Finally, numerical examples are provided and the influence of different parameters on ground motion is discussed.

Scattering of SH-waves on triangular hill joined bysemi-cylindrical canyon

Scattering of SH-waves on the triangular hill joined by semi-cylindrical canyon in half-space is studied using the method of complex function and moving coordinates. The model being studied is divided into two domains. The wave functions satisfying the required condition at each wedge are constructed in each equation. The equations are solved with Fourier expansion. Numerical results are provided to discuss the influence of scattering of SH-waves.

On New Dispersive SH-Waves Propagating in Piezoelectromagnetic Plate

This theoretical work discovers four new dispersive shear-horizontal (SH) waves propagating in the transversely isotropic piezoelectromagnetic plate of class 6 mm. In this work, the following mechanical, electrical, and magnetic boundary conditions at both the upper and lower free surfaces of the piezoelectromagnetic plate are utilized: the mechanically free surface, continuity of both the electrical and magnetic potentials, and continuity of both the electrical and magnetic inductions. The solutions for the new SH-wave velocities (dispersion relations) are found in explicit forms and then graphically studied. The graphical investigation has soundly illuminated several interesting peculiarities that were also discussed. The piezoelectromagnetic materials, also known as the magnetoelectroelastic media, are famous as smart materials because the electrical subsystem of the materials can interact with the magnetic subsystem via the mechanical subsystem, and vice versa. Therefore, it is very important to know the wave characteristics of such (composite) materials because of possible constitution of new technical devices with a high level of integration. It is obvious that the plate waves can be preferable for further miniaturization of the technical devices and used for the nondestructive testing and evaluation of thin piezoelectromagnetic films.

SCATTERING OF SH-WAVES BY AN INTERACTING INTERFACE LINEAR CRACK AND A CIRCULAR CAVITY NEAR BIMATERIAL INTERFACE

An analytical method is developed for scattering of SH-waves and dynamic stressconcentration by an interacting interface crack and a circular cavity near bimaterial interface.Asuitable Green’s function is contructed,which is the fundamental solution of the displacement fieldfor an elastic half space with a circular cavity impacted by an out-plane harmonic line source loadingat the horizontal surface.First,the bimaterial media is divided into two parts along the horizontalinterface,one is an elastic half space with a circular cavity and the other is a complete half space.Then the problem is solved according to the procedure of combination and by the Green’s functionmethod.The horizontal surfaces of the two half spaces are loaded with undetermined anti-plane forcesin order to satisfy continuity conditions at the linking section,or with some forces to recover cracks bymeans of crack-division technique.A series of Fredholm integral equations of first kind for determiningthe unknown forces can be set up through continuity conditions as expressed in terms of the Green’sfunction.Moreover,some expressions are given in this paper,such as dynamic stress intensity factor(DSIF)at the tip of the interface crack and dynamic stress concentration factor(DSCF)around thecircular cavity edge.Numerical examples are provided to show the influences of the wave numbers,the geometrical location of the interface crack and the circular cavity,and parameter combinations ofdifferent media upon DSIF and DSCF.

Time-history responses on the surface by regularly distributed enormous embedded cavities:Incident SH-waves

The time-history responses of the surface were obtained for a linear elastic half-plane including regularly distributed enormous embedded circular cavities subjected to propagating obliquely incident plane SH-waves. An advanced numerical approach named half-plane time-domain boundary element method（BEM）, which only located the meshes around the cavities, was used to create the model. By establishing the modified boundary integral equation（BIE）independently for each cavity and forming the matrices, the final coupled equation was solved step-by-step in the timedomain to obtain the boundary values. The responses were developed for a half-plane with 512 cavities. The amplification patterns were also obtained to illustrate the frequencydomain responses for some cases. According to the results,the presence of enormous cavities affects the scattering and diffraction of the waves arrived to the surface. The introduced method can be recommended for geotechnical/mechanical engineers to model structures in the fields of earthquake engineering and composite materials.

SCATTERING OF PLANE SH-WAVES ON SEMI-CANYON TOPOGRAPHY OF ARBITRARY SHAPE WITH LINGING IN ANISOTROPIC MEDIA

The purpose of this paper is to use the conforma mapping method[1]to analyzeand evaluate the ground displacement and scattering of incident SH-waves, on thesurface of semi-canyon topography of arbitrary shape with lining in anisotropic media.The problem to be solved can be reduced to the solution of an infinite algebraicequation set by using the method of full-space expansion of Fourier progression Usingthe mapping function and scattering theory to solve problems due to semi-canyon topography with lining is just like mapping the semi-cylindrical canyon of arbitraryshape into a cylindrical canyon in full-space.Moreover,it is far practical inengineering practice.From the computational examples,it is obvious that the variation of displacement amplitudes on the surface near the canyon topography is rather sharp. especially when the freqencies of incident SH-waves increase.

SCATTERING FAR FIELD SOLUTION OF SH-WAVE BY MOVABLE RIGID CYLINDRICAL INTERFACE INCLUSION

The Green's function is used to solve the scattering far fieldsolution of SH-wave by a mov- able rigid cylindrical interfaceinclusion in a linear elastic body. First, a suitable Green'sfunction is devel- oped, which is the fundamental displacementsolution of an elastic half space with a movable rigid half-cylin-drical inclusion impacted by out-of-plane harmonic line source loadedat any point of its horizontal surface.

Investigation of SH-Wave Fundamental Modes in Piezoelectromagnetic Plate: Electrically Closed and Magnetically Closed Boundary Conditions

It is theoretically considered the propagation (first evidence) of new dispersive shear-horizontal (SH) acoustic waves in the piezoelectromagnetic (magnetoelectroelastic) composite plates. The studied two-phase composites (BaTiO3-CoFe2O4 and PZT-5H-Terfenol-D) possess the piezoelectric phase (BaTiO3, PZT-5H) and the piezomagnetic phase (CoFe2O4, Terfenol-D). The mechanical, electrical, and magnetic boundary conditions applied to both the upper and lower free surfaces of the plate are as follows: the mechanically free, electrically closed, and magnetically closed surfaces. As a result, the fundamental modes of two new dispersive SH-waves recently discovered in book [Zakharenko, A.A. (2012) ISBN: 978-3-659-30943-4] were numerically calculated. It was found that for large values of normalized plate thickness kd (k and d are the wavenumber and plate half-thickness, respectively) the velocities of both the new dispersive SH-waves can approach the nondispersive SH-SAW velocity of the piezoelectric exchange surface Melkumyan (PEESM) wave. It was also discussed that for small values of kd, the experimental study of the new dispersive SH-waves can be preferable in comparison with the nondispersive PEESM wave. The obtained results can be constructive for creation of various technical devices based on (non)dispersive SH-waves and two-phase smart materials. The new dispersive SH-waves propagating in the plates can be also employed for nondestructive testing and evaluation. Also, it is obvious that the plates can be used in technical devices instead of the corresponding bulk samples for further miniaturization.

Scattering of SH-wave by multiple circular cavities in half space

In this paper,an analytic method is developed to address steady SH-wave scattering and perform dynamic analysis of multiple circular cavities in half space.The scattered wave function used for scattering of SH-waves by multiple circular cavities,which automatically satisfies the stress-free condition at the horizontal surface,is constructed by applying the symmetry of the SH-wave scattering and the method of multi-polar coordinates system.Applying this scattered wave function and method of moving coordinates,the original problem can be transformed to the problem of SH-wave scattering by multiple circular cavities in the full space.Finally,the solution of the problem can be reduced to a series of algebraic equations and solved numerically by truncating the infinite algebraic equations to the finite ones.Numerical examples are provided for case with two cavities to show the effect of wave number,and the distances between the centers of the cavities and from the centers to the ground surface on the dynamic stress concentration around the cavity impacted by incident steady SH-wave.

Interaction of a circular cavity and a beeline crack in right-angle plane impacted by SH-wave

The problem of scattering of SH-wave by a circular cavity and an arbitrary beeline crack in right-angle plane was investigated using the methods of Green's function,complex variables and muti-polar coordinates.Firstly,we constructed a suitable Green's function,which is an essential solution to the displacement field for the elastic right-angle plane possessing a circular cavity while bearing out-of-plane harmonic line source load at arbitrary point.Secondly,based on the method of crack-division,integration for solution was established,then expressions of displacement and stress were obtained while crack and circular cavities were both in existence.Finally,the dynamic stress concentration factor around the circular cavity and the dynamic stress intensity factor at crack tip were discussed to the cases of different parameters in numerical examples.Calculation results show that the crack produces adverse engineering influence on both of the dynamic stress concentration factor and the dynamic stress intensity factor.

Dynamic stress concentration and scattering of SH-wave by interface elliptic cylindrical cavity

In this paper,the dynamic stress concentration and scattering of SH-waves by bi-material structures that possess an interface elliptic cavity are investigated.First,by using the complex function method,the Green's function is constructed.This yields the solution of the displacement field for an elastic half space with a semi-elliptic canyon impacted by an anti-plane harmonic line source loading on the horizontal surface.Then,the problem is divided into an upper and lower half space along the horizontal interface,regarded as a harmony model.In order to satisfy the integral continuity condition, the unknown anti-plane forces are applied to the interface.The integral equations with unknown forces can be established through the continuity condition,and after transformation,the algebraic equations are solved numerically.Finally,the distribution of the dynamic stress concentration factor(DSCF)around the elliptic cavity is given and the effect of different parameters on DSCF is discussed.

SCATTERING OF CIRCULAR CAVITY IN RIGHT-ANGLE PLANAR SPACE TO STEADY SH-WAVE

Complex function method and multi-polar coordinate transformation technology are used here to study scattering of circular cavity in right-angle planar space to SH-wave with out-of-plane loading on the horizontal straight boundary. At first, Green function of right-angle planar space which has no circular cavity is constructed; then the scattering solution which satisfies the free stress conditions of the two right-angle boundaries with the circular cavity existing in the space is formulated. Therefore, the total displacement field can be constructed using overlapping principle. An infinite algebraic equations of unknown coefficients existing in the scattering solution field can be gained using multi-polar coordinate and the free stress condition at the boundary of the circular cavity. It can be solved by using limit items in the infinite series which can give a high computation precision. An example is given to illustrate the variations of the tangential stress at the boundary of the circular cavity due to different dimensionless wave numbers, the location of the circular cavity, the loading center and the distributing range of the out-of-plane loading. The results show the efficiency and effectiveness of the mothod introduced here.

SCATTERING OF PLANE SH-WAVE BY A CYLINDRICAL HILL OF ARBITRARY SHAPE

The problems of scattering of plane SH-wave by a cylindrical hill of arbitrary shape is studied based on the methods of conjunction and division of solution zone. The scattering wave function is given by using the complex variable and conformal mapping methods. The conjunction boundary conditions are satisfied. Furthermore appling orthogonal function expanding technique, the problems can finally be summarized into the solution of a series of infinite algebraic equations. At last, numerical results of surface displacements of a cylindrical arc hill and of a semi-ellipse hill are obtained. And those computational results are compared with the results of finite element method (FEM).

Scattering of SH-wave from interface cylindrical elastic inclusion with a semicircular disconnected curve

Scattering of SH wave from an interface cylindrical elastic inclusion with a semicircular disconnected curve is investigated. The solution of dynamic stress concentration factor is given using the Green＇s function and the method of complex variable functions. First, the space is divided into upper and lower parts along the interface. In the lower half space, a suitable Green＇s function for the problem is constructed. It is an essential solution of the displacement field for an elastic half space with a semi-cylindrical hill of cylindrical elastic inclusion while bearing out-plane harmonic line source load at the horizontal surface. Thus, the semicircular disconnected curve can be constructed when the two parts are bonded and continuous on the interface loading the undetermined anti-plane forces on the horizontal surfaces. Also, the expressions of displacement and stress fields are obtained in this situation. Finally, examples and results of dynamic stress concentration factor are given. Influences of the cylindrical inclusion and the difference parameters of the two mediators are discussed.

Global SH-wave propagation in a 2D whole Moon model using the parallel hybrid PSM/FDM method

We present numerical modeling of SH-wave propagation for the recently proposed whole Moon model and try to improve our understanding of lunar seismic wave propagation. We use a hybrid PSM/FDM method on staggered grids to solve the wave equations and implement the calculation on a parallel PC cluster to improve the computing efficiency. Features of global SH-wave propagation are firstly discussed for a 100-km shallow and900-km deep moonquakes, respectively. Effects of frequency range and lateral variation of crust thickness are then investigated with various models. Our synthetic waveforms are finally compared with observed Apollo data to show the features of wave propagation that were produced by our model and those not reproduced by our models. Our numerical modeling show that the low-velocity upper crust plays significant role in the development of reverberating wave trains. Increasing frequency enhances the strength and duration of the reverberations.Surface multiples dominate wavefields for shallow event.Core–mantle reflections can be clearly identified for deep event at low frequency. The layered whole Moon model and the low-velocity upper crust produce the reverberating wave trains following each phases consistent with observation. However, more realistic Moon model should be considered in order to explain the strong and slow decay scattering between various phases shown on observation data.

Surface Effects on the Scatter of SH-Wave by a Shallow Buried the Elliptical Hole

The methods of complex function, multi-polar coordinate system, and conformal mapping are used to solve dynamic stress concentration factor. The surface elasticity theory is applied to obtain the stress boundary conditions on the surface. The effects of frequency and the ration of the major and minor axis of the ellipse on the dynamic stress concentration factor around the elliptical nano-hole are discussed in detail. When the size of elliptical hole shrinks to nanometers, the numerical results show that the surface effect has a significant effect on the scattering of SH-wave.

Global SH-wavefield calculation for a two-dimensional whole-Earth model with the parallel hybrid PSM/FDM algorithm

We present a parallel hybrid algorithm based on pseudospectral method （PSM） and finite difference method （FDM） for two-dimensional （2-D） global SH- wavefield simulation. The whole-Earth model is taken as a cross section of spherical Earth, and corresponding wave equations are defined in 2-D cylindrical coordinates. Spatial derivatives in the wave equations are approximated with efficient and high accuracy PSM in the lateral and high-order FDM in the radial direction on staggered grids. This algorithm allows us to divide the whole-Earth into sub-domains in radial direction and implement efficient parallel computing on PC cluster, while retains high accuracy and efficiency of PSM in lateral direction. A transformation of moment tensor between 3-D spherical Earth and our 2-D model was proposed to give corre- sponding moment tensor components used in 2-D modeling. Comparison of modeling results with those obtained by direct solution method shows very good accuracy of our algorithm. We also demonstrate its feasibility with a lateral heterogeneous whole-Earth model with localized velocity perturbation.

Measurement of Surface SH-Wave Velocities Generated on the Surface of Japanese Cypress Column

Polyetherimide resin wedge transducers were used to generate a shear wave that was obliquely incident relative to the surface of a Japanese cypress column for measuring the surface SH-wave velocity. As the inter-transducer distance increased, the propagation time increased and the am-plitude became smaller. The propagation time and the amplitude were significantly correlated with the inter-transducer distance. The SH-wave velocity ranged from 1270 m/s to 1496 m/s. Surface SH-wave velocity was lower in the central part of the column and higher in the outer part. Velocity was negatively correlated with moisture content at 1% of significance level. These results suggest the accomplishment of the first target for applying the surface SH-wave acoustoelastic technique to nondestructive evaluation of drying stress in wood.