Arxiv:2104.14394V1 [Physics.Class-Ph] 29 Apr 2021 Metry Is Generally Undesired and Diﬃcult to Grasp, in Sition Between Reﬂection and Trapping Mechanisms [57]

Home , Seismic metamaterial, Wave Motion (journal)

Selective mode conversion and rainbow trapping via graded elastic waveguides

Jacopo M. De Ponti1, Luca Iorio1, Emanuele Riva2, Raﬀaele Ardito1, Francesco Braghin2, Alberto Corigliano1 1 Dept. of Civil and Environmental Engineering, Politecnico di Milano, Piazza Leonardo da Vinci, 32, 20133 Milano, Italy 2 Dept. of Mechanical Engineering, Politecnico di Milano, Via Giuseppe La Masa, 1, 20156 Milano, Italy

We experimentally achieve wave mode conversion and rainbow trapping in an elastic waveguide loaded with an array of resonators. Rainbow trapping is a phenomenon that induces wave confinement as a result of a spatial variation of the wave velocity, here promoted by gently varying the length of consecutive resonators. By breaking the geometrical symmetry of the waveguide, we combine the wave speed reduction with a reflection mechanism that mode-converts flexural waves impinging on the array into torsional waves travelling along opposite directions. The framework presented herein may open new opportunities in the context of wave manipulation through the realization of novel structural components with concurrent wave conversion and energy trapping capabilities.

I. INTRODUCTION ial and flexural vibrations [37, 38]. As reported by the aforementioned examples, the additional complexity induced by the symmetry break often reflects in relevant The emergence of new mechanisms to manipulate behaviours that are relatively unexplored in the field of the propagation of waves has attracted growing interest mechanics. Other approaches to induce modal coupling across different realms of physics, with multiple realiza- leverage nonlinear dynamics which, however, can be of- tions in quantum [1], electromagnetic [2, 3], acoustic [4–6] ten unpractical due to the large amplitudes required to and elastic systems [7–17]. In the context of mechanics, activate sufficiently strong nonlinear interactions [39, 40]. a number of works have recently investigated the emer- A recent line of work employs a graded array of res- gence of non-trivial topological phases in elastic struc- onators embedded in a host structure to manipulate tures, in analogy with relevant behaviours previously ob- wave propagation by taking advantage of the resonator- served in quantum physics [18]. Notable examples exten- structure interaction. This modulation strategy pro- sively explored in mechanics include defect-immune and motes a wavenumber transformation that, in turn, ac- scattering-free waveguides, which have been conceived tivates a spatial decrease of the wave velocity at dis- in analogy to the quantum Hall (QH) [19, 20], quan- tinct frequencies, which is blueprint of the so called rain- tum spin Hall (QSH) [21–24], and quantum valley Hall bow effect [41–45]. Originally discovered in electromag- (QVH) effects [25–28]. Other approaches to topology- netic systems in non-uniform, linearly tapered, planar based design leverage 1D structures augmented by a vir- waveguides with cores of negative index material [46], tual dimension in parameters space, to access topologi- the rainbow effect has been pursued within different re- cal properties typically attributed to 2D systems to pur- search fields and through numerous physical platforms, sue pumping of elastic waves [29–34] and nonreciprocity among which acoustic systems [47–49], water waves [50], [35, 36], to name a few. In other words, elastic waveg- and fluid loaded elastic plates [51]. Similar configura- uides are excellent candidates to explore new physics, es- tions have been combined with deep elastic substrates pecially due to the abundance of supported modes with to induce conversion between Rayleigh (R) into Shear distinct polarizations and coupling among them, which (S) or Pressure (P ) waves [53–56]. Despite the under- can be relatively simple to be established as compared to lying physics is driven by the spatial variation of the the electromagnetic counterparts. While the existence of wave speed, the existence of distinct wave modes and multiple modes that can hybridize in presence of asym- the nature of the coupling among them delineate a tran- arXiv:2104.14394v1 [physics.class-ph] 29 Apr 2021 metry is generally undesired and difficult to grasp, in sition between reflection and trapping mechanisms [57]. many cases the break of symmetry creates new opportu- The former is induced by resonances and leads to wave nities by lifting the accidental degeneracies. Mechanical scattering at the boundary of the first Brillouin zone. systems with broken symmetries along the thickness di- The latter is instead promoted by the coupling between rection can be employed for the nucleation of a double crossing wave modes, typical of the second Brillouin zone, Dirac cone dispersion that features coupling between oth- and hereafter employed to achieve concurrent wave mode- erwise degenerate states, which is the key to emulate the conversion and trapping. While the majority of the stud- spin orbital coupling in QSH-based waveguides [24]. Un- ies are focused on a single phenomenon, the coexistence dulated structures have been shown to exhibit frequency of multiple functionalities is not common, especially in gaps and wave directionality due to coupling between ax- 2 the field of mechanics. Motivated by this and by the opportunity to simulta- neously achieve trapping and mode-conversion, in this paper we consider an elastic waveguide loaded by a graded array of resonators with variable length. This gradual variation promotes rainbow trapping of flexural waves throughout the array. The considered implementation embodies a broken symmetry with respect to the shear center of the waveguide which, in turn, activates a mode-conversion between impinging flexural waves and torsional waves traveling along the opposite direction through a phenomenon that is known as mode locking in mechanics [58]. The numerical and experimental results presented herein demonstrate how a simple structure can be used to confine and mode convert elastic FIG. 1: (a) Graphical representation of the waveguide, waves, expanding the range of possibilities in the con- which is made of a homogeneous beam rigidly text of wave manipulation and control, with implications connected to a graded array of resonators that activates of technological relevance for applications involving me- the rainbow effect. Thanks to the broken geometrical chanical vibrations, such as nondestructive evaluation, symmetry, the waveguide operates as a converter ultrasonic imaging, and energy harvesting. + between input A0 and output T0− wave modes. (b) The paper is organized as follows. In section II, a simpli- Idealization of the conversion mechanism: the system fied model for the analysis of symmetry broken waveg- can be interpreted as two distinct elastic waveguides uides is presented and employed to tailor a suitable supporting A0 and T0 wave modes, which are coupled graded profile of locally unbalanced resonators. The the- using an array of resonators. (c) Schematic of the oretical aspects of the trapping and conversion mecha- reference cell, that consists in identical resonator pairs nisms are unfolded by way of a dispersion analysis, which separated by a non null distance ξ, along wih the is preparatory for numerical end experimental analysis employed reference system. pursued in section III. Concluding remarks are presented in section IV. in Figure 1(a) and further wave confinement due to multiple conversion between T A wave modes. 0 − 0 II. SIMULTANEOUS TRAPPING AND MODE For simplicity, we assume that the further wave modes CONVERSION: SIMPLIFIED MODELING AND supported by the waveguide are orthogonal to the excita- DESIGN tion mechanism and weakly coupled to the wave modes of relevance A0 and T0. To elucidate the design paradigm and conditions for As such, the functional design and conceptual analysis of selective trapping and reflection, we consider the waveg- the grading profile are carried out considering the elas- uide displayed in Figure 1. The implementation consists todynamic equations for an Euler-Bernoulli beam domi- in a slender beam, which is equipped with a set of at- nated by torsional-flexural coupling [60]: tachments of variable length linearly increasing along the main dimension of the waveguide, to activate the rainbow 2 + ∞ 2 effect for the asymmetric Lamb A0 (input) wave mode. EIyw,xxxx + ρAw,tt = mω0[Ψr (xr + na) + r=1 n= Each side of the beam accommodates an identical graded X X−∞ array of elements which, in turn, is rigidly connected to b + l w (x + na) + ( 1)r φ (x + na)]δ (x x na) the host structure. In addition to the graded profile, the − r − 2 r r − r − two arrays are spatially shifted by a distance ξ, to locally 2 + ∞ b + l break the geometrical symmetry with respect to the shear GJ φ ρJ w = ( 1)rmω2 [Ψ (x + na) + t ,xx − p ,tt − − 0 2 r r centre of the cross-section. Such a tailored symmetry- r=1 n= X X−∞ break drives a spatially growing coupling between trans- r b + l verse motion w and rotation φ of the waveguide, which w (xr + na) + ( 1) φr (xr + na)]δ (x xr na) + − − 2 − − are inherently linked to the asymmetric Lamb A0 and r b + l torsional T − wave modes. This coupling activates a con- 2 0 Ψr,tt (xr + na) + ω0[( 1) φ (xr + na) . version mechanism that can be interpreted in the light of − 2 modal interaction between distinct waveguides, accord- w (xr + na) + Ψr (xr + na)] = 0. ing to the schematic in Figure 1(b). The schematic il- − (1) lustrates how consecutive resonators locally trigger the where ( ) denotes the partial derivative ∂( )/∂x. w · ,x · transformation between A0 and T0 wave modes, induc- and φ are the vertical displacement and the rotation of th ing a reflected wave T0− exiting from the array, as shown the beam, Ψr is the displacement of the n resonator 3

(a) (b)

FIG. 2: (a) Dispersion surface ω (µ, lr) /2π. For ease of visualization, only the first three dispersion surfaces are displayed. The colour is proportional to the polarization factor p, which allows to discriminate between rotation φ and vertical displacement. Due to the broken symmetry and local resonance, the second and third dispersion surfaces are separated by a frequency gap. (b) First and second dispersion surfaces, mainly associated to flexural and torsional wave modes, respectively. The isofrequency line highlighted in black represents the operating frequency of the structure, while the coloured version is reported in the inset. (c) Group velocity cg upon varying the length of the resonators lr for the first two dispersion bands. In the figure, the group velocity line associated to the target dispersion branch is highlighted with darker dots. Such a target branch corresponds to the dispersion tangent to the isofrequency line displayed in (b). The group velocity cg (lr) computed for the target frequency is superimposed to the coloured plot with black dots. (d) The dispersion relation ω (µ) upon varying lr. The schematic illustrate the conversion mechanism that undergoes in correspondence of the target resonator. The system is excited + through the A0 wave mode in correspondence to the red star (i), which is wavenumber transformed to the cyan star (ii), due to the grading. As the wave reaches the zero group velocity region within the second Brillouin Zone, it is backscattered to the cyan star (iii) by a quantity equal to the reciprocal lattice vector µm, and wavenumber transformed to the back traveling T0− wave mode marked with the blue star (iv). Black dots represent the numerical solution computed in COMSOL Multiphysics. in the absolute reference frame. a = 40 mm is the the width and thickness of the host waveguide, respec- lattice size, EIy, ρA, GJt, and ρJp are the elastic and tively. For simplicity, the dynamic contribution of the inertial parameters associated to flexural and torsional resonator pair is approximated in terms of equivalent motion of the beam, consistently with the unit cell and bending stiffness k, participating mass m, and resonance reference frame displayed in Figure 1(c). E = 70 GPa frequency ω0 = k/m, which are dependent upon the and ρ = 2710 kg/m3 are the Young’s modulus and resonator geometry c, l, and h. The Dirac delta func- the material density; b = 7 mm and h = 2 mm are tion δ (x x nap) locally accounts for the presence of − r − 4 the resonators, placed at a distance xr = a/2 ξ/2 from gap, the dispersion relation exhibits wave modes charac- the left boundary of the lattice. Additional± details on terized by a balanced torsional and flexural motion, espe- the simplified model are reported in the supplementary cially in correspondence of the zero group velocity region material [59]. (a∂ω/∂µ = 0); this is of paramount importance for the We investigate the dispersion properties ω (µ, l) of the interplay between wave conversion and energy trapping waveguide, where l is considered as a free param- mechanisms, which will be discussed in the remainder of eter and µ = κxa is the normalized wavenumber. this section. To this end, we consider Ansatz for the displacement The functional design of the graded profile of resonators j(κ x ωt) w (x, t) =w ˆ (x) e− x − and for the rotation φ (x, t) = is accomplished following the general guidelines provided P ˆ j(κxx ωt) jnκmx in prior works [41, 43]. That is, rainbow trapping is here- φ (x) e− − wherew ˆ (x) = p= P wˆpe− , P − after pursued targeting a group velocity decrease along φˆ (x) = φˆ e jnκmx embody the x-periodicity of p= P p − P the beam at an operating frequency f = 2.12 kHz. To the medium and− κ = 2π/a is the modulation wavenum- P m that end, l is denoted as the relevant parameter linearly ber. As such, the transverse and torsional motions are varied along the main dimension of the beam. Such approximated in terms of p = P,..., +P plane wave − a variation yields a local wavenumber distribution µ (l) components. P = 3 is found to be sufficient for an accu- highlighted with black dots in Figure 2(b), where the dis- rate description of the dynamic behavior at the operating persion analysis is now limited to the first two branches. frequency region. Harmonic motion is also assumed for ˆ jωt The coloured version of the isofrequency line µ (l) is re- the resonators Ψr (xr) = Ψr (xr) e . The application of ported in the inset, illustrating that a variation of l the Plane Wave Expansion Method (PWEM), whose for- not only activates wave modes characterized by differ- mulation is detailed in the supplementary material [59], ent wavenumbers µ, but also promotes a transformation yields the following eigenvalue problem: between distinct polarizations (coloured dots). Consis- tently, the coloured group velocity profile cg (ω(µ), lr) = a∂ω/∂µ, is evaluated by finite difference for the entire K (κ , l ) ηˆ = ω2Mηˆ (2) x r wavenumber-parameter space and represented in Figure where K and M are the (2P + 2) (2P + 2) stiff- 2(c) for the first two dispersion branches. For ease of × T visualization, the second dimension l is eliminated and ˆ ˆ ˆ ness and mass matrices and ηˆ = wˆ, φ, Ψ1, Ψ2 substituted with arrows, to better illustrate the group accommodates the vector coefficients forh the distincti velocity profile cg (ω(µ)) in response to a variation of wave modes and resonators pair. The solution of the l. Here, the target dispersion curve, i.e. the dispersion eigenvalue problem ω (κx, l) is displayed in Figure 2(a) branch that touches the isofrequency line with zero group for the first three dispersion bands. In the figure, velocity, is represented with darker (coloured) dots. As the nature of the motion is discriminated through a the length of the resonator l is modified in space, the dis- color scale proportional to the polarization factor p = persion properties at the operating frequency naturally a wˆ (x) 2/ a wˆ (x) 2 + b a φˆ (x) 2 , which can be follow the black curve µ (l) displayed in Figure 2(b). This | 0 | | 0 | | 0 | modification is accompanied by a change in the group interpreted as a measure of the coupling between waves R R R velocity c (l) highlighted with black dots in Figure 2(c) characterized by distinct polarizations. Thus, a transi- g and, in turn, terminate in the zero group velocity region tion from blue to red denotes a transformation from pure away from the edges of the Brillouin zone, which is the rotation to a transverse motion. Some considerations fol- key factor to achieve rainbow trapping [57]. The analy- low. (i) The lone beam features an accidental degener- sis shows that the wave speed transformation is accom- acy, which corresponds to crossing flexural and rotational panied with a change in the polarization from a purely dispersion curves when l 0. As the length of the res- flexural mode that terminates into a region in which the onator increases, the geometrical→ symmetry of the cross- wave is characterized by a mixed torsional-flexural mo- section is locally broken, and a frequency gap emerges tion and c 0. from the otherwise degenerate states; this phenomenon g ≈ is generally known as mode locking in mechanics, and is In addition to the rainbow effect, the tailored broken hereafter employed to tailor selective rainbow trapping symmetry of the array is responsible for a reflection and reflection. (ii) Interestingly, when the attachments mechanism that mode-converts the impinging flexural length is smaller than a threshold (i.e. l 12 mm), the waves into torsional waves traveling along the oppo- bandgap formation mechanism is dominated≈ by purely site direction. This mechanism is hereafter discussed geometrical reasons; in contrast, for sufficiently high l assuming that most of the conversion takes place in values, the resonance frequency of the attachment lies correspondence of the zero group velocity region, where in the neighborhood of the gap, and the nature of the the coupling is maximum. As such, among the possible coupling is driven by a combination of broken symme- dispersion branches, we focus on the (target) curve try and local resonance; practically, the combination of highlighted with darker (coloured) dots in Figures 2(c-d) the two mechanisms yields a flattening of the dispersion and we assume that a smooth variation of l is provided curves that is beneficial in terms of wave velocity de- to reach the target curve. Now, consider an input wave + crease and trapping. (iii) In the neighborhood of the A0 , represented in figure 2(d) with the red star (i). 5

(a)

(b) (c) (d)

(e) (f) (g)

+ FIG. 3: (a) Numerical displacement field measured at different time instants for an input wave mode A0 , along with the corresponding dispersion relations for a waveguide endowed with absorbing boundaries. The diagrams show that + the input wave is reflected and converted from A0 to T0− in time. (b-d) The spectrograms display the temporal evolution of the spectral content for (b) a waveguide with the graded array of resonator; (c) a waveguide without the initial part of the grading, i.e. with only resonators 7 9; (d) a homogeneous waveguide, whithout array of resonators. The excitation signal is shown alongside the− spectrogram. The input wave is confined, trapped and delayed only if the complete array of resonators is present. (e-g) Corresponding out-of-plane displacement field filtered to separate the contribution of the A0 wave mode.

Such a wave impinges on the graded array of resonators Due to the conservation of the crystal momentum, part and, due to the smooth variation of l in space, the wave of the traveling wave is back-scattered to µ = µi µm mode experiences a wavenumber transformation that and, once again, is wavenumber-transformed from− the drives the transfer of spectral content to the cyan star cyan star (iii) to the blue star (iv) as a result of the (ii), where the group velocity at the operating frequency linear decrease of l, which produces a back-propagating is nullified and the normalized wavenumber is µ = µi. torsional wave mode T0− exiting the graded array. Part 6 of the energy is instead back-scattered as a flexural wave tion. In addition, the numerical dispersion is reduced to A− with spectral content µ = µi and delayed in time, dˆ(κ , ω) by taking the Root Mean Square (RMS) value 0 − x according with the wave speed reduction described by along κy which yields the diagram displayed in the fig- the dispersion analysis. ure. As expected, the spectral content in the waveg- To conclude the design of the graded array, the numerical uide is initially located in the positive half of the recip- + dispersion relation is evaluated numerically via COM- rocal space and in correspondence of the A0 dispersion SOL multiphysics environment spanning the parameter branch. When the flexural wave reaches the array, the space l [0, 22] mm and reported with black circles in wave is trapped and converted into a component travel- ∈ Figure 2(d) for the target dispersion curve. According to ing along the opposite direction, corresponding to the T0− the analysis, the necessary array of resonators to reach and the A0− dispersion branches exiting the graded array. near-zero group velocity is characterized by initial and The concept is elucidated in the spectrogram displayed final (target) length of l1 = 8 mm and l7 = 18.3 mm, in Figure 3(b), which is obtained by smoothly varying respectively, which are distributed along L = 7a mm for the position of the Gaussian function spanning the range a number of N = 7 cells. The array is then continued t [0, 20] ms. Also the dependence from frequency is 0 ∈ to a final length of l9 = 21.2 mm to prevent wave eliminated by taking the RMS value, which results in the propagation through the array for frequencies in the amplitude dˆ(κx, t0). The same analysis is performed on neighborhood of the target mode. a beam without the initial part of the array (see Figure 3(c)), responsible for the wave speed reduction, and on a homogeneous beam without graded array of resonators (reported in Figure 3(d)) equipped with clamps in corre- III. THEORY MEETS EXPERIMENTS spondence of resonators N = 8, 9.

The discussion is now focused on the transient analysis of wave propagation simulated in Abaqus environment, to corroborate the theoretical claims. The implementation consists in a homogeneous beam that serves as input domain for right traveling waves, which is fol- lowed by the graded array of resonators. Undesired reflections are avoided by way of absorbing boundaries applied to the left and right ends of the beam [61]. A transverse force with central frequency f0 = 2.12 kHz, width ∆f = 0.14 kHz and time duration of 15 ms is employed to provide excitation and confine the energy content in the neighborhood of the operating frequency. Additional details on the numerical methods are reported in the supplementary material [59]. Numerical results for consecutive time instants t0 = FIG. 4: Experimental setup employed to measure wave 2, 6, 10, 12 ms are displayed in Figure 3(a). As expected, propagation. The waveguide is excited using an at the beginning of the time simulation, the motion is electrodynamic shaker (i), while the opposite side is dominated by flexural waves. As time elapses, part of the clamped (ii). Wave propagation is measured on the energy is transferred to torsional wave modes, starting surface of the waveguide through a 3D Polytec Scanner from the snapshot captured in Figure 3(a)II at t0 = 6 ms, Laser Doppler Vibrometer (iii), which is able to where the homogeneous trait is dominated by transverse separate the 3D velocity field. The top view of the waves and the graded array exhibits rotation, due to the graded array is also provided (iv). broken geometrical symmetry. In contrast, the displacement field displayed in Figures 3(a)III-IV is characterized A qualitative comparison of Figures 3(b-d) reveals by mixed wave motion within the entire spatial domain. that: (i) the presence of the graded array shifts the These considerations are further confirmed by the nu- backward propagating wave A0− by 5 ms, as a result of ˆ 2 2 2 merical dispersion d (κx, κy, ω) = √wˆ +u ˆ +v ˆ , which the tailored trapping mechanism, while the backward is illustrated alongside the displacement fields in Figure propagating A0− component is not delayed in a relevant 3(a), wherew ˆ =w ˆ (κx, κy, ω),v ˆ =w ˆ (κx, κy, ω), and way in case the initial part of the array is not considered; uˆ =w ˆ (κx, κy, ω) are the 3D Fourier-transformed dis- (ii) the wave speed reduction is accompanied by a + placement fields of the homogeneous part of the waveg- transfer of energy between A0 and T0− wave modes, uide. To evaluate the numerical dispersion, the time his- which occurs approximately 3 ms before the A0− wave tories w (x, y, t), v (x, y, t), and u (x, y, t) are windowed in mode back-propagates; (iii) the amount of conversion the neighborhood of the probed time instants t0 by way is dramatically reduced in case the first part of the (t t )2/2c2 of a Gaussian function g(t) = e− − 0 where c is a array is not present. To further confirm that the graded parameter that controls the shape of the Gaussian func- array operates as delay line and mode converter for 7

(a)

(b) (c)

FIG. 5: (a) Experimental velocity field and corresponding dispersion relation measured at different time instants. (b) Experimental and (c) numerical spectrogram, illustrating the evolution of the spectral content in time. Experimental and numerical results are characterized by a comparable temporal delay for the A0− wave mode, while different conversion amplitudes are observed, due to the presence of imperfections and misalignments that enhance the torsional coupling.

flexural/torsional waves, the displacement field w (x, y, t) to avoid excessive geometrical deformation. The wave- is filtered in wavenumber domain for κx values outside field is measured on the surface of the system through a 40 rad/m and 100 rad/m. The filtered displacement field Polytec 3D Scanner Laser Doppler Vibrometer (SLDV), allows to isolate and graphically visualize the A0 mode which is able to separate the 3D velocity field in both evolution in time, that is displayed in Figure 3(e-g) for space and time. The excitation signal is synchronously the three reference configurations. As expected, only if started with the acquisition which, in turn, is averaged in the graded array is present, the numerical wavefield is time to decrease the noise. Additional information on the altered in a way that the energy slowly vanishes and experimental methods are reported in the supplementary is mode-converted to the counter-propagating T0− wave material [59]. The experimental results in terms of the mode. The non-converted energy is back-reflected as 3D displacement field w (x, y, t), v (x, y, t), and u (x, y, t), a A0− mode and shifted in time, consistently with the along with the corresponding dispersion dˆ(κx, ω) cen- above discussion. tered at different time instants t0 are displayed in Figure 5(a), which are similar to the ones displayed in Figure + 3(a). After the input A0 mode is injected in the struc- The numerical results are validated through an exper- ture, part of the energy is mode-converted in a torsional imental analysis in the transient regime, which is per- T0± mode and part is delayed in time. In contrast to the formed on a structure identical to the one employed in results displayed in Figure 3(a), the propagating waves the simulations, except for the boundary conditions that cannot escape the system, and the amount of dissipa- are hereafter discussed. The schematic of the experimen- tion is negligible. We also observe that, in comparison tal setup is shown in Figure 4. The system is rigidly con- to the numerical case-study, the mode-conversion takes nected to a LDS v406 electrodynamic shaker, to provide place in shorter time, and the associated amplitude is excitation, while the opposite end is clamped to ground, 8 greater. This is attributed to: (i) the absence of ab- waves into torsional waves traveling along the opposite sorbing boundaries, which facilitates the accumulation of direction. Such a conversion is accompanied with trap- energy, especially that transferred to the torsional wave ping of flexural waves that are delayed while propagat- mode characterized by faster propagation speed and trig- ing through the array. The complex interplay between gers the formation of a standing mode; (ii) the presence trapping and conversion is explained in terms of disper- of misalignment and imperfections that facilitate the con- sion analysis, where the dynamic behavior is dominated version between wave modes that populate the dispersion by a local symmetry break and the smooth variation of at the operating frequency. the resonator length in space, which locally activates the The same considerations apply for the experimental and energy transfer between distinct wave modes and glob- numerical spectrograms in Figure 5(b,c). The figures il- ally determines a wave speed reduction. The concept + lustrate a similar amount of delay between incident A0 can be generalized to multi-dimensional or multi-physics and reflected A0−,T0− wave modes, whereas in the exper- systems and may open new opportunities in the context imental diagram in Figure 5(b) we observe an extra con- of wave manipulation and control in metamaterial struc- version amplitude between A0 and T0 waves that is trig- tures with concurrent wave conversion and trapping ca- gered by the undesired curvature of the beam, due to the pabilities. manufacturing process that is not present in the numerical analysis. Also, consistently with the simulation performed with absorbing boundaries, (3(b)) the converted A0− wave mode is delayed by the same amount as for the experimental data displayed in 5(b), which further con- firms that the rainbow trapping occurs as expected. We + finally notice that an additional T0 component is mea- ACKNOWLEDGEMENTS sured; this is attributed to the finite length of the system, which allows for spurious edge reflections of the T0− waves exiting the array. The support of the H2020 FET-proactive project MetaVEH under grant agreement No. 952039 is acknowl- edged. We also gratefully acknowledge the Italian Min- CONCLUSIONS istry of Education, University and Research for the support provided through the Project “Department of Ex- In this manuscript we have experimentally realized a cellence LIS4.0—Lightweight and Smart Structures for graded structure that is capable of converting flexural Industry 4.0.”

[1] Hasan, M. Zahid, and Charles L. Kane. Colloquium: look.” Applied Mechanics Reviews 66.4 (2014). topological insulators. Reviews of modern physics 82.4 [10] C . Sugino, A. Erturk, Analysis of multifunctional piezo- (2010): 3045. electric metastructures for low-frequency bandgap forma- [2] Lu, Ling, John D. Joannopoulos, and Marin Soljaˇcić. tion and energy harvesting, J. Phys. D: Appl. Phys. 51 ”Topological photonics.” Nature photonics 8.11 (2014): 215103 (2018). 821-829. [11] E. G. Williams, P. Roux, M. Rupin, W. A. Kuperman, [3] Khanikaev, Alexander B., et al. ”Photonic topological Theory of multiresonant metamaterials for A0 Lamb insulators.” Nature materials 12.3 (2013): 233-239. waves, Phys. Rev. B 91, 104307 (2015). [4] Cummer, Steven A., Johan Christensen, and Andrea Alù. [12] C. Sugino, S. Leadenham, M. Ruzzene, A. Erturk, On ”Controlling sound with acoustic metamaterials.” Nature the mechanism of bandgap formation in locally resonant Reviews Materials 1.3 (2016): 1-13. finite elastic metamaterials, Journal of Applied Physics [5] Ma, Guancong, and Ping Sheng. ”Acoustic metamateri- 120, 134501 (2016). als: From local resonances to broad horizons.” Science [13] S. D. M. Adams, R. V. Craster, S. Guenneau, Guided and advances 2.2 (2016): e1501595. standing Bloch waves in periodic elastic strips, Waves in [6] Yang, Zhaoju, et al. ”Topological acoustics.” Physical Random and Complex Media 19 (2) (2009), 321—346. review letters 114.11 (2015): 114301. [14] V. Pagneux, A. Maurel, Lamb wave propagation in in- [7] Nassar, Hussein, et al. ”Nonreciprocity in acoustic and homogeneous elastic waveguides, Proc. R. Soc. A 462 elastic materials.” Nature Reviews Materials 5.9 (2020): (2006), 1315–1339 667-685. [15] A. Khelif, B. Djafari-Rouhani, J. O. Vasseur, [8] Wang, Yan-Feng, et al. ”Tunable and active phononic P. A. Deymier, Transmission and dispersion relations of crystals and metamaterials.” Applied Mechanics Reviews perfect and defect-containing waveguide structures in 72.4 (2020). phononic band gap materials, Phys. Rev. B 68, 024302 [9] Hussein, Mahmoud I., Michael J. Leamy, and Massimo (2003). Ruzzene. ”Dynamics of phononic materials and struc- [16] Y. Pennec, B. D. Rouhani, C. Li, J. M. Escalante, tures: Historical origins, recent progress, and future out- A. Martinez, S. Benchabane, V. Laude, N. Papanikolaou, 9

Band gaps and cavity modes in dual phononic and pho- servation of nonreciprocal band gaps in a space-time- tonic strip waveguides, AIP Advances 1, 041901 (2011). modulated beam using a shunted piezoelectric array. [17] R. V. Craster, S. Guenneau, S. D. M. Adams, Mecha- Physical Review Applied, 13(3), 031001. nism for slow waves near cutoff frequencies in periodic [36] Attarzadeh, M. A., Callanan, J., Nouh, M. (2020). waveguides, Phys. Rev. B 79, 045129 (2009). Experimental observation of nonreciprocal waves in a [18] Huber, Sebastian D. ”Topological mechanics.” Nature resonant metamaterial beam. Physical Review Applied, Physics 12.7 (2016): 621-623. 13(2), 021001. [19] Wang, Pai, Ling Lu, and Katia Bertoldi. ”Topologi- [37] G. Trainiti, J. J. Rimoli, M. Ruzzene, Wave propagation cal phononic crystals with one-way elastic edge waves.” in periodically undulated beams and plates, International Physical review letters 115.10 (2015): 104302. Journal of Solids and Structures 75-76 (2015) 260-276. [20] Chen, H., et al. ”Mechanical quantum hall effect in time- [38] G. Trainiti, J. J. Rimoli, M. Ruzzene, Wave propagation modulated elastic materials.” Physical Review Applied in undulated structural lattices, International Journal of 11.4 (2019): 044029. Solids and Structures 97–98, (2016) 431-444. [21] Süsstrunk,Roman, and Sebastian D. Huber. ”Observa- [39] W. Jiao, S. Gonella, Intermodal and Subwavelength En- tion of phononic helical edge states in a mechanical topo- ergy Trapping in Nonlinear Metamaterial Waveguides, logical insulator.” Science 349.6243 (2015): 47-50. Phys. Rev. Applied 10, 024006 (2018). [22] Chen, Hui, et al. ”Elastic quantum spin Hall effect in [40] R. Ganesh, S. Gonella, From Modal Mixing to Tun- kagome lattices.” Physical Review B 98.9 (2018): 094302. able Functional Switches in Nonlinear Phononic Crystals, [23] Chaunsali, Rajesh, Chun-Wei Chen, and Jinkyu Yang. Phys. Rev. Lett. 114, 054302 (2015). ”Subwavelength and directional control of flexural waves [41] J. M. D. Ponti, A. Colombi, R. Ardito, F. Braghin, in zone-folding induced topological plates.” Physical Re- A. Corigliano, R. V. Craster, Graded elastic metasurface view B 97.5 (2018): 054307. for enhanced energy harvesting, New J. Phys. 22 (2019) [24] Miniaci, M., Pal, R. K., Morvan, B., Ruzzene, M. (2018). 013013. Experimental observation of topologically protected he- [42] J. M. De Ponti, A. Colombi, E. Riva, R. Ardito, lical edge modes in patterned elastic plates. Physical Re- F. Braghin, A. Corigliano, R. V. Craster, Experimen- view X, 8(3), 031074. tal investigation of amplification, via a mechanical delay- [25] Riva, E., Quadrelli, D. E., Cazzulani, G., Braghin, F. line, in a rainbow-based metamaterial for energy harvest- (2018). Tunable in-plane topologically protected edge ing, Appl. Phys. Lett. 117, 143902 (2020). waves in continuum Kagome lattices. Journal of Applied [43] J. M. De Ponti, Graded Elastic Metamaterials for Energy Physics, 124(16), 164903. Harvesting, Springer International Publishing (2021). [26] Vila, J., Pal, R. K., Ruzzene, M. (2017). Observation of [44] M. Alshaqaq, A. Erturk, Graded multifunctional piezo- topological valley modes in an elastic hexagonal lattice. electric metastructures for wideband vibration attenu- Physical Review B, 96(13), 134307. ation and energy harvesting, Smart Mater. Struct. 30 [27] Liu, T. W., Semperlotti, F. (2018). Tunable acoustic val- 015029 (2021). ley–hall edge states in reconfigurable phononic elastic [45] G. J. Chaplain, J. M. De Ponti, G. Aguzzi, A. Colombi, waveguides. Physical Review Applied, 9(1), 014001. R. V. Craster, Topological Rainbow Trapping for Elastic [28] Liu, T. W., Semperlotti, F. (2019). Experimental evi- Energy Harvesting in Graded Su-Schrieffer-Heeger Sys- dence of robust acoustic valley Hall edge states in a non- tems, Phys. Rev. Applied 14, 054035 (2020). resonant topological elastic waveguide. Physical Review [46] K. L. Tsakmakidis, A. D. Boardman, O. Hess, ‘Trapped Applied, 11(1), 014040. rainbow’ storage of light in metamaterials, Nature vol- [29] Rosa, M. I., Pal, R. K., Arruda, J. R., Ruzzene, M. ume 450, pages397–401(2007). (2019). Edge states and topological pumping in spa- [47] J. Zhu, Y. Chen, X. Zhu, F. J. Garcia-Vidal, X. Yin, tially modulated elastic lattices. Physical review letters, W. Zhang, and X. Zhang, “Acoustic rainbow trapping,” 123(3), 034301. Sci. Rep. 3, 1728 (2013). [30] Y. Xia, A. Erturk, M. Ruzzene, Topological Edge States [48] V. Romero-Garc´ıa,R. Picó,A. Cebrecos, V. J. Sánchez- in Quasiperiodic Locally Resonant Metastructures, Phys. Morcillo, and K. Staliunas, Enhancement of sound in Rev. Applied 13, 014023 (2020). chirped sonic crystals, Appl. Phys. Lett. 102, 091906 [31] Riva, E., Rosa, M. I., Ruzzene, M. (2020). Edge states (2013). and topological pumping in stiffness-modulated elastic [49] A. Cebrecos, R. Picó,V. J. Sánchez-Morcillo, K. Stal- plates. Physical Review B, 101(9), 094307. iunas, V. Romero-Garc´ıa,L. M. Garcia-Raffi, Enhance- [32] Grinberg, I. H., Lin, M., Harris, C., Benalcazar, W. A., ment of sound by soft reflections in exponentially chirped Peterson, C. W., Hughes, T. L., Bahl, G. (2020). Robust crystals, AIP Advances 4, 124402 (2014). temporal pumping in a magneto-mechanical topological [50] L. G. Bennetts, M. A. Peter, R. V. Craster, Graded res- insulator. Nature communications, 11(1), 1-9. onator arrays for spatial frequency separation and ampli- [33] Riva, E., Casieri, V., Resta, F., Braghin, F. (2020). fication of water waves, J. Fluid Mech. 854 R4 (2018). Adiabatic pumping via avoided crossings in stiffness- [51] E. A. Skelton, R. V. Craster, A. Colombi, D. J. Colquitt, modulated quasiperiodic beams. Physical Review B, The multi-physics metawedge: graded arrays on fluid- 102(1), 014305. loaded elastic plates and the mechanical analogues of [34] Cheng, W., Prodan, E., Prodan, C. (2020). Experi- rainbow trapping and mode conversion, New J. Phys. mental Demonstration of Dynamic Topological Pumping 20 053017 (2018). across Incommensurate Bilayered Acoustic Metamateri- [52] A. Deitmar, S. Echterhoff, Principles of Harmonic Anal- als. Physical Review Letters, 125(22), 224301. ysis, Springer, ISBN 978-3-319-05791-0 (2014). [35] Marconi, J., Riva, E., Di Ronco, M., Cazzulani, G., [53] A. Colombi, D. Colquitt, P. Roux, S. Guenneau, Braghin, F., Ruzzene, M. (2020). Experimental ob- R. V. Craster, A seismic metamaterial: The resonant 10

metawedge, Sci. Rep. 6 (1) (2016) 1–6. cations for energy harvesting, New J. Phys. 22 063024 [54] A. Colombi, V. Ageeva, R. J. Smith, A. Clare, R. Pa- (2020). tel, M. Clark, D. Colquitt, P. Roux, S. Guenneau, R. V. [58] B. R. Mace, E. Manconi, Wave motion and dispersion Craster, Enhanced sensing and conversion of ultrasonic phenomena: Veering, locking and strong coupling effects, Rayleigh waves by elastic metasurfaces, Sci. Rep. 7 (1) The Journal of the Acoustical Society of America 131, (2017) 6750. 1015 (2012). [55] D. Colquitt, A. Colombi, R. Craster, P. Roux, S. Guen- [59] See Supplemental Material at xxxx for more details neau, Seismic metasurfaces: Sub-wavelength resonators on the simulation procedures, experimental setup and and rayleigh wave interaction, Journal of the Mechanics methodology and Physics of Solids 99 (2017) 379–393. [60] J. Achenbach, Wave propagation in elastic solids, Vol. 16, [56] G. J. Chaplain, J. M. De Ponti, A. Colombi, R. Fuentes- Elsevier, 2012. Dominguez, P. Dryburg, D. Pieris, R. J. Smith, A. Clare, [61] P. Rajagopal, M. Drozdz, E. A. Skelton, M. J. Lowe, M. Clark, R. V. Craster, Tailored elastic surface to body R. V. Craster, On the use of absorbing layers to simulate wave Umklapp conversion, Nat. Commun. 11 (2020) 1–6. the propagation of elastic waves in unbounded isotropic [57] G. J. Chaplain, D. Pajer, J. M. De Ponti, R. V. Craster, media using commercially available finite element pack- Delineating rainbow reflection and trapping with appli- ages, NDT & E Int. 51 (2012) 30–40. Selective mode conversion and rainbow trapping via graded elastic waveguides: Supplementary Material

To aid insight into the underlying physics of selective mode conversion and rainbow trapping, we give a detailed description of the analytical procedures, numerical models, sample manufacturing and experimental methods.

I. SUPPLEMENTARY NOTE 1: NUMERICAL PROCEDURE FOR THE PLANE WAVE EXPANSION METHOD

Consider the elastodynamic equation for an Euler Bernoulli beam subjected to ﬂexural-torsional coupling:

2 + ∞ 2 EIyw,xxxx + ρAw,tt = mω0 Ψr (xr + na) + r=1 n= X X−∞ b + l w (x + na) + ( 1)r φ (x + na) δ (x x na) − r − 2 r r − r − 2 + ∞ b + l GJ φ ρJ w = ( 1)rmω2 Ψ (x + na) + (S1) t ,xx − p ,tt − − 0 2 r r r=1 n= X X−∞ b + l w (x + na) + ( 1)r φ (x + na) δ (x x na) − r − 2 r r − r − b + l Ψ (x + na) + ω2 ( 1)r φ (x + na) w (x + na) + Ψ (x + na) = 0 r,tt r 0 − 2 r − r r r

Where EIy, ρA are the flexural stiffness and the linear mass density of the beam, Iy is the second moment of area 3 Iy = bh /12. b = 7 mm and h = 2 mm are the width and the thickness of the beam. E = 70 GPa is the Young’s 3 Modulus, ρ = 2710 Kg/m is the material density and A = bh is the cross-section area. GJt and ρJp are the rotary stiffness and inertia. G = E/2(1 + ν) is the shear modulus (ν = 0.33), Jt is the torsional constant, and Jp is the polar 3 3 second moment of area. The resonator stiffness and mass are k = 3EIres/l and m = ρchl, where Ires = ch /12 is the second moment of area of the resonator. l, h = 2 mm and c = 5 mm are the length, thickness, and width of the resonator. ω0 = k/m is the bending resonance frequency of the attachments. In the attempt to evaluate the dispersion properties ω (κx, l) of the undelying medium, Ansatz are considered p j(κ x ωt) j(κ x ωt) for the displacement w (x, t) =w ˆ (x) e− x − and for the rotation φ (x, t) = φˆ (x) e− x − wherew ˆ (x) = P P jnκmx ˆ ˆ jnκmx p= P wˆpe− , φ (x) = p= P φpe− reflect the periodicity of the medium and κm = 2π/a is the mod- arXiv:2104.14394v1 [physics.class-ph] 29 Apr 2021 − − ˆ jωt ulationP wavenumber. HarmonicP motion is assumed for the resonator displacement Ψr (xr) = Ψr (xr) e . Due to the periodic nature of the underlying medium, the relation between displacements and rotations between consecutive units writes:

jκna jκna jκna w (xr + na) = w (xr) e− φ (xr + na) = φ (xr) e− Ψr (xr + na) = Ψr (xr) e− (S2)

plugging Equation S2 into Equation S1 gives: 2

2 + ∞ 2 EIyw,xxxx + ρAw,tt = mω0 Ψr (xr) + r=1 n= X X−∞ r b + l jκna w (x ) + ( 1) φ (x ) e− δ (x x na) − r − 2 r r − r − 2 + ∞ b + l GJ φ ρJ w = ( 1)rmω2 Ψ (x ) + (S3) t ,xx − p ,tt − − 0 2 r r r=1 n= X X−∞ r b + l jκna w (x ) + ( 1) φ (x ) e− δ (x x na) − r − 2 r r − r − b + l Ψ (x ) + ω2 ( 1)r φ (x ) w (x ) + Ψ (x ) = 0 r,tt r 0 − 2 r − r r r The right hand side of Equation S3 is non-null if na = x xr, which allows to eliminate the dependence on n in the exponential term. Equation S3 is conveniently rewritten− as:

2 + ∞ 2 EIyw,xxxx + ρAw,tt = mω0 Ψr (xr) + r=1 n= X X−∞ r b + l jκ(x x ) w (x ) + ( 1) φ (x ) e− − r δ (x x na) − r − 2 r r − r − 2 + ∞ b + l GJ φ ρJ w = ( 1)rmω2 Ψ (x ) + (S4) t ,xx − p ,tt − − 0 2 r r r=1 n= X X−∞ r b + l jκ(x x ) w (x ) + ( 1) φ (x ) e− − r δ (x x na) − r − 2 r r − r − b + l Ψ (x ) + ω2 ( 1)r φ (x ) w (x ) + Ψ (x ) = 0 r,tt r 0 − 2 r − r r r

jnκm(x xr ) The term n∞= δ (x xr na) = 1/a n∞= e − can be be written in terms of exponential functions, where 1/a are the−∞ Fourier− coeﬃcients− of the expansion.−∞ We are now ready to enforce the Bloch-wave solution for the displacementsP and rotations. We get to: P

2 2 + ∞ 4 j(κ+pκ )x 2 ∞ j(κ+pκ )x mω0 ∞ jnκ (x x ) EI (κ + pκ ) wˆ e− m ω ρA wˆ e− m = Ψ e m − r + y m p − p a r p= p= r=1 n= " X−∞ X−∞ X X−∞

∞ j(κ+pκ )x jnκ (x x ) r b + l ∞ j(κ+pκ )x jnκ (x x ) jκ(x x ) wˆ e− m r e m − r + ( 1) φˆ e− m r e m − r e− − r − p − 2 p p= p= # X−∞ X−∞ 2 2 + ∞ 2 j(κ+pκ )x 2 ∞ j(κ+pκ )x mω0 b + l r ∞ jnκ (x x ) GJ (κ + pκ ) φˆ e− m ρJ ω φˆ e− m = ( 1) Ψ e m − r + t m p − p p − a 2 − r p= p= r=1 n= " X−∞ X−∞ X X−∞

∞ j(κ+pκ )x jnκ (x x ) r b + l ∞ j(κ+pκ )x jnκ (x x ) jκ(x x ) wˆ e− m r e m − r + ( 1) φˆ e− m r e m − r e− − r − p − 2 p p= p= # X−∞ X−∞ b + l ∞ ∞ ω2Ψˆ + ω2 ( 1)r φˆ ej(κ+pκm)x wˆ ej(κ+pκm)x + Ψˆ = 0 − r 0 − 2 p − p r " p= p= # −∞ −∞ X X (S5) 3

for ease of notation, the resonator displacement at a distance xr is replaced with Ψr (xr) = Ψr. We now multiply by an orthogonal function ejsκmx and integrate over the unit cell spatial domain D [0, a]. The equations are simpliﬁed as: ∈

2 2 4 2 mω0 j(κ+sκ )x ∞ j[p s]κ x r b + l ∞ j[p s]κ x EI (κ + sκ ) wˆ ω ρAwˆ Ψ e m r wˆ e− − m r + ( 1) φˆ e− − m r = 0 y m s − s − a r − p − 2 p r=1" p= p= # X X−∞ X−∞

2 2 2 2 mω0 b + l r j(κ+sκ )x ∞ j[p s]κ x GJ (κ + sκ ) φˆ ρJ ω φˆ + ( 1) Ψ e m r wˆ e− − m r + t m s − p s a 2 − r − p r=1 " p= X X−∞

r b + l ∞ j[p s]κ x +( 1) φˆ e− − m r = 0 − 2 p p= # X−∞

b + l ∞ ∞ ω2Ψˆ + ω2 ( 1)r φˆ ej(κ+pκm)x wˆ ej(κ+pκm)x + Ψˆ = 0 − r 0 − 2 p − p r " p= p= # −∞ −∞ X X (S6) Equation S6 is truncated to the ﬁrst s [ P,P ], p [ P,P ] plane wave components and is written in a matrix form, making explicit the equations for the∈ consecutive− ∈ resonators:−

K wˆ ω2M wˆ + Λ Ψˆ + Λ Ψˆ + C wˆ + C φˆ = 0 1 − 1 11 1 12 2 11 12 2 K2φˆ ω M2φˆ + Λ21Ψˆ 1 + Λ22Ψˆ 2 + C21wˆ + C22φˆ = 0  − (S7)  ˆ 2 ˆ 2 ˆ Γ11wˆ + Γ12φ + ω0Ψ1 ω Ψ1 = 0  − 2 2 Γ21wˆ + Γ22φˆ + ω Ψˆ 2 ω Ψˆ 2 = 0  0 −    where the displacement and rotation vectors wˆ and φˆ accommodate the expansion coeﬃcients wˆ = [w P , . . . , wP ], − φˆ = [φ P , . . . , φP ]. The above matrices are deﬁned as: − 4

...... 4 2 K1 =  EIy (κ + sκm)  K2 =  GJt (κ + sκm)  . .  ..  ..         . . . . mω2 mω2 b + l 0 2 j[s p]κmx1 0 2 r j[s p]κmx1 C11 = ... r=1 e − ... C12 = ... r=1( 1) e − ... a . − a 2 − .  .   .   P   P      . . mω2 b + l 0 2 j[s p]κmx1 C21 = C12 C22 = ... r=1 e − ... a 2 .  .   P    (S8) . . 2 . 2 . mω0 j(κ+sκ )x mω0 j(κ+sκ )x Λ11 = e m 1  Λ12 = e m 2  − a . − a .  .   .          . . 2 . 2 . mω0 b + l j(κ+sκ )x mω0 b + l j(κ+sκ )x Λ21 = e m 1  Λ22 = e m 2  − a 2 . a 2 .  .   .          b + l Γ = ω2 ... e j(κ+sκm)x1 ... Γ = ω2 ... e j(κ+sκm)x1 ... 11 − 0 − 12 − 0 2 − b + l Γ = ω2 ... e j(κ+sκm)x2 ... Γ = ω2 ... e j(κ+sκm)x2 ... 21 − 0 − 22 0 2 − and M1 = ρAI and M2 = ρJpI are the mass and inertia matrices. I is the (2P + 1) (2P + 1) identity matrix. The × T vector coeﬃcients wˆ, φˆ and the resonator displacements Ψ1,2 are stored in ηˆ = wˆ, φˆ, Ψˆ 1, Ψˆ 2 for a more compact representation of the dynamic equations, yielding the dispersion relation: h i

2 K (κx, lr) ηˆ = ω Mηˆ (S9)

where K and M are the (2P + 2) (2P + 2) wavenumber dependent stiﬀness and mass matrices of the unit cell. ×

II. SUPPLEMENTARY NOTE 2: NUMERICAL METHODS

The numerical models employed for the transient analysis are based on the commercial software Abaqus and, specifically, a finite element discretisation is accomplished through full 3D stress quadratic elements (C3D20). The analysis is performed via implicit time integration based on the Hilber-Hughes-Taylor operator [1], an extension of the Newmark β-method with a constant time increment dt = 0.02 ms. The system is forced through an imposed pressure 2 p(t) applied to a rectangular surface S0 = 42.4 mm . The input signal is a finite burst p(t) = p0 sin (2πfc)w(n) with amplitude p0 = 1 MPa, where w(n) is a Hann window with time duration 15 ms. As a result, the input signal is charaterized by central frequency fc = 2.12 kHz and a spectral content of width ∆f = 0.14 kHz. In the numerical part of the paper, an infinite waveguide model is employed to avoid spurious edge reflections. This is achieved by imposing absorbing boundary conditions at both ends of the strip. The absorbing boundaries are implemented in Abaqus using the ALID (Absorbing Layers using Increasing Damping) method, as widely used in the 6 elastic wave community, adopting a cubic law function for mass proportional Rayleigh damping with cM max = 10 . A detailed description is provided in [2]. The absorbing boundary is discretized using 100 finite elements along x, covering a total length of approximately 3λmax, where λmax is the maximum wavelength travelling in the system. 5

III. SUPPLEMENTARY NOTE 3: SAMPLE MANUFACTURING AND EXPERIMENTAL METHODS

The system is manufactured from a 2 mm thick aluminium plate through a laser cutting technology. A Bystronic Fibra (3000 W) laser machine is used to cut the slab with a position tolerance of 0.1 mm. A top view of the system is shown in Fig. S1, which is made of an elastic waveguide 850 mm long, 7 mm± wide and 2 mm thick. The system includes an initial portion 350 mm long that serves as an input domain. The ﬁnal shape of the specimen includes the array of resonators, thereby avoiding additional operations such as soldering or mechanical connections that would have led to imperfections in the prototype. The array is composed of 9 unit cell of size a = 40 mm, each one containing two identical resonators shifted by a distance ξ = 20 mm. The geometrical dimensions of the resonators for each unit cell are stored in Tab. S1 and numbered from 1 to 9.

FIG. S1: Top view of the system, made of a homogeneous beam rigidly connected to a graded array of resonators.

CELL NUMBER: x y z { | } ~ li [mm]: 8.1 10.0 11.6 13.0 14.8 16.3 18.0 19.7 21.2 ci [mm]: 5 hi [mm]: 2

TABLE S1: Geometric parameters of the graded array of resonators. Each unit cell contains two identical resonators of length li, width ci, and thickness hi, shifted by a constant distance of ξ = 20 mm.

The experimental tests are performed according to the following setup. The waveguide is mechanically joined to a LDS v406 electrodynamic shaker at the left boundary, to provide an out-of-plane input excitation and avoid the excitation of torsional and longitudinal wave modes. The opposite end is clamped, to avoid undesired deformations of the beam. A reﬂective powder is employed to cover the beam, to increase the quality of the measurements, which are performed through a Polytec 3D scanner laser Doppler vibrometer (SLDV), that is able to measure and separate the 3D velocity ﬁeld u(x, y), v(x, y), and w(x, y) on the top surface of the beam. Finally, the excitation signal, provided through a KEYSIGHT 33500B waveform generator, synchronously starts with the measurement system and consists in a input function V (t) = V0w(n)sin(2πfc) with amplitude V0 = 2 V, Hann window w(n), central frequency fc = 2.12 kHz and time duration 15 ms.

[1] H. M. Hilber, T. J. R. Hughes, R. L. Taylor, Improved Numerical Dissipation for Time Integration Algorithms in Structural Dynamics, Earthquake Engineering and Structural Dynamics, vol. 5, pp. 283–292, 1977. [2] P. Rajagopal, M. Drozdz, E. A. Skelton, M. J. Lowe, R. V. Craster, On the use of absorbing layers to simulate the propagation of elastic waves in unbounded isotropic media using commercially available ﬁnite element packages, NDT & E Int. 51 (2012) 30–40.