1. Introduction
Lightweight technology has emerged as an undisputed trend topic in recent years [
1,
2,
3,
4,
5]. Due to the traditional mass law, lightweight materials and designs have a major drawback: they worsen vibrational and acoustic performance. Conventional solutions to suppress unwanted mechanical vibrations and noise often require heavy and bulky constructions, contradicting the recent trend toward the introduction of lightweight materials and designs [
6,
7].
The concept of vibroacoustic metamaterial based on the local resonance effect has shown a promising solution to achieve lightweight and good vibroacoustic performance [
8,
9,
10]. The working principle of VAMM is mainly based on their microstructure, which consists of repeated elements called unit cells (UCs). Each UC is composed of a spring-mass-damper system called a local resonator. The resonators could either be attached to the host structure or integrated into the host structure. In the vicinity of the resonator’s resonance frequency, VAMM forms a so-called stop band. A stop band is a frequency range in which no free propagation of elastic waves is possible [
11,
12,
13]. Only undamped systems have perfect stop bands. All real-world structures possess a certain degree of damping. In real damped systems, a stop band is, therefore, typically characterized as a frequency band where vibrations are significantly reduced. The interaction between the local resonators and the traveling wave leads to the formation of a stop band [
14]. To achieve a stop band, the distances between the resonant elements must be smaller than half the wavelengths of the waves to be affected [
13,
15].
Since VAMMs usually have periodic configurations, wave propagation techniques such as Bloch waves, Floquet models, and transfer matrix methods are commonly used to derive the dispersion relations from the periodic symmetry of the metamaterial UCs to generate the wave dispersion curves [
16,
17,
18,
19]. With the help of dispersion curves, the stop bands of VAMMs could be predicted. Even though VAMMs do not require periodicity, they are often implemented for manufacturing reasons or modeling purposes [
19,
20,
21]. The prediction of the stop band behavior with the methods mentioned earlier is purely based on the dynamic behavior of the UC model, representing infinite periodic structures, and not taking into consideration the real boundary conditions [
19,
22]. Additionally, the UC model may not take into account the design complexities or non-uniform variations in material properties or geometries within the periodic structure. Consequently, the responses of actual metamaterials often deviate from the predictions of the UC dispersion curves [
19]. This deviation can result in the VAMM not performing as expected and can limit the effectiveness and applicability of the material. In addition, the Bloch approach becomes complex when damping is included, leading to complex frequencies and wavenumbers in the dispersion analysis [
23,
24]. It is, therefore, of enormous interest to researchers to investigate the dynamic behaviors of finite VAMMs using other approaches.
As the performance of VAMM depends primarily on the properties of the resonators, any uncertainty or variability in the resonator parameters, such as mass, spring stiffness, and the damping ratio, has an impact on the vibroacoustic performance of the VAMM [
25,
26]. In the context of industrial applications, the material and manufacturing tolerances can also introduce disturbances that affect the periodicity [
27,
28]. In [
29,
30,
31], the influence of uncertainties in periodic structures has been investigated. In order to put the VAMM concept into real-world applications, it is important to quantify the impact of uncertainties in factors such as resonator parameters, the spatial distance between resonators, real-boundary conditions, and non-periodicity on the stop band of the finite VAMM. Studying the impact of these factors on the finite VAMM performance is crucial in identifying optimal design and manufacturing strategies of finite VAMM that can minimize the impact of parameter variability and ensure consistent performance.
The stopband width of a metamaterial can serve as an important performance indicator, particularly when examining the influences of parameter uncertainty. Compared to other approaches mentioned above, the analysis of the frequency response function (FRF) emerges as the ideal method for evaluating uncertainties, as it presents fewer limitations compared to other techniques. However, predicting the width of the stop band solely based on the FRF analysis poses a significant challenge. One potential approach is to compare the FRFs of the VAMM with its host structure. By examining the intersection between these two curves in the vicinity of the resonance frequency, the upper and lower limits of the stop band can be defined (see
Figure 1a). Another possible criterion for determining the stop band could be a specific reduction in decibels (dB) value (e.g., a 20 dB reduction) of the VAMM FRF in relation to that of the host structure (see
Figure 1a).
However, identifying the start and end points of the stop band can be challenging due to the sensitivity of this method to the FRF curves. Even slight modifications of the resonator parameters can lead to significant variations in the width of the stop band. For instance,
Figure 1b displays the FRFs of the host structure and the VAMM when the resonator’s mass is reduced by 21%. According to the theory, a decrease in the mass of the resonators should result in a smaller stop band. However, upon examining the intersection points in
Figure 1b, it is observed that the predicted stop band using this method is wider, contradicting the theoretical expectation. Moreover, when considering a reduction of 20 dB, two separate stop bands can be identified instead of a single, smaller continuous stop band. Introducing uncertainties into the system further complicates the analysis. Therefore, a robust approach is necessary to overcome these limitations and accurately determine the boundaries of the stop band.
In the context of developing an automated stop band identification process for VAMMs with the presence of uncertainties, it is imperative to clearly outline the definition of the stop band. The automated determination of stop bands becomes vital for the optimization and evaluation of VAMM performance when considering parameter uncertainties. Therefore, this overall research aims to provide an effective and robust algorithm for accurately identifying the stop band boundaries, thus enabling automated processes for uncertainty assessment. One potentially valuable tool for investigating vibroacoustic problems is the utilization of active structural intensity (STI) [
32].
Recently, the utilization of a power-based approach that relies on an active STI in the scope of VAMM has demonstrated promising results in evaluating the existence of a stop band. H. Al Ba’ba’a et al. presented a numerical investigation of vibrational power flow in one-dimensional (1D) metamaterials, demonstrating that the STI analysis can accurately predict power flow and stop band frequencies, providing an alternative to traditional wave dispersion analysis and negative effective mass methods, while also enabling an accurate depiction of the performance of dissipative metamaterials and facilitating the optimization of actual metamaterial designs [
19]. In addition, an active STI can also be used to predict the occurrence of stop bands in 1D dissipative phononic structures with viscoelastic components, providing a clear distinction between Bragg scattering effects and wave attenuation due to material damping [
24]. Most recently, the STI approach has been extended to a two-dimensional (2D) metamaterial plate, presenting a detailed assessment of the experimental measurements of the different STI representations and comparing the experimentally determined power flow pattern with the mechanical energy amplitude to find the metamaterial plate’s stop band, which was validated against the dispersion curve predictions [
33].
Despite the limited number of published papers on the utilization of the STI approach in the field of VAMM, the above analyses have shed light on its application. These analyses include numerical investigations into the STI approach in 1D locally resonant metamaterials and 1D dissipative phononic structures, as well as experimental evaluations of the STI method in 2D locally resonant elastic metamaterials. However, as far as what is currently known, there has been no numerical study on the utilization of the STI method in 2D VAMM plates. Furthermore, the current state of research lacks a precise definition of the stop band limits for finite VAMM, which poses a challenge in conducting the uncertainty analysis and automating the detection of stop bands. This paper aims to address these gaps by providing a clear definition of the stop band in terms of a percentage, thereby facilitating the automation and accurate localization of stop bands in finite VAMM.
This work is mainly divided into two parts. The first part of this work utilizes the active STI to develop a new power-based approach that quantifies the power loss in a 2D VAMM plate using the finite element method. In this work, a clear definition of the stop band in terms of percentage has been established, which can be leveraged for a variety of applications based on the intended purpose and desired degree of vibration attenuation of the VAMM structures. In order to facilitate the discussion and validate the results, the stop band predicted using the power-based method is compared with the unit cell dispersion analysis method (UCDA), the negative effective mass method (), and the FRF.
In the second part of this work, the developed power-based approach is then applied in uncertainty analysis to create a robust algorithm that accurately evaluates and determines the impact of parameter uncertainties on the stop band behavior of a finite VAMM plate. The current study concentrates solely on the mass parameter uncertainty to demonstrate the applicability and robustness of the developed algorithm for uncertainty analysis, with the intention of extending its use to examine the effects of other parameter uncertainties on stop band behavior.
2. Structural Intensity
STI describes the magnitude and direction of the power flow in a structure excited by an external force. STI is defined as the vibration power flow per unit area transmitted in an elastic structure subjected to a dynamic loading [
33,
34]. Similar to airborne sound intensity, STI is a vector resulting from the product of the stress tensor,
S, and the velocity vector,
, as shown in Equation (
1) [
35].
In other words, STI is the work conducted by a unit area in the direction of velocity [
36]. This means that STI is a vector quantity that indicates not only the flow of energy across a unit area over time but also the direction in which the energy is transported. The negative sign in Equation (
1) comes from the fact that
and
are tensile stresses, which are distinguished from compressive forces by the sign [
37].
The STI in Equation (
1) shows energy transport as a function of time and is referred to as instantaneous STI [
38]. This means that the STI varies depending on its location in the elastic medium and over time. However, when dealing with harmonic structure-borne sound waves, the structural intensity in the frequency domain,
, is of interest and it could be derived from time-averaged
(e.g., over one period) [
39]. The STI at a steady state can be calculated in the frequency domain using a Fourier transform in Equation (
1) [
33,
40], which yields Equation (
2), where
is a frequency-dependent stress tensor and
is a complex conjugate of the frequency-dependent velocity.
The underline denotes complex quantities and the asterisk (*) conjugates complex quantities. The subscript letter,
s in
, refers to steady-state STI. The factor 1/2 in Equation (
2) is derived from the time-averaging process of the function. Thus, the STI is the product of the peak values of the stresses and the vibrational velocities. According to Equation (
2), the complex STI is composed of a real component known as active STI,
(Equation (
3)), and an imaginary component known as reactive STI or reactive power,
(Equation (4)) [
39].
The active STI,
, describes the energy flow from the source to the sink (moving wave) on a time-average basis [
39]. On the other hand, the reactive STI,
, describes the amount of energy that constantly oscillates in a structure (standing wave), allowing conclusions to be derived about the amplitude distribution of the natural mode of vibration. This presented work utilizes the active STI to quantify the vibrational energy flow and energy transmission pathways within a locally resonant finite VAMM structure.
3. STI Calculation in Thin Structures
For thin-walled structures, such as thin plates, it can be assumed that the energy transport across the plate thickness is negligible
[
41]. As discussed by Noiseux [
42], the STI calculation can be simplified from 3D to 2D shell structures by performing integration over the thickness of the shell, thereby reducing it into
x and
y components only. The integrated intensity is now equal to the net energy flow per unit width (length) of the mid-surface of the shell with the unit of W/m, as shown by Equation (
5) [
43], where
,
,
,
,
,
,
, and
denote bending moments, torsional moments, in-plane shear forces, out-of-plane shear forces, and in-plane axial forces per unit width of the plate, respectively. Also,
,
,
, and
,
represent translational velocities and rotational velocities. As a result, it is possible to specify the STI integrated over the shell thickness as a function of the internal forces and moments and, therefore, easily compute using FEM simulations [
44].
The first two summands in the equation describe the energy transport of the in-plane waves by the membrane forces
,
, and
. The other summands represent the energy transport of the out-of-plane waves caused by the shear forces (
and
), the bending moments (
and
), and the torsional moment (
). These so-called in-plane waves consist of transverse and extensional waves. The quantities from Equation (
5) are defined positively according to
Figure 2.
4. Energy Balance in Elastic Medium
STI in an elastic medium can be described using the first law of thermodynamics [
41]. The law of conservation of energy is represented by Equation (
6) [
46]. The law states that the total energy of a closed system is constant, implying that the sum of all energy changes in a closed system is zero.
where
e is the energy density inside the control volume,
is the displacement vector of any particle at the boundary of the control volume,
is the input power density (or the energy supplied per unit volume per unit time),
is the dissipated power density (or the energy dissipated per unit volume per unit time),
is the vector perpendicular (normal) to the surface of the control volume for a given point on the surface, and
S is the stress tensor.
Although the law was developed for thermodynamic processes, it can be used to represent the energy balance in an elastic medium using a control volume approach. The rate of change of total energy within the surface enclosing the volume is equivalent to the energy flow through a closed surface [
46]. The first summand on the left side of Equation (
6) describes the change in energy density within the control volume. This change in energy density is due to the work conducted by the stress on the medium. In terms of stress and velocity, the local outflow of the energy density from the control volume is known as STI, as described in Equation (
7) [
46].
By substituting Equation (
7) into Equation (
6), the equation can be rewritten as
where
is the input power (in-flow of the energy into the control volume) and
is the dissipated power density (out-flow of the energy from the control volume). Equation (
8) represents the energy relationship for all elastic media and is valid for steady-state or transient analysis. When dealing with harmonic structure-borne sound waves, the change in energy over time,
, is zero for steady-state (time-averaged value) vibrational energy propagation and, hence, Equation (
8) becomes
where the angled brackets
reflect the time-averaging. The correlation between time-averaged STI and the consideration of STI in the frequency domain is as follows:
8. Power-Based Method
As this work primarily focuses on waves propagating in the x-direction only, the power flow can be summed across the width of the plate, transforming the 2D plate model into a 1D beam-like model, revealing a similar power flow pattern to a 1D system, like a beam.
Figure 12 shows the power flow along the length of the plate for a frequency of 131 Hz, which is within the desired stop band frequency range. The vertical spike in power corresponds to the input power contributed by harmonic force excitation,
, and the decline of power across the length of the bare plate is due to the material loss factor of the plate. The power flow of the VAMM plate in
Figure 12 also illustrates a step-like decrease in power, which suggests that the energy is being dissipated at every column of resonators. It could be seen that in the region before the resonators and after the force excitation point (between
m and
m), the power in the VAMM plate is higher than the bare plate. This accumulation of power could be understood as a compensatory effect. Such an increase compensates for the local resonance band gap effect, which blocks the transmission of energy beyond the resonator region [
33].
Figure 13a,b denote the graphical representations of power flow along the length of the VAMM plate, calculated for frequencies ranging from 0 to 200 Hz. The power,
, is normalized to the injected force,
. In
Figure 13b, the horizontal dark blue region starting at 0.9 m indicates that the wave is effectively blocked as the energy is absorbed by the resonators. The frequency range of this horizontal dark blue region is consistent with the theoretical stop band and the drop in the FRF graph.
The power loss in the resonator region can be determined by using Equation (
11) by defining the contour line around the region with resonators, as illustrated in
Figure 14. The vectors normal to the closed area (top, bottom, left, and right) are represented by
,
,
, and
, respectively. Since no power can flow beyond the clamped and free edges [
41,
51], the injected power may only flow over the left or right vertical lines.
Therefore, Equation (
11) can be rewritten as follows:
where
, and
denote the power flow into and out of the resonator region, respectively. Therefore, by substituting the values of
and
into Equation (
23), the power loss basically equals the difference between power flowing into and out of the resonator region, (
=
−
). To express the power loss in terms of percentage, the power inflow and outflow,
and
, are normalized with respect to the power entering the resonator zone,
.
In order to quantify power loss in the VAMM plate, relative to the host structure, the power loss calculation is repeated for the host structure (i.e., bare plate) in a similar manner, and then the power flows in the VAMM plate are normalized with respect to the power flows in the host structure,
and
, respectively (Equation (
24)), to yield the relative power loss caused by the resonators.
The STI method only yields meaningful results in the presence of damping components, such as internal material damping or dissipative resonators. In the absence of any damping, the active STI,
, becomes zero. For validation purposes, a so-called undamped system is achieved by setting the material damping of the VAMM plate and Lehr’s damping ratio of the resonator to very low values instead of zeros (
).
Figure 15a,b shows the power loss in the VAMM plate,
(green line), the host structure,
(blue line), and the relative power loss,
(orange line) for undamped and damped VAMM plates. It could be seen that approximately 40% of the power (blue line) is dissipated within the resonator region of the host plate due to material damping alone. Therefore, the normalization process is important to yield the relative power loss,
, induced by the resonators only.
Quantitative analysis of the vibration attenuation capabilities of a VAMM structure requires establishing a clear definition of the stop band in percentage terms. This step is crucial to accurately evaluate the structure’s vibroacoustic performance. The percentage value denotes the magnitude of the vibration attenuation achieved. Two threshold percentages were chosen in this work, namely, 90 and 99%, as shown in
Figure 15c. The stop bands predicted based on these threshold values are compared with the FRF graphs, as illustrated in
Figure 15d. The stop band gaps align well with the dip in the FRF of the VAMM plate. A relative power loss of 99% indicates that the resonators in VAMM structures absorbed 99% of the vibrational energy relative to the host plate. This definition of the stop band can be applied to various applications, depending on the intended purpose and the desired degree of vibration attenuation for the VAMM structures.
10. Parameter Uncertainty Analysis
The uncertainty quantification in this work is based on a one-factor-at-a-time (OAT) method, which means that only one parameter, uncertainty, is investigated. Here, the uncertainty in the resonator mass and its influence on the stopband width is of interest.
Figure 17 provides an overview of the simulation flow using a power-based approach in uncertainty analysis.
In order to create uncertainties in the model parameter, a probability distribution function is required, from which the samples can be drawn. The probability density function of the uncertain parameter depends on the source of its uncertainty. Experimental characterization of uncertainty is needed to appropriately define the distribution function. In the absence of such detailed characterization, a normal or Gaussian distribution is considered to characterize the parameter uncertainty [
52].
The uncertainty in the parameter is described by a normal distribution function and quantified by the mean, , of the distribution, which denotes the nominal value of the uncertain parameter, and its standard deviation, , which is proportional to the uncertainty of the parameter. The uncertainty level can be expressed by the coefficient of variation (CV), which is also known as the normalized standard deviation. The CV is defined as the ratio of the standard deviation to the mean, .
In order to accurately represent the output distributions, in this case, the stopband width corresponding to the uncertainty of the input parameter, the sampling technique, and the number of samples are crucial for the accuracy of the uncertainty analysis. Instead of using the random sampling method used in a classical Monte Carlo simulation, the Latin hypercube sampling (LHS) method is chosen, which is based on stratified sampling. Compared to pure random sampling, the LHS method eases the computational burden by sampling more effectively, thereby reducing the number of samples needed [
52].
Based on LHS, masses for 25 resonators will be selected from the normal distribution of the mass described by its mean and standard deviation, N (, ). In this study, the maximum CV value of 0.5 is chosen in order to evaluate the impact of larger uncertainty. During sampling, if the chosen mass has a value less than 0 g, the mass of the resonator is set to 0 g. The purpose of this very high maximum standard deviation is to find a general correlation between the variability of mass and the width of the stop band. It may not necessarily represent real manufacturing uncertainties.
During the Monte Carlo simulations, the CV is increased at a step size of 0.01, dividing the CV into 50 intervals . To achieve an accurate output probability distribution, it is necessary to have a considerably high number of samples for each CV value. Taking into consideration the computational burden of performing the STI analysis, 40 simulations are performed for each CV value. This means that the STI of the VAMM plate is computed 40 times for each CV before increasing the CV by 0.01. For each simulation, the probability distribution of the mass of the sample is calculated. The probability distribution follows a normal distribution and, therefore, has a slightly different mean and standard deviation compared to the distribution function from which these samples are drawn. The probability density function of the 25-resonator mass is of particular interest since it demonstrates the actual variation among them.
The power loss, , is computed, and the stopband width is determined for the two threshold values, , and . The stopband widths for all simulations are plotted as functions of mass uncertainty, expressed as the CV of the input parameter, . In the case of an ideal damped finite VAMM plate where there is no variation in the mass of the resonators, a single stop band is present. However, with increasing uncertainty, there is the possibility that two or more stop bands will be identified. Since the influence of uncertainty is evaluated in terms of stopband width, the stopband widths of two or more stop bands are added to obtain the total stopband width. The larger deviation of the mass of the resonators affects the natural frequencies of the individual resonators and the mass ratio of the resonators to the host structure. This leads to the formation of more than one stop band.
Figure 18 shows the stopband width as a function of CV for two different threshold values (a)
and (b)
. Inconsistencies in the mass of the resonators shift the resonance frequency of individual resonators, causing all resonators not to work on the same frequency anymore. This widens the stop band at the expense of the stop band depth. It is noticeable that the maximum CV for this simulation is above 0.5. Although, the maximum specified CV for the normal distribution is 0.5, the actual CV of each sample for the 25 resonators may differ from this value. The higher the specified CV, the greater the likelihood that the sample’s CV will deviate greatly from this value. The 95% confidence interval (CI) is indicated by the two blue lines, which signifies that 95% of the time, the stopband width will lie within this range. It is useful to examine the fluctuation of the mean stopband width values in relation to their corresponding standard deviations, in order to estimate the maximum achievable mean stopband width given an acceptable maximum variation in resonator masses.
In
Figure 18a, the mean width of the stop band increases and reaches a peak at a CV value of 0.3 before decreasing again. However, the mean width always remains higher than the ideal stopband width (27 Hz). Similarly, in
Figure 18b, the mean width of the stop band reaches a maximum CV value of 0.1 before decreasing below the ideal stopband width (18 Hz), eventually nearing zero. This indicates that the VAMM performance has significantly decreased, and it cannot achieve the desired stop band with
power loss. Moreover, the comparison of the confidence intervals in
Figure 18a,b reveals that the CI
in
Figure 18a is generally wider than that in
Figure 18b.
By investigating the variation of mean stopband width values with respect to their corresponding standard deviations, the maximum allowable variation in resonator masses can be estimated while still achieving the required mean stopband width. This can be useful in the design and optimization of VAMM structures to achieve the desired degree of vibration attenuation, where the stopband width is an important performance parameter.
Based on the
method (
Figure 18b), for the VAMM structure where great vibration reduction is the goal, the tolerances on the mass fabrication process must be tightened to ensure that the resonator masses are more consistent, thereby improving the reliability of the VAMM performance. On the other hand, the
method (
Figure 18a) indicates that larger fabrication tolerances could be positively exploited to achieve a wider stop band range while still fulfilling the desired performance level at lower manufacturing costs.
11. Conclusions
This study presents a novel methodology for quantifying the power flow in a VAMM plate using the active STI integrated with FEM. The study was conducted on a rectangular Kirchhoff–Love plate model made of steel and excited by a point force. The resonators attached to the middle section of the plate were tuned to create a stop band in the frequency range of around 130 Hz. By utilizing the active STI, the relative power loss caused by the resonators in comparison to the host structure was calculated, predicting the stop band of the VAMM plate. Validation of the STI results with other approaches, such as the negative effective mass method, the UCDA method, and the frequency response function method for both undamped and damped systems, has shown that the STI method is a valuable tool that can be used to reliably predict the stop band of a finite VAMM structure. Even though this power-based approach is computationally intensive, it is a powerful and reliable tool used for predicting the stop band of real, complex 2D structures, considering damping, real boundary conditions, non-uniformities in the structure, and not depending on periodicity. This work also establishes the definition of a stop band in terms of percentage, which can be leveraged for a variety of applications based on the intended purpose and desired degree of vibration attenuation of the VAMM structure.
In the second part of this work, a robust algorithm based on relative power loss in VAMM was presented to evaluate and determine the impact of parameter uncertainties on the stop band behavior of a finite VAMM. Based on the method, an increase in mass uncertainty up to a CV of 0.1 has a positive effect on the stopband width. However, if the mass uncertainty increases beyond this value, the performance of VAMM deteriorates significantly. On the other hand, the method demonstrates that the variability in the mass uncertainty can be utilized in widening the width of the stop band, where the focus is on reducing vibrations by 90% compared to the host structure. In short, the relationship between the parameter uncertainties and the stopband width can be employed in the design and optimization process of VAMM, such as in determining the maximum allowable uncertainty while still achieving the desired level of vibration reduction. In order to study the impacts of other parameter uncertainties, such as spring stiffness, damping ratio, and the position of resonators on the stop band behavior, this proposed uncertainty quantification algorithm can be easily implemented. The findings can be useful in industrial applications for determining whether to loosen or tighten the manufacturing tolerances based on the sensitivity of the stopband width to parameter variations. In the future, the developed framework will be implemented on a more complex VAMM structure.