scispace - formally typeset
Search or ask a question
Journal ArticleDOI

Frequency filtering in disordered granular chains

28 Jun 2014-Acta Mechanica (Springer)-Vol. 225, Iss: 8, pp 2385-2407
TL;DR: In this article, disorder-induced frequency filtering is studied for one-dimensional systems composed of random, pre-stressed masses interacting through both linear and nonlinear (Hertzian) repulsive forces.
Abstract: The study of disorder-induced frequency filtering is presented for one-dimensional systems composed of random, pre-stressed masses interacting through both linear and nonlinear (Hertzian) repulsive forces. An ensemble of such systems is driven at a specified frequency, and the spectral content of the propagated disturbance is examined as a function of distance from the source. It is shown that the transmitted signal contains only low-frequency components, and the attenuation is dependent on the magnitude of disorder, the input frequency, and the contact model. It is found that increased disorder leads to a narrower bandwidth of transmitted frequencies at a given distance from the source and that lower input frequencies exhibit less sensitivity to the arrangement of the masses. Comparison of the nonlinear and linear contact models reveals qualitatively similar filtering behavior; however, it is observed that the nonlinear chain produces transmission spectrums with a greater density at the lowest frequencies. In addition, it is shown that random masses sampled from normal, uniform, and binary distributions produce quantitatively indistinguishable filtering behavior, suggesting that knowledge of only the distribution’s first two moments is sufficient to characterize the bulk signal transmission behavior. Finally, we examine the wave number evolution of random chains constrained to move between fixed end-particles and present a transfer matrix theory in wave number space, and an argument for the observed filtering based on the spatial localization of the higher-frequency normal modes.

Content maybe subject to copyright    Report

Citations
More filters
Journal ArticleDOI
TL;DR: In this paper, experimental data on a model soil in a cubical cell are compared with both discrete element (DEM) simulations and continuum analyses and the results show that 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.
Abstract: 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.

43 citations


Cites background from "Frequency filtering in disordered g..."

  • ...This confirms that frequency filtering occurs as a function of the wave travelling through a granular system as noted by [49,50] and others....

    [...]

Journal ArticleDOI
TL;DR: This work studies the spreading of initially localized excitations in one-dimensional disordered granular crystals to investigate localization phenomena in strongly nonlinear systems, which it is demonstrated to differ fundamentally from localization in linear and weakly non linear systems.
Abstract: 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.

40 citations

Journal ArticleDOI
TL;DR: In this paper, a discrete element method (DEM) model of a face-centred cubic packing of uniform spheres was created to simulate bender element tests to investigate this test from a fundamental perspective.
Abstract: While bender element testing is now well-established as a laboratory technique to determine soil stiffness, a robust technique to interpret the data remains elusive. A discrete element method (DEM) model of a face-centred cubic packing of uniform spheres was created to simulate bender element tests to investigate this test from a fundamental perspective. During the DEM simulations transmitter and receiver signals were recorded, analogous to the data available in laboratory tests, and these macro-scale data were supplemented with particle scale measurements (forces, stresses and displacements). A range of approaches previously applied in experimental and numerical studies were used to analyse the resulting data in both the time and frequency domains. The shortcomings in these approaches are clear from the differences in the resultant shear stiffness values and the frequency-dependent nature of the values. The particle-scale data enabled visualization of the passage of the wave through the sample, and it was found not to be possible to precisely link the arrival of the shear wave at the receiver and any of the previously proposed characteristic points along the signal recorded at the receiver. The most reliable determination of the shear wave velocity was obtained by applying a two-dimensional fast Fourier transform (2D FFT) to the data describing the velocity of the particles lying between the transmitter and receiver elements. Use of the DEM model and this 2D FFT approach facilitated the sensitivity of the system response to small variations in the interparticle force–displacement law (the contact model) to be established.

38 citations

Book ChapterDOI
01 Jan 1972

32 citations

Journal ArticleDOI
TL;DR: In this article, discrete element method (DEM) simulations of planar compression wave propagation were performed to generate the data for the study, and the assembly stiffness and mass matrices were extracted from the DEM model and these data were used in an eigenmode analysis that provided significant insight into the observed overall dynamic response.
Abstract: Laboratory geophysics tests including bender elements and acoustic emission measure the speed of propagation of stress or sound waves in granular materials to derive elastic stiffness parameters. This contribution builds on earlier studies to assess whether the received signal characteristics can provide additional information about either the material’s behaviour or the nature of the material itself. Specifically it considers the maximum frequency that the material can transmit; it also assesses whether there is a simple link between the spectrum of the received signal and the natural frequencies of the sample. Discrete element method (DEM) simulations of planar compression wave propagation were performed to generate the data for the study. Restricting consideration to uniform (monodisperse) spheres, the material fabric was varied by considering face-centred cubic lattice packings as well as random configurations with different packing densities. Supplemental analyses, in addition to the DEM simulations, were used to develop a more comprehensive understanding of the system dynamics. The assembly stiffness and mass matrices were extracted from the DEM model and these data were used in an eigenmode analysis that provided significant insight into the observed overall dynamic response. The close agreement of the wave velocities estimated using eigenmode analysis with the DEM results confirms that DEM wave propagation simulations can reliably be used to extract material stiffness data. The data show that increasing either stress or density allows higher frequencies to propagate through the media, but the low-pass wavelength is a function of packing density rather than stress level. Prior research which had hypothesised that there is a simple link between the spectrum of the received signal and the natural sample frequencies was not substantiated.

24 citations


Cites background from "Frequency filtering in disordered g..."

  • ...The FCC sample consisted of 3200 particles (4 × 4 × 200 layers) and so is equivalent to that considered by Mouraille et al. [19] and Mouraille and Luding [12]; it was created by considering the lattice geometry of the packing....

    [...]

  • ...Lawney and Luding [15] examined a 1-...

    [...]

  • ...Santamarina and Aloufi [10] and Santamarina et al. [11] related the maximum transmitted frequency ( flow−pass) and the associated wavelength (λlow−pass) to particle size, while Mouraille and Luding [12] related λlow−pass to the layer spacing....

    [...]

  • ...Santamarina et al. [11] and Santamarina and Aloufi [10] assumed the particle diameter to be an internal scale ( α) of granular materials where λlow−pass = 2α, while Mouraille and Luding [12] took α to be the layer distance for a FCC sample, i.e. α = √2R....

    [...]

  • ...D chain model, Lawney and Luding [15] showed that the low-frequency eigenmodes are not affected by small random variations in particle mass. Somfai et al. [30] considered a 2-...

    [...]

References
More filters
Journal ArticleDOI
TL;DR: In this paper, a universal relation between the wavefront speed and the force amplitude is derived, valid for all the granular structures studied, including hexagonal packing, face-centered cubic and bodycentered cubic packings.
Abstract: Plane wave propagation in periodic ordered granular media comprising of elastic spherical particles is investigated. The spheres are under zero precompression and are assumed to interact via the Hertzian contact potential. Various two- and three-dimensional granular structures such as hexagonal packing (2D and 3D), face-centered cubic and body-centered cubic packings are considered in the present study, with the plane impact either normal or oblique to the granular system. For the normal impact case, 1D chains equivalent to the 2D and 3D structures are obtained. A universal relation between the wavefront speed and the force amplitude is derived, valid for all the granular structures studied. In the angular impact case, the shear component of the amplitude of the particle velocity is found to initially decay exponentially and further in a series of linear regimes. By employing simpler models, semi-analytical predictions are obtained for the decay of shearing effect.

25 citations


"Frequency filtering in disordered g..." refers background in this paper

  • ...The interest in non-linear, realistic wavepropagation phenomena is documented by the most recent publications on contact dissipation mechanisms in one-dimensional [31,36,38] and higher dimensional idealized granular systems [7,24,31]....

    [...]

Journal ArticleDOI
TL;DR: A theory for the optimal distribution of particle masses which maximize the energy transfer along the chain and it is found that for epsilon=const, the mass distribution is a monotonous function which does not depend on the value of ePSilon.
Abstract: The transmission of kinetic energy through chains of inelastically colliding spheres is investigated for the case of constant coefficient of restitution epsilon=const and impact-velocity-dependent coefficient epsilon(v) for viscoelastic particles. We derive a theory for the optimal distribution of particle masses which maximize the energy transfer along the chain and check it numerically. We found that for epsilon=const, the mass distribution is a monotonous function which does not depend on the value of epsilon. In contrast, for epsilon(v) the mass distribution reveals a pronounced maximum, depending on the particle properties and on the chain length. The system investigated demonstrates that even for small and simple systems, the velocity dependence of the coefficient of restitution may lead to new effects with respect to the same systems under the simplifying approximation epsilon=const.

23 citations


"Frequency filtering in disordered g..." refers background in this paper

  • ...One may imagine applications exploiting disorder effects, as could be pursued in the spirit of [5,10,13,14,32,40,47,49] where such systems were engineered to produce a desired output....

    [...]

  • ...Studies have employed smoothly varying mass distributions [40], “decoration” (e....

    [...]

Journal ArticleDOI
TL;DR: In this article, a mixed numerical/analytical correction algorithm was developed to estimate the velocity of the particles during pulse propagation, which is in turn able to correctly predict the momentum and kinetic energy along the chain for several tapering configurations.
Abstract: We study momentum and energy propagation in 1D tapered chains of spherical granules which interact according to a Hertz potential. In this work we apply the binary collision approximation, which is based on the assumption that transfer of energy along the chain occurs via a succession of two-particle collisions. Although the binary theory correctly captures the trends of increase or decrease of kinetic energy and momentum, the actual values of these quantities are not in good quantitative agreement with those obtained by numerically integrating the full equations of motion. To address this difficulty we have developed a mixed numerical/analytical correction algorithm to provide an improved estimate of the velocity of the particles during pulse propagation. With this corrected velocity we are in turn able to correctly predict the momentum and kinetic energy along the chain for several tapering configurations, specifically for forward linear, forward exponential, backward linear and backward exponential tapering.

16 citations


"Frequency filtering in disordered g..." refers background in this paper

  • ..., deliberate insertion of different sized masses) [10,13,14,27] tapering [5,32,47,49] and controlled variation of the particle material [3,15]....

    [...]

Journal ArticleDOI
TL;DR: In this article, the authors numerically solved the elastic multiple impact problem in a linear chain of balls using the impulse momentum and Hertz contact theory based methods and developed an energy absorption scheme using an arrangement of balls with different sizes.
Abstract: In this paper, we numerically solved the elastic multiple impact problem in a linear chain of balls using the impulse momentum and Hertz contact theory based methods. The first method depends on a parameter called the impulse correlation ratio (ICR) while the latter one depends on the material properties and geometries of the colliding bodies. We compared the post impact velocities from the two methods and we found that the two methods yield similar solutions for elastic collisions. Then, we develop an energy absorption scheme using an arrangement of balls with different sizes. The impulse momentum method is used and the ICRs are experimentally estimated. We use numerical and experimental analyses to demonstrate that one can significantly improve the energy absorption by placing small balls in specific locations in a linear chain of large balls.

16 citations


"Frequency filtering in disordered g..." refers background in this paper

  • ...One may imagine applications exploiting disorder effects, as could be pursued in the spirit of [5,10,13,14,32,40,47,49] where such systems were engineered to produce a desired output....

    [...]

  • ..., deliberate insertion of different sized masses) [10,13,14,27] tapering [5,32,47,49] and controlled variation of the particle material [3,15]....

    [...]

Journal ArticleDOI
TL;DR: In this paper, the average strain theory was used to analyze the behavior of a dense granular aggregate made by identical, elastic spheres, uni-axially compressed at constant pressure.
Abstract: We analyze the behavior of a dense granular aggregate made by identical, elastic spheres, uni-axially compressed at constant pressure. Our goal is to predict the evolution of the effective moduli along the loading path when small perturbations are applied to stressed states. The analytical model is based upon the average strain theory. We show that the moduli in the anisotropic state normalized with the corresponding initial isotropic value, are captured by a so crude model. Numerical simulations support this result.

15 citations