Abstract
Breaking reciprocity has recently gained significant attention due to its broad range of applications in engineering systems. Here, we introduce the first experimental demonstration of a broadband mechanical beam waveguide, which can be reconfigured to represent wave nonreciprocity. This is achieved by using spatiotemporal stiffness modulation with piezoelectric patches in a closed-loop controller. Using a combination of analytical methods, numerical simulations, and experimental measurements, we show that contrary to the conventional shunted piezoelectrics or nonlinearity based methods, our setup is stable, less complicated, reconfigurable, and precise over a broad range of frequencies. Our reconfigurable nonreciprocal system has potential applications in phononic logic, wave diodes, energy trapping, and localization.
1. Introduction
Reciprocity is a fundamental property of various physical systems, where the transmission of a physical quantity, such as waves, between 2 points in space is symmetrical. Breaking this reciprocity offers an enhanced control over wave signal transmission and has recently become of interest in many branches of physics such as optics Miri et al. (2017); Sounas and Alù (2017), electronics Dobson (1995), thermodynamics Torrent et al. (2018); Nakai and Nagaosa (2019), electromagnetism Mahmoud et al. (2015); Caloz et al. (2018), acoustics Nassar et al. (2020); Fleury et al. (2014), and classical mechanics Sugino et al. (2020); Marconi et al. (2020). In mechanical systems, breaking wave reciprocity has numerous applications in trapping waves for efficient energy harvesting devices and designing transistors in mechanical logic circuits Cummer et al. (2016); Cullen (1958); Hadad et al. (2016); Felsen and Whitman (1970); Auld et al. (1968).
Breaking reciprocity in mechanical systems can be achieved either by (i) passively employing nonlinear elements in an asymmetric structure or (ii) actively varying material properties in space and time periodically Casimir (1945). Passive methods of breaking reciprocity require high wave amplitudes, and as a result are impractical in compact devices Sugino et al. (2020). Additionally, since passively designed structures cannot be reconfigured or reprogrammed due to their static design, creating robust devices for broad frequency ranges arises additional complexities Darabi et al. (2019); Fronk et al. (2019); Wu et al. (2018); Boechler et al. (2011). Active metamaterials, on the other hand, are reprogrammable and tunable by leveraging the active spatiotemporal modulations. Such active metamaterials offer an effective platform for breaking reciprocity Riva et al. (2019); Trainiti and Ruzzene (2016); Nassar et al. (2017), and as a result, have found many applications in (i) increasing the width of bandgaps Airoldi and Ruzzene (2011), (ii) focusing or redirecting wave propagation Celli et al. (2017); Darabi et al. (2018a); Zareei et al. (2018), (iii) changing the amplitude and phase of transmitted and reflected waves Chen et al. (2018), and (iv) one-way wave blocking and cloaking Ning et al. (2019); Darabi et al. (2018b).
In active metamaterials, active elements such as magnetoelastic Chen et al. (2016); Korivand et al. (2021), photosensitive Yannopoulos and Trunov (2009), or piezoelectric Airoldi and Ruzzene (2011); Mehrvarz et al. (2019) are used to modulate physical properties. In the photosensitive and magnetoelastic materials, the stiffness can be varied by changing the magnetic field Chen et al. (2016); Korivand et al. (2021) and temperature Yannopoulos and Trunov (2009), respectively. In piezoelectric transducers, however, the stiffness is modulated by connecting these patches to shunted circuits Airoldi and Ruzzene (2011); Marconi et al. (2020) with negative capacitance Tateo et al. (2015). Although shunted piezoelectric patches are precise Darabi et al. (2020); Sugino et al. (2020); Trainiti et al. (2019); Chen et al. (2019) and able to function in a wide frequency range, a dramatic change in the equivalent Young’s modulus only happens when the system operates very close to the unstable zones of the circuit Airoldi and Ruzzene (2011). As such, the shunted piezoelectric patches are prone to an error where a small variation of the applied negative capacitance can make the system unstable and result in a large deviation from desired values Tateo et al. (2015). Additionally, implementing a system with shunted piezoelectric patches is extremely challenging since each piezoelectric element is separately connected to its own independent circuit and this requires working with too many elements and connections.
Here, we present a closed-loop feedback control system connected to two sets of parallel piezoelectric patches 1 bonded on a host beam to modulate the beam’s stiffness. The controller here continually measures the voltage from one set of PZTs and applies the required analog signals to the other set, and as such, changes the beam’s effective stiffness Li et al. (2017); Sugino et al. (2018). Contrary to shunted PZTs with negative capacitance circuits, a closed-loop stiffness modulation system is stable and precise with reduced complexity and has functionality over a broad range of frequencies. As a result, stiffness modulation with a closed-loop control system offers a robust platform for breaking reciprocity. Here, we first analytically/numerically show that spatiotemporal (or spatial) stiffness modulation in beam results in directional (or total) bandgaps. We then experimentally exhibit the bandgaps using a closed-loop controlled system. We show that different behaviors (i.e., wave transmission, blockage, and nonreciprocal transmission) are all attainable on the same reconfigurable beam setup by only tuning the controller’s parameters.
2. Model
We consider elastic wave propagation in a thin beam with linear mass density ρ
l
, rigidity E0, second moment of inertial I0, and stiffness D0 = E0I0 where the displacement w is governed by the Euler-Bernoulli equation This straightforward equation can clearly show the concept of nonreciprocity.
Assuming a spatiotemporal modulation in stiffness with spatial modulation wavelength λ
m
and temporal modulation period T
m
, the stiffness is obtained as Nonreciprocal wave transmission: (a) Schematic of a beam structure with a spatiotemporally modulated section in the middle. The spatiotemporal modulation results in a nonreciprocal wave propagation allowing wave to propagate from right to left and blocking waves from left to right. (b) The stiffness modulation diagram for the structure is shown in space and time. The spatial/temporal modulation wavelength/period are λ
m
and T
m
.
The band diagram of a spatiotemporally modulated beam is the solution of equation (2) and can be found for different modulation parameters k
m
, ω
m
, and α
m
. Achieving different behaviors is possible by selecting proper modulation functions and simply changing the controller’s parameters accordingly (see Supplementary Information for the effects of density and stiffness modulation with or without a phase difference between them.) To avoid nonlinear effects, we assume a small fixed value of α
m
= 0.15 in our modulation. Initially assuming only a spatial modulation, that is, k
m
≠ 0, ω
m
= 0, the band diagram exhibits a full bandgap as shown in Figure 2(a). Next, by introducing and varying temporal modulation, that is, ω
m
≠ 0, we find that the bandgap becomes asymmetrical (Figure 2(b)). Furthermore, increasing the temporal modulation results in a fully asymmetrical bandgap, which has a nonreciprocal wave behavior with unidirectional wave propagation (Figure 2(c)). In the case of only temporal modulation, that is, k
m
= 0, ω
m
≠ 0, we find that the bandgaps are shifted to the wavenumbers where certain wavelengths are not allowed in the system (see Figure 2(d)) Trainiti et al. (2019). (A continuous changing of the modulation parameters (i.e. k
m
and ω
m
) can provide a better insight into spatiotemporal modulation, shown in supplementary videos A and B). Band diagrams of an aluminum beam (E = 69.9 GPa and ρ = 2700 kg.m−3) with 1 mm thickness and 2.5 cm width. (a) Complete bandgap for only space modulated beam α
m
= 0.15, k
m
= 80 m−1, and ω
m
= 0. (b) Directional and complete bandgaps in the beam with α
m
= 0.15, k
m
= 80 m−1, and ω
m
= 400 Hz. (c) Directional bandgaps in the beam with α
m
= 0.15, k
m
= 80 m−1, and ω
m
= 800 Hz. (d) Wavenumber bandgap with α
m
= 0.15, k
m
= 0, and ω
m
= 800 Hz (for more information on the effect of k
m
and ω
m
on band diagram and bandgaps, readers are referred to supplementary videos A and B).
Next, we use a closed-loop circuit with two parallel PZT patches (one actuator and one sensor) on both sides of the beam to implement spatiotemporal modulation in stiffness. Based on the Euler-Bernoulli theory, the stiffness in a beam is defined M = −Dw
xx
, where M is the moment and w
xx
is the curvature of the beam. When a PZT actuator and a PZT sensor are attached on each side of a beam extending between x
l
and x
r
, the stiffness relation is modified to Jalili (2010)
3. Numerical results
We first numerically test the effect of spatiotemporal modulation of the PZTs on an aluminum beam using Finite Element simulations in COMSOL Multiphysics. Finding the best results which have the closest results to the theoretical results is the goal of this simulation. We consider an aluminum beam with the thickness of 1 mm, width of 2.54 mm, and length of 1.5 m. We further assume 12 pairs of PZT patches perfectly bonded on both sides of the beam. Each PZT covers an area of 21 mm × 21 mm, thickness 0.55 mm, and we have a distance of 5 mm in-between them (see Figure 3(a)). A few of the parameters are general and not dependent on the experimental setup, parameters such as frequency of time modulation and the phase of modulation that can be tuned in the feedback controller. On the other hand, for the other parameters, such as frequency of space modulation we are limited to the experimental setup because it is related to the length of PZTs and the distance between them. So, we choose these numbers based on an experimentally feasible setup. A low reflection boundary is used on both sides of the beam to avoid reflecting the wave. The closed-loop parameters are selected so that the stiffness of the piezoelectric actuators vary by 15%, i.e. α
m
= 0.15. The closed-loop feedback controller based on equation (5) is used to modulate the PZT (see Figure 3(c)). The phase difference between subsequent PZTs is set to ϕ = 2π/3, which results in three PZT pairs per spatial wavelength λ
m
(Figure 3(b)). It is to be noted that the equivalent stiffness of each PZT pair, contrary to the shunted circuits, can vary continuously based on the continuous modulation signal (see Figure 3(b)). To test the setup, the wave is initiated on one side of the structure, and measurements are done on both sides of the spatiotemporally modulated section. Finally, a sweep over a range of wave frequencies is done to obtain the transfer function (i) in the absence of any modulation, in the presence of (ii) spatial modulation (i.e., ω
m
= 0, k
m
≠ 0), or (iii) spatiotemporal modulation (i.e., ω
m
≠ 0, k
m
≠ 0). (a) The schematic of the spatiotemporally modulated beam with PZT actuators (blue) and sensors (red) bonded on both sides of the beam. (b) The stiffness modulation coefficient α
m
of three consecutive PZT patches over time. Note that the phase difference between consequent patches results in the spatial modulation. (c) The schematic of the closed-loop feedback controller to implement the active modulation. The PZT sensor’s voltage is used in equation (5) to obtain actuating voltage, which is applied to the actuator PZT. (d) The numerical transfer function (ratio between transmitted wave to incident wave) without any modulation (black line) and with spatial modulation with k
m
= 80.5 m−1 for waves propagating left-to-right (red line) and waves propagating right-to-left (blue line). Since there is no time modulation, the results are identical for both LR and RL propagations. The yellow boxes indicate the bandgaps. (e) Numerically calculated transmission ratio for spatiotemporally modulated beam with f
m
= 400 Hz and k
m
= 80.5 m−1 for waves propagating RL (blue line) and LR (red line). The red boxes indicate the nonreciprocal bandgaps.
In Figure 3(d)–(e), we plot the transfer function, that is, the ratio of transmitted wave amplitude to incident wave amplitude, in the frequency range 2 − 10 kHz for the spatial/spatiotemporal modulation of the beam. Here, the transmission ratio below r t ≤ 0.1 is considered blockage because the energy level, which is proportional to the deflection square, drops more than 0.01 Khelif and Adibi (2015). Note that this transmission level corresponds to −20dB which is commonly considered as bandgap Yi et al. (2017); Zouari et al. (2018); Tang and Cheng (2017). First, in the absence of modulation, no bandgap is observed over the frequency range of interest (black line in Figure 3(d)). In the presence of spatial modulation, however, the transfer function changes and results in identical transfer functions for left-to-right (LR) and right-to-left (RL) wave of propagation (see red and blue lines in Figure 3(d), respectively). Additionally, the results show two identical bandgaps for opposite propagation direction (LR and RL) bounded between 2.5 − 4.5 kHz, marked with yellow boxes in Figure 3(d), and 7–8.5 kHz. Interestingly, the first bandgap contains a localized mode at 3.16 kHz, which dominates the transfer function response, and results in a spike in the transmission ratio in the bandgap. Next, we run the spatiotemporal modulated simulations (i.e., ω m ≠ 0, k m ≠ 0) and report the results in Figure 3(e) for the wave propagation in LR and RL directions. This figure shows the presence of three directional bandgaps (marked with red rectangles), in which the magnitude of the wave for LR is different from RL. As demonstrated, for two of these bandgaps (2.65–2.86 kHz, and 6.46–7.48 kHz), waves only travel for LR, while for the (7.97–8.69 kHz) waves propagate in the opposite direction. Furthermore, bandgaps move to lower frequencies for the RL and higher frequencies for the LR, respectively. Nonetheless, this change is not visible for the first bandgap due to the presence of localized mode in this frequency range. See Supplementary Information for the band structure diagram of FEM.
4. Experimental results
Next, to show the feasibility of having nonreciprocity with the proposed structure, we implement the time-periodic stiffness modulation of the elastic waveguide beam using an array of PZT actuators and sensors attached on a flexible aluminum beam controlled using closed-loop circuits (Figure 4). The aluminum beam has a rectangular cross-section similar to the numerical simulation, and layers of Butyl rubber at its boundaries damp the wave reflection on the structure. The waves are produced using Macro Fiber Composite (MFC) actuators that are bonded on two ends of the beam, and a laser vibrometer (Polytec CLV-2544) measures the system’s response at a distance of 5 mm before and after the modulated section (Figure 4(a)). The voltages of PZT sensors are measured using a DS1006 R&D controller board. The measured voltages are used in the closed-loop circuits in Simulink-MATLAB with reconfigurable parameters of equation (5) (see Supplementary Information for more details). As can be seen in Figure 4(a), ControlDesk software with DS1006 R&D controller board applies the control voltages to the piezodrive amplifier, and the outputs are applied to the PZT actuators. See Figure 4(b) for the experimental setup. The experimental results for the transformation ratio without any modulation are shown in Figure 4(c) with a solid black line. Introducing spatial modulation with k
m
= 80.5 m−1, the transmission ratio changes and results in two identical bandgaps for RL and LR wave propagation direction at frequencies 2.86–3.33 kHz and 3.18–3.6 kHz (shown with yellow boxes in Figure 4(c)). Note that in the analytical results, this bandgap is continuous from 2.86 to 3.6 kHz without any break in the middle. However, in the experiments similar to the numerical simulations, this bandgap breaks into two separate bandgaps due to an internal localized mode of the system. Introducing a time modulation in addition to the spatial modulation with f
m
= 400 Hz results in asymmetrical bandgaps that depend on the direction of wave propagation. Figure 4(d) shows the presence of four directional bandgaps (red boxes), in which, same as the numerical simulations, the magnitude of the wave for LR is different from RL. Interestingly, two of the bandgaps (specifically 2.84–2.9 kHz, and 3.16–3.27 kHz) allow only LR wave propagation (RL nonreciprocal bandgaps), while on the contrary, the other two bandgaps (i.e., 3.01–3.13 kHz and 3.57–3.6 kHz) allow only RL waves propagation (LR nonreciprocal bandgaps). Moreover, similar to the numerical results, bandgaps move to lower/higher frequencies for the RL/LR bandgaps. It should be noted that two of the directional bandgaps are observed only because the bandgaps are in the vicinity of the interior localized mode. Additionally, we expect the experimental result to defer from the numerics due to the following sources of error, (i) In the numerical simulations, the contacts between the PZTs and beam are assumed to be perfect contacts; however, there is an adhesive layer between the PZTs and aluminum beam that changes the effective thickness of the system but also changes the corresponding density of the system Rabinovitch and Vinson (2002). Additionally, (ii) our experiments include structural damping and nonlinear effects, which are ignored in the simulations and theory. Lastly, (iii) including the equipment’s error (e.g., signal generator, piezodrive amplifier, and controller) would further defer our experimental result from the numerical simulation. (a) Schematic of the experimental setup showing the beam, absorbing patches, PZTs sensors and actuators, piezodrives, MFC actuators, controller, and the laser. MFCs on both sides are used to introduce the wave to the system. The wave amplitude is measured before/after the first/last PZTs using a laser vibrometer. (b) Experimental setup. (c) The experimental measurements of transmission ratio without any modulation (solid black line) and with spatial modulation with k
m
= 80.5 m−1 for wave propagating LR (red) and RL (blue). Note that the transmission ratio for RL and LR wave are identical. The bandgaps are indicated with the yellow box. (d) Transmission ratio in presence of spatiotemporal modulation with k
m
= 80.5 m−1 and f
m
= 400 Hz for LR (red) and RL (blue) waves. The nonreciprocal bandgaps are identified with the red box.
5. Conclusion
In summary, we experimentally implemented a beam structure to achieve total/nonreciprocal bandgaps. For the first time, an active closed-loop feedback control system was employed using an array of PZTs bonded on both sides of an aluminum host layer to implement the stiffness modulation. At first, we investigated the changing effects of spatiotemporal parameters on the bandgaps. Then, we numerically simulated spatial/spatiotemporal stiffness modulation using the closed-loop system and validated our analytical results. Finally, we experimentally showed that actively controlled stiffness modulation using closed-loop circuits is reconfigurable and can be used for both spatial and spatiotemporal modulation and has potential applications in waveguides, diodes, phononic logic circuits, or energy localization.
Supplemental Material
Supplemental Material - Experimental demonstration of broadband reconfigurable mechanical nonreciprocity
Supplemental Material for Experimental demonstration of broadband reconfigurable mechanical nonreciprocity by Amin Mehrvarz, Mohammad Javad Khodaei, Amir Darabi, Ahmad Zareei, and Nader Jalili for Journal of Vibration and Control.
Supplemental Material
Supplemental Material
Supplemental Material
Supplemental Material - Experimental demonstration of broadband reconfigurable mechanical nonreciprocity
Supplemental Material for Experimental demonstration of broadband reconfigurable mechanical nonreciprocity by Amin Mehrvarz, Mohammad Javad Khodaei, Amir Darabi, Ahmad Zareei, and Nader Jalili for Journal of Vibration and Control.
Footnotes
Declaration of conflicting interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The authors received no financial support for the research, authorship, and/or publication of this article.
Supplemental Material
Supplemental material for this article is available online.
Note
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
