Next Article in Journal
High-Stability Ti3C2-QDs/ZnIn2S4/Ti(IV) Flower-like Heterojunction for Boosted Photocatalytic Hydrogen Evolution
Previous Article in Journal
Non-Restoring Array Divider Using Optimized CAS Cells Based on Quantum-Dot Cellular Automata with Minimized Latency and Power Dissipation for Quantum Computing
Previous Article in Special Issue
Multiscale Innovative Materials and Structures (MIMS)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Free Vibrations of Non-Classically Damped Locally Resonant Metamaterial Plates

by
Andrea Francesco Russillo
1,
Giuseppe Failla
1,*,
Ada Amendola
2 and
Raimondo Luciano
3
1
Department of Civil, Energy, Environmental and Materials Engineering (DICEAM), University of Reggio Calabria, Via Graziella, 89124 Reggio Calabria, Italy
2
Department of Civil Engineering, University of Salerno, 84084 Fisciano, Italy
3
Department of Civil Engineering, Parthenope University of Naples, 80133 Naples, Italy
*
Author to whom correspondence should be addressed.
Nanomaterials 2022, 12(3), 541; https://doi.org/10.3390/nano12030541
Submission received: 27 December 2021 / Accepted: 28 January 2022 / Published: 5 February 2022
(This article belongs to the Special Issue Multiscale Innovative Materials and Structures)

Abstract

:
In this paper, the focus is on the free vibrations of locally resonant metamaterial plates with viscously damped resonators. Upon formulating a dynamic-stiffness model where the resonators are represented via pertinent reaction forces depending on the deflections of the attachment points, the complex eigenvalues are calculated by a contour-integral algorithm introduced in the literature for general nonlinear eigenvalue problems. The interest in the proposed approach is twofold. The dynamic-stiffness model involves a limited number of generalised coordinates compared to the nodal degrees of freedom of a standard finite-element model, and the contour-integral algorithm proves successful in evaluating all complex eigenvalues, without missing any one, with remarkable computational efficiency. Numerical results are presented for Lévy plates, but are readily extendible to other plate theories. Finally, an ad hoc dynamic-stiffness approach is formulated to calculate the frequency response of the plate under arbitrarily placed loads, which is of particular interest to investigate its elastic wave attenuation properties.

1. Introduction

There is a considerable body of recent literature on locally resonant metamaterial plates (LRMPs), i.e., plates engineered with a periodic array of small resonators. Indeed, these plates exhibit the inherent attenuation properties of elastic waves that make them ideally suitable for several applications in dynamics. Typically, the resonators may be distributed over the external surface of homogeneous or composite plates [1,2,3,4,5,6,7,8,9,10] or, alternatively, they may be embedded within the core matrix of sandwich plates [11,12,13].
Computational models of LRMPs generally rely on standard finite-element (FE) analysis. The model involves a standard FE for the plate, which is coupled with concentrated mass–spring or mass–spring–damper subsystems representing the resonators [6,8]. In some cases, the resonators are modelled by FEs as well [12,13].
A very appealing approach for modelling bare plates, i.e., plates without resonators, is a dynamic-stiffness approach, as it ensures a very accurate description of the plate dynamics with a limited number of generalised coordinates compared to a standard FE model. For instance, in a comprehensive treatment for plates of various types, Banerjee and coworkers [14,15,16,17,18,19,20,21] developed dynamic-stiffness models where a few coefficients of appropriate Fourier series expansions proved capable of representing, with remarkable accuracy, the plate dynamic response. In this context, natural frequencies of the undamped modes were calculated by the powerful Wittrick–Williams algorithm [17,22], without missing any one and including multiple ones. As for LRMPs, a pertinent dynamic-stiffness model and an ad hoc formulation of the Wittrick–Williams algorithm were recently proposed by Russillo et al. [23].
A main limitation in the existing dynamic-stiffness models for LRMPs is that no damping is considered [17,22]. Yet, damping is often intentionally introduced into LRMPs, for instance by endowing the resonators with viscous dampers [24]. However, dealing with the free vibrations of dynamic-stiffness models for damped LRMPs poses severe computational challenges, because the eigenvalues are complex and the Wittrick–Williams algorithm, i.e., the classical algorithm used in dynamic-stiffness approaches, no longer applies. Moreover, a further difficulty is that the complex eigenvalues are very close to each other, as a result of local resonance. In view of these complications, no existing study deals with the free vibrations of dynamic-stiffness models for damped LRMPs, to the best of the authors’ knowledge.
This paper addresses the free vibrations of LRMPs with viscously damped resonators, focusing on Lévy plates. Two main novelties are introduced. A dynamic-stiffness formulation is adopted, which builds on the formulation proposed by Banerjee and coworkers for bare plates [14,15,16,17,18,19,20,21] and considers the resonators via pertinent reaction forces, expressed in terms of the deflections of the attachment points by a frequency-dependent stiffness obtained from the resonator motion equations. In this way, the size of the dynamic-stiffness model depends only on the number of generalised coordinates of the plate and does not include any degrees of freedom of the resonators, with a significant computational advantage as, in general, LRMPs involve a very large number of resonators. Next, the free-vibration problem is tackled by a contour-integral algorithm [25,26,27,28], introduced a decade ago for nonlinear eigenvalue problems and used by the authors, very recently, for the free-vibration problem of damped locally resonant sandwich beams [29]. Comparing the proposed dynamic-stiffness approach to a standard FE one implemented in Abaqus demonstrates that: (1) the size of the dynamic-stiffness model is much smaller than the size of the FE model required to attain the same accuracy; (2) the contour-integral algorithm is capable of calculating accurately all the complex eigenvalues, without missing any one; (3) calculating the eigenvalues by the dynamic-stiffness approach via the contour-integral algorithm is computationally more efficient than by the FE method via standard eigensolvers. The formulation presented here for Lévy plates is readily extendible to Kirchhoff ones with arbitrary boundary conditions (BCs). Finally, the paper proposes a dynamic-stiffness approach to calculate the frequency response of the LRMP under arbitrarily placed concentrated loads, for the purpose of investigating the elastic wave attenuation properties of the plate.
The paper is organised as follows. The dynamic-stiffness model of Lévy LRMPs is described in Section 2. Details on the contour-integral algorithm are illustrated in Section 3, and the dynamic-stiffness approach to calculate the frequency response is described in Section 4. A numerical example is discussed in Section 5.

2. Dynamic-Stiffness Model

Consider the equation of motion of a Kirchhoff plate strip in Figure 1 where the rotational inertia is taken into account [14]:
D 4 w x 4 + 2 4 w x 2 y 2 + 4 w y 4 + ρ h 2 w t 2 1 12 ρ h 3 4 w x 2 t 2 + 4 w y 2 t 2 = 0
w ( x , y , t ) being the transverse deflection, h the thickness, ρ the volumetric mass density, and D the bending rigidity of the plate. A solution that enforces simply supported BCs at y = 0 and y = L is considered:
w ( x , y , t ) = m = 1 W m ( x ) sin ( α m y ) e i ω t
where α m = m π L for m = 1 , 2 , , . Replacing Equation (2) in Equation (1) yields:
d 4 W m d x 4 + ρ h 3 ω 2 12 D 2 α m 2 d 2 W m d x 2 + α m 4 ρ h ω 2 D ρ h 3 ω 2 12 D α m 2 W m = 0
Introduce the vectors collecting the generalised displacements and forces along the two unconstrained edges of the strip:
u m s = W 1 m Φ 1 m W 2 m Φ 2 m = W m ( 0 ) Φ y m ( 0 ) W m ( b ) Φ y m ( b ) ; f m s = V 1 m M 1 m V 2 m M 2 m = V x m ( 0 ) M x m ( 0 ) V x m ( b ) M x m ( b )
where Φ y m ( x ) = d W m / d x is the rotation in the x z plane and V x m ( x ) and M x m ( x ) are the shear force and bending moment per unit length along the edge of the strip. Making use of the Kirchhoff plate equations providing V x m ( x ) and M x m as functions of W m ( x ) , the following matrix relation is obtained [14]:
f m s = D m s ( ω ) u m s
where D m s is the dynamic-stiffness matrix of the single strip for m = 1 , 2 , , .
The dynamic-stiffness matrix of a plate consisting of n e 1 strips ( n e is number of lines) is assembled in a finite-element fashion, and the following equation of equilibrium holds:
f m = D m ( ω ) u m
for m = 1 , 2 , , , where f m = f m , 1 T f m , n e T T and u m = u m , 1 T u m , n e T T . Truncating the Fourier series (2) up to N terms, Equation (6) can be equivalently written in matrix form:
f 1 f N f = D 1 ( ω ) 0 0 D N ( ω ) D ( ω ) u 1 u N u
It is noticed that the displacement/force q i ( y ) along the i th line, either at the boundary or between two adjacent strips of the assembled plate, can be expanded into a sine Fourier series:
q i ( y ) = m = 1 Q i m sin ( α m y )
where the coefficients Q i m are given as:
Q i m = 2 L 0 L q i ( y ) sin ( α m y ) d y
Now, consider a resonator placed on the ith line at y = y j . Resonators attached along the ith line exert a force per unit length:
f i ( y ) = k e q ( ω ) j = 1 n s w i ( y j ) δ ( y y j )
where δ is Dirac’s delta function, n s is the number of resonators attached along the ith line, and k e q ( ω ) is the frequency-dependent stiffness of the resonator, which can be obtained from its dynamic-stiffness matrix, as demonstrated by the authors in [23,29,30,31].
By means of Equation (8), the displacement of the attachment point is:
w i ( y j ) = = 1 W i sin ( α y j )
Replacing Equation (11) in Equation (10) yields:
f i ( y ) = k e q ( ω ) j = 1 n s = 1 W i sin ( α y j ) δ ( y y j )
Applying the transformation in Equation (9) to Equation (12) gives the Fourier series coefficients associated with the reaction forces exerted by the resonators:
F i m = 2 k e q ( ω ) L j = 1 n s = 1 W i sin ( α y j ) sin ( α m y j )
Equation (13) is readily written in matrix form as:
F i 1 F i 2 = D res ( ω ) W i 1 W i 2
matrix D res being defined as follows:
( D res ) m ( ω ) = 2 k e q ( ω ) L j = 1 n s sin ( α m y j ) sin ( α y j )
Matrix D res in Equation (15) is the dynamic-stiffness matrix associated with the resonators. Now, consider Equation (15), and define the following diagonal matrix with n e submatrices:
D ¯ m ( ω ) = Γ m ( ω ) 0 0 Γ m ( ω ) Γ m ( ω ) = ( D res ) m ( ω ) 0 0 0
Notice that the ith block refers to the ith line of the assembled plate. Using Equation (16) and truncating the Fourier series (13) up to N terms, the following relation is readily obtained, which mirrors Equation (7):
f 1 f 2 f N f = D ¯ 11 ( ω ) D ¯ 12 ( ω ) D ¯ 1 N D ¯ 21 ( ω ) D ¯ 22 ( ω ) D ¯ 2 N D ¯ N 1 ( ω ) D ¯ N 2 ( ω ) D ¯ N N D r ( ω ) u 1 u 2 u N u
Finally, replacing Equation (17) for f in Equation (7) leads to the following nonlinear eigenproblem for the free-vibration response:
D ˜ ( ω ) u = D ( ω ) + D r ( ω ) u = 0
The matrix D ˜ ( ω ) is the dynamic-stiffness matrix of the Lévy LRMP.

3. Contour-Integral Algorithm

In order to compute the eigenvalues of the Lévy LRMP, it is necessary to solve the nonlinear eigenproblem given by Equation (18). The presence of viscous dampers within the resonators makes the system non-classically damped and, consequently, the eigenvalues are complex. In this case, the usual algorithm employed to solve the eigenproblems associated with dynamic-stiffness models, that is the Wittrick–Williams algorithm, is no longer applicable. Therefore, the recently introduced contour-integral algorithm formulated by Asakura and coworkers [26] is used here, for the first time, to solve this challenging problem. The application of the algorithm follows these steps:
  • Select a circle on the complex plane Γ = γ 0 + ρ 0 e i θ with centre γ 0 and radius ρ 0 with 0 θ 2 π ;
  • Compute two complex random source matrices U and V having dimensions n 0 × L 0 , n 0 being the size of the dynamic-stiffness matrix D ˜ ( ω ) in Equation (18) and L 0 the number of source vectors collected in U and V ;
  • Compute the shifted and scaled moments M k using the N 0 -point trapezoidal rule:
    S k = 1 N 0 j = 0 N 0 1 ω j γ 0 ρ 0 k + 1 D ˜ ( ω j ) 1 V , k = 0 , 1 , , 2 K 1 M k = U H S k
    with K the maximum moment degree considered for the moment and U H the Hermitian transpose of U ;
  • Construct the Hankel matrices H ^ K L 0 and H ^ K L 0 < C K L 0 × K L 0 such that:
    H ^ K L 0 = [ M i + j 2 ] i , j = 1 K H ^ K L 0 < = [ M i + j 1 ] i , j = 1 K ;
  • Perform the singular-value decomposition of H ^ K L 0 ;
  • Omit small singular-value components σ i < ϵ · max i σ i ; set the number m ˜ of remaining singular value components ( m ˜ < K L 0 ) ; construct H ^ m ˜ and H ^ m ˜ < extracting the principal submatrix with maximum index m ˜ from H ^ K L 0 and H ^ K L 0 < , i.e.,
    H ^ m ˜ = H ^ K L 0 ( 1 : m ˜ , 1 : m ˜ ) ; H ^ m ˜ < = H ^ K L 0 < ( 1 : m ˜ , 1 : m ˜ ) ;
  • Compute the eigenvalues ζ j of the linear pencil:
    H ^ m ˜ < = ζ H ^ m ˜ ;
  • Calculate the eigenvalues:
    ω j = γ 0 + ρ 0 ζ j , j = 1 , , m ˜ .
The algorithm allows calculating all the eigenvalues ω j of the nonlinear eigenproblem (18) that fall within the selected circle Γ , including multiple roots [25,26,27].

4. Frequency Response

The frequency response of a Lévy LRMP can be computed by considering the equilibrium equation:
f = D ˜ ( ω ) u
A distributed force p i ( y ) and moment m i ( y ) acting along the ith line, either at the boundary or between two adjacent strips of the assembled plate, can be expanded into a sine Fourier series by Equation (8), whose coefficients V i m and M i m are:
V i m = 2 L 0 L f i ( y ) sin ( α m y ) d y M i m = 2 L 0 L m i ( y ) sin ( α m y ) d y
The coefficients in Equation (20) are collected in the subvector f m , i = V i m , M i m T of the force vector f m in Equation (6) (represented in Figure 2), which in turn is collected in the global force vector f in Equation (19). Once the vector f is built, the Fourier coefficients of the frequency response are given as:
u = D ˜ 1 ( ω ) f
where u = u 1 T u N T T collects displacement vectors associated with the terms of the Fourier series, i.e., u m = u m , 1 T u m , n e T T with m = 1 , , N and u m , i = W i m Φ i m T . The displacement u i ( y ) and rotation ϕ i ( y ) along the i th line are simply given as:
u i ( y ) = m = 1 W i m sin ( α m y ) ϕ i ( y ) = m = 1 Φ i m sin ( α m y )
where W i m and Φ i m are respectively the coefficients of the Fourier series expansion of the displacement and rotation along the ith line.
Notice that concentrated loads can readily be handled within the framework above: for this, two strips separated by a line passing through the load application point can be considered, and the load can be modelled as a standard 1D Dirac’s delta.

5. Numerical Applications

To validate the proposed method, consider a simply supported locally resonant steel plate coupled with 2-DOF viscously damped resonators and having dimensions 0.20 m × 0.20 m × 0.01 m . The plate and resonator parameters are: Young’s modulus E = 200 GPa , Poisson’s ratio ν = 0.3 , volumetric mass density ρ = 7750 kg / m 3 , k 1 = k 2 = 10 kN / m , m 1 = m 2 = 0.01 kg , c 1 = c 2 = 0.05 Ns m 1 .
The proposed dynamic-stiffness approach is implemented by applying the contour-integral algorithm to solve the eigenvalue problem involving the dynamic-stiffness matrix in Equation (18), where the frequency-dependent stiffness of the resonators is calculated as in [30]. The parameters of the contour-integral algorithm are N 0 = 36 , L 0 = 60 , K = 15 .
The first 150 eigenvalues are computed using the proposed method and, in order to assess the accuracy, modelling the plate in Figure 3 with the FE code Abaqus using a mesh of 200 × 200 S4R5 elements. Some eigenvalues are reported in Table 1, while the complete list of the first 150 is given in Table A1, Table A2, Table A3 in Appendix A, for conciseness. The calculated eigenvalues are in excellent agreement and the maximum relative error computed is equal to ϵ = 3.65 % for the real part and ϵ = 0.78 % for the imaginary part. Furthermore, it is seen that the contour-integral algorithm is capable of capturing eigenvalues very close to each other, which may differ even by a few digits. An excellent agreement is found in terms of mode shapes, as shown in Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8. Notice that the mode shapes are complex, as they are associated with complex eigenvalues; see Table 1 and Table A1, Table A2, Table A3 in Appendix A; in particular, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8 show the dimensionless real parts, as obtained upon dividing all components by the one with the maximum absolute value.
As for the size of the two models, the proposed dynamic-stiffness model involves a 144 × 144 dynamic-stiffness matrix (corresponding to eight strips with nine generalised coordinates each), while the FE model involves 11,8791 × 11,8791 stiffness and mass matrices (corresponding to 40,401 nodes). It is seen that the eigenvalues and the maximum relative error do not change appreciably by increasing the sizes of the two models. Moreover, no significant differences are found in the eigenvalues if the parameters of the contour-integral algorithm are changed. As for the computational effort, it is noteworthy that the contour-integral algorithm is implemented in an in-house MATLAB code. For the parameters assumed, i.e., N 0 = 36 , L 0 = 60 , K = 15 , the computation time required to calculate the eigenvalues in Table 1 and Table A1, Table A2, Table A3 of Appendix A is 44.4 s. On the other hand, the corresponding eigenvalues of the finite-element model in ABAQUS are calculated by the default Lanczos eigensolver, the computation time of which is 111.2 s. Therefore, it can be concluded that the dynamic-stiffness model in conjunction with the contour-integral algorithm provides very accurate results compared to a standard FE approach, requiring a very limited size of the model and computational effort.
Finally, the dynamic-stiffness approach proposed in Section 4 is applied to calculate the transmittance of the plate in Figure 3, considering two sets of BCs: (i) all edges are simply supported; (ii) two edges are simply supported and two are free. Specifically, a unit harmonic concentrated load is applied at ( x 0 , y 0 ) = ( 0.025 , 0.10 ) , and the transmittance is evaluated as the ratio of the deflection at ( x 1 , y 1 ) = ( 0.175 , 0.10 ) to the deflection at the load application point ( x 0 , y 0 ) . For completeness, Figure 9 includes the band gaps of the corresponding infinite LRMP, as computed by the standard FE approach [32]. The transmittance varies significantly with the BCs and, within the band gaps of the corresponding infinite plate, better wave attenuation properties of the finite plate are found when two edges of the plate are free and two simply supported, compared to the case where all the edges are simply supported.

6. Discussion

This paper has addressed the free vibrations of non-classically damped locally resonant metamaterial plates, focusing on Lévy plates. The main novelties of the proposed approach can be summarized as: (1) the formulation of a reduced-order dynamic-stiffness model where multi-degree-of-freedom resonators are represented via frequency-dependent reaction forces, involving the deflection of the attachment points only and no resonator DOFs; (2) the application of the contour-integral algorithm to calculate the complex eigenvalues; (3) a dynamic-stiffness approach to the calculation of the frequency response under arbitrarily placed concentrated loads, which is of interest to investigate the elastic wave attenuation properties of the plate. Considering a simply supported plate coupled with 2-DOF viscously damped resonators, it has been demonstrated that the proposed approach provides very accurate eigenvalues and mode shapes, with a very limited size of the model and computational costs compared to a standard FE approach implemented in Abaqus. It is noteworthy that the proposed approach requires only that the dynamic-stiffness matrix of the bare plate is available. As such, it is readily extendible to locally resonant metamaterial Kirchhoff plates with arbitrary BCs, using the appropriate dynamic-stiffness matrix formulated in [18,19,21].
It is important to remark that the LRMPs under study are of interest not only as structural/mechanical components at the macroscale, but also at much smaller scales. For instance, the concept of periodic arrays of resonators coupled with a primary host system was proposed at the microscale, e.g.,: Reference [33] investigated numerically and experimentally the formation of locally resonant band gaps in a 2D surface phononic crystal with inverted conical pillars; specifically, the inverted conical pillars were deposited on a semi-infinite lithium–niobate substrate and arranged in a honeycomb lattice array for applications in low-frequency guiding, acoustic wave isolation, acoustic absorbers, and acoustic filters. Moreover, numerical and experimental investigations demonstrated the existence of band gaps at a multi-GHz frequency range in a pillar-based hypersonic 2D phononic crystal with nanoscale dimensions [34]; in this case, the fabricated phononic crystal consisted of a periodic array of nanopillars arranged according to the triangular lattice structure of the crystal. An interesting discussion on the potential applications of locally resonant phononic nanostructures with other pertinent references may be found in the work by Guo et al. [35]. These recent studies demonstrate the interest in LRMPs at various scales, for which the proposed approach may represent a valuable analysis tool.

Author Contributions

Conceptualisation, A.F.R., G.F. and R.L.; methodology, G.F. and R.L.; software, A.F.R.; validation, A.F.R. and A.A.; data curation, A.F.R. and A.A.; writing—original draft preparation, G.F. and R.L.; visualisation, A.F.R.; supervision, G.F. and R.L.; project administration, G.F.; funding acquisition, G.F. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge financial support from the Italian Ministry of Education, University and Research (MIUR) under the P.R.I.N. 2017 National Grant “Multiscale Innovative Materials and Structure” (Project Code 2017J4EAYB).

Data Availability Statement

The raw/processed data required to reproduce the above findings cannot be shared at this time as the data also forms part of an ongoing study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BCBoundary condition
DOFDegree of freedom
FEFinite element
LRMPLocally resonant metamaterial plate

Appendix A. Eigenvalues

Table A1. Complex eigenvalues (from 1 to 50) of the LRMP in Figure 3 with 2-DOF resonators, computed by means of the proposed dynamic-stiffness approach (DS) and the FE method (FE).
Table A1. Complex eigenvalues (from 1 to 50) of the LRMP in Figure 3 with 2-DOF resonators, computed by means of the proposed dynamic-stiffness approach (DS) and the FE method (FE).
ModeEigenvalue ( · 10 4 ) (DS)Eigenvalue ( · 10 4 ) (FE)
1 0.06099453 ± 0.00009052 i 0.06097390 ± 0.00009040 i
2 0.06167468 ± 0.00009470 i 0.06167134 ± 0.00009468 i
3 0.06167470 ± 0.00009470 i 0.06167133 ± 0.00009468 i
4 0.06175305 ± 0.00009518 i 0.06175079 ± 0.00009517 i
5 0.06177108 ± 0.00009529 i 0.06176983 ± 0.00009528 i
6 0.06177113 ± 0.00009529 i 0.06176980 ± 0.00009528 i
7 0.06178420 ± 0.00009537 i 0.06178295 ± 0.00009537 i
8 0.06178422 ± 0.00009537 i 0.06178293 ± 0.00009537 i
9 0.06179204 ± 0.00009542 i 0.06179124 ± 0.00009542 i
10 0.06179212 ± 0.00009542 i 0.06179117 ± 0.00009542 i
11 0.06179329 ± 0.00009543 i 0.06179227 ± 0.00009542 i
12 0.06179513 ± 0.00009544 i 0.06179431 ± 0.00009544 i
13 0.06179520 ± 0.00009544 i 0.06179425 ± 0.00009544 i
14 0.06179803 ± 0.00009546 i 0.06179724 ± 0.00009545 i
15 0.06179807 ± 0.00009546 i 0.06179761 ± 0.00009546 i
16 0.06179834 ± 0.00009546 i 0.06179774 ± 0.00009546 i
17 0.06179849 ± 0.00009546 i 0.06179721 ± 0.00009545 i
18 0.06179928 ± 0.00009547 i 0.06179866 ± 0.00009546 i
19 0.06179941 ± 0.00009547 i 0.06179855 ± 0.00009546 i
20 0.06180004 ± 0.00009547 i 0.06179963 ± 0.00009547 i
21 0.06180033 ± 0.00009547 i 0.06179971 ± 0.00009547 i
22 0.06180043 ± 0.00009547 i 0.06179929 ± 0.00009547 i
23 0.06180064 ± 0.00009547 i 0.06180020 ± 0.00009547 i
24 0.06180090 ± 0.00009548 i 0.06180034 ± 0.00009547 i
25 0.06180100 ± 0.00009548 i 0.06180053 ± 0.00009547 i
26 0.06180122 ± 0.00009548 i 0.06180054 ± 0.00009547 i
27 0.06180123 ± 0.00009548 i 0.06180059 ± 0.00009547 i
28 0.06180128 ± 0.00009548 i 0.06179999 ± 0.00009547 i
29 0.06180145 ± 0.00009548 i 0.06180092 ± 0.00009548 i
30 0.06180154 ± 0.00009548 i 0.06180095 ± 0.00009548 i
31 0.06180164 ± 0.00009548 i 0.06180109 ± 0.00009548 i
32 0.06180172 ± 0.00009548 i 0.06180121 ± 0.00009548 i
33 0.06180187 ± 0.00009548 i 0.06180123 ± 0.00009548 i
34 0.06180188 ± 0.00009548 i 0.06180127 ± 0.00009548 i
35 0.06180195 ± 0.00009548 i 0.06180133 ± 0.00009548 i
36 0.06180196 ± 0.00009549 i 0.06180134 ± 0.00009548 i
37 0.06180201 ± 0.00009548 i 0.06180141 ± 0.00009548 i
38 0.06180210 ± 0.00009548 i 0.06180079 ± 0.00009548 i
39 0.06180220 ± 0.00009548 i 0.06180159 ± 0.00009548 i
40 0.06180223 ± 0.00009549 i 0.06180160 ± 0.00009548 i
41 0.06180227 ± 0.00009549 i 0.06180166 ± 0.00009548 i
42 0.06180231 ± 0.00009549 i 0.06180178 ± 0.00009548 i
43 0.06180240 ± 0.00009549 i 0.06180179 ± 0.00009548 i
44 0.06180245 ± 0.00009549 i 0.06180181 ± 0.00009548 i
45 0.06180246 ± 0.00009549 i 0.06180184 ± 0.00009548 i
46 0.06180251 ± 0.00009549 i 0.06180187 ± 0.00009548 i
47 0.06180254 ± 0.00009549 i 0.06180193 ± 0.00009548 i
48 0.06180261 ± 0.00009549 i 0.06180195 ± 0.00009548 i
49 0.06180263 ± 0.00009549 i 0.06180200 ± 0.00009548 i
50 0.06180264 ± 0.00009549 i 0.06180201 ± 0.00009548 i
Table A2. Complex eigenvalues (from 51 to 100) of the LRMP in Figure 3 with 2-DOF resonators, computed by means of the proposed dynamic-stiffness approach (DS) and the FE method (FE).
Table A2. Complex eigenvalues (from 51 to 100) of the LRMP in Figure 3 with 2-DOF resonators, computed by means of the proposed dynamic-stiffness approach (DS) and the FE method (FE).
ModeEigenvalue ( · 10 4 ) (DS)Eigenvalue ( · 10 4 ) (FE)
51 0.06180269 ± 0.00009550 i 0.06180204 ± 0.00009548 i
52 0.06180273 ± 0.00009549 i 0.06180205 ± 0.00009548 i
53 0.06180277 ± 0.00009549 i 0.06180207 ± 0.00009548 i
54 0.06180278 ± 0.00009551 i 0.06180211 ± 0.00009548 i
55 0.06180279 ± 0.00009549 i 0.06180213 ± 0.00009548 i
56 0.06180280 ± 0.00009549 i 0.06180214 ± 0.00009548 i
57 0.06180281 ± 0.00009548 i 0.06180217 ± 0.00009548 i
58 0.06180287 ± 0.00009549 i 0.06180220 ± 0.00009548 i
59 0.06180289 ± 0.00009553 i 0.06180158 ± 0.00009548 i
60 0.06180290 ± 0.00009549 i 0.06180224 ± 0.00009548 i
61 0.06180290 ± 0.00009546 i 0.06180227 ± 0.00009548 i
62 0.06180294 ± 0.00009549 i 0.06180230 ± 0.00009549 i
63 0.06180297 ± 0.00009548 i 0.06180222 ± 0.00009548 i
64 0.16086702 ± 0.00063822 i 0.16083957 ± 0.00063771 i
65 0.16167059 ± 0.00065236 i 0.16166713 ± 0.00065230 i
66 0.16167061 ± 0.00065236 i 0.16166711 ± 0.00065230 i
67 0.16175137 ± 0.00065368 i 0.16174909 ± 0.00065365 i
68 0.16176966 ± 0.00065398 i 0.16176840 ± 0.00065396 i
69 0.16176970 ± 0.00065398 i 0.16176838 ± 0.00065396 i
70 0.16178289 ± 0.00065420 i 0.16178165 ± 0.00065418 i
71 0.16178291 ± 0.00065420 i 0.16178162 ± 0.00065418 i
72 0.16179077 ± 0.00065433 i 0.16178997 ± 0.00065431 i
73 0.16179085 ± 0.00065433 i 0.16179100 ± 0.00065433 i
74 0.16179202 ± 0.00065435 i 0.16178991 ± 0.00065431 i
75 0.16179387 ± 0.00065438 i 0.16179304 ± 0.00065436 i
76 0.16179394 ± 0.00065438 i 0.16179299 ± 0.00065436 i
77 0.16179677 ± 0.00065442 i 0.16179599 ± 0.00065441 i
78 0.16179681 ± 0.00065442 i 0.16179637 ± 0.00065442 i
79 0.16179708 ± 0.00065443 i 0.16179649 ± 0.00065442 i
80 0.16179723 ± 0.00065443 i 0.16179730 ± 0.00065443 i
81 0.16179803 ± 0.00065444 i 0.16179742 ± 0.00065443 i
82 0.16179816 ± 0.00065445 i 0.16179596 ± 0.00065441 i
83 0.16179879 ± 0.00065446 i 0.16179838 ± 0.00065445 i
84 0.16179908 ± 0.00065446 i 0.16179845 ± 0.00065445 i
85 0.16179918 ± 0.00065446 i 0.16179873 ± 0.00065445 i
86 0.16179939 ± 0.00065447 i 0.16179895 ± 0.00065446 i
87 0.16179965 ± 0.00065447 i 0.16179910 ± 0.00065446 i
88 0.16179975 ± 0.00065447 i 0.16179928 ± 0.00065446 i
89 0.16179997 ± 0.00065447 i 0.16179929 ± 0.00065446 i
90 0.16179997 ± 0.00065447 i 0.16179932 ± 0.00065446 i
91 0.16180003 ± 0.00065448 i 0.16179955 ± 0.00065447 i
92 0.16180020 ± 0.00065448 i 0.16179805 ± 0.00065444 i
93 0.16180029 ± 0.00065448 i 0.16179969 ± 0.00065447 i
94 0.16180039 ± 0.00065448 i 0.16179984 ± 0.00065447 i
95 0.16180047 ± 0.00065448 i 0.16179996 ± 0.00065447 i
96 0.16180062 ± 0.00065448 i 0.16179999 ± 0.00065447 i
97 0.16180063 ± 0.00065448 i 0.16180003 ± 0.00065448 i
98 0.16180070 ± 0.00065449 i 0.16180008 ± 0.00065448 i
99 0.16180070 ± 0.00065448 i 0.16180010 ± 0.00065448 i
100 0.16180076 ± 0.00065449 i 0.16180016 ± 0.00065448 i
Table A3. Complex eigenvalues (from 101 to 150) of the LRMP in Figure 3 with 2-DOF resonators, computed by means of the proposed dynamic-stiffness approach (DS) and the FE method (FE).
Table A3. Complex eigenvalues (from 101 to 150) of the LRMP in Figure 3 with 2-DOF resonators, computed by means of the proposed dynamic-stiffness approach (DS) and the FE method (FE).
ModeEigenvalue ( · 10 4 ) (DS)Eigenvalue ( · 10 4 ) (FE)
101 0.16180085 ± 0.00065449 i 0.16180033 ± 0.00065448 i
102 0.16180095 ± 0.00065449 i 0.16180034 ± 0.00065448 i
103 0.16180098 ± 0.00065449 i 0.16180035 ± 0.00065448 i
104 0.16180103 ± 0.00065449 i 0.16180041 ± 0.00065448 i
105 0.16180106 ± 0.00065449 i 0.16180053 ± 0.00065448 i
106 0.16180119 ± 0.00065449 i 0.16180054 ± 0.00065448 i
107 0.16180121 ± 0.00065449 i 0.16180056 ± 0.00065448 i
108 0.16180125 ± 0.00065450 i 0.16180060 ± 0.00065448 i
109 0.16180126 ± 0.00065450 i 0.16180062 ± 0.00065449 i
110 0.16180129 ± 0.00065449 i 0.16180068 ± 0.00065449 i
111 0.16180136 ± 0.00065450 i 0.16180069 ± 0.00065449 i
112 0.16180136 ± 0.00065450 i 0.16180074 ± 0.00065449 i
113 0.16180139 ± 0.00065450 i 0.16180075 ± 0.00065449 i
114 0.16180146 ± 0.00065450 i 0.16180078 ± 0.00065449 i
115 0.16180148 ± 0.00065450 i 0.16180079 ± 0.00065449 i
116 0.16180152 ± 0.00065450 i 0.16180081 ± 0.00065449 i
117 0.16180153 ± 0.00065450 i 0.16180085 ± 0.00065449 i
118 0.16180155 ± 0.00065450 i 0.16180087 ± 0.00065449 i
119 0.16180156 ± 0.00065450 i 0.16180089 ± 0.00065449 i
120 0.16180161 ± 0.00065450 i 0.16180092 ± 0.00065449 i
121 0.16180162 ± 0.00065450 i 0.16180096 ± 0.00065449 i
122 0.16180164 ± 0.00065450 i 0.16180098 ± 0.00065449 i
123 0.16180165 ± 0.00065450 i 0.16180099 ± 0.00065449 i
124 0.16180169 ± 0.00065450 i 0.16180102 ± 0.00065449 i
125 0.16180172 ± 0.00065450 i 0.16180105 ± 0.00065449 i
126 0.16180228 ± 0.00065458 i 0.16179967 ± 0.00065447 i
127 0.38637283 ± 0.00015015 i 0.38188243 ± 0.00015078 i
128 0.94979965 ± 0.00013164 i 0.93953762 ± 0.00013174 i
129 0.94979965 ± 0.00013164 i 0.93953552 ± 0.00013173 i
130 1.51582776 ± 0.00012963 i 1.48995104 ± 0.00012971 i
131 1.89304723 ± 0.00012909 i 1.87023849 ± 0.00012922 i
132 1.89304723 ± 0.00012909 i 1.87044411 ± 0.00012920 i
133 2.45835312 ± 0.00012860 i 2.41032791 ± 0.00012871 i
134 2.45835312 ± 0.00012860 i 2.41037235 ± 0.00012874 i
135 3.21092421 ± 0.00012816 i 3.15927753 ± 0.00012834 i
136 3.21092421 ± 0.00012816 i 3.15939426 ± 0.00012844 i
137 3.39884491 ± 0.00012807 i 3.31429443 ± 0.00012829 i
138 3.77441335 ± 0.00012790 i 3.68814316 ± 0.00012808 i
139 3.77441335 ± 0.00012790 i 3.68894501 ± 0.00012809 i
140 4.71172216 ± 0.00012751 i 4.57325706 ± 0.00012775 i
141 4.71172217 ± 0.00012751 i 4.57372169 ± 0.00012786 i
142 4.89890560 ± 0.00012743 i 4.79083290 ± 0.00012797 i
143 4.89890560 ± 0.00012744 i 4.79105072 ± 0.00012793 i
144 5.45989728 ± 0.00012722 i 5.30883064 ± 0.00012762 i
145 5.45989728 ± 0.00012722 i 5.30934525 ± 0.00012751 i
146 6.02005043 ± 0.00012701 i 5.80780907 ± 0.00012723 i
147 6.39302010 ± 0.00012688 i 6.17382736 ± 0.00012755 i
148 6.39302011 ± 0.00012688 i 6.17507515 ± 0.00012727 i
149 6.95177672 ± 0.00012667 i 6.74670029 ± 0.00012756 i
150 6.95177673 ± 0.00012668 i 6.74753723 ± 0.00012768 i

References

  1. Xiao, Y.; Wen, J.; Wen, X. Sound transmission loss of metamaterial-based thin plates with multiple subwavelength arrays of attached resonators. J. Sound Vib. 2012, 331, 5408–5423. [Google Scholar] [CrossRef]
  2. Xiao, Y.; Wen, J.; Yu, D.; Wen, X. Flexural wave propagation in beams with periodically attached vibration absorbers: Band-gap behavior and band formation mechanisms. J. Sound Vib. 2013, 332, 867–893. [Google Scholar] [CrossRef]
  3. Xiao, Y.; Wen, J.; Huang, L.; Wen, X. Analysis and experimental realization of locally resonant phononic plates carrying a periodic array of beam-like resonators. J. Phys. Appl. Phys. 2013, 47, 045307. [Google Scholar] [CrossRef]
  4. Peng, H.; Pai, P.F.; Deng, H. Acoustic multi-stopband metamaterial plates design for broadband elastic wave absorption and vibration suppression. Int. J. Mech. Sci. 2015, 103, 104–114. [Google Scholar] [CrossRef]
  5. He, Z.C.; Xiao, X.; Li, E. Design for structural vibration suppression in laminate acoustic metamaterials. Compos. Part B Eng. 2017, 131, 237–252. [Google Scholar] [CrossRef] [Green Version]
  6. Wang, T.; Sheng, M.; Guo, Z.; Qin, Q. Acoustic characteristics of damped metamaterial plate with parallel attached resonators. Arch. Mech. 2017, 69, 29–52. [Google Scholar]
  7. Gao, P.; Climente, A.; Sánchez-Dehesa, J.; Wu, L. Single-phase metamaterial plates for broadband vibration suppression at low frequencies. J. Sound Vib. 2019, 444, 108–126. [Google Scholar] [CrossRef]
  8. Miranda, E., Jr.; Nobrega, E.; Ferreira, A.; Dos Santos, J. Flexural wave band gaps in a multi-resonator elastic metamaterial plate using Kirchhoff-Love theory. Mech. Syst. Signal Process. 2019, 116, 480–504. [Google Scholar] [CrossRef]
  9. Xiao, X.; He, Z.C.; Li, E.; Cheng, A.G. Design multi-stopband laminate acoustic metamaterials for structural-acoustic coupled system. Mech. Syst. Signal Process. 2019, 115, 418–433. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, T. Tunable band gaps in an inertant metamaterial plate with two-degree-of-freedom local resonance. Phys. Lett. A 2020, 384, 126420. [Google Scholar] [CrossRef]
  11. Song, Y.; Wen, J.; Tian, H.; Lu, X.; Li, Z.; Feng, L. Vibration and sound properties of metamaterial sandwich panels with periodically attached resonators: Simulation and experiment study. J. Sound Vib. 2020, 489, 115644. [Google Scholar] [CrossRef]
  12. Li, J.; Fan, X.; Li, F. Numerical and experimental study of a sandwich-like metamaterial plate for vibration suppression. Compos. Struct. 2020, 238, 111969. [Google Scholar] [CrossRef]
  13. Wang, Q.; Li, J.; Zhang, Y.; Xue, Y.; Li, F. Bandgap properties in metamaterial sandwich plate with periodically embedded plate-type resonators. Mech. Syst. Signal Process. 2021, 151, 107375. [Google Scholar] [CrossRef]
  14. Boscolo, M.; Banerjee, J. Dynamic stiffness elements and their applications for plates using first order shear deformation theory. Comput. Struct. 2011, 89, 395–410. [Google Scholar] [CrossRef]
  15. Boscolo, M.; Banerjee, J. Dynamic stiffness formulation for composite Mindlin plates for exact modal analysis of structures. Part I: Theory. Comput. Struct. 2012, 96, 61–73. [Google Scholar] [CrossRef]
  16. Boscolo, M.; Banerjee, J.R. Dynamic stiffness formulation for composite Mindlin plates for exact modal analysis of structures. Part II: Results and applications. Comput. Struct. 2012, 96–97, 74–83. [Google Scholar] [CrossRef] [Green Version]
  17. Banerjee, J.; Papkov, S.; Liu, X.; Kennedy, D. Dynamic stiffness matrix of a rectangular plate for the general case. J. Sound Vib. 2015, 342, 177–199. [Google Scholar] [CrossRef]
  18. Liu, X.; Banerjee, J. An exact spectral-dynamic stiffness method for free flexural vibration analysis of orthotropic composite plate assemblies—Part I: Theory. Compos. Struct. 2015, 132, 1274–1287. [Google Scholar] [CrossRef] [Green Version]
  19. Liu, X.; Banerjee, J.R. An exact spectral-dynamic stiffness method for free flexural vibration analysis of orthotropic composite plate assemblies—Part II: Applications. Compos. Struct. 2015, 132, 1288–1302. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, X.; Banerjee, J. Free vibration analysis for plates with arbitrary boundary conditions using a novel spectral-dynamic stiffness method. Comput. Struct. 2016, 164, 108–126. [Google Scholar] [CrossRef] [Green Version]
  21. Liu, X.; Kassem, H.; Banerjee, J. An exact spectral dynamic stiffness theory for composite plate-like structures with arbitrary non-uniform elastic supports, mass attachments and coupling constraints. Compos. Struct. 2016, 142, 140–154. [Google Scholar] [CrossRef] [Green Version]
  22. Wittrick, W.; Williams, F. A general algorithm for computing natural frequencies of elastic structures. Q. J. Mech. Appl. Math. 1971, 24, 263–284. [Google Scholar] [CrossRef]
  23. Russillo, A.F.; Failla, G. On the free vibrations of locally-resonant structures. Comput. Struct. 2020, 241, 106356. [Google Scholar] [CrossRef]
  24. Hussein, M.I.; Leamy, M.J.; Ruzzene, M. Dynamics of phononic materials and structures: Historical origins, recent progress, and future outlook. Appl. Mech. Rev. 2014, 66, 040802. [Google Scholar] [CrossRef]
  25. Sakurai, T.; Sugiura, H. A projection method for generalized eigenvalue problems using numerical integration. J. Comput. Appl. Math. 2003, 159, 119–128. [Google Scholar] [CrossRef]
  26. Asakura, J.; Sakurai, T.; Tadano, H.; Ikegami, T.; Kimura, K. A numerical method for nonlinear eigenvalue problems using contour integrals. JSIAM Lett. 2009, 1, 52–55. [Google Scholar] [CrossRef] [Green Version]
  27. Ikegami, T.; Sakurai, T.; Nagashima, U. A filter diagonalization for generalized eigenvalue problems based on the Sakurai–Sugiura projection method. J. Comput. Appl. Math. 2010, 233, 1927–1936. [Google Scholar] [CrossRef] [Green Version]
  28. Sakurai, T.; Futamura, Y.; Tadano, H. Efficient Parameter Estimation and Implementation of a Contour Integral-Based Eigensolver. J. Algorithms Comput. Technol. 2013, 7, 249–269. [Google Scholar] [CrossRef]
  29. Russillo, A.F.; Failla, G.; Fraternali, F. Free and forced vibrations of damped locally-resonant sandwich beams. Eur. J. Mech.—A/Solids 2021, 86, 104188. [Google Scholar] [CrossRef]
  30. Failla, G.; Santoro, R.; Burlon, A.; Russillo, A.F. An exact approach to the dynamics of locally-resonant beams. Mech. Res. Commun. 2020, 103, 103460. [Google Scholar] [CrossRef]
  31. Russillo, A.F.; Failla, G. A novel reduced-order dynamic-stiffness formulation for locally resonant metamaterial plates. Compos. Struct. 2022, 280, 114811. [Google Scholar] [CrossRef]
  32. Åberg, M.; Gudmundson, P. The usage of standard finite element codes for computation of dispersion relations in materials with periodic microstructure. J. Acoust. Soc. Am. 1997, 102, 2007–2013. [Google Scholar] [CrossRef]
  33. Hsu, J.C.; Lin, F.S. Measurement of locally resonant band gaps in a surface phononic crystal with inverted conical pillars. Jpn. J. Appl. Phys. 2018, 57, 07LB01. [Google Scholar] [CrossRef]
  34. Yudistira, D.; Boes, A.; Graczykowski, B.; Alzina, F.; Yeo, L.Y.; Sotomayor Torres, C.M.; Mitchell, A. Nanoscale pillar hypersonic surface phononic crystals. Phys. Rev. B 2016, 94, 094304. [Google Scholar] [CrossRef] [Green Version]
  35. Guo, Y.; Schubert, M.; Dekorsy, T. Finite element analysis of surface modes in phononic crystal waveguides. J. Appl. Phys. 2016, 119, 124302. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Generalised displacements and forces for a single plate strip with simply supported edges parallel with the x-axis.
Figure 1. Generalised displacements and forces for a single plate strip with simply supported edges parallel with the x-axis.
Nanomaterials 12 00541 g001
Figure 2. Plate consisting of n e 1 strips: (a) displacements and (b) generalised forces.
Figure 2. Plate consisting of n e 1 strips: (a) displacements and (b) generalised forces.
Nanomaterials 12 00541 g002
Figure 3. Simply supported LRMP with mass–spring–damper subsystems modelling 2-DOF resonators.
Figure 3. Simply supported LRMP with mass–spring–damper subsystems modelling 2-DOF resonators.
Nanomaterials 12 00541 g003
Figure 4. Simply supported LRMP in Figure 3, real part of Mode 127: (a) dynamic-stiffness approach and (b) finite-element method.
Figure 4. Simply supported LRMP in Figure 3, real part of Mode 127: (a) dynamic-stiffness approach and (b) finite-element method.
Nanomaterials 12 00541 g004
Figure 5. Simply supported LRMP in Figure 3, real part of Mode 130: (a) dynamic-stiffness approach and (b) finite-element method.
Figure 5. Simply supported LRMP in Figure 3, real part of Mode 130: (a) dynamic-stiffness approach and (b) finite-element method.
Nanomaterials 12 00541 g005
Figure 6. Simply supported LRMP in Figure 3, real part of Mode 144: (a) dynamic-stiffness approach and (b) finite-element method.
Figure 6. Simply supported LRMP in Figure 3, real part of Mode 144: (a) dynamic-stiffness approach and (b) finite-element method.
Nanomaterials 12 00541 g006
Figure 7. Simply supported LRMP in Figure 3, real part of Mode 150: (a) dynamic-stiffness approach and (b) finite-element method.
Figure 7. Simply supported LRMP in Figure 3, real part of Mode 150: (a) dynamic-stiffness approach and (b) finite-element method.
Nanomaterials 12 00541 g007
Figure 8. Simply supported LRMP in Figure 3, view along the line x = y of the real part of the modes: (a) 127, (b) 130, (c) 144, and (d) 150; dynamic-stiffness approach (black continuous line) and finite-element method in ABAQUS (black dots).
Figure 8. Simply supported LRMP in Figure 3, view along the line x = y of the real part of the modes: (a) 127, (b) 130, (c) 144, and (d) 150; dynamic-stiffness approach (black continuous line) and finite-element method in ABAQUS (black dots).
Nanomaterials 12 00541 g008
Figure 9. LRMP in Figure 3, amplitude ratio of the deflection frequency response at ( 0.175 , 0.10 ) to the deflection frequency response at the application point ( 0.025 , 0.10 ) of a unit harmonic force, calculated by the dynamic-stiffness approach (black continuous line) and FE model in ABAQUS (black dots), considering: (a) all edges simply supported; (b) two edges simply supported and two edges free.
Figure 9. LRMP in Figure 3, amplitude ratio of the deflection frequency response at ( 0.175 , 0.10 ) to the deflection frequency response at the application point ( 0.025 , 0.10 ) of a unit harmonic force, calculated by the dynamic-stiffness approach (black continuous line) and FE model in ABAQUS (black dots), considering: (a) all edges simply supported; (b) two edges simply supported and two edges free.
Nanomaterials 12 00541 g009
Table 1. Complex eigenvalues (from 1 to 150) of the LRMP in Figure 3 with 2-DOF resonators, computed by means of the proposed dynamic-stiffness approach (DS) and the FE method (FE).
Table 1. Complex eigenvalues (from 1 to 150) of the LRMP in Figure 3 with 2-DOF resonators, computed by means of the proposed dynamic-stiffness approach (DS) and the FE method (FE).
ModeEigenvalue ( · 10 4 ) (DS)Eigenvalue ( · 10 4 ) (FE)
1 0.06099453 ± 0.00009052 i 0.06097390 ± 0.00009040 i
2 0.06167468 ± 0.00009470 i 0.06167134 ± 0.00009468 i
3 0.06167470 ± 0.00009470 i 0.06167133 ± 0.00009468 i
4 0.06175305 ± 0.00009518 i 0.06175079 ± 0.00009517 i
5 0.06177108 ± 0.00009529 i 0.06176983 ± 0.00009528 i
6 0.06177113 ± 0.00009529 i 0.06176980 ± 0.00009528 i
7 0.06178420 ± 0.00009537 i 0.06178295 ± 0.00009537 i
8 0.06178422 ± 0.00009537 i 0.06178293 ± 0.00009537 i
9 0.06179204 ± 0.00009542 i 0.06179124 ± 0.00009542 i
10 0.06179212 ± 0.00009542 i 0.06179117 ± 0.00009542 i
11 0.06179329 ± 0.00009543 i 0.06179227 ± 0.00009542 i
12 0.06179513 ± 0.00009544 i 0.06179431 ± 0.00009544 i
13 0.06179520 ± 0.00009544 i 0.06179425 ± 0.00009544 i
14 0.06179803 ± 0.00009546 i 0.06179724 ± 0.00009545 i
15 0.06179807 ± 0.00009546 i 0.06179761 ± 0.00009546 i
90 0.16179997 ± 0.00065447 i 0.16179932 ± 0.00065446 i
91 0.16180003 ± 0.00065448 i 0.16179955 ± 0.00065447 i
92 0.16180020 ± 0.00065448 i 0.16179805 ± 0.00065444 i
93 0.16180029 ± 0.00065448 i 0.16179969 ± 0.00065447 i
94 0.16180039 ± 0.00065448 i 0.16179984 ± 0.00065447 i
95 0.16180047 ± 0.00065448 i 0.16179996 ± 0.00065447 i
96 0.16180062 ± 0.00065448 i 0.16179999 ± 0.00065447 i
97 0.16180063 ± 0.00065448 i 0.16180003 ± 0.00065448 i
98 0.16180070 ± 0.00065449 i 0.16180008 ± 0.00065448 i
99 0.16180070 ± 0.00065448 i 0.16180010 ± 0.00065448 i
100 0.16180076 ± 0.00065449 i 0.16180016 ± 0.00065448 i
135 3.21092421 ± 0.00012816 i 3.15927753 ± 0.00012834 i
136 3.21092421 ± 0.00012816 i 3.15939426 ± 0.00012844 i
137 3.39884491 ± 0.00012807 i 3.31429443 ± 0.00012829 i
138 3.77441335 ± 0.00012790 i 3.68814316 ± 0.00012808 i
139 3.77441335 ± 0.00012790 i 3.68894501 ± 0.00012809 i
140 4.71172216 ± 0.00012751 i 4.57325706 ± 0.00012775 i
141 4.71172217 ± 0.00012751 i 4.57372169 ± 0.00012786 i
142 4.89890560 ± 0.00012743 i 4.79083290 ± 0.00012797 i
143 4.89890560 ± 0.00012744 i 4.79105072 ± 0.00012793 i
144 5.45989728 ± 0.00012722 i 5.30883064 ± 0.00012762 i
145 5.45989728 ± 0.00012722 i 5.30934525 ± 0.00012751 i
146 6.02005043 ± 0.00012701 i 5.80780907 ± 0.00012723 i
147 6.39302010 ± 0.00012688 i 6.17382736 ± 0.00012755 i
148 6.39302011 ± 0.00012688 i 6.17507515 ± 0.00012727 i
149 6.95177672 ± 0.00012667 i 6.74670029 ± 0.00012756 i
150 6.95177673 ± 0.00012668 i 6.74753723 ± 0.00012768 i
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Russillo, A.F.; Failla, G.; Amendola, A.; Luciano, R. On the Free Vibrations of Non-Classically Damped Locally Resonant Metamaterial Plates. Nanomaterials 2022, 12, 541. https://doi.org/10.3390/nano12030541

AMA Style

Russillo AF, Failla G, Amendola A, Luciano R. On the Free Vibrations of Non-Classically Damped Locally Resonant Metamaterial Plates. Nanomaterials. 2022; 12(3):541. https://doi.org/10.3390/nano12030541

Chicago/Turabian Style

Russillo, Andrea Francesco, Giuseppe Failla, Ada Amendola, and Raimondo Luciano. 2022. "On the Free Vibrations of Non-Classically Damped Locally Resonant Metamaterial Plates" Nanomaterials 12, no. 3: 541. https://doi.org/10.3390/nano12030541

APA Style

Russillo, A. F., Failla, G., Amendola, A., & Luciano, R. (2022). On the Free Vibrations of Non-Classically Damped Locally Resonant Metamaterial Plates. Nanomaterials, 12(3), 541. https://doi.org/10.3390/nano12030541

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop