Next Article in Journal
Biocompatibility and Osteogenic Potential of Calcium Silicate-Based Cement Combined with Enamel Matrix Derivative: Effects on Human Bone Marrow-Derived Stem Cells
Next Article in Special Issue
Active Vibration Control of Piezoelectric Sandwich Plates
Previous Article in Journal
Analysis of Multiphase Flow and Heat and Mass Transfer for Ammonium Chloride Crystallization of the High-Pressure Heat Exchanger in Hydrogenation Unit
Previous Article in Special Issue
Structural Lattice Topology and Material Optimization for Battery Protection in Electric Vehicles Subjected to Ground Impact Using Artificial Neural Networks and Genetic Algorithms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transverse Vibration of Viscoelastic Sandwich Structures: Finite Element Modeling and Experimental Study

1
School of Mechanical and Electronic Engineering, Jingdezhen Ceramic University, Jingdezhen 333001, China
2
Department of Mechanical Engineering, Tsinghua University, Beijing 100084, China
*
Author to whom correspondence should be addressed.
Materials 2021, 14(24), 7751; https://doi.org/10.3390/ma14247751
Submission received: 19 October 2021 / Revised: 20 November 2021 / Accepted: 8 December 2021 / Published: 15 December 2021
(This article belongs to the Special Issue Sandwich Composites: Design, Simulation and Applications)

Abstract

:
In the present work, the nonlinear vibration behavior of elastic-viscoelastic-elastic sandwich (EVES) beams is studied. A finite element (FE) equation taking intoaccount the transverse compression deformation of the viscoelastic core for the EVES beams is derived. In order toaccurately characterize the frequency-dependent feature of the viscoelastic materials layer, athird-order seven-parameter Biot model isused. A 2-node 8-DOF element is established to discretize the EVES beams. The experimental testing onEVES beams validates the numerical predication of the FE model. Numerical and analytical investigations are carried on a series of EVES beams with different thicknesses. The results indicate that the presented FE model has better accuracy in predicting the natural frequency of the sandwich beams, and in predicting damping, the accuracy is related to the thickness of each layer. The results of this paper have important reference values for the design and optimization of the viscoelastic sandwich structure.

1. Introduction

Sandwich structure is a special type of composite structure. It has the ability to significantly reduce weight while maintaining mechanical properties. That is, it has typical characteristics of light weight, high rigidity, and high strength [1].The weight reduction in the sandwich structure brings many benefits, including increased stroke, greater load, and reduced fuel consumption. All of these have a positive impact on costs and reduce the impact on the environment. In addition to the high stiffness-to-weight ratio, sandwich composites can also achieve different functions by choosing different core materials. For example, the use of foam [2,3,4], functionally graded materials [5,6,7], and viscoelastic materials [8,9,10] for the core material can achieve the effects of flame retardant, heat insulation, vibration, and noise reduction. Therefore, sandwich composite structures are more and more widely used in aerospace, automotive, marine, naval industries, and civil engineering.
Figure 1 shows a typical elastic-viscoelastic-elastic sandwich (EVES) beam structure. When the base beam vibrates, the viscoelastic layer will produce the corresponding deformation to dissipate vibration energy. However, what deformation of the viscoelastic layer dissipates energy is a question. There are two main assumptions for the energy dissipation mode of the viscoelastic layer: shear and compression energy dissipation assumptions. When the relative movement of the constraining layer and the base beam layer is parallel to the neutral plane of the sandwich beam, the viscoelastic core layer will produce shear deformation to dissipate vibration energy. When the relative movement of the constrained layer and the base beam layer is perpendicular to the neutral plane of the beam, the viscoelastic core layer will produce compression-tension deformation to dissipate vibration energy. Simply put, the former considers the viscoelastic layer to be incompressible, and the latter considers the viscoelastic layer to be compressible. The former dissipates vibration energy through shear deformation, and the latter dissipates vibration energy through compression deformation.
Most studies adopt the shear energy dissipation assumption. They believe that the transverse displacement of all points on the same cross-section of EVES structure is the same, so there is no compression/tensile deformation in the viscoelastic layer. The relevant research work can be traced back to the 1950s. Ross Kerwin Ungar’s (RKU) model is considered the first analytical model of EVES structure [11,12,13,14]. Later, Ditaranto extended Kerwin’s work and proposed a sixth-order linear homogeneous differential equation to study the free vibration of sandwich beams. However, the calculation results of the model show that the structural loss factor does not depend on the boundary conditions, which is obviously not in line with reality [15]. Mead and Markus modified Ditaranto’s model and proposed the famous Mead-Markus model, which can be applied to various boundary conditions [16,17]. The above-mentioned RKU model and Mead-Markus model were based on the shear assumption in the processing of the displacement domains of the three different layers of the EVEC beam structure. Later, many scholars such as Johnson et al. [18], Galucio et al. [19], Kumar et al. [20], Bilasse et al. [21], Kpeky et al. [22], Hamdaoui et al. [23], and Huang et al. [24,25] developed the FE method for viscoelastic sandwich structure, which were still based on the assumption of shear energy dissipation. Recently, some scholars have also carried out related researches on viscoelastic sandwich structures based on the shear hypothesis. Karmi et al. [26] carried out a dynamic analysis of sandwich beams with a viscoelastic core subjected to a moving load, and the mathematical formulation was based on shear deformation theory. Mario et al. [27] studied the fractional viscoelastic characterization of sandwich beams under time-varying loading. The fractional viscoelasticity is incorporated in the analytical model of a three-layer laminated beam. Maleki-Bigdeli et al. [28] presented an analytical model for the bending problem of a viscoelastic sandwich plate. The generalized Maxwell model was used to describe the viscoelastic response. Amanieh et al. [29] studied the nonlinear vibrationof a sandwich plate. The four-parameter fractional derivative model was used to determine the frequencydependence of the viscoelastic core. Garbowski et al. [30,31] studied the transverse shear stiffness in sandwich plates.
The above-referred papers deal with the analysis of the sandwich structures under theshear energy dissipation assumption. On the other hand, some scholars have found the compression damping of EVES structure. In the 1970s and 1980s, Douglas and Yang [32,33] proved the existence of compression damping through experiments on EVES beams. In addition, based on the assumption of compression energy dissipation, an analytical model (Douglas-Yang model) of the EVES beam was established, ignoring the shear deformation of the viscoelastic layer. They believed that compression damping is the main form in a narrow frequency band centered on the compression resonance frequency of the viscoelastic layer. However, their work did not attract people’s attention at the time. Later, Sisemore et al. [34,35] conducted a more in-depth study on the cantilever EVES beam structure. They not only proved the existence of compression damping experimentally but also established an analytical model of the EVES beam structure based on the assumption of compression energy dissipation. The Sisemore model assumed that the viscoelastic layer was a compressible spring, and the vibration energy was dissipated through the transverse compression of the spring, and there was no shear effect. They proved that the analytical model can predict the resonance frequency of the structure well, but it cannot predict the loss factor well Funari et al. [36] presented a nonlinear approach to investigate the behavior of composite sandwich structures. The core was considered compressible.
It can be seen from the above literature that the research based on the shear energy dissipation assumption of the viscoelasticis relatively in-depth. However, the research based on the assumption of compression energy dissipation is still in its infancy. There are only a few analytical models, and FE models are rarely seen. In engineering practice, EVES structures that need to be designed or analyzed often have complex geometric shapes and boundary conditions and bear complex loads. The analytical method cannot effectively solve these problems. The FE method has shown its excellent applicability and has become the most important method for studying the vibration problems of EVES structures. The study of the FE compression model of the EVES structures has great engineering significance.
In order to make up for this deficiency, in this work, a FE model for EVES beam is established based on the compression energy dissipation assumption. A sandwich composite element with two-node eight degrees of freedom is established. In the process of modeling, it is considered that the energy dissipation is caused by the compression deformation of the viscoelastic layer, and the transverse displacement of the viscoelastic layer is a linear interpolation between the transverse displacement of the constraining layer and base beam layer. The Biot model is used to characterize the frequency-dependent characteristics of viscoelastic materials. Finally, the FE model is verified by comparing with the analytical solutions and experimental values, and some useful conclusions are obtained.

2. FE Model for an EVES Beam

2.1. Fundamental Assumptions

The present analysis is based on the following assumptions:
  • The structural damping is only caused by the transverse compression/tensile deformation of the viscoelastic sandwich layer;
  • The constraint layer and the base beam are regarded as Euler-Bernoulli beams;
  • Considering the compression deformation of the viscoelastic layer perpendicular to the neutral plane of the beam, it is considered that the base beam layer, the damping layer, and the constraint layer have different deflection functions, and the transverse displacement of the viscoelastic layer is the linear interpolation of the two surface layers;
  • The materials of each layer are firmly pasted, and there is no relative sliding between the layers;
  • The influence of the moment of inertia of each layer is ignored.

2.2. Kinematics

The geometric and deformation relationship of the EVES beam is shown in Figure 2.
According to the above assumptions, for the viscoelastic core layer, the transverse displacement of the mid-plane and the compressive strain can be respectively written as:
y 2 = 1 2 ( y 1 + y 3 ) ,   ε 2 = w 1 w 3
where y 2 and ε 2 are the transverse displacement and the compressive strain of the viscoelastic layer, respectively. y 1 and y 3 are the transverse displacement of the beam and the constraining layer, respectively.

2.3. Degrees of Freedom and Shape Functions

The elementof the EVES beam is established, as shown in Figure 3. It is a sandwich composite element with two nodes. The length and width of the element are l e and b , respectively. The thickness of the base beam, viscoelastic layer, and the constraining layer of the element are h 1 , h 2 , and h 3 , respectively. Each node of the element has four degrees of freedom (DOF), representing the transverse displacement and the rotation angle of the base beam layer and the constraining layer, respectively.
The node displacement vector of the element is given by
Δ e = { y 1 i θ 1 i y 3 i θ 3 i y 1 j θ 1 j y 3 j θ 3 j } T
In the element, the transverse displacement y 1 and y 3 , the rotation angle θ 1 and θ 3 of the base beam layer and the constraining layer are expressed in the nodal displacements by FE shape functions
y 1 = N 1 Δ e ,   θ 1 = N 2 Δ e ,   y 3 = N 3 Δ e ,   θ 3 = N 4 Δ e
N 1 = [ 2 ( x l e ) 3 3 ( x l e ) 2 + 1 x 3 l e 2 2 ( x 2 l e ) + x 0 0 2 ( x l e ) 3 + 3 ( x l e ) 2 x 3 l e 2 x 2 l e 0 0 ]     N 2 = [ 6 ( x 2 l e 3 ) 6 ( x l e 2 ) 3 ( x l e ) 2 4 ( x l e ) + 1 0 0 6 ( x 2 l e 3 ) + 6 ( x l e 2 ) 3 ( x l e ) 2 2 ( x l e ) 0 0 ] N 3 = [ 0 0 2 ( x l e ) 3 3 ( x l e ) 2 + 1 x 3 l e 2 2 ( x 2 l e ) + x 0 0 2 ( x l e ) 3 + 3 ( x l e ) 2 x 3 l e 2 x 2 l e ]      N 4 = [ 0 0 6 ( x 2 l e 3 ) 6 ( x l e 2 ) 3 ( x l e ) 2 4 ( x l e ) + 1 0 0 6 ( x 2 l e 3 ) + 6 ( x l e 2 ) 3 ( x l e ) 2 2 ( x l e ) ]
where the shape functions are given by substituting Equation (3) into Equation (1) gives
y 2 = N 5 Δ e ,   ε 2 = N 6 Δ e
where
N 5 = 1 2 ( N 1 + N 3 ) ,   N 6 = N 1 N 3

2.4. Energy Terms

2.4.1. The Strain Energy

The potential energy from the bending of the base beam layer and the constraining layer are respectively written as
U 1 = 1 2 E 1 I 1 0 l e ( 2 y 1 x 2 ) 2 d x = 1 2 Δ e T K 1 e Δ e ,   U 3 = 1 2 E 3 I 3 0 l e ( 2 y 3 x 2 ) 2 d x = 1 2 Δ e T K 3 e Δ e
where U 1 and U 3 , E 1 and E 3 , I 1 and I 3 and the potential energy from bending, the elastic modulus, and the moment of inertia of the base beam layer and the constraining layer, respectively. K 1 e and K 3 e are the stiffness matrix of the base beam and the constraining layer, respectively. Applying the shape functions, they can be expressed in the form
K 1 e = E 1 I 1 0 l e [ 2 N 1 x 2 ] T [ 2 N 1 x 2 ] d x ,   K 3 e = E 3 I 3 0 l e [ 2 N 3 x 2 ] T [ 2 N 3 x 2 ] d x
The viscoelastic layer is transversely compressed or stretched to dissipate energy. Therefore its potential energy from compression or stretching can be expressed as
U 2 = 1 2 E 2 b h 2 0 l e ( y 1 y 3 ) 2 d x = 1 2 Δ e T K 2 e Δ e
where U 2 , E 2 , and K 2 e are the potential energy, the elastic modulus, and the stiffness matrix of the viscoelastic layer, respectively. Applying the shape functions, and the K 2 e can be obtained as
K 2 e = E 2 b h 2 0 l e ( N 6 T N 6 ) d x = G v 2 b ( 1 + ν 2 ) h 2 0 l e ( N 6 T N 6 ) d x = G v K v v e K v v e = 2 b ( 1 + ν 2 ) h 2 0 l e ( N 6 T N 6 ) d x                     
where G v , v 2 , and K v v e are the shear modulus, the Poisson’s ratio and the viscous stiffness matrix of the viscoelastic core layer.
Obviously, the total stiffness matrix of the element is the sum of the stiffness matrix of each layer. It can be written as
K e = K 1 e + K 3 e K e e + K 2 e
where the sum of the first two terms on the right side of the equation is the elastic stiffness matrix K e e .
The total potential energy of the element is the sum of the potential energy of each layer. It can be written as
U = U 1 + U 2 + U 3

2.4.2. The Kinetic Energy

The kinetic energy from the transverse vibration of the element is given by
T 1 = 1 2 ρ 1 A 1 0 l e ( y 1 t ) 2 d x = 1 2 Δ ˙ e T M 1 e Δ ˙ e T 2 = 1 2 ρ 2 A 2 0 l e ( y 2 t ) 2 d x = 1 2 Δ ˙ e T M 2 e Δ ˙ e T 3 = 1 2 ρ 3 A 3 0 l e ( y 3 t ) 2 d x = 1 2 Δ ˙ e T M 3 e Δ ˙ e
where T 1 , T 2 , and T 3 , M 1 e , M 2 e , and M 3 e are the kinetic energy and mass matrices of the base beam layer, the viscoelastic layer, and the constraining layer, respectively. These mass matrices are written as
M 1 e = ρ 1 A 1 0 l e N 1 T N 1 d x ,   M 2 e = ρ 2 A 2 0 l e N 5 T N 5 d x ,   M 3 e = ρ 3 A 3 0 l e N 3 T N 3 d x
where ρ 1 , ρ 2 and ρ 3 , A 1 , A 2 ,   and A 3 are the density and the cross-sectional area of the base beam layer, the viscoelastic layer, and the constraining layer, respectively.
Obviously, the total mass matrix of the element is the sum of the mass matrix of each layer.
M e = M 1 e + M 2 e + M 3 e
The total kinetic energy of the element is the sum of the kinetic energy of each layer.
T = T 1 + T 2 + T 3

2.5. Equation of Motion

The FE dynamic equation of viscoelastic composite structure can be expressed in Laplace domain as [15,16]
( s 2 M e + K e e + s G ˜ ( s ) K v v e ) X ( s ) = F e ( s )
where x ( s ) is the displacement vector, F e ( s ) is the excitation vector.   s G ˜ ( s ) = G * ( s ) is the complex shear modulus of the viscoelastic material, which is frequency dependent.
The frequency-dependent characteristics of viscoelastic materials can be expressed by the Biot model
s G ˜ ( s ) = G [ 1 + i = 1 N a i s s + b i ]
where G is the steady-state value of the shear modulus of the viscoelastic material, N is the number of series items, { a i , b i } are positive constants.
Substituting the Biot model Equation (18) into Equation (17) and then introducing an auxiliary coordinate Z ^ k ( s ) = a k s + b k X ( s ) , ( k = 1 ,   2 ,   3 N ) , the Formula (17) can be extended as
M - q .. + D - q . + K - q = f -
where
M - = [ M e 0 0 0 0 0 0 0 0 ] ,   D - = [ 0 0 0 0 a 1 b 1 Λ 0 0 0 a N b N Λ ] K - = [ K e e + K ˜ ( 1 + K = 1 N a k ) a 1 R a N R a 1 R T a 1 Λ 0 a N R T 0 a N Λ ] ,   q = { x Z 1 Z N } ,   f - = { F e 0 0 }
where Λ v is a diagonal matrix, which is composed of the positive eigenvalues of the viscous stiffness matrix K v v e , and R v is a matrix with the corresponding orthogonal eigenvector as the column, Z j = R v T Z ^ j ,   ( j = 1 ,   2 ,   , N ) .
Assembling all of the elements, the overall dynamic equation of the viscoelastic sandwich beam structure can be obtained as
M x ¨ + D x ˙ + K x = F
where M ,   D and K   are the overall mass matrix, damping matrix, and stiffness matrix, respectively, F is the external force of the viscoelastic sandwich beam.
The numerical strategy used to extract the transverse vibration and damping characteristics was described in detail in reference [37].

3. Experiment

The experiment is designed to test the response of an EVES beam. Figure 4 shows the experimental setup. In the experiment, use a self-made lightweight, soft hammer with a rubber head to lightly hit the end of the EVES beam to excite vibration. A lightweight acceleration sensor is pasted at the free end to measure the vibration signal, which is collected and adjusted by the LMS SCADAS III 316 data acquisition system, and then analyzed by the LMS test. Lab vibration and noise test analysis system to obtain the natural frequencies and damping ratios.
The material and mechanical properties are shown in Table 1. The viscoelastic layer material is ZN-1, and its frequency-dependent Biot model parameters are given in reference [24]. The FE compression model is used to calculate the first three order natural frequencies and damping ratios of the EVES beam. A total of 30 finite elements are used for discreting it. The experimental and calculation results are listed in Table 2.
It can be seen from Table 2 that in the prediction of natural frequency, the calculation results are in suitable agreement with the experimental results. The prediction error is less than 4%. However, in the prediction of damping, the FE model shows deficiencies. It is only adequate for the first damping ratio, and for the second and third damping ratios, it shows a poor predictor of damping for the tested beam.

4. Numerical and Analytical Investigations

In order to study the compressive damping of EVES beams, Sisemore [35] carried out experiments on cantilever EVES beams with different thicknesses. They established an analytical model to study the transverse vibration. The thickness of the base beam layer of these beams is the same. The constraining layer and the viscoelastic layer have different thicknesses. The viscoelastic material is EAR-C1002. Table 3 shows the material parameters of the EVES beams, and Table 4 shows the thickness parameters.
Reference [38] shows that at the temperature of 23.9 °C, and the frequency range of f = 0 ,   ,   2000 Hz, the storage modulus and loss factor of the complex shear modulus of the viscoelastic material layer (EAR-C1002) are, respectively,
G v = ( 44.4 17.6 / a ) N / mm 2
η v = 1.643 0.6025 z 2 0.2557 × 10 20 z 16 + 0.1260 × 10 9 z 8 0.1959 × 10 4 z 4
where a = 0.4 + 0.0003 f , z = 0.05 + 0.000475 f .
The FE method presented in this paper is used to calculate the natural frequencies and damping ratios corresponding to the first two modes of the 9 beams listed in Table 2. The Biot model is used to describe the frequency-dependent characteristics of viscoelastic materials. Based on Equations (22) and (23), the Biot model parameters (Equation (18)) of EAR-C1002 can be fitted according to the method introduced in reference [24]. Table 5 lists these parameters.
Taking the experimental data in the reference [35] as the standard, the calculation results of the FE model presented in this paper are compared with the calculation results of the Sisemore model [35] and Mead-Markus model [16]. The last two models are both analytical models, in which the Sisemore modelisthe compression model, and the Mead-Markus model is the shear model. When calculating with the FE method in this paper, the EVES beams are discretized into 31 elements. The calculation results of the first two natural frequencies and damping ratios are presented in Table 6 and Table 7, respectively.
It can be seen from Table 6 that the accuracy of the compressional model is better than that of the Mead-Markus model in predicting the first two natural frequencies. Most of the prediction errors of the 9 beams are within 10%, the error range is 0.73–14%, and the average error is 7%. In comparison, the error of the Mead-Markus model is mostly larger than that of the compression model. The error range is 13–51%, and the average error is 28.5%. The error range of the Sisemore model is 0.1–11%, and the average error is 6%. The accuracy of the compressional model is close to the Sisemoremodel. This is because they are based on the same assumption of compressive energy dissipation.
According to Euler-Bernoulli beam theory, the first two natural frequencies of the base beam without additional damping layer and constraint layer is given by
f n = A n 2 π L 2 E I ρ A
where f n ( n = 1 ,   2 ) is the first two natural frequencies of the base beam (Hz), E is the Young’s modulus, I is the moment of inertia of section, ρ is the density, A is the cross-sectional area, A n = { 3.52 ,   22.4 } is modal factor. Substituting the material parameters of the base beam in Table 1 into Equation (24), the first two natural frequencies of the base beam are 53 and 333 Hz. With the addition of the viscoelastic layer and the constraining layer, the stiffness and mass increase accordingly. It can be seen from Equation (24) that the increase in additional stiffness will cause the natural frequency of the system to increase, but the increase in additional mass will lead to the decrease in the natural frequency. It can be seen from Table 1 that the density of the viscoelastic layer is about 50% of the base beam. However, the stiffness of the viscoelastic layer is only 0.01% of it. Obviously, the effect of the additional mass on the natural frequency far exceeds the additional stiffness. Therefore, after adding a viscoelastic layer, the natural frequency of the structure will decrease. It can be seen from Table 4 that the natural frequency values obtained by the experiment are all smaller than the theoretical value of the natural frequency of the base beam, which is in line with the above analysis. The predicted value of the natural frequency of the FE compression model is very close to the experimental results. On the contrary, the natural frequencies of most sandwich beams predicted by the shear model are higher than those of the base beam, and the error is larger than that of the compression model. Therefore, it is obvious that under the geometric conditions, the shear model will overestimate the natural frequency of the EVES beams, while the compression model has higher accuracy.
Table 7 presents a comparison of the predicted and experimental damping ratios of the nine EVES beams. The results indicate that all three models are not accurate enough to predict the damping ratio. However, some laws can still be seen from it. In the prediction of the damping ratio corresponding to the first-order natural frequency of the beams, the compression damping model only predicts beams 3, 6, and 9 relatively accurately. The common point of the three beams is that their constraint layers have the maximum thickness. This shows that when the stiffness of the constrained layer increases, the prediction accuracy of the compression damping model will improve. This is because when the stiffness of the constraint layer increases, it helps to make the viscoelastic layer produce transverse tension and compression deformation, which is closer to the assumption of the compression damping model. When the constraining layer is relatively thin, its stiffness is small, and at this time, it will produce bending deformation under the action of the viscoelastic layer, which is contrary to the assumption of the compression damping model. Therefore, in beams 3, 6, and 9, compression damping accounts for the main proportion. In addition, it can be seen that in these three beams, when the viscoelastic layer is thicker, the compression damping is more obvious, and the prediction accuracy of the FE compression model for the damping ratio is higher. Similarly, in beams 1, 4, and 7, the thickness of the constrained layer is the smallest, and the prediction accuracy of the FE compression model is the worst. In addition, the average prediction error of theMead-Markus model for the first-order damping ratio exceeds 130%. Obviously, the shear model overestimates the first-order damping ratio of the EVES beams. This overestimation may have a bad influence on the design and optimization of the EVES beam structures.
In the prediction of the second-order damping ratio, the compression models for all beams are distorted. However, in comparison, the Mead-Markus model is more adequate. This shows that in the second-order mode, the shear damping accounts for the main proportion, and in this case, the compression model is not applicable. This can be explained by the vibration mode of the cantilever beam. In the second-order vibration mode, the constraining layer and the foundation beam are more likely to produce relative axial displacement. In this case, the shear deformation of the viscoelastic layer dissipates more vibration energy, which means that the shear damping is dominant. This also shows that the vibration mode is also a factor that determines the damping mechanism of the EVES beam. This also shows the complexity of the damping mode of the viscoelastic sandwich beam, which is related not only to the thickness of each layer of the structure but also to the vibration mode.

5. Conclusions

A FE compressional model of EVES beams is established to study its transverse vibration and damping characteristics. The FE is a three-layer two-node beam element with four degrees of freedom at each node. It is considered that the viscoelastic material is compressible, and the damping of the EVES beam is only caused by the transverse compression/tension of the viscoelastic layer.Biot model is used to describe the frequency-dependent characteristics of viscoelastic materials. Nine numerical examples are introduced to compare the calculation results of the FE model in this paper with those of a compressional analytical model and a shear analytical model. Finally, an EVES beam is experimentally studied. Some conclusions are obtained:
  • The transverse compressional vibration is a broadband existence in EVES beams;
  • The FE compressional model presented inthis paper is in suitable agreement with the analytical model. It is a relatively accurate and simple method for predicting the natural frequencies of EVES beams. However, it is quite poor for damping predictions because it ignores shear damping;
  • The geometry and the vibration mode determine the damping mechanisms. The shear model and the compressional model have different applicable conditions, which are related to the thickness of each layer and the vibration mode of the structure. In general, the shear energy dissipation assumption is applicable to EVES thin-walled beams, and the compression energy dissipation assumption is applicable to EVES beams with relatively thick constraining layer and base beams;
  • When the base beam and the constrained layer are relatively thick, and the viscoelastic layer is relatively thin, neither the current shear model nor the compression model can accurately predict the damping. A model including shear damping and compressional damping should greatly improve the accuracy of damping prediction, which needs further research.
The EVES structure studied in this paper is in a simple environment. In fact, they sometimes work in a state of motion or a special environment. When they are moving in a straight line or rotating, their dynamic characteristics need further research. In addition, when the structure is in strong airflow, flutter may occur, and buckling may occur in an aero-thermal environment. These special behaviors also need to be further studied. The work proposedin this paper can be used in these studies.

Author Contributions

Z.Y., Data Curation; F.C. and X.W., Proofreading of the grammar; Z.Y. and J.P., Designed the experiment; Z.H., Project administration; X.W., Supervision; F.C., Validation; J.P. and Z.H., Writing—original draft; Z.H., Writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Science Foundation of China (11862007), Science and Technology Projects of Jiangxi Education Department of China, Jingdezhen Science and Technology Bureau Project, and Innovation and entrepreneurship training Project (S202110408057) for College Students of Jingdezhen Ceramic University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ma, Q.J.; Rejab, R.; Siregar, J.P.; Guan, Z.W. A review of the recent trends on core structures and impact response of sandwich panels. J. Compos. Mater. 2021, 55, 2513–2555. [Google Scholar] [CrossRef]
  2. Bethke, C.; Weber, L.; Goedderz, D.; Standau, T.; Altstädt, V. Fire behavior of flame retarded sandwich structures containing PET foam cores and epoxy face sheets. Polym. Compos. 2020, 41, 5195–5208. [Google Scholar] [CrossRef]
  3. Burlayenko, V.N.; Sadowski, T. Simulations of Post-impact Skin/core Debond Growth in Sandwich Plates Under Impulsive Loading. J. Appl. Nonlinear Dyn. 2014, 3, 369–379. [Google Scholar] [CrossRef]
  4. Funari, M.F.; Spadea, S.; Lonetti, P.; Loureno, P.B. On the elastic and mixed-mode fracture properties of PVC foam. Theor. Appl. Fract. Mech. 2021, 112, 102924. [Google Scholar] [CrossRef]
  5. Liu, Y.; Qin, Z.; Chu, F. Nonlinear forced vibrations of FGM sandwich cylindrical shells with porosities on an elastic substrate. Nonlinear Dynam. 2021, 104, 1007–1021. [Google Scholar] [CrossRef]
  6. Liu, Y.; Qin, Z.; Chu, F. Nonlinear dynamic responses of sandwich functionally graded porous cylindrical shells embedded in elastic media under 1:1 internal resonance. Appl. Math. Mech. 2021, 42, 805–818. [Google Scholar] [CrossRef]
  7. Koutoati, K.; Mohri, F.; Daya, E.M.; Carrera, E. A finite element approach for the static and vibration analyses of functionally graded material viscoelastic sandwich beams with nonlinear material behavior. Compos. Struct. 2021, 274, 114315. [Google Scholar] [CrossRef]
  8. Sun, C.T.; Sankar, B.V.; Rao, V.S. Damping and vibration control of unidirectional composite laminates using add-on viscoelastic materials. J. Sound Vib. 1990, 139, 277–287. [Google Scholar] [CrossRef]
  9. Hamdaoui, M.; Ledi, K.S.; Robin, G.; Daya, E.M. Identification of frequency-dependent viscoelastic damped structures using an adjoint method. J. Sound Vib. 2019, 453, 237–252. [Google Scholar] [CrossRef]
  10. Li, W.; He, Y.; Xu, Z.; Zhang, Z. A reduced passive constrained layer damping finite element model based on the modified improved reduced system method. J.Sandw. Struct. Mater. 2019, 21, 758–783. [Google Scholar] [CrossRef]
  11. Kerwin, E.M.J. Damping of flexural waves by a constrained viscoelastic layer. J. Acoust. Soc. Am. 1959, 31, 952–962. [Google Scholar] [CrossRef]
  12. Ross, D.; Ungar, E.E.; Kerwin, E.M. Damping of plate flexural vibrations by means of viscoelastic laminae. Struct. Damping 1959, 3, 44–87. [Google Scholar]
  13. Ungar, E.E. Loss factors of viscoelastically damped beam structures. J. Acoust. Soc. Am. 1962, 34, 1082–1089. [Google Scholar] [CrossRef]
  14. Ungar, E.E.; Kerwin, E.M.J. Loss factors of viscoelastic systems in terms of energy concepts. J. Acoust. Soc. Am. 1962, 34, 954–957. [Google Scholar] [CrossRef]
  15. Ditarant, R.A. Theory of vibratory bending for elastic and viscoelastic layered finite length beams. Int. J. Appl. Mech. 1965, 32, 881–886. [Google Scholar] [CrossRef]
  16. Mead, D.J.; Markus, S. The forced vibration of a three-layer, damped sandwich beam with arbitrary boundary conditions. J. Sound Vib. 1969, 10, 163–175. [Google Scholar] [CrossRef]
  17. Mead, D.J.; Markus, S. Loss factors and resonant frequencies of encastré damped sandwich beams. J. Sound Vib. 1970, 12, 99–112. [Google Scholar] [CrossRef]
  18. Johnson, C.D.; Kienholz, D.A. Finite Element Prediction of Damping in Structures with Constrained Layers. AIAA J. 1982, 20, 1284–1290. [Google Scholar] [CrossRef]
  19. Galucio, A.C.; Deu, J.F.; Ohayon, R. Finite element formulation of viscoelastic sandwich beamsusing fractional derivative operators. Comput. Mech. 2004, 33, 282–291. [Google Scholar] [CrossRef] [Green Version]
  20. Kumar, N.; Singh, S.P. Vibration and damping characteristics of beams with active constrainedlayer treatments under parametric variations. Mater.Des. 2009, 30, 4162–4174. [Google Scholar] [CrossRef]
  21. Bilasse, M.; Daya, E.M.; Azrar, L. Linear and nonlinear vibrations analysis of viscoelastic sandwich beams. J. Sound Vib. 2010, 329, 4950–4969. [Google Scholar] [CrossRef]
  22. Kpeky, F.; Boudaoud, H.; Abed-Meraim, F.; Daya, E.M. Modeling of viscoelastic sandwich beams using solid–shell finiteelements. Compos. Struct. 2015, 133, 105–116. [Google Scholar] [CrossRef] [Green Version]
  23. Hamdaoui, M.; Akoussan, K.; Daya, E.M. Comparison of non-linear eigensolvers for modal analysis of frequencydependent laminated visco-elastic sandwich plates. Finite Elem. Anal. Des. 2016, 121, 75–85. [Google Scholar] [CrossRef]
  24. Huang, Z.; Wang, X.; Wu, N.; Chu, F.; Luo, J. A Finite Element Model for the Vibration Analysis of Sandwich Beam with Frequency-DependentViscoelastic Material Core. Materials 2019, 12, 3390. [Google Scholar] [CrossRef] [Green Version]
  25. Huang, Z.; Wang, X.; Wu, N.; Chu, F.; Luo, J. TheFinite Element Modeling and Experimental Study of Sandwich Plates with Frequency-DependentViscoelastic Material Model. Materials 2020, 13, 2296. [Google Scholar] [CrossRef] [PubMed]
  26. Karmi, Y.; Khadri, Y.; Tekili, S.; Daouadji, A.; Daya, E.M. Dynamic Analysis of Composite Sandwich Beams with a Frequency-Dependent Viscoelastic Core under the Action of a Moving Load. Mech. Compos. Mater. 2021, 56, 755–768. [Google Scholar] [CrossRef]
  27. Paola, M.D.; Galuppi, L.; Carfagni, G.R. Fractional viscoelastic characterization of laminated glass beams under time-varying loading. Int. J. Mech. Sci. 2021, 196, 106274. [Google Scholar] [CrossRef]
  28. Maleki-Bigdeli, M.A.; Sheikhi, S.; Baghani, M. Development of an analytical framework for viscoelastic corrugated-core sandwich plates and validation against FEM. Meccanica 2021, 56, 2103–2120. [Google Scholar] [CrossRef]
  29. Amanieh, H.R.T.; Roknizadeh, S.A.S.; Reza, A. Nonlinear Vibration Analysis of Viscoelastic Smart Sandwich Plates Through the use of Fractional Derivative Zener Model. Int. J. Struct. Stab. Dyn. 2021, 21, 2150061. [Google Scholar] [CrossRef]
  30. Garbowski, T.; Gajewski, T.; Grabski, J.K. Torsional and Transversal Stifness of Orthotropic Sandwich Panels. Materials 2020, 13, 5016. [Google Scholar] [CrossRef]
  31. Garbowski, T.; Gajewski, T.; Grabski, J.K. Role of Transverse Shear Modulus in the Performance of Corrugated Materials. Materials 2020, 13, 3791. [Google Scholar] [CrossRef]
  32. Douglas, B.E.; Yang, J. Transverse compressional damping in the vibratory response of elastic-viscoelastic-elastic beams. AIAA J. 1978, 16, 925–930. [Google Scholar] [CrossRef]
  33. Douglas, B.E. Compressional damping in three-layer beams incorporating nearly incompressible viscoelastic cores. J. Sound Vib. 1986, 104, 343–347. [Google Scholar] [CrossRef]
  34. Sisemore, C.L.; Smaili, A.A.; Darvennes, C.M. Experimental measurement of compressional damping in an elastic-viscoelastic-elastic sandwich beam. Proc. Am. Soc. Mech. Eng. Noise Control Acoust. Div. 1999, 26, 223–227. [Google Scholar]
  35. Sisemore, C.L.; Darvennes, C.M. Transverse vibration of elastic-viscoelastic-elastic sandwich beams: Compression-experimental and analytical study. J. Sound Vib. 2002, 252, 155–167. [Google Scholar] [CrossRef]
  36. Funari, M.F.; Greco, F.; Lonetti, P. Sandwich panels under interfacial debonding mechanisms. Compos. Struct. 2018, 203, 310–320. [Google Scholar] [CrossRef]
  37. Huang, Z.; Qin, Z.; Chu, F. Damping mechanism of elastic–viscoelastic–elastic sandwich structures. Compos. Struct. 2016, 153, 96–107. [Google Scholar] [CrossRef]
  38. Barkanov, E. Transient response analysis of structures made from viscoelastic materials. Intern. J. Numer. Methods Eng. 1999, 44, 393–403. [Google Scholar] [CrossRef]
Figure 1. Elastic-viscoelastic-elastic sandwich beam structure.
Figure 1. Elastic-viscoelastic-elastic sandwich beam structure.
Materials 14 07751 g001
Figure 2. The geometric and deformation relationship of the EVES beam.
Figure 2. The geometric and deformation relationship of the EVES beam.
Materials 14 07751 g002
Figure 3. The element of the EVES beam.
Figure 3. The element of the EVES beam.
Materials 14 07751 g003
Figure 4. Experimental setup.
Figure 4. Experimental setup.
Materials 14 07751 g004
Table 1. Material parameters of experimental specimen.
Table 1. Material parameters of experimental specimen.
Material PropertiesConstraining LayerBase BeamViscoelastic Layer (ZN-1)
Elastic modulus (GPa)699699Biotmodel [24]
Density (kg/m3)270027001010
Poisson’s ratio0.30.30.3
Thickness (mm)4.914.854.91
Length (mm)290290290
Width (mm)252525
Table 2. Comparison of experimental results and FE calculations.
Table 2. Comparison of experimental results and FE calculations.
OrderExperimental ResultFE Model This Paper
Natural Frequency (Hz)Damping RatioNatural Frequency (Hz)Error (%)Damping RatioError (%)
152.50.135051.61.70.12656.3
2315.50.0846279.81.10.012485.3
3842.00.0546747.53.360.000998.4
Table 3. Structure and material parameters of cantilever EVESbeams.
Table 3. Structure and material parameters of cantilever EVESbeams.
Material PropertiesConstraining LayerBase BeamViscoelastic Material Layer(EAR-C1002)
Elastic modulus (GPa)7171Frequency dependent
Density (kg/m3)271027101280
Poisson’s ratio0.30.30.3
Length (mm)314314314
Width (mm)25.425.425.4
Table 4. Thickness of each layer of cantilever EVES beams.
Table 4. Thickness of each layer of cantilever EVES beams.
NumberThickness of Base Beam (mm)Thickness of the Viscoelastic Layer (mm)Thickness of the Constraining Layer (mm)
16.3500.3811.588
26.3500.3813.175
36.3500.3816.350
46.3503.0481.588
56.3503.0483.175
66.3503.0486.350
76.3506.3501.588
86.3506.3503.175
96.3506.3506.350
Table 5. Fitted parameters of the Biotmodel for EAR-C1002 at 23.9 °C.
Table 5. Fitted parameters of the Biotmodel for EAR-C1002 at 23.9 °C.
k = 1 k = 2 k = 3
G4 × 105
ak8.22441.1116 × 1034.8334 × 102
bk2.2936 × 1051.7267 × 1065.9245 × 106
Table 6. Natural frequencies of the first two modes for cantilever EVES beams: comparison of experimental and calculated results.
Table 6. Natural frequencies of the first two modes for cantilever EVES beams: comparison of experimental and calculated results.
BeamExperiment [35]Sisemore Model [35]Mead-Markus Model [16]Compressional Model
Frequency
(Hz)
Frequency
(Hz)
Error
(%)
Frequency
(Hz)
Error
(%)
Frequency
(Hz)
Error
(%)
First natural frequency
147.347.30.155.71847.81.0
244.545.21.556.62743.81.6
344.248.59.766.95148.69.9
443.044.02.448.41344.53.4
540.942.33.647.61741.20.73
640.744.18.555.43645.612
739.640.72.845.01441.24.0
838.139.53.644.41738.51.0
937.841.08.5523841.710
Second natural frequency
13292979.8376142999.0
2319284114233327414
3350313115285130114
42932766.0357222785.1
52842666.2364282579.5
63052934.0419372808.2
72702555.3344282584.4
82622495.0344322418.0
92862763.6388362599.4
Table 7. Damping ratios of the first two modes for cantilever EVES beams: comparison of experimental and calculated results.
Table 7. Damping ratios of the first two modes for cantilever EVES beams: comparison of experimental and calculated results.
BeamExperiment [24]Sisemore Model [24]Mead-Markus Model [11]Compressional Model
Damping RatioDamping RatioError (%)Damping RatioError (%)Damping RatioError (%)
Damping ratio corresponding to the first natural frequency
10.02230.00122000.04621100.014635
20.02670.001517000.07681900.019726
30.02150.0151420.09143300.019012
40.02140.0002110000.05761700.028835
50.02260.00238800.06631900.032644
60.02260.0203110.06631800.027120
70.02290.000375000.06631900.030834
80.02240.00287000.07022100.032545
90.02310.022340.05961600.02425
Damping ratio corresponding to the second natural frequency
10.05139 × 10−5570000.0392240.05517
20.08260.0011740000.0654210.07657
30.10840.007314000.0809250.080126
40.09140.0002460000.080970.08794
50.10550.001569000.0980180.130323
60.10270.007413000.1243240.147443
70.09770.0005190000.1274230.128731
80.09670.001853000.1205770.121826
90.10700.007214000.1707120.135026
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, Z.; Pan, J.; Yang, Z.; Wang, X.; Chu, F. Transverse Vibration of Viscoelastic Sandwich Structures: Finite Element Modeling and Experimental Study. Materials 2021, 14, 7751. https://doi.org/10.3390/ma14247751

AMA Style

Huang Z, Pan J, Yang Z, Wang X, Chu F. Transverse Vibration of Viscoelastic Sandwich Structures: Finite Element Modeling and Experimental Study. Materials. 2021; 14(24):7751. https://doi.org/10.3390/ma14247751

Chicago/Turabian Style

Huang, Zhicheng, Jinbo Pan, Ziheng Yang, Xingguo Wang, and Fulei Chu. 2021. "Transverse Vibration of Viscoelastic Sandwich Structures: Finite Element Modeling and Experimental Study" Materials 14, no. 24: 7751. https://doi.org/10.3390/ma14247751

APA Style

Huang, Z., Pan, J., Yang, Z., Wang, X., & Chu, F. (2021). Transverse Vibration of Viscoelastic Sandwich Structures: Finite Element Modeling and Experimental Study. Materials, 14(24), 7751. https://doi.org/10.3390/ma14247751

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