This paper presents an analogical study between electromagnetic and elastic wave fields,with a one-to-one correspondence principle established regarding the basic wave equations,the physical quantities and the differe...This paper presents an analogical study between electromagnetic and elastic wave fields,with a one-to-one correspondence principle established regarding the basic wave equations,the physical quantities and the differential operations.Using the electromagnetic-to-elastic substitution,the analogous relations of the conservation laws of energy and momentum are investigated between these two physical fields.Moreover,the energy-based and momentum-based reciprocity theorems for an elastic wave are also derived in the time-harmonic state,which describe the interaction between two elastic wave systems from the perspectives of energy and momentum,respectively.The theoretical results obtained in this analysis can not only improve our understanding of the similarities of these two linear systems,but also find potential applications in relevant fields such as medical imaging,non-destructive evaluation,acoustic microscopy,seismology and exploratory geophysics.展开更多
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...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.展开更多
Incident particles in the Klein tunnel phenomenon in quantum mechanics can pass a very high potential barrier.Introducing the concept of tunneling into the analysis of phononic crystals can broaden the application pro...Incident particles in the Klein tunnel phenomenon in quantum mechanics can pass a very high potential barrier.Introducing the concept of tunneling into the analysis of phononic crystals can broaden the application prospects.In this study,the structure of the unit cell is designed,and the low frequency(<1 k Hz)valley locked waveguide is realized through the creation of a phononic crystal plate with a topological phase transition interface.The defect immunity of the topological waveguide is verified,that is,the wave can propagate along the original path in the cases of impurities and disorder.Then,the tunneling phenomenon is introduced into the topological valley-locked waveguide to analyze the wave propagation,and its potential applications(such as signal separators and logic gates)are further explored by designing phononic crystal plates.This research has broad application prospects in information processing and vibration control,and potential applications in other directions are also worth exploring.展开更多
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...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.展开更多
The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-eleme...The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.展开更多
In this work,the three-dimensional(3 D)propagation behaviors in the nonlinear phononic crystal and elastic wave metamaterial with initial stresses are investigated.The analytical solutions of the fundamental wave and ...In this work,the three-dimensional(3 D)propagation behaviors in the nonlinear phononic crystal and elastic wave metamaterial with initial stresses are investigated.The analytical solutions of the fundamental wave and second harmonic with the quasilongitudinal(qP)and quasi-shear(qS_(1) and qS_(2))modes are derived.Based on the transfer and stiffness matrices,band gaps with initial stresses are obtained by the Bloch theorem.The transmission coefficients are calculated to support the band gap property,and the tunability of the nonreciprocal transmission by the initial stress is discussed.This work is expected to provide a way to tune the nonreciprocal transmission with vector characteristics.展开更多
The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data...The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.展开更多
In this study,the wave propagation properties of lattice metamaterials with Koch fractal structures are investigated in terms of band structures and directional wave propagation.The analytical models of lattice metama...In this study,the wave propagation properties of lattice metamaterials with Koch fractal structures are investigated in terms of band structures and directional wave propagation.The analytical models of lattice metamaterials are established using the finite element method,and the dispersion relation is solved using the Bloch’s theorem.The band structures of the lattice metamaterials with different numbers of iterations are studied,and the group velocities at a selected frequency are calculated to analyze the directional wave propagation characteristics.Furthermore,dynamic responses of the finite structures are calculated using commercial finite element software to verify the band gaps and directional wave propagation behaviors in the lattice metamaterials.The results show that multiple and low band gaps are present in the lattice materials with various geometric parameters of the Koch fractal,and the position of the lowest band gap decreases as the number of iterations increases.The results indicate the potential applications of lattice metamaterials with Koch fractals for vibration isolation and multi-functional design.展开更多
This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial...This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.展开更多
Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coeff...Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.展开更多
It is an important subject to probe the structure in the medium by various kinds of detection methods in the geotechnical engineering. Based on the propagation theory of elastic wave in half-space layered medium, the ...It is an important subject to probe the structure in the medium by various kinds of detection methods in the geotechnical engineering. Based on the propagation theory of elastic wave in half-space layered medium, the propagation characteristics of elastic wave in layered medium with different elastic parameters are discussed using dynamic analysis of finite element method. It is known that the S-wave velocity, density and thickness of layer are related to the properties of the elastic wave including waveform characteristics, spectral characteristics and time-frequency characteristics. We pay special attention to the structure with low velocity interlayer. The impact imaging method is applied to the grouting construction of the immersed tube tunnel. Data acquisition and analytical method are introduced in detail. The grouting effects can be qualitatively evaluated by comparing the characteristics of elastic wave before grouting with those after grouting. Finally, a quantitative evaluation is obtained according to the relationship between energy response of elastic wave and impedance ratio.展开更多
Due to their potential properties unlike traditional materials and structures,elastic wave metamaterials have received significant interests in recent years.With the coupling between the acoustic and vibration,their m...Due to their potential properties unlike traditional materials and structures,elastic wave metamaterials have received significant interests in recent years.With the coupling between the acoustic and vibration,their mechanical characteristics can be tuned by the active feedback control system at low frequency ranges in which the traditional passive control is limited.This work illustrates that the superior performances of the effective mass density and sound pressure level(SPL)of an elastic wave metamaterial can be significantly changed by the active control,in which the periodic array of local resonators and orthogonal stiffeners are included.Significantly,based on the locally resonant mechanism,the negative density occurs over a frequency range.Due to the effects of lattice constant,structural damping and other parameters,the SPL with the function of fluid-solid coupling are illustrated and discussed.展开更多
Mesh-free finite difference(FD)methods can improve the geometric flexibility of modeling without the need for lattice mapping or complex meshing process.Radial-basisfunction-generated FD is among the most commonly use...Mesh-free finite difference(FD)methods can improve the geometric flexibility of modeling without the need for lattice mapping or complex meshing process.Radial-basisfunction-generated FD is among the most commonly used mesh-free FD methods and can accurately simulate seismic wave propagation in the non-rectangular computational domain.In this paper,we propose a perfectly matched layer(PML)boundary condition for a meshfree FD solution of the elastic wave equation,which can be applied to the boundaries of the non-rectangular velocity model.The performance of the PML is,however,severely reduced for near-grazing incident waves and low-frequency waves.We thus also propose the complexfrequency-shifted PML(CFS-PML)boundary condition for a mesh-free FD solution of the elastic wave equation.For two PML boundary conditions,we derive unsplit time-domain expressions by constructing auxiliary differential equations,both of which require less memory and are easy for programming.Numerical experiments demonstrate that these two PML boundary conditions effectively eliminate artificial boundary reflections in mesh-free FD simulations.When compared with the PML boundary condition,the CFS-PML boundary condition results in better absorption for near-grazing incident waves and evanescent waves.展开更多
A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding cr...A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding criterion. Based on the HLLCE, a third-order cell-centered Lagrangian scheme is built for one-dimensional elastic-plastic problems. A number of numerical experiments are carried out. The numerical results show that the proposed third-order scheme achieves the desired order of accuracy. The third-order scheme is used to the numerical solution of the problems with elastic shock waves and elastic rarefaction waves. The numerical results are compared with a reference solution and the results obtained by other authors. The comparison shows that the pre- sented high-order scheme is convergent, stable, and essentially non-oscillatory. Moreover, the HLLCE is more efficient than the two-rarefaction Riemann solver with elastic waves （TRRSE）展开更多
Using the active feedback control system on the elastic wave metamaterial,this research concentrates on the sound transmission with the dynamic effective model.The metamaterial is subjected to an incident pressure and...Using the active feedback control system on the elastic wave metamaterial,this research concentrates on the sound transmission with the dynamic effective model.The metamaterial is subjected to an incident pressure and immersed in the external mean flow.The elastic wave metamaterial consists of double plates and the upper and lower four-link mechanisms are attached inside.The vertical resonator is attached by the active feedback control system and connected with two four-link mechanisms.Based on the dynamic equivalent method,the metamaterial is equivalent as a single-layer plate by the dynamic effective parameter.With the coupling between the fluid and structure,the expression of the sound transmission loss(STL)is derived.This research shows the influence of effective mass density on sound transmission properties,and the STL in both modes can be tuned by the acceleration and displacement feedback constants.In addition,the dynamic response and the STL are also changed obviously by different values of structural damping,incident angle(i.e.,the elevation and azimuth angles)and Mach number of the external fluid with the mean flow property.The results for sound transmission by two methods are compared,i.e.,the virtual work principle for double plates and the dynamic equivalent method corresponding to a single one.This paper is expected to be helpful for understanding the sound transmission properties of both pure single-and double-plate models.展开更多
Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condit...Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.展开更多
To understand the evolution of stress-induced elastic wave anisotropy,three triaxial experiments were performed on sandstone specimens with bedding orientations parallel,perpendicular,and oblique to the maximum princi...To understand the evolution of stress-induced elastic wave anisotropy,three triaxial experiments were performed on sandstone specimens with bedding orientations parallel,perpendicular,and oblique to the maximum principal stress.P-wave velocities along 64 different directions on each specimen were monitored frequently to understand the anisotropy change at various stress levels by fitting Thomsen’s anisotropy equation.The results show that the elastic wave anisotropy is very sensitive to mechanical loading.Under hydrostatic loading,the magnitude of anisotropy is reduced in all three specimens.However,under deviatoric stress loading,the evolution of anisotropic characteristics(magnitude and orientation of the symmetry axis)is bedding orientation dependent.Anisotropy reversal occurs in specimens with bedding normal/oblique to the maximum principal stress.P-wave anisotropyε0 is linearly related to volumetric strain Sv and dilatancy,indicating that stress-induced redistribution of microcracks has a significant effect on P-wave velocity anisotropy.The closure of initial cracks and pores aligned in the bedding direction contributes to the decrease of the anisotropy.However,opening of new cracks,aligned in the maximum principal direction,accounts for the increase of the anisotropy.The experimental results provide some insights into the microstructural behavior under loading and provide an experimental basis for seismic data interpretation and parameter selection in engineering applications.展开更多
It is essential to assess the evolution of soil fabric as it has an important role in the mechanical responses of soils during complex loading conditions.This contribution carries out the physical experiments using th...It is essential to assess the evolution of soil fabric as it has an important role in the mechanical responses of soils during complex loading conditions.This contribution carries out the physical experiments using three granular materials in the laboratory.The variations of compression and shear wave velocities(Vp and Vs)are investigated during load-unload cycles under dry and drained conditions.Supplementary discrete element method(DEM)simulations are performed to understand the evolution of soil fabric during the equivalent load-unload cycles using spherical particles.Vp and Vs are not always reversible even though the stress state regains its isotropic condition after unload,indicating that Vp and Vs are governed by not only the stress state but also the fabric change.The variations of Vp/Vs are density-and stress-dependent;a higher level of stress ratio(s01/s03)threshold is observed for denser packings to trigger a significant change in wave velocity ratio(Vp/Vs)for experimental results using spherical glass beads and simulation data using spherical particles.Considering the particle shape,a higher s01/s03 threshold is found for more angular particles than rounded particles.The DEM result reveals that Vp/Vs of spherical particles can be correlated linearly with the evolution of fabric ratio(Fver/Fhor)during loadunload in a pre-peak range under dry and drained conditions.展开更多
Using a stiffness matrix method, we in- vestigate the propagation behaviors of elastic waves in one-dimensional （1D） piezoelectric/piezomagnetic （PE/PM） phononic crystals （PCs） with line defects by calculating e...Using a stiffness matrix method, we in- vestigate the propagation behaviors of elastic waves in one-dimensional （1D） piezoelectric/piezomagnetic （PE/PM） phononic crystals （PCs） with line defects by calculating energy reflection/transmittion coefficients of quasi-pressure and quasi-shear waves. Line defects are created by the re- placement of PE or PM constituent layer. The defect modes existing in the first gap are considered and the influences on defect modes of the material properties and volume fraction of the defect layers, the type of incident waves, the location of defect layer and the number of structural layers are discussed in detail. Numerical results indicate that defect modes are the most obvious when the defect layers are inserted in the middle of the perfect PCs; the types of incidence wave and material properties of the defect layers have important effects on the numbers, the location of frequencies and the peaks of defect modes, and the defect modes are strongly de- pendent on volume fraction of the defect layers. We hope this paper will be found useful for the design of PE/PM acoustic filters or acoustic transducer with PCs structures.展开更多
Basedonthe first-order Biot-equation with simplified coefficients,astaggered irregu- lar-grid finite difference method(FDM)is developed to simulate elastic wave propagation in 3-D heterogeneous anisotropic porous medi...Basedonthe first-order Biot-equation with simplified coefficients,astaggered irregu- lar-grid finite difference method(FDM)is developed to simulate elastic wave propagation in 3-D heterogeneous anisotropic porous media.The ‘slow’P wave in porous media wave simulation is highly dispersive.Finer grids are needed to get a precise wavefield calculation for models with curved interface and complex geometric structure.Fine grids in a global model greatly increase computation costs of regular grids scheme.Irregular fine or coarse grids in local fields not only cost less computing time than the conventional velocity-stress FDM,but also give a more accu- rate wavefield description.A dispersion analysis of the irregular-grid finite difference operator has confirmed the stability and high efficiency.The absorbing boundary condition is used to elimi- nate artificial reflections.Numerical examples show that this new irregular-grid finite difference method is of higher performance than conventional methods using regular rectangular grids in simulating elastic wave propagation in heterogeneous anisotropic porous media.展开更多
基金funded by the National Natural Science Foundation of China(Grant No.12404507)the Natural Science Research of Jiangsu Higher Education Institutions of China(Grant No.24KJB140013)the Scientific Startup Foundation of Nanjing Normal University(Grant No.184080H201B49).
文摘This paper presents an analogical study between electromagnetic and elastic wave fields,with a one-to-one correspondence principle established regarding the basic wave equations,the physical quantities and the differential operations.Using the electromagnetic-to-elastic substitution,the analogous relations of the conservation laws of energy and momentum are investigated between these two physical fields.Moreover,the energy-based and momentum-based reciprocity theorems for an elastic wave are also derived in the time-harmonic state,which describe the interaction between two elastic wave systems from the perspectives of energy and momentum,respectively.The theoretical results obtained in this analysis can not only improve our understanding of the similarities of these two linear systems,but also find potential applications in relevant fields such as medical imaging,non-destructive evaluation,acoustic microscopy,seismology and exploratory geophysics.
基金supported by the National Key R&D Program of China(No.2021YFE0110900)the National Natural Science Foundation of China(Nos.U22B2078 and 11991033)。
文摘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.
基金supported by the National Natural Science Foundation of China(No.12172297)the Open Foundation of State Key Laboratory of Structural Analysis for Industrial Equipment of China(No.GZ22106)the Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University of China(No.CX2023055)。
文摘Incident particles in the Klein tunnel phenomenon in quantum mechanics can pass a very high potential barrier.Introducing the concept of tunneling into the analysis of phononic crystals can broaden the application prospects.In this study,the structure of the unit cell is designed,and the low frequency(<1 k Hz)valley locked waveguide is realized through the creation of a phononic crystal plate with a topological phase transition interface.The defect immunity of the topological waveguide is verified,that is,the wave can propagate along the original path in the cases of impurities and disorder.Then,the tunneling phenomenon is introduced into the topological valley-locked waveguide to analyze the wave propagation,and its potential applications(such as signal separators and logic gates)are further explored by designing phononic crystal plates.This research has broad application prospects in information processing and vibration control,and potential applications in other directions are also worth exploring.
基金supported by the National Natural Science Foundation of China(No.12134002)。
文摘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.
基金sponsored by the National Natural Science Foundation of China Research(Grant No.41274138)the Science Foundation of China University of Petroleum(Beijing)(No.KYJJ2012-05-02)
文摘The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.
基金Project supported by the National Natural Science Foundation of China(Nos.11922209,11991031 and 12021002)。
文摘In this work,the three-dimensional(3 D)propagation behaviors in the nonlinear phononic crystal and elastic wave metamaterial with initial stresses are investigated.The analytical solutions of the fundamental wave and second harmonic with the quasilongitudinal(qP)and quasi-shear(qS_(1) and qS_(2))modes are derived.Based on the transfer and stiffness matrices,band gaps with initial stresses are obtained by the Bloch theorem.The transmission coefficients are calculated to support the band gap property,and the tunability of the nonreciprocal transmission by the initial stress is discussed.This work is expected to provide a way to tune the nonreciprocal transmission with vector characteristics.
基金supported by the National Science and Technology Major Project of China(Grant No. 2011ZX05004-003,2011ZX05014-006-006)the National Key Basic Research Program of China(Grant No. 2013CB228602)the Natural Science Foundation of China(Grant No. 40974066)
文摘The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.
基金Funding for this work has been provided by the National Natural Science Foundation of China(Nos.11872313 and 11502202)National Key R&D Program of China(2017YFB1102801)Fundamental Research Funds for the Central Universities and Seed Foundation of Innovation and Creation for Graduate Students in Northwestern Polytechnical University(CX2020107).
文摘In this study,the wave propagation properties of lattice metamaterials with Koch fractal structures are investigated in terms of band structures and directional wave propagation.The analytical models of lattice metamaterials are established using the finite element method,and the dispersion relation is solved using the Bloch’s theorem.The band structures of the lattice metamaterials with different numbers of iterations are studied,and the group velocities at a selected frequency are calculated to analyze the directional wave propagation characteristics.Furthermore,dynamic responses of the finite structures are calculated using commercial finite element software to verify the band gaps and directional wave propagation behaviors in the lattice metamaterials.The results show that multiple and low band gaps are present in the lattice materials with various geometric parameters of the Koch fractal,and the position of the lowest band gap decreases as the number of iterations increases.The results indicate the potential applications of lattice metamaterials with Koch fractals for vibration isolation and multi-functional design.
文摘This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.
文摘Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.
基金the National Basic Research Program (973) of China(No.2011CB013505)
文摘It is an important subject to probe the structure in the medium by various kinds of detection methods in the geotechnical engineering. Based on the propagation theory of elastic wave in half-space layered medium, the propagation characteristics of elastic wave in layered medium with different elastic parameters are discussed using dynamic analysis of finite element method. It is known that the S-wave velocity, density and thickness of layer are related to the properties of the elastic wave including waveform characteristics, spectral characteristics and time-frequency characteristics. We pay special attention to the structure with low velocity interlayer. The impact imaging method is applied to the grouting construction of the immersed tube tunnel. Data acquisition and analytical method are introduced in detail. The grouting effects can be qualitatively evaluated by comparing the characteristics of elastic wave before grouting with those after grouting. Finally, a quantitative evaluation is obtained according to the relationship between energy response of elastic wave and impedance ratio.
基金the supports by the National Natural Science Foundation of China(Grants 11922209,11991031 and 12021002)for this research work.
文摘Due to their potential properties unlike traditional materials and structures,elastic wave metamaterials have received significant interests in recent years.With the coupling between the acoustic and vibration,their mechanical characteristics can be tuned by the active feedback control system at low frequency ranges in which the traditional passive control is limited.This work illustrates that the superior performances of the effective mass density and sound pressure level(SPL)of an elastic wave metamaterial can be significantly changed by the active control,in which the periodic array of local resonators and orthogonal stiffeners are included.Significantly,based on the locally resonant mechanism,the negative density occurs over a frequency range.Due to the effects of lattice constant,structural damping and other parameters,the SPL with the function of fluid-solid coupling are illustrated and discussed.
基金supported by the National Science and Technology Major Project(2016ZX05006-002)the National Natural Science Foundation of China(Nos.41874153,41504097)
文摘Mesh-free finite difference(FD)methods can improve the geometric flexibility of modeling without the need for lattice mapping or complex meshing process.Radial-basisfunction-generated FD is among the most commonly used mesh-free FD methods and can accurately simulate seismic wave propagation in the non-rectangular computational domain.In this paper,we propose a perfectly matched layer(PML)boundary condition for a meshfree FD solution of the elastic wave equation,which can be applied to the boundaries of the non-rectangular velocity model.The performance of the PML is,however,severely reduced for near-grazing incident waves and low-frequency waves.We thus also propose the complexfrequency-shifted PML(CFS-PML)boundary condition for a mesh-free FD solution of the elastic wave equation.For two PML boundary conditions,we derive unsplit time-domain expressions by constructing auxiliary differential equations,both of which require less memory and are easy for programming.Numerical experiments demonstrate that these two PML boundary conditions effectively eliminate artificial boundary reflections in mesh-free FD simulations.When compared with the PML boundary condition,the CFS-PML boundary condition results in better absorption for near-grazing incident waves and evanescent waves.
基金Project supported by the National Natural Science Foundation of China(Nos.11172050 and11672047)the Science and Technology Foundation of China Academy of Engineering Physics(No.2013A0202011)
文摘A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding criterion. Based on the HLLCE, a third-order cell-centered Lagrangian scheme is built for one-dimensional elastic-plastic problems. A number of numerical experiments are carried out. The numerical results show that the proposed third-order scheme achieves the desired order of accuracy. The third-order scheme is used to the numerical solution of the problems with elastic shock waves and elastic rarefaction waves. The numerical results are compared with a reference solution and the results obtained by other authors. The comparison shows that the pre- sented high-order scheme is convergent, stable, and essentially non-oscillatory. Moreover, the HLLCE is more efficient than the two-rarefaction Riemann solver with elastic waves （TRRSE）
基金The authors wish to express gratitude for the support provided by the National Natural Science Foundation of China(Grant Nos.11922209,11991031 and 12021002).
文摘Using the active feedback control system on the elastic wave metamaterial,this research concentrates on the sound transmission with the dynamic effective model.The metamaterial is subjected to an incident pressure and immersed in the external mean flow.The elastic wave metamaterial consists of double plates and the upper and lower four-link mechanisms are attached inside.The vertical resonator is attached by the active feedback control system and connected with two four-link mechanisms.Based on the dynamic equivalent method,the metamaterial is equivalent as a single-layer plate by the dynamic effective parameter.With the coupling between the fluid and structure,the expression of the sound transmission loss(STL)is derived.This research shows the influence of effective mass density on sound transmission properties,and the STL in both modes can be tuned by the acceleration and displacement feedback constants.In addition,the dynamic response and the STL are also changed obviously by different values of structural damping,incident angle(i.e.,the elevation and azimuth angles)and Mach number of the external fluid with the mean flow property.The results for sound transmission by two methods are compared,i.e.,the virtual work principle for double plates and the dynamic equivalent method corresponding to a single one.This paper is expected to be helpful for understanding the sound transmission properties of both pure single-and double-plate models.
基金supported by the National Natural Science Foundation of China(No.41474110)
文摘Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.
基金The research was partially supported by the National Natural Science Foundation of China(Grant Nos.41902297,41872210)the Natural Science Foundation of Hubei Province(Grant No.2018CFB292)Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences(Grant No.Z017006).
文摘To understand the evolution of stress-induced elastic wave anisotropy,three triaxial experiments were performed on sandstone specimens with bedding orientations parallel,perpendicular,and oblique to the maximum principal stress.P-wave velocities along 64 different directions on each specimen were monitored frequently to understand the anisotropy change at various stress levels by fitting Thomsen’s anisotropy equation.The results show that the elastic wave anisotropy is very sensitive to mechanical loading.Under hydrostatic loading,the magnitude of anisotropy is reduced in all three specimens.However,under deviatoric stress loading,the evolution of anisotropic characteristics(magnitude and orientation of the symmetry axis)is bedding orientation dependent.Anisotropy reversal occurs in specimens with bedding normal/oblique to the maximum principal stress.P-wave anisotropyε0 is linearly related to volumetric strain Sv and dilatancy,indicating that stress-induced redistribution of microcracks has a significant effect on P-wave velocity anisotropy.The closure of initial cracks and pores aligned in the bedding direction contributes to the decrease of the anisotropy.However,opening of new cracks,aligned in the maximum principal direction,accounts for the increase of the anisotropy.The experimental results provide some insights into the microstructural behavior under loading and provide an experimental basis for seismic data interpretation and parameter selection in engineering applications.
文摘It is essential to assess the evolution of soil fabric as it has an important role in the mechanical responses of soils during complex loading conditions.This contribution carries out the physical experiments using three granular materials in the laboratory.The variations of compression and shear wave velocities(Vp and Vs)are investigated during load-unload cycles under dry and drained conditions.Supplementary discrete element method(DEM)simulations are performed to understand the evolution of soil fabric during the equivalent load-unload cycles using spherical particles.Vp and Vs are not always reversible even though the stress state regains its isotropic condition after unload,indicating that Vp and Vs are governed by not only the stress state but also the fabric change.The variations of Vp/Vs are density-and stress-dependent;a higher level of stress ratio(s01/s03)threshold is observed for denser packings to trigger a significant change in wave velocity ratio(Vp/Vs)for experimental results using spherical glass beads and simulation data using spherical particles.Considering the particle shape,a higher s01/s03 threshold is found for more angular particles than rounded particles.The DEM result reveals that Vp/Vs of spherical particles can be correlated linearly with the evolution of fabric ratio(Fver/Fhor)during loadunload in a pre-peak range under dry and drained conditions.
基金supported by the National Natural Science Foundation of China(11102122)
文摘Using a stiffness matrix method, we in- vestigate the propagation behaviors of elastic waves in one-dimensional （1D） piezoelectric/piezomagnetic （PE/PM） phononic crystals （PCs） with line defects by calculating energy reflection/transmittion coefficients of quasi-pressure and quasi-shear waves. Line defects are created by the re- placement of PE or PM constituent layer. The defect modes existing in the first gap are considered and the influences on defect modes of the material properties and volume fraction of the defect layers, the type of incident waves, the location of defect layer and the number of structural layers are discussed in detail. Numerical results indicate that defect modes are the most obvious when the defect layers are inserted in the middle of the perfect PCs; the types of incidence wave and material properties of the defect layers have important effects on the numbers, the location of frequencies and the peaks of defect modes, and the defect modes are strongly de- pendent on volume fraction of the defect layers. We hope this paper will be found useful for the design of PE/PM acoustic filters or acoustic transducer with PCs structures.
文摘Basedonthe first-order Biot-equation with simplified coefficients,astaggered irregu- lar-grid finite difference method(FDM)is developed to simulate elastic wave propagation in 3-D heterogeneous anisotropic porous media.The ‘slow’P wave in porous media wave simulation is highly dispersive.Finer grids are needed to get a precise wavefield calculation for models with curved interface and complex geometric structure.Fine grids in a global model greatly increase computation costs of regular grids scheme.Irregular fine or coarse grids in local fields not only cost less computing time than the conventional velocity-stress FDM,but also give a more accu- rate wavefield description.A dispersion analysis of the irregular-grid finite difference operator has confirmed the stability and high efficiency.The absorbing boundary condition is used to elimi- nate artificial reflections.Numerical examples show that this new irregular-grid finite difference method is of higher performance than conventional methods using regular rectangular grids in simulating elastic wave propagation in heterogeneous anisotropic porous media.
