Numerical Assessment of the Performance of Elastic Cloaks for Transient Flexural Waves

Abstract

A relevant application of transformation elastodynamics has shown that flexural waves in a Kirchhoff-Love plate can be diverted and channeled to cloak a region of the ambient space. To achieve the goal, an orthotropic meta-structural plate should be employed. However, the corresponding mathematical transformation leads to the presence of an unwanted strong compressive prestress, likely beyond the buckling threshold of the structure, with a set of in-plane body forces to warrant equilibrium. In addition, the plate must possess, at the same time, high bending stiffnesses, but a null twisting rigidity. With the aim of estimating the performance of cloaks modelled with approximate parameters, an in-house finite element code, based on a subparametric technique, is implemented to deal with the cloaking of transient waves in orthotropic thin plates. The tool allows us to explore the sensitivity of specific stiffness parameters that may be difficult to match in a real cloak design. In addition, the finite element code is extended to investigate a meta-plate interacting with a Winkler foundation, to confirm how the subgrade modulus should transform in the cloak region.

Introduction

State of the Art and Research Challenges

The control of elastic waves to cloak a region of the ambient space has been shown achievable by transformation elastodynamics, which provides the mechanical properties of the material surrounding the region (Milton et al., 2006; Brun et al., 2009; Norris and Shuvalov, 2011). A relevant application of this broad area is the control of transverse waves in plates, for which solutions have been proposed mainly based on two approaches: a “passive” one, where the features of the cloak are achieved by a given microstructure (Farhat et al., 2009a; Farhat et al., 2009b; Brun et al., 2014; Colquitt et al., 2014; Climente et al., 2016; Zareei and Alam, 2017; Liu and Zhu, 2019; Misseroni et al., 2019) and an “active” one, in which tunable quantities depending on the actual mechanical input are employed for the same goal (see e.g., Futhazar et al., 2015; O’Neill et al., 2015; Chen et al., 2016; Ning et al., 2020). Experiments have been also proposed to validate some of the previous theoretical proposals (Stenger et al., 2012; Misseroni et al., 2016; Darabi et al., 2018; Misseroni et al., 2019). Most of these investigations show that this technique can be feasible and likely to be employed in centimetre-size real-life systems with wave frequencies in the order of a few kHz.

Within the set of attempts based on the “passive” approach, Colquitt et al. (2014) proposed a transformation to conceive a Kirchhoff-Love plate theory for cloaking flexural waves. The resulting governing equations for the thin plate in the transformed domain involve the presence of variable (in space) bending and torsional stiffnesses, and density, in addition to the presence of in-plane body forces and prestresses in the plate. In the same paper, the general framework was specialised to the case of a square cloak composed of four trapezoidal elements (see Figure 1) embedded in an isotropic, homogeneous domain. For this geometry, a set of relationships was established to provide explicit expressions of the quantities concerned in each part of the cloak. The broadband effectiveness of the meta-material plate cloak was then assessed by means of numerical tests.

FIGURE 1

Colquitt et al. (2014): (A) initial, undeformed domain where the plate is isotropic and homogeneous; (B) transformed domain where the four trapezoids composing the cloak are numbered. The principal directions of orthotropy are represented with thin lines (the local axes at a generic point P are labelled x and y as depicted). The perimeter of trapezoid is marked with a dashed red line.

The set of listed properties demonstrates that the design and engineering of an effective cloak for flexural waves based on transformation elastodynamics of Kirchhoff-Love plates is an exceptional challenge. Actually, the presence of severe compressive prestresses leads to the conclusion that it is impossible to construct a stable, ideal cloak. Therefore, an approach based on carefully considered approximations are required to achieve the goal. In the experimental test conducted by Misseroni et al. (2016), for instance, the design of the approximated cloak has focused mainly on the match of bending stiffnesses along the directions of orthotropy.

Background

The definition of the flexural cloak is based on a transformation, say , that maps a two-dimensional subdomain occupied by a homogeneous, isotropic Kirchhoff-Love plate, to the domain where the meta-plate is to be constructed, i.e., .¹ The map is such that the transformed equation in the latter is still the governing equation for flexural waves supported by a Kirchhoff-Love plate. In our case, is a square deprived of a tiny hole in the centre (Figure 1A) that is mapped to a collection of four trapezoids (Figure 1B), such that . The length of one side of the external boundary of both and is equal to whereas that of is equal to , where ε is a small parameter.²

In both configurations, the exterior domain is occupied by the initial, homogeneous structure that possesses constant thickness H, mass density per unit volume and bending stiffness , where E is the Young’s modulus of the material and ν its Poisson’s ratio. When free-standing, the governing equation of the plate under harmonic vibrations takes the formwhere is the transverse displacement, is the biharmonic operator in and is the mass density per unit area. For the relevant case of a plate resting on a homogeneous substrate modelled as a bed of independent linear springs (Winkler model), a term is to be added on the left-hand side of Eq. 1, where is the subgrade coefficient. We emphasise that in Eq. 1, and are independent of the position.

The function is assumed to be invertible, a requisite that is met if its evaluation in the inner boundary of the cloak is non-singular as in the case under study, After having recalled that uppercase indices refer to points belonging to and that a comma indicates partial differentiation, the gradient of the transformation , its Jacobian and the tensor³are instrumental in defining the properties of the cloak as shown by Colquitt et al. (2014).

In order to have a complete overview on the established theory of meta-plate cloaks, we remind that the local governing equation of a Kirchhoff-Love plate takes the form (Timoshenko and Woinowsky-Krieger, 1959; Lekhnitskii et al., 1968)where σ is the mass density per unit area, represents the set of moments per unit length whilst and denote membrane forces and in-plane body forces per unit length, respectively, in equilibrium through fulfillment of the condition

The next step is to substitute the constitutive equations that connect to the generalised curvature as , where is the stiffness tensor which may depend on the position. The substitution in Eq. 2 yieldsan equation used by Colquitt et al. (2014) to identify the different terms corresponding to the transformation and confirm that the flexural displacement of a Kirchhoff-Love plate under an arbitrary coordinate mapping may be interpreted as a generalised plate.

The identification of terms shows that the stiffness tensor of the transformed plate corresponds tothe in-plane body forces toand the transformed density per unit area turns out to be

When the Winkler foundation is involved, the subgrade coefficient transforms similarly to the density, i.e.,a quantity that is now dependent on the position.

Focusing now on the square cloak transformation sketched in Figure 1 and, in particular, on trapezoid , ,and , whereand is the orthonormal basis associated with both co-ordinate systems and ; a is the half-length of the side of the square hole to be cloaked and is the “cloak thickness”. This information defines the six independent components of the stiffness tensor , the prestress and the in-plane body forces in the four subdomains. For trapezoid , all these quantities are reported in Appendix A. For the other parts of the cloak, they can be easily inferred by taking into account the relevant symmetries. Note that α is dimensionless while [L].

For the same trapezoid , to reach the orthotropic principal directions at a point P, the angle (anticlockwise in the domain , see Figure 1B) is now introduced. In particular, this angle identifies the principal axis x (respectively y) moving from (respectively ). Its value iswhere

In the new local system , the constitutive relationships assume the formwhere the new stiffness moduli contained therein are provided in terms of in Appendix A.

Interestingly, as shown by Colquitt et al. (2014), the twisting stiffness brought about by the transformation vanishes in all points of the transformed domain. This is a particular feature of the cloak under investigation that must be properly addressed in a meaningful design. It is worth recalling that in an isotropic plate with compact cross section, .

The cloak is subjected to traction-free boundary conditions in the inner free boundary and to continuity conditions at the interface between itself and the outer homogeneous plate.

Numerical Model and Implementation in a In-House Finite Element Code

With the aim of performing numerical simulations via the Finite Element Method (FEM) of the transient dynamic behaviour of a square cloak, a weak form of the governing equations is now obtained. The formulation can be derived by multiplying Eq. 4 by the test function and using the Galerkin method (Zienkiewicz and Taylor, 2005; Reddy, 1993). After the application of the divergence theorem to the generic sub-domain of the entire plate and use of the equilibrium Eq. 3, the final expression becomeswhere the term involving is now taken into account. In expression Eq. 11, and are the directional derivatives along the normal unit vector and tangent unit vector of the boundary of ; , are the bending and twisting moments per unit length along , respectively, whereas , describe the in-plane internal actions along the normal and tangent unit vectors; is the specific shear force on the boundary and is the suitable Sobolev space. It should be noted that in the standard procedure leading to the weak form Eq. 11, the terms and naturally disappear by integrating by part. The analysed plate is composed of various subdomains with different stiffness properties represented by the terms : the uncloaked region possesses a homogeneous stiffness whereas the cloak is described by the parameters given in Eq. 5.

Geometric and analytical approximations are performed to achieve the discretised formulation by dividing the entire domain into a finite subset of elements, i.e., , where is the number of the elements, and substituting unknown and test functions with the linear combinationsrespectively, where , , are shape functions defined on the spatial domain and indices range within a suitable interval as specified later. For each Finite Element, the final discretised formulation can be written aswhere and are the number of the shape functions that discretise w and v, respectively, and all the terms related to the in-plane prestress have been neglected.

A parametric representation is implemented in the code: the shape functions are defined in the local coordinate system of the master element , which is mapped into each element of the mesh using appropriate coordinate transformations. Since an isoparametric representation is not possible for thin-plate elements unless a previous distortion of the mesh (see for example Petera and Pittman, 1994), a subparametric representation has been adopted, so the bilinear shape functions that are used for the transformation of the master element are different from the hermitian ones used for the analytical approximation of the unknown.

The structured mesh employed in the numerical analyses is composed of quadrangular elements with four nodes in the corner of each element whose generalised displacements areleading to () 16 degrees of freedom for the element (Imbert, 1979; Petyt, 1990; Reddy, 1993; Zienkiewicz and Taylor, 2005).

Finally, after substitution and integration of the aforementioned shape functions in Eq. 1 and the assembly operation over all the Finite Elements of the mesh, the set of algebraic equations used for the spatial approximation of this problem can be written aswhere and are the mass matrix and the stiffness matrix, respectively; is the vector of the degrees of freedom whereas is the vector of generalised external loads.

The approximation in the time domain is performed with the single-step time-integration algorithm Generalized-α method (Chung and Hulbert, 1993; Chung and Hulbert, 1994) which is an evolution of the well-known Newmark method that is frequently used for transient dynamical analysis in mechanics (see Imbert, 1979; Hughes, 1987; Petyt, 1990; Zienkiewicz and Taylor, 2005). According to this method, Eq. 15 is written aswhere a generic damping matrix comes into play and the time-domain discretisationis performed with respect to the time-domain .

The Taylor’s expansions of the solution and its first time derivative are given bywhere and are positive real parameters that define the integration scheme. Eq. 18 are substituted into Eq. 16 calculated in the intermediate time value , thus obtainingwhere and are two additional parameters that will be specified later. The following implicit scheme can be used to achieve the solution , namelywhere . Equations 18 and 20 are employed to calculate the complete solution at each time step .

The Generalized-α method is completely described by the four parameters , , , that define the stability and the convergence properties of the time-integration scheme. They can be written as a function of the asymptotic spectral radius which is related to the numerical damping in the high-frequency limit as

It is worth mentioning that Eq. 15 can be profitably employed to perform time-harmonic simulations once a harmonic form for and is introduced.

The mass matrix can generally be computed in different ways: the consistent mass matrix is obtained directly from the weak Eq. 13, so the inertia contribution is assigned to each node via the shape function; the lumped mass matrix is a diagonal matrix obtained condensing the mass of each element in its nodes via equilibrium. In the ensuing simulations we considered a mixed mass matrix , where [see Hughes (1987) and Zienkiewicz and Taylor (2005)].

We implemented this numerical model in a Finite Element code, written using the software Matlab. The program deals with anisotropic and inhomogeneous thin elastic plates with different possible shapes of the domain and various kind of boundary and loading conditions. Using this code, it is also possible to perform static and modal analyses.

Numerical Simulations and Performance of the Cloak

Using the code described above, we studied the transient propagation of flexural waves in a simply-supported square plate in which a square hole is cloaked with respect to out-of-plane vibration generated by a point force varying sinusoidally. With reference to Figure 2, the homogeneous plate in is a steel one ( GPa, ,  kg/m) with thickness  mm and bending stiffness  Nm. The outer side of the cloak is such that  m while  m,  m and  m.⁴ When a Winkler-like soil comes into play in the analysis, the stiffness of the elastic substratum is represented by  N/m, qualitatively corresponding to that of a silicone material (Bigoni et al., 2008; Gei, 2008).

FIGURE 2

The boundary conditions on the outer perimeter are those of a simply supported plate, namely null transversal displacement , null rotation and null bending moment . The distance between the centre of the plate and the centre of the cloaked square hole is equal to 2.5 m. A harmonic transversal point force is applied at the centre of the plate where we assumed  N and ω varying from 50 to 600 rad/s. The information of the point load is contained in the vector load in the right-hand-side of Eqs. 15 and 16.

The mesh that we used for the transient analysis is composed of () 83,340 quadrangular elements with 335,864 total degrees of freedom; the elements are rectangular and regular outside the area of the cloak whereas the mesh is refined and composed of trapezoidal elements inside the area of the cloak. Each trapezoid of the cloak is composed of 1,080 Finite Elements.

We performed a transient dynamical analysis with the Generalized-α method (see above) for which we chose to neglect numerical damping, therefore the method reduced to a second order and unconditionally stable algorithm with which is equivalent to the constant-average acceleration method (Reddy, 1993). This choice of the parameters of the numerical algorithm leads to a stable scheme for any time interval, so we considered a time domain  ms with  ms. We assumed homogeneous initial condition, so and . The initial value of (necessary in the first step) has been calculated assuming the initial condition on and and solving Eq. 16 for . In the time-domain integration we neglected the contribution of the damping matrix, i.e., , and we assumed .

For all case studies, the displacement maps and the other relevant quantities that are illustrated are computed at a time t compatible with a propagation of the wavefront located downline of the region of the cloak, but not significantly disturbed by secondary waves reflected at the boundary of the plate.

Performance of the Cloak

The displacement map of a TTWP cloak computed at  ms for  rads is displayed in Figure 3B, to be compared with that reported in Figure 3A for a plate with an uncloaked hole. In the latter diagram, a wake is evident just on the right-hand side of the hole as well as a perturbed wave pattern on the left-hand side due to reflection caused by the hole. Those two features have been eliminated by the cloak that is able to regularise the wavefront emerging from the device into the expected circular pattern and minimise back scattering. A cloak following the transformation with prestress (Eqs. 5 and 6) would have produced a circular wavefront matching that of a homogeneous plate as presented by Colquitt et al. (2014).

FIGURE 3

In order to appraise more in detail the quality of the cloaking, Figure 4 reports the out-of-plane displacement w for the four sections sketched in Figure 2; the first three sections are transverse with respect to the propagation of the wave, at a distance from the centre of the hole of 1.25 m (Section no. 1), 1.75 m (Section no. 2) and 3 m (Section no. 3); the fourth one runs longitudinally from a point at 0.25 m just outside the right-hand boundary of the cloak to the external side of the plate. In all plots, the red line, which describes the displacements of the TTWP cloak, follows quite well the dashed line representing the response of the square, hole-less, homogeneous plate . Figure 4D shows that the approximate cloak is able to reproduce quite satisfactorily the phase of the wave of the homogeneous plate/perfect cloak at points in the region just after the cloaked region, to confirm the remark on Figure 3B that the TTWP cloak is able to reconstruct the circular wavefront.

FIGURE 4

In order to measure quantitatively the quality of the cloak we adopt the index (Colquitt et al., 2014)that is computed along the relevant section ς that is described by the co-ordinate r. Perfection corresponds to . In Eq. 22, are the displacements computed by the numerical model for the non-homogeneous case under investigation and corresponds to the solution for the homogeneous plate. With reference to the four sections displayed in Figure 4, it turns out that the quality indices are those reported in Table 1, showing that the cloaked solutions are in average 5.95 times more effective than the uncloaked ones, therefore a good result is in any case obtained.

The effect of an amount of twisting stiffness on the displacement map of a cloak TTWP is displayed in Figure 5 where several increasing values of the parameter are analysed. The plot on the far left is the one reported in Figure 3B. The detrimental effect of this parameter is readily observed. A quality factor can also be computed for the displayed cases. For Section no. 2, increases from 0.1241 (second panel) to 1.348 (far right panel) showing a constant worsening in the response of the meta-plate. In Appendix B, a strategy to minimise the twisting stiffness of a fibre-composite micro-structured plate is illustrated.

FIGURE 5

A similar analysis can be reiterated for the coupling stiffness . This term plays an unexpected important role as revealed by the simulations reported in Figure 6 where TTWP meta-plates are implemented, but possessing a whose point-wise value corresponds to 90, 50, and 10 of the theoretical one. Unexpected because, on the one hand, in a rational microstructured plate design where the principal bending stiffnesses and are matched first, the coupling term simply results as an outcome of the choices made previously; on the other hand, the mathematical transformation generally predicts values of this parameter that are remarkably larger than those that can be reasonably reachable for an orthotropic plate (see Appendix B, where the list of theoretical stiffnesses are reported for two selected cases). The three panels displayed in the figure clearly demonstrate the relevance of the coupling bending stiffness in assuring an accurate behaviour of the cloak.

FIGURE 6

The last picture of this subsection (Figure 7) refers to a TTWP cloak subjected to point-loads pulsating at different ω, to show that the approximate model is behaving well for a quite broad range of frequencies. The time of computation t differs from one picture to the other in order to compare a similar displacement pattern. For all displayed cases, the meta-plate is able to reconstruct satisfactorily the circular wavefront.

FIGURE 7

Flexural Cloaks on Winkler Substrates

In real applications, cloaks are likely grounded to a substrate which also provides the support of the object to hide. The cloaking features of a plate resting on a Winkler foundation is analysed in Figure 8 for  rads at  ms. In part b), the TTWP cloak rests on a bed of springs whose subgrade modulus is constant (). This assumption clearly worsens the performance of the device with respect to the uncloaked hole case (Figure 8A). Conversely, the position dependent stiffness , whose value in follows Eq. 7, guarantees a satisfactory behaviour of the device, as it can be noted that the wavefront reconstructs correctly after the transformed domain. Similarly to the simply-supported plate, solutions along two sections (i.e. no. 2 and no. 4 in Figure 2A) are studied in Figure 9. The quality index for the two sections are , and , , respectively. Again, the behaviour of the perfect cloak, both in phase and amplitude, is captured remarkably well by the approximate one with position dependent subgrade modulus. This shows that a correct design of a flexural cloak resting on a substrate requires the presence of a variable subgrade modulus that follows the provision of the theoretical transformation.

FIGURE 8

FIGURE 9

As a final remark of this section, one must note that the employed geometrical, elastic and loading conditions correspond to a chosen prototype geometry, but the scale of the dimensions, and the load intensity and frequency can be changed to obtain an identical behaviour on a different scaled structure subjected to different load conditions.

Concluding Remarks

The engineering of a meta-structural cloak for elastic flexural waves based on transformation elastodynamics is an exceptional challenge which theoretically would require the implementation of unfeasible compressive prestresses and in-plane body forces to warrant equilibrium. Therefore, reasonable assumptions should be adopted that are however grounded on the following main points: 1) the meta-plate has a locally orthotropic response with 2) an almost vanishing twisting stiffness, 3) spatially-varying bending stiffnesses and 4) density. In addition, the structure invariably would interact with a substrate.

With the aim of dealing with a fully open simulation tool, a FE code is developed and implemented with the specific purpose of studying transient wave propagation in locally orthotropic Kirchhoff-Love plates. A subparametric technique is adopted for spatial discretization, whereas the approximation in the time domain is performed with the single-step time-integration algorithm Generalized-α method.

The performance of a prototype meta-plate square cloak based on the approximate assumption of null prestress is then assessed parametrically by focusing specially on the role of the twisting stiffness and the coupling bending stiffness. The reason is that while the local principal bending rigidities can be matched quite easily in a microstructured plate, for those two parameters the process is much more difficult.

The simulations show that for an effective cloaking response, the twisting stiffness should be lower than 10 of the bending one for the homogeneous plate, while for the coupling bending stiffness a departure of 10 from the theoretical value already worsens with some evidence the performance of the meta-plate.

A second contribution of this work consists in the extension of the general theory of thin-plate cloaks to comprise interaction of the meta-structure with an elastic substrate via Winkler-foundation model. The conclusion is that the subgrade modulus transforms similarly to the mass density of the plate. The numerical simulations confirm this finding and clearly show that a constant modulus beneath the cloak jeopardizes dramatically the functionality of the structure.

The FE tool can be further expanded to embed inelastic responses of plate and substrate, thus enabling simulations of approximate supported meta-plates subjected to transient waves arising from seismic shocks.

References

• 1

  BigoniD.GeiM.MovchanA. B. (2008). Dynamics of a prestressed stiff layer on an elastic half space: filtering and band gap characteristics of periodic structural models derived from long-wave asymptotics. J. Mech. Phys. Solid.56, 2494–2520. 10.1016/j.jmps.2008.02.007

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  BrunM.ColquittD. J.JonesI. S.MovchanA. B.MovchanN. V. (2014). Transformation cloaking and radial approximations for flexural waves in elastic plates. New J. Phys.16, 093020. 10.1088/1367-2630/16/9/093020

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BrunM.GuenneauS.MovchanA. B. (2009). Achieving control of in-plane elastic waves. Appl. Phys. Lett.94, 061903. 10.1063/1.3068491

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  ChenY.HuJ.HuangG. (2016). A design of active elastic metamaterials for control of flexural waves using the transformation method. J. Intell. Mater. Syst. Struct.27, 1337–1347. 10.1177/1045389x15590273

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  ChungJ.HulbertG. M. (1993). A time integration algorithm for structural dynamics with improved numerical dissipation: the generalized-α method. J. Appl. Mech.60, 371–375. 10.1115/1.2900803

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  ChungJ.HulbertG. M. (1994). A family of single-step houbolt time integration algorithms for structural dynamics. Comput. Methods Appl. Mech. Eng.118, 1–11. 10.1016/0045-7825(94)90103-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  ClimenteA.TorrentD.Sánchez-DehesaJ. (2016). Analysis of flexural wave cloaks. AIP Adv.6, 121704. 10.1063/1.4968611

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  ColquittD. J.BrunM.GeiM.MovchanA. B.MovchanN. V.JonesI. S. (2014). Transformation elastodynamics and cloaking for flexural waves. J. Mech. Phys. Solid.72, 131–143. 10.1016/j.jmps.2014.07.014

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  DarabiA.ZareeiA.AlamM.-R.LeamyM. J. (2018). Experimental demonstration of an ultrabroadband nonlinear cloak for flexural waves. Phys. Rev. Lett.121, 174301. 10.1103/physrevlett.121.174301

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  FarhatM.GuenneauS.EnochS. (2009a). Ultrabroadband elastic cloaking in thin plates. Phys. Rev. Lett.103, 024301. 10.1103/PhysRevLett.103.024301

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  FarhatM.GuenneauS.EnochS.MovchanA. B. (2009b). Cloaking bending waves propagating in thin elastic plates. Phys. Rev. B79, 033102. 10.1103/PhysRevB.79.033102

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  FuthazarG.ParnellW. J.NorrisA. N. (2015). Active cloaking of flexural waves in thin plates. J. Sound Vib.356, 1–19. 10.1016/j.jsv.2015.06.023

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  GeiM. (2008). Elastic waves guided by a material interface. Eur. J. Mech. Solid.27, 328–345. 10.1016/j.euromechsol.2007.10.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  HalpinJ. C.KardosJ. L. (1976). Halpin-tsai equations: a review. Polym. Eng. Sci.16, 344–352. 10.1002/pen.760160512

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  HughesT. (1987). The finite element method - linear static and dynamic finite element analysis. Upper Saddle River, NJ: Prentice-Hall.

  • Google Scholar
  • View reference in article

• 16

  ImbertF. (1979). Analyse des structures par elements finis. Paris, France: Cepadues editions.

  • Google Scholar
  • View reference in article

• 17

  KawA. (2006). Mechanics of composite materials. Boca Raton, FL: Taylor & Francis Group.

  • Google Scholar
  • View reference in article

• 18

  KohnR. V.ShenH.VogeliusM. S.WeinsteinM. I. (2008). Cloaking via change of variables in electric impedance tomography. Inverse Probl.24, 015016. 10.1088/0266-5611/24/1/015016

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  LekhnitskiiS.TsaiS.CheronT. (1968). Anisotropic plates. New York, NY: Gordon and Breach.

  • Google Scholar
  • View reference in article

• 20

  LiuM.ZhuW. D. (2019). Nonlinear transformation-based broadband cloaking for flexural waves in elastic thin plates. J. Sound Vib.445, 270–287. 10.1016/j.jsv.2018.12.025

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  MiltonG. W.BrianeM.WillisJ. R. (2006). On cloaking for elasticity and physical equations with a transformation invariant form. New J. Phys.8, 248. 10.1088/1367-2630/8/10/248

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  MisseroniD.ColquittD. J.MovchanA. B.MovchanN. V.JonesI. S. (2016). Cymatics for the cloaking of flexural vibrations in a structures plate. Sci. Rep.6, 23929. 10.1038/srep23929

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  MisseroniD.MovchanA. B.BigoniD. (2019). Omnidirectional flexural invisibility of multiple interacting voids in vibrating elastic plates. Proc. R. Soc. A475, 2229. 10.1098/rspa.2019.0283

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  MrozA. (2011). Stability analysis of a plane, rectangular, boron-epoxy laminated plate basing on strength properties determined by different methods. Mech. Mech. Eng.15, 161–181.

  • Google Scholar
  • View reference in article

• 25

  NingL.WangY.-Z.WangY.-S. (2020). Active control cloak of the elastic wave metamaterial. Int. J. Solid Struct.202, 126–135. 10.1016/j.ijsolstr.2020.06.009

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  NorrisA. N.ShuvalovA. L. (2011). Elastic cloaking theory. Wave Motion.48, 525–538. 10.1016/j.wavemoti.2011.03.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  O’NeillJ.SelsilO.McPhedranR. C.MovchanA. B.MovchanN. V. (2015). Active cloaking of inclusions for flexural waves in thin elastic plates. Q. J. Mech. Appl. Math.68, 263–288. 10.1093/qjmam/hbv007

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  PeteraJ.PittmanJ. F. T. (1994). Isoparametric hermite elements. Int. J. Numer. Meth. Engng.37, 3489–3519. 10.1002/nme.1620372006

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  PetytM. (1990). An introduction to the finite element method. Cambridge, UK: Cambridge University Press.

  • Google Scholar
  • View reference in article

• 30

  ReddyJ. (1993). An introduction to the finite element method. 2nd Edn. New York, NY: McGraw-Hill.

  • Google Scholar
  • View reference in article

• 31

  StengerN.WilhelmM.WegenerM. (2012). Experiments on elastic cloaking in thin plates. Phys. Rev. Lett.108, 014301. 10.1103/PhysRevLett.108.014301

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  TimoshenkoS.Woinowsky-KriegerS. (1959). Theory of plates and shells. New York, NY: McGraw-Hill.

  • Google Scholar
  • View reference in article

• 33

  ZareeiA.AlamM.-R. (2017). Broadband cloaking of flexural waves. Phys. Rev. E95, 063002. 10.1103/PhysRevE.95.063002

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  ZienkiewiczO.TaylorR. (2005). The finite element method. 6th Edn. Amsterdam, Netherlands: Elsevier.

  • Google Scholar
  • View reference in article