Normal form computation of nonlinear dispersion relationship for locally resonant metamaterial

Abstract

This article is devoted to the application of the parametrisation method for invariant manifold with a complex normal form style (CNF), for the derivation of high-order approximations of underdamped nonlinear dispersion relationships for periodic structures, more specifically by considering the case of a locally resonant metamaterial chain incorporating damping and various nonlinear stiffnesses. Two different strategies are proposed to solve the problem. In the first one, Bloch's assumption is first applied to the equations of motion, and then the nonlinear change of coordinates provided by the complex normal form style in the parametrisation method is applied. This direct procedure, which applies first the wave dependency to the original physical coordinates of the problem, is referred to as CNF-BP (for CNF applied with Bloch's assumption on physical coordinates). In the second strategy, the nonlinear change of coordinates provided by the parametrisation method, which relates the physical coordinates to the so-called normal coordinates, is first applied. Then the periodic assumption is used, thus imposing a Bloch wave ansatz on the normal coordinates. This method will be referred to as CNF-PN (for CNF with a periodic assumption on normal coordinates). In the conservative case, the CNF-PN strategy exhibits superior capability in capturing complex wave propagation phenomena, whereas the CNF-BP strategy encounters limitations in handling non-fundamental harmonics and the nonlinear interactions between host oscillators. For underdamped systems, the CNF-PN is rigorously validated and systematically compared against numerical techniques, a classical analytical perturbation technique (the method of multiple scales), and direct numerical time integration of annular chain structures.

Figures (9)

• Figure 1: Schematic diagram of the infinite locally resonant metamaterial chain with damping and various nonlinear stiffnesses.
• Figure 2: Nonlinear dispersion characteristics of a periodic chain with quadratic nonlinearity, $\beta_2$, in attachments: (a) Nonlinear dispersion manifold in the space of $k-\omega_{nl}-\log_{10}{u_n}$ calculated via EPMC-HBM; (b)-(d) Comparison of CNF-BP, CNF-PN, and EPMC-HBM in predicting nonlinear dispersion spectrum at different wave amplitudes, $u_n = 5\times10^{-3},1.5\times10^{-2}$, and $10^{-1}$. The red, blue, and cyan planes in (a) represent the cross-section of $u_n = 5\times10^{-3},1.5\times10^{-2}$, and $10^{-1}$. The acoustic branch in (d) is omitted. The used parameters are: $c_1 = c_2= \beta_1 =\gamma_1=\gamma_2=0$ and $\beta_2 = 2$.
• Figure 3: Backbone curves and invariant manifolds of (a)-(b) the acoustic branch (AB) at $k = 0.8\pi$ and (c)-(d) the optic branch (OB) at $k=0.3\pi$. The used parameters are: $c_1 = c_2= \beta_1 =\gamma_1=\gamma_2=0$ and $\beta_2 = 2$.
• Figure 4: Nonlinear dispersion characteristics of a periodic chain with cubic nonlinearity, $\gamma_1$, between host lattices: (a) Nonlinear dispersion manifold in the space of $k-\omega_{nl}-\log_{10}{u_n}$ calculated via EPMC-HBM; (b) Comparison of CNF-BP, CNF-PN, and EPMC-HBM in predicting nonlinear dispersion spectrum at wave amplitude, $\max{(u_n)} = 10^{-0.75}$. The red plane in (a) represents the cross-section of $u_n = 5\times10^{-0.75}$. The used parameters are: $c_1 = c_2= \beta_1 =\beta_2=\gamma_2=0$ and $\gamma_1 = 3$.
• Figure 5: Backbone curves and invariant manifolds of (a)-(b) the acoustic branch (AB) at $k = 0.4\pi$ and (c)-(d) optic branch (OB) at $k=0.6\pi$. The used parameters are: $c_1 = c_2= \beta_1 =\beta_2=\gamma_2=0$ and $\gamma_1 = 3$.