Time domain analysis of locally resonant elastic metamaterials under impact

The microstructure of a material can be engineered to achieve unique properties not found in nature. Microstructured materials, also known as metamaterials (MMs), can exhibit properties utilizing local resonance and dynamics of their heterogeneous microstructure that are activated below the traditional Bragg limit. In this study, the linear dynamic response of a low-frequency resonant ceramic MM slab is analyzed using the Finite Element Method (FEM) in the time domain. The MM is compared to monolithic slabs and other microstructured designs in terms of stress wave mitigation, peak load retardation, and energy transfer. Simulations are conducted using various boundary conditions and domain sizes to evaluate their influence on the performance. Potential graded slab designs and material damping effects are also discussed and are both shown to reduce the energy transmitted from the impact surface to the opposing surface significantly. The results showed that the MM slabs had superior performance in reducing the peak stress wave and reducing the transfer of energy. This study demonstrates that resonant ceramic MMs are a promising material design with unique and tunable properties that can be used for stress wave mitigation and structural protection applications.

💡 Research Summary

This paper presents a comprehensive time‑domain finite‑element investigation of a low‑frequency locally resonant elastic metamaterial (MM) slab and compares its dynamic response to that of conventional monolithic slabs and several alternative microstructured designs. The authors focus on an H‑shaped ceramic unit cell (10 mm × 10 mm) made of alumina (ρ = 3985 kg m⁻³, E = 300 GPa, ν = 0.27) whose resonator is connected to the surrounding frame by narrow bridges, creating a resonant frequency well below 5 kHz and a band‑gap in the 20–40 kHz range. Three reference cells—monolithic, equal‑gap, and square‑cavity—share the same outer dimensions and mass, allowing a fair comparison of geometry‑induced effects.

Numerical simulations are performed in LS‑DYNA with a uniform 0.025 mm quadrilateral mesh. Impact loading is modeled by a cold‑rolled steel projectile (ρ = 7980 kg m⁻³, E = 200 GPa, ν = 0.3) striking the slab at 30 m s⁻¹, which generates an initial normal traction of about 600 MPa. Contact is handled by a mortar algorithm; no internal contacts within the unit cells are observed. The authors extract spatial averages of stress σ and particle velocity v at each interface between cells, compute the traction T = σ·n, and define the instantaneous energy flux P = T·v. Integrating P over time and interface area yields the transmitted energy E. Two performance metrics are introduced: (i) the energy‑transfer ratio, defined as the ratio of energy at the last interface to that at the first, and (ii) the “half‑time,” the time required for the transmitted energy to reach half of its peak value, which quantifies wave slowdown.

A parametric study on slab length is carried out using slabs composed of 4, 8, 16, and 32 unit cells (total thicknesses of 40 mm to 320 mm). Results show a clear trend: increasing the number of cells reduces the energy‑transfer ratio and increases the half‑time, with convergence observed beyond 16 cells, indicating an optimal slab thickness for the chosen geometry.

Boundary conditions are examined by applying prescribed traction (PT), prescribed velocity (PV), and a realistic projectile‑impact (PI) scenario. Under PT, the low mechanical impedance of the MM slab actually leads to higher energy transmission compared with the monolithic slab, whereas PV—representative of a free‑surface condition—produces a substantial reduction in transmitted energy. Consequently, the PV case is adopted for the remainder of the study as the most representative of an impact event.

Two graded‑slab configurations are designed by linearly varying the geometric parameters of the H‑cells through the thickness. Compared with a uniform‑property MM slab of identical mass and thickness, the graded slabs achieve an additional 10–15 % reduction in transmitted energy and exhibit a gradual decrease in wave speed, confirming the “wave‑slow‑down” effect predicted by effective‑medium theory.

Finally, material damping is introduced by assigning a loss factor of 0.02 to the alumina. The damped MM slab shows a further 30 % drop in the energy‑transfer ratio and a 1.5‑fold increase in half‑time relative to the undamped case, demonstrating a synergistic benefit of combining intrinsic material loss with local resonance.

Overall, the study demonstrates that locally resonant H‑shaped ceramic metamaterial slabs outperform conventional monolithic slabs in three key aspects: (1) peak stress reduction (up to ~40 % lower than monolithic), (2) transmitted energy mitigation (up to ~35 % lower), and (3) wave attenuation speed (half‑time up to twice as long). The authors also show that performance can be further enhanced through graded designs and modest material damping. Because the H‑cell geometry is compatible with additive manufacturing of ceramics, the proposed MM slabs are viable for practical applications such as blast and impact protection, structural health monitoring, and high‑performance vibration isolation.