Wave Propagation In Elastic Solids Research Articles

Planar piezoelectric metamaterials: Sound transmission and applicable frequency range in oblique incidence

Frequency limit and applicable frequency range of planar piezoelectric metamaterials connected to external circuits have not been well defined in estimating the sound transmission loss. This article extends the classical transfer matrix method for use in evaluating the sound transmission of thin-plate piezoelectric metamaterials in oblique incidence. Using the Kirchhoff thin plate theory, a modified transfer matrix method that takes factors of external circuits into consideration is developed for attenuation and control of acoustic waves. Several vibro-acoustic analytical models are compared, including the Kirchhoff thin plate theory, the Reissner–Mindlin thick plate theory and the theory of wave propagation in elastic solids. These theories are used to determine the dispersion relation, coincidence and transition frequency of thin and thick plate theories in analysing piezoelectric acoustic metamaterials, alongside a validation using the finite element method. Then acoustic properties of piezoelectric plates connected to passive external circuits are studied parametrically with both dimensional and dimensionless variables based on an equivalent Kirchhoff plate approximation. The findings show that external electrical impedance alone can be used to adjust the resonance frequency over a broad range and thereby control sound transmission loss. This provides considerable flexibility in modifying the acoustic properties of the piezoelectric metamaterial in comparison to traditional, fixed-structure metamaterials. The study indicates a straightforward and powerful analytical approach for the optimization of acoustic insulation using thin-plate piezoelectric metamaterials.

Homogenization and analysis of planar piezoelectric metamaterials: A comprehensive study

This article presents a comprehensive investigation into the homogenization and analysis of planar piezoelectric metamaterials. The classical transfer matrix method is extended by using homogenization methods to evaluate the sound transmission of piezoelectric metamaterials in both normal and oblique incidences. This has been validated numerically using the finite element method and experimentally in the solid medium propagation in normal incidence. Several vibro-acoustic analytical models are compared in oblique incidence, including the Kirchhoff thin plate theory, the Reissner-Mindlin thick plate theory, and the theory of wave propagation in elastic solids. These theories are used to determine the dispersion relation, coincidence, and transition frequency of thin and thick plate theories in analyzing piezoelectric metamaterials. Additionally, acoustic properties of piezoelectric layers connected to external circuits are parametrically studied, exploring both dimensional and dimensionless variables. Results indicate that significant control over the resonance frequency and sound transmission can be enabled by adjusting the external electrical impedance, flat-layer structures, and piezoelectric materials. This demonstrates excellent tunability and compactness of planar piezoelectric metamaterials for space-sensitive applications. The study indicates a straightforward and powerful analytical approach for the optimization of acoustic insulation using planar piezoelectric metamaterials.

Wave Polarization Control in Anisotropic Locally Resonant Materials

Elastic wave propagation in solids can be controlled and manipulated by properly designed metamaterials. In particular, polarization conversion can be obtained by using anisotropic materials. In this paper, we propose a three-component locally resonant material with non-symmetrically coated inclusions, and we study the effect of the anisotropic equivalent mass on band gap formation and the polarization conversion of elastic waves. The equivalent frequency-dependent mass tensor is obtained through the two-scale homogenization approach. The study of the eigenvalues of the mass tensor enables to predict band gaps and polarization bands, as well as identifying a priori the effect of different geometric and material parameters, thus opening the way to metamaterial optimization.

Dirichlet and Neumann boundary conditions in a lattice Boltzmann method for elastodynamics

AbstractRecently, Murthy et al. (Commun Comput Phys 2:23, 2017. http://dx.doi.org/10.4208/cicp.OA-2016-0259 ) and Escande et al. (Lattice Boltzmann method for wave propagation in elastic solids with a regular lattice: theoretical analysis and validation, 2020. arXiv.doi:1048550/ARXIV.2009.06404. arXiv:2009.06404) adopted the Lattice Boltzmann Method (LBM) to model the linear elastodynamic behaviour of isotropic solids. The LBM is attractive as an elastodynamic solver because it can be parallelised readily and lends itself to finely discretised simulations of dynamic effects in continua, allowing transient phenomena such as wave propagation to be modeled efficiently. This work proposes simple local boundary rules which approximate the behaviour of Dirichlet and Neumann boundary conditions with an LBM for elastic solids. The boundary rules are shown to be consistent with the target boundary values in the first order. An empirical convergence study is performed for the transient tension loading of a rectangular plate, with a Finite Element (FE) simulation being used as a reference. Additionally, we compare results produced by the LBM for the sudden loading of a stationary crack with an analytical solution from Freund (Dynamic fracture mechanics. Cambridge Monographs on Mechanics. Cambridge University Press, Cambridge, 1990. https://doi.org/10.1017/CBO9780511546761).

Vibrational loss analysis of a new type of phononic crystal with a tungsten block embedded inside a rubber matrix

Vibration and noise pollution is one of the main pollution in modern society. In order to obtain good vibration damping effect, more effective phononic crystals should be designed. Based on the theory of elastic wave propagation in solids, three types of phononic crystals are designed in this paper using tungsten blocks embedded in a rubber matrix, which are cup-shaped phononic crystal, solid cylindrical phononic crystal, and hollow cylindrical phononic crystal. Firstly, the band gap characteristics and vibration losses of the three phononic crystals are analyzed by using the finite element method. Secondly, the physical mechanism of band gap formation is explored by vibration modes. Finally, the cup-shaped phononic crystal was introduced into the core layer of the sandwich plate to form the cup-shaped phononic crystal sandwich plate, and its vibration damping performance was analyzed. The results show that the three phononic crystals can form three band gaps in the range of 0–800 Hz, and the first low-frequency band gap starts at about 140 Hz and is all wider than 200 Hz. It is noteworthy that the average loss of vibration transmission of the three types of phononic crystals is more than 69 dB, which possesses stronger damping capability and wider low- and middle-frequency band gaps than the flat plate type phononic crystals. The vibration direction of the phononic crystal is at an angle of 90° to the wave vector, which prevents the propagation of elastic waves. The phononic crystal embedded in the sandwich panel plate is more in line with the actual demand of vibration damping and obtains good vibration damping effect. The research in this paper can provide more feasibility for phononic crystal damping.

Modeling wave propagation in elastic solids via high-order accurate implicit-mesh discontinuous Galerkin methods

A high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial discretization, which enables boundary and interface conditions to be enforced with high-order accuracy on curved embedded geometries. High-order accuracy is achieved via high-order quadrature rules for implicitly-defined domains and boundaries, whilst a cell-merging strategy addresses the presence of small cut cells. The framework is used to discretize the governing equations of elastodynamics, written using a first-order hyperbolic momentum-strain formulation, and an exact Riemann solver is employed to compute the numerical flux at the interface between dissimilar materials with general anisotropic properties. The space-discretized equations are then advanced in time using explicit high-order Runge–Kutta algorithms. Several two- and three-dimensional numerical tests including dynamic adaptive mesh refinement are presented to demonstrate the high-order accuracy and the capability of the method in the elastodynamic analysis of single- and bi-phases solids containing complex geometries.

A Lagrangian DG-Method for Wave Propagation in a Cracked Solid with Frictional Contact Interfaces

We developed a discontinuous Galerkin (DG) numerical scheme for wave propagation in elastic solids with frictional contact interfaces. This type of numerical scheme is useful in investigations of wave propagation in elastic solids with micro-cracks (cracked solid) that involve modeling the damage in brittle materials or architected meta-materials. Only processes with mild loading that do not trigger crack fracture extension or the nucleation of new fractures are considered. The main focus lies on the contact conditions at crack surfaces, including provisions for crack opening and closure and stick-and-slip with Coulomb friction. The proposed numerical algorithm uses the leapfrog scheme for the time discretization and an augmented Lagrangian algorithm to solve the associated non-linear problems. For efficient parallelization, a Galerkin discontinuous method was chosen for the space discretization. The frictional interfaces (micro-cracks), where the numerical flux is obtained by solving non-linear and non-smooth variational problems, concerns only a limited number the degrees of freedom, implying a small additional computational cost compared to classical bulk DG schemes. The numerical method was tested through two model problems with analytical solutions. The proposed Lagrangian approach of the nonlinear interfaces had excellent results (stability and high accuracy) and only required a reasonable additional amount of computational effort. To illustrate the method, we conclude with some numerical simulations on the blast propagation in a cracked material and in a meta-material designed for shock dissipation.

Experimental observation of transverse and longitudinal wave propagation in a metamaterial periodically arrayed with nonlinear resonators

Elastic wave propagation in solids composed of elastic metamaterials always exists as both transverse and longitudinal waves simultaneously. By using the Hamilton principle, Bloch theorem, and finite-element method, an approach is proposed to explore the band energy structures of the metamaterial, this is observed in an experiment on an ultralow frequency metamaterial periodically positioned nonlinear resonators and this means the broadband isolation of both longitudinal and transverse waves was verified. Comparison of the propagation of transverse and longitudinal waves in the ultralow frequency metamaterials, the Bragg band-gaps (BBGs) and local resonance band-gaps (LRBGs) differ with both the frequency of the bandwidth and the depth of the band-gap.

An analysis of axisymmetric Sezawa waves in elastic solids

The wave propagation in elastic solids covered by a thin layer has received significant attention due to the existence of Sezawa waves in many applications such as medical imaging. With a Helmholtz decomposition in cylindrical coordinates and subsequent solutions with Bessel functions, it is found that the velocity of such Sezawa waves is the same as the one in Cartesian coordinates, but the displacement will be decaying along the radius with eventual conversion to plane waves. The decaying with radius exhibits a strong contrast to the uniform displacement in the Cartesian formulation, and the asymptotic approximation is accurate in the range about one wavelength away from the origin. The displacement components in the vicinity of origin are naturally given in Bessel functions which can be singular, making it more suitable to analyze waves excited by a point source with solutions from cylindrical coordinates. This is particularly important in extracting vital wave properties and reconstructing the waveform in the vicinity of source of excitation with measurement data from the outer region.

A non-local method in peridynamic theory for simulating elastic wave propagation in solids

In order to simulate the elastic wave propagation in unbounded solids, this paper introduces a non-local method to construct absorbing boundary conditions in the bond-based peridynamic theory. To construct non-local absorbing boundary conditions in the form of absorbing layers with increasing damping, we first re-derived the particle motion equilibrium equation of the bond-based peridynamic theory, based on the principle of virtual work and Lagrange's equation on the premise of considering material damping. We also give the expression for the inter-particle damping force density which satisfies the balance laws of linear momentum and angular momentum. This paper also provides a process and algorithm for realizing the non-local absorbing boundary conditions in the bond-based peridynamic theory, and adopts the Fortran language for computer program compilation. Finally, we simulate the elastic wave propagation in semi-unbounded and unbounded solids as numerical examples. The results show that such non-local absorbing boundary conditions can stably absorb and attenuate incident elastic waves. Since it is derived in the time domain, this method does not require frequency-domain operations, such as Fourier and Laplace transformations, nor does it require wave field splitting.

Normalized SPH without boundary deficiency and its application to transient solid mechanics problems

Smoothed particle hydrodynamics (SPH) method is a powerful tool for modeling solid mechanics problems, especially for large deformation problems. However, it suffers from boundary deficiency and difficulty of boundary condition treatment. In this work, a normalized SPH method is proposed to overcome these problems. The method is based on a newly developed normalized particle approximation. To derive this particle approximation, a normalized kernel approximation which is accurate for derivatives of linear functions everywhere in a problem domain is constructed, and all integral terms of the normalized kernel approximation including boundary terms are discretized by particle summations. The normalized particle approximation is free of matrix inversion, consequently attractive in computational stability and simplicity compared with other corrective particle approximations. Its approximation accuracy is demonstrated by calculating derivatives of test functions. Based on this particle approximation, the formulation of the normalized SPH method for transient solid mechanics problems is derived. Moreover, a direct method of treating traction boundary conditions is presented by making use of the boundary term of the normalized particle approximation. The accuracy and capability of the normalized SPH method are validated by the calculation of elastic wave propagation in solids and compared with commonly used SPH method.

Simulation of crack induced nonlinear elasticity using the combined finite-discrete element method

Numerical simulation of nonlinear elastic wave propagation in solids with cracks is indispensable for decoding the complicated mechanisms associated with the nonlinear ultrasonic techniques in Non-Destructive Testing (NDT). Here, we introduce a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), to explicitly simulate the crack induced nonlinear elasticity in solids with both horizontal and inclined cracks. In the FDEM model, the solid is discretized into finite elements to capture the wave propagation in the bulk material, and the finite elements along the two sides of the crack also behave as discrete elements to track the normal and tangential interactions between crack surfaces. The simulation results show that for cracked models, nonlinear elasticity is generated only when the excitation amplitude is large enough to trigger the contact between crack surfaces, and the nonlinear behavior is very sensitive to the crack surface contact. The simulations reveal the influence of normal and tangential contact on the nonlinear elasticity generation. Moreover, the results demonstrate the capabilities of FDEM for decoding the causality of nonlinear elasticity in cracked solid and its potential to assist in Non-Destructive Testing (NDT).

Lattice Boltzmann Method for Wave Propagation in Elastic Solids

Lattice Boltzmann Method for Wave Propagation in Elastic Solids

Spectral BEM for the Analysis of Wave Propagation and Fracture Mechanics

This paper presents a spectral boundary element formulation for analysis of structures subjected to dynamic loading. Two types of spectral elements based on Lobatto polynomials and Legendre polynomials are used. Two-dimensional analyses of elastic wave propagation in solids with and without cracks are carried out in the Laplace frequency domain with both conventional BEM and the spectral BEM. By imposing the requirement of the same level of accuracy, it was found that the use of spectral elements, compared with conventional quadratic elements, reduced the total number of nodes required for modeling high-frequency wave propagation. Benchmark examples included a simple one-dimensional bar for which analytical solution is available and a more complex crack problem where stress intensity factors were evaluated. Special crack tip elements are developed for the first time for the spectral elements to accurately model the crack tip fields. Although more integration points were used for the integrals associated with spectral elements than the conventional quadratic elements, shorter computation times were achieved through the application of the spectral BEM. This indicates that the spectral BEM is a more efficient method for the numerical modeling of structural health monitoring (SHM) processes, in which high-frequency waves are commonly used to detect damage, such as cracks, in structures.

Two dimensional modeling of elastic wave propagation in solids containing cracks with rough surfaces and friction – Part II: Numerical implementation

Our study aims at the creation of a numerical toolbox that describes wave propagation in samples containing internal contacts (e.g. cracks, delaminations, debondings, imperfect intergranular joints) of known geometry with postulated contact interaction laws including friction. The code consists of two entities: the contact model and the solid mechanics module. Part I of the paper concerns an in-depth description of a constitutive model for realistic contacts or cracks that takes into account the roughness of the contact faces and the associated effects of friction and hysteresis. In the crack model, three different contact states can be recognized: contact loss, total sliding and partial slip. Normal (clapping) interactions between the crack faces are implemented using a quadratic stress-displacement relation, whereas tangential (friction) interactions were introduced using the Coulomb friction law for the total sliding case, and the Method of Memory Diagrams (MMD) in case of partial slip. In the present part of the paper, we integrate the developed crack model into finite element software in order to simulate elastic wave propagation in a solid material containing internal contacts or cracks. We therefore implemented the comprehensive crack model in MATLAB® and introduced it in the Structural Mechanics Module of COMSOL Multiphysics®. The potential of the approach for ultrasound based inspection of solids with cracks showing acoustic nonlinearity is demonstrated by means of an example of shear wave propagation in an aluminum sample containing a single crack with rough surfaces and friction.

Two dimensional modeling of elastic wave propagation in solids containing cracks with rough surfaces and friction – Part I: Theoretical background

Our study aims at the creation of a numerical toolbox that describes wave propagation in samples containing internal contacts (e.g. cracks, delaminations, debondings, imperfect intergranular joints) of known geometry with postulated contact interaction laws including friction. The code consists of two entities: the contact model and the solid mechanics module. Part I of the paper concerns the modeling of internal contacts (called cracks for brevity), while part II is related to the integration of the developed contact model into a solid mechanics module that allows the description of wave propagation processes. The contact model is used to produce normal and tangential load-displacement relationships, which in turn are used by the solid mechanics module as boundary conditions at the internal contacts. Due to friction, the tangential reaction curve is hysteretic and memory-dependent. In addition, it depends on the normal reaction curve. An essential feature of the proposed contact model is that it takes into account the roughness of the contact faces. On one hand, accounting for roughness makes the contact model more complicated since it gives rise to a partial slip regime when some parts on the contact area experience slip and some do not. On the other hand, as we will show, the concept of contact surfaces covered by asperities receding under load makes it possible to formulate a consistent contact model that provides nonlinear load-displacement relationships for any value of the drive displacements and their histories. This is a strong advantage, since this way, the displacement-driven model allows for a simple explicit procedure of data exchange with the solid mechanics module, while more traditional flat-surface contacts driven by loads generate a complex iterative procedure. More specifically, the proposed contact model is based on the previously developed method of memory diagrams that allows one to automatically obtain memory-dependent solutions to frictional contact problems in the particular case of partial slip. Here we extend the solution onto cases of total sliding and contact loss which is possible while using the displacement-driven formulation. The method requires the knowledge of the normal contact response obtained in our case as a result of statistical consideration of roughness of contact faces.

High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: Viscous heat-conducting fluids and elastic solids

This paper is concerned with the numerical solution of the unified first order hyperbolic formulation of continuum mechanics recently proposed by Peshkov & Romenski, denoted as HPR model. In that framework, the viscous stresses are computed from the so-called distortion tensor A, which is one of the primary state variables. A very important key feature of the model is its ability to describe at the same time the behavior of inviscid and viscous compressible Newtonian and non-Newtonian fluids with heat conduction, as well as the behavior of elastic and visco-plastic solids. This is achieved via a stiff source term that accounts for strain relaxation in the evolution equations of A. Also heat conduction is included via a first order hyperbolic evolution equation of the thermal impulse, from which the heat flux is computed. The governing PDE system is hyperbolic and fully consistent with the principles of thermodynamics. It is also fundamentally different from first order Maxwell-Cattaneo-type relaxation models based on extended irreversible thermodynamics. The connection between the HPR model and the classical hyperbolic-parabolic Navier-Stokes-Fourier theory is established via a formal asymptotic analysis in the stiff relaxation limit. From a numerical point of view, the governing partial differential equations are very challenging, since they form a large nonlinear hyperbolic PDE system that includes stiff source terms and non-conservative products. We apply the successful family of one-step ADER-WENO finite volume and ADER discontinuous Galerkin finite element schemes in the stiff relaxation limit, and compare the numerical results with exact or numerical reference solutions obtained for the Euler and Navier-Stokes equations. To show the universality of the model, the paper is rounded-off with an application to wave propagation in elastic solids.

Modeling granular soils liquefaction using coupled lattice Boltzmann method and discrete element method

In this paper, a novel coupled pore-scale model of pore-fluid interacting with discrete particles is presented for modeling liquefaction of saturated granular soil. A microscale idealization of the solid phase is achieved using the discrete element method (DEM) while the fluid phase is modeled at a pore-scale using the lattice Boltzmann method (LBM). The fluid forces applied on the particles are calculated based on the momentum exchange between the fluid and particles. The presented model is based on a first principles formulation in which pore-pressure develops due to actual changes in pore space as particles׳ rearrangement occurs during shaking. The proposed approach is used to model the response of a saturated soil deposit subjected to low and large amplitude seismic excitations. Results of conducted simulations show that at low amplitude shaking, the input motion propagates following the theory of wave propagation in elastic solids. The deposit response to the strong input motion indicates that liquefaction took place and it was due to reduction in void space during shaking that led to buildup in pore-fluid pressure. Soil liquefaction was associated with soil stiffness degradation and significant loss of interparticle contacts. Simulation results also indicate that the level of shaking-induced shear strains and associated volumetric strains play a major role in the onset of liquefaction and the rate of pore-pressure buildup.

Complex wavenumber Fourier analysis of the B-spline based finite element method

We present the results of one-dimensional complex wavenumber Fourier analysis of the B-spline variant of Finite Element Method (FEM). Generally, numerical results of elastic wave propagation in solids obtained by FEM are polluted by dispersion and attenuation. It was shown for the higher-order B-spline based FEM, that the optical modes did not occur in the case of infinite domains, unlike the higher-order Lagrangian and Hermitian finite elements, and also the dispersion errors are smaller. The paper’s main focus is on the wave propagation through B-spline multi-patch/segment discretization with the C0 connection of B-spline segments and, chiefly, to the determining of dispersion and attenuation dependences. The numerical approach employed leads to substantial minimization of dispersion errors. Furthermore, the errors decrease in line with the increasing order of the B-spline elements/segments, with the local refinement, and also by the particular choice of the positions of control points through the optimizing procedure.

Coupled method of finite and dynamic infinite elements for simulating wave propagation in elastic solids involving infinite domains

This paper deals with the coupled method of finite and dynamic infinite elements for simulating wave propagation in elastic and viscoelastic solids involving infinite domains. This method can be used to simultaneously simulate material complexities in the near field and the infinite extent of the far field. Based on the governing equations of wave motion in two-dimensional and three-dimensional elastic/viscoelastic solids, the mass and stiffness matrices of the dynamic infinite element have been derived. The proposed two-dimensional dynamic infinite element can be used to simulate both the P-wave and the SV-wave propagation within the element, while the proposed three-dimensional dynamic infinite element can be used to simultaneously simulate the Rayleigh wave, P-wave and S-wave propagation within the element. The related simulation results have demonstrated that the coupled method of finite and dynamic infinite elements can be accurately used to simulate, both physically and computationally, wave propagation in elastic/viscoelastic solids involving infinite domains. Thus, this method provides an advanced scientific tool for dealing with both scientific and engineering problems involving infinite domains.