Finite Element Studies of Solitary Waves in Granular Chains

Abstract

Solitary wave propagation in a monodisperse metallic granular chain was simulated using the finite element method. The model was built to address a discrepancy between numerical and experimental results from Lazaridi and Nesterenko.¹ In their work, solitary waves were generated in a chain of particles through impact of a piston, and results were quantified by comparing the chains' reactions to a rigid wall. Their numerical calculations resulted in a solitary wave with a force amplitude of 83 N, while it was measured experimentally to be 71 N. In the present work, the configuration of the granular chain and piston was duplicated from Lazaridi and Nesterenko.¹ Qualitatively similar solitary waves were produced, and von Mises stress values indicated that localized plastic deformation is possible, even at low piston impact velocities. These results show that localized plastic deformation was a likely source of dissipation in experiments performed by Lazaridi and Nesterenko. Solitary wave response was investigated in the same metallic granular chain-piston system using LS-DYNA. A power-law hardening material model was used to show that localized plastic deformation is present in the metallic granular chain, even for an impact velocity of 0.5 m/s. This loss due to plastic deformation was quantified via impulse, and it was shown that the loss scales nearly linearly with impact velocity. Therefore, metallic grains may not be suitable for devices that require high amplitude solitary waves. There would be too much energy lost to plastic deformation. The response of an aluminum oxide granular chain was subsequently compared to that of a steel chain because ceramics are inherently elastic. It was shown that solitary waves travel faster and the initial peak is slightly lower when compared to a steel chain. The response of granular chains to impulse loading was investigated as a function of material properties. COMSOL Multiphysics was used to study the effect of density and elastic modulus on a granular chain with fixed Poisson's ratio. Solitary wave velocity and amplitude increased with elastic modulus. Increasing density caused a decrease in wave velocity and an increase in amplitude. In addition, higher density granular chains exhibited a decrease in the number of solitary waves in their respective solitary wave trains. LS-DYNA was then used to explore the response of a variety of ceramic and metallic granular chains. Density, elastic modulus, and Poisson's ratio were all set to representative values for the respective material. It was shown that solitary wave development and decay occur at different rates for different materials. In addition, the kinetic energy decay of the impactor was slower for glass compared with tungsten. Finally, it was shown that a single solitary wave with no train could be produced by impacting a high density, high modulus chain such as tungsten with a glass piston, which has relatively low density and elastic modulus. increasing impact velocity for this case resulted in a single high-amplitude solitary wave with no train.