Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic half-space (Ⅱ): Numerical results and discussion

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.

Diffraction of plane SV waves by a cavity in poroelastic half-space

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.

Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic half-space (Ⅰ): Formulation

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.

Diffraction of plane P waves around an alluvial valley in poroelastic half-space

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 ILLUSTRATION CALCULATIONS OF SECOND ORDER EFFECTS IN ELASTIC HALF-SPACE ACTED UPON BY A NON-UNIFORM SHEAR LOAD

This paper is a continuation of [1]. An example is discussed in derail to illustrate the second order effects. Numerical calculations for the second order elastic material for the z-direction displacement and the stress t(rz) are carried out. It is found that the second order effect is to reduce z-direction displacement and to decrease t(rz)inside the circle but to increase its value outside the circle.

THE GENERAL SOLUTION OF SECOND ORDER EFFECTS IN ELASTIC HALF-SPACE ACTED UPON BY A NON-UNIFORM SHEAR LOAD

This paper is a continuation of [1]. A closed form solution to the second order elasticity problem, when an isotropic compressible elastic half-space undergoes a deformation owing to a non-uniformly distributed shear load, is presented. The method of integral transform is employed to determine the solutions.

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.

Scattering and diffraction of plane P-waves in a 2-D elastic half-space Ⅱ：shallow arbitrary shaped canyon

Scattering and Diffraction of elastic in-plane P-and SV-waves by a surface topography such as an elastic canyon at the surface of a half-space is a classical problem which has been studied by earthquake engineers and strong-motion seismologists for over forty years. The case of out-of-plane SH waves on the same elastic canyon that is semi-circular in shape on the half-space surface is the first such problem that was solved by analytic closed form solutions over forty years ago by Trifunac. The corresponding case of in-plane P-and SV-waves on the same circular canyon is a much more complicated problem because, the in-plane P-and SV-scattered waves have different wave speeds and together they must have zero normal and shear stresses at the half-space surface. It is not until recently in 2014 that analytic solution for such problem is found by the author in the work of Lee and Liu. This paper uses the technique of Lee and Liu of defining these stress-free scattered waves to solve the problem of the scattered and diffraction of these in-plane waves on an almost-circular surface canyon that is arbitrary in shape.

STEADY-STATE RESPONSE OF A TIMOSHENKO BEAM ON AN ELASTIC HALF-SPACE UNDER A MOVING LOAD

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.

Frequency domain fundamental solutions for a poroelastic half-space

In frequency domain, the fundamental solutions for a poroelastic half-space are re-derived in the context of Biot＇s theory. Based on Biot＇s theory, the governing field equations for the dynamic poroelasicity are established in terms of solid displacement and pore pressure. A method of potentials in cylindrical coordinate system is proposed to decouple the homogeneous Biot＇s wave equations into four scalar Helmholtz equations, and the general solutions to these scalar wave equations are obtained. After that, spectral Green＇s functions for a poroelastic full-space are found through a decomposition of solid displacement, pore pressure, and body force fields. Mirror-image technique is then applied to construct the half-space fundamental solutions.Finally, transient responses of the half-space to buried point forces are examined.

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.

Reflection and transmission of elastic waves at five types of possible interfaces between two dipolar gradient elastic half-spaces

Reflection and transmission of an incident plane wave at five types of possible interfaces between two dipolar gradient elastic solids are studied in this paper. First, the explicit expressions of monopolar tractions and dipolar tractions are derived from the postulated function of strain energy density. Then, the displacements, the normal derivative of displacements, monopolar tractions, and dipolar tractions are used to create the nontraditional interface conditions. There are five types of possible interfaces based on all possible combinations of the displacements and the normal derivative of displacements. These interfacial conditions with consideration of microstructure effects are used to determine the amplitude ratio of the reflection and transmission waves with respect to the incident wave. Further, the energy ratios of the reflection and transmission waves with respect to the incident wave are calculated. Some numerical results of the reflection and transmission coefficients are given in terms of energy flux ratio for five types of possible interfaces. The influences of the five types of possible interfaces on the energy partition between the refection waves and the transmission waves are discussed, and the concept of double channels of energy transfer is first proposed to explain the different influences of five types of interfaces.

Propagation of wave at the boundary surface of transversely isotropic thermoelastic material with voids and isotropic elastic half-space

The purpose of this research is to study the effect of voids on the surface wave propagation in a layer of a transversely isotropic thermoelastic material with voids lying over an isotropic elastic half-space. The frequency equation is derived after developing a mathematical model for welded and smooth contact boundary conditions. The dispersion curves giving the phase velocity and attenuation coefficient via wave number are plotted graphically to depict the effects of voids and anisotropy for welded contact boundary conditions. The specific loss and amplitudes of the volume fraction field, the normal stress, and the temperature change for welded contact are obtained and shown graphically for a particular model to depict the voids and anisotropy effects. Some special cases are also deduced from the present investigation.

A Note on the Effect of Negative Poisson’s Ratio on the Deformation of a Poroelastic Half-Space by Surface Loads

The aim of this note is to study the effect of negative Poisson’s ratio on the quasi-static deformation of a poroelastic half-space with anisotropic permeability and compressible fluid and solid constituents by surface loads. Two particular cases considered are: two-dimensional normal strip loading and axisymmetric normal disc loading. It is found that a negative Poisson’s ratio makes the Mandel-Cryer effect more prominent. It also results in an increase in the magnitude of the surface settlement.

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

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

Microwave-induced thermoacoustic elastic imaging:A simulation study

Microwave-induced thermoacoustic imaging(MTI)has the advantages of high resolution,high contrast,non-ionization,and non-invasive.Recently,MTI was used in the¯eld of breast cancer screening.In this paper,based on the¯nite element method(FEM)and COMSOL Multiphysics software,a three-dimensional breast cancer model suitable for exploring the MTI process is proposed to investigate the in°uence of Young's modulus(YM)of breast cancer tissue on MTI.It is found that the process of electromagnetic heating and initial pressure generation of the entire breast tissue is earlier in time than the thermal expansion process.Besides,compared with normal breast tissue,tumor tissue has a greater temperature rise,displacement,and pressure rise.In particular,YM of the tumor is related to the speed of thermal expansion.In particular,the larger the YM of the tumor is,the higher the heating and contraction frequency is,and the greater the maximum pressure is.Di®erent Young's moduli correspond to di®erent thermoacoustic signal spectra.In MTI,this study can be used to judge di®erent degrees of breast cancer based on elastic imaging.In addition,this study is helpful in exploring the possibility of microwave-induced thermoacoustic elastic imaging(MTAE).

Mathematical framework of nonlinear elastic waves propagating in pre-stressed media

Acoustic nonlinearity holds potential as a method for assessing material stress.Analogous to the acoustoelastic effect,where the velocity of elastic waves is influenced by third-order elastic constants,the propagation of nonlinear acoustic waves in pre-stressed materials would be influenced by higher-order elastic constants.Despite this,there has been a notable absence of research exploring this phenomenon.Consequently,this paper aims to establish a theoretical framework for governing the propagation of nonlinear acoustic waves in pre-stressed materials.It delves into the impact of pre-stress on higher-order material parameters,and specifically examines the propagation of one-dimensional acoustic waves within the contexts of the uniaxial stress and the biaxial stress.This paper establishes a theoretical foundation for exploring the application of nonlinear ultrasonic techniques to measure pre-stress in materials.

Reconfigurable mechanism-based metamaterials for ternary-coded elastic wave polarizers and programmable refraction control

Elastic metamaterials with unusual elastic properties offer unprecedented ways to modulate the polarization and propagation of elastic waves.However,most of them rely on the resonant structural components,and thus are frequency-dependent and unchangeable.Here,we present a reconfigurable 2D mechanism-based metamaterial which possesses transformable and frequency-independent elastic properties.Based on the proposed mechanism-based metamaterial,interesting functionalities,such as ternarycoded elastic wave polarizer and programmable refraction,are demonstrated.Particularly,unique ternary-coded polarizers,with 1-trit polarization filtering and 2-trit polarization separating of longitudinal and transverse waves,are first achieved.Then,the strong anisotropy of the proposed metamaterial is harnessed to realize positive-negative bi-refraction,only-positive refraction,and only-negative refraction.Finally,the wave functions with detailed microstructures are numerically verified.

Publisher Correction to:Highly Elastic,Bioresorbable Polymeric Materials for Stretchable,Transient Electronic Systems

Due to typesetting mistake,Hanul Min was missed to be denoted as a corresponding author in the article.The type-setter apologizes for this.The original article has been corrected.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License,which permits use,sharing,adaptation,distribution and reproduction in any medium or format,as long as you give appropriate credit to the original author(s)and the source,provide a link to the Creative Commons licence,and indicate if changes were made.The images or other third party material in this article are included in the article’s Creative Commons licence,unless indicated otherwise in a credit line to the material.If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use,you will need to obtain permission directly from the copyright holder.

A novel relationship between elastic modulus and void ratio associated with principal stress for coral calcareous sand

Elastic moduli,e.g.shear modulus G and bulk modulus K,are important parameters of geotechnical materials,which are not only the indices for the evaluation of the deformation ability of soils but also the important basic parameters for the development of the constitutive models of geotechnical materials.In this study,a series of triaxial loading-unloading-reloading shear tests and isotropic loading-unloadingreloading tests are conducted to study several typical mechanical properties of coral calcareous sand(CCS),and the void ratio evolution during loading,unloading and reloading.The test results show that the stress-strain curves during multiple unloading processes are almost parallel,and their slopes are much greater than the deformation modulus at the initial stage of loading.The relationship between the confining pressure and the volumetric strain can be defined approximately by a hyperbolic equation under the condition of monotonic loading of confining pressure.Under the condition of confining pressure unloading,the evolution of void ratio is linear in the e-lnp0 plane,and these lines are a series of almost parallel lines if there are multiple processes of unloading.Based on the experimental results,it is found that the modified Hardin formulae for the elastic modulus estimation have a significant deviation from the tested values for CCS.Based on the experimental results,it is proposed that the elastic modulus of soils should be determined by the intersection line of two spatial surfaces in the G/K-e-p’/pa space(pa:atmosphere pressure).“Ye formulation”is further proposed for the estimation of the elastic modulus of CCS.This new estimation formulation for soil elastic modulus would provide a new method to accurately describe the mechanical behavior of granular soils.