Research Article  Open Access
Vibroacoustical Analysis of MultipleLayered Structures with Viscoelastic Damping Cores
Abstract
This paper presents a modeling technique to study the vibroacoustics of multiplelayered viscoelastic laminated beams using the Biot damping model. In this work, a complete simulation procedure for studying the structural acoustics of the system using a hybrid numerical model is presented. The boundary element method (BEM) was used to model the acoustical cavity, whereas the finite element method (FEM) was the basis for vibration analysis of the multiplelayered beam structure. Through the proposed procedure, the analysis can easily be extended to another complex geometry with arbitrary boundary conditions. The nonlinear behavior of viscoelastic damping materials was represented by the Biot damping model taking into account the effects of frequency, temperature, and different damping materials for individual layers. The curvefitting procedure used to obtain the Biot constants for different damping materials for each temperature is explained. The results from structural vibration analysis for selected beams agree with published closedform results, and results for the radiated noise for a sample beam structure obtained using a commercial BEM software are compared with the acoustical results of the same beam by using the Biot damping model.
1. Introduction
The traditional designs of freelayer, constrainedlayer or sandwichlayer, damping treatment using viscoelastic materials have been around for over forty years. Recent improvements in the understanding and application of the damping principles, together with advances in materials science and manufacturing, have led to many successful applications and the development of patch damping and multiplelayered damping structures. The key point in any design is to recognize that the damping material must be applied in such a way that it is significantly strained whenever the structure is deformed in the vibration mode under investigation. Numerous researchers have successfully implemented the passive constrained layer (PCL) and active constrained layer (ACL) systems. In 1959, Kerwin [1] and Ross et al. [2] presented a general analysis of viscoelastic material structure. The damping was attributed to the extension and shear deformations of the viscoelastic layers. Ditaranto [3] developed sixthorder equations of motion in terms of axial displacements and developed a closedform solution. Mead and Markus [4] extended the sixthorder equations of motion for transverse displacement to include various boundary conditions. A paper by Rao [5] presented the equations of motion of viscoelastic sandwich beams with various boundary conditions using the energy method. The equations were solved numerically, and a practical design guideline was presented. Similar to Rao’s theory, Cottle [6] used Hamilton’s principle to derive equations of motion. The damping could also be increased by adding passive standoff layer (PSOL) and slotted standoff layer (SSOL) to the layered systems. Falugi [7] and Parin et al. [8] conducted theoretical and experimental work on a fourlayered panel and a fivelayered beam with PSOL treatment. Rogers and Parin [9] and Yellin et al. [10] have performed experimental investigations and demonstrated that PSOL treatment increased the damping significantly in aeronautical structures and beams. Yellin and colleagues [11, 12] also developed normalized equations of motion for beam, fully treated with PSOL using nonideal standoff layer assumption. The equations were solved using the method of distributed transfer functions [13].
In addition to the closedform analytical approach, many researchers have used the finite element method (FEM), the most popular numerical modeling method in building the numerical model of the multiple layers system. In 2000, Chen and Chan [14] studied four different types of integral FEM models with the viscoelastic cores. The numerical stability and accuracy as well as for convergence issue of these four different FEM models were demonstrated by comparing the numerical results with those from experiments. Lesieutre and Lee [15] proposed a 3node, 10 DOF FEM model for the threelayer ACL damping beam. This FEM model is advantageous in active control application due to its features of nonshear locking and adaptability to segmented constraining layers.
Other than the Hamilton and FEM numerical methods in building the models, other researchers have proposed many irregular modeling techniques for the numerical representation of continuous/discontinuous systems. Kung and Singh [16] calculated the natural frequencies and loss factor using the RayleighRitz energy method and modalstrain energy technique in modeling a 3layer patch damping structure. These approximate modeling methods were also extended to rectangular damping patch of plates and shells with viscoelastic cores. Zhang and Sainsbury [17] combined the Gerlerkin orthogonal function with the traditional finite element method and successfully applied to the vibration analysis of the damped sandwich plates.
While the FEM is used widely in the modeling of the structure, many researchers sought for proper mathematical models to represent the damping behavior of the viscoelastic material, as well as incorporating the damping model in commercial FEM software packages. Currently, many FEM commercial software incorporate damping models based mostly on viscous/hysteretic damping. Some allow incorporation of damping energy dissipation in the time domain using the Prony series. None of these damping models, however, is suitable to capture the damping behavior in the frequency domain, which is the most important issue in predicting the vibroacoustical response of complex structural systems. The drawback of these damping models raised considerable interest and motivation in the development of damping models of viscoelastic material in the frequency domain compatible with FEM software. These damping models can be classified as derivative type and integral type.
The “Fractional Derivative” is essentially the representative damping model in the derivative form family proposed by Bagley and Torvik [18] in 1983. This damping model not only described the material properties of viscoelastic damping but so established the closedform equation compatible with the FEM technique. Compared with the other integralform models, the fractional derivative is only able to capture the relatively weak frequencydependent information; however, it was an important milestone in the area of damping research.
Lesieutre et al. [20] mathematically modeled the relaxation behavior of viscoelastic material in terms of augmenting thermodynamic field (ATF) in 1989. Initially, introducing a single augment field, this damping model provided the ability to represent the lightdamping behavior, with the application of a 1D viscoelastic structure. In the subsequent research, using a series of augment fields, the ATF model is able to model the damping material of higher loss factor with the weak frequency dependence. Remedying the limitation of 1D application, Lesieutre and Lee [15] proposed an anelastic displacement field (ADF) technique in 1996 and successfully extended its application from the 1D problem to the 3layer sandwich beam and 3D problems.
As far as the minioscillator damping models are concerned, the complex shear modulus which is a function of both frequency and temperature can be expressed by a series of minioscillation perturbations. Biot [21] first proposed the firstorder relaxation function with the introduction of the “dissipative variables” into the dynamic equations using the theory of irreversible thermodynamics. In 2007, Zhang and Zheng [22] utilized the Biot model to describe the dynamic behavior of a viscoelastic structure. The dimension reduction technique and nonlinear curvefitting procedure were discussed in the paper. McTavish [23] developed another minioscillator damping model called “GHM” by the usage of Secondorder relaxation function. Compared with the Biot model, the GHM model has a more complicated expression and also requires better performance of the computational tool.
The popularity of these integralform damping models in recent years brought two research interests: nonlinear curve fitting and dimension reduction. The advanced curvefitting techniques in the damping models guarantee the accuracy of the numerical representation of the actual shear modulus data from the experiment. The dimension reduction technique increases the computational efficiency due to the additional orders of equation in order to gain the frequency independence of the frequencyform damping model.
Zhang et al. [24] converted the nonlinear curvefitting problem in frequency domain with respect to the GHM parameters into the constrained nonlinear optimization problem. The efficiency and correctness were demonstrated for a commercial viscoelastic material.
Park et al. [25] examined the GHM damping model with the application to the FEM method associated with the Guyan reduction technique. The numerical example in this research leads to an FEM model applied to the GHM dynamic equation quantitatively without increasing the number of order.
Hao and Rao [26] carried out the optimum design of a threelayer sandwich beam for the vibration analysis in 2005. In this research, the numerical model is a comprehensive formulation for a threelayer unsymmetrical sandwich beam with two different damping materials adjacent to each other. The criterion of the optimization is to minimize the mass of the structure while maximizing the system damping. In 2008, Lee [27] published the semicoupled vibroacoustical analysis and optimization of a simply supported threelayer sandwich beam. The modal superposition method was used to investigate the vibration problem with the fractional derivative damping model. The interior acoustical problem was studied by BEM numerical technique, and the optimization problem was established through the appropriate sizing parameters of the sandwich beam.
The objective of this paper is to extend the previous work by the authors [28] on the vibration analysis of a multiplelayered beam structure incorporating the Biot damping model to solve the acoustic problem to predict the radiated noise. In this paper, we present a complete numerical procedure for the vibroacoustical analysis and design for a multiplelayer laminated damping beam. Results obtained from the proposed vibration analysis are compared with the previous closedform results to show the validity of this approach. The radiated noise spectrum at selected field point shows good agreement between the 2D BEM acoustical analysis and the result without system damping calculated by commercial software for a sample viscoelastic damping structure. The acoustical solution is demonstrated, and the correlation between sound pressure level (SPL) and the loss factor is also highlighted.
2. FEM Modeling and the Biot Dynamic Equation
The FEM modeling procedure and the establishment of the Biot dynamic equation will be discussed in this section. The structure chosen for illustration is a sevenlayer viscoelastic sandwich beam. The elastic beam and the constrained damping layer are the two fundamental components in this FEMmodeling technique. The concept of transfer matrix is used to convert the local coordinates to the global coordinates in order to assemble and construct the complete model of the sandwich damping structure with arbitrary number of layers. The Biot viscoelastic damping model will be used to describe the damping behavior. Through the use of the FEM, the structure is discretized which will enable the use of the Biot damping model for different damping layers in the structure. The reader is referred to the nomenclature for the definition of different variables used in the derivation.
2.1. FEM Modeling of Component I: The Elastic Layer
Figure 1 shows the elastic layer in the FEM model, containing 2 nodes and 6 degrees of freedom (DOF). The element displacements of each node can be expressed as follows:
The stiffness matrix can be derived based on the following energy method: as the shape functions are the following: in which : the local coordinate, , , : longitudinal length of elastic layer, : crosssectional area of the elastic layer, : Young’s modulus of the elastic layer, and : moment of inertia of elastic layer.
Similarly, the element mass matrix can be expressed as:
2.2. FEM Modeling of Fundamental Component II: The Constrained Damping Layer
The FEM model of the constrained layout containing the damping layer sandwiched between two outer layers is shown in Figure 2. This Figure illustrates each element consisting of 2 nodes and 8 DOF, where the nodal displacement vector is as follows:
Through the introduction of transfer matrix, in which each means the following vector: ; the element elastic stiffness and the element viscoelastic stiffness matrix for this 3layer component, respectively, are the following: where : crosssectional area of the damping layer, : longterm shear modulus of the damping layer, and : correction factor of the shear strain energy, for the rectangular crosssection, .
Also, the element mass matrix for this 3layer component is where
2.3. FEM Modeling of a SevenLayer Constrained Damping Beam
The sevenlayer sandwich beam consists of seven alternating layers—four elastic layers and three damping layers. Figure 3 shows the FEM model of a sevenlayer sandwich beam containing 2 nodes and 10 DOF, and the node displacement vector is as follows:
The transfer matrix to obtain the element stiffness and the mass matrix when the 1st, 3rd, 5th, and 7th layers are elastic are follows:
Similarly, the element stiffness and the mass matrix for the 2nd, 4th, and 6th layers of the constrained damping layer can be derived through the transfer matrix: where the notation means
Based on the above equations and design parameters of each layer, the element mass/stiffness/damping matrix of the sevenlayer sandwich damping beam can be expressed as follows:
Thus, the element matrices can be assembled to obtain the global mass/stiffness/damping matrix and can be applied to the boundary condition through the conventional FEM technique. Taking into the consideration of the viscoelastic damping properties, the global matrices need to be manipulated as a portion of the Biot dynamic equation.
2.4. Introduction of the Biot Dynamic Equation
To consider the vibration problem numerically, the dynamic equation discretized by FEM technique needs to be expressed by the following secondorder ordinary differential equation (ODE) form:
The Biot viscoelastic damping model numerically represents the complex shear modulus with a series of minioscillator perturbing terms: in which is the longterm shear moduli; and are the Biot constants. These parameters are positive and can be determined by nonlinear curve fitting from the experimental data. The curvefitting procedure will be discussed in Section 3.
Substituting the Biot damping model into (15), the dynamic equation with terms of the Biot parameters for the first viscoelastic material and terms for the second viscoelastic material can be developed as follows: where and are the eigenvector and diagonal eigenvalue matrices, respectively, from the damping matrix . Additionally, , , and denote terms of the Biot parameters and the dissipative coordinates, respectively, for first viscoelastic material.
Similarly, , , and denote terms of the Biot parameters and the dissipation coordinates, respectively, for second viscoelastic material. A detailed derivation can be found in the previous publication [28].
3. Parametric Determination of the Biot Damping Model
A curvefitting technique is used to provide the accurate Biot constants to the dynamic equation and to establish the dynamic characteristics of the viscoelastic materials. In this section, the nonlinear curvefitting procedure for the complex shear modulus in the frequency domain is converted into a nonlinear constrained optimization problem.
The complex shear modulus with the Biot damping model form can be broken into real and imaginary parts separately:
The Biot parameters—, , and —are estimated from experimental data with the certain fitting frequency range, on real part and imaginary parts separately. Generally speaking, one set of the Biot parameters needs to be determined for each ambient temperature independently. In (18), is the number of the Biot perturbing items, defining the capability of this numerical approximation. As the Biot terms () are increased, the relative error between the experimental data and the curvefitting result reduces.
Assuming ; ; ; ; ; with the constraint condition: ; , the target equation of the optimization problem is the following:
In the target equation (19), stands for the complex shear modulus from the experimental data with interested points (larger than the number of unknowns). The 3 M ISD110/112 viscoelastic polymer is selected in this simulation. The experimental data is obtained by the Arrhenius empirical equation from [19]. With a specific fitting range at a particular temperature, the complex shear modulus can be synthesized from one set of the Arrhenius coefficients. The number of terms () in (18) needs to be determined to ensure the precision of this approximation. The curve fitting of the experimental data is accomplished using the commercial software package Auto2fit on the real and imaginary parts simultaneously. Using the Biot terms equal to six and four with respect to two commercial damping materials 3 M ISD110 and 112, respectively, the results are shown in Table 1 for ambient temperature () equal to 45°C and frequency range of 500 Hz.

Figures 4(a) and 4(b) show the comparison between the Arrhenius data and curvefitting data for the real and imaginary parts, respectively. Figure 5 shows the relative error in the fitting range.
(a)
(b)
As shown in Figures 4(a) and 4(b), the Biot parametric determination technique estimates the dynamic properties of 3 M ISD110/112 at 45°C with almost zero error. The constants determined using the above procedure along with the FEM model of sandwich beam will now be incorporated to solve the complete Biot dynamic equation using the decoupling transformation technique.
4. Decoupling Transformation and Dynamic Solution
In this section, the algorithm used to obtain the frequency response function (FRF) will be discussed with respect to the vibroacoustical problem for a multiplelayer viscoelastic damping structure. In this research, the damping matrix in (15) does not have a proportional relationship with the mass and stiffness matrix. Thus, a decoupling transformation is needed to construct the firstorder state equation by introducing the auxiliary equation as follows: where
Here, is the number of DOF in the , , and matrices, the DOF of and matrices is .
Firstly, the free vibration of (20) will be considered. Assuming , the following form of solution is obtained: or where matrix stands for complex conjugate eigenvalues including the natural frequencies and loss factors information:
It must be noted that zero items will appear in the eigenvalue matrix if the damping matrix is not fully ranked. The mode shape vector for the vector can be extracted from the eigenvector matrix with respect to the vector :
In addition, (22b) can be numerically solved by using mathematical software package such as MATLAB or Mathematica.
Secondly, the forced vibration solution of (20) in the time domain will be discussed. Assuming , the variable substitution can be made by assuming , converting the statespace equation from the time space to the modal space. By left multiplying of with the substitution of , we get:
The diagonal modal mass and stiffness matrix are:
Then rewrite the equation with the diagonal mass and stiffness matrices
The FRF in the frequency domain can be easily determined through the complex conjugate eigenvalue matrix , eigenvector matrix _{,} and the modal mass matrix . The modal scaling factor matrix can be calculated through the following:
Thus, FRF can be established through the modal parameters, being expressed in partial fraction form in terms of the residue vector and system poles as follows:
The system velocity can now be obtained from the above equation by a simple Fourier transformation. By doing so, the vibration problem can be extended to an FRFbased acoustical problem and the combination of these two analyses is the particle velocities information calculated by the following:
5. Acoustical Boundary Element Method (BEM) Analysis
5.1. Introduction of Acoustical BEM Theory
In Section 4, the vibration problem of the multiplelayer sandwich beam is solved through the timedomain dynamic ordinary differential equation of the Biot damping model with numerical analysis by the FEM technique. The vibration problem can be extended to the acoustical problem by the semicoupled method: the vibration will induce a change in sound pressure, yet the sound pressure will not cause the vibration. In this section, the acoustical interior problem will be numerically solved by 2D boundary element method (BEM) technique [29] in a bounded fluid domain as shown in Figure 6.
The sound pressure distribution () of the timeharmonic wave in the domain satisfies the governing differential equation, well known as the Helmholtz equation, associated with the boundary conditions on boundary as follows:
Here, is equal to , which means that the wave number is equal to the radiant frequency over the speed of sound; , , stand for the normal velocity, density of the fluid (normally the air), and acoustical impedance of the fluid , respectively.
In this work, the link between the vibration and the acoustics analysis is the normal velocity at the acoustical boundaries. Recalling the dynamic solution of the decoupling transformation, the particle velocity in the time domain at each node can be calculated through (30) if the multiplelayer sandwich beam is discretized by the FEM; alternatively the FRF, the complex ratio between the output and input response in the frequency domain, can be determined through (29). Once the input signal is given, the particle velocity of the system displacement versus frequency relationship can be conveniently obtained through the FRF.
To solve the governing differential equation (31) in the bounded fluid domain , the Helmholtz Equation can be transformed into the integral equation, converting the 2D area integration to the 1D curve integration around the area:
in which : geometrydependent coefficient, normally when is in the domain and when is on the smooth boundary , : sound pressure at source point , : is the field point and for the 2D BEM problem the Euclidian distance between and , and the Secondtype Henkel function, : normal vector pointing away to the fluid domain .
By discretizing the boundary into a series of curvelinear elements through the introduction of the shape functions, the integral equation can be calculated numerically by solving the following linear matrix: where comes from the terms of and , is derived from , and the vector and include sound pressure and particle velocity values, both unknowns and known from the boundary condition.
Thus, each set of node velocities due to the force input results in one set of solutions on the sound pressure by BEM discussed in this section. In sum, through the proposed acoustical BEM, it is possible to compute the timeharmonic sound pressure distribution corresponding to each single frequency point in the frequency spectrum.
5.2. Calculation Details in This BEM Analysis
For this particular acoustical BEM interior problem, the boundary of acoustical cavity is discretized as 18 quadratic equally spaced boundary elements. The quadratic curvilinear element has three nodes, and the interpolation between each node represents the geometry of each element. The shape functions are as folows: with respect to the following element coordinates: where and are the coordinates at each nodal point, and stands for the local coordinate between −1 and 1 on a master element.
When the sevenlayered sandwich beam ( m) is simply supported at the bottom of the acoustical cavity, the sound pressure level at the field point ( m, m) is calculated through this proposed method, and the calculation results are presented in Section 6. Figure 7 demonstrates the detailed layout of this 2D acoustical cavity problem. The anechoic boundary condition is applied on the inside of the acoustical cavity, and the thickness of the multiplelayered beam is neglected.
6. Numerical Results and Discussion
6.1. Design Parameter of Sandwich Beam and Vibration Analysis Result
In Figure 8, a sevenlayer sandwich beam with viscoelastic cores is shown, with the design parameters listed in Table 2.

The data presented in Table 2 are used to predict the vibration performance of the system using the numerical simulation method presented in this paper, and the results are compared with the closedform solution of Hao [19]. The curvefitting results for the damping material 3M ISD110 at 45°C discussed earlier are selected for the shear modulus of the viscoelastic layers in this example. The results are shown in Table 3. It shows that the simulation presented in this paper conforms to the closedform solution. This validates the analysis methodology proposed in the paper.

6.2. FrequencySpectrum Analysis under the Arbitrary Input
Figure 9 shows the transverse velocity of the middle node (node number 7) with a 10 N step input in the frequency domain vertically applied at the middle (node number 7) of the simplysupported sevenlayer sandwich beam with the same design parameters as the previous example. The same curvefitting results of 3 M ISD110 at the ambient temperature of 45°C for the shear modulus are used in this example. This pivotal result is the demonstration of extending the vibration to the acoustical problem in the frequency domain when an arbitrary force is applied on the structure.
6.3. Acoustical BEM Results
Figure 10 illustrates the contour plot ( Hz) of SPL when the sevenlayer sandwich beam (using the same design parameters as before) is subjected to a 10 N step input in the frequency domain at the middle node.
The interpolation of each elements result in Figures 10 and 11 shows the continuous sound pressure distribution in the acoustical cavity with an anechoic boundary condition. Figure 12 extracts the frequency spectrum of SPL at the filed point (0.5, 0.4 m) indicated by red dot in Figure 10. From the results of Figure 12, it can be found that the dominant contribution is due to the peak value of the first flexible vibration mode, which is in agreement with the frequencyspectrum analysis of the vibration problem.
6.4. Validation Using a BEM Commercial Software
In this section, a hybrid FEMBEM model of a beam without the viscoelastic damping was developed using the commercial software packages ANSYS ADPL and LMS Virtual. Lab Acoustics. The harmonic vibration analysis is conducted in ANSYS APDL module, and the frequency spectrum of field point SPL was calculated in Virtual.Lab Acoustics module for comparison with the SPL frequency spectrum presented in Section 5. The analysis sequence consists of the following steps.(a)Build the FEM model and apply appropriate boundary conditions in ANSYS ADPL. The 8node element SOLID45 (element size =10 mm for each layer) was used to build the 3D sevenlayer model. The design parameters are identical with the parameters in Tables 1 and 2 for the comparison and the geometry boundary conditions are simply supported. A 10 N force at each frequency is applied at the middle nodes.(b)Conduct the harmonic vibration analysis in ANSYS ADPL. The harmonic analysis is used to calculate the nodal displacements for a forced vibration problem in the frequency domain. The frequency range is 0–200 Hz with a 2 Hz for step size, and the full method is being utilized in this analysis. The comparison of system frequencies between ANSYS modal results and calculation results by the Biot dynamic equation is shown in Table 4. The results show that the 3D model built in ANSYS APDL has good correlation with the FEM model. (c)Prepare the BEM mesh in LMS Virtual. Lab PreAcoustics module. It converts from a solid FEM model to a skin mesh that the BEM analysis requires. The BEM mesh, can be seen as a wrap around the structural mesh and usually the BEM mesh is coarser.(d)Calculate the sound pressure in LMS Virtual.lab Acoustics module. Both acoustical and structural meshes are imported to VL Acoustics. The nodal displacement at each vibration mode calculated in ANSYS APDL is also imported and meshmapped to the acoustical skin mesh as the vibration boundary condition. The location of field plane and field point is consistent with the 2D BEM analysis in this research. The acoustical pressure is solved over the frequency range from 2 to 200 Hz.

As shown in Figure 13, the peak frequency from the 2D BEM calculation matches with the first dominant SPL peak obtained from the VL Acoustic result without the damping. Comparing the two results, it is clear that the introduction of viscoelastic damping not only causes almost a 20 dB reduction in the first peak SPL but also attenuates the sound at other peaks as well. This proves that the use of viscoelastic damping material will greatly attenuate the vibroacoustical response of the structure.
6.5. Acoustical Performance for a Combination of Several Viscoelastic Materials at Different Temperatures
The temperature is a significant external factor affecting the performance of viscoelastic damping material in a mechanical system. With an increase in temperature, the loss factor approaches its best performance towards the transition region and then decreases afterwards. In this example, the objective is to study the effects of both 3 M ISD110 material (that has a better damping performance) and the 3 M ISD112 over the chosen temperature between 40 and 60 degree Celsius. It is of interest to study the effect of the combination of these two materials on the damping of the structure.
To introduce the different viscoelastic materials, the sevenlayer sandwich beam (with the same parameters as in the previous example) is redesigned incorporating both damping materials (3 M ISD110 and ISD112). This system is compared to an identical structure with only one damping material (either 3 M ISD110 or ISD112). In the system including two viscoelastic materials, the outer damping layers (2nd and 6th) are 3 M ISD112 and the inner damping layer (5th) is the 3 M ISD110. The simply supported boundary condition is examined in this numerical example, and the temperature range is from 40 to 60 degree Celsius. The acoustical response is also calculated with the step input in the frequency domain (equivalent to impulse input in the time domain). Table 5 shows the first order natural frequency, the system loss factor, and the corresponding peak value (dB) of the sound pressure level over the temperature range with the simply supported boundary condition applied to the FEM model.

It can be seen that for the same damping material, as the ambient temperature is increased, the value of SPL increases while the loss factor decreases.
7. Conclusions
A framework for conducting vibroacoustical analysis for multiplelayer beam structures containing different types of viscoelastic materials is presented in this paper. Several observations and conclusions can be drawn from the results of this research.(1)The vibration section of the proposed analysis consists of FEM model of multiplelayered damping beam incorporating the Biot damping model. The FEM model of the beam structure can be extended to more complicated damping structures using the same procedure. The nonlinear curvefitting technique accurately estimates the Biot constants. The Biot damping model can then be solved using the decoupling transformation to yield the frequencyspectrum analysis.(2)The Biot damping model is also capable of improving a structure’s damping performance by adding new features such as different viscoelastic materials and the variation of operating temperature. The result obtained through the procedure of vibration analysis discussed in this paper compares well to the closedform solution from a previous work. The first peak from the frequency spectrum is the predominant cause of the vibration issue in this damping structure.(3)The direct boundary element method of analysis for acoustical cavity applied under anechoic boundaries was chosen as the basis for predicting the particle velocity from the frequencyspectrum analysis. The acoustical result validates the frequencyspectrum result from vibration analysis and has good agreement with the predicted SPL spectrum of the identical sandwich beam without damping calculated by commercial software.
Nomenclature
Mass matrix  
, :  Elastic stiffness/viscous stiffness matrix 
:  Damping matrix 
:  Displacement vector 
:  Force vector 
, :  Coefficient matrix of state equation 
:  Dissipation coordinate vector 
, :  Number of minioscillators for first/second type of viscoelastic material 
:  Laplace variable 
:  Time 
:  Young’s modulus 
:  Shear modulus 
, :  FEM shape function of longitudinal/transverse deflection 
:  Number of DOF 
:  Density of material 
:  Thickness of layer 
:  Length of beam 
:  Eigenvector matrix 
:  Eigenvalue matrix 
, :  Biot constants 
:  Sound pressure 
:  Velocity vector 
:  Nodal normal component of boundary velocity. 
Disclosure
The authors (D. Rao and F. Lin) hereby declare that they do not have any direct or indirect financial relation leading to any conflict of interests with the commercial identities (BEM software, FEM software, Auto2fit, MATLAB, and Mathematic) mentioned in the text of their paper.
References
 E. M. Kerwin, “Damping of flexural waves by a constrained viscoelastic layer,” Journal of the Acoustical Society of America, vol. 31, pp. 952–962, 1959. View at: Google Scholar
 D. Ross, E. E. Ungar, and E. M. Kerwin, “Damping of plate flexural vibration by means of viscoelastic laminae,” in Structural Dampinga Colloquium on Structural Damping Held at the ASME Annual Meeting, pp. 49–87, 1959. View at: Google Scholar
 R. A. Ditaranto, “Thery of vibratory bending for elastic and viscoelastic layered finitelength beams,” Journal of Applied Mechanics, vol. 32, pp. 881–886, 1965. View at: Google Scholar
 D. J. Mead and S. Markus, “The forced vibration of a threelayer, damped sandwich beam with arbitrary boundary conditions,” Journal of Sound and Vibration, vol. 10, no. 2, pp. 163–175, 1969. View at: Google Scholar
 D. K. Rao, “Frequency and loss factors of sandwich beams under various boundary conditions,” Journal of Mechanical Engineering Science, vol. 20, no. 2, pp. 271–282, 1978. View at: Google Scholar
 E. T. Cottle, Damping of layered beams with mixed boundary conditions [M.S. thesis], Air Force Institute of Technology, 1990.
 M. Falugi, “Analysis of a fivelayer, viscoelastic, constrainedlayer beam,” in Proceedings of the Damping Workshop Paper CCB, 1991. View at: Google Scholar
 M. Parin, L. C. Rogers, and M. Falugi, “Practical stand off damping treatment for sheet metal,” in Proceedings of the Damping Workshop Paper IBA, 1989. View at: Google Scholar
 L. C. Rogers and M. Parin, “Experimental results for standoff passive vibration damping treatment,” in Proceedings of the Smart Structures and Materials 1995: Passive Damping, pp. 374–383, March 1995. View at: Google Scholar
 J. M. Yellin, I. Y. Shen, P. G. Reinhall, and P. Y. H. Huang, “Experimental investigation of a passive standoff layer damping treatment applied to an EulerBernoulli beam,” in Proceedings of the 1999 Smart Structures and Materials—Passive Damping and Isolation, vol. 3672 of Proceedings of SPIE, pp. 228–233, March 1999. View at: Google Scholar
 J. M. Yellin, I. Y. Shen, P. G. Reinhall, and P. Y. H. Huang, “An analytical and experimental analysis for a onedimensional passive standoff layer damping treatment,” Journal of Vibration and Acoustics, vol. 122, no. 4, pp. 440–447, 2000. View at: Google Scholar
 J. M. Yellin and I. Y. Shen, “An analytical model for a passive standoff layer damping treatment applied to an EulerBernoulli beam,” in Smart Structures and Materials 2002: Damping and Isolation, Proceedings of SPIE, pp. 349–357, June 1998. View at: Publisher Site  Google Scholar
 B. Yang and C. A. Tan, “The transfer functions of one dimensional distributed parameter systems,” Journal of Applied Mechanics, vol. 116, pp. 341–349, 1959. View at: Google Scholar
 Q. Chen and Y. W. Chan, “Integral finite element method for dynamical analysis of elasticviscoelastic composite structures,” Computers and Structures, vol. 74, no. 1, pp. 51–64, 2000. View at: Publisher Site  Google Scholar
 G. A. Lesieutre and U. Lee, “A finite element for beams having segmented active constrained layers with frequencydependent viscoelastics,” Smart Materials and Structures, vol. 5, no. 5, pp. 615–627, 1996. View at: Publisher Site  Google Scholar
 S. W. Kung and R. Singh, “Vibration analysis of beams with multiple constrained layer damping patches,” Journal of Sound and Vibration, vol. 212, no. 5, pp. 781–805, 1998. View at: Google Scholar
 Q. J. Zhang and M. G. Sainsbury, “The Galerkin element method applied to the vibration of rectangular damped sandwich plates,” Computers and Structures, vol. 74, no. 6, pp. 717–730, 2000. View at: Publisher Site  Google Scholar
 R. L. Bagley and P. J. Torvik, “Fractional calculusa different approach to the analysis of viscoelastically damped structures,” AIAA Journal, vol. 21, no. 5, pp. 741–748, 1983. View at: Google Scholar
 M. Hao, Vibration analysis of constrained layered beams with multiple damping layers [Ph.D. thesis], Michigan Tech, 2005.
 G. A. Lesieutre, E. Bianchini, and A. Maiani, “Finite element modeling of onedimensional viscoelastic structures using anelastic displacement fields,” Journal of Guidance, Control, and Dynamics, vol. 19, no. 3, pp. 520–527, 1996. View at: Google Scholar
 M. A. Biot, “Variational principles in irreversible thermodynamics with application to viscoelasticity,” Physical Review, vol. 97, no. 6, pp. 1463–1469, 1955. View at: Publisher Site  Google Scholar
 J. Zhang and G. T. Zheng, “The Biot model and its application in viscoelastic composite structures,” Journal of Vibration and Acoustics, vol. 129, no. 5, pp. 533–540, 2007. View at: Publisher Site  Google Scholar
 D. J. McTavish, “Shock response of a damped linear structure using GHM finite elements,” in Proceedings of the 44th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, pp. 1681–1689, April 2003, Paper 1591. View at: Google Scholar
 L. Zhang, H. P. Du, Y. M. Shi, and X. Z. Shi, “Parametric determination for GHM of ZN1 viscoelastic material,” Rare Metal Materials and Engineering, vol. 31, no. 2, pp. 91–95, 2002. View at: Google Scholar
 C. H. Park, D. J. Inman, and M. J. Lam, “Model reduction of viscoelastic finite element models,” Journal of Sound and Vibration, vol. 219, no. 4, pp. 619–637, 1999. View at: Google Scholar
 M. Hao and M. D. Rao, “Vibration and damping analysis of a sandwich beam containing a viscoelastic constraining layer,” Journal of Composite Materials, vol. 39, no. 18, pp. 1621–1643, 2005. View at: Publisher Site  Google Scholar
 D. H. Lee, “Optimal placement of constrainedlayer damping for reduction of interior noise,” AIAA Journal, vol. 46, no. 1, pp. 75–83, 2008. View at: Publisher Site  Google Scholar
 F. Lin and M. D. Rao, “Vibration analysis of a multiplelayered viscoelastic structure using the biot damping model,” AIAA Journal, vol. 48, no. 3, pp. 624–634, 2010. View at: Publisher Site  Google Scholar
 T. W. Wu, Boundary Element Acoustics: Fundamentals and Computer Codes, WIT Press, Ashurst, UK, 2000.
Copyright
Copyright © 2013 Fei Lin and Mohan D. Rao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.