Universal power-law decay of the impulse energy in granular protectors

Hybrid Composite Material Reinforced with Carbon Nanolaminates for Gradient Stiffness: Preparation and Characterization

Currently, the procurement of lightweight, tough, and impact resistant materials is garnering significant industrial interest. New hybrid materials can be developed on the basis of the numerous naturally found materials with gradient properties found in nature. However, previous studies on granular materials demonstrate the possibility of capturing the energy generated by an impact within the material itself, thus deconstructing the initial impulse into a series of weaker impulses, dissipating the energy through various mechanisms, and gradually releasing undissipated energy. This work focuses on two production methods: spin coating for creating a granular material with composition and property gradients (an acrylonitrile–butadiene–styrene (ABS) polymer matrix reinforced by carbon nanolaminates at 0.10%, 0.25%, and 0.50%) and 3D printing for generating viscoelastic layers. The aim of this research was to obtain a hybrid material from which better behaviour against shocks and impacts and increased energy dissipation capacity could be expected when the granular material and viscoelastic layers were combined. Nondestructive tests were employed for the morphological characterization of the nanoreinforcement and testing reinforcement homogeneity within the matrix. Furthermore, the Voronoï tessellation method was used as a mathematical method to supplement the results. Finally, mechanical compression tests were performed to reveal additional mechanical properties of the material that had not been specified by the manufacturer of the 3D printing filaments.

Time scale disparity yielding acoustic nonreciprocity in a two-dimensional granular-elastic solid interface with asymmetry

We study nonreciprocal wave transmission across the interface of two dissimilar granular media separated by an elastic solid medium. Specifically, a left, larger-scale and a right smaller-scale granular media composed of two-dimensional, initially uncompressed hexagonally packed granules are interfacing with an intermediate linearly elastic solid, modeled either as a thin elastic plate or a linear Euler-Bernoulli beam. The granular media are modeled by discrete elements and the elastic solid by finite elements assuming a plane stress approximation for the thin plate. Accounting for the combined effects of Hertzian, frictional and rotational interactions in the granular media, as well as the highly discontinuous interfacial effects between the (discrete) granular media and the (continuous) intermediate elastic solid, the nonlinear acoustics of the integrated system is computationally studied subject to a half-sine shock excitation applied to a boundary granule of either the left or right granular medium. The highly discontinuous and nonlinear interaction forces coupling the granular media to the elastic solid are accurately computed through an algorithm with interrelated iteration and interpolation at successive adaptive time steps. Numerical convergence is ensured by monitoring the (linearized) eigenvalues of a nonlinear map of interface forces at each (variable) time step. Due to the strong nonlinearity and hierarchical asymmetry of the left and right granular media, time scale disparity occurs in the response of the interface which breaks acoustic reciprocity. Specifically, depending on the location and intensity of the applied shock, propagating wavefronts are excited in the granular media, which, in turn, excite either (slow) low-frequency vibrations or (fast) high-frequency acoustics in the intermediate elastic medium. This scale disparity is due to the size disparity of the left and right granular media, which yields drastically different wave speeds in the resulting propagating wavefronts. As a result, the continuum part of the interface responds with either low-frequency vibrations-when the shock is applied to the larger-scale granular medium, or high-frequency waves-when the shock is applied to the smaller-scale granular medium. This provides the fundamental mechanism for breaking reciprocity in the interface. The nonreciprocal interfacial acoustics studied here apply to a broad class of asymmetric hybrid (discrete-continuum) nonlinear systems and can inform predictive designs of highly effective granular shock protectors or granular acoustic diodes.

Piezoelectric Impact Energy Harvester Based on the Composite Spherical Particle Chain for Self-Powered Sensors

In this study, a novel piezoelectric energy harvester (PEH) based on the array composite spherical particle chain was constructed and explored in detail through simulation and experimental verification. The power test of the PEH based on array composite particle chains in the self-powered system was realized. Firstly, the model of PEH based on the composite spherical particle chain was constructed to theoretically realize the collection, transformation, and storage of impact energy, and the advantages of a composite particle chain in the field of piezoelectric energy harvesting were verified. Secondly, an experimental system was established to test the performance of the PEH, including the stability of the system under a continuous impact load, the power adjustment under different resistances, and the influence of the number of particle chains on the energy harvesting efficiency. Finally, a self-powered supply system was established with the PEH composed of three composite particle chains to realize the power supply of the microelectronic components. This paper presents a method of collecting impact energy based on particle chain structure, and lays an experimental foundation for the application of a composite particle chain in the field of piezoelectric energy harvesting.

Wireless Module for Nondestructive Testing/Structural Health Monitoring Applications Based on Solitary Waves

In recent years, there has been an increasing interest in the use of highly nonlinear solitary waves (HNSWs) for nondestructive evaluation and structural health monitoring applications. HNSWs are mechanical waves that can form and travel in highly nonlinear systems, such as granular particles in Hertzian contact. The easiest setup consists of a built-in transducer in drypoint contact with the structure or material to be inspected/monitored. The transducer is made of a monoperiodic array of spherical particles that enables the excitation and detection of the solitary waves. The transducer is wired to a data acquisition system that controls the functionality of the transducer and stores the time series for post-processing. In this paper, the design and testing of a wireless unit that enables the remote control of a transducer without the need to connect it to sophisticated test equipment are presented. Comparative tests and analyses between the measurements obtained with the newly designed wireless unit and the conventional wired configuration are provided. The results are corroborated by an analytical model that predicts the dynamic interaction between solitary waves and materials with different modulus. The advantages and limitations of the proposed wireless platform are given along with some suggestions for future developments.

Nonlinear dynamics of coupled transverse-rotational waves in granular chains

The nonlinear dynamics of coupled waves in one-dimensional granular chains with and without a substrate is theoretically studied accounting for quadratic nonlinearity. The multiple time scale method is used to derive the nonlinear dispersion relations for infinite granular chains and to obtain the wave solutions for semi-infinite systems. It is shown that the sum frequency and difference frequency components of the coupled transverse-rotational waves are generated due to their nonlinear interactions with the longitudinal wave. Nonlinear resonances are not present in the chain with no substrate where these frequency components have low amplitudes and exhibit beating oscillations. In the chain positioned on a substrate two types of nonlinear resonances are predicted. At resonance, the fundamental frequency wave amplitudes decrease and the generated frequency component amplitudes increase along the chain, accompanied by the oscillations due to the wave number asynchronism. The results confirm the possibility of a highly efficient energy transfer between the waves of different frequencies, which could find applications in the design of acoustic devices for energy transfer and energy rectification.

Waves in strongly nonlinear discrete systems

The paper presents the main steps in the development of the strongly nonlinear wave dynamics of discrete systems. The initial motivation was prompted by the challenges in the design of barriers to mitigate high-amplitude compression pulses caused by impact or explosion. But this area poses a fundamental mathematical and physical problem and should be considered as a natural step in developing strongly nonlinear wave dynamics. Strong nonlinearity results in a highly tunable behaviour and allows design of systems with properties ranging from a weakly nonlinear regime, similar to the classical case of the Fermi-Pasta-Ulam lattice, or to a non-classical case of sonic vacuum. Strongly nonlinear systems support periodic waves and one of the fascinating results was a discovery of a strongly nonlinear solitary wave in sonic vacuum (a limiting case of a periodic wave) with properties very different from the Korteweg de Vries solitary wave. Shock-like oscillating and monotonous stationary stress waves can also be supported if the system is dissipative. The paper discusses the main theoretical and experimental results, focusing on travelling waves and possible future developments in the area of strongly nonlinear metamaterials.This article is part of the theme issue 'Nonlinear energy transfer in dynamical and acoustical systems'.

Small nanoparticles, surface geometry and contact forces

In this molecular dynamics study, we examine the local surface geometric effects of the normal impact force between two approximately spherical nanoparticles that collide in a vacuum. Three types of surface geometries-(i) crystal facets, (ii) sharp edges, and (iii) amorphous surfaces of small nanoparticles with radii R<10 nm-are considered. The impact forces are compared with their macroscopic counterparts described by nonlinear contact forces based on Hertz contact mechanics. In our simulations, edge and amorphous surface contacts with weak surface energy reveal that the average impact forces are in excellent agreement with the Hertz contact force. On the other hand, facet collisions show a linearly increasing force with increasing compression. Our results suggest that the nearly spherical nanoparticles are likely to enable some nonlinear dynamic phenomena, such as breathers and solitary waves observed in granular materials, both originating from the nonlinear contact force.

Nonlinear coherent structures in granular crystals

The study of granular crystals, which are nonlinear metamaterials that consist of closely packed arrays of particles that interact elastically, is a vibrant area of research that combines ideas from disciplines such as materials science, nonlinear dynamics, and condensed-matter physics. Granular crystals exploit geometrical nonlinearities in their constitutive microstructure to produce properties (such as tunability and energy localization) that are not conventional to engineering materials and linear devices. In this topical review, we focus on recent experimental, computational, and theoretical results on nonlinear coherent structures in granular crystals. Such structures-which include traveling solitary waves, dispersive shock waves, and discrete breathers-have fascinating dynamics, including a diversity of both transient features and robust, long-lived patterns that emerge from broad classes of initial data. In our review, we primarily discuss phenomena in one-dimensional crystals, as most research to date has focused on such scenarios, but we also present some extensions to two-dimensional settings. Throughout the review, we highlight open problems and discuss a variety of potential engineering applications that arise from the rich dynamic response of granular crystals.

Energy decay in a tapped granular column: Can a one-dimensional toy model provide insight into fully three-dimensional systems?

The decay of energy within particulate media subjected to an impulse is an issue of significant scientific interest, but also one with numerous important practical applications. In this paper, we study the dynamics of a granular system exposed to energetic impulses in the form of discrete taps from a solid surface. By considering a one-dimensional toy system, we develop a simple theory, which successfully describes the energy decay within the system following exposure to an impulse. We then extend this theory so as to make it applicable also to more realistic, three-dimensional granular systems, assessing the validity of the model through direct comparison with discrete particle method simulations. The theoretical form presented possesses several notable consequences; in particular, it is demonstrated that for suitably large systems, effects due to the bounding walls may be entirely neglected. We also establish the existence of a threshold system size above which a granular bed may be considered fully three dimensional.

Emergence of dispersive shocks and rarefaction waves in power-law contact models

In the present work, motivated by generalized forms of the Hertzian dynamics associated with granular crystals, we consider the possibility of such models to give rise to both dispersive shock and rarefaction waves. Depending on the value p of the nonlinearity exponent, we find that both of these possibilities are realizable. We use a quasicontinuum approximation of a generalized inviscid Burgers model in order to predict the solution profile up to times near the formation of the dispersive shock, as well as to estimate when it will occur. Beyond that time threshold, oscillations associated with the highly dispersive nature of the underlying model emerge, which cannot be captured by the quasicontinuum approximation. Our analytical characterization of the above features is complemented by systematic numerical computations.

Elastic Vector Solitons in Soft Architected Materials

We demonstrate experimentally, numerically, and analytically that soft architected materials can support the propagation of elastic vector solitons. More specifically, we focus on structures comprising a network of squares connected by thin and highly deformable ligaments and investigate the propagation of planar nonlinear elastic waves. We find that for sufficiently large amplitudes two components-one translational and one rotational-are coupled together and copropagate without dispersion. Our results not only show that soft architected materials offer a new and rich platform to study the propagation of nonlinear waves, but also open avenues for the design of a new generation of smart systems that take advantage of nonlinearities to control and manipulate the propagation of large amplitude vibrations.

Wave propagation in one-dimensional microscopic granular chains

We employ noncontact optical techniques to generate and measure stress waves in uncompressed, one-dimensional microscopic granular chains, and support our experiments with discrete numerical simulations. We show that the wave propagation through dry particles (150 μm radius) is highly nonlinear and it is significantly influenced by the presence of defects (e.g., surface roughness, interparticle gaps, and misalignment). We derive an analytical relation between the group velocity and gap size, and define bounds for the formation of highly nonlinear solitary waves as a function of gap size and axial misalignment.

Stress Wave Propagation in Two-dimensional Buckyball Lattice

Orderly arrayed granular crystals exhibit extraordinary capability to tune stress wave propagation. Granular system of higher dimension renders many more stress wave patterns, showing its great potential for physical and engineering applications. At nanoscale, one-dimensionally arranged buckyball (C60) system has shown the ability to support solitary wave. In this paper, stress wave behaviors of two-dimensional buckyball (C60) lattice are investigated based on square close packing and hexagonal close packing. We show that the square close packed system supports highly directional Nesterenko solitary waves along initially excited chains and hexagonal close packed system tends to distribute the impulse and dissipates impact exponentially. Results of numerical calculations based on a two-dimensional nonlinear spring model are in a good agreement with the results of molecular dynamics simulations. This work enhances the understanding of wave properties and allows manipulations of nanoscale lattice and novel design of shock mitigation and nanoscale energy harvesting devices.

Nonlinear vibrational-state excitation and piezoelectric energy conversion in harmonically driven granular chains

This article explores the excitation of different vibrational states in a spatially extended dynamical system through theory and experiment. As a prototypical example, we consider a one-dimensional packing of spherical particles (a so-called granular chain) that is subject to harmonic boundary excitation. The combination of the multimodal nature of the system and the strong coupling between the particles due to the nonlinear Hertzian contact force leads to broad regions in frequency where different vibrational states are possible. In certain parametric regions, we demonstrate that the nonlinear Schrödinger equation predicts the corresponding modes fairly well. The electromechanical model we apply predicts accurately the conversion from the obtained mechanical energy to the electrical energy observed in experiments.

Scattering of waves by impurities in precompressed granular chains

We study scattering of waves by impurities in strongly precompressed granular chains. We explore the linear scattering of plane waves and identify a closed-form expression for the reflection and transmission coefficients for the scattering of the waves from both a single impurity and a double impurity. For single-impurity chains, we show that, within the transmission band of the host granular chain, high-frequency waves are strongly attenuated (such that the transmission coefficient vanishes as the wavenumber k→±π), whereas low-frequency waves are well-transmitted through the impurity. For double-impurity chains, we identify a resonance-enabling full transmission at a particular frequency-in a manner that is analogous to the Ramsauer-Townsend (RT) resonance from quantum physics. We also demonstrate that one can tune the frequency of the RT resonance to any value in the pass band of the host chain. We corroborate our theoretical predictions both numerically and experimentally, and we directly observe almost complete transmission for frequencies close to the RT resonance frequency. Finally, we show how this RT resonance can lead to the existence of reflectionless modes in granular chains (including disordered ones) with multiple double impurities.

Universal relations for solitary waves in granular crystals under shocks with finite rise and decay times

We focus on solitary waves generated in arrays of lightly contacting spherical elastic granules by shock forces of steep rise and slow decay durations and establish a priori: (i) whether the peak value of the resulting solitary wave would be greater than, equal to, or less than the peak value of the input shock force; (ii) the magnitude of the peak value of the solitary waves; (iii) the magnitude of the linear momentum in each solitary wave; (iv) the magnitude of the linear momentum added to the remaining granules if the first granule is ejected; and (v) a quantitative estimate of the effect of the granules' radius, density, and stiffness on force amplification or mitigation. We have supported the analytical results by direct numerical simulations.

Solitary Wave in One-dimensional Buckyball System at Nanoscale

We have studied the stress wave propagation in one-dimensional (1-D) nanoscopic buckyball (C60) system by molecular dynamics (MD) simulation and quantitative modeling. Simulation results have shown that solitary waves are generated and propagating in the buckyball system through impacting one buckyball at one end of the buckyball chain. We have found the solitary wave behaviors are closely dependent on the initial temperature and impacting speed of the buckyball chain. There are almost no dispersion and dissipation of the solitary waves (stationary solitary wave) for relatively low temperature and high impacting speed. While for relatively high temperature and low impacting speed the profile of the solitary waves is highly distorted and dissipated after propagating several tens of buckyballs. A phase diagram is proposed to describe the effect of the temperature and impacting speed on the solitary wave behaviors in buckyball system. In order to quantitatively describe the wave behavior in buckyball system, a simple nonlinear-spring model is established, which can describe the MD simulation results at low temperature very well. The results presented in this work may lay a solid step towards the further understanding and manipulation of stress wave propagation and impact energy mitigation at nanoscale.

Energy transport in one-dimensional disordered granular solids

We investigate the energy transport in one-dimensional disordered granular solids by extensive numerical simulations. In particular, we consider the case of a polydisperse granular chain composed of spherical beads of the same material and with radii taken from a random distribution. We start by examining the linear case, in which it is known that the energy transport strongly depends on the type of initial conditions. Thus, we consider two sets of initial conditions: an initial displacement and an initial momentum excitation of a single bead. After establishing the regime of sufficiently strong disorder, we focus our study on the role of nonlinearity for both sets of initial conditions. By increasing the initial excitation amplitudes we are able to identify three distinct dynamical regimes with different energy transport properties: a near linear, a weakly nonlinear, and a highly nonlinear regime. Although energy spreading is found to be increasing for higher nonlinearities, in the weakly nonlinear regime no clear asymptotic behavior of the spreading is found. In this regime, we additionally find that energy, initially trapped in a localized region, can be eventually detrapped and this has a direct influence on the fluctuations of the energy spreading. We also demonstrate that in the highly nonlinear regime, the differences in energy transport between the two sets of initial conditions vanish. Actually, in this regime the energy is almost ballistically transported through shocklike excitations.

Superdiffusive transport and energy localization in disordered granular crystals

We study the spreading of initially localized excitations in one-dimensional disordered granular crystals. We thereby investigate localization phenomena in strongly nonlinear systems, which we demonstrate to differ fundamentally from localization in linear and weakly nonlinear systems. We conduct a thorough comparison of wave dynamics in chains with three different types of disorder-an uncorrelated (Anderson-like) disorder and two types of correlated disorders (which are produced by random dimer arrangements)-and for two types of initial conditions (displacement excitations and velocity excitations). We find for strongly precompressed (i.e., weakly nonlinear) chains that the dynamics depend strongly on the type of initial condition. In particular, for displacement excitations, the long-time asymptotic behavior of the second moment m̃(2) of the energy has oscillations that depend on the type of disorder, with a complex trend that differs markedly from a power law and which is particularly evident for an Anderson-like disorder. By contrast, for velocity excitations, we find that a standard scaling m̃(2)∼t(γ) (for some constant γ) applies for all three types of disorder. For weakly precompressed (i.e., strongly nonlinear) chains, m̃(2) and the inverse participation ratio P(-1) satisfy scaling relations m̃(2)∼t(γ) and P(-1)∼t(-η), and the dynamics is superdiffusive for all of the cases that we consider. Additionally, when precompression is strong, the inverse participation ratio decreases slowly (with η<0.1) for all three types of disorder, and the dynamics leads to a partial localization around the core and the leading edge of a propagating wave packet. For an Anderson-like disorder, displacement perturbations lead to localization of energy primarily in the core, and velocity perturbations cause the energy to be divided between the core and the leading edge. This localization phenomenon does not occur in the sonic-vacuum regime, which yields the surprising result that the energy is no longer contained in strongly nonlinear waves but instead is spread across many sites. In this regime, the exponents are very similar (roughly γ≈1.7 and η≈1) for all three types of disorder and for both types of initial conditions.

Conical wave propagation and diffraction in two-dimensional hexagonally packed granular lattices

Linear and nonlinear mechanisms for conical wave propagation in two-dimensional lattices are explored in the realm of phononic crystals. As a prototypical example, a statically compressed granular lattice of spherical particles arranged in a hexagonal packing configuration is analyzed. Upon identifying the dispersion relation of the underlying linear problem, the resulting diffraction properties are considered. Analysis both via a heuristic argument for the linear propagation of a wave packet and via asymptotic analysis leading to the derivation of a Dirac system suggests the occurrence of conical diffraction. This analysis is valid for strong precompression, i.e., near the linear regime. For weak precompression, conical wave propagation is still possible, but the resulting expanding circular wave front is of a nonoscillatory nature, resulting from the complex interplay among the discreteness, nonlinearity, and geometry of the packing. The transition between these two types of propagation is explored.

Mechanisms for impulsive energy dissipation and small-scale effects in microgranular media

We study impulse response in one-dimensional homogeneous microgranular chains on a linear elastic substrate. Microgranular interactions are analytically described by the Schwarz contact model which includes nonlinear compressive as well as snap-to and from-contact adhesive effects forming a hysteretic loop in the force deformation relationship. We observe complex transient dynamics, including disintegration of solitary pulses, local clustering, and low-to-high-frequency energy transfers resulting in enhanced energy dissipation. We study in detail the underlying dynamics of cluster formation in the impulsively loaded medium and relate enhanced energy dissipation to the rate of cluster formation. These unusual and interesting dynamical phenomena are shown to be robust over a range of physically feasible conditions and are solely scale effects since they are attributed to surface forces, which have no effect at the macroscale. We establish a universal relation between the reclustering rate and the effective damping in these systems. Our findings demonstrate that scale effects generating new nonlinear features can drastically affect the dynamics and acoustics of microgranular materials.

Nonlinear low-to-high-frequency energy cascades in diatomic granular crystals

We study wave propagation in strongly nonlinear one-dimensional diatomic granular crystals under an impact load. Depending on the mass ratio of the "light" to "heavy" beads, this system exhibits rich wave dynamics from highly localized traveling waves to highly dispersive waves featuring strong attenuation. We demonstrate experimentally the nonlinear resonant and antiresonant interactions of particles, and we verify that the nonlinear resonance results in strong wave attenuation, leading to highly efficient nonlinear energy cascading without relying on material damping. In this process, mechanical energy is transferred from low to high frequencies, while propagating waves emerge in both ordered and chaotic waveforms via a distinctive spatial cascading. This energy transfer mechanism from lower to higher frequencies and wave numbers is of particular significance toward the design of novel nonlinear acoustic metamaterials with inherently passive energy redistribution properties.

Highly nonlinear wave propagation in elastic woodpile periodic structures

In the present work, we experimentally implement, numerically compute with, and theoretically analyze a configuration in the form of a single column woodpile periodic structure. Our main finding is that a Hertzian, locally resonant, woodpile lattice offers a test bed for the formation of genuinely traveling waves composed of a strongly localized solitary wave on top of a small amplitude oscillatory tail. This type of wave, called a nanopteron, is not only motivated theoretically and numerically, but is also visualized experimentally by means of a laser Doppler vibrometer. This system can also be useful for manipulating stress waves at will, for example, to achieve strong attenuation and modulation of high-amplitude impacts without relying on damping in the system.

A quasi-unidimensional granular chain to attenuate impact

We study via numerical simulations a granular chain not only with decreasing radii (forward tapering) in geometric progression, but also decorated with grains positioned on the top and bottom of the chain, without altering its original length. The decorating grains act as an auxiliary chain which traps part of the energy and linear momentum, both propagating as a pulse due to disturbances produced at the end of the chain. Thus, this configuration optimizes the impact attenuation and aligns the chain naturally (by symmetry) which facilitates the construction of the experimental setup. Furtheremore, since the decorating grain radii increase along the chain, this new type of chain is necessarily short in order to avoid precompression. Nevertheless, even with short chains, it is possible to mitigate impacts almost completely.

Wave tailoring by precompression in confined granular systems

We present a granular system whose response under an impact load can be varied from rapidly decaying to almost constant amplitude waves by an external regulator. The system consists of a granular chain of larger spheres surrounded by small spheres, confined in a hollow cylindrical tube and supporting wave propagation along the axis of the cylinder. We demonstrate using numerical simulations that the response can be controlled by applying radial precompression. These observations are then complemented by an asymptotic analysis, which shows that the decay in the leading wave is due to energy leakage to the oscillating small beads in the tail of the wave. This system has potential applications in systems requiring tuning of elastic waves.

Characterizing traveling-wave collisions in granular chains starting from integrable limits: the case of the Korteweg-de Vries equation and the Toda lattice

Our aim in the present work is to develop approximations for the collisional dynamics of traveling waves in the context of granular chains in the presence of precompression. To that effect, we aim to quantify approximations of the relevant Hertzian FPU-type lattice through both the Korteweg-de Vries (KdV) equation and the Toda lattice. Using the availability in such settings of both one-soliton and two-soliton solutions in explicit analytical form, we initialize such coherent structures in the granular chain and observe the proximity of the resulting evolution to the underlying integrable (KdV or Toda) model. While the KdV offers the possibility to accurately capture collisions of solitary waves propagating in the same direction, the Toda lattice enables capturing both copropagating and counterpropagating soliton collisions. The error in the approximation is quantified numerically and connections to bounds established in the mathematical literature are also given.

Traveling waves in 2D hexagonal granular crystal lattices

This study describes the dynamic response of a two-dimensional hexagonal packing of uncompressed stainless steel spheres excited by localized impulsive loadings. The dynamics of the system are modeled using the Hertzian normal contact law. After the initial impact strikes the system, a characteristic wave structure emerges and continuously decays as it propagates through the lattice. Using an extension of the binary collision approximation for one-dimensional chains, we predict its decay rate, which compares well with numerical simulations and experimental data. While the hexagonal lattice does not support constant speed traveling waves, we provide scaling relations that characterize the directional power law decay of the wave velocity for various angles of impact. Lastly, we discuss the effects of weak disorder on the directional amplitude decay rates.

Damped-driven granular chains: an ideal playground for dark breathers and multibreathers

By applying an out-of-phase actuation at the boundaries of a uniform chain of granular particles, we demonstrate experimentally that time-periodic and spatially localized structures with a nonzero background (so-called dark breathers) emerge for a wide range of parameter values and initial conditions. We demonstrate a remarkable control over the number of breathers within the multibreather pattern that can be "dialed in" by varying the frequency or amplitude of the actuation. The values of the frequency (or amplitude) where the transition between different multibreather states occurs are predicted accurately by the proposed theoretical model, which is numerically shown to support exact dark breather and multibreather solutions. Moreover, we visualize detailed temporal and spatial profiles of breathers and, especially, of multibreathers using a full-field probing technology and enable a systematic favorable comparison among theory, computation, and experiments. A detailed bifurcation analysis reveals that the dark and multibreather families are connected in a "snaking" pattern, providing a roadmap for the identification of such fundamental states and their bistability in the laboratory.

Directed ratchet transport in granular chains

Directed-ratchet transport (DRT) in a one-dimensional lattice of spherical beads, which serves as a prototype for granular chains, is investigated. We consider a system where the trajectory of the central bead is prescribed by a biharmonic forcing function with broken time-reversal symmetry. By comparing the mean integrated force of beads equidistant from the forcing bead, two distinct types of directed transport can be observed-spatial and temporal DRT. Based on the value of the frequency of the forcing function relative to the cutoff frequency, the system can be categorized by the presence and magnitude of each type of DRT. Furthermore, we investigate and quantify how varying additional parameters such as the biharmonic weight affects DRT velocity and magnitude. Finally, friction is introduced into the system and is found to significantly inhibit spatial DRT. In fact, for sufficiently low forcing frequencies, the friction may even induce a switching of the DRT direction.

Frequency bands of strongly nonlinear homogeneous granular systems

Recent numerical studies on an infinite number of identical spherical beads in Hertzian contact showed the presence of frequency bands [Jayaprakash, Starosvetsky, Vakakis, Peeters, and Kerschen, Nonlinear Dyn. 63, 359 (2011)]. These bands, denoted here as propagation and attenuation bands (PBs and ABs), are typically present in linear or weakly nonlinear periodic media; however, their counterparts are not intuitive in essentially nonlinear periodic media where there is a complete lack of classical linear acoustics, i.e., in "sonic vacua." Here, we study the effects of PBs and ABs on the forced dynamics of ordered, uncompressed granular systems. Through numerical and experimental techniques, we find that the dynamics of these systems depends critically on the frequency and amplitude of the applied harmonic excitation. For fixed forcing amplitude, at lower frequencies, the oscillations are large in amplitude and governed by strongly nonlinear and nonsmooth dynamics, indicating PB behavior. At higher frequencies the dynamics is weakly nonlinear and smooth, in the form of compressed low-amplitude oscillations, indicating AB behavior. At the boundary between the PB and the AB large-amplitude oscillations due to resonance occur, giving rise to collisions between beads and chaotic dynamics; this renders the forced dynamics sensitive to initial and forcing conditions, and hence unpredictable. Finally, we study asymptotically the near field standing wave dynamics occurring for high frequencies, well inside the AB.

Interaction of traveling waves with mass-with-mass defects within a Hertzian chain

We study the dynamic response of a granular chain of particles with a resonant inclusion (i.e., a particle attached to a harmonic oscillator, or a mass-with-mass defect). We focus on the response of granular chains excited by an impulse, with no static precompression. We find that the presence of the harmonic oscillator can be used to tune the transmitted and reflected energy of a mechanical pulse by adjusting the ratio between the harmonic resonator mass and the bead mass. Furthermore, we find that this system has the capability of asymptotically trapping energy, a feature that is not present in granular chains containing other types of defects. Finally, we study the limits of low and high resonator mass, and the structure of the reflected and transmitted pulses.

Dark breathers in granular crystals

We present a study of the existence, stability, and bifurcation structure of families of dark breathers in a one-dimensional uniform chain of spherical beads under static load. A defocusing nonlinear Schrödinger equation (NLS) is derived for frequencies that are close to the edge of the phonon band and is used to construct targeted initial conditions for numerical computations. Salient features of the system include the existence of large amplitude solutions that emerge from the small amplitude solutions described by the NLS equation, and the presence of a nonlinear instability that, to the best of the authors' knowledge, has not been observed in classical Fermi-Pasta-Ulam lattices. Finally, it is also demonstrated that these dark breathers can be detected in a physically realistic experimental settings by merely actuating the ends of an initially at rest chain of beads and inducing destructive interference between their signals.

A comparative study on three different transducers for the measurement of nonlinear solitary waves

In the last decade there has been an increasing interest in the use of highly- and weakly- nonlinear solitary waves in engineering and physics. Nonlinear solitary waves can form and travel in nonlinear systems such as one-dimensional chains of particles, where they are conventionally generated by the mechanical impact of a striker and are measured either by using thin transducers embedded in between two half-particles or by a force sensor placed at the chain's base. These waves have a constant spatial wavelength and their speed, amplitude, and duration can be tuned by modifying the particles' material or size, or the velocity of the striker. In this paper we propose two alternative sensing configurations for the measurements of solitary waves propagating in a chain of spherical particles. One configuration uses piezo rods placed in the chain while the other exploits the magnetostrictive property of ferromagnetic materials. The accuracy of these two sensing systems on the measurement of the solitary wave's characteristics is assessed by comparing experimental data to the numerical prediction of a discrete particle model and to the experimental measurements obtained by means of a conventional transducer. The results show very good agreement and the advantages and limitations of the new sensors are discussed.

Wave propagation in square granular crystals with spherical interstitial intruders

We investigate the propagation and scattering of highly nonlinear waves in granular systems composed of spheres in contact arranged in a square packing, and study how the presence of small and light spherical interstitial defects, also referred to as intruders, affects the wave propagation. The effects of a single defect are investigated experimentally and compared to numerical simulations, showing very good quantitative agreement. Transmitted and scattered waves are formed, whose characteristics depend on the material properties of the defect in relation to the properties of the particles in the lattice. Experiments and numerical simulations reveal that stiffer defects are more efficient at redistributing energy outside the impacted chain and soft defects induce a localization of the energy at the defect. Finally, the effects of the presence of two defects, placed diagonally or aligned in the square packing are also investigated, as well as how their interaction depends on their relative positions.

Breathers in periodic granular chains with multiple band gaps

We consider the localized nonlinear breathing modes that emerge in heterogeneous granular configurations of two materials with a periodicity of three and four beads. We examine as characteristic examples chains with one steel and two aluminum beads, as well as ones with one steel and three aluminum beads, although we briefly touch upon other configurations as well, illustrating their similar characteristics. We analyze the higher order gaps that emerge in such settings and explore the intrinsic localized modes that bifurcate from the edge of the upper bands. A generic surprising feature of such states is that they appear to be more robust than their counterparts bifurcating from the edges of the lower bands. Direct numerical simulations, using driving of the system at suitable frequencies through an actuator or taking advantage of the modulational instabilities of extended band edge states in the system, illustrate the spontaneous formation of localized modes within the corresponding nearest gaps. In these numerical experiments, we also account for the presence of dissipation and analyze its potential role.

Site-Specific Quantification of Bone Quality Using Highly Nonlinear Solitary Waves

Osteoporosis is a well recognized problem affecting millions of individuals worldwide. The ability to diagnose problems in an effective, efficient, and affordable manner and identify individuals at risk is essential. Site-specific assessment of bone mechanical properties is necessary, not only in the process of fracture risk assessment, but may also be desirable for other applications, such as making intraoperative decisions during spine and joint replacement surgeries. The present study evaluates the use of a one-dimensional granular crystal sensor to measure the elastic properties of bone at selected locations via direct mechanical contact. The granular crystal is composed of a tightly packed chain of particles that interact according to the Hertzian contact law. Such chains represent one of the simplest systems to generate and propagate highly nonlinear acoustic signals in the form of compact solitary waves. First, we investigated the sensitivity of the sensor to known variations in bone density using a synthetic cancellous bone substitute, representing clinical bone quality ranging from healthy to osteoporotic. Once the relationship between the signal response and known bone properties was established, the sensor was used to assess the bone quality of ten human cadaveric specimens. The efficacy and accuracy of the sensor was then investigated by comparing the sensor measurements with the bone mineral density (BMD) obtained using dual-energy x-ray absorptiometry (DEXA). The results indicate that the proposed technique is capable of detecting differences in bone quality. The ability to measure site-specific properties without exposure to radiation has the potential to be further developed for clinical applications.

Wave propagation in random granular chains

The influence of randomness on wave propagation in one-dimensional chains of spherical granular media is investigated. The interaction between the elastic spheres is modeled using the classical Hertzian contact law. Randomness is introduced in the discrete model using random distributions of particle mass, Young's modulus, or radius. Of particular interest in this study is the quantification of the attenuation in the amplitude of the impulse associated with various levels of randomness: two distinct regimes of decay are observed, characterized by an exponential or a power law, respectively. The responses are normalized to represent a vast array of material parameters and impact conditions. The virial theorem is applied to investigate the transfer from potential to kinetic energy components in the system for different levels of randomness. The level of attenuation in the two decay regimes is compared for the three different sources of randomness and it is found that randomness in radius leads to the maximum rate of decay in the exponential regime of wave propagation.

Defect modes in one-dimensional granular crystals

We study the vibrational spectra of one-dimensional statically compressed granular crystals (arrays of elastic particles in contact) containing light-mass defects. We focus on the prototypical settings of one or two spherical defects (particles of smaller radii) interspersed in a chain of larger uniform spherical particles. We present a systematic measurement, using continuous noise, of the near-linear frequency spectrum within the spatial vicinity of the defect(s). Using this technique, we identify the frequencies of the localized defect modes as a function of the defect size and the position of the defects relative to each other. We also compare the experimentally determined frequencies with those obtained by numerical eigenanalysis and by analytical expressions based on few-site considerations. These approximate analytical expressions, based on normal-mode analysis, are found to be in excellent agreement with numerics for a wide range of mass ratios. We also observe that the experimentally measured frequencies of the localized defect modes are uniformly upshifted, compared to the numerically and theoretically predicted values.

Highly nonlinear solitary wave propagation in Y-shaped granular crystals with variable branch angles

We study the propagation of highly nonlinear waves in a branched (Y-shaped) granular crystal composed of chains of spherical particles of different materials, arranged at variable branch angles. We experimentally test the dynamic behavior of a solitary pulse, or of a train of solitary waves, crossing the Y-junction interface, and splitting between the two branches. We describe the dependence of the split pulses' speed and amplitude on the branch angles. Analytic predictions based on the quasiparticle model and numerical simulations based on Hertzian interactions between the particles are found to be in excellent agreement with the experimental data.

Mesoscopic approach to granular crystal dynamics

We present a mesoscopic approach to granular crystal dynamics, which comprises a three-dimensional finite-element model and a one-dimensional regularized contact model. The approach investigates the role of vibrational-energy trapping effects in the dynamic behavior of one-dimensional chains of particles in contact (i.e., granular crystals), under small to moderate impact velocities. The only inputs of the models are the geometry and the elastic material properties of the individual particles that form the system. We present detailed verification results and validate the model comparing its predictions with experimental data. This approach provides a physically sound, first-principles description of dissipative losses in granular systems.

Sustained strong fluctuations in a nonlinear chain at acoustic vacuum: beyond equilibrium

Here we consider dynamical problems as in linear response theory but for purely nonlinear systems where acoustic propagation is prohibited by the potential, e.g., the case of an alignment of elastic grains confined between walls. Our simulations suggest that in the absence of acoustic propagation, the system relaxes using only solitary waves and the eventual state does not resemble an equilibrium state. Further, the studies reveal that multiple perturbations could give rise to hot and cold spots in these systems. We first use particle dynamics based simulations to understand how one of the two unequal colliding solitary waves in the chain can gain energy. Specifically, we find that for head-on collisions the smaller wave gains energy, whereas when a more energetic wave overtakes a less energetic wave, the latter gains energy. The balance between the rate at which the solitary waves break down and the rate at which they grow eventually makes it possible for the system to reach a peculiar equilibriumlike phase that is characteristic of these purely nonlinear systems. The study of the features and the robustness of the fluctuations in time has been addressed next. A particular characteristic of this equilibriumlike or quasiequilibrium phase is that very large energy fluctuations are possible--and by very large, we mean that the energy can vary between zero and several times the average energy per grain. We argue that the magnitude of the fluctuations depend on the nature of the nonlinearity in the potential energy function and the feature that any energy must eventually travel as a compact solitary wave in these systems where the solitary wave energies may vary widely. In closing we address whether these fluctuations are peculiar to one dimension or can exist in higher dimensions. The study hence raises the following intriguing possibility. Are there physical or biological systems where these kinds of nonlinear forces exist, and if so, can such large fluctuations actually be seen? Implications of the study are briefly discussed.

Highly nonlinear solitary waves in chains of ellipsoidal particles

We study the dynamic response of a one-dimensional chain of ellipsoidal particles excited by a single compressive impulse. We detail the Hertzian contact theory describing the interaction between two ellipsoidal particles under compression, and use it to model the dynamic response of the system. We observe the formation of highly nonlinear solitary waves in the chain, and we also study their propagation properties. We measure experimentally the traveling pulse amplitude (force), the solitary wave speed, and the solitary wave width. We compare these results with theoretical predictions in the long wavelength approximation, and with numerical results obtained with a discrete particle model and with finite element simulations. We also study the propagation of highly nonlinear solitary waves in the chain with particles arranged in different configurations to show the effects of the particle's geometry on the wave propagation characteristics and dissipation. We find very good agreement between experiment, theory, and simulations for all the ranges of impact velocity and particle arrangements investigated.

Laser-based excitation of nonlinear solitary waves in a chain of particles

Highly nonlinear solitary waves (HNSWs) are stress waves that can form and travel in highly nonlinear systems. They are characterized by a constant spatial wavelength and by a tunable propagation speed, dependent on the wave amplitude. Conventionally, HNSW's are generated in one-dimensional chains of spherical particles by means of a mechanical impact. In this paper, we demonstrate that short-duration laser pulses can be used to generate HNSW's, and we characterize their propagating properties in terms of shape, speed, and duration. We compare the waves' characteristics with theoretical predictions, finding excellent agreement. In addition a simplified formulation is given to estimate the dynamic contact force generated by laser pulses onto the chain.

Bifurcation-based acoustic switching and rectification

Switches and rectification devices are fundamental components used for controlling the flow of energy in numerous applications. Thermal and acoustic rectifiers have been proposed for use in biomedical ultrasound applications, thermal computers, energy- saving and -harvesting materials, and direction-dependent insulating materials. In all these systems the transition between transmission states is smooth with increasing signal amplitudes. This limits their effectiveness as switching and logic devices, and reduces their sensitivity to external conditions as sensors. Here we overcome these limitations by demonstrating a new mechanism for tunable rectification that uses bifurcations and chaos. This mechanism has a sharp transition between states, which can lead to phononic switching and sensing. We present an experimental demonstration of this mechanism, applied in a mechanical energy rectifier operating at variable sonic frequencies. The rectifier is a granular crystal, composed of a statically compressed one-dimensional array of particles in contact, containing a light mass defect near a boundary. As a result of the defect, vibrations at selected frequencies cause bifurcations and a subsequent jump to quasiperiodic and chaotic states with broadband frequency content. We use this combination of frequency filtering and asymmetrically excited bifurcations to obtain rectification ratios greater than 10(4). We envisage this mechanism to enable the design of advanced photonic, thermal and acoustic materials and devices.

Actuators for the generation of highly nonlinear solitary waves

In this paper we present the design of two actuators for the generation of highly nonlinear solitary waves (HNSWs), which are mechanical waves that can form and travel in highly nonlinear systems. These waves are characterized by a constant spatial wavelength and by a tunable propagation speed, dependent on the wave amplitude. To date, the simplest and widely adopted method to generate HNSWs is by impacting a striker onto a chain of beads of equal size and mass. This operation is conducted manually and it might be impracticable if repetition rates higher than 0.1 Hz are necessary. It is known that the HNSWs' properties, such as amplitude, duration, and speed can be modified by changing the size or the material of the particles, the velocity of the striker, and/or the precompression on the chain. To address the limitations associated with the manual generation of HNSWs we designed, built, and tested two actuators. The first actuator consists of a chain of particles wrapped by an electromagnet that induces static precompression on the chain. This design allows for the generation of solitary waves with controlled properties. The second actuator consists of a chain surmounted by an electromagnet that lifts and releases a striker. This actuator permits the remote and noncontact generation of solitary waves. The performance of both actuators is evaluated by comparing the experimental HNSWs to theoretical predictions, based on the long wavelength approximation.

Intrinsic energy localization through discrete gap breathers in one-dimensional diatomic granular crystals

We present a systematic study of the existence and stability of discrete breathers that are spatially localized in the bulk of a one-dimensional chain of compressed elastic beads that interact via Hertzian contact. The chain is diatomic, consisting of a periodic arrangement of heavy and light spherical particles. We examine two families of discrete gap breathers: (1) an unstable discrete gap breather that is centered on a heavy particle and characterized by a symmetric spatial energy profile and (2) a potentially stable discrete gap breather that is centered on a light particle and is characterized by an asymmetric spatial energy profile. We investigate their existence, structure, and stability throughout the band gap of the linear spectrum and classify them into four regimes: a regime near the lower optical band edge of the linear spectrum, a moderately discrete regime, a strongly discrete regime that lies deep within the band gap of the linearized version of the system, and a regime near the upper acoustic band edge. We contrast discrete breathers in anharmonic Fermi-Pasta-Ulam (FPU)-type diatomic chains with those in diatomic granular crystals, which have a tensionless interaction potential between adjacent particles, and note that the asymmetric nature of the tensionless interaction potential can lead to hybrid bulk-surface localized solutions.

Highly nonlinear pulse splitting and recombination in a two-dimensional granular network

The propagation of highly nonlinear signals in a branched two-dimensional granular system was investigated experimentally and numerically for a system composed of chains of spherical beads of different materials. The system studied consists of a double Y-shaped guide in which high- and low-modulus/mass chains of spheres are arranged in various geometries. We observed the transformation of a single or a train of solitary pulses crossing the interface between branches. We report fast splitting of the initial pulse, rapid chaotization of the signal and impulse redirection and bending. Pulse and energy trapping was also observed in the branches. Numerical analysis based on Hertzian interaction between the particles and the side walls of the guide was found in agreement with the experimental data, except for nonsymmetric arrangements of particles excited by a large mass striker.

Nonlinear waves in disordered diatomic granular chains

We investigate the propagation and scattering of highly nonlinear waves in disordered granular chains composed of diatomic (two-mass) units of spheres that interact via Hertzian contact. Using ideas from statistical mechanics, we consider each diatomic unit to be a "spin," so that a granular chain can be viewed as a spin chain composed of units that are each oriented in one of two possible ways. Experiments and numerical simulations both reveal the existence of two different mechanisms of wave propagation: in low-disorder chains, we observe the propagation of a solitary pulse with exponentially decaying amplitude. Beyond a critical level of disorder, the wave amplitude instead decays as a power law, and the wave transmission becomes insensitive to the level of disorder. We characterize the spatiotemporal structure of the wave in both propagation regimes and propose a simple theoretical interpretation for a transition between the two regimes. Our investigation suggests that an elastic spin chain can be used as a model system to investigate the role of heterogeneities in the propagation of highly nonlinear waves.