Abstract
We discuss the propagation of Rayleigh waves at the boundary of a semi-infinite elastic lattice connected to a system of gyroscopic spinners. We present the derivation of the analytical solution of the equations governing the system when the lattice is subjected to a force acting on the boundary. We show that the analytical results are in excellent agreement with the outcomes of independent finite element simulations. In addition, we investigate the influence of the load direction, frequency and gyroscopic properties of the model on the dynamic behavior of the micro-structured medium. The main result is that the response of the forced discrete system is not symmetric with respect to the point of application of the force when the effect of the gyroscopic spinners is taken into account. Accordingly, the gyroscopic lattice represents an important example of a non-reciprocal medium. Hence, it can be used in practical applications to split the energy coming from an external source into different contributions, propagating in different directions.
1 Introduction
According to their original definition, Rayleigh waves are a class of elastic waves that propagate on the surface of an infinite homogenous isotropic solid, and they are confined within a superficial region whose thickness is comparable with their wavelength (). These waves are well known by seismologists, since they are usually detected after the occurrence of an earthquake. They are also observed in other common urban activities, such as construction and demolition works and vehicular traffic, in addition to industrial processes and technologies, like mining exploration, non-destructive testing and design of electronic instruments. In the literature, they have been studied in depth especially with reference to continuous media (see, for instance, the classical treatizes by , and ).
Surface waves traveling on the free boundaries of periodic media are usually referred to as Rayleigh-Bloch waves. They are of great importance in problems concerning the dynamic propagation of cracks in discrete systems, as discussed by and for uniform media and in successive works () for non-uniform lattices. Similar localized phenomena may play a substantial role in non-uniform crack propagation as evidenced in (), where lattice dissimilarity has been shown to promote or diminish localized deformations around the faces of the crack. Trapped modes associated with Rayleigh-Bloch waves in systems incorporating periodic gratings or periodic arrays of resonators were analyzed by , , and for scalar problems, by for vector systems and by and for plates.
Recently, there has been an increasing interest in the design and fabrication of elastic media with unusual properties, often referred to as metamaterials (see the recent works by , , , , , , and amongst others). In this paper, we focus the attention on an elastic metamaterial consisting of a triangular lattice connected to a system of gyroscopic spinners. The presence of gyroscopic spinners breaks the time-reversal symmetry of the system and makes the medium non-reciprocal, as proved by . We observe that the gyroscopic effect plays the role of magnetic bias () or angular momentum () in linear non-reciprocal electromagnetic metamaterials and of circulating fluids in acoustic linear circulators ().
Firstly introduced by and later developed by , a gyroscopic elastic lattice is a tunable system, whose dispersive properties can be varied by changing the spin and precession rates of the spinners. This type of metamaterial can be utilized to force waves to propagate along a line, whose direction is defined by the geometry of the medium (). Gyroscopic spinners can also be employed to design topological insulators, where waves travel in one direction and are immune to backscattering (; ; ; ; ; ; ; ; ). Furthermore, systems with gyroscopic spinners can be used to design coatings to hide the presence of objects in a continuous or discrete medium (; ). In a recent work by , the dispersion analysis of Rayleigh waves in a semi-infinite triangular gyroscopic lattice has been carried out. The non-symmetry of the eigenmodes of the lattice’s particles at the free boundary for positive and negative values of the wave number has been linked to the non-symmetry of the system’s response to an applied force, determined by means of finite element simulations. In addition, a comparison with an effective gyroscopic continuum discussed in that paper has corroborated the results for the discrete system when low values of the wave number are considered.
In this paper, we present for the first time the analytical derivation of the displacement field in a semi-infinite elastic lattice incorporating gyroscopic spinners, focusing the attention on Rayleigh waves. Conversely, the main objective of previous papers (, ) was related to the study of topologically-protected waveforms in lattices incorporating sub-domains with different values of the parameters of the gyroscopic spinners. The analytical results of the present paper are also verified with an independent finite element code. The displacement field produced by a point load acting on the boundary of the medium and the calculation of the energy flow demonstrate that the considered system is non-reciprocal, as the response of the forced lattice is not symmetric with respect to the point of application of the concentrated load. The influence of different physical quantities on the behavior of the gyroscopic system is also investigated through a detailed parametric analysis.
While in () the wave field produced by a force on the boundary was determined numerically by using a finite element code, here the response of the system is calculated by means of a novel analytical formulation. The latter has many advantages. First, it does not require Adaptive Absorbing Layers (AAL) to prevent wave reflections at the boundaries of the computational domain; in fact, AAL are frequency-dependent and, as a consequence, their characteristic parameters need to be tuned manually every time the frequency is changed. The analytical formulation does not require the introduction of fictitious boundaries as in finite element codes. Second, the analytical formulation allows one to perform a parametric analysis quickly and efficiently, a task that is not straightforward in many finite element packages and that is performed here. Third, analytical results are necessary to validate the outcomes of finite element models. For these reasons, it is envisaged that the proposed analytical approach can be useful to the interested reader to tackle similar dynamic problems in discrete elastic systems.
2 Materials and Methods
The material under consideration consists of a semi-infinite two-dimensional triangular lattice of masses m, linked by elastic springs of stiffness γ, length L and whose inertia is negligible in comparison with m. The planar and three-dimensional representations of the lattice are shown in Figures 1A,C, respectively. In addition, each mass is attached to a gyroscopic spinner (see Figure 1B), whose configuration is described by the Euler angles θ, ϕ and ψ, denoted as the nutation, precession and spin angles, respectively. We assume that the nutation angle θ is small (), together with the displacement of the mass attached to the spinner. Each gyroscopic spinner has length l, moments of inertia about the axis of revolution and about the other two principal axes passing through the spinner’s base. Firstly derived by and , the gyricity is an independent parameter, representing the sum of the initial precession and spin rates, that remains constant throughout the motion: . We also introduce the effective gyricity and the effective mass, that both include the inertial contribution of the spinners.
FIGURE 1
In practice, such a lattice can be realized by constructing a triangular array of masses (represented, for example, by spheres) connected by thin elastic rods. At the bottom part of each mass a cylindrical hole can be drilled, where the tip of the gyroscopic spinner can be inserted. The connection needs to be frictionless, so that the spinning motion of the gyroscope is not transmitted to the mass, which can only move in the x1- or x2‐direction without rotating. The spinning motion of the gyroscopic spinner can be applied by using an electric motor; in this way, its spin rate can remain constant during the motion. Each gyroscopic spinner is pinned at the base and its axis is parallel to the ‐direction in the initial configuration.
The boundary of the semi-infinite lattice is subjected to an oscillatory force of amplitude and prescribed radian frequency (see Figure 1A), which generates Rayleigh waves localized at the boundary and bulk waves propagating inside the medium. The displacement of a generic particle, whose position is identified by the multi-index , is denoted by and is a function of time t.
We introduce the following normalizations:where the tilde symbol indicates a dimensionless quantity. In the following, the tilde is omitted for ease of notation and it is assumed that all the appearing quantities are dimensionless.
2.1 Governing Equations of the Forced Lattice in the Transient Regime
The equations of motion of a mass within the bulk of the lattice are given in normalized form as ()wherewith
Further, the dot denotes the time derivative, the vectors () areas shown in Figure 1A, and the matrix is
The governing equations for a mass belonging to the boundary of the semi-infinite lattice () arewhereis the Kronecker delta.
In addition, we assume that the lattice is at rest initially, namely:
2.2 Alternative Representations of the Governing Equations in the Transient Regime
We introduce the integer and the transient solution in the formwhere the time-dependent complex displacement amplitude satisfies the conditions of zero initial velocity and displacement. In addition, note that for
In terms of these complex displacement amplitudes, Eqs. 2,7 become:andwhere is the identity matrix. In the last equation, the complex representation for the applied force was used.
2.2.1 Laplace and Fourier Transformed Equations for the Complex Displacement Amplitudes
Next, we apply the Laplace transform in time t to Eqs. 12,13 and we use the fact that the displacement amplitudes satisfy zero initial conditions. After this, we apply the discrete Fourier transform with respect to . In what follows, denotes the Laplace and discrete Fourier transform of the solution, defined as
Here, s and are the Laplace and discrete Fourier transform parameters, respectively. The Fourier transform parameter will also be referred to as the (normalized) wave number.
For the equations in the bulk, we obtain that the complex displacement amplitude satisfiesfor , where the factor appearing throughout due to the application of the discrete Fourier transform has been cancelled (see Eq. 10). Similar conversions for Eq. 7 describing the forced boundary at lead to
2.3 Analysis of the Forced Problem in the Steady-State Regime
The transition to the steady-state regime (i.e. when ) is made by multiplying Eqs. 15,16 by the Laplace transform parameter s and taking the limit as . In this limit, we define the discrete Fourier transform of the displacement amplitude in the steady-state regime as
From Eqs. 15,16, these transformed amplitudes then satisfyfor andfor . Here, , and in what follows “” will represent a small imaginary term with positive imaginary part. The presence of this term enables one to determine how singular points of the Fourier transformed solution approach the real axis in the complex plane in passing to the steady-state regime. In this limit, the location of these singular points in the complex plane provides a way of ascertaining the direction of waves propagating in the medium relative to the load position. This is linked to the causality principle (see Section 3.3.2 in () for more details).
2.4 Solution to the Transformed Steady-State Forced Problem
The transformed amplitudes are sought in the form
Here, is such that for , since we are considering waves that decay into the bulk of the lattice.
The solution in Eq. 20 is inserted into Eq. 18 in order to find the eigensolutions for waves in the bulk; using Eq. 19 it is then possible to satisfy conditions at and find the response to the applied load. The full derivation is reported in the Supplementary Material, where it is shown that the solution has the formwhere , and are the matrices specified in Supplementary Equations S8, S7, S11, respectively.
Note, the right-hand side of Eq. 21 defines a -periodic ( -periodic) function for real when is even (odd).
If the gyricity is zero and the load acts in the horizontal (vertical) direction, the solution in Eq. 21 describes a vector function whose first component has even (odd) real and imaginary parts as functions of , whereas the second component has odd (even) real and imaginary parts. When gyricity is non-zero, these components have neither even nor odd real and imaginary parts. Such properties can be observed in Figure 2, where we report the horizontal and vertical components of the Fourier transform of the steady-state solution , when (parts (A) and (B)) and (parts (C) and (D)) for a load acting in the horizontal direction. By further inspection of the analytical solution in Eq. 21, we pose the attention on the numerator in (), belonging to and (see Supplementary Equations S9, S11), where the presence of gyricity induces a competition between terms and the amplitude of the components depends on the sign of the wave number. This feature of the solution in Eq. 21 for is responsible for the symmetry breaking in the considered problem (see also ()).
FIGURE 2
2.5 The Displacements in the Forced Lattice and Associated Wave Phenomena
The displacements in the lattice can be found from inverting the discrete Fourier transform, taking into account the periodicity properties of Eq. 21. Then, the displacements of the nodes in the lattice are given bywhere
Here, we note that , where and the complex time-dependent amplitude were introduced in Section 2.2.
2.5.1 The Rayleigh Waves
Rayleigh waves carry some of the energy produced by the point load in both directions along the boundary. The lack of symmetry in the integral kernel of Eq. 23 with respect to the zero wave number results in a disparity between amplitudes of waves outgoing from the source to the left and to the right of the load.
The Rayleigh waves are defined by the degenerate values of the wave number of for a given frequency , and these singular points are represented by the simple poles of Eq. 21. One can check that there exist no other singular points of this function for (the special case will be discussed later). On the dispersion diagram, these simple poles correspond to the intersections of the line with the curve , where is given byand represents the dispersion curve of the system associated with Rayleigh waves.1 The dispersion curve is plotted in Figure 3.
FIGURE 3
In Figure 3 we limit our attention to the intervals of periodicity for Eq. 21 mentioned in Section 2.4. Since is present in Eq. 21, as discussed in Section 2.3, the degenerate wave numbers of Eq. 21 are perturbed and are slightly shifted from the real -axis in the complex plane. The causality principle enables the new locations of these points to be determined using information concerning the group velocity (see also ()). In particular, if at the previously mentioned intersection points () the frequency and the corresponding wave number define waves propagating to the right (left) of the load. Additionally, when the frequency is , the associated degenerate value of Eq. 21 is located at (). Consequently, the dispersion diagram in Figure 3 predicts that the degenerate values required later for the computation of outgoing waves from the source are at and .
Using the above information, we can calculate the form of the Rayleigh waves produced by the point load by determining the simple poles of Eq. 23 and employing the residue theorem to compute Eq. 23 at a point located far from the application point of the oscillatory force. With this approach, it is possible to show that for where the term on the right-hand side is the Rayleigh wave produced to the right of the point load. On the other hand, we havewhere the right-hand side defines the Rayleigh wave traveling away from the load to the left. Here,
It can be verified that when
Hence, the resulting dynamic response at the boundary of the semi-infinite gyro-elastic lattice to the far left and right of the point load is different. Moreover, these outwardly-propagating boundary waves cause the nodes to follow elliptical trajectories, as shown by the eigenmode analysis developed by .
2.5.2 Bulk Wave Radiation
Waves are also radiated along the rows of the lattice in the bulk. As , the amplitudes of such waves decrease. These waves are again attributed to the singular points of Eq. 21 associated with , which also define Rayleigh waves along the boundary. The form of the waves depends on the row of the lattice considered and it is obtained following the same procedure outlined in Section 2.5.1.
For with even, the displacements to the far right of the lattice behave aswhere the term on the right-hand side describes a wave traveling along the row defined by to the right. On the other hand, we havewhere the term on the right-hand side represents a wave propagating to the left along the row . We note that for the wave numbers and the frequency , the functions , , are complex with modulus less than unity. Hence, the matrix in Eq. 21 ensures that these waves have a decreasing amplitude for increasing .
For rows in the lattice defined by odd , we recall that Eq. 21 is a -periodic function in . Hence, the inversion formula for the discrete Fourier transform is as presented in the second equation in Eq. 23. In this case, a slightly modified procedure is required to compute the waves radiated in the bulk, taking into account the singular points of the function Eq. 21 on the interval . For odd , the displacements for arewhereas for we have
We point out that the amplitudes of the waves along the odd rows of the lattice, specified in Eqs. 32,33, take into account the contributions from the wave numbers (see Figure 3).
There also exist preferential directions for energy radiation in the bulk. When the gyricity is zero, along these specific lines in the lattice, the nodal displacements decay slowly as for , as shown by . These directions can be identified by determining when , involved in Eq. 21, is complex with , for the frequency equal to . As discussed by , this effect is purely attributed to the lattice’s micro-structure and is not found in the analogous continuous model. When the gyricity is introduced, preferential directions for wave propagation can remain and increasing the gyricity causes the associated displacement amplitudes to decrease.
2.5.3 Resonant Modes
Next, we discuss the resonant case when the steady-state solution to the considered problem does not exist. We show this by investigating the derived solution in Eq. 21 along the boundary, where .
When , for either horizontal or vertical loading at the lattice boundary, the solution in Eq. 21 is bounded for and admits the following asymptotic representation:whereand denotes a constant vector with non-zero entries depending on the gyricity , the frequency and the load amplitude . Clearly, when . In this degenerate case, the asymptote near takes the formthat indicates the appearance of other singular points near in the solution in Eq. 21. In the above, if the load is horizontal (vertical) the vector ().
We note that the singular points in this asymptote approach each other as the Laplace transform parameter , with their limiting location being the real axis at . An illustration showing the typical behavior of the first component of Eq. 21 in this case is given in Figure 4A. There, the kernel has a simple pole for positive wave numbers but not for negative wave numbers, which implies that the structure does not support Rayleigh waves propagating to the left of the load on the boundary. The second component of Eq. 21 possesses the same singular features. Further to Eq. 36, in this situation can be written in the form:where represents the term that is singular at of Eq. 21 and is a bounded function for . Upon applying the inverse of the discrete Fourier transform, while the last two terms in the right-hand side of Eq. 37 give bounded contributions to the displacement, the first term leads to a function that is singular for and . Hence, the displacement at the location of the load is unbounded for .
FIGURE 4
The asymptotic representation in Eq. 36 is also in agreement with the dispersion diagram for , shown in Figure 4B. The associated horizontal line on this diagram intersects the curves along which (), represented by gray lines. These additional curves are useful in characterizing the main contribution of the integral contained in Eq. 23 in describing the lattice far-field and are connected with the appearance of waveforms and preferential directions in the bulk lattice (see also ). At on the dispersion diagram, the group velocity at the associated intersection point is zero (see Figure 4B, where at the point the group velocity is zero). Physically, such a point represents a resonance mode of the considered system, as discussed in (, Section 3.3.5). The energy associated with this mode is unable to leave the location where the force is applied. Thus, the energy density at this location becomes unbounded as . In other words, in the case the steady-state solution does not exist.
2.6 Determination of Energy Flow in the Steady-State Regime
In this section, we analyze the energy carried by the Rayleigh waves and by the waves radiated into the bulk of the gyroscopic lattice by calculating the energy flow through the boundary of a sufficiently large region of the medium, as shown in Figure 5. The considered region is the rectangle S enclosed by the half-plane boundary and by the segments (), indicated in Figure 5.
FIGURE 5
In the steady-state regime, the rate at which energy is introduced into the system through the action of an oscillating point load applied at the node , having the vector amplitude and frequency , can be computed using the formula:where the overline denotes the complex conjugate and is defined in Eq. 23. The rate at which energy flows through each segment is given by the classical formula introduced by Brillouin (see Chapter V in ):for . In Eq. 39, represents the elastic force supplied by the nodes outside the considered region and connected to the node indicated by , positioned in the immediate proximity of , within S. The set of the nodes with multi-indices located in the vicinity of , inside S, is denoted by . For the conservation of energy law, the following equality is satisfied:
In Section 3, the above formulae for the energy flow will be used to show quantitatively that the presence of gyroscopic spinners breaks the symmetry of energy propagation with respect to a vertical line passing through the application point of the load.
2.7 Finite Element Model
The results of the analytical formulation illustrated in the previous sections will be verified with a finite element model built in the commercial software Comsol Multiphysics (version 5.4).
The numerical model consists of massless truss elements and point masses inserted at the lattice’s nodes. The effect of gyricity is simulated by imposing at each node a force that is proportional to the velocity, as in the first term on the right-hand side of the governing Eq. 2. Here the computational domain is finite, while the analytical treatment developed above is for a semi-infinite medium. In particular, the lattice has dimensions . Adaptive Absorbing Layers (AAL) are introduced close to the top and vertical boundaries, in order to prevent waves from being reflected at these boundaries. AAL are created by assigning to the links located inside those regions a complex elastic modulus. A point load is applied to the middle point of the bottom boundary of the domain. The numerical simulations are performed in the time-harmonic regime.
3 Results
3.1 The Displacement Field in the Semi-Infinite Gyro-Elastic Lattice Loaded on the Boundary
In this section, we show the response of the semi-infinite lattice with embedded gyroscopic spinners to a point force applied on the boundary. In particular, we investigate how the gyricity affects the behavior of Rayleigh waves traveling along the boundary and of the waves propagating into the bulk of the medium.
Figure 6 shows the total displacement amplitude calculated at each node of the discrete system, produced by an oscillating force applied on the boundary. The results are based on the solution in Eq. 22, derived in Section 2. The force has unit amplitude and frequency , and it acts in the horizontal (vertical) direction in parts (A) and (C) (parts (B) and (D)).
FIGURE 6
In parts (A) and (B) of Figure 6 the effective gyricity is set equal to zero. It is apparent that when the response of the system is symmetric with respect to the point of application of the force. This in agreement with the classical theory on Rayleigh wave propagation in elastic lattices without gyroscopic effects. We also observe that in the vicinity of the load there exist evanescent modes, which induce localized deformations in the neighboring links. In the bulk there are preferential directions along which waves are radiated from the source, as described in Section 2.5.2. The region between these rays show outwardly propagating oscillations that decay like . Comparing parts (A) and (B), we note that in (B) there exist waves with larger amplitudes. Moreover, by looking at the wavelengths, we observe that the bulk waves generated by a horizontal force are of the shear type, while the vertical force induces bulk waves of the pressure type, as expected.
Figures 6C,D illustrate the total displacement amplitude field produced by a horizontal and vertical force, respectively, when the effective gyricity is . As in Figures 6A,B, the force has unit amplitude and frequency . The main difference between Figures 6A–D is that when the effective gyricity is non-zero the response of the system ceases to be symmetric with respect to the vertical line passing through the point of application of the force. In part (C) we note that the amplitude of Rayleigh waves propagating to the right of the force is much larger compared to that of the surface waves traveling to the left. In addition, the force appears to activate waves in its vicinity having significant amplitude and propagating at 120° to the positive horizontal axis. On the other hand, the displacement field produced by a vertical force (see part (D)) shows that the symmetry is still broken, but the amplitudes of the waves traveling to the right and to the left of the force along the boundary are now comparable.
We also point out that the non-symmetric displacement field in Figure 6 can be inverted by changing the sign of the effective gyricity .
In Figure 7 we present the total displacement amplitude fields computed by using the finite element model developed in Comsol Multiphysics and described in Section 2.7. The values of the parameters are the same as those considered in Figure 6. Comparing Figures 6, 7, we observe that the numerical and analytical results show an excellent agreement. This confirms the validity and accuracy of the analytical treatment discussed in Section 2.
FIGURE 7
3.2 The Energy Flow
Here, we report the analytical results concerning the energy flow furnished by the external force, referred to as , and the energy flow passing through each segment in Figure 5, denoted as (). The formulae for and are given in Section 2.6.
As in Section 3.1, two different values of the effective gyricity are taken, namely and . Moreover, two directions of the point load are considered, i.e. horizontal and vertical. The values of and for all the examined cases are summarized in Table 1.
TABLE 1
| Horizontal force | Vertical force | Horizontal force | Vertical force | |
|---|---|---|---|---|
| 0.214 | 0.322 | 0.119 | 0.227 | |
| 0.029 | 0.104 | 0.028 | 0.083 | |
| 0.078 | 0.057 | 0.030 | 0.031 | |
| 0.078 | 0.057 | 0.058 | 0.022 | |
| 0.029 | 0.104 | 0.003 | 0.091 | |
| 0.214 | 0.322 | 0.119 | 0.227 | |
Values of energy flows in the micro-structured lattice due to a horizontal or vertical force, when the effective gyricity is either or .
By looking at Table 1, we notice that when the effective gyricity is zero , implying that the response of the system is symmetric with respect to the application point of the force. The above equality ceases to hold when . We also note that the energy partition between vertical and horizontal boundaries depends on the gyricity and on the direction of the force. Furthermore, for any of the four cases considered in Table 1, the energy balance is satisfied, since .
3.3 Parametric Analysis
In order to assess how the response of the micro-structured medium is affected by different physical quantities, we perform a parametric analysis where we vary the direction of the force, the radian frequency of the force and the effective gyricity of the spinners. The behavior of the medium is evaluated quantitatively by calculating the percentages of the energy flows passing through the boundaries () (see Figure 5) with respect to the energy input due to the external source.
In Figure 8 we show how the energy introduced into the system by the external oscillating force is divided into two parts, propagating in opposite directions relative to the position of the point force. In particular, the circles (squares) indicate the percentages of the energy flowing to the right (left) of the force with respect to the input energy flow. In Figure 8 it is assumed that the force acts in the horizontal direction. The five diagrams correspond to ; in each diagram, several values of the frequency of the point force are considered.
FIGURE 8
In Figure 9 the values of the effective gyricity and of the frequency of the external source are identical to those considered in Figure 8, but the outcomes are obtained by applying a concentrated oscillating load acting in the vertical direction. In both Figures 8, 9, when the incoming energy is split into two equal contributions that propagate to the left and to the right of the force, both along the boundary and inside the bulk. Conversely, when we observe a non-symmetrical distribution of energy in the system. We also notice that even when the symmetry in the energy flow partition is retrieved for some specific values of the frequency, that change with the value of the effective gyricity. The crosses in both Figures 8, 9 represent the sums of the energy flow components traveling inside the system with respect to the input energy; it can be seen that, for every case considered, there is balance between input and output energy flows.
FIGURE 9
The insets in Figure 8 present the color maps of the displacement fields, computed at given gyricities and frequencies of the external force. In part (A), where , the wave pattern is clearly symmetric. When gyricity is introduced, as in part (B), the symmetry of the displacement field is broken, even if the frequency of the external force remains the same. In part (C), the wave pattern for a frequency near the intersection between the two sets of data is shown; in this case, the energy is split into two approximately equal parts traveling to the left and to the right of the load, but on each side the energy amounts propagating into the bulk and on the boundary are different. Here, it is apparent that the energy flowing into the bulk in the left-hand part of the lattice is approximately equal to the amount of energy carried in the right-hand part, where the energy is mainly concentrated on the boundary.
In Figures 10, 11 we focus the attention on Rayleigh waves, showing how the contributions of the energy flows corresponding to surface waves traveling to the right () and to the left () of the point load vary with the effective gyricity of the spinners and with the frequency of the external force. In Figure 10 (Figure 11) the force acts in the horizontal (vertical) direction. It is apparent that the direction of the force, the gyricity and the frequency all influence the response of the system and, in particular, the relative amount of total and surface energy propagating to the left and to the right of the applied force.
FIGURE 10
FIGURE 11
The insets in Figure 10 illustrate, for different values of gyricity and frequency of the excitation, the amplitudes of the displacement components and of the lattice’s nodes at the boundary and the displacement magnitude . The symmetry of the displacement field in the lattice without gyroscopic spinners, shown in Figure 10A, is broken when gyricity is incorporated into the system, as illustrated by the non-symmetric displacement profiles of Figures 10B,E.
Now, we consider a special case, where the surface waves propagating to the left of the point load exhibit a negligibly small amplitude. This can be obtained by taking the effective gyricity , the frequency of the external force and applying a horizontally-acting force (see also Figures 8, 10). The displacement field for this choice of the parameters, presented in Figure 12, clearly shows that the energy coming from the external source propagates into the bulk and practically only to the right of the force along the boundary. A similar effect could also be achieved in the resonant case , but there the steady-state solution does not exist (see Section 2.5.3). The energy flow partition for the system in the neighborhood of is investigated in Figure 8C, which shows that the energy flux distributions vary rapidly in the vicinity of .
FIGURE 12
Another interesting case is represented by the scenario where almost all the energy propagates along the boundary. This situation is shown in Figure 13, where , and the point force acts in the vertical direction (see also Figures 9, 11). From the displacement field in this figure, it is apparent that most of the energy is localized at the boundary of the medium. Moreover, in this case the discrete system shows negligible preferential directionality (see also Figure 9E).
FIGURE 13
4 Discussion
The diagrams of the displacement amplitude fields in Figures 6, 7, as well as the values of the energy flows in Table 1, show that the gyricity is capable of breaking the symmetry in the energy propagation of both Rayleigh and bulk waves propagating from the external source. This is a consequence of the non-reciprocity of the gyro-elastic lattice.2
Examining the outcomes of the parametric analysis presented in Section 3.3, in particular Figures 8, 9, it is apparent that the distribution of energy in the system strongly depends on both the effective gyricity of the spinners and the radian frequency of the external oscillating force. Generally, at low frequencies it is easier to partition the energy flow into two significantly different contributions, propagating in opposite directions. Nonetheless, for any given value of the effective gyricity, it is possible to find one or more values of at which the energy is split into two equal parts. The lowest value of this frequency increases as the effective gyricity is increased. Comparing Figures 8, 9, we observe that the energy symmetry breaking is more evident when the force acts in the horizontal direction. Of course, for an inclined direction of the force, the response of the system can be determined from the principle of superposition by summing the effects due to the horizontally- and vertically-acting forces, since the considered problem is linear.
The investigation of the individual components of the energy flows associated with Rayleigh waves, presented in Figures 10, 11 for a horizontally- and vertically-acting point force respectively, reveals that the propagation of energy along the boundary of the medium is affected significantly by the direction of the force. In the case when , the energy flow along the boundary is the same in both directions and it decreases (increases) as the frequency of the external source acting in the horizontal (vertical) direction is increased. On the other hand, when , the energy flowing to the right (left) is larger than in the opposite direction if the direction of the force is horizontal (vertical). Independently of the direction of the force, the difference between the energy flows in the two opposite directions generally decreases as the frequency is increased.
The parametric analysis discussed in Section 3.3 has been helpful in identifying special cases, where Rayleigh waves propagate only in one direction (see Figure 12) or where bulk waves have very small amplitudes compared with those of surface waves (see Figure 13).
The analytical formulation also made it possible to analyze resonant regimes where the steady-state solution cannot be reached (see Section 2.5.3).
The capability of the considered micro-structured system in creating preferential directionality in the propagation of waves can be exploited to design novel energy splitters, where the desired amount of energy propagating in a specific prescribed direction can be varied by changing the effective gyricity of the spinners. The tunability of the proposed model is an essential tool, that can be used in many practical applications where it is required to vary the output depending on the contingent needs. Important examples include electronic instruments converting mechanical energy into electric energy (and vice versa) and elastic filters that can be utilized both for protection and energy harvesting.
Funding
MN and MB gratefully acknowledge the financial support of the EU H2020 grant MSCA-IF-2016-747334-CAT-FFLAP. The program “Mobilità dei giovani ricercatori” financed by Regione Autonoma della Sardegna is also acknowledged (MB and VP).
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
All the authors contributed to the discussion of the results and the revision of the manuscript. MN had a leading role in the analytical development of the problem, GC in the numerical computations, VP in the physical interpretation of the model and MB in the conceptualization of the work.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmats.2020.602960/full#supplementary-material.
Footnotes
1.^The dispersion analysis for Rayleigh waves in a semi-infinite gyroscopic triangular lattice is discussed in detail by , where it is explicitly shown that the dispersion diagram remains symmetric.
2.^A formal proof of the non-reciprocity of the gyroscopic medium has been presented by .
References
1
AchenbachJ. D. (1973). “Wave propagation in elastic solids.” in North-holland series in applied mathematics and mechanics. Amsterdam, London: North-Holland Publishing Company, 16.
2
Al Ba’ba’aH.YuK.WangQ. (2020). Elastically-supported lattices for tunable mechanical topological insulators. Extreme Mechanics Letters. 38, 100758. 10.1016/j.eml.2020.100758
3
AntonakakisT.CrasterR. V.GuenneauS.SkeltonE. A. (2014). An asymptotic theory for waves guided by diffraction gratings or along microstructured surfaces. Proc. Math. Phys. Eng. Sci. 470, 20130467. 10.1098/rspa.2013.0467
4
ArmaniniC.Dal CorsoF.MisseroniD.BigoniD. (2017). From the elastica compass to the elastica catapult: an essay on the mechanics of soft robot arm. Proc. Math. Phys. Eng. Sci. 473, 20160870. 10.1098/rspa.2016.0870
5
BacigalupoA.De BellisM. L.GneccoG. (2019). Complex frequency band structure of periodic thermo-diffusive materials by Floquet–Bloch theory. Acta Mech. 230, 3339–3363. 10.1007/s00707-019-02416-9
6
BiL.HuJ.JiangP.KimD. H.DionneG. F.KimerlingL. C.et al (2011). On-chip optical isolation in monolithically integrated non-reciprocal optical resonators. Nat. Photon. 5, 758–762. 10.1038/nphoton.2011.270
7
BordigaG.CabrasL.PiccolroazA.BigoniD. (2018). Prestress tuning of negative refraction and wave channeling from flexural sources. Appl. Phys. Lett. 114. 10.1063/1.5084258
8
BrillouinL. (1953). Wave propagation in periodic structuresElectric filters and crystal lattices. Mineola, NY: Dover Publications, Inc.
9
BrunM.JonesI. S.MovchanA. B. (2012). Vortex-type elastic structured media and dynamic shielding. Proc. Roy. Soc. Lond. 468, 3027–3046. 10.1098/rspa.2012.0165
10
CartaG.JonesI. S.MovchanN. V.MovchanA. B.NievesM. J. (2017). Deflecting elastic prism” and unidirectional localisation for waves in chiral elastic systems. Sci. Rep. 7, 26. 10.1038/s41598-017-00054-6
11
CartaG.BrunM.MovchanA. B.MovchanN. V.JonesI. S. (2014). Dispersion properties of vortex-type monatomic lattices. Int. J. Solid Struct. 51, 2213–2225. 10.1016/j.ijsolstr.2014.02.026
12
CartaG.ColquittD. J.MovchanA. B.MovchanN. V.JonesI. S. (2020). Chiral flexural waves in structured plates: directional localisation and control. J. Mech. Phys. Solid. 137, 103866. 10.1016/j.jmps.2020.103866
13
CartaG.ColquittD. J.MovchanA. B.MovchanN. V.JonesI. S. (2019). One-way interfacial waves in a flexural plate with chiral double resonators. Philosophical Transactions of the Royal Society A. 378, 20190350. 10.1098/rsta.2019.0350
14
CartaG.NievesM. J.JonesI. S.MovchanN. V.MovchanA. B. (2018). Elastic chiral waveguides with gyro-hinges. Q. J. Mech. Appl. Math. 71, 157–185. 10.1093/qjmam/hby001
15
ColquittD. J.CrasterR. V.AntonakakisT.GuenneauS. (2015). Rayleigh-Bloch waves along elastic diffraction gratings. Proc. Math. Phys. Eng. Sci. 471, 20140465. 10.1098/rspa.2014.0465
16
FilipovE. T.TachiT.PaulinoG. H.WeitzD. A. (2015). Origami tubes assembled into stiff, yet reconfigurable structures and metamaterials. Proc. Natl. Acad. Sci. U.S.A. 112, 12321–12326. 10.1073/pnas.1509465112
17
FleuryR.SounasD. L.SieckC. F.HabermanM. R.AlùA. (2014). Sound isolation and giant linear nonreciprocity in a compact acoustic circulator. Science. 343, 516–519. 10.1126/science.1246957
18
GarauM.CartaG.NievesM. J.JonesI. S.MovchanN. V.MovchanA. B. (2018). Interfacial waveforms in chiral lattices with gyroscopic spinners. Proc. Roy. Soc. Lond. 474, 20180132. 10.1098/rspa.2018.0132
19
GarauM.NievesM. J.CartaG.BrunM. (2019). Transient response of a gyro-elastic structured medium: unidirectional waveforms and cloaking. Int. J. Eng. Sci. 143, 115–141. 10.1016/j.ijengsci.2019.05.007
20
GraffK. F. (1975). Wave motion in elastic solids. Oxford, United Kingdom: Oxford University Press.
21
HaslingerS. G.MovchanN. V.MovchanA. B.JonesI. S.CrasterR. V. (2017). Controlling flexural waves in semi-infinite platonic crystals with resonator-type scatterers. Q. J. Mech. Appl. Math. 70, 216–247. 10.1093/qjmam/hbx005
22
LeeC. H.LiG.JinG.LiuY.ZhangX. (2018). Topological dynamics of gyroscopic and Floquet lattices from Newton’s laws. Phys. Rev. B. 97, 085110. 10.1103/PhysRevB.97.085110
23
LintonC. M.McIverM. (2002). The existence of Rayleigh-Bloch surface waves. J. Fluid Mech. 470, 85–90. 10.1017/S0022112002002227
24
MarderM.GrossS. (1995). Origin of crack tip instabilities. J. Mech. Phys. Solid. 43, 1–48. 10.1016/0022-5096(94)00060-I
25
MisseroniD.ColquittD. J.MovchanA. B.MovchanN. V.JonesI. S. (2016). Cymatics for the cloaking of flexural vibrations in a structured plate. Sci. Rep. 6, 23929. 10.1038/srep23929
26
MitchellN. P.NashL. M.IrvineW. T. M. (2018). Tunable band topology in gyroscopic lattices. Phys. Rev. B. 98, 174301. 10.1103/PhysRevB.98.174301
27
MorvaridiM.CartaG.BrunM. (2018). Platonic crystal with low-frequency locally-resonant spiral structures: wave trapping, transmission amplification, shielding and edge waves. J. Mech. Phys. Solid. 121, 496–516. 10.1016/j.jmps.2018.08.017
28
NashL. M.KlecknerD.ReadA.VitelliV.TurnerA. M.IrvineW. T. (2015). Topological mechanics of gyroscopic metamaterials. Proc. Natl. Acad. Sci. U.S.A. 112, 14495–14500. 10.1073/pnas.1507413112
29
NievesM. J.CartaG.JonesI. S.MovchanN. V.MovchanA. B. (2018). Vibrations and elastic waves in chiral multi-structures. J. Mech. Phys. Solid. 121, 387–408. 10.1016/j.jmps.2018.07.020
30
NievesM. J.CartaG.PagneuxV.BrunM. (2020). Rayleigh waves in micro-structured elastic systems: non-reciprocity and energy symmetry breaking. Int. J. Eng. Sci. 156, 103365. 10.1016/j.ijengsci.2020.103365
31
NievesM. J.MovchanA. B.JonesI. S.MishurisG. S. (2013). Propagation of Slepyan’s crack in a non-uniform elastic lattice. J. Mech. Phys. Solid. 61, 1464–1488. 10.1016/j.jmps.2012.12.006
32
PiccolroazA.GorbushinN. A.MishurisG. S.NievesM. J. (2020). Dynamic phenomena and crack propagation in dissimilar elastic lattices. Int. J. Eng. Sci. 149, 103208. 10.1016/j.ijengsci.2019.103208
33
PorterR.EvansD. V. (2005). Embedded Rayleigh-Bloch surface waves along periodic rectangular arrays. Wave Motion. 43, 29–50. 10.1016/j.wavemoti.2005.05.005
34
PorterR.EvansD. V. (1999). Rayleigh-Bloch surface waves along periodic gratings and their connection with trapped modes in waveguides. J. Fluid Mech. 386, 233–258. 10.1017/S0022112099004425
35
RayleighJ. W. S. (1885). On waves propagated along the plane surface of an elastic solid. Proc. Lond. Math. Soc. 17, 4–11. 10.1112/plms/s1-17.1.4
36
SlepyanL. I. (2002). Models and phenomena in fracture Mechanics. Berlin, Germany: Springer.
37
SlepyanL. I. (2010). Wave radiation in lattice fracture. Acoust Phys. 56, 962–971. 10.1134/S1063771010060217
38
SounasD. L.CalozC.AlùA. (2013). Giant non-reciprocity at the subwavelength scale using angular momentum-biased metamaterials. Nat. Commun. 4, 2407. 10.1038/ncomms3407
39
SüsstrunkR.HuberS. D. (2015). Observation of phononic helical edge states in a mechanical topological insulator. Science349, 47–50. 10.1126/science.aab0239
40
TallaricoD.HannemaG.MiniaciM.BergaminiA.ZempA.Van DammeB. (2020). Superelement modelling of elastic metamaterials: complex dispersive properties of three-dimensional structured beams and plates. J. Sound Vib. 484, 115499. 10.1016/j.jsv.2020.115499
41
ViktorovI. A. (1967). Rayleigh and lamb waves. New York: Plenum Press.
42
WangP.LuL.BertoldiK. (2015). Topological phononic crystals with one-way elastic edge waves. Phys. Rev. Lett. 115, 104302. 10.1098/rspa.2018.0132
43
WenzelM.BursiO. S.AntoniadisI. (2020). Optimal finite locally resonant metafoundations enhanced with nonlinear negative stiffness elements for seismic protection of large storage tanks. J. Sound Vib. 483, 115488. 10.1016/j.jsv.2020.115488
Summary
Keywords
Rayleigh waves, elastic lattice, gyroscopic spinners, dispersion properties, energy flow, non-reciprocity, energy symmetry breaking
Citation
Nieves MJ, Carta G, Pagneux V and Brun M (2021) Directional Control of Rayleigh Wave Propagation in an Elastic Lattice by Gyroscopic Effects. Front. Mater. 7:602960. doi: 10.3389/fmats.2020.602960
Received
04 September 2020
Accepted
05 November 2020
Published
27 January 2021
Volume
7 - 2020
Edited by
Andrea Bacigalupo, University of Genoa, Italy
Reviewed by
Andrea Colombi, ETH Zürich, Switzerland
Enrico Radi, University of Modena and Reggio Emilia, Italy
Updates
Copyright
© 2021 Nieves, Carta, Pagneux and Brun.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: G. Carta, giorgio_carta@unica.it
This article was submitted to Mechanics of Materials, a section of the journal Frontiers in Materials
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.