Continuum model of two-dimensional crystal lattice of metamaterials

The continuum models of two-dimensional crystal lattice of metamaterial are investigated in this paper. The equivalent classical continuum requires the introduction of frequency-dependent mass density that becomes negative or infinite near the resonance frequency. In order to avoid the frequency-dependent mass density and make the dispersive characteristic of the elastic waves propagating in the equivalent continuum approximating that of lattice wave in two-dimensional crystal lattice of metamaterial, three kinds of continuum models, namely, the multiple displacements continuum model, the strain gradient continuum model and the nonlocal strain gradient continuum model, are investigated. It is found that the nonlocal gradient continuum model may better represent the dispersive relation and the bandgap characteristics of the metamaterial by the appropriate selection of nonlocal parameters.     

Microstructure continuum modeling of an elastic metamaterial

Elastic metamaterials have unusual microstructures that can make them exhibit unusual dynamic behavior. For instance, if treated as classical elastic solids, these materials may have frequency-dependent effective mass densities which may become negative in certain frequency range. In this study, an approach for developing microstructure continuum models to represent elastic metamaterials was presented. Subsequently, this continuum model was used to study wave propagation and band gaps in elastic metamaterial with resonators. In contrast to the use of the conventional continuum theory with which the effective mass density would become frequency-dependent and negative, the main advantage of the microstructure continuum model is that the local microstructural deformation/motion is accounted for with the introduction of additional kinematic variables. Moreover, the material constants of the microstructure continuum model are explicitly expressed in term of the properties of the host medium and the resonator. The accuracy of the microstructure continuum model was evaluated by comparing dispersion curves of harmonic waves to those obtained by the finite element analysis based on the exact geometry of the elastic metamaterial.     

Ultrasonic wave characteristics of a monolayer graphene on silicon substrate

The present work deals with an ultrasonic type of wave propagation characteristics of monolayer graphene on silicon (Si) substrate. An atomistic model of a hybrid lattice involving a hexagonal lattice of graphene and surface atoms of diamond lattice of Si is developed to identify the carbon-silicon bond stiffness. Properties of this hybrid lattice model is then mapped into a nonlocal continuum framework. Equivalent force constant due to Si substrate is obtained by minimizing the total potential energy of the system. For this equilibrium configuration, the nonlocal governing equations are derived to analyze the ultrasonic wave dispersion based on spectral analysis. From the present analysis we show that the silicon substrate affects only the flexural wave mode. The frequency band gap of flexural mode is also significantly affected by this substrate. The results also show that, the silicon substrate adds cushioning effect to the graphene and it makes the graphene more stable. The analysis also show that the frequency bang gap relations of in-plane (longitudinal and lateral) and out-of-plane (flexural) wave modes depends not only on the y-direction wavenumber but also on nonlocal scaling parameter. In the nonlocal analysis, at higher values of the y-directional wavenumber, a decrease in the frequency band gap is observed for all the three fundamental wave modes in the graphene–silicon system. The atoms movement in the graphene due to the wave propagation are also captured for all the tree fundamental wave modes. The results presented in this work are qualitatively different from those obtained based on the local analysis and thus, are important for the development of graphene based nanodevices such as strain sensor, mass and pressure sensors, atomic dust detectors and enhancer of surface image resolution that make use of the ultrasonic wave dispersion properties of graphene.     

Equivalent continuum modeling for non‐linear wave propagation in jointed media

This paper describes a methodology to build equivalent continuum (ECC) models to represent jointed media for non‐linear wave propagation. In the present work, we study the response of randomly jointed rock media both to quasi‐static and dynamic loadings. We have found that, unlike in the quasi‐static case where no relaxation is observed, during dynamic loading, the deviatoric stress drops in the plastic wave after reaching its peak. Such stress relaxation phenomenon is calculated in jointed media, although the model used for intact rock blocks is rate‐independent. This is a strong indication that EC models for jointed rock should include some rate‐dependence when used in wave propagation simulations. Copyright © 2010 John Wiley & Sons, Ltd.     

A micropolar continuum model for vibrating grid frameworks

A principle of converting a discrete system of grid frameworks to an equivalent micropolar continuum model is given with the degree of approximation taken into consideration. A micropolar continuum is then defined in the form of higher order extension. In order to supplement defects of previous theories, a complex-valued micropolar continuum model is constructed for grid frameworks vibrating with an arbitrary frequency by means of the variational principle related to the average energy of the system. An analysis of wave propagation reveals the existence of high frequency waves. The accuracy of solutions is also investigated.     

On the modelling of dynamic behavior of periodic lattice structures

Summary. The aim of this contribution is to propose and apply a new approach to the formulation of mathematical models for the analysis of dynamic behavior of dense periodic lattice structures (space or plane trusses) of an arbitrary form. The modelling approach is carried out on two levels. First, we formulate a discrete model, represented by the system of finite difference equations with respect to the spatial coordinates. The obtained equations describe both low- and high-frequency wave propagation problems. Second, two continuum models are derived directly from the finite difference equations and represented respectively by the second- and the fourth-order PDEs with constant coefficients. These models have a physical sense provided that the considerations are restricted to the long wave propagation phenomena. The proposed approach is applied to the vibration analysis for a certain plane lattice structure. Special attention is given to the determination of the range of applicability of the continuum models.     

Remarks on tension instability of Eulerian and Lagrangian corrected smooth particle hydrodynamics (CSPH) methods

The paper discusses the problem of tension instability of particle‐based methods such as smooth particle hydrodynamics (SPH) or corrected SPH (CSPH). It is shown that tension instability is a property of a continuum where the stress tensor is isotropic and the value of the pressure is a function of the density or volume ratio. The paper will show that, for this material model, the non‐linear continuum equations fail to satisfy the stability condition in the presence of tension. Consequently, any discretization of this continuum will result in negative eigenvalues in the tangent stiffness matrix that will lead to instabilities in the time integration process. An important exception is the 1‐D case where the continuum becomes stable but SPH or CSPH can still exhibit negative eigenvalues. The paper will show that these negative eigenvalues can be eliminated if a Lagrangian formulation is used whereby all derivatives are referred to a fixed reference configuration. The resulting formulation maintains the momentum preservation properties of its Eulerian equivalent. Finally a simple 1‐D wave propagation example will be used to demonstrate that a stable solution can be obtained using Lagrangian CSPH without the need for any artificial viscosity. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical simulation of acoustic band gaps in homogenized elastic composites

The dispersion of acoustic or elastodynamic waves in elastic composites are studied using the homogenized model. We consider heterogeneous periodic structures consisting of soft but heavy inclusions embedded in a stiffer matrix. By virtue of the asymptotic homogenization technique in conjunction with an appropriate scaling of the elasticity coefficients in the inclusions, the limit model exhibits the band gaps in wave propagation due to the negative effective mass. This phenomenon can be revealed by studying guided waves in discrete mass-spring structures with scale-dependent parameters. The main purpose of the paper is to justify the applicability of the homogenized model of the heterogeneous elastic continuum for prediction of the band gaps in structures featured by a finite scale of heterogeneities. We show the band gaps numerical identification and discus aspects of anisotropy, microstructure geometry and material contrast between the constituents in the context of the long wave dispersion.     

Magnetoactive Acoustic Metamaterials

Acoustic metamaterials with negative constitutive parameters (modulus and/or mass density) have shown great potential in diverse applications ranging from sonic cloaking, abnormal refraction and superlensing, to noise canceling. In conventional acoustic metamaterials, the negative constitutive parameters are engineered via tailored structures with fixed geometries; therefore, the relationships between constitutive parameters and acoustic frequencies are typically fixed to form a 2D phase space once the structures are fabricated. Here, by means of a model system of magnetoactive lattice structures, stimuli‐responsive acoustic metamaterials are demonstrated to be able to extend the 2D phase space to 3D through rapidly and repeatedly switching signs of constitutive parameters with remote magnetic fields. It is shown for the first time that effective modulus can be reversibly switched between positive and negative within controlled frequency regimes through lattice buckling modulated by theoretically predicted magnetic fields. The magnetically triggered negative‐modulus and cavity‐induced negative density are integrated to achieve flexible switching between single‐negative and double‐negative. This strategy opens promising avenues for remote, rapid, and reversible modulation of acoustic transportation, refraction, imaging, and focusing in subwavelength regimes.     

Continuum model of acoustic metamaterials with diatomic crystal lattice

The continuum model of one-dimensional acoustic metamaterial with diatomic crystal lattice is studied in this article. First, the dispersive relation of lattice wave in one-dimensional diatomic crystal lattice of metamaterial is established and compared with that of the classic material. Then, the continuum model of the acoustic metamaterial leads to the classical continuum model, the stain gradient continuum model, and the nonlocal gradient continuum model based on different assumptions. The dispersive curves, which correspond to the three kinds of models, are shown graphically and compared with that of discrete crystal lattice of metamaterial. The disadvantage of the classic continuum model and the strain gradient continuum model are discussed. The nonlocal gradient continuum model is derived based on the nonlocal assumption of continuous displacement field. The stability of dispersive curves is guaranteed. However, the actual prediction effects are still dependent upon the appropriate selection of the nonlocal parameter in the nonlocal gradient continuum model.     

Study of wave propagation in nanowires with surface effects by using a high-order continuum theory

A high-order continuum model is developed to study wave propagation in nanowires. By using the model, heterogeneous nanostructure effects can be captured especially for high wave frequency cases. Surface stress effects are also included by using the incremental deformation approach. The governing equations of motion in the nanowire are derived including both the strain-independent and strain-dependent surface stresses. For simplicity and clarity, specific attention will be paid to the effects of strain-independent surface stress in this study. The accuracy of the proposed model is validated by comparing dispersion curves of longitudinal wave propagation from the current model with those from the exact solution. By conducting a reduced formulation, the results predicted by the current model will be compared with those based on existed high-order models to show capability of the current model. Numerical simulations are then conducted to study both longitudinal and flexural wave propagation in nanowires. The surface stress effects upon both longitudinal and flexural wave propagation in nanowires are demonstrated, from which the size dependent wave information in nanowires can be observed. Some new physical wave phenomena related to the surface stress effects are discussed.     

Anomalies in stiffness and damping of a 2D discrete viscoelastic system due to negative stiffness components

The recent development of using negative stiffness inclusions to achieve extreme overall stiffness and mechanical damping of composite materials reveals a new avenue for constructing high performance materials. One of the negative stiffness sources can be obtained from phase transforming materials in the vicinity of their phase transition, as suggested by the Landau theory. To understand the underlying mechanism from a microscopic viewpoint, we theoretically analyze a 2D, nested triangular lattice cell with pre-chosen elements containing negative stiffness to demonstrate anomalies in overall stiffness and damping. Combining with current knowledge from continuum models, based on the composite theory, such as the Voigt, Reuss, and Hashin-Shtrikman model, we further explore the stability of the system with Lyapunov's indirect stability theorem. The evolution of the microstructure in terms of the discrete system is discussed. A potential application of the results presented here is to develop special thin films with unusual in-plane mechanical properties.     

Micro-mechanical modelling of granular material. Part 2: Plane wave propagation in infinite media

Summary This paper discusses the propagation of plane body waves through a second-gradient micropolar elastic continuum. In an accompanying paper, this macroscopic constitutive law has been derived from the micro-level particle characteristics, which are the inter-particle stiffness, the particle size and the package density. As a result of incorporating the micro-scale effects, the body waves propagate in a dispersive manner, where dispersion becomes more prominent when the wavelength of the generated body waves reaches the order of magnitude of the particle size. After successively deriving the equations of motion and the dispersion relations for plane body wave propagation, the compressional wave properties for the second-gradient micro-polar model are compared to those for the Born-Karman lattice structure. Furthermore, distinguished features of the second-gradient micro-polar model are exhibited by comparing the dispersion relations of the coupled propagation of the shear wave and the micro-rotational wave with those of more simple constitutive models. The paper ends with a parameter study, where the effect by the translational particle contact stiffness and the rotational particle contact stiffness is examined.     

Wave propagation in presence of oscillators on the free surface

This article is devoted to wave propagation in presence of a periodic distribution of oscillators (3D sprung-damped mass) on the surface of an elastic medium. The case of an elastic half space loaded by quasi-periodic surface forces is studied at first. The equivalent boundary conditions at the macro-scale are derived at the first and second order, by assuming the existence of a boundary layer and using a two-dimensional homogenisation method. When the force distribution results from attached oscillators, it is shown that boundary conditions can be expressed in terms of equivalent surface impedance at the first order, with local and non-local correctors at the second order. These results are used to study wave refraction by an oscillator layer. The main phenomena—atypical redistribution of mode and mode conversion, frequency range of efficiency, characteristic time of response—are identified. The first correctors are determined analytically for plane waves of oblique incidence. Finally the validity range of the modelling is discussed.     

Enhanced continua and discrete lattices for modelling granular assemblies

This article discusses the derivation of continuum models that can be used for modelling the inhomogeneous mechanical behaviour of granular assemblies. These so-called kinematically enhanced models are of the strain-gradient type and of the strain-gradient micro-polar type, and are derived by means of homogenizing the micro-structural interactions between discrete particles. By analysis of the body wave dispersion curves, the enhanced continuum models are compared to corresponding discrete lattice models. Accordingly, it can be examined up to which deformation level the continuum models are able to accurately describe the discrete particle behaviour. Further, the boundary conditions for the enhanced continuum models are formulated, and their stability is considered. It is demonstrated how to use the body wave dispersion relations for the assessment of stability.     

Micromechanics of seismic wave propagation in granular materials

In this study experimental data on a model soil in a cubical cell are compared with both discrete element (DEM) simulations and continuum analyses. The experiments and simulations used point source transmitters and receivers to evaluate the shear and compression wave velocities of the samples, from which some of the elastic moduli can be deduced. Complex responses to perturbations generated by the bender/extender piezoceramic elements in the experiments were compared to those found by the controlled movement of the particles in the DEM simulations. The generally satisfactory agreement between experimental observations and DEM simulations can be seen as a validation and support the use of DEM to investigate the influence of grain interaction on wave propagation. Frequency domain analyses that considered filtering of the higher frequency components of the inserted signal, the ratio of the input and received signals in the frequency domain and sample resonance provided useful insight into the system response. Frequency domain analysis and analytical continuum solutions for cube vibration show that the testing configuration excited some, but not all, of the system’s resonant frequencies. The particle scale data available from DEM enabled analysis of the energy dissipation during propagation of the wave. Frequency domain analysis at the particle scale revealed that the higher frequency content reduces with increasing distance from the point of excitation.     

Generation of chaotic dissipative solitons in active ring resonator with one-dimensional periodic ferromagnetic miscrostructure

Generation of chaotic dissipative solitons has been observed in an active ring resonator with one-dimensional periodic ferromagnetic microstructure comprising a single-crystalline film of yttrium iron garnet (YIG) with a lattice of grooves oriented perpendicular to the direction of propagation of a magnetostatic surface wave (MSSW). A quasi-periodic train of chaotic dissipative solitons was generated in this system under the conditions of three-magnon processes of MSSW decay due to the passive synchronization (PS) of spin wave self-modulation in a frequency band corresponding to the first bandgap. The PS onset was caused by the saturable absorption of microwave signals within the bandgap of the MSSW transmission line.     

A partition of unity‐based multiscale approach for modelling fracture in piezoelectric ceramics

The development of models for a priori assessment of the reliability of micro electromechanical systems is of crucial importance for the further development of such devices. In this contribution a partition of unity‐based cohesive zone finite element model is employed to mimic crack nucleation and propagation in a piezoelectric continuum. A multiscale framework to appropriately represent the influence of the microstructure on the response of a miniaturized component is proposed. It is illustrated that using the proposed multiscale method a representative volume element exists. Numerical simulations are performed to demonstrate the constitutive homogenization framework. Copyright © 2009 John Wiley & Sons, Ltd.     

Timoshenko-beam effects on transverse wave propagation in carbon nanotubes

This paper studies effects of rotary inertia and shear deformation on transverse wave propagation in individual carbon nanotubes (CNTs) within terahertz range. Detailed results are demonstrated for transverse wave speeds of doublewall CNTs, based on Timoshenko-beam model and Euler-beam model, respectively. The present models predict some terahertz critical frequencies at which the number of wave speeds changes. The effects of rotary inertia and shear deformation are negligible and transverse wave propagation can be described satisfactorily by the existing single-Euler-beam model only when the frequency is far below the lowest critical frequency. When the frequency is below but close to the lowest critical frequency, rotary inertia and shear deformation come to significantly affect the wave speed. Furthermore, when the frequency is higher than the lowest critical frequency, more than one wave speed exists and transverse waves of given frequency could propagate at various speeds that are considerably different than the speed predicted by the single-Euler-beam model. In particular, rotary inertia and shear deformation have a significant effect on both the wave speeds and the critical frequencies especially for CNTs of larger radii. Hence, terahertz transverse wave propagation in CNTs should be better modeled by Timoshenko-beam model, instead of Euler-beam model.