Logo
ScienceStack
Table of Contents
Fetching ...
Sign In

Phase-Space Analysis of Elastic Vector Solitons in Flexible Mechanical Metamaterials

M. H. Duong, M. J. Reynolds

TL;DR

This work derives a revised continuum model for a nonlinear, flexible mechanical metamaterial built from a lattice of rotating squares, by applying a Galilean transformation to a discretely coupled displacement-rotation system. The authors obtain a coupled pair of PDEs for displacement and rotation, and analyze traveling-wave solitons in phase space by reducing the problem to a nonlinear ODE with an integrating-factor approach, revealing conditions for homoclinic orbits. The soliton solutions are compared to Deng et al., demonstrating a $1/9$ improvement in velocity-amplitude agreement with numerical simulations and establishing a close match to full numerics when using the sech$^2$ ansatz. The study highlights the importance of higher-order terms in the continuum limit and points to potential fully third-order regimes (e.g., cnoidal waves) for future work. Overall, this phase-space analysis provides a rigorous framework for understanding vector solitons in elastic metamaterials and improves predictive accuracy over prior continuum models.

Abstract

The purpose of this paper is to propose a revised continuum model from the discrete system introduced in [Deng et.al., PRL, 2017] . Using a Galilean transformation, we obtain an equation governing the soliton solutions in the phase plane - a second-order nonlinear ODE related to the Klein-Gordon equation with quadratic nonlinearity. These admit the well-known $\mathrm{sech}^2$ solutions, which we employ as an ansätz following [Deng et.al., PRL, 2017]. The resulting analysis yields soliton amplitudes and velocities that agree closely with numerical simulations, achieving an improvement of exactly 1/9 relative to the benchmark reported by the Harvard group.

Phase-Space Analysis of Elastic Vector Solitons in Flexible Mechanical Metamaterials

TL;DR

This work derives a revised continuum model for a nonlinear, flexible mechanical metamaterial built from a lattice of rotating squares, by applying a Galilean transformation to a discretely coupled displacement-rotation system. The authors obtain a coupled pair of PDEs for displacement and rotation, and analyze traveling-wave solitons in phase space by reducing the problem to a nonlinear ODE with an integrating-factor approach, revealing conditions for homoclinic orbits. The soliton solutions are compared to Deng et al., demonstrating a improvement in velocity-amplitude agreement with numerical simulations and establishing a close match to full numerics when using the sech ansatz. The study highlights the importance of higher-order terms in the continuum limit and points to potential fully third-order regimes (e.g., cnoidal waves) for future work. Overall, this phase-space analysis provides a rigorous framework for understanding vector solitons in elastic metamaterials and improves predictive accuracy over prior continuum models.

Abstract

The purpose of this paper is to propose a revised continuum model from the discrete system introduced in [Deng et.al., PRL, 2017] . Using a Galilean transformation, we obtain an equation governing the soliton solutions in the phase plane - a second-order nonlinear ODE related to the Klein-Gordon equation with quadratic nonlinearity. These admit the well-known solutions, which we employ as an ansätz following [Deng et.al., PRL, 2017]. The resulting analysis yields soliton amplitudes and velocities that agree closely with numerical simulations, achieving an improvement of exactly 1/9 relative to the benchmark reported by the Harvard group.
Paper Structure (10 sections, 4 theorems, 103 equations, 16 figures)

This paper contains 10 sections, 4 theorems, 103 equations, 16 figures.

Table of Contents

  1. Introduction
  2. Discrete model
  3. Derivation of the continuum limit
  4. Continuum equation for the displacement
  5. Continuum equation for the rotation
  6. Soliton solutions
  7. Phase space dynamics
  8. Numerical investigations
  9. A direct comparison with deng2017elastic
  10. Appendix: detailed derivation of the discrete model

Key Result

Lemma 1

For $\sigma\ne -3,-2,-1$ we have

Figures (16)

  • Figure 1: Full numerics
  • Figure 2: $c^2=0.625$
  • Figure 3: $c^2=0.630$
  • Figure 4: $c^2=0.635$
  • Figure 5: $c^2=0.640$
  • ...and 11 more figures

Theorems & Definitions (7)

  • Lemma 1
  • proof
  • Theorem 1: Reduction to a First-Order ODE using the Integrating Factor method
  • proof
  • Theorem 2: Equilibrium Points and Stability
  • proof
  • Theorem 3: Phase Portrait of the Reduced System